摘 要：針對降雨入滲條件下非飽和土堤的穩(wěn)定性，文章采用Fortran軟件結(jié)合強度折減系數(shù)原理以某段非飽和土堤為研究對象，考慮水力滲透系數(shù)特征曲線、修正的Mohr-coulomb破壞準(zhǔn)則，建立了一個考慮非飽和滲流的應(yīng)力應(yīng)變模型。根據(jù)瞬時安全系數(shù)和抗滑薄弱區(qū)分布綜合評判土堤的抗滑穩(wěn)定性，得到降雨作用下土堤穩(wěn)定性的若干規(guī)律。從分析結(jié)果來看：隨著降雨的持續(xù)入滲，土體抗剪強度大幅度下降，抗滑塑性區(qū)始于堤腳并不斷擴大，當(dāng)降雨強度超過一定范圍時，土堤穩(wěn)定性與降雨強度的關(guān)系不再明顯。

關(guān)鍵詞：降雨入滲；非飽和土堤；體積含水率；強度折減

降雨開始后，雨水滲透浸潤線不斷向土堤下部擴展，并在滲透路徑較短的堤腳處率先形成飽和區(qū)域，堤腳處開始形成塑性區(qū)，圖2（a）。但坡腳節(jié)點位移增長緩慢；在17h后，上游流域地面土體大致飽和，地面徑流使河道水位抬高，加之土堤上部雨水不斷向下部浸潤補充，堤腳極易形成超飽和區(qū)，并且這種超飽和區(qū)開始向內(nèi)部擴展，孔隙水壓力增大，基質(zhì)吸力減少，節(jié)點位移迅速增加，塑性區(qū)向上延伸，如圖2（b）；在50h后，雨水浸潤線接近洪水浸潤線，土體飽和趨于穩(wěn)定，坡腳位移雖有增加但增長緩慢，此時土堤內(nèi)部的塑性區(qū)接近堤頂，圖2（c），折減系數(shù)達(dá)到最小1.43，并從圖中可以清晰看到土堤邊坡的抗滑薄弱面始于堤腳大致呈圓弧狀向堤頂發(fā)展，這與瑞典分條法的假定是一致的。

（a） t=17h

（b） t=45h

（c） t=60h

圖2 降雨時間對土堤薄弱區(qū)發(fā)展的影響

土堤穩(wěn)定性不僅與降雨時長密切相關(guān)，同時又受降雨強度的影響。為研究降雨強度對土堤穩(wěn)定性的影響，分別取R=0.42mm/h（小雨），1.24mm/h（中雨），2.49mm/h（大雨），5.83mm/h（暴雨），11.66mm/h（大暴雨），根據(jù)式（1）（2）改變模型透水邊界的初始入滲流量值。為控制河道內(nèi)洪水水位的影響，假設(shè)各不同降雨強度的水位在同時刻點相同并隨時間逐級增加，計算各工況下的折減系數(shù)F與降雨強度R關(guān)系如圖3。當(dāng)降雨強度小于土體滲透系數(shù)k時，土堤表面的入滲速率取決于降雨強度R，因此不同降雨強度的曲線差異較顯著；而當(dāng)降雨強度繼續(xù)增大并大于土體滲透系數(shù)k時，土堤表層土的飽和度并不會繼續(xù)增加，而是以大致相同的速度向下推移，向下推移的速度取決于土體的飽和滲透系數(shù)，因此不同降雨強度的折減系數(shù)曲線差異較小。

為驗證文章所采用方法的有效性，取降雨強度R=11.66mm/h，在已分析的抗滑薄弱面的基礎(chǔ)上，采用基于極限平衡理論的瑞典分條法計算各對應(yīng)工況的穩(wěn)定系數(shù)如表2，對比發(fā)現(xiàn)兩種方法吻合較好，最大差值為8.9%。不難發(fā)現(xiàn)傳統(tǒng)方法較有限元方法偏大，主要是因為有限元方法考慮了土堤模型的初始應(yīng)力、土體材料的應(yīng)力應(yīng)變關(guān)系及坡腳集中滲流所產(chǎn)生的滲流力等因素。

表2 不同方法安全系數(shù)對比

3 結(jié)束語

文章基于有限元理論，改進(jìn)了傳統(tǒng)極限平衡理論對土堤穩(wěn)定分析的不足，建立一個考慮土-水耦合情況下的滲流模型，分析降雨強度對土堤穩(wěn)定分析的影響規(guī)律。

（1）降雨入滲作用下土堤非飽和區(qū)體積含水率增加，引起土的抗剪強度降低，從而影響土堤邊坡的穩(wěn)定性。

（2）降雨條件下坡腳處易產(chǎn)生高飽和區(qū)，在滲透力作用下形成不穩(wěn)定的抗滑塑性區(qū)，并向內(nèi)部擴展。因此，在土堤施工和運行階段，要做好土堤的排水工程措施。

（3）降雨強度越大，對土堤穩(wěn)定越不利，但當(dāng)降雨強度繼續(xù)增大超過土體滲透系數(shù)時，這種差異不再明顯，而是更多的取決于降雨歷時。

參考文獻(xiàn)

[1]周建，余嘉澍.防洪堤穩(wěn)定性的研究[J].水利學(xué)報，2002，33（7）：98-103.

[2]曾理，錢曉慧，祁小輝.霧化雨入滲條件下高拱壩下游邊坡穩(wěn)定性分析[J].人民黃河，2011，33（6）：129-131.

[3]陳善雄，譚新，柳治國.降雨條件下土質(zhì)邊坡穩(wěn)定性預(yù)測預(yù)報方法[J].巖土力學(xué)，2002，23（增刊1）：31-36.

[4]魏保義.用非飽和滲流理論研究降雨入滲的特性及對邊坡穩(wěn)定的影響[D].南京：河海大學(xué)，2001.

[5]徐晗，朱以文，蔡元奇，等.降雨入滲條件下非飽和土邊坡穩(wěn)定分析[J].巖土力學(xué)，2005，26：1957-1962.

[6]Cho S E， Lee S R. Instability of unsaturated soil slopes due to infiltration[J].Computers and Geotechnics，2001.

[7]錢家歡.土力學(xué)[M].南京：河海大學(xué)出版社，1995：17-21.

[8]吳宏偉，陳守義，龐宇威.雨水入滲對非飽和土坡穩(wěn)定性影響的參數(shù)研究[J].巖土力學(xué)，1999，01：001-013.

[9]Zienkiewice O C，Humpheson Cand Lewis RW.Associated and Non-associated Visco-plasticity and Plasticity in Soil Mechanics[J].Geotechnique，1975，25（4）：671-689.

[10]張曉岸，郭小剛，王晶晶，等.降雨入滲對邊坡穩(wěn)定性影響的強度折減法[J].工程建設(shè)與設(shè)計，2014，06：88-92.

作者簡介：孫文濤（1991-），男，漢族，河南商丘人，碩士研究生，華北水利水電大學(xué)，研究方向：水工結(jié)構(gòu)數(shù)值分析及模型試驗。