黃軍朋，張景科，王南，李卷強，郭青林，趙林毅

(1. 蘭州大學 西部災(zāi)害與環(huán)境力學教育部重點實驗室；土木工程與力學學院，蘭州 730000；2. 敦煌研究院，甘肅 敦煌 736200；3.國家古代壁畫與土遺址保護工程技術(shù)研究中心，甘肅 敦煌 736200)

錨固技術(shù)是解決土遺址變形開裂病害的主要方法之一。土遺址加固工程中，錨固桿材先后出現(xiàn)過薄壁鋼管[1]、白蠟桿[2]、楠竹[3]、楠竹加筋復合錨桿[4]等。多年實踐證明，金屬桿材及水泥系錨固漿液與遺址體的兼容性較差，不符合中國文物保護的理念和原則?；趥鹘y(tǒng)材料與工藝的科學挖掘，作為西北干旱半干旱地區(qū)夯筑土遺址中常用的加筋體，天然木材成為較理想的錨固桿材。

目前，土遺址加固中木錨桿主要為天然白蠟桿，區(qū)別于其他巖土錨固領(lǐng)域中規(guī)則、均質(zhì)、順直的金屬桿材，其具有通長變徑、表面形態(tài)不均一、近似圓臺狀、強度各向異性等特點。由于土遺址用木錨桿本身結(jié)構(gòu)變徑的特殊性，其剪應(yīng)力傳遞與分布規(guī)律、破壞形式等錨固特性不同于其他巖土錨固領(lǐng)域中鋼筋、鋼絞線等標準規(guī)則桿體。對于標準規(guī)則錨桿研究較為深入，如Phillips[5]提出剪應(yīng)力按指數(shù)分布；蔣忠信[6]提出剪應(yīng)力按三參數(shù)高斯曲線分布；張季如等[7]提出剪應(yīng)力按雙曲函數(shù)分布等。以上錨固理論在實踐中均已得到映證。而土遺址中變徑木錨桿錨固機制的研究遠滯后于工程應(yīng)用，在實際應(yīng)用中多采用剪應(yīng)力均勻分布的理想形式設(shè)計，這與試驗結(jié)果不符[8]。因此，有必要對土遺址用變徑木錨桿的力學傳遞機制、界面剪應(yīng)力分布規(guī)律和錨固效果的影響因素等進行更為深入的研究。

夯筑土遺址具有夯筑分層的土體結(jié)構(gòu)致使土體力學特征不均一，土遺址建造具有明顯的地域特征，在遺址本體上進行錨固試驗不符合土遺址保護的理念，以上問題極大地制約土遺址錨固機理的研究。此外，室內(nèi)外拉拔試驗受界面應(yīng)變測試元件限制，實驗值較為離散。然而，數(shù)值方法在錨固機制及錨固效果研究方面具有顯著的優(yōu)勢。在錨固機制研究方面，戰(zhàn)玉寶等[9]、龐有師等[10]采用有限元法，李英勇等[11]、葉紅等[12]采用有限差分法模擬研究錨索錨固段剪應(yīng)力分布和傳遞規(guī)律并進行錨固參數(shù)分析；在錨固效果研究方面，賈金青等[13]、林杭等[14]、Murphy等[15]運用FLAC3D分別模擬研究錨索支護深基坑、邊坡、礦井等錨固工程的錨固效果，王連國等[16]采用有限元法模擬錨注支護對深部軟巖巷道的錨固效果。以上數(shù)值模擬研究主要是標準規(guī)則桿材的錨固機制及其應(yīng)用于基坑、邊坡、礦井、巷道等工程，而對土遺址加固領(lǐng)域中變徑桿材數(shù)值模擬研究還是空白。

鑒于此，本文基于前期室內(nèi)試驗物理模型的研究成果[8]，利用FLAC3D模擬研究變徑木錨桿錨固系統(tǒng)的錨固性能、界面力學傳遞規(guī)律，并對室內(nèi)試驗結(jié)果與數(shù)值模擬計算結(jié)果進行對比和分析，進而研究錨孔直徑、錨桿直徑和錨固長度對變徑木錨桿錨固效果的影響。

1 錨固試驗案例簡介

2 變徑木錨桿錨固機理的數(shù)值模擬

2.1 計算模型建立及參數(shù)

錨固體的計算模型選取層高0.1 m、半徑14.5 cm、共5層夯土(夯層1～5)的圓柱體，木錨桿布置在模型的中心(圖3)。模型由11 040個單元，13 034個節(jié)點組成，變徑木錨桿采用實體單元模擬，直徑由上至下線性增大，錨固長度0.5 m。夯土采用8個節(jié)點的六面體單元模擬，網(wǎng)格劃分以錨固體為中心向四周呈放射狀由密變疏(圖4)?？卓诘藉^固末端為z軸負方向，孔口位置為坐標原點。夯層1上部為內(nèi)半徑4.5 cm、外半徑14.5 cm、高2 cm的圓形護層鋼板。計算模型的邊界條件為底部和側(cè)面均施加法向約束(圖3)。

圖3 模型及邊界條件Fig.3 Boundary condition

圖4 網(wǎng)格劃分Fig.4 Mesh of geometry

選擇文獻[8]中M3木錨桿錨固系統(tǒng)作為研究對象，室內(nèi)拉拔試驗對M3錨固系統(tǒng)進行單級2 kN加載以獲取極限抗拔力。M3桿長1 m，錨固長度0.5 m，起始端半徑1.6 cm，末端半徑2.6 cm。由于木錨桿天然結(jié)構(gòu)變徑的特殊性，其直徑變化非線性，按理想條件假定直徑線性變化， 則錨固長度0.5 m時，錨固段頂端半徑2.1 cm，末端半徑2.6 cm。錨固漿體為基于水硬石灰和石英砂的漿液，水硬石灰、石英砂按照質(zhì)量1∶1配比，水灰比為0.33。由于遺址土具有夯筑分層的性質(zhì)，遺址土表層土風化程度不同，土樣顆粒不均勻，底層夯筑次數(shù)多于頂層，均致使同一遺址土不同夯層土樣的物理力學性質(zhì)有差異[17]，總體上是下部夯土強度高于上部夯土強度。實測M3錨固系統(tǒng)材料的物理力學參數(shù)見表1，界面力學參數(shù)見表2。

表1 材料的物理力學參數(shù)[8,18]Table 1 Physico-mechanical parameters of materials

表2界面力學參數(shù)
Table2Interfacialmechanicalparameters

接觸面ks/(Pa·m-1)kn/(Pa·m-1)с/kPaφ/(°)σt/MPa接觸面16.579×10106.579×1010500.024.00.5接觸面22.500×1082.500×108212.030.00.4

2.2 屈服破壞準則

采用目前巖土工程中應(yīng)用最為廣泛的非關(guān)聯(lián)流動剪切屈服Mohr-Coulomb準則，其表達式為[19]

式中：I1、J2分別為應(yīng)力張量第1不變量和應(yīng)力偏量第2不變量；θσ為應(yīng)力羅德角；c、φ分別為粘聚力和內(nèi)摩擦角。

2.3 數(shù)值模擬結(jié)果分析

圖5 荷載位移關(guān)系曲線Fig.5 Load displacement curves

2.3.2 界面剪應(yīng)力分布規(guī)律 綜合FLAC3D數(shù)值計算下木錨桿漿體界面剪應(yīng)力分布云圖(圖6)和木錨桿漿體界面剪應(yīng)力沿錨桿軸向分布圖(圖7)可知：剪應(yīng)力沿錨固段呈不均勻分布，在木錨桿頂端和末端的0.1 m范圍內(nèi)量值較大，中部量值相對較小。其剪應(yīng)力分布不同于其他錨固工程中規(guī)則桿體剪應(yīng)力分布的理論模型[5-7]。當荷載P≤6 kN時，界面剪應(yīng)力呈U型曲線，在錨桿頂端的量值最大，且沿錨桿軸向向中部呈減小的趨勢，在距錨桿錨固體頂端0.4 m處剪應(yīng)力呈增大的趨勢，此時錨固界面處于彈性黏結(jié)狀態(tài)。當6 kN 8 kN時，錨固段進入滑移破壞階段，頂端和末端界面剪應(yīng)力繼續(xù)增大，而中間錨固段界面剪應(yīng)力增幅較小。

圖6 不同荷載下界面剪應(yīng)力分布云圖Fig.6 Distributions of the interfacial shear stress under different loads(unit: Pa)

圖7 界面剪應(yīng)力分布圖Fig.7 Distributions of the interfacial shear stres

鑒于木錨桿近似圓臺狀，其頂端直徑小于末端直徑，隨著荷載及位移的增加，變徑誘發(fā)的剪脹作用愈發(fā)明顯，在錨桿末端0.1 m范圍內(nèi)界面剪應(yīng)力呈現(xiàn)增大的趨勢。當上部界面粘結(jié)力不足以抵抗拉拔荷載而發(fā)生剪切破壞時，粘結(jié)力主要由錨固段末端提供，木錨桿變徑的特點在一定程度上提高了錨固系統(tǒng)的抗拔力。

2.3.3 土體漿體應(yīng)力場、位移場 不同荷載作用下土體漿體的應(yīng)力云圖(圖8)表明：漿體土體所受應(yīng)力以木錨桿為中心呈軸對稱分布；土體上部為壓應(yīng)力區(qū)，下部為拉應(yīng)力區(qū)；土體中拉、壓應(yīng)力區(qū)分界面(szz=0 Pa等值線)可用一拋物面來描述，且拉拔荷載越大，拋物分界面越靠近錨固段末端。木錨桿受拉時，上部及中下部漿體處于受拉狀態(tài)；由于木錨桿直徑漸大，荷載在木錨桿末端對漿體施加壓應(yīng)力，漿體與周圍土體產(chǎn)生剪切抗力，以此提供錨桿所需的承載力，說明變徑木錨桿同時具有拉力型和壓力型錨桿的特征。

圖8 土體漿體應(yīng)力云圖Fig.8 Stress distributions of soil and slurr

不同荷載作用下漿體土體的位移云圖(圖9)表明：土體漿體的位移以木錨桿為中心呈軸對稱分布，且位移隨著荷載的增加而增大。由于液壓油缸對上部鋼板的作用力傳遞到上部土體，在土體下部，木錨桿末端直徑較大向四周擠壓漿體土體，使得土體位移等值線呈拋物線狀，即在縱向兩端土體位移小于中部土體位移。

2.3.4 剪脹現(xiàn)象 數(shù)值模擬結(jié)果表明變徑木錨桿對漿體土體具有剪脹作用。由圖9可知：由于木錨桿本身變徑的結(jié)構(gòu)特點，錨桿末端直徑相對較大，而向四周擠壓土體使土體位移增大，土體位移等值線呈拋物線狀，在縱向兩端位移小于中間位移；越靠近木錨桿漿體土體橫向位移越大，徑向擠壓土體，體現(xiàn)變徑木錨桿徑向的剪脹作用，這與拉拔試驗中剪脹作用導致PVC管破裂一致[8]。

圖9 土體漿體位移云圖Fig.9 Displacement distributions of soil and slurr

3 錨固參數(shù)分析

木錨桿極限抗拔力是土遺址錨固研究的基礎(chǔ)。變徑木錨桿錨固作用機理復雜，影響錨固效果的因素較多，其中木錨桿錨固系統(tǒng)本身的結(jié)構(gòu)特征是影響錨固效果的主要因素。

3.1 錨孔直徑

圖10 不同錨孔直徑荷載位移曲線Fig.10 Load displacement curves of different

不同錨孔直徑界面剪應(yīng)力分布圖(圖11)表明：1)界面剪應(yīng)力呈不均勻分布，木錨桿兩端界面剪應(yīng)力顯著大于中部界面剪應(yīng)力，末端界面剪應(yīng)力增大主要是木錨桿本身變徑的結(jié)構(gòu)特征所引起。2)在荷載一定時，界面剪應(yīng)力與錨孔直徑密切相關(guān)。錨孔直徑越小，剪應(yīng)力峰值越大；錨孔直徑越大，剪應(yīng)力峰值越小，剪應(yīng)力分布曲線越平緩。

圖11 不同孔徑界面剪應(yīng)力分布圖Fig.11 Distributions of the interfacial shear stress under different anchor hole diameter

3.2 錨桿直徑

圖12 不同錨桿直徑荷載位移曲線Fig.12 Load displacement curves of different anchor diameters

不同桿徑界面剪應(yīng)力分布圖(圖13)表明：1)界面剪應(yīng)力沿錨桿軸向呈不均勻分布，木錨桿直徑逐漸增大而引起錨桿末端界面剪應(yīng)力呈增大的趨勢；2)桿徑越小其頂端界面剪應(yīng)力隨著拉拔荷載增大首先達到界面容許粘結(jié)力而產(chǎn)生破壞；3)隨著桿徑及荷載的增大，木錨桿頂端界面剪應(yīng)力峰值減小而末端界面剪應(yīng)力增大，說明抗拔力由錨桿頂端界面粘結(jié)力提供轉(zhuǎn)向由錨桿末端界面粘結(jié)力提供，木錨桿變徑的特點能充分發(fā)揮漿體的抗剪強度。

圖13 不同桿徑界面剪應(yīng)力分布圖Fig.13 Distributions of the interfacial shear stress under different anchor diameter

3.3 錨固長度

圖14 不同錨固長度荷載位移曲線Fig.14 Load displacement curves of different anchorage lengths

不同錨固長度界面剪應(yīng)力的分布圖(圖15)表明：1)界面剪應(yīng)力沿著錨桿軸向分布不均，主要集中在木錨桿頂端和末端的0.1 m范圍內(nèi)，且頂端界面剪應(yīng)力大于末端界面剪應(yīng)力，木錨桿變徑的特征引起木錨桿末端界面剪應(yīng)力呈增大的趨勢；2)錨固長度越小，木錨桿末端的直徑相對較小，隨著荷載的增大，首先達到界面容許粘結(jié)力，從而界面剪應(yīng)力相對較大，越先發(fā)生滑移破壞；3)隨著荷載及錨固長度的增大，錨固段中部界面剪應(yīng)力增值較小，剪應(yīng)力的分布特征沒有太大的改變。

圖15 不同錨固長度界面剪應(yīng)力分布圖Fig.15 Distributions of the interfacial shear stress under different anchorage length

3.4 極限抗拔力計算

綜合錨固參數(shù)分析結(jié)果可知：木錨桿抗拔力隨錨孔直徑φ、錨固長度h增大呈線性增長。對極限抗拔力Pu與錨孔直徑、錨固長度的乘積(φh)采用二次多項式擬合(圖16)，可得木錨桿錨固系統(tǒng)極限抗拔力計算式為

Pu=1 593.69(φh)2+156.87φh

圖16 極限抗拔力擬合曲線Fig.16 Fitted carve for the ultimate pullout capacit

從式中可看出，木錨桿極限抗拔力隨著錨孔直徑與錨固長度的乘積(φh)增加而增大。而實際工程中，錨固長度增加至一定長度后，木錨桿的抗拔力增加甚少。鑒于文物的特殊性，錨孔直徑、錨固長度應(yīng)根據(jù)工程實際來選擇。

4 結(jié)論

2)木錨桿漿體界面剪應(yīng)力沿錨固段分布不均，主要集中在錨固段頂端和末端的0.1 m范圍內(nèi)，末端界面剪應(yīng)力呈增大的趨勢主要與木錨桿變徑的結(jié)構(gòu)特征有關(guān)，其變徑的特征在一定程度上提高了抗拔力。

3)變徑木錨桿同時具有拉力型和壓力型錨桿的特征，且其徑向具有剪脹作用。

4)模擬試驗所考察的3個參數(shù)中，錨孔直徑、錨固長度對木錨桿抗拔力有顯著影響，而錨桿直徑對抗拔力的影響較小。木錨桿抗拔力隨著錨孔直徑、錨固長度增加而提高。土遺址中小體量危巖體的加固，適當增加錨固長度可以有效提高抗拔力。綜合模擬試驗錨固參數(shù)分析，提出木錨桿極限抗拔力建議計算公式。

[1] 李最雄，王旭東.古代土建筑遺址保護加固研究的新進展[J].敦煌研究，1997(4)：167-172.

LI Z X，WANG X D. The new progress of protection and reinforcement for the ancient earthen architecture site [J]. Dunhuang Research，1997(4)：167-172. (in Chinese)

[2] 李最雄. 絲綢之路古遺址保護[M]. 北京：科學出版社，2010.

LI Z X. Conservation of ancient sites on the Silk Road [M]. Beijing：Science Press，2010. (in Chinese)

[3] 孫滿利，王旭東，李最雄，等.木質(zhì)錨桿加固生土遺址研究[J].巖土工程學報，2006，28(12)：2156-2159.

SUN M L，WANG X D，LI Z X，et al. Study on immature earthen sites reinforced with wood anchor [J]. Chinese Journal of Geotechnical Engineering，2006，28(12)： 2156-2159. (in Chinese)

[4] 張虎元，王曉東，王旭東，等.楠竹加筋復合錨桿內(nèi)部界面黏結(jié)滑移模型[J].巖土力學，2011，32(3)：789-796.

ZHANG H Y，WANG X D，WANG X D，et al. Bond-slip model for bamboo-steel cable composite anchor [J]. Rock and Soil Mechanics，2011，32(3)：789-796. (in Chinese)

[5] PHILLIPS S H E. Factors affecting the design of anchorages in rock [R]. London：Cementation Research Ltd.，1970：1-10.

[6] 蔣忠信.拉力型錨索錨固段剪應(yīng)力分布的高斯曲線模式[J].巖土工程學報，2001，23(6)：696-699.

JIANG Z X. A gauss curve model on shear stress along anchoring section of anchoring rope of extensional force type [J]. Chinese Journal of Geotechnical Engineering，2001，23(6)：696-699. (in Chinese)

[7] 張季如，唐保付.錨桿荷載傳遞機理分析的雙曲函數(shù)模型[J].巖土工程學報，2002，24(2)：188-192.

ZHANG J R，TANG B F. Hyperbolic function model to analyze load transfer mechanism on bolts [J]. Chinese Journal of Geotechnical Engineering，2002，24(2)：188-192. (in Chinese)

[8] 張景科，陳鵬飛，諶文武，等.夯土-水硬石灰與石英砂漿液-木錨桿錨固系統(tǒng)性能[J].土木建筑與環(huán)境工程，2015，37(3)：86-93.

ZHANG J K，CHEN P F，CHEN W W，et al. Anchor performance of wood bolt fully grouted by natural hydraulic lime and quartz sand slurry [J]. Journal of Civil，Architectural & Environmental Engineering，2015，37(3)：86-93. (in Chinese)

[9] 戰(zhàn)玉寶，畢宣可，尤春安. 預(yù)應(yīng)力錨索錨固段應(yīng)力分布影響因素分析[J].土木工程學報，2007，40(6)：49-53.

ZHAN Y B，BI X K，YOU C A. Analysis of factors influencing the stress distribution in prestressed cables [J]. China Civil Engineering Journal，2007，40(6)：49-53. (in Chinese)

[10] 龐有師，劉漢龍，陳育民.可回收式錨桿拉拔試驗的數(shù)值模擬與影響因素分析[J].解放軍理工大學學報(自然科學版)，2009，10(2)：170-174.

PANG Y S，LIU H L，CHEN Y M. Numerical simulation of removable anchor pullout test and influence factors [J]. Journal of PLA University of Science and Technology (Natural Science Edition)，2009，10(2)：170-174. (in Chinese)

[11] 李英勇，張頂立，王松根.預(yù)應(yīng)力錨索錨固作用機理的數(shù)值模擬研究[J].巖土力學，2006，27(Sup)：921-925.

LI Y Y，ZHANG D L，WANG S G. Research on numerical simulation of anchoring mechanism of prestressed anchor cables [J]. Rock and Soil Mechanics，2006，27(Sup)：921-925. (in Chinese)

[12] 葉紅，盛建龍，陳燕平.壓力型錨索孔周邊巖體應(yīng)力的數(shù)值模擬分析[J].武漢理工大學學報，2012，34(7)：115-118.

YE H，SHENG J L，CHEN Y P. Numerical simulation analysis of rock stress surrounding pressure type cable hole [J]. Journal of Wuhan University of Technology，2012，34(7)：115-118. (in Chinese)

[13] 賈金青，鄭衛(wèi)鋒，陳國周.預(yù)應(yīng)力錨桿柔性支護技術(shù)的數(shù)值分析[J].巖石力學與工程學報，2005，24(21)： 3978-3982.

JIA J Q，ZHENG W F，CHEN G Z. Numerical analysis of prestressed anchor flexible retaining method [J]. Chinese Journal of Rock Mechanics and Engineering，2005，24(21)：3978-3982. (in Chinese)

[14] 林杭，曹平，周正義. FLAC3D模擬全長注漿錨桿的作用效果[J]. 巖土力學，2005，26(Sup)：167-170.

LIN H，CAO P，ZHOU Z Y. Simulating effects of wholly grouted anchor rods with FLAC3D [J]. Rock and Soil Mechanics，2005，26(Sup)：167-170. (in Chinese)

[15] MURPHY M M，ESTERHUIZEN G S，TULU I B. Lateral reinforcement of fully grouted roof bolts in a FLAC3D simulated coal mine entry [C] //48th US Rock Mechanics / Geomechanics Symposium,American Rock Mechanics Association (ARMA)，Minneapolis，Minnesota,2014：1945-1953.

[16] 王連國，繆協(xié)興，董健濤，等.深部軟巖巷道錨注支護數(shù)值模擬研究[J].巖土力學，2005，26(6)：983-985.

WANG L G，MIAO X X，DONG J T，et al. Numerical simulation research of bolt-grouting support in deep soft roadway [J]. Rock and Soil Mechanics，2005，26(6)：983-985. (in Chinese)

[17] 趙海英，李最雄，韓文峰，等.西北干旱區(qū)土遺址的主要病害及成因[J].巖石力學與工程學報，2003，22(Sup2)：2875-2880.

ZHAO H Y，LI Z X，HAN W F，et al. Main diseases and their causes of earthen ruins in arid region of northwestern China [J]. Chinese Journal of Rock Mechanics and Engineering，2003，22(Sup2)：2875-2880. (in Chinese)

[18] 張景科，郭青林，李最雄，等.土遺址錨固機理初探[M].蘭州：蘭州大學出版社，2014.

ZHANG J K，GUO Q L，LI Z X，et al. Preliminary study on anchoring mechanism of earthen sites [M]. Lanzhou：Lanzhou University Press，2014. (in Chinese)

[19] Itasca Consulting Group Inc. FLAC3D (fast Lagrangian analysis of continua in three-dimensions)，version2.1，User's Mannual [R]. Minneapolis：Itasca Consulting Group Inc.，2002.