王斌 張力

(山西省聞合高速公路黃河大橋建設(shè)管理處，山西運(yùn)城 044300)

邊坡穩(wěn)定性是邊坡工程中最重要的控制因素和工程技術(shù)人員最關(guān)心的問題之一。然而實(shí)際工程中，由于邊坡坡體材料本身的空間變異性，以及取樣和實(shí)驗(yàn)過程中人和環(huán)境所造成的誤差，使得實(shí)際試驗(yàn)參數(shù)具有明顯的不確定性和隨機(jī)性。目前，設(shè)計(jì)規(guī)范仍然采用確定性的設(shè)計(jì)思想，以此進(jìn)行的設(shè)計(jì)必然與實(shí)際結(jié)構(gòu)存在著較大的誤差。

本文通過對(duì)影響邊坡穩(wěn)定的有關(guān)因素的分析得出不同參數(shù)取值的隨機(jī)性對(duì)邊坡穩(wěn)定安全系數(shù)的影響。

1 基本原理

邊坡穩(wěn)定分析的常用方法包括極限平衡法、特征曲線法、極限分析法和數(shù)值分析法等。其中數(shù)值分析法克服了傳統(tǒng)方法中對(duì)于滑裂面形狀、土條間作用力等假設(shè)的缺陷且考慮了材料的應(yīng)力應(yīng)變關(guān)系使得計(jì)算結(jié)果的精確性得到了保證。

邊坡工程數(shù)值分析常采用強(qiáng)度折減法評(píng)價(jià)邊坡的安全穩(wěn)定性。Bishop等將土質(zhì)邊坡穩(wěn)定安全系數(shù)定義為沿整個(gè)滑動(dòng)面的抗剪強(qiáng)度與實(shí)際抗剪強(qiáng)度之比，這也是實(shí)際工程中采用的計(jì)算方法。強(qiáng)度折減法的基本思想與此一致，均為邊坡的強(qiáng)度儲(chǔ)備安全度。其基本原理是將土體參數(shù)c，φ同時(shí)除以一個(gè)折減系數(shù)F，得到一組新的 c′，φ′值，即:

然后作為新的材料參數(shù)代入數(shù)值軟件中進(jìn)行試算，依據(jù)某種破壞判據(jù)判定系統(tǒng)是否進(jìn)入極限平衡狀態(tài)，經(jīng)過多次迭代，邊坡出現(xiàn)失穩(wěn)時(shí)對(duì)應(yīng)的F便為所求安全系數(shù)。

邊坡失穩(wěn)判據(jù)的選擇對(duì)強(qiáng)度折減法的實(shí)現(xiàn)是很重要的，目前的邊坡失穩(wěn)判據(jù)主要有以下幾類:

1)收斂性判據(jù)。

采用迭代不收斂作為邊坡失穩(wěn)的標(biāo)志。

2)突變性判據(jù)。

以計(jì)算過程中關(guān)鍵部位發(fā)生位移突變，土體無限移動(dòng)作為邊坡失穩(wěn)的標(biāo)志。

3)塑性變形判據(jù)。

以廣義塑性應(yīng)變或者等效塑性應(yīng)變從坡腳到坡頂貫通作為邊坡破壞的標(biāo)志。本文綜合上述三種判據(jù)評(píng)價(jià)邊坡的穩(wěn)定性。

在進(jìn)行邊坡穩(wěn)定性分析時(shí)，通常將土體假設(shè)成理想彈塑性體，本文選用摩爾—庫侖屈服準(zhǔn)則進(jìn)行穩(wěn)定性分析。

2 安全系數(shù)參數(shù)敏感性分析

文中采用1987年澳大利亞計(jì)算機(jī)應(yīng)用協(xié)會(huì)(ACADS)對(duì)澳大利亞所使用的邊坡穩(wěn)定分析程序進(jìn)行的調(diào)查評(píng)比考核題中的1(a)作為算例:一均質(zhì)邊坡，其尺寸見圖1，其中土體重度 γ=20 kN/m3，彈性模量 E=10 MPa，泊松比 μ =0.25，粘聚力 c=3.0 kPa，內(nèi)摩擦角 φ =20°(參考結(jié)果:簡(jiǎn)化 Bishop條分法為0.993，Janbu 法為 0.978，裁判答案 1.000)。

圖1 邊坡尺寸

根據(jù)張魯渝、鄭穎人等建議，建模時(shí)采用坡腳至左端邊界距離為坡高的1.5倍，坡頂?shù)接叶诉吔缇嚯x為坡高的2.5倍，上下邊界距離為坡高的2倍。

邊坡問題為平面應(yīng)變問題，根據(jù)上述建議建立如圖2所示的準(zhǔn)三維模型，極差分析數(shù)據(jù)見表1。

模型采用2 421個(gè)網(wǎng)格點(diǎn)，9 642個(gè)單元格。邊界條件取為:左右邊界x，z向受約束，下邊界x，y，z向受約束，其余為自由邊界。采用摩爾—庫侖模型計(jì)算出的安全系數(shù)為1.007 81，與裁判答案基本一致。最大不平衡力采樣記錄圖見圖3，坡肩x軸位移采樣記錄圖見圖4，水平方向位移圖見圖5，最大位移出現(xiàn)在坡腳，坡頂坡腳位移貫通，剪應(yīng)變率云圖見圖6。

圖2 準(zhǔn)三維模型

表1 極差分析數(shù)據(jù)表

邊坡問題屬于多因素多水平問題，本文選擇土體的剪脹角、彈性模量和泊松比作為3個(gè)因素，各個(gè)因素選擇三個(gè)水平，通過構(gòu)造正交試驗(yàn)進(jìn)行參數(shù)敏感性分析(見表2)。根據(jù)工程中常用的取值:對(duì)彈性模量E分別取10 MPa，15 MPa，20 MPa，對(duì)泊松比μ 分別取 0.2，0.25，0.3，剪脹角在 0 ＜ ψ ＜ φ 范圍內(nèi)分別取 10°，15°，20°。

圖3 最大不平衡力采樣記錄圖

圖4 坡肩x軸位移采樣記錄圖

由于 F0.05(2，2)=19 且 Fψ=96.695 82 ＞ 19，F(xiàn)μ=25.651 04＞19，故因素ψ，μ作用顯著，E的作用不明顯。這與極差分析的結(jié)果一致，且安全系數(shù)隨著剪脹角的增大逐漸增大。

表2 正交試驗(yàn)方差分析表

圖5 水平方向位移圖

圖6 剪應(yīng)變率云圖

3 結(jié)語

1)與傳統(tǒng)方法極限平衡法相比，運(yùn)用數(shù)值模擬手段基于強(qiáng)度折減法對(duì)邊坡進(jìn)行穩(wěn)定性分析，可得到邊坡坡體內(nèi)各單元的應(yīng)力應(yīng)變圖，能夠更詳細(xì)地了解邊坡內(nèi)應(yīng)力應(yīng)變的變化規(guī)律，為邊坡后期治理方案的選擇及設(shè)計(jì)提供可靠的理論依據(jù)。

2)通過運(yùn)用正交試驗(yàn)方法得出各參數(shù)對(duì)邊坡穩(wěn)定安全系數(shù)的影響程度，能夠?yàn)楣こ谭治稣咛峁┮欢ǖ膮⒖家罁?jù)。

[1]吳 明，博旭東，劉 歡.邊坡穩(wěn)定分析中的強(qiáng)度折減法[J］.土工基礎(chǔ)，2006，20(1):49-52.

[2]張魯渝，鄭穎人，趙尚毅，等.有限元強(qiáng)度折減系數(shù)法計(jì)算土坡穩(wěn)定安全系數(shù)的精度研究[J］.水利學(xué)報(bào)，2003(1):21-27.

[3]欒茂田，武亞軍，年廷凱.強(qiáng)度折減有限元法中邊坡失穩(wěn)的塑性區(qū)判據(jù)及其應(yīng)用[J］.防災(zāi)減災(zāi)工程學(xué)報(bào)，2003，24(4):644-652.

[4]周翠英，劉祚秋，董立國，等.邊坡變形破壞過程的大變形有限元分析[J］.巖土力學(xué)，2003，24(4):644-652.

[5]關(guān)立軍.基于強(qiáng)度折減的土坡穩(wěn)定分析方法的研究[D］.大連:大連理工大學(xué)，2003.

[6]李旭東.FLAC3D在邊坡穩(wěn)定分析中的應(yīng)用[J］.中國水運(yùn)，2008，8(4):77-79.

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