極限平衡法與強度折減法計算穩(wěn)定安全系數(shù)的對比分析

魏蕓，袁琳

（中交第一航務(wù)工程勘察設(shè)計院有限公司，天津 300222）

0 引言

邊坡穩(wěn)定性問題一直以來都是巖土工程界的一項重要研究內(nèi)容，二維巖土邊坡穩(wěn)定分析方法及應(yīng)用相對較成熟，在工程上應(yīng)用較為廣泛的是極限平衡法。對于純巖土三維邊坡穩(wěn)定問題，由于極限平衡分析涉及眾多的未知內(nèi)力和繁復(fù)的矢量運算，至今尚無較好的實用方法[1]。有不少學(xué)者對三維邊坡穩(wěn)定問題從基本方程到分析方法上進(jìn)行了研究[2-4]，可見，純巖土邊坡的算法已相當(dāng)復(fù)雜，涉及到結(jié)構(gòu)的巖土邊坡穩(wěn)定分析算法就更為復(fù)雜。隨著有限元技術(shù)及算法的進(jìn)步，運用三維有限元分析樁基穩(wěn)定性的論文也比較多[5-9]，但是，得出來的安全系數(shù)到底可靠度如何，在工程上需要用經(jīng)驗去標(biāo)定。

目前用于邊坡穩(wěn)定分析的方法有：工程類比法、圖解法、極限平衡法、有限元法等。其中應(yīng)用較為廣泛的是建立在剛體極限平衡理論上的剛體極限平衡法和以有限元為代表的數(shù)值分析方法即通常的強度折減系數(shù)法。極限平衡法是將滑體視為剛性體，不考慮其本身的變形，大多數(shù)的破壞都簡化為平面問題。根據(jù)土體沿著假想滑動面上的極限平衡條件進(jìn)行分析，其破壞遵從摩爾-庫倫定律，并認(rèn)為當(dāng)邊坡的穩(wěn)定系數(shù)等于1時，即為邊坡失穩(wěn)的極限平衡狀態(tài)。極限平衡法是工程實踐中應(yīng)用最早且經(jīng)驗積累最多的方法，缺點是必須通過試算，且絕大多數(shù)是平面問題，對于復(fù)雜的三維問題的算法有不少學(xué)者在研究，但推廣度不夠。目前多數(shù)情況下只能求助于大型的有限元軟件，將強度折減一個系數(shù)，該系數(shù)由小到大逐步增加，折減后的抗剪強度逐漸減小，直到發(fā)生整體失穩(wěn)之前的那個折減系數(shù)值，即為土坡的安全系數(shù)。

在實際應(yīng)用中發(fā)現(xiàn)，強度折減法與極限平衡法計算的安全系數(shù)有時相去甚遠(yuǎn)，有必要對該方面的問題進(jìn)行系統(tǒng)的分析總結(jié)，為較準(zhǔn)確地評估復(fù)雜三維問題的安全性奠定基礎(chǔ)。

1 研究內(nèi)容及研究方案

極限平衡法在工程實踐中使用時間長，經(jīng)驗豐富，而強度折減法是三維有限元中普遍采用的方法，但經(jīng)驗缺乏。本研究的目的是以極限平衡法計算的安全系數(shù)作為基準(zhǔn)，標(biāo)定有限元強度折減法安全系數(shù)的計算值，為將來復(fù)雜三維有限元計算安全系數(shù)值的評估應(yīng)用提供參考。主要通過以下對比分析展開研究：

1）不同土質(zhì)的純巖土邊坡采用極限平衡法、二維強度折減法和三維強度折減法計算安全系數(shù)并進(jìn)行差異分析。

2）有樁基結(jié)構(gòu)時極限平衡法、二維強度折減法和三維強度折減法計算安全系數(shù)并進(jìn)行差異分析。

2 有限元計算

2.1 無樁時純土坡安全系數(shù)計算

以下分別針對純土質(zhì)邊坡采用上述3種方法計算，土體分層及參數(shù)見表1（參照天津港某碼頭工程地質(zhì)資料選取），幾何模型見圖1。

極限平衡法、二維強度折減法、三維強度折減法3種方法安全系數(shù)計算結(jié)果依次為1.224、1.261、1.278，以極限平衡法為基準(zhǔn)1，則二維強度折減法和三維強度折減法的安全系數(shù)分別為1.030和1.044。極限平衡法將滑動體視為剛體，滑弧形狀是圓弧形，而其它兩種強度折減法，滑弧位置不固定，隨邊坡坡率及土層分布變化而變化，且三維分析中出現(xiàn)主滑弧影響帶，即滑弧外一定厚度范圍的土體在抗滑過程中發(fā)揮了一部分抗滑力。對于土體而言，存在主滑弧影響帶是比較符合實際的，而極限平衡法無法考慮這部分土的抗滑貢獻(xiàn)力。

表1 土層分布及其參數(shù)Table 1 Distribution of soil layer and its parameters

圖1 幾何模型（m）Fig.1 Geometric model(m)

2.2 有樁時土坡安全系數(shù)計算

當(dāng)土坡穩(wěn)定性不滿足時，通常采用阻滑樁提高土坡的穩(wěn)定性。選擇帶阻滑樁的邊坡作為分析對象，分別用二維等效方法和三維方法計算。如果邊坡太過復(fù)雜，本身在向二維簡化過程中就會有失真，用本身就已失真的模型去分析算法對結(jié)果的影響就會失真更多。因此，這里選擇最簡單的沿長度方向斷面不變的邊坡作為分析對象。

極限平衡法雖然經(jīng)驗豐富，但應(yīng)用于三維計算的成熟的商業(yè)軟件還沒有。本研究的目的在于，通過簡單邊坡在考慮阻滑樁受剪承載力時，隨著樁數(shù)的變化，二維計算結(jié)果與三維計算結(jié)果的差異，分析其原因，從而評估在復(fù)雜三維空間計算中安全系數(shù)的可靠性。

樁型選擇：混凝土管樁，外徑1 000 mm，壁厚130 mm，參照GB 13476—2009《先張法預(yù)應(yīng)力混凝土管樁》附錄C表C.1管樁的抗剪性能，外徑為1 000 mm，型號為A的抗裂剪力為695 kN，假定斷裂面與正截面夾角45毅（面積最大，應(yīng)力最?。?，則每根樁相應(yīng)的抗裂剪應(yīng)力約為1 390 kPa。設(shè)計以下3種計算工況，每隔14 m布置2根樁、7根樁、14根樁，見圖2所示。計算時樁底標(biāo)高插入土層淤1-2黏土中4 m。

圖2 各工況樁數(shù)設(shè)置Fig.2 Number of piles in each working condition

樁按彈性體考慮，用三維強度折減法計算以上3種工況，得到的安全系數(shù)分別為：1.352、1.389和1.578。在強度折減法中，由于彈性材料沒有強度參數(shù)，不會參與強度折減，也即樁不會破壞，因此，只要有樁，破壞時滑弧必定會在樁端以下。然而，實際上，當(dāng)樁數(shù)或樁的抗剪能力不足時，樁會被截斷。

樁體用修正摩爾庫侖模型，可以真實地考慮樁的抗剪作用。分別計算3種工況下的安全系數(shù)為 1.347、1.391、1.600。

用極限平衡法（二維）計算3種工況的安全系數(shù)分別為1.254、1.369和1.369。

用二維強度折減法考慮樁的抗剪能力計算3種工況下安全系數(shù)分別為：1.355、1.397和1.397。

3 分析

各種工況安全系數(shù)匯總見表2。

表2 安全系數(shù)匯總表Table 2 Summary of safety factors

在三維應(yīng)力應(yīng)變有限元分析中，為方便提取樁體內(nèi)力，常常將樁設(shè)為彈性材料，因此三維模型中用彈性樁計算的安全系數(shù)與二維極限平衡法計算結(jié)果的關(guān)系對評價三維模型計算安全系數(shù)有重要意義。將各種工況下二維極限平衡法計算的安全系數(shù)定為基準(zhǔn)1，則三維彈性樁計算的相對安全系數(shù)見表3。

表3 相對安全系數(shù)表Table 3 Relative safety factor

4 討論

根據(jù)上述計算結(jié)果，對比分析各種算法。

1）極限平衡法與強度折減法（二維）相比：純巖土邊坡，強度折減法計算的安全系數(shù)比極限平衡法大3%；有抗滑樁時，強度折減法計算的安全系數(shù)比極限平衡法大，且樁的抗剪能力越高，安全系數(shù)偏大得越多，當(dāng)抗剪能力達(dá)到某一數(shù)值時，通過增加樁的抗剪能力來提高抗滑穩(wěn)定性不再有效果。

2）三維強度折減法與二維強度折減法相比：純巖土邊坡三維強度折減法計算結(jié)果比二維強度折減法大1.4%；有抗滑樁但樁數(shù)較少（2根樁、7根樁）時，二維計算結(jié)果比三維大0.5%左右；有抗滑樁但樁數(shù)較多（14根樁）時，三維計算結(jié)果比二維大14.55%。究其原因是因為二維算法是將樁的抗剪強度均勻折算到土里，但實際上當(dāng)樁數(shù)不足時，樁增強的土的范圍有限，當(dāng)多根樁及其周圍土破壞時，其余范圍的土的強度來不及發(fā)揮作用就會被一起帶動而失穩(wěn)，因此用二維算法將樁的強度均勻折算有時得出的安全系數(shù)失真且偏大。

3）三維計算中樁為彈性和考慮其實際抗剪強度時計算結(jié)果的對比：抗滑樁較少時（2根樁），彈性樁算出的安全系數(shù)偏大0.45%；抗滑樁較多時（14根樁），樁考慮抗剪強度時計算的安全系數(shù)偏大1.4%；當(dāng)抗滑樁數(shù)為7根時，兩種方法算出的安全系數(shù)基本相等。由于彈性樁不會受剪破壞，如考慮樁的抗剪強度計算的安全系數(shù)大于彈性樁計算的安全系數(shù)時，說明邊坡的破壞不受樁的抗剪強度控制（即破壞面不會穿過樁身），此時，用彈性樁計算更為合理。

4）三維彈性樁模型與二維極限平衡法計算結(jié)果的對比：樁數(shù)較少（2根）時，由于滑弧位置不同，三維計算結(jié)果偏大7.8%；當(dāng)樁數(shù)從7根增加到14根時，滑弧位置基本相同，樁數(shù)越多，三維計算結(jié)果越大，在14根樁時偏大15.3%。

5 結(jié)語

本文旨在通過各種算法的比較，評估三維彈性樁計算安全系數(shù)的可靠性，通過分析，可以得到以下結(jié)論：

1）極限平衡法計算安全系數(shù)時，抗滑力只與滑弧經(jīng)過的土單元的強度有關(guān)，而強度折減法求得的滑弧是一個滑動帶，影響范圍更大，滑動帶周圍的土對滑弧及安全系數(shù)都有貢獻(xiàn)，所以任何情況下，極限平衡法求得的安全系數(shù)均最小。

2）同等條件下，當(dāng)樁數(shù)偏少時，二維強度折減法比三維強度折減法計算的安全系數(shù)偏大。因為二維等效方法假定土的強度全部等比例發(fā)揮。實際上，在三維空間中，樁及其周圍的土體承擔(dān)的荷載較大，發(fā)揮的抗滑強度也較大，相比之下，離樁較遠(yuǎn)的土容易被動滑動，其強度不一定能夠完全發(fā)揮。因此，在樁數(shù)不足時，慎用二維等效算法。

3）當(dāng)抗滑樁足夠即樁不會受剪破壞時，三維彈性樁強度折減法求得的安全系數(shù)更真實。在樁數(shù)較多時，考慮樁的抗剪強度比直接用彈性樁計算安全系數(shù)更大，說明群樁在破壞時不一定每一根都能夠達(dá)到最大的抗剪能力，可能是幾根樁破壞帶動土體滑動，其它樁隨即一起滑動。

4）當(dāng)樁數(shù)較少時，例如14延米只有2根樁時，實際的破壞很可能是樁受剪破壞，此時，三維彈性樁計算出的安全系數(shù)大概率是失真的。

因此，在三維邊坡穩(wěn)定分析中，應(yīng)首先判斷邊坡最危險滑弧是否通過樁身。具體方法是分別采用彈性樁和考慮抗剪強度的樁兩種方法計算安全系數(shù)，當(dāng)前者大于后者，說明滑弧可能通過樁身，應(yīng)補樁；當(dāng)后者大于前者，說明有部分樁功能不能發(fā)揮，應(yīng)減少樁的布置；當(dāng)兩者幾乎相等時，可認(rèn)為樁的抗滑功能發(fā)揮到最大，此時的安全系數(shù)應(yīng)該是邊坡的真實安全系數(shù)（在土體參數(shù)真實的情況下）。同時為了評估樁在抗滑中發(fā)揮的作用，還應(yīng)計算在不考慮樁時結(jié)構(gòu)的安全系數(shù)。