土體含水率與根系埋深對銀杏樹苗抗拉拔性能的影響

傅胤榕 仉文崗

摘 要：中國是受熱帶臺風災害影響最為嚴重的國家之一，由于強風暴雨作用引起的林木根系拔出及傾倒破壞致使交通中斷、車輛受損、人員傷亡等次生災害與日俱增。為進一步研究土體物理特性對林木抗拉拔性能的影響，采用自主設計的根系橫向加載試驗系統(tǒng)，通過對銀杏樹苗施加定速率橫向加載的拉拔試驗，模擬外力荷載下根系的傾倒過程。試驗設置了5組梯度含水率和4種埋深工況分析土體含水率、根系埋深對拔出過程中抗拉拔力的影響。研究表明：在粉質砂土中，黏聚力根系位移抗拉拔力曲線中抗拉拔力的變化可以概括為平緩上升、陡峭上升和緩慢上升達到峰值3個階段;最大抗拉拔力隨含水率提升呈現(xiàn)先升后降的總體趨勢，與土體黏聚力變化的總體趨勢相同，當含水率趨近飽和含水率時，最大抗拉拔力衰減顯著;當含水率一定時，根系的最大抗拉拔力隨著根系埋深的增加，在試驗設置的埋深區(qū)間內(nèi)呈線性趨勢增長。

關鍵詞：含水率;根系;抗拉拔力;埋深;模型試驗

中圖分類號：TU431 ? 文獻標志碼：A ? 文章編號：2096-6717（2019）05-0042-07

Abstract：China is one of the countries most seriously affected by the tropical typhoon disasters. Root pull-out and toppling damage caused by severe storms and rains occur frequently. Secondary disasters such as traffic disruption， vehicle damage and casualties are increasing. In order to investigate the effect of soil physical properties on the anti-overturning performance of trees， the pulling test of Ginkgo biloba seedlings under constant lateral loading rate was carried out through a self-designed lateral root loading test system to simulate the dumping process of roots under external load. Five sets of gradient water content and four buried depths were set to analyze the influence of soil moisture content and the root depth on the pull-out resistance during the pull-out process. The research indicates that in silty sand， the change of the anti-overturning force in the root displacement-drawing force curve can be summarized as three stages： gentle rise， steep rise and slowly rising to reach the peak. The anti-overturning force with the increase of moisture content showed the first rise. The overall trend of the post-fall is the same as the overall trend of the change of soil cohesion. When the moisture content approximates the saturated moisture content， the maximum anti-overturning force is significantly reduced. When the moisture content is constant， the maximum anti-overturning force of the root system increases linearly with the depth of the root system in the buried depth range set by the test.

Keywords：moisture content; root system; anti-overturning force; burial depth; model tests

? 隨著城市化進程的不斷推進和綠色生態(tài)文明的發(fā)展，城市綠化覆蓋率不斷提升，綠化樹在城市建設中起到凈化城區(qū)空氣、降低環(huán)境噪音、美化城市道路等多種作用。近年來，東南沿海省份陸續(xù)遭受“莫蘭蒂”、“蓮花”等臺風入境，遭受強風和暴雨，部分城市行道樹倒伏率接近30%，個別林木樹種倒伏率達到59%，總受損率達96%以上[1]。由于林木傾倒導致城市交通道路阻斷、車輛及公共基礎設施受損，甚至造成嚴重的人員傷亡，因此，強風、暴雨引起的林木傾倒破壞現(xiàn)象不容輕視。

學者們對林木倒伏和根土相互作用方面進行了大量研究，取得了一系列成果。在根系研究方面，向師慶等[2]根據(jù)對北京地區(qū)主要林木根系的挖掘研究，將林木根系歸納為垂直根型、斜生根型、水平根型等主要根型;Yang等[3]建立了根系錨固強度模型，模擬了根系連續(xù)錨固破壞過程，提出根組分對樹木錨地的等級影響如下：主根>迎風淺根>垂直淺根>迎風垂根>其他影響較小的組分，揭示了根系個體破壞與樹木總體傾倒破壞之間的根本聯(lián)系;王可鈞等[4]重點分析了林木根系的抗拉性能和生長方向與土質坡體穩(wěn)定性的關系，針對根系固坡的力學特性進行了簡析。不少學者對現(xiàn)場試驗獲得的相關數(shù)據(jù)進行分析，Coutts[5]使用絞車試驗對云杉進行了現(xiàn)場拉拔試驗，測量包括抗拉拔力、傾斜角、土壤和根系位移等方面，表明根的數(shù)量和尺寸分布的重要性;劉秀萍[6]采用野外試驗與室內(nèi)試驗相結合的方式，構建了刺槐、油松樹種根系形態(tài)分布模型，通過對根系與黃土接觸面單元摩阻特性的試驗研究，揭示了影響根土接觸面摩阻特性的因素。在單根參數(shù)方面，鄭力文等[7]通過對油松單根施加拔出荷載進行拉拔實驗，分析土壤含水率、干密度、加載速率等因素對拔出過程中根土界面摩擦性能的影響;夏振堯等[8]、管世烽等[9]通過干密度改變和多梯度含水率室內(nèi)試驗，測定了麥冬和多花木藍根系根土界面的抗拉拔力特性，為根土界面摩擦特性的進一步研究提供依據(jù);呂春娟等[10]選取了5 種常見喬木根系進行室內(nèi)單根拉伸試驗，發(fā)現(xiàn)了喬木根系極限抗拉力和直徑之間存在的關系，揭示了植被根系的固土力學機制;劉小光等[11]通過直接拉拔試驗研究了加載速率、根系埋深和不同樹種對根系與土體摩擦特性的影響，揭示了根系拉拔破壞過程中位移抗拉拔力關系曲線的階段性變化。此外，通過有限元軟件對根系形態(tài)的模擬，Rahardjo等[12]探究了土體抗剪強度參數(shù)和根系幾何形狀與林木傾倒最大風力的關系，研究了粗粒土改良對樹木穩(wěn)定性的影響;Fourcaud等[13]采用了有限元法對飽和軟黏土和土壤中的根系拔出進行二維模擬，模擬了樹根破壞時板滑移面形狀和旋轉軸位置，估算了覆土過程中土體和根系的應力分布，確定了不同生根方式下不同根系在樹根錨固中發(fā)揮的作用。學者們針對根系形態(tài)理論和小部分根段的抗拉、摩擦特性兩個方面開展了大量研究，然而，林木傾倒與根系整體的關聯(lián)性大，局部根段研究無法較好地反映其根系對林木抗拉拔性能的作用，目前，針對全根段抗拉拔性能的研究較少，需進一步探索研究。

為更好地分析林木受橫向拉拔作用下的傾倒破壞過程，研究不同因素的改變對根系抗拉拔性能的影響作用，從而為強風暴雨作用下林木傾倒防治提供一定的理論依據(jù)。本文著眼于強風、暴雨作用后土壤在降雨作用下含水率升高、林木受橫向風力作用的情況，通過室內(nèi)模擬絞車試驗，分析了在外力荷載下林木根系傾倒的破壞過程，通過改變含水率以及根系埋深，探究不同含水率及根系埋深對林木傾倒破壞的影響。

1 模型試驗

1. 1 試驗裝置

開展橫向拉拔荷載室內(nèi)試驗的裝置是自主設計的根系橫向加載試驗系統(tǒng)，該裝置由拉拔力試驗機、滑輪傳動軸、模型箱和含水率監(jiān)測裝置4部分組成。拉拔力試驗機的主要技術規(guī)格為：1）橫梁移動速度v=1～10 mm/s;2）荷重精度≤0.5%，荷載上限100 kN;3）最大行程1 000 mm。滑輪傳動軸為介于拉拔力試驗機與模型箱之間的T字型鋼架，鋼架上帶有可固定滑輪，拉拔力試驗機通過向上移動吊臂帶動鋼索，經(jīng)滑輪傳動軸改變力的方向后對樹木施加橫向拉拔力。吊臂上搭載壓力傳感器，并通過TM2101pro測控系統(tǒng)軟件記錄拉拔距離。模型箱尺寸長×寬×高為600 mm×600 mm×400 mm，側面軸線上，距離底面100、200、300 mm處分別打平行雙排孔，用于兩針型含水率傳感器探頭插入測量，含水率數(shù)據(jù)記錄處理采用配套的CYY-3055型8通道數(shù)據(jù)采集儀。模型箱鋼制底板試驗機鉚釘固定形成整體剛度。裝置整體形態(tài)如圖1所示。

1.2 試樣材料與方案

試樣材料：模型試驗樹木根系選用銀杏樹苗，為深根性樹種，主根明顯并起主導作用，側根較少且處從屬地位。選取幾何形態(tài)相近的樹苗，去除其細小分支，保留垂直主根備用。選取土壤為重慶地區(qū)粉質砂土，原狀土通過現(xiàn)場取土實驗計算得出原土樣含水率為13.3 %，干密度為1.65 g/cm3，在模型試驗中重塑土樣制備也以此為基準。篩前工樣如圖2所示。受模型箱大小限制，選取的銀杏根系根莖在28～32 mm之間，根條形態(tài)如圖3所示。

拉拔試驗方案：根據(jù)模型試驗需求，稱取適量土壤，依據(jù)對應試驗組的含水率需求進行土樣配制。用環(huán)刀取部分配制好的重塑土樣，開展直接剪切試驗。隨后，將測定后符合試驗要求的土壤采用稱重分層裝填、分層擊實的形式，使用夯錘將試驗土體每10 cm裝填擊實1次，共計4層。裝填至一定高度土層時，于模型箱中央埋入樹苗，后繼續(xù)裝填。裝填完成后，將含水率傳感器探頭插入模型箱側面鉆孔中，進行含水率校對及變化監(jiān)測，采用伸縮鋼圈將鋼索與樹根主干連接固定。裝填完成后，試樣靜置2 h以確保傳感器讀數(shù)精確，待各監(jiān)測點含水率讀數(shù)與試驗所需含水率數(shù)值相差<0.5%時，認定該土樣含水率符合試驗條件，隨后啟動試驗機施加外力荷載，以恒定速率8 mm/s進行拉拔試驗，直至采集系統(tǒng)記錄到外力荷載達到峰值后停止拉拔試驗，輸出根系拉拔破壞全過程位移抗拉拔力（S-F）曲線。通過控制變量法改變含水率或根系埋深，進行多組拉拔破壞試驗。模型試驗工況中根系埋深選取為15、20、25、30 cm，含水率控制區(qū)間為7%～22%。為減小試驗離散型偏差對試驗結果影響，每組工況重復進行多次試驗。試驗結束后土體破壞面具體形態(tài)如圖4所示。

2 結果與分析

2.1 含水率對土體抗剪強度的影響

土樣的直剪實驗根據(jù)《土工試驗規(guī)程》[14]開展，分別取含水率ω為7.0%、10.0%、13.0%、16.0%、19.0%和22.0%，共計6組土樣，每組土樣4個，根據(jù)環(huán)刀容積及所需的干密度，制備土樣所需的濕土量按式（1）計算。m0=（1+0.01ω0）ρdV

（1）式中：m0為濕土質量;ω0為濕土含水率;ρd為試樣的干密度;V為試樣體積（環(huán)刀容積）。

對同一含水率下的4組土樣分別施加100、200、300、400 kPa法向應力進行快剪試驗，剪切速率控制在0.8 mm/min。根據(jù)量力環(huán)系數(shù)的轉換關系，得到不同垂直壓力作用下的抗剪強度，以垂直壓力為橫坐標，剪切強度為縱坐標，繪制剪切強度與垂直壓力的關系曲線，將圖上各點采取線性擬合，直線的傾角即為土的內(nèi)摩擦角φ，截距為土的黏聚力c，從而得到不同含水率下土樣的c、φ值，土樣的抗剪強度按式（2）計算。τf=c+σ tan φ

（2）式中：τf為土體的抗剪強度;c為土體的黏聚力;σ為法向應力;φ為土體的內(nèi)摩擦角。

根據(jù)試驗得到不同含水率下土體受各級垂直壓力作用下的法向應力及推算出的土體強度指標，見表1。其中，x為法向應力，y為剪應力。

由圖5～圖6可知，黏聚力隨含水率的變化曲線分為兩部分：當含水率ω<13%時，土體黏聚力隨著含水率的增大顯著增加;當含水率ω>13%時，隨著含水率的進一步增大，黏聚力降低。在ω=22%時，黏聚力急劇減小，接近于0，這表明該含水率已接近土體的飽和含水率。含水率對內(nèi)摩擦角變化的影響要小于黏聚力，內(nèi)摩擦角隨含水率增大變化不明顯，在28.8°～32.9°間變化。土體的抗剪強度與黏聚力的變化規(guī)律類似，也隨著含水率的增加先增大，隨后逐漸減小。林鴻州等[15]、黃琨等[16]研究的非飽和土抗剪強度規(guī)律也得到了類似結論。

2.2 根系埋深對根系抗拉拔力的影響

分析S-F曲線可以得知（圖7），在控制拔出速率恒定的情況下，改變根系埋深進行銀杏根系模擬絞車試驗得到的S-F曲線存在著明顯的階段特性。由圖7可以看出，以30 cm埋深為例，在達到最大抗拉拔力對應的峰值位移前，抗拉拔力曲線呈現(xiàn)平緩上升（AB段左側）、陡峭上升（AB段與CD段之間）和緩慢上升達到峰值（CD段右側）3個階段。分析這3個階段可得出：在平緩上升階段，可以觀測到該階段土體表面并未發(fā)生明顯破壞，但樹干部分發(fā)生彎曲，抗拉拔力主要由樹干受橫向位移而產(chǎn)生彈性形變的抵抗力提供;隨著橫向荷載進一步增大，根系達到最大靜摩擦力后，開始與土體之間發(fā)生錯動，土體下方根系朝模擬迎風側剪出，土體上部枝干擠壓模擬背風側土壤，因此，抗拉拔力激增，形成了陡峭上升階段;繼續(xù)施加橫向位移，大多數(shù)根系周圍土體已發(fā)生錯動，達到抗拉拔力峰值，隨后土體剪破導致林木徹底破壞失穩(wěn)，發(fā)生傾倒破壞。

[6] 劉秀萍. 林木根系固土有限元數(shù)值模擬[D]. 北京： 北京林業(yè)大學， 2008.

LIU X P. Finite element method numerical simulation of forest roots reinforcement[D]. Beijing： Beijing Forestry University， 2008.（in Chinese）

[7] 鄭力文， 劉小光， 涂志華， 等. 土壤含水率與干密度對油松根-土界面摩擦性能的影響[J]. 中國水土保持科學， 2014， 12（6）： 36-41.

ZHENG L W， LIU X G， TU Z H， et al. Effects of moisture content and dry density of soil on the friction characteristics between Pinus tabulaeformis root and soil [J]. Science of Soil and Water Conservation， 2014， 12（6）： 36-41.（in Chinese）

[8] 夏振堯， 管世烽， 牛鵬輝， 等. 麥冬和多花木藍根系抗拉拔特性試驗研究[J]. 水土保持通報， 2015， 35（6）： 110-113.

XIA Z Y， GUAN S F， NIU P H， et al. Experimental study on pull-out resistive properties of ophiopogon japonicus and indigofera amblyantha roots [J]. Bulletin of Soil and Water Conservation， 2015， 35（6）： 110-113.（in Chinese）

[9] 管世烽， 夏振堯， 張倫， 等. 水平荷載作用下多花木藍根系拉拔試驗研究[J]. 長江科學院院報， 2016， 33（6）： 24-28.

GUAN S F， XIA Z Y， ZHANG L， et al. Pull-out test of indigofera amblyantha craib root under horizontal load [J]. Journal of Yangtze River Scientific Research Institute， 2016， 33（6）： 24-28.（in Chinese）

[10] 呂春娟， 陳麗華， 周碩， 等. 不同喬木根系的抗拉力學特性[J]. 農(nóng)業(yè)工程學報， 2011， 27（Sup1）： 329-335.

LYU C J， CHEN L H， ZHOU S， et al. Root mechanical characteristics of different tree species [J]. Transactions of the Chinese Society of Agricultural Engineering， 2011， 27（Sup1）： 329-335.（in Chinese）

[11] 劉小光， 趙紅華， 冀曉東， 等. 油松和落葉松根與土界面摩擦特性[J]. 摩擦學學報， 2012， 32（6）： 550-556.

LIU X G， ZHAO H H， JI X D， et al. Friction characteristics of root-soil interface of pinus tabulaeformis and larix gmelinii [J]. Tribology， 2012， 32（6）： 550-556.（in Chinese）

[12] RAHARDJO H， HARNAS F R， LEONG E C， et al. Tree stability in an improved soil to withstand wind loading [J]. Urban Forestry & Urban Greening， 2009， 8（4）： 237-247.

[13] FOURCAUD T， JI J N， ZHANG Z Q， et al. Understanding the impact of root morphology on overturning mechanisms： a modelling approach[J]. Annals of Botany， 2007， 101（8）： 1267-1280.

[14] 土工試驗規(guī)程：SL 237—1999[S]. 北京：中國水利水電出版社，1999.

Code of geotechnical testin： SL 237-1999 [S]. Beijing： China Water & Power Press， 1999. （in Chinese）

[15] 林鴻州， 李廣信， 于玉貞， 等. 基質吸力對非飽和土抗剪強度的影響[J]. 巖土力學， 2007， 28（9）： 1931-1936.

LIN H Z， LI G X， YU Y Z， et al. Influence of matric suction on shear strength behavior of unsaturated soils[J]. Rock and Soil Mechanics， 2007， 28（9）： 1931-1936.（in Chinese）

[16] 黃琨， 萬軍偉， 陳剛， 等. 非飽和土的抗剪強度與含水率關系的試驗研究 [J]. 巖土力學， 2012， 33（9）： 2600-2604.

HUANG K， WAN J W， CHEN G， et al. Testing study of relationship between water content and shear strength of unsaturated soils [J]. Rock and Soil Mechanics， 2012， 33（9）： 2600-2604. （in Chinese）

（編輯 胡玥）