[論文レビュー] An efficient 146-line 3D sensitivity analysis code of stress-based topology optimization written in MATLAB

Hao Deng, Praveen S. Vulimiri+1 more|arXiv (Cornell University)|Apr 2, 2021

Topology Optimization in Engineering参考文献 68被引用数 6

ひとこと要約

本稿では、p-ノルムの全応力測度とアドジョイント感度解析を用いた、効率的な3次元応力ベースのトポロジー最適化のための146行のMATLABコードを提示する。解析的微分を有限差分と比較して検証した結果、感度計算が正確に行えることが確認され、移動漸近法を用いて体積制約付きの応力最小化が複雑な3次元構造物において効果的に実行可能である。

ABSTRACT

This paper presents an efficient and compact MATLAB code for three-dimensional stress-based sensitivity analysis. The 146 lines code includes the finite element analysis and p-norm stress sensitivity analysis based on the adjoint method. The 3D sensitivity analysis for p-norm global stress measure is derived and explained in detail accompanied by corresponding MATLAB code. The correctness of the analytical sensitivity is verified by comparison with finite difference approximation. The nonlinear optimization solver is chosen as the Method of moving asymptotes (MMA). Three typical volume-constrained stress minimization problems are presented to verify the effectiveness of sensitivity analysis code. The MATLAB code presented in this paper can be extended to resolve different stress related 3D topology optimization problems. The complete program for sensitivity analysis is given in the Appendix and is intended for educational purposes only.

研究の動機と目的

• 3次元応力ベースのトポロジー最適化のための、簡潔で効率的かつ再利用可能なMATLAB実装を開発すること。
• アドジョイント法を用いてp-ノルムの全応力測度の解析的感度解析を導出し、実装すること。
• 有限差分近似との比較により、解析的感度の正確さを検証すること。
• 感度コードが3次元で体積制約付き応力最小化問題を効果的に解くことの有効性を示すこと。
• 計算力学およびトポロジー最適化分野の研究者や学生が利用できる教育的ツールを提供すること。

提案手法

• コードは3次元で線形六面体要素を用いた有限要素解析を実行する。
• 要素応力のp-ノルムを計算し、最適化のための全応力測度を形成する。
• 設計変数に関する応力測度の微分を計算するために、アドジョイントベースの感度解析を採用する。
• 感度の微分を解析的に導出し、146行のコードに直接実装する。
• 非線形最適化器として、移動漸近法（MMA）を用いて応力最小化問題を解く。
• コードは明確さと拡張性を考慮して構造化されており、他の忤力関連3次元トポロジー最適化問題への応用が可能である。

実験結果

リサーチクエスチョン

• RQ1 MATLABで、3次元応力ベースのトポロジー最適化のためのコンパクトで効率的な感度解析コードをどのように実装できるか？
• RQ2 p-ノルムの全応力測度に対する解析的感度の正確さは、有限差分近似と比較してどの程度か？
• RQ3 提案された感度コードは、体積制約付き3次元応力最小化問題を効果的に解けるか？
• RQ4 アドジョイント法は、3次元トポロジー最適化における感度の効率的計算をどのように可能にするか？
• RQ5 146行のこのコードは、応力ベースのトポロジー最適化の再利用可能な教育的・研究的ツールとしてどの程度活用可能か？

主な発見

• アドジョイント法を用いて計算された解析的感度は、有限差分近似と非常に良好に一致しており、正しさが確認された。
• 146行のMATLABコードは、1つのコンパクトな実装で3次元有限要素解析とp-ノルム応力感度計算を効果的に実行できた。
• 3つのベンチマーク問題（体積制約付き応力最小化）が効果的に解かれた。これにより、コードの強固さと効率性が示された。
• コードは教育的用途に適しており、他の忤力関連3次元トポロジー最適化問題への拡張も可能である。
• 実装は高いパフォーマンスと明確さを維持しており、計算力学分野の研究者や学生にとってアクセスしやすい。
• p-ノルムの全応力測度の使用により、3次元で滑らかで良好に条件付けられた最適化問題が実現された。