Shear strength of reinforced soil interface under normal cyclic loading and its prediction
-
摘要: 为了研究法向动荷载作用下筋土界面的剪切应力变化规律,采用大型动态直剪仪,对相对密实度为75%的砾石和土工格栅的界面进行了剪切试验,研究了4种法向初始应力(20,40,60,80 kPa)和4种法向荷载振动幅值(10,20,30,40 kPa)对筋土界面循环剪切特性的影响。在试验的基础上,建立了法向循环荷载作用下筋土界面在峰值前和残余阶段的剪切应力-法向应力表达式。同时,结合应力时间差规律,并考虑法向初始应力和荷载振幅等的影响,提出了界面剪切应力-剪切位移的表达式。将两种预测表达式与试验结果进行了对比,均具有较好的吻合度,验证了方法的正确性。Abstract: In order to study the variation of shear stress at the interface between reinforced soil and geogrid under normal cyclic loading, a large-scale dynamic direct shear apparatus is used to conduct the shear tests on the interface between gravel and geogrid with relative density of 75%. The effects of four initial normal stresses (20, 40, 60 and 80 kPa) and four normal load vibration amplitudes (10, 20, 30 and 40 kPa) on the interfacial cyclic shear characteristics are studied. On the basis of the tests, the expression for the normal stress at the interface before peak and at the residual stage under normal cyclic loading is established. Considering the influences of the initial normal stress and load vibration amplitudes, the expression for the interfacial shear stress-shear displacement is proposed based on the relative phase shift law. The results by the two proposed prediction expressions are in good agreement with the experimental ones, which verifies the correctness of the proposed method.
-
0. 引言
土工合成材料已经在土木和水利领域得到了广泛使用,而土工合成材料与土体界面的相互作用往往是影响工程安全性的关键因素[1-2]。直剪试验是评估界面力学特性的重要方式,但相关的研究主要针对法向静荷载的工况[3],较少考虑法向动荷载的作用。刘飞禹等[4]、徐超等[5]、王家全等[6]从室内试验和数值模拟等角度探索了法向循环荷载下筋土界面的剪切行为,并指出了动力响应的复杂性。
在界面剪切理论研究方面,已有学者提出了一些描述界面力学特征的理论模型。Basudhar等[7]基于邓肯-张模型和Kondner[8]对模型参数的研究,用双曲线模型描述了土–土工布界面直剪试验达到峰值前的规律。王军等[9]将Esterhuizen等[10]提出的位移软化模型用于拟合达到峰值后的砂土–土工格栅界面剪切应力–位移的非线性关系,取得了较好的拟合效果。Toufigh等[11]对碳纤维增强聚合物–回填土界面进行了一系列剪切试验,并在此基础上建立了筋土界面的塑性本构模型。Zornberg等[12]发展了一套小位移下的土–土工合成材料界面模型,并将其应用于实际道路工程设计中。陈榕等[13]采用罚函数法模拟界面纵肋作用,而在横肋上应用黏聚力模型,发现新的方法比传统界面模拟方法更准确。杜常博等[14]为了描述土工合成材料拉拔试验的应变软化、硬化现象,分别提出了双线性弹塑性硬化模型和三线性应力–位移软化模型,并与试验结果、几种经典界面模型进行了对比。Cen等[15]基于土工膜与砂砾在不同剪切幅值、剪切速率下的界面循环剪切试验,建立了筋土界面相互作用的计算模型。
由上述分析可知,在筋土界面单调剪切特性和循环剪切特性的计算模型方面,已有了一些研究。然而以往的理论分析大多是基于法向静荷载的条件,缺乏法向动荷载下的界面强度预测研究。本文在试验研究的基础上,进行了法向循环荷载下筋土界面剪切应力–法向应力和剪切应力–剪切位移的相关理论研究,并将所得结果与试验进行了对比验证。
1. 试验方案及界面剪切行为
1.1 试验材料与方案
试验所用仪器为微机控制电液伺服动态直剪仪,相关介绍见文献[16]。试验用土样为相对密实度为75%的砾石:相对质量密度GS为2.62;最大、最小孔隙比分别为0.66,0.47。级配分布特征:平均粒径D50为5.03 mm,有效粒径D10为3.19 mm,连续粒径D30为4.30 mm,限制粒径D60为5.41 mm,不均匀系数Cu为1.70,曲率系数Cc为1.07。试验用加筋材料为聚丙烯塑料土工格栅,单位面积的质量为330 g/m2,其网孔尺寸为25 mm×25 mm,横、纵肋的厚度均为5 mm,横、纵向的极限延伸率分别为13%和15%,极限抗拉强度为20 kN/m。
法向荷载加载方式见图 1,选取的试验采用相同的法向加载频率(0.5 Hz)、正弦波形、水平剪切速率(1 mm/min)和最大剪切位移(90 mm),而法向初始应力σi和振幅σa不相同,法向荷载分别为20~50,40~70,60~90,80~110 kPa和60~70,60~80,60~90,60~100 kPa。
1.2 界面剪切行为
法向循环荷载作用(初始应力40 kPa,振幅30 kPa)下的筋–土界面剪切特性如图 2所示。整体上表现出剪切应力逐渐增加至峰值,而后软化至残余状态的特征,且其值在一定范围内动态变化。剪切应力与法向应力相似,随时间以接近正弦波形循环变化,但二者的变化相位不一致,剪切应力峰值滞后于法向应力峰值。根据图 2(c)可看出两个方向的应力发展路径,随着法向应力的循环变化,剪切应力也在循环变化。
2. 剪切应力–法向应力规律预测
从应力比(剪切应力/法向应力)与法向应力的相对相位关系(图 2(b))可以看出,应力比的波形相较于剪切应力更接近正弦波形,且剪切应力的相邻峰、谷值间隔均约为半个周期,这与Dang等[17-18]进行的混凝土界面剪切试验现象类似。基于这一规律,可建立剪切应力与法向应力的关系式为
τσn=c1+c2sin[2π f(t+Δt)], (1) σn=(σi+0.5σa)+0.5σasin(2π ft), (2) TR=ΔtT, (3) 式中,c1,c2为表达式的经验系数,△t为应力比与法向应力的时间差,σi为初始应力,σa为振幅,f为频率,TR为相对时间差,以周期为单位,T为一个循环周期所经历的时间。
图 3(a)为选取的相同振幅(30 kPa)、不同初始应力(20,40,60,80 kPa)作用下峰值阶段的应力变化曲线,从图 3中可以看出,不管初始应力值如何改变,应力比与法向应力均存在约0.5个周期的相对时间差。同时,随着初始应力的增加,最大应力比在减小,而最小应力比在增大。将图中4个最大应力比和4个最小应力比取平均值,可计算得到c1和c2,并用对数函数进行拟合,得到峰值阶段的系数与初始应力的表达式。
图 3(c)为选取的与图 3(a)相同条件下残余阶段最后4个循环周期的应力变化曲线,不管初始应力值如何改变,应力比与法向应力均存在约0.5个周期的相对时间差。随着初始应力的增加,最大应力比在减小,而最小应力比在增大,计算可得到残余阶段的系数与初始应力的表达式。
由图 4可见,在相同初始应力(60 kPa)、不同振幅(10,20,30,40 kPa)条件下,无论是峰值还是残余阶段,应力比与法向应力均存在约0.5个周期的相对时间差。随着振幅的增加,最大应力比在增大,而最小应力比在减小。用线性函数拟合,可得到经验系数c1和c2随振幅变化的表达式。
根据上述讨论与计算,得到峰值和残余阶段剪切应力–法向应力的表达式如下:
τ={c1+c2sin[2π f(t+1)]}σn。 (4) 将剪切试验在峰值和残余阶段的剪切应力随法向应力变化曲线与根据式(4)计算得到的预测曲线进行了对比,部分数据如图 5所示。结果表明,计算结果与试验结果吻合较好,验证了本文所建立的表达式的正确性。
3. 剪切应力–剪切位移规律预测
3.1 剪切应力–剪切位移关系分析
图 6(a)是剪切位移为10,20,40,80 mm时的法向应力与应力比随时间变化的曲线,可见各阶段的数据均满足相位差为0.5个周期的规律,且应力和应力比都可看作稳定的循环变化。图 6(b)为应力比随剪切位移变化的曲线,可见应力比在一定范围内随位移变化,上、下界控制应力比的变化区间,可拟合上、下界曲线得到应力比的上下限变化规律。
上、下界剪切应力–剪切位移表达式的建立思路是:先建立剪切位移与时间的表达式,并将该表达式代入式(4),再将式(4)中的系数c1和c2表达成剪切位移的函数,即可得到用剪切位移表达的上、下界剪切应力–剪切位移表达式。以剪切位移零点为时间起点,可知t=δ/ν,其中,t为时间参数,δ为剪切位移,ν为剪切速率。由于c1和c2受到剪切位移的影响,可将其表达为剪切位移的函数:
c1=f1(δ) ,c2=f2(δ) 。} (5) 将式(5)代入式(4)得
τ={f1(δ)+f2(δ)sin[2πf(δv+1)]}[(σi+0.5σa)+0.5σasin(2πfδv)]。 (6) 由图 6(b)可以发现上、下界曲线分峰前、峰后两个阶段,峰前类似于应力–应变的邓肯-张模型:
Rf=δa+bδ。 (7) 峰后类似于文献[8]提出的非线性模型,对峰后剪切应力进归一化处理得到残余系数:
S=Rp−RbRp−Rr。 (8) 又由于上、下界峰后剪切位移与残余系数形成S形的曲线,且残余系数在0~1变化,可用下式描述二者间的关系:
S=1−exp(−A(δ−δp)B。 (9) 由此可得到峰后应力比与位移的关系:
Rb=Rp−(Rp−Rr)(1−exp(−A(δ−δp)B) , (10) 式中,Rf,Rb为峰前和峰后的边界曲线应力比,S为残余系数,a,b,A,B为经验参数,受法向荷载等的影响,δp为峰值应力比对应的剪切位移。
在确定了应力比的内部变化形式和边界曲线后,可得到f1(δ)和f2(δ)的表达式为
f1(δ)={Rf,u+Rf,l2 (δ≤δp)Rb,u+Rb,l2(δ>δp) , (11) f2(δ)={Rf,u+Rf,l2 (δ≤δp)Rb,u+Rb,l2 (δ>δp), (12) 式中,Rf,u,Rf,l,Rb,u,Rb,l为上界峰前应力比、下界峰前应力比、上界峰后应力比、下界峰后应力比。
通过以上步骤的推导,可得到式(6)中未确定的参数,进而得出剪切应力τ与剪切位移δ的关系式。
3.2 参数确定
若能确定应力比峰值前参数a和b、峰值后参数A和B,即可得到应力比偏距函数f1(δ)和振幅函数f2(δ),进而推出剪切应力–剪切位移的关系式。法向应力受初始应力和振幅这两个影响因素控制,需要分别考虑这两种工况下的参数。
(1)不同初始应力
根据式(7),参照邓肯-张非线性模型分析方法,δ/Rf与δ形成线性关系,通过绘制并拟合两者的关系曲线,可得到线性方程的两个参数:截距a和斜率b。通过比较发现,若用0~100%的应力比水平所相应的点进行拟合,在峰值处与试验结果存在一定误差,故而本文采用文献[11]的方法,取70%~95%应力比水平绘制δ/Rf−δ曲线,以减小结果误差,所得的参数见表 1。总的来说,曲线上界参数随初始应力增加而增大,而曲线下界参数随初始应力变化有增有减,数值比较接近。
表 1 不同初始应力峰前参数Table 1. Pre-peak parameters of different initial stresses法向应力/kPa 应力比 参数a 参数b 20~50 上界 1.4169 0.4599 下界 2.6122 0.9694 40~70 上界 1.8774 0.5513 下界 3.0433 0.8234 60~90 上界 2.0269 0.6465 下界 2.8182 0.8959 80~110 上界 2.4041 0.6821 下界 3.0726 0.8712 对于应力比曲线上、下界峰后表达式,主要受峰值应力比、残余应力比和曲线形态控制,峰值和残余应力比统计见表 2。依据式(9)绘制上、下界峰后剪切位移与残余系数的关系曲线,如图 7,拟合不同初始应力下的试验数据,发现可以较好的描述上下界的非线性变化规律。不同初始应力下参数A为0.02559,B为1.73186,相关系数R2为0.92。
表 2 应力比汇总Table 2. Summary of stress ratio法向应力/kPa 应力比 峰值应力比 残余应力比 20~50 上界 1.5374 1.1121 下界 0.7728 0.5624 40~70 上界 1.3343 1.0260 下界 0.8645 0.6537 60~90 上界 1.2507 0.9214 下界 0.9137 0.6628 80~110 上界 1.2059 0.8806 下界 0.9347 0.6736 (2)不同法向振幅
利用式(7),通过对不同法向振幅下的直剪试验数据的峰前δ/Rf与δ进行回归分析,计算可得参数值(表 3)。总体来说,应力比曲线下界参数随初始应力增加而较小,而曲线上界参数随初始应力变化有增有减,数值比较接近。不同振幅下峰后参数A为0.006,B为2.25831,相关系数R2为0.93(图 8)。
表 3 不同振幅峰前参数Table 3. Pre-peak model parameters of different amplitudes法向应力/kPa 应力比 参数a 参数b 60~70 上界 2.1310 0.6618 下界 2.3306 0.7647 60~80 上界 1.9142 0.6565 下界 2.2515 0.8247 60~90 上界 2.0269 0.6465 下界 2.8182 0.8959 60~100 上界 2.1722 0.6518 下界 2.9624 0.9790 3.3 预测验证
图 9为采用本文所提方法的预测结果与试验结果的对比验证图。从图中可以看出,基于本文所提方法得到的拟合曲线既能较好地反映应力变化的整体趋势,也能代表界面循环剪切应力的变化规律,表明预测结果与试验结果符合的较好,验证了本文所提方法的正确性。
4. 结论
本文根据法向循环荷载作用下筋土界面试验结果,对界面剪切行为进行了预测研究,主要得到3点结论。
(1)根据应力比与法向应力这两种循环波形的相位差特点,建立了一种新型的峰值和残余阶段的剪切应力–法向应力关系式,可以反映峰值和残余阶段剪切应力的变化规律。
(2)基于应力比循环变化特点和应力比–位移上、下界曲线模型,提出了剪切应力–剪切位移关系表达式,可以描述法向循环荷载下筋土界面的软化规律。
(3)两种预测研究的参数均考虑了法向初始应力和振幅的影响,并得到了应力水平对表达式参数的影响规律,为后续研究打下了基础。
-
表 1 不同初始应力峰前参数
Table 1 Pre-peak parameters of different initial stresses
法向应力/kPa 应力比 参数a 参数b 20~50 上界 1.4169 0.4599 下界 2.6122 0.9694 40~70 上界 1.8774 0.5513 下界 3.0433 0.8234 60~90 上界 2.0269 0.6465 下界 2.8182 0.8959 80~110 上界 2.4041 0.6821 下界 3.0726 0.8712 表 2 应力比汇总
Table 2 Summary of stress ratio
法向应力/kPa 应力比 峰值应力比 残余应力比 20~50 上界 1.5374 1.1121 下界 0.7728 0.5624 40~70 上界 1.3343 1.0260 下界 0.8645 0.6537 60~90 上界 1.2507 0.9214 下界 0.9137 0.6628 80~110 上界 1.2059 0.8806 下界 0.9347 0.6736 表 3 不同振幅峰前参数
Table 3 Pre-peak model parameters of different amplitudes
法向应力/kPa 应力比 参数a 参数b 60~70 上界 2.1310 0.6618 下界 2.3306 0.7647 60~80 上界 1.9142 0.6565 下界 2.2515 0.8247 60~90 上界 2.0269 0.6465 下界 2.8182 0.8959 60~100 上界 2.1722 0.6518 下界 2.9624 0.9790 -
[1] WANG J, LIU F Y, WANG P, et al. Particle size effects on coarse soil-geogrid interface response in cyclic and post-cyclic direct shear tests[J]. Geotextiles and Geomembranes, 2016, 44(6): 854–861. doi: 10.1016/j.geotexmem.2016.06.011
[2] FERREIRA F B, VIEIRA C S, LOPES M L, et al. HDPE geogrid-residual soil interaction under monotonic and cyclic pullout loading[J]. Geosynthetics International, 2020, 27(1): 79–96. doi: 10.1680/jgein.19.00057
[3] LASHKARI A, JAMALI V. Global and local sand–geosynthetic interface behaviour[J]. Géotechnique, 2021, 71(4): 346–367. doi: 10.1680/jgeot.19.P.109
[4] 刘飞禹, 朱晨, 王军. 剪切速率和法向加载频率对筋土界面剪切特性的影响[J]. 岩土工程学报, 2021, 43(5): 832–840. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18614.shtml LIU Fei-yu, ZHU Chen, WANG Jun. Influences of shear rate and loading frequency on shear behavior of geogrid-soil interfaces[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(5): 832–840. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18614.shtml
[5] 徐超, 石志龙. 循环荷载作用下筋土界面抗剪特性的试验研究[J]. 岩土力学, 2011, 32(3): 655–660. doi: 10.3969/j.issn.1000-7598.2011.03.003 XU Chao, SHI Zhi-long. Experimental research on shearing resistance property of sand-geogrid interface under cyclic load[J]. Rock and Soil Mechanics, 2011, 32(3): 655–660. (in Chinese) doi: 10.3969/j.issn.1000-7598.2011.03.003
[6] 王家全, 王宇帆, 黄世斌, 等. 循环荷载作用下土工格栅剪切特性的颗粒流细观分析[J]. 水利学报, 2014, 45(9): 1082–1090. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201409009.htm WANG Jia-quan, WANG Yu-fan, HUANG Shi-bin, et al. The particle flow mesoscopic analysis of geogrid shear properties under cyclic loading[J]. Journal of Hydraulic Engineering, 2014, 45(9): 1082–1090. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201409009.htm
[7] ANUBHAV, BASUDHAR P K. Modeling of soil-woven geotextile interface behavior from direct shear test results[J]. Geotextiles and Geomembranes, 2010, 28(4): 403–408. doi: 10.1016/j.geotexmem.2009.12.005
[8] KONDNER R L. Hyperbolic stress-strain response: cohesive soils[J]. Journal of the Soil Mechanics and Foundations Division, 1963, 89(1): 115–143. doi: 10.1061/JSFEAQ.0000479
[9] 王军, 林旭, 符洪涛. 砂土-格栅筋土界面特性的本构模型研究[J]. 岩土力学, 2014, 35(增刊2): 75–84. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2014S2011.htm WANG Jun, LIN Xu, FU Hong-tao. Study of constitutive model of sand-geogrid interface behavior in geogrid/geotextile reinforced soil[J]. Rock and Soil Mechanics, 2014, 35(S2): 75–84. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2014S2011.htm
[10] ESTERHUIZEN J J B, FILZ G M, DUNCAN J M. Constitutive behavior of geosynthetic interfaces[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2001, 127(10): 834–840. doi: 10.1061/(ASCE)1090-0241(2001)127:10(834)
[11] TOUFIGH V, DESAI C S, SAADATMANESH H, et al. Constitutive modeling and testing of interface between backfill soil and fiber-reinforced polymer[J]. International Journal of Geomechanics, 2014, 14(3): 04014009. doi: 10.1061/(ASCE)GM.1943-5622.0000298
[12] ZORNBERG J G, ROODI G H, GUPTA R. Stiffness of soil–geosynthetic composite under small displacements: I. model development[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2017, 143(10): 04017075. doi: 10.1061/(ASCE)GT.1943-5606.0001768
[13] 陈榕, 李博, 郝冬雪, 等. 基于黏聚力模型的土工格栅筋土界面作用模拟方法[J]. 岩土工程学报, 2020, 42(5): 934–940. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18206.shtml CHEN Rong, LI Bo, HAO Dong-xue, et al. Simulation for interaction between geogrids and soil by cohesive zone model[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(5): 934–940. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18206.shtml
[14] 杜常博, 易富. 筋土拉拔界面弹塑性模型的全过程分析[J]. 煤炭学报, 2020, 45(12): 4062–4073. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202012010.htm DU Chang-bo, YI Fu. Full-range analysis of elastic-plastic model of pull-out interface between geosynthetics and soil[J]. Journal of China Coal Society, 2020, 45(12): 4062–4073. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202012010.htm
[15] CEN W J, BAUER E, WEN L S, et al. Experimental investigations and constitutive modeling of cyclic interface shearing between HDPE geomembrane and sandy gravel[J]. Geotextiles and Geomembranes, 2019, 47(2): 269–279. doi: 10.1016/j.geotexmem.2018.12.013
[16] YING M J, LIU F Y, WANG J, et al. Coupling effects of particle shape and cyclic shear history on shear properties of coarse-grained soil-geogrid interface[J]. Transportation Geotechnics, 2021, 27: 100504. doi: 10.1016/j.trgeo.2020.100504
[17] DANG W G, KONIETZKY H, FRÜHWIRT T. Direct shear behavior of planar joints under cyclic normal load conditions: effect of different cyclic normal force amplitudes[J]. Rock Mechanics and Rock Engineering, 2017, 50(11): 3101–3107. doi: 10.1007/s00603-017-1270-7
[18] DANG W G, KONIETZKY H, CHANG L F, et al. Velocity-frequency-amplitude-dependent frictional resistance of planar joints under dynamic normal load (DNL) conditions[J]. Tunnelling and Underground Space Technology, 2018, 79: 27–34. doi: 10.1016/j.tust.2018.04.038
-
期刊类型引用(3)
1. 姚文杰,李明宝,郑俊杰,郭亿辉. 筋土峰值拉拔力的影响因素研究. 森林工程. 2024(01): 201-206 . 百度学术
2. 侯娟,楚辰玺,孙银玉,文平,马艳霞,张吾渝. 法向循环荷载下格栅加筋盐渍土剪切性能研究. 岩石力学与工程学报. 2023(S1): 3832-3841 . 百度学术
3. 张琬,陈建峰. 加筋土挡墙墙趾界面剪切刚度计算模型研究. 水资源与水工程学报. 2023(04): 152-158 . 百度学术
其他类型引用(5)