🔢 mathematics

A robust mixed finite element formulation for third medium contact

本文介绍了一种用于第三类接触的鲁棒混合有限元公式，该公式利用辅助场稳定化结合惩罚耦合来正则化问题，并证明了其在涉及大变形和自接触的各种插值方案及基准场景中的有效性。

原作者： M. Vorwerk, J. Schröder, P. Wriggers

发布于 2026-06-29

📖 1 分钟阅读🧠 深度阅读

原作者： M. Vorwerk, J. Schröder, P. Wriggers

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下，你正试图在计算机上模拟两个柔软、有弹性的物体（比如橡胶块）如何相互碰撞。在工程模拟领域，这在技术上是极其困难的。

旧问题：“硬壁”法
传统上，计算机程序将接触处理为两个人试图穿过一堵实心墙。软件必须不断地询问：“它们接触了吗？它们重叠了吗？停！退回去！”它必须寻找碰撞的确切时刻并执行严格的规则。如果物体在进行复杂的滑动、扭转或折叠（比如橡胶圈回弹并缠绕自身），计算机就会陷入混乱，数学模型会崩溃，模拟也会中断。这就像是在试图驱赶一群正在迷宫中奔跑的猫。

新思路：“魔力海绵”
这篇论文介绍了一个巧妙的技巧，称为第三介质接触（Third Medium Contact）。它不再将物体之间的空间视为真空或硬壁，而是用一种“虚拟介质”填满这个间隙——想象一下，这就像是一个超软、隐形的魔力海绵。

接触前： 海绵非常柔软，物体可以穿过它而感觉不到任何阻力。这就像是在温暖的空气中移动。
接触期间： 随着物体互相挤压，海绵被压缩了。你挤压得越厉害，它就越难以被压缩。最终，它会变得像岩石一样坚硬，将物体推开。
结果： 在计算机眼中，物体从未真正“接触”；它们只是在压缩彼此之间的海绵。这避免了所有棘手的“搜索与停止”规则，使模拟过程更加平滑且不易崩溃。

新问题：“海绵坍塌”
这个海绵构想有一个隐患。如果物体过度挤压海绵，计算机的数学网格（mesh）会发生扭曲、变形，并最终坍塌，导致模拟失败。这就像试图挤压一块海绵，直到它变成一张扁平、无用的薄饼。

解决方案：“幽灵向导”
作者提出了一种新的稳定方法来防止海绵坍塌。他们引入了第二个不可见的场，称之为**“幽灵向导”**（技术上称为辅助场）。

以下是类比：
想象海绵是由微小的、独立的瓷砖组成的。如果你挤压海绵，某些瓷砖可能会被压扁，而它们的邻居却保持原样。这会产生一个参差不齐、不稳定的混乱局面。

幽灵向导就像是一层铺在这些瓷砖之上的、光滑且有弹性的薄片。
瓷砖与这张薄片松散地连接在一起。如果某块瓷砖相对于其邻居挤压得太过分，薄片就会将其拉回原位。
至关重要的是，这张薄片不需要是复杂的高科技材料。它可以是简单的、低分辨率的（使用“低阶”数学），这使得计算机计算既快速又廉价。

为什么这很重要
论文通过几个场景测试了这个想法：

两个方块： 将两个方块推在一起。新方法允许它们靠得非常近而不会导致数学崩溃，即使“海绵”是由非常简单的、块状的数学瓷砖组成的。
盒子： 一个在盒子内部弯曲并与其自身接触的框架。旧方法在这里经常失败，因为“海绵”变得过于扭曲。新的“幽灵向导”保持了海绵的稳定，使框架能够自然地折叠并相互接触。
C形结构： 一个弯曲的梁，其两端弯曲直到相互接触。新方法能够平滑地处理扭转和转向，其结果与更复杂、更旧的方法相匹配，但没有那些令人头疼的问题。

核心结论
作者发现，如果你让“幽灵向导”是连续的（在各部分之间平滑连接），它就像一个强大的稳定器，即使物体被挤压成极端形状，也能防止模拟崩溃。如果你让引导是“不连续的”（按每个瓷砖进行拆分），虽然计算起来更容易，但无法像前者那样有效阻止坍塌。

简而言之，这篇论文为工程师们提供了一种稳健、防崩溃的方法，通过用智能的“魔力海绵”和简单的“幽灵向导”取代硬性的接触规则，来模拟物体相互摩擦、滑动和折叠的过程，从而防止数学模型崩塌。

技术摘要：一种用于第三介质接触的鲁棒混合有限元公式化方法

问题陈述
经典的有限元接触公式依赖于对接触界面的显式离散化，通过惩罚法、拉格朗日乘子法或障碍法来施加约束。尽管这些方法已经成熟，但在涉及大滑动、自接触、严重网格畸变或接触边界先验未知（如拓扑优化）的情景中，面临着显著的算法负担。“第三介质接触”（Third Medium Contact, TMC）方法通过在潜在接触体之间的空间内填充一种高度顺应性的虚拟介质来解决这一问题，从而取代了显式约束。该介质在接触前保持柔软，但在受压时会剧烈硬化，通过连续体传递接触力。

然而，TMC 的一个核心挑战是稳定性。该介质必须足够柔软，以避免接触前的伪力；同时又必须足够稳定，以承受极端的压缩和剪切而不发生网格畸变。现有的稳定策略通常需要高阶有限元（以计算位移的二阶导数），或者引入增加系统规模的全局辅助场，亦或是依赖专门的离散化方法。因此，需要一种能够与简单、高效的一阶有限元兼容，且避免直接计算位移二阶导数的鲁棒稳定机制。

方法论
作者提出了一种引入辅助场稳定化的混合有限元公式化方法，用于处理第三介质。该方法的核心包括：

引入辅助场： 在第三介质中引入一个类似于变形梯度（deformation-gradient-like）的场，记作 $\Theta$ ，并将其视为一个独立的未知量，与物理位移场 $u$ 解耦。
弱耦合： 辅助场 $\Theta$ 通过惩罚项（ $W^p_{tm} \propto |\Theta - F|$ ）与物理变形梯度 $F$ 进行弱耦合。这确保了 $\Theta$ 能够追踪物理变形。
梯度正则化： 该公式化方法并非惩罚物理变形梯度的梯度（ $\nabla F$ ）——这需要计算位移的二阶导数——而是惩罚辅助场的梯度（ $\nabla \Theta$ ）。正则化能量定义为 $W^r_{tm} \propto |\nabla \Theta|^2$ 。
混合离散化： 该公式化允许对 $u$ $u$ 和 $\Theta$ $Θ$ 进行独立的插值。作者研究了 $\Theta$ $Θ$ 的连续与单元内不连续近似形式。
- 连续 $\Theta$ ： 节点值在单元间共享，从而产生单元间的耦合。
- 不连续 $\Theta$ ： 自由度局限于每个单元内部。这些自由度可以通过静态凝聚（static condensation）被消除，从而确保全局系统仅保留位移自由度。

主要贡献

低阶兼容性： 该公式化方法实现了梯度型正则化，而无需计算位移的二阶导数，使其能够兼容标准的低阶（线性）和高阶（二次）有限元。
稳定机制： 本文阐明了对于低阶位移单元（其中 $F$ 为单元内常数），使用连续插值的 $\Theta$ 是至关重要的。它将 $F$ 的单元间跳跃转化为 $\Theta$ 的非零梯度，从而激活正则化。相比之下，不连续的 $\Theta$ 场无法产生这些单元间梯度，从而削弱了低阶单元的稳定效果。
静态凝聚： 对于不连续辅助场，该方法证明了额外的未知量可以被局部凝聚，从而在不增加全局系统规模的前提下保留局部稳定性。

数值结果
该公式化方法通过了以下基准测试问题的评估：

双块压缩： 结果表明，第三介质缩放参数 $\gamma$ 控制着残余间隙（接触精度）与非线性收敛鲁棒性之间的权衡。较小的 $\gamma$ 会产生更紧密的接触，但需要更多的牛顿迭代次数。该方法在组合压缩和剪切条件下仍保持稳定。
箱体内的自接触：
- 插值对比： 连续辅助场（ $T_u^1 T_\Theta^1$ , $T_u^2 T_\Theta^2$ ）提供了鲁棒的接触闭合。混合对 $T_u^2 T_\Theta^1$ （二次位移，线性辅助）在提供与 $T_u^2 T_\Theta^2$ 相当的结果的同时，减少了辅助自由度。
- 连续与不连续对比： 不连续辅助场表现出正则化效果降低的问题，特别是对于线性单元，其梯度项会消失或无法实现单元间的耦合。
- 与变形梯度平均法的比较： 与 Faltus 等人 [18] 的方法相比，所提出的混合公式化方法所需的载荷步和牛顿迭代次数显著减少，展示了更优越的非线性鲁棒性。虽然混合公式化方法的代数系统规模更大，但其改进的收敛性补偿了成本；不过在极端精细网格（ $>10^6$ 个未知数）下，辅助自由度的成本会变得占据主导地位。
C 形箱体： 该公式化成功捕捉了大旋转、渐进式自接触以及严重的第三介质压缩。其力-位移响应与经典的节点-段接触公式高度吻合，并表现出平滑的接触过渡。
3D 扩展： 一个三维自接触示例证实了该方法在平面应变设置之外的适用性。

意义与主张
本文声称，所提出的辅助场稳定化提供了一种鲁棒且兼容低阶的替代方案，可用于传统的接触算法及现有的 TMC 正则化方法。通过避免直接评估位移的二阶导数，该方法使得使用简单、计算廉价的有限元（线性及二次）来处理涉及大变形和自接触的复杂接触问题成为可能。

作者强调，该方法的有效性本质上取决于辅助场的连续性；连续插值对于为低阶位移单元生成基于梯度的稳定机制是必不可少的。该公式化特别适用于接触边界未知或不断演化的应用场景，例如接触辅助机构、软体机器人以及拓扑优化，在这些领域，基于连续体的处理方式消除了对显式接触搜索算法和活跃集策略的需求。研究结论指出，虽然该方法是以更大的局部代数系统规模来换取更强的全局鲁棒性，但它为在不使用高阶单元或专门网格技术的情况下模拟剧烈接触场景提供了一条可行的路径。

类似论文