紫花苜蓿對黃土邊坡淺層破壞防護時間效應的數(shù)值分析

毛正君，張瑾鴿，畢銀麗，孫偉博，安 寧

紫花苜蓿對黃土邊坡淺層破壞防護時間效應的數(shù)值分析

毛正君1,2,3，張瑾鴿1，畢銀麗1，孫偉博4，安 寧1

（1. 西安科技大學地質與環(huán)境學院，西安 710054；2. 西安科技大學煤炭綠色開采地質研究院，西安 710054；3. 陜西省煤炭綠色開發(fā)地質保障重點實驗室，西安 710054；4. 西安科技大學能源學院，西安 710054）

為探究紫花苜蓿對黃土邊坡淺層破壞防護的時間效應，該研究考慮降雨入滲條件并以含不同生長期紫花苜蓿的黃土邊坡為例，基于室內土工試驗獲取參數(shù)，分別采用含植物根系的無限邊坡模型和數(shù)值模擬兩種方法進行穩(wěn)定性分析。結果表明：降雨作用使得邊坡潛在滑動面位置由坡體內部轉移到降雨最大入滲深度處；隨著坡面植被的生長，邊坡淺層土體的強度逐漸提高，潛在滑動面又逐漸轉移到邊坡內部。根據(jù)數(shù)值模擬強度折減法，邊坡穩(wěn)定系數(shù)隨著植被生長不斷增大。對比無限邊坡模型與數(shù)值模擬方法計算的邊坡穩(wěn)定系數(shù)，植物根系生長到1.0 m且不論在天然狀態(tài)或降雨條件下黃土邊坡均處于穩(wěn)定狀態(tài)，植物根系在0～0.8 m且在降雨條件下時處于不穩(wěn)定狀態(tài)，隨著邊坡植物生長期的延長，邊坡穩(wěn)定系數(shù)總體變化趨勢一致，即先降低再逐漸增大。草本植物隨著生長期的延長，對降雨條件下黃土邊坡坡面侵蝕有明顯的抑制作用，從而草本植物根系生長能夠提高黃土邊坡穩(wěn)定性。當草本植物生長時間達到150 d時，能夠有效防治黃土坡面侵蝕，并提高黃土邊坡穩(wěn)定性，因此草本植物生長的前5個月為關鍵期。研究成果對于黃土邊坡淺層破壞防護以及實現(xiàn)黃土地區(qū)可持續(xù)發(fā)展具有重要的理論意義和應用價值。

邊坡穩(wěn)定性；數(shù)值模擬；黃土；紫花苜蓿；邊坡淺層破壞；時間效應；無限邊坡模型

0 引 言

黃土是一種碎散的多孔多相顆粒介質，由于其獨特的物質來源、特有的生成環(huán)境及應力-應變歷史而具有特殊的微結構[1]。當黃土遇水浸濕或遭受強震作用時，其結構會出現(xiàn)崩潰性破壞，導致強度突發(fā)性完全喪失，從而誘發(fā)嚴重的地質災害或者工程病害[2]。黃土邊坡在降雨作用下往往先出現(xiàn)邊坡坡面侵蝕，隨著降雨次數(shù)的增加將會導致邊坡局部破壞，甚至影響邊坡的整體穩(wěn)定性[3]。黃土梁峁地形切割強烈，不利于降水入滲及地下水賦存，因此降雨直接入滲的影響深度是有限的[4]；表層黃土軟化、質量增加，造成數(shù)量多、規(guī)模小、多發(fā)于30°～50°黃土斜坡上的邊坡淺層破壞[5]。中國黃土高原作為世界上黃土沉積厚度最大地區(qū)之一，黃土邊坡淺層破壞數(shù)目每年高達數(shù)百個[6-7]。黃土邊坡淺層破壞一般發(fā)生在坡體表層或坡面下2 m的范圍內，其破壞類型有坡面侵蝕、風化剝落、溜坍、錯落、淺層滑坡等[8]，可造成黃土地區(qū)土地可持續(xù)生產(chǎn)力的損失和破壞，是當前黃土地區(qū)研究的重點。

自20世紀60年代末以來，計算土力學一直是熱點研究領域[9]。近年來各種數(shù)值模擬計算技術的迅猛發(fā)展，為土體破壞現(xiàn)象的模擬提供了有力的手段[10]。通過合理的對根系加固作用進行數(shù)值模擬，進而對根系加固土體的機理進行分析，能為植物防護技術的應用提供強有力的理論依據(jù)[11]。周群華等[12]采用有限元數(shù)值模擬方法，對植物根系錨固邊坡的作用及其力學機理進行分析，發(fā)現(xiàn)根系固土是因為根系增加了土體的抗剪強度，從而增強了邊坡土體的抗滑力。Temgoua等[13]采用SIMULIA軟件建立三維數(shù)值模型，發(fā)現(xiàn)整個邊坡的穩(wěn)定性主要取決于根系的深度和根系提供的附加黏聚力。Li等[14]采用數(shù)值模擬的方法，將植物根系結構概化，分析了植物根系空間布局對邊坡穩(wěn)定性的影響。黃建坤等[15]基于有限元軟件，模擬了含草邊坡和無草邊坡的應力和應變規(guī)律，通過邊坡安全系數(shù)反映了根系的固土效應。由于高密度的細根能夠和土體顆粒緊密地纏繞在一起，草本植物能有效地提高根土復合體的穩(wěn)定性[16]，是防治黃土邊坡淺層破壞土木工程解決方案的一種生態(tài)替代方案[17]。紫花苜蓿屬多年生草本植物，具有繁殖能力強、適應性強、經(jīng)濟效益高等特點[18]，能夠有效提高土體的抗沖刷能力，并且根系的固土能力強，具有顯著的水土保持效果，是中國黃土高原地區(qū)水土保持的先鋒植物之一[19]。

為了探究紫花苜蓿對黃土邊坡淺層破壞防護的時間效應，本文采用套管人工種植紫花苜蓿的方式，獲取含紫花苜蓿根系的原狀土樣進行試驗，并分別利用數(shù)值模擬和無限邊坡模型，對不同生長期紫花苜蓿天然狀態(tài)和降雨條件下黃土邊坡的不利情況進行穩(wěn)定性分析，探究不同生長期紫花苜蓿對黃土邊坡淺層破壞防護作用。研究成果對于黃土邊坡淺層破壞防護以及實現(xiàn)黃土地區(qū)可持續(xù)發(fā)展具有重要的理論意義和應用價值。

1 材料與方法

1.1 試驗材料與試樣制備

試驗用土取自陜西省西安市臨潼區(qū)，地理坐標為109°20′25″E，34°25′46″N。試驗用土為離石黃土（Q2黃土），顏色為褐黃色，天然密度為1.67g/cm3，天然含水率為16.3%，塑限為17.4%，塑性指數(shù)為19.8%，液限為37.2%，液限指數(shù)為?0.06，不均勻系數(shù)（C）和曲率系數(shù)（C）分別為8.36和0.78，為不均勻級配土。

試驗所用紫花苜蓿種子來自內蒙古赤峰市林西縣草原站，采用人工種植培養(yǎng)的方式，獲得含不同生長期紫花苜蓿根系的黃土試樣。方法如下：首先，將試驗用土取回后經(jīng)自然風干，挑出其中的礫石與垃圾，把土塊砸碎后過2 mm篩，之后置于陰涼、干燥處保存；其次，在種植箱內鋪一層5 cm左右的土，將切割好的內徑為67.8 mm，外徑為75 mm的PVC種植管插入種植箱內，并在PVC種植管中裝滿黃土，在PVC種植管上端澆水浸潤土體后靜置一段時間，繼續(xù)填滿黃土并澆水浸潤，直至黃土裝滿PVC種植管澆水后不再下沉；再次，在每根PVC種植管中撒上數(shù)量大致相同的紫花苜蓿種子，生長期為60與90 d的紫花苜蓿（根長約15～20 cm）種植于20 cm長的PVC管中，生長期為120與150 d的紫花苜蓿（根長約30～40 cm）種植于40 cm長的PVC管中；最后，為了保證PVC種植管中黃土的含水率大致相同，在養(yǎng)護過程中定時給PVC管內澆入相同質量的水。

為了滿足三軸試驗中直徑為39.1 mm、高為80 mm的標準三軸樣的要求，且便于對比，對于兩組根土復合體土樣，取樣位置分別選取20 cm長PVC管的中部和40 cm長PVC管上半部分20 cm管的中部。

1.2 根土復合體力學指標求解

本文考慮邊坡表層1 m內完全飽和，當在強度相對較高且埋深較淺的根系分布范圍內已然發(fā)生了破壞，那么在強度較低且埋深較大的根系分布范圍內必然更容易發(fā)生破壞。此外，所取三軸樣即根土復合體土樣被看做一個整體進行試驗。因此，考慮植物不同生長期為變量，根系分布對于本文試驗抗剪指標測定的影響忽略。采用中國南京土壤儀器廠生產(chǎn)的ZJ型應變控制式直剪儀，進行直剪快剪試驗測定天然素黃土的抗剪強度指標黏聚力c和內摩擦角。采用中國南京土壤儀器廠生產(chǎn)的TSZ-1型三軸儀，進行不固結不排水（UU）三軸試驗測定天然素黃土的彈性模量和泊松比，固結不排水（CU）三軸試驗測定飽和素黃土和含不同生長期紫花苜蓿根系黃土試樣的抗剪強度指標黏聚力c、內摩擦角、彈性模量和泊松比。

根據(jù)直剪試驗剪切應力與剪切位移關系曲線，選取峰值強度作出不同垂直壓力下所對應的抗剪強度坐標點，并通過擬合求解抗剪強度指標。本文根據(jù)三軸試驗應力應變曲線畫出了應力莫爾圓，但各個應力圓之間無明顯規(guī)律，難以繪制強度包絡線，因此選用最小二乘法求解抗剪強度指標。

根據(jù)三軸試驗得出的應力-應變關系曲線可以求得土的變形模量和泊松比，它們不是常量，而是隨應力水平變化的參數(shù)。在常規(guī)三軸試驗中，當試樣只承受豎向應力增量Δ1時，除了產(chǎn)生豎向應變Δ1之外，同時還產(chǎn)生側向應變Δ3（側向膨脹）。側向應變增量Δ3與豎向應變增量Δ1的比值可用彈性力學中的泊松比表示。

在不允許有側向變形的情況下，土樣受到軸向壓力增量Δ1將會引起側向壓力的相應增量Δ3，二者比值（Δ3/Δ1）稱為土的側壓力系數(shù)或靜止土壓力系數(shù)0：

側壓力系數(shù)與泊松比的關系為

1.3 含植物根系的無限邊坡模型

現(xiàn)實中的天然邊坡和人工邊坡都是有限長的，但對于一些具有平行于坡面的淺層滑動面的長坡，滑動面的上部和坡腳處的滑動面彎曲部分、兩側邊界的約束都可以忽略，可以視為二維的無限邊坡[20]。無限邊坡穩(wěn)定性分析模型常常應用于淺層滑坡的穩(wěn)定性分析[21]和坡面穩(wěn)定性分析[22]。已有研究認為植物根系擁有防止淺層滑坡和地表侵蝕的潛能[23]，因此無限邊坡穩(wěn)定性分析模型適用于含植物根系邊坡穩(wěn)定性分析[24-26]。此外，根土復合體與素土的內摩擦角基本一致，僅黏聚力發(fā)生了變化[27-28]。黃土地區(qū)地下水位埋深較大[29-30]，降雨入滲深度較淺[5]，本文假定不存在地下水滲流作用，且降水入滲影響范圍內的土體完全飽和，邊坡穩(wěn)定系數(shù)可以表示為

式中F是邊坡穩(wěn)定系數(shù)；是邊坡土體的重度，kN/m3；是滑動面深度，m；是坡面與水平地面的夾角，(°)；是內摩擦角，(°)；c是根土復合體的黏聚力，kPa；h是植物根系分布深度，m；c是素土的黏聚力，kPa。

紫花苜蓿根系分布深度較淺[31]，本文假定含紫花苜蓿根系的邊坡土體厚度為1 m。已有研究[5,32]表明黃土地區(qū)降水入滲深度約為1～3 m，本文僅考慮一次降雨事件造成的影響，因此將降水入滲影響深度假定為1 m。本文僅考慮含紫花苜蓿根系的淺層邊坡土體穩(wěn)定性，由于邊坡穩(wěn)定系數(shù)與滑動面深度成反比，因此邊坡穩(wěn)定系數(shù)最小值對應的滑動面深度=1.0 m。將黏聚力、內摩擦角、重度以及滑動面深度=1.0 m代入式（1），發(fā)現(xiàn)坡度在49°～51°范圍內邊坡穩(wěn)定系數(shù)F最小，因此本文假定坡度=50°。

1.4 含植物根系的邊坡穩(wěn)定性數(shù)值模擬分析

與現(xiàn)場試驗和模型試驗方法相比，數(shù)值模擬需要的人力和財力更少，研究周期更短，靈活性和便利性更好[33]。數(shù)值模擬是用數(shù)值方法定量地描述一個物理系統(tǒng)的演化過程，允許分析復雜的因素，如非線性行為，包含復雜的不連續(xù)模式、動態(tài)載荷、非均勻性、各向異性、流體流動和力學行為的耦合、材料軟化、時間依賴的行為等，還可以不受限制地設置任意時間尺度或長度的數(shù)值，可以快速、有效地改變測試條件和組件[34]。強度折減法是指采用降低巖土強度的方法來計算巖土工程的穩(wěn)定系數(shù)，通過不斷降低邊坡巖土體的抗剪強度參數(shù)來逐步逼近極限破壞狀態(tài)，利用數(shù)值模擬計算程序自動計算得到滑動破壞面和強度折減穩(wěn)定系數(shù)；運用強度折減法分析巖土體穩(wěn)定性具有不需要假定滑動面的形狀和位置、不需進行條分和離散的優(yōu)勢[8]。

1.4.1 工況設定

含不同生長期紫花苜蓿根系邊坡穩(wěn)定性數(shù)值模擬分析，共設置6種工況，見表1。植被的適宜播種時間是秋季和春季，黃土地區(qū)的雨季一般在7－9月，為了使植被護坡能夠在同年中起到預期效果，春季是進行播種的最佳時機。假定春季3月進行播種，則到雨季來臨時植物的生長周期可能是150 d，因此本次試驗中采樣周期選擇60、90、120和150 d共4個采樣周期。本文所選草本植物紫花苜蓿，屬多年生深根系草本植物，主根入土深，在播種當年根系深度可達2 m左右[35]，為考慮野外邊坡紫花苜蓿的實際生長情況，反映植物根系隨生長期的變化特征，假定植物生長期為60、90、120、150 d時，根系深度分別為0.4、0.6、0.8、1.0 m。

表1 工況設定

1.4.2 模型建立

本文基于數(shù)值模擬軟件FLAC3D，進行含不同生長期紫花苜蓿根系黃土邊坡淺層破壞防護效應分析。根據(jù)基本假設，含紫花苜蓿根系的邊坡土體厚度隨生長時間分別為0.4、0.6、0.8和1.0 m，降水入滲影響深度為1 m，坡度=50°。為了達到理想計算精度，坡腳到左端邊界的距離設置為坡高12 m的1.5倍即18 m，坡肩到右端邊界的距離為坡高12 m的2.5倍即30 m，且上下邊界總高等于坡高12 m的2倍即24 m[36]，建立邊坡幾何模型，如圖1所示。模型網(wǎng)格劃分后共包含87 720個六面體單元和96 938個節(jié)點。圖1中①分區(qū)為受到降雨作用和植物根系共同影響的土體，①分區(qū)隨著植物根系的發(fā)育而增大；①分區(qū)與②分區(qū)的分界線為植物根系的發(fā)育深度，②分區(qū)隨著植物根系的發(fā)育而減小，為僅受到降雨作用影響的土體；③分區(qū)為僅受到降雨影響的土體；④分區(qū)為天然狀態(tài)土體。屈服準則采用Mohr-Coulomb準則，并將素土和根-土復合體都視為均質材料進行模擬分析。模型的側面和底面設置為固定邊界，模型的上表面設置為自由邊界，模型在計算過程中僅受重力荷載的作用。在模型中設置了34個監(jiān)測點（各監(jiān)測點位由上而下依次為邊坡表面、①分區(qū)中點、①②分區(qū)界面、②分區(qū)中點、②④分區(qū)界面和坡面豎直向下2 m處；斷面6在邊坡內部，故放置4個點位），用于監(jiān)測計算過程中的水平位移和豎向位移的變化，如圖2所示。

注：①～④為4個分區(qū)?！　ote: ①-④are 4 partitions.

注：圓點表示每個斷面下各監(jiān)測點（層位），從上到下依次為1～6號點，記為1#～6#。

1.4.3 參數(shù)取值

重度通過環(huán)刀法測得，抗拉強度σ參考文獻[37]獲得，剪脹角參考文獻[38]取值為0°。不同狀態(tài)下黃土物理力學指標參數(shù)取值如表2所示，不同工況下各分區(qū)參數(shù)選取如表3所示。

表2 不同狀態(tài)下黃土物理力學指標參數(shù)取值

注：素黃土與含紫花苜蓿根系黃土的重度值略有波動，此處將其看作定值。：重度；：內摩擦角；c：根土復合體的黏聚力；c：素土的黏聚力；：彈性模量；σ：抗拉強度。

Note: The gravity values of plain loess and root loess containing alfalfa fluctuated slightly, which was regarded as a fixed value here.: unit weight of soil;: internal friction angle;c: cohesion of root soil complex;c: cohesion of plain soil;: elastic modulus;σ:tensile strength.

表3 不同工況下各分區(qū)參數(shù)選取

注：A、B、C、D、E、F指代的為表2中六種狀態(tài)；工況6時植物根系深度為1 m，此時沒有②分區(qū)。

Note: A, B, C, D, E, F refer to the six states in table 2; under condition 6, the root depth was 1 m, and there was no ②partition.

2 結果與分析

2.1 數(shù)值模擬分析結果

2.1.1 自重作用下邊坡監(jiān)測點水平和豎向位移

本文工況2～工況5的計算結果不收斂，工況2～工況5在計算至450時步時塑性區(qū)都發(fā)生貫通，本文以此作為計算結束的標準，此時邊坡正處于極限平衡狀態(tài)。限于篇幅，圖3和圖4僅展示了模型的監(jiān)測斷面1和坡面層位1的位移情況，能夠反映模型的整體位移變化規(guī)律。從工況1的監(jiān)測點水平位移（圖3a）和豎向位移（圖4a）可以看出：監(jiān)測點水平位移隨時步增加，位移量逐漸增大，且都為負值，其中邊坡坡肩處（監(jiān)測斷面5）位移值最大；豎向位移各點位移量出現(xiàn)波動，但量值均很?。贿吰缕旅媾c素黃土邊坡整體的位移無明顯差異。工況2～工況5監(jiān)測點的水平位移均為負值，在邊坡坡腳處監(jiān)測斷面1-1#點的水平位移量最大，從邊坡坡腳到邊坡坡肩各監(jiān)測斷面水平位移量總體上越來越??；邊坡坡腳處由于受到擠壓隆起飽水層或根土復合體層與無根系飽水層豎向位移量大且為正值，其余監(jiān)測點豎向位移量較小；邊坡根土復合體層與無根系飽水層位移量無明顯差異，但兩者與天然狀態(tài)土體位移量具有明顯差異，表明邊坡破壞面可能處于邊坡飽水層與天然土體的交界處。但在工況5時，邊坡根土復合體層與天然狀態(tài)土體位移量具有明顯差異，無根系飽水層位移量與天然狀態(tài)土體位移量差異較小，表明邊坡破壞面可能處于邊坡飽水層。工況6監(jiān)測點的水平位移均為負值，在靠近邊坡坡腳處監(jiān)測斷面2的水平位移量最大，從邊坡靠近坡腳到邊坡坡肩各監(jiān)測斷面水平位移量總體上越來越小，邊坡坡腳處監(jiān)測斷面1的水平位移量最小；邊坡坡腳處監(jiān)測斷面1的豎向位移為正值，其余監(jiān)測斷面豎向位移均為負值，且越靠近邊坡坡肩其絕對值越大；坡面表層與素黃土邊坡整體的位移無明顯差異。

2.1.2 自重作用下邊坡水平位移與塑性區(qū)

由圖5和圖6可知，工況2～工況5在自重作用下邊坡的最大位移出現(xiàn)在邊坡坡面；隨著計算進行，首先在邊坡坡面以下1 m處即降水最大入滲深度處形成了平行于坡面的剪切屈服塑性區(qū)，之后邊坡坡面中上部位置的少量土體開始因拉伸而產(chǎn)生剪切屈服，并隨后剪切屈服區(qū)域沿邊坡坡面逐漸向邊坡坡肩方向移動。工況6在自重作用下邊坡的最大位移發(fā)生在邊坡坡底面靠近邊坡坡腳處，邊坡整體有向左滑動的趨勢，邊坡坡底面靠近邊坡坡腳處由于擠壓出現(xiàn)隆起，此時邊坡整體處于穩(wěn)定狀態(tài)；從塑性區(qū)分布來看，此時僅在邊坡坡底面靠近邊坡坡腳處有少部分網(wǎng)格發(fā)生剪切屈服。在圖5中隨著紫花苜蓿生長期的延長，邊坡的最大水平位移依次為0.193、0.194、0.203、0.261、0.394 cm，都出現(xiàn)在坡面位置，當紫花苜蓿生長至150 d即工況6時，最大水平位移僅出現(xiàn)在坡腳位置，坡面水平位移明顯減小，且相對于整體邊坡而言該水平位移量相對較小。

2.1.3 強度折減法計算下的最大剪應變增量和塑性區(qū)

工況1采用強度折減法計算的邊坡穩(wěn)定系數(shù)為1.457，此時邊坡處于穩(wěn)定狀態(tài)。由圖7a可以看出，邊坡的潛在滑動面在坡體內部較深的位置，沿著潛在滑動面出現(xiàn)大量的剪切屈服區(qū)域，即圖8a中紅色區(qū)域，邊坡發(fā)生失穩(wěn)破壞時將從坡腳位置剪出。工況2～工況5的邊坡穩(wěn)定系數(shù)分別為0.600、0.637、0.660、0.695，此時邊坡可能會發(fā)生失穩(wěn)破壞。從圖7和圖8可以看出，在邊坡坡腳和邊坡坡肩處均出現(xiàn)剪切屈服區(qū)域，邊坡潛在滑動面發(fā)育在降雨最大入滲深度處，邊坡發(fā)生失穩(wěn)破壞時將從坡腳位置剪出。工況6的邊坡穩(wěn)定系數(shù)為1.367，此時邊坡處于穩(wěn)定狀態(tài)。從圖7f和圖8f可以看出，此時邊坡失穩(wěn)破壞可能有兩個潛在滑動面，一個在降雨入滲能夠影響的最大深度，即坡面以下1 m處，另一個在坡體內部較深位置。沿著坡體深部潛在滑動面出現(xiàn)大量的剪切屈服區(qū)域，同時塑性區(qū)已經(jīng)貫通；降雨最大入滲深度的淺部潛在滑動面還未貫通；對于這兩個潛在滑動面，邊坡發(fā)生失穩(wěn)破壞時都將從坡腳位置剪出。隨著紫花苜蓿生長期的延長，邊坡的最大剪應變增量依次為0.033、0.156、0.209、0.362、3.829、0.082，當紫花苜蓿的生長期達到150 d時，邊坡的最大剪應變增量驟然下降，此時紫花苜蓿對降雨條件下黃土邊坡坡面侵蝕具有明顯的抑制作用。

注：圖例中每條曲線的第一個數(shù)字代表斷面的編號，第二個數(shù)字代表監(jiān)測點（層位）序號。下同。

圖4 自重作用下各工況監(jiān)測點豎向位移變化曲線

圖5 自重作用下各工況450時步的水平位移云圖

注：圖中紅色區(qū)域為剪切屈服區(qū)。

圖7 基于強度折減法的各工況最大剪應變增量云圖

圖8 基于強度折減法的各工況塑性區(qū)分布圖

由數(shù)值模擬結果發(fā)現(xiàn)：工況1時黃土邊坡的潛在滑動面出現(xiàn)在坡體內部較深位置處；當處于工況2～工況5時，黃土邊坡潛在滑動面出現(xiàn)在降雨入滲能夠影響到的最大深度處，即坡面以下1 m處，表明降雨作用使得黃土邊坡可能發(fā)生失穩(wěn)破壞的潛在滑動面位置由邊坡坡體內部轉移到邊坡坡面；隨著植物生長至工況6時，黃土邊坡內部又開始發(fā)育潛在滑動面。在工況2～工況5時，黃土邊坡的破壞過程都是首先在邊坡坡面以下1 m處，即降水最大入滲深度處，形成了平行于坡面的剪切屈服塑性區(qū)，之后邊坡坡面中上部位置的少量土體開始因拉伸而產(chǎn)生剪切屈服，并隨后剪切屈服區(qū)域沿邊坡坡面逐漸向邊坡坡肩方向移動，且邊坡的最大位移都出現(xiàn)在邊坡坡腳處。因此，黃土邊坡在天然狀態(tài)下是穩(wěn)定的，而在降雨條件下邊坡坡面土體強度迅速降低，重力作用下坡面土體可能沿著降雨入滲的最大影響深度處發(fā)生失穩(wěn)破壞。隨著植被的生長，植物根系的發(fā)育，植物根系對邊坡土體的加固作用顯著增強，增大了邊坡的抗剪強度，降低了邊坡的坡面侵蝕作用，阻止了邊坡的淺層破壞，最終提高了邊坡的整體穩(wěn)定性。

2.2 無限邊坡模型結果

將已測得的黏聚力、內摩擦角、重度（表2）、坡面與水平地面的夾角=50°以及滑面深度=1.0 m代入無限邊坡模型，得出無限邊坡模型計算的6種工況下的邊坡穩(wěn)定系數(shù)F值，如表4所示。當邊坡處于工況1時邊坡穩(wěn)定系數(shù)大于1，表明邊坡處于穩(wěn)定狀態(tài)。當邊坡處于工況2～工況5即受到降雨影響時，邊坡穩(wěn)定系數(shù)迅速減小，處于不穩(wěn)定狀態(tài)。但隨著邊坡坡面紫花苜蓿生長期的延長，邊坡穩(wěn)定系數(shù)逐漸增大。當達到工況6即紫花苜蓿生長150 d時，邊坡處于穩(wěn)定狀態(tài)。

2.3 數(shù)值模擬與無限邊坡模型結果對比

將含不同生長期紫花苜蓿根系的土體的物理力學參數(shù)值帶入FLAC3D數(shù)值模擬軟件與無限邊坡模型中，探究植物根系對黃土邊坡淺層破壞的防護效應。將無限邊坡模型計算的邊坡穩(wěn)定系數(shù)與數(shù)值模擬強度折減法的計算結果進行對比（表4），可以看出：黃土邊坡在工況1與工況6時處于穩(wěn)定狀態(tài)，在工況2～工況5時均處于不穩(wěn)定狀態(tài)；采用數(shù)值模擬與分析模型兩種方法計算結果總體一致。但是，在工況1條件下兩種方法的計算結果相差較大，這是由于兩種方法計算過程中的潛在滑動面位置不同所引起的，數(shù)值模擬強度折減法計算的潛在滑動面位置在坡體內部較深處（圖7），而無限邊坡模型計算的潛在滑動面假定為坡面以下1 m處。在工況2～工況5條件下，兩種方法的計算結果基本一致，這是由于降雨工況下強度折減法計算的潛在滑動面位置與無限邊坡模型假定的坡面以下1 m處潛在滑動面位置基本一致。對于工況6，數(shù)值模擬計算結果表明坡體的潛在滑動面開始從淺部向深部轉移，強度折減法計算的邊坡穩(wěn)定系數(shù)為深部潛在滑動面對應的邊坡穩(wěn)定系數(shù)，無限邊坡模型仍然以1 m作為潛在滑動面進行計算對應的邊坡穩(wěn)定系數(shù)，此時兩種計算方法結果數(shù)值上很接近，說明當植物根系發(fā)育至1 m時，即便是降雨入滲深度也達到1 m，含植物根系的黃土邊坡整體處于穩(wěn)定狀態(tài)，且黃土邊坡淺層破壞也不會發(fā)生。

數(shù)值模擬與無限邊坡模型兩種方法的結果表明，黃土邊坡在降雨條件下，其邊坡坡面土體強度迅速降低，邊坡發(fā)生失穩(wěn)破壞，但隨著植物根系的生長發(fā)育，其加固作用顯著增強，并增大了黃土邊坡的抗剪強度，降低了黃土邊坡坡面侵蝕作用，阻止了邊坡的淺層局部破壞，最終提高了黃土邊坡的整體穩(wěn)定性。

表4 數(shù)值模擬與分析模型邊坡穩(wěn)定系數(shù)對比

3 討 論

3.1 根系固土時間效應及其數(shù)值分析

本文與?wita?a等[39-42]的研究結果一致，都表明了植物根系對邊坡穩(wěn)定性具有加強作用。張興玲等[43]發(fā)現(xiàn)植物根系的生長指標與生長時間之間呈現(xiàn)顯著相關關系，并通過對生長期為1～4 a的4種不同灌木進行現(xiàn)場拉拔試驗，指出根系平均抗拔力隨生長時間的增加呈顯著增加趨勢。Zhang等[44]通過進行原位不排水抗剪強度測試，指出根系數(shù)量、單位體積根截面積與土體抗剪強度增強值均極顯著相關，并且根土復合體抗剪強度隨樹齡的增大而增大。Valizade等[45]通過直剪試驗，發(fā)現(xiàn)隨著植物齡期的增長，植物根系與周圍土壤顆粒之間的結合更好，從而增強了土體的抗剪強度。本文在不同生長期植物根系對邊坡穩(wěn)定性的影響方面，結論與以上研究一致。

本文采用數(shù)值模擬方法，研究了不同生長期紫花苜蓿根系對黃土邊坡的護坡效果，驗證了紫花苜蓿隨著生長期的延長，坡面土體的強度逐漸提高，邊坡穩(wěn)定系數(shù)逐漸增大，這與已有文獻[46-47]結論一致。將根土復合體層看做一個整體，通過測定根土復合體試樣的抗剪強度指標，這與將根系與土分別考慮進行數(shù)值模擬計算有所不同，但所的結論一致。通過分別測定多花木蘭根系與土體的物理力學參數(shù)，嵇曉雷等[46]采用點對點錨桿單元，模擬不同生長期多花木蘭根系；鄭明新等[47]采用桁架單元模擬多花木蘭根系。本文發(fā)現(xiàn)了紫花苜蓿生長的前5個月為關鍵期，當紫花苜蓿生長時間達到150 d時，紫花苜蓿對降雨條件下黃土邊坡坡面侵蝕有明顯的抑制作用。

總之，植物根系對邊坡穩(wěn)定性具有加強作用；隨著植物生長期的延長，根系固坡的作用也在逐步加強；數(shù)值模擬方法能夠有效分析不同生長期植物根系對邊坡淺層破壞的防護效果。今后將進一步開展模型試驗或現(xiàn)場監(jiān)測分析不同生長期植物根系對邊坡淺層破壞的防護效果。

3.2 數(shù)值模型中的參數(shù)取值合理性分析

本文通過三軸壓縮試驗測定了含紫花苜蓿根系黃土試樣的黏聚力、內摩擦角、彈性模量和泊松比4個參數(shù)。由于植物根系的生長不可避免會存在上下部根密度、根直徑的差異，三軸試驗試樣取樣處的根土復合體強度對于埋深較大的根系分布范圍內土樣確實相對較高。本文考慮降雨強度的最不利因素，即邊坡表層1 m內完全飽和，同時在相對較高的土體強度下邊坡坡面已然都發(fā)生了破壞，在強度較低的埋深較大的根系分布范圍內發(fā)生破壞更容易。另外，本文模型計算是針對邊坡整體的穩(wěn)定性，文章研究內容重點關注的是邊坡表層1 m，相對整體邊坡的幾何尺寸，邊坡表層僅1 m幾何尺寸相差較大，因此對于整體邊坡的穩(wěn)定性分析并沒有太大影響。若再將邊坡表層1 m進行細分，不僅會增加數(shù)值模擬模型建立的復雜度，數(shù)值模擬模型的建立并非越精細越好，適當簡化模型能夠更加清楚地反映計算結果的規(guī)律[48-49]，并且節(jié)約了時間成本。

泊松比是側向應變增量與豎向應變增量的比值。曾超峰等[50]通過開展基坑抽水模型試驗與三維數(shù)值分析，引用其他文獻的材料參數(shù)，將圍護結構、內隔墻、降水井的泊松比統(tǒng)一設定為0.2，探究內隔墻對于限制開挖前抽水引發(fā)基坑變形的有效性。常成等[51]開發(fā)了一個概念性的熱-孔隙流體-應力應變耦合的二維平面應變有限元模型，將圍巖、侵入巖體和花崗巖基底的泊松比統(tǒng)一定為0.25，并模擬了斑巖成礦系統(tǒng)的溫度、應力、應變、孔隙流體壓力的演化及巖漿侵入體周圍的破裂形成過程。楊衛(wèi)波等[52]為了獲得熱力耦合作用下相變混凝土能量樁的熱-力學特性，建立了其三維數(shù)值模型，引用其他文獻的材料參數(shù)，將泊松比設定為定值，比較了傳統(tǒng)和相變混凝土能量樁熱-力學特性的差異。本文泊松比參數(shù)為三軸試驗實測獲得，隨著植物生長期的延長，植物根系逐漸生長發(fā)育，根系固土作用逐漸增強，試樣的側向應變增量逐漸減小，導致泊松比的逐漸減小。左學龍[53]采用GDS非飽和土試驗儀，針對四類不同含根量的根-土復合體試樣，分別進行了常規(guī)三軸壓縮排水試驗和等向固結試驗等，基于鄧肯-張模型和K-G模型理論，探索了含根量對根-土復合體力學特性的影響規(guī)律，發(fā)現(xiàn)隨著含根量的增大，根土復合體的泊松比減小。而其他論文多將泊松比設置為一定值，如鄭明新等[47]將不同生長期多花木蘭根系的泊松比參數(shù)統(tǒng)一設置為0.2，采用二維有限元法模擬試驗區(qū)不同生長期多花木蘭根系對邊坡穩(wěn)定性的增強作用。

3.3 黃土邊坡淺層植物防護的時間效應

中國西北黃土高原地處干旱-半干旱地區(qū)，地質生態(tài)環(huán)境十分脆弱，水土流失嚴重，滑坡、崩塌災害頻發(fā)，給當?shù)氐墓まr(nóng)業(yè)生產(chǎn)以及人民的生命財產(chǎn)安全帶來巨大的威脅[54]。植物護坡是恢復黃土邊坡生態(tài)、減少土壤侵蝕和防治黃土邊坡地質災害的有效手段[55]。

在開展植被護坡措施的同時，研究人員已經(jīng)關注到植物護坡的時間效應問題，如張興玲等[43,56]為研究植物護坡的時間效應，對生長期為1～4 a的四種不同灌木進行現(xiàn)場拉拔試驗，指出根系平均抗拔力隨生長時間的增加呈顯著增加趨勢；陳潮[57]從時間尺度，對護坡灌木根系與土體之間相互作用力學機理展開研究，發(fā)現(xiàn)邊坡淺層穩(wěn)定性隨著時間的推移穩(wěn)步提升，且植被防護效果越來越好；劉治興等[58]通過對素土及植物生長期61、95、125 d的根土復合體抗剪強度進行現(xiàn)場測定和物理模型預測，發(fā)現(xiàn)抗剪強度隨植物生長階段呈逐漸增大趨勢。

當紫花苜蓿從120 d生長至150 d時，邊坡穩(wěn)定系數(shù)從0.695增大至1.367。究其原因，首先工況5在坡面處分布有3層，分別為飽和根土復合體層、飽和素黃土層和天然素黃土層，工況6僅有飽和根土復合體層和天然素黃土層。當邊坡處于工況5時，由于飽和素黃土的強度較低，在整個黃土邊坡中成為了軟弱層帶，軟弱層帶通常是邊（斜）坡失穩(wěn)的控滑帶[59]，因此其邊坡穩(wěn)定系數(shù)較低。當處于工況6時，植物根系生長至1 m，根系的加固作用大大增強了邊坡軟弱層帶的強度，提高了黃土邊坡的穩(wěn)定系數(shù)。

黃土地區(qū)的雨季一般在7月–9月，雨季坡面侵蝕現(xiàn)象多發(fā)。根據(jù)本文的研究結果表明：當生長時間達到150 d時，草本植物根系能夠有效防治黃土坡面侵蝕，并保持黃土邊坡穩(wěn)定性。并且隨著植物生長期的延長，植物根系對邊坡土體的固土效果更加顯著[46]。因此，草本植物生長的前5個月為關鍵期，假定對裸露黃土邊坡于春季（3月或4月）進行播種，則同年雨季來臨之時，坡面植被能夠起到防治降雨條件下的坡面侵蝕的作用，若為其他時間播種，在前5個月時需要進行一定的人為保護，如鋪麥草、塑料布等，對裸露坡面進行保護。坡面侵蝕普遍發(fā)育面積大，采用傳統(tǒng)工程措施進行防治需要龐大的資金支持且容易造成資源浪費，相比而言造價低廉的植被種植技術具有廣闊的推廣前景，這一技術在黃土邊坡淺層破壞防治的同時，還兼顧了景觀效應和生態(tài)可持續(xù)發(fā)展。

4 結 論

本文考慮降雨入滲條件并以含不同生長期紫花苜蓿的黃土邊坡為例，分別采用無限邊坡模型和數(shù)值模擬兩種方法進行了穩(wěn)定性分析，探究了紫花苜蓿對黃土邊坡淺層破壞防護的時間效應。本文得出以下結論：

1）不同工況下黃土邊坡淺層破壞過程都是首先在邊坡坡面以下1 m處，即降水最大入滲深度處，形成了平行于坡面的剪切屈服塑性區(qū)，之后邊坡坡面中上部位置的少量土體開始因拉伸而產(chǎn)生剪切屈服，并隨后剪切屈服區(qū)域沿邊坡坡面逐漸向邊坡坡肩方向移動，且邊坡的最大位移都出現(xiàn)在邊坡坡腳處。

2）黃土邊坡在天然條件下處于穩(wěn)定狀態(tài)，潛在滑動面在坡體內部較深位置，降雨作用使得邊坡潛在滑動面位置由坡體內部轉移到坡體降雨最大入滲深度處，在植物根系的加固作用下，坡面位置的穩(wěn)定性提高，此時邊坡潛在滑動面又轉為內部。

3）草本植物隨著生長期的延長，坡面土體的強度逐漸提高。數(shù)值模擬強度折減法計算的邊坡穩(wěn)定系數(shù)隨著植被生長不斷增大，表明植被對降雨條件下黃土邊坡坡面侵蝕有明顯的抑制作用，在坡面上種植植物能夠有效的防止黃土邊坡在雨季降雨事件時發(fā)生黃土邊坡淺層破壞。

4）通過對比無限邊坡模型與數(shù)值模擬方法計算的邊坡穩(wěn)定系數(shù)，發(fā)現(xiàn)隨著植物生長期的延長，采用數(shù)值模擬與無限邊坡模型兩種方法計算的邊坡穩(wěn)定系數(shù)，總體變化趨勢一致，表現(xiàn)出草本植物對黃土邊坡淺層破壞防護具有時間效應；黃土邊坡在天然狀態(tài)和降雨條件下植物生長期為150 d時穩(wěn)定系數(shù)大于1，處于穩(wěn)定狀態(tài)，在降雨條件下植物生長期為0～120 d時穩(wěn)定系數(shù)小于1，都處于欠穩(wěn)定狀態(tài)。

5）當草本植物生長時間達到150 d時，植物根系生長能夠有效防治黃土坡面侵蝕，并提高黃土邊坡穩(wěn)定性，因此草本植物生長的前5個月為關鍵期。

[1] 王蘭民，蒲小武，吳志堅，等. 地震和降雨耦合作用下黃土邊坡失穩(wěn)滑移的振動臺試驗研究[J]. 巖石力學與工程學報，2017，36(增刊2)：3873-3883.

Wang Lanmin, Pu Xiaowu, Wu Zhijian, et al. The shaking table test of the instability sliding of loess slope under the coupling effects of earthquake and rainfall[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(Supp.2): 3873-3883. (in Chinese with English abstract)

[2] 孫軍杰，王蘭民，龍鵬偉，等. 地震與降雨耦合作用下區(qū)域滑坡災害評價方法[J]. 巖石力學與工程學報，2011，30(4)：752-760.

Sun Junjie, Wang Lanmin, Long Pengwei, et al. An assessment method for regional susceptibility of landslides under coupling condition of earthquake and rainfall[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(4): 752-760. (in Chinese with English abstract)

[3] 王鮮，倪萬魁，劉海松，等. 黃土邊坡坡面降雨沖刷破壞機制模擬試驗[J]. 中國地質災害與防治學報，2015，26(3)：41-45.

Wang Xian, Ni Wankui, Liu Haisong, et al. Simulation experiment on loess slope rainfall erosion mechanism[J]. The Chinese Journal of Geological Hazard and Control, 2015, 26(3): 41-45. (in Chinese with English abstract)

[4] Tu X B, Kwong A K L, Dai F C, et al. Field monitoring of rainfall infiltration in a loess slope and analysis of failure mechanism of rainfall-induced landslides[J]. Engineering Geology, 2009, 105(1/2): 134-150.

[5] 張茂省，李同錄. 黃土滑坡誘發(fā)因素及其形成機理研究[J]. 工程地質學報，2011，19(4)：530-540.

Zhang Maosheng, Li Tonglu. Triggering factors and forming mechanism of loess landslides[J]. Journal of Engineering Geology, 2011, 19(4): 530-540. (in Chinese with English abstract)

[6] 董凌勃，海旭瑩，汪曉珍，等. 黃土高原退耕還草地植物群落動態(tài)對生態(tài)系統(tǒng)碳儲量的影響[J]. 生態(tài)學報，2020，40(23)：8559-8569.

Dong Linbo, Hai Xuying, Wang Xiaozhen, et al. Effects of plant community dynamics on ecosystem carbon stocks since returning farmlands to grasslands on the Loess Plateau[J]. Acta Ecologica Sinica, 2020, 40(23): 8559-8569. (in Chinese with English abstract)

[7] 邱海軍，崔鵬，曹明明，等. 基于最大熵原理的黃土丘陵區(qū)地質災害規(guī)模頻率分布研究[J]. 巖土力學，2014，35(12)：3541-3549，3555.

Qiu Haijun, Cui Peng, Cao Mingming, et al. Distribution of frequency-size of geological disaster base on principle of maximum entropy in loess plateau[J]. Rock and Soil Mechanics, 2014, 35(12): 3541-3549, 3555. (in Chinese with English abstract)

[8] 殷躍平，胡時友，石勝偉，等. 滑坡防治技術指南[M]. 北京：地質出版社，2018：267-268.

[9] 姚仰平，張丙印，朱俊高. 土的基本特性、本構關系及數(shù)值模擬研究綜述[J]. 土木工程學報，2012，45(3)：127-150.

Yao Yangping, Zhang Bingyin, Zhu Jungao. Behaviors, constitutive models and numerical simulation of soils[J]. China Civil Engineering Journal, 2012, 45(3): 127-150. (in Chinese with English abstract)

[10] 張丙印，溫彥鋒，朱本珍，等. 土工構筑物和邊坡工程發(fā)展綜述：作用機理與數(shù)值模擬方法[J]. 土木工程學報，2016，49(8)：1-15，35.

Zhang Bingyin, Wen Yanfeng, Zhu Benzhen, et al. Action mechanism and numerical simulation methods of soil structures and slope engineering: the state-of-the-art[J]. China Civil Engineering Journal, 2016, 49(8): 1-15, 35. (in Chinese with English abstract)

[11] 姚成成. 土質淺層邊坡的草本植物防護加固效應分析[D]. 成都：成都理工大學，2017.

Yao Chengcheng. Analysis of Herb Protection and Reinforcement Effect of Soil Shallow Slope[D]. Chengdu: Chengdu University of Technology, 2017. (in Chinese with English abstract)

[12] 周群華，鄧衛(wèi)東. 植物根系固坡的有限元數(shù)值模擬分析[J]. 公路，2007(12)：132-136.

Zhou Qunhua, Deng Weidong. Finite element numerical simulation analysis of slope reinforcement by plant roots[J]. Highway, 2007(12): 132-136. (in Chinese with English abstract)

[13] Temgoua A G T, Kokutse N K, Kavazovic Z, et al. 3D model applied to analyze the mechanical stability of real-world forested hillslopes prone to landslides[J]. Ecological Engineering, 2017, 106: 609-619.

[14] Li Y P, Wang Y Q, Ma C, et al. Influence of the spatial layout of plant roots on slope stability[J]. Ecological Engineering, 2016, 91: 477-486.

[15] 黃建坤，王學林，及金楠，等. 基于漸近均勻化理論的黃土高原草本植物固土效果模擬[J]. 農(nóng)業(yè)工程學報，2020，36(9)：168-176.

Huang Jiankun, Wang Xuelin, Ji Jinnan, et al. Numerical simulation of root reinforcement for herbs in Loess Plateau based on asymptotic homogenization theory[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(9): 168-176. (in Chinese with English abstract)

[16] Fattet M, Fu Y, Ghestem M, et al. Effects of vegetation type on soil resistance to erosion: Relationship between aggregate stability and shear strength[J]. Catena, 2011, 87(1): 60-69.

[17] Stokes A, Douglas G B, Fourcaud T, et al. Ecological mitigation of hillslope instability: Ten key issues facing researchers and practitioners[J]. Plant Soil, 2014, 377(1): 1-23.

[18] 楊吉華，張光燦，劉霞，等. 紫花苜蓿保持水土效益的研究[J]. 土壤侵蝕與水土保持學報，1997，3(2)：91-96.

Yang Jihua, Zhang Guangcan, Liu Xia, et al. Study on benefits of alfalfa conservation soil and water[J]. Journal of Soil Erosion and Soil and Water Conservation, 1997, 3(2): 91-96. (in Chinese with English abstract)

[19] 趙祥，朱慧森. 牧草科學種植與利用[M]. 北京：中國社會出版社，2006.

[20] 李廣信. 高等土力學[M]. 北京：清華大學出版社，2004.

[21] 雷小芹，劉恩龍，何思明，等. 非飽和堆積土邊坡降雨-滲流潛蝕耦合過程模擬[J]. 工程科學與技術，2021，53(2)：28-37.

Lei Xiaoqin, Liu Enlong, He Siming, et al. Modelling of coupled process of rainfall infiltration and internal erosion within unsaturated deposited soil slopes[J]. Advanced Engineering Sciences, 2021, 53(2): 28-37. (in Chinese with English abstract)

[22] 金保明，林鵬，李光敦. 基于運動波壤中流理論與無限邊坡穩(wěn)定分析理論的滑坡預測模型[J]. 工程科學與技術，2021，53(4)：140-148.

Jin Baoming, Lin Peng, Lee Kwantun. Integrating kinematic subsurface-wave approximation with infinite-slope stability analysis for landslide predictions[J]. Advanced Engineering Sciences, 2021, 53(4): 140-148. (in Chinese with English abstract)

[23] 吳宏偉. 大氣-植被-土體相互作用:理論與機理[J]. 巖土工程學報，2017，39(1)：1-47.

Wu Hongwei. Atmosphere-plant-soil interactions: Theories and mechanisms[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(1): 1-47. (in Chinese with English abstract)

[24] Arnone E, Caracciolo D, Noto L V, et al. Modeling the hydrological and mechanical effect of roots on shallow landslides[J]. Water Resources Research, 2016, 52(11): 8590-8612.

[25] Hayati E, Abdi E, Saravi M M, et al. Soil water dynamics under different forest vegetation cover: Implications for hillslope stability[J]. Earth Surface Processes and Landforms, 2018, 43(10): 2106-2120.

[26] Leung F T Y , Yan W M, Hau B C H, et al. Root systems of native shrubs and trees in Hong Kong and their effects on enhancing slope stability[J]. Catena, 2015, 125: 102-110.

[27] Waldron L J. The shear resistance of root-permeated homogeneous and stratified soil[J]. Soil Science Society of America Journal, 1977, 41(5): 843-849

[28] Wu T H, McKinnell W P, Swanston D N. Strength of tree roots and landslides on Price of Wales Island, Alaska[J]. Canadian Geotechnical Journal, 1979, 16(1): 19-33.

[29] 李國良，邵生俊，靳寶成，等. 黃土隧道地基的濕陷性問題研究[J]. 鐵道工程學報，2015，32(12)：12-16，50.

Li Guoliang, Shao Shengjun, Jin Baocheng, et al. Research on the Problems of collapsibility of loess tunnel foundation[J]. Journal of Railway Engineering Society, 2015, 32(12): 12-16, 50. (in Chinese with English abstract)

[30] 邵生俊，楊春鳴，焦陽陽，等. 濕陷性黃土隧道的工程性質分析[J]. 巖土工程學報，2013，35(9)：1580-1590.

Shao Shengjun, Yang Chunming, Jiao Yangyang, et al. Engineering properties of collapsible loess tunnel[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(9): 1580-1590. (in Chinese with English abstract)

[31] 白文明，左強，李保國. 烏蘭布和沙區(qū)紫花苜蓿根系吸水模型[J]. 植物生態(tài)學報，2001，25(4)：431-437.

Bai Wenming, Zuo Qiang, Li Baoguo. A simulation model for water uptake by alfalfa roots in the Wulanbuhe sandy desert[J]. Acta Phytoecologica Sinica, 2001, 25(4): 431-437. (in Chinese with English abstract)

[32] Tu X B, Kwong A K L, Dai F C, et al. Field monitoring of rainfall infiltration in a loess slope and analysis of failure mechanism of rainfall-induced landslides[J]. Engineering Geology, 2009, 105(1/2): 134-150.

[33] Mao Z J, Wang X K, An N, et al. Water leakage susceptible areas in loess multi-arch tunnel operation under the lateral recharge conditions[J]. Environmental Earth Sciences, 2020, 79(15): 1-32.

[34] Li Z Q, Oyediran I A, Tang C, et al. FEM application to loess slope excavation and support: Case study of Dong Loutian coal bunker, Shuozhou, China[J]. Bulletin of Engineering Geology and the Environment, 2014, 73(4): 1013-1023.

[35] 趙秋. 北方低產(chǎn)土壤實用綠肥作物栽培與利用[M]. 天津：天津科技翻譯出版公司，2010.

[36] 張魯渝，鄭穎人，趙尚毅，等. 有限元強度折減系數(shù)法計算土坡穩(wěn)定安全系數(shù)的精度研究[J]. 水利學報，2003(1)：21-27.

Zhang Luyu, Zheng Yingren, Zhao Shangyi, et al. The feasibility study of strength-reduction method with FEM for calculating safety factors of soil slope stability[J]. Journal of Hydraulic Engineering, 2003(1): 21-27. (in Chinese with English abstract)

[37] 趙勇，李國良，喻渝. 黃土隧道工程[M]. 北京：中國鐵道出版社，2011：32.

[38] Yin Y P, Li B, Wang W P. Dynamic analysis of the stabilized Wangjiayan landslide in the Wenchuan Ms 8.0 earthquake and aftershocks[J]. Landslides, 2015, 12(3): 537-547.

[39] ?wita?a B M, Wu W. Numerical modelling of rainfall- induced instability of vegetated slopes[J]. Géotechnique, 2018, 68(6): 481-491.

[40] ?wita?a B M. Numerical simulations of triaxial tests on soil-root composites and extension to practical problem: rainfall-induced landslide[J]. International Journal of Geomechanics, 2020, 20(11): 04020206.

[41] Park J, Kim I, Kang J K. Root Reinforcement effect on cover slopes of solid waste landfill in soil bioengineering[J]. Sustainability, 2021, 13(7): 3991.

[42] Eab K H, Likitlersuang S, Takahashi A. Laboratory and modelling investigation of root-reinforced system for slope stabilization[J]. Soils and Foundations, 2015, 55(5): 1270-1281.

[43] 張興玲，胡夏嵩，李國榮，等. 青藏高原東北部黃土區(qū)灌木幼林根系護坡的時間效應[J]. 農(nóng)業(yè)工程學報，2012，28(4)：136-141.

Zhang Xingling, Hu Xiasong, Li Guorong, et al. Time effect of young shrub roots on slope protection of loess area in Northeast Qinghai-Tibetan plateau[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(4): 136-141. (in Chinese with English abstract)

[44] Zhang D, Cheng J H, Liu Y, et al. Spatio-temporal dynamic architecture of living brush mattress:root system and soil shear strength in riverbanks[J]. Forests, 2018, 9(8): 493.

[45] Valizade N, Tabarsa A. Laboratory investigation of plant root reinforcement on the mechanical behaviour and collapse potential of loess soil[J]. European Journal of Environmental and Civil Engineering, 2020, 26(4): 1475-1491.

[46] 嵇曉雷，尹鵬. 不同生長時間根系對邊坡穩(wěn)定性的影響[J]. 中南林業(yè)科技大學學報，2020，40(12)：144-150，178.

Ji Xiaolei, Yin Peng. Influence of different root growth stages on the slope anti-erosion performance[J]. Journal of Central South University of Forestry & Technology, 2020, 40(12): 144-150, 178. (in Chinese with English abstract)

[47] 鄭明新，黃鋼，彭晶. 不同生長期多花木蘭根系抗拉拔特性及其根系邊坡的穩(wěn)定性[J]. 農(nóng)業(yè)工程學報，2018，34(20)：175-182.

Zheng Mingxin, Huang Gang, Peng Jing. Tensile-pullout properties of roots of Magnolia multiflora in different growth stages and stability of slope with its root[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(20): 175-182. (in Chinese with English abstract)

[48] 朱家奇，朱晟. 單元網(wǎng)格密度及畸變度對面板壩有限元計算精度的影響[J]. 水電能源科學，2022，40(6)：88-92.

Zhu Jiaqi, Zhu Sheng. Influence of element mesh density and distortion degree on finite element calculation accuracy of face rockfill dams[J]. Water Resources and Power, 2022, 40(6): 88-92. (in Chinese with English abstract)

[49] 汪必升，李毅波，廖雅詩，等. 基于擴展有限元模型的動態(tài)應力強度因子計算[J]. 中國機械工程，2019，30(11)：1294-1301.

Wang Bisheng, Li Yibo, Liao Yashi, et al. Calculation of dynamic stress intensity factors based on XFFM model[J]. China Mechanical Engineering, 2019, 30(11): 1294-1301. (in Chinese with English abstract)

[50] 曾超峰，王碩，宋偉煒，等. 內隔墻對開挖前抽水引發(fā)軟土區(qū)地鐵深基坑變形的控制效果[J]. 巖石力學與工程學報，2021，40(6)：1277-1286.

Zeng Chaofeng, Wang Shuo, Song Weiwei, et al. Control effect of cross wall on metro foundation pit deformation induced by pre-excavation dewatering in soft soils[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(6): 1277-1286. (in Chinese with English abstract)

[51] 常成，羅綱. 斑巖礦床侵入體頂部破裂系統(tǒng)形成的力學機制：多場耦合數(shù)值模擬的啟示[J]. 地球物理學報，2022，65(8)：3006-3024.

Chang Cheng, Luo Gang. The fracture forming mechanism at the top of the intrusion body in porphyry deposit: The insights from a multi-physical field coupled numerical model[J]. Chinese Journal of Geophysics, 2022, 65(8): 3006-3024. (in Chinese with English abstract)

[52] 楊衛(wèi)波，楊彬彬，汪峰. 相變混凝土能量樁熱-力學特性的數(shù)值模擬與試驗驗證[J]. 農(nóng)業(yè)工程學報，2021，37(2)：268-277.

Yang Weibo, Yang Binbin, Wang Feng. Numerical simulation and experimental validation of the thermo-mechanical characteristics of phase change concrete energy pile[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(2): 268-277. (in Chinese with English abstract)

[53] 左學龍. 基于彈性本構關系的根-土復合體力學特性試驗研究[D]. 長沙：中南林業(yè)科技大學，2016.

Zuo Xuelong. Experimental Study on Mechanical Properties of Root Soil Composite Based on Elastic Constitutive Relationship[D]. Changsha: Central South University of Forestry Science and Technology, 2016. (in Chinese with English abstract)

[54] 彭建兵，王啟耀，門玉明，等. 黃土高原滑坡災害[M]. 北京：科學出版社，2019.

[55] Li Q, Liu G B, Zhang Z, et al. Effect of root architecture on structural stability and erodibility of topsoils during concentrated flow in hilly Loess Plateau[J]. Chinese Geographical Science, 2015, 25(6): 757-764.

[56] 張興玲，胡夏嵩，李國榮，等. 寒旱環(huán)境草本植物根系護坡的時間尺度效應[J]. 水文地質工程地質，2009，36(4)：117-120.

Zhang Xingling, Hu Xiasong, Li Guorong, et al. Time scale effect of herb roots on slope protection in a cold, arid and semi-arid environment[J]. Hydrogeology and Engineering Geology, 2009, 36(4): 117-120. (in Chinese with English abstract)

[57] 陳潮. 護坡灌木根—土力學作用時間尺度效應研究[D]. 成都：西南交通大學，2017.

Chen Chao. Study on Time Scale Effect of Shrub Root Soil Mechanical Action on Slope Protection[D]. Chengdu: Southwest Jiaotong University, 2017. (in Chinese with English abstract)

[58] 劉治興，楊建英，楊陽，等. 高速公路不同植物防護邊坡根土復合體抗剪能力研究[J]. 生態(tài)環(huán)境學報，2015，24(4)：631-637.

Liu Zhixing, Yang Jianying, Yang Yang, et al. Research on the root-soil composite anti-shearing strength on expressway slope in different plant protection measures[J]. Ecology and Environmental Sciences, 2015, 24(4): 631-637. (in Chinese with English abstract)

[59] 崔圣華，裴向軍，黃潤秋. 大光包滑坡啟動機制：強震過程滑帶非協(xié)調變形與巖體動力致?lián)p[J]. 巖石力學與工程學報，2019，38(2)：237-253.

Cui Shenghua, Pei Xiangjun, Huang Runqiu. An initiation model of DGB landslide: Non-coordinated deformation inducing rock damage in sliding zone during strong seismic shaking[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38 (2): 237-253.

Numerical analysis of protection time effect on planting alfalfa in loess slope with shallow failure

Mao Zhengjun1,2,3, Zhang Jinge1, Bi Yinli1, Sun Weibo4, An Ning1

(1.,,710054,; 2.,,710054,; 3.,710054,; 4.,,710054,)

The special microstructure of loess has often posed a great risk on the slope in the loess areas. The loess structure can suddenly collapse, due to the complete loss of strength, particularly when the loess is wetted by water or subjected to a strong earthquake. The loess slope can usually be triggered by the slope surface erosion under the action of rainfall. The local failure of the slope body can often occur under heavy rainfall. There is a serious threat to the overall stability of the slope. For instance, hundreds of slope shallow failures have occurred in the Chinese Loess Plateau every year, as one of the regions with the thickest loess deposition in the world. The loess slope shallow failure can cause the loss and destruction of land sustainable productivity in the loess region. Fortunately, the slope protection with vegetation can be expected to alleviate the current slope shallow failure. A reasonable numerical simulation has been used to analyze the mechanism of root reinforcement soil, providing a strong theoretical basis for the application of slope protection with vegetation. This study aims to explore the protection time effect on planting alfalfa on the loess slope with shallow failure. The example was set as the loess slope with the alfalfa in different growth periods (0, 60, 90, 120, and 150 days). The rainfall infiltration was also considered in this case. A series of tests were carried out using the indoor planting of PVC pipe soil samples, including the quick direct shear tests, Consolidated Undrained (CU), and Unconsolidated Undrained (UU) triaxial tests. The physical and mechanical parameters were measured for the root-soil composite under different conditions. The slope stability was analyzed using the infinite slope model with the plant roots and the numerical simulation. The results show that three stages were divided for the loess slope under different conditions in the process of shallow failure. Firstly, a shear yield plastic zone parallel to the slope surface was formed at the maximum infiltration depth of precipitation. Secondly, there was a small amount of soil in the middle and upper position of the slope surface, due to the tension action. Thirdly, the shear yield zone gradually moved along the slope surface toward the slope shoulder. The loess slope was in a stable state under natural conditions. The potential sliding surface was at a deep position inside the slope. Therefore, the rainfall infiltration was attributed to the potential sliding surface of the slope transferring from the inside of the slope to the maximum infiltration depth of rainfall. The soil strength gradually increased in the shallow portion of the slope with the growth of vegetation on the slope. The potential sliding surface of the slope gradually shifted to the inside of the slope. The slope stability coefficient increased with the growth of vegetation, according to the strength reduction of numerical simulation. A comparison was made to determine the slope stability coefficient from the infinite slope model and the numerical simulation. The loess slope was in a stable state, when the plant roots grew to 1.0 m, either under natural or rainfall conditions, whereas, in an unstable state, when the plant roots grew to 0-0.8 m under rainfall conditions. The slope stability coefficient went down and then gradually increased with the growth of vegetation. The overall trend was all the same. The herbaceous plants were used to significantly inhibit the loess slope surface erosion under rainfall conditions, particularly with the extension of the growing period. Consequently, the herbaceous roots can be expected to improve the stability of the loess slope. Once the growth time of herbaceous plants reached 150 days, the loess slope surface erosion was effectively prevented for the better stability of the loess slope. Therefore, the first five months of herbaceous plant growth can be the key period to improving the stability of the loess slope in this case. The finding can provide the theoretical reference and application for the shallow failure protection of the loess slope in the sustainable development of the loess areas.

slope stability; numerical simulation; loess;alfalfa; shallow slope failure; time effect; infinite slope model

10.11975/j.issn.1002-6819.2022.15.008

S157.2

A

1002-6819(2022)-15-0072-12

毛正君，張瑾鴿，畢銀麗，等. 紫花苜蓿對黃土邊坡淺層破壞防護時間效應的數(shù)值分析[J]. 農(nóng)業(yè)工程學報，2022，38(15)：72-83.doi：10.11975/j.issn.1002-6819.2022.15.008 http://www.tcsae.org

Mao Zhengjun, Zhang Jinge, Bi Yinli, et al. Numerical analysis of protection time effect on planting alfalfa in loess slope with shallow failure[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(15): 72-83. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2022.15.008 http://www.tcsae.org

2022-06-15

2022-07-30

陜西省重點研發(fā)計劃項目（2020SF-379）；寧夏回族自治區(qū)重點研發(fā)計劃項目（2022BEG03059，2022BEG03052）

毛正君，博士，副教授，研究方向為地質環(huán)境保護與國土空間生態(tài)修復。Email：zhengjun_mao@163.com。