蘇永華 李誠誠

摘 ? 要：研究了間歇性強降雨作用下邊坡穩(wěn)定性. 間歇性強降雨下濕潤鋒以上土體將會經(jīng)受干濕循環(huán)作用，考慮干濕循環(huán)對邊坡土體滲透特性及土水特性的影響，改進了傳統(tǒng)Green-Ampt入滲（GA）模型;考慮干濕循環(huán)對土體強度的劣化作用，建立間歇性強降雨下邊坡穩(wěn)定性表達(dá)式.對于不考慮間歇性強降雨作用的邊坡，改進GA模型可簡化為傳統(tǒng)GA模型，說明傳統(tǒng)GA模型是改進GA模型的一個特例;利用改進GA模型和穩(wěn)定性評價方法得到的結(jié)果與工程實際情況基本吻合，證明了該方法的可靠性;根據(jù)降雨入滲速率的變化特征，降雨過程可分為穩(wěn)定階段、持減階段和突變階段等3個階段;傳統(tǒng)GA模型與改進GA模型相比在失穩(wěn)深度和時間上有較大延緩;僅改進穩(wěn)定評價方法的傳統(tǒng)GA模型與改進GA模型相比在失穩(wěn)深度上一致，在失穩(wěn)時間上有延緩.

關(guān)鍵詞：間歇性強降雨;Green-Ampt入滲模型;滲透特性;土水特性;強度衰減;穩(wěn)定性

中圖分類號：P642.22 ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼：A

Slope Stability Analysis Based on Green-Ampt Infiltration

Model under Intermittent Heavy Rainfall

SU Yonghua？覮，LI Chengcheng

（College of Civil Engineering，Hunan University，Changsha 410082，China）

Abstract：The stability of slopes under intermittent heavy rainfall is studied. The soil above the wetting front is subjected to dry-wet cycles under intermittent heavy rainfall. Considering the influence of dry-wet cycles on the soil permeability and soil-water characteristics of slope，the traditional Green-Ampt infiltration （GA） model is improved. The expression of slope stability coefficient under intermittent heavy rainfall is established considering the deterioration of dry-wet cycles on soil strength. For slopes without considering intermittent heavy rainfall，the improved GA model can be simplified to the traditional GA model，which indicates that the traditional GA model is a special case of the improved GA model. The results obtained by using the improved GA model and the stability evaluation method are basically consistent with the actual engineering situation，which proves the reliability of the method. According to the variation characteristics of rainfall infiltration rate，the rainfall process can be divided into three stages： stable stage，continuous reduction stage and mutation stage. The slope infiltration depth and instability time of traditional GA are obviously behind that of the improved GA model. Compared with the improved GA model，the traditional GA model which only improves the stability evaluation method is consistent in the instability depth and delayed in the instability time.

Key words：intermittent heavy rainfall;Green-Ampt infiltration model;permeability;soil-water characteristics;strength deterioration;stability

隨著我國經(jīng)濟的高速發(fā)展，高陡型路塹邊坡日益增多，加之極端天氣多發(fā)，邊坡穩(wěn)定性問題日益突出.強降雨是誘發(fā)邊坡失穩(wěn)的最主要、最常見因素之一. 據(jù)統(tǒng)計，1980年以來，我國大陸所發(fā)生的大型災(zāi)害性滑坡約50%由強降雨直接觸發(fā)，造成大量人員傷亡，經(jīng)濟損失巨大[1]. 因此，建立行之有效的強降雨作用下邊坡失穩(wěn)預(yù)測預(yù)報方法，對防災(zāi)減災(zāi)具有非常重要的理論意義和應(yīng)用價值.

基于降雨入滲模型的邊坡穩(wěn)定分析是評價降雨時邊坡穩(wěn)定性的有效方法.基于毛細(xì)管理論提出的Green-Ampt土壤水分入滲（GA）模型[2]，因其參數(shù)少、物理意義明確、精度高、適用性強、求解簡單而在降雨型滑坡研究中逐漸得到應(yīng)用與發(fā)展. 但傳統(tǒng)GA模型僅適用于均質(zhì)且初始含水率均勻分布的土壤水分入滲過程，因此，國內(nèi)外不少學(xué)者對該模型的適用條件和范圍以及主要參數(shù)的確定等問題進行了大量研究與改進，從而使該模型不斷得到完善和發(fā)展. 常金源等[3]以GA模型為基礎(chǔ)，在考慮動水壓力的基礎(chǔ)上分析了降雨入滲條件下淺層邊坡穩(wěn)定性.石振明等[4]提出適用于多層土邊坡降雨入滲的GA模型，并采用強度折減法分析了不同降雨強度和歷時下邊坡穩(wěn)定. Sung[5]基于初始含水率非均勻分布的GA模型提出了淺層基巖邊坡穩(wěn)定性分析方法. Loáiciga等[6]基于GA模型與徑流波動方程建立了入滲-徑流耦合模型，并探討了該模型下的邊坡穩(wěn)定性. Yao等[7]提出了GA模型和分層假定下的邊坡穩(wěn)定性分析方法.

由于間歇性強降雨將導(dǎo)致坡內(nèi)土體處于飽和-非飽和干濕循環(huán)過程，因此有必要用發(fā)展動態(tài)的觀點研究間歇性降雨對邊坡穩(wěn)定性的影響.然而目前GA模型應(yīng)用于邊坡穩(wěn)定性研究時主要考慮持續(xù)性降雨下邊坡穩(wěn)定性，鮮有利用GA模型研究間歇性降雨下邊坡穩(wěn)定性問題的報道，這限制了GA模型在邊坡穩(wěn)定性分析中的進一步應(yīng)用.

本文首先提出一種適用于間歇性強降雨的入滲模型，該模型以 GA 模型為基礎(chǔ)，考慮了間歇性強降雨引發(fā)的干濕循環(huán)效應(yīng)對飽和層土體滲透特性及土水特性的影響;然后考慮干濕循環(huán)對土體強度的衰減作用，利用雙強度折減法分析邊坡穩(wěn)定性;最后分析了本文模型與GA模型的關(guān)系，并通過與實際工程情況對比，驗證了本文計算方法的可靠性.

1 ? 降雨入滲模型

Mein等[8]將GA 模型應(yīng)用到降雨入滲情況，認(rèn)為降雨入滲過程應(yīng)分為降雨強度控制階段和土體入滲能力控制階段.如圖1所示，假設(shè)強降雨強度為q，飽和體積含水率為θs，初始體積含水率為θi，邊坡傾角為β. 降雨初期，坡表土體處于非飽和狀態(tài)，土體入滲能力大于降雨強度，降雨全部滲入土體，此時，垂直于坡面方向的降雨入滲速率為：

以每次強降雨濕潤鋒所達(dá)到的最大深度邊界線，對土體進行分層.記第1次強降雨時濕潤鋒所達(dá)到的最大深度為hs，m1，第2次強降雨時濕潤鋒所達(dá)到的最大深度為hs，m2，第N次強降雨時濕潤鋒所達(dá)到的最大深度為hs，m N，則土層1的深度范圍為h∈[0，hs，m1]，土層2的深度范圍為h∈（hs，m1，hs，m2]，土層N的深度范圍為h∈（h s，m（N-1），hs，m N]，其余土層為初始土層，其具體分層如圖3所示. 由于干濕循環(huán)由強降雨引起，故第N次強降雨時，地表土體滲透特性參數(shù)與土水特性參數(shù)采用土體經(jīng)歷的n次干濕循環(huán)之后相應(yīng)的參數(shù)，對土層1：n = N - 1，對于土層2：n = N - 2，其余依次類推，且n ≥ 0.當(dāng)n = 0時，表示該層土體為未受降雨影響的初始土層，即此次降雨前雨水入滲深度未到達(dá)該土層.

假設(shè)降雨強度控制階段轉(zhuǎn)為土體入滲能力控制階段發(fā)生在h < hs，m1范圍內(nèi)，令hs，N、hsp，N和tp，N分別為第N次強降雨時濕潤鋒深度、臨界時刻濕潤鋒深度和臨界時刻，故第N次強降雨時土體入滲能力控制階段的降雨入滲速率式（2）可修改為：

式（26）反映了某次強降雨下濕潤鋒深度隨降雨歷時增長的動態(tài)變化規(guī)律，考慮了間接性強降雨對邊坡土體滲透特性及土水特性的影響，推廣了GA模型的應(yīng)用范圍.

2 ? 邊坡穩(wěn)定性分析

目前邊坡穩(wěn)定性分析方法主要有極限平衡法[12]和強度折減法[13]. 降雨入滲條件下，非飽和土邊坡多發(fā)生平行于邊坡表面的淺層滑坡，破壞面往往發(fā)生在濕潤鋒或相對隔水層處.本文以無限長邊坡體為研究對象，假定濕潤鋒以上土體處于飽和狀態(tài)，并以濕潤鋒處作為潛在滑動面，在考慮干濕循環(huán)對土體力學(xué)特性的衰減作用下，利用雙強度折減法評價邊坡穩(wěn)定性隨強降雨濕潤鋒深度、強降雨次數(shù)的動態(tài)變化過程，其計算簡圖如圖4所示.

3 ? 模型驗證與分析

3.1 ? 與傳統(tǒng)GA模型解析解比較

傳統(tǒng)GA模型不考慮間歇性降雨對邊坡土體滲透特性及土水特性的影響，即令改進GA模型強降雨次數(shù)N = 1，則式（26）可化簡為：

3.2 ? 案例計算與分析

廣西百隆路某段右側(cè)高10 m、坡率1 ∶ 1.5的新修路塹邊坡在經(jīng)歷第1個雨季時發(fā)生淺層坍滑.土體設(shè)計參數(shù)見表1.土體土水特征曲線函數(shù)表達(dá)式、有效黏聚力c′、有效內(nèi)摩擦角φ′和孔隙比與干濕循環(huán)次數(shù)n滿足的函數(shù)關(guān)系分別如式（35）～（39）所示.假設(shè)當(dāng)?shù)仄骄鶜鉁?5 ℃，每次強降雨歷時為4 h，降雨強度 q = 0.000 8 cm/s.

不同干濕循環(huán)次數(shù)下，土體脫濕曲線和吸濕曲線的偏移系數(shù)見表2.圖5所示為n = 0和n = 10土水特征曲線.同文獻(xiàn)[11]所反映的規(guī)律相同，該土體不同干濕循環(huán)下土水特征曲線均存在明顯的滯后現(xiàn)象，且隨著干濕循環(huán)次數(shù)的增加，滯回效應(yīng)明顯減小，并趨于穩(wěn)定;若干次干濕循環(huán)后土體土水特征曲線中脫濕曲線的偏移量大于吸濕曲線的偏移量.

圖6所示為不同干濕循環(huán)次數(shù)下，土體吸濕過程濕潤鋒處基質(zhì)吸力與土層飽和導(dǎo)水率變化曲線.可知，土體吸濕時濕潤鋒處基質(zhì)吸力隨著干濕循環(huán)次數(shù)的增加而減小，且首次循環(huán)降幅最大，多次循環(huán)后趨于穩(wěn)定.飽和導(dǎo)水率隨著干濕循環(huán)次數(shù)的增加而增大，但增大速率逐漸放緩.土體吸濕時濕潤鋒處基質(zhì)吸力和飽和導(dǎo)水率的變化皆表明了隨著干濕循環(huán)次數(shù)的增加，土體逐漸松散，土體持水能力逐漸降低，符合相關(guān)試驗研究和自然條件下土體實際滲流特性[18].

根據(jù)降雨入滲速率的變化特征，將降雨過程分為穩(wěn)定階段、持減階段和突變階段等3個階段，如圖7所示. 但當(dāng)N = 1時，只有穩(wěn)定階段、持減階段2個階段;不同降雨次數(shù)下每個階段開始時間和持續(xù)時間也并不相同.

階段的入滲速率與強降雨次數(shù)無關(guān)，保持為qcos β.某次降雨下，該階段濕潤鋒深度與降雨歷時呈正比，邊坡穩(wěn)定性系數(shù)隨降雨歷時的增長快速下降.不同強降雨次數(shù)下，該階段持續(xù)時間隨強降雨次數(shù)的增加而增加，其原因在于干濕循環(huán)過程中，飽和導(dǎo)水率的變化占主導(dǎo)位置，淺層土體入滲能力隨干濕循環(huán)次數(shù)的增加而增強;邊坡穩(wěn)定性系數(shù)隨強降雨次數(shù)的增加而減小，其原因在于干濕循環(huán)降低了土體的抗剪強度.

持減階段是指降雨入滲速率隨降雨歷時的增長連續(xù)穩(wěn)定下降階段.假定第一次降雨時其濕潤鋒所達(dá)到的最大深度為hs，m1，某次降雨下（N > 1），濕潤鋒到達(dá)hs，m1深度的降雨歷時為ts，m N，該階段的持續(xù)時間段為tp，N ～ ts，m N;濕潤鋒深度隨降雨歷時的增加而增加，但增加速度逐漸變慢;邊坡穩(wěn)定性系數(shù)隨降雨歷時的增長而下降，但下降速度亦逐漸放緩.不同強降雨次數(shù)下，該階段的持續(xù)時間隨強降雨次數(shù)的增加而減小，其原因在于干濕循環(huán)效應(yīng)使得土層飽和導(dǎo)水率不斷增大，而基質(zhì)吸力變化并不明顯，故上層土層內(nèi)的滲透速率增大，降雨穿越上層土層的時間變短，而且隨著強降雨次數(shù)的增加，淺層土體會更多地按雨強控制階段的入滲速率進行滲透;該階段的入滲速率、濕潤鋒深度隨強降雨次數(shù)的增加而增加;邊坡穩(wěn)定性系數(shù)隨強降雨次數(shù)的增加而減小.

突變階段是指ts，m N至降雨結(jié)束階段.該階段的入滲速率呈階梯式突變，其原因在于不同土層土體所經(jīng)歷的干濕循環(huán)次數(shù)并不相同，土體埋深越深其經(jīng)歷的干濕循環(huán)次數(shù)就越少，其飽和導(dǎo)水率小于上層土體，基質(zhì)吸力大于上層土體，但本案例中飽和導(dǎo)水率的變化占主導(dǎo)地位，故其入滲速率會在土層接觸面發(fā)生向下突變.某次降雨下，突變點之間的入滲速率逐漸降低;濕潤鋒深度隨降雨歷時增加而增加，但增加速率減緩較快;邊坡穩(wěn)定性系數(shù)隨降雨歷時的增大而減小，但在土層接觸面發(fā)生突變而變大，其原因在于深部土體所受干濕循環(huán)較少，土體抗剪強度隨深度逐漸提高.不同強降雨次數(shù)下，入滲速率的降低速率隨強降雨次數(shù)的增加而增大，該階段結(jié)束時刻，降雨次數(shù)越多，入滲速率越小，其原因在于強降雨下濕潤鋒最終達(dá)到的深度會逐漸增加，且其土層為原狀土體，其飽和導(dǎo)水率與基質(zhì)吸力相同，但深度加深;濕潤鋒深度隨強降雨次數(shù)的增加而增加;邊坡穩(wěn)定性系數(shù)隨強降雨次數(shù)的增加而減小.

如圖11所示，本文模型邊坡失穩(wěn)時，降雨次數(shù)N=5，降雨歷時t=5 956 s，滑動面埋深hs=70.74 cm，說明邊坡在前4次強降雨下能保持穩(wěn)定，但間歇性降雨引發(fā)的干濕循環(huán)效應(yīng)導(dǎo)致土體土水保持能力和強度降低，穩(wěn)定性系數(shù)越來越小，第5次強降雨時發(fā)生淺層滑坡，這與廣西百隆路修建過程中的邊坡實際破壞情況基本吻合[19-20]，證實本文方法的可靠性.

一般而言，邊坡穩(wěn)定性系數(shù)不會隨著降雨歷時的增加而增加，故本文方法下邊坡穩(wěn)定性系數(shù)應(yīng)取整個降雨過程穩(wěn)定性系數(shù)的最小值. 由圖12可知，不同強降雨次數(shù)下，邊坡穩(wěn)定性系數(shù)的最小值皆出現(xiàn)在第二土層內(nèi)，解釋了降雨條件下非飽和土質(zhì)邊坡多發(fā)生淺層破壞的原因.本文方法同時考慮了干濕循環(huán)作用下邊坡土體風(fēng)化層逐步加深和土體強度的逐步衰減，說明了本文方法的全面性.同時，由以上分析可知，初始狀態(tài)下均質(zhì)土邊坡在經(jīng)歷若干次降雨后將變成多層非均質(zhì)土坡，邊坡土體滲透特性參數(shù)、土水特征參數(shù)及強度參數(shù)非均勻分布，且表現(xiàn)為沿深度方向發(fā)生變化，此變化過程及變化結(jié)果符合實際情況，從側(cè)面證明了本文所提模型的合理性.

3.3 ? 不同模型對比

仍以上述案例為例，分別以傳統(tǒng)GA模型+極限平衡法、傳統(tǒng)GA模型+雙強度折減法、改進GA模型+雙強度折減法3種模型計算邊坡失穩(wěn)時間及失穩(wěn)深度.每次計算中降雨強度相同，傳統(tǒng)GA模型+極限平衡法的降雨歷時為5 × 104 s，其余模型每次降雨歷時為5 × 104 s除以對應(yīng)降雨次數(shù). 其計算結(jié)果如表3所示，其中計算模型一、二、三分別與傳統(tǒng)GA模型+極限平衡法、傳統(tǒng)GA模型+雙強度折減法和改進GA模型+雙強度折減法相對應(yīng)，表中計算結(jié)果（1，400 093，400）含義：1表示邊坡失穩(wěn)時的降雨次數(shù)N=1，400 093表示失穩(wěn)時的降雨歷時t=400 093 s，400表示失穩(wěn)時濕潤鋒深度h = 400 cm.

由表3可知，相同降雨量下，持續(xù)性降雨誘發(fā)邊坡失穩(wěn)的時間和滑坡體量遠(yuǎn)大于間歇性降雨的結(jié)果，這在某種程度上證實了“大雨大滑，小雨小滑”. 僅利用GA模型結(jié)合極限平衡法評價強降雨下邊坡的穩(wěn)定性，而忽略前期間歇性強降雨的影響會使預(yù)測結(jié)果在時間和深度上有較大延遲，進而導(dǎo)致邊坡失穩(wěn)得不到及時有效的預(yù)測，并且在防治時造成不必要的資源浪費.因此降雨入滲時，有必要考慮間歇性強降雨作用，改進模型更有利于準(zhǔn)確預(yù)測入滲深度與邊坡穩(wěn)定性分析.模型二與模型三相比在失穩(wěn)的深度、體量上基本相同，但模型二在失穩(wěn)時間上有較大延遲，兩種模型所預(yù)測邊坡失穩(wěn)的時間和深度皆隨降雨次數(shù)的增加而減小.模型三考慮了間歇性降雨對土體滲透特性、土水特性和土體強度的影響，相比于前人研究成果更加符合實際.

4 ? 討 ? 論

1）本文所提模型的入滲速率和穩(wěn)定性系數(shù)存在突變階段的根本原因在于GA模型假定濕潤鋒面至坡面間的土體體積含水率為飽和含水率，間歇性強降雨作用下濕潤鋒面面上面下土體所經(jīng)歷的干濕循環(huán)次數(shù)不同，故在每次強降雨濕潤鋒所達(dá)到的最大深度附近土體在滲透特性和土水特性有較大不同，進而產(chǎn)生突變.事實上，濕潤鋒面至坡面間土體并非完全飽和，而是存在飽和層、非飽和過渡區(qū)、傳導(dǎo)層等[21].間歇性降雨下基于分層假定的邊坡穩(wěn)定性分析有待進一步研究.

2）土體入滲能力隨強降雨次數(shù)增長而變化的過程中，飽和導(dǎo)水率的變化或基質(zhì)吸力的變化占主導(dǎo)地位將對降雨入滲速率、濕潤鋒深度等隨降雨歷時、干濕循環(huán)次數(shù)變化而變化的物理量產(chǎn)生不同的影響，限于篇幅本文不再分析，具體影響可結(jié)合工程實例具體分析，分析方法同本文所提案例分析方法相同.

3）本文對前人優(yōu)秀研究成果加以綜合利用，以

提供一種可用于間歇性強降雨下邊坡降雨入滲計算及邊坡穩(wěn)定性分析方法.為使論文在表達(dá)上簡潔有效，在公式推導(dǎo)過程中用函數(shù)關(guān)系表達(dá)孔隙比、土體強度參數(shù)隨干濕循環(huán)次數(shù)變化規(guī)律.這種函數(shù)關(guān)系是真實存在的，只是函數(shù)表達(dá)式不一定是連續(xù)函數(shù)，在實際工程或試驗中此函數(shù)更可能符合離散函數(shù)，利用本文方法進行間歇性強降雨入滲及邊坡穩(wěn)定性分析時，可利用工程或試驗實測結(jié)果進行計算.

5 ? 結(jié) ? 論

1）基于傳統(tǒng)GA模型，本文推導(dǎo)了間歇性強降雨下考慮土體滲透特性及土水特性變化的退化GA模型，根據(jù)降雨入滲速率的特征，將降雨歷程劃分為穩(wěn)定階段、持減階段、突變階段，并有機結(jié)合極限平衡法和強度折減法進行邊坡穩(wěn)定性分析.

2）某次降雨下，濕潤鋒深度隨降雨歷時的增長持續(xù)向下發(fā)展，但發(fā)展速率逐漸放緩;同一降雨歷時下，濕潤鋒深度隨著強降雨次數(shù)的增加而增長，增長速率亦逐漸放緩.

3）某次降雨下隨降雨歷時的增加，穩(wěn)定階段、持減階段邊坡穩(wěn)定性系數(shù)逐漸減小，突變階段邊坡穩(wěn)定性系數(shù)發(fā)生向上突變而增大;隨著強降雨次數(shù)的增加，邊坡穩(wěn)定性系數(shù)逐漸減小.

4）利用傳統(tǒng)GA模型結(jié)合極限平衡法評價邊坡穩(wěn)定性時在失穩(wěn)深度和時間上有較大延緩;傳統(tǒng)GA模型結(jié)合雙強度折減法評價邊坡穩(wěn)定性時能很好地反映失穩(wěn)深度，但在失穩(wěn)時間上有較大延遲;改進GA模型結(jié)合雙強度折減法，既考慮了干濕循環(huán)對土體滲透特性和土水特性的影響，又考慮了其對土體抗剪強度的衰減作用，更加符合實際.

參考文獻(xiàn)

[1] ? ?黃潤秋. 20世紀(jì)以來中國的大型滑坡及其發(fā)生機制[J].巖石力學(xué)與工程學(xué)報，2007，26（3）：433—454.

HUANG R Q. Large-scale landslides and their sliding mechanisms in China since the 20th century[J]. Chinese Journal of Rock Mechanics and Engineering，2007，26（3）：433—454. （In Chinese）

[2] ? ?HEBER G W，AMPT G A . Studies on soil physics [J]. The Journal of Agricultural Science，1911，4（1）：1—24.

[3] ? ?常金源，包含，伍法權(quán)，等. 降雨條件下淺層滑坡穩(wěn)定性探討[J]. 巖土力學(xué)，2015，36（4）：995—1001.

CHANG J Y，BAO H，WU F Q，et al. Discussion on stability of shallow landslide under rainfall [J]. Rock and Soil Mechanics，2015，36（4）：995—1001. （In Chinese）

[4] ? ?石振明，沈丹祎，彭銘，等. 考慮多層非飽和土降雨入滲的邊坡穩(wěn)定性分析[J]. 水利學(xué)報，2016，47（8）：977—985.

SHI Z M，SHEN D Y，PENG M，et al. Slope stability analysis by considering rainfall infiltration in multi-layered unsaturated soils [J]. Journal of Hydraulic Engineering，2016，47（8）：977—985. （In Chinese）

[5] ? ?SUNG E C. Prediction of shallow landslide by surficial stability analysis considering rainfall infiltration [J]. Engineering Geology，2017，231：126—138.

[6] ? ?LO？魣ICIGA H A，JOHNSON J ?M. Infiltration on sloping terrain and its role on runoff generation and slope stability [J]. Journal of Hydrology，2018，561：584—597.

[7] ? ?YAO W M，LI C D，ZHAN H B，et al. Time-dependent slope stability during intense rainfall with stratified soil water content [J]. Bulletin of Engineering Geology and the Environment，2019，2：1—15.

[8] ? ?MEIN R G，LARSON C L . Modeling infiltration during a steady rain [J]. Water Resources Research，1973，9（2）：384—394.

[9] ? ?江強強，劉路路，焦玉勇，等. 干濕循環(huán)下滑帶土強度特性與微觀結(jié)構(gòu)[J]. 巖土力學(xué)，2019，40（3）：1005—1012+1022.

JIANG Q Q，LIU L L，JIAO Y Y，et al. Strength properties and microstructure characteristics of slip zone soil under wetting-drying cycles [J]. Rock and Soil Mechanics，2019，40（3）：1005—1012+1022. （In Chinese）

[10] ?劉昌鑫，潘健，鄧羽松，等. 干濕循環(huán)對崩崗?fù)馏w穩(wěn)定性的影響 [J]. 水土保持學(xué)報，2016，30（6）：253—258.

LIU C X，PAN J，DENG Y S，et al. The influence of wet-dry cycles on collapsing gully soil stability [J]. Bulletin of Soil and Water Conservation，2016，30（6）：253—258. （In Chinese）

[11] ?張俊然，許強，孫德安. 多次干濕循環(huán)后土-水特征曲線的模擬[J]. 巖土力學(xué)，2014，35（3）：689—695.

ZHANG J R，XU Q，SUN D A. Simulation of soil-water characteristic curves during drying and wetting cycles [J]. Rock and Soil Mechanics，2014，35（3）：689—695. （In Chinese）

[12] ?陳昌富，楊宇. 邊坡穩(wěn)定性分析水平條分法及其進化計算[J].湖南大學(xué)學(xué)報（自然科學(xué)版），2004，31（3）：72—75.

CHEN C F，YANG Y. Horizontal slice method of slope stability analysis based on evolutionary algorithm [J]. Journal of Hunan University （Natural Sciences），2004，31（3）：72—75. （In Chinese）

[13] ?蘇永華，張航，蹇宜霖. 基于滑帶剛度退化的古滑坡自穩(wěn)能力研究[J]. 湖南大學(xué)學(xué)報（自然科學(xué)版），2018，45（7）：77—84.

SU Y H，ZHANG H，JIAN Y L. Research on self stabilizing ability of ancient landslide based on stiffness degradation of sliding zone [J]. Journal of Hunan University （Natural Sciences），2018，45（7）：77—84. （In Chinese）

[14] ?涂義亮，劉新榮，鐘祖良，等. 干濕循環(huán)下粉質(zhì)黏土強度及變形特性試驗研究[J].巖土力學(xué)，2017，38（12）：3581—3589.

TU Y L，LIU X R，ZHONG Z L，et al. Experimental study on strength and deformation characteristics of silty clay during wetting-drying cycles [J]. Rock and Soil Mechanics，2017，38（12）：3581—3589.（In Chinese）

[15] ?CHEN L ，YOUNG M H . Green-Ampt infiltration model for sloping ? ?surfaces [J]. Water Resources Research，2006，42（7）：1—9.

[16] ?黨發(fā)寧，劉海偉，王學(xué)武，等. 基于有效孔隙比的黏性土滲透系數(shù)經(jīng)驗公式研究[J]. 巖石力學(xué)與工程學(xué)報，2015，34（9）：1909—1917.

DANG F N，LIU H W，WANG X W，et al. Empirical formulas of permeability of clay based on effective pore ratio [J]. Chinese Journal of Rock Mechanics and Engineering，2015，34（9）：1909—1917. （In Chinese）

[17] ?SL237—1999 ? 土工試驗規(guī)程[S].北京：中國水利水電出版社，1999：68—69.

SL237—1999 ? Specification of soil test [S]. Beijing：China Water Power Press，1999：68—69. （In Chinese）

[18] ?孫德安，黃丁俊. 干濕循環(huán)下南陽膨脹土的土水和變形特性[J].巖土力學(xué)，2015，36（S1）：115—119.

SUN D A，HUANG D J. Soil-water and deformation characteristics of Nanyang expansive soil after wetting-drying cycles [J]. Rock and Soil Mechanics，2015，36（S1）：115—119. （In Chinese）

[19] ?楊和平，肖杰，程斌，等. 開挖膨脹土邊坡坍滑的演化規(guī)律[J]. 公路交通科技，2013，30（7）：18—24.

YANG H P，XIAO J，CHENG B，et al. Landslide evolution regularity of excavated expansive soil slope [J]. Journal of Highway and Transportation Research and Development，2013，30（7）：18—24. （In Chinese）

[20] ?王興正.模擬干濕循環(huán)及低應(yīng)力條件的膨脹土抗剪強度試驗研究[D]. 長沙：長沙理工大學(xué)交通運輸工程學(xué)院，2014：45—47.

WANG X Z. The research of shear strength tests of expansive soil with dry-wet cycle simulation and lower stress condition [D].Changsha：School of Traffic & Transportation Engineering，Changsha University of Science & Technology，2014：45—47. （In Chinese）

[21] ?BODMAN G B，COLMAN E A. Moisture and energy conditions during downward entry of water into soils[J]. Soil Science Society of America Journal，1944，8（2）：166—182.