降雨及滲流條件下層狀砂質(zhì)板巖邊坡變形破壞模式研究

毋遠(yuǎn)召 馬文禮 魏占璽 董順德

摘 要：為增強(qiáng)邊坡在惡劣條件下的穩(wěn)定性，針對(duì)降雨及滲流條件下的層狀砂質(zhì)板巖邊坡進(jìn)行研究，采用邊坡模型試驗(yàn)和有限元模擬方法，研究了邊坡的變形破壞形態(tài)和位移變化情況。結(jié)果表明，邊坡的變形破壞是一個(gè)持續(xù)發(fā)展的過(guò)程。主要從坡體中部開(kāi)始，形成一個(gè)D字狀的位移曲線，隨著降雨時(shí)間的增加，最終坡體完全被破壞;當(dāng)降雨強(qiáng)度分別為600 mm/d和800 mm/d時(shí)，邊坡中部出現(xiàn)塌落，其降雨時(shí)間分別為16 800 s和13 300 s，說(shuō)明了強(qiáng)降雨對(duì)邊坡失穩(wěn)有被更大的加速作用;通過(guò)對(duì)邊坡模型試驗(yàn)和有限元模擬試驗(yàn)的結(jié)果對(duì)比表明，層狀砂質(zhì)板巖邊坡破壞形態(tài)和位移變化趨勢(shì)基本保持一致，驗(yàn)證了模型試驗(yàn)監(jiān)測(cè)數(shù)據(jù)的準(zhǔn)確性和有效性。研究結(jié)果為降雨及滲流條件下層狀砂質(zhì)板巖邊坡變形破壞機(jī)理提供了實(shí)驗(yàn)依據(jù)和數(shù)值模擬參考，可為類似砂質(zhì)板巖邊坡滲流問(wèn)題研究提供借鑒。

關(guān)鍵詞：巖土力學(xué);降雨;邊坡;穩(wěn)定性;有限元模擬

中圖分類號(hào)：P642?文獻(xiàn)標(biāo)識(shí)碼：A

doi： 10.7535/hbgykj.2020yx04004

文章編號(hào)：1008-1534（2020）04-0230-07

Abstract：In order to enhance the stability of side slopes in harsh conditions， the research on layered sandy slate slopes under rainfall and seepage conditions was carried on， and the slope model test and finite element simulation method were used to study the deformation failure form and displacement change of the slope. The results show that the deformation and failure is a continuous development process. It mainly starts from the middle of the slope and forms a D-shaped displacement curve. With the increase of rainfall time， the slope is completely destroyed. When the rainfall intensity is 600 mm/d and 800 mm/d， the rainfall time required for the collapse in the middle of the slope is 16 800 s and 13 300 s respectively， which shows that the continuous rainfall has a greater acceleration effect on the slope instability. The comparative analysis of the slope model test and the finite element simulation shows that the failure pattern and displacement change trend of the layered sandy slate slope are basically consistent， which also verifies the accuracy and effectiveness of the model test monitoring data. The research results not only provide experimental basis and numerical simulation foundation for the deformation failure mechanism of underlying sandy slate slope under the condition of rainfall and seepage， but also provide reliable reference for similar research on seepage problems of sandy slate slope.

Keywords：geotechnical mechanics; rainfall; side slope; stability; finite element simulation

在公路、隧道、鐵路城建等基建項(xiàng)目的建設(shè)過(guò)程中會(huì)遇到各種各樣的地質(zhì)條件[1]。

為了保證道路及建設(shè)工程的安全性，需要對(duì)各種惡劣條件下的邊坡性能進(jìn)行研究。在持續(xù)降雨的條件下，邊坡會(huì)出現(xiàn)失穩(wěn)的現(xiàn)象，輕則影響工程的質(zhì)量和施工進(jìn)度，重則導(dǎo)致人員的傷亡[2-3]。

針對(duì)邊坡的變形破壞有很多研究[4-6]。楊世平等[7]結(jié)合歷史資料，利用數(shù)值模擬的方法研究了水庫(kù)右岸趾板邊坡的變形機(jī)理;何忠明等[8]結(jié)合非飽和抗剪強(qiáng)度理論，分析了強(qiáng)降雨條件下邊坡的安全系數(shù);劉杰等[9]基于有限元計(jì)算理論對(duì)影響邊坡降雨入滲深度的因素進(jìn)行了分析。孫萍等[10]模擬降雨系統(tǒng)，研究了不同降雨條件下天然黃土邊坡的變形破壞模式;王磊等[11]針對(duì)均質(zhì)黃土邊坡開(kāi)展了人工降雨現(xiàn)場(chǎng)試驗(yàn)，探究了邊坡變形開(kāi)裂的特征;AMINI等[12]針對(duì)巖石邊坡的傾覆失效進(jìn)行了研究，并建立了具有撓曲傾覆失穩(wěn)潛力的巖石邊坡理論模型;PARK等[13]通過(guò)單軸壓縮試驗(yàn)對(duì)巖石的微觀結(jié)構(gòu)進(jìn)行數(shù)值模擬分析，并驗(yàn)證了巖石的力學(xué)性能;LI等[14]研究了運(yùn)動(dòng)學(xué)在滑坡之間相互作用的影響機(jī)理。

降雨對(duì)于邊坡有很大影響，大部分滑坡現(xiàn)象的產(chǎn)生緣于長(zhǎng)時(shí)間降雨。針對(duì)工程實(shí)際中降雨條件下的層狀砂質(zhì)板巖邊坡的模型試驗(yàn)尚未有文獻(xiàn)報(bào)道。筆者結(jié)合邊坡模型和有限元模擬軟件，分析降雨及滲流條件下層狀砂質(zhì)板巖邊坡的變形破壞模式，以期能夠?yàn)樘岣哌吰碌姆€(wěn)定性能提供解決方案。

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

以層狀砂質(zhì)板巖邊坡為分析對(duì)象，在降雨和滲流條件下，巖體會(huì)承受以下集中力的作用：流體對(duì)邊坡裂隙壁的拉力T、流體的阻力F、裂隙內(nèi)的靜水壓力V、邊坡坡體上的靜水壓力U以及邊坡的自重G[15-16]。邊坡的受力分析如圖1所示。

在降雨和滲流條件下，此時(shí)流體對(duì)邊坡裂隙的推力T增加，同時(shí)黏聚力c減小，使得邊坡容易失穩(wěn)[17]。

2?邊坡模型試驗(yàn)設(shè)計(jì)

選擇邊坡模型試驗(yàn)其材料需要滿足以下幾個(gè)方面的要求：1）模型的物理力學(xué)性能需要與邊坡的巖石結(jié)構(gòu)相似;2）試驗(yàn)過(guò)程中需要保持模型的性能穩(wěn)定，不能被空氣中的溫度和濕度影響其力學(xué)性能;3）模型制作簡(jiǎn)單，成本低，材料易得[18-19]。

筆者選擇某高速公路其中一段層狀砂質(zhì)板巖邊坡作為研究對(duì)象，其原狀巖土的物理性能如表1所示。

為了能夠簡(jiǎn)單、準(zhǔn)確的模擬高速路層狀砂質(zhì)板巖邊坡的性質(zhì)，邊坡模型試驗(yàn)簡(jiǎn)化后的物理力學(xué)性能如表2所示。

為了能夠更好地了解降雨及滲流條件下邊坡的變形破壞情況，設(shè)置人工降雨模擬系統(tǒng)。降雨采用人工控制調(diào)節(jié)模式，降雨的有效面積為1.5 m×2.5 m，雨量的高度設(shè)置為3 m，持續(xù)時(shí)間可以手動(dòng)控制，降雨的強(qiáng)度根據(jù)當(dāng)?shù)赝杲涤炅窟M(jìn)行模擬。降雨系統(tǒng)的主要框架如圖2所示。

此降雨系統(tǒng)主要由供水系統(tǒng)、降雨管路和控水閥門(mén)組成，人工模擬降雨的情況。模型箱上部的降雨噴頭如圖3所示。

本研究的邊坡模型砌筑過(guò)程如下：1）選擇合適的模型箱，鋪設(shè)一定角度的墊層;2）待墊層養(yǎng)護(hù)達(dá)到一定的強(qiáng)度后，選擇合適的砌筑方式砌筑邊坡;3）養(yǎng)護(hù)邊坡模型，直至其達(dá)到一定的強(qiáng)度，再在對(duì)應(yīng)的位置進(jìn)行打孔、安裝。邊坡模型的成型過(guò)程如圖4所示。

3?邊坡模型監(jiān)測(cè)數(shù)據(jù)分析

本研究模擬邊坡坡面傾角25°的層狀砂質(zhì)板巖邊坡的破壞情況，采用層鋪法砌筑邊坡模型，在邊坡的上、中、下部鋪設(shè)相應(yīng)的位移傳感監(jiān)測(cè)系統(tǒng)。邊坡模型試驗(yàn)前需要對(duì)降雨系統(tǒng)進(jìn)行檢測(cè)，確認(rèn)邊坡模型能正常工作，監(jiān)測(cè)模擬降雨強(qiáng)度分別為600，?800 mm/d?時(shí)的數(shù)據(jù)的變化情況。

邊坡坡面傾角25°層狀砂質(zhì)板巖邊坡在降雨強(qiáng)度為600 mm/d時(shí)的邊坡模型水平、豎向位移圖如圖5和圖6所示。

從圖5和圖6可以看出。

1）在降雨初期，邊坡模型的水平和豎向位移沒(méi)有變化，說(shuō)明降雨初期邊坡處于穩(wěn)定的狀態(tài)。隨著降雨時(shí)間的逐漸增加，邊坡的含水率也逐漸增加。當(dāng)坡體某一區(qū)域的含水率達(dá)到飽和狀態(tài)時(shí)，則在滲流的作用下邊坡會(huì)出現(xiàn)明顯的位移。

2）當(dāng)降雨時(shí)間達(dá)到一定時(shí)長(zhǎng)時(shí)，邊坡會(huì)出現(xiàn)滲流現(xiàn)象。此時(shí)邊坡的位移會(huì)進(jìn)一步增加，中部的水平、豎向位移出現(xiàn)得更快，最終形成一個(gè)D字狀的位移曲線。

3）隨著降雨時(shí)間的繼續(xù)增加，邊坡大范圍地出現(xiàn)含水率飽和的現(xiàn)象，此時(shí)的邊坡位移變化情況繼續(xù)增加。當(dāng)降雨時(shí)間為16 800 s時(shí)，邊坡中部位移最大，此時(shí)坡體中部的標(biāo)志牌滑落，邊坡破壞。

4）整個(gè)降雨過(guò)程中，邊坡的下部位移基本保持不變。當(dāng)邊坡中部破壞后，隨著降雨的持續(xù)，上部的邊坡也慢慢的塌落，最終整體邊坡垮塌。

邊坡坡面傾角為25°層狀砂質(zhì)板巖邊坡在降雨強(qiáng)度為800 mm/d時(shí)的位移變形情況如圖7和圖8所示。

從圖7和圖8可以看出。

1）在降雨初期邊坡模型的水平和豎向位移基本沒(méi)有變化，說(shuō)明降雨初期邊坡處于穩(wěn)定的狀態(tài)，與降雨強(qiáng)度600 mm/d時(shí)邊坡的變化情況基本一致。隨著降雨時(shí)間的逐漸增加，邊坡的含水率也逐漸增加。當(dāng)坡體某一區(qū)域的含水率達(dá)到飽和狀態(tài)時(shí)，在滲流的作用下邊坡會(huì)出現(xiàn)明顯的位移。

2）隨著降雨時(shí)間的持續(xù)增加，邊坡中部開(kāi)始出現(xiàn)滑動(dòng)的現(xiàn)象。當(dāng)降雨時(shí)間達(dá)到13 300 s時(shí)，邊坡中部的豎向位移突然大幅度的增加，坡體破壞，中部標(biāo)志牌也脫落。

對(duì)比邊坡模型分別在600，800 mm/d降水強(qiáng)度下其水平、豎向位移情況可以發(fā)現(xiàn)：邊坡的變形破壞是一個(gè)逐漸發(fā)展的過(guò)程。隨著降雨時(shí)間的增加，邊坡內(nèi)部的土體慢慢達(dá)到飽和狀態(tài)。坡體水平方向的位移變化相對(duì)較緩，當(dāng)降雨時(shí)間持續(xù)到一定的時(shí)長(zhǎng)后，邊坡中部出現(xiàn)破壞，其豎向位移快速增加，最終使得整個(gè)邊坡變形破壞。對(duì)比600，?800 mm/d?降水強(qiáng)度對(duì)邊坡的影響可以看出：降水強(qiáng)度越大，邊坡模型變形破壞所需要的時(shí)間越短，邊坡越容易破壞。

4?降雨誘發(fā)邊坡變形破壞有限元模擬研究

為了能夠驗(yàn)證層狀砂質(zhì)板巖邊坡模型在降雨條件下的監(jiān)測(cè)試驗(yàn)結(jié)果，筆者結(jié)合有限元模擬方法分析降雨強(qiáng)度為800 mm/d時(shí)邊坡的變形破壞模式。對(duì)比試驗(yàn)監(jiān)測(cè)數(shù)據(jù)和模擬結(jié)果，研究降雨及滲流條件下邊坡的破壞機(jī)理。

4.1?有限元模型建立

層狀砂質(zhì)板巖邊坡中主要成分的巖土性質(zhì)在物理性能上類似于混凝土，在模擬邊坡變形破壞過(guò)程中，筆者采用ANSYS軟件中solid164單元模擬邊坡的巖體。在有限元分析中，網(wǎng)格劃分得越密集，其模擬的結(jié)果越準(zhǔn)確，但是這對(duì)計(jì)算機(jī)的計(jì)算性能提出了很高的要求。

筆者將巖體的網(wǎng)格劃分為0.5 m[20]。通過(guò)?LS-DYNA?進(jìn)行建模并導(dǎo)出k文件，根據(jù)模型的特征對(duì)相關(guān)的參數(shù)進(jìn)行修改。取坡頂、坡腳和坡中的相應(yīng)位置作為監(jiān)測(cè)點(diǎn)，分析坡體在降雨及滲流作用下的位移變化情況。在實(shí)際模擬邊坡變形破壞的過(guò)程中，忽略邊坡坡面上的位移，邊坡與模型框架之間水平向和側(cè)向分別取2個(gè)和1個(gè)約束。模擬的邊坡模型約束及網(wǎng)格劃分情況如圖9所示。

4.2?降雨強(qiáng)度為800 mm/d時(shí)邊坡模擬結(jié)果分析

當(dāng)降雨強(qiáng)度為800 mm/d時(shí)，分別取邊坡底、坡中、坡頂?shù)臄?shù)值模擬數(shù)據(jù)進(jìn)行分析。邊坡X方向上的模擬位移結(jié)果如圖10所示，其中A1114，B3392，C6260分別表示的為邊坡的坡頂、坡中、坡底的數(shù)據(jù)。可以看出：邊坡的位移曲線呈現(xiàn)拋物線的形狀，其中主要是以邊坡中部的變化為主，邊坡下部基本無(wú)變化。

模型試驗(yàn)的水平方向位移曲線如圖11所示。將邊坡模擬試驗(yàn)結(jié)果與實(shí)際的模型位移結(jié)果對(duì)比可以發(fā)現(xiàn)：兩者曲線的變化趨勢(shì)基本一致，都是邊坡中部變化最為明顯，底部基本無(wú)變化。隨著降雨時(shí)間的增加，邊坡模型試驗(yàn)的最終監(jiān)測(cè)結(jié)果顯示邊坡出現(xiàn)標(biāo)志牌脫落，邊坡中部開(kāi)始塌落，直至最終邊坡變形破壞。而邊坡模擬的結(jié)果顯示隨著降雨時(shí)間的增加邊坡最終位移不再變化，邊坡徹底破壞。邊坡模型試驗(yàn)和有限元模擬結(jié)果能夠很好地吻合。

圖12所示為降雨強(qiáng)度為800 mm/d時(shí)邊坡試驗(yàn)與豎直方向的應(yīng)變?cè)茍D?？梢钥闯觯哼吰履M的變形形態(tài)與實(shí)際邊坡模型的變形形態(tài)基本一致，兩者相似度很高，最終的破壞形態(tài)也基本類似。這也證明了選擇的有限元數(shù)值模擬方法的合理性和有效性。

5?結(jié)?語(yǔ)

筆者以某公路旁的層狀砂質(zhì)板巖邊坡作為研究對(duì)象，通過(guò)模型試驗(yàn)和有限元數(shù)值模擬分析了其在降雨及滲流條件下的變形破壞模式，得出如下結(jié)論。

1）持續(xù)降雨條件下，層狀砂質(zhì)板巖邊坡由于水滲流到巖石土體中，土體的黏聚力減小，邊坡容易失穩(wěn)。當(dāng)邊坡開(kāi)始緩慢地出現(xiàn)位移現(xiàn)象時(shí)，證明邊坡已達(dá)到失穩(wěn)的臨界狀態(tài)，此時(shí)邊坡的滑移面減小，加大了流體對(duì)邊坡裂隙的拉力，導(dǎo)致其加速失穩(wěn)，最終出現(xiàn)崩塌現(xiàn)象。

2）通過(guò)計(jì)算原始邊坡的尺寸，設(shè)計(jì)了一套包含降雨系統(tǒng)的縮尺邊坡模型。記錄了邊坡在降雨過(guò)程中不同位置的位移變化情況，分析了層狀巖質(zhì)邊坡在降雨及滲流作用下的持續(xù)變形過(guò)程。降雨開(kāi)始時(shí)無(wú)變化，當(dāng)降雨時(shí)間達(dá)到一定時(shí)長(zhǎng)時(shí)，邊坡會(huì)出現(xiàn)滲流現(xiàn)象。此時(shí)邊坡的位移會(huì)進(jìn)一步的增加，其中部的水平、豎向位移出現(xiàn)得更快，形成一個(gè)D字狀的位移曲線，直至坡體破壞。

3）邊坡的位移曲線呈拋物線形狀，水平位移的變化幅度相比豎向位移更大，表明邊坡的變形破壞主要以水平位移變化為主。當(dāng)降雨強(qiáng)度為600 mm/d且降雨時(shí)間持續(xù)16 800 s時(shí)，邊坡中部位移最大，此時(shí)坡體中部的標(biāo)志牌滑落，邊坡破壞;當(dāng)降雨強(qiáng)度為800 mm/d且降雨時(shí)間持續(xù)13 300 s時(shí)，邊坡中部的豎向位移突然大幅度增加，坡體破壞，中部標(biāo)志牌脫落。

4）利用ANSYS有限元模擬的邊坡變形破壞與實(shí)際邊坡模型的變形形態(tài)基本一致，最終的破壞形態(tài)基本類似。

模型試驗(yàn)中的邊坡穩(wěn)定性是基于表現(xiàn)位移進(jìn)行判斷的，實(shí)際工程中邊坡的穩(wěn)定性問(wèn)題往往始于內(nèi)部，因此，在今后的研究中應(yīng)從表現(xiàn)位移向微觀破壞機(jī)理深入研究，以便進(jìn)一步掌握巖質(zhì)邊坡的變形規(guī)律。

參考文獻(xiàn)/References：

[1]劉震濤，尚彥軍，邵鵬，等.降雨與開(kāi)挖聯(lián)合作用下邊坡位移矢量及速率變化分析——以韓江高陂水利樞紐右岸尾水渠邊坡為例[J].工程地質(zhì)學(xué)報(bào)，2020，28（1）：122-131.

LIU Zhentao， SHANG Yanjun， SHAO Peng， et al. Analysis on variations of displacement rates under influences of rainfall and excavation-Taking right bank slope of Gaobei key water control project as case study[J].Journal of Engineering Geology，2020，28（1）：122-131.

[2]劉暢，季凡凡，鄭剛，等.降雨對(duì)軟土基坑支護(hù)結(jié)構(gòu)影響實(shí)測(cè)及機(jī)理研究[J].巖土工程學(xué)報(bào)，2020，42（3）：447-456.

LIU Chang， JI Fanfan， ZHENG Gang， et al. Measurement and mechanism of influences of rainfall on supporting structures of foundation pits in soft soils[J].Chinese Journal of Geotechnical Engineering， 2020，42（3）：447-456.

[3]張碩，裴向軍，黃潤(rùn)秋，等.黃土填方邊坡降雨入滲特征及變形破壞模式的模型試驗(yàn)[J].中國(guó)公路學(xué)報(bào)，2019，32（9）：32-41.

ZHANG Shuo， PEI Xiangjun， HUANG Runqiu， et al. Model test on seepage characteristics and deformation failure modes of loess fill slope under rainfall [J]. China Journal of Highway and Transport， 2019，32 （9）： 32-41.

[4]韓廷文，李軍，康曉鵬.降雨條件下的節(jié)理化露天礦高邊坡穩(wěn)定性分析[J].科學(xué)技術(shù)與工程，2019，19（19）：75-80.

HAN Tingwen， LI Jun， KANG Xiaopeng. Stability analysis of high slope of jointed open-pit mine under rainfall conditions [J].Science Technology and Engineering， 2019，19 （19）： 75-80.

[5]SUN Dongmei， LI Xiaomin， FENG Ping， et al. Stability analysis of unsaturated soil slope during rainfall infiltration using coupled liquid-gas-solid three-phase model [J]. Water Science and Engineering，2016，9（3）：193-194.

[6]劉廣寧，李聰，盧波，等.降雨誘發(fā)全-強(qiáng)風(fēng)化巖邊坡淺層失穩(wěn)模型試驗(yàn)研究[J/OL].長(zhǎng)江科學(xué)院院報(bào)， http：//kns.cnki.net/kcms/detail/42.1171.TV.20191104.1013.002.html， 2019-11-05.

LIU Guangning， LI Cong， LU Bo， et al. Model test of shallow failure of fully strong weathered rock slope induced by rainfall [J].Journal of Yangtze River Scientific Research Institute，?http：?//kns.cnki.net/kcms/detail/42.1171.TV.20191104.1013.002.html，2019-11-05.

[7]楊世平，陳泉.四川涼山州某水庫(kù)右岸趾板邊坡的變形特征及機(jī)制研究[J].三峽大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，41（sup1）：103-107.

YANG Shiping， CHEN Quan. Deformation characteristics and mechanism of toe-slab slope on right bank of a Liangshan reservoir in Sichuan province[J].Journal of China Three Gorges University （Natural Sciences）， 2019，41 （sup1）： 103-107.

[8]何忠明，王保林，胡慶國(guó).強(qiáng)降雨條件下含軟弱夾層粘性土坡穩(wěn)定性分析[J]. 土木建筑與環(huán)境工程， 2018，40（5）：109-116.

HE Zhongming， WANG Baolin， HU Qingguo. Stability analysis of cohesive soil slopes with weak interlayer under heavy rainfall[J].Journal of Civil，Architectural & Environmental Engineering， 2018，40 （5）： 109-116.

[9]劉杰，曾鈴，付宏淵，等.土質(zhì)邊坡降雨入滲深度及飽和區(qū)變化規(guī)律[J].中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，50（2）：452-459.

LIU Jie， ZENG Ling， FU Hongyuan， et al. Variation law of rainfall infiltration depth and saturation zone of soil slope[J]. Journal of Central South University（Science and Technology）， 2019，50 （2）： 452-459.

[10]孫萍，王剛，李榮建，等.降雨條件下黃土邊坡現(xiàn)場(chǎng)試驗(yàn)研究[J].工程地質(zhì)學(xué)報(bào)，2019，27（2）：466-476.

SUN Ping， WANG Gang， LI Rongjian， et al. Study on field test of loess slope?under?the artificial rainfall condition[J].Journal of Engineering Geology， 2019，27 （2）： 466-476.

[11]王磊，李榮建，潘俊義，等.隔離邊界條件下連續(xù)降雨誘發(fā)黃土邊坡開(kāi)裂試驗(yàn)研究[J].天津大學(xué)學(xué)報(bào)（自然科學(xué)與工程技術(shù)版），2019，52（11）：1163-1170.

WANG Lei， LI Rongjian， PAN Junyi， et al. Experimental study on loess slope cracking induced by continuous rainfall?under?isolated boundary[J]. Journal of Tianjin University （Science and Technology）， 2019，52 （11）： 1163-1170 .

[12]AMINI M， MAJDI A， VESHADI M A. Stability analysis of rock slopes against block-flexure toppling failure[J].Rock?Mechanics?and Rock Engineering，2012，45（4）：519-532.

[13]PARK J W， PARK C， SONG J W， et al. Polygonal grain-based distinct element modeling for mechanical behavior of brittle rock[J]. International Journal for Numerical and Analytical Methods in Geomechanics， 2017，41（6）：880-898.

[14]LI Longqi，JU Nengpan，ZHANG Shuai，et al. Seismic wave propagation characteristic and its effects on the failure of steep jointed anti-dip rock slope[J]. Landslides，2019，16（1）： 105-123.

[15]龍安發(fā)，陳開(kāi)圣，季永新.不同降雨強(qiáng)度下紅黏土邊坡干濕循環(huán)試驗(yàn)研究[J].巖土工程學(xué)報(bào)，2019，41（sup2）：193-196.

LONG Anfa， CHEN Kaisheng， JI Yongxin. Experimental study on wetting-drying cycles of red clay slopes under different rainfall intensities[J].Chinese Journal of Geotechnical Engineering， 2019，41 （sup2）： 193-196.

[16]金亞兵，黃健，梁軍，等.降雨—變形耦合預(yù)警模型在荔景山莊邊坡的應(yīng)用[J].廣西大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，44（3）：834-844.

JIN Yabing， HUANG Jian， LIANG Jun， et al. Application of rainfall-deformation coupling early warning model in Lijing mountain slope [J].Journal of Guangxi University （Natural Science Edition）， 2019，44 （3）： 834-844.

[17]劉新榮，鄧志云，劉永權(quán)，等.地震作用下水平層狀巖質(zhì)邊坡累積損傷與破壞模式研究[J].巖土力學(xué)，2019，40（7）：2507-2516.

LIU Xinrong， DENG Zhiyun， LIU Yongquan， et al. Study of cumulative damage and failure mode of horizontal layered rock slope subjected to seismic loads[J]. Rock and Soil Mechanics， 2019，40 （7）： 2507-2516.

[18]李玉瑞，吳紅剛，趙金，等.模擬降雨作用樁錨-加筋土組合結(jié)構(gòu)加固邊坡研究[J].防災(zāi)減災(zāi)工程學(xué)報(bào)，2019，39（3）：516-523.

LI Yurui， WU Honggang， ZHAO Jin， et al. Study on simulated rainfall effect on pile-anchor-reinforced soil combined structure for slope reinforcement[J].Journal of Disaster Prevention and Mitigation Engineering， 2019，39 （3）： 516-523.

[19]梁樹(shù)，謝強(qiáng)，郭永春，等.降雨入滲條件下膨脹土基坑邊坡離心試驗(yàn)[J].地質(zhì)科技情報(bào)，2019，38（2）：249-255.

LIANG Shu， XIE Qiang， GUO Yongchun， et al. Centrifugal test of the expansive soil foundation slope under rainfall infiltration[J]. Geological Science and Technology Information， 2019，38 （2）： 249-255.

[20]許方領(lǐng)，阮懷寧，黃肖，等.考慮降雨入滲的非飽和土邊坡流固耦合數(shù)值分析[J].河南科學(xué)，2019，37（6）：955-961.

XU Fangling， RUAN Huaining， HUANG Xiao， et al. Fluid solid coupling numerical analysis of unsaturated soil slope considering rainfall infiltration[J]. Henan Science， 2019，37 （6）： 955-961.