4. 接触算法详解与选择策略

在 Ansys Mechanical 中，接触算法决定了法向接触力如何被施加、接触约束如何被满足，以及求解器在迭代中如何判断“接触已收敛”。同一个装配模型，换一种算法，可能从“反复不收敛”变成“十几步收敛”，也可能从“位移合理”变成“穿透明显”。理解四种主流算法——罚函数法、拉格朗日乘子法、增广拉格朗日法、MPC——以及梁连接的定位，是复杂装配仿真中最具工程价值的一项判断能力。

一、四种算法的本质差异

接触问题的核心矛盾是：两个面不能互相穿透，但求解器又必须用有限元方程表达这个“不能”。不同算法的区别，就在于它们如何处理这个约束。

罚函数法（Pure Penalty）把接触约束转化为“弹簧”处理。当从面节点侵入目标面时，求解器施加一个与侵入量成正比的力：F = k·x，其中 k 是法向接触刚度（Normal Stiffness）。侵入量越小，需要越大的 k；k 越大，约束越“硬”，但方程条件数越差。罚函数法不引入额外自由度，方程规模小、迭代快，但必然存在一定穿透，且穿透量取决于 k 与接触压力的比值。它适合接触状态明确、以法向压力传递为主、对微小穿透不敏感的场景，例如大多数螺栓预紧后的承压面、过盈配合的近似分析。

拉格朗日乘子法（Lagrange Multiplier）引入额外的乘子自由度，把接触力作为独立未知量直接求解，理论上可实现零穿透。代价是方程规模增大、可能出现负对角项，导致求解器需要更多迭代，且对初始接触状态敏感。它适合穿透必须严格控制的场合，例如密封面接触、薄壁件压合、需要精确提取接触反力的分析。

增广拉格朗日法（Augmented Lagrange）是前两者的折中：在罚函数基础上叠加拉格朗日乘子迭代，先用适中的刚度快速逼近，再通过乘子修正逐步消除穿透。Ansys Mechanical 中它常作为默认法向算法，兼顾收敛性与精度。对于大变形、接触状态频繁切换、摩擦参与的复杂装配，增广拉格朗日通常比纯罚函数更稳健，比纯拉格朗日更容易收敛。

MPC（多点约束）不通过“面—面”接触力，而是用约束方程把从面节点的位移与目标面节点位移绑定，例如把从面节点位移表示为目标面若干节点位移的加权组合。它不产生穿透，也不依赖接触刚度，但要求绑定关系在求解前确定，因此适合“始终接触、不允许分离”的连接，如绑定接触、远端位移耦合、刚性连接。若接触可能分离或滑移，MPC 会给出错误约束。

梁连接（Beam Connection）不是接触算法，而是一种以无质量梁单元传递力和力矩的连接方式，常用于模拟螺栓、销钉、点焊或结构件之间的柔性连接。它把两个面或两个点通过梁单元耦合，可传递拉压、剪切、弯曲和扭转，刚度由用户定义的梁截面与材料决定。它适合需要体现连接柔度、又不想建立完整螺栓实体和接触对的装配。

二、算法选择不是“越精确越好”

很多初学者会直接选拉格朗日乘子法，认为“零穿透”一定最好。实际工程中，这往往导致收敛困难。原因在于：拉格朗日乘子法把接触力作为独立未知量，当接触状态在迭代中反复切换（开—闭—滑移）时，乘子会剧烈振荡，残差曲线出现平台甚至发散。此时改用增广拉格朗日法，让罚刚度先“兜住”接触，再逐步修正，通常能在更少迭代内稳定下来。

反过来，如果分析目标是提取密封面的接触压力分布，穿透量必须控制在极低水平，纯罚函数即使把法向刚度调到 1e6 N/mm³ 仍可能残留可见穿透，这时拉格朗日乘子法或增广拉格朗日法配合较严格的收敛准则更合适。判断标准不是“哪个算法高级”，而是“这个装配的物理行为是否允许近似”。

三、复杂装配中的组合策略

真实装配很少只用一种算法。一个典型的减速器箱体装配可能同时包含：

• 箱体结合面：绑定接触，用 MPC 或绑定约束，避免无意义的接触迭代；
• 轴承外圈与箱体孔：过盈配合，用罚函数法配合摩擦，关注接触压力；
• 齿轮啮合面：增广拉格朗日法，因为接触状态随转动频繁切换；
• 螺栓连接：若只关心被连接件的夹紧效果，可用梁连接替代完整螺栓实体，显著降低规模；若需评估螺纹牙受力，则必须建实体并采用增广拉格朗日法。

这种“分区选算法”的思路，比全局统一设置更贴近工程实际。在 Mechanical 中，每个接触对都可以独立设置法向算法、切向算法、刚度更新方式和摩擦系数，正是为这种组合策略服务的。

四、关键参数与操作要点

选定算法后，真正影响收敛的是参数。以增广拉格朗日法为例，在接触对细节设置中需要关注：

法向刚度因子（Normal Stiffness Factor）：默认 1.0，可先降到 0.1 帮助收敛，再逐步恢复
刚度更新（Update Stiffness）：Each Iteration / Each Substeps / Never
穿透容差（Penetration Tolerance）：增广拉格朗日下用于判断乘子修正是否停止
摩擦系数（Friction Coefficient）：0.1～0.2 为常见钢—钢干摩擦示例值
摩擦算法：Pure Penalty 或 Augmented Lagrange，与法向算法独立

一个可复现的调试顺序是：先用增广拉格朗日 + 刚度因子 0.1 + 每子步更新刚度，确认模型能收敛；再逐步把刚度因子提高到 1.0，观察接触压力与穿透量是否稳定；若仍不收敛，检查是否接触对定义过宽、是否存在初始穿透、是否刚体位移未约束。对于梁连接，则需检查梁截面方向、端点耦合范围以及是否与周围接触对产生过约束。

五、选择策略速查