[論文レビュー] A simple Python code for computing effective properties of 2D and 3D representative volume element under periodic boundary conditions

Fan Ye, Hu Wang|arXiv (Cornell University)|Mar 11, 2017

Composite Material Mechanics参考文献 22被引用数 38

ひとこと要約

本論文では、周期的境界条件の下で2次元および3次元代表体積要素（RVE）の有効な機械的性質を計算するためのシンプルで拡張可能なPythonコードを提示する。この手法により、複合材料における効率的なマルチスケール解析が可能となり、3次元のケースにおいて数値結果と理論的・実験的ベンチマークとの強い一致が確認された。

ABSTRACT

Multiscale optimization is an attractive research field recently. For the most of optimization tools, design parameters should be updated during a close loop. Therefore, a simple Python code is programmed to obtain effective properties of Representative Volume Element (RVE) under Periodic Boundary Conditions (PBCs). It can compute the mechanical properties of a composite with a periodic structure, in two or three dimensions. The computation method is based on the Asymptotic Homogenization Theory (AHT). With simple modifications, the basic Python code may be extended to the computation of the effective properties of more complex microstructure. Moreover, the code provides a convenient platform upon the optimization for the material and geometric composite design. The user may experiment with various algorithms and tackle a wide range of problems. To verify the effectiveness and reliability of the code, a three-dimensional case is employed to illuminate the code. Finally numerical results obtained by the code agree well with the available theoretical and experimental results

研究の動機と目的

• 2次元および3次元複合材料のRVEの有効性を計算するためのユーザーフレンドリーで拡張可能なPython実装を開発すること。
• 周期的境界条件と漸近的均質化を統合することで、効率的なマルチスケール最適化を可能にすること。
• さまざまなアルゴリズムや材料の微細構造をテストするための柔軟な計算プラットフォームを提供すること。
• 3次元ケースにおける理論的および実験的結果と比較することで、コードの正確性を検証すること。

提案手法

• この手法は、不均質な微細構造を有効なマクロ特性に均質化することができる漸近的均質化理論（AHT）に基づいている。
• 無限大の周期的複合材料を模擬するために周期的境界条件（PBC）が適用され、代表的かつ一貫性のある境界応答が保証される。
• 構造的グリッド上での標準的なFEM離散化を用いて、平衡方程式の弱形式を有限要素法で解く。
• 有効な剛性テンソルは、ユニットセルのシミュレーションから得られた応力およびひずみ場の体積平均を用いて計算される。
• 実装はモジュラーであり、単純な周期的配置を超える複雑な微細構造に対しても容易に修正可能である。
• コードは純粋なPythonで記述されており、最小限の依存関係を有するため、最適化ワークフローへの統合に向けたポータビリティと拡張性が向上する。

実験結果

リサーチクエスチョン

• RQ1 2次元および3次元RVEの周期的境界条件の下での有効性を計算するためのシンプルで効率的なPython実装は、どのように設計できるか？
• RQ2 理論的および実験的データと比較した場合、コードは有効な機械的性質をどの程度正確に予測できるか？
• RQ3 コードは、最小限のコード修正でより複雑な微細構造をモデル化するために容易に拡張可能か？
• RQ4 PythonベースのFEMフレームワークにおける漸近的均質化理論の使用は、マルチスケール材料設計および最適化をどのように支援するか？

主な発見

• 3次元ケースにおける数値結果は、利用可能な理論的および実験的データと優れた一致を示しており、コードの信頼性が確認された。
• 実装は、2次元および3次元の両領域において、周期的複合材料微細構造の有効剛性テンソルを正常に計算できた。
• コードは強固でスケーラブルであり、ユーザーがさまざまな微細構造的配置および境界条件で実験を行うことが可能である。
• モジュラー設計により、より複雑な材料系や高度な最適化アルゴリズムへの容易な拡張が可能である。
• 周期的境界条件の使用により、無限大の周期的複合材料に対して物理的に一貫性があり代表的な均質化特性が保証された。
• コードのシンプルさと依存関係の少なさのおかげで、大規模なマルチスケールシミュレーションおよび設計パイプラインへの統合に適している。