艾國棟，李云青，資麗君，梁昌奇

（1.湖南省礦產(chǎn)資源調(diào)查所，湖南 郴州 423000；2.中南大學有色金屬成礦預測與地質(zhì)環(huán)境監(jiān)測教育部重點實驗室，湖南 長沙 410083；3.湖南省有色資源與地質(zhì)災害探查湖南省重點實驗室，湖南 長沙 410083；4.中南大學地球科學與信息物理學院，湖南 長沙 410083）

0 引言

邊坡失穩(wěn)產(chǎn)生滑坡是一個典型的巖土工程大變形問題［1］，其造成的破壞后果通常取決于邊坡失穩(wěn)后滑坡的運動速度、滑動距離和堆積形態(tài)等［2-5］?；峦ǔ＞哂型话l(fā)性、滑速快和動能大等特點［6］，嚴重威脅著其影響范圍內(nèi)的建筑物和人民生命財產(chǎn)安全，如發(fā)生在我國的灑勒山滑坡［7］和深圳渣土場滑坡［8-9］，造成了嚴重的經(jīng)濟損失和不良社會影響。特別是對于體積超過105萬方的滑坡，在重力勢能的作用下，其滑動距離通?？蛇_其坡高的10 倍以上［10-11］，超過數(shù)公里，因而影響范圍巨大。因此，為準確定量評估邊坡風險，科學合理地分析邊坡失穩(wěn)后的狀態(tài)，即預測滑體的運動范圍和堆積過程，對實際滑坡防災減災，具有十分重要的理論和現(xiàn)實意義［12-13］。

傳統(tǒng)的邊坡穩(wěn)定性分析方法如極限平衡法和有限元法只能求解邊坡破壞前的狀態(tài)，不能有效處理滑坡大變形問題。物質(zhì)點法在每個計算時間步都采用規(guī)則的背景網(wǎng)格，計算時網(wǎng)格不會扭曲變形，是求解與材料變形有關(guān)問題的有效方法，近年來在巖土工程領域中得到了廣泛的應用。例如，王雙等［14］采用物質(zhì)點強度折減法研究了邊坡失穩(wěn)后對滑動路徑上建筑物的影響程度，并通過算例驗證了MPM 預測失穩(wěn)后滑坡堆積形態(tài)和滑動距離的能力；孫玉進等［15］利用物質(zhì)點法計算大變形的優(yōu)勢模擬了某高邊坡從開始變形滑動到穩(wěn)定的全過程，證明該方法可以較好地再現(xiàn)邊坡滑動全過程；Zhao 等［16］將物質(zhì)點法應用于流體力學模型，模擬了不同縱橫比條件下的潰壩流動問題，驗證了物質(zhì)點法的準確性和穩(wěn)定性；Troncone 等［17］用物質(zhì)點法分析了1986 年發(fā)生在意大利南部的森尼斯滑坡的破壞過程；Yerro 等［18］采用物質(zhì)點法模擬了華盛頓奧索鎮(zhèn)滑坡的滑動過程；Xu 等［19］使用物質(zhì)點法、DP 本構(gòu)模型和流變關(guān)系，結(jié)合ArcGIS 軟件對中國紅石巖滑坡進行了數(shù)值模擬；Conte 等［20］用物質(zhì)點法對1986 年發(fā)生在意大利南部Senise 滑坡的徑流過程進行了二維分析；Coelho等［21］和Dong 等［22］使用物質(zhì)點法研究了剛性基礎上的滑坡流動問題，并用物質(zhì)點法和深度平均法對海底滑坡的實際滑動過程進行了反演分析，再現(xiàn)了從陡坡到中等坡度的滑坡滑動過程；Li 等［23］使用物質(zhì)點法對王家?guī)r滑坡破壞后的滑動過程進行了數(shù)值模擬；Huang 等［24］使采用物質(zhì)點法和DP 本構(gòu)模型對粘性土坡和無粘性土坡的大變形破壞過程進行了數(shù)值分析。

上述研究證明了物質(zhì)點法在模擬邊坡滑坡大變形方面具有較強可行性，但鮮有研究從工程實際出發(fā)，研究不同工況下的邊坡大變形機制。為此，本文以湖南省郴州市某滑坡為例，基于物質(zhì)點法分析該邊坡開挖前后在天然工況和暴雨工況下的大變形特征，并與傳統(tǒng)邊坡穩(wěn)定性分析方法作對比。

1 物質(zhì)點法

物質(zhì)點法是1994 年Sulsky 等［25］提出的一種用于計算大變形問題的數(shù)值計算方法。其基本原理是將連續(xù)體離散成一系列質(zhì)點，每個質(zhì)點代表一塊區(qū)域并攜帶了該區(qū)域的所有物質(zhì)信息，如：質(zhì)量、密度、速度、加速度、應力、應變等，而網(wǎng)格在空間中固定，不攜帶任何參數(shù)，僅用于動量方程的求解和空間導數(shù)的計算。該方法具有拉格朗日質(zhì)點和歐拉網(wǎng)格的雙重優(yōu)勢，通過背景網(wǎng)格與物質(zhì)點間的交替映射完成數(shù)值計算，在每個時間步都采用未變形的背景網(wǎng)格，避免了因網(wǎng)格畸變和單元扭曲帶來的數(shù)值求解困難，可以很好的解決與邊坡大變形有關(guān)的問題。

物質(zhì)點法在每個時間步內(nèi)的計算可分為3 個階段，即初始化階段、拉格朗日計算階段和映射階段［26］（如圖1 所示）。在初始化階段，物質(zhì)點攜帶的所有物質(zhì)信息通過線性形函數(shù)被映射到背景網(wǎng)格上，為拉格朗日計算階段做準備；在拉格朗日計算階段，在背景網(wǎng)格上計算動量方程和導數(shù)，并通過計算獲得新的節(jié)點速度、加速度和位移等信息；最后在映射階段，更新所有物質(zhì)點狀態(tài)信息，如物質(zhì)店的位移、速度等物理量在該時間步計算結(jié)束后，背景網(wǎng)格恢復至初始化階段時的形態(tài)，并在下一個時間步重復上述三個階段的計算。

圖1 物質(zhì)點法計算流程Fig.1 Calculation process of material point method

2 滑坡基本情況介紹

該滑坡位于郴州市第七完全小學北側(cè)，地貌為剝蝕殘丘地貌，總體地勢北高南低，地勢起伏變化較大。場地內(nèi)及其周邊未見有明顯活動的斷裂構(gòu)造通過，且新構(gòu)造運動不明顯，處于相對穩(wěn)定狀態(tài)，無地表水系，周邊未發(fā)現(xiàn)崩塌、泥石流、地下采空區(qū)等不良地質(zhì)作用。地層主要為第四系坡殘積成因粉質(zhì)粘土，二疊系下統(tǒng)當沖組全風化泥頁巖以及全-強風化炭質(zhì)頁巖，節(jié)理裂隙極發(fā)育，具體計算參數(shù)見表1。

表1 邊坡各土層土體參數(shù)Table 1 Soil parameters of each soil layer of slope

該滑坡目前處于蠕滑變形階段，局部已經(jīng)處于變形加速階段。地表水下滲后富集在粉質(zhì)粘土內(nèi)部滑動面，水量不斷富集、運移，該粉質(zhì)粘土的抗剪強度不斷降低，形成軟弱滑帶，使粉質(zhì)粘土層內(nèi)部產(chǎn)生滑動。該滑坡在天然狀態(tài)下處于基本穩(wěn)定狀態(tài)，但在長期降雨或暴雨時，地表水不斷入滲滑體，造成滑體自重增加，滑面抗剪強度降低，可能發(fā)生更大的位移，嚴重威脅滑坡體下方居民生命、財產(chǎn)安全。

因此本文用有限元強度折減法［27-28］計算了該邊坡開挖前后在天然工況和暴雨工況下的穩(wěn)定性，并用物質(zhì)點法［29］計算了開挖前后在天然工況和暴雨工況下的滑動距離，以分析邊坡失穩(wěn)后對滑動路徑上建筑物的威脅程度。

3 穩(wěn)定性分析

利用有限元分析軟件MIDAS GTS NX，分析該邊坡開挖前后在天然工況和暴雨工況下的穩(wěn)定性。在構(gòu)建有限元模型時，根據(jù)現(xiàn)場實際情況考慮了不同土層材料。同時，邊坡模型左右邊界在水平方向進行了固定，底邊界在水平方向和垂直方向進行了固定。

其中，局部計算結(jié)果通過在有限元計算過程中對邊坡模型進行全局搜索滑面得到，而整體計算結(jié)果通過在搜索滑面時排除前述安全系數(shù)小的局部區(qū)域得到。開挖前邊坡在天然工況下的計算得到的等效塑性應變云圖見圖2。由圖2 可知，該邊坡首先在坡腳處形成較小的等效塑性應變貫通區(qū)，安全系數(shù)為1.08，隨后，在邊坡上部產(chǎn)生較大的等效塑性應變貫通區(qū)，安全系數(shù)為1.61。表明：邊坡整體處于穩(wěn)定狀態(tài)，但在坡腳處可能會發(fā)生局部破壞。

圖2 天然工況下開挖前邊坡等效塑性應變云圖Fig.2 Equivalent plastic strain nephogram of slope before excavation under natural condition

開挖前邊坡在暴雨工況下的塑性應變云圖見圖3。與天然工況下的等效塑性應變云圖類似，邊坡首先在坡腳處形成等效塑性應變貫通區(qū)，安全系數(shù)為0.80，隨后，在邊坡上部形成較大的等效塑性應變貫通區(qū)，安全系數(shù)為1.22。表明：在暴雨工況下，邊坡處于不穩(wěn)定狀態(tài)。

圖3 暴雨工況下開挖前邊坡等效塑性應變云圖Fig.3 Equivalent plastic strain nephogram of slope before excavation under rainstorm condition

開挖后邊坡在天然工況和暴雨工況下的等效塑性應變云圖見圖4，安全系數(shù)分別為0.60 和0.50，處于不穩(wěn)定狀態(tài)。這是因為開挖使得邊坡前緣形成臨空面，邊坡抗滑力減小，下滑力不變，進而導致邊坡失穩(wěn)［30］。

圖4 開挖暴雨工況邊坡安全系數(shù)Fig.4 Slope safety factor after excavation under rainstorm condition

4 大變形分析

由前面的分析可知，開挖前邊坡在暴雨工況下處于基本穩(wěn)定狀態(tài)，開挖后處于不穩(wěn)定狀態(tài)，因此本節(jié)使用物質(zhì)點法［31］計算開挖前后邊坡在天然工況和暴雨工況下的滑動距離，以分析其對滑動路徑上建筑物的威脅程度。

4.1 計算參數(shù)

本研究中，邊坡長260 m，高87 m，物質(zhì)點大小為2.0 m，共離散成6668 個物質(zhì)點，如圖5 所示。為了簡化問題，本研究僅考慮平面應變情況。為了更好的還原邊坡現(xiàn)場的真實狀態(tài)，邊坡底部選用固定邊界條件以模擬基巖，兩側(cè)和前后使用對稱邊坡邊界條件以模擬平面應變狀態(tài)下土體的各種狀態(tài)。本次計算選用8 節(jié)點六面體單元，單元尺寸4 m，采用非關(guān)聯(lián)流動準則和服從Drucker-Prager（DP）屈服準則的彈塑性本構(gòu)模型。計算時間步長為0.465×10-3s，總計算時間20 s。

圖5 邊坡離散示意Fig.5 Schematic diagram of slope dispersion

4.2 大變形模擬

開挖前邊坡在天然和暴雨工況下的滑動距離計算云圖如圖6 所示，開挖后邊坡在天然和暴雨工況下的滑動距離計算云圖如圖7 所示。由圖6 可知，開挖前邊坡在天然工況和暴雨工況下的滑動距離分別為1.37 和3.75 m。由圖7 可知，開挖后邊坡在天然工況和暴雨工況下的滑動距離分別為20.11 和24.12 m。由此可知，開挖會使邊坡形成不穩(wěn)定臨空面，使邊坡抗滑力減小，穩(wěn)定性降低；在雨水的長期作用或暴雨情況下，地表水向邊坡土體滲透，滑體自重增加，滑面抗剪強度降低，其穩(wěn)定性進一步降低，失穩(wěn)時滑動的距離也進一步增大［32-34］；并且，邊坡對于在滑動距離范圍內(nèi)的建筑物威脅程度較大，計算結(jié)果可為提高邊坡安全性防護提供理論參考。

圖6 開挖前邊坡滑動距離云圖Fig.6 Contour plot of slope sliding distance before excavation

圖7 開挖后邊坡滑動距離云圖Fig.7 Contour plot of slope sliding distance after excavation

由穩(wěn)定性分析可知，在天然狀態(tài)下，開挖前邊坡的整體穩(wěn)定性系數(shù)為1.61，坡腳的局部穩(wěn)定性系數(shù)為1.08，表明：邊坡不會滑動，坡腳處于基本穩(wěn)定狀態(tài)。由大變形分析結(jié)果可知，天然狀態(tài)下邊坡滑動距離為1.37 m，這是因為在本次計算中視土體在空間上為各項同性，且土體是經(jīng)過室內(nèi)試驗測試獲得的，在測試前土體性質(zhì)已經(jīng)擾動，原有性質(zhì)發(fā)生了改變，所以安全系數(shù)計算結(jié)果存在一定的誤差。開挖前暴雨工況和開挖后邊坡安全系數(shù)均＜1，所以邊坡會發(fā)生滑動。需要說明的是，由于本文進行的有限元模擬假定土體發(fā)生小變形，因此該方法不用于計算滑坡的大變形過程，而用于進行邊坡的初始破壞階段安全系數(shù)的計算。在物質(zhì)點法計算過程中，物質(zhì)點法可以計算土體的大變形過程，從而用于本文邊坡滑動過程的模擬。

5 結(jié)論

（1）開挖前邊坡在天然工況下的整體安全系數(shù)為1.61，處于穩(wěn)定狀態(tài)，局部安全系數(shù)為1.08，處于基本穩(wěn)定狀態(tài)；在暴雨工況下的整體安全系數(shù)1.22，處于基本穩(wěn)定狀態(tài)，局部安全系數(shù)為0.80，處于失穩(wěn)狀態(tài)。表明：在天然工況下，該邊坡處于基本穩(wěn)定狀態(tài)，在暴雨工況下，該邊坡處于失穩(wěn)狀態(tài)。

（2）開挖后邊坡在天然工況和暴雨工況下的安全系數(shù)分別為0.60 和0.50，處于失穩(wěn)狀態(tài)，且失穩(wěn)后的滑動距離分別為20.11 和24.12 m。表明：開挖和降雨對邊坡滑動距離有顯著影響，使失穩(wěn)后的邊坡滑動距離增大。

（3）物質(zhì)點法在計算邊坡在不同工況下的滑動距離的同時，還可以評估邊坡失穩(wěn)對滑動路徑上建筑物的威脅程度，對提高邊坡安全性防護具有重要意義。