劉 永，張國(guó)和*，楊建華，董 鵬，毛吉成，甘小迎

(中國(guó)電建集團(tuán)貴州電力設(shè)計(jì)研究院有限公司，貴州 貴陽(yáng))

引言

隨著近代科學(xué)技術(shù)的快速發(fā)展，計(jì)算機(jī)科學(xué)與技術(shù)在地質(zhì)工程領(lǐng)域得到了極大的推廣發(fā)展與應(yīng)用，基于計(jì)算機(jī)有限元數(shù)值模擬軟件的計(jì)算分析法在滑坡穩(wěn)定性分析中不斷發(fā)展運(yùn)用，并成為一種新的穩(wěn)定性分析方法[1]，同時(shí)運(yùn)用數(shù)值模擬方法對(duì)滑坡的穩(wěn)定性及治理方案進(jìn)行分析研究將成為未來(lái)研究分析的主要發(fā)展趨勢(shì)[2]。Midas GTS / NX 作為一種巖土有限元分析軟件，其具有強(qiáng)大的分析計(jì)算能力，廣泛適用于巖土工程、樁基工程、地下工程、水利工程等領(lǐng)域[3]，有限元數(shù)值模擬方法已經(jīng)成為滑坡坡穩(wěn)定性分析中的重要方法，這類方法不但能考慮邊坡巖土體本身的變形對(duì)邊坡穩(wěn)定性的影響，而且能給出邊坡巖土體中應(yīng)力應(yīng)變分布，分析邊坡破壞的發(fā)生和發(fā)展過(guò)程等，具有傳統(tǒng)的極限平衡法所不可比擬的優(yōu)勢(shì)[4]。

文章將借助Midas GTS / NX 巖土有限元分析軟件，對(duì)滑坡的穩(wěn)定性及治理方案可行性進(jìn)行研究，對(duì)滑坡在不同工況下的穩(wěn)定性進(jìn)行模擬計(jì)算[5]，從滑坡穩(wěn)定性系數(shù)，位移及最大剪應(yīng)變等方面對(duì)滑坡穩(wěn)定性進(jìn)行分析，并借助Midas GTS / NX 對(duì)該滑坡治理方案進(jìn)行模擬分析，對(duì)比采取治理方案前模擬結(jié)果，對(duì)治理方案的可行性進(jìn)行進(jìn)一步分析研究。

1 滑坡概況

該滑坡位于邊坡中上部，兩側(cè)沖溝發(fā)育，下部受人類工程活動(dòng)開挖嚴(yán)重，滑坡上部坡度在12°～22°之間，中部和下部坡度在43°～60°之間，滑坡平面展布上呈現(xiàn)不規(guī)則半橢圓狀，縱向長(zhǎng)約120 m，中部寬約220 m，前緣寬約260 m，主滑方向約為240°?；抡w平均厚度達(dá)24.9 m；滑坡平面面積約2.66×104m2，滑體體積約為69.46×104m3，屬中型滑坡。

根據(jù)現(xiàn)場(chǎng)調(diào)查及鉆探資料，滑體土厚度在7.5～33.3 m 范圍之間，滑體中間厚，前緣較薄，主要由第四系中更新統(tǒng)離石黃土（Qp2eol）和上更新統(tǒng)馬蘭黃土（Qp3eol）及少量紅黏土組成，土體覆蓋于侏羅系砂巖（J2y）之上。

2 模型建立

根據(jù)滑坡主剖面地層狀況和相關(guān)地形特征，按照實(shí)際尺寸確定模擬剖面的模型大小。依據(jù)勘察資料，將滑坡的地層劃分為馬蘭黃土、離石黃土、紅黏土、滑帶土、滑體土以及滑床六大部分，并對(duì)相應(yīng)地層進(jìn)行合理的網(wǎng)格劃分，從而得到該滑坡穩(wěn)定性分析模擬的網(wǎng)格模型圖（見(jiàn)圖1）。

圖1 模型網(wǎng)格劃分圖

3 模擬工況及材料賦值

由于滑坡所在區(qū)域降雨時(shí)間分布極為不均，雨季降雨量較大并常伴有暴雨極端天氣，故滑坡模擬計(jì)算采用天然工況和飽和工況兩種狀態(tài)。

通過(guò)對(duì)地質(zhì)勘察取樣進(jìn)行室內(nèi)試驗(yàn)數(shù)據(jù)測(cè)定，并結(jié)合相關(guān)工程類比和簡(jiǎn)化模型綜合得到各地層參數(shù)（見(jiàn)表1），作為用于數(shù)值模擬計(jì)算的相關(guān)賦值計(jì)算參數(shù)。

表1 材料參數(shù)賦值

4 基于Midas GTS / NX 的滑坡穩(wěn)定性分析

滑坡穩(wěn)定性狀態(tài)評(píng)判按照《滑坡防治工程勘查規(guī)范》(GBT32864-2016），數(shù)值模擬穩(wěn)定性分析計(jì)算采用天然工況和飽和工況兩種狀態(tài)。

4.1 天然工況

在天然工況下，由有限元數(shù)值模擬計(jì)算得出剖面的穩(wěn)定系數(shù)為1.024 3，處于欠穩(wěn)定狀態(tài)。

從天然工況滑坡的位移云圖（見(jiàn)圖2），我們可以較清楚地看到位移最大出現(xiàn)在滑坡后緣附近，這主要是因?yàn)榛潞缶壚鄯e變形，在滑坡后部形成較多裂縫導(dǎo)致。而滑坡前部坡腳產(chǎn)生的位移，是由于該段滑體比較薄，容易發(fā)生應(yīng)力集中[6]，這與天然狀態(tài)下坡體破壞特征相符合。

圖2 天然工況位移云圖

從剪應(yīng)變?cè)茍D（見(jiàn)圖3）可知，滑坡后部滑面為切割土體形成，而中前部滑面為滑體與基巖接觸面，滑面目前尚未完全貫通。

圖3 天然工況剪應(yīng)變?cè)茍D

4.2 飽和工況

在飽和工況下，據(jù)有限元數(shù)值模擬計(jì)算得出剖面的穩(wěn)定系數(shù)為0.857 6，穩(wěn)定性較天然狀態(tài)明顯降低，滑坡處于不穩(wěn)定狀態(tài)。

從飽和工況下滑坡的位移云圖（見(jiàn)圖4），可以看到位移與天然工況相比，位移量顯著增加，滑坡已經(jīng)發(fā)生變形破壞，整個(gè)滑坡已處于不穩(wěn)定狀態(tài)。

圖4 飽和工況位移云圖

從剪應(yīng)變?cè)茍D（見(jiàn)圖5）可以明顯的看出，剪應(yīng)變?cè)陲柡凸r下擴(kuò)展并貫通了整個(gè)滑坡體，形成了貫通的滑面，說(shuō)明在降水持續(xù)入滲過(guò)程中，黃土內(nèi)部與基巖接觸面附近產(chǎn)生了貫通滑坡體的軟弱帶，隨著軟弱帶土體抗剪強(qiáng)度大幅度降低，最終滑坡發(fā)生剪切破壞貫通滑面，造成滑坡發(fā)生滑動(dòng)。

圖5 飽和工況剪應(yīng)變?cè)茍D

綜合分析得出該滑坡天然狀態(tài)下處于欠穩(wěn)定狀態(tài)，在降雨等其他因素作用影響下，滑坡將處于不穩(wěn)定狀態(tài)，隨時(shí)可能發(fā)生滑動(dòng)，需要及時(shí)采取相關(guān)治理措施進(jìn)行加固治理。

5 基于Midas GTS / NX 的滑坡治理方案可行性分析

5.1 擬采用治理方案

根據(jù)滑坡的地形地貌以及工程地質(zhì)特征，并結(jié)合滑坡穩(wěn)定性分析，設(shè)計(jì)治理方案為：抗滑樁+擋土板+平整坡面及削坡+截排水

（1） 抗滑樁

根據(jù)該滑坡工程地質(zhì)特征狀況，在滑坡前部設(shè)置一排抗滑樁，抗滑樁規(guī)格2.5×3 m，樁長(zhǎng)21 m。

（2） 擋土板

由于坡體整體結(jié)構(gòu)較為疏松，為防止相鄰抗滑樁之間土體受坡體自身擠壓作用溜出，在相鄰抗滑樁間設(shè)計(jì)擋土板。

（3） 平整坡面及削坡

根據(jù)《建筑邊坡工程技術(shù)規(guī)范》，按照設(shè)計(jì)要求，從上到下對(duì)滑坡進(jìn)行削坡處理，對(duì)坡體抗滑樁以上坡體按照坡比1:1 進(jìn)行削坡，每級(jí)高度8 m，護(hù)坡道平臺(tái)寬3 m，并對(duì)坡體裂縫及落水洞進(jìn)行回填夯實(shí)處理，確保坡面完整性。

（4） 排水工程

根據(jù)研究，隨著水體入滲坡體，土體容重隨之增大，坡體下滑力增加，穩(wěn)定性降低；另一方面由于黃土整體結(jié)構(gòu)疏松多孔，垂直節(jié)理發(fā)育，坡面徑流、匯水等侵蝕破壞坡面，并沿土體裂隙下滲，軟化土體形成軟弱帶[7]，做好排水是有效治理滑坡的關(guān)鍵。

5.2 治理方案可行性分析

治理方案均嚴(yán)格按照邊坡治理規(guī)范設(shè)計(jì)，在天然狀態(tài)下設(shè)計(jì)方案安全系數(shù)不小于1.3，暴雨?duì)顟B(tài)下不小于1.15，為了檢驗(yàn)治理方案的極限穩(wěn)定性，模擬選擇在暴雨飽和工況下的治理方案的滑坡穩(wěn)定性數(shù)值模擬。材料參數(shù)采取飽和狀態(tài)下的參數(shù)進(jìn)行模擬計(jì)算分析。

在飽和工況下，治理方案模擬結(jié)果穩(wěn)定性系數(shù)為1.183，與治理前飽和工況相比，穩(wěn)定性得到了顯著的提高，滿足飽和工況安全系數(shù)要求，證明該方案能夠達(dá)到滑坡治理穩(wěn)定性要求。

從治理方案的總位移云圖（見(jiàn)圖6）我們可以看出，采取治理方案之后，最大位移雖出現(xiàn)在坡體后緣，但整體位移量較采取方案前大幅度降低，證明治理方案對(duì)該滑坡治理效果明顯。

圖6 治理方案位移云圖

從治理方案剪應(yīng)變?cè)茍D（見(jiàn)圖7）可以看出，采取治理方案后，滑坡最大剪應(yīng)變主要出現(xiàn)在抗滑樁樁后附近，與采取方案之前相比，剪應(yīng)變值明顯減小，并且不再出現(xiàn)貫通，證明該治理方案加固效果較好，治理方案可行。

圖7 治理方案剪應(yīng)變?cè)茍D

6 結(jié)論

本研究主要借助于Midas GTS / NX 有限元模擬軟件，對(duì)滑坡在天然工況和飽和工況分別模擬分析，得出天然狀態(tài)下滑坡的穩(wěn)定系數(shù)為1.024 3，飽和工況下的穩(wěn)定系數(shù)為0.857 6。從模擬結(jié)果分析來(lái)看，滑坡從天然狀態(tài)過(guò)渡到飽和狀態(tài)的過(guò)程中，滑坡整體位移量逐漸增大，并且坡體內(nèi)部逐漸形成貫通的滑面；在采取治理方案后，受前部抗滑樁支撐，滑坡位移量得到明顯控制，且不再形成貫通滑面，滑坡穩(wěn)定性系數(shù)增大，能夠滿足規(guī)范治理安全系數(shù)要求，該治理方案是可行的。通過(guò)Midas GTS / NX 有限元模擬軟件能夠較快速準(zhǔn)確對(duì)滑坡不同狀態(tài)下的穩(wěn)定性進(jìn)行快速分析判斷，并且能夠有效的對(duì)治理方案的治理效果進(jìn)行比對(duì)分析，為相關(guān)工程提供一定的借鑒。