伍韻瑩，王志鵬，孫立宇

（1.廣東省水利電力勘測(cè)設(shè)計(jì)研究院，廣州510635；2.中國(guó)市政工程西北設(shè)計(jì)研究院深圳分院，廣東 深圳518048；3.長(zhǎng)春工程學(xué)院水利與環(huán)境工程學(xué)院，長(zhǎng)春130021）

0 引言

鑒于邊坡穩(wěn)定性問(wèn)題涉及到工程中的許多領(lǐng)域，近年來(lái)邊坡分析方法也越來(lái)越多，特別是隨著巖土材料的彈塑性有限元計(jì)算技術(shù)的發(fā)展，應(yīng)用有限元法分析邊坡穩(wěn)定性得到了普遍推廣。

目前，強(qiáng)度折減彈塑性有限元法是在邊坡穩(wěn)定分析中適用廣泛、前景良好的一種數(shù)值分析方法。它是將強(qiáng)度折減技術(shù)與彈塑性有限元方法相結(jié)合，對(duì)邊坡的穩(wěn)定性進(jìn)行分析，求得邊坡的最小穩(wěn)定安全系數(shù)。但是，影響計(jì)算結(jié)果計(jì)算精度的因素很多，除了有限元法引入的誤差外，還依賴(lài)于所選用的屈服準(zhǔn)則；同時(shí)，在通常的有限元分析程序中，材料參數(shù)是需要事先給定，為了得到邊坡的最小穩(wěn)定安全系數(shù)，需要多次修改給定值，計(jì)算量大，即使采用二分法一般也需要試算多次。

對(duì)于巖土材料來(lái)說(shuō)，與其他屈服準(zhǔn)則相比，Mohr－Coulomb屈服準(zhǔn)則（簡(jiǎn)稱(chēng) M－C準(zhǔn)則）應(yīng)用最廣、時(shí)間最長(zhǎng)；但是，M－C屈服準(zhǔn)則假定材料的破壞與中主應(yīng)力無(wú)關(guān)，而且在主應(yīng)力空間的π平面上屈服面是一個(gè)等邊不等角的六邊形，屈服面存在尖角，在數(shù)值計(jì)算中收斂性較差，計(jì)算困難；Drucker－Prager屈服準(zhǔn)則（簡(jiǎn)稱(chēng)D－P準(zhǔn)則）在π平面上的屈服面不存在尖角，在數(shù)值計(jì)算中收斂性較好，但是需要采用合適的參數(shù)，否則計(jì)算精度較低，可能造成計(jì)算的邊坡最小穩(wěn)定安全系數(shù)較大。

基于此，本文采用ABAQUS軟件中擴(kuò)展的DP模型，結(jié)合 M－C內(nèi)切圓屈服準(zhǔn)則參數(shù)，并使用ABAQUS軟件中的場(chǎng)變量對(duì)材料強(qiáng)度參數(shù)進(jìn)行連續(xù)折減，整個(gè)分析過(guò)程一次完成，既大大減少了計(jì)算工作量，也為邊坡的穩(wěn)定性分析提供較為精確的結(jié)果。

1 數(shù)值分析方法

ABAQUS軟件中提供了擴(kuò)展的D－P模型，擴(kuò)展的D－P模型的屈服面在π平面上不是圓形，屈服面在子午面上包括線性模型、雙曲線模型和指數(shù)模型。擴(kuò)展的D－P模型可以模擬土、巖石等摩擦材料在單調(diào)加載下的材料特性，同時(shí)也可以考慮材料的剪脹性。

經(jīng)典的D－P模型的屈服準(zhǔn)則的表達(dá)式如下：

式中：I1—— 應(yīng)力張量的第一不變量；

J2——應(yīng)力偏量的第二不變量。

α、k與材料參數(shù)相關(guān)。

在關(guān)聯(lián)流動(dòng)法則下推導(dǎo)平面應(yīng)變狀態(tài)時(shí)與MC準(zhǔn)則匹配的D－P準(zhǔn)則，即可得到M－C內(nèi)切圓屈服準(zhǔn)則（簡(jiǎn)稱(chēng) DP1）的參數(shù)關(guān)系如下［1－2］：

式中：c、θ——M－C 準(zhǔn)則中材料的凝聚力和內(nèi)摩擦角。

本文采用 ABAQUS軟件中線性的D－P模型［3］，屈服準(zhǔn)則的表達(dá)式如下：

式中：β——材料的摩擦角；

d——材料的凝聚力，其值與輸入的硬化參數(shù)有關(guān)，可以分別進(jìn)行計(jì)算。

當(dāng)硬化參數(shù)由單軸壓縮試驗(yàn)參數(shù)σc（屈服應(yīng)力）定義時(shí)：

當(dāng)硬化參數(shù)由單軸拉伸試驗(yàn)參數(shù)σt（屈服應(yīng)力）定義時(shí)：

式中：K——三軸拉伸與壓縮屈服應(yīng)力之比。

聯(lián)立式（1）～ 式（4）可得：

式（6）即為ABAQUS軟件中在關(guān)聯(lián)流動(dòng)法則平面應(yīng)變狀態(tài)時(shí)的參數(shù)轉(zhuǎn)換式。

在ABAQUS軟件中，所謂的場(chǎng)變量只是一個(gè)中間量，材料參數(shù)可以定義為場(chǎng)變量的函數(shù)，以便在問(wèn)題的分析過(guò)程中對(duì)模型中各部分乃至全部區(qū)域進(jìn)行參數(shù)的調(diào)整。例如溫度場(chǎng)是其中最簡(jiǎn)單的一個(gè)場(chǎng)變量，此溫度場(chǎng)只是一個(gè)變量場(chǎng)，不代表真實(shí)溫度，只起到帶動(dòng)材料參數(shù)變化的作用，計(jì)算原理詳見(jiàn)文獻(xiàn)［4］。

抗剪強(qiáng)度折減系數(shù)的概念是由Zienkiewicz等［5］1975年在土工彈塑性有限元數(shù)值分析一文中首次提出的，由此所確定的強(qiáng)度儲(chǔ)備安全系數(shù)與Bishop［6］在極限平衡法中所給的穩(wěn)定安全系數(shù)在概念上是一致的。邊坡的安全系數(shù)定義為把強(qiáng)度參數(shù)減小到邊坡臨界破壞時(shí)的強(qiáng)度參數(shù)折減系數(shù)，強(qiáng)度參數(shù)如下式進(jìn)行折減：

彈性模量E、泊松比μ在計(jì)算中假設(shè)為定值，不隨c、θ值的改變而變化。

2 失穩(wěn)判據(jù)

采用強(qiáng)度折減有限元法分析邊坡穩(wěn)定性的一個(gè)關(guān)鍵問(wèn)題就是如何根據(jù)有限元計(jì)算結(jié)果來(lái)判別邊坡是否處于整體破壞狀態(tài)。目前，邊坡失穩(wěn)破壞的判據(jù)主要有以下幾點(diǎn)［7］。

2.1 數(shù)值收斂判據(jù)

根據(jù)有限元解的收斂性確定失穩(wěn)狀態(tài)，即在給定的非線性迭代次數(shù)極限值條件下，最大位移或不平衡力的殘差值不能滿足所要求的收斂條件，則認(rèn)為邊坡在所給定的強(qiáng)度折減系數(shù)下失穩(wěn)破壞。

2.2 位移判據(jù)

根據(jù)計(jì)算域內(nèi)某一部位的位移與折減系數(shù)之間關(guān)系曲線的變化特征確定失穩(wěn)狀態(tài)，如當(dāng)折減系數(shù)增大到某一特定值時(shí)，某一部位位移突然迅速增大，則認(rèn)為邊坡發(fā)生失穩(wěn)。

2.3 塑性區(qū)判據(jù)

通過(guò)域內(nèi)廣義剪應(yīng)變等某些物理量的變化和分布來(lái)判斷，如當(dāng)域內(nèi)的塑性區(qū)（或某一幅值的廣義剪應(yīng)變）連通時(shí)，則判斷邊坡發(fā)生破壞。理論上，邊坡的變形過(guò)程總是伴隨著一些物理量的出現(xiàn)和發(fā)展，如塑性區(qū)、塑性應(yīng)變、廣義剪應(yīng)變和應(yīng)力水平等，當(dāng)這些物理量達(dá)到一定的值時(shí)，邊坡失穩(wěn)。

本文主要針對(duì)位移判據(jù)和塑性區(qū)判據(jù)這兩種判據(jù)進(jìn)行分析，數(shù)值收斂判據(jù)僅起輔助效果，進(jìn)行簡(jiǎn)單闡述。

3 算例分析

為了便于討論，這里選用趙尚毅等［8－9］所采用的算例邊坡作為計(jì)算對(duì)象，如圖1所示。選用與文獻(xiàn)［8－9］相同的計(jì)算參數(shù)，坡高H＝20m，坡角β＝45°，土的重度γ＝20kN／m3。對(duì)于土體，采用服從M－C屈服準(zhǔn)則與關(guān)聯(lián)流動(dòng)法則的理想彈塑性本構(gòu)模型以及服從DP1準(zhǔn)則與關(guān)聯(lián)流動(dòng)法則的理想彈塑性本構(gòu)模型進(jìn)行計(jì)算，土的黏聚力與內(nèi)摩擦角分別為c＝42kPa和θ＝17°。這里假定彈性模量E＝100MPa，泊松比v＝0.3。

邊坡兩側(cè)為水平約束邊界，底部為完全固定邊界。采用四節(jié)點(diǎn)四邊形單元進(jìn)行剖分，如圖1所示（X軸以水平向右為正，Y軸以豎直向上為正），有限元計(jì)算模型共含有1 062個(gè)節(jié)點(diǎn)和986個(gè)單元。

圖1 計(jì)算模型／m

在重力荷載作用下，隨著土體強(qiáng)度的折減，邊坡內(nèi)特征點(diǎn)的位移變化和塑性應(yīng)變發(fā)展的最終情況如圖2、3所示。

圖2 特征點(diǎn)位移變化曲線圖

從圖2中可以看出，M－C準(zhǔn)則與DP1準(zhǔn)則兩者計(jì)算結(jié)果基本上是一致的；但是DP1準(zhǔn)則計(jì)算時(shí)，數(shù)值的收斂性較好，能夠計(jì)算到邊坡失穩(wěn)的工況，而且還能夠繼續(xù)計(jì)算下去，此時(shí)數(shù)值收斂判據(jù)失真，采用位移判據(jù)進(jìn)行分析，在F＝1.23時(shí)位移曲線出現(xiàn)了明顯的轉(zhuǎn)折點(diǎn)，可以認(rèn)為其最小安全系數(shù)為1.23；M－C準(zhǔn)則的數(shù)值收斂性不好，得到的位移曲線沒(méi)有出現(xiàn)明顯的轉(zhuǎn)折點(diǎn)，位移已經(jīng)基本上不再變化，很大程度上是迫于數(shù)值的不收斂，此時(shí)得到的最小安全系數(shù)也為1.23，與DP1的計(jì)算結(jié)果是相同的。

圖3中PEEQ表示的是邊坡在整個(gè)變形過(guò)程中塑性應(yīng)變的累積結(jié)果。根據(jù)本文所建議的塑性應(yīng)變與塑性區(qū)發(fā)展的失穩(wěn)判據(jù)，隨著土體強(qiáng)度參數(shù)的不斷折減，土體塑性區(qū)發(fā)展區(qū)域逐漸增大。對(duì)于M－C準(zhǔn)則來(lái)說(shuō)，當(dāng)F＝1.23時(shí)，邊坡塑性區(qū)基本貫通，坡腳塑性應(yīng)變已達(dá)3.67%。對(duì)于DP1準(zhǔn)則來(lái)說(shuō)，當(dāng)F＝1.23時(shí)，邊坡塑性區(qū)基本貫通，坡腳塑性應(yīng)變已達(dá)2.00%。當(dāng)F＝1.24時(shí)，塑性區(qū)基本貫通至坡頂，邊坡已經(jīng)失穩(wěn)，認(rèn)為其最小安全系數(shù)均為1.23。從上述計(jì)算結(jié)果來(lái)看，2個(gè)材料模型在采用不同的失穩(wěn)判據(jù)得到的計(jì)算結(jié)果基本一致。

圖3 邊坡塑性應(yīng)變最終分布圖

4 結(jié)語(yǔ)

（1）在平面應(yīng)變問(wèn)題中，采用DP1準(zhǔn)則進(jìn)行模擬計(jì)算，能得到很好的結(jié)果，克服了M－C準(zhǔn)則在數(shù)值分析時(shí)存在的不足，減少了數(shù)值計(jì)算的困難，在一定程度上，避免了收斂性的問(wèn)題，同時(shí)充分吸取了M－C準(zhǔn)則在描述巖土強(qiáng)度特性時(shí)的有效性，完全可以在巖土工程中獲得廣泛的應(yīng)用。

（2）在ABAQUS中采用場(chǎng)變量進(jìn)行強(qiáng)度參數(shù)的折減，能起到很好的效果，不需要每算一次改一次土體強(qiáng)度參數(shù)，大大減少了計(jì)算工作量，而且所得計(jì)算結(jié)果合理，表明了該方法的可行性。

（3）當(dāng)采用位移判據(jù)時(shí)應(yīng)選取多個(gè)位移特征點(diǎn)，以避免在位移變化曲線上找不到明顯的轉(zhuǎn)折點(diǎn)。

（4）當(dāng)采用強(qiáng)度折減彈塑性有限元法進(jìn)行計(jì)算分析時(shí)，數(shù)值收斂判據(jù)在應(yīng)用上存在局限性，不具有廣泛的適用性，盡管位移判據(jù)和塑性區(qū)判據(jù)在計(jì)算過(guò)程中也受到很多因素的影響，但是這些因素不足以改變邊坡破壞的本質(zhì)趨勢(shì)，可以聯(lián)合采用位移判據(jù)和塑性區(qū)判據(jù)來(lái)對(duì)邊坡的失穩(wěn)進(jìn)行分析，以便得到更合理適用的結(jié)果。

［1］鄧楚鍵，何國(guó)杰，鄭穎人.基于 M－C準(zhǔn)則的D－P系列準(zhǔn)則在巖土工程中的應(yīng)用研究［J］.巖土工程學(xué)報(bào)，2006，28（6）：735－739.

［2］張艷山，潘玉珍.基于ABAQUS對(duì)系列準(zhǔn)則的討論［J］.水電能源科學(xué)，2010，28（11）：70－72.

［3］蘇凱，伍鶴皋，李沖.Mohr—Coulomb等面積圓屈服準(zhǔn)則在圍巖穩(wěn)定分析中的應(yīng)［J］.中國(guó)農(nóng)村水利水電，2008，11：81－86.

［4］姚燕稚，陳國(guó)興.基于場(chǎng)變量控制強(qiáng)度參數(shù)折減的有限元邊坡穩(wěn)定分析［J］.湖南大學(xué)學(xué)報(bào)，2008，35（11）：124－129.

［5］Zienkiewicz O C，Humpheson C，Lewis R W.A ssoci－ated and non－associated visco－plasticity and plasticity in soil mechanics［J］.Geostechnique，1975，25（4）：671－689.

［6］Bishop A W.The use of the slip circle in the stability analysis of slopes［J］.Geostechnique，1955（5）：7－17.

［7］張陽(yáng)陽(yáng)，李宗坤，李艷.基于ABAQUS的強(qiáng)度折減法邊坡失穩(wěn)判據(jù)研究［J］.浙江水利水電專(zhuān)科學(xué)校學(xué)報(bào)，2009，21（1）：13－15.

［8］趙尚毅，鄭穎人，張玉芳.極限分析有限元法講座—Ⅱ有限元強(qiáng)度折減法中邊坡失穩(wěn)的判據(jù)探討［J］.巖土力學(xué)，2005，26（2）：332－336.

［9］劉金龍，欒茂田，趙少飛，等.關(guān)于強(qiáng)度折減有限元法中邊坡失 穩(wěn) 判 據(jù)的 討 論 ［J］.巖 土 力 學(xué)，2005，26（8）：1345－1348.