Influences of morphology and hierarchy of roots on mechanical characteristics of root-soil composites
-
摘要: 植物根系的结构特征包括形态和层次结构,对根土复合体力学特性影响显著。通过开展黑麦草生长参数测试及直剪试验,获得了不同时期根土复合体的力学特征;采用自主研发的MechRoot程序建立了符合黑麦草根系结构特征的根土复合体模型,研究黑麦草不同形态和层次结构根系在直剪过程中的轴力水平以及占比情况,阐明了根系形态与层次结构对根土复合体力学特性影响及固土作用机理。研究表明:植物根系能明显提高土体的强度,根土复合体的抗剪强度增加主要是由黏聚力增加引起,增幅最大达4.99 kPa;随根系形态复杂程度增加,剪切过程中根系能够调动更大范围土体抵抗剪切变形,根系周围的剪切带和塑性区分布范围不断增加,并向根系周围集中,根土复合体的抗剪强度提高;根土复合体剪切过程中,各层次根系发挥作用逐渐变化,一级根、二级根、三级根的轴力水平分别为3.87,1.50,0.15 N,随根系生长发展,二级根、三级根参与到根系固土作用的比例不断增大,最大贡献分别为43.69%,13.80%。Abstract: The structural characteristics of plant roots, including morphology and topological structure, have significant effects on the mechanical characteristics of root-soil composites. The mechanical behaviors of the root-soil composites at different time are obtained by carrying out the Lolium perenne root system morphology parameter tests and direct shear tests. By using the self-developed MechRoot program, a numerical model for the root-soil composites with more realistic root structural characteristics of Lolium perenne is established, and the axial force levels and proportions of roots with different shapes and topological structures in the process of direct shear are studied, which illustrates the influences of morphology and topological structures of roots on the mechanical characteristics of root-soil composites and the mechanism of soil consolidation. The results show that the plant roots can strengthen soil significantly, and the increase of shear strength of root-soil composites is mainly caused by the increase of cohesion, with the maximum increase of 4.99 kPa. During the shearing process, with the increase of the morphological complexity of roots, the more range of soil can participate in the resistance of shear by the roots, the shear zone and plastic zone around the roots increase and are concentrated around the roots, and the shear strength of the root-soil composites increases. At the same time, the effects of the roots at several levels change gradually. The axial forces of the primary roots, secondary roots and fibrous roots are 3.87, 1.50 and 0.15 N. With the growth of the roots, the proportion of the secondary roots and fibrous roots participating in soil fixation of the roots increases continuously, with the maximum contributions being 43.69% and 13.80%, respectively.
-
0. 引言
对于饱和砂土地震液化触发判别问题,应用最广泛的Seed简化判别方法主要建立在包含液化宏观震害现象、地震动加速度记录和场地原位测试指标三要素的地震液化实例数据和液化分界线的基础上[1-3]。以基于原位剪切波速测试指标的液化实例为例,图1给出了地震液化和非液化实例,以及区分液化与否的分界线[4]。由于强震液化事件彻底破坏易液化土体原有结构性,导致了消除液化土层既往沉积作用的“时钟重置”效应,使得震后获得的场地原位测试指标与真正控制场地抗液化性能的震前土体工程力学特性不对应。考虑到针对地震液化场地的原位测试一般发生在震后数天或数月甚至更长的时间里,图1中地震液化实例数据的剪切波速和相应的液化分界线可能存在与测试时间相关的系统误差。因此,研究震后易液化场地饱和砂土土层的剪切波速(或其它原位测试指标)随时间的变化规律,对正确认识现有基于原位测试指标的液化判别简化方法的可靠性并对震后测试指标进行合理修正具有重要的科学意义。
土体在自然沉积过程中或者经受外荷载剧烈扰动后,其刚度和强度等工程力学参数在恒定应力状态下随时间变化的现象被称为“时间效应”,且该变化可以通过原位测试指标(例如,SPT、CPT和Vs)体现[5-7]。例如,Andrus等[8]通过静力触探与剪切波速的经验关系获得了原位砂土剪切波速和完全扰动后剪切波速的比值MEVR,初步揭示了地质时间尺度上的砂土时间效应(见图2);周燕国等[9]结合剪切波速(剪切模量)的Hardin公式[10]和静力触探测试,提出了砂土时间效应的表征指标AI,并利用不同地质年代砂土实测数据对其表征能力进行了检验。然而,上述研究主要面向地质沉积年代和土性测试指标的经验关系,关注数万年到数十万年甚至更长的地质时间尺度,并没有考虑地质沉积年代中可能受到的各种扰动(如强震事件)对土体工程力学特性的影响。以图2为例,深色椭圆圈内的数据并不符合图中拟合斜线的趋势,而这个时间范围(10-1~102 a)正是实际工程服役寿命或者高烈度区强震重现周期的范围。Pavlenko等[11]研究1995年神户地震中场地刚度演化问题,发现液化使得场地刚度急剧降低,但在震后较长一段时间内无法恢复到地震液化前的值。可见,对于实际工程,震后较短的时间尺度内场地土体的时间效应更具有工程意义,是影响易液化场地在余震中的地震响应和再液化可能性的重要因素。
针对上述问题,本文对新西兰基督城的REHS强震台站在2010—2011年间的若干强震记录的尾波进行水平与竖向谱比(HVSR)分析,获得了该台站所在的易液化场地在地震前后小应变剪切刚度随时间的发展规律。在此基础上,提出了综合考虑主固结和次固结作用的震后饱和砂土小应变刚度计算模型,并合理预测了REHS台站场地震后刚度恢复变化过程。
1. 场地条件与强震台站记录分析
HVSR方法通过对强震台站的三向加速度记录进行分析得到场地卓越频率以及场地土层结构分布,相比于其它场地频率分析方法,HVSR方法相对简单,能够基本消除震源和传播路径的影响,得到的场地卓越频率更可靠[12-13]。场地的卓越频率与场地平均剪切波速的关系可表示为[14]
vsavg=4Hf。 (1) 式中,
vsavg 为场地平均剪切波速,H为场地土层厚度,f为场地卓越频率。根据弹性理论,土体小应变剪切模量表达为
Gmax=ρ⋅v2s。 (2) 本文针对新西兰基督城CBD区的REHS台站(见图3),通过触发液化的强震本身和之后发生余震的加速度记录的尾波信号,获得场地卓越频率,根据式(1),(2)分析从2010年以来3次强震事件之后场地刚度随时间变化规律。REHS台站所处场地的上覆软弱土层厚度约为20 m,存在较厚的饱和粉质砂土层和洁净砂土层,下卧刚度较大的砂卵石层(见图4),是较典型的易液化场地[15]。
本文采用Geopsy软件进行HVSR分析,该软件能够多窗口同时计算H/V谱比和上覆土层的卓越频率f,并计算它们的平均值和方差。H/V谱比曲线的峰值对应的频率即为该尾波窗口时间对应的卓越频率。图5是尾波窗口选取和计算结果示意图,其中黑色实线代表所有窗口的平均值,两条虚线分别代表它的上、下限,竖向灰色条状区域代表了平均峰值所在的频率段。采用地震动记录的尾波进行HVSR分析获得场地卓越频率需满足一定的条件,即首先要保证主震后的选取的余震强度足够小,一是为了避免该余震事件对场地刚度恢复过程产生额外的显著影响,二是为了尾波段对应场地土体处于小应变弹性应变响应状态;其次,为保证数据处理结果的可靠性,均要求H/V谱比值大于2.0,而且卓越频率的误差在±0.3 Hz以内。
2. 震后场地小应变刚度演化规律
图6给出了REHS台站场地2010年至2011年间3次强震后场地卓越频率随时间的发展情况。考虑到每次强震事件之后大部分余震在台站的地震动强度都很小,可以忽略其对震后场地的进一步扰动作用。每张图上方的黑色虚线代表2010年9月3日Mw7.2地震事件之前场地的卓越频率(fi=2.20 Hz),由2010年大震之前的小震记录进行HVSR分析得到;右上角红色虚线表示强震事件后场地的稳定值,可视为下次强震前的场地卓越频率初始值。由图6可见,每次强震后场地刚度恢复至相对稳定的状态所需要的时间基本一致(为数天到数十天),而且场地稳定后的卓越频率均值接近但低于最初的水平(即2.20 Hz)。
由图6(b)可见,20110221 Mw6.2地震事件中场地卓越频率受到强余震地震历史的影响,因此图7仅将图6(a)和(c)的归一化场地卓越频率(f/fi)随时间发展规律归纳在一起。由图7可见,震后场地刚度恢复速率基本一致,呈现对数形式的增长,这与PENG等[16-17]的研究发现类似。而基督城场地小应变刚度在震后10~100 d会有显著的恢复。根据Downes等[18]的研究,2010年距上次基督城发生强震液化的时间间隔接近100 a,可以认为基督城典型易液化场地的新近沉积时间(或液化的重现周期)为100 a左右,所以尽管REHS台站场地的砂土层刚度在液化后的数月内会有较大的恢复,但仍小于震前水平。需要指出的是,在基督城的其它强震台站(如CBGS和CCCC台站)场地也观测到了类似规律[19]。另外,注意到图7中同一时间对应的频率数据存在一定离散性,原因包括:首先,HVSR方法不能完全消除震源及传播路径的影响;其次,地震动加速度信号(尤其是小震事件)信噪比水平不同;再次,地下水位的季节性升降等环境变化带来的可能影响。
3. 饱和砂土震后小应变刚度计算模型
需要注意,Andrus等[8]将砂土的MEVR参数用于剪切模量的时间效应修正,这种修正主要反映数十万年甚至更久的由于颗粒胶结引起的时间效应,不能合理反映砂土液化后短期内刚度的变化规律(即孔压消散与次固结引起的时间效应),而大部分震后实例调查的原位测试均是在这一时间段获得的;其次,这种修正对特定地震和场地数据的依赖性较大,数据十分有限,规律的一般性不强。本节基于饱和砂土的主固结与次固结过程,提出了考虑有效应力变化和砂土颗粒间接触力均质化(包括孔隙比变化)这两个主要因素的震后饱和砂土小应变刚度计算模型。
3.1 主固结影响函数
首先考虑由孔压消散和孔隙变化引起的剪切模量增量,可采用Hardin公式表示土体小应变剪切模量与土体孔隙比及有效应力的关系:
Gmax=A⋅F(e)(σ′m)n, (3) 式中,Gmax为砂土小应变剪切模量,A为结构性参数,F(e)为孔隙比方程,
σ′m 为平均有效应力,n约等于0.5。平均固结度U通常定义为
U=SctSc∞, (4) 式中,Sct为地基某时刻的主固结变形,Sc∞为地基的最终主固结变形。
根据太沙基一维固结理论,平均固结度U为
U=1−∞∑m=12M2e−M2Tv。 (5) 式中,竖向固结时间因子Tv为
Tv=kvtH2γwmv, (6) 式中,kv为土的竖向渗透系数,t为固结时间,H为土层厚度,
γw 为水的重度,mv为土体体积压缩模量。设砂土震后初始孔隙比为e0,主固结完成后孔隙比为ec∞,固结中土体孔隙比为ect,则由式(4)可以得到土体平均固结度:
U=e0−ecte0−ec∞。 (7) 超静孔压比与砂土有效应力的关系如下:
1−ru(t)=[σ′m,tσ′m,i]2, (8) 式中,
ru(t) 为t时刻砂土的超静孔压比,σ′m,t 为t时刻砂土的有效应力,σ′m,i 为砂土的初始最大有效应力。考虑到在主固结阶段砂土剪切模量增长主要由有效应力增加和孔隙比减小引起的,而结构性变化的影响很小。因此本文假设砂土主固结过程中结构性参数A不变,结合式(3),(7),(8)得到砂土震后剪切模量随时间的变化函数如下:
Gpc(t)GR=F(e0−Ue0+Uec∞)F(ec∞)[1−p0σ′(1−U)]0.5, (9) 式中,Gpc(t)为任意固结时刻的砂土小应变剪切模量,GR为砂土在参考时间时的小应变剪切模量,一般参考时间tR大于主固结时间tp。
3.2 次固结影响函数
砂土次固结过程中,颗粒间接触力均质化(包括孔隙比减小)也是引起砂土刚度随时间增长的一个重要因素[20-21]。本文采用Howie等[22]基于室内试验提出的经验公式加以描述:
Gsc(t)GR=1+ΔGGR⋅lg(ttR), (10) 式中,t为固结时间,tR为参考时间,Gsc(t)为任意固结时刻的小应变剪切模量,Gsc,R为在参考时间的小应变剪切模量,
ΔG /Gsc,R为剪切模量增长速率,对于特定的土体需根据室内固结试验确定。3.3 饱和砂土震后小应变刚度计算模型
在震后一定的时间段,砂土剪切模量变化受上述两个因素共同影响,本文假设这两种机理相互独立,结合式(9),(10)得
G(t)GR=1+ΔGpc(t)GR+ΔGsc(t)GR。 (11) 式(11)结合
G=ρv2s 和vs=4Hf 可得到场地卓越频率随随时间发展的函数:f(t)fR=√G(t)GR。 (12) 根据现有围绕基督城场地的研究成果,式(12)的参数取值如下:该场地上覆软弱土层厚20 m,主要由砂土和粉土组成,排水条件良好[15];根据CPT测试的Ic和渗透系数的经验关系[23]取土层平均渗透系数kv=2.5×10-5 m/s;土层的平均压缩模量由动力触探锤击数与压缩模量的经验关系[24]得到,取Es=15.0 MPa;土层初始平均超静孔压比取p0=0.8。另外,假设震后场地刚度需要经历震前同样的沉积时间才能恢复至2010年9月3日强震液化事件之前的震前值,则参考时间取为tR=100 a。
根据式(12)计算得到基督城REHS台站场地液化后基频随时间的发展规律如图8中的曲线所示,图中上、下两条边界线代表由于参数取值带来的标准差。由图8可见,该计算模型对图6(a),(c)中的两次强震数据的预测效果较好,能够考虑震后一定时间内饱和砂土刚度的恢复规律。
4. 结论
为研究震后饱和砂土场地原位测试指标变化规律及其与震前值的差异,本文针对新西兰基督城的典型易液化场地,开展了从实际场地卓越频率监测分析到砂土小应变刚度计算模型的研究。主要结论如下:
(1)针对基督城系列强震,对REHS强震台站记录信号进行HVSR分析,揭示了震后场地刚度整体呈现指数形式的增长规律。场地小应变刚度会在震后经过数天的快速恢复后进入相对稳定的缓慢上升期。对于震前沉积时间较久的场地,场地刚度在震后短时间内难以恢复到震前水平,现有基于震后原位调查获得的液化实例数据应当根据震后原位测试时间进行适当的修正。
(2)假设震后场地刚度需要经历震前同样的沉积时间才能恢复至震前值,基于主固结与次固结过程提出了具有物理意义的饱和砂土震后小应变刚度计算模型,可估算基督城典型易液化场地在强震后平均小应变刚度随时间的发展规律。该计算模型为将特定场地的震后原位测试指标(如剪切波速)修正到对应的震前值提供了一种有效手段。
注意到本文所揭示的震后场地小应变刚度演化规律和所建议的计算模型及参数主要针对新西兰基督城场地,有必要进一步研究更多的其它地震和场地刚度演化规律,并获取各类易液化场地和土体参数,才能为将液化实例中的震后原位测试指标修正到震前原状土的对应值提供定量依据,提高当前基于液化实例调查的地震液化简化判别方法的准确性。
-
表 1 试验土体基本物理性质指标
Table 1 Basic physical properties of test soils
土粒相对质量密度Gs /(g·cm-3) 孔隙比 液限/% 塑限/% 渗透系数k/(cm·s-1) 2.58 1.63 0.82 48 23 2.5×10-6 表 2 不同层次根系的生长数量指标
Table 2 Growth indices of root system
根系指标 时间/d 数量/条 长度/mm 直径/mm 一级根 7 3 98.160 — 10 3 175.320 0.28 14 3 244.060 0.52 21 3 278.100 0.65 28 3 285.980 0.68 二级根 7 11 7.640 — 10 21 15.260 0.16 14 33 17.520 0.36 21 36 17.600 0.47 28 42 17.700 0.50 三级根 7 0 — — 10 45 1.660 0.08 14 91 4.700 0.14 21 134 5.350 0.16 28 168 5.400 0.16 注:表中根系直径为距离根顶100 mm处根系的平均值;二级根和三级根的长度为距离根顶80~120 mm范围内的平均值。 表 3 根系的力学参数指标
Table 3 Mechanical parameters of roots
根径/mm 抗拉力/N 弹性模量/MPa 0.21 0.50 73.88 0.32 1.43 55.72 0.46 2.74 52.82 0.56 3.52 40.21 0.67 5.88 28.79 0.79 7.56 31.11 0.84 9.12 25.41 0.91 11.87 23.98 1.02 15.64 20.59 表 4 室内直剪试验结果
Table 4 Results of direct shear test
生长时间/d c/kPa φ/(°) 0 13.60 25.5 7 13.82 25.3 10 15.32 25.0 14 17.26 24.2 21 18.21 24.8 28 18.59 24.5 表 5 黑麦草的关键生长参数
Table 5 Key growth parameters of Ryegrass
生长参数 s/(mm·d-1) lb/mm la/mm ln/mm n /(rad) c/(mm·d-1) 0th
(种子)30 0 0 0 3 1 1st
(一级根)60 32.2 42 15.3 15 0.218 2nd
(二级根)10 3.2 3.3 2.8 10 0.146 3rd
(三级根)2.0 3.3 2.1 — — — 0.0925 表 6 不同生长时间的根系模型生物量
Table 6 Growth parameters of root model
不同生长时间/d 单元数 根系总长/mm 根系总面积/mm2 根系总体积/mm3 10 589 1445.236 1089.63 116.59 14 1143 2441.007 2090.78 252.34 28 2172 4456.762 3777.45 452.01 表 7 土体参数
Table 7 Soil parameters
天然重度/(kN·m-3) 内摩擦角/(°) 黏聚力/kPa 杨氏模量/Pa 泊松比 14.9 25.5 13.6 2.72×106m2 0.22 表 8 剪切面上各层次根系轴力水平统计表
Table 8 Statistics of axial force of roots at different levels on shear plane
不同生长时间/d 一级根数量/条 一级根平均轴力/N 二级根数量/条 二级根平均轴力/N 三级根数量/条 三级根平均轴力/N 10 3 2.46 3 0.60 6 0.002 14 3 3.87 6 0.78 17 0.100 28 3 3.47 6 1.50 18 0.150 表 9 不同形态根系的根土复合体抗剪强度统计表
Table 9 Statistics of shear strength of root-soil composites with different root forms
生长时间/d 根土复合体抗剪强度 试验结果/kPa 模拟结果/kPa 误差/% 10 15.58 15.84 1.60 14 17.19 18.94 10.18 28 17.98 20.46 13.93 -
[1] 徐华, 李天斌, 周雄华, 等. 高寒地区JYC生态基材护坡现场试验及测试研究[J]. 岩土工程学报, 2009, 31(5): 799–804. doi: 10.3321/j.issn:1000-4548.2009.05.026 XU Hua, LI Tian-bin, ZHOU Xiong-hua, et al. Field tests on JYC ecological base material for slope protection in high-cold areas[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(5): 799–804. (in Chinese) doi: 10.3321/j.issn:1000-4548.2009.05.026
[2] 李天斌, 徐华, 周雄华, 等. 高寒高海拔地区岩质陡边坡JYC生态基材护坡技术[J]. 岩石力学与工程学报, 2008, 27(11): 2332–2339. doi: 10.3321/j.issn:1000-6915.2008.11.022 LI Tian-bin, XU Hua, ZHOU Xiong-hua, et al. Protection techniques of steep rock slope with jyc ecological base materials in high-cold and high-altitude area[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(11): 2332–2339. (in Chinese) doi: 10.3321/j.issn:1000-6915.2008.11.022
[3] XU H, LI T B, CHEN J N, et al. Characteristics and applications of ecological soil substrate for rocky slope vegetation in cold and high-altitude areas[J]. Science of the Total Environment, 2017, 609: 446–455. doi: 10.1016/j.scitotenv.2017.07.156
[4] 周云艳, 陈建平, 王晓梅. 植物根系固土护坡机理的研究进展及展望[J]. 生态环境学报, 2012, 21(6): 1171–1177. https://www.cnki.com.cn/Article/CJFDTOTAL-TRYJ201206030.htm ZHOU Yun-yan, CHEN Jian-ping, WANG Xiao-mei. Progress of study on soil reinforcement mechanisms by root and its expectation[J]. Ecology and Environmental Sciences, 2012, 21(6): 1171–1177. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TRYJ201206030.htm
[5] GONZALEZ-OLLAURI A, MICKOVSKI S B. Plant-soil reinforcement response under different soil hydrological regimes[J]. Geoderma, 2017, 285: 141–150. doi: 10.1016/j.geoderma.2016.10.002
[6] KIM J H, FOURCAUD T, JOURDAN C, et al. Vegetation as a driver of temporal variations in slope stability: the impact of hydrological processes[J]. Geophysical Research Letters, 2017, 44(10): 4897–4907. doi: 10.1002/2017GL073174
[7] NG C W W, WOON K X, LEUNG A K, et al. Experimental investigation of induced suction distribution in a grass-covered soil[J]. Ecological Engineering, 2013, 52(2): 219–223.
[8] 周成, 路永珍, 黄月华. 香根草加固不同含水率膨胀土的侧限膨胀和直剪试验[J]. 岩土工程学报, 2016, 38(增刊2): 30–35. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract16653.shtml ZHOU Cheng, LU Yong-zhen, HUANG Yue-hua. Oedometer expansion and direct shear tests on vetiver root-reinforced expansive soil with different water contents[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(S2): 30–35. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract16653.shtml
[9] 刘亚斌, 胡夏嵩, 余冬梅, 等. 西宁盆地黄土区2种灌木植物根–土界面微观结构特征及摩擦特性试验[J]. 岩石力学与工程学报, 2018, 37(5): 1270–1280. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201805021.htm LIU Ya-bin, HU Xia-song, YU Dong-mei, et al. Microstructural features and friction characteristics of the interface of shrub roots and soil in loess area of Xining Basin[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(5): 1270–1280. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201805021.htm
[10] 王一冰, 吴美苏, 周成, 等. 组合根系加固坡土的直剪试验及数值模拟[J]. 岩土工程学报, 2020, 42(增刊1): 177–182. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18472.shtml WANG Yi-bing, WU Mei-su, ZHOU Cheng, et al. Direct shear tests and numerical simulation on slope soils reinforced by composite roots[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(S1): 177–182. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18472.shtml
[11] 曾红艳, 吴美苏, 周成, 等. 根系与植筋带固土护坡的力学机理试验研究[J]. 岩土工程学报, 2020, 42(增刊2): 151–156. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18539.shtml ZENG Hong-yan, WU Mei-su, ZHOU Cheng, et al. Experimental study on reinforcement mechanism of vegetated slopes with root system and vertical geotextile belts[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(S2): 151–156. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18539.shtml
[12] 孔纲强, 文磊, 刘汉龙, 等. 植物根系分布形态及含根复合土强度特性试验[J]. 岩土力学, 2019, 40(10): 3717–3723. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910004.htm KONG Gang-qiang, WEN Lei, LIU Han-long, et al. Strength properties of root compound soil and morphological observation of plant root[J]. Rock and Soil Mechanics, 2019, 40(10): 3717–3723. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910004.htm
[13] MICKOVSKI S B, STOKES A, VAN BEEK R, et al. Simulation of direct shear tests on rooted and non-rooted soil using finite element analysis[J]. Ecological Engineering, 2011, 37(10): 1523–1532. doi: 10.1016/j.ecoleng.2011.06.001
[14] BERNTSON G M. Modelling root architecture: are there tradeoffs between efficiency and potential of resource acquisition? [J]. New Phytologist, 1994, 127(3): 483–493. doi: 10.1111/j.1469-8137.1994.tb03966.x
[15] XU H, WANG X Y, LIU C N, et al. A 3D root system morphological and mechanical model based on L-Systems and its application to estimate the shear strength of root-soil composites[J]. Soil and Tillage Research, 2021, 212: 105074. doi: 10.1016/j.still.2021.105074
[16] JEWELL R A, WROTH C P. Direct shear tests on reinforced sand[J]. Géotechnique, 1987, 37(1): 53–68. doi: 10.1680/geot.1987.37.1.53
[17] LEITNER D, KLEPSCH S, KNIEß A, et al. The algorithmic beauty of plant roots-an L-System model for dynamic root growth simulation[J]. Mathematical and Computer Modelling of Dynamical Systems, 2010, 16(6): 575–587. doi: 10.1080/13873954.2010.491360
[18] PAGÈS L, VERCAMBRE G, DROUET J L, et al. Root Typ: a generic model to depict and analyse the root system architecture[J]. Plant and Soil, 2004, 258(1): 103–119. doi: 10.1023/B%3APLSO.0000016540.47134.03
[19] MATTIA C, BISCHETTI G B, GENTILE F. Biotechnical characteristics of root systems of typical mediterranean species[J]. Plant and Soil, 2005, 278(1/2): 23–32.
[20] 郝郑芳. 高速公路切方边坡防护草本植物根系固土能力研究[D]. 雅安: 四川农业大学, 2014. HAO Zheng-fang. Research of Capability of Soil Conservation of Herb Plant root Used in Cutted Slope of Highway[D]. Yaan: Sichuan Agricultural University, 2014. (in Chinese)
[21] 沈庆双. 草本植物加固边坡的试验探究[D]. 北京: 中国地质大学(北京), 2018. SHEN Qing-shuang. Experimental Investigation of Slope Reinforcement by Herbaceous Plants[D]. Beijing: China University of Geosciences, 2018. (in Chinese)
[22] 田佳, 曹兵, 金楠, 等. 花棒根-土复合体直剪试验的有限元数值模拟与验证[J]. 农业工程学报, 2015, 31(16): 152–158. https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201516021.htm TIAN Jia, CAO Bing, JI Jin-nan, et al. Numerical simulation and validation test of direct shear test for root-soil composite of Hedysarum scoparium using finite element method[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(16): 152–158. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201516021.htm
[23] 闫海燕. 香根草根土复合体力学性能研究[D]. 重庆: 重庆交通大学, 2013. YAN Hai-yan. Research on Mechanical Properties of Vetiver Root-Soil Composite[D]. Chongqing: Chongqing Jiaotong University, 2013. (in Chinese)
[24] 卜宗举. 植被根系浅层加筋作用对边坡稳定性的影响[J]. 北京交通大学学报, 2016, 40(3): 55–60. https://www.cnki.com.cn/Article/CJFDTOTAL-BFJT201603010.htm BU Zong-ju. Effect of shallow layer of vegetation root on slope stability[J]. Journal of Beijing Jiaotong University, 2016, 40(3): 55–60. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BFJT201603010.htm
[25] WU T H. Investigation of landslides on Prince of Wales Island, Alaska, Geotechnical Engr. Report No 5, dept. of Civil Engr[M]. Columbus: Ohio State University, 1976: 94.
-
期刊类型引用(2)
1. 刘志文,王媛,董琪,高山. 不连续界面渗流的堤防防渗膜防渗效果模拟研究. 人民黄河. 2024(06): 48-53+67 . 百度学术
2. 张升,兰鹏,苏晶晶,熊海斌. 基于PINNs算法的地下水渗流模型求解及参数反演. 岩土工程学报. 2023(02): 376-383+443-444 . 本站查看
其他类型引用(0)