[論文レビュー] 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要素について、従来は2つの三角形に分割する必要があったが、その分解を回避することで統合要素を削減する。
• 必要に応じて局所的グリッド細分化を実施し、要素タイプの一貫性を保ちながら計算負荷を低減する。
• 統合点の分布が従来手法よりもより滑らかであることを保証し、数値的安定性を向上させる。

実験結果

リサーチクエスチョン

• RQ1 どのようにして、トリムされた幾何形状の等参画的解析における統合の複雑さを低減できるか？
• RQ2 統合点のより均一な分布は、解の正確性と条件数の改善に寄与するか？
• RQ3 タイプBのトリム要素について、三角形分解を排除することで計算効率にどのような影響があるか？
• RQ4 解の正確性を損なわずに、統合要素および統合点の数をどの程度まで削減できるか？
• RQ5 提案手法は、従来の手法と比較して、誤差および剛性行列の条件数の観点でどの程度優れているか？

主な発見

• 提案手法は、[13]の手法と比較して、特にタイプAおよびB要素において、統合要素および統合点の数を最大50％まで削減した。
• 細分化されたグリッドにおいて、剛性行列の条件数が約50％低減され、数値的安定性が向上した。
• $20\times20$グリッドにおいて、数値解の$L^2$誤差は[13]の手法と比較してほぼ半減した。
• 統合点の分布がより滑らかであり、解の正確性の向上と条件数の低減に寄与した。
• タイプB要素では、三角形分解を回避したため、統合要素が半分に削減され、関連する数値積分の複雑さが解消された。
• 4つのベンチマーク例を通じて、複数の突出部や穴を有する複雑なトリムドドメインにおいても、高い正確性と効率性を維持した。