An Investigation of Stabilization Scaling in Finite-Strain Virtual Element Methods for Hyperelasticity

本文提出了一种无需子网格、仅作用于核空间的有限应变虚拟元稳定化方法，通过解耦偏量与体积通道并采用与牛顿切线能匹配的缩放策略，有效解决了传统方法在近乎不可压缩超弹性问题中因引入体变代理项而导致的剪切模态过度刚化及网格敏感性问题。

要点

这是一篇关于计算机模拟材料变形（特别是橡胶、生物组织等软材料）的学术论文。为了让你轻松理解，我们可以把这篇论文的研究内容想象成是在解决一个"如何给形状奇怪的积木搭建更稳固的脚手架"的问题。

1. 背景：什么是“虚拟单元法”？

想象你正在用电脑模拟一块橡皮泥被拉伸、扭曲的过程。为了计算，我们需要把这块橡皮泥切成很多小块（网格）。

• 传统方法：通常把橡皮泥切成标准的三角形或正方形。
• 虚拟单元法 (VEM)：这是一种更高级的方法，它允许你把橡皮泥切成任意形状的多边形（比如五边形、六边形，甚至像不规则的石头）。这就像是用乐高积木搭房子，但允许你使用各种奇形怪状的积木块，这让模拟更灵活、更适应复杂的形状。

2. 核心问题：当积木块“太软”时会发生什么？

这篇论文关注的是超弹性材料（比如橡胶），它们在受力时会发生巨大的变形，而且几乎不可压缩（体积不变，像水一样，你压它，它只会往旁边鼓出来，不会变小）。

在计算机模拟中，为了计算方便，我们只计算每个积木块“平均”的变形（就像只看积木块的中心点）。但是，积木块内部其实还有很多看不见的、细微的变形模式（比如像手风琴一样折叠，或者像沙漏一样扭曲）。

• 问题所在：传统的计算方法只关注“平均变形”，忽略了这些“看不见的模式”。如果这些模式没有受到控制，电脑模拟就会出错，要么积木块乱飞（不稳定），要么变得像石头一样硬（锁死），导致算不出真实的柔软变形。
• 现有的“补丁”：以前的工程师为了修补这个漏洞，会在每个积木块内部再强行切出很多小三角形（子网格），然后给这些看不见的模式加上一层“胶水”（稳定项）来固定它们。
  • 缺点：这层“胶水”有时候太粘了。特别是当材料变得非常像水（不可压缩）时，这层胶水会错误地把“体积变化”的力加到了“剪切变形”上。结果就是：明明材料很软，电脑却算出它硬得像铁，而且这种错误会随着网格切得越细越明显（这就是所谓的“体积锁死”）。

3. 这篇论文的解决方案：一种“聪明”的新胶水

作者提出了一种全新的、更聪明的“胶水”（稳定化方法），它不需要在内部切小三角形，而且能分清“剪切”和“体积”两种力。

我们可以用三个生动的比喻来理解它的创新点：

比喻一：不再依赖“内部地图”（无子网格）

• 旧方法：就像你要检查一个房间是否稳固，必须进屋把房间分成很多小格子，在每个小格子里放传感器。这很麻烦，而且如果你把房间形状改一下，传感器就得重新摆。
• 新方法：作者说，不需要进屋。只要站在门口，看看房间四个角的顶点（自由度）有没有偏离它“应该”在的位置，就能知道房间内部是否稳固。
  • 好处：不需要在内部切小三角形，计算更简单，也不受内部怎么切的影响。

比喻二：分清“推挤”和“扭曲”（解耦）

• 旧方法：就像你给一个弹簧上油，不管弹簧是被拉长（体积变化）还是被扭弯（剪切变形），你都倒同样多的油。当弹簧很难被拉长（不可压缩）时，油倒多了，弹簧反而被粘住扭不动了。
• 新方法：作者设计了一个双通道阀门。
  • 通道 A（剪切）：专门负责处理“扭曲”。无论材料多硬，这个通道只根据材料的“剪切刚度”（像橡胶的韧性）来调节力度。
  • 通道 B（体积）：专门负责处理“推挤”。当材料变得几乎不可压缩时，这个通道可以自动关闭或限制力度。
  • 效果：这样就不会因为材料“太硬”（不可压缩）而把“扭曲”的通道也堵死了。

比喻三：给不同形状的积木“量身定制”（方向感知）

• 旧方法：给所有形状的积木用同样强度的胶水。如果积木被拉得很长（像一根面条），用同样的胶水可能会导致它要么太软（塌了），要么太硬（断了）。
• 新方法：作者让胶水变得“眼力好”。它能看出积木是圆的还是长的。
  • 如果积木是长条形，胶水就会在长方向上给多一点支撑，在短方向上给少一点，就像给长面条加了一根筷子支撑，既稳固又不会把它压扁。
  • 这种方向感让模拟在形状很怪（比如被压扁的六边形）的网格上也能非常精准。

4. 实验结果：真的好用吗？

作者做了一个经典的测试：库克膜（Cook's Membrane）。

• 场景：想象一块梯形的橡胶板，一端被死死夹住，另一端被用力向下拉。这会产生巨大的扭曲和剪切。
• 挑战：当橡胶变得非常像水（不可压缩）时，旧方法算出来的橡胶板几乎拉不动（太硬了），而且网格越细，结果越奇怪。
• 新方法的胜利：使用作者的新“胶水”，无论网格是整齐的正方形，还是乱七八糟的六边形、多边形，算出来的橡胶板都能平滑地、真实地被拉弯，而且随着网格变细，结果越来越稳定，不再出现“假硬”的现象。

总结

这篇论文的核心思想是：在模拟软材料的大变形时，不要盲目地给所有变形模式加“胶水”。

作者发明了一种不需要内部切分、能分清“剪切”和“体积”、且能根据形状自动调整力度的新方法。

• 以前：像用大锤砸钉子，不管钉子多细，都用同样的力气，容易把钉子砸弯（锁死）。
• 现在：像用精密的螺丝刀，根据钉子的形状和材质，精准地施加合适的力。

这使得计算机模拟在处理橡胶、生物组织、软体机器人等复杂材料时，更加精准、快速且不容易出错。

技术摘要

这篇论文题为《有限应变虚拟元方法中双曲弹性稳定化缩放的调查》（AN INVESTIGATION OF STABILIZATION SCALING IN FINITE-STRAIN VIRTUAL ELEMENT METHODS FOR HYPERELASTICITY），由保罗·阿基拉·F·埃纳贝（Paulo Akira F. Enabe）和罗德里戈·普罗瓦西（Rodrigo Provasi）撰写。文章深入探讨了在有限应变超弹性问题中，虚拟元方法（VEM）稳定化项的设计及其缩放行为，特别是针对近不可压缩（nearly incompressible） regime 下的鲁棒性问题。

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

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

• 虚拟元方法（VEM）的机制：VEM 是一种在任意多边形/多面体网格上求解偏微分方程的伽辽金框架。在低阶（$k=1$）设置中，离散场由边界自由度唯一确定，并通过多项式投影算子（$\Pi^\nabla_E$）计算一致性贡献（consistency contribution）。
• 稳定化的必要性：由于一致性贡献仅依赖于投影分量，投影核（$\ker(\Pi^\nabla_E)$）中的未解析模式（missing modes）在一致性能量中是“不可见”的。因此，必须引入**稳定化项（stabilization）**来控制这些模式，防止零能模式（hourglass modes）并保证离散问题的适定性。
• 现有方法的缺陷：
  1. 依赖子网格（Submesh-dependence）：现有的非线性（超弹性）VEM 稳定化通常通过在辅助子三角剖分（sub-triangulation）上积分非线性代理能量（surrogate energy）来定义。这使得稳定化项依赖于任意内部剖分的几何形状，而非单元本身。
  2. 体积代理污染（Volumetric Proxy Contamination）：在接近不可压缩极限（泊松比 $\nu \to 1/2$）时，现有的稳定化方案常通过修正的拉梅参数（modified Lamé parameters）进行缩放。这导致体积项（bulk term）的代理量被错误地注入到剪切型（shear-type）稳定化惩罚中。
  3. 后果：这种“体积代理泄漏”会导致对等体积（isochoric）未解析模式的过度惩罚，人为地增加刚度，从而在有限应变模拟中引发**体积锁定（volumetric locking）**现象，特别是在近不可压缩材料和扭曲网格上。

2. 方法论 (Methodology)

作者提出了一种**无子网格（submesh-free）、仅针对核（kernel-only）且偏量/体积解耦（deviatoric/volumetric decoupled）**的新型稳定化策略。

• 谱分析框架（Spectral Framework）：

  • 文章首先建立了稳定化的谱分析框架，将 VEM 的稳定性要求重新表述为稳定化双线性形式与目标牛顿切线能量（Newton tangent energy）在投影核上的谱等价性（spectral equivalence）。
  • 利用广义瑞利商（generalized Rayleigh quotients）和广义特征值，提供了一种与基无关的诊断工具，用于评估稳定化项如何为未解析模式分配能量。

• 新型稳定化构造：

  • 基于残差：稳定化项直接基于投影残差 $r(u_h) = u_h - \Pi^\nabla_E u_h$ 在自由度（顶点）处的值构建，无需内部子三角剖分。
  • 偏量/体积解耦：
    • 偏量项（Deviatoric term）：仅由剪切模量 $\mu_E$ 缩放。为了适应各向异性多边形，引入了基于几何主方向（principal directions）的有界各向异性权重（anisotropy weights），使稳定化刚度能沿特定方向重新分配。
    • 体积项（Volumetric term）：由独立的体积模量 $\kappa_E$ 缩放，定义为边界法向残差的惩罚。在 $\nu \to 1/2$ 时，$\kappa_E$ 可以被截断或设为零，从而完全消除核上的体积惩罚，避免残留压缩性。
  • 解析线性化：该构造具有闭式残差和切线，支持牛顿法的一致线性化（consistent linearization）。

• 理论证明：

  • 在标准多边形正则性假设下，证明了提出的偏量稳定化在 $\ker(\Pi^\nabla_E)$ 上与 $\mu_E |\cdot|^2_{1,E}$ 是一致等价的。
  • 等价常数独立于网格尺寸 $h_E$ 和泊松比 $\nu$，确保了近不可压缩极限下的鲁棒性。

3. 主要贡献 (Key Contributions)

1. 谱视角的引入：将有限应变 VEM 的稳定化要求形式化为核上的谱等价条件，提供了基于广义特征值的精确稳定性度量标准。
2. 新型解耦稳定化方案：提出了一种无子网格、仅依赖核残差的稳定化方法。该方法严格分离了偏量和体积通道，避免了体积代理对剪切惩罚的污染。
3. 各向异性感知缩放：引入了基于几何特征的方向性权重，使稳定化能在各向异性或扭曲多边形上受控地重新分配刚度，解决了标量参数在扭曲网格上缩放失效的问题。
4. 理论缩放结果：证明了偏量稳定化在 $\nu \to 1/2$ 时保持剪切缩放（shear-scaled），且常数不随 $\nu$ 恶化，从理论上消除了近不可压缩锁定。
5. 诊断与验证：
   • 设计了单单元等体积核模式诊断，直接揭示了经典方法中体积代理泄漏的机制。
   • 通过**Cook 膜（Cook's membrane）**基准测试，验证了新方法在近不可压缩大变形下的鲁棒性。

4. 数值实验结果 (Results)

• 单单元诊断（Single-element diagnostic）：

  • 测试了一个等体积的核模式。结果显示，经典稳定化方法在 $\nu \to 1/2$ 时，归一化后的稳定化能量显著增加（表明体积代理导致剪切模式过刚）。
  • 相比之下，提出的解耦稳定化能量保持恒定，仅随剪切模量 $\mu$ 变化，证明了其纯剪切缩放特性。

• 模态分析（Modal analysis）：

  • 在投影核上限制稳定化算子进行特征值分析。结果显示，提出的方法能精确分离偏量和体积模式：偏量特征值保持恒定，而体积特征值随 $\kappa_E$ 变化。
  • 对于各向异性网格，特征值比率与设计的主轴权重比率高度一致，证明了方向性刚度分配的有效性。

• Cook 膜基准测试（Cook's membrane benchmark）：

  • 在 $\nu = 0.499$（近不可压缩）条件下，对四种不同类型的网格（规则四边形、扭曲四边形、高度扭曲四边形、Voronoi 多边形）进行了大变形剪切加载测试。
  • 经典方法：表现出明显的体积锁定特征。粗网格下位移被严重低估（锁定），且随着网格细化收敛缓慢，不同网格族之间的结果差异巨大。
  • 提出的方法：在所有网格族和细化级别上均表现出一致的收敛行为。粗网格响应虽有波动，但迅速收敛到一个共同的细网格解（约 8.5），消除了锁定现象，证明了其在近不可压缩 regime 下的卓越鲁棒性。

5. 意义与结论 (Significance)

• 理论意义：文章澄清了 VEM 稳定化在有限应变下的本质——即作为未解析模式的能量分配机制。通过谱等价性分析，明确了稳定化必须独立于体积代理，并在近不可压缩极限下保持剪切缩放。
• 实践意义：提出的解耦稳定化方案无需内部子三角剖分，计算效率高，且能一致线性化。它有效解决了多边形网格在超弹性大变形模拟中常见的锁定问题，特别是针对近不可压缩材料和高度扭曲网格。
• 未来方向：该方法可推广至高阶 VEM、三维多面体以及更复杂的超弹性模型。此外，将谱诊断集成到自动化参数选择中也是未来的潜在研究方向。

总结：该论文通过理论分析和数值实验，有力地证明了传统的基于子网格和代理能量的稳定化方法在近不可压缩有限应变问题中存在根本性的缺陷（体积代理泄漏），并提出了一种基于核残差、偏量/体积解耦的新型稳定化策略，显著提高了 VEM 在复杂几何和材料非线性问题中的鲁棒性和精度。