傅胤榕 仉文岗

摘 要：中国是受热带台风灾害影响最为严重的国家之一，由于强风暴雨作用引起的林木根系拔出及倾倒破坏致使交通中断、车辆受损、人员伤亡等次生灾害与日俱增。为进一步研究土体物理特性对林木抗拉拔性能的影响，采用自主设计的根系横向加载试验系统，通过对银杏树苗施加定速率横向加载的拉拔试验，模拟外力荷载下根系的倾倒过程。试验设置了5组梯度含水率和4种埋深工况分析土体含水率、根系埋深对拔出过程中抗拉拔力的影响。研究表明：在粉质砂土中，黏聚力根系位移抗拉拔力曲线中抗拉拔力的变化可以概括为平缓上升、陡峭上升和缓慢上升达到峰值3个阶段;最大抗拉拔力随含水率提升呈现先升后降的总体趋势，与土体黏聚力变化的总体趋势相同，当含水率趋近饱和含水率时，最大抗拉拔力衰减显著;当含水率一定时，根系的最大抗拉拔力随着根系埋深的增加，在试验设置的埋深区间内呈线性趋势增长。

关键词：含水率;根系;抗拉拔力;埋深;模型试验

中图分类号：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

  隨着城市化进程的不断推进和绿色生态文明的发展，城市绿化覆盖率不断提升，绿化树在城市建设中起到净化城区空气、降低环境噪音、美化城市道路等多种作用。近年来，东南沿海省份陆续遭受“莫兰蒂”、“莲花”等台风入境，遭受强风和暴雨，部分城市行道树倒伏率接近30%，个别林木树种倒伏率达到59%，总受损率达96%以上[1]。由于林木倾倒导致城市交通道路阻断、车辆及公共基础设施受损，甚至造成严重的人员伤亡，因此，强风、暴雨引起的林木倾倒破坏现象不容轻视。

学者们对林木倒伏和根土相互作用方面进行了大量研究，取得了一系列成果。在根系研究方面，向师庆等[2]根据对北京地区主要林木根系的挖掘研究，将林木根系归纳为垂直根型、斜生根型、水平根型等主要根型;Yang等[3]建立了根系锚固强度模型，模拟了根系连续锚固破坏过程，提出根组分对树木锚地的等级影响如下：主根>迎风浅根>垂直浅根>迎风垂根>其他影响较小的组分，揭示了根系个体破坏与树木总体倾倒破坏之间的根本联系;王可钧等[4]重点分析了林木根系的抗拉性能和生长方向与土质坡体稳定性的关系，针对根系固坡的力学特性进行了简析。不少学者对现场试验获得的相关数据进行分析，Coutts[5]使用绞车试验对云杉进行了现场拉拔试验，测量包括抗拉拔力、倾斜角、土壤和根系位移等方面，表明根的数量和尺寸分布的重要性;刘秀萍[6]采用野外试验与室内试验相结合的方式，构建了刺槐、油松树种根系形态分布模型，通过对根系与黄土接触面单元摩阻特性的试验研究，揭示了影响根土接触面摩阻特性的因素。在单根参数方面，郑力文等[7]通过对油松单根施加拔出荷载进行拉拔实验，分析土壤含水率、干密度、加载速率等因素对拔出过程中根土界面摩擦性能的影响;夏振尧等[8]、管世烽等[9]通过干密度改变和多梯度含水率室内试验，测定了麦冬和多花木蓝根系根土界面的抗拉拔力特性，为根土界面摩擦特性的进一步研究提供依据;吕春娟等[10]选取了5 种常见乔木根系进行室内单根拉伸试验，发现了乔木根系极限抗拉力和直径之间存在的关系，揭示了植被根系的固土力学机制;刘小光等[11]通过直接拉拔试验研究了加载速率、根系埋深和不同树种对根系与土体摩擦特性的影响，揭示了根系拉拔破坏过程中位移抗拉拔力关系曲线的阶段性变化。此外，通过有限元软件对根系形态的模拟，Rahardjo等[12]探究了土体抗剪强度参数和根系几何形状与林木倾倒最大风力的关系，研究了粗粒土改良对树木稳定性的影响;Fourcaud等[13]采用了有限元法对饱和软黏土和土壤中的根系拔出进行二维模拟，模拟了树根破坏时板滑移面形状和旋转轴位置，估算了覆土过程中土体和根系的应力分布，确定了不同生根方式下不同根系在树根锚固中发挥的作用。学者们针对根系形态理论和小部分根段的抗拉、摩擦特性两个方面开展了大量研究，然而，林木倾倒与根系整体的关联性大，局部根段研究无法较好地反映其根系对林木抗拉拔性能的作用，目前，针对全根段抗拉拔性能的研究较少，需进一步探索研究。

为更好地分析林木受横向拉拔作用下的倾倒破坏过程，研究不同因素的改变对根系抗拉拔性能的影响作用，从而为强风暴雨作用下林木倾倒防治提供一定的理论依据。本文着眼于强风、暴雨作用后土壤在降雨作用下含水率升高、林木受横向风力作用的情况，通过室内模拟绞车试验，分析了在外力荷载下林木根系倾倒的破坏过程，通过改变含水率以及根系埋深，探究不同含水率及根系埋深对林木倾倒破坏的影响。

1 模型试验

1. 1 试验装置

开展横向拉拔荷载室内试验的装置是自主设计的根系横向加载试验系统，该装置由拉拔力试验机、滑轮传动轴、模型箱和含水率监测装置4部分组成。拉拔力试验机的主要技术规格为：1）横梁移动速度v=1～10 mm/s;2）荷重精度≤0.5%，荷载上限100 kN;3）最大行程1 000 mm。滑轮传动轴为介于拉拔力试验机与模型箱之间的T字型钢架，钢架上带有可固定滑轮，拉拔力试验机通过向上移动吊臂带动钢索，经滑轮传动轴改变力的方向后对树木施加横向拉拔力。吊臂上搭载压力传感器，并通过TM2101pro测控系统软件记录拉拔距离。模型箱尺寸长×宽×高为600 mm×600 mm×400 mm，侧面轴线上，距离底面100、200、300 mm处分别打平行双排孔，用于两针型含水率传感器探头插入测量，含水率数据记录处理采用配套的CYY-3055型8通道数据采集仪。模型箱钢制底板试验机铆钉固定形成整体刚度。装置整体形态如图1所示。

1.2 试样材料与方案

试样材料：模型试验树木根系选用银杏树苗，为深根性树种，主根明显并起主导作用，侧根较少且处从属地位。选取几何形态相近的树苗，去除其细小分支，保留垂直主根备用。选取土壤为重庆地区粉质砂土，原状土通过现场取土实验计算得出原土样含水率为13.3 %，干密度为1.65 g/cm3，在模型试验中重塑土样制备也以此为基准。筛前工样如图2所示。受模型箱大小限制，选取的银杏根系根茎在28～32 mm之间，根条形态如图3所示。

拉拔试验方案：根据模型试验需求，称取适量土壤，依据对应试验组的含水率需求进行土样配制。用环刀取部分配制好的重塑土样，开展直接剪切试验。随后，将测定后符合试验要求的土壤采用称重分层装填、分层击实的形式，使用夯锤将试验土体每10 cm装填击实1次，共计4层。装填至一定高度土层时，于模型箱中央埋入树苗，后继续装填。装填完成后，将含水率传感器探头插入模型箱侧面钻孔中，进行含水率校对及变化监测，采用伸缩钢圈将钢索与树根主干连接固定。装填完成后，试样静置2 h以确保传感器读数精确，待各监测点含水率读数与试验所需含水率数值相差<0.5%时，认定该土样含水率符合试验条件，随后启动试验机施加外力荷载，以恒定速率8 mm/s进行拉拔试验，直至采集系统记录到外力荷载达到峰值后停止拉拔试验，输出根系拉拔破坏全过程位移抗拉拔力（S-F）曲线。通过控制变量法改变含水率或根系埋深，进行多组拉拔破坏试验。模型试验工况中根系埋深选取为15、20、25、30 cm，含水率控制区间为7%～22%。為减小试验离散型偏差对试验结果影响，每组工况重复进行多次试验。试验结束后土体破坏面具体形态如图4所示。

2 结果与分析

2.1 含水率对土体抗剪强度的影响

土样的直剪实验根据《土工试验规程》[14]开展，分别取含水率ω为7.0%、10.0%、13.0%、16.0%、19.0%和22.0%，共计6组土样，每组土样4个，根据环刀容积及所需的干密度，制备土样所需的湿土量按式（1）计算。m0=（1+0.01ω0）ρdV

（1）式中：m0为湿土质量;ω0为湿土含水率;ρd为试样的干密度;V为试样体积（环刀容积）。

对同一含水率下的4组土样分别施加100、200、300、400 kPa法向应力进行快剪试验，剪切速率控制在0.8 mm/min。根据量力环系数的转换关系，得到不同垂直压力作用下的抗剪强度，以垂直压力为横坐标，剪切强度为纵坐标，绘制剪切强度与垂直压力的关系曲线，将图上各点采取线性拟合，直线的倾角即为土的内摩擦角φ，截距为土的黏聚力c，从而得到不同含水率下土样的c、φ值，土样的抗剪强度按式（2）计算。τf=c+σ tan φ

（2）式中：τf为土体的抗剪强度;c为土体的黏聚力;σ为法向应力;φ为土体的内摩擦角。

根据试验得到不同含水率下土体受各级垂直压力作用下的法向应力及推算出的土体强度指标，见表1。其中，x为法向应力，y为剪应力。

由图5～图6可知，黏聚力随含水率的变化曲线分为两部分：当含水率ω<13%时，土体黏聚力随着含水率的增大显著增加;当含水率ω>13%时，随着含水率的进一步增大，黏聚力降低。在ω=22%时，黏聚力急剧减小，接近于0，这表明该含水率已接近土体的饱和含水率。含水率对内摩擦角变化的影响要小于黏聚力，内摩擦角随含水率增大变化不明显，在28.8°～32.9°间变化。土体的抗剪强度与黏聚力的变化规律类似，也随着含水率的增加先增大，随后逐渐减小。林鸿州等[15]、黄琨等[16]研究的非饱和土抗剪强度规律也得到了类似结论。

2.2 根系埋深对根系抗拉拔力的影响

分析S-F曲线可以得知（图7），在控制拔出速率恒定的情况下，改变根系埋深进行银杏根系模拟绞车试验得到的S-F曲线存在着明显的阶段特性。由图7可以看出，以30 cm埋深为例，在达到最大抗拉拔力对应的峰值位移前，抗拉拔力曲线呈现平缓上升（AB段左侧）、陡峭上升（AB段与CD段之间）和缓慢上升达到峰值（CD段右侧）3个阶段。分析这3个阶段可得出：在平缓上升阶段，可以观测到该阶段土体表面并未发生明显破坏，但树干部分发生弯曲，抗拉拔力主要由树干受横向位移而产生弹性形变的抵抗力提供;随着横向荷载进一步增大，根系达到最大静摩擦力后，开始与土体之间发生错动，土体下方根系朝模拟迎风侧剪出，土体上部枝干挤压模拟背风侧土壤，因此，抗拉拔力激增，形成了陡峭上升阶段;继续施加横向位移，大多数根系周围土体已发生错动，达到抗拉拔力峰值，随后土体剪破导致林木彻底破坏失稳，发生倾倒破坏。

[6] 刘秀萍. 林木根系固土有限元数值模拟[D]. 北京： 北京林业大学， 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]. 农业工程学报， 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] 土工試验规程：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）

（编辑 胡玥）