何忠明，蔡中心，曹平，劉建華，周蓮君

(1. 長沙理工大學 交通運輸工程學院，湖南 長沙，410004；2. 中南大學 資源與安全工程學院，湖南 長沙，410083；3. 中南大學 地球科學與信息物理學院，湖南 長沙，410083)

全長黏結(jié)式錨桿在邊坡工程中廣泛應用[1-3]。以往一般采用理論分析的方法對錨桿加固工程進行分析，這種方法中將錨桿的受力和變形視為一維情況。即只承受軸向力和軸向變形，而不發(fā)生橫向變形。但在加固節(jié)理邊坡過程中，錨桿往往受到橫向剪切力的作用，一些學者指出其重要性與軸向力作用相當[4-6]。另外，由于錨桿和巖土體相互作用的復雜性[7-9]，使得理論分析模型不能完全反映錨桿的特征，因此，一些研究者采用數(shù)值模擬的方法進行研究[10-12]，但這些研究中較少考慮錨桿的橫向作用特征。在此，本文作者在前人研究的基礎(chǔ)上，采用錨桿數(shù)值計算單元，利用FLAC3D[13]建立層狀邊坡的計算模型，分析錨桿在同時受到軸向和橫向作用情況下邊坡的加固效應，以便為工程實際加固提供參考。

1 數(shù)值模型

利用拉格朗日差分軟件FLAC3D建立數(shù)值分析模型，如圖1所示。對于節(jié)理的模擬，采用低強度彈塑性夾層單元，節(jié)理傾角為40°，厚度為0.1 m；節(jié)理以外的巖體仍視為均質(zhì)體。利用自編的ANSYS-FLAC3D接口程序，按照平面應變建立數(shù)值模型，模型共36 253個單元，11 977個節(jié)點。自然坡角為10°，開挖后形成的坡角為 75°。邊界條件為下部固定約束，左右兩側(cè)水平約束，上部為自由邊界；初始應力場按自重應力考慮；計算收斂準則為不平衡力比率(節(jié)點平均內(nèi)力與最大不平衡力的比值)滿足10-7的求解要求。邊坡計算模型見圖1，計算參數(shù)見表1。錨桿結(jié)構(gòu)單元參數(shù)根據(jù)實際工程中的錨桿參數(shù)取值：錨桿長為32 m，傾角為15°，間距為2.5 m；彈性模量為200 GPa，泊松比為0.25，截面積為314 mm2，鉆孔周長為189.6 mm；切向黏結(jié)力為1.75×105N/m，切向摩擦角為30°，黏結(jié)剛度為1.0×109N/m2，法向黏結(jié)力為1.75×108N/m，法向剛度為1.0×109N/m2。邊坡和錨桿設(shè)置順序為邊開挖邊支護，錨桿布置方式如圖2所示，從上到下共布置18排錨桿。

圖1 數(shù)值模擬模型Fig.1 Numerical simulation model

圖2 錨桿布置圖Fig.2 Location of bolts

表1 計算參數(shù)Table 1 Calculation parameters

為了分析邊坡巖土體的變形特征，設(shè)置相應位移監(jiān)測點，監(jiān)測點具體布置位置如圖3所示。沿坡頂向下每隔 5 m設(shè) 1點，共布設(shè) 10個監(jiān)測點，分別為P01～P10，并在每個記錄點位置沿自然坡傾向的方向，向坡體內(nèi)側(cè)每隔2 m布設(shè)1點，每條監(jiān)測線共設(shè)26個記錄點。

圖3 監(jiān)測點位置Fig.3 Location of monitoring points

2 數(shù)值計算結(jié)果與討論

2.1 監(jiān)測點水平位移

圖4所示為邊坡在錨固后的水平位移分布情況。從圖4可以看出：邊坡的破壞形式為明顯的直線型滑動破壞，滑動面位于第3條結(jié)構(gòu)面處，并且滑塊的位移從上到下逐漸增大，水平位移最大處位于第3條結(jié)構(gòu)面處，即剪出口位置(見圖5)，并且剪出口位置處布置有監(jiān)測點P09。層狀邊坡錨固后監(jiān)測點的水平位移如圖6所示，從圖6可以看出：監(jiān)測點從上到下的位移曲線沿自然坡的變化趨勢從均勻變化過渡到尖點突變變化形態(tài)，其中：P01～P04的位移曲線比較均勻，但從P05開始曲線形態(tài)發(fā)生變化，出現(xiàn)多個臺階。這是由于此時監(jiān)測線經(jīng)過節(jié)理面，位移傳遞受到阻隔；另外，P07和P09監(jiān)測線位移突變現(xiàn)象明顯，這是由于這些監(jiān)測點離邊坡剪出口位置較近，受到的影響較大。以上說明節(jié)理邊坡錨桿加固后，沿節(jié)理面仍然會發(fā)生一定的位錯變形，表明錨桿的橫向抗變形能力較弱。各監(jiān)測點位移沿自然坡傾向位置逐漸減小，說明越往巖體內(nèi)部受到的開挖擾動越小，邊坡的水平位移也越小，這與實際情況相符。

2.2 監(jiān)測點總位移

圖4 錨固后邊坡水平位移云圖Fig.5 Horizontal displacement of slope after reinforcement

圖5 邊坡剪應變增量云圖Fig.5 Shear strain increment of slope

圖6 錨固后監(jiān)測點水平位移Fig.6 Horizontal displacement of monitoring points after reinforcement

圖7 錨固后監(jiān)測點總位移Fig.7 Total displacement of monitoring points after reinforcement

監(jiān)測點總位移在逆自然坡傾向方向上的變化曲線如圖7所示。從圖7可見：同一測試點的總位移與開挖坡體相比有所減?。粶y試點的總位移被節(jié)理面分為3個區(qū)域，每個區(qū)域內(nèi)測試點的總位移從坡面至坡里側(cè)方向逐漸減小，但同一區(qū)域內(nèi)的測試點總位移變化不大。這說明錨桿加固邊坡后，測試點的平動位移仍然占其總位移的較大部分(監(jiān)測點的位移等于巖體的平動位移與變形之和)，即巖體沿節(jié)理面還會由一定的位錯變形，其加固屬于柔性加固方式，但與開挖后測試點的總位移曲線相比，節(jié)理兩側(cè)測試點的總位移突變值變小，說明錨桿對節(jié)理的上下盤巖體起到了有效的拉結(jié)作用，抑制了兩盤之間的較大變形位移。但其對滑動變形的控制作用是被動的，即只有上下盤巖體發(fā)生一定的位錯后，錨桿的拉結(jié)力才能發(fā)揮出來。

2.3 錨桿軸力

圖8 錨桿軸力分布Fig.8 Distribution of bolt axial stress

圖9 錨桿軸力分布云圖Fig.9 Distribution of bolt axial stress

選擇第1，5，10，15和18排錨桿作為軸力監(jiān)測的典型錨桿，圖8所示為其軸力分布曲線(其中，拉為正，壓為負)。從圖8可以看出：由于節(jié)理的存在，全長黏結(jié)式錨桿的軸力分布為多峰值曲線，峰值均出現(xiàn)在節(jié)理面位置，這是由于錨桿跟隨節(jié)理面的錯動產(chǎn)生了較大軸向變形。位于節(jié)理面之間的錨桿曲線呈現(xiàn)先下降再上升的趨勢，如通過第1條節(jié)理后，錨桿的軸力迅速衰減，也符合以上理論分析的結(jié)果。但是經(jīng)過一段距離后，由于第2條節(jié)理的出現(xiàn)，其軸力曲線又迅速攀升至峰值。并且錨桿軸力的最大峰值處可表征該節(jié)理面的錯動最大，對于邊坡的穩(wěn)定性影響最大，如圖9所示。同時，錨桿軸力最大位置可表征滑動面位置；在以往工程設(shè)計中，都假設(shè)錨桿軸力分布模式為均勻分布，但是從圖9可看出：錨桿的峰值軸力遠大于平均軸力。因此，按照均勻分布模式計算是不合理的。

3 結(jié)論

(1) 邊坡的破壞形式為明顯的直線型滑動破壞，滑塊的位移從上到下逐漸增大，水平位移最大處位于剪出口位置；監(jiān)測點從上到下的位移曲線沿自然坡的變化趨勢從均勻變化過渡到尖點突變變化形態(tài)，監(jiān)測線經(jīng)過節(jié)理面，位移傳遞收到阻隔。

(2) 節(jié)理邊坡錨桿加固后，沿節(jié)理面仍然會發(fā)生一定的位錯變形，表明錨桿的橫向抗變形能力較弱。各監(jiān)測點位移沿自然坡傾向位置逐漸減小，越往巖體內(nèi)部受到的開挖擾動越小，邊坡的水平位移也越小，符合實際情況。

(3) 由于存在節(jié)理，全長黏結(jié)式錨桿的軸力分布為多峰值曲線，峰值均出現(xiàn)在節(jié)理面位置；位于節(jié)理面之間的錨桿曲線呈現(xiàn)先下降再上升的趨勢，如通過第1條節(jié)理后，錨桿的軸力迅速衰減。但是經(jīng)過一段距離后，由于第2條節(jié)理的出現(xiàn)，其軸力曲線又迅速攀升至峰值；同時，錨桿軸力最大位置可表征滑動面位置。

[1] 楊慶, 朱訓國, 欒茂田. 全長注漿巖石錨桿雙曲線模型的建立及錨固效應的參數(shù)分析[J]. 巖石力學與工程學報, 2007,26(4): 692-698.YANG Qing, ZHU Xun-guo, LUAN Mao-tian. Development of hyperbolic model for fully grouting rock bolt and parameters analysis for anchoring effect[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(4): 692-698.

[2] 鄒金鋒, 李亮, 楊小禮, 等. 基于損傷理論的全長式錨桿荷載傳遞機理研究[J]. 鐵道學報, 2007, 29(6): 84-88.ZOU Jin-feng, LI Liang, YANG Xiao-li, et al. Study on load transfer mechanism for span-type anchor based on the damage theory[J]. Journal of the China Railway Society, 2007, 29(6):84-88.

[3] 曹國金, 姜弘道, 熊紅梅. 一種確定拉力型錨桿支護長度的方法[J]. 巖石力學與工程學報, 2003, 22(7): 1141-1145.CAO Guo-jin, JIANG Hong-dao, XIONG Hong-mei.Calculation method of supportin length for stretched bolts[J].Chinese Journal of Rock Mechanics and Engineering, 2003,22(7): 1141-1145.

[4] 楊松林, 徐衛(wèi)亞, 黃啟平. 節(jié)理剪切過程中錨桿的變形分析[J]. 巖石力學與工程學報, 2004, 23(9): 3268-3273.YANG Song-lin, XU Wei-ya, HUANG Qi-ping. Analysis on the bolt deformation as result of joint shear displacement[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(9):3268-3273.

[5] 劉波, 李先煒, 陶龍光. 錨拉支架中錨桿橫向效應分析[J]. 巖土工程學報, 1998, 20(4): 36-39.LIU Bo, LI Xian-wei, TAO Long-guang. Analysis on lateral behavior of bolts in roof truss system[J]. Chinese Journal of Geotechnical Engineering, 1998, 20(4): 36-39.

[6] LI C, Stillborg B. Analytical models for rock bolts[J]. Int J Rock Mech Sci and Geomech Abstr, 1999, 36(8): 1013-1029.

[7] 彭文祥, 趙明華, 袁海平, 等. 基于拉格朗日差分法的全長注漿錨桿支護參數(shù)優(yōu)化[J]. 中南大學學報: 自然科學版, 2006,37(5): 1002-1007.PENG Wen-xiang, ZHAO Ming-hua, YUAN Hai-ping, et al.Parameters analysis of grouted bolts by Lagrangian difference method[J]. Journal of Central South University: Science and Technology, 2006, 37(5): 1002-1007.

[8] Samuel Frimpong, LI Ying. Stress loading of the cable shovel boom under in-situ digging conditions[J]. Engineering Failure Analysis, 2007, 14(4): 702-715.

[9] Soneji B B, Jangid R S. Influence of soil-structure interaction on the response of seismically isolated cable-stayed bridge[J]. Soil Dynamics and Earthquake Engineering, 2008, 28(4): 245-257.

[10] Guo R, Thompson P. Influences of changes in mechanical properties of an overcored sample on the far-field stress calculation[J]. International Journal of Rock Mechanics &Mining Sciences, 2002(39): 1153-1166.

[11] 林杭, 曹平. 錨桿長度對邊坡穩(wěn)定性影響的數(shù)值分析[J]. 巖土工程學報, 2009, 31(3): 470-474.LIN Hang, CAO Ping. Numerical analysis for the effect of cable length to the stability of slope[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(3): 470-474.

[12] 林杭, 曹平, 周正義. FLAC3D模擬全長注漿錨桿的作用效果[J]. 巖土力學, 2005, 26(s2): 167-170.LIN Hang, CAO Ping, ZHOU Zheng-yi. Simulating effects of wholly grouted anchor rods with FLAC3D[J]. Rock and Soil Mechanics, 2005, 26(s2): 167-170.

[13] Itasca Consulting Group. User’s Guide[Z]. Minnesota: Itasce consulting group, 2002.