張貴生 杜衍慶 練象平 路遙

【摘? ? 要】：筋-土界面作用機理對研究加筋土結構性能具有重要意義?；隍炞C的有限元數(shù)值方法模擬筋-土界面拉拔過程，研究界面循環(huán)正應力初值、幅值和頻率對筋-土界面拉拔特性的影響，結果表明：增加筋-土界面循環(huán)正應力頻率或幅值，導致筋-土界面相互嵌固咬合和摩擦作用力降低。因此，高頻重載作用下加筋土結構的穩(wěn)定性應引起重視。

【關鍵詞】：循環(huán)動載；筋-土界面；拉拔試驗；加筋土結構

【中圖分類號】：TU472.34【文獻標志碼】：C【文章編號】：1008-3197（2023）03-04-05

【DOI編碼】：10.3969/j.issn.1008-3197.2023.03.002

Numerical Analysis of Pull-out Characteristics of Reinforcement-soil Interface

under Cyclic Dynamic Loading

ZHANG Guisheng1，DU Yanqing1， LIAN Xiangping1， LU Yao2

（1. Tianjin Municipal Engineering Design & Research Institute Co. Ltd.， Tianjin 300392， China;

2. School of Civil and Transportation Engineering， Hebei University of Technology， Tianjin 300401 ， China）

【Abstract】：The interaction mechanism of reinforcement-soil interface is of great significance to study the structural performance of reinforced soil. Based on the verified finite element numerical method， the pull-out process of the reinforcement-soil interface was simulated， and the influence of the initial value， amplitude and frequency of the interface cyclic normal stress on the pull-out characteristics of the reinforcement-soil interface was studied. The results showed that when the frequency and amplitude of interface cyclic normal stress was increased， mutual embedding， occlusion and friction between reinforcement and soil interface decreased. The stability of reinforced soil structure under high frequency and heavy load should be paid attention to.

【Key words】：cyclic dynamic load ; reinforcement -soil interface; pull-out test； reinforced soil structure

加筋土通過向填料中加入筋材來增加自身強度，筋材與填料土體的相互作用決定加筋效果；因此，開展筋-土界面相互作用研究具有重要意義。目前，拉拔試驗或直剪試驗是國內(nèi)外研究筋-土界面相互作用機理的重要試驗手段，在筋材與填土類型、界面參數(shù)等對作用機理的影響方面取得了較為豐碩的成果[1～2]；但受模型尺寸與條件、標準選擇等因素影響，試驗結果在全面反映筋-土界面相互作用機理方面有待進一步完善。

現(xiàn)有大量的研究結果表明[3～4]，基于有限元或離散元等數(shù)值計算方法，通過設置合理的筋-土界面本構模型，可有效對筋-土界面特性開展多因素、多尺度的變參數(shù)分析并研究筋材拉拔力及位移、應力和應變的發(fā)生、發(fā)展與變化規(guī)律。如Jensen等[5]、Wang等[6]和易富等[7]采用PFC離散元軟件模擬筋-土間的直剪試驗，研究了直剪過程中筋-土界面特性及筋材周邊土體的變化規(guī)律，重點研究了筋-土界面間力的傳導機理，在格柵加筋效果的最佳孔徑比等方面取得了較好的成果。

既有研究主要針對靜載作用下筋-土界面特性，考慮到實際中大量加筋土用于道路支擋結構中，筋-土界面或筋材往往承受交通循環(huán)荷載作用，但動載作用下筋-土界面特性的研究相對較少。本文基于有限元數(shù)值方法，模擬研究拉拔試驗中循環(huán)正應力為循環(huán)荷載作用下的筋-土界面作用特性。

1 筋-土界面拉拔特性數(shù)值模擬

1.1 模型的建立

采用ABAQUS有限元軟件平面應變狀態(tài)來模擬拉拔過程，建模尺寸與既有實際箱子尺寸一致，長×寬為800 mm×420 mm，格柵筋材長度和厚度分別1 040 、1 mm，網(wǎng)格劃分尺寸設定0.04 m。見圖1。

模擬拉拔試驗筋材從左側加載拉拔，土體底部邊界限制位移，左側即加載一側土體受箱子邊界限制水平方向位移。

填料土體采用摩爾-庫倫屈服準則的彈塑性本構、筋材采用彈性本構模擬，采用罰函數(shù)模擬筋-土界面摩擦接觸屬性，根據(jù)室內(nèi)拉拔試驗測定格柵和填土參數(shù)。見表1。

通過3個分析步完成，即初始分析步、靜載施加分析步和動載施加分析步。靜載分析步施加重力荷載和恒定荷載；動載分析步建立時間-荷載關系，施加循環(huán)荷載。通過設置格柵和上下土體的接觸屬性完成筋-土相互作用的構建，筋-土界面摩擦系數(shù)設置為0.9。

1.2 模擬方案

為分析動載作用下的筋-土界面特性，筋-土界面正應力為動載循環(huán)作用，即重點研究界面正應力初始值s0、幅值An和頻率fn變化時，筋材拉拔特性和筋材位移沿長度分布的特征。見表2。

1.3 模擬驗證

為了驗證數(shù)值模擬方法的準確性，選取既有拉拔試驗結果[8]進行對比分析，以界面正應力sn=10 kPa為例，對比分析拉拔過程中不同加載端位移時筋材沿長度方向位移變化。不同加載端位移的數(shù)值模擬和拉拔試驗測得的筋材位移量相差較小，表明模擬筋-土界面拉拔試驗的數(shù)值方法是可靠的。見圖2。

選取sn為50～75 kPa、頻率f為0.1 Hz對比分析界面正應力變化時剪應力與加載端位移的變化。數(shù)值模擬結果較試驗結果曲線走勢更為順滑、曲線波動范圍更小，可準確完成整個循環(huán)荷載的施加，較好反映試驗數(shù)據(jù)結果走勢，基于本文數(shù)值模擬開展界面動載作用下的筋-土界面特性研究是可行性的。見圖3。

2 模擬結果分析

2.1 界面循環(huán)正應力初值對拉拔特性的影響

選取3種界面循環(huán)正應力進行分析，即An、fn和拉拔速率相同，僅改變s0，模擬所得筋材拉力與加載端位移。初始階段，拉拔力均隨加載端位移增加而增加，由于界面正應力為循環(huán)荷載，拉拔力同加載端位移關系呈波動變化，隨拉拔進行，拉拔力逐漸達到峰值，后續(xù)拉拔力達到相對穩(wěn)定狀態(tài)，此后加載端位移持續(xù)增加，而拉拔力保持穩(wěn)定，終至筋材被拔出破壞（20 mm），模擬拉拔過程結束；而且，隨界面循環(huán)正應力初值增大，拉拔力峰值呈遞增趨勢且進入穩(wěn)定波動狀態(tài)時加載端位移隨之增加，即達到穩(wěn)定狀態(tài)的時間后移。見圖4和圖5。

分析界面循環(huán)正應力初值對拉拔過程中筋材不同位置處的位移變化影響，相同界面循環(huán)正應力初值下，隨加載端位移增加，測點1、2和3相繼發(fā)生移動且隨著距加載端越遠，出現(xiàn)初始位移時加載端位移亦越大。另外，隨界面循環(huán)正應力初值的增大，初始位移的滯后性表現(xiàn)愈加明顯且在相同加載端位移時，測點位移量呈減小趨勢，反映出筋材拔出的阻抗作用增大，筋-土間相互作用得到提升。見圖6。

對比同一測點處s0為10、15、20 kPa時的位移量，s0為10 kPa時，位移值最大，隨初值增加位移呈遞減趨勢，這是因為在幅值為5 kPa即界面正應力較小時，拉力沿筋材長度方向傳遞較完全且界面正應力對位移沿筋材傳遞影響較大。

2.2 界面循環(huán)正應力頻率對拉拔特性的影響

選取s0為15 kPa、An為5 kPa、fn分別為0.1、0.2、0.5 Hz，研究界面循環(huán)正應力頻率對筋材拉力與拉拔位移的影響。隨著fn增大，峰值拉拔力呈衰減趨勢且至峰值穩(wěn)定狀態(tài)時的加載端位移量略微減小，即達到峰值穩(wěn)定狀態(tài)時間提前，主要原因為隨著fn增大，上覆荷載震動周期減小，整體振動強度增加，導致土體顆粒易發(fā)生重新排列，使得筋材拔出的阻抗作用減弱，從而呈現(xiàn)出峰值拉拔力產(chǎn)生減小的情況。見圖7和圖8。

隨著測點位置距加載端距離的增加，位移量呈遞減趨勢且表現(xiàn)出明顯的滯后性。相同加載端位移時，測點位移量隨頻率增大而增加，主要原因是隨著頻率增大，填土顆粒發(fā)生更多次的重新排列組合，導致筋-土間的嵌固咬合和摩擦作用力降低，埋置土體內(nèi)部筋材拔出難度降低，其位移量也隨之增加。表明增加界面正應力作用頻率導致筋-土界面作用降低。見圖9。

2.3 界面循環(huán)正應力幅值對拉拔特性的影響

選取s0為15 kPa、fn為0.1 Hz、An分別為2.5、5、10 kPa，研究幅值對筋材拉力與拉拔位移關系的影響。相同加載端位移時，隨著An不斷增加，拉拔力減小，主要是因為界面正應力最小值隨幅值的增加而減小，由于循環(huán)速度較快，筋材拔出過程中表現(xiàn)出兩種不同狀態(tài)，即界面正應力最小值時筋材持續(xù)拔出，但在界面正應力較大時筋材不會移動，故而界面正應力最小值決定了筋材拔出時筋-土界面阻抗作用的大小，呈現(xiàn)出峰值拉拔力隨幅值增大而減小的情況。An為10 kPa時，加載端拉拔力-加載端位移關系曲線中上升階段存在波動起伏狀況，主要原因是隨著幅值增大，界面循環(huán)正應力最大差值增大，拉拔過程中不同界面正應力時土體產(chǎn)生的塑性變形程度不同，即在拉拔階段時筋-土界面接觸部位易產(chǎn)生部分土體凸起，也同樣存在土體凸起部位消失情況，從而在拉拔力方面表現(xiàn)出驟增和驟降現(xiàn)象，即筋材拉拔過程中拉拔力-加載端位移曲線的異常波動情況。見圖10和圖11。

相同加載端位移時，測點位移隨幅值的增加而增加；各測點初始移動時的加載端位移量減小且隨距加載端距離增加而愈加明顯，表明增加界面循環(huán)正應力的幅值，促進了拉力沿筋材長度方向的傳遞，表現(xiàn)出相同加載端位移時筋材各測點變形量有所增加。見圖12。

3 結論與展望

1）隨著筋-土界面循環(huán)正應力初值增加，拉拔過程中拉拔力峰值呈遞增趨勢，筋材沿長度距加載端越遠，測點位移量呈逐步遞減趨勢。

2）隨界面循環(huán)正應力頻率或幅值的增加，拉拔力峰值呈遞減趨勢，距加載端不同距離處的位移量呈遞增趨勢。

3）增加界面循環(huán)正應力頻率，導致筋-土界面相互嵌固咬合和摩擦作用力降低；因此，在高頻荷載作用下，加筋土結構的穩(wěn)定性要引起重視。

筋材水平拉拔力作用模式是影響筋-土界面特性的重要因素，未來可通過改變筋材水平拉拔力的大小和頻率，進一步深入研究筋-土相互作用特性。

參考文獻：

[1]凌天清，周濱，吳春波， 等. 筋土界面摩擦特性影響因素分析[J]. 交通運輸工程學報，2009，9（5）：7-12.

[2]萬亮，楊和平. 格柵加筋膨脹土拉拔試驗測試值的數(shù)值模擬[J]. 公路交通科技，2020，37（2）：31-39.

[3]張孟喜，陳高峰，朱引，等. H-V加筋飽和砂土性狀的三軸試驗研究[J]. 巖土力學，2010，31（5）：1345-1351.

[4]周健，王家全，孔祥利，等. 砂土顆粒與土工合成材料接觸界面細觀研究[J]. 巖土工程學報，2010，32（1）：61-67.

[5]Jensen R P，EdilT B，Bosscher P J，et al. Effect of particle shape on interface behavior of DEM-simulated granular materials[J]. International Journal of Geomechanics，2001，1（1）：1-19.

[6]Wang Z J，Jacobs F，Ziegler M. Visualization of load transfer behaviour between geogrid and sand using PFC2D[J]. Geotextiles and Geomembranes，2014，42（2）：83-92.

[7]易富，張利陽. 土工格柵拉拔試驗模型優(yōu)化及加筋特性研究[J]. 安全與環(huán)境學報，2020， 20（1）：73-81.

[8]徐超，石志龍. 循環(huán)荷載作用下筋土界面抗剪特性的試驗研究[J]. 巖土力學，2011， 32（3）：655-660.