[論文レビュー] Overhang penalization in additive manufacturing via phase field structural topology optimization with anisotropic energies

Harald Garcke, Kei Fong Lam+2 more|arXiv (Cornell University)|Nov 28, 2021

Topology Optimization in Engineering被引用数 3

ひとこと要約

本稿では、積層造形におけるオーバーハングをペナルティ化する非等方的エネルギー関数を用いたフェーズフィールドトポロジー最適化法を提案する。これにより、外部補助具を必要としない自己支持構造が保証される。凸および非凸の非等方的性質を組み込むことで、部分微分法と漸近解析を用いてオーバーハング角度制約を強制し、特に強い非等方性下でより急勾配で自己支持性を有する構造が得られる。数値シミュレーションにより、キャンチレバー梁およびブリッジ設計の事例が検証された。

ABSTRACT

A phase field approach for structural topology optimization with application to additive manufacturing is analyzed. The main novelty is the penalization of overhangs (regions of the design that require underlying support structures during construction) with anisotropic energy functionals. Convex and non-convex examples are provided, with the latter showcasing oscillatory behavior along the object boundary termed the dripping effect in the literature. We provide a rigorous mathematical analysis for the structural topology optimization problem with convex and non-continuously-differentiable anisotropies, deriving the first order necessary optimality condition using subdifferential calculus. Via formally matched asymptotic expansions we connect our approach with previous works in the literature based on a sharp interface shape optimization description. Finally, we present several numerical results to demonstrate the advantages of our proposed approach in penalizing overhang developments.

研究の動機と目的

• 支持構造を必要とする積層造形におけるオーバーハング問題に取り組み、コストおよび後処理の複雑さを低減する。
• 最小オーバーハング角度制約を強制することで、トポロジー最適化中にオーバーハングをペナルティ化するための数学的枠組みを構築する。
• ビルド方向に応じてオーバーハング領域をペナルティ化することで、自己支持性を有する幾何形状を促進する非等方的エネルギー関数を導入する。
• 非連続的に微分可能でない非等方的性質に対し、部分微分法を用いた厳密な数学的解析を実施する。
• 形式的一致漸近展開を用いてフェーズフィールドモデルとシャープインターフェース形状最適化の間の関係を確立する。

提案手法

• 設計をフェーズフィールド変数 ϕ ∈ [0,1] で表現するフェーズフィールドアプローチを定式化し、ϕ = 1 が材料、ϕ = 0 が空隙を示す。
• オーバーハング角度が臨界閾値 ψ 未満の領域をペナルティ化する非等方的エネルギー関数 γ(ν) を導入する。ここで ν は界面における単位法線ベクトルである。
• パラメータ α を用いた凸正則化非等方的性質（6.2）を用い、ペナルティの強度および方向性を制御する。α = 1 の場合、等方的挙動に対応する。
• 非凸および非連続的に微分可能な非等方的性質に対処するため、部分微分法を用いて一次最適性条件を導出する。
• フェーズフィールドモデルとシャープインターフェース形状最適化の間の関係を検証するため、形式的一致漸近展開を実施する。
• 適応的時間刻みと質量保存制約を備えた有限要素法を用いて、得られたPDE制約付き最適化問題を数値的に解く。

実験結果

リサーチクエスチョン

• RQ1 どのように非等方的エネルギー関数を数学的に定式化することで、積層造形におけるオーバーハングをペナルティ化しつつ構造的整合性を保てるか？
• RQ2 凸非等方性と非凸非等方性の違いが、トポロジー最適化における振動境界層（ドリップ効果）の形成に与える影響は何か？
• RQ3 提案されたフェーズフィールドモデルにおける非等方的エネルギーは、古典的シャープインターフェース形状最適化モデルとどのように関係しているか？
• RQ4 非等方性パラメータ α の選択が、キャンチレバーおよびブリッジ設計における自己支持構造の幾何形状に及ぼす影響はどの程度か？
• RQ5 本手法は外部補助具を必要としない自己支持構造を生成できるか？また、それらは等方的または標準的なトポロジー最適化結果と比較してどのように異なるか？

主な発見

• 提案手法は、特に強い非等方性（α = 0.2）下で、キャンチレバー梁において内部ブリッジの傾斜を強く促進し、オーバーハングを効果的にペナルティ化した。この場合、内部構造はほぼ垂直に近づく。
• α = 0.2 の場合、内部接続部材を有しない構造に収束することが示され、自己支持構造への強い傾向が裏付けられた。
• 数値結果から、α が小さい値のとき、界面層（フェーズフィールド遷移領域）が厚くなることが確認された。これは漸近解析の予測と一致し、特定の法線方向に γ(ν) の値が高くなるためである。
• ブリッジ設計においては、非等方的ケース（α = 0.7, 0.5）が等方的ケース（α = 1）と比較して、下部に鋭い角を有する形状を示し、オーバーハング制約への適合性が向上していることが明らかになった。
• ε を 1/(32π) から 1/(64π) に減少させた結果、界面層がより薄くなることが確認され、フェーズフィールドモデルがシャープインターフェース極限に漸近収束することを裏付けた。
• 非凸非等方性ケースにおいてドリップ効果が効果的に抑制されたが、境界付近では振動的挙動が観察され、先行研究と整合的であった。