[論文レビュー] Automated shape and thickness optimization for non-matching isogeometric shells using free-form deformation

Han Zhao, David Kamensky+2 more|arXiv (Cornell University)|Aug 2, 2023

Advanced Numerical Analysis TechniquesEngineering参考文献 53被引用数 3

ひとこと要約

本稿では、NURBSに基づくアイソジオメトリック解析（IGA）を用いて、非一致するアイソジオメトリックシェル構造の自動形状および厚さ最適化のための統一的自由形式変形（FFD）ベースのフレームワークを提案する。FFDをラグランジュ抽出およびFEniCSベースのコード生成と統合することで、パッチ界面を越えて滑らかで連続的な設計更新が可能となり、自動的解析的勾配を備えた効率的な勾配ベース最適化が実現され、eVTOL機の翼最適化において最大95.09％の内部エネルギー低減を達成した。

ABSTRACT

Isogeometric analysis (IGA) has emerged as a promising approach in the field of structural optimization, benefiting from the seamless integration between the computer-aided design (CAD) geometry and the analysis model by employing non-uniform rational B-splines (NURBS) as basis functions. However, structural optimization for real-world CAD geometries consisting of multiple non-matching NURBS patches remains a challenging task. In this work, we propose a unified formulation for shape and thickness optimization of separately-parametrized shell structures by adopting the free-form deformation (FFD) technique, so that continuity with respect to design variables is preserved at patch intersections during optimization. Shell patches are modeled with isogeometric Kirchhoff--Love theory and coupled using a penalty-based method in the analysis. We use Lagrange extraction to link the control points associated with the B-spline FFD block and shell patches, and we perform IGA using the same extraction matrices by taking advantage of existing finite element assembly procedures in the FEniCS partial differential equation (PDE) solution library. Moreover, we enable automated analytical derivative computation by leveraging advanced code generation in FEniCS, thereby facilitating efficient gradient-based optimization algorithms. The framework is validated using a collection of benchmark problems, demonstrating its applications to shape and thickness optimization of aircraft wings with complex shell layouts.

研究の動機と目的

• 複数の非一致NURBSパッチから構成される複雑なシェル構造における形状および厚さ最適化の課題に対処すること。
• 最適化中にパッチ界面で幾何的連続性を保証し、形状の不連続性を回避すること。
• NURBS基底関数を用いたアイソジオメトリック解析（IGA）により、CAD幾何と解析モデルのシームレス統合を実現すること。
• 効率的な勾配ベース最適化を支援するため、解析的感度計算を自動化すること。
• 手動メッシュ生成を必要とせず、実世界の航空機翼構造における設計－解析－最適化ワークフローを簡素化すること。

提案手法

• 自由形式変形（FFD）を用いて、最適化中に非一致NURBSパッチ間で連続性を保つ形状および厚さの変動をパラメータ化する。
• ラグランジュ抽出を用いて、FFDブロックの制御点をシェルパッチの制御点にマッピングし、一貫性のある幾何更新を可能にする。
• アイソジオメトリック解析における非一致キルヒホッフ＝ローヴェシェルフパッチ間の連続性を弱く強制するためのペナルティベースの結合法を適用する。
• FEniCS PDEライブラリ内の既存の有限要素アセンブリ手続きを活用し、同じ抽出行列を用いてIGAを実行する。
• FEniCSにおける高度なコード生成を活用して解析的導出計算を自動化し、効率的な勾配ベース最適化を実現する。
• 正則化係数を用いて形状と厚さの更新をバランスさせ、設計における悪条件化や振動を防止する。

実験結果

リサーチクエスチョン

• RQ1 別々にパarameter化された非一致NURBSパッチからなるシェル構造に対して、形状と厚さ最適化をどのように統合できるか？
• RQ2 適合メッシュを必要とせずに、自由形式変形（FFD）は最適化中にパッチ界面で幾何的連続性を維持できるか？
• RQ3 ラグランジュ抽出とFEniCSベースのコード生成は、アイソジオメトリック最適化における感度計算をどの程度自動化できるか？
• RQ4 正則化係数は、複雑な翼構造における最適化の収束性と設計品質にどのような影響を与えるか？
• RQ5 提案されたフレームワークは、実世界の航空機翼幾何形状において顕著な構造的性能向上（例：内部エネルギー低減）を達成できるか？

主な発見

• 提案されたFFDベースの最適化フレームワークは、正則化係数$10^{-5}$を用いたeVTOL翼において最大95.09％の内部エネルギー低減を達成し、顕著な性能向上を示した。
• 正則化係数$10^{-3}$を用いた場合、83.23％の内部エネルギー低減が達成され、形状と厚さの両方の更新が顕著に寄与するバランスの取れた設計が得られた。
• 高い正則化係数（$ackslashlambda = 1$および$0.1$）を用いた最適化では、形状変更に対する高いペナルティが加わったため、厚さの更新が支配的となり、形状変更は最小限に抑えられた。
• 低正則化係数（$10^{-4}$および$10^{-5}$）では、顕著な形状変更と高いエネルギー低減が得られたが、翼端部で振動が生じ、悪条件化の兆候を示した。
• フレームワークは最適化全体を通じてパッチ界面での幾何的連続性を正常に維持し、独立したパッチパラメータ化に起因する不連続性を回避した。
• FEniCS上に構築されたオープンソースPythonライブラリGOLDFISHにより、再現可能で自動的かつ効率的な多分野最適設計が、複雑なシェル構造に対して可能になった。