[论文解读] Enhancing SfePy with Isogeometric Analysis

Robert Cimrman|arXiv (Cornell University)|Dec 19, 2014

Advanced Numerical Analysis Techniques参考文献 5被引用 5

一句话总结

本文通过引入等几何分析（IGA）扩展了开源 SfePy 有限元框架，实现了更高阶的光滑性与直接的 CAD 兼容几何处理。通过基于 Bézier 提取的 NURBS 近似，IGA 在保持更少自由度的情况下实现了比 FEM 更精确的解，已在二维拉普拉斯与泊松问题及三维弹性问题上得到验证，且网格生成开销极低。

ABSTRACT

In the paper a recent enhancement to the open source package SfePy (Simple Finite Elements in Python, http://sfepy.org) is introduced, namely the addition of another numerical discretization scheme, the isogeometric analysis, to the original implementation based on the nowadays standard and well-established numerical solution technique, the finite element method. The isogeometric removes the need of the solution domain approximation by a piece-wise polygonal domain covered by the finite element mesh, and allows approximation of unknown fields with a higher smoothness then the finite element method, which can be advantageous in many applications. Basic numerical examples illustrating the implementation and use of the isogeometric analysis in SfePy are shown.

研究动机与目标

• 将等几何分析（IGA）集成到现有的 SfePy 有限元框架中，以支持更高阶光滑性与直接的 CAD 几何输入。
• 解决高精度模拟中对全局 C^p-连续近似的需求，特别是针对从头算电子结构计算。
• 通过二维与三维基准 PDE 问题，评估 IGA 在 SfePy 中的可行性与性能表现。
• 通过 Bézier 提取技术实现单一计算框架内从 FEM 到 IGA 的无缝过渡，实现基函数的高效处理。

提出的方法

• 采用基于 NURBS 的几何与近似方法实现 IGA，从而在整个计算域内获得平滑的高阶基函数。
• 利用 Bézier 提取算子将 NURBS 基函数映射为与 SfePy 现有 FEM 基础设施兼容的形式。
• 通过 NURBS 形状函数与高斯积分对 PDE 的弱形式进行组装，积分在参数化单元上进行。
• 通过在参数域选定区域施加本质（狄利克雷）约束来应用边界条件。
• 该框架支持标量与矢量值 PDE，包括扩散与弹性问题。
• 采用直接求解器与牛顿法求解非线性问题，网格与几何通过 Bézier 控制点定义。

实验结果

研究问题

• RQ1 能否在不显著改变架构的前提下，将等几何分析有效集成到类似 SfePy 的现有 FEM 框架中？
• RQ2 在简单区域上的标准 PDE 问题中，IGA 与 FEM 在自由度数量与解的精度方面有何差异？
• RQ3 在有限元环境中支持基于 NURBS 的离散化面临哪些计算与实现挑战？
• RQ4 与 FEM 相比，IGA 在保持解的精度前提下，能在多大程度上减少自由度数量？
• RQ5 IGA 是否能原生支持复杂边界条件与变系数 PDE，且在 SfePy 框架中实现？

主要发现

• SfePy 中的 IGA 实现成功支持了标量与矢量 PDE，包括拉普拉斯方程、泊松方程与线性弹性方程。
• 在二维拉普拉斯问题中，IGA 解仅使用了 20 个自由度，而使用结构化网格的 FEM 解则需要 1363 个自由度。
• IGA 解实现了平滑的、全局 C^p-连续的近似，无需进行域三角剖分或网格生成。
• 该方法与基于 NURBS 的 CAD 式几何描述具有原生兼容性，消除了对独立网格生成的需求。
• 在三维弹性问题示例中，曲面上空间变化的位移边界条件被准确施加。
• Bézier 提取技术实现了 IGA 与现有 SfePy FEM 流程的相对无缝集成，同时保持了模块化与可扩展性。