Four-field mixed finite elements for incompressible nonlinear elasticity

本文提出了一种基于位移、位移梯度、第一 Piola-Kirchhoff 应力和压力四场混合公式的有限元方法，用于求解不可压缩非线性弹性问题，该方法采用不连续位移场且无需稳定化，并建立了相应的适定性理论、先验误差估计及高效后处理技术，在二维和三维数值实验中均展现出最优收敛率与卓越的鲁棒性。

要点

这篇论文介绍了一种新的数学计算方法，用来模拟那些**“怎么压都不变形体积”（不可压缩）的“橡皮泥”或“生物组织”**（非线性弹性体）在受力时的复杂行为。

为了让你更容易理解，我们可以把这篇论文的核心内容想象成**“如何用最聪明的方法给橡皮泥建模”**。

1. 背景：为什么这很难？（橡皮泥的烦恼）

想象你手里有一块超硬的橡皮泥（比如橡胶或人体组织）。

• 特性：无论你怎么捏、怎么拉，它的体积永远不变（不可压缩）。
• 挑战：如果你用传统的计算机方法（有限元法）去模拟它，就像试图用一张完全不能拉伸的网去包裹一个正在变形的球。
  • 如果网太紧（为了模拟体积不变），计算机就会“卡死”，算不出结果，或者算出一些很假的形状（比如网格乱成一团，或者出现奇怪的棋盘格花纹）。
  • 以前的方法（像 CSFEM）就像是一个挑剔的厨师：在二维（平面）时做得很好，但一到了三维（立体），就需要加很多“调料”（稳定化参数）或者用特殊的、很难找的“厨具”（非标准单元），否则就会翻车。

2. 核心创新：四场混合法（四个角色的交响乐）

这篇论文提出了一种新的**“四场混合”方法。你可以把它想象成一个四人乐队**，每个人负责演奏不同的乐器，共同完成一首曲子（模拟物理过程）：

1. 位移 (Displacement)：橡皮泥**“移到了哪里”**。
2. 位移梯度 (Displacement Gradient)：橡皮泥**“被拉伸或扭曲的程度”**（局部变形）。
3. 应力 (Stress)：橡皮泥**“内部有多大的力在对抗变形”**。
4. 压力 (Pressure)：橡皮泥**“为了保持体积不变而产生的内部挤压力”**。

以前的做法：通常是先算出“位移”，然后强行推导其他三个量。这就像只让主唱唱歌，其他乐器手只能跟着猜，容易走调。
这篇论文的做法：让这四个角色同时独立演奏，互相配合。这样不仅更准确，而且不需要额外的“调料”（稳定化参数）来防止走调。

3. 最大的突破：让“位移”变得“不连续”

这是这篇论文最反直觉但也最聪明的地方。

• 传统思维：想象橡皮泥是连续的，像一块完整的布。在数学上，这意味着相邻的网格点必须严丝合缝地连在一起（连续位移）。
• 这篇论文的思维：作者说，“嘿，我们允许这块布在网格之间稍微断开一点点（不连续位移）”。
  • 比喻：想象你在拼乐高。传统方法要求每一块乐高必须完美咬合，不能有一点缝隙。但这篇论文允许乐高块之间有微小的缝隙，只要它们整体拼出来的形状是对的就行。
  • 好处：这种“允许断开”的策略，反而让计算变得极其简单和稳定。它不需要那些复杂的“特殊乐高块”（特殊单元），用标准的、随处可见的乐高块（标准有限元）就能搞定。
  • 结果：无论是在二维还是三维，它都不需要加任何“稳定剂”，就能算得稳稳当当。

4. 后期修补：把“断开”的乐高拼回去（后处理）

你可能会问：“既然允许断开，那算出来的橡皮泥看起来是不是破破烂烂的？”

• 答案：是的，直接算出来的结果，网格之间确实会有小缝隙。
• 解决方案：作者发明了一个**“快速修补术”**（后处理技术）。
  • 这就像是你先快速搭好一个大概的模型（允许有缝隙），然后花很少的时间，用一种简单的算法把缝隙抹平，让它看起来像一块完美的、连续的橡皮泥。
  • 这个修补过程非常快，而且能让结果变得超级精准，甚至比那些一开始就追求完美的传统方法还要好。

5. 实验结果：为什么它更牛？

作者做了很多测试（比如吹气球、拉橡皮膜、拉伸带孔的方块），把他们的新方法 (DDFEM) 和 旧方法 (CSFEM) 做对比：

• 三维大变形：旧方法在三维大变形时经常算崩（出现负体积、棋盘格噪音），而新方法稳如泰山。
• 精度：新方法算出来的应力分布更平滑，没有那些奇怪的“棋盘格”假象。
• 通用性：旧方法在二维和三维要用不同的“配方”，而新方法一套配方通吃，不管是 2D 还是 3D，都不用改。

总结：这篇论文说了什么？

简单来说，这篇论文发明了一种更聪明、更鲁棒、更通用的方法来模拟**“怎么捏都不变体积的软物体”**。

• 它的秘诀：允许计算过程中的“小断开”（不连续位移），通过四个变量同时求解，最后再花一点点时间“修补”成完美的连续体。
• 它的优势：不需要复杂的特殊设置，不需要额外的稳定参数，在三维大变形模拟中表现极佳，且能避免各种数值计算的“假象”。

这就好比以前我们要模拟橡皮泥变形，需要请一位挑剔的大师（旧方法），还得给他配特殊的工具；现在，我们只需要一位灵活的工匠（新方法），用普通的工具，稍微允许一点“不完美”的中间过程，最后却能做出更完美、更真实的作品。

技术摘要

这篇论文提出了一种用于不可压缩非线性弹性力学的稳定四场混合有限元方法。该方法基于位移、位移梯度、第一 Piola-Kirchhoff 应力和压力这四个独立场的混合公式，旨在解决传统方法在处理大变形不可压缩材料时遇到的体积自锁（volumetric locking）和数值不稳定问题。

以下是该论文的详细技术总结：

1. 研究背景与问题 (Problem)

• 背景：不可压缩非线性固体力学在生物组织力学、橡胶弹性和软材料设计等领域至关重要。由于解析解难以获得，有限元方法（FEM）是主要的数值工具。
• 挑战：
  • 传统的基于位移的单场方法在处理不可压缩材料时容易发生严重的体积自锁。
  • 现有的混合方法（如两场或三场公式）虽然能缓解自锁，但在大变形下可能面临稳定性下降的问题。
  • 现有的四场混合方法（如兼容应变混合有限元方法 CSFEM）虽然在 2D 和 3D 中表现稳健，但存在显著局限性：
    • 需要针对不同维度（2D/3D）使用不同的单元对。
    • 在 3D 情况下，需要非标准的有限元空间（如特殊的 Nédélec 元素）和额外的稳定化项（stabilisation），这增加了实现的复杂性和参数调节的难度。
    • 位移场通常是连续的，导致在粗网格上可能产生非物理的间隙或重叠。

2. 方法论 (Methodology)

作者提出了一种新的四场混合公式，其核心特征如下：

• 四场变量：
  1. 位移 ($u$)：在不连续的 Lagrange 空间中进行近似（这是与 CSFEM 的关键区别）。
  2. 位移梯度 ($K = \nabla u$)：在 $H(\text{div})$ 相容空间（如 BDM 元素）中近似。
  3. 第一 Piola-Kirchhoff 应力 ($P$)：在 $H(\text{div})$ 相容空间中近似。
  4. 压力 ($p$)：在连续 Lagrange 空间中近似。
• 弱形式推导：
  • 基于 Hu-Washizu 变分原理。
  • 通过对位移梯度方程进行分部积分，将位移边界条件弱施加，而牵引力边界条件强施加。
  • 推导了基于 Newton-Raphson 方法的线性化方案，用于求解非线性方程组。
• 稳定单元对：
  • 识别并证明了在 2D 和 3D 中均稳定的有限元对，无需针对不同维度更换单元类型。
  • 主要使用的稳定对包括：
    • 低阶：$P_0 d_1 d_1 P_1$（位移不连续，梯度/应力为 BDM1，压力为 $P_1$）。
    • 高阶：$P_1 d_2 d_2 P_2$。
  • 无需稳定化：该方法在 2D 和 3D 中均不需要额外的稳定化参数或项。
• 后处理技术：
  • 由于位移场在单元间不连续，为了获得物理上合理的连续位移场，作者提出了一种高效的后处理校正技术。
  • 通过求解一个定义在更高阶连续空间上的椭圆问题，利用计算出的位移梯度 $K_h$ 重构全局连续的位移场 $\tilde{u}_h$。这一步计算成本极低（只需解一个对称正定系统）。

3. 理论分析 (Theoretical Analysis)

• 适定性：证明了线性化连续问题的适定性（存在唯一解），基于双鞍点问题的抽象理论，验证了算子的有界性、强制性以及关键的 inf-sup 条件。
• 先验误差估计：推导了离散解的先验误差估计。
  • 证明了位移、位移梯度、应力和压力的收敛阶。
  • 证明了后处理后的位移场具有超收敛性（在 $L^2$ 范数下达到 $O(h^{k+2})$，在 $H^1$ 范数下达到 $O(h^{k+1})$）。

4. 数值实验结果 (Results)

作者在 2D 和 3D 中进行了广泛的数值实验，并与现有的 CSFEM 方法进行了对比：

• 测试算例：
  1. 径向膨胀（Radial Inflation）：2D 圆柱壳和 3D 空心球。
  2. Cook 膜（Cook's Membrane）：经典的弯曲和剪切组合问题，包含奇点。
  3. 带孔块体拉伸（Stretching of Perforated Blocks）：包含大变形和复杂几何。
• 主要发现：
  • 收敛性：提出的方法（DDFEM）在所有场变量上均达到了理论预期的最优收敛阶，甚至在某些情况下（如应力场的 $H(\text{div})$ 误差）表现出超收敛性。
  • 3D 稳健性：在 3D 问题中，该方法无需稳定化即可保持稳健，而 CSFEM 在 3D 中通常需要稳定化且对网格质量敏感。
  • 消除数值伪影：
    • CSFEM 在某些情况下（特别是低阶或粗网格）会出现**棋盘格（checkerboarding）**现象（应力和压力的非物理振荡）。
    • 提出的方法（特别是使用 BDM 元素）产生的应力和压力场非常平滑，无棋盘格伪影。
  • 体积保持：在带孔块体拉伸的大变形测试中，CSFEM 的某些单元对出现了非物理的负雅可比行列式（导致网格翻转），而提出的 DDFEM 方法能更好地保持体积，雅可比行列式始终接近 1。
  • 后处理效果：位移校正技术显著提高了位移场的精度，使其满足边界条件并消除了单元间的不连续，且计算代价很小。

5. 主要贡献与意义 (Contributions & Significance)

• 统一性：提出了一套适用于 2D 和 3D 的统一四场混合公式，无需针对不同维度设计不同的单元对。
• 无需稳定化：消除了对稳定化参数和复杂稳定化项的依赖，简化了实现过程，提高了方法的鲁棒性。
• 标准单元：仅使用标准的有限元空间（Lagrange 和 BDM 元素），避免了 CSFEM 所需的特殊 Nédélec 元素，使得在现有 FEM 框架中更容易实现。
• 高精度与鲁棒性：通过不连续位移场结合高效的后处理，既保留了混合方法的稳定性，又获得了高精度的连续位移解，有效避免了体积自锁和棋盘格现象。
• 应用前景：该方法为生物力学（如心脏电 - 机械耦合、肿瘤生长）和软材料的大变形模拟提供了一种可扩展、稳健的数值工具。

总结：这篇论文通过引入不连续位移场和特定的四场混合公式，成功克服了现有不可压缩非线性弹性混合有限元方法在 3D 实现中的复杂性和不稳定性问题，提供了一种无需稳定化、基于标准单元且高精度的数值解决方案。