抗滑樁加固邊坡三維數(shù)值分析中的幾個問題

年廷凱 ，徐海洋，劉紅帥

（1. 大連理工大學(xué) 土木水利學(xué)院 海岸和近海工程國家重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116024；

2. 成都理工大學(xué) 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國家重點(diǎn)實(shí)驗(yàn)室，成都 610059；3. 中國地震局工程力學(xué)研究所，哈爾濱 150080）

1 引 言

抗滑樁是一種大截面、側(cè)向受荷的排樁或樁群，它穿過滑體錨入滑床一定深度，借助與樁周巖土的共同作用，將滑坡推力傳遞到穩(wěn)定地層，其抗滑機(jī)制體現(xiàn)于樁、滑體與滑床3 者間的相互協(xié)調(diào)工作[1]。幾十年來，許多學(xué)者針對抗滑樁開展了大量的研究，綜合現(xiàn)有的文獻(xiàn)成果，大體上可分為3 類：（1）理論分析，通過建立相應(yīng)的數(shù)學(xué)或物理模型，研究土體側(cè)向移動對抗滑樁的影響[2－9]；（2）數(shù)值模擬，主要利用有限元或有限差分法，建立二維或三維數(shù)值計(jì)算模型，研究樁-土相互作用下樁身變形與內(nèi)力分布，樁周土體應(yīng)力和變形情況[10－16]；（3）試驗(yàn)研究，通過室內(nèi)小比尺物理模型與離心試驗(yàn)等[17－18]，直觀研究樁-土相互作用及抗滑樁工作性能。

在這3 類研究工作中，數(shù)值分析方法得到了快速發(fā)展和廣泛關(guān)注，特別是基于強(qiáng)度折減技術(shù)的有限元或有限差分法[10－17]，給抗滑樁加固邊坡的數(shù)值模擬注入了新的活力。然而，通過綜合分析文獻(xiàn)資料，筆者注意到抗滑樁數(shù)值分析中仍有幾個問題有待明晰，如抗滑樁計(jì)算模型尺度、樁底接觸模式、加固位置、適宜樁長、樁間距與樁徑比、樁頭約束條件及臨界滑動面等。為此，利用考慮樁-土-邊坡相互作用的強(qiáng)度折減有限元方法[19]，結(jié)合典型算例邊坡，開展了抗滑樁加固邊坡的三維數(shù)值分析，著重研究了抗滑樁-邊坡體系的計(jì)算模型尺度、設(shè)樁位置、樁間距與樁徑比、樁長與樁底接觸模式等因素對邊坡安全系數(shù)及臨界滑動面的影響，以及不同樁頭約束下抗滑樁內(nèi)力分布等，以期為抗滑樁工程設(shè)計(jì)及規(guī)范修訂提供參考。

2 強(qiáng)度折減有限元法

利用有限元程序進(jìn)行強(qiáng)度折減計(jì)算，邊坡土體采用服從Mohr-Coulomb 破壞準(zhǔn)則與非關(guān)聯(lián)流動法則的理想彈塑性本構(gòu)模型?？够瑯侗灰暈槔硐刖€彈性體，樁-土間采用接觸相互作用，摩擦系數(shù)取為0.3，約為tan(0.75φ)。以數(shù)值迭代不收斂并結(jié)合坡面特征點(diǎn)位移陡增作為邊坡失穩(wěn)判據(jù)[19]，即在此時(shí)的強(qiáng)度折減系數(shù)為邊坡穩(wěn)定安全系數(shù)。強(qiáng)度折減計(jì)算中，折減后的強(qiáng)度參數(shù)表達(dá)式為

式中：c′和φ′為土體的實(shí)際抗剪強(qiáng)度參數(shù)；mc′ 和mφ′為土體發(fā)揮的抗剪強(qiáng)度參數(shù)；SRF 為強(qiáng)度折減系數(shù)。

圖1 抗滑樁-邊坡模型的側(cè)面、平面圖 Fig.1 Side view and plan view for slope-pile model

表1 抗滑樁-邊坡土性質(zhì)參數(shù) Table 1 Properties of soil and pile

3 抗滑樁-邊坡計(jì)算模型尺度問題

針對當(dāng)前抗滑樁加固三維邊坡數(shù)值計(jì)算中樁-邊坡體系計(jì)算模型尺度選擇的問題，這里建立單樁取半、單樁、雙樁取半、雙樁、單樁加雙樁取半5種尺度有限元計(jì)算模型，其整體網(wǎng)格劃分及樁與周圍土體的局部網(wǎng)格型式如圖2 所示。其中數(shù)字①～⑤表示抗滑樁-邊坡計(jì)算模型尺度（寬度）不同，依次代表單樁取半（計(jì)算模型寬度為0.5S，樁取半樁）、單樁（計(jì)算模型寬度為1S，取全樁）、雙樁取半（計(jì)算模型寬度為1S，取兩個半樁）、雙樁（計(jì)算模型寬度為2S，取兩個全樁）、單樁加雙樁取半（計(jì)算模型寬度為2S，取兩個半樁加一個全樁）。為了減小接觸面單元的病態(tài)條件，抗滑樁周圍土體采用較細(xì)的網(wǎng)格形式，土與樁的單元形式皆采用六面體減縮積分。所有模型采用統(tǒng)一的邊界條件，即邊坡兩側(cè)面Z 方向水平位移約束，前后兩面X 方向水平位移約束，底面X、Y、Z 方向全約束。

利用強(qiáng)度折減有限元法計(jì)算加固前邊坡（模型寬度Z =3.0～16.0 m）的安全系數(shù)為1.158～1.162，與Cai 等[10]、Won 等[12]和Wei 等[13]所得的安全系數(shù)1.13～1.20 基本一致，微小的差別可能源于有限元網(wǎng)格、數(shù)值算法、失穩(wěn)判據(jù)等因素?？紤]采用抗滑樁加固方案，有關(guān)參數(shù)同前，分別計(jì)算了5 種尺度有限元模型下的邊坡安全系數(shù)，所得結(jié)果列于表2。

圖2 5 種尺度有限元計(jì)算模型及網(wǎng)格形式 Fig.2 Five computational models and the FE meshes

表2 5 種尺度有限元計(jì)算模型下的邊坡安全系數(shù) Table 2 FOS of slope with a row of piles under five models

分析表2 可知，當(dāng)樁間距與樁徑之比S/D 和樁長(Lp)一定時(shí)，5 種尺度計(jì)算模型所得邊坡的安全系數(shù)基本相同，相對誤差不超過3%，說明在計(jì)算樁-邊坡體系安全系數(shù)時(shí)，有限元計(jì)算模型尺度取0.5S、1S（兩個半樁或全樁）和2S（兩個全樁或兩個半樁加一個全樁），所得安全系數(shù)并無差異；所不同的是，在網(wǎng)格大小劃分一致的情況下，樁-邊坡體系計(jì)算模型尺度取0.5S 時(shí)計(jì)算工作量將大大減少，這在很大程度上提高了三維數(shù)值計(jì)算的效率。

為進(jìn)一步驗(yàn)證上述結(jié)論的正確性，將抗滑樁的彈性模量調(diào)整為60 000 MPa，土的變形模量調(diào)整為200 MPa，計(jì)算了樁間距S 為3.0 m，樁徑D 為0.8 m和1.0 m 時(shí)的樁-邊坡體系安全系數(shù)，結(jié)果列于表3。

表3 修改材料參數(shù)后樁-邊坡體系的安全系數(shù) Table 3 FOS of slope with piles after modified parameters

分析表3 同樣可以得出前述結(jié)論，即5 種尺度計(jì)算模型在相同參數(shù)條件下所得邊坡安全系數(shù)相同。特別地，當(dāng)樁徑為D =0.8 m 時(shí)，加固邊坡的安全系數(shù)為1.730，與Wei 等[13]利用有限差分法所得結(jié)果1.72 相一致，從而驗(yàn)證了上述方法的可行性。

此外，分析表3 還可發(fā)現(xiàn)，不同的樁底接觸模式（接觸與固定模式）對抗滑樁加固邊坡的安全系數(shù)有一定的影響，固定模式下的安全系數(shù)明顯大于接觸模式；對比表2、3 可見，提高樁與土的彈性模量，邊坡的安全系數(shù)也在一定程度上得到改善，樁的抗彎剛度增加，一定程度上提高了抗滑效果。

4 最佳設(shè)樁位置問題

利用前述邊坡算例，取樁間距為S =3.0、4.0、5.0 m，改變抗滑樁在坡體中的位置，使Lx從0 變化到15.0 m，所得安全系數(shù)隨Lx與坡面水平長度L比值的變化曲線如圖3 所示。分析這3 條曲線，當(dāng)設(shè)樁于Lx=7.5 m 處，即坡中位置，樁-邊坡體系獲得最大的安全系數(shù)；當(dāng)設(shè)樁于坡頂或坡腳處，安全系數(shù)略高于無樁狀態(tài)，但提高幅度不顯著。進(jìn)一步采用二次拋物線對3 條曲線進(jìn)行擬合，得到相應(yīng)的拋物線方程，如圖3 中所示，其對稱軸位于Lx/L=0.5處，說明在該處邊坡加固效果最好；當(dāng)樁間距比S/D較小時(shí)（圖中S/D =3）拋物線開口較小，不同設(shè)樁位置處的安全系數(shù)明顯高于樁間距比S/D =4 和5 時(shí)的情況（拋物線趨于平緩），說明樁間距比影響著加固邊坡的安全系數(shù)。

圖4 顯示抗滑樁不同加固位置時(shí)的邊坡臨界滑動面。對比分析可見，當(dāng)抗滑樁位于邊坡中部時(shí)，其滑裂面以樁為界被分為前、后兩個部分，起到了最有效的阻滑作用，這與無樁狀態(tài)下臨界滑裂面反映出的邊坡中部出現(xiàn)高應(yīng)變區(qū)是一致的；而設(shè)樁于坡頂與坡腳附近時(shí)，滑裂面基本為一條圓弧形的曲面，抗滑樁的阻滑作用并不明顯。

圖3 安全系數(shù)隨樁加固位置的變化曲線 Fig.3 Effects of pile location on FOS

圖4 不同設(shè)樁位置時(shí)邊坡的臨界滑裂面 Fig.4 Critical slip surface at the different pile location

5 合理樁間距與樁徑比問題

仍采用前述算例參數(shù)，改變樁間距S 值，得到不同樁間距與樁徑比S/D（以下簡稱間距比）條件下邊坡的安全系數(shù)變化曲線，如圖5 所示。隨著間距比S/D 增大，安全系數(shù)逐漸減小并趨于無樁狀態(tài)；當(dāng)S/D≤6 時(shí)，邊坡安全系數(shù)隨間距比S/D 的增大而大幅減小，這說明抗滑樁的群樁效應(yīng)在快速減弱；當(dāng)S/D>6 時(shí)，隨著S/D 增加，邊坡安全系數(shù)小幅下降并趨于無樁狀態(tài)，這反應(yīng)群樁效應(yīng)已不明顯，只有單樁在起阻滑作用[16]。

圖5 安全系數(shù)隨樁間距比變化曲線 Fig.5 Effect of ratio of pile spacing to pile diameter on FOS

圖6 顯示不同樁間距比S/D 條件下的邊坡臨界滑裂面。其中圖6(a)、(c)、(e)和(g)顯示過樁中心線的截面，圖6(b)、(d)、(f)和(h)顯示兩樁之間的中心截面，其中z=0 為過樁形心的剖面，z=1、2、3、4 m分別為不同S/D 時(shí)過兩樁中心的剖面。分析圖6(b)可見，當(dāng)S/D =2 時(shí)，滑裂面以樁為界被分為前、后2 個部分，此時(shí)抗滑樁發(fā)揮擋墻的功效；當(dāng)S/D =4時(shí)，群樁效應(yīng)已發(fā)揮，樁間有土拱產(chǎn)生，如圖6(d)；而圖6(f)、(h)表明，當(dāng)S/D≥6 邊坡破壞時(shí)，其臨界滑裂面為一條連續(xù)貫通的圓弧形曲面，且滑裂面的深度基本相同，說明樁與樁間的相互作用減弱，土拱效應(yīng)已接近消失，這與圖5 結(jié)果是一致的。

變化樁徑D 從0.6 m 到2.0 m，考察安全系數(shù)隨樁徑的變化關(guān)系，所得曲線如圖7 所示。當(dāng)樁間距為S =3.0 m 時(shí)，隨著樁徑的增大，安全系數(shù)呈線性增長。間距比S/D 不變時(shí)，樁徑在0.6～2.0 m 之間變化（不考慮彈性模量的變化）所得安全系數(shù)基本相同，說明抗滑樁加固邊坡工程，根據(jù)預(yù)期的安全系數(shù)，樁間距與樁徑有多種可供選擇方案。

6 合理樁長問題

采用算例參數(shù)，考察樁長Lp從4.3 m 到15.3 m變化時(shí)邊坡的安全系數(shù)，所得結(jié)果繪于圖8。

圖6 不同間距比條件下的邊坡臨界滑裂面 Fig.6 Effect of ratio of pile spacing to pile diameter on CSS

圖7 安全系數(shù)隨樁徑變化曲線 Fig.7 Effect of pile diameter on FOS

圖8 安全系數(shù)隨樁長的變化曲線 Fig.8 Effect of pile length on FOS

分析圖8 可見，隨著樁長增加，邊坡安全系數(shù)不斷增大；當(dāng)樁長4.3≤Lp≤7.3 m 時(shí)，安全系數(shù)呈線性變化關(guān)系；而當(dāng)Lp＞7.3 m 時(shí)，曲線斜率突然變大，邊坡安全系數(shù)顯著提高；但當(dāng)樁長小于4.3 m時(shí)無法改善邊坡的穩(wěn)定性，因?yàn)榇藭r(shí)樁長小于臨界滑動面深度，這從圖9 中臨界滑動面深度hc=4.3 m 可反映出。進(jìn)一步分析圖9，當(dāng)樁長Lp≤7.3 m 時(shí)，塑性應(yīng)變剪切帶最深處基本都通過樁底，說明加固效果不理想；當(dāng)Lp>7.3 m 后，塑性應(yīng)變剪切帶沿樁周發(fā)展，而沒有擴(kuò)展至樁底，由此表明，抗滑樁已起到明顯的抗滑效果；這與圖8 顯示的當(dāng)樁長Lp≥7.3 m 時(shí)安全系數(shù)顯著增加（曲線斜率明顯增大）的變化趨勢是一致的，由此可將此時(shí)的樁長 7.3 m作為設(shè)計(jì)樁長，相當(dāng)于錨固深度為2/5樁長（臨界滑動面深度為hc=4.3 m），這與抗滑樁設(shè)計(jì)手冊[20]中建議的土質(zhì)邊坡抗滑樁錨固深度約為1/3～1/2 樁長是一致的。

7 樁頭約束效應(yīng)問題

考慮4 種樁頭約束條件[5,10]，即自由樁頭（可任意轉(zhuǎn)動和位移）、無轉(zhuǎn)動樁頭（可無轉(zhuǎn)動地發(fā)生位移）、鉸接樁頭（無位移地轉(zhuǎn)動）和固定樁頭（無轉(zhuǎn)動也無位移），利用算例邊坡中的樁土參數(shù)，分別計(jì)算樁底接觸模式（Lp=13.0 m）與固定模式（Lp= 15.3 m）情況下的樁身彎矩、剪力及樁側(cè)位移、土壓力分布，所得固定模式結(jié)果如圖10 所示。

從圖可見，鉸結(jié)和固定樁頭條件下樁身的彎矩較小，在8.0 m 以下（嵌巖段）樁身彎矩分布更小，表明約束樁頭位移，能有效減小樁身彎矩，并使之分布趨于合理，充分發(fā)揮抗滑潛力；而不轉(zhuǎn)動和自由樁頭條件下樁身彎矩相對較大。對比2 種不同樁底接觸模式下的彎矩或剪力圖可見，樁底接觸模式的選取對樁身最大彎矩和剪力的影響較小，而對其沿樁身的分布模式有一定影響；樁底固定模式下的樁身位移和樁側(cè)土壓力，總體上小于樁底接觸模式。

圖9 不同樁長時(shí)邊坡的臨界滑裂面 Fig.9 Critical slip surface at different pile lengths

圖10 樁底固定模式下的樁身內(nèi)力及變位 Fig.10 Behavior of pile with fixed mode under four boundary conditions

8 結(jié) 論

（1）三維有限元數(shù)值計(jì)算中，單樁取半、單樁、雙樁取半、雙樁、單樁加雙樁取半5 種尺度計(jì)算模型所得邊坡的安全系數(shù)并無差異；所不同的是，在網(wǎng)格大小劃分一致的情況下，樁-邊坡體系計(jì)算模型尺度取0.5S 時(shí)的計(jì)算工作量大大減少，效率倍增；

（2）抗滑樁加固于邊坡中部可獲得最大的安全系數(shù)，而位于坡頂或坡腳處所得安全系數(shù)略高于無樁狀態(tài)，總體上其安全系數(shù)與設(shè)樁位置的變化曲線近似為一拋物曲線；

（3）抗滑樁加固邊坡的安全系數(shù)隨樁間距與樁徑之比S/D 的增加而減小，最優(yōu)間距比宜為S/D =2～6，此時(shí)樁間存在土拱效應(yīng)；間距比S/D 不變情況下，增加樁徑安全系數(shù)基本不變；而樁間距S不變，增加樁徑D，則安全系數(shù)呈線性增加；均質(zhì)土坡中抗滑樁設(shè)計(jì)樁長宜為2.5倍臨界滑動面深度。

（4）本文所得結(jié)論適用于直線型均質(zhì)土坡，對于其他幾何形狀和荷載條件的復(fù)雜邊坡，仍需開展進(jìn)一步的研究工作。

[1] 陳立新, 王士川. 抗滑樁的彈塑性理論分析[J]. 工業(yè)建筑, 1997, 27(7): 28－33. CHEN Li-xin, WANG Shi-chuan. The analysis of elastic-plastic design theory of anti-slide piles[J]. Industrial Construction, 1997, 27(7): 28－33.

[2] ITO T, MATSUI T. Methods to estimate lateral force acting on stabilizing piles[J]. Soils and Foundations, 1975, 15(4): 43－59.

[3] ITO T, MATSUI T, HONG W P. Design method for stabilizing piles against landslide—one row of piles[J]. Soils and Foundations, 1981, 21(1): 21－37.

[4] ITO T, MATSUI T, HONG W P. Extended design method for multi-row stabilizing piles against landslide[J]. Soils and Foundations, 1982, 22(1): 1－13.

[5] HASSIOTIS S, CHAMEAU J L, GUNARATNE M. Design method for stabilization of slopes with piles[J]. Journal of Geotechnical and Geoenvironmental Engineering, 1997, 123(4): 314－323.

[6] 周德培, 肖世國, 夏雄. 邊坡工程中抗滑樁合理樁間距的探討[J]. 巖土工程學(xué)報(bào), 2004, 26(1): 132－135. ZHOU De-pei, XIAO Shi-guo, XIA Xiong. Discussion on rational spacing between adjacent anti-slide piles in some cutting slope projects[J]. Chinese Journal of Geotechnical Engineering, 2004, 26(1): 132－135.

[7] 張建勛, 陳福全, 簡洪鈺. 被動樁中土拱效應(yīng)問題的數(shù)值分析[J]. 巖土力學(xué), 2004, 25(2): 174－178. ZHANG Jian-xun, CHEN Fu-quan, JIAN Hong-yu. Numerical analysis of soil arching effects in passive piles[J]. Rock and Soil Mechanics, 2004, 25(2): 174－178.

[8] 呂慶, 孫紅月, 尚岳全. 抗滑樁樁后土拱形狀及影響因素[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2010, 42(4): 629－633. Lü Qing, SUN Hong-yue, SHANG Yue-quan. Shape of soil arch behind anti-slide piles and its major influence factors[J]. Journal of Harbin Institute of Technology, 2010, 42(4): 629－633.

[9] 趙明華, 廖彬彬, 劉思思. 基于拱效應(yīng)的邊坡抗滑樁樁間距計(jì)算[J]. 巖土力學(xué), 2010, 31(4): 1211－1216. ZHAO Ming-hua, LIAO Bin-bin, LIU Si-si. Calculation of anti-slide piles spacing based on soil arching effect[J]. Rock and Soil Mechanics, 2010, 31(4): 1211－1216.

[10] CAI F, UGAI K. Numerical analysis of the stability of a slope reinforced with piles[J]. Soils and Foundations, 2000, 40(1): 73－84.

[11] LIANG R, ZENG S. Numerical study of soil arching mechanism in drilled shafts for slope stabilization[J]. Soils and Foundations, 2002, 42(2): 83－92.

[12] WON JINOH, YOU KWANGHO, JEONG SANGSEOM, et al. Coupled effects in stability analysis of pile-slope systems[J]. Computers and Geotechnics, 2005, 32(4): 304－315.

[13] WEI W B, CHENG Y M. Strength reduction analysis for slope reinforced with one row of piles[J]. Computers and Geotechnics, 2009, 36(7): 1176－1185.

[14] 高長勝, 陳生水, 楊守華, 等. 基于強(qiáng)度折減有限單元法的抗滑樁加固邊坡特性分析[J]. 水利與建筑工程學(xué)報(bào), 2010, 8(4): 119－122. GAO Chang-sheng, CHEN Sheng-shui, YANG Shou-hua, et al. Analysis for slope stability reinforced with anti-slide piles based on strength reduction method with3D FEM[J]. Journal of Water Resources and Architectural Engineering, 2010, 8(4): 119－122.

[15] DANIEL PRADEL, JASON GARNER, ANNIE ON LEI KWOK. Design of drilled shafts to enhance slope stability[C]//Proceedings of 2010 Earth Retention Conference. Los Angeles: [s. n.], 2010: 920－927.

[16] KOURKOULIS R, GELAGOTI F, ANASTASOPOULOS I, et al. Slope stabilizing piles and pile-groups: parametric study and design insights[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2011, 137(7).doi: 10.1061/(ASCE)GT.1943-5606. 0000479.

[17] 楊明, 姚令侃, 王廣軍. 樁間土拱效應(yīng)離心模型試驗(yàn)及數(shù)值模擬研究[J]. 巖土力學(xué), 2008, 29(3): 817－822. YANG Ming, YAO Ling-kan, WANG Guang-jun. Study of centrifuge model tests and numerical simulation on soil arching in space of piles[J]. Rock and Soil Mechanics, 2008, 29(3): 817－822.

[18] 高長勝, 魏汝龍, 陳生水. 抗滑樁加固邊坡變形破壞特性離心模型試驗(yàn)研究[J]. 巖土工程學(xué)報(bào), 2009, 31(1): 145－148. GAO Chang-sheng, WEI Ru-long, CHEN Sheng-shui. Centrifugal model tests on deformation of slopes reinforced with piles[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(1): 145－148.

[19] 年廷凱. 樁-土-邊坡相互作用數(shù)值分析及阻滑樁簡化設(shè)計(jì)方法研究[D]. 大連: 大連理工大學(xué), 2005.

[20] 鐵道部第二勘測設(shè)計(jì)院. 抗滑樁設(shè)計(jì)與計(jì)算[M]. 北京: 中國鐵道出版社, 1983.