[论文解读] Nitsche methods for constrained problems in mechanics

Tom Gustafsson, Antti Hannukainen+2 more|arXiv (Cornell University)|Mar 5, 2026

Advanced Numerical Methods in Computational Mathematics被引用 0

一句话总结

该论文提出了一种通用的、基于最小化的尼茨赫框架，用于在力学中对等式和不等式约束进行强制 enforcing，并在各种接触和不等式问题上给出数值收敛证据。

ABSTRACT

We present guidelines for deriving new Nitsche Finite Element Methods to enforce equality and inequality constraints that act on the value of the unknown mechanical quality. We first formulate the problem as a stabilized finite element method for the saddle point formulation where a Lagrange multiplier enforces the underlying constraint. The Nitsche method is then presented in a general minimization form, suitable for nonlinear finite element methods and allowing straightforward computational implementation with automatic differentation. To validate these ideas, we present Nitsche formulations for a range of problems in solid mechanics and give numerical evidence of the convergence rates of the Nitsche method.

研究动机与目标

• 提供将尼茨赫有限元方法导出以作用于固体力学未知量的等式与不等式约束的指南。
• 将尼茨赫方法重新表述为稳定化的拉格朗日乘子问题并作为适用于非线性有限元实现的能量最小化问题。
• 展示针对具体问题（如膜/板接触与不等式边界约束）的新颖尼茨赫形式。
• 提供数值收敛速率证据并讨论实现细节包括自动微分。

提出的方法

• 将尼茨赫方法重新表述为 constrained minimization 的稳定化拉格朗日乘子方法。
• 制定包含约束 β(u_h) ≥ 0 及相应乘子 λ(u_h) 的一般尼茨赫最小化泛函。
• 通过强形式定义连续乘子 λ，并给出与网格和材料参数相关的稳定化 γ(h,κ)。
• 推导包含一个正分布项 (β/γ)+ 和一个二次稳定化项的通用尼茨赫能量 J_h(u_h)。
• 将通用形式具体化为等式与不等式约束，展示如何在元素级消去拉格朗日乘子，从而获得纯主变量的公式。
• 给出代表性问题（如膜障碍、两体接触、板接触、Kirchhoff板等）的显式公式，并借助自动微分的牛顿法求解。

Figure 1 : Numerical solution for two membranes in contact using the Nitsche method.

实验结果

研究问题

• RQ1 如何将尼茨赫方法推广以强制执行固体力学中形式为 β(u_h) ≥ 0 的任意不等式约束？
• RQ2 在选择稳定化参数 γ 以确保受限问题的稳定性与最优收敛性方面有哪些指南？
• RQ3 如何在离散层面消去拉格朗日乘子以得到可实现的非线性问题的原变量公式？
• RQ4 对代表性接触与不等式问题，所提出的尼茨赫公式的收敛行为如何？
• RQ5 与惩罚法相比，所提出的方法在稳定性、条件性与收敛方面有何差异？

主要发现

• 当正确考虑拉格朗日乘子残量时，通用尼茨赫框架可对等式与不等式约束给出稳定的公式。
• 当残量尺度与连续 Sobolev 范数相匹配且稳定化 γ 结合网格和材料参数时，可达到最优收敛速率。
• 拉格朗日乘子可以显式投影到非负约束集合，从而得到适用于如 Signorini-type 问题的实用尼茨赫实现。
• 数值演示显示：对两膜接触使用线性单元时在 H1 上的线性收敛，以及对 Kirchhoff 板的不等式条件时的二次收敛，验证所提出的方法。
• 该方法通过一个可用牛顿法和自动微分求解的一般最小化问题实现，并提供公开源码包以便复现。
• 框架通过多类问题类来示例：膜障碍、两体接触、膜-固体接触、板接触，以及带不等式边界条件的 Kirchhoff 板。

Figure 2 : Two membrane contact problem convergence rate follows the theoretical linear convergence in the $H^{1}$ norm for linear finite element basis.