李 濤，付宏淵，周功科，莫 凱，曾 鈴

（1.長沙理工大學交通運輸工程學院，湖南 長沙 410004；2.浙江通衢交通建設監(jiān)理咨詢有限公司，浙江 衢州 324000；3.長沙理工大學土木與建筑學院，湖南 長沙 410004）

公路路堤穩(wěn)定性問題歷來受到人們的關注，盡管在路堤變形規(guī)律、穩(wěn)定性、填筑材料與施工工藝等方面開展了大量的工作［1～3］，對于路堤而言，填料含水量控制著它的密度和強度，含水量變大會使路堤填料的重度與孔隙水壓力增加，抗剪強度降低，最終導致路堤變形與失穩(wěn)，許多學者對此進行了深入的研究。王瑞鋼［4］等采用等效粘聚力的概念，利用延伸的MOHRCOULOMB破壞準則對路堤邊坡進行彈塑性有限元分析，指出在降雨作用下，路堤邊坡會發(fā)生局部失穩(wěn)問題，降雨重現(xiàn)期、土參數(shù)φb與各向異性滲透系數(shù)比對路堤邊坡滲流穩(wěn)定均有較大的影響。陳曉斌［5］等對降雨入滲作用下粗粒土路堤變形及穩(wěn)定性進行了計算分析，結果表明:當單次入滲深度為2m時，路堤產(chǎn)生的附加沉降可達1.30cm，且對路堤穩(wěn)定性的影響也比較明顯。李聰［6］在降雨條件下進行了壓實度不同的兩組路堤室內(nèi)模型試驗，并指出突然降雨或長期暴雨對路堤土壓力和變形影響顯著，入滲深度與填土類型和粒徑有關。李汝成［7］等通過現(xiàn)場試驗與數(shù)值分析相結合，揭示了降雨對泥巖-土混填路堤邊坡穩(wěn)定影響規(guī)律，并指出裂縫的出現(xiàn)對路堤穩(wěn)定性安全系數(shù)的影響顯著。

現(xiàn)有研究成果表明，人們盡管認識到降雨入滲引起路堤內(nèi)部含水率的增加是導致路堤變形與失穩(wěn)的主要因素，但是對于填方路堤而言，其飽和-非飽和滲流過程受降雨量、降雨時間、表面徑流量、填料性質(zhì)與壓實度等多種因素的影響［8～10］，其滲流過程較邊坡而言更為復雜。目前，公路路堤降雨入滲方面主要以數(shù)值模擬與模型實驗相結合對其進行分析，且研究內(nèi)容主要側重于降雨入滲引起路堤變形與穩(wěn)定性兩方面，而較少對雨水在路堤內(nèi)部的入滲過程進行探討。因此，對路堤降雨入滲過程進行研究，不僅能更加全面地揭示路堤內(nèi)部的雨水入滲規(guī)律，還能對路堤變形與穩(wěn)定性研究提供理論基礎，以此來更加深入地揭示降雨引起路堤邊坡失穩(wěn)的本質(zhì)原因和外在因素。圖1為路堤降雨入滲流向示意圖，當降雨強度大于滲透系數(shù)時，降雨量將由坡面入滲量與徑流量組成，而降雨強度小于滲透系數(shù)時，邊坡表面則不會形成坡面徑流。

在前人研究的基礎上，本文采用數(shù)值模擬分析方法，結合實際降雨資料，對廣西六寨至河池高速公路K19+133路堤斷面降雨入滲過程中路堤內(nèi)部滲流場以及穩(wěn)定性進行研究，以期獲得路堤內(nèi)部暫態(tài)飽和區(qū)的變化規(guī)律及路堤失穩(wěn)機理，為降雨影響下的粗粒土路堤邊坡設計與治理提供參考。

圖1 路堤降雨流向示意圖Fig.1 Sketch map of rainfall flow direction of an embankment

1 飽和-非飽和滲流理論

對于公路路堤而言，經(jīng)過碾壓密實后的路堤一般處于非飽和狀態(tài)，其降雨入滲屬于典型的飽和-非飽和過程，采用經(jīng)典的飽和土理論無法真實的反映路堤土的入滲過程。此外，由于路堤填筑碾壓過程造成了其內(nèi)部土體的各向異性，表現(xiàn)為豎向滲透系數(shù)小于水平方向［11］，即Ky/Kx＜1。在本文的分析中，降雨條件下路堤內(nèi)部飽和-非飽和滲流過程的控制偏微分方程如下［4］:

式中:kx、ky——土體水平和豎直方向滲透系數(shù)(m/s)；

Q——邊界降雨量(m3)；

m2w——比水容重，為體積含水量和基質(zhì)吸力的偏導數(shù)的負值，=

降雨入滲作用下的非飽和滲流計算中，飽和區(qū)滲透系數(shù)取值為飽和滲透系數(shù)，其值一般通過現(xiàn)場取樣再經(jīng)室內(nèi)或現(xiàn)場實驗獲得；而非飽和滲透系數(shù)的取值則與土體含水率、基質(zhì)吸力存在密切的關系。由于降雨入滲過程實際上是飽和-非飽和入滲過程，降雨過程中路堤中某一區(qū)域的滲透系數(shù)隨基質(zhì)吸力與含水量的變化而時刻變化，本文采用Van Genuchten模型進行擬合，模型表達式如下［12］:

式中:θ，θr，θs——體積含水率、殘余含水率、飽和含水率(cm3/cm3)；

H——負壓(kPa)；

Ks——飽和滲透系數(shù)(m/s)；

α，m，n——土水特征曲線形狀參數(shù)。

2 數(shù)值計算模型及方案確定

2.1 數(shù)值計算模型

在建廣西六寨至河池高速公路K19+133斷面路面寬60m，填方高度為12m。路堤中使用典型的土-石(崩解完成后的炭質(zhì)泥巖)混填粗粒土填料，由于路堤兩側對稱，故選取路堤橫斷面的一半作為研究對象。路堤典型剖面二維有限元網(wǎng)格如圖2所示。計算模型單元數(shù)量為2997，節(jié)點數(shù)量為3114，為了保證計算的精度，通過設置輔助線的方法將網(wǎng)格劃分為四邊形單元。

圖2 路堤典型剖面二維有限元網(wǎng)格Fig.2 Typical profile of the embankment of two-dimensional finite element mesh

滲流邊界條件:模型底部、兩側鉛直位置、路堤路面設置為不透水邊界，路堤坡面、坡腳水平面設置為單位流量自由溢出邊界。當程序運行時，自動判斷降雨強度與土體滲透性的關系。如果降雨強度小于路堤坡面、坡腳水平面表層土體的滲透性，按流量邊界處理，大小為降雨強度；如果雨強大于路堤坡面、坡腳水平面表層土體的滲透性，一部分雨水沿坡面流失，會在坡面形成一薄層水膜，此時可按定水頭邊界處理。由于水膜很薄，計算中取水頭值等于地表高程。

滲流初始條件:根據(jù)路堤坡面及坡腳多個位置鉆孔取樣資料，設置如圖2所示的初始地下水位進行穩(wěn)態(tài)流計算，以此計算結果作為邊坡初始滲流場。

滲流水力學模型采用各向異性達西滲流模型。

滲流計算水力學參數(shù):根據(jù)現(xiàn)場獲取多組原狀土樣進行室內(nèi)飽和滲流試驗，測得壓實度為93%的炭質(zhì)泥巖路堤填料飽和滲透系數(shù)為2.0×10－7～7.0×10－7m/s，路堤內(nèi)部非飽和區(qū)非飽和滲透系數(shù)采用式(2)、(3)進行擬合，當k=4.0×10－7m/s時滲透系數(shù)及含水率隨基質(zhì)吸力變化見圖3。各向異性滲透系數(shù)比Ky/Kx=0.5，飽和體積含水率為0.15。通過調(diào)取當?shù)貧庀筚Y料可知［13］，廣西六寨地區(qū)有歷史記錄以來最大日降雨、最大3日降雨、最長 21日降雨分別為 295.3、357.0、585.8mm，經(jīng)換算為單位降雨強度分別為:q=3.14×10－6m/s、1.37 ×10－6m/s、3.23 ×10－7m/s。

圖3 體積含水率、滲透系數(shù)與基質(zhì)吸力關系曲線圖Fig.3 Volatile water content and permeability coefficient versus matrix suction

2.2 數(shù)值分析方案

為獲得在不同降雨情況下路堤邊坡的滲流狀態(tài)，本文擬根據(jù)實際降雨資料及巖土體室內(nèi)外試驗數(shù)據(jù)對炭質(zhì)泥巖粗粒土路堤邊坡在降雨入滲作用下的滲流過程進行探討，數(shù)值計算模擬方案見表1。

表1 數(shù)值計算模擬方案Table 1 Numerical simulation schemes

3 路堤邊坡暫態(tài)飽和區(qū)變化規(guī)律分析

圖4為計算方案2條件下路堤邊坡在降雨過程中及降雨停止后暫態(tài)飽和區(qū)及地下水位隨時間的變化過程(限于篇幅，未將3個計算方案濕潤鋒變化過程一一列出)。由圖4可知:路堤邊坡隨著降雨的持續(xù)，其表面逐漸會出現(xiàn)一定深度的暫態(tài)飽和區(qū)，暫態(tài)飽和區(qū)面積隨著降雨歷時的增長而逐漸擴大，并與地下水位面連通，致使水位線迅速上升，降雨停止后，飽和區(qū)由路堤邊坡坡面從上往下開始消散。

通過對比路堤各位置暫態(tài)飽和區(qū)濕潤鋒變化規(guī)律可以看出，在降雨歷時相同的情況下，坡面入滲深度比坡腳外側的入滲深度要??；在坡腳附近所形成的暫態(tài)飽和區(qū)最先與地下水位面連通，連通后路堤邊坡坡面以下水位線出現(xiàn)大幅上升，而路堤中部地下水位線則變化較??；降雨停止后，路堤上部暫態(tài)飽和區(qū)沿坡面向下逐漸消散，而其它區(qū)域暫態(tài)飽和區(qū)面積仍出現(xiàn)短時間持續(xù)增長，在一定時間后才開始消散，暫態(tài)飽和區(qū)面積消散相對于降雨停止時間而言具有明顯的滯后特征。其主要原因是由于降雨強度要大于邊坡的入滲量，邊坡表面未入滲的雨水以坡面徑流的形式向坡腳匯集，在匯集過程中，路堤邊坡中、下部入滲量得到坡面徑流雨水的持續(xù)供給，使其暫態(tài)飽和區(qū)面積得以短時間繼續(xù)增大。

圖4 降雨及降雨后暫態(tài)飽和區(qū)及地下水位隨時間的變化Fig.4 Changes in transient saturated area and groundwater level with time during and after rainfall

圖5為3種降雨計算方案暫態(tài)飽和區(qū)面積隨降雨時間的變化規(guī)律，表2為暫態(tài)飽和區(qū)特征值對比。由圖5可知，在降雨階段，暫態(tài)飽和區(qū)的形成及擴展與降雨強度有直接關系。降雨歷時相同，所形成暫態(tài)飽和區(qū)的面積與降雨強度成正比。降雨停止后，由于沒有雨水補給，暫態(tài)飽和區(qū)開始消散，消散速率具有逐漸減緩的趨勢。由表2可知，暫態(tài)飽和區(qū)形成時間以及與地下水位面相連通的時間都與降雨強度密切相關，降雨強度越大所需要的時間越短。由此可以推斷:降雨強度越大，雨水在路堤內(nèi)部擴展的速度亦越大。暫態(tài)飽和區(qū)最大面積則受降雨歷時與降雨強度的共同影響。以方案3為例，在降雨強度最小，降雨歷時最長的條件下，所形成的暫態(tài)飽和區(qū)面積均比方案1、方案2大。

圖5 暫態(tài)飽和區(qū)面積大小隨時間的變化過程Fig.5 Process of changes in the transient saturated area with time

表2 暫態(tài)飽和區(qū)特征值Table 2 Characteristic value of the transient saturated area

4 降雨入滲條件下路堤邊坡穩(wěn)定性分析

4.1 路堤邊坡失穩(wěn)機理分析

降雨入滲是導致邊坡失穩(wěn)的主要原因［5，14］。為進一步揭示降雨條件下路堤邊坡的失穩(wěn)機理，在路堤邊坡表面分別設置A、B、C(圖2)3個特征點，研究路堤在降雨條件下基質(zhì)吸力與體積含水率的變化規(guī)律。圖6為整個降雨過程及降雨停止后一段時間內(nèi)特征點A、B、C的孔隙水壓力與體積含水率變化過程曲線(限于篇幅，只列出方案3的變化過程曲線)。由圖6可知，降雨入滲使路堤內(nèi)部非飽和區(qū)基質(zhì)吸力逐漸降低，最終出現(xiàn)正孔隙水壓力。隨著降雨的持續(xù)，孔隙水壓力有不斷增大的趨勢，在靠近坡腳處孔隙水壓力變化幅度最大。降雨停止后，3個特征點的孔隙水壓力均表現(xiàn)迅速下降，并重新進入非飽和狀態(tài)。同時，路堤在受到降雨入滲的影響后，其體積含水率由初始含水率逐漸變化為飽和含水率，降雨停止一段時間后，隨著暫態(tài)飽和區(qū)的消散，體積含水率也相應降低。

圖6 降雨過程中孔隙水壓力和體積含水率的變化Fig.6 Changes in pore water pressure and volatile water content during rainfall process

一般，巖土體的抗剪強度與含水特征、基質(zhì)吸力大小有密切的關系，降雨條件下飽和-非飽和抗剪強度因雨水入滲引起滲流場變化而時刻變化。因此，確定不同時刻巖土體內(nèi)含水率及基質(zhì)吸力的大小對于進行降雨條件下路堤邊坡穩(wěn)定性的研究顯得至關重要。為此，加拿大學者Fredlund基于摩爾-庫倫公式提出了非飽和巖土體抗剪強度計算方法［15］:

式中:σ——作用于坡體底面上的總法向應力；

c'、φ'——有效粘聚力和有效內(nèi)摩擦角；

φb——因基質(zhì)吸力增加引起抗剪強度增加的曲線傾角；

(σ－ua)——破壞面上的凈應力；

(ua－uw)——破壞面上的基質(zhì)吸力。

從(4)式可以看出，抗剪強度是基質(zhì)吸力的線性函數(shù)。在降雨階段，孔隙水壓力增大引起基質(zhì)吸力最終喪失，同時c'、φ'也相應地減小，從而導致非飽和抗剪強度降低，增大了路堤邊坡發(fā)生失穩(wěn)的可能性。

4.2 降雨過程中路堤邊坡過程穩(wěn)定性分析

(1)穩(wěn)定性計算方法

將式(4)按照ua=0，φb為常數(shù)進行整理得:

式中:c″——非飽和等效粘聚力。

式(5)為本文推導的考慮非飽和抗剪強度理論中基質(zhì)吸力存在的修正Mohr-Coulomb破壞準則。

考慮降雨入滲條件下路堤邊坡穩(wěn)定性計算是在二維飽和-非飽和計算軟件Seep/W及FLAC3D軟件平臺上實現(xiàn)的，首先建立二維有限元模型進行二維飽和-非飽和滲流分析，然后基于二維模型建立三維等效模型，通過在FLAC3D中編寫用戶子程序，將降雨入滲過程中某一時刻的孔隙水壓力導入FLAC3D計算模型。并設置孔隙水壓力大于0的區(qū)域土體重度為γsat，對于孔隙水壓力小于0的區(qū)域則將土體重度在γd與γsat進行線性插值賦值。最終運用強度折減法對該時刻的邊坡穩(wěn)定性進行計算。

(2)計算模型及計算參數(shù)

圖7為計算方案2降雨36h、72h時刻的三維計算結果。應力計算邊界條件采用鉛直邊界上施加水平約束、水平底部邊界上施加鉛直約束。本構模型采用摩爾庫侖彈塑性本構模型。安全系數(shù)計算方法采用強度折減法［16］。計算中所用粘聚力c按式(10)進行設置，土體重度根據(jù)孔隙水壓力的大小在γd與γsat進行線性插值賦值。其余相關物理力學根據(jù)筆者進行的室內(nèi)試驗與業(yè)主工地試驗室提供參數(shù)進行綜合取值(表3)。

圖7 方案2降雨36h、72h時刻的三維計算結果Fig.7 3D model diagram of the embankment after rainfall 36 h and 72 h

表3 路堤邊坡土石混填粗粒土物理力學參數(shù)表Table 3 Physical＆mechanical parameters of rock and soil of the embankment slope

(3)路堤失穩(wěn)模式及安全系數(shù)變化過程分析

圖8為與圖7對應時刻的路堤在計算方案二降雨影響下的塑性區(qū)分布圖。由圖8可知，塑性區(qū)主要分布在路堤邊坡表層，隨著降雨的持續(xù)，邊坡塑性區(qū)逐漸向內(nèi)部延伸并持續(xù)擴大，表現(xiàn)出與暫態(tài)飽和區(qū)發(fā)展規(guī)律一致的擴展規(guī)律。在實際工程中則表現(xiàn)為路堤邊坡在降雨影響下表層土體的局部滑移與整體失穩(wěn)。引起上述塑性區(qū)發(fā)展變化規(guī)律的主要因素有:邊坡表層處于臨空面，其自穩(wěn)能力不強；由于暫態(tài)飽和區(qū)的形成引起土體有效應力減小、非飽和等效粘聚力(c″=(c'－uwtanφb))因為水壓力的存在而降低，最終導致非飽和抗剪強度降低；再者，由于雨水的入滲，邊坡表層及地下水位以下區(qū)域巖土體重度增大，引起邊坡坡面下滑力增大。表4所列為計算方案2路堤邊坡安全系數(shù)與降雨時間的關系，在降雨過程中安全系數(shù)持續(xù)下降，降雨停止后安全系數(shù)緩慢增大。

圖8 降雨36 h和72 h塑性區(qū)分布Fig.8 Distribution of the plastic zone after rainfall 36 h and 72 h

表4 不同降雨時刻路堤穩(wěn)定性系數(shù)(計算方案2)Table 4 Safety coefficient of slope stability in different rainfall times

由此可見，在降雨過程中路堤邊坡穩(wěn)定性持續(xù)下降，降雨停止后，隨著暫態(tài)飽和區(qū)的消散、基質(zhì)吸力的恢復、含水率的降低，路堤邊坡穩(wěn)定性也將得到一定程度的提高。由特征點孔隙水壓力與含水率變化曲線、暫態(tài)飽和區(qū)發(fā)展規(guī)律以及非飽和抗剪強度理論可以推斷:路堤坡腳處抗剪強度降低幅度最大，且在巖土體重度增大的影響下路堤邊坡失穩(wěn)模式為整體滑移與表面局部坍塌相結合，邊坡發(fā)生整體滑移的剪出口位置最有可能出現(xiàn)在坡腳附近。

5 結論

(1)路堤降雨入滲過程實際上是暫態(tài)飽和區(qū)的形成與發(fā)展過程；暫態(tài)飽和區(qū)形成以及與地下水位面連通所需的時間均取決于降雨強度的大小，但擴展范圍大小則受到降雨強度與降雨歷時的共同影響。

(2)降雨過程中，路堤邊坡坡腳附近暫態(tài)飽和區(qū)率先與地下水位面連通，降雨入滲僅會引起路堤邊坡坡面以下水位線的大幅上升，對路堤中部地下水位線影響較小。降雨停止后，路堤上部暫態(tài)飽和區(qū)迅速消散，而路堤中下部暫態(tài)飽和區(qū)面積則繼續(xù)增大一段時間后才逐漸減小，表現(xiàn)出明顯的滯后特點，但總體而言暫態(tài)飽和區(qū)的面積呈下降趨勢。

(3)降雨入滲影響區(qū)域內(nèi)基質(zhì)吸力的喪失與體積含水率的增大，是引起路堤邊坡穩(wěn)定性降低的主要影響因素，路堤邊坡失穩(wěn)模式表現(xiàn)為整體滑移與表面局部坍塌相結合，降雨過程中安全系數(shù)逐漸降低，降雨停止后安全系數(shù)得到恢復。

［1］董云，閻宗嶺.土石混填路基沉降變形特征的二維力學模型試驗研究［J］.巖土工程學報，2007，29(6):943 – 947.［DONG Y，YAN Z L.2D mechanical model tests on settlement of rock-soil filled roadbed［J］. Chinese Journal of Geotechnical Engineering，2007，29(6):943－947.(in Chinese)］

［2］李亮，劉興旺，趙煉恒.降雨入滲對路基穩(wěn)定性影響因素分析［J］.鐵道科學與工程學報，2007，4(2):19－23.［LI L，LIU X W，ZHAO L H.Factors analysis of subgrade stability due to rainfall infiltration［J］.Journal of Railway Science and Engineering，2007，4(2):19 －23.(in Chinese)］

［3］鄒維列.路堤壓實及路基土工復合材料排水層的研究［D］.武漢:武漢大學，2004.［ZOU W L.Compaction of road embankment and subgrade geocomposite drainage layers［D］.Wuhan:Wuhan University，2004.(in Chinese)］

［4］王瑞鋼，閆澍旺，鄧衛(wèi)東.降雨作用下高填土路堤邊坡的滲流穩(wěn)定分析［J］.中國公路學報，2004，17(4):25－30.［WANG R G，YAN S W，DENG W D.Analysis of seepage stability of higll-filled embankment slope due to rainfall infiltration［J］.China Journal of Highway and Transport，2004，17(4):25 － 30.(in Chinese)］

［5］陳曉斌，徐望國，劉小平.降雨入滲對粗粒土路堤變形與穩(wěn)定性的影響［J］.中南大學學報:自然科學版，2011，42(3):766 －771.［CHEN X B，XU W G，LIU X P.Effect of rainfall infiltration on additional settlement and stability debasement of granular soil fillings embankment［J］.Journal of Central South University:Science and Technology，2011，42(3):766－771.(in Chinese)］

［6］李聰.路堤土降雨入滲響應的模型試驗研究［J］.武漢理工大學學報:交通科學與工程版，2010，34(3):550 － 553.［LI C.Experimental Study on the Response of Road Embankment to Rainfall Infiltration［J］.Journal of Wuhan University of Technology:Transportation Science＆ Engineering，2010，34(3):550－553.(in Chinese)］

［7］李汝成，王復明.降雨入滲對泥巖-土混填路堤穩(wěn)定性的研究［J］.巖石力學與工程學報，2008，27(11):2260－2266.［LI R C，WANG F M.Effect of rainfall infiltration on stability of mudstone-soil mixture embankment［J］.Chinese Journal of Rock Mechanics and Engineering，2008，27(11):2260 － 2266.(in Chinese)］

［8］邱路陽，劉毓氚，李大勇.高填方殘積土路堤降雨滑塌機理與治理對策［J］.巖土力學，2007，28(10):2161－2166.［QIU L Y，LIU Y C，LI D Y.A case study of failure mechanism of residual soil filled embankment［J］.Rock and Soil Mechanics，2007，28(10):2161－2166.(in Chinese)］

［9］劉俊新，劉育田，胡啟軍.非飽和地表徑流-滲流和流固體耦合條件下降雨入滲對路堤邊坡穩(wěn)定性研究［J］.巖土力學，2010，31(3):903 －910.［LIU J X，LIU Y T，HU Q J.Stability of embankment slope subjected to rainfall infiltration considering both runoff-underground seepage and fluid-solid coupling［J］.Rock and Soil Mechanics，2010，31(3):903 －910.(in Chinese)］

［10］沈水進，孫紅月，尚岳全，等.降雨作用下路堤邊坡的沖刷-滲透耦合分析［J］.巖石力學與工程學報，2011，30(12):2456 －2462.［SHEN S J，SUN H Y，SHANG Y Q，et al.Scouring-penetration coupling analysis of embankment slopeunder rainfallaction［J］.Chinese Journal of Rock Mechanics and Engineering，2011，30(12):2456 －2462.(in Chinese)］

［11］劉新喜，夏元友，蔡俊杰，等.降雨入滲下強風化軟巖高填方路堤邊坡穩(wěn)定性研究［J］.巖土力學，2007，28(8):1705 － 1709.［LIU X X，XIA Y Y，CAI J J， et al. Study on stability on high-filled embankment slope of highly weathered soft rock under rainfall infiltra-tion［J］.Rock and Soil Mechanics，2007，28(8):1705 －1708.(in Chinese)］

［12］Van Genuchten M T.A closed form equation for prediction the hydraulic conductivity of unsaturated soisl［J］.Soil Sci.Soc.Am.，J.，1980，44:892 －898.

［13］韋力行.廣西壯族自治區(qū)分區(qū)降雨量統(tǒng)計報告［R］.2011.［WEI L X.Guangxi Zhuang Autonomous Region rainfall statistics report［R］.2011.(in Chinese)］

［14］戚國慶，黃潤秋，速寶玉，等.巖質(zhì)邊坡降雨入滲過程的數(shù)值模擬［J］.巖石力學與工程學報，2003，22(4):625－629.［QI G Q，HUANG R Q，SU B Y，et al.Numeric simulation on rainfall infiltration on rock slope［J］.Chinese Journal of Rock Mechanics and Engineering，2003， 22(4):625 － 629. (in Chinese)］

［15］劉建華，查旭東，付宏淵，等.考慮降雨入滲條件下巖質(zhì)邊坡穩(wěn)定性分析［J］.公路交通科技，2009，26(10):33－37.［LIU Ji H，ZHA X D，F(xiàn)U H Y，et al.Analysis of Rock Slope Stability under Rainfall Infiltration Condition［J］.Journal of Highway and Transportation and Research and Development，2009，26(10):33－37.(in Chinese)］

［16］Matsui T，San K C.Finite Element Slope Stability Analysis by Shear Strength Reduction Technique［J］.Soil and foundations，JSSMFE，1992，32(1):59 －70.