DOI： 10.11817/j.issn.1672-7207.2020.07.023

全长黏结GFRP抗浮锚杆荷载分布函数模型研究

白晓宇^{1, 2},张明义^{1, 2},匡政¹,王永洪^{1, 2},闫楠^{1, 2}

（1. 青岛理工大学 土木工程学院，山东 青岛，266033；

2. 山东省高等学校蓝色经济区工程建设与安全协同创新中心，山东 青岛，266033）

摘要:基于荷载传递法理论与Kelvin问题的位移解，进一步推导全长黏结GFRP(glass fiber reinforced plastics)抗浮锚杆在轴向拉拔荷载作用下轴力沿锚固深度的分布函数。为验证该理论应用于GFRP抗浮锚杆的合理性，借助植入式光纤光栅传感技术，对2根同型号GFRP抗浮锚杆进行现场拉拔破坏性试验。研究结果表明：根据荷载传递法与Kelvin位移解得到锚杆轴力与剪应力分布函数曲线形式与试验结果相近，说明该理论合理；孔口锚固体的开裂导致锚杆轴力及剪应力分布曲线试验值主要分布范围比理论值的更大；由于试验过程中岩土体位移较小，锚杆的剪应力分布曲线形式较理论值呈现“矮胖”的特点。此外，对锚杆轴力、剪应力理论分布函数曲线进行修正后的结果与试验结果吻合度显著提高。

关键词:GFRP抗浮锚杆；轴力；剪应力；荷载传递法；Kelvin位移解；拉拔试验

中图分类号:TU473      文献标志码:A            开放科学(资源服务)标识码(OSID)

文章编号:1672-7207（2020）07-1977-12

Load distribution function model of full-length bond GFRP anti-floating anchor

BAI　Xiaoyu^{1, 2}, ZHANG　Mingyi^{1, 2}, KUANG　Zheng¹, WANG　Yonghong^{1, 2}, YAN　Nan^{1, 2}

(1. College of Civil Engineering, Qingdao University of Technology, Qingdao 266033, China；

2. Cooperative Innovation Center of Engineering Construction and Safety in Shandong Blue Economic Zone, Qingdao 266033, China)

Abstract: Based on the load-transfer method and Kelvin’s displacement solution, the distribution function for the axial force of full-length grouted anchor under the pullout load was deduced. In order to prove the rationality that the present theory applies to the GFRP(glass fiber reinforced plastics) anti-floating anchor, an in-site pullout test was conducted with embedded bare fiber sensing technology, for the two same type GFRP anti-floating anchors. The results show that the curve form of the anchor axial force and shear stress, on the basis of the load-transfer method and Kelvin’s displacement solution, is similar to the test results, which indicates the rationality of the theory mentioned above. Due to the concrete cracking near the anchor hole port, the main distribution range for the experimental curve of the axial force and shear stress is deeper than that of the theoretical value. Because of a smaller displacement of the rock-soil body during the test process, the form for the distribution curve of the shear stress presents a “short and fat” property, compared with the curve form of the theoretical value. Additionally, the test-theory conformity is improved obviously after corrected the theoretical distribution curve of the axial force and shear stress.

Key words: GFRP anti-floating anchor; axial force; shear stress; load-transfer method; Kelvin’s displacement solution; pullout test

抗浮锚杆作为一种建(构)筑物基础抗浮措施，具有单点受力小、施工简单、成本低等优点，受到工程界的广泛青睐^[1-4]。由于抗浮锚杆的工作环境比普通锚杆的更为恶劣，对传统钢筋锚杆的腐蚀相当严重，大大降低了抗浮锚杆的服役寿命，新型非金属抗浮锚杆应运而生。近年来，玻璃纤维增强聚合物(glass fiber reinforced plastics，GFRP)锚杆因其抗拉性能好、质量小、抗腐蚀性强和抗电磁干扰能力优良等特点^[5-8]被越来越多的应用于复杂地质环境下建(构)筑物的抗浮工程中^[9-11]。GFRP材料作为一种新兴的建筑材料，其工程使用性备受瞩目，国内外专家学者对GFRP抗浮锚杆的材料与力学特性进行了大量试验与理论研究。LEE等^[12]对锚固在高强混凝土与普通混凝土中GFRP筋进行拉拔试验，探究了GFRP筋材在不同强度混凝土之间黏结强度的变化规律。黄生文等^[13]也进行了GFRP锚杆与不同强度水泥砂浆间拉拔试验，探究水泥砂浆强度对GFRP锚杆的锚固特性及黏结强度的影响。白晓宇等^[14]对GFRP抗浮锚杆施加稳定拉拔荷载，探究GFRP抗浮锚杆在长期荷载作用下的蠕变特性。孙晓云等^[15]提出改进的指-幂函数模型，预测了抗浮锚杆拉拔过程中的荷载-位移曲线，具有较好的拟合效果。FAVA等^[16]通过大量拉拔试验探究锚杆直径不同对GFRP锚杆与混凝土间黏结特性的影响，并且采用非线性有限元数值模拟了拉拔试验过程中GFRP锚杆与混凝土锚固界面处破裂面的发展过程，模拟结果与试验结果较为相近。以上研究成果均较好地反映了GFRP筋作为抗浮锚杆的力学性能与工作特性。对于剪应力沿锚固深度的分布模式，钟志彬等^[17]利用荷载传递法与Kelvin位移解进行预测，但轴力沿锚杆长度的分布函数尚未可知。本文作者基于荷载传递理论与Kelvin问题的位移解，在文献[17]的基础上，对全长黏结GFRP抗浮锚杆的剪应力和轴力沿锚固深度的分布函数进行更深入的研究，将2根全长黏结的GFRP抗浮锚杆现场拉拔破坏性试验结果与函数计算结果进行对比，并分析二者之间误差，据此对理论解进行修正。

1 基本理论分析

通常所说的锚杆是由锚杆杆体与锚固体共同组成，拉拔荷载通过二者组成的锚杆传递至周围岩土体中；此外，GFRP锚杆杆体的弹性模量与锚固体的比较相近，二者协同工作效果较好。鉴于上述原因，尽管GFRP锚杆杆体与锚固体为2种不同的材料，本文将二者视为整体，称之为锚杆，受力分析时共同考虑。

假设全长黏结GFRP抗浮锚杆与周围岩土体为均匀弹性体，在工作状态下将二者视为无相对滑移的整体。将锚杆视为微型抗拔桩，考虑到其荷载传递特性与抗拔桩类似，根据荷载传递理论^[18]，将锚杆分解成多个以非线性弹簧与周围岩土层相连的弹性单元，弹性单元体受力情况如图1所示。其中：P为锚杆所受到的拉拔荷载；P(x)为深度x时单元体的轴力；τ(x)为深度x时单元体所受摩阻力；dx为单元体长度；ds(x)为深度x时单元体的弹性变形量；s(x)为深度x时锚杆的弹性变形量。

图1　荷载传递理论模型

Fig. 1　Load-transfer theory model

对于锚杆单元体，由静力平衡条件得：

(1)

根据胡克定律，单元体的弹性变形量为

(2)

式中：r为锚杆半径；E_a为锚杆综合弹性模量，；E_G为锚杆杆体弹性模量；E_c为锚固体弹性模量；r_G为锚杆杆体半径。

对式(2)左右两边同时求导，代入式(1)得

(3)

可见：求解s(x)是解决问题的关键。由于锚杆受力情况与Kelvin位移解^[19]假设相似，因此本文借助Kelvin位移解求解锚杆周围岩土体的竖向位移。由于锚杆与岩土体变形协调，岩土体的竖向位移亦是锚杆的弹性变形量。

Kelvin位移解的力学模型如图2所示。假设不考虑锚杆自身的质量，将x轴与锚杆锚固方向重合，锚杆顶端即为坐标原点O。当一集中力P作用于坐标原点O时，空间内某一点A(x，y，z)的竖直方向位移为

(4)

式中：μ为岩土体的泊松比；E为岩土体的弹性模量；。

图2　Kelvin 模型示意图

Fig. 2　Schematic diagram of Kelvin model

由于杆体与x轴重合，锚杆上各单元体坐标y=z=0，且锚杆实际受拉方向与图2模型中的集中力相反，则式(4)可变为

(5)

根据图1所示锚杆模型，锚杆在深度为x处的轴力为

(6)

将深度为x处的锚杆杆体所在截面视为Kelvin模型中xOy平面，则该截面的锚杆轴力P(x)即为图2中的集中力P，因此可对式(5)和式(6)进行联立：

(7)

式中：L_a为锚杆锚固长度。

将式(7)代入式(3)，则

(8)

将式(8)展开后化简得到二阶变系数齐次常微分方程：

(9)

式中：。

假设锚杆为半无限长，则有边界条件为，，代入式(9)得后解得

(10)

在锚杆顶端，亦有边界条件：，代入式(10)得

(11)

为提高锚固效果，通常抗浮锚杆的锚固长度较长(不小于3.0 m)，故式(11)中计算值很小，可以忽略不计，则式(11)简化为

(12)

将式(12)代入式(10)中得到的GFRP抗浮锚杆剪应力沿深度分布的表达式为

(13)

最后，将式(13)代入式(6)中即可得到GFRP抗浮锚杆轴力沿深度分布的函数为

(14)

2 现场试验

2.1　试验方案

为探究GFRP抗浮锚杆在实际工作情况下剪应力与轴力的分布情况，对编号为G28-1和G28-2的2根全长黏结GFRP抗浮锚杆进行破坏性拉拔试验。为提高试验的准确性，本试验使用植入式光纤光栅传感器^[20]，对GFRP抗浮锚杆杆体上剪应力与轴力进行测量。传感器植入的具体制作过程为：首先将微型FBG(fiber Bragg grating)传感器串沿杆体长度方向放置于1根锚杆所需的玻璃纤维束中心，然后将其与玻璃纤维束一起浸泡于基体材料(环氧树脂)中浇筑成型，再经过挤拉成型、缠绕螺纹后定型固化等工序制备得到本次试验使用的植入FBG传感器的GFRP抗浮锚杆。FBG传感器的安装布置方式如图3所示。

图3　FBG传感器布置示意图

Fig. 3　Schematic diagram installation of FBG sensors

本试验采用低温敏型FBG传感器，温度变化对本次试验影响不大，在此忽略。在拉拔过程中，随着锚杆受到向上的拉拔荷载而发生伸长形变，FBG传感器的反射或透射峰波长也随之发生改变，使用光纤光栅解调仪对波长变化进行记录后，利用式(15)进行计算后即可得到杆体的应变^[21]：

(15)

式中：△ε为GFRP抗浮锚杆杆体应变；△λ_B为FBG传感器中心波长变化量；K_ε为FBG传感器应变灵敏度系数。

根据式(16)可进一步计算得到杆体轴力F为

(16)

由于前文推导过程中假设GFRP锚杆杆体与锚固体为整体，故本文将杆体轴力视为锚杆轴力，根据锚杆轴力即利用式(17)对应关系计算锚杆与岩土体之间的剪应力：

(17)

式中：τ为锚杆剪应力，MPa；△L为相邻的FBG传感器间距离，mm。

2.2　试验场地及锚杆参数

在青岛市崂山区某建筑工地进行拉拔试验，试验场地地质组成成分主要为中风化花岗岩，岩体节理裂隙发育，主要呈块状结构，具体力学参数如表1所示。

本试验GFRP锚杆杆体基体材料为环氧树脂，其体积分数接近25%，玻璃纤维体积分数约为75%，总长8 m，其中地下锚固段为6.5 m，地面自由加载段1.5 m，其他基本力学参数如表2所示。

2.3　试验流程

1) 钻孔。为充分发挥杆体材料抗拉性能，按照施工要求，钻机成孔过程中，锚孔的深度大于杆体锚固段长度0.5 m，即锚孔深度为7 m，孔径为120 mm，此外，为避免群锚效应^[22-23]对试验结果产生影响，2根试验锚杆间距定为3 m。需要说明的是，由于锚杆受拉，锚杆杆体底端以下的纯锚固砂浆部分对于抗拉作用不大，因此，进行计算时锚杆的锚固长度不考虑该部分，即锚固长度取6.5 m。

2) 灌浆及养护。锚孔开挖完毕后，将试验锚杆竖直伸入孔内指定深度，灌入M30水泥砂浆，固化后其弹性模量取20 GPa，锚杆综合弹性模量为21.143 GPa。注浆完毕后的锚杆如图4所示。

表1　试验场地力学参数

Table 1　Mechanical parameters of test site

表2　试验GFRP锚杆基本参数

Table 2　Basic parameters of test GFRP anchors

图4　锚杆注浆

Fig. 4　Anchor grouting

3) 加载装置安装。待灌浆体养护28 d后，在地表以上锚杆段安装如图5所示的试验加载装置，该试验装置的荷载传递机理为：将锚具焊接在钢套管外表面，千斤顶施加的拉拔荷载经垫板传递至焊接锚具，再由焊接锚具通过钢套管和其内部的结构胶传递至锚杆杆体，锚杆杆体最终将该拉拔力经由锚固体传递至周围岩土体。

4) 试验加载制度。试验采用逐级加载法，以0.2 kN/s左右的速率、每级40 kN逐级匀速加载，相邻2级荷载加载时间间隔15 min，直至锚杆发生破坏。每级荷载加载完成后，迅速对百分表显示的相对位移进行测读，此后每间隔5 min测读1次，直至15 min 后进行下级荷载加载。根据GB 50007—2011 “建筑地基基础设计规范”确定GFRP抗浮锚杆的极限抗拔承载力及最大滑移量。

图5　试验加载装置示意图

Fig. 5　Schematic diagram of loading device

3 试验结果及分析

3.1　试验结果

2根锚杆锚固体均发生开裂，锚固体顶面出现破碎。各试验锚杆的破坏荷载及最终破坏形式如表3所示。

表3　锚杆破坏情况

Table 3　Failure condition of test anchors

由表3可知：试验过程中锚固体与锚杆杆体之间的黏结作用力大于锚杆杆体本身的材料强度，因此锚杆材料性能得到充分发挥，产生了工程应用中的理想破坏形式—杆体被拉断；此外，GFRP锚杆可以提供较大的抗拔力，证明其应用于建筑物抗浮结构的可行性。

根据式(16)和(17)计算得到的全长黏结的GFRP锚杆轴力与剪应力沿其锚固深度的分布曲线，并且将其与基于上述理论推导结果进行对比，对比结果如图6和图7所示。

图6　不同荷载下全长黏结锚杆轴力沿锚固深度分布

Fig. 6　Distribution of full length bond anchor axial force along anchorage depth under different loads

图7　不同荷载下全长黏结锚杆剪应力沿锚固深度分布图

Fig. 7　Distribution of full length bond anchor shear stress along anchorage depth under different loads

3.2　误差描述及分析

由图6和图7可知：试验锚杆G28-1和G28-2试验数据的分布情况较为一致，说明试验较为成功，未出现较大偏差。总的来说，在不同荷载作用下，根据计算得到的理论曲线与根据FBG传感器测得锚杆轴力、剪应力沿锚固深度变化的曲线趋势基本一致，亦证实了基于Kelvin位移解及荷载传递法计算得到的GFRP锚杆轴力和剪应力分布函数模型是科学有效的且具有实用意义。但是，无论是轴力还是剪应力，其理论曲线与试验结果之间存在明显误差，主要包括以下2点：

1) 试验锚杆轴力、剪应力试验曲线的主要分布范围较理论值更广，即锚杆开始有效抵抗拉拔荷载的深度比理想条件更深；且拉拔荷载越大，试验曲线分布范围越广。

2) 试验锚杆剪应力曲线主要分布范围比理论曲线大且试验曲线的高度较高，即理论剪应力曲线线型“高瘦”，试验曲线线型“矮胖”。

为了给实际应用提供相关依据，分析图6和图7中轴力及剪应力沿锚固深度变化理论曲线与实测曲线出现误差的原因。

出现误差1) 的原因主要有：本文理论推导中将锚杆杆体与锚固体视为整体，二者在拉拔荷载作用下共同变形，无脱开现象，外锚固段锚杆杆体始终与锚固体充分作用，作为一个整体将自由段所受拉拔荷载向周围岩土层传递，造成锚杆锚固孔口附近剪应力迅速上升，而轴力则随之迅速下降，在距离孔口不远处即将至较低水平。但在试验过程中可以明显看到锚杆锚固孔口处在拉拔过程中出现开裂现象(图8)，二者之间的相互作用效率降低，对拉拔荷载的抵消作用严重降低，导致该部分锚杆杆体与锚固体组成的整体不再满足理论假设的条件；此时，更深处的锚杆与锚固体之间尚未产生脱开现象，二者组成整体的力学性能相对完整，可继续视为理论假设中的锚杆整体，因此理论推导中的锚杆上表面位置将在实际试验中逐渐下移，从更深部位开始有效抵消锚杆所受拉拔力。

拉拔荷载通过继续向下转移，由埋深更大处的锚杆及锚固体向周围岩土层传递，因此，由FBG传感器实测的GFRP锚杆轴力、剪应力沿埋深分布曲线的分布区域比理论解分布区域要深。

此外，随着施加荷载的增加锚固孔口处的开裂更为明显，表明随着荷载增加，锚固体与锚杆杆体间无效锚固范围增大，理论假设中的锚固体顶面继续下移，导致拉拔荷载持续下移，只能由锚固更深的锚固体向周围岩土层传递，GFRP抗浮锚杆杆体的轴力和剪应力的零值点深度也随之下移。

图8　锚固体开裂

Fig. 8　Concrete cracking

出现误差2) 的原因主要有：本文假设中，锚杆与岩土体协调变形，二者弹性位移相同，但实际试验中，荷载作用下岩土体位移很小，远低于锚杆位移，因此，可视为只有锚杆产生变形，为限制岩土体位移，由于加载支墩的存在，将理想条件下锚杆受力情况转换为图9所示模型。

图9　锚杆等效受力模型

Fig. 9　Equivalent force model of anchor

从图9可以看出：为限制岩土体位移，锚杆受拉时，在岩土体上施加额外的向下压力T_e，额外压力与锚杆所受剪应力的合力与拉拔荷载相等，由式(18)表示：

(18)

其中：d为锚杆杆体直径。根据锚杆系统静力平衡条件可知，对于式(18)，在拉拔荷载P一定的情况下，压力T_e的出现必然导致剪应力τ的降低。

此外，函数推导时，假设锚杆与岩土体是均匀材料，则锚杆及岩土体二者接触面位置在拉拔荷载作用下相互作用，将拉拔荷载由锚杆传递至岩土体中。然而，由于在实验过程中施工及养护工艺等缺陷，锚固砂浆并非均匀分布于锚孔内，锚杆与锚固岩土体并非始终牢固地黏结在一起，当拉拔荷载传递至这些截面时无法通过锚杆-岩土体界面横向传递，只能继续向下传递，导致需要更长的锚杆锚固范围来抵消拉拔荷载，因此，试验剪应力曲线主要分布范围比理论曲线的大。

3.3　讨论

令式(6)右边的深度x取值为0，则式(6)变为

(19)

由式(19)可知：锚杆所受拉拔荷载P和x轴与剪应力分布曲线所围图形的面积成正相关，在锚杆半径一定的情况下，相同拉拔荷载P对应相同的剪应力曲线与x轴所围的面积。因此，当试验剪应力峰值因岩土体未产生明显位移而降低时，剪应力曲线分布范围必然变大，从而使得该曲线变得“矮胖”。由图7也可以看出：实际剪应力曲线与x轴所围图形的面积与理论曲线相似，间接证明了本次试验测量结果的准确性。

4 锚杆荷载分布函数修正

4.1　轴力修正

对于锚杆轴力分布函数的修正，观察图6发现，可以考虑在式(14)中的自变量x前引入修正系数，从而使轴力分布曲线的主要分布范围增大，因此，从试验和理论轴力分布曲线接近x轴的位置入手，对理论轴力分布曲线进行修正。GFRP锚杆试验轴力分布函数接近x轴(2根试验锚杆的平均轴力P＜1 kN)的位置如表4所示，锚杆理论轴力分布函数与x轴的距离为0.71 m。

表4　试验锚杆轴力曲线与x轴距离

Table 4　Distance of test anchor axial force curves to x-axis

由表4可知：锚杆理论轴力分布曲线接近x轴的位置为试验值的0.23~0.40倍，因此，取其中间值0.315作为理论轴力分布函数的修正系数，则式(14)修正为

(20)

根据式(20)得到的GFRP锚杆轴力分布曲线如图10所示。

图10表明，对锚杆轴力分布函数的修正具有较好的效果，使得按式(20)求得的轴力分布函数与试验数据比较接近。

4.2　剪应力修正

对于剪应力，需解决理论曲线较“高”和较“瘦”的问题，其中理论曲线较“瘦”问题的解决方案与4.1节轴力修正方法相同，即在式(13)自变量x前引入修正系数0.315，使得曲线主要分布范围增大。

对于理论曲线较“高”的问题，拟从剪应力峰值点的理论与实际误差入手进行修正，本试验中，不同荷载作用下锚杆的剪应力理论峰值与实际峰值之间的误差如表5所示，其中，试验剪应力峰值取2根锚杆试验剪应力峰值的平均值。

由表5可知：GFRP锚杆剪应力理论峰值与试验峰值的比值在2.3~3.2之间，且随拉拔荷载的增大，该比值有逐渐增大的趋势。

对于式(13)所示理论剪应力分布曲线，需引入适当的折减系数，用以降低其与实际试验结果在取值上的误差，根据表5，本文取折减系数为1/2.75，则式(13)变为

(21)

根据式(21)得到的GFRP锚杆剪应力分布曲线如图11所示。

图10　函数修正后不同荷载下锚杆轴力沿深度分布

Fig. 10　Depth distribution of anchor axial force along different loads after correction function

图11　函数修正后不同荷载下锚杆剪应力沿深度分布

Fig. 11　Depth distribution of bolt shear stress along different loads after correction function

由图11可知：修正后的理论剪应力分布函数在取值上与实际试验结果吻合度大大提升。

4.3　轴力和剪应力分布函数修正系数适用性验证

白晓宇等^[24]利用三重光纤光栅传感器串同步测得GFRP抗浮锚杆杆体和锚固体-岩土体界面的轴力及剪应力分布形式。本文基于文献[24]中锚杆G6-01对上述剪应力和轴力修正系数进行验证。根据文献[24]相关试验计算参数，按照式(20)和式(21)计算得到的GFRP锚杆剪应力和轴力分布曲线与试验结果的对比分别如图12和图13所示。

表5　锚杆剪应力峰值误差

Table 5　Error for maximum value of anchor shear stress

图12　锚杆G6-01^[24]剪应力分布情况

Fig. 12　Distribution for shear stress of anchor G6-01^[24]

图13　锚杆G6-01^[24]轴力分布情况

Fig. 13　Distribution for axial force of the anchor G6-01^[24]

由图12和图13可以看出：按照式(20)和式(21)计算得到的全长黏结式GFRP抗浮锚杆的轴力和剪应力分布函数曲线与文献[24]中的试验结果吻合度较高，说明式(20)和式(21)中的修正系数取值适用性较高。结合本文试验结果可知：当GFRP抗浮锚杆锚固于硬度较高的中风化岩土层中时，利用上述计算方法预测锚杆在拉拔荷载作用下的荷载沿锚固深度的分布情况较为准确。

5 结论

1) 将荷载传递法及半无限空间内Kelvin位移解的相关理论融合进锚杆力学模型中，推导了GFRP抗浮锚杆在受荷状态下的轴力与剪应力沿锚固深度的分布函数。

2) 借助植入式裸光纤光栅传感技术，对2根全长黏结式GFRP抗浮锚杆进行现场拉拔破坏性试验，依据试验结果描绘的锚杆轴力、剪应力分布曲线与理论推导结果形式相近，表明本文理论推导方法合理。

3) 试验过程中孔口附近锚固体的开裂导致GFRP锚杆轴力、剪应力试验分布曲线的主要分布范围较理论值更大，即锚杆开始有效抵抗拉拔荷载的深度大于理想条件的深度。

4) 理论假设中锚杆与岩土体协调变形，但试验过程中岩土体的位移非常小，因此，GFRP锚杆实际剪应力分布曲线较“矮胖”。

5) 修正了理想条件下全长黏结GFRP抗浮锚杆的轴力和剪应力沿锚固深度分布函数，修正后的计算结果与本文试验结果吻合度显著提高。

参考文献：

[1] 贾金青, 宋二祥. 滨海大型地下工程抗浮锚杆的设计与试验研究[J]. 岩土工程学报, 2002, 24(6): 769-771.

JIA Jinqing, SONG Erxiang. The design and test on anti-floating anchorage of large substructure in coastal region[J]. Chinese Journal of Geotechnical Engineering, 2002, 24(6): 769-771.

[2] KOU Hailei, GUO Wei, ZHANG Mingyi. Pullout performance of GFRP anti-floating anchor in weathered soil[J]. Tunnelling and Underground Space Technology, 2015, 49: 408-416.

[3] ASHRAFI H, BAZLI M, VATANI OSKOUEI A, et al. Effect of sequential exposure to UV radiation and water vapor condensation and extreme temperatures on the mechanical properties of GFRP bars[J]. Journal of Composites for Construction, 2018, 22(1): 04017047.

[4] 匡政, 白晓宇, 张明义, 等. 考虑锚固体不均匀及杆体脱黏效应的GFRP抗浮锚杆杆体荷载分布函数[J]. 岩石力学与工程学报, 2019, 38(6): 1158-1171.

KUANG Zheng, BAI Xiaoyu, ZHANG Mingyi, et al. Load distribution function of GFRP anti-floating anchors considering the anchorage body unevenness and the anchor debonding effect[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(6): 1158-1171.

[5] MALVAR L J. Tensile and bond properties of GFRP reinforcing bars[J]. ACI Materials Journal, 1995, 92(3): 276-285.

[6] American Concrete Institute. Guide for the design and construction of concrete reinforced with FRP bars[R]. Farmington Hills, Michigan: American Concrete Institute, 2006.

[7] 刘汉东, 于新政, 李国维. GFRP锚杆拉伸力学性能试验研究[J]. 岩石力学与工程学报, 2005, 24(20): 3719-3723.

LIU Handong, YU Xinzheng, LI Guowei. Experimental study on tensile mechanical properties of glass fiber reinforced plastic rebar[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(20): 3719-3723.

[8] 白晓宇, 张明义, 寇海磊. 基于裸光纤光栅传感技术GFRP抗浮锚杆荷载传递机制的原位试验研究[J]. 工程力学, 2015, 32(8): 172-181.

BAI Xiaoyu, ZHANG Mingyi, KOU Hailei. Field experimental study of load transfer mechanism of GFRP anti-floating anchors based on embedded bare fiber Bragg grating sensing technology[J]. Engineering Mechanics, 2015, 32(8): 172-181.

[9] 杨蒙, 张彦忠, 谭跃虎, 等. GFRP锚杆在抗浮工程中的应用研究[J]. 四川建筑科学研究, 2013, 39(3): 112-115.

YANG Meng, ZHANG Yanzhong, TAN Yuehu, et al. Research on the application of GFRP anchor in anti-floating engineering[J]. Sichuan Building Science, 2013, 39(3): 112-115.

[10] PENDHARI S S, KANT T, DESAI Y M. Application of polymer composites in civil construction: a general review[J]. Composite Structures, 2008, 84(2): 114-124.

[11] 白晓宇, 张明义, 匡政, 等. 风化岩地基中玻璃纤维增强聚合物抗浮锚杆体系受力特性数值模拟[J]. 中南大学学报(自然科学版), 2019, 50(6): 1466-1473.

BAI Xiaoyu, ZHANG Mingyi, KUANG Zheng, et al. Numerical simulation on bearing behavior of glass fiber reinforcement polymer anti-floating anchor system in weathered rock foundation[J]. Journal of Central South University(Science and Technology), 2019, 50(6): 1466-1473.

[12] LEE J Y, KIM T Y, KIM T J, et al. Interfacial bond strength of glass fiber reinforced polymer bars in high-strength concrete[J]. Composites Part B: Engineering, 2008, 39(2): 258-270.

[13] 黄生文, 邱贤辉, 罗文兴. GFRP锚杆锚固特性研究[J]. 长沙理工大学学报(自然科学版), 2009, 6(3): 33-39.

HUANG Shengwen, QIU Xianhui, LUO Wenxing. The bond feature researching of the GFRP (glass fiber reinforced polymer) bolt[J]. Journal of Changsha University of Science & Technology(Natural Science), 2009, 6(3): 33-39.

[14] 白晓宇, 张明义, 张舜泉. 全长黏结螺纹玻璃纤维增强聚合物抗浮锚杆蠕变试验研究[J]. 岩石力学与工程学报, 2015, 34(4): 804-813.

BAI Xiaoyu, ZHANG Mingyi, ZHANG Shunquan. Creep testing on anti-floating anchors of full-length bonding thread glass fiber reinforced polymer(GFRP)[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(4): 804-813.

[15] 孙晓云, 张涛, 王明明, 等. 基于改进指-幂混合函数模型的锚杆承载力预测方法研究[J]. 岩石力学与工程学报, 2015, 34(8): 1641-1649.

SUN Xiaoyun, ZHANG Tao, WANG Mingming, et al. A revised model for predicting the bearing capacity of rock bolts based on mixed exponential and power unction[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(8): 1641-1649.

[16] FAVA G, CARVELLI V, PISANI M A. Remarks on bond of GFRP rebars and concrete[J]. Composites Part B: Engineering, 2016, 93: 210-220.

[17] 钟志彬, 吕蕾, 邓荣贵. 考虑轴力分布的全长粘结锚杆受力分析[J]. 防灾减灾工程学报, 2013, 33(3): 311-315.

ZHONG Zhibin, LU Lei, DENG Ronggui. Mechanical analysis of wholly grouted bolt considering axial distribution along the bars[J]. Journal of Disaster Prevention and Mitigation Engineering, 2013, 33(3): 311-315.

[18] SEED H B, REESE L C. The action of soft clay along friction piles[J]. American Society of Civil Engineers Transactions, 1955, 122(1): 731-754.

[19] 王光钦. 弹性力学[M]．北京: 中国铁道出版社, 2008: 178-179.

WANG Guangqin. Elastic mechanics[M]. Beijing: China Railway Publishing House, 2008: 178-179.

[20] 白晓宇, 张明义, 匡政, 等. 光纤光栅传感技术在GFRP抗浮锚杆现场拉拔试验中的应用[J]. 岩土力学, 2018, 39(10): 3891-3899.

BAI Xiaoyu, ZHANG Mingyi, KUANG Zheng, et al. Application of fiber grating sensing technology in pull-out test on glass fiber reinforced polymer anti-floating anchor[J]. Rock and Soil Mechanics, 2018, 39(10): 3891-3899.

[21] 邱松, 顾浩, 曹进捷, 等. FBG传感技术在混凝土预制桩水平载荷试验中的应用[J]. 岩土工程学报, 2011, 33(S2): 105-107.

QIU Song, GU Hao, CAO Jinjie, et al. Application of FBG sensing technology in lateral load tests on precast concrete piles[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(S2): 105-107.

[22] GHALY A, HANNA A. Model investigation of the performance of single anchors and groups of anchors[J]. Canadian Geotechnical Journal, 1994, 31(2): 273-284.

[23] 张昕, 乐金朝, 刘汉东. 砂土中群锚锚周土体变形特性模型试验研究[J]. 岩土力学, 2016, 37(S1): 240-248.

ZHANG Xin, YUE Jinchao, LIU Handong. Experimental study of soil deformation around group anchors in sand[J]. Rock and Soil Mechanics, 2016, 37(S1): 240-248.

[24] 白晓宇, 张明义, 朱磊, 等. 全长黏结GFRP抗浮锚杆界面剪切特性试验研究[J]. 岩石力学与工程学报, 2018, 37(6): 1407-1418.

BAI Xiaoyu, ZHANG Mingyi, ZHU Lei, et al. Experimental study on shear characteristics of interface of full-bonding glass fiber reinforced polymer anti-floating anchors[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(6): 1407-1418.

（编辑 赵俊）

收稿日期： 2019 -07 -24; 修回日期： 2019 -09 -30

基金项目(Foundation item)：国家自然科学基金资助项目(51708316，51778312，51809146)；中国博士后科学基金资助项目(2018M632641)；山东省重点研发计划项目(2017GSF16107，2018GSF117008)；山东省自然科学基金资助项目(ZR2016EEQ08，ZR2017PEE006)；山东省高等学校科技计划项目(J16LG02)；青岛市博士后应用研究项目(2018101) (Projects(51708316, 51778312, 51809146) supported by the National Natural Science Foundation of China; Project(2018M632641) supported by China Postdoctoral Science Foundation; Projects(2017GSF16107, 2018GSF117008) supported by the Key Research and Development Program of Shandong Province; Projects(ZR2016EEQ08, ZR2017PEE006) supported by the Natural Science Foundation of Shandong Province; Project(J16LG02) supported by the Higher Educational Science and Technology Program of Shandong Province; Project(2018101) supported by the Applied Research Programs of Qingdao Postdoctoral)

通信作者：张明义，博士，教授，从事土力学及地基基础研究；E-mail: zmy58@163.com