[論文レビュー] Remeshing-Free Graph-Based Finite Element Method for Ductile and Brittle Fracture

Avirup Mandal, Parag Chaudhuri+1 more|arXiv (Cornell University)|Mar 27, 2021

Structural Response to Dynamic Loads被引用数 1

ひとこと要約

本論文は、四面体ボリュメトリックメッシュから導かれるグラフ上で破壊をモデル化することにより、延性破壊および破壊破壊の両方を再メッシュ不要でシミュレートするグラフベースの有限要素法を提案する。超弾性ひずみエネルギーをエッジ長のみの関数として表現し、破壊変数を用いて亀列の進行を追跡することで、システム行列のサイズを固定したままにし、再メッシュなしで高いスケーラビリティと安定性を実現する。

ABSTRACT

Fracture produces new mesh fragments that introduce additional degrees of freedom in the system dynamics. Existing finite element method (FEM) based solutions suffer from an explosion in computational cost as the system matrix size increases. We solve this problem by presenting a graph-based FEM model for fracture simulation that is remeshing-free and easily scales to high-resolution meshes. Our algorithm models fracture on the graph induced in a volumetric mesh with tetrahedral elements. We relabel the edges of the graph using a computed damage variable to initialize and propagate fracture. We prove that non-linear, hyper-elastic strain energy is expressible entirely in terms of the edge lengths of the induced graph. This allows us to reformulate the system dynamics for the relabeled graph without changing the size of system dynamics matrix and thus prevents the computational cost from blowing up. The fractured surface has to be reconstructed explicitly only for visualization purposes. We simulate standard laboratory experiments from structural mechanics and compare the results with corresponding real-world experiments. We fracture objects made of a variety of brittle and ductile materials, and show that our technique offers stability and speed that is unmatched in current literature.

研究の動機と目的

• 亀列の形成に伴い自由度が増加することによる従来のFEMベースの破壊シミュレーションにおけるスケーラビリティおよび計算コストの問題を解決すること。
• 破壊進行中に再メッシュの不安定性および高コストを抱える既存のFEMおよびXFEM手法の課題を克服すること。
• 元の計算メッシュを維持することで、再メッシュなしに安定的かつ高解像度の破壊シミュレーションを可能にすること。
• 非局所破壊基準を用いることで、1つの統一フレームワーク内で破壊的および延性的破壊を両立させること。
• エッジベースの破壊変数を用いたグラフ上でのシステムダイナミクスの再定式化により、システム行列サイズを一定に保ちながら効率的なシミュレーションを実現すること。

提案手法

• 四面体ボリュメトリックメッシュからグラフを構築し、メッシュの頂点をグラフのノード、メッシュのエッジをグラフのエッジとする。
• 導かれたグラフにおけるエッジ長のみに依存する形で、システムの超弾性ひずみエネルギーを完全に再定式化し、ノード位置に依存しないようにする。
• 材料の劣化を表すエッジベースの破壊変数を導入し、非局所破壊基準に基づいて更新することで亀列の拡散を制御する。
• グラフベースのひずみエネルギー定式化を用いてノード力の計算とシステムダイナミクスの更新を実行し、システム行列サイズを変更しない。
• 亀列が進行しても元の計算メッシュを維持し、再メッシュを一切行わない。
• 破壊表面の再構築は可視化のための処理に限定し、計算メッシュは変更されない。

実験結果

リサーチクエスチョン

• RQ1 グラフベースのFEM定式化は、亀列進行中でもシステム行列サイズを維持できるか。再メッシュおよび関連計算コストを回避できるか。
• RQ2 3次元ボリュメトリックFEMにおける超弾性ひずみエネルギーは、グラフ表現におけるエッジ長のみで完全に表現可能か。
• RQ3 非局所破壊基準をどのようにグラフベースモデルに統合することで、現実的な亀列の拡散および進行を実現できるか。
• RQ4 提案手法は、高解像度メッシュ上でも破壊的および延性的破壊の両方を安定的かつスケーラブルにシミュレートできるか。
• RQ5 再メッシュ不要なアプローチは、既存のFEMおよびXFEM技術と比較して計算効率および数値的安定性の面で優れているか。

主な発見

• 破壊シミュレーションの全過程においてシステム行列サイズが一定のまま維持され、高解像度メッシュへの効率的スケーリングが可能である。
• 再メッシュや再メッシュのオーバーヘッドなしに、破壊的および延性的破壊の両方を安定的かつ正確にシミュレートできる。
• グラフベースの定式化により、エッジ長のみを用いた超弾性ひずみエネルギーの正確な表現が可能となり、一貫性のあるシステムダイナミクスが実現される。
• 本手法は構造力学分野の標準的実験を成功裏に再現でき、実世界の物理的挙動と密接に一致する結果を示した。
• 特に複数または複雑な亀列パターンを伴うシナリオにおいて、既存の文献と比較して類似しない速度と安定性を示した。
• 破壊表面の再構築は計算と分離されており、シミュレーション性能に影響を与えることなく、可視化時のみの処理で十分である。