[論文レビュー] Volumetric expressions of the shape gradient of the compliance in structural shape optimization

Matteo Giacomini, Olivier Pantz+1 more|arXiv (Cornell University)|Jan 20, 2017

Topology Optimization in Engineering参考文献 31被引用数 8

ひとこと要約

本稿では、線形弾性における剛性の形状勾配について、純粋な変位および二重混合有限要素法の両形式を用いて体積積分表現を導出する。数値シミュレーションにより、メッシュ変形下でも表面ベースの手法に比べて、体積表現がより精度高く安定した降下方向を提供することを示し、メッシュ歪みや勾配近似誤差に起因する形状最適化の課題を浮き彫りにする。

ABSTRACT

In this article, we consider the problem of optimal design of a compliant structure under a volume constraint, within the framework of linear elasticity. We introduce the pure displacement and the dual mixed formulations of the linear elasticity problem and we compute the volumetric expressions of the shape gradient of the compliance by means of the velocity method. A preliminary qualitative comparison of the two expressions of the shape gradient is performed through some numerical simulations using the Boundary Variation Algorithm.

研究の動機と目的

• 線形弾性における剛性最小化のための体積表現形状勾配を導出すること。境界のみに依存する表現を避けることを目的とする。
• 線形弾性問題の純粋な変位および二重混合形式から導出される体積形状勾配の数値的性能を比較すること。
• 離散化された体積形状勾配を用いた境界変動アルゴリズム（BVA）によって計算された降下方向の安定性および精度を調査すること。
• 形状最適化における数値的不安定要因、例えば負の方向微分が得られているにもかかわらず降下でない行動を特定すること。
• 構造形状最適化の耐障害性を向上させるために、認証済み降下アルゴリズムおよびレベルセット法の有用性を示唆すること。

提案手法

• 変分的枠組み内での速度法を用いて、境界積分を体積積分に変換するドメイン写像を介して、剛性汎関数の形状勾配を導出する。
• 変位を主変数とする線形弾性の純粋な変位形式を適用し、対応する体積形状勾配表現を導出する。
• ヘリンジャー＝ライスナー形式およびスティュレス対称性を弱く強制する変種を含む、二つの二重混合形式を導入し、両者とも体積形状勾配表現を導出する。
• 境界変動アルゴリズム（BVA）を用いて、計算された形状勾配に従い、降下方向に従ってドメインを繰り返し変形する。メッシュ更新は降下方向によって駆動される。
• 事後誤差推定器および認証済み降下アルゴリズム（CDA）を用いて、離散形状勾配における非降下行動を検出し、是正する。
• 初期形状が異なる2次元テストケースにおける数値シミュレーションを実施し、純粋な変位形式と二重混合形式の収束性および安定性を比較する。

実験結果

リサーチクエスチョン

• RQ1 線形弾性における剛性の形状勾配を、表面積分を回避する完全な体積表現としてどのように表現できるか？
• RQ2 純粋な変位形式と二重混合形式から導出される体積形状勾配の、数値的利点と欠点は何か？
• RQ3 境界変動アルゴリズム（BVA）が、負の方向微分が得られているにもかかわらず、真の降下方向を生成しないのはなぜか？
• RQ4 メッシュ歪みや速度場の正則性が低いことが、形状最適化アルゴリズムの収束性にどの程度影響を及えるか？
• RQ5 離散化された体積形状勾配を用いる場合、認証済み降下アルゴリズムは最適化の信頼性を向上させ得るか？

主な発見

• 純粋な変位形式および二重混合形式から導出された体積形状勾配表現は、連続的状況下で古典的な表面ベース形状勾配と数学的に同等である。
• 数値的シミュレーションの結果、特に初期段階の反復で、二重混合形式が純粋な変位形式に比べ、より安定的かつ正確な降下方向を提供することが示された。
• ⟨dₕL(Ω), θₕ⟩ < 0 という負の方向微分が得られているにもかかわらず、ペナルティ付き剛性汎関数 L(Ω) が増加する場合（例：表1cの10.64 → 10.66）があり、θₕ が真の降下方向でないことが示された。
• 純粋な変位形式を用いたテストケース（c）では、30反復で剛性が10.64から4.66に減少したが、振動的挙動を示し、降下方向の不安定性が確認された。
• テストケース（d）において二重混合形式は、30反復で剛性を11.01から3.37に減少させ、収束性が良好であったが、依然として振動が観察された。
• 本研究は、標準的なBVAが勾配近似の悪さやメッシュ歪みの影響により失敗する可能性があるため、離散形状勾配の降下性を検証するための認証済み降下アルゴリズムの必要性を確認した。