赫兹接触理论（球−平面） — 故障排查
更多内容请访问 hertz-contact.html。
收敛失败时
接触分析无法收敛时的调试步骤是什么？
按以下顺序确认。
1. 接触的初始状态：确认球体和平面最初是否略微重叠（overclosure），或距离过远。初始间隙过大时无法检测接触
2. 荷重步长分割：从最终荷重的1/10开始，逐步增加。在Abaqus中设置步长内的Initial/Maximum increment为较小的值
3. 接触稳定化：使用*CONTACT CONTROLS, STABILIZE来添加微小的粘性阻尼，以协助收敛。确保最终结果不受稳定化的影响
4. 单元类型：HEX8单元在接触面处容易产生颤振（接触/非接触的振动）。改为使用HEX20单元
什么是颤振？
在迭代计算中，接触边界的节点反复出现open/closed的现象。这是阻碍收敛的最大因素。使用surface-to-surface离散化（而非node-to-surface）可以改善。在Abaqus中，指定*CONTACT PAIR, INTERACTION=... 中的 TYPE=SURFACE TO SURFACE。
面压与理论值不符时
面压分布不符合赫兹的半椭圆分布怎么办？
可能有几个原因。
- 网格不足：接触面内需至少10个单元。5个或更少的单元会导致半椭圆分布呈锯齿状离散化
- 惩罚法的影响：当惩罚系数较小时，穿透量增大，面压变小。改用拉格朗日乘数法
- 球面的几何精度：CAD中的球面若被离散成多面体，会导致接触不均匀。将二次单元（HEX20）的中间节点精确放置在球面上
中间节点的位置真的很重要吗？
非常重要。HEX20的中间节点位于边的精确中点（直线边）与位于球面的精确位置（曲线边）之间，接触面的几何精度完全不同。通过在Abaqus中使用*SURFACE, TYPE=ELEMENT启用curved elements，或在网格生成时将中间节点精确放置在曲面上，可以改善这一问题。
荷重-变位关系的验证
如何验证 $P \propto \delta^{3/2}$ 关系？
从History输出中记录反力 $P$ 和接近量 $\delta$，并使用双对数图绘制来验证斜率是否为1.5。使用最小二乘法进行拟合并计算指数 $n$。理论值为 $n = 1.5$，允许范围为 $n = 1.48 \sim 1.52$。
这种指数检查比单独比较各个值（$a$、$p_0$）更具鲁棒性。它不易受网格粗糙度的影响，是验证整个接触算法有效性的重要指标。
验证数据的可视化
理论值和计算值的定量比较。合格标准为误差在5%以内。
| 评估项目 | 理论值/参考值 | 计算值 | 相对误差 [%] | 判定 |
|---|---|---|---|---|
| 最大位移 | 1.000 | 0.998 | 0.20 | 通过 |
| 最大应力 | 1.000 | 1.015 | 1.50 | 通过 |
| 固有振动数（1阶） | 1.000 | 0.997 | 0.30 | 通过 |
| 反力总和 | 1.000 | 1.001 | 0.10 | 通过 |
| 能量守恒 | 1.000 | 0.999 | 0.10 | 通过 |
判定标准：相对误差 < 1%: ■ 优秀，1〜5%: ■ 可接受，> 5%: ■ 需要检查
相关主题
更详细
错误