A higher order polytopal method for contact mechanics with Tresca friction

这是一份关于论文《A higher order polyhedral method for contact mechanics with Tresca friction》（具有 Tresca 摩擦的接触力学高阶多面体方法）的详细技术总结。

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

• 应用场景：该研究针对多孔介质中的接触力学问题，特别是涉及裂缝网络（Fracture Networks）的情况。这类问题在地下应用（如二氧化碳封存）中至关重要，因为注入气体引起的压力变化可能导致断层重新激活，引发岩石变形、泄漏甚至地震活动。
• 物理模型：
  • 几何结构：将裂缝建模为嵌入在三维基质中的共维数 1（codimension-1）的平面网络。
  • 力学行为：考虑线性各向同性弹性本构关系。
  • 接触条件：在基质 - 裂缝界面处应用 Tresca 摩擦模型。该模型规定了切向牵引力的屈服阈值（与法向牵引力无关），比库仑摩擦（Coulomb friction）在数学上更易处理，但能捕捉关键的摩擦行为。
  • 挑战：裂缝网络的几何复杂性（角点、交点）以及接触问题的非光滑性（不等式约束、非穿透条件、摩擦非线性）使得数值离散极具挑战性。传统的有限元方法（FEM）在处理非单纯形网格（如多面体）时往往面临困难。

2. 方法论 (Methodology)

本文提出并分析了一种基于 离散 de Rham (Discrete de Rham, DDR) 复形的高阶多面体格式。

• 混合变分公式：
  • 采用混合形式，引入位移场 $u$ 和拉格朗日乘子 $\lambda$（代表接触面上的牵引力）。
  • 拉格朗日乘子用于强制接触条件（非穿透和摩擦约束）。
• 离散空间构造：
  • 位移空间 ($X_{grad,h}$)：基于 DDR 复形的高阶版本。自由度（DOFs）附着在每个单元（体）、面、边和顶点上。
    • 允许位移在裂缝处不连续（通过裂缝两侧的独立自由度实现），但在非裂缝网格内部保持连续性。
    • 利用这些自由度可以重构出分片二次（piecewise quadratic）的位移近似场。
  • 拉格朗日乘子空间 ($M_h$)：在裂缝面上定义为分片常数向量。
• 算子重构：
  • 定义了局部和全局的梯度重构算子 ($G^1$) 和位移重构算子 ($\Upsilon^2$)，用于计算应变和应力。
  • 定义了裂缝面上的位移跳跃算子 $[[\cdot]]_\sigma$。
• 离散方案：
  • 包含一个稳定化双线性形式 $S_h$，用于控制单元内部重构与自由度之间的差异。
  • 通过半光滑牛顿算法（Semi-smooth Newton algorithm）求解由此产生的非线性方程组。

3. 主要理论贡献 (Key Contributions)

• 离散 Korn 不等式：证明了适应于 DDR 类型空间且考虑裂缝存在的离散 Korn 不等式。证明过程借鉴了低阶方法，但因引入了二次面位移重构而存在显著差异。
• 离散 inf-sup 条件：建立了一个基于非标准 $H^{-1/2}$ 范数的离散 inf-sup 条件，连接了离散拉格朗日乘子空间和离散位移空间，确保了混合格式的稳定性。
• 锁死自由（Locking-free）特性：
  • 设计了一种特殊的插值算子（修正了标准 DDR 插值算子的单元值），使其满足 Fortin 性质（与连续和离散散度算子交换）。
  • 关键成果：证明了该方法在准不可压缩极限（泊松比 $\nu \to 0.5$）下，对原始变量（位移）的误差估计是鲁棒的，即不会出现体积锁死（volumetric locking）现象。
• 误差分析：
  • 提供了完整的先验误差分析。
  • 对于 $H^{3/2+r}$ 正则性的解，证明了能量范数下的收敛率为 $O(h^{1/2+r})$。
  • 在无摩擦情况下 ($g=0$)，理论允许 $r$ 接近 1；在有摩擦情况下 ($g \neq 0$)，理论限制 $r \le 1/2$，但数值实验表明高阶精度依然保持。

4. 数值实验结果 (Results)

作者通过一系列数值实验验证了理论分析：

• 收敛性测试：
  • 在摩擦无接触 ($g=0$) 和 Tresca 摩擦 ($g \neq 0$) 两种情况下，使用制造解（Manufactured Solution）进行测试。
  • 在笛卡尔、四面体和六面体（Hexa-cut）网格上，位移梯度的收敛率接近 2 阶（优于理论预测的 1.5 阶），拉格朗日乘子在 $L^2$ 范数下收敛率为 1 阶。
  • 高阶方案在达到相同精度时，所需的自由度数量显著少于低阶方案（如参考文献 [32] 中的方法）。
• 不可压缩性测试：
  • 测试了拉梅系数 $L$ 从 1 增加到 $10^6$（对应 $\nu \to 0.5$）的情况。
  • 结果显示误差幅度保持不变，证实了方法的锁死自由特性。
  • 刚度矩阵的条件数随 $1/(0.5-\nu)$ 几乎线性增长，表明数值稳定性良好。
• 复杂裂缝网络 (DFM 模型)：
  • 在一个包含相交裂缝的 3D 离散裂缝 - 基质模型中，展示了高阶方法在捕捉接触状态（粘着/滑动）和接触过渡区方面的优越性。
  • 高阶方法产生的法向牵引力分布更平滑，且收敛所需的牛顿迭代次数更少。

5. 意义与影响 (Significance)

• 几何灵活性：该方法能够处理地质应用中常见的非单纯形、非结构化多面体网格，克服了传统有限元在处理复杂裂缝网络时的局限性。
• 高阶精度与效率：作为首个应用于接触力学的 DDR 方法，它展示了高阶多项式重构带来的显著精度提升和计算效率（更少的自由度达到同等精度）。
• 鲁棒性：解决了接触力学中常见的不可压缩材料锁死问题，这对于模拟岩石力学等准不可压缩材料至关重要。
• 理论突破：建立了适用于裂缝网络的高阶离散 inf-sup 条件和 Korn 不等式，为后续更复杂的摩擦模型（如库仑摩擦）和更广泛的接触力学问题提供了坚实的理论基础。

总结：这篇论文成功地将高阶离散 de Rham 格式扩展到了具有 Tresca 摩擦的接触力学问题中，不仅提供了严格的数学证明（稳定性、收敛性、锁死自由），还通过广泛的数值实验证明了其在处理复杂裂缝网络和不可压缩材料时的实用性和优越性。