樁板墻結構加固邊坡的穩定性分析

羅 渝， 許 強， 何思明， 何盡川

（1.中國科學院山地災害與地表過程重點實驗室，四川成都610041；2.中國科學院水利部成都山地災害與環境研究所，四川成都610041；3.成都理工大學地質災害防治與地質環境保護國家重點實驗室，四川成都610059；4.中國市政工程西南設計研究總院，四川成都610081）

樁板墻結構加固邊坡的穩定性分析

羅 渝1，2，3， 許 強3， 何思明1，2， 何盡川4

（1.中國科學院山地災害與地表過程重點實驗室，四川成都610041；2.中國科學院水利部成都山地災害與環境研究所，四川成都610041；3.成都理工大學地質災害防治與地質環境保護國家重點實驗室，四川成都610059；4.中國市政工程西南設計研究總院，四川成都610081）

為探討樁板墻加固邊坡的穩定性，根據多塊體滑移理論，采用極限分析上限定理，導出了為保證邊坡穩定樁板墻結構中抗滑樁應提供的抗力的計算公式，提出了樁板墻結構加固邊坡穩定性的一種新的分析方法.用該方法對算例進行了分析，并與旋轉破壞模型獲得的結果進行了比較.研究結果表明：樁板墻結構中，抗滑樁應提供的抗力隨邊坡坡角增大而增大，隨坡體材料抗剪強度增大而減小；邊坡潛在滑移面埋深隨坡體抗剪強度增大變淺；按新方法得到的潛在滑移面比基于旋轉破壞模型獲得的潛在滑移面更危險.

邊坡穩定性；極限分析；多塊體滑移理論；樁板墻

樁板墻是一種新型的邊坡加固結構，代替了高大的重力式擋墻，具有良好的加固效果［1］.樁板墻結構由抗滑樁和擋土板兩部分組成，結構的主要受力構件是抗滑樁.通過抗滑樁提供的抗力提高邊坡穩定性，確保邊坡不發生破壞.擋土板的主要作用則是將抗滑樁連接起來，同時攔截抗滑樁間的松散土體，把樁間土壓力和推力傳遞到抗滑樁上，并經由抗滑樁將這些力傳到穩定的地層中.由此可見，樁板墻結構加固邊坡穩定與否主要由抗滑樁決定.

目前用于邊坡穩定性分析的方法較多，極限分析方法在近幾十年來得到了許多學者的關注，有較多研究［2-13］.在抗滑樁結構加固邊坡穩定性的研究方面，Ausilio等提出了利用極限分析上限定理計算抗滑樁加固邊坡穩定性的計算方法［4］；羅渝等利用極限分析上限定理，導出了地震作用下樁錨組合結構加固邊坡時抗滑樁需提供的抗力的計算公式和加固邊坡屈服加速度的計算公式，研究了地震作用下樁錨組合結構加固邊坡的穩定性［11］.然而，對抗滑樁加固邊坡穩定性的上述研究都建立在假設邊坡的潛在滑移面形狀為對數螺旋面的基礎上，而實際潛在滑移面不一定是對數螺旋面，因此，與實際情況可能存在一定差距.

本文基于多塊體滑移理論，用一系列直線將樁板墻結構加固邊坡的滑體土劃分成若干個三角形土體，這一系列三角形土體底邊組合而成的復合直線作為邊坡的潛在滑移面；采用極限分析上限定理，分析各塊體的外力功率和內能耗散，進而導出樁板墻結構中抗滑樁應提供抗力的計算公式；采用數學優化方法，確定邊坡的潛在滑移面和作用在抗滑樁上的抗力，分析樁板墻結構加固邊坡的穩定性，為樁板墻結構設計提供理論依據.

1 分析方法

1.1 基于多塊體滑移理論的速度場計算

假設樁板墻結構加固邊坡潛在滑移面的剪出口位于坡腳.基于多塊體滑移理論［5，14］，以邊坡頂點為這一系列直線的起點，假設這些直線將滑移面劃分為n段，因此，邊坡滑體土被分為n個三角形塊體.把各三角形土體看作剛體，沿滑面產生平移運動.設Li為任意一個三角形土體的底邊長，di為腰長，θi為頂角，αi為底角，ξ為邊坡坡度，β為以抗滑樁為底邊的三角形土體的頂角（見圖1）.這樣，當θi和αi確定后，樁板墻加固邊坡的潛在滑移面分布即可確定.

對于樁板墻加固邊坡而言，要對其穩定性進行分析，首先要構造并計算多塊體滑移破壞模式下的速度場［14］（如圖2）.設第1塊土體ACD沿直線CD的速度為v1，與滑移面的夾角為φ，φ為邊坡土體的內摩擦角［5］.同時，假設土體ABC和樁板墻結構中的抗滑樁是在第1塊土體的推動下運動的，因此，具有與該土體一樣的運動速度v1.以vi代表土體i的速度，與滑移面的夾角為φ；以vi，i＋1代表土體i與i＋1之間的相對速度，與塊體界面的夾角也為φ［5，13］.由此，可推導出：

圖1 樁板墻加固邊坡的多塊體滑移破壞模式Fig.1 Mulit-wedge translation failure mechanism for stability analysis of slopes reinforced with sheet pile wall

圖2 樁板墻結構加固邊坡的多塊體速度場Fig.2 Velocity field of the mulit-wedge of a slope reinforced with sheet pile wall

i＋1

根據式（1）和（2），任意三角形土體的速度都可以表示為v1（第1塊土體ACD的速度）的函數.

1.2 樁板墻加固邊坡的外力功率和內能耗散

1.2.1 外力功率計算

根據極限分析上限定理，對于圖1所示的多塊體破壞模式，外力功率包括土體重力所做的外力功率Ps和抗滑樁所做的功率Pf.

首先，計算土體自重力G：

式中：G0為三角形土體ABC的重力；Gi為三角形土體i的重力.

式中：γ為土體的重度；a和d0分別為三角形土體ABC的AB和AC邊邊長.

根據圖1，由幾何關系可得

由此，土體重力所做的功率：

如圖2，三角形土體i的速度vi與x軸負半軸的夾角ηi可由式（11）計算：

抗滑樁所做的功率

式中：F為樁板墻結構作用在抗滑樁上的側向力，即為保證邊坡穩定，抗滑樁所需提供的水平抗力.

需要說明的是，抗滑樁所做的功率是阻止邊坡發生滑移破壞，而土體重力所做的功率是促使邊坡發生滑移破壞，這2個外力的功率對邊坡起的作用是不同的.因此，若土體重力所做的功率為正，抗滑樁所做的功率就為負.整個機構的外力功率

1.2.2 內能耗散率的計算

內能耗散率D為沿滑移面Li和相鄰三角形土體交界面Ri上的能量耗散之和，

其中，

1.3 樁板墻中抗滑樁所需提供抗力及邊坡潛在滑移面的確定

根據極限分析上限定理，任意機動容許的破壞機制，內能損耗率不小于外力功率［13］.由此，將式（10）和（12）代入式（13）得到外力功率，建立其與內能耗散（式（14））的平衡，使外力功率與內能耗散相等，從而求出

綜上所述可見，式（16）中參數β、ηi、di、vi、vi，i＋1、Gi和Li均為αi和θi的函數，其他參數均為常數.也就是說，樁板墻結構中，抗滑樁所需提供的抗力實際上是αi和θi的函數.因此，可以采用等分的思想，將θi定義為一個常數組，θi＝（π－ξ－β）／n.以減少未知數，簡化計算.

這樣，抗滑樁所需提供抗力F就轉化為只是未知參數αi的函數.對于邊坡問題而言，在所有可能的滑動面中，抗滑樁所需提供的使邊坡達到穩定的抗力值越大，相對應的滑面應越危險［4］，即最危險的滑面就是抗滑樁所需提供抗力最大的滑面.因此，問題的實質是求解F＝max［F（α1，α2，…，αn）］.

可以采用數學優化的方法求解，求出滿足式（17）的αi后，代入式（16），即可得樁板墻加固邊坡時，為滿足邊坡穩定，抗滑樁所需提供的最小抗力.

2 算例分析

用本文提出的方法，以圖3所示的樁板墻加固邊坡為實例，用Mathematic數學軟件求解，研究土體抗剪強度對樁板墻結構中抗滑樁所需提供的抗力及相應潛在滑移面的影響.

基于多塊體滑移理論的樁板墻加固邊坡穩定性分析，劃分的三角形土體數越多，計算結果精度應越高，但達到條一定數目時，計算結果將基本穩定［13］.對于本算例，取三角形土體條分數n＝15.

2.1 土體抗剪強度的影響

這里，用無量綱土體抗剪強度參數λcφ描述土體的黏聚力c和內摩擦力φ，定義

圖3 樁板墻加固邊坡穩定性分析算例Fig.3 Stability analysis example of a slope reinforced with sheet pile wall

同時，用無量綱抗力系數K表示樁板墻結構中抗滑樁所需提供的抗力F，定義

圖5為邊坡坡度為60°，不同土體抗剪強度參數λcφ時潛在滑移面.可見，隨土體抗剪強度參數λcφ增大，潛在滑移面埋深變淺.也就是說，坡體材料抗剪強度參數c、φ越大，潛在滑移面埋深越淺.

圖4 不同坡度下K與土體抗剪強度參數λcφ的關系Fig.4 Relationship between parameters K and λcφfor different values of slope gradient ξ

圖5 邊坡坡度ξ＝60°時，不同土體抗剪強度參數λcφ的潛在滑移面Fig.5 Potential sliding surfaces of a slope when ξ＝60°and λcφhas different values

2.2 與旋轉破壞模式計算結果的比較

為驗證方法的正確性，將本文方法的結果與旋轉破壞模式［3］的結果進行比較，見圖6和圖7.

圖6 邊坡坡度ξ＝60°時，2種方法K與λcφ關系的比較Fig.6 Comparison of relationships between K and λcφobtained by two methods when ξ＝60°

圖7 邊坡坡度ξ＝60°，λcφ＝0.05時，2種方法潛在滑移面的比較Fig.7 Comparison of potential sliding surfaces obtained by the two methods when ξ＝60°and λcφ＝0.05

可見，（1）本文方法得到的抗力值參數K比基于旋轉破壞模式的K值大，這說明本文方法得到的潛在滑移面比旋轉破壞模式的潛在滑移面更危險.這可以理解為，旋轉破壞模式假設潛在滑移面為對數螺旋面，而本文方法并未假設滑移面的形狀，因而可以是任何可能的形狀，當然也包括對數螺旋面，因此可以搜索到比對數螺旋面更危險的滑移面.（2）本文方法得到的潛在滑移面比旋轉破壞模式下的潛在滑移面埋深淺.

3 結 論

（1）基于多塊體滑移理論，提出了樁板墻結構加固邊坡穩定性分析的一種新方法.

（2）樁板墻結構中抗滑樁所需提供的抗力隨邊坡坡角增大而增大，隨坡體材料抗剪強度（參數）c、φ增大而減小；邊坡潛在滑移面的埋深隨土體抗剪強度參數λcφ增大而變淺.

（3）本文方法獲得的潛在滑移面比旋轉破壞模式下的潛在滑移面更危險，驗證了本文方法的正確性.

［1］ 李海光.新型支擋結構設計與工程實例［M］.北京：人民交通出版社，2004：1-9.

［2］ MICHALOWSKI R L.Slope stability analysis：a kinematical approach［J］.Geotechnique，1995，45（2）：283-293.

［3］ CHEN W F.極限分析與土體塑性［M］.詹世斌，譯.北京：人民交通出版社，1995：28-60.

［4］ AUSILIO E，CONTE E，DUENTE G.Stability analysis of slopes reinforced with piles［J］.Computers and Geotechnics，2001，28（8）：591-611.

［5］ DONALD I，CHEN Z Y.Slope stability analysis by the upper bound approach：fundamentals and methods［J］.Canadian Geotechnical Journal，1997，34（6）：853-862.

［6］ 劉小麗，周德培.有軟弱夾層巖體邊坡的穩定性評價［J］.西南交通大學學報，2002，47（4）：382-386.LIU Xiaoli，ZHOU Depei.Stability evaluation of rock mass slope with weak intercalated layers［J］.Journal of Southwest Jiaotong University，2002，47（4）：382-386.

［7］ JEONG S，KIM B，WON J，et al.Uncoupled analysis of stabilizing piles in weathered slopes［J］.Computers and Geotechnics，2003，30（8）：671-682.

［8］ LI X P，WU Y，HE S M.Seismic stability analysis of gravity retaining walls［J］.Soil Dynamics and Earthquake Engineering，2010，30（10）：875-878.

［9］ HE S M，OUYANG C J，LUO Y.Seismic stability analysis of soil nail reinforced slope using kinematic approach of limit analysis［J］.Environmental Earth Sciences，2012，66（1）：319-326.

［10］ LI X P，PEI X J，GUTIERREZ M，et al.Optimal locationofpilesinslopestabilizationbylimit analysis［J］.Acta Geotechnica，2012，7（3）：253-259.

［11］ 羅渝，何思明，歐陽朝軍，等.地震作用下樁錨組合結構加固邊坡穩定性分析［J］.四川大學學報：工程科學版，2010，42（增刊1）：93-99.LUO Yu，HE Siming，OUYANG Chaojun，et al.Stability analysis of pile and anchor composite structure reinforced slope under the earthquake loading［J］.Journal of Sichuan University：Engineering Science Edition，2010，42（Sup.1）：93-99.

［12］ 何思明，張曉曦，歐陽朝軍.基于非線性破壞準則超前支護樁加固高切坡的靜動穩定分析［J］.工程力學，2011，28（12）：119-125.HE Siming，ZHANG Xiaoxi，OUYANG Chaojun.Static-dynamic stability of pre-piled for high cutting slope with nonlinear failure criterion［J］.Engineering Mechanics，2011，28（12）：119-125.

［13］ 羅渝，許強，何思明，等.重力式擋墻加固邊坡的動力穩定性分析［J］.工程力學，2013，30（9）：132-138.LUO Yu，XU Qiang，HE Siming，et al.Stability analysis of slopes reinforced with gravity walls under the earthquake loading［J］.Engineering Mechanics，2013，30（9）：132-138.

［14］ 陳祖煜.土力學經典問題的極限分析上、下限解［J］.巖土工程學報，2002，24（1）：1-11.CHEN Zuyu.Limit analysis for the classic problems of soil mechanics［J］.Chinese Journal of Geotechnical Engineering，2002，24（1）：1-11.

（中、英文編輯：付國彬）

Stability Analysis of Slopes Reinforced with Sheet Pile Wall

LUO Yu1，2，3， XU Qiang3， HE Siming1，2， HE Jinchuan4
（1.Key Laboratory of Mountain Hazards and Earth Surface Process，CAS，Chengdu 610041，China；2.Institute of Mountain Hazards and Environment，CAS，Chengdu 610041，China；3.State Key Laboratory of Geohazard Prevention and Geoenvironment Protection，Chengdu University of Technology，Chengdu 610059，China；4.Southwest Municipal Engineering Design＆Research Institute of China，Chengdu 610081，China）

In order to analyze the stability of a slope reinforced with sheet pile wall，based on the multi-wedge translation mechanism the upper bound method was applied to derive the expressions to calculate needed lateral force acting on piles to guarantee the stability of a slope.A new method was presented to analyze the stability of slopes reinforced with sheet pile wall.In addition，an example was given to illustrate the presented method，and the results were compared with those obtained using the rotation failure mechanism.The research results show that needed lateral force acting on piles increases as slope angle increases，and decrease as soil strength increases.The potential sliding surface becomes increasingly shallow as soil strength increases，and the potential sliding surface obtained using the presented method is more critical than that obtained using the rotation failure mechanism.

slope stability；limit analysis；multi-wedge translation mechanism；sheet pile wall

U418.52

：A

0258-2724（2014）06-0967-05

10.3969／j.issn.0258-2724.2014.06.006

2013-09-28

國家自然科學基金資助項目（41401004）；地質災害與地質環境保護國家重點實驗室開放基金資助項目（SKLGP2012K024）

羅渝（1981－），女，助理研究員，博士，研究方向為地質災害形成機理及其防治結構，E-mail：ly＠imde.ac.cn

何思明（1968－），男，研究員，博士，研究方向為地質災害形成機理及其防治結構，E-mail：hsm112003＠yahoo.com.cn

羅渝，許強，何思明，等.樁板墻結構加固邊坡的穩定性分析［J］.西南交通大學學報，2014，49（6）：967-971.