Hertz-Mindlin 颗粒接触模型：原理、公式与农业工程应用

Hertz-Mindlin 颗粒接触模型：原理、公式与农业工程应用

 

核心速览：Hertz-Mindlin 是离散元法（DEM）模拟颗粒 - 颗粒 / 颗粒 - 壁面接触的经典非线性弹簧 - 阻尼模型，核心是结合Hertz 法向弹性理论与Mindlin 切向粘滑理论，精准描述散体物料（如土壤、秸秆、种子、肥料）的接触力学行为，是农业智能装备（捡拾器、排种器、旋耕机等）仿真的核心工具。

 

---

 

一、核心理论基础

 

模型由两部分耦合而成，分别描述接触的 ** 法向（垂直接触面）与切向（平行接触面）** 行为，适配 DEM 的显式动力学计算逻辑。

 

表格

 

理论模块	核心来源	适用场景	核心假设
Hertz 法向理论	赫兹（1882）弹性接触理论	颗粒碰撞、法向挤压变形	光滑弹性球体、小变形、无摩擦
Mindlin 切向理论	Mindlin-Deresiewicz（1949/1953）	颗粒滑移、剪切、摩擦行为	考虑切向弹性位移、粘滑（stick-slip）机制

 

---

 

二、接触力计算核心公式

 

模型将接触力分解为弹性力（弹簧项）与阻尼力（阻尼项），并通过库仑摩擦准则限制切向力上限，完整公式如下：

 

1. 基础等效参数（必算）

 

接触前需先计算两接触体的等效参数，统一描述材料与几何特性：

 

• 等效杨氏模量 E∗：
   
  E∗1​=E1​1−ν12​​+E2​1−ν22​​
   
  （E1​/E2​：材料杨氏模量；ν1​/ν2​：泊松比）
• 等效半径 R∗：
   
  R∗1​=R1​1​+R2​1​
   
  （R1​/R2​：颗粒 / 壁面曲率半径，壁面可视为 R→∞）
• 等效质量 m∗：
   
  m∗1​=m1​1​+m2​1​

 

2. 法向接触力 Fn​

 

包含法向弹性力与法向阻尼力，由法向重叠量 δn​（接触后两物体的法向穿透量）与法向相对速度 vn,rel​ 决定：

 

Fn​=弹性力 Fn,elastic​34​E∗R∗​δn3/2​​​+阻尼力 Fn,damping​265​​β2m∗E∗R∗​δn​​​​

 

 

• 弹性力：非线性项，与 δn3/2​ 成正比，描述法向挤压变形；
• 阻尼力：与恢复系数 β（由恢复系数 e 计算：β=lne/(lne)2+π2​）相关，耗散碰撞能量，模拟冲击衰减。

 

3. 切向接触力 Ft​

 

包含切向弹性力、切向阻尼力与库仑摩擦限制，核心是粘滑机制（颗粒接触时先粘滞滑移，达到摩擦极限后滑动）：

 

Ft​=弹性力 Ft,elastic​−kt​δt​​​−阻尼力 Ft,damping​265​​β2m∗G∗R∗​δn​​​​

 

同时满足库仑摩擦准则：

 

∣Ft​∣≤μ∣Fn​∣

 

 

• 切向刚度 kt​=8G∗R∗δn​​（G∗：等效剪切模量，1/G∗=(1−ν12​)/2G1​+(1−ν22​)/2G2​）；
• δt​：切向重叠量（切向相对位移累积量）；
• μ：摩擦系数（决定颗粒滑移的临界切向力）。

 

---

 

三、关键物理参数与意义

 

表格

 

参数	物理意义	农业仿真影响
E∗	材料法向刚度，决定法向变形难易	土壤越硬 / 种子越硬，E∗ 越大，变形越小
R∗	接触几何等效半径	颗粒粒径越大，R∗ 越大，接触力非线性越显著
μ	摩擦系数，控制切向滑移阈值	秸秆与弹齿间 μ 大，易粘附漏捡；土壤与旋耕刀 μ 大，耕作阻力大
e	恢复系数，控制能量耗散	种子碰撞排种盘 e 大，弹跳小，排种更均匀；土壤碰撞 e 小，破碎更充分
kt​/kn​	切向 / 法向刚度比	比值越大，切向抗滑移能力越强，颗粒流动性越差

 

---

 

四、农业工程核心应用（适配你的研究方向）

 

模型是农业智能装备离散元仿真的标配，核心应用场景如下：

 

1. 土壤 - 机具互作仿真

 

• 旋耕机、深松机：模拟土壤颗粒与旋耕刀 / 深松铲的接触，计算耕作阻力、土壤破碎效果，优化刀具结构（如旋耕刀曲率、铲刃角度）；
• 播种机开沟器：模拟土壤与开沟器的接触，优化开沟深度与土壤回填效果，减少漏播 / 重播。

 

2. 散体物料输送与捡拾

 

• 圆捆打捆机捡拾器：模拟秸秆颗粒与弹齿 / 轨道的接触，优化弹齿运动轨迹，降低漏捡率与机械冲击（对应你之前关注的捡拾器结构）；
• 气力 / 机械排种器：模拟种子颗粒与排种盘 / 窝眼的接触，优化排种参数（如转速、窝眼尺寸），提升排种均匀性与合格率。

 

3. 物料粘结与粘附模拟（模型变体）

 

• Hertz-Mindlin with Bonding：用于秸秆、稻株等粘结性物料，模拟秸秆粉碎、稻株破碎过程，计算断裂力与能量消耗；
• Hertz-Mindlin with JKR：用于湿润土壤、粘性肥料，考虑颗粒间粘附力，提升仿真精度。

 

---

 

五、模型局限性与注意事项

 

1. 核心假设限制：默认颗粒为光滑弹性球体，实际农业物料（如秸秆、不规则土壤颗粒）多为非球形、非光滑，需通过颗粒形状修正（如簇团模型）适配；
2. 小变形假设：大变形场景（如秸秆碾压、土壤压实）需结合弹塑性模型扩展；
3. 参数标定关键：农业物料参数（如土壤 E∗、秸秆摩擦系数）需通过 ** 室内试验（堆积角、剪切试验）** 标定，否则仿真结果误差极大。

 

---

 

六、与你研究的结合建议

 

1. 捡拾器优化：用 Hertz-Mindlin 模型模拟秸秆颗粒与弹齿的接触，通过调整弹齿轨道曲率（改变 δn​）与弹齿材料（改变 E∗），降低漏捡率与磨损；
2. 种子排种仿真：针对排种器，对比不同 μ（种子 - 排种盘材料）下的排种效果，优化窝眼设计；
3. 土壤耕作仿真：结合 JKR 变体模型，模拟湿润土壤与旋耕部件的互作，精准计算耕作阻力，提升装备能效。