[論文レビュー] A 55-line code for large-scale parallel topology optimization in 2D and 3D

Abhinav Gupta, Rajib Chowdhury+2 more|arXiv (Cornell University)|Dec 15, 2020

Topology Optimization in Engineering参考文献 71被引用数 9

ひとこと要約

本稿では、FEniCS有限要素フレームワークとPETScを用いた並列線形代数により、大規模な2次元および3次元トポロジー最適化を効率的に行う55行のPythonコードを提示する。SIMPに基づく最適化を、ユーグリッド距離行列を用いたベクトル化フィルタリングにより実現し、複雑なメッシュにおいて高い性能とスケーラビリティを達成する。複数の荷重ケース、境界条件、パassing要素をサポートしており、教育的および大規模な研究用途に最適である。

ABSTRACT

This paper presents a 55-line code written in python for 2D and 3D topology optimization (TO) based on the open-source finite element computing software (FEniCS), equipped with various finite element tools and solvers. PETSc is used as the linear algebra back-end, which results in significantly less computational time than standard python libraries. The code is designed based on the popular solid isotropic material with penalization (SIMP) methodology. Extensions to multiple load cases, different boundary conditions, and incorporation of passive elements are also presented. Thus, this implementation is the most compact implementation of SIMP based topology optimization for 3D as well as 2D problems. Utilizing the concept of Euclidean distance matrix to vectorize the computation of the weight matrix for the filter, we have achieved a substantial reduction in the computational time and have also made it possible for the code to work with complex ground structure configurations. We have also presented the code's extension to large-scale topology optimization problems with support for parallel computations on complex structural configuration, which could help students and researchers explore novel insights into the TO problem with dense meshes. Appendix-A contains the complete code, and the website: \url{https://github.com/iitrabhi/topo-fenics} also contains the complete code.

研究の動機と目的

• 学生や研究者にとってアクセス可能でありながら、大規模な問題を扱えるようにする、コンactで読みやすく拡張可能なトポロジー最適化コードの開発。
• 商業用FEAソフトウェアに依存せずに、オープンソースツールを用いて効率的な大規模2次元および3次元トポロジー最適化を可能にすること。
• 特に複雑なグランド構造に対して、計算時間の短縮を図るために、フィルタ計算をベクトル化されたユーグリッド距離行列を用いて行うこと。
• HPCシステムにおけるスケーラビリティを実現するため、PETScによる並列化を可能にすること。
• 最小限のコードベース内で、古典的なSIMP法を複数の荷重ケース、さまざまな境界条件、およびパassing要素へと拡張すること。

提案手法

• FEniCSを用いて有限要素のアセンブリと解法を実装し、SOLID ISOTROPIC MATERIAL WITH PENALIZATION (SIMP) 法を実装する。
• 線形代数のバックエンドとしてPETScを用い、行列演算の高速化と並列計算の実現を図る。
• フィルタ重みの計算を効率的に行うために、ベクトル化されたユーグリッド距離行列を計算する。
• フィルタ操作は行列-ベクトル乗算として実装され、感度 = (D @ (密度 * 感度)) / (密度 * 距離合計) で表され、Dは距離行列である。
• 体積制約を満たすために、ラインサーチを用いた最適性基準法を密度更新に用いる。
• メッシュと荷重マーカーのモジュラーな設定により、複数の荷重ケース、ディリクレ境界条件、およびパassing要素をサポートする。

実験結果

リサーチクエスチョン

• RQ1 55行という極めてコンactなコードで、オープンソースツールを用いてスケーラブルで高性能な2次元および3次元トポロジー最適化を実現できるか？
• RQ2 従来のループベースの方法と比較して、ベクトル化されたユーグリッド距離行列による計算は、フィルタの効率をどのように向上させるか？
• RQ3 最小限のコードベースが、複雑な構造的配置、複数の荷重ケース、並列実行をどの程度サポートできるか？
• RQ4 このようなコンactな実装が、HPCシステム上で大規模メッシュにスケーリングする際、正確性と収束性を維持できるか？
• RQ5 FEniCSとPETScの統合は、教育的明快さと生産レベルのパフォーマンスを両立できるか？

主な発見

• 55行のコードは、複数の荷重ケース、境界条件、パassing要素をサポートする完全な2次元および3次元SIMPベースのトポロジー最適化を実装している。
• ベクトル化されたユーグリッド距離行列の使用により、フィルタ計算時間は顕著に短縮され、複雑なグランド構造の効率的処理が可能になった。
• PETScの並列線形代数バックエンドのおかげで、大規模メッシュにおける計算が顕著に高速化され、スケーラブルな計算が実現した。
• コードは高い読みやすさとモジュラリティを維持しており、教育的用途や迅速なプロトタイピングに適している。
• FEniCSのプラットフォーム非依存性とオープンソーススタックのおかげで、Windows、Linux、macOSの各環境でポータブルに動作する。
• 体積制約が保持されたまま、目的関数と密度の反復履歴から示されるように、最適設計に収束した。