[論文レビュー] Meshless Methods for Large Deformation Elastodynamics

Matthias Röthlin, Hagen Klippel+1 more|arXiv (Cornell University)|Jul 3, 2018

Fluid Dynamics Simulations and Interactions参考文献 78被引用数 6

ひとこと要約

本論文は、大変形弾性動力学におけるメッシュレス法の包括的ベンチマークを提示し、特にSPH、RKPM、CSPM、EFGの定式化に焦点を当てる。一貫性のあるトータル・ラグランジュアンフレームワークと一貫性のある微分を導入し、1000%に達する極限的大変形下での安定性、精度、耐障害性の再現可能な比較を可能にするオープンソースコードを提供する。

ABSTRACT

Meshless methods are a promising candidate to reliably simulate materials undergoing large deformations. Unlike mesh based methods like the FEM, meshless methods are not limited in the amount of deformation they can reproduce since there are no mesh regularity constraints to consider. However, other numerical issues like zero energy modes, the tensile instability and disorder of the discretization points due to the deformation may impose limits on the deformations possible. It is thus worthwhile to benchmark a wide array of these methods since a proper review to this end has been missing from the literature so far. In the interest of reproducibility, the complete source code of all methods considered is made public.

研究の動機と目的

• メッシュレス法における大変形弾性動力学の分野において、体系的な比較と再現可能なベンチマークが不足している問題を解決すること。
• 一貫したトータル・ラグランジュアン定式化を用いて、1000%に達する極限的大変形下でのSPH、RKPM、CSPM、EFGの安定性、精度、耐障害性を評価すること。
• すべての手法の一貫性があり、オープンソースで利用可能な実装を開発・公開し、再現可能性を確保するとともに、今後の研究を促進すること。
• 一貫性のある微分定式化および安定化技術が、大変形領域における数値的性能に与える影響を調査すること。
• 今後のメッシュレス法の開発および検証のための基準フレームワークを提供すること。

提案手法

• 大変形における材料点の分離問題を回避するため、トータル・ラグランジュアンフレームワークを用いて弾性動力学方程式の弱形式を定式化する。
• 再現核粒子法（RKPM）を用いて一貫性のある微分を適用し、勾配近似における線形完全性と精度の向上を実現する。
• ひずみおよび応力計算に適した一貫性のある形状関数の微分を導出し、∇ᵤⱼε = 2·Fᵀ·dⱼを介して正確な力の計算を可能にする。
• ∇ᵤᵢεおよびσを用いた一貫性のある力の計算を実装し、fⱼ = -2·vᵢ·F·σ·dⱼを導出することで、運動量およびエネルギー保存を確保する。
• 時計回りのモード（hourglass modes）および引張不安定性を軽減するため、節点統合に安定化項を適用する。
• 大変形下での粒子の分布品質の維持とクラスタリングや無秩序の防止のため、rSPHに類似した粒子正規化を適用する。

実験結果

リサーチクエスチョン

• RQ1 SPH、RKPM、CSPM、EFGといった異なるメッシュレス法は、1000%を超える大変形下において、精度および安定性の面でどのように性能を発揮するか？
• RQ2 一貫性のあるトータル・ラグランジュアン定式化と一貫性のある微分を用いることで、数値的安定性および収束性にどのような影響を与えるか？
• RQ3 人工応力や粒子正規化といった安定化技術は、弾性動力学におけるメッシュレス法の耐障害性にどのような影響を与えるか？
• RQ4 統一された実装フレームワークは、複数のメッシュレス法間での公平で再現可能なベンチマークを可能にするか？
• RQ5 既存のメッシュレス法には、固体動力学における極限的大変形、接触、トポロジーの変化を処理する上でどのような限界があるか？

主な発見

• 一貫性のあるトータル・ラグランジュアン定式化と一貫性のある微分を用いることで、標準的なSPH定式化と比較して、精度および安定性が顕著に向上する。
• RKPMおよびCSPMは、1000%に達する大変形を処理する上で、歪みの固定（locking）が低減され、収束性が優れている。
• RKPMによる一貫性のある微分の導入により、不適切なゼロエネルギーモードが排除され、動的シミュレーションにおけるエネルギー保存が向上する。
• rSPHに類似した粒子正規化は、極限的大変形下でも粒子の分布品質を効果的に維持し、数値的不安定性を防止する。
• GitHub上のオープンソース実装により、完全な再現可能性が達成され、今後のメッシュレス法開発の基準ベンチマークとして機能する。
• EFGおよびTLEDに基づく手法は、ソフト tissue や衝撃シミュレーションにおいて有望であるが、引張不安定性を回避するための注意深い安定化が求められる。