网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于双界面滑移耦合的水泥土锚杆荷载传递模型  PDF

  • 张根宝 1,2
  • 陈昌富 3
  • 徐长节 4
  • 陈海军 1
1. 湖南城市学院,土木工程学院,湖南 益阳 413000; 2. 湖南城市学院,城市地下基础设施结构安全与防灾湖南省工程研究中心,湖南 益阳 413000; 3. 湖南大学 土木工程学院,长沙 410082; 4. 江西省地下空间技术开发工程 研究中心,南昌 33001

中图分类号: TU476

最近更新:2022-11-09

DOI:10.11835/j.issn.2096-6717.2021.244

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

水泥土锚杆的承载过程伴随筋体—黏结体界面和黏结体—岩土体界面的黏结强度调动,界面剪应力的径向传播机制受到不同界面附近材料的应力和变形条件影响,对水泥土锚杆进行荷载传递分析需考虑双界面的剪切变形耦合。结合锚杆受力变形分析常用的荷载传递法和剪切位移法,将二者分别用于界面剪应力在锚杆轴向和径向引起的变形分析计算,考虑锚杆双界面的剪应力与剪切变形耦合,建立基于界面特性测试的水泥土锚杆双界面滑移耦合荷载传递分析模型。通过水泥土的材料性质试验及单元体尺度和模型试验尺度的水泥土锚杆拉拔试验,获得了模型的计算参数取值,并验证了模型对锚杆拉拔响应的预测能力。

锚杆锚固体系广泛应用于岩土工程中,随着材料组分、现场工况的迭代更新,各种新型锚杆逐渐涌现,相应的荷载传递机理和设计计算理论也随之成为研究热[

1]。锚杆材料组分中,筋体材料已从传统的钢筋和钢绞线等发展到抗腐蚀性能更优的纤维增强树脂材[2]和应用于遗址保护的竹木材[3]等,黏结材料也从传统的水泥砂浆拓展到树[4]以及水泥[5]等。现场工况方面,地基层状特性、地基土弱化特性、岩体离层作用、拉剪共同作用等复杂工况得到研究者较多关注,并有针对性地提出了具体的锚杆承载特性分析模[6-9]。值得说明的是,虽然锚杆的组分材料和服役工况得到了进一步发展,但其承载过程分析框架变化不大,基本沿着锚固界面黏结—滑移和径向剪应力传播两条路径分析锚固体的轴力和轴向位移等响应,前者常被统称为荷载传递[10-11],后者则被称作剪切位移[12-13]。近年来,也出现了基于小孔扩张理论的锚杆承载过程分析,但尚停留在设计锚固力评估和锚固界面强度预测等极限状态分[14]。此外,光纤光栅等测试技术的发展也促进了基于现场监测数据的锚杆承载过程的实时分析,即在受力变形分析模型中所采用的锚固界面模型可以根据锚杆轴力等监测数据进行动态调[15]。综合来看,荷载传递法和剪切位移法是基于锚杆材料和界面力学模型的解析方法,相比基于小孔扩张和光纤光栅监测的分析方法,可以实现锚杆承载全过程的力学响应预测,从而在各参数物理意义明确的条件下深入揭示锚杆荷载传递机理。

虽然锚杆承载涉及黏结体和筋体(T-B)以及黏结体与岩土体(B-M)两个界面,如图1所示,但锚杆锚固体系的设计计算一般基于最不利工况,故只需关注最薄弱界面及其相关的材料变[

16]。通常情况下,B-M界面强度弱于T-B界面,锚杆承载过程分析框架多为基于B-M界面的单界面分析,并假设T-B界面不发生相对滑[17-18]。在确定了最薄弱界面后,需要进一步确定该界面两侧的材料变形才能得到界面的荷载响应,如对土层灌浆锚杆进行受力变形分析时,一般会考虑筋体与黏结体作为锚固体的轴向拉伸变形及土体的剪切变[6,8];岩层灌浆锚杆分析则考虑注浆体的剪切变形和岩体的压剪变[13]。然而,对于水泥土锚杆,其黏结体为水泥搅拌土,强度变形特性要逊于水泥砂浆,应用地层则多为软土层,强度特性也劣于一般黏土层,其承载过程同时伴随T-B界面和B-M界面的剪切变形和界面强度调[5]。因此,对水泥土锚杆进行受力变形分析时,需考虑T-B和B-M两个界面共同剪切滑移下的协同承载,并着重分析黏结体的剪切变形和地层的压剪变形。需要指出的是,在剪应力传播过程中,T-B-M三者在界面附近会产生相对位移,相对位移大小决定了界面剪应力的大小,同时,界面剪应力又决定了界面附近材料的剪切变形大小并进一步影响界面相对位移的大小,因此,界面剪切强度的发挥过程涉及界面剪应力与界面相对位移之间的相互耦[19]

图1  锚杆双界面组成示意图

Fig. 1  Schematic of dual-interface composition of anchor

(b) 横截面

笔者综合荷载传递法和剪切位移法在锚杆承载分析中的优势,提出一种可以考虑T-B和B-M双界面剪切—滑移耦合的锚杆受力变形分析方法,并通过开展T-B和B-M界面黏结—滑移特性测试及水泥土中加筋体拉拔模型试验,对该方法中荷载传递分析的有效性和拉拔响应预测的准确性进行验证。

1 基本假设

1.1 边界条件与变形模式

由于水泥土锚杆的破坏一般由界面失效引起,即界面强度一般都小于材料抗拉强度和抗剪强[

5]。为了简化分析,进行锚杆受力变形分析时,假设加筋体一直处于线弹性拉伸阶段,水泥土和土体则处于线弹性剪切阶段,且不考虑剪切模量在径向和轴向的分布变化,即分别采用弹性模量Et以及剪切模量GbGm刻画钢筋及水泥土和土体的变形特性。

考虑到工程锚杆的自由段地层可能发生轴向变形,且锚杆现场拉拔试验多将反力施加于锚固体直径以外的地层表面,在受力变形分析模型中,采用均布应力边界条件来近似模拟锚杆自由段地层和现场拉拔试验的反力对锚固段地层表面的约束状态(图1)。

锚杆锚固体在双界面承载下的变形模式如图2所示。其主要变形特征为锚杆单元横截面上的轴向位移从剪应力影响边界到筋体轴线沿径向逐渐增大,并在两个界面上都发生相对位移。根据剪切位移法在桩基沉降计算应用中的研究结[

12-13],界面剪应力的影响半径可取黏结体直径的5~10倍。

图2  锚杆单元的变形模式示意图

Fig. 2  Schematic diagram of deformation mode of elementary anchor

1.2 界面剪应力的径向分布

筋体在拉拔力和T-B界面剪应力的作用下发生拉伸变形,受力分析单元如图3(a)所示。黏结体在T-B界面和B-M界面剪应力作用下发生剪切变形,受力分析单元如图3(b)所示,双界面间剪应力径向分布为传统剪滞模型中的计算解[

12],表达式为

τbr=τtrtr (1)

式中:τbr为黏结体内剪应力的径向分布;τt为T-B界面剪应力;rt为筋体半径。

图3  锚杆受力分析单元

Fig. 3  Force analysis elements of anchor

(a) 筋体 (b) 黏结体 (c) 地层基体

在B-M界面剪应力影响半径内,地层基体主要发生剪切变形,但由于拉拔端地层表面受到压应力作用,地层剪应力径向分布需考虑正应力轴向分布的影响,参考岩石锚杆变形分析中采用的改进剪滞模[

13],其分布函数可表示为

τmr=τbAr1r-Brr (2a)

式中:τmr为地层基体内剪应力的径向分布;τb为B-M界面剪应力;ArBr为系数,可用黏结体半径和剪应力影响半径表示为

Ar=1+(rt+tb)2(rt+tb+bm)2rt+tb (2b)
Br=1(rt+tb+bm)2rt+tb (2c)

1.3 双界面的界面剪切模型

笔者曾进行过不同含水量下砂浆—土体界面及不同配比和龄期下水泥土—加筋体界面的承载特性测试,发现在土体含水量较高时,砂浆—土体界面特性可用理想弹塑性界面剪切模型刻[

20],加筋体—水泥土界面特性则可简化为三折线的黏结—滑移模[21-22]。考虑到水泥土与水泥砂浆在材料配比及凝结硬化过程等方面具有较强的相似性,且广泛应用水泥土锚杆的软土地层一般具有较高的含水量,故水泥土锚杆B-M界面剪切行为采用理想弹塑性模型,T-B界面剪切行为则采用三折线模型刻画,如图4所示。

图4  水泥土锚杆双界面剪切模型

Fig. 4  Dual interface bond-slip models of soil mixing anchor

(a) B-M界面 (b) T-B界面

2 公式推导

由于拉拔力直接作用于锚杆筋体上,在受力分析时以筋体为主要对象,列出图3(a)中筋体单元拉伸应力—应变关系及静力平衡方程,可以得到

Px=πrt2Etut(x)x (3)
Pxx=2πrtτtx (4)

式中:Px为筋体在坐标x处的轴力;ut(x)为筋体在坐标x处的位移;τtx为坐标x处的T-B界面剪应力;rt为筋体半径;Et为筋体弹性模量。

联立式(3)式(4),可得

2ut(x)x2=2rtEtτtx (5)

为方便公式推导,统一采用切线法来表示T-B界面和B-M界面的界面剪切模型,其示意图如图5所示,相应的表达式为

τtx=ktδtx+ct (6a)
τmx=kmδmx+cm (6b)

式中:τmx为坐标x处B-M界面的剪应力;δtxδmx分别为坐标x处T-B界面、B-M界面的剪切位移;ktct以及kmcm分别为T-B和B-M界面模型上对应δtxδmx的切线斜率以及该切线在界面剪应力轴上的截距。

图5  界面模型的切线法表示示意图

Fig. 5  Schematic diagram of tangential form of interface model

由于T-B界面和B-M界面剪应力属于黏结体剪应力分布的边界条件,根据式(1)可得

τmx=τtxrtrt+tb (6c)

假设T-B界面和B-M界面黏结—滑移模型的函数表达式分别为

τt(x)=f(δt(x)) (7a)
τm(x)=g(δm(x)) (7b)

锚杆上任一点xi对应的T-B界面模型切线参数kticti及B-M界面模型切线参数kmicmi可分别表示为

kti=f '(δt(xi)) (8a)
kmi=g'(δm(xi)) (8b)
cti=τtxi-δt(xi)f '(δt(xi)) (9a)
cmi=τmxi-δm(xi)g'(δm(xi)) (9b)

根据图2中的变形模式,T-B界面剪切位移δtx)和B-M界面剪切位移δmx)可分别表示为

δtx=utx-ubtx (10a)
δmx=ubmx-umx (10b)

式中:ubtxubmx分别为坐标x处黏结体在T-B界面和B-M界面的轴向位移;umx为坐标x处地层基体在B-M界面的轴向位移。

黏结体和地层基体的轴向位移来源于剪切变形在径向的累积,不同半径处的轴向位移差可以采用剪滞模型计算,因而黏结体在T-B界面和B-M界面的位移差,及地层基体在B-M界面剪应力影响半径内的累积位移可分别表示为

ubtx-ubmx=rtrt+tbτb(r)Gbdr (11a)
umx=rt+tbrt+tb+bmτm(r)Gmdr (11b)

式中:GbGm分别为黏结体和地层基体的剪切模量;tb为黏结体与筋体的半径差;bm为B-M界面剪应力影响半径与黏结体的半径差。

联立式(1)式(11a),可得

ubtx-ubmx=Crτt (12a)

式中:Cr为与尺寸参数和黏结体剪切模量相关的系数,表达式为

Cr=rtGblnrt+tbrt (12b)

联系式(2)和(11b),可得

umx=Drτt (13a)

式中:Dr为与尺寸参数和地层基体剪切模量相关的系数,表达式为

Dr=rtGm1+rt+tb2rt+tb+bm2lnrt+tb+bmrt+tb-0.51-rt+tb2rt+tb+bm2 (13b)

联立式(10)~式(13),可得

δtx=utx-δmx-(Cr+Dr)τt (14)

式(6)代入式(14),可以得到锚杆T-B界面剪应力与筋体位移之间的关系式

τtx=Kuutx+Cu (15a)

需要说明的是,由于τt - ut关系是进行荷载传递法t-z分析的关键,此处沿用荷载传递法中的表达形式,其中,KuCu分别为

Ku=ktkmkm+ktrtrt+tb+ktkmCr+ktkmDr (15b)
Cu=ctkm+cmktkm+ktrtrt+tb+ktkmCr+ktkmDr (15c)

式(15)代入式(5)可得锚杆双界面受力变形分析的控制微分方程。

2ut(x)x2=Auutx+Bu (16a)

式中:AuBu为待定系数,可分别表示为

Au=2rtEtKu (16b)
Bu=2rtEtCu (16c)

需要注意的是,AuBu的表达式中包含T-B界面和B-M界面的界面模型参数ktctkmcm,这4个参数需要根据界面相对位移δtxδmx来确定,而界面相对位移的确定则需要筋体位移utx。因此,虽然式(16)在形式上为非齐次线性二阶偏微分显性方程,但本质上为高度非线性的二阶偏微分隐性方程,其解析或半解析求解涉及复杂的公式推导,计算难度较大。考虑到微分方程待定系数的确定涉及迭代计算,可以采用数值计算方法对方程进行求解。已有研究者将有限差分方法引入岩石锚杆的受力变形分析计算中,计算效果较[

10-11],故沿用该方法对方程式(16)进行有限差分数值求解。

3 模型解算

数值计算的基础是对锚杆的离散化处理,如图6所示,锚杆黏结段被n个结点均匀分成n-1个单元微段,每个单元长度为La/n-1,从自由端到拉拔端对结点进行编号,相应的第i个结点的坐标可表示为

xi=i-1n-1La (17)

图6  锚杆黏结段的离散化示意图

Fig. 6  Schematic diagram of discretization of anchor bonding segment

式(16)转换成对应的有限差分形式,表达式为

utxi+1-utxixi+1-xi-utxi-utxi-1xi-xi-1xi+1-xi-12=     Auiutxi+Bui (18)

式中:AuiBui分别为AuBu对应坐标xi的取值。

按照结点位移顺序对式(18)进行整理,可得

1xi+1-xiutxi+1-1xi+1-xi+1xi-xi-1+xi+1-xi-12Auiutxi+1xi-xi-1utxi-1-xi+1-xi-12Bui=0 (19)

相应地,锚杆黏结段自由端(结点x1)和受荷端(结点xn)的边界条件可分别表示为

1x2-x1utx2-1x2-x1+x2-x12Au1utx1-x2-x12Bu1=0 (20)
Ppreπrt2Et-1xn-xn-1+xn-xn-12Aunutxn+       1xn-xn-1utxn-1-xn-xn-12Bun=0 (21)

式中:Ppre为作用在受荷端的拉拔力的预测值。

按照图7所示的求解流程图,根据式(18)~式(21)编制MatLab程序对方程进行迭代求解,其具体步骤如下:

图7  锚杆荷载传递方程的求解流程

Fig. 7  Flow chart of solving load-transfer equation of anchor

1)根据离散分段数,由式(17)确定每个结点的坐标xi;将材料参数和尺寸参数代入式(12b)式(13b),确定参数CrDr的取值,需要说明的是,这两个参数取值与结点坐标无关。

2)根据拉拔力P0假定锚杆黏结段自由端的筋体位移utx1,将其代入式(15),确定τtx1-δtx1的关系,联合T-B界面黏结—滑移模型的函数表达式,即式(6a),确定结点x1处对应的T-B界面剪应力τtx1和界面剪切位移δtx1;将T-B界面剪切位移δtx1代入式(14),确定B-M界面剪切位移δmx1;将界面剪切位移δtx1δmx1代入式(8)式(9),确定切线参数kt1ct1km1cm1,并代入式(15)式(16)确定参数Au1Bu1;将utx1Au1Bu1代入式(20),确定utx2

3)根据utx2的取值,重复步骤2),确定参数Au2Bu2;将utx1utx2Au2Bu2的值代入式(19),求得utx3的取值。

4)重复步骤2)和步骤3),通过迭代式(19)求得其他结点的参数AuiBui及位移utxi的取值。

5)将utxn-1utxnAunBun的值代入式(21),求得Ppre的取值;按式(22)评估拉拔力预测值Ppre与实际值P0的误差。

Ppre/P0-1<0.01 (22)

式(22)成立,则输出所有结点位移utxi;否则,更新步骤2)中utx1的假定值,重复步骤2)~步骤5),直到式(22)得到满足。

6)将求得的结点位移utxi、界面剪应力τtxi代入式(5),确定锚杆筋体的轴力分布。

4 计算参数的试验确定

4.1 试验方案

为了对提出的荷载传递分析框架的有效性进行验证,分别开展单元体尺度和大尺度水泥土的钢管拉拔试验,两个尺度的试样尺寸参数如表1所示。

表1  两个尺度的水泥土锚杆拉拔试验的试样尺寸参数
Table 1  Dimension parameters of testing specimens of soil mixing anchor pull-out test in two scales
尺度黏结长度/mm钢管外径/mm钢管内径/mm试样直径/mm

筋体/水泥土

直径比

单元体 100 40 30 200 0.20
模型 800 40 30 250 0.16

为了消除试样套筒的边界效应,试样的筋体与水泥土直径比均不高于0.2[

21-22]。两种尺度的试验采用相同的水泥土(20%水泥掺入比,65%含水量以及28 d养护龄期)和筋体(直径40 mm的光面钢管),前者作为锚固界面特性测试用以确定界面模型参数,后者作为模型试验用以验证分析框架对拉拔响应预测的准确性。考虑到水泥土锚杆B-M界面滑移量测可控性难以保证,水泥土内剪切带的变形发展同样适用该变形模式,试验未进行B-M界面变形量测,而是将T-B界面剪切带与非剪切带形成的变形差异界面作为B-M替代界面处理。此外,为了确定水泥土材料的力学参数,进行了水泥土的无侧限压缩试验。上述试验现场如图8所示,具体试验步骤可参阅笔者前期研究报[21-22],此处不再赘述。

(a) 锚杆单元体拉拔试验

(b) 水泥土无侧限压缩试验

  

(c)  水泥土锚杆拉拔模型试验

图8 水泥土锚杆试验现场

Fig. 8 Laboratory view of soil mixing anchor tests

4.2 筋体与水泥土力学参数

水泥土无侧限压缩试验测得3个平行试样的压力—位移曲线如图9所示,根据曲线峰值点和线性段斜率计算水泥土的强度和变形参数。需要说明的是,模型计算中水泥土弹性模量取用平行试样试验结果的均值,水泥土剪切模量根据弹性模量和泊松比推算求[

23];试验用钢管为商用采购品,其钢材弹性模量在产品出厂参数中给出,依据截面等效刚度原则,钢管弹性模量可由钢材弹性模量确定,具体参数取值如表2所示。

图9  水泥土平行试样的压力—位移曲线

Fig. 9  Compression-displacement curves of parallel specimens of cemented soils

表2  筋体和水泥土力学参数
Table 2  Mechanical parameters of tested reinforcement and cemented soils
材料屈服强度/MPa拉伸模量Et/ GPa抗压强度/ MPa弹性模量Em/ MPa剪切模量Gm/ MPa
钢管 335 91.87 NA NA NA
水泥土 NA NA 3.35 269.05 126.91

注:  NA表示未对该性质进行测试。

4.3 锚固界面模型参数

水泥土锚杆单元体拉拔试验测得的界面黏结—滑移曲线如图10所示。可以发现,两个平行试样的曲线较为接近,表明试样的一致性控制较好,可取用均值曲线刻画界面特性。考虑到上述公式推导中界面模型已采用切线法表示,故采用三折线模型对均值曲线进行拟合,如图10所示,获得的界面模型参数取值如表3所示。

图10  界面模型拟合均值曲线

Fig. 10  The fitted mean curve of interface model

表3  界面模型参数取值
Table 3  Parameter values for interface model
τtf /kPaδtf /mmτtr /kPaδtr /mm
336.22 0.99 178.87 3.05

5 模型试验结果对比

上述水泥土锚杆试验中所用的锚杆试样未引入岩土体,即锚固界面只涉及T-B界面。考虑到锚杆拉拔过程中水泥土在T-B界面附近会形成一定厚度的剪切带,为了匹配双界面荷载传递分析框架的变形模式,将剪切带内水泥土看作图2中的黏结体,将剪切带外的水泥土看作岩土体,剪切带内外边界上不发生滑移。结合界面和水泥土特性测试所获得的力学参数及模型试验的尺寸参数,应用本文提出的荷载传递分析框架对水泥土锚杆模型试验的拉拔力学响应进行预测。将预测结果与量测结果进行对比,如图11所示。

(a)  拉拔力—位移响应

(b)  轴力分布响应

图11 锚杆力学响应预测—量测对比

Fig. 11 Comparison between predictions and measurements for mechanical response of soil ;mixing anchor

图11(a)可以看出,拉拔力—钢管加载端相对位移曲线的计算结果与量测数据的整体走势比较吻合,同一拉拔力对应的相对位移量测值比计算值偏大。产生这种偏差的原因在于,直接采用无侧限压缩试验数据作为水泥土压缩模量可能会高估水泥土的变形模量。值得说明的是,拉拔力—钢管自由端位移计算曲线与量测数据更为接近,考虑到该曲线本质上是界面黏结—滑移模型在更大尺度下的宏观表现,表明单元拉拔试验测得的界面黏结特性可以较准确地反映模型试验中的实际界面行为。由图11(b)可以看出,根据模型计算得到的钢管轴力分布与量测数据吻合良好,在拉拔力较小时,轴力计算值略低于量测值;在拉拔力接近峰值和在峰值后下降时,轴力计算值则略高于量测值。由于较小的拉拔力对应界面的最大剪应力尚未达到界面黏结强度,峰值附近的拉拔力对应界面的最大剪应力已经达到界面黏结强度且部分处于界面软化阶段,因此可以判断,在弹性段采用界面黏结—滑移模型会高估界面实际剪应力,在软化段则会低估界面实际剪应力。整体来看,轴力分布的计算值和量测值差异不大,说明采用锚杆单元拉拔的界面测试结果来刻画模型试验中的界面黏结特性是可以接受的。

表4比较了模型试验中钢管极限拉拔力及对应相对位移的计算值和量测值。可以看出,基于界面单元测试结果提出的水泥土锚杆双界面荷载传递分析模型可以较准确地预测筋体的极限拉拔力,预测偏差低于5%。此外,达到极限拉拔力所需的加载端相对位移的计算偏差低于15%,说明所采用的改进剪滞模型对水泥土实际变形响应的刻画效果大致可以接受,但仍有进一步优化的空间。

表4  锚杆极限拉拔力及对应位移的计算值与量测值对比
Table 4  Comparison between calculated and measured values for ultimate pullout resistance and corresponding displacement of anchor
数值极限拉拔力/kN加载端相对位移/mm
量测值 30.00 4.18
计算值 30.79 4.65
偏差/% 2.6 11.2

6 结论

提出一种可以考虑界面滑移耦合的水泥土锚杆双界面承载受力变形分析方法。通过开展水泥土锚杆单元体界面测试和水泥土力学试验确定该方法中所需计算参数,并进一步开展水泥土锚杆拉拔模型试验,对该方法对锚杆拉拔响应的预测效果进行了对比分析。主要结论如下:

1)锚固界面强度的发挥程度由界面两侧材料的相对位移决定,由于剪应力在界面两侧的连续性,界面两侧材料的变形存在明显的耦合效应。

2)界面剪应力在径向的传播与分布会受到正应力轴向分布的影响,计算剪切变形时应对剪滞模型进行修正。

3)基于切线法的界面剪切模型可对理想弹塑性和应变软化等不同锚固界面的黏结—滑移特性进行统一表达,从而有效提高锚杆荷载传递模型的解算效率。

4)水泥土锚杆双界面承载分析方法可以有效连接基于锚杆拉拔的单元体界面测试结果和大尺度模型试验结果,从而建立从单元体尺度锚固界面特性到大尺度锚杆拉拔响应的分析框架和预测模型。

参考文献

1

王树仁, 澳)保罗·黑根(Paul C. Hagan), (澳)曹晨. 岩体锚固与岩土工程新进展[M]. 北京: 清华大学出版社, 2016: 1-18. [百度学术] 

WANG S R, HAGAN P C, CAO C. Advances in rock-support and geotechnical engineering [M]. Beijing: Tsinghua University Press, 2016: 1-18. (in Chinese) [百度学术] 

2

匡政, 张明义, 白晓宇, . 风化岩地基GFRP抗浮锚杆力学与变形特性现场试验[J]. 岩土工程学报, 2019, 41(10): 1882-1892. [百度学术] 

KUANG Z, ZHANG M Y, BAI X Y, et al. Field tests on mechanics and deformation properties of GFRP anti-floating anchors in decomposed rock foundation [J]. Chinese Journal of Geotechnical Engineering, 2019, 41(10): 1882-1892. (in Chinese) [百度学术] 

3

李金峰, 张景科, 王南, . 夯筑土遗址木锚杆四锚系统锚固机制物理模型试验研究[J]. 岩石力学与工程学报, 2019, 38(11): 2321-2331. [百度学术] 

LI J F, ZHANG J K, WANG N, et al. Physical model test on working mechanisms of an anchor group system with four wood bolts in rammed earthen sites [J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(11): 2321-2331. (in Chinese) [百度学术] 

4

MA S Q, NEMCIK J, AZIZ N. An analytical model of fully grouted rock bolts subjected to tensile load [J]. Construction and Building Materials, 2013, 49: 519-526. [百度学术] 

5

加筋水泥土桩锚技术规程: CECS 147—2016 [S]. 北京: 中国计划出版社, 2016. [百度学术] 

Technical specification for soil mass with reinforced cement soil pile and anchors: CECS 147-2016 [S]. Beijing: China Planning Press, 2016. (in Chinese) [百度学术] 

6

黄明华, 李嘉成, 赵明华, . 层状地基中锚杆拉拔荷载传递非线性分析[J]. 中国公路学报, 2019, 32(1): 12-20, 56. [百度学术] 

HUANG M H, LI J C, ZHAO M H, et al. Nonlinear analysis on load transfer mechanism of bolts in layered ground [J]. China Journal of Highway and Transport, 2019, 32(1): 12-20, 56. (in Chinese) [百度学术] 

7

何栋梁, 杨伟军, 林杭, . 拉—剪状态下锚固节理剪切强度与破坏模式分析[J]. 公路交通科技, 2019, 36(4): 95-101. [百度学术] 

HE D L, YANG W J, LIN H, et al. Analysis on shear strength and failure mode of anchored joint under tensile-shear load [J]. Journal of Highway and Transportation Research and Development, 2019, 36(4): 95-101. (in Chinese) [百度学术] 

8

郭锐剑, 谌文武, 段建, . 考虑界面软化特性的土层锚杆拉拔受力分析[J]. 中南大学学报(自然科学版), 2012, 43(10): 4003-4009. [百度学术] 

GUO R J, CHEN W W, DUAN J, et al. Pullout mechanical analysis of soil anchor based on softening behavior of interface [J]. Journal of Central South University (Science and Technology), 2012, 43(10): 4003-4009. (in Chinese) [百度学术] 

9

黄明华, 李嘉成, 赵明华, . 岩体离层作用下锚杆受力特性全历程分析[J]. 岩石力学与工程学报, 2017, 36(9): 2177-2184. [百度学术] 

HUANG M H, LI J C, ZHAO M H, et al. Full-range analysis on mechanical characteristics of bolts in rock mass with bedding separation [J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(9): 2177-2184. (in Chinese) [百度学术] 

10

MA S Q, NEMCIK J, AZIZ N, et al. Analytical model for rock bolts reaching free end slip [J]. Construction and Building Materials, 2014, 57: 30-37. [百度学术] 

11

MA S Q, ZHAO Z Y, NIE W, et al. A numerical model of fully grouted bolts considering the tri-linear shear bond-slip model [J]. Tunnelling and Underground Space Technology, 2016, 54: 73-80. [百度学术] 

12

杨嵘昌, 宰金珉. 广义剪切位移法分析桩—土—承台非线性共同作用原理[J]. 岩土工程学报, 1994(6): 103-116. [百度学术] 

YANG R C, ZAI J M. The principle of the generalized shear displacement method for analyzing the nonlinear interaction of pile-soil-pilecap [J]. Chinese Journal of Geotechnical Engineering, 1994(6): 103-116. (in Chinese) [百度学术] 

13

CAI Y, ESAKI T, JIANG Y J. A rock bolt and rock mass interaction model [J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(7): 1055-1067. [百度学术] 

14

YIN J H, HONG C Y, ZHOU W H. Simplified analytical method for calculating the maximum shear stress of nail-soil interface [J]. International Journal of Geomechanics, 2012, 12(3): 309-317. [百度学术] 

15

FABRIS C, SCHWEIGER H F, PULKO B, et al. Numerical simulation of a ground anchor pullout test monitored with fiber optic sensors [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2021, 147(2): 04020163. [百度学术] 

16

岩土锚杆与喷射混凝土支护工程技术规范: GB 50086—2015 [S]. 北京: 中国建筑工业出版社. [百度学术] 

Technical code for engineering of ground anchorages and shotcrete support: GB 50086-2015 [S]. Beijing: China Architecture & Building Press. (in Chinese) [百度学术] 

17

宋洋, 王贺平, 许怀玉, . 全长黏结锚杆锚—浆界面应力的分布规律[J]. 土木与环境工程学报(中英文), 2019, 41(4): 44-50. [百度学术] 

SONG Y, WANG H P, XU H Y, et al. Analysis of stress distribution patten along the anchor-slurry interface of full-length bond anchor [J]. Journal of Civil and Environmental Engineering, 2019, 41(4): 44-50. (in Chinese) [百度学术] 

18

任非凡, 徐超, 谌文武. 多界面复合锚杆荷载传递机制的数值模拟[J]. 同济大学学报(自然科学版), 2011, 39(12): 1753-1759. [百度学术] 

REN F F, XU C, CHEN W W. Numerical simulation on load transfer mechanism of multi-interface composite rockbolt [J]. Journal of Tongji University (Natural Science), 2011, 39(12): 1753-1759. (in Chinese) [百度学术] 

19

陈昌富, 梁晓腾, 张根宝. 考虑土体变形与界面剪切耦合作用锚杆拉拔响应[J]. 科学技术与工程, 2018, 18(23): 276-280. [百度学术] 

CHEN C F, LIANG X T, ZHANG G B. Pullout response of anchor under coupled effect of soil deformation and interface shear [J]. Science Technology and Engineering, 2018, 18(23): 276-280. (in Chinese) [百度学术] 

20

ZHANG G B, CHEN C F, ZORNBERG J G, et al. Interface creep behavior of grouted anchors in clayey soils: Effect of soil moisture condition [J]. Acta Geotechnica, 2020, 15(8): 2159-2177. [百度学术] 

21

CHEN C F, ZHANG G B, ZORNBERG J G, et al. Interface behavior of tensioned bars embedded in cement-soil mixtures [J]. Construction and Building Materials, 2018, 186: 840-853. [百度学术] 

22

CHEN C, ZHANG G, ZORNBERG J G, et al. Interface bond behavior of tensioned glass fiber-reinforced polymer (GFRP) tendons embedded in cemented soils [J]. Construction and Building Materials, 2020, 263: 120132. [百度学术] 

23

陈甦, 彭建忠, 韩静云, . 水泥土强度的试件形状和尺寸效应试验研究[J]. 岩土工程学报, 2002, 24(5): 580-583. [百度学术] 

CHEN S, PENG J Z, HAN J Y, et al. Experimental research on the shape and size effect of cement-stabilized soil specimen strength [J]. Chinese Journal of Geotechnical Engineering, 2002, 24(5): 580-583. (in Chinese) [百度学术]