閱讀 3063 次 開挖邊坡穩(wěn)定有限元研究

開挖邊坡穩(wěn)定有限元研究

高志輝，蔣和洋

(北京工業(yè)大學建筑工程學院，北京100022)

1.邊坡穩(wěn)定計算方法簡述

    邊坡穩(wěn)定性問題是長久以來巖土工程界一直關注卻叉懸而未決的問題．經(jīng)濟合理地設計邊坡及評價邊坡的穩(wěn)定性．對國民經(jīng)濟建設有著重要意義。目前．邊坡穩(wěn)定分析的方法已有數(shù)十種之多，主要分為兩大類：(I)經(jīng)典方法極限平衡法、極限分析法、滑移線法等；(2)數(shù)值分析方法：有限元法、離散元法、邊界元法等。

    經(jīng)典方法把破壞分析和變形分析分開計算，破壞分析往往基于理想剛塑性假定．而，變形分析卻基于巖土介質(zhì)的平均線彈性假定。如目前應用最廣泛的極限平衡法，它通過極限狀態(tài)時土體外荷載與內(nèi)部強度所提供的抗力之間的靜力平衡關系．計算土體的安全系數(shù)。這種破壞分析將士體視為理想剛塑性體．不考慮變形協(xié)調(diào)與材料的本構(gòu)差系，只假定破壞瞬間的變形機制雖然它有原理簡單．計算方便．直接給出安全系數(shù)等優(yōu)點，但它有較大局限性，對于復雜一些的情況就不適用。有限元法在邊坡分析中有著特殊的優(yōu)勢，它不受幾何形狀的不規(guī)則、邊界條件的多樣性和材料的不均勻性的限制。更重要的是有限元法可以分析變形與應力變化的過程，而不是僅僅進行臨近破壞時的分析，所以有限元法是一種把破壞分析與變形分析相結(jié)臺同時進行的方法，而且隨著工程實踐解決問題的深度和廣度的需要．有限元方法將具有更J'泛而良好的發(fā)展前景。

    與經(jīng)典方法相比．采用有限元法進行土體穩(wěn)定分析的優(yōu)勢可以歸納如下。

    (1)可以對具有復雜的地貌、地質(zhì)的邊坡進行穩(wěn)定性分析：

    (2)可以得q極限狀態(tài)下的失效形式．確定潛在滑動面的大致位置，因而不需要事先假定破壞面的形狀或位置這對于指導施工設計是很重要的；

    (3)由于有限元法引入變形協(xié)調(diào)的車構(gòu)關系，因此也不必引入假定條件就可以得到不同工作狀態(tài)下土體的真實受力狀態(tài)，全面了解應力、變形的分布狀況，使方法保持了嚴密的理論體系；

    (4)可以了解土體隨強度的惡化而呈現(xiàn)出的漸進失穩(wěn)過程，了解土體的薄弱部位；

    (5)可以考慮不同的施工工序?qū)ν馏w穩(wěn)定最終安全度的影啊。天然邊坡、開挖邊坡、填筑邊坡的安全系數(shù)往往是不同的。

2.安全系數(shù)及臨界狀態(tài)的評定

    利用有限元法強大的計算功能．得到工程實用的安全系數(shù)，首先需要解決安全系數(shù)的計算方法和臨界(破壞)狀態(tài)的評定標準問題。

    在邊坡工程中．安全系數(shù)是衡量土體穩(wěn)定的定量指標，安全系數(shù)的概念在工程界有兩種不同的定義方式：一是將荷載乘以，確定對應于結(jié)構(gòu)址于極限狀態(tài)的F值，即為結(jié)構(gòu)的超載安全系數(shù)，實際上是求結(jié)構(gòu)的極限荷載；另一種是基于強度儲備定義的，通常的做法是將材料的強度陳以F，并將F逐漸增大或減小，直到有限元計算結(jié)果表明結(jié)構(gòu)已達到臨界狀態(tài)，對應的，即為結(jié)構(gòu)的安全系數(shù)對于同一體系．這兩種不同概念的安全系數(shù)是不同的。對于土體，目前多趨向于采用強度儲備意義的安全系數(shù)。Duncan指出安全系數(shù)可以定義為使土體剛好達到臨界破壞狀態(tài)時，對土體剪切強度進行折減的程度。

    目前，將有限元計算的成果轉(zhuǎn)化為工程上適用的安全系數(shù)主要有兩類方法一類是建立在滑裂面應力分析基礎上，一類建立在強度折減基礎上。

    建立在滑裂面應力分析基礎上的有限元法，是在邊坡中定義一個潛在的滑裂面，根據(jù)有限元計算得出的應力分布，計算滑裂面上各點的應力水平，然后根據(jù)加權(quán)平均的原則定義安全系數(shù)。

    建立在強度折減基礎上的邊坡穩(wěn)老分析的基本原理是假定外荷載不變，利用公式(1)和公式(2)折減土體的強度指標c和p，然后對土坡進行彈塑性有限元分析，通過不斷地增加折減系數(shù)f．反復分析土坡，直至邊坡達到臨界破壞，此時的折減系數(shù)即為安全系數(shù)F。

    這種強度折剪技術特別適合甩有限元方法實現(xiàn)．能得到隨著材料強度的降低而呈現(xiàn)出的漸進破壞過程，而且大量數(shù)值結(jié)果表明．強度折減有限元法計算得到的安全系數(shù)與極限平衡法幾乎接近。

    有限元邊坡穩(wěn)定分析的另一個關鍵問題是臨界狀態(tài)(或破壞狀態(tài))的判定．目前主要有兩類方法。

    (1)基于解的收斂性．即以力和位移的不收斂作為邊坡失穩(wěn)破壞的標志

    (2)基于塑性區(qū)的開展和分布，以廣義塑性啦變或者廣義剪應變從坡腳到坡頂貫通作為邊坡破壞的標志之一。

3.在ABAQUS上的實現(xiàn)

    有限元計算中導致解不收斂的原因很多。而且在 ABAQUS有限元軟件中，收斂準則(如迭代次數(shù))的定義古有許多人為的因素，通常很難得到一個穩(wěn)定的解而眾多文獻認為，邊坡土體中塑性應變的發(fā)生與發(fā)展表明了土體屈服或破壞的發(fā)生和發(fā)晨的程度，即不可恢復變形的發(fā)展程度，因此筆者建議根據(jù)ABAQUS的等效塑性應變的開展情況來評判邊坡的穩(wěn)定，即認為當邊坡中等效塑性應變上下貫通時．邊坡失穩(wěn)破壞。

    ABAQUS分析程序?qū)Σ牧蠀��?shù)進行改變方法，不適合實現(xiàn)強度的自動折剪．因此需要手動進行折剪．并建議采用：分法求解方案．一般都能夠很陜地得到一個收斂解分析表明．彈性模量對土體塑性應變的開腱無太大影響．而膨脹角影響較大，隨著膨脹角的增大，安全系數(shù)也隨之提高。

4.算例分析

    分別針對天然垂直邊坡與開挖邊坡及填筑邊坡．進行數(shù)值計算．通過分析驗證了所建議分析方法的合理性與實用性。

    4.1天然邊坡

    如圖l所示均質(zhì)直立邊坡(文獻[8])．坡高為4 5 in幾何參數(shù)見圖。土體本構(gòu)摸型采用彈一理想塑性模型．屈服準則為Molu—Coulomb準則，土體的材料參數(shù)如下：土的重度y=20 kN·m。；粘結(jié)力f；20 kPa；內(nèi)摩擦角f=10，剪脹角為口=10。；彈性模量E；10MPR:泊松比H =0 3；刪壓力系數(shù)k=0 5邊坡右側(cè)和坡腳下左傭均為豎向滑動水f約束邊界條件．底面為全部固定邊界條件。取1／2對稱結(jié)構(gòu)計算．有限元網(wǎng)格劃分如圖I所示。采用八節(jié)點四邊形減縮積分單兀有限元計算模型中共含有270個單元、883個節(jié)點對于這一天然邊坡。外荷載僅為重力。

    安全系數(shù)及邊坡臨界滑動面的計算結(jié)果如圖2所示可見折減(安全)系教為1 055時，邊坡內(nèi)等效塑性應變貫通。

    4.2開挖邊坡

    土體材料參數(shù)與以上天然邊坡同．但分三次開挖．每次開挖深度為15 rH，需時1天，開挖寬度為6 m計算得到的安全系數(shù)和臨界滑動面位置加圖3所示

從圖2，圖3的對比中可以看出，開挖邊坡在第二次開挖時．土坡中便已出現(xiàn)一部分等效塑性應變．這說明開挖工序影響了土坡中塑性變形的分布而且可以看出，開挖邊坡破壞時的等效塑性應變分布與天然邊坡相似；但開挖邊坡比同幾何形式下的天然邊坡更安全。

5.結(jié)論

    通過以上分析，可以得出如下的結(jié)論。

    (1)基于ABAQUS的有限元強度折減法不需要作任何假定，可以分析各種復雜形狀的邊坡，而且能考慮邊坡的施工過程的影響，計算結(jié)果可靠。

    (2)采用等效塑性應變的臨界狀態(tài)破壞準則能夠有效評定邊坡的穩(wěn)定，邊坡的破壞過程可通過ABAQUS軟件的后處理功能動態(tài)明確顯示出來，具有工程使用價值。

    (3)開挖等施工工序?qū)�邊坡的安全穩(wěn)定性有一定影響。 在相同材料參數(shù)和幾何形式下，其所得到的安全系數(shù)比天然邊坡要大．因此進行邊坡穩(wěn)定分析建模時上需要充分考慮。

參考文獻：

[1]鄭穎人邊坡穩(wěn)定性理論及其局限性[J1后勤工程學院學報2005(1)：15—19

[2]趙尚毅，鄭潁人,用有限元強度折減法求邊坡穩(wěn)定安全系數(shù)[J]巖±工程學報，2002．24(3)343—346

[3]蔣和洋，唐立民非線性有限元的若干基本問題[J]力學進展．l988(1)

[4]連陣營，韓國城．孔憲京強度折減有限元法研究開挖邊坡的穩(wěn)定性[J巖土工程學報2001 23(4)407—411

(本文來源：陜西省土木建筑學會  文徑網(wǎng)絡：尚雯瀟 尹維維 編輯  文徑 審核)