[論文レビュー] A fast and robust discrete FFT-based solver for computational homogenization

A. Finel|arXiv (Cornell University)|May 18, 2024

Composite Material Mechanics被引用数 4

ひとこと要約

本稿では、任意の弾性対比、特に無限大の対比下でも数学的安定性を保証する四面体ステンシルを用いた、計算均質化のための新規離散FFTベースのスケーラーを提示する。この手法は物理的でない振動の発生を防ぎ、完全な異方性および空間的に変化する材料関数を扱える。基本的な反復スキームを用いても、高対比下でも計算コストの増加なしに高速収束を達成する。

ABSTRACT

We propose a new discrete FFT-based method for computational homogenization of micromechanics on a regular grid that is simple, fast and robust. The discretization scheme is based on a tetrahedral stencil that displays three crucial properties. First, and most importantly, the Fourier representation of the associated Green operator is defined for any finite q-vector generated by the periodic boundary conditions and that does not belong to the Reciprocal Lattice of the discrete grids. As shown in the paper, this property guaranties that, for any elastic contrats, even infinite, mechanical equilibrium is always mathematically stable, i.e. free of any unphysical patterns, such as oscillations, ringing or checkerboarding, a property which is not shared by the original Moulinec-Suquet method \cite{moulinec1994fast,moulinec1998numerical} nor by the rotated scheme proposed by Willot \cite{willot2015fourier}. Second, the components of tensorial quantities are all defined on the same location, which permits the use of any elastic anisotropy and any spatial variation of the material fields. Third, convergence to equilibrium using the simplest iterative scheme, the "basic scheme", is fast and the number of iterates stabilizes at high contrasts, so that infinite contrast is obtained without additional computational cost.

研究の動機と目的

• 正方格子上での計算均質化に適した、安定的で高速かつ頑健な離散FFTベースのスケーラーの開発。
• 高対比、特に無限大の弾性対比下で一般的に見られる不安定性（振動、チェッカーボードパターンなど）の問題を解決すること。
• 無限大の対比下でも、逆格子に属さないすべてのqベクトルに対して数学的安定性を保証すること。
• すべてのテンソル成分を同じグリッド位置に配置することにより、任意の弾性異方性および空間的に変化する材料関数を扱えるようにすること。
• 最も単純な反復スキームを用いても、対比の変化に関係なく安定した反復回数で収束を達成すること。

提案手法

• 離散化に四面体ステンシルを採用し、周期的境界条件によって生成されるすべてのqベクトル（逆格子に属するものを除く）に対して、グリーン作用素のフーリエ表現が定義されることを保証する。
• このステンシル設計により、解の数学的安定性が保証され、リッピングやチェッカーボードパターンのような物理的でないパターンが排除される。
• 応力・ひずみなどのテンソリアル量のすべての成分を同じグリッド点に配置することで、異方性および空間的に変化する材料特性の一貫した取り扱いが可能になる。
• Lippmann-Schwinger方程式の解法に基本的な反復スキーム（固定点反復）を用い、高対比下でも収束特性が安定する。
• 無限大の対比下でも、追加の数値安定化やフィルタリングの必要がなく、本手法はその場で安定性を維持する。
• 離散的定式化により、すべての関連するqベクトルに対してグリーン作用素がフーリエ空間で適切に定義され、逆行列可能であることが保証され、あらゆる材料対比にわたる頑健性が実現される。

実験結果

リサーチクエスチョン

• RQ1 任意の弾性対比、特に無限大の対比下でも物理的でない振動を発生させず、かつ条件を満たす離散FFTベースのスケーラーを設計可能か？
• RQ2 すべてのテンソル成分を同じグリッド上に一貫して定義することで、任意の異方性および空間的に変化する材料関数を扱えるか？
• RQ3 材料対比が無限大に近づくにつれて、基本的な反復スキームの収束速度が安定的かつ高速に保たれるか？
• RQ4 高度な反復スキームやフィルタリング技術を用いなくても、頑健性と正確性を維持できるか？
• RQ5 逆格子に属さないすべてのqベクトルに対して、グリーン作用素のフーリエ表現が適切に定義されており、数学的安定性が保証されるか？

主な発見

• 本手法は、逆格子に属さないすべてのqベクトルに対してグリーン作用素のフーリエ表現が定義されることを保証することで、すべての弾性対比（無限大を含む）で数学的安定性を達成する。
• オリジナルのMoulinec-Suquet法やWillotの回転スキームとは異なり、振動、リッピング、チェッカーボードパターンなどの物理的でないパターンは完全に存在しない。
• 基本的な反復スキームを用いた収束に必要な反復回数は、対比の増加に関係なく安定しており、無限大の対比下でも増加しない。
• すべてのテンソル成分を同じグリッドに配置することで、任意の弾性異方性および空間的に変化する材料関数を扱える。
• 高度な反復技術や追加のフィルタリングを必要とせず、高速な収束を達成するため、シンプルかつ効率的である。
• 本手法は、すべてのテストされた材料対比で頑健であり、数値アーチファクトを発生させることなく、正確性を維持する。