横向超声振动辅助切削的减阻机制与土壤扰动：基于标定离散元模型的参数效应、耦合优化与接触层面起源研究

《Agriculture》：Effects of Transverse Ultrasonic Vibration on Drag Reduction and Soil Disturbance in Soil Tillage: A DEM Investigation

　　

土壤耕作是能耗密集型作业，土壤耕作部件（soil-engaging tools）的切削阻力（cutting resistance）是其最大的单项损耗。超声振动辅助切削（ultrasonic vibration-assisted cutting，UVAC）能够降低该阻力，但定量研究集中于纵向振动（longitudinal vibration）；横向振动（transverse vibration）通过循环侧向挤压与释放对土壤施加载荷，其减阻行为在很大程度上仍未得到探索。研究人员建立了横向超声振动辅助切削工具–土壤系统的标定离散元模型（discrete element method，DEM），用以确定振幅（amplitude）、频率（frequency）与切削速度（cutting speed）如何支配减阻（drag reduction），并识别该效应在接触层面的起源。20 kHz台架试验（bench tests）以额外的作业条件补充了模拟设计空间，并在八个工况——五个振动条件与三个基线——下证实了预测力，平均相对误差为8.7%，基线误差在6%以内；全程使用同一套标定参数。在30–70 μm、25–45 kHz与0.03–0.07 m·s^?1范围内，振幅起主导作用，使减阻率从44.10%提升至68.18%，而频率主要通过其二次项（quadratic term）及其与振幅的交互作用（interaction）起作用，切削速度则进入同一耦合关系。Box–Behnken优化将最优点确定在70 μm、25 kHz与0.07 m·s^?1处，离散元验证的减阻率为69.25%，而模型预测值为69.53%。接触分析（contact analysis）表明，振动使荷载由法向压缩（normal compression）向切向剪切（tangential shear）转移，并将键断裂（bond failure）由集中断裂（concentrated rupture）转化为渐进降解（progressive degradation）。因此，在含水率9%的砂壤土（sandy loam）中，振幅是横向超声振动辅助切削的主要设计变量，频率与速度起调节作用。

土壤耕作是作物生产中能量消耗最大的田间作业，其能耗主要由土壤耕作部件的切削阻力构成，直接影响燃油消耗、CO₂排放与机具磨损。超声振动辅助切削虽已在精密加工领域证实可显著降低切削力，但农业耕作领域的定量研究集中于纵向振动模式；横向振动以循环侧向挤压–释放方式加载土壤，其减阻规律、参数耦合机制与微观接触起源尚未定量建立。Qiao等前期已标定超声振动–刀具–土壤体系的离散元参数，但未比较振动方向，也未系统考察振幅、频率与切削速度的交互作用。上述空白导致农业超声振动工具的参数选择仍主要依赖经验，因此有必要开展横向振动减阻的系统性定量研究。研究人员以标定离散元模型为核心手段，结合单因素模拟、响应面优化、台架验证与接触层面机理分析，得出振幅为砂壤土耕作主导设计变量、横向振动使荷载由压缩向剪切转移且破坏由集中断裂向渐进降解转变等核心结论，为节能型超声耕作机具设计提供了理论依据与数值指导。

关键技术与方法方面，研究人员采用离散元法（EDEM 2024）开展模拟：以Hertz–Mindlin with Bonding接触与键合模型描述颗粒间力学行为，键合球粒直径4 mm、约17万个颗粒生成床高0.11 m的砂壤土试样；通过单因素模拟（振幅30–70 μm、频率25–45 kHz、切削速度0.03–0.07 m·s^?1各五水平）筛选参数范围，并以15次运行的Box–Behnken设计拟合二次多项式回归模型，经Design-Expert 13.0进行方差分析与数值优化；土壤–土壤接触参数与键合参数直接采用前期标定结果，土壤–刀具接触系数在试验前调整并保持不变。土壤样本于2026年3月采集自青岛农业大学试验场（36.3193° N，120.3960° E）0–15 cm土层五点混合样（2.5 kg），风干过2 mm筛后调湿至9%含水率，分五层装填压实；通过20 kHz楔形钛合金刀片台架切割试验与S型拉压传感器测力对模拟进行独立验证。

3.1 横向振动下的宏观减阻
通过单因素模拟与台架试验对比得出：振幅增大使Z轴切削力终值由约50 N降至29 N，减阻率由44.10%单调升至68.18%；频率在25–45 kHz内各水平减阻率介于51.45%–52.69%，波动仅1.24个百分点；切削速度影响较弱，最大减阻率（53.07%）出现于0.04 m·s^?1，高速时趋近常规切削基线。台架试验在五个振动条件与三个常规基线共八个工况下的平均相对误差为8.68%（3.00%–13.56%），基线误差在6%以内，证实标定参数在物理切削条件下的适用性。

3.2 减阻机理
3.2.1 切削力动态与颗粒速度响应：通过力–时间/深度曲线与颗粒速度场分析得出，振幅是减阻强度的直接控制因素，且25–45 kHz下颗粒运动基本一致；切削速度通过"相同深度累积振动循环数（fd/v_c）"调节减阻——低速时同一土体承受更多循环，高速区范围更大，高速时减阻效果衰减。3.2.2 接触力重分布与结构弱化：通过颗粒间法向/切向接触力与键状态、接触数演化分析得出，振动使法向接触力持续低于基线（压缩减少）、切向接触力持续高于基线（剪切增强），振幅越大转移越充分；常规切削呈现最快键断裂，横向振动则保留更多完整键，破坏模式由集中破坏转向渐进降解；但振幅达70 μm时键数下降呈非单调；切削速度对键存活数无可分辨影响（差异小于0.3%）。

3.3 多因素交互作用与最优参数组合
3.3.1 回归模型与方差分析：通过二次多项式回归与方差分析得出，模型极显著（F=173.38，p<0.0001，R²=0.9968，R²_adj=0.9911，R²_pred=0.9489）；振幅主效应（p<0.0001）与振幅–速度交互（p=0.0056）极显著，振幅–频率交互（p=0.0375）与频率二次项（p=0.0380）显著，频率线性主效应（p=0.0571）与速度线性主效应（p=0.1773）不显著。3.3.2 响应面分析：通过三因素响应曲面得出，响应沿频率与速度轴方向平坦、沿振幅方向陡峭，证实振幅为主导变量；振幅–速度正交互表明高速时单位振幅增量带来更大减阻增幅。3.3.3 最优参数组合与验证：通过数值优化与验证模拟得出，最优点位于探索范围边界（70 μm、25 kHz、0.07 m·s^?1），预测减阻率69.53%，验证值69.25%，相对误差0.40%。

3.4 激励方向：横向与纵向
通过在最优工况（70 μm、25 kHz、0.07 m·s^?1）下仅改变激励方向的对照模拟得出：横向振动减阻69.25%，纵向振动仅16.88%，横向模式降低的切削阻力为纵向模式的4.1倍，说明接触荷载重分配比缩短净切削行程是更高效的减阻路径。

讨论部分，研究人员对结果给出物理诠释：颗粒集合对周期激励的响应由激励位移振幅与特征颗粒间距之比控制，振幅通过逐步破坏与重建颗粒间接触将加载模式由压缩主导转为剪切主导；频率的作用体现为控制接触循环速率，在25–45 kHz内时间平均接触力分布对其无线性依赖；切削速度通过工具与颗粒网相互作用时间调节，高速时冲击主导的接触状态放大位移幅值对接触力重分配的作用。与纵向振动研究的比较表明，横向模式在工具尺度保持连续接触、在颗粒尺度产生循环加卸载，其四倍减阻优势仅在土壤状态、刀具几何与减阻率定义完全匹配的条件下成立。研究局限包括：均匀键合球粒与单一含水率（9%）简化了土体；刚性工具边界运动、波传播与热效应由标定接触系数有效承担；台架试验频率（20 kHz）与模拟范围（25–45 kHz）存在错位，未来需拓展全频率范围台架试验、田间速度、纵横复合模式与土壤粉碎质量响应。

研究结论部分翻译如下：本研究表明，土壤耕作是能耗密集型作业，土壤耕作部件切削阻力是最大单项损耗，超声振动辅助切削可降低该阻力，但定量研究多集中于纵向振动；横向振动通过循环侧向挤压与释放加载土壤，其减阻行为尚未充分探索。在30–70 μm、25–45 kHz与0.03–0.07 m·s^?1范围内，振幅是主导因素，减阻率由44.10%单调升至68.18%，频率线性主效应不显著（p=0.0571），速度线性效应不显著（p=0.1773）；横向振动减阻的本质在于接触力学状态的转化——荷载由法向压缩主导转向切向剪切主导，接触网松动、键破坏模式由集中破坏转向渐进降解；回归模型拟合与预测精度高（R²=0.9968，R²_pred=0.9489），最优参数组合（70 μm、25 kHz、0.07 m·s^?1）位于探索范围边界，预测减阻率69.53%，离散元验证69.25%（相对误差0.40%）；20 kHz台架试验证实标定参数在物理切削条件下适用，频率范围内结论依赖已验证模型。研究同时指出，总能耗与田间性能尚未评估，需在未来工作中完善。