[论文解读] An Interface-enriched Generalized Finite Element Method for Levelset-based Topology Optimization

Sanne J. van den Boom, Jian Zhang+2 more|arXiv (Cornell University)|Mar 30, 2020

Topology Optimization in Engineering参考文献 74被引用 24

一句话总结

该论文提出一种界面增强型广义有限元法（IGFEM）与径向基函数（RBF）参数化水平集拓扑优化相结合的方法，实现无需重剖分的平滑、非像素化设计。该方法可实现精确的边界分辨率并直接施加本质边界条件，从而生成具有最小数值伪影的清晰黑白拓扑结构，已在基准柔度最小化问题上得到验证。

ABSTRACT

During design optimization, a smooth description of the geometry is important, especially for problems that are sensitive to the way interfaces are resolved, e.g., wave propagation or fluid-structure interaction. A levelset description of the boundary, when combined with an enriched finite element formulation, offers a smoother description of the design than traditional density-based methods. However, existing enriched methods have drawbacks, including ill-conditioning and difficulties in prescribing essential boundary conditions. In this work we introduce a new enriched topology optimization methodology that overcomes the aforementioned drawbacks; boundaries are resolved accurately by means of the Interface-enriched Generalized Finite Element Method (IGFEM), coupled to a levelset function constructed by radial basis functions. The enriched method used in this new approach to topology optimization has the same level of accuracy in the analysis as standard the finite element method with matching meshes, but without the need for remeshing. We derive the analytical sensitivities and we discuss the behavior of the optimization process in detail. We establish that IGFEM-based levelset topology optimization generates correct topologies for well-known compliance minimization problems.

研究动机与目标

• 解决基于密度的拓扑优化方法的局限性，如因网格依赖的材料界面导致的灰度伪影和阶梯效应。
• 克服现有富集FEM方法（如X/GFEM）的缺点，包括病态条件、边界条件弱约束以及界面附近应力过度估计的问题。
• 开发一种鲁棒的、无需重剖分的拓扑优化框架，在保持高分析精度的同时实现精确的边界追踪。
• 通过解耦设计与分析网格，实现本质边界条件的直接施加，并改善收敛行为。
• 证明IGFEM在基于水平集的拓扑优化中用于柔度最小化问题的可行性与有效性，为更复杂应用提供概念验证。

提出的方法

• 采用界面增强型广义有限元法（IGFEM），仅在被材料界面切割的单元中对有限元形函数进行富集，确保富集函数在原始网格节点处为零。
• 使用径向基函数（RBFs）对水平集函数进行参数化，将设计空间与固定分析网格解耦，从而实现平滑的边界演化。
• 通过RBF插值节点水平集值构建水平集函数作为有符号距离场，实现平滑且连续的边界表示。
• 应用IGFEM求解控制方程的弱形式（如弹性或导热方程），其中富集函数仅局限于被切割的单元。
• 通过保留原始网格节点的物理解释，直接在非匹配边界上强约束本质边界条件。
• 利用伴随法推导目标函数关于RBF控制点的解析灵敏度，实现高效的基于梯度的优化。

实验结果

研究问题

• RQ1 基于IGFEM的水平集拓扑优化能否在无需后处理或重剖分的情况下生成清晰、非像素化的设计？
• RQ2 基于RBF的设计网格与固定分析网格的解耦如何影响拓扑优化中的收敛性和数值稳定性？
• RQ3 IGFEM在多大程度上缓解了富集FEM中的常见问题，如病态条件、界面附近应力过度估计以及混合单元中的收敛性差？
• RQ4 当设计空间未与分析网格对齐时，该方法能否准确分辨细小特征并避免锯齿状伪影？
• RQ5 在基准问题中，该方法与标准SIMP、Ersatz及离散水平集方法相比，在性能和解质量方面表现如何？

主要发现

• 基于IGFEM的水平集拓扑优化方法由于采用平滑且贴合几何的水平集表示，可生成无需后处理的非像素化、黑白分明设计。
• 即使界面切割单元，该方法在无需重剖分的情况下仍能达到与匹配网格标准FEM相当的精度。
• 已推导并以低计算开销计算出解析灵敏度，支持高效的基于梯度的优化。
• 设计中的锯齿状伪影由粗略的分析网格和近似误差引起，可通过使用RBF限制设计复杂度或细化网格加以缓解。
• 该方法在标准柔度最小化基准问题（包括MBB梁和3D悬臂梁）中成功生成了正确拓扑，且收敛稳定。
• RBF的使用起到了滤波作用，减少了数值伪影，并通过允许每次迭代中更大的边界移动，使优化器更快收敛。