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

《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）表征土壤强度，其核心理念是将施加应力控制在土壤既往承受过的范围内以避免永久变形。然而，传统方法普遍依赖施应力的对数转换，这一操作存在显著缺陷。本研究针对比利时黄土母质底土样品开展限制性单轴压缩试验，采用应变控制方式加载至约800 kPa，直接在线性尺度下计算土壤压缩性（定义为单位应力增量引起的应变增量）。结果显示，在低应力区间，土壤压缩性通常随应力增加而下降；约88%的样品在压缩性曲线中出现局部最小值，该点对应从应变硬化向应变软化的转变，可视为预压应力的物理等价指标，所有试验的中位值为75 kPa。约12%的样品在整个应力区间内压缩性持续下降。此外，58个出现局部最小值的样品中，有33个在更高应力水平出现第二个局部最小值。研究人员强烈建议在压缩试验数据分析中采用线性尺度的应变与应力。本研究表明，土壤未必存在单一的明确预压应力。附加分析证实，对仅含约11个观测值的数据集应用样条插值（spline interpolation），可可靠估算观测点之间的应变值。

研究背景与意义

随着农业机械大型化导致接地压强持续增加，底土压实风险已成为威胁农田生产力与生态系统服务的关键问题。传统土壤力学沿用Casagrande于1936年提出的半对数法（应变-对数应力）确定预压应力σ_p，该方法源自饱和均质黏土的固结理论，后被推广至非饱和结构性农业土壤。然而，对数转换会人为制造弹性材料的“假预压应力”，导致物理意义缺失。同时，农业土壤具有多级结构体（团聚体-大结构单元）与生物扰动形成的复杂异质性，其力学响应与非结构性沉积土存在本质差异。现有模型过度依赖单一σ_p阈值，但田间观测常发现预测值与实测土壤行为不符。为此，本研究提出基于线性尺度的物理机制分析方法，重新评估预压应力的存在性与判定逻辑。

关键技术方法

研究人员在比利时Gembloux的侵蚀棕壤（Eutric Cambisol）长期试验田采集40–44 cm底土层原状土芯，共72个样品，质地为粉壤土（粘粒0.27、粉粒0.71、砂粒0.02 kg·kg^-1）。样品经毛细管饱和后，通过张力台控制至6个基质势（-30至-300 hPa）。采用应变控制式限制性单轴压缩试验，以1 mm·min^-1速率加载至约800 kPa，连续采集应变-应力数据。为验证稀疏数据的可靠性，额外构建仅含11个应力点的简化数据集。数据分析采用样条插值（spline interpolation）估算5 kPa间隔的应变值，通过计算应变-应力曲线的斜率获取土壤压缩性（C = dε/dσ），并识别压缩性曲线的局部最小值作为屈服应力（σ_yield）。橡胶对照试验用于验证对数转换的物理误导性。

研究结果

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

线性尺度下，土壤压缩性随应力增加呈非线性变化。约88%样品（58个）出现首个压缩性局部最小值（中位σ_yield= 75 kPa），标志应变硬化转为应变软化；其中57%（33个）在更高应力处（中位190 kPa）出现第二个局部最小值。约12%样品（8个）在全应力区间压缩性持续下降，无明确转折点。基质势对σ_yield无系统性影响，但在-30 hPa下出现未预期的高值。

3.2 稀疏数据集的压缩性估算

基于11个应力点的样条插值可高精度重建全数据集的应变曲线（平均误差<0.0005 m·m^-1）。对具单一明显最小值的样品，简化数据与全数据集测得的σ_yield高度一致；对具双最小值的样品，插值会部分平滑化曲线细节，但仍能反映基本形态。

3.3 橡胶对照试验

橡胶芯的应变-应力曲线在线性尺度下显示压缩性持续下降，符合其熵弹性物理特性。但在半对数坐标中，传统方法（Gompertz模型、最大曲率法）错误识别出约4 kPa的“预压应力”，证实对数转换会赋予完全弹性材料虚假的力学阈值。

讨论与结论

对数转换的物理缺陷

对数变换将线性弹性响应扭曲为S形曲线，迫使任意土壤数据拟合出“预压应力”，导致橡胶等无塑性变形的材料被误判。这与Wesley等学者对e-log(p)法的批判一致，说明传统方法缺乏物理基础。

土壤屈服行为的复杂性

压缩性局部最小值对应结构体破坏引发的应变软化，但其数量与位置受土壤多级结构控制。单最小值反映历史荷载统一破坏了特定尺度团聚体；双最小值揭示不同强度结构单元的分批屈服；持续下降则表明无主导结构破坏阈值，此时预压应力概念失效。这解释了为何传统σ_p常与土壤属性关联微弱。

方法论革新意义

线性尺度分析法消除了对数转换的人为偏差，使σ_yield与物理过程直接对应。样条插值证明，即使经典试验仅获取少量数据点，仍可可靠识别屈服应力。这为重新分析历史压缩试验数据提供了可行路径。

结论

本研究证实，基于线性尺度的压缩性分析能透明揭示土壤力学行为。多数底土存在至少一个压缩性局部最小值，但12%样品无明确转折点，50%以上具多重转折点，证明单一预压应力概念不适用于所有农业土壤。研究人员呼吁在农田交通承载力评估中放弃单一阈值假设，采用物理机制驱动的线性分析法，并建议对既有压缩试验数据进行重新评估。本研究发表于《Soil Science Society of America Journal》。