拉格朗日乘数法 — CAE术语解释

最典型的应用是"接触问题"。例如，两个零件相互接触时，零件之间不能相互穿透，这就是"非穿透条件"。这个条件是约束条件，拉格朗日乘数法用于严格满足这个条件。乘数可以理解为接触面产生的"接触压力"。

从变分原理看比较容易理解。弹性体的势能Π(u)在位移u处取最小值，但接触面Γ_c上有约束条件g(u) ≤ 0。拉格朗日函数为

完全正确，这也是接触问题成为非线性问题的原因之一。求解器在迭代计算中判断每个接触点是"已接触"还是"分离"状态，然后寻找满足条件的λ和u。Abaqus/Standard使用的"拉格朗日乘数法"和Ansys Mechanical使用的"增广拉格朗日法"都是基于这个原理的不同实现。

关键的权衡在于约束的严格性和收敛性。纯拉格朗日乘数法可以在数值误差范围内严格满足非穿透条件，但整体刚性矩阵会出现零对角元素，方程求解会变得较困难。增广拉格朗日法使用这样的函数

这取决于模型的总体尺寸和单元大小。例如，10mm的零件、0.5mm的单元，允许穿透量通常被控制在1e-3 mm到1e-5 mm。在大多数实际的结构分析中，这个量级的误差对应力结果的影响可以忽略不计。Ansys的默认设置采用了这个增广拉格朗日法。

有两种情况。第一种是包含流体的密闭容器密封分析等，几何间隙绝对不能允许穿透。第二种是需要精确处理过度约束产生的刚体运动的情况。例如，轴承中球与滚道的接触，增广拉格朗日法会产生微小的滑移，影响旋转摩擦的计算。Abaqus/Standard的"拉格朗日乘数法（默认）"正是为这类情况设计的。

有几种对策。首先检查接触面的定义。单元表面应该光滑，主面和从属面的网格尺寸不应该相差太大（理想1:1，可接受1:10）。其次调整接触刚性（惩罚刚性）。Ansys中是"FKN"参数，Abaqus中是"接触刚性缩放因子"。试试从默认的1.0改为0.1或10.0，收敛往往会改善。

确实是权衡。经验规则是，穿透量不超过最小单元尺寸的1%，对应力分布的影响就很小。先改为0.1运行，检查软件输出的穿透监控结果（Abaqus的`.msg`文件、Ansys的求解器输出）。如果穿透在允许范围内且收敛，就可以了。如果穿透太大，就试试0.5等中间值。

最少要检查三项。1. **接触压力（CPRESS）的连续性和特异点**：压力分布是否物理合理（尖锐峰值可能是网格依赖）。2. **接触状态（CSTATUS）**：接触/分离区域是否符合预期。3. **穿透量/间隙量（COPEN）**：是否在允许范围内。特别是Abaqus，要检查"接触约束"的满足程度。这些都是拉格朗日乘数法是否正确运作的直接证据。

术语和默认设置差异很大。**Abaqus/Standard**的面-面接触默认使用"拉格朗日乘数法"，这接近纯拉格朗日乘数法，严格处理约束。而**Ansys Mechanical**的默认是"增广拉格朗日法"，在软件中有时也简称为"拉格朗日法"。如果想在Ansys中得到接近纯拉格朗日乘数法的效果，需要选择"纯惩罚法"并设置非常高的惩罚刚性，比较绕。

COMSOL的接触基础是"惩罚法"，但也有"Nitsche法"这个选项。Nitsche法不引入额外的乘数变量，而是在弱形式中直接施加约束，稳定性和精度的平衡很好。它特别适合多物理场问题，比如流体-结构耦合（FSI）中的移动界面"滑动网格"条件。不过COMSOL的接触算法调优选项没有Ansys或Abaqus那么丰富。

"更好"应该改成"更合适"。**Abaqus/Standard**在汽车螺栓紧固、橡胶密封等工业界要求严格接触分析的领域积累丰富，算法稳健。**Ansys Mechanical**的增广拉格朗日法收敛性好，适合复杂大变形接触（如超弹性材料的滑动）。**COMSOL**灵活定义热接触阻力、电接触阻力等特殊物理场耦合接触，更适合研发。

反而是拉格朗日乘数法的"优势"引起的现象。这个方法力求严格满足约束，如果一个节点同时被多个主面约束，就会产生约束重复（过度约束），导致方程组数值奇异。惩罚法因为有某种"柔性"，反而不易出现这个错误。

首先定位出错的节点。常见的两个原因：1. **凹角处的自接触**：在锐角凹陷处，一个从属节点可能同时被相邻的两个主面段判定为接触。对策是分割接触面定义，或在几何上添加微小倒角。2. **对称面与接触面的交线**：对称边界条件和接触条件在同一节点上约束。对策是从接触面定义中排除靠近对称面的节点（Abaqus有`Adjust`选项）。

完全正确，物理上不合理。这是数值计算不稳定的信号，通常是迭代过程中的临时现象。主要原因是接触刚性太低或荷载增量太大，导致接触状态急剧变化。对策是更细致地分割荷载步（使用Ansys的`Auto Time Stepping`、Abaqus的`Stabilization`），或将接触刚性（FKN等）上调一个档级再计算。如果还是不行，检查模型本身（例如初始穿透）是否有问题。

好的观察。乘数作为接触面节点的附加自由度，如果接触面节点很多（例如10万节点），会产生不可忽略的内存和计算成本。对策有两个。第一，定义接触面时只选择真正可能接触的区域。第二，重新考虑求解器设置。如果Abaqus/Standard使用的是直接法求解器，疏矩阵直接法（例如`PARALLEL=ALL`）对乘数引入的零对角元素有较好的容错性。如果使用迭代法求解器（`ITERATIVE`），需要仔细调优预处理选项。