[论文解读] A finite element implementation of surface elasticity at finite strains using the deal.II library

Andrew McBride, Ali Javili+2 more|arXiv (Cornell University)|Jun 2, 2015

Elasticity and Material Modeling参考文献 1被引用 4

一句话总结

本文提出了一种基于开源 deal.II 库的完全非线性有限元实现，用于在有限大变形下模拟表面弹性，采用总拉格朗日格式，并从体积域和表面域精确计算切线贡献。该实现能够准确模拟有限大变形范围内的表面效应，通过基准测试验证，并为纳米尺度应用提供完整的源代码。

ABSTRACT

The potentially significant role of the surface of an elastic body in the overall response of the continuum can be described using the mature theory of surface elasticity. The objective of this contribution is to detail the finite element approximation of the underlying governing equations (both in the volume and on its surface) and their solution using the open-source finite element library deal.II. The fully-nonlinear (geometric and material) setting is considered. The nonlinear problem is solved using a Newton--Raphson procedure wherein the tangent contributions from the volume and surface are computed exactly. The finite element formulation is implemented within the total Lagrangian framework and a Bubnov-Galerkin spatial discretization of the volume and the surface employed. The surface is assumed material. A map between the degrees of freedom on the surface and on the boundary of the volume is used to allocate the contribution from the surface to the global system matrix and residual vector. The deal.II library greatly facilitates the computation of the various surface operators, allowing the numerical implementation to closely match the theory developed in a companion paper. Key features of the theory and the numerical implementation are elucidated using a series of benchmark example problems. The full, documented source code is provided.

研究动机与目标

• 开发有限大变形范围内表面弹性的有限元实现，解决现有数值方法在无限小变形限制下的局限性。
• 采用总拉格朗日框架与 Bubnov–Galerkin 空间离散化，提供完全非线性、几何与材料一致的公式化方法。
• 通过精确计算切线矩阵，整合表面与体积贡献，实现对纳米结构中表面效应的精确建模。
• 使实现开源并提供完整文档，通过 TBB 支持并行计算，提升复杂几何结构下的性能。
• 通过利用 deal.II 的共维数一能力，弥合理论表面弹性与实际数值模拟之间的鸿沟，实现对三维空间中曲面流形的建模。

提出的方法

• 有限元公式基于总拉格朗日框架，采用 Bubnov–Galerkin 方法对体积域和表面域进行空间离散化。
• 通过自由度映射将表面贡献映射到体积边界，实现表面与体积自由度之间的耦合。
• 采用牛顿-拉夫森格式求解非线性系统，精确计算来自体积和表面弱形式的切线贡献。
• 利用 deal.II 库的共维数一例程，处理嵌入三维空间中的二维曲面流形的几何与算子。
• 表面建模为具有类 neo-Hookean 自由能的材料膜，并扩展以包含表面张力，实现流体行为。
• 通过共享内存范式下的 TBB 库实现并行化，分布式内存扩展方案正在研究中。

实验结果

研究问题

• RQ1 如何在一致的有限元公式下，准确实现有限大变形范围内的表面弹性？
• RQ2 精确切线矩阵计算在确保非线性表面弹性问题中的收敛性与稳定性方面起到何种作用？
• RQ3 如何利用 deal.II 库的共维数一功能，对三维连续体中的复杂表面几何进行建模？
• RQ4 在网格细化过程中，表面与体积比增加带来的数值挑战是什么，如何缓解？
• RQ5 该实现的并行化程度如何，能否高效模拟具有复杂表面拓扑的现实纳米结构？

主要发现

• 该实现成功以高精度捕捉有限应变下的表面效应，通过涉及表面张力和曲率依赖变形的基准问题得到验证。
• 牛顿-拉夫森格式中采用精确切线贡献，确保了非线性迭代的二次收敛性与鲁棒性。
• deal.II 的共维数一特性能够精确高效地计算复杂曲面上的表面算子，即使在三维空间中的嵌入曲面也适用。
• 完整且文档齐全的源代码公开可用，有利于结果复现，并可轻松扩展至新材料模型或应用场景。
• 通过 TBB 实现的并行化显著提升了大规模问题的性能，未来在分布式环境中的扩展潜力巨大。
• 该框架能够模拟具有复杂表面几何的现实纳米结构，证明其在学术基准之外的实际应用价值。