降雨入滲對土質路塹邊坡的穩定性影響分析

姚志勇

摘要：為了研究土質路塹邊坡在降雨入滲條件下的穩定性，文章以廣西桂三高速公路高架橋右側某土質路塹邊坡為研究對象，采用數值分析軟件，建立相應的數值分析模型，描述了土質路塹邊坡在降雨入滲條件下的滲流場演變規律，揭示了滲流作用下土質路塹邊坡的穩定性變化規律。研究結果表明：降雨入滲條件下，邊坡土體體積含水率隨降雨入滲量的不斷累積而逐漸增加;受滲流演變規律的影響，塑性區在邊坡內部形成，并隨降雨歷時的增加，逐漸向邊坡內部拓展;研究歷時內，邊坡安全系數從1.56減小到了1.12，邊坡的穩定性受降雨影響較大。

關鍵詞：降雨入滲;路塹邊坡;滲流作用;穩定性

中圖分類號：U416.14 文獻標識碼：A DOI：10.13282/j.cnki.wccst.2019.10.017

文章編號：1673-4874（2019）10-0058-03

0引言

在道路建設過程中，受路線整體規劃的要求，常常會對原有土基進行開挖，形成大量的土質路塹邊坡工程。在降雨入滲條件下，該類邊坡土體強度會逐漸劣化，導致邊坡失穩事件時有發生。目前國內外學者針對土質路塹邊坡的穩定性已經進行了相關研究：張漣英等在對依托工程進行現場調研的基礎上，提出了紅黏土邊坡的穩定性計算方法，并利用該方法對紅黏土邊坡進行了穩定性計算;張坤等對黃土路塹邊坡土體的水熱變化規律進行了監測，根據監測結果，分析了干濕循環和凍融循環條件下路塹邊坡內的體積含水率及溫度變化規律;喬翔根據典型的高路塹邊坡失穩案例，分析了誘發高路塹邊坡發生滑坡破壞的影響因素，并提出了相應的加固方案;侯志峰等以典型的黃土路塹邊坡為研究對象，采用數值軟件，對黃土路塹邊坡進行了數值模擬，揭示了黃土路塹邊坡的變形演變規律。從上述研究中可以看出，目前關于路塹邊坡穩定性的研究主要集中在穩定性計算及加固技術這兩個方面，較少有學者從滲流演變角度出發，考慮降雨入滲對土質路塹邊坡穩定性的影響。

鑒于此，本文針對目前國內外關于降雨入滲條件下土質路塹邊坡穩定性影響研究的不足，以廣西桂三高速公路高架橋右側某土質路塹邊坡為研究對象，采用Geo-studio軟件，建立相應的數值分析模型，通過數值模擬土質路塹邊坡在降雨入滲條件下的孔隙水壓力、體積含水率等滲流演變特征，揭示降雨入滲對土質路塹邊坡的穩定性影響機制，為工程實踐提供理論參考。

1土質路塹邊坡降雨入滲機理及穩定性分析理論

1.1滲流邊界

雨水在土質路塹邊坡表面入滲的情況如圖1所示。根據圖1所示的雨水入滲特性，可將滲流邊界條件分為3種。

第一種邊界條件又稱作水頭邊界條件。通過在邊界上指定水頭分布，在降雨強度大于土體滲透系數時，雨水在邊坡表面入滲的條件為第一種邊界條件，如圖1中BC段中B對應的階段所示，具體表達式為：

h|=H（x，y，z，t） （1）

第二種邊界條件又稱為流量邊界條件。該情況下，當降雨強度小于土體滲透系數，對應于圖1中的A階段，具體表達式為：

式中：q——一流過邊界單位面積上的流量，等同于v;l、l、l。指外法線n與坐標之間的方向余弦。

第三種邊界條件又稱為混合邊界條件。該情況下，土質路塹邊坡表面水頭差與邊坡內部交換的流量呈函數關系，雨水在邊坡表面產生徑流。如圖1中BC段中C對應于該階段，其中線D則表示了土體的入滲能力。具體表達式為：

式中：α——正常數，α與β均是此種邊界中各點的已知數。

1.2 非飽和抗剪強度理論

土質路塹邊坡的非飽和抗剪強度計算方法采用雙應力變量（υ-μ，μ-μ）理論基礎上發展而來的非飽和抗剪強度公式：

τ=c+（σ-μ）tan+（μ-μ）tan （5）式中：μ——粗粒土土體的進氣值，由于考慮邊坡孔隙與大氣相連μ=0;

ψ——基質吸力對內摩擦角的貢獻，按經驗取值ψ=14°;

μ——孔隙水壓力。

2 土質路塹邊坡計算參數及數值模型

2.1計算參數

本文計算的土質路塹邊坡位于廣西桂三高速公路高架橋右側。對路塹邊坡進行現場取樣，在室內進行變水頭滲透試驗，測得邊坡土體滲透系數為：1.35×10cm/s;通過直剪試驗測得土樣的粘聚力為27kPa，內摩擦角為21°。此外，參考相關文獻，確定土樣的物理力學參數如表1所示。

根據當地歷年降雨情況，考慮本文主要研究滲流影響下的路塹邊坡穩定性，選取本次數值模擬的降雨強度為9.1×10mm/S.降雨時長為24h。通過室內試驗可得土樣飽和體積含水率和殘余體積含水率分別為0.18、0.10，在Geosudio軟件的seep模塊中，利用Von Genuchten模型對土樣的土水特征曲線進行擬合，如圖2和圖3所示。

2.2 數值計算模型的建立

根據現場裯研資料可知，土質路塹邊坡總高度達20m，分兩級開挖，第一級邊坡坡比為：1：1.75，第二級邊坡坡比為1：1.5。為了分析降雨入滲條件下土質路塹邊坡的滲流演變規律，選取了兩個監測截面。模型邊界設定為：邊坡表面及坡頂設定為降雨邊界，模型底部及兩側邊界設定為不透水邊界。見圖4。

3數值模擬結果

本算例采用摩爾庫倫模型，在降雨入滲條件下對土質路塹邊坡的滲流特性進行計算，降雨24h時邊坡監測截面1和截面2的體積含水率變化規律如圖5所示。

從圖5中可以看出：在截面1處，高程在19-20m內，體積含水率基本達到飽和狀態，高程在13-19m內體積含水率隨高程的減小而逐漸減小;在截面2處，高程在30-28m內，體積含水率基本達到飽和狀態，高程在13-28m內，體積含水率隨高程的減小而逐漸減小。高程在5-10m內，兩個監測截面體積含水率隨高程減小而逐漸增加;高程在0-5m內，兩個監測截面體積含水率維持在飽和狀態。這主要是因為在高程5m處，土質路塹邊坡內存在地下水位，因此，高程5m內的體積含水率基本都維持在飽和狀態，受地下水位的影響，地下水位以上的土體在基質吸力的作用下，對地下水具有毛細作用。因此，在一定高程內，土體體積含水率隨高程的增加逐漸減小;受降雨入滲的影響，邊坡表面土體基本達到飽和狀態，下部土體受入滲雨水的影響，體積含水率也有所增加。

將數值模擬的滲流特性通過FiSh語言導入到FLAC3D軟件中，計算得到的土質路塹邊坡的塑性區如圖6昕示.

從圖6中可以看出：降雨入滲條件下，雨水在路塹邊坡內部滲流，使得邊坡內部出現塑性區，隨著降雨歷時的增加，塑性區面積逐漸增大，在第一邊坡坡腳處的塑性區面積也逐漸向邊坡內部拓展，直至形成貫通區，導致不能再繼續承受外部荷載或自重力，最終出現沿塑性區發育的失穩現象。

運用強度折減法計算的邊坡安全系數變化規律如圖7所示。

從圖7中可以看出：土質路塹邊坡在降雨入滲條件下，安全系數從最初的1.56減小到了1.12，減小了28%，平均每小時減小0.018。可見，在降雨入滲作用下，雨水從路塹邊坡表面入滲后，在邊坡內部進行滲流，形成滲流場。隨著降雨入滲量的不斷累積，滲流場逐漸向邊坡內部演變，體積含水率逐漸增加，土體抗剪強度逐漸降低，導致路塹邊坡的整體穩定性逐漸降低。

4 結語

本文從滲流演變的角度，深入分析了土質路塹邊坡在降雨入滲條件下的滲流特性，根據滲流演變規律，采用數值軟件，對滲流影響下的土質路塹邊坡的穩定性進行了計算，揭示了土質路塹邊坡在降雨入滲條件下的穩定性變化規律，得到以下主要結論：

（1）降雨入滲條件下，邊坡表面土體體積含水率隨降雨入滲量的不斷累積而逐漸增加至飽和狀態，邊坡內部土體體積含水率受滲流場演變規律的影響。

（2）受滲流演變規律的影響，塑性區在邊坡內部形成，并隨降雨歷時的增加，塑性區面積逐漸擴展，直至貫通。

（3）土質路塹邊坡在降雨入滲條件下，安全系數隨降雨歷時的增加而逐漸減小。