姚 云

(懷化市農(nóng)村公路建設(shè)辦公室，湖南 懷化 418000)

基于物質(zhì)點(diǎn)法考慮剪脹特性的土質(zhì)邊坡穩(wěn)定性分析

姚 云

(懷化市農(nóng)村公路建設(shè)辦公室，湖南 懷化 418000)

自然界中的土體其體積變形往往具有剪脹特性，為了探究土體的剪脹特性對土質(zhì)邊坡穩(wěn)定性的影響，采用物質(zhì)點(diǎn)法模擬了具有不同剪脹角大小的土質(zhì)邊坡的滑動破壞全過程。計算結(jié)果表明，隨著土體剪脹角的逐漸增大，土質(zhì)邊坡的滑動破壞區(qū)域逐漸減小，并且其最終的滑動距離也相應(yīng)地減小，表明增大土體的剪脹角能夠有效地提高土質(zhì)邊坡的穩(wěn)定性，在實際工程應(yīng)用中具有一定的指導(dǎo)意義。最后初步探究了采用最終滑動距離作為土質(zhì)邊坡穩(wěn)定性指標(biāo)的補(bǔ)充的可行性。

土質(zhì)邊坡；剪脹特性；物質(zhì)點(diǎn)法；大變形破壞

0 引言

在地震或強(qiáng)降雨等自然因素的作用下邊坡可能發(fā)生失穩(wěn)破壞，產(chǎn)生大面積滑塌，造成嚴(yán)重的經(jīng)濟(jì)損失。針對邊坡的穩(wěn)定性分析和破壞機(jī)理探究引起了廣大學(xué)者極大的關(guān)注，人們對此展開了大量的研究，并且取得了很多有益的成果。

近些年，隨著計算機(jī)運(yùn)算水平的提升，數(shù)值模擬方法，尤其是有限單元法在巖土工程領(lǐng)域具有廣泛的應(yīng)用。然而，由于邊坡在滑動破壞的過程中會產(chǎn)生很大的變形，采用傳統(tǒng)的有限單元法處理該問題時，在計算的過程中會出現(xiàn)網(wǎng)格畸變的問題，導(dǎo)致計算結(jié)果不收斂。物質(zhì)點(diǎn)法作為一種新型的無網(wǎng)格粒子法，能夠有效地避免計算過程中可能出現(xiàn)的網(wǎng)格畸變的問題，能夠模擬連續(xù)介質(zhì)任意大變形問題。物質(zhì)點(diǎn)法由Sulsky等[1]首先提出，采用歐拉和拉格朗日雙重描述，即采用質(zhì)點(diǎn)離散材料區(qū)域，用背景網(wǎng)格計算空間導(dǎo)數(shù)和求解動量方程。根據(jù)現(xiàn)有的資料，物質(zhì)點(diǎn)法被大量應(yīng)用于邊坡穩(wěn)定性分析。Andersen等[2]采用物質(zhì)點(diǎn)法分析了邊坡失效過程。廉燕平等[3]應(yīng)用自適應(yīng)物質(zhì)點(diǎn)有限元法分析了邊坡的失效問題，計算結(jié)果表明物質(zhì)點(diǎn)能夠較好地模擬邊坡失穩(wěn)破壞的全過程。王雙等[4]提出物質(zhì)點(diǎn)強(qiáng)度折減法計算了不同安全系數(shù)下邊坡內(nèi)部等效塑性變形的分布特點(diǎn)。

然而，上述研究在采用物質(zhì)點(diǎn)法探究邊坡滑動破壞機(jī)理時并沒有考慮土體剪脹特性對其安全系數(shù)的影響。Sloan等[5]指出在計算邊坡安全系數(shù)時有必要考慮土體的剪脹特性，其對于邊坡內(nèi)部滑動面的形成過程具有顯著的影響。張培文等[6]采用強(qiáng)度折減法計算邊坡安全系數(shù)時考慮了土體的剪脹特性，計算結(jié)果表明，隨著土體剪脹角的逐漸增大邊坡的安全系數(shù)逐漸增大。

本文考慮土體的剪脹特性，采用物質(zhì)點(diǎn)法分別模擬了不同剪脹角下的邊坡滑動破壞過程。計算結(jié)果土體的剪脹特性對于邊坡內(nèi)部的等效應(yīng)變的分布具有顯著的影響，并且從能量的角度對該現(xiàn)象進(jìn)行了一定程度的解釋。

1 物質(zhì)點(diǎn)法

物質(zhì)點(diǎn)法最主要的特點(diǎn)在于其采用歐拉和拉格朗日雙重描述。物質(zhì)點(diǎn)法采用質(zhì)點(diǎn)積分，將積分轉(zhuǎn)化為被積函數(shù)在各質(zhì)點(diǎn)處的值與該質(zhì)點(diǎn)所代表的體積之積的和。物質(zhì)點(diǎn)法的背景網(wǎng)格旨在為每個計算步內(nèi)和物體固連，在每個計算步結(jié)束后丟棄已經(jīng)變形的背景網(wǎng)格。由于質(zhì)點(diǎn)已經(jīng)攜帶了物質(zhì)體的所有物質(zhì)信息，在下一個時刻，可以通過將指點(diǎn)的信息映射到新的背景網(wǎng)格上來求得網(wǎng)格的信息，因此在物質(zhì)點(diǎn)法中不需要記錄每個時刻網(wǎng)格節(jié)點(diǎn)的任何信息，如圖1所示。

在顯示物質(zhì)點(diǎn)法求解過程，應(yīng)力更新主要有兩種方法，一是在每個計算時間步開始時進(jìn)行，稱為USF方法，二是在每個計算步結(jié)束后進(jìn)行，稱為USL。Naim等指出USF格式具有較好的能量守恒性。 物質(zhì)點(diǎn)法的求解過程主要包含以下幾方面：

1)將每個質(zhì)點(diǎn)的質(zhì)量和動量映射到網(wǎng)格上，計算背景網(wǎng)格節(jié)點(diǎn)的質(zhì)量和動量；

2)對節(jié)點(diǎn)動量施加本質(zhì)邊界條件；

3)由網(wǎng)格節(jié)點(diǎn)的動量計算節(jié)點(diǎn)速度，由此計算每個質(zhì)點(diǎn)的應(yīng)變增量和旋量增量，并對各支點(diǎn)的密度和應(yīng)力進(jìn)行更新；

4)計算背景網(wǎng)格節(jié)點(diǎn)力；

5)在背景網(wǎng)格節(jié)點(diǎn)上對動量方程進(jìn)行積分；

6)將背景網(wǎng)格節(jié)點(diǎn)的速度變化量和位置變化量映射回相應(yīng)的節(jié)點(diǎn)；

7)生成新的規(guī)則的背景網(wǎng)格，開始下一步的計算。

在整個計算過程中，計算信息在物質(zhì)點(diǎn)和背網(wǎng)格間相互映射，能夠有效地避免網(wǎng)格畸變和歐拉對流項處理。

圖1 物質(zhì)點(diǎn)法計算過程示意圖

2 邊坡穩(wěn)定性分析

2.1 考慮剪脹特性的土體力學(xué)模型

邊坡的穩(wěn)定性分析一直是巖土領(lǐng)域的一個熱點(diǎn)的研究課題，通常采用安全系數(shù)來評價其穩(wěn)定性。在具體求解邊坡安全系數(shù)的過程中，很多學(xué)者忽略了土體剪脹特性對土體力學(xué)特性的影響。為了考慮土體的剪脹特性，本文采用非關(guān)聯(lián)的莫爾-庫侖彈塑性模型描述土體的力學(xué)特性。

莫爾-庫倫彈塑性模型塑性勢函數(shù)和屈服方程有如下形式:

(1)

(2)

其中，φ為內(nèi)摩擦角;c為黏聚力;ψ為膨脹角。

由于本文中采用張雄等[7]研發(fā)的MPM3 D程序，其內(nèi)部只包含了Drucker-Prager彈塑性模型，在實際計算的過程需要將莫爾-庫倫彈塑性模型的參數(shù)進(jìn)行轉(zhuǎn)化。

Drucker-Prager彈塑性模型的屈服函數(shù)為：

(3)

其中，J2為偏應(yīng)力張量的第2不變量；σm為球應(yīng)力；α為摩擦系數(shù)；k為純剪切狀態(tài)時的屈服應(yīng)力。

材料參數(shù)α、k和c、φ之間的關(guān)系如下：

(4)

(5)

在張雄等所給出的Drucker-Prager彈塑性模型中，剪切式函數(shù)ψs具有如下表達(dá)式：

ψs=βσm

(6)

其中，β為考慮剪脹特性的力學(xué)參數(shù)，其與土體的剪脹角ψ之間具有如下關(guān)系。

(7)

2.2 計算模型

在實際工程應(yīng)用中，由于人工開挖等因素的影響，道路兩旁的邊坡中的土體往往具有超固結(jié)性，在自身重力或外部荷載的作用下其體積會發(fā)生剪脹，進(jìn)而對邊坡的穩(wěn)定性產(chǎn)生影響。為了探究不同剪脹角所對應(yīng)的土體內(nèi)部等效塑性變形的分布，本文假定土體的剪脹角為0°、5°、10°、15°、30°和45°這6種情況針對均質(zhì)土質(zhì)邊坡進(jìn)行數(shù)值模擬，計算模型如圖2所示。

圖2 土質(zhì)邊坡計算模型(單位：m)

計算模型中的物質(zhì)點(diǎn)數(shù)目為11 968個，計算模型的參數(shù)如表1所示。

表1 模型計算參數(shù)

每個計算步的時間為0.1 s，總的計算時間為15 s。

2.3 計算結(jié)果及分析

不考慮土體剪脹角時，土質(zhì)邊坡在自身重力的作用下發(fā)生滑動破壞，如圖3所示，可以看出，土質(zhì)邊坡坡腳處首先發(fā)生屈服，產(chǎn)生較大的等效塑性變形，并且沿著滑動破壞面向上擴(kuò)展，當(dāng)滑動破壞面完全貫通后，土質(zhì)邊坡發(fā)生大變形破壞。由于物質(zhì)點(diǎn)法屬于無網(wǎng)格法，在計算模型產(chǎn)生較大的體積變形時，不需要考慮網(wǎng)格畸變所帶來的計算結(jié)果發(fā)生不收斂的問題，因此能夠采用物質(zhì)點(diǎn)法準(zhǔn)確模擬土質(zhì)邊坡的滑動破壞全過程，為后面探究土體剪脹特性對邊坡穩(wěn)定性的影響打下了良好的基礎(chǔ)。

圖3 不考慮剪脹特性時土質(zhì)邊坡滑動破壞全過程(等效塑性變形)

不同剪脹角下土體邊坡最終的滑動形態(tài)如圖4所示，可以看出，隨著剪脹角的逐漸增大，該邊坡的最終滑動距離逐漸減小，另外，該邊坡滑動區(qū)的范圍也逐漸的減小。該計算結(jié)果表明隨著土體剪脹角的逐漸增大，土質(zhì)邊坡的穩(wěn)定性會逐漸增大，該結(jié)論與增大土體剪脹角會導(dǎo)致土體的安全系數(shù)增大的結(jié)論保持一致。

本文計算的土質(zhì)邊坡最終滑動距離隨土體剪脹角的變化曲線如圖5所示，可以看出，隨著剪脹角的增大，土質(zhì)邊坡的最終滑動距離逐漸減小，并且具有較好的線性規(guī)律。

圖4 土體剪脹特性對土質(zhì)邊坡最終的破壞形態(tài)的影響(等效塑性變形)

圖5 土質(zhì)邊坡最終的滑動距離隨土體剪脹角的變化規(guī)律

值得注意的是，在實際工程應(yīng)用通常采用安全系數(shù)作為邊坡穩(wěn)定性大小的評判指標(biāo)。然而安全系數(shù)主要考慮邊坡不發(fā)生滑動破壞的情況。事實上，對于一些安全性要求并不是很高的邊坡，在治理的過程中只需要保證邊坡最終的滑動距離遠(yuǎn)離附近的建筑物，因此可以采用物質(zhì)點(diǎn)法等數(shù)值計算手段計算得到邊坡的最終滑動距離，根據(jù)其與附近建筑物的距離進(jìn)行相應(yīng)的治理，進(jìn)而節(jié)約治理成本。

不同剪脹角下土質(zhì)邊坡在滑動破壞的過程中總動能和總內(nèi)能隨著計算時間的變化規(guī)律如圖6、圖7所示。根據(jù)圖6可以看出，隨著計算時間的增大土質(zhì)邊坡內(nèi)部的總動能先逐漸增大，達(dá)到峰值后逐漸減小。并且隨著土體剪脹角的增大土質(zhì)邊坡總動能峰值逐漸增大。根據(jù)圖7可以看出隨著隨著計算時間的增大土質(zhì)邊坡的總內(nèi)能逐漸增大，最終趨于穩(wěn)定，并且隨著土體剪脹角的增大，土質(zhì)邊坡的總內(nèi)能逐漸增大。

圖6 不同剪脹角下土質(zhì)邊坡滑動破壞過程中總動能隨時間的變化規(guī)律

圖7 不同剪脹角下土質(zhì)邊坡滑動破壞過程中總內(nèi)能隨時間變化規(guī)律

3 結(jié)論

本文采用物質(zhì)點(diǎn)法對土質(zhì)邊坡開展了一系列的數(shù)值試驗，重點(diǎn)探究了土體剪脹特性對邊坡穩(wěn)定性的影響，主要結(jié)論有：

1)物質(zhì)點(diǎn)法能夠較好地模擬邊坡的滑動破壞全過程，能夠用來探究不同影響因素對邊坡穩(wěn)定性的影響。

2)土體剪脹角對邊坡的滑動機(jī)理具有顯著的影響，表現(xiàn)為隨著剪脹角的增大，土質(zhì)邊坡的最終滑動距離逐漸減小，土質(zhì)邊坡的穩(wěn)定性逐漸增大。

3)邊坡最終的滑動距離可以作為邊坡穩(wěn)定性分析的新的參考指標(biāo)，考慮邊坡的失穩(wěn)破壞對臨近建筑物的影響，從而縮減邊坡治理成本。

[1] Sulsky D, Schreyer H L. Axisymmetric form of the material point method with applications to upsetting and Taylor impact problems[J]. Computer Methods in Applied Mechanics & Engineering, 1996, 139(1):409-429.

[2] Andersen S, Andersen L. Modelling of landslides withthe material-point method[J]. Computational Geosciences, 2010, 14(1):137-147.

[3] 廉艷平, 張帆, 劉巖,等. 物質(zhì)點(diǎn)法的理論和應(yīng)用[J]. 力學(xué)進(jìn)展, 2013, 43(2):237-264.

[4] 王雙, 李小春, 石露,等. 物質(zhì)點(diǎn)強(qiáng)度折減法及其在邊坡中的應(yīng)用[J]. 巖土力學(xué), 2016, 37(9):2672-2678.

[5] Tschuchnigg F, Schweiger H F, Sloan S W. Slope stability analysis by means of finite element limit analysis and finite element strength reduction techniques. Part I: Numerical studies considering non-associated plasticity[J]. Computers & Geotechnics, 2015, 70:169-177.

[6] 張培文, 陳祖煜. 剪脹角對求解邊坡穩(wěn)定的安全系數(shù)的影響[J]. 巖土力學(xué), 2004, 25(11):1757-1760.

[7] 張雄, 宋康祖, 陸明萬. 無網(wǎng)格法研究進(jìn)展及其應(yīng)用[J]. 計算力學(xué)學(xué)報, 2003, 20(6):730-742.

1008-844X(2017)02-0082-04

U 416.1+4

A