強(qiáng)降雨條件下非飽和土邊坡穩(wěn)定性研究

華月新

(惠州市華禹水利水電工程勘測設(shè)計(jì)有限公司，廣東 惠州 516003)

滑坡是指在一定地形地質(zhì)條件下，由于外界條件發(fā)生變化，導(dǎo)致原狀邊坡土體應(yīng)力場發(fā)生改變，破壞了原有的力學(xué)平衡條件，使邊坡在自重和其他荷載聯(lián)合作用下，沿滑裂面突發(fā)性地向下移動(dòng)[1-3]。我國幅員遼闊，約2/3為山地，受地質(zhì)構(gòu)造和氣候條件的影響，每年發(fā)生滑坡災(zāi)害數(shù)以萬計(jì)，西南地區(qū)尤為嚴(yán)重，給人民群眾的生命財(cái)產(chǎn)安全造成巨大損失，降雨條件下邊坡穩(wěn)定性分析已成為工程界研究的熱點(diǎn)問題[4]。近年來，隨著我國基礎(chǔ)設(shè)施建設(shè)的快步推進(jìn)，其施工過程往往不可避免地對原狀邊坡產(chǎn)生擾動(dòng)，在遭遇降雨等不利因素影響時(shí)，極易導(dǎo)致滑坡災(zāi)害的發(fā)生，因此有必要對降雨條件下的邊坡穩(wěn)定性分析開展研究[5-7]。本文以均質(zhì)土坡為例，考慮降雨入滲對邊坡應(yīng)力場的影響，采用有限元強(qiáng)度折減法，對邊坡進(jìn)行穩(wěn)定性分析，對不同降雨條件下邊坡安全系數(shù)的變化規(guī)律進(jìn)行研究，研究結(jié)果可為邊坡抗滑穩(wěn)定分析與設(shè)計(jì)提供參考。

1 計(jì)算理論與方法

1.1 流固耦合模型

利用Galerkin 有限元方法，離散節(jié)點(diǎn)孔壓和位移，可得流固耦合控制方程[8]：

(1)

1.2 飽和非飽和理論

降雨在非飽和土中滲流滿足質(zhì)量守恒方程，以達(dá)西定律為基礎(chǔ)，結(jié)合滲流過程中流體質(zhì)量守恒，推導(dǎo)出以基質(zhì)勢為因變量的滲流控制方程[9]：

(2)

(3)

式中：kx、ky分別為X、Y方向的滲透系數(shù)；h為基質(zhì)勢；ρw為水的密度；g為重力加速度；mw為比水容重。

1.3 強(qiáng)度折減理論

強(qiáng)度折減法主要思想是將原始土體的強(qiáng)度參數(shù)c和φ按照一定比例進(jìn)行折減，將折減后的強(qiáng)度參數(shù)作為新的材料輸入進(jìn)行計(jì)算。以此類推，直到邊坡塑性區(qū)貫通或坡頂相對位移急劇增大時(shí)停止，此時(shí)折減系數(shù)即為邊坡安全系數(shù)Ks?；谀枎靷惸Ｐ蛷?qiáng)度折減公式為[10]：

tanφ′=tanφ/k

(4)

c′=c/k

(5)

式中：k為強(qiáng)度折減系數(shù)；c、φ為土體實(shí)際抗剪強(qiáng)度指標(biāo)；c′、φ′為折減后的強(qiáng)度參數(shù)。

2 降雨入滲條件下非飽和土邊坡穩(wěn)定性分析

2.1 有限元計(jì)算模型

強(qiáng)降雨入滲條件下，邊坡土體快速由非飽和狀態(tài)進(jìn)入飽和狀態(tài)，相應(yīng)區(qū)域孔隙水壓力會(huì)急劇升高，同時(shí)抗剪強(qiáng)度會(huì)有所降低，瞬態(tài)滲流場和應(yīng)力場耦合作用下，邊坡更容易發(fā)生失穩(wěn)破壞。本文針對以均值土坡為例，對不同降雨條件下邊坡安全系數(shù)的變化規(guī)律進(jìn)行研究。邊坡有限元模型見圖1，A、B、C分別為坡頂、坡中和坡腳處3個(gè)位移觀測點(diǎn)。邊坡坡角45°，坡腳長15 m, 坡頂長25 m, 高度25 m，初始水位位于坡腳處。模型共劃分1 585個(gè)單元和1 687個(gè)節(jié)點(diǎn)，采用ABAQUS孔壓應(yīng)力耦合CPE4P單元模擬。初始應(yīng)力場為地下水和自重引起的應(yīng)力場。土體彈性模量為200 MPa，泊松比0.3，密度為1 300 kg/m3。凝聚力為19 kPa，摩擦角為27°，滲透系數(shù)為0.44E-6 m/s。土體孔壓和滲透折減系數(shù)隨飽和度變化見圖2。

圖1 二維邊坡模型

圖2 孔壓和滲透系數(shù)隨飽和度變化曲線

2.2 不同降雨持時(shí)對邊坡穩(wěn)定性影響

為研究不同降雨持時(shí)對邊坡穩(wěn)定性影響，取雨強(qiáng)為20 mm/h，分別設(shè)置降雨持時(shí)0、10、20、30、40、50、60和70 h共8種工況，用有限元強(qiáng)度折減法對不同工況下邊坡穩(wěn)定性進(jìn)行分析。圖3為降雨持時(shí)分別為0、30和60 h 3種工況下的邊坡孔隙水壓力分布云圖。從圖3中可以看出，3種工況下，邊坡孔隙水壓力最大值分別為150、157.5和167.7 kPa，最小值分別為-250、-59.7和-53.5 kPa。隨著降雨持時(shí)的增加，邊坡孔隙水壓力逐漸增大，其中邊坡上部負(fù)孔隙水壓力在降雨初期變化顯著，隨著降雨時(shí)間增加則趨于穩(wěn)定。

圖3 孔隙水壓力分布云圖(不同降雨持時(shí))

圖4為3種工況下的邊坡飽和度分布云圖。從圖4中可以看出，3種工況下，邊坡飽和度最小值分別為0.08、0.091和0.101。隨著降雨持時(shí)的增加，邊坡土體飽和度呈上升趨勢。對比3種工況下飽和度云圖可以發(fā)現(xiàn)，隨著降雨持時(shí)增加，邊坡下部土體逐漸趨于飽和，上部非飽和土體部分區(qū)域逐漸縮小。

圖4 邊坡飽和度分布云圖(不同降雨持時(shí))

圖5為邊坡觀測點(diǎn)水平向位移隨降雨持時(shí)變化曲線。從圖5中可以看出，其水平向位移隨觀測點(diǎn)高程增加而增大。在降雨0～5 h以內(nèi)，A、B兩觀測點(diǎn)位移隨降雨持時(shí)變化較為陡峭。這是由于在降雨初期，邊坡上部土體孔隙水壓力變化較大；隨著降雨持時(shí)增加，邊坡上部土體孔隙水壓力逐漸趨于穩(wěn)定。當(dāng)降雨持時(shí)達(dá)到30 h左右時(shí)，3個(gè)觀測點(diǎn)水平向位移變化速率均發(fā)生明顯變化。圖6為邊坡安全系數(shù)隨降雨持時(shí)變化曲線。從圖6中可以看出，在降雨0～10 h以內(nèi)，邊坡安全系數(shù)降低較快。當(dāng)降雨持時(shí)大于10 h時(shí)，邊坡安全系數(shù)隨降雨持時(shí)增加線性變化。

圖5 觀測點(diǎn)水平向位移隨降雨持時(shí)變化曲線

圖6 邊坡安全系數(shù)隨降雨持時(shí)變化曲線

為研究不同降雨強(qiáng)度對邊坡穩(wěn)定性影響，設(shè)置降雨持時(shí)30 h，分別設(shè)置降雨強(qiáng)度5、10、15、20、25、30和35 mm/h 7種工況，對邊坡穩(wěn)定性進(jìn)行分析。圖7為降雨強(qiáng)度分別為5、20和35 mm/h 3種工況下的邊坡孔隙水壓力分布云圖。從圖7中可以看出，3種工況下，邊坡孔隙水壓力最大值分別為152.7、157.5和162.9 kPa，最小值分別為-206.7、-59.7和-50.6 kPa。隨著降雨強(qiáng)度的增加，邊坡孔隙水壓力逐漸增大，其中邊坡上部負(fù)孔隙水壓力隨著降雨強(qiáng)度增加變化速率由大到小。

圖7 孔隙水壓力分布云圖(不同降雨強(qiáng)度)

圖8為3種工況下的邊坡飽和度分布云圖。從圖8中可以看出，3種工況下，邊坡飽和度最小值分別為0.08、0.091和0.106。隨著降雨持時(shí)的增加，邊坡土體飽和度呈上升趨勢。對比3種工況下飽和度云圖可以發(fā)現(xiàn)，隨著降雨強(qiáng)度增加，邊坡下部土體逐漸趨于飽和，上部非飽和土體部分區(qū)域逐漸縮小。

圖8 邊坡飽和度分布云圖(不同降雨強(qiáng)度)

圖9為不同降雨強(qiáng)度下邊坡觀測點(diǎn)水平向位移隨降雨持時(shí)變化曲線。從圖9中可以看出，其水平向位移隨降雨強(qiáng)度增加而增大。對比3種工況下位移變化速率可知，邊坡位移變化速率隨降雨強(qiáng)度增加而不斷增大。圖10為邊坡安全系數(shù)隨降雨強(qiáng)度變化曲線。總體來看，邊坡安全系數(shù)隨降雨強(qiáng)度增加而逐漸減小，兩者呈線性變化。

圖9 不同降雨強(qiáng)度下邊坡觀測點(diǎn)水平向位移隨降雨持時(shí)變化曲線

圖10 邊坡安全系數(shù)隨降雨強(qiáng)度變化曲線

3 結(jié) 語

本文基于非飽和土體理論，并結(jié)合有限元強(qiáng)度折減法，對不同降雨持時(shí)和降雨強(qiáng)度作用下的邊坡穩(wěn)定性進(jìn)行研究。研究結(jié)果表明，在強(qiáng)降雨條件下，邊坡上部負(fù)孔隙水壓力在降雨初期變化顯著，隨著降雨時(shí)間增加則趨于穩(wěn)定。在降雨0～5 h以內(nèi)，邊坡土體上部位移速率較大。在降雨0～10 h以內(nèi)，邊坡安全系數(shù)降低較快。當(dāng)降雨持時(shí)大于10 h時(shí)，邊坡安全系數(shù)隨降雨持時(shí)增加線性變化。隨著降雨強(qiáng)度增加，邊坡位移變化速率隨逐漸增大。總體來看，邊坡安全系數(shù)隨降雨強(qiáng)度增加而逐漸減小，兩者呈線性變化。