抗滑樁加固邊坡的穩(wěn)定性分析及最優(yōu)樁位的確定

彭文哲 趙明華 肖堯 楊超煒

摘 ? 要：為合理分析抗滑樁加固后的邊坡穩(wěn)定性及確定最危險(xiǎn)滑動(dòng)面和最優(yōu)樁位、樁長，基于有限元極限分析軟件OptumG2建立抗滑樁加固前后的邊坡數(shù)值模型，通過與已有研究中的邊坡算例進(jìn)行對比，驗(yàn)證了本文數(shù)值分析方法的合理性;然后，基于算例驗(yàn)證的對比結(jié)果，討論了OptumG2的2種分析類型（重力乘數(shù)與強(qiáng)度折減）的區(qū)別并認(rèn)為基于強(qiáng)度折減極限分析得出的安全系數(shù)偏于安全;最后，探討了邊坡坡度及抗剪強(qiáng)度（內(nèi)摩擦角與粘聚力）對安全系數(shù)、最優(yōu)樁位及樁長的影響，分析了抗滑樁加固后的最危險(xiǎn)滑動(dòng)面變化規(guī)律，根據(jù)參數(shù)分析結(jié)果擬合出安全系數(shù)公式，得出一些規(guī)律性的結(jié)論并總結(jié)了4種常見的滑動(dòng)面形式及其形成條件，可為后續(xù)邊坡穩(wěn)定性分析理論研究提供參考，具有一定的理論及工程應(yīng)用價(jià)值.

關(guān)鍵詞：抗滑樁;邊坡穩(wěn)定性;有限元極限分析;最優(yōu)樁位

中圖分類號(hào)：TU473 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼：A

Abstract：This paper is to analyze the stability with anti-slide pile reinforced slope and determine the critical failure surface， optimal pile position and pile length. Firstly， a numerical model of slope was established based on the finite element limit analysis software OptumG2. Secondly， two cases were employed to verify the rationality of the numerical analysis. Then， the differences between two different analysis types （gravity multiplier and strength reduction） of OptumG2 were discussed based on the comparison results and the safety factor based on the strength reduction limit analysis was biased towards safety. Finally， the effects of slope angle and soil shear strength （internal friction angle and cohesion） on safety factor， optimal pile position and pile length were investigated. In addition， the variation law of the critical failure surface of anti-slide pile reinforced slope was analyzed and the fitting equation of safety factor was derived based on the parametric study results. Some regular conclusions were presented， and four common sliding surface forms and their forming conditions were summarized. The results can provide reference for the subsequent theoretical research of slope stability analysis.

Key words：anti-slide pile;slope stability;finite element limit analysis;optimal pile position

邊坡穩(wěn)定性分析是邊坡處治工程設(shè)計(jì)的依據(jù)，作為最直觀評價(jià)標(biāo)準(zhǔn)之一，邊坡穩(wěn)定性安全系數(shù)的分析方法已日趨完善[1-7]. 然而，自然狀態(tài)下邊坡的安全系數(shù)往往難以達(dá)到一些重要工程的要求，因此，需要對邊坡采取一定的加固措施改善其穩(wěn)定性. 抗滑樁因其抗滑能力強(qiáng)、適用范圍廣及施工便捷等優(yōu)點(diǎn)而被廣泛使用[8-9]，相比抗滑樁加固前的邊坡，加固后的邊坡穩(wěn)定性分析更為復(fù)雜，已成為當(dāng)前重要研究課題之一.

目前國內(nèi)外諸多學(xué)者對經(jīng)抗滑樁加固的邊坡穩(wěn)定性研究主要可分為試驗(yàn)研究[10-14]、理論分析[15-19]及數(shù)值模擬[20-24]三個(gè)方面.

孫書偉等[10]基于室內(nèi)試驗(yàn)探究剛度較大抗滑樁的抗滑機(jī)制，認(rèn)為樁后土體壓裂破壞是抗滑樁變形的主要原因;戴自航等[11]依托福建省某大型滑坡治理工程進(jìn)行大型推樁現(xiàn)場試驗(yàn)，對抗滑樁分別施加矩形、三角形及拋物線形水平分布荷載，揭示不同分布形式滑坡推力作用下的抗滑樁受力及變形特性;劉洪佳等[12]基于抗滑樁模型試驗(yàn)實(shí)測數(shù)據(jù)，分析滑坡推力分布形式、樁前土抗力的變化規(guī)律及抗滑樁-邊坡系統(tǒng)的失穩(wěn)模式，認(rèn)為導(dǎo)致失穩(wěn)的主要原因是樁前坡體達(dá)到屈服狀態(tài)，樁頂位移過大，從而導(dǎo)致樁體發(fā)生折斷破壞;傅翔等[13]通過室內(nèi)推樁模型試驗(yàn)，探究抗滑樁的傾覆破壞模式及相應(yīng)的變形特征;章為民等[14]為研究滑動(dòng)面位置（深度）對抗滑樁的抗滑機(jī)理及其失穩(wěn)破壞模式的影響，引入無量綱參數(shù)（錨固比），提出4種抗滑及破壞模式及其界限判別條件，并基于極限平衡理論推導(dǎo)出相應(yīng)的極限抗滑力理論解，最終通過模型試驗(yàn)對此進(jìn)行驗(yàn)證. 然而由于試驗(yàn)難度較大，且難以量化邊坡穩(wěn)定性，國內(nèi)外不少學(xué)者采用理論方法分析其穩(wěn)定性. Hassiotis等[15]基于摩擦圓法對邊坡進(jìn)行穩(wěn)定性分析，并引入考慮抗滑樁的穩(wěn)定數(shù)來確定新的最危險(xiǎn)滑動(dòng)面和安全系數(shù);Ausilio等[16]采用動(dòng)態(tài)極限分析方法探討了抗滑樁對邊坡穩(wěn)定性的影響，并提出抗滑樁最優(yōu)位置的確定方法;李家平等[17]采用瑞典法及簡化Bishop法探討了抗滑樁對邊坡的加固效果，并認(rèn)為抗滑樁最優(yōu)樁位應(yīng)為最危險(xiǎn)滑動(dòng)面最低處;年廷凱等[18]基于強(qiáng)度折減理論及極限分析上限定理，確定了土坡的安全系數(shù)及潛在破壞模式，并在建立抗滑樁加固土坡的極限平衡狀態(tài)方程的同時(shí)探討了抗滑樁的最優(yōu)樁位問題;梁冠亭等[19]基于Morgenstern-Price法及自適應(yīng)遺傳優(yōu)化算法，提出了抗滑樁加固邊坡后的穩(wěn)定性分析優(yōu)化模型，搜索非圓弧最危險(xiǎn)滑動(dòng)面并探討樁位對邊坡穩(wěn)定性的影響.

然而，以上理論研究均需對邊坡加固前后的最危險(xiǎn)滑動(dòng)面形狀及位置進(jìn)行假定，難以全面反映邊坡-抗滑樁系統(tǒng)的耦合效應(yīng). 考慮到數(shù)值分析方法能有效地彌補(bǔ)這一不足，雷文杰等[20]基于有限元強(qiáng)度折減法對抗滑樁加固的邊坡進(jìn)行數(shù)值模擬，探討樁位對邊坡的穩(wěn)定安全系數(shù)、最危險(xiǎn)滑動(dòng)面位置及形狀的影響;李榮建等[21]自行開發(fā)強(qiáng)度折減有限元程序，用以考慮抗滑樁及基質(zhì)吸力對穩(wěn)定性的影響;胡新麗等[22]以張桓候廟園區(qū)東側(cè)滑坡為例，通過有限元軟件Geostudio，探討不同樁位下的穩(wěn)定性變化規(guī)律;陳樂求等[23]采用Fortan95語言編制有限元程序，用以分析抗滑樁加固后的邊坡穩(wěn)定性，并探討了抗滑樁支護(hù)參數(shù)對安全系數(shù)和滑動(dòng)面的影響;王聰聰?shù)萚24]基于有限差分軟件FLAC3D分析抗滑樁對邊坡的加固效果，探討樁位、樁長及樁體彈性模量等因素對邊坡穩(wěn)定系數(shù)及臨界滑移面的影響，并提出了有效嵌固深度的概念.

綜上可知，目前對抗滑樁加固后邊坡穩(wěn)定性分析已取得了一定的成果，但仍存在以下不足：①理論研究多為極限平衡法，需在事先假定滑動(dòng)面形式再進(jìn)一步搜索最危險(xiǎn)滑動(dòng)面，然而搜索最危險(xiǎn)滑動(dòng)面較為繁瑣，其計(jì)算效率主要取決于假定的滑動(dòng)面形式是否合理;②未能綜合考慮邊坡坡度、內(nèi)摩擦角及粘聚力等邊坡參數(shù)對最優(yōu)樁位及樁長的影響;③未能全面考慮抗滑樁樁位及樁長對最危險(xiǎn)滑動(dòng)面的變化規(guī)律.

基于此，本文首先采用有限元極限分析軟件OptumG2建立抗滑樁加固前后的數(shù)值模型，并通過算例驗(yàn)證了運(yùn)用本文數(shù)值分析方法的可行性;然后，探討邊坡坡度及抗剪強(qiáng)度對安全系數(shù)、最優(yōu)樁位及樁長的影響;最后，分析了抗滑樁加固后的最危險(xiǎn)滑動(dòng)面變化規(guī)律，以期為邊坡處治工程設(shè)計(jì)提供一定的參考.

1 ? OptumG2計(jì)算邊坡穩(wěn)定性的基本原理

如上所述，在抗滑樁-滑坡系統(tǒng)穩(wěn)定性分析中，經(jīng)典極限分析方法適用于工程設(shè)計(jì)，但需事先知道破壞面的位置;有限元法適應(yīng)性廣，但無法量化穩(wěn)定性，求解其安全系數(shù). 二維有限元極限分析軟件OptumG2與傳統(tǒng)的極限分析方法相比，有限元極限分析不需要在計(jì)算之前假定邊坡破壞形式，而是通過計(jì)算自動(dòng)找出最危險(xiǎn)滑動(dòng)面;而相比有限元法，有限元極限分析只需計(jì)算極限狀態(tài)下的邊坡穩(wěn)定性，因而可以在不降低精度的前提下提高計(jì)算效率，兼有極限分析法與有限元分析法的優(yōu)勢，尤其適用于巖土工程中滑坡及隧道等工程的穩(wěn)定性分析. 考慮邊坡穩(wěn)定性安全系數(shù)有2種方式：1）重力乘數(shù)法. 不斷增大土體重度，并定義其安全系數(shù)FSγ為邊坡剛好失穩(wěn)時(shí)的土體重度γcr和土體實(shí)際重度γ之比，即FSγ = γcr /γ;2）強(qiáng)度折減法. 同比例降低土體強(qiáng)度指標(biāo)c及tanφ，并定義其安全系數(shù)FSS為土體實(shí)際強(qiáng)度指標(biāo)（c、tanφ）和邊坡剛好失穩(wěn)時(shí)的土體強(qiáng)度指標(biāo)（ccr、tanφcr）之比，即FSS = c/ccr = tanφ/tanφcr. 總體而言，基于強(qiáng)度折減的安全系數(shù)比基于重力增大的安全系數(shù)更保守，即強(qiáng)度折減法算得的安全系數(shù)更小. 2種安全系數(shù)計(jì)算方式各有其優(yōu)點(diǎn)，一般認(rèn)為實(shí)際邊坡破壞源于土體強(qiáng)度的降低而非重力的增大，基于強(qiáng)度折減的安全系數(shù)更合理一些;而某些情況（現(xiàn)場或模型試驗(yàn)）下，由于坡頂堆載而引起邊坡的破壞，此時(shí)重力增大的安全系數(shù)可能更符合實(shí)際. OptumG2軟件判斷邊坡破壞的方式為：滑面塑性區(qū)貫通或滑面上位移與應(yīng)變將產(chǎn)生突變，產(chǎn)生很大的且無限制的塑性流動(dòng).

2 ? 邊坡數(shù)值模型的建立及驗(yàn)證

2.1 ? 算例1（無抗滑樁加固的邊坡）

為驗(yàn)證本文數(shù)值分析方法用于未加固邊坡穩(wěn)定性分析的合理性，現(xiàn)引入文獻(xiàn)[2]中的邊坡算例進(jìn)行對比. 假定某均勻土坡的破壞模式為坡腳破壞，數(shù)值模型中的土坡高度H、邊坡坡度β、土體重度γ、內(nèi)摩擦角φ和粘聚力c等參數(shù)見圖1. 根據(jù)圖1的尺寸及參數(shù)建立邊坡數(shù)值模型，見圖2. 土坡采用Mohr-Coulomb材料模擬. 邊界條件選取為標(biāo)準(zhǔn)邊界條件，即模型左右兩端約束橫向位移，模型底部約束兩個(gè)方向的位移. 值得注意的是，OptumG2可采用重力增大及強(qiáng)度折減2種方式計(jì)算安全系數(shù). 通過“工況階段管理器”界面選擇分析類型“極限分析”或“強(qiáng)度折減”可以實(shí)現(xiàn).

采用強(qiáng)度折減極限分析來分析時(shí)，時(shí)間范圍設(shè)置為“長期”，單元類型選為上限或下限（分別用以計(jì)算上限解或下限解），單元數(shù)量設(shè)置為10 000. 采用自適應(yīng)網(wǎng)格，初始單元數(shù)量設(shè)置為1 000，進(jìn)行3次自適應(yīng)迭代，控制變量選為剪切耗散. 在采用重力乘數(shù)極限分析來分析邊坡穩(wěn)定性時(shí)，需將乘數(shù)設(shè)置為“重力”，其余設(shè)置與強(qiáng)度折減極限分析均保持一致.

表1為本文數(shù)值分析結(jié)果與各參考文獻(xiàn)的對比結(jié)果. 對比可知：1）基于強(qiáng)度折減的安全系數(shù)與各參考文獻(xiàn)值吻合較好，且與Bishop法及動(dòng)態(tài)極限分析理論得出的安全系數(shù)非常接近，說明采用強(qiáng)度折減極限分析對邊坡穩(wěn)定性進(jìn)行評價(jià)是合理的;2）總體而言，基于強(qiáng)度折減的安全系數(shù)相對于基于重力增大的安全系數(shù)更小一些，考慮到現(xiàn)有研究及工程多采用強(qiáng)度折減法分析邊坡穩(wěn)定性，因而，采用強(qiáng)度折減極限分析得到的安全系數(shù)偏于安全.

2.2 ? 算例2（經(jīng)抗滑樁加固的邊坡）

為驗(yàn)證本文數(shù)值分析方法用于抗滑樁加固的邊坡穩(wěn)定性分析的合理性，現(xiàn)引入文獻(xiàn)[15]中邊坡算例進(jìn)行對比. 該邊坡參數(shù)、抗滑樁樁長h及樁徑D等如圖3所示.

根據(jù)圖3的尺寸及參數(shù)建立邊坡數(shù)值模型，如圖2所示. 模型建立的過程與算例1幾乎一致，唯一不同的是，需在離坡腳13.7 m處設(shè)置抗滑樁，抗滑樁采用線彈性模擬，如圖4中灰色部分，抗滑樁基本參數(shù)選為：樁長h=25 m，樁徑D=0.91 m，彈性模量Ep =30 000 MPa，泊松比ν=0.2，樁體重度γp = 24 kN/m3. 表2為基于本文數(shù)值分析的抗滑樁加固前后的安全系數(shù)與各參考文獻(xiàn)的對比結(jié)果.

對比可知：1） 抗滑樁可以有效地提高邊坡的穩(wěn)定性;2） 本文基于強(qiáng)度折減極限分析的安全系數(shù)與各理論及數(shù)值分析推導(dǎo)結(jié)果相差不大，說明本文數(shù)值方法用以分析抗滑樁前后的邊坡穩(wěn)定性是合理的.

基于不同邊坡參數(shù)，改變樁位及樁長，求得相應(yīng)的穩(wěn)定性安全系數(shù)，如圖7～圖9所示. 由圖可知：

1）邊坡剪切強(qiáng)度對加固效果的影響.不同粘聚力下的安全系數(shù)隨樁長及樁位變化規(guī)律基本一致，而內(nèi)摩擦角對抗滑樁加固效果的影響較大，當(dāng)內(nèi)摩擦角較小時(shí)（φ=5°），抗滑樁對邊坡的加固效果不明顯，應(yīng)嘗試考慮其他加固方法;隨著內(nèi)摩擦角的增大，最優(yōu)樁位逐漸由中部向中上部移動(dòng).

2）樁長與最優(yōu)樁位的影響. 樁長較小時(shí)，抗滑樁靠近坡腳的安全系數(shù)最大，樁長較長時(shí)，抗滑樁位于邊坡中部的安全系數(shù)最大;樁位相同時(shí)，安全系數(shù)隨樁長增大而逐漸增大，但樁長達(dá)到一定程度時(shí)，繼續(xù)增加樁長，安全系數(shù)將不再增加，因而工程實(shí)際中需結(jié)合邊坡安全系數(shù)要求、工程造價(jià)及施工難度確定樁位.

3）樁位對樁長-安全系數(shù)的影響. 當(dāng)抗滑樁靠近坡頂或坡腳時(shí)，樁長對于安全系數(shù)的影響很小;當(dāng)抗滑樁處于邊坡中部時(shí)，樁長的變化對于安全系數(shù)的影響最為顯著，尤其是邊坡較陡時(shí)（β=45°），靠近坡腳處抗滑樁的加固效果不明顯，當(dāng)抗滑樁樁長較小時(shí)（0～10 m），安全系數(shù)基本未發(fā)生變化.

4 ? 結(jié) ? 論

本文基于有限元極限分析軟件OptumG2建立邊坡數(shù)值模型，并通過算例驗(yàn)證其合理性，進(jìn)而探討各參數(shù)對邊坡穩(wěn)定性的影響，得出以下結(jié)論：

1）本文數(shù)值方法用以分析邊坡穩(wěn)定性是可行

的，基于強(qiáng)度折減極限分析的安全系數(shù)偏于安全. 數(shù)值分析得出的剪切耗散較大區(qū)域可視為最危險(xiǎn)滑動(dòng)面.

2）粘聚力對安全系數(shù)隨樁長及樁位變化規(guī)律的影響不明顯. 內(nèi)摩擦角對抗滑樁加固效果的影響較大. 內(nèi)摩擦角較小時(shí)，抗滑樁對邊坡的加固效果不明顯，可嘗試考慮其他加固方法;隨著內(nèi)摩擦角的增大，最優(yōu)樁位逐漸由中部向中上部移動(dòng).

3）樁長較小時(shí)，最優(yōu)樁位靠近坡腳;樁長較長時(shí)，最優(yōu)樁位處于邊坡中部. 當(dāng)抗滑樁樁長較小時(shí)（0～10 m），安全系數(shù)基本未發(fā)生變化，隨著樁長的增加，安全系數(shù)逐漸增大，但樁長達(dá)到一定程度時(shí)，繼續(xù)增加樁長，安全系數(shù)將不再增加.

4）當(dāng)抗滑樁靠近坡頂或坡腳時(shí)，樁長對于安全系數(shù)的影響很小;當(dāng)抗滑樁處于邊坡中部時(shí)，樁長的變化對于安全系數(shù)的影響最為顯著. 尤其是邊坡較陡時(shí)（β = 45°），靠近坡腳處抗滑樁的加固效果不明顯.

5）當(dāng)樁長較小時(shí)，滑動(dòng)面大多未穿過抗滑樁而是繞過抗滑樁底部;隨著樁長的不斷增大，最危險(xiǎn)滑動(dòng)面逐漸向邊坡內(nèi)部移動(dòng);最終，根據(jù)樁位的不同，將出現(xiàn)3種最危險(xiǎn)滑動(dòng)面情況.

參考文獻(xiàn)

[1] ? ?肖專文，張奇志，梁力，等. 遺傳進(jìn)化算法在邊坡穩(wěn)定性分析中的應(yīng)用[J]. 巖土工程學(xué)報(bào)，1998，20（1）：44—46.

XIAO Z W，ZHANG Q Z，LIANG L，et al. Application of genetic evolutionary algorithm for slope stability analysis[J]. Chinese Journal of Geotechnical Engineering，1998，20（1）：44—46. （In Chinese）

[2] ? ?CAO J，ZAMAN M M. Analytical method for analysis of slope stability[J]. International Journal for Numerical and Analytical Methods in Geomechanics， 1999， 23（5）：439—449.

[3] ? ?王國強(qiáng)， 吳道祥， 劉洋， 等. 巢湖鳳凰山滑坡形成機(jī)制和穩(wěn)定性分析[J]. 巖土工程學(xué)報(bào)， 2002， 24（5）：645—648.

WANG G Q， WU D X， LIU Y， et al. Analysis of the stability of mountain Phoenix Landslide in Chaohu city[J]. Chinese Journal of Geotechnical Engineering， 2002， 24（5）：645—648. （In Chinese）

[4] ? ?馬建勛， 賴志生， 蔡慶娥， 等. 基于強(qiáng)度折減法的邊坡穩(wěn)定性三維有限元分析[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2004， 23（16）：2690—2693.

MA J X， LAI Z S， CAI Q E， et al. 3D FEM analysis of slope stability based on strength reduction method[J]. Chinese Journal of Rock Mechanics and Engineering， 2004， 23（16）：2690—2693. （In Chinese）

[5] ? ?郭明偉， 李春光， 王水林. 基于有限元應(yīng)力的三維邊坡穩(wěn)定性分析[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2012， 31（12）：2494—2500.

GUO M W， LI C G， WANG S L. Three-dimensional slope stability analysis based on finite element stress[J]. Chinese Journal of Rock Mechanics and Engineering， 2012， 31（12）：2494—2500. （In Chinese）

[6] ? ?周芬， 郭奧飛， 杜運(yùn)興. 邊坡穩(wěn)定可靠度分析的新型四階矩法[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2016， 43（5）：113—119.

ZHOU F， GUO A F， DU Y X. A new fourth-order moment method for slope stability reliability analysis[J]. Journal of Hunan University（Natural Sciences）， 2016， 43（5）：113—119. （In Chinese）

[7] ? ?李亮， 李彥軍， 趙煉恒， 等. 土石混合邊坡塊石隨機(jī)生成方法與穩(wěn)定性分析[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2017， 44（7）：170—178.

LI L， LI Y J， ZHAO L H， et al. Random generation method and stability analysis of earth-rock mixed slope block stone[J]. Journal of Hunan University（Natural Sciences）， 2017， 44（7）：170—178. （In Chinese）

[8] ? ?戴自航， 沈蒲生. 抗滑樁內(nèi)力計(jì)算懸臂樁法的改進(jìn)[J].湖南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2003，30（3）：81—85.

DAI Z H， SHEN P S. Improvements on calculation of internal forces of cantilev er anti-sliding piles[J]. Journal of Hunan University（Natural Sciences）， 2003，30（3）：81—85. （In Chinese）

[9] ? ?陳昌富， 杜翠翠， 張根寶. 基于雙參數(shù)法剛性抗滑樁嵌固深度可靠性設(shè)計(jì)[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2014， 41（2）：40—46.

CHEN C F， DU C C， ZHANG G B. Reliability design of the embedded depth of the rigid anti-slide piles based on bi-paramater method[J]. Journal of Hunan University（Natural Sciences）， 2014， 41（2）：40—46. （In Chinese）

[10] ?孫書偉， 朱本珍， 馬惠民， 等. 微型樁群與普通抗滑樁抗滑特性的對比試驗(yàn)研究[J]. 巖土工程學(xué)報(bào)， 2009， 31（10）：1564—1570.

SUN S W， ZHU B Z， MA H M， et al. Model tests on anti-sliding mechanism of micropile groups and anti-sliding piles[J]. Chinese Journal of Geotechnical Engineering， 2009， 31（10）：1564—1570. （In Chinese）

[11] ?戴自航， 張曉詠， 鄒盛堂， 等. 現(xiàn)場模擬水平分布式滑坡推力的抗滑樁試驗(yàn)研究[J]. 巖土工程學(xué)報(bào)， 2010， 32（10）：1513—1518.

DAI Z H， ZHANG X Y， ZOU S T， et al. Field modeling of laterally distributed landslide thrusts over anti-slide piles[J]. Chinese Journal of Geotechnical Engineering， 2010， 32（10）：1513—1518. （In Chinese）

[12] ?劉洪佳， 門玉明， 李尋昌， 等. 懸臂式抗滑樁模型試驗(yàn)研究[J]. 巖土力學(xué)， 2012， 33（10）：2960—2966.

LIU H J， MEN Y M， LI X C， et al. Study of model test on cantilever anti-slide pile[J]. Rock and Soil Mechanics， 2012， 33（10）：2960—2966. （In Chinese）

[13] ?傅翔， 謝強(qiáng)， 張永興， 等. 全埋式抗滑樁傾覆破壞的室內(nèi)模型試驗(yàn)研究[J]. 巖土力學(xué)， 2014， 35（8）：2205—2211，2218.

FU X， XIE Q， ZHANG Y X， et al. Model experimental study of toppling failure of fully-embedded anti-sliding pile[J]. Rock and Soil Mechanics， 2014， 35（8）：2205—2211，2218. （In Chinese）

[14] ?章為民， 顧行文， 任國峰， 等. 剛性抗滑樁的樁土相互作用機(jī)理與極限阻滑力研究[J]. 土木工程學(xué)報(bào)， 2017， 50（10）：86—94.

ZHANG W M， GU X W， REN G F， et al. Pile-soil interaction mechanism and ultimate lateral load on rigid stabilizing pile[J]. China Civil Engineering Journal， 2017，50（10）：86—94. （In Chinese）

[15] ?HASSIOTIS S，CHAMEAU J L，GUNARATNE M. Design method for stabilization of slopes with piles[J]. Journal of Geotechnical & Geoenvironmental Engineering，1997，125（10）：314—323.

[16] ?AUSILIO E，CONTE E，DENTE G. Stability analysis of slopes reinforced with piles[J]. Computers & Geotechnics，2001，28（8）：591—611.

[17] ?李家平，賴允瑾，李永盛，等. 利用抗滑樁加固臨水岸坡的邊坡穩(wěn)定性分析[J]. 水利水運(yùn)工程學(xué)報(bào)，2005，27（2）：53—58.

LI J P，LAI Y J，LI Y S，et al. Stabilization analysis of river-side slope reinforced by stabilizing piles[J]. Hydro-science and Engineering，2005，27（2）：53—58. （In Chinese）

[18] ?年廷凱，欒茂田，楊慶. 阻滑樁加固土坡穩(wěn)定性分析與樁基的簡化設(shè)計(jì)[J]. 巖石力學(xué)與工程學(xué)報(bào)，2005，24（19）：29—35.

NIAN T K，LUAN M T，YANG Q. Stability analysis of slopes with stabilizing piles and their simplified design[J]. Chinese Journal of Rock Mechanics and Engineering，2005，24（19）：29—35. （In Chinese）

[19] ?梁冠亭，陳昌富，朱劍鋒，等. 基于M-P法的抗滑樁支護(hù)邊坡穩(wěn)定性分析[J]. 巖土力學(xué)，2015，36（2）：451—456，469.

LIANG G T，CHEN C F，ZHU J F，et al. Stability analysis of pile stabilized slope based on Morgenstern-Price method[J]. Rock and Soil Mechanics，2015，36（2）：451—456，469. （In Chinese）

[20] ?雷文杰，鄭穎人，馮夏庭. 滑坡治理中抗滑樁樁位分析[J]. 巖土力學(xué)，2006，27（6）：950—954.

LEI W J，ZHEN Y R，F(xiàn)ENG X T. Analysis of pile location on landslide control[J]. Rock and Soil Mechanics，2006，27（6）：950—954. （In Chinese）

[21] ?李榮建，于玉貞，李廣信. 抗滑樁加固非飽和土邊坡三維穩(wěn)定性分析[J]. 巖土力學(xué)，2008，29（4）：968—972.

LI R J，YU Y Z，LI G X. 3D global stability analysis of unsaturated soil slope reinforced with piles[J]. Rock and Soil Mechanics，2008，29（4）：968—972. （In Chinese）

[22] ?胡新麗，張永忠，李長冬，等. 庫水位波動(dòng)條件下不同樁位抗滑樁抗滑穩(wěn)定性研究[J]. 巖土力學(xué)，2011，32（12）：3679—3684.

HU X L，ZHANG Y Z，LI C D，et al. Anti-sliding stability of different pile positions under reservoir water level fluctuation[J]. Rock and Soil Mechanics，2011，32（12）：3679—3684. （In Chinese）

[23] ?陳樂求，楊恒山，林杭. 抗滑樁加固邊坡穩(wěn)定性及影響因素的有限元分析[J]. 中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，42（2）：490—494.

CHEN L Q，YANG H S，LIN H. Finite element analysis for slope stability and its influencing factors with pile reinforcement[J]. Journal of Central South University （Science and Technology），2011，42（2）：490—494. （In Chinese）

[24] ?王聰聰，李江騰，廖峻，等. 抗滑樁加固邊坡穩(wěn)定性分析及其優(yōu)化[J]. 中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，46（1）：231—237. （In Chinese）

WANG C C，LI J T，LIAO J，et al. Stability analysis of slope reinforced with piles and optimization[J]. Journal of Central South University （Science and Technology），2015，46（1）：231—237.

[25] ?BISHOP A W. The use of slip circle in the stability analysis of earth slopes[J]. Geotechnique，1955，5（1）：7—17.

[26] ?HUANG S L，YAMASAKI K. Slope failure analysis using local minimum factor-of-safety approach[J]. Journal of Geotechnical Engineering，ASCE 1993，119（12）：1974—1987.

[27] ?HULL T S，POULOS H G. Discussion of “Design method for stabilization of slopes with piles” by T. S. Hull and H. G. Poulos[J]. Journal of Geotechnical & Geoenvironmental Engineering，1999，125（10）：911—913.

[28] ?LI R P，CAI A G，SHENG C G. Stability analysis of pile-slope interaction by finite element method[J]. Journal of China Three Gorges University （Natural Sciences），2007，27（5）：396—401.

[29] ?年廷凱，欒茂田，楊慶，等. 基于強(qiáng)度折減彈塑性有限元法的抗滑樁加固邊坡穩(wěn)定性分析[J]. 巖土力學(xué)，2007，28（S1）：558—562.

NIAN T K，LUAN M T，YANG Q，et al. Stability analysis of slope reinforced with piles by using strength reduction FEM[J]. Rock and Soil Mechanics，2007，28（S1）：558—562. （In Chinese）