[论文解读] An improved isogeometric analysis method for trimmed geometries

Jinlan Xu, Ningning Sun+3 more|arXiv (Cornell University)|Jul 2, 2017

Advanced Numerical Analysis Techniques参考文献 18被引用 3

一句话总结

该论文提出了一种改进的等几何分析（IGA）方法，用于处理修剪几何体，通过参数变换将修剪单元映射到参考单元，从而减少积分单元和积分点数量，相较于先前方法，数值误差和刚度矩阵的条件数降低50%，同时显著提升了计算效率和解的正则性。

ABSTRACT

Trimming techniques are efficient ways to generate complex geometries in Computer-Aided Design(CAD). In this paper, an improved isogeometric analysis(IGA) method for trimmed geometries is proposed. We will show that the proposed method reduces the numerical error of physical solution by 50% for simple trimmed geometries, and the condition number of stiffness matrix is also decreased. Furthermore, the number of integration elements and integration points involved in the solving process can be significantly reduced compared to previous approaches, drastically improving the computational efficiency for IGA problems on the trimmed geometry. Several examples are illustrated to show the effectiveness of the proposed approach.

研究动机与目标

• 解决修剪几何体等几何分析中的高计算成本和数值不准确性问题。
• 减少IGA中修剪曲面所需积分单元和积分点的数量。
• 改善修剪单元上积分点分布的正则性。
• 提高解的精度并降低刚度矩阵的条件数。
• 为复杂修剪CAD模型的IGA提供更高效、更鲁棒的框架。

提出的方法

• 通过参数变换将每个修剪单元（类型A、B、C）映射到参考单位正方形，以简化积分过程。
• 从物理修剪单元到参数参考域应用双射映射，保持几何形状并支持标准积分。
• 利用由控制点和节点向量导出的变换矩阵，将计算域映射到规则网格。
• 通过避免对B型单元进行三角剖分，减少积分单元数量，此前B型单元需拆分为两个三角形。
• 仅在必要时进行局部加密，保持单元类型一致性，降低计算负载。
• 确保积分点分布比以往方法更规则，从而提升数值稳定性。

实验结果

研究问题

• RQ1 如何降低修剪几何体等几何分析中的积分复杂度？
• RQ2 更规则的积分点分布能否提升解的精度和条件数？
• RQ3 消除三角剖分对B型修剪单元的计算效率有何影响？
• RQ4 在不牺牲解精度的前提下，积分单元和积分点的数量最多可减少多少？
• RQ5 与现有方法相比，所提方法在误差和刚度矩阵条件数方面表现如何？

主要发现

• 与文献[13]中的方法相比，该方法将积分单元和积分点数量减少了最多50%，尤其在A型和B型单元中效果显著。
• 在精细网格上，刚度矩阵的条件数降低了约50%，显著提升了数值稳定性。
• 在20×20网格上，数值解的$L^2$误差相比文献[13]的方法几乎减少了一半。
• 积分点分布更加规则，有助于提升解的精度并减少条件数问题。
• 对于B型单元，该方法避免了三角剖分，使积分单元数量减少一半，并消除了相关的高斯积分复杂性。
• 在包含多个凸起或孔洞的复杂修剪区域中，该方法仍能保持高精度和高效性，四个基准算例验证了这一点。