NAFEMS 接触基准 — V&V结果总结报告

严格来说，"接触"是边界条件的一种。由两个物体互不穿透的"非穿透条件"和接触面压力始终为压缩（不产生拉力）的"压力条件"定义，还有接触面上是否存在摩擦力。NAFEMS基准的重点是这些条件如何在数值上得到满足。

在简单例子中，用接触点间隙

完全同意。有摩擦时，需要定义接触面上剪力与滑动速度、位移的关系。NAFEMS基准问题中常采用经典库仑摩擦模型

根本区别在于非穿透条件的处理方式。惩罚法用小虚拟弹簧防止穿透，穿透量与反力成正比。计算相对稳定，但原则上允许微小穿透。而拉格朗日乘数法引入附加自由度（乘数）严格满足条件，但方程规模变大，接触状态判定困难。看NAFEMS结果，在许多情况下，使用扩展拉格朗日法（两者结合）的高精度设置更接近参考解。

与线性弹性问题的简单收敛不同。接触点和接触边界依赖于网格离散化，粗网格的结果可能波动很大。NAFEMS的"LE10: 赫兹接触"问题报告称，要得到目标接触压力分布，接触区域至少需要8个以上的单元。误差不是线性递减，而是取决于算法和网格的匹配。

主要是为了速度，但间接影响精度。全搜索检查所有节点和面组合（

首先检查"接触刚度（惩罚刚度）"或"拉格朗日乘数允许穿透量"。Ansys中对应`FKN`，Abaqus中对应"接触刚度缩放因子"。参考NAFEMS结果，默认值（如Abaqus的1.0）在粗网格下会导致过大穿透，需要根据问题调整为0.1到10之间。但设置过高会导致数值病态（矩阵过硬）。

NAFEMS指南通常推荐"面对面"接触。这样接触压力分布计算更光滑。"节点对面"计算成本低，但从面粗糙时易发生"接触锁定"（过刚现象）。主面应设为更刚硬或网格较粗的那个。

非常有效，这是NAFEMS很多问题采用的基本策略。特别是初始接触（有间隙或干涉），先用小荷载稳定接触状态，再加本荷载。Abaqus/Standard有"步"功能，Ansys Mechanical有"分步加载"或"自动步"。一步加全荷载常导致接触状态剧变而不收敛。

默认设置下有差异，但适当调优后都能接近参考解。特点上，Abaqus/Standard的接触算法（默认扩展拉格朗日法）在业界广受信任，非常鲁棒。Ansys Mechanical默认惩罚法，但可精细控制穿透量。COMSOL的"接触对"设置直观，但大变形、大滑动问题设置难度相对较高。

"鲁棒"是指网格质量、初始间隙设置稍差，也能收敛并返回物理合理答案的倾向。Abaqus/Standard对接触状态判定和方程更新极其谨慎，自动调步能力强。但相比Ansys的显式动力学（LS-DYNA），计算时间通常较长。

CalculiX（Abaqus风格）和Code_Aster官方通过了许多NAFEMS基准，核心算法与高价软件相当。但前处理（网格生成）、后处理以及"收敛困难时的自动处理能力"和用户支持不如商业软件。即用户需更深入理解理论和求解器细节，初学者门槛较高。

首先怀疑"初始过干涉"。部件初始就陷入，求解器产生大反力而发散。对策1: 检查模型，设置间隙。对策2: 多数软件有"初始接触调整"功能（Abaqus的`*CONTACT INTERFERENCE`、Ansys的"接触工具"），初阶段解除过干涉。NAFEMS指南强调初始几何定义的重要性。

网格粗是可能，但另一主因是"接触条件离散化"。使用节点对面接触时最易出现。对策：1)改为面对面接触，2)细化网格（特别是接触区），3)后处理中对压力"节点平滑化"。NAFEMS LE11问题专门讨论如何抑制压力振动的网格要求。

摩擦引入强非线性。分步对策：第一，确认无摩擦时（

接触分析中，数%的不平衡多是设置错误信号。可能是接触点离荷载传递路径太远（接触未建立），或惩罚刚度过低导致大穿透产生"间隙"。首先检查接触状态可视化，看目标面是否接触。NAFEMS验证通常要求力平衡误差在0.5%内。