基于线性尺度应变-应力数据评估土壤预压应力：亟需物理驱动的方法

《Soil Science Society of America Journal》：Soil precompression stress assessed from strain-stress data in linear scale: The need for a physically based approach

　　

摘要：前沿土壤压实模型以预压应力（precompression stress, σp）表征土壤强度，其核心假设为施加应力不应超过土壤既往承受过的荷载，以避免产生永久变形。然而，传统σp量化普遍依赖施加油门的对数变换，存在显著缺陷。本研究针对比利时黄土母质底土样品开展试验，采用应变控制式加载至约800 kPa，实施限制单向压缩测试。研究人员将应变与应力均按线性尺度处理，以单位应力增量引起的应变增量定义土壤压缩性（compressibility, C）。结果显示，低应力下土壤压缩性随应力升高而降低；约88%的样品中，压缩性出现局部最小值，标志着由应变硬化向应变软化的转变，该转折点可视为σp的物理等效指标，所有测试的σp中位数为75 kPa。约12%的样品在全应力范围内压缩性持续下降，无明确阈值。此外，58个出现压缩性局部最小值的样品中，33个在高应力区出现第二个局部最小值。研究强烈建议在压缩测试数据分析中采用线性尺度的应变与应力，结果表明土壤未必存在单一的明确预压应力。附加分析表明，对仅含约11个观测值的数据集采用样条插值（spline interpolation），可可靠估计观测点间的应变值。

论文解读：线性尺度下土壤预压应力的物理重构与农业压实风险评估的新启示

研究背景与必要性

随着农业机械大型化导致接地压强持续升高，底土压实已成为威胁农田生产力与生态系统服务的关键问题。传统土壤力学沿用Casagrande于1936年提出的半对数法（应变-对数应力曲线）确定预压应力（σ_p），将其定义为弹性变形与塑性变形的分界阈值。然而，该方法最初针对饱和均质黏土层开发，应用于非饱和、结构性复杂的农业土壤时存在根本性缺陷：Keller等已通过理论证明，理想线弹性材料在半对数坐标下也会人为产生“预压应力”，与物理事实相悖。橡胶对照实验进一步证实，完全弹性的橡胶在半对数坐标下同样被误判为有σ_p，揭示了传统方法的系统性偏差。此外，大量研究表明σ_p与土壤理化性质的关联性极弱，难以支撑精准的田间交通风险评估。因此，亟需建立基于物理机制、摆脱对数变换误导的土壤力学分析方法。

技术方法概述

研究人员依托比利时Gembloux长期定位试验田，采集12个小区40–44 cm深度的底土样品，共获得72个原状土芯（100 cm3）。样品经毛细管饱和后，通过张力台与真空陶瓷板分别平衡至-30、-50、-75、-100、-160、-300 hPa共6个基质势，每个小区对应一个基质势。采用应变控制式限制单向压缩测试，以1 mm min?1速率加载至约800 kPa，连续采集56组应变-应力数据。核心分析流程包括：将应变与应力统一置于线性尺度，以5 kPa间隔通过样条插值补全数据；计算压缩性C为应变对应力的导数（C=dε/dσ），采用二次多项式滑动回归与点对点线性回归两种方法验证；定义压缩性局部最小值对应的应力为屈服应力（σ_yield）；为验证稀疏数据的适用性，特意构建仅含11个观测值的简化数据集，对比全数据与简化数据集的σ_yield差异，以评估插值可靠性。

研究结果

3.1 应变-应力关系与压缩性特征

样品初始容重平均为1.53 g cm?3，孔隙度为0.43 m3 m?3，测试前饱和度介于0.776–0.858之间。应变随应力升高而增加，最大应变介于0.1–0.2之间。压缩性C随应力变化呈现三类模式：第一类（占比约88%）表现为低应力下C持续下降，达局部最小值后转为上升，对应σ_yield中位数为75 kPa；第二类（占比约12%）在全应力区间C单调下降，无局部最小值；第三类（占第一类的57%，即总样本的50%）出现两个局部最小值，第一个中位数75 kPa，第二个中位数190 kPa。基质势对σ_yield无显著影响，仅-30 hPa下略高，可能与个别样品特性有关。

3.2 稀疏数据集的可靠性验证

对仅含11个观测点的简化数据集进行样条插值，其应变估计值与实测值的平均误差小于0.0005 m m?1，精度极高。简化数据集能准确识别第一类样品的单一σ_yield，但对第二类双最小值样品，插值会平滑曲线，使两个极小值合并为一个中间值，或弱化其区分度。31个样品的全数据与简化数据集σ_yield呈显著正相关，证实传统稀疏采样数据可通过插值有效提取线性尺度下的力学特征。

3.3 橡胶对照的证伪实验

橡胶样品的应变-应力关系呈非线性，压缩性C随应力升高持续下降，符合弹性材料的熵弹性特征。但在半对数坐标下，其曲线呈典型的S型，传统Gompertz模型与最大曲率法均给出约4 kPa的“预压应力”。这一结果直接证明传统半对数法的物理谬误——它无法区分真实的结构破坏阈值与坐标变换带来的数学假象。

讨论与结论总结

传统半对数法的根本问题在于将经验公式误作物理定律，导致农业土壤σ_p估算普遍存在偏差。基于阻力概念（resistance concept，即单位作用产生的响应，此处为应力/应变）的线性尺度分析，将压缩性定义为C=dε/dσ，其局部最小值对应着土壤结构的屈服：低应力下颗粒重排致密化（应变硬化，压缩性下降），超过历史荷载后原生结构单元破坏（应变软化，压缩性上升）。但研究首次发现，相当比例的土壤并无单一阈值，而是呈现多阶屈服或持续硬化，这与土壤结构的层级性（大团聚体由更小单元构成，强度各异）直接相关，也解释了为何传统模型常低估实际压实风险。

结论部分明确指出：线性尺度分析是评估土壤力学行为的透明物理途径，阻力概念应成为未来田间交通临界应力研究的基础。约88%的样品存在由应变硬化向软化转变的σ_yield，但其中过半存在多阶屈服，且12%的样品无明确阈值，证明“单一预压应力”概念在结构性土壤中并不普适。稀疏数据经样条插值可有效识别σ_yield，建议重新分析历史压缩试验数据。该研究发表于《Soil Science Society of America Journal》，为农业土壤压实风险评估提供了去伪存真的物理框架，推动模型从经验拟合转向机理驱动。