接触解析的收敛失败

"接触"这种状态的变化本身是不连续的。接触力由间隙函数

拉格朗日乘数法确实禁止贯穿约束，但接触状态（活跃/非活跃）的判定在每个增量步骤中变化的"状态非线性"仍然存在。例如，Abaqus/Standard使用增强拉格朗日法时，需要重复进行接触压力修正的内迭代，这可能发散。特别是在接触表面出现"角"或"边缘"处面压分布急剧变化的情况下，一个步骤中接触集合可能大幅改变，这是原因所在。

摩擦引入了"粘滞"和"滑动"这样的另一种状态变化。在库伦摩擦模型中，摩擦力

许多商业软件采用"预测子-修正子"策略。在预测步骤中固定接触条件求解变位，在修正步骤中检查接触状态（贯穿或拉伸），更新刚性矩阵或残差力。COMSOL的"连续接触"功能是将接触边界近似为光滑函数使其可微，以稳定化这个迭代的尝试。但如果初始贫穿很大，预测步骤时已经出现物理上无法进行的变形，修正就追不上。

影响很大。简单的"节点对表面"容易发生"口袋效应"贯穿，残差力会振动。而"表面对表面"算法（Abaqus和MSC Marc的标准）能更光滑地处理接触压力或力的分布，收敛性更好。但计算成本会提高。Ansys的"基于MPC"或"纯惩罚"也是算法选择，需要根据问题切换。

主要有两个标准。一是Newton迭代的收敛速率。例如迭代5次内残差力范数无法收敛到规定值（默认0.5%）。二是接触状态的急剧变化量。Abaqus/Standard中，如果一个增量步骤内接触状态变化的单元比例超过一定值（默认约30%），就会将步长减半重新计算。重复这个过程会达到最小增量大小（例1e-8），导致解析停止。

首先是"排除初始贯穿"。检查CAD装配的间隙（Clearance）或干涉（Interference），使其与设计相符。可以用Ansys SpaceCllaim或CATIA的"间隙分析"功能。其次是"接触面网格整合性"。接触面之间的网格大小或位置差异过大，力的传递会不稳定。经验法则是，从属侧网格应比主侧网格细约1.5倍。

"调整"是数值稳定化的便利处理。Abaqus的"Contact Interference Fit"或Ansys的"Offset"功能，将微小的初始贯穿或间隙"调整"为零，使求解器从一开始就能正确判定接触状态。但是如果硬生生调整较大的干涉（例如0.1mm以上），会产生意外的初应力，要注意。它只是为了消除网格微小不一致（约网格大小的1e-3倍）的工具。

"阶段施加荷载"是基本原则。例如螺栓接合部的分析，先只施加螺栓轴力建立接触，再加上横向荷载。此外，位移控制往往比力控制收敛性好。原因是求解器能直接控制位移，可以防止接触面"逃逸"的现象。Ansys的"稳定化（Damping）"功能也是在最初几个步骤中加入人工阻尼抑制振动，实务中常用。

经验上从接触面周围最硬部材的杨氏模量和网格尺寸估算。一个参考是

/Standard是隐式方法，如前所述通过迭代严格满足接触条件。/Explicit是显式方法，每个时间步根据当前几何信息直接计算接触力。无迭代过程，对复杂接触状态变化（如多零件相互碰撞）的收敛性优势明显。但稳定时间步会很小，准静解析需要质量缩放等技巧。汽车碰撞分析用Explicit，螺栓接合分析用Standard比较合适。

"Program Controlled"由软件内部判断，多数情况下默认是"Pure Penalty"。这计算快但为了将贯穿量保持在允许值内需要高接触刚性，成为振动原因。"Augmented Lagrange"在惩罚项上加拉格朗日乘数修正，指定允许贯穿量（默认为接触工具厚度的0.1%）同时能保低接触刚性。收敛性改善，特别是非线性材料或大变形接触问题中，许多实务工作者从一开始就选择Augmented Lagrange。

"接触对"是离散接触算法，类似Abaqus。"连续接触"将间隙函数近似为光滑函数（如逻辑函数），使接触条件连续化和可微化。这使Newton法迭代非常稳定，收敛半径扩大。代价是严格的"接触/非接触"边界模糊，总会有微小贯穿或拉伸。对微细机构的初期行为研究连续接触有效，但对接触压力分布的精密评估需要接触对。

基本算法（惩罚法、表面对表面）已实现。CalculiX（Abaqus风格输入）或Code_Aster（Salome-Meca环境）也能进行简单的面压接触分析。但收敛性高级功能（Ansys的"稳定化"、Abaqus的"自动接触稳定化"）或高级搜索算法、错误发生时的详细诊断消息不足。因此难题中需要用户自己细调网格或荷载条件，经验很重要。商业软件则将这些"知识"黑箱化、自动化了，占比很大。

这是Newton迭代未在默认上限（通常16次）内收敛的错误。首先检查.msg或.dat文件，看残差在哪次迭代中振动，接触状态是否频繁切换。具体对策有：(1)进一步减小初始增量大小（如从1e-3到1e-5），(2)将接触刚性缩放因子（惩罚刚性）从默认1.0降至0.1，(3)将接触定义从"小滑动"改为"有限滑动"（大变形情况）。

是数值上超过了允许贯穿量。Augmented Lagrange法设置"许可贯穿量"，接触刚性太低或荷载变化太急剧，导致计算上产生了超过这个许可量的贯穿。对策是提高接触刚性（Normal Stiffness）或用"Interface Treatment"中"Adjust to Touch"调整初始状态。另外网格太粗可能无法表现接触面曲率，导致局部深贯穿，细化接触区域网格也有效。

这叫"接触振动"或"接触噪音"，一种数值不稳定现象。每个增量步接触状态（活跃/非活跃）在单元单位上频繁变化。表面上收敛了但物理上不自然。原因多是接触面网格非整合且接触刚性太高。对策有：(1)整合网格，(2)使用"表面对表面"接触，(3)稍降接触刚性，(4)Abaqus中应用"光滑化"选项。这现象的接触压力积分值（总力）不差，但局部应力不可信。

摩擦系数的绝对值不是问题，而是"粘滞/滑动"判定的不稳定。对策是体系的。第一，最初几步关闭摩擦建立接触状态，然后用*MODEL CHANGE（Abaqus）或分离分析设置（Ansys）启用摩擦。第二，降低摩擦刚性（Tangential Stiffness）（如默认的0.01倍）。第三，改摩擦模型为"弹性滑动"或微滑模型，缓和严格库伦摩擦的不连续性。都试了还不行，可能那个接触面本质上"粘滞"状态物理上不能成立（一直滑动）。