龍飛 謝杰 張油軍

廣西大學(xué)土木建筑工程學(xué)院，廣西 南寧 530004

基于強(qiáng)度折減法的邊坡穩(wěn)定性分析

龍飛 謝杰 張油軍

廣西大學(xué)土木建筑工程學(xué)院，廣西 南寧 530004

基于ABAQUS軟件在巖土分析計算中的應(yīng)用，采用強(qiáng)度折減法討論了忽略剪脹角的均質(zhì)邊坡和考慮軟弱夾下臥層二種工況下某邊坡的安全系數(shù)。結(jié)果表明：隨著折減系數(shù)的不斷增大并達(dá)到某一數(shù)值時，邊坡內(nèi)塑性應(yīng)變自坡底向坡頂貫通，邊坡達(dá)到極限狀態(tài)，此時的折減系數(shù)即為安全系數(shù)；有限元強(qiáng)度折減法對邊坡穩(wěn)定性分析具有良好的適用性；采用強(qiáng)度折減法無須預(yù)先假設(shè)滑動面位置。

ABAQUS；強(qiáng)度折減；邊坡穩(wěn)定；安全系數(shù)

Abaqus; strength reduction ;Slope stability; safety factor

前言

邊坡是天然斜坡和人工邊坡的總稱。邊坡失穩(wěn)是生態(tài)環(huán)境和工程建設(shè)中經(jīng)常遇到的廣泛且嚴(yán)重的地質(zhì)災(zāi)害之一，給世界各國人民的生命財產(chǎn)和經(jīng)濟(jì)建設(shè)帶來了巨大的損失。因此有必要對邊坡穩(wěn)定性進(jìn)行分析，從而指導(dǎo)工程建設(shè)。

邊坡穩(wěn)定分析一直是巖土工程中的重要研究領(lǐng)域。目前邊坡穩(wěn)定的分析方法主要可以分為兩大類，即極限平衡法和有限元分析方法[1]。在極限平衡分析方法中，以安全系數(shù)來評價邊坡的穩(wěn)定性，其原理簡單，物理意義明確，是最重要、最常用和最直觀的穩(wěn)定性評價指標(biāo)。而在以往的有限元分析中，是沒有辦法直接評價邊坡穩(wěn)定性的，通常都是根據(jù)邊坡的位移場、應(yīng)力場、塑性區(qū)等來間接地評價，或者根據(jù)有限元計算出應(yīng)力分布之后再利用極限平衡分析方法計算出一個安全系數(shù)指標(biāo)，其計算結(jié)果往往難以理解，也很難被工程技術(shù)人員所接受。強(qiáng)度折減法的出現(xiàn)改變了這一格局，它直接通過有限元分析獲得一個安全系數(shù)，不僅保持了有限元在模擬復(fù)雜問題上的優(yōu)點，而且概念明確，結(jié)果直觀，在工程中得到了運(yùn)來越多的應(yīng)用。

1、強(qiáng)度折減法的基本原理

強(qiáng)度折減法最早由Zienkiewicz等提出，后被許多學(xué)者廣泛采用。他們提出了一個抗剪強(qiáng)度折減的概念，其定義為：在外荷載保持不變的情況下，邊坡內(nèi)土體所能提供的最大抗剪強(qiáng)度與外荷載在邊坡內(nèi)所產(chǎn)生的實際剪應(yīng)力之比。在極限狀況下，外荷載所產(chǎn)生的實際剪應(yīng)力與抵御外荷載所發(fā)揮的最低抗剪強(qiáng)度即按照實際強(qiáng)度指標(biāo)折減后所確定的、實際中得以發(fā)揮的抗剪強(qiáng)度相等[2]。當(dāng)假定邊坡內(nèi)所有土體抗剪強(qiáng)度的發(fā)揮程度相同時，這種抗剪強(qiáng)度折減系數(shù)相當(dāng)于傳統(tǒng)意義上的邊坡整體穩(wěn)定安全系數(shù)FS，又稱為強(qiáng)度儲備安全系數(shù)，與極限平衡法中所給出的穩(wěn)定安全系數(shù)在概念上是一致的。強(qiáng)度折減法已經(jīng)有許多的學(xué)者運(yùn)用到邊坡的穩(wěn)定分析之中[3-5]

折減后的抗剪強(qiáng)度參數(shù)可分別表達(dá)為：

2、強(qiáng)度折減法在Abaqus中的實現(xiàn)

3、模型概況及參數(shù)

有一高H=10m，坡腳β=45°的土坡，土體容重γ=20KN/m3，黏聚力c=12.38kPa，摩擦角=20°，分別計算忽略剪脹角的均質(zhì)邊坡（工況1）和考慮軟弱夾下臥層（工況2）。邊坡尺寸模型見圖1、圖2所示, 在本模擬中通過分段直線來模擬，2種工況的具體計算參數(shù)見表1[6]。若按極限平衡法分析，本算例的土坡穩(wěn)定安全系數(shù)為1.0。

圖1 忽略剪脹角均質(zhì)邊坡尺寸模型（工況1）

圖2 軟弱夾下臥層邊坡尺寸模型（工況2）

表1 計算參數(shù)表

以忽略剪脹角均質(zhì)邊坡（工況1）為例，網(wǎng)格劃分如圖3所示,采用的是有限差分網(wǎng)格，單元選取CPE4（四節(jié)點平面應(yīng)變單元）。約束條件為:模型的左側(cè)和右側(cè)約束水平向X位移,底部約束水平X和豎直向Y位移,其它為自由約束。在本文分析中，從眾多模型中選取Mohr-Koulomb模型對邊坡進(jìn)行分析。

圖3 邊坡網(wǎng)格劃分圖（工況1）

4、結(jié)果分析

4.1 忽略剪脹角均質(zhì)邊坡

4.1.1 安全系數(shù)

輸出X方向位移與安全系數(shù)Fs的關(guān)系曲線如圖4，由圖可見，若以數(shù)值計算不收斂作為土坡穩(wěn)定的評價標(biāo)準(zhǔn)，對應(yīng)的安全系數(shù)Fs為1.06，即Fs=1.06。同時注意到頂部節(jié)點水平位移有一個明顯的拐點，若以位移的拐點作為評價標(biāo)準(zhǔn)，則安全系數(shù)為Fs=0.98。這兩個數(shù)值與極限平衡分析方法給出的Fs=1.0相比都比較接近，說明本例都是可行的。

將分析步t=0.294和t=0.320的積分點上的等效塑性應(yīng)變圖分別繪制于圖5和圖6中。從圖5和圖6這兩個圖清楚地表明了土坡失穩(wěn)的過程，即一開始是土坡坡腳出現(xiàn)屈服，然后向上延伸，直到分析步t=0.320時出現(xiàn)了塑性區(qū)的貫通現(xiàn)象，對應(yīng)的安全系數(shù)為Fs=0.98（見圖7），這和前面兩種判斷標(biāo)準(zhǔn)得到的安全系數(shù)差不多，尤其是位移拐點方法并無差別。初步分析的原因是塑性區(qū)貫通后位移自然快速增加，而計算并不一定不收斂。

這條等效塑性帶貫穿這個邊坡，且經(jīng)過坡腳。這也是邊坡即將破壞的重要標(biāo)志。通過這個事例的計算說明Abaqus能夠處理復(fù)雜的邊坡穩(wěn)定性問題。為我們利用數(shù)值解法研究邊坡穩(wěn)定性提供了一個重要途徑。

圖4 位移-安全系數(shù)關(guān)系曲線

圖5 t=0.294的塑性區(qū)

圖6 t=0.320的塑性區(qū)

4.1.2 滑動面

在邊坡穩(wěn)定分析中，還需要確定滑動面的位置。將計算終止時的位移等值線云圖繪制于圖8中。由圖8可以很清楚的判斷出滑動面的位置，與極限平衡分析法中的一樣，呈大致的圓弧狀，并且通過坡角點。

圖7 t=0.320的安全系數(shù)（Fs=Fv=0.98）

圖8 滑坡位移等值線云圖

4.2 軟弱夾下臥層

現(xiàn)模擬的工況為邊坡以下的局部圖層設(shè)為軟弱土層，將邊坡向左和向右均延伸20米。邊坡模型尺寸見圖2及軟弱土層的參數(shù)見表1所示。邊界條件及計算方法同工況1一致，工況2的網(wǎng)格劃分如圖9所示。下面將對模型的計算結(jié)果進(jìn)行分析，如下：

圖9 邊坡網(wǎng)格劃分圖（工況2）

4.2.1 安全系數(shù)

圖10是土層頂部標(biāo)記節(jié)點(見圖9)的X方向位移與安全系數(shù)Fs的關(guān)系曲線（如圖10），由圖10可見，若以數(shù)值計算不收斂作為土坡穩(wěn)定的評價標(biāo)準(zhǔn)，對應(yīng)的安全系數(shù)Fs為0.87，即Fs=0.87。同時注意到頂部節(jié)點水平位移有一個非常明顯的拐點，若以位移的拐點作為評價標(biāo)準(zhǔn)，則安全系數(shù)為Fs=0.86。這兩種方式得到的安全系數(shù)都十分接近，并且也符合邊坡底部夾軟弱層的邊坡的一般規(guī)律相吻合。根據(jù)所假定的條件，該邊坡將發(fā)生滑坡現(xiàn)象，此時邊坡不安全。

圖10 位移與安全系數(shù)關(guān)系曲線

圖11 用增量位移判斷滑動面

4.2.2 滑動面

在邊坡穩(wěn)定分析中，還需要確定滑動面的位置。將計算終止時的位移等值線云圖繪制于圖11中。圖11給出了計算結(jié)束后的位移等值線云圖。計算結(jié)果很好地模擬了有軟弱下臥層的滑動趨勢，即有部分滑動面沿著軟弱下臥層水平方向，整個滑動面呈兩端圓弧中間直線的復(fù)合狀，與一般規(guī)律吻合，這說明強(qiáng)度折減法無須預(yù)先假設(shè)滑動面的位置。由圖可以很清楚的判斷出滑動面的位置，與極限平衡分析法中的一樣，呈大致的圓弧狀，并且通過坡角點。

5、結(jié)語

通過以上分析和計算可以得出以下結(jié)論：

5.1 與極限平衡法相比，有限單元法不需要作任何假定，計算結(jié)果可靠。有限元法能分析各種復(fù)雜形狀、多種材料組成的邊坡，并且不需要事先假定滑動面。

5.2 建立在強(qiáng)度縮小有限元分析基礎(chǔ)上的邊坡穩(wěn)定分析理論，折減系數(shù)本身就是傳統(tǒng)意義上的邊坡穩(wěn)定系數(shù)，通過折減土體材料的強(qiáng)度來分析邊坡的穩(wěn)定性，直到收斂條件崩潰為止，此時的強(qiáng)度折減系數(shù)即為所求的邊坡安全系數(shù)。

5.3 采用建立在強(qiáng)度縮小有限元分析基礎(chǔ)上的邊坡穩(wěn)定分析的基本原理，利用Abaqus提供的非線性彈塑性模型進(jìn)行邊坡穩(wěn)定分析計算是可行的。

[1]劉文平,鄭穎人，劉元雪.邊坡穩(wěn)定性理論及其局限性[J].后勤工程學(xué)院學(xué)報，2005，1

[2]郭愛斌，萬智. 基于強(qiáng)度折減法的紅粘土邊坡穩(wěn)定分析研究[J].湖南交通科技，2010，36(1):36～37

[3]趙尚毅，鄭穎人，張玉芬.極限分析有限元法講座-Ⅱ有限元強(qiáng)度折減法中邊坡失穩(wěn)的判據(jù)討論[J].巖土力學(xué)，2005，26(2)：332～336

[4]奕茂田，武亞軍，年廷凱.強(qiáng)度折減法中邊坡失穩(wěn)的塑性區(qū)判據(jù)及其應(yīng)用[J].防災(zāi)減災(zāi)工程學(xué)報，2003,23(3):1～8

[5]陳祖煜.土質(zhì)邊坡穩(wěn)定分析—原理–方法–程序[M]. 北京:中國水利水電出版社, 2003

[6]費康，張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京：中國水利水電出版社，2010.393～403

Analyze the Slope Stability Based On Strength Reduction Fem

Based on the great geotechnical analyzing function of Abaqus，through strength reduction fem，Discussing the safe factor of the slope under dilatancy ignored uniform slope angle and consider the weak two folder under the conditions underlying layer. The results showed that: With the increasing efficiency factor and achieve a certain value, the slope of the plastic strain from its base to the top of the hill through the slope to reach the limit state, when the reduction factor is the safety factor;Limited strength reduction of the slope stability analysis has a good applicability; The slip surface location should not be assumed when using strength reduction fem.

10.3969/j.issn.1001-8972.2012.06.023

龍飛，男，碩士研究生，巖土工程方向。