陳 賓， 張 亮， 劉正才

(1.湘潭大學(xué) 能源工程學(xué)院，湖南 湘潭 411105；2.湘潭大學(xué) 土木工程與力學(xué)學(xué)院，湖南 湘潭 411105)

基于應(yīng)變軟化模型的樁錨支護過程穩(wěn)定性分析*

陳 賓1*， 張 亮2， 劉正才2

(1.湘潭大學(xué) 能源工程學(xué)院，湖南 湘潭 411105；2.湘潭大學(xué) 土木工程與力學(xué)學(xué)院，湖南 湘潭 411105)

為了研究深基坑在樁錨支護過程中的穩(wěn)定性變化規(guī)律，利用FLAC3D軟件對深基坑進行數(shù)值模擬，采用應(yīng)變軟化模型對基坑在開挖過程中的支護樁變形、錨索受力、穩(wěn)定性以及塑性區(qū)變化進行計算，結(jié)合實測數(shù)據(jù)進行計算分析，并與Mohr-Coulomb模型進行對比.研究結(jié)果表明：FLAC3D應(yīng)變軟化模型數(shù)值模擬結(jié)果與實際監(jiān)測結(jié)果基本吻合，且優(yōu)于Mohr-Coulomb模型，并能較好地反映軟巖土體的變形特性；隨著基坑開挖的進行，支護樁身最大水平位移點位置隨著開挖步的不斷增加呈現(xiàn)出逐漸下降趨勢；在基坑開挖過程中，錨索錨固力逐漸增大，并沿索體呈現(xiàn)出先增大后減小的變化趨勢，且最大值出現(xiàn)在剪切滑動面位置處；開挖過程中，基坑張拉、剪切塑性區(qū)域主要集中分布在開挖面及開挖面3～7 m深度內(nèi).

應(yīng)變軟化；樁錨支護；變形特征；穩(wěn)定性；塑性區(qū)

基坑開挖過程中，由于巖土體的應(yīng)力迅速釋放，使巖土體產(chǎn)生一定位移，向基坑內(nèi)側(cè)移動，從而影響了基坑穩(wěn)定性.為提高基坑穩(wěn)定性，需采取一定支護措施，而樁錨支護結(jié)構(gòu)因其能夠合理發(fā)揮巖土體自身的承載能力，且能充分發(fā)揮樁錨支護體系與巖土體之間的黏結(jié)作用，被廣泛應(yīng)用于深基坑支護工程中[1～8].目前，對于樁錨支護的深基坑穩(wěn)定性研究，一般采用簡化模型，而在實際工程中，由于樁錨結(jié)構(gòu)與巖土體之間涉及到復(fù)雜的力學(xué)機制，若采用簡化模型對樁錨支護結(jié)構(gòu)進行分析往往不能有效地反映巖土體與樁錨支護結(jié)構(gòu)之間作用的真實狀態(tài)，并且無法有效合理的對基坑支護進行設(shè)計，從而造成工程浪費.此外，在開挖過程中，許多巖土體顯現(xiàn)出開挖變形的非線性與巖土體力學(xué)強度非線性的特性，且在軟巖土體的開挖中尤為明顯，其原因主要是巖土體在應(yīng)力達到強度峰值之后，巖體強度隨著變形繼續(xù)增加，會逐漸降低到一個較低值，這種由于巖體變形所引起的物理力學(xué)性能劣化的行為稱為“應(yīng)變軟化”[2～5].近年來，隨著計算機的發(fā)展，推動了數(shù)值分析方法在巖土工程中的廣泛應(yīng)用，如蔡海波[6]等利用規(guī)范分析方法以及數(shù)值模擬分析手段，計算并分析了深基坑樁錨支護體系的開挖變形與受力情況，驗證樁錨支護結(jié)構(gòu)的合理性.劉自由[7]基于應(yīng)變軟化模型，分析樁錨支護參數(shù)對基坑穩(wěn)定性的影響.鄒志強[8]等將數(shù)值模擬結(jié)果和實際監(jiān)測數(shù)據(jù)進行對比分析，得出應(yīng)變軟化模型的適應(yīng)性.本文基于以上研究成果，以某深基坑工程為背景，利用FLAC3D軟件，采用應(yīng)變軟化模型分析開挖過程中基坑在樁錨支護過程中的變形、錨索受力、穩(wěn)定性變化情況，并與Mohr-Coulomb模型進行對比，結(jié)合實測數(shù)據(jù)進行分析，探討基坑在開挖過程中的穩(wěn)定性變化規(guī)律.

1 工程概況

湖南省長沙市公安局高新區(qū)分局業(yè)務(wù)辦公用房深基坑工程，位于高新區(qū)長川路以北，東方紅路以東.場地西邊較平整，地勢總體呈東北高西南低，標(biāo)高為54.71～65.30 m.項目場地設(shè)計標(biāo)高56.24 m，基坑底設(shè)計標(biāo)高44.96 m.場地平整開挖將在北側(cè)及東側(cè)形成4～12 m高邊坡，且大部分區(qū)域，下部基坑與邊坡以組合形式出現(xiàn)，局部開挖高度超過14.4 m.

本文主要研究深基坑北側(cè)的樁錨支護在基坑開挖中的變形.場地北側(cè)基坑與邊坡呈組合形式出現(xiàn)，幾乎無放坡空間，需垂直開挖，下部基坑設(shè)計開挖深度10.1 m，上部邊坡設(shè)計開挖4.3 m，支護高度為14.4 m，該段邊坡坡頂為新建7層居民住宅，建筑基礎(chǔ)距離紅線最短10.2 m，基坑(邊坡)一旦變形過大或失穩(wěn)，后果十分嚴(yán)重.根據(jù)基坑支護設(shè)計，北側(cè)基坑支護形式采用支護樁+預(yù)應(yīng)力錨索支護.基坑北側(cè)長43.5 m，基坑支護高度16 m，設(shè)計樁徑1 m，樁間距2 m，嵌固深度5.0 m，采用三排錨索，錨索豎向間距4.5 m、2.5 m，傾角15°，冠梁、腰梁均采用混凝土澆筑.基坑開挖后，樁間采用50 mm厚掛網(wǎng)噴射混凝土支護.基坑回填后，上部4.5 m高邊坡部分支護樁設(shè)200 mm厚鋼筋混凝土擋土板，支護結(jié)構(gòu)剖面圖如圖1所示.將深基坑開挖與支護分為以下六個工況進行：

工況1：在冠梁上設(shè)置錨索并施加預(yù)應(yīng)力值；

工況2：基坑開挖距離地面5 m處；

工況3：距離地面4.5 m處設(shè)置錨索并施加預(yù)應(yīng)力值；

工況4：基坑開挖距離地面10 m處；

工況5：距離地面7 m處設(shè)置錨索并施加預(yù)應(yīng)力值；

工況6：基坑開挖距離地面16 m處.

根據(jù)工程地質(zhì)勘察資料報告所提供地層資料，場地北側(cè)主要為平江穹褶斷裂和潭寧凹褶斷裂兩構(gòu)造單元接觸處,基底巖層為冷家溪群板巖.巖土體力學(xué)參數(shù)如表1所示,錨索計算參數(shù)如表2所示.

表1 巖土體力學(xué)參數(shù)

表2 錨索設(shè)計參數(shù)

2 模型建立與參數(shù)選取

2.1 計算模型

根據(jù)上述工程案例，采用FLAC3D軟件建立如圖2所示數(shù)值單元模型.其中，模型整體長度為80 m，高度為45 m，基坑開挖深度為16 m，寬度取錨桿橫向布置間距為2 m，錨索加固單元與樁單元取FLAC3D中的實體單元cable單元與pile單元，基坑模型底部為固定約束，模型兩側(cè)為法向約束，其上部在自然狀態(tài)下為自由邊界，模型共劃分為2 000個單元，3 198個節(jié)點.

2.2 材料本構(gòu)模型及參數(shù)選取

(1) 材料本構(gòu)模型

數(shù)值模擬中巖土體本構(gòu)模型選取應(yīng)變軟化模型，屈服準(zhǔn)則采用Mohr-Coulomb強度屈服準(zhǔn)則.基于應(yīng)變軟化特性的數(shù)值模型實際為Mohr-Coulomb模型的一種特殊的演化形式[9～11]，該模型與Mohr-Coulomb模型的主要區(qū)別在于塑性應(yīng)變產(chǎn)生后，剪脹角、黏聚力、拉伸強度、摩擦角都可能發(fā)生變化，在此后的每個時間增量步內(nèi)，增量硬化參數(shù)會對總的塑性拉應(yīng)變和剪應(yīng)變進行校驗，使得應(yīng)變軟化模型的參數(shù)自動調(diào)節(jié)為與自定義方程相適應(yīng).

在主應(yīng)力空間σ1-σ3平面內(nèi)，根據(jù)摩爾-庫倫屈服準(zhǔn)則[12，13]，剪切失效包絡(luò)線fs=0，表示為

(1)

式中：c為粘聚力；φ為內(nèi)摩擦角；Nφ=(1+sinφ)/(1-sinφ).

拉伸失效包絡(luò)線ft=0可表示為

ft=σ1-σt.

(2)

數(shù)值計算中須對每個單元體設(shè)置2個硬化參數(shù)ks和kt，由此來作為剪切應(yīng)變與拉應(yīng)變的增量度量的和，單元體中的剪切增量與拉伸硬化增量，可表示為[14，15]

(3)

式中：Δεmps為體塑性剪切應(yīng)變增量,Δεmps=(Δε3ps+Δεmps)2/3；Δε1ps與Δε3ps分別為第1和第3主應(yīng)力方向的塑性剪切應(yīng)變增量.

塑性拉伸應(yīng)變增量則將拉伸硬化增量表示為

(4)

(2) 計算參數(shù)選取

計算參數(shù)的選取對于深基坑樁錨支護在開挖過程中的變形分析至關(guān)重要，并且直接影響到計算結(jié)果的可靠性.對于巖土體的參數(shù)，如表1所示.錨索單元參數(shù)：鉆孔直徑150mm，彈性模量195GPa，內(nèi)摩擦角為25°，泊松比為0.25，黏結(jié)力為18.0kPa，黏結(jié)剛度為1.0GN/m2.支護樁單元參數(shù)：彈性模量為30.0GPa，泊松比為0.35，切向摩擦角28°，切向黏結(jié)力為0.185MN/m，法向黏結(jié)力0.185GN/m，法向剛度為1.5GN/m2，黏結(jié)剛度為1.0GN/m2.基于巖土體的應(yīng)變軟化原理，材料發(fā)生塑性變形后，其定義的強度參數(shù)cp和φp與原始的ci和φi之間的對應(yīng)關(guān)系為：cp=wcci；φp=φi-wφ.其中，wc與wφ分別為描述粘聚力與內(nèi)摩擦角的強度變化因子.wc和wφ與塑性應(yīng)變εp的關(guān)系[11]如表3所示.

表3 wc和wφ與塑性應(yīng)變εp的關(guān)系

3 計算結(jié)果分析

3.1 支護樁頂位移分析

選取支護樁頂?shù)奶卣鼽cA作為基坑開挖過程當(dāng)中水平位移與鉛直位移變化的研究對象，采用應(yīng)變軟化模型與Mohr-Coulomb理想彈塑性模型對該問題進行分析,并結(jié)合實際監(jiān)測數(shù)據(jù)進行對比，從而對基坑開挖過程中的變化規(guī)律進行分析，結(jié)果如圖3、圖4所示.

由圖3可知，隨著工況的進行，特征點A水平位移整體上逐漸增大；開挖完成時，應(yīng)變軟化模型計算所得最大水平位移為16.1 mm，彈塑性模型計算所得最大水平位移為10.5 mm，實際監(jiān)測值為17.5 mm.由此可知：應(yīng)變軟化模型計算結(jié)果明顯大于彈塑性模型計算結(jié)果.由圖4可知，基坑開挖過程中，特征點A的鉛直位移逐漸增大，且應(yīng)變軟化模型計算結(jié)果明顯大于彈塑性模型.在基坑開挖進行至第三工況時，特征點對應(yīng)的水平位移與鉛直位移減小，這是由于基坑在開挖過程中引起樁錨結(jié)構(gòu)的動態(tài)響應(yīng)，錨索發(fā)揮了拉拔效應(yīng)，導(dǎo)致特征點向基坑土體內(nèi)部移動.由圖3和圖4可以看出，基坑開挖過程中，數(shù)值模擬計算得到的位移值與現(xiàn)場實際監(jiān)測的位移值變化趨勢一致，且應(yīng)變軟化模型所得結(jié)果比Mohr-Coulomb模型所得結(jié)果更加接近實際監(jiān)測結(jié)果，說明應(yīng)變軟化模型相對于Mohr-Coulomb模型而言，能夠較好地反映軟巖土體的變形特性.

3.2 支護樁深層水平位移分析

選取支護樁身作為研究對象，并采用數(shù)值模擬所得數(shù)據(jù)與實測數(shù)據(jù)進行對比，從而對支護樁身隨基坑開挖過程中的水平位移變化情況進行分析，結(jié)果如圖5所示.

圖5為基坑在開挖過程中通過測斜儀測出不同工況下樁身水平位移實際監(jiān)測變化值并與開挖完成時的數(shù)值模擬結(jié)果進行對比.由圖5可知，基坑開挖過程中，支護樁身水平位移呈現(xiàn)中間大、兩頭小的變形規(guī)律，主要由于隨著開挖步的不斷增加以及錨索預(yù)應(yīng)力值的施加，錨索與支護樁發(fā)生結(jié)構(gòu)動態(tài)響應(yīng)，錨索扼制住支護樁身的變形，使得樁體的最大水平位移點已不處于支護樁頂處，且呈現(xiàn)逐漸下降的趨勢.此外，由圖5可見，基坑開挖完成時，應(yīng)變軟化模型下的支護樁頂水平位移為16.1 mm，Mohr-Coulomb模型下的支護樁頂水平位移為10.5 mm，實際監(jiān)測數(shù)據(jù)為18.3 mm，應(yīng)變軟化模型計算所得結(jié)果明顯大于彈塑性模型且更接近實測數(shù)據(jù)值.

3.3 錨索的受力分析

為了研究錨索在基坑開挖過程中的受力變化情況，選取各層錨索作為研究對象，采用應(yīng)變軟化模型對在不同工況下，每層錨索軸力分布情況以及基坑開挖完成時各層錨索軸力分布情況進行分析計算，如圖6～9所示.

由圖6～8可知，隨著基坑開挖工況不斷增加，錨索逐漸發(fā)揮效應(yīng)，其錨固力值逐漸增加；錨索錨固力沿索體的分布形態(tài)呈現(xiàn)出先增大后減小的變化趨勢，錨索錨固力向錨索錨固端的深部傳遞，且錨固端遠端軸力值越來越??；隨著工況的進行，錨索的最大軸力點逐漸向后移動，說明隨著基坑的開挖，其影響范圍逐漸擴大，基坑滑動面逐步向坑后移動.由圖9可知，基坑開挖完成時，各層錨索的軸力變化形態(tài)基本保持一致；各層錨索軸力的最大值位置處有所變化，第二排錨索軸力最大值位于錨索長度方向9 m處，為285.5 kN，而第三排錨索軸力最大值位于錨索長度方向6 m處，為380.4 kN，這是由于其軸力最大值處潛在于滑動面.開挖完成時，錨索尾部的錨固力值相差不大，由此也說明錨索尾部變形較小.

3.4 基坑開挖過程中整體穩(wěn)定性分析

為研究基坑在不同工況下的整體穩(wěn)定性的變化規(guī)律，采用強度折減法[16]對基坑在不同工況下的穩(wěn)定性進行計算并分析，結(jié)果如表4所示.從表4可見：在整個的基坑開挖過程中，基坑安全系數(shù)整體呈現(xiàn)下降趨勢.基坑開挖進行至工況3與工況5時，錨索設(shè)置并施加預(yù)應(yīng)力值，錨索發(fā)揮拉拔效應(yīng)，從而使得基坑的整體穩(wěn)定性提高，此時基坑安全系數(shù)增大.隨著工況的進行，基坑開挖完成時較開挖前降低顯著.工況5與工況6進行對比發(fā)現(xiàn)：相對于整個開挖過程而言，基坑的整體安全系數(shù)降低尤為顯著，主要由于在工況6下并未進行相應(yīng)的支護措施，而只是進行基坑的垂直開挖，說明此步開挖對基坑的穩(wěn)定性影響最大，與支護樁點的位移分析相吻合.由表4可知，開挖完成時，應(yīng)變軟化模型計算得到基坑整體安全系數(shù)為1.21，Mohr-Coulomb模型計算得到基坑整體安全性系數(shù)為1.40，基坑開挖應(yīng)變軟化模型計算所得基坑安全系數(shù)要比Mohr-Coulomb模型要小,應(yīng)變軟化模型計算下的基坑相比Mohr-Coulomb模型計算下的基坑是偏于危險的.

表4 各個工況下基坑安全系數(shù)的變化情況

3.5 塑性區(qū)發(fā)展規(guī)律分析

圖10和圖11為深基坑在樁錨支護過程中由于巖土體的應(yīng)變軟化特性所引起的張拉、剪切塑性區(qū)分布及發(fā)展變化云圖.由圖10和圖11可以看出，在基坑開挖過程中，塑性區(qū)主要分布于基坑開挖面及開挖面3～6 m深度內(nèi)，且基坑并未出現(xiàn)整體滑裂帶，也并未產(chǎn)生深層的塑性剪切帶；隨著基坑開挖的進行，基坑土體塑性區(qū)域明顯增大，但并未存在潛在的滑裂面，塑性區(qū)亦未向深部發(fā)展.對比圖5與圖6可知，隨著基坑的開挖進行，巖土體的應(yīng)力迅速釋放，剪切塑性破壞區(qū)繼續(xù)發(fā)展的同時，在基坑局部區(qū)域出現(xiàn)了明顯的張拉剪切破壞區(qū).

4 結(jié) 論

(1) FLAC3D應(yīng)變軟化模型數(shù)值模擬結(jié)果與實際監(jiān)測結(jié)果基本吻合，且優(yōu)于Mohr-Coulomb模型，并能較好地反映軟巖土體的變形特性，更加合理地反映軟巖基坑開挖的變形情況，并做出合理地支護措施，對相類似的工程實踐具有指導(dǎo)意義.

(2) 在基坑開挖過程中，樁身最大水平位移點隨著開挖步的不斷增加呈現(xiàn)逐漸下降趨勢；錨索錨固力逐漸增大，且錨索錨固力沿索體的分布形態(tài)呈現(xiàn)出先增大后減小的變化規(guī)律；基坑開挖完成時，各層錨索的軸力變化形態(tài)基本一致，且最大錨固力值出現(xiàn)在剪切滑動面位置處，錨索尾部的錨固力相差不大.

(3) 由于受到基坑巖土體應(yīng)變軟化特性以及開挖擾動等各種因素的影響，基坑開挖過程中的塑性區(qū)域主要集中分布在基坑開挖面及開挖面3～7 m深度內(nèi).

[1] 張帆，盛謙，朱澤奇，等．三峽花崗巖峰后力學(xué)特性及應(yīng)變軟化模型研究[J].巖石力學(xué)與工程學(xué)報，2008，27(增 1): 2 651-2 655．

[2] 楊超，崔新明，徐水平.軟巖應(yīng)變軟化數(shù)值模型的建立與研究[J].巖土力學(xué)，2002，23(6):695-697，701．

[3] 王水林，吳振君，李春光，等． 應(yīng)變軟化模擬與圓形隧道襯砌分析[J].巖土力學(xué)，2010，31(6):1 929-1 936．

[4] 盧敦華,王星華.基于應(yīng)變軟化模型的深基坑土釘加固分析[J].解放軍理工大學(xué)學(xué)報, 2009, 10(1): 61-66.

[5] 周勇,王濤,呂慶,等.基于FLAC3D巖石應(yīng)變軟化模型的研究[J].長江科學(xué)院院報,2012,29(5)：51-56.

[6] 蔡海波,吳順川，周喻.既有基坑延深開挖穩(wěn)定性評價與支護方案確定[J].巖土力學(xué), 2011, 32(11): 3 306-3 311.

[7] 劉自由.基坑應(yīng)變軟化的樁錨支護數(shù)值計算[J].中南大學(xué)學(xué)報(自然科學(xué)版), 2012, 43(7): 2 832-2 837.

[8] 鄒志強,萬宇. 基坑樁錨聯(lián)合支護分析及彈塑性驗算[J].建筑科學(xué), 2013 ,29(7):24-29.

[9] 周斌，張可能，柳群義.考慮應(yīng)變軟化特征的層狀巖土體三軸壓縮數(shù)值試驗分析[J].中南大學(xué)學(xué)報(自然科學(xué)版), 2009, 40(3): 828-832.

[10] 徐志英.巖石力學(xué)[M].北京：水利水電出版社，1993.

[11] Itasca Consulting Group . User's guid [M].Minnesota: Itasce Consulting Group, 2002.

[12] 劉繼國,曾亞武. FLAC3D在深基坑開挖與支護數(shù)值模擬中的應(yīng)用[J].巖土力學(xué), 2006 ,27(3):505-508.

[13] 張嘎, 張建民. 基于瑞典條分法的應(yīng)變軟化邊坡穩(wěn)定性評價方法[J].巖石力學(xué), 2007, 28(1): 12-16.

[14] 李廣信. 高等土力學(xué)[M]. 北京: 清華大學(xué)出版社, 2004: 41-52.

[15] 何滿潮. 軟巖工程力學(xué)[M]. 北京: 科學(xué)出版社, 2003: 1-3.

[16] 趙尚毅,鄭穎人,時衛(wèi)民,等.用有限元強度折減法求邊坡穩(wěn)定安全系數(shù)[J] .巖土工程學(xué)報,2002,24(3):343-346.

責(zé)任編輯：羅 聯(lián)

Process Stability Analysis of Pile-Anchor Supporting Based on Strain Softening Model

CHENBin1*，ZHANGLiang2，LIUZheng-cai2

(1.School of Energy Engineering, Xiangtan University, Xiangtan 411105; 2. School of Civil Engineering and Mechanics, Xiangtan University, Xiangtan 411105 China)

In order to study the change law of stability during deep foundation pit pile-anchor supporting process，deformation of pile, anchor cable force, changes of safety coefficient and plastic zone distribution were calculated and analyzed by taking FLAC3Dsoftware based on strain softening model during the deep foundation pit excavation process, and compared with the Mohr-Coulomb model, and analyzes the results combining with the measured data.The results show that the strain softening model results are in substantially the same with the actual monitoring results，which is better than Mohr-Coulomb model, and the strain softening model can reflect the strain softening properties of soft rock better.During foundation pit excavation process, the maximum horizontal displacement of supporting pile showed a gradual downward trend with the increasing of excavation step. Along with the excavation, the overall anchor force increases gradually,which causes the anchoring force to first increase and then decrease, and the maximum anchor force appeared in the shear slip surface. The distribution of slope shear and tension plastic zone mainly concentrate on the foundation pit surface layer within 3～7 m.

strain softening; pile-anchor supporting; deformation characteristic; stability; plastic zone

2014-12-26

湖南省自然科學(xué)基金項目(14JJ7038)；湖南省教育廳項目(12C0377)

陳賓(1977— ),男，河南 駐馬店人，博士，副教授.E-mail: chenbinhnxt@126.com

TU457

A

1000-5900(2015)01-0031-06