用于三维不可压缩线弹性分析的有效重叠有限元

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：An effective overlapping finite element for three-dimensional incompressible linear elasticity

　　

本文继续探讨用于不可压缩介质分析的重叠有限元（Overlapping Finite Elements, OFE）。此前，研究人员已考察了离散化于四边形单元上的二维纳维 - 斯托克斯（Navier-Stokes）方程的重叠单元解法 [36]。基于该研究成果，研究人员将工作推进至三维分析，重点关注离散化于四面体单元上的线弹性问题。研究再次证实这些单元是稳定的（即它们满足 inf-sup 条件），从而排除了出现虚假压力解的可能性。然而，真正贡献于重叠单元有效性的关键在于其所需的自由度（Degrees of Freedom, DoF's）更少，从而导致矩阵规模更小、带宽更紧凑，且相较于其他常用于不可压缩介质分析的经典有限元，其总求解时间更短。正如文中所示，所提出的重叠单元特别优于广泛使用的泰勒 - 胡德（Taylor-Hood）单元。这种有效性的提升至关重要，因为在三维分析中，要获得精确结果通常需要的自由度数及求解时间非常高。

有限元方法多年来已成功应用于求解固体力学问题（包括线性和非线性问题）。尽管该领域已有许多卓越进展，但现有研究焦点集中在与传统有限元方法形成对比的无网格法（Meshfree methods），特别是有限球法（Method of Finite Spheres, MFS）。虽然 MFS 能提供精确解，但其基函数的非多项式性质导致数值积分过程极其昂贵，尤其在三维情形下，这阻碍了其在工程实践中的广泛应用。为了解决传统有限元对高质量网格的依赖与无网格法高昂计算成本之间的矛盾，重叠有限元（OFE）方法应运而生。该方法旨在结合两者的优点：既使用多项式基函数以简化数值积分，又能在低质量或畸变网格上保持良好性能。早期的 OFE 研究曾尝试混合使用 MFS 基函数，但随后的发展彻底转向了完全基于多项式的插值方案，显著提高了计算效率并扩展至三维问题。先前的研究已证实 OFE 在二维不可压缩流体动力学（Navier-Stokes 方程）中的有效性，满足了 Babuska-Brezzi 或 inf-sup 稳定性条件。

本研究由 Williams L. Nicomedes 和 Klaus-Jürgen Bathe 开展，旨在将 OFE 方法进一步拓展至三维不可压缩线弹性介质的分析中。研究人员利用最初在四面体单元上引入的重叠有限元概念，分别构建用于近似位移场和压力场的离散化空间，旨在建立一种既稳定又无虚假解的数值格式。研究的核心目标是验证该方法在三维情形下的稳定性、精度以及相对于传统有限元（如 Taylor-Hood 单元）在计算效率上的显著优势。研究结论表明，所提出的四面体重叠单元对不仅满足 inf-sup 稳定性条件，有效避免了体积锁死（locking）现象，而且在达到同等精度的前提下，所需的自由度数量远少于传统单元，生成的线性代数方程组规模更小且矩阵带宽更窄，从而大幅降低了存储需求和处理时间。这一发现对于解决三维大规模计算中资源消耗巨大的难题具有重要意义。该论文发表于《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》。

为开展此项研究，研究人员主要采用了以下关键技术方法：首先，基于重叠有限元范式，构建了用于三维域离散化的四面体单元网格，其中计算域被划分为四面体单元集合；其次，在单元内部构造了完全多项式性质的基函数，用于近似位移场（采用一阶 Lagrange 基函数结合局部富集空间）和压力场（采用常数或低阶多项式），摒弃了早期研究中使用的 Shepard 函数或非多项式基函数，以确保数值积分的高效性；再次，应用混合有限元公式（Mixed Formulation），针对不可压缩介质的特性，分别定义位移和压力的插值空间，并通过 inf-sup 测试验证其稳定性；最后，利用立方体、圆柱体及环形体等几何模型，通过对比不同网格畸变程度下的数值解与解析解，评估方法的精度、收敛性及对网格质量的鲁棒性。研究未涉及具体的生物样本队列或化学试剂操作，纯粹为计算力学领域的数值方法学研究。

研究结果部分详细展示了重叠有限元在三维不可压缩线弹性问题中的表现。
在“数值属性”部分，研究人员通过立方体模型验证了 OFE 离散化的可解性、带宽及所需自由度数量。结果表明，与传统有限元相比，OFE 方法生成的稀疏矩阵具有更紧凑的带宽，且达到相同精度所需的自由度显著减少，这意味着求解线性方程组的操作次数更少，计算效率更高。
在“稳定性”部分，研究重点考察了离散格式的 inf-sup 稳定性。通过理论分析与数值测试证明，所提出的 OFE1/OFE0 单元对（即位移场采用一阶近似，压力场采用零阶近似）满足 inf-sup 条件。这一结果至关重要，因为它从理论上保证了数值解的唯一存在性，并有效排除了不可压缩介质计算中常见的体积锁死效应及虚假压力振荡模式。
在“误差分析”部分，研究人员以圆柱体为模型，对比了 OFE1/OFE0 解与传统 P2/P1（二次位移/一次压力）Taylor-Hood 单元解的精度，特别是在网格发生畸变时的表现。结果显示，即使在网格质量较差的情况下，OFE 方法仍能保持较高的预测精度，其误差收敛率与理论预期一致，且对网格畸变表现出极强的不敏感性，优于传统单元。
在“附加实例：环形问题”部分，研究进一步考察了一个复杂的环形几何体（由两个同心圆柱截取四分之一得到）。分析再次证实，当应用于具有复杂边界和非结构化网格的问题时，重叠有限元方法依然能够给出精确且稳定的位移与压力分布，进一步佐证了该方法在处理实际工程几何形状时的有效性和鲁棒性。

讨论与结论部分总结指出，本研究成功探索并验证了一种专为三维不可压缩线弹性介质分析设计的新型重叠有限元。该单元基于完全多项式的基函数，极大地简化了弱形式的数值积分过程。其最显著的优势在于极高的计算效率：相较于广泛使用的 Taylor-Hood 等传统有限元，该重叠单元在保持甚至提高精度的同时，大幅减少了自由度数量和矩阵带宽。这种特性使得该方法在处理大规模三维工程问题时具有优异的可扩展性，能够显著降低对存储容量和计算能力的需求。此外，该方法对网格畸变的高度不敏感性，降低了对网格生成质量的严苛要求，为复杂几何域的自动化网格划分提供了更有力的工具。综上所述，重叠有限元方法为解决三维不可压缩介质力学分析中的计算瓶颈提供了一条高效、稳定且精确的新途径。