🔢 mathematics

A Hybridized Staggered Discontinuous Galerkin--Mixed Finite Element Method for Strain Gradient Elasticity

本文提出了一种用于应变梯度弹性学的混合交错间断伽辽金—混合有限元法，该方法通过将四阶问题重构为具有强对称性强制执行的一阶系统，并采用一种最小化全局耦合的混合化方案，从而实现了参数鲁棒稳定性与最优收敛性。

原作者： Bohan Yang, Eric T. Chung

发布于 2026-09-15

📖 1 分钟阅读🧠 深度阅读

原作者： Bohan Yang, Eric T. Chung

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

想象一个材料如何弯曲和拉伸的规则取决于观察尺度大小的世界。一个多世纪以来，工程师们一直依赖于经典弹性理论，这是一套将材料视为连续、无缝织物的法则。这些法则在应用于桥梁、建筑和大型机械时表现得非常完美。然而，当应用于现代技术中的微观结构（例如微芯片中的细微导线或医疗传感器中的薄梁）时，这些法则就会失效。在这些尺度下，材料会表现出“尺寸效应”，即仅仅因为体积微小而表现出不同的行为。一根细丝可能会比经典物理学预测的扭转得更厉害，或者一根微型梁可能会以意想不到的方式抵抗弯曲。为了解释这一现象，科学家们开发了应变梯度弹性理论，这是一种更先进的理论，它在方程中增加了一个特定的“长度尺度”。这个额外的参数解释了材料的内部结构，使得数学模型能够捕捉到那些神秘的尺寸相关行为。

然而，挑战在于这些高级方程在计算机上极难求解。它们是四阶问题，这意味着涉及材料如何变形以及这种变形如何随空间变化的复杂关系。用于求解这些问题的传统方法在处理近乎不可压缩（如橡胶）的材料或长度尺度非常小时经常会遇到困难，导致结果不准确，或者在物体边缘产生“数值边界层”，即计算机在边缘附近产生的虚假误差。研究人员尝试过各种方法，但许多方法需要对参数进行精细调节，或者无法提供计算总应力的直接途径，而总应力是工程师设计安全微型设备的至关重要的量。

在一项新的研究中，研究人员 Bohan Yang 和 Eric T. Chung 提出了一种解决该问题的新方法。他们开发了一种混合方法，结合了两种强大的数学技术：交错间断伽辽金法（staggered discontinuous Galerkin method）和混合有限元法（mixed finite element method）。他们的做法并非仅仅尝试求解材料的位移，而是引入了几个新变量来表示材料内的不同类型的应力。他们将总应力以及内部力的缩放版本作为独立的未知数进行处理。通过这种方式，他们将困难的四阶问题重新构建为一个由更简单的一阶方程组成的系统。这种转变使他们能够使用特定的数学空间，这些空间能够自然地处理问题的复杂性，而无需人工修复或惩罚项。

他们的核心创新在于如何组织计算。他们将材料的几何形状划分为三角形网格，然后在其中创建第二个重叠网格。在这个双重结构上，他们执行不同的连续性规则，使数学片段能够以保持应力物理对称性的方式组合在一起。这种设计确保了无论材料的长度尺度多么小，或者材料多么接近不可压缩，该方法都能保持稳定和准确。他们工作的一个关键特征是，他们可以从数学上证明该方法的稳健性：即当参数变化时，准确度不会下降，而这正是其他方法常见的失效点。此外，他们还引入了一种“混合化”方案，允许他们在每个小三角形内局部求解大多数变量，仅留下极少数全局变量进行统一求解。这极大地降低了计算成本，使该方法在效率上足以投入实际应用。

研究人员通过数值实验测试了他们的方法，模拟了具有平滑解以及具有强边界层（即物体边缘附近行为的剧烈变化）的情景。在每种情况下，他们的方法都达到了预期的准确率，并在广泛的材料属性（包括极端的刚度和不可压缩性）下保持了稳定性。与在长度尺度极小时会遇到困难的传统位移法不同，他们的方法在不产生人工误差的情况下保持了最优准确度。他们还证明了该方法能直接且准确地计算总应力，而这在其他技术中通常难以获得。虽然目前的工作侧重于具有特定边界条件的二维形状，但这一方法的成功为模拟复杂的微观尺度材料提供了一条充满前景的路径。这项研究证实，通过重新思考这些方程的结构和求解方式，是有可能克服应变梯度弹性长期存在的难题，从而为下一代微工程提供可靠的工具。

问题陈述

本文研究了应变梯度弹性（SGE）模型（特别是 Aifantis 及其合作者提出的公式化模型）的数值近似。该模型是一个由材料特征长度参数 $\iota \in (0, 1]$ 以及拉梅常数 $\lambda$ 和 $\mu$ 控制的四阶奇异摄动问题。其控制方程 $-(1 - \iota^2\Delta)(\mu\Delta u + (\lambda + \mu)\nabla\text{div} u) = f$ 在离散化方面面临两个主要挑战：

高阶连续性： 由于这是一个四阶问题，标准的有限元方法需要 $H^2$ 连续性。虽然存在 $C^1$ 连续单元，但其构造复杂且计算成本高昂。非连续方法往往难以保证稳定性，或者需要仔细调整惩罚参数。
参数鲁棒性： 该问题相对于 $\iota$ 是奇异摄动的。当 $\iota \to 0$ 时，模型退化为经典弹性力学，但高阶狄利克雷边界条件（ $\partial_n u = 0$ ）无法由经典解满足，从而导致一个厚度为 $O(\iota)$ 的物理边界层。在离方面通过强行施加该条件的位移类方法会产生由网格尺寸 $h$ 控制的“数值边界层”，导致在 $\iota \ll 1$ 时收敛阶降至 $O(h^{1/2})$ 。此外，当 $\lambda \to \infty$ （不可压缩极限）时，方法可能会出现体积自锁现象。

方法论

作者提出了一种杂化交错间断伽辽金（SDG）–混合有限元方法。该方法包含以下步骤：

一阶重构： 通过引入总应力 $\sigma$ 以及两个辅助变量：缩放后的柯西应力 $\tilde{\sigma}$ 和缩放后的超应力 $\chi$ ，将四阶方程重构为一个一阶系统。这产生了一个包含位移 $u$ 、总应力 $\sigma$ 、缩放柯西应力 $\tilde{\sigma}$ 和缩放超应力 $\chi$ 的四场系统。对称性在所有三个应力张量上被强行施加。
离散空间：
- 位移 $u$ 和总应力 $\sigma$ 使用交错间断伽辽金（SDG）空间进行离散。这些空间利用原网格的 Alfeld 分裂来创建一个对偶网格，在原网格边上强制执行 $u$ 的连续性，并在对偶边上强制执行 $\sigma$ 的法向连续性。
- 缩放应力 $\tilde{\sigma}$ 和 $\chi$ 分别使用 Raviart–Thomas (RT) 对和间断伽辽金空间进行离散，并施加对称性。
杂化： 为了减小全局系统规模，总应力和缩放超应力的单元间法向连续性被放宽，并通过定义在网格骨架（对偶边）上的拉格朗日乘子进行弱实现。这允许进行局部静态凝聚，其中内部自由度被局部消除，仅留下作为全局耦合未知数的乘子。
边界条件： 高阶狄利克雷边界条件被自然地纳入变分形式中，而不是强行施加在离散空间上，从而避免了数值边界层。

核心贡献与理论结果

论文提供了支持所提方法的严谨理论分析：

稳定性： 作者建立了对称子空间上的 inf-sup 条件以及离散 Korn 不等式。这些工具被用于证明该方案的代数稳定性。
收敛性：
- 在充分正则性的条件下，该方法在位移和应力方面均达到最优收敛阶。
- 至关重要的是，误差估计在 $\iota$ 和 $\lambda$ 方面都是一致的。即使在 $\iota \ll 1$ 时，该方法仍能保持最优收敛阶，避免了位移类方法中出现的性能退化。
- 分析明确追踪了对解的正则性和参数的依赖关系。
杂化： 论文开发了杂化方案，证明了两个乘子（近似于 $u$ 的迹和 $\iota\tilde{\sigma}$ ）满足它们各自的 inf-sup 条件，从而确保了简化后的全局系统的稳定性。
不可压缩极限： 针对 $\lambda \to \infty$ （不可压缩极限）的情况提供了一种改进方法，涉及应力的分解，以处理实现良定性所需的全局归一化约束。

数值结果

报告了两组数值实验以验证理论：

光滑解： 一个具有光滑、无散解的测试用例证实了位移在 $L^2$ 和离散 $H^1$ 范数下的最优收敛阶，以及总应力的最优收敛阶。结果展示了在广泛的 $\lambda$ 范围（从 1 到 $10^8$ ）和 $\iota$ 范围（从 1 到 $10^{-8}$ ）内的鲁棒性。
边界层： 一个具有强边界层的测试用例（其中 SGE 问题的精确解在边界附近与经典解不同）证实了即使在 $\iota$ 非常小（ $10^{-6}, 10^{-8}$ ）的情况下，该方法仍能保持最优收敛阶和参数鲁棒性。

意义

该论文声称提出了首个得到严谨分析支持的应变梯度弹性四场混合方案。其意义在于：

避免数值边界层： 通过将高阶边界条件以变分方式纳入，该方法在 $\iota$ 很小时仍能保持最优收敛阶，解决了位移类方法中已知的一个难题。
参数鲁棒性： 该方案在长度参数 $\iota$ 和拉梅常数 $\lambda$ 方面都是稳定且一致收敛的，同时解决了奇异摄动和不可压缩性问题。
直接应力近似： 与位移类方法不同，这种混合公式直接提供了总应力的近似，而总应力是一个具有实际意义的量，无需通过高阶导数的后处理获得。
计算效率： 杂化版本将全局系统减少到仅在网格骨架上的两个乘子，与非杂化的四场系统相比，显著降低了计算成本。

目前该工作局限于二维及齐次双向夹紧边界条件，向三维及一般边界条件的扩展被列为未来的工作。

类似论文