鄭明新，黃 鋼,，彭 晶

?

不同生長期多花木蘭根系抗拉拔特性及其根系邊坡的穩(wěn)定性

鄭明新1，黃 鋼1,2※，彭 晶2

（1. 華東交通大學(xué)土木建筑學(xué)院，南昌 330013； 2. 黃岡師范學(xué)院建筑工程學(xué)院，黃岡 438000）

為研究不同生長期灌木植物根系抗拉拔特性，并分析其對護(hù)坡效應(yīng)的影響，對生長期分別為3個月、6個月、9個月、12個月、15個月和18個月的典型護(hù)坡灌木植物多花木蘭根系進(jìn)行室內(nèi)抗拉試驗(yàn)和現(xiàn)場抗拔試驗(yàn)，并采用二維有限元法模擬試驗(yàn)區(qū)不同生長期多花木蘭根系邊坡穩(wěn)定性的增強(qiáng)作用。結(jié)果表明：受不同生長期的影響，多花木蘭單根的抗拉力隨著根徑的增加呈冪函數(shù)增大，抗拉強(qiáng)度隨著根徑的增加呈冪函數(shù)減小，單根直徑一定時(shí)，抗拉力和抗拉強(qiáng)度隨著生長期的增長而增大；受不同生長期的影響，多花木蘭植物根系的抗拔力隨著地徑、地下生物量和側(cè)根數(shù)的增加呈指數(shù)函數(shù)增大，多花木蘭根系的平均生物指標(biāo)和平均抗拔力的增大幅度隨著生長期的增長而減小，其中平均抗拔力增大幅度由102.78%減小到31.55%；邊坡的安全系數(shù)隨著植物生長期的增長而增大，等效塑性應(yīng)變和最大位移隨著植物生長期的增長而減小，其中最大等效塑性應(yīng)變較裸坡減小了24%～72%,最大水平位移的減少程度約為最大豎向位移的一半。

邊坡穩(wěn)定性；數(shù)值模擬；植物；根系；生長期；抗拉強(qiáng)度；抗拔力；生物指標(biāo)

0 引 言

隨著生態(tài)護(hù)坡技術(shù)的不斷發(fā)展，植物根系固土護(hù)坡效果已引起了國內(nèi)外學(xué)者的高度重視[1-2]。根系的護(hù)坡效果與植物根系抗拉拔特性緊密相關(guān)，根系的抗拉強(qiáng)度與土-根粘結(jié)強(qiáng)度是研究植被護(hù)坡力學(xué)機(jī)理的重要力學(xué)指標(biāo)。國內(nèi)外學(xué)者對植物根系抗拉拔方面試驗(yàn)進(jìn)行了大量研究[3]。在抗拉特性方面，李可等[4-9]開展單根抗拉試驗(yàn)，對草本和灌木植物不同根徑根系的抗拉力、抗拉強(qiáng)度進(jìn)行了研究，得出單根抗拉力隨根徑增大而減小的結(jié)論；李光瑩等[10]、Yan等[11]通過分析認(rèn)為，根系受土體運(yùn)動而產(chǎn)生的抗拉力是植物根系的主要固土力學(xué)效應(yīng)，故反映根系受力潛能的單根抗拉強(qiáng)度，可作為評判根系固土能力的一個有效指標(biāo)；Roering等[12]分析了根系化學(xué)成分的差異對根系抗拉強(qiáng)度的影響；Vergani等[13]發(fā)現(xiàn)植物根系的抗拉強(qiáng)度與環(huán)境因素有關(guān)；馮國建等[14]研究了根系抗拉強(qiáng)度和根系分布特征研究；萬娟等[15]探討了不同坡度邊坡灌木根系的抗拉強(qiáng)度與邊坡穩(wěn)定性。在抗拔特性方面，李紹才等[16]通過抗拉試驗(yàn)對巖石邊坡植被護(hù)坡工程坡面巖體根系互作的力學(xué)特性進(jìn)行了初步探討，研究了抗拉強(qiáng)度與復(fù)合體含根量等因子間的關(guān)系；Hu等[17]研究了根系纖維含量對抗拉力的影響；陳麗華等[18]通過對邊坡單株木本進(jìn)行拉拔試驗(yàn)，研究了根系的分布形態(tài)對抗拔能力影響；劉旭菲等[19-22]通過對根系的室內(nèi)外試驗(yàn)，研究了根-土摩擦力的影響因素；劉小光[23]研究了根系與土壤之間的錨固性能的影響因素；曹云生等[24]研究了根-土摩擦力的影響因素，根的抗拔力隨著根系根徑和埋深的增加而增大。這些研究成果中有關(guān)不同生長期的灌木植物根系的抗拉拔特性及護(hù)坡效果的研究較少。然而植物護(hù)坡是一個長期的、復(fù)雜的、動態(tài)的過程，根系護(hù)坡坡效果與植物的生長有很大關(guān)系[25-27]。為了探討生態(tài)護(hù)坡的可持續(xù)能力，需要對生態(tài)護(hù)坡中植物生長期對邊坡穩(wěn)定性影響進(jìn)行研究。

本文依據(jù)對試驗(yàn)區(qū)內(nèi)的灌木植物多花木蘭（）進(jìn)行室外現(xiàn)場抗拔試驗(yàn)和室內(nèi)抗拉試驗(yàn)，探討灌木植物多花木蘭根系抗拉力、抗拉強(qiáng)度和抗拔力與不同生長期關(guān)系，分析不同生長期的多花木蘭根系抗拉拔能力強(qiáng)弱變化。并采用二維有限元法，模擬和對比試驗(yàn)區(qū)黏土裸坡與不同生長期多花木蘭根系邊坡的安全系數(shù)和等效塑性區(qū)應(yīng)變、最大位移，從而評價(jià)試驗(yàn)區(qū)不同生長期多花木蘭根系護(hù)坡貢獻(xiàn)。這對于采用多花木蘭類護(hù)坡，實(shí)現(xiàn)基礎(chǔ)設(shè)施建設(shè)與生態(tài)環(huán)境保護(hù)之間的協(xié)調(diào)發(fā)展具有重要的學(xué)術(shù)價(jià)值和現(xiàn)實(shí)意義。

1 試驗(yàn)材料與方法

為保證護(hù)坡植物的生長適應(yīng)能力和水土保持能力，在黃岡師范學(xué)院周邊邊坡上建立試驗(yàn)區(qū)，選取優(yōu)良的護(hù)坡物種多花木蘭，種植在坡比為1∶1.5邊坡上。在該邊坡上，自植物發(fā)芽開始，每隔3個月（3個月、6個月、9個月、12個月、15個月和18個月）分別選取多花木蘭根系做現(xiàn)場抗拔試驗(yàn)和室內(nèi)抗拉試驗(yàn)。

抗拉試驗(yàn)方法：在試驗(yàn)區(qū)邊坡上隨機(jī)選取10株多花木蘭根系用水清洗掉根系上的泥土和側(cè)根，從主根上選取較直的10 cm段作為抗拉試驗(yàn)材料。為了保證單根新鮮、完整，把剪好單根放入保鮮袋里面儲存，并在24 h內(nèi)采用電子拉力儀進(jìn)行室內(nèi)抗拉試驗(yàn)，拉伸加載速度控制為0.02 mm/s，測得最大拉力,用游標(biāo)卡尺測量根斷裂處的平均根徑。

2 結(jié)果及分析

2.1 不同生長期多花木蘭根系抗拉特性分析

2.1.1 不同生長期多花木蘭根系抗拉特性分析

多花木蘭根系生長期、單根根徑與抗拉力和抗拉強(qiáng)度與的關(guān)系見圖1。

圖1 不同生長期根系的抗拉力、抗拉強(qiáng)度與根徑關(guān)系

由圖1可以看出，在單根生長期相同的情況下，單根抗拉力隨著根徑的增大而增加，兩者存在明顯的正向關(guān)系；單根抗拉強(qiáng)度隨著根徑的增大而減小，兩者存在明顯的負(fù)向關(guān)系。

相同生長期的單根根徑與抗拉力、抗拉強(qiáng)度之間具有較好的冪函數(shù)關(guān)系：

式中為抗拉力或抗拉強(qiáng)度，為單根根徑；，為擬合參數(shù)。不同生長時(shí)間根系擬合結(jié)果及相關(guān)系數(shù)見表1，相關(guān)系數(shù)均較高。

表1 不同生長期抗拉力、抗拉強(qiáng)度與根徑的擬合參數(shù)及相關(guān)系數(shù)

由表1分析可知，抗拉力與根徑擬合參數(shù)1在1.7左右變化，1隨著生長時(shí)間的增加而增大，從而同一根徑的根系抗拉力隨著生長時(shí)間的增加而增大；抗拉強(qiáng)度與根徑的擬合參數(shù)2在?0.45左右變化，2隨著生長期的增加而增大，從而可知，同一根徑一定時(shí)，根系抗拉強(qiáng)度隨著生長期的增加而增大，說明隨著植物的生長發(fā)育，根系的抗拉能力越來越強(qiáng)，這是因?yàn)殡S著植物生長發(fā)育，植物根系內(nèi)部的纖維含量逐漸增多，木質(zhì)化程度增加，提高根系的抗拉斷裂極限能力，這與劉治興[28]通過對不同邊坡防護(hù)根系的研究發(fā)現(xiàn)，隨著植物生長期的增長其根系平均抗拉力逐漸增加的結(jié)論一致。

2.1.2 根系抗拔特性

蔡曉松解釋說，有些客戶接到品牌訂單，覺得難以達(dá)到要求，公司就把其小樣拿過來，在實(shí)驗(yàn)室里試樣，得到工廠認(rèn)可后，再將配方和工藝交給客戶，并派技術(shù)人員前往指導(dǎo)。如果在售后中碰到問題，品牌商可以將產(chǎn)品拿到實(shí)驗(yàn)室，進(jìn)行分析，尋找解決方案，再到工廠實(shí)地解決。

根系抗拔力的大小反映了根-土相互作用的強(qiáng)度，因此抗拔力愈大，其根系愈有利于邊坡土體的穩(wěn)定性。根系的抗拔力與地徑、地下生物量、側(cè)根數(shù)等生物指標(biāo)有關(guān)[19]。

不同生長期多花木蘭抗拔力與根系生物指標(biāo)的關(guān)系見圖2。由圖2可以看出，隨著生長期的增長，地徑、地下生物量、側(cè)根數(shù)都增大，抗拔力增大。生長期3～6個月內(nèi)的生物指標(biāo)和抗拔力增幅度最大，之后生長期每增加3個月平均地徑增大幅度由102.78%減小到9.80%，平均地下生物量增大幅度由117.55%減小到19.70%，平均側(cè)根數(shù)增加幅度由104.35%減小到6.67%，平均抗拔力增大幅度由102.78%減小到31.55%。根系生物指標(biāo)的增加幅度隨著生長期的增長而減小。說明根系的生長發(fā)育程度與抗拔力之間存在明顯的正向關(guān)系，多花木蘭根系和土壤的粘結(jié)作用隨著生長不斷加強(qiáng)，但是生長半年后發(fā)育增快幅度有變緩趨勢。這是由于隨著根系的生長發(fā)育，根系中側(cè)根的數(shù)量、根系的根徑和延伸的深度隨之增加，根系對四周土體產(chǎn)生膨壓作用，使根和土之間的接觸更加緊密,增大了根和土之間的接觸面積，增強(qiáng)了根系和土壤之間的粘液、有機(jī)膠質(zhì)和摩擦等產(chǎn)生的粘結(jié)強(qiáng)度，從而根系的固持土壤的能力在不斷增強(qiáng)。

圖2 不同生長期根系的抗拔力與地徑、地下生物量和側(cè)根數(shù)的關(guān)系

圖3 邊坡幾何模型圖

圖4 邊坡網(wǎng)格劃分

2.2 多花木蘭根系邊坡穩(wěn)定性數(shù)值分析

目前，邊坡穩(wěn)定性的分析方法主要分為極限平衡法和有限元法兩大類。極限平衡法是基于摩爾-庫倫強(qiáng)度準(zhǔn)則，必需事先假定滑動面，計(jì)算式不考慮土體應(yīng)力應(yīng)變的關(guān)系，通過引入若干假定而轉(zhuǎn)化為靜定問題進(jìn)行計(jì)算。而強(qiáng)度折減有限元法可求解非線性問題，對巖土工程有較好的適用性，分析折減之后抗剪強(qiáng)度參數(shù)來計(jì)算邊坡安全系數(shù)，從而使得計(jì)算更加方便快捷[29-30]。以下使用ABAQUS軟件來模擬不同生長期多花木蘭根系邊坡的穩(wěn)定性。

2.2.1 數(shù)值模型建立

1）邊坡幾何模型與邊界條件

采用二維平面應(yīng)變模型來分析多花木蘭根系對邊坡土體穩(wěn)定作用。結(jié)合種植多花木蘭植物的試驗(yàn)區(qū)邊坡地形條件和根系密度（間隔0.5 m），建立簡化試驗(yàn)區(qū)邊坡平面幾何模型如圖3所示，模型左右邊界施加水平方向的位移約束，底部邊界施加水平和豎向約束。邊坡計(jì)算模型網(wǎng)格如圖4所示，劃分3 650個單元，2 146個節(jié)點(diǎn)，土體采用Mohr-Coulomb破壞準(zhǔn)則。

2）多花木蘭根系模型

多花木蘭根系屬于直根型灌木，垂直主根優(yōu)勢不明顯，斜向側(cè)根較為發(fā)達(dá)，與垂直主根之間的，角度在20°～65°之間。真實(shí)狀態(tài)下的多花木蘭根系在土體里面盤根錯節(jié)，數(shù)值分析時(shí)把每個多花木蘭根系全簡化成直線段，垂直主根和斜向側(cè)根均按等徑處理，側(cè)根的直徑簡化成主根地徑的一半，垂直主根和斜向側(cè)根之間的角度為45°，側(cè)根間距取0.1 m，對稱布置，最上面的側(cè)根與邊坡距離0.1 m。模擬計(jì)算時(shí)假定不考慮土體和多花木蘭根系界面發(fā)生錯動，把垂直主根和斜向側(cè)根使用約束命令嵌固土體里[31]，根系和土體能自動實(shí)現(xiàn)位移協(xié)調(diào)。多花木蘭根系為受拉強(qiáng)度遠(yuǎn)大于受壓強(qiáng)度柔性材料，采用桁架單元模擬多花木蘭根系。

2.2.2 數(shù)值模型工況及計(jì)算參數(shù)

本次模擬試驗(yàn)采用的飽和黏土裸坡以及植物生長期分別為3個月、6個月、9個月、12個月、15個月、18個月的飽和黏土邊坡的7個工況的模型，模擬工況如表2所示。

表2 模擬工況

飽和黏土邊坡土體的參數(shù)由土工試驗(yàn)測得，含水率為16.90%，彈性模量為11.23 MPa，黏聚力為12.68 kPa，內(nèi)摩擦角為21.35°，干密度為1.85 kg/m3，模型中根系的物理參數(shù)由試驗(yàn)測得，具體數(shù)據(jù)見表3。

表3 不同生長期根系參數(shù)表

2.2.3 計(jì)算結(jié)果及分析

1）邊坡安全系數(shù)分析

通過有限元計(jì)算得到各工況下邊坡安全系數(shù)分別為1.190、1.203、1.221、1.246、1.345、1.405和1.468。

由圖5可知，6種不同生長期多花木蘭的植物邊坡的安全系數(shù)大于裸坡，并隨著多花木蘭植物生長期增加，其安全系數(shù)逐漸增大。工況2～7的安全系數(shù)分別較工況1增幅1.09%、2.61%、4.71%、13.03%、18.07%和23.36%，環(huán)比增幅1.09%、1.50%、2.05%、7.95%、4.46%和4.48%。說明多花木蘭根的存在能過提高邊坡的穩(wěn)定性能,且隨著根系的生長發(fā)育，主根越長，側(cè)根越多，根系越發(fā)達(dá)，根系護(hù)坡能力越強(qiáng)，邊坡的安全系數(shù)越高。

2）邊坡塑性應(yīng)變

基于試驗(yàn)所測多花木蘭根系的增強(qiáng)邊坡穩(wěn)定性分析，各工況邊坡失穩(wěn)時(shí)的等效塑性應(yīng)變云圖和最大等效性應(yīng)變分別如圖6、表4所示。由此可知：

圖5 邊坡的安全系數(shù)

圖6 各工況下邊坡失穩(wěn)時(shí)等效塑性應(yīng)變云圖

表4 各工況下邊坡失穩(wěn)時(shí)最大等效塑性應(yīng)變

裸坡失穩(wěn)時(shí)最大等效塑性應(yīng)變出現(xiàn)在坡腳，為0.228 6，且主要集中坡面以下4 m之內(nèi)，其破壞屬于淺層滑坡。不同生長期的根系邊坡，最大等效塑性應(yīng)變比裸坡減少了24%～72%，并且隨著根系生長期的增加而減小，表明多花木蘭根系的存在能有效減少邊坡等效塑性應(yīng)變。

通過對裸坡與多花木蘭根系邊坡的滑動面對比發(fā)現(xiàn)，隨著邊坡根系生長期的增加，塑性區(qū)逐漸向邊坡的深層移動，塑性變形集中區(qū)域逐漸遠(yuǎn)離坡腳，并出現(xiàn)范圍越來越大隆起“延性”破壞。說明灌木多花木蘭植物根系能加固土體，控制邊坡淺層土體的變形，使土體變形區(qū)域向深層移動，增加淺層土體的穩(wěn)定性。

3）最大水平和豎向位移

基于試驗(yàn)所測多花木蘭根系的增強(qiáng)邊坡穩(wěn)定性分析，1～7工況的最大水平位移和豎向位移如表5所示。由表5可知，6種不同生長期多花木蘭根系均能有效減少邊坡的最大水平位移和豎向位移，并且隨著生長期的增大，最大水平位移和豎向位移的減少程度越來越高，最大水平位移的減少程度是最大豎向位移的一半左右，表明根系的生長發(fā)育減小邊坡的變形。這是由于剛度較大的多花木蘭根系與邊坡周邊土體的共同作用，增大邊坡整體剛度，約束了邊坡的水平位移和豎向位移，約束效應(yīng)并隨著根系的生長發(fā)育而增強(qiáng)，因根系剛度主要對水平向應(yīng)力起作用，所以植入根系后的邊坡水平位移減小幅度要大于豎向位移減小幅度。

表5 各工況下最大水平位移和豎向位移

3 討 論

本文研究發(fā)現(xiàn)多花木蘭植物的單根抗拉力隨著根徑的增大呈增大趨勢，抗拉強(qiáng)度隨著根徑的增大呈減小趨勢。這一研究結(jié)論與大多數(shù)研究者有關(guān)單根抗拉特性與根徑的關(guān)系的研究結(jié)果一致[4-6,9-10,15]；不同之處在于，不同生長期條件下植物單根的抗拉力增加、抗拉強(qiáng)度減少的函數(shù)不同，在直徑一定的情況下，生長期越長的根系抗拉力和抗拉強(qiáng)度越高。該結(jié)論需要結(jié)合植物根系在不同生長期的微觀結(jié)構(gòu)和纖維素含量進(jìn)行進(jìn)一步分析研究。由圖1可知多花木蘭植物的單根抗拉力和抗拉強(qiáng)度增大隨著生長期的增長而增大，這一結(jié)論與苑淑娟[31]關(guān)于2種生長期灌木檸條單根抗拉力和抗拉強(qiáng)度的研究結(jié)論一致，與劉艷琦等[7]關(guān)于灌木檸條、沙棘、沙柳的研究結(jié)論不一致。分析原因，可能與植物的種類、生長的立地條件不同有關(guān)，這些影響因素需要更進(jìn)一步的研究。

已有研究得到植物根系的抗拔力受到植物生物指標(biāo)的影響[22-24]。本文研究發(fā)現(xiàn)多花木蘭植物根系抗拔力在不同生長期條件下隨著生物指標(biāo)的增加呈增大趨勢，多花木蘭根系的平均生物指標(biāo)和平均抗拔力的增大幅度隨著生長期的增長而減小。這與李紹才等[16]研究發(fā)現(xiàn)巖石風(fēng)化條件一定時(shí)灌木抗拔力隨生物指標(biāo)的增加而增大的結(jié)論一致。但是，本文關(guān)于多花木蘭根系抗拔力與側(cè)根數(shù)呈冪函數(shù)關(guān)系的結(jié)論與李國榮等[32]研究寒旱地區(qū)灌木根系得出抗拔力呈線性相關(guān)結(jié)論不一致。分析原因，可能與植物的種類、土壤含水等水文條件不同有關(guān)，因?yàn)閺埮d玲等[22]研究發(fā)現(xiàn)根系抗拔力隨著土壤含水量有關(guān)，關(guān)于這些需要更進(jìn)一步的研究。

本文采用數(shù)值分析的研究發(fā)現(xiàn)結(jié)果與萬娟等[15]植物根系理論分析結(jié)果和瞿文斌等[20]、陳潮[30]有限元分析結(jié)果相一致。本文采用數(shù)值模擬中根系簡化模型反映了根系的特征和根土的相互作用，比瞿文斌等[20]更接近實(shí)際情況，但與實(shí)際還有一定差距，因根系在土中的分布特征、根系與土的相互作用非常復(fù)雜，這些將在以后做進(jìn)一步的研究。

4 結(jié) 論

1）相同生長期多花木蘭植物單根的抗拔力隨著直徑的增大而增大，抗拉強(qiáng)度隨著根徑的增大而減小，抗拉力、抗拉強(qiáng)度與根徑之間存在顯著冪函數(shù)關(guān)系；單根直徑一定時(shí)，抗拉力和抗拉強(qiáng)度增大隨著生長期的增長而增大，反映出根系的生長發(fā)育能提高根系的抗拉斷裂極限能力。

不同生長期多花木蘭植物根系的抗拔力隨著生物指標(biāo)的增大而增大，抗拔力與地徑、地下生物量和側(cè)根數(shù)三者之間存在顯著的指數(shù)函數(shù)關(guān)系。多花木蘭根系的平均生物指標(biāo)和平均抗拔力的增大幅度隨著生長期的增長而減小，其中平均抗拔力增大幅度由102.78%減小到31.55%，反映出根系的生長發(fā)育能提高根-土的粘結(jié)作用。

2）植物根系的存在能增強(qiáng)邊坡的穩(wěn)定性。邊坡的安全系數(shù)隨著植物生長期的增長而增大，等效塑性應(yīng)變和最大位移隨著植物生長期的增長而減小，其中最大等效塑性應(yīng)變較裸坡減小了24%～72%，最大水平位移的減少程度約為最大豎向位移的一半。

[1] 夏振堯，張倫，陳毅，等. 香根草根系與土壤接觸特性及抗拔模型研究[J]. 水生態(tài)學(xué)雜志，2016，37(4)：36－41.

Xia Zhenyao, Zhang Lun, Chen Yi, et al. Vetiver root-soil contact characteristics and pull-out model[J]. Journal of Hydroecology, 2016, 37(4): 36－41. (in Chinese with English abstract)

[2] 譚儒蛟，楊旭朝，胡瑞林. 反傾巖體邊坡變形機(jī)制與穩(wěn)定性評價(jià)研究綜述[J]. 巖土力學(xué)，2009，30(增刊2)：479－484.

Tan Rujiao, Yang Xuzhao, Hu Ruilin. Review of deformation mechanism and stability analysis of anti-dipped rock slopes[J]. Rock Mechanics, 2009, 30(Supp.2): 479－484. (in Chinese with English abstract)

[3] 楊旸，字淑慧，余建新，等. 植物根系固土機(jī)理及模型研究進(jìn)展[J].云南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2014，29(5)：759－765.

Yang Yang, Zi Shuhui, Yu Jianxin, et al.Advances in research on mechanism and models of plant roots reinforcement[J]. Journal of Yunnan Agricultural University, 2014, 29(5): 759－765.(in Chinese with English abstract)

[4] 李可，朱海麗，宋路，等. 青藏高原兩種典型植物根系抗拉特性與其微觀結(jié)構(gòu)上的關(guān)系[J]. 水土保持研究，2018，25(2)：240－249.

Li Ke, Zhu Haili, Song Lu, et al. Relationship between tensile properties and microstructure of two typical plant roots in the qinghai-tibei plateau[J]. Soil and Water Conservation Research, 2018, 25(2): 240－249. (in Chinese with English abstract)

[5] 朱海麗，胡夏嵩，毛小青，等. 青藏高原黃土區(qū)護(hù)坡灌木植物根系力學(xué)特性研究[J]. 巖石力學(xué)與工程學(xué)報(bào)，2008(增刊2)：3445－3452.

Zhu Haili, Hu Xiasong, Mao Xiaoqing, et al. Study on mechanical characteristcs of shrub roots for slope protection in loess area of tibetan plateau[J].Chinese Journal of Rock Mechanics and Engineering, 2008(Supp.2): 3445－3452. (in Chinese with English abstract)

[6] 陸桂紅，歐國強(qiáng)，楊順. 泥石流源區(qū)典型草本植物根系抗拉試驗(yàn)[J]. 山地學(xué)報(bào)，2014，32(6)：725－731.

Lu Guihong, Ou Guoqiang, Yang Shun. Root tensile test of typical grass in debris flow triggering area[J]. Mount Ain Research, 2014, 32(6): 725－731. (in Chinese with English abstract)

[7] 劉艷琦，格日樂，阿如旱，等. 2個生長時(shí)期5種植物單根抗拉力學(xué)特性比較[J]. 內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版，2017，38(6)：25－30.

Liu Yanqi, Gerile, Aruhan, et al.Comparison of single root tensile resistance characteristics of five plants in two growth periods[J]. Journal of Inner Mongolia Agricultural University: Natural Science Edition, 2017, 38(6): 25－30. (in Chinese with English abstract)

[8] 嵇曉雷. 基于植被根系分布形態(tài)的生態(tài)邊坡穩(wěn)定性研究[D]. 南京：南京林業(yè)大學(xué)，2013.

Ji Xiaolei. A Roots Distribution-based Study on the Stability of Ecological Slope[D]. Nanjing: Nanjing Forestry University, 2013. (in Chinese with English abstract)

[9] 嵇曉雷，楊平，王磊. 狗牙根根系抗拉強(qiáng)度試驗(yàn)[J].湖北農(nóng)業(yè)科學(xué)，2015，54(23)：5941－5948.

Ji Xiaolei, Yang Ping, Wang Lei. Experiments on root tensile strengths of bermuda grass[J].Hubei Agricultural Sciences, 2015, 54(23): 5941－5948. (in Chinese with English abstract)

[10] 李光瑩，杜江濤，余冬梅，等. 堿脅迫條件下草本植物根系力學(xué)強(qiáng)度試驗(yàn)研究[J]. 工程地質(zhì)學(xué)報(bào)，2016，24(4)：585－595.

Li Guangying, Fu Jiangtao, Yu Dongmei, et al.Mechanical strength characteristics of herbaceous plant roots under alkali stressing[J]. Journal of Engineering Geology, 2016, 24(4): 585－595. (in Chinese with English abstract)

[11] Yan Z, Song Y, Jiang P, et al. Mechanical analysis of interaction between plant roots and rock and soil mass in slope vegetation[J]. Applied Mathematics and Mechanics, 2010, 31(5): 617－622.

[12] Roering J J, Schmidt K M, Stock J D, et al. Shallow landsliding, root reinforcement, and the spatial distribution[J]. Canadian Geotechnical Journal, 2003, 40(2): 237－253.

[13] Vergani C, Chiaradia E, Bischetti G. Variability in the tensile resistance of roots in Alpine forest tree species[J]. Ecological Enginerring, 2012, 46(9): 43－56.

[14] 馮國建，王世通. 護(hù)坡植物根系分布特征及抗拉強(qiáng)度研究[J]. 重慶師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，30(2)：115－118.

Feng Guojian, Wang Shitong. Study on distributiaon characteristic and tensile strength of slope eco-engineering plant root[J]. Journal of Chongqing Normal University: Natural Science, 2013, 30(2): 115－118. (in Chinese with English abstract)

[15] 萬娟，肖衡林，何俊. 多花木蘭根系抗拉特性及邊坡穩(wěn)定性分析[J]. 華中科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2014(8)：110－116.

Wan Juan, Xiao Henglin, He Jun. Tensile Properties of magnolia multiflora’s roots and slope stability analysis[J]. Journal of Huazhong University of Science and Technology: Natural Science Edition, 2014(8): 110－116. (in Chinese with English abstract)

[16] 李紹才，孫海龍，楊志榮，等. 護(hù)坡植物根系與巖體相互作用的力學(xué)特性[J]. 巖石力學(xué)與工程學(xué)報(bào)，2006，25(10)：2051－2057.

Li Shaocai, Sun Hailong, Yang Zhirong, et al. Mechanical characteristics of interaction between root system of plants and rock for slope protection[J]. Chinese Journal of Rock and Engineering, 2006, 25(10): 2051－2057. (in Chinese with English abstract)

[17] Hu X, Brierley G, Zhu H, et al. An exploratory analysis of vegetation strategies to reduce shallow landslide activity on loess hillslopes, Northeast qinghai-tibet plateau[J]. Journal of Mountain Science, 2013, 10(4): 668－686.

[18] 陳麗華，余新曉，張東升，等. 整株林木垂向抗拉試驗(yàn)研究[J]. 資源科學(xué)，2004，26(增刊)：39－41.

Chen Linhua, Yu Xinxiao, Zhang Dongsheng, et al.Experimental study on vertically tensile strength of whole tree[J]. Resources Science, 2004, 26(Supp.): 39－41. (in Chinese with English abstract)

[19] 劉旭菲. 多因素作用植被對岸坡穩(wěn)定影響的試驗(yàn)及數(shù)值模擬研究[D]. 天津：天津大學(xué)，2015.

Liu Xufei. Experimental and numerical investigation of slope stability with multiple factors based on vegetation Protection[D]. Tianjin: Tianjin University, 2015. (in Chinese with English abstract)

[20] 瞿文斌，及金楠，陳麗華，等. 黃土高原植物根系增強(qiáng)土體抗剪強(qiáng)度的模型與試驗(yàn)研究[J]. 北京林業(yè)大學(xué)學(xué)報(bào)，2017，39(12)：79－86.

Qu Wenbin, Ji Jinnan, Chen Lihua, et al.Research on model and test of reinforcing shear strength by vegetation roots in the Loess Plateau of northern China[J]. Journal of Beijing Forestry University, 2017, 39(12): 79－87. (in Chinese with English abstract)

[21] 蔣坤云，陳麗華，蓋小剛，等. 華北護(hù)坡闊葉樹種根系抗拉性能與其微觀結(jié)構(gòu)的關(guān)系[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29(3)：115－123.

Jiang Kunyun, Chen Lihua, Gai Xiaogang, et al. Relationship be-tween tensile properties and microstruc-tures of three different broad-leaf tree roots in North China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(3): 115－123. (in Chinese with English abstract)

[22] 張興玲，胡夏嵩，毛小青，等. 青藏高原東北部黃土區(qū)護(hù)坡灌木檸條錦雞兒根系拉拔摩擦試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2011(增刊2)：3739－3745.

Zhang Xingling, Hu Xiasong, Mao Xiaoqing, et al. Research on pull-out friction test of shrub caraganak or shinskii roots for slope protection in loess area of noetheast qinghai-tibetan plateau[J]. Journal of Rock Mechanics and Engineering, 2011(Supp.2): 3739－3745. (in Chinese with English abstract)

[23] 劉小光. 林木根系與土壤摩擦錨固性能研究[D]. 北京：北京林業(yè)大學(xué)，2013.

Liu Xiaoguang. Study on Firction and Anchorage Characteristics Between Root System and Soil[D]. Beijing: Beijing Forestry University, 2013. (in Chinese with English abstract)

[24] 曹云生，陳麗華，劉小光. 植物根土界面摩擦力的影響因素研究[J]. 摩擦學(xué)學(xué)報(bào)，2014(5)：482－488.

Cao Yunsheng, Chen Lihua, Liu Xiaoguang.The lnfluence factors of plant root-soil interface friction[J]. Tribology, 2014(5): 482－488. (in Chinese with English abstract)

[25] 吳宏偉. 大氣–植被–土體相互作用：理論與機(jī)理[J]. 巖土工程學(xué)報(bào)，2017(1)：1－40.

Charles Wangwai. Atmosphere-plant-soil interactions: theories and mechanisms[J]. Chinese Journal of Geotechnical Engineering, 2017(1): 1－40. (in Chinese with English abstract)

[26] Mao Z, Bourrier F, Stokes A, et al. Three dimensional modelling of slope stability in heterogeneous montane forest ecosystems[J]. Engineering Division, 2014, 273: 11－22.

[27] Fan C C, Laiy F. Influence of the spatial layout of vegetation on the stability of slopes[J]. Plant and Soil, 2014, 377(1/2): 83－95.

[28] 劉治興. 高速公路邊坡植物不同生長期防護(hù)效果[D]. 北京：北京林業(yè)大學(xué)，2016.

Liu Zhixing. Study on the Protective Effect of Highway Slope Plant in Different Growth Stages[D]. Beijing: Beijing Forestry University, 2016. (in Chinese with English abstract)

[29] 閆俊維，王曙光，陳靜瑜. ABAQUS強(qiáng)度折減法在邊坡穩(wěn)定性分析中的應(yīng)用[J]. 西部交通科技，2015(8)：21－24.

Yan Junwei, Wang Shuguang, Chen Jun. Applications of ABAQUS strength reduction method in slope stability analysis[J]. Western China Communications Science & Technology, 2015(8): 21－24. (in Chinese with English abstract)

[30] 陳潮. 護(hù)坡灌木根-土力學(xué)作用時(shí)間尺度效應(yīng)研究[D]. 成都：西南交通大學(xué)，2017：1－54.

Chen Chao. Study on the Effect of Shrub Roots-soil in Bio-slope[D]. Chengdu: South Jiaotong University, 2017: 1－54. (in Chinese with English abstract)

[31] 苑淑娟. 4種植物單根抗拉力學(xué)特性的研究[D]. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2010.

Yuan Shujuan. Study on Mechanical Properties of Single Root Tensions of Four Plants[D]. Hohhot: Inner Mongolia Agricultural University, 2010. (in Chinese with English abstract)

[32] 李國榮，胡夏嵩，毛小青，等. 寒旱環(huán)境黃土區(qū)灌木根系護(hù)坡力學(xué)效應(yīng)研究[J]. 水文地質(zhì)工程地質(zhì)，2008(1)：94－97.

Li Guorong, Hu Xiasong, Mao Xiaoqing, et al. A study of the mechanical effects of shrub roots for slope protection in frigid and arid-semiarid loes area[J]. Hydrogeology and Engineering Geology, 2008(1): 94－97. (in Chinese with English abstract)

Tensile-pullout properties of roots ofin different growth stages and stability of slope with its root

Zheng Mingxin1, Huang Gang1,2※, Peng Jing2

(1.,330013,; 2.438000,)

To study the root tensile-pullout characteristics of shrubs in different growth periods and analyze their effects on slope protection. A indoor tensile test and a field pull-out test using the roots ofwhich is a typical slope-protection shrubs with a growth period of three months, six months, nine months, twelve months, fifteen months and eighteen months were carried out, respectively. The indoor tensile test and field pull-out test of’s roots were conducted at a constant displacement rate of 0.02 mm/s and 10 mm/min, respectively. Under the different growth periods, root mechanical effect indexes and physiological indexes were determined, respectively. Root mechanical effect indexes include the single root tensile resistance, single root tensile strength and root pull-resistance. The physiological indexes include the ground diameter, subterranean biomass and lateral root quality. The two-dimensional finite element method was used to simulate the enhancement of the slope stability in terms of the’s roots at different growth stages. The results are as follows: Under different growth periods, the tensile resistance of the’s single root increases with the increasing of root diameter according to a good power function, and the tensile strength decreases with the increasing of root diameter according to a good power function, when the single root diameter is constant, both the tensile resistance and the tensile strength increase with the growth period. Under the influence of different growth periods, the pull-resistance increases with the increasing of ground diameter, subterranean biomass and lateral root quantity according to a significant exponential function, the increase rate of the average biological index and the average pull-resistance of the’s roots are the largest in the growth period of three to six months. After six months, the average biological index and the average pull-resistance of the’s roots decrease with the growth period increasing, the increase rates of the ground diameter decrease from 102.78% to 9.8%, the increase rates of subterranean biomass decrease from 117.55% to 19.70% and the increase rates of lateral root quantity decrease from 104.35% to 6.67%, the increase rates of pull-resistance decrease from 102.78% to 31.55%.The safety factor of’s roots slope is greater than that of the bare slope, the increase rates of the safety factor are from 1.09% to 23.36% and the safety factor of the slope increases with the’s growth period increasing, which indicates that Magnolia magnolia’s roots can improve the stability of slope. The equivalent plastic strain and the maximum displacement decreases with the’s growth period increasing, among them, the maximum equivalent plastic strain is reduced by 24% to 72% compared with that of the bare slope. The’s roots can effectively reduce the slope’s maximum horizontal displacement and vertical displacement at different growth stages, the decreasing proportions of maximum horizontal displacement and maximum vertical displacement is increase with the growth period increasing and the reduction of the maximum horizontal displacement are about half of the maximum vertical displacement, it indicates that the growth and development of roots can effectively reduce the displacement of slopes. This research has important research value and practical significance for the coordinated development of infrastructure construction and ecological environment protection using shrub slopes such as.

slope stability; numerical simulation; vegetation; roots; growth period; tensile resistance; pull-resistance; biological indicator

10.11975/j.issn.1002-6819.2018.20.022

TU411.3

A

1002-6819(2018)-20-0175-08

2018-03-29

2018-06-05

國家自然科學(xué)基金項(xiàng)目（51568022）；江西省自然科學(xué)基金項(xiàng)目（20171BAB206056）

鄭明新，博士，博士生導(dǎo)師，長期從事道路邊坡穩(wěn)定與路基病害防治、環(huán)境巖土保護(hù)的教學(xué)和研究。Email：492001473@qq.com

黃 鋼，講師，博士生，主要從事環(huán)境巖土保護(hù)的研究。Email：627823405@qq.com

鄭明新，黃 鋼，彭 晶. 不同生長期多花木蘭根系抗拉拔特性及其根系邊坡的穩(wěn)定性[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(20)：175－182. doi：10.11975/j.issn.1002-6819.2018.20.022 http://www.tcsae.org

Zheng Mingxin, Huang Gang, Peng Jing.Tensile-pullout properties of roots ofin 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) doi：10.11975/j.issn.1002-6819.2018.20.022 http://www.tcsae.org