全 文 :第42卷 第8期
2014年 8月
华 中 科 技 大 学 学 报 (自 然 科 学 版)
J.Huazhong Univ.of Sci.&Tech.(Natural Science Edition)
Vol.42No.8
Aug. 2014
DOI:10.13245/j.hust.140822
收稿日期 2013-11-16.
作者简介 万 娟(1977-),女,讲师;肖衡林(通信作者),教授,E-mail:7531450@qq.com.
基金项目 国家自然科学基金资助项目 (51178166);2012年湖北省科技厅创新群体项目 (2012FFA035);教育部
新世纪优秀人才支持计划资助项目 (NCET-11-0962).
多花木兰根系抗拉特性及边坡稳定性分析
万 娟1,2 肖衡林2 何 俊2
(1华中科技大学土木工程与力学学院,湖北 武汉430074;
2湖北工业大学土木工程与建筑学院,湖北 武汉430068)
摘要 以护坡典型植物多花木兰根系为对象进行室内抗拉试验;在假定滑动面平行于坡面的前提下,推导生
态护坡中坡体稳定系数的计算公式;研究不同坡度下根系的抗拉力、抗拉强度、弹性模量及边坡的稳定性.研
究表明:多花木兰根系抗拉力随着直径的增大呈幂函数增加;根系抗拉强度为21.15~76.08MPa,随着直径
的增加呈幂函数减小;坡度越大,同一直径根系的抗拉力与抗拉强度越大.多花木兰根系弹性模量为0.64~
2.33GPa,随着直径的增加其弹性模量逐渐减小;根系直径为1.2~2.2mm时,弹性模量减小的幅度大,其他
区间减小的幅度较小.坡体的稳定系数与根系的抗拉力呈正相关关系,与边坡的坡度呈负相关关系,并受根的
生长方向和密度、植物的种植密度以及原状土的性质等影响.坡体的稳定系数随根系直径的增加呈线性增加.
关键词 护坡;根系;抗拉力;抗拉强度;弹性模量;稳定系数
中图分类号 TU411.3 文献标志码 A 文章编号 1671-4512(2014)08-0109-05
Tensile properties of Magnolia Multiflora′s roots
and slope stability analysis
Wan Juan1,2 Xiao Henglin2 He Jun2
(1School o f Civil Engineering and Mechanics,Huazhong University of Science and Technology,Wuhan 430074,
China;2School o f Civil Engineering and Architecture,Hubei University of Technology,Wuhan 430068,China)
Abstract The tensile tests indoors were conducted targeted roots of Magnolia Multiflora,a typical
slope plant.On the assumption that the sliding surface was paralel to the slope,the formula of slope
stability coefficient was deduced.The tension,tensile strength and modulus of roots and the stability
of the slope were researched on different gradients.The results show that the tensile force of Magnoli-
a Multiflora′s roots increase with diameter increasing according to a power function.Their tensile
strength is between 21.15to 76.08MPa,which decreases with diameter increasing according to a
power function.The larger the gradient,the greater are the roots′tensile force and tensile strength
with the same diameter.Roots'modulus is between 0.64to 2.33GPa,which reduce gradualy along
with the increase of roots′diameter.The rate of the decrease is sharp when the diameter is between
1.2and 2.2mm,that is slower in other cases.The slope stability coefficient is positively associated
with roots′tensile force,negative correlation with the slope gradient,and influenced by many factors,
such as roots′growth direction and density,plants′planting density and soil′s properties.The slope
stability coefficient increases linearly with the increase of root diameter.
Key words slope protection;roots;tensile force;tensile strength;modulus;stability coefficient
公路工程等基础设施的建设产生大量裸露的 边坡,破坏了原有的生态环境,导致水土流失,甚
至滑坡失稳.植被作为一种护坡材料,用于减轻坡
面的不稳定性和侵蚀[1].植物根系与土之间的相
互作用,主要有准黏聚力原理、摩擦加筋原理与等
效围压原理[2-3].
根系固土护坡的试验研究主要集中在根系的
抗拉强度及根土复合体的抗剪强度两类.抗拉强
度的测量可通过原位拉拔试验[4-7]与实验室抗
拉[8-11]来实现.
根系的抗拉强度受到植物的种类、根系的立
地环境等诸多因素的影响,故需大量的数据来充
实.对于不同生长期、不同坡度下根系的抗拉特性
及其对边坡稳定性的定量分析较少.
本研究选用典型的护坡灌木多花木兰的根系
进行室内抗拉试验,分析其根系在不同的坡度下
的抗拉力、抗拉强度及本构关系曲线.目的是探明
在外部环境变化的条件下根系的抗拉特性及其对
边坡稳定性的影响,揭示根系的护坡机理.
1 试验材料与方法
根系的抗拉特性是坡体稳定的重要参数,本
研究选取灌木多花木兰根系为研究对象.
试验地为湖北工业大学生态护坡基地.基地
为不同坡度(0°,37°,45°与53°)草灌混播试验地.
为了防止土体的下滑,边坡都采用三维土工网固
定[12-13].每隔一段时间从四个不同的坡度选取上
部高度大致相等的多花木兰为试样,挖出的根系
用水冲掉其上的泥土,选取根系中较直的部分,用
剪刀剪取约10cm左右贴好标签,立即放入保鲜
袋,在24h内试验.
试验设备采用济南华衡试验设备有限公司的
微机控制电子万能试验机 WDW-2,试验量程为2
kN,多花木兰根系试验时的拉伸速度取1mm/
min.
为避免根系在夹具处滑脱,实验时根系在夹
具处增加橡皮垫以增大摩擦.根系拉断后用游标
卡尺测量根断裂处直径D.试验可以得出根系的
应力应变曲线及抗拉力T.选取应力 -应变曲线中
直线段的斜率作为根系的弹性模量E.根系抗拉
强度为
p=4T/(πD2). (1)
2 试验结果与分析
为了分析多花木兰根系的抗拉特性,对不同
坡度的根系进行抗拉试验,试验结果见表1.
表1 多花木兰根系抗拉试验结果
编号 θ/(°) D/mm T/N p/MPa E/GPa
1
2
3
4
5
0
0.93
1.52
1.65
2.16
3.48
24
58.9
65.6
78.8
201.2
35.33
32.46
30.68
21.50
21.15
1.77
1.47
1.28
0.87
0.64
6
7
8
9
10
37
0.72
1.39
1.63
2.02
2.78
17.4
58.3
65.3
75.8
140.7
42.74
38.42
31.29
23.65
23.18
1.89
1.67
1.30
0.94
0.77
11
12
13
14
15
45
0.57
1.31
1.42
1.87
2.48
13.4
58.1
63.2
74.6
127.4
52.51
43.11
39.91
27.16
26.37
1.94
1.72
1.52
1.12
0.88
16
17
18
19
20
53
0.45
1.21
1.28
1.41
2.32
12.1
54.9
59.7
64.5
118.4
76.08
47.74
46.39
41.31
28.01
2.33
1.82
1.64
1.35
0.96
2.1 根系抗拉力、抗拉强度分析
为了分析影响抗拉力、抗拉强度的因素,作出
不同坡度下抗拉力、抗拉强度与直径的关系曲线,
如图1和图2所示(R2 为相关系数).
1—θ=0°,T=28.013 D1.552 1,R2=0.981 3;2—θ=
37°,T=30.334 D1.496 4,R2=0.978 5;3—θ=45°,T=
33.587 D1.502 1,R2 = 0.979 3;4—θ=53°,T =
39.37 D1.412 1,R2=0.993 4.
图1 抗拉力随直径与坡度的变化
由图1和图2可知:多花木兰根系直径为
0.45~3.48mm,抗拉力为12.1~201.2N,拟合
曲线相关系数均在0.95以上.在同一坡度下,抗
拉力随着直径的增加而增大.抗拉力与直径存在
着幂函数关系,即T=a1Db1,其中:b1 大小在1.5
·011· 华 中 科 技 大 学 学 报 (自 然 科 学 版) 第42卷
1—θ=0°,p=35.668 D-0.447 9,R2=0.813 7;2—θ=
37°,p=38.623 D-0.503 6,R2=0.837 8;3—θ=45°,
p=42.764 D-0.497 9,R2=0.838 9;4—θ=53°,p=
50.128 D-0.587 9,R2=0.962 8.
图2 抗拉强度随直径与坡度的变化
左右变化;a1 随着坡度的增加而增大.多花木兰
根系抗拉强度为21.15~76.08MPa,与 HPB235
钢筋的抗拉强度235MPa相比,多花木兰根系的
抗拉强度约为钢筋的1/10~1/4,拟合曲线相关
系数均在0.8以上.在同一坡度下,抗拉强度随着
直径的增加而减小.抗拉强度与直径存在着幂函
数关系,即p=a2Db2,其中:b2 大小在-0.5左右
变化;a2 随着坡度的增加而增大.同一直径的根
系坡度越大,抗拉力与抗拉强度越大.
2.2 根系应力 -应变关系与弹性模量分析
图3给出在平地上不同的多花木兰根系的应
力 -应变曲线,图中1~5为根系编号.分析可知:
多花木兰根系受拉初期表现出较好的弹性,随着
应力的增加曲线呈锯齿状上升,这时根系表现出
一定的塑性变形.
此次试验根系的破坏表现出三种形式:1号
根系在达到应力峰值后直接在根系的细端被拉
断,断口参差不齐;2号与3号根系应力瞬间下降
后仍有较小拉力,原因是根系已大部分损伤,但没
有完全拉断;4号与5号根系应力 -应变曲线呈锯
齿状下降,根系逐步损伤.后两种形式最后阶段应
变不断增大,而力基本上保持不变,说明应力经过
下降阶段的振动,根系与夹具之间出现了一定的
松动,根系相对于夹具出现了滑移.可见,随着根
系的生长其塑形越来越明显.其他坡度的根系的
曲线有类似的规律.
把不同坡度的多花木兰根系的弹性模量绘于
图4.可以看出:在同一坡度下,伴随着根系的成
图3 平地上多花木兰根系应力 -应变曲线
1—θ=0°;2—θ=37°;3—θ=45°;4—θ=53°.
图4 不同坡度的多花木兰根系弹性模量随直径
的变化
长,其弹性模量减小,即抵抗变形的能力减小.根
系的弹性模量在0.64~2.33GPa之间,其弹性
模量远大于土体的弹性模量.各个坡度的根系的
曲线有类似的规律.平坡上生长的根系,其直径D
为1.52~2.16mm弹性模量减小的幅度较大,而
D<1.52mm或D>2.16mm时下降幅度相对平
缓;37°坡度上生长的根系,D=1.39~2.02mm
弹性模量减小的幅度较大,而D<1.39mm 或
D>2.02mm时下降幅度相对平缓;45°坡度上生
长的根系,D为1.31~1.87mm弹性模量减小的
幅度较大,而D<1.31mm或D>1.87mm时下
降幅度相对平缓;53°坡度上生长的根系,D 为
1.21~1.41mm 弹性模量减小的幅度较大,而
D<1.21或D>1.41mm时下降幅度相对平缓.
可见:根系的弹性模量减小的幅度与直径的大小
相关,在D为1.2~2.2mm区间内,随着直径的
增大弹性模量减小的速度快;直径偏小或偏大时,
随着直径的增加弹性模量减小的速度相对较慢,
·111·第8期 万娟,等:多花木兰根系抗拉特性及边坡稳定性分析
原因可能是直径在1.2~2.2mm区间内正是根
系的木质素急剧增加的时期.
3 根系的护坡稳定性分析
3.1 理论计算
为了明确根系的抗拉强度、分布特征对边坡
稳定性的影响,取滑坡体进行分析.在斜坡中,往
往产生与坡面或坡底面平行的压致拉裂面[14],这
里取滑动面平行于坡面的滑坡体为研究对象,其
受力图如图5所示,图中:G 为滑坡体的重力,θ
为边坡的倾角,N 为下部土体对滑坡体的支撑
力,Cf和为土体的黏聚力和内摩擦角,A 为滑
动面的面积.
图5 滑坡体的受力图
无植物根系时,滑坡体沿坡面的下滑力为
Gsinθ,抗滑力为下部土体对上部滑动土体的黏
结力与摩擦力之和,即F=CfA+Gcosθtanφ.无
根系时滑坡体的稳定系数为
K1 =F/(Gsinθ), (2)
若存在与滑动面夹角为α的根系,则滑坡体
发生滑移,0°<α≤90°时,根系受拉;90°<α<180°
时,根系受压.
取一根系进行分析,当0°<α≤90°时,根系受
拉,土体滑移不仅要克服自身的黏聚力与内摩擦
力,还要克服根系的抗拉力T.抗拉力产生了两个
效果,一是产生了垂直滑动面的向下分力Tsinα,
使土体间的摩擦力增大了Tsinαtan;二是产生
平行于滑动面的向上分力Tcosα,也提高了滑坡
体的稳定性.此时滑坡体的稳定系数为
K′2 =F/(Gsinθ)+[T(cosα+
sinαtan)]/(Gsinθ). (3)
当α>90°时,根系在滑坡体滑动时受压,其抗
拉强度未得到发挥,对提高滑坡体的稳定性无作
用.若根系较好地锚固在稳定岩土层上,则坡体滑
动还须克服根系自身的抗剪力.若这部分的根系
较细小,则对稳定性的提高可忽略.
若单株植物质量为G′,该植物有n个根系与
滑动面的夹角为0°<αi≤90°(i=1,2,…,n)时,则
坡体的稳定系数为
K″2 [= CfA+(G+G′)cosθtan+
∑
n
i=1
(Ticosαi)+tan∑
n
i=1
(Tisinαi ])
[(G+G′)sinθ]. (4)
取沿着坡面向上为x 轴,平行于坡面为y
轴.设坡面上同种护坡植物,沿x轴间距为Sx,沿
y轴间距为Sy;滑坡体厚度为 H,土体的容重为
γ.
假设每株植物根系分布相同,则植物的存在
产生两个方面的影响:一方面根系的存在提高了
坡体的稳定系数;另一方面植物的质量使坡体的
稳定性降低.于是坡体的稳定系数为
K2 {= CfSxSy+
∑
n
i=1
(Ticosαi)+tan∑
n
i=1
(Tisinαi[ ]})
[(SxSyγH +G′)sinθ]+tan/tanθ. (5)
可见:坡体的稳定系数与根系的抗拉力呈正相关
关系,与边坡的坡度呈负相关关系,还与根的生长
方向和密度、植物的种植密度以及原状土的性质
等密切相关.
3.2 实例分析
为了分析护坡基地中根系对坡体的稳定性影
响,测得土的容重γ=18kN/m3,Cf=2.5kPa,
=8°,滑坡体厚度取 H=0.1cm,多花木兰间距
取Sx=50cm,Sy=50cm.由于多花木兰质量远
小于土体质量,因此计算时忽略不计.为了简化计
算,多花木兰根系与竖直方向夹角均取75°,有效
根系均取5根.相关数据代入式(5)计算坡体的稳
定系数.坡体的稳定系数K 与直径D 拟合曲线的
关系式与相关系数R2 列于表2中.
表2 不同坡体稳定系数K与根系直径D 的拟合曲线
θ/(°) K R2
37 0.818 9 D+2.211 1 0.857 4
45 0.487 9 D+2.026 3 0.968 9
53 0.327 8 D+1.823 9 0.992 2
可以看出,在同一坡度下,随着根系的生长稳
定系数逐渐增加,多花木兰的后期护坡效果好于
前期护坡效果.稳定系数与根系直径具有高度的
线性相关性,相关系数均大于0.85,且坡度越小
稳定系数增加的幅度越大.在相同的直径下,随着
坡度的增加稳定系数降低,即坡越陡越不稳定.
4 结论
a.多花木兰根系抗拉力随着直径的增加而
·211· 华 中 科 技 大 学 学 报 (自 然 科 学 版) 第42卷
增大,抗拉力与直径存在着幂函数关系.根系的抗
拉强度为21.15~76.08MPa,与HPB235的钢筋
的抗拉强度235MPa相比,多花木兰根系的抗拉
强度约为钢筋的1/10~1/4.抗拉强度随着直径
的增加而减小,抗拉强度与直径存在着幂函数关
系.坡度越大,同一直径根系的抗拉力与抗拉强度
越大.
b.多花木兰根系的受拉破坏过程是一个从弹
性过渡到弹塑性的过程,随着直径的增加其塑形
越来越明显.根系的弹性模量为0.64~2.33
GPa,其弹性模量远大于土体的弹性模量.伴随着
根系的生长其弹性模量逐渐减小.根系的弹性模
量减小的幅度与直径的大小相关,当D=1.2~
2.2mm时,随着直径的增大弹性模量减小的速
度快;当直径偏小或偏大时,随着直径的增加弹性
模量减小的速度相对较慢.
c.在假定滑动面平行于坡面的前提下,推导
出生态护坡的坡体稳定系数的计算公式.植物的
存在产生两个方面的影响:一方面根系的存在提
高了坡体的稳定系数;另一方面植物的质量使坡
体的稳定性降低.坡体的稳定系数与根系的抗拉
力呈正相关关系,与边坡的坡度呈负相关关系,还
与根的生长方向和密度、植物的种植密度以及原
状土的性质等密切相关.
d.在同一坡度下,随着根系的生长坡体稳定
系数逐渐增加,说明多花木兰的后期护坡效果好
于前期护坡效果.坡体稳定系数与根系直径有着
高度的线性相关性,且坡度越小稳定系数增加的
幅度越大.在相同的根系直径下,随着坡度的增加
稳定系数降低,即坡越陡越不稳定.
参 考 文 献
[1]周德培,张俊云.植被护坡工程技术[M].北京:人
民交通出版社,2003.
[2]Waldron L J.The shear resistance of root-permeated
homogeneous and stratified soil[J].Soil Science Soc
Amer,1977,41:843-849.
[3]封金财,王建华.植物根的存在对边坡稳定性的作用
[J].华东交通大学学报,2003:20(5):42-45.
[4]Docker B B,Hubble T C T.Quantifying root-rein-
forcement of river bank soils by four Australian tree
species[J].Geomorphology,2008,100(3-4):401-
418.
[5]李臻,余芹芹,杨占风,等.西宁盆地两种灌木植物原
位拉拔试验及其护坡效应[J].水土保持研究,2011,
18(3):206-209.
[6]Riestenberg M M,Sovonick D S.The role of woody
vegetation in stabilizing slopes in the Cincinnati area,
Ohio[J].Geological Society of American Buletin,
1983,94(4):506-518.
[7]李国荣,胡夏嵩,毛小青,等.寒旱环境黄土区灌木根
系护坡力学效应研究[J].水文地质工程地质,
2008(1):94-97.
[8]朱海丽,胡夏嵩,毛小青,等.青藏高原黄土区护坡灌
木植物根系力学特性研究[J].岩石力学与工程学
报,2008,27(增2):3445-3452.
[9]Yang Yonghong,Liu Shuzhen,Wang Chenghua,et
al.A study of tensile strength tests of Arborous spe-
cies root system in forest engineering technique of
shalow landslide[J].Wuhan University Journal of
Natural Sciences,2006,11(4):892-896.
[10]朱清科,陈丽华,张东升,等.贡嘎山森林生态系统
根系固土力学机制研究[J].北京林业大学学报,
2002,24(4):64-67.
[11]程洪,张新全.草本植物根系网固土原理的力学试
验探究[J].水土保持通报,2002,22(5):20-23.
[12]刘晓路、肖衡林.浅谈三维土工网垫边坡防护体系
的建立[J].路基工程,2006,(6):114-117.
[13]肖成志,孙建诚,李雨润,等.三维土工网垫植草护
坡防坡面径流冲刷的机制分析[J].岩土力学,
2011,32(2):453-458.
[14]陈丽华,余新晓,宋维峰,等.林木根系固土力学机
制[M].北京:科学出版社,2008.
·311·第8期 万娟,等:多花木兰根系抗拉特性及边坡稳定性分析