基于堆积试验的黏壤土仿真物理参数标定

为获取南方黏壤土的离散元仿真模型精准接触参数,构建土壤离散元仿真模型,基于土壤堆积试验,结合试验测定和EDEM软件推荐的参数构建土壤仿真模型,以休止角为响应值,采用Design Expert软件依次设计Plackett-Burman试验、最陡爬坡试验和Box-Behnken试验完成土壤仿真物理参数标定及优化。进一步通过成穴装置成穴的仿真试验与土槽试验的对比分析,验证黏壤土仿真模型的精准性。成穴试验的仿真与土槽试验结果表明,装置在仿真土壤中运动规律与土槽试验中差异较小,成型穴孔的开口纵长和有效深度的误差分别为3.98%和1.87%,模拟仿真土壤的物理力学特性与实际土壤一致,表明黏壤土的离散元仿真参数系统标定研究的方法准确可行。研究构建了南方黏壤土精准的离散元仿真模型,为该类型土壤与触土部件相互作用的动力学研究提供技术支持。     

煤矸石振陷变形特性试验研究

动荷载作用下累积轴向残余应变是煤矸石路基发生振陷变形的主要原因,利用自燃煤矸石原料制作圆柱形试件进行动三轴试验,研究了不同相对密度Dr、动荷载幅值σ_(dmax)条件下煤矸石累积轴向残余应变ε_(pa)和动模量E_d特性,建立了考虑相对密度、动荷载幅值二因素的累积轴向残余应变和动模量公式。试验结果表明:动应力小于临界动应力σ_(dcr)时,累积轴向残余应变与动荷载振次对数值log N呈线性关系,可用公式ε_(pa)=a+b log N进行计算;动应力大于临界动应力σ_(dcr)时,累积轴向残余应变呈加速增长形式;动模量随累积轴向残余应变增长而降低,呈双曲线形式;初始动模量随相对密度增加而增大。研究结果为煤矸石路基振陷变形分析提供了试验、理论依据。     

不同含水率黏重黑土与触土部件互作的离散元仿真参数标定

为了获得可用于东北地区黏重黑土与触土部件相互作用的离散元仿真模拟参数,该文利用EDEM中Hertz-Mindlin with JKR Cohesion接触模型对不同含水率的东北地区黏重黑土进行相关参数标定,针对含水率在10%～20%的实际作业环境,分别配置含水率为12.46%±1%和17.15%±1%的2种黏重黑土,以土壤颗粒间的滚动摩擦系数、恢复系数、表面能参数及静摩擦系数为标定对象,并以土壤颗粒的仿真堆积角为响应值,基于Box-Behnken的响应面优化方法得到堆积角回归模型,并对回归模型进行寻优,得到2种含水率的模型参数优化解,并给出了模型参数范围。测定了4种含水率下黏重黑土对3种触土部件材料(65Mn、UHMW-PE和PTFE)的静摩擦系数,并以此为基础分别对65Mn(典型铁基材料)和PTFE(典型低表面能材料)板进行斜面试验,以含水率为17.15%±1%的黏重黑土为试验对象,分别搭建斜面物理试验平台和仿真模型,以土壤颗粒与触土部件材料之间的滚动摩擦系数、恢复系数、表面能参数及静摩擦系数为标定对象,以仿真得到的土球在65Mn和PTFE板上的滚动距离为响应值,基于响应面优化法得到滚动距离回归模型,以实测的滚动距离为目标对回归模型进行寻优,得到黏重黑土对2类典型触土部件材料接触模型参数的优化解。研究结果表明,标定优化后的土壤模型能够近似代替真实的东北地区黏重黑土进行仿真,可利用标定后的参数进行黏重黑土与触土部件间的离散元仿真,可为东北黏重黑土作业条件下的农业机械触土部件仿生减阻设计与优化提供基础数据。     

基于离散元非线性弹塑性接触模型的免耕土壤参数标定

为了进一步提升基于离散元法对免耕作业机具与土壤互作关系研究的准确性,以华北麦玉两熟区免耕壤土为研究对象,基于离散元软件（Discrete Element Modeling,DEM）扩展的The Edinburgh Elasto-Plastic Adhesion非线性弹塑性接触模型开展常年免耕农田土壤离散元仿真模型参数标定研究。应用Plackett-Burman 设计敏感性分析试验,选择对土壤压板试验沉陷量影响显著的关键参数（即土壤颗粒恢复系数、土壤颗粒-颗粒静摩擦系数、土壤颗粒表面能和土壤颗粒-颗粒塑性变形比）,以土壤压板试验沉陷量为评价指标,应用Box-Behnken 试验建立了沉陷量与4个显著性参数的二次多项式回归模型,以物理试验得到的实际沉陷量6.2 mm为目标值,对显著性参数进行寻优,得到最优组合为：土壤颗粒恢复系数为0.47,土壤颗粒-颗粒间静摩擦系数为0.62,土壤颗粒表面能为6.12 J/m2,土壤颗粒塑性变形比为0.41。最后通过标定优化的参数进行土壤颗粒应力传递试验,与室内试验结果对比发现,仿真的土壤颗粒接触传力特性与实际土壤压实过程中应力传递差别较小,误差范围在9.21%以内,并对比分析仿真和实测的土壤应力传递特性曲线的拟合情况,两者之间的决定系数R2为0.91,表明两条曲线的拟合相似度较高,验证了免耕土壤参数标定的准确可靠,参数标定后的离散元免耕土壤模型精确度高。研究可为快速构建免耕模式下的土壤离散元模型及免耕农机装备优化提供技术支撑。     

面向夹持采收的油菜薹夹段茎秆离散元参数标定与优化

为提高离散元法对指导油菜薹有序采收装备设计与优化的准确性和可靠性,该研究以双行垄作移栽的“农大1号”双低甘蓝型油菜机械化适收期油菜薹夹段茎秆为对象,测定其本征参数、表面接触参数以及破碎力学参数,利用EDEM仿真软件Hertz-Mindlin无滑移模型和Hertz-Mindlin with bonding粘结模型建立夹段茎秆堆积仿真标定模型和破碎仿真标定模型。采用逐步调整仿真参数使仿真试验值与物理试验值逼近的方法,利用夹段茎秆堆积仿真标定模型,以休止角仿真试验值与实际物理试验值的相对误差为目标,完成夹段茎秆表面接触参数的标定与优化;利用破碎仿真标定模型,以轴向压缩和弯曲仿真试验与实际物理试验的最大轴向压缩力和最大弯曲力的相对误差为目标,利用标定后的表面接触参数完成夹段茎秆粘结参数的标定与优化。最后,利用夹段茎秆的径向压缩与剪切、内芯与表皮拉伸的破碎仿真力学试验和有序采收EDEM-Recurdyn耦合仿真试验验证标定后的表面接触参数和粘结参数。结果表明,仿真与实测试验的破碎力学参数相对误差在5% 以内,且仿真与实测的“时间-载荷”曲线变化趋势一致,低速、中速和高速档的有序采收仿真试验结果与实际物理试验结果相对误差在7.0% 以内。研究结果表明,采用离散元仿真方法研究油菜薹采收过程具有可行性,标定结果可用于指导油菜薹机械化生产。     

土壤-水稻前茬离散元仿真参数标定及其旋耕轨迹分析

在进行稻油轮作区域水稻收获后旋耕作业中的前茬运动、受力仿真分析时,由于缺乏可靠的离散元仿真参数难以准确计算土壤-水稻前茬运动规律和相互间受力。该研究采用EDEM软件对土壤-水稻前茬混合物内部各物质的物理与接触参数进行离散元仿真标定。首先测定了土壤-水稻前茬混合物内各物质物理与相互间接触参数,以堆积角作为评价指标,利用仿真试验分析进行了显著因素筛选、最优水平取值范围确定及二阶回归模型下最佳参数寻优对参数进行了标定,得到土壤-水稻前茬混合物模型极显著参数的最优组合为土壤-土壤恢复系数0.407、土壤-土壤动摩擦系数0.123、土壤-水稻前茬静摩擦系数0.596、土壤表面能1.860 J,此时与试验堆积角误差为0.58%。并进行了土壤-水稻前茬间的直剪试验与仿真,得到了稳定后不同垂直应力下试验与仿真的直剪切应力误差分别为：5.4%（50 kPa）,4.1%（100 kPa）,3.1%（150 kPa）。最后在最优参数组合下,开展了基于MBD-DEM的旋耕土壤-水稻前茬轨迹分析试验,分析得到前茬跟随旋耕刀及直接被掩埋两种埋覆场景,仿真与土槽试验轨迹的Spearman秩相关系数分别为0.91,0.84,验证了仿真模型标定参数的准确性。研究结果可为水稻前茬旋耕轨迹调控过程的离散元仿真分析提供可靠的接触参数。     

小麦籽粒振动筛分黏弹塑性接触模型构建及其参数标定

为更好地模拟小麦籽粒(麦粒)振动筛分过程,该文在现有离散元接触模型基础上,通过EDEM软件应用程序编程接口,构建了一种麦粒黏弹塑性接触模型。该接触模型法向方向通过将Kuwabara and Kono非线性黏弹性接触模型中的黏性耗散项引入Thornton滞回接触模型进行构建;切向方向采用简化Thornton切向接触模型;滚动摩擦力矩计算同Hertz-Mindlin(no slip)接触模型。法向模型参数标定采用单轴加载-卸载试验、碰撞试验,分别构建了模拟-试验接触力误差平方和与麦粒屈服重叠量、恢复系数与法向阻尼系数的二阶回归方程,得到麦粒屈服重叠量为7.63×10-6 m,麦粒-麦粒/钢板法向阻尼系数分别为190.68和306.65。切向模型参数利用旋转鼓试验进行标定,得到最佳麦粒-麦粒/钢板静摩擦系数组合为0.40和0.44。最后利用振动筛分试验对所标定参数进行验证,模拟与试验所得筛下麦粒质量分数最大误差为8.97%,模拟中筛下物分布规律与试验结果无显著性差异,表明所建立的接触模型及标定的参数能够很好地模拟麦粒振动筛分过程。该文也可为其他农业物料黏弹塑性接触模型构建及其参数标定提供参考。     

葡萄藤防寒土与清土部件相互作用的离散元仿真参数标定

为系统地研究中国北方地区沙壤土质地的葡萄藤防寒土及其与清土机清土部件常用材料(Q235钢、橡胶)相互作用的离散元仿真参数,以构建准确的土壤离散元仿真模型,该文选用整合延迟弹性模型(hysteretic spring contact model,HSCM)和线性粘附模型(liner cohesion model,LCM)作为土壤颗粒间的接触模型;基于土壤堆积试验,以土壤颗粒间恢复系数、静摩擦系数、滚动摩擦系数和土壤粘附能量密度为因素,以土壤堆积角为指标,利用EDEM进行通用旋转中心组合模拟试验,采用Design-Expert软件对试验数据进行回归分析,以实测的土壤堆积角作为优化目标值,获得土壤颗粒间的最佳接触参数组合;利用土壤屈服试验获得HSCM模型参数;基于斜面滑动法原理,利用倾斜板试验台测得土壤与Q235钢和橡胶之间的静摩擦系数,并以此为基础,采用土壤滑落试验,以滑动摩擦角为响应值,对土壤颗粒与Q235钢和橡胶之间的恢复系数和滚动摩擦系数进行寻优,得到最优解参数组合。为验证标定优化的离散元模型参数的准确性,采用刮土板土槽试验和仿真试验进行对比分析,获得刮土板在土槽试验和仿真试验中的水平前进阻力分别为228.36 N和213.79 N,两者之间的相对误差为6.38%,表明仿真模型中土壤的物理力学特性与实际土壤基本一致,验证了葡萄藤防寒土离散元仿真参数标定结果和研究方法准确可靠。研究结果可为基于离散元法研制适用于北方地区沙壤土质地的葡萄藤防寒土清土机提供理论基础和技术支撑。     

基于准静态压缩的果实黏弹塑性模型

目前对果实黏弹塑性特征与模型进行表征比较困难。基于果实等速压缩变形特征，提出了一表征其黏弹塑性的流变模型。该模型包括3个非线性弹簧元件、一个黏性元件以及一个滑块元件。结合压缩应力—应变试验曲线特征与模型特性提出了模型参数的数值确定方法。采用梨果实全果肉、带核果肉两种试样进行试验，获得了其黏弹塑性模型参数。结果表明试验与理论模型结果非常吻合，果核对其黏弹塑性有较大影响。采用本模型结合试验结果能较准确简便地确定类似特性果实的流变特征参数。研究结果为果品黏弹塑性流变特性分析以及模型表征提供了一种新的方法。     

新疆棉田耕后土壤模型离散元参数标定

为提高耕后分层施肥开沟覆土过程离散元仿真模拟的准确性,采用EDEM离散元软件对分层施肥作业土壤的堆积和滑落过程进行仿真模拟,来标定土壤接触参数。通过通用旋转中心组合试验,采用 Design-Expert 软件对试验数据进行回归分析,以实测土壤休止角、土壤与65 Mn钢滑动摩擦角为优化目标,获得最优的离散元接触参数组合为：土壤间恢复系数0.48、土壤间滚动摩擦系数0.56、土壤间静摩擦系数0.24、土壤与65 Mn钢间恢复系数0.5、土壤与65 Mn钢间滚动摩擦系数0.1、土壤与65 Mn钢间静摩擦系数0.31。为验证标定优化的离散元模型参数的准确性,对土壤堆积试验和滑落试验进行仿真试验与实际试验对比,两者相对误差分别为1.7%和2.5%;并在最优标定参数组合条件下,采用离散法仿真模拟分层施肥装置的开沟覆土过程,获得分层施肥装置5、6和7 km/h作业速度下,仿真试验和田间试验的工作阻力相对误差分别为10.2%、7.95%、7.04%,误差在可接受范围内。仿真试验和田间试验开沟覆土效果基本一致,验证了土壤参数标定的准确可靠,可为后期分层施肥装置减阻研究提供理论基础和技术支持。     

有机质浸染砂水泥土的力学特性及本构关系

为研究海南省海湾地区分布的有机质浸染砂水泥土的力学特性,该文首先对有机质浸染砂水泥土和标准砂水泥土进行无侧限抗压强度对比试验(采用配合比均为熟石灰掺入比7.5%、水泥掺入比20%、水灰比0.45),定量分析养护龄期对其无侧限抗压强度及试样破坏形式的影响。然后对有机质浸染砂水泥土进行单轴抗压试验,获得了水泥土材料的应力-应变全过程曲线、刚度变化规律以及改进邓肯-张本构模型。结果表明:1)有机质浸染砂水泥土试样的破坏类型为塑性剪切破坏和脆性剪切破坏;2)有机质浸染砂水泥土抗压强度随着养护龄期的增加基本呈指数形式增长,但在养护龄期14 d后,增长速度逐渐降低并趋于稳定;3)随着养护龄期增长,水泥土刚度增加。在加载初期,水泥土切线模量随着轴向应变增加而增大,呈现刚度硬化现象;4)基于单轴抗压试验得到应力-应变全过程曲线,可分为2个阶段:塑性阶段、软化阶段;5)通过对应力-应变全曲线的描述,得到了修正的邓肯-张模型,确定模型参数后,与实测数据对比发现,该修正模型可以模拟有机质浸染砂水泥土的应力-应变关系。     

Mathematical simulation and calculation of the soil compaction under dynamic loads

The deformation and compaction of loamy sandy soddy-podzolic soils under linear dynamic changes in the compressive stresses and in the course of the soil creeping were studied in field experiments. The rheological properties of these soils occurring in the viscoelastic state were described by a first-order differential equation relating the compressive stresses, the rates of their changes, and the velocities of the relative vertical compressive deformation. Regression equations were derived for the viscoelastic properties of the studied soil as functions of its density, moisture, and linear compaction velocity. Methods were proposed for the calculation of indices of the stress-strain state and the compaction of soils under specified conditions of changes in their compressive stresses with time and in the course of the soil creeping after the initial linear increase in load. Corresponding computer programs were developed. The effect of the main factors due to the linear increase in the compressive loads and in the course of the soil creeping on the rheological properties, the stress-strain state, and the density of soils was quantitatively estimated. The calculation showed that the values of the soil deformation and the density under compressive stresses lower than the ultimate strength were stabilized with time, and the properties of the viscoelastic soil approached elastic ones.     

木质素纤维改良膨胀土的冻融特性及微观机理

高寒区复杂环境场是渠基膨胀土劣化的主要原因，渠基膨胀土改良是保障渠系工程安全运行的重要手段。为探索冻融循环下膨胀土的纤维改良效果，以北疆某引调水工程为背景，开展木质素纤维改良膨胀土的系列室内试验，探明冻融循环下改良膨胀土的体变特性、抗压强度特性、抗拉强度特性及其作用机理。结果表明：相较于未改良试样，冻融循环下改良试样的体积变化率减小了24%～37%；在15次冻融循环后，未改良试样抗剪强度衰减率在42.2%以上，而改良试样抗剪强度衰减率则显著低于29.1%；反复冻融作用并未显著影响试样的压应力-应变关系曲线特征，而是造成试样拉应力-应变关系由“线弹性增长-强化-软化”特征发展为“线弹性增长-软化”特征；木质素纤维的掺入可以有效改善冻融循环下膨胀土的体变特性、抗剪强度特性、抗压强度特性、抗拉强度特性，且当掺量和长度分别为2%和1 mm时，木质素纤维改良膨胀土冻融特性的作用效果最为显著。同时，冻融循环下改良试样内部结构损伤程度显著低于未改良试样，原因在于木质素纤维的存在限制了冻融循环下土颗粒间的错动、重分布，在一定程度上保证了膨胀土内部结构完整性。     

养护龄期对改良膨胀土无侧限抗压强度试验的影响

以湖北省宜昌市某公路膨胀土为研究对象,分别采用不同掺量的石灰、水泥、粉煤灰、风化砂对膨胀土进行改良。在经过7,14和28d标准条件养生后,进行无侧限抗压强度试验。试验结果表明,上述4种材料均能有效地提高改良膨胀土的无侧限抗压强度,在最初的14d内,水泥改良膨胀土的强度增长较为明显,后期随着养护龄期的增长,无侧限抗压强度增长速率较缓慢;石灰改良膨胀土的无侧限抗压强度值随着养护龄期的增加,一直保持近似直线形的增长趋势;粉煤灰改良膨胀土的无侧限抗压强度亦随着养护龄期的增长而呈直线增长,但其值小于石灰改良膨胀土;养护龄期对风化砂改良膨胀土的无侧限抗压强度影响很小。     

湖南邵阳地区高液限红黏土干缩裂隙演化过程的量化分析

红黏土对环境湿度变化非常敏感,在干燥环境中极易开裂,纵横交错的裂缝网络损害了土壤结构的完整性,很容易诱发红黏土边坡失稳和崩塌,导致农田水利设施的破坏,甚至加剧整个生态环境的水旱灾害。为探究红黏土裂隙的演化规律,该研究采用自制试验装置和三维应变测量系统开展了自然湿热条件下的红黏土泥浆样干燥试验,通过采集土体水分和土表位移、应变和裂隙的变化,定量分析脱湿过程中土体表面裂隙形态和应变场的演变特征,并进一步探讨水分变化对裂隙形态和应变场的影响。结果表明：1）土样表面干缩裂隙的演化一共经历了6个阶段,后阶段裂隙分割前阶段裂隙围成的区域,且不同阶段裂隙的交叉角接近90o;2）裂隙产生初始,裂隙尖端处拉应变约为0.5%,土表面大部分区域处于受拉状态;随着裂隙进一步发展,裂隙周边土体逐渐由拉应变状态向压应变状态转变;当所有裂隙发育完成,裂隙周边土体处于压应变状态;3）裂隙演化阶段与界限含水率有关,当泥浆土样的含水率接近液限时（67.7%）,土体表面裂隙开始发育,裂隙迅速张开和延伸;当土的含水率达到塑限时（28.3%）,裂隙发展的速率逐渐变缓;当含水率小于缩限时（18.8%）,裂隙变化已经很小,裂隙发育接近完成;4）在裂隙演化过程中,早期裂隙的发展持续时间和裂隙宽度均超过后期裂隙;土表不同位置的位移和应变均不相同,土块中心竖向收缩大于边缘竖向收缩,而土块中心位移及应变均小于土块边缘,研究可为红黏土开裂引发的工程地质灾害的预防及治理提供参考。     

装配式地下粮仓钢板-混凝土组合仓壁轴压受力性能分析

为了研究装配式地下粮仓钢板-混凝土组合仓壁的轴压力学性能,该研究在仓壁试件轴压试验的基础上,对仓壁试件进行非线性有限元分析,模拟试件加载的全过程,进一步分析试件及其组件在加载过程中的受力性能及工作机理,并对钢板强度、混凝土强度、距厚比等不同的参数影响规律进行了分析。结果表明:有限元模拟结果与试验结果吻合较好,在最大试验载荷5 000 kN时模拟轴向压缩总变形与试验平均轴向压缩总变形相对差值为4.2%,建立的有限元模型适用;在弹性阶段仓壁预制块部位混凝土和接头部位止水钢板分别承担了79.7%和50.9%的荷载,峰值荷载主要由混凝土和传力钢板决定,达到峰值荷载后接头部位止水钢板承担更大的荷载,增强止水钢板可以改善试件的延性;相较钢板强度、距厚比和止水钢板厚度,混凝土强度对仓壁试件的初始刚度和峰值荷载影响最大,混凝土强度、钢板强度、距厚比和止水钢板厚度对峰值荷载回归得到的回归系数值分别为0.910、0.154、-0.005和0.301;止水钢板的强度和厚度较小时,试件易发生脆性破坏;结合设计参数分析中得到的荷载-变形曲线,提出2种荷载-变形模型曲线,进一步提出装配式钢板-混凝土组合仓壁轴压峰值荷载简化计算式,得到的计算峰值荷载与有限元峰值荷载相对差值均不超过9%,计算结果具有较高的精度。研究结果可为装配式地下粮仓仓壁的工程设计提供指导和参考。     

基于流变学法研究容重和含水率对土壤结构力学稳定性的影响

土壤结构力学稳定性不仅与土壤质量和肥力密切相关,而且还与农业器具设计和农业水土工程建设紧密联系。该研究以黄土高原广泛分布的塿土和黄绵土为研究对象,采用振幅扫描试验模拟振荡荷载过程,研究土壤容重和含水率对土壤结构力学稳定性的影响。结果表明：1）随土壤容重的增加,土粒间接触点增多,使得剪切强度参数：线性黏弹区的剪切应力、屈服点的剪切应力、最大剪切应力及屈服点的储能模量均增加,土壤结构强度增强;其在较大剪切应变下的弹性和总的弹性（屈服点的剪切应变和积分z）先增大后降低,表明土壤颗粒存在一个最稳定的排列组合方式。同时,剪切强度参数对土壤容重的响应更为敏感。2）随含水率的增加,土壤颗粒间黏聚力和摩擦力降低,剪切强度参数：线性黏弹区的剪切应力、屈服点的剪切应力、最大剪切应力及屈服点的储能模量均降低,土壤结构强度降低;当剪切应变在线性黏弹区时,弹性随含水率的增加而增加,而较大的剪切应变下的弹性和总的弹性随含水率增加均降低。3）对比2种土壤,因塿土的黏粒、有机质、阳离子交换量、比表面积等较高,增加了土粒间的胶结强度,使得塿土的弹性和剪切强度较高,而黄绵土结构更具脆性。该研究结果表明基于流变仪中的振幅扫描测试所获取的流变学参数能够定量表征土壤细观结构力学稳定性,为进一步深入认识土壤微观力学特性提供了丰富的评价参数。     

An expression for wind erosion rate from dry and dense soil surface

Wind erosion is a physical process predominated by airborne grains rather than the wind itself. The soil deformation is either elastic–plastic or fully plastic during its collisions with solid grains, corresponding to dense and loose soils, respectively. Only the fully plastic deformation was previously taken into account in the physically-based wind erosion models. The impact erosion of dry and dense soil surface can be well quantified by the contact mechanics and strength theories. In this study, with the help of Owen's uniform saltation model, a simple expression for wind erosion rate is derived based upon the eroded volume per impact when soil deforms in the elastic–plastic manner. The new expression, a simple analytical function of several conventional variables and parameters of wind and soil, describes the effects of saltation flux and soil susceptibility. Its validity is verified by some field and laboratory experiments.     

添加砾石对崩岗岩土无侧限抗压强度的影响

崩岗是中国南方红壤地区常见的一种土壤侵蚀类型，该研究借鉴砾石影响降雨入渗、产流产沙和力学性质的研究成果，对砾石质量分数、形状和直径三因素进行正交设计，研究不同组合砾石对崩岗岩土无侧限抗压强度的影响。崩岗4层土体较好的处理分别为：淋溶层A：A3B1C1，质量分数15%，直径2～4 mm，圆砾；黏化层Bt：A3B2C1，质量分数15%，直径5～7 mm，圆砾；淀积层B：A3B2C1，质量分数15%，直径5～7 mm，圆砾；母质层C：A3B2C1，质量分数15%，直径5～7 mm，圆砾。4层土体的轴向应力随轴向应变均呈急剧上升、急剧下降、减速衰减和衰减稳定4个阶段，但砾石复合土高于未加砾石土。4层土体在较好处理下的无侧限抗压强度分别比未加砾石土提高59.56%、71.70%、49.51%和83.64%，且二者呈线性递增函数关系（R2＝0.99），添加砾石的土柱在受压时破坏程度较小。本研究将为崩岗侵蚀预防和分层治理提供一定的理论依据。