Dynamic brittle fracture using Lip-field approach in an explicit dynamics context

技術要約：明示的動力学におけるリップフィールド・アプローチを用いた動的脆性破壊

問題提起
本論文は、衝撃や爆発などの高歪速度荷重を受ける材料における、動的な脆性破壊の数値シミュレーションについて扱っている。古典的な破壊力学（例：グリフィスのモデル）は、亀裂の発生や分岐の扱いに苦慮するが、現代的なアプローチとして、凝集ゾーンモデル（CZM）や拡散損傷モデルが開発されてきた。しかし、局所的な拡散損傷モデルは病理的なメッシュ感度（mesh sensitivity）に陥るため、長さスケールを導入しメッシュに依存しない結果を得るための正則化が必要となる。フェーズフィールド法や厚いレベルセット（TLS）法といった既存の正則化手法は、多くの場合、計算コストの高い非局所的なエネルギー汎関数を伴う。ここで取り組む具体的な課題は、リップフィールド・アプローチ（損傷場に対してリプシッツ制約を課すことで長さスケールを導入する変分法）を、これまでの準静的および一次元動的な適用から、二次元動的破壊および明示的な時間積分フレームワークへと拡張することである。

手法
著者らは、増分ポテンシャルの最小化に基づく動的損傷力学の変分定式化を提案している。手法の核となる内容は以下の通りである：

1. 変分定式化： 力学的問題は、運動エネルギー、自由エネルギー、および散逸ポテンシャルを含む増分ポテンシャル $F$ の最小化によって定義される。状態変数は変位 ($u$) とスカラー損傷変数 ($d$) である。構成モデルは、特定のエネルギー低減関数 $g(d)$ および軟化関数 $h(d)$ を持つ一般化標準材料を利用する。
2. リップフィールド正則化： エネルギー汎関数に勾配項を加えるフェーズフィールド法とは異なり、リップフィールド・アプローチはエネルギー汎関数を局所的（非正則化）に保つ。その代わりに、損傷場に対してリプシッツ制約を課し、損傷勾配が臨界値 ($1/l$) によって制限されるようにする。これにより、損傷場を特定のリプシッツ空間 $\mathcal{L}$ に制限し、局所的なエネルギー密度を変更することなく、実質的に長さスケール $l$ を導入する。
3. 明示的交互計算スキーム（Explicit Staggered Scheme）： 著者らは、明示的なニューマーク時間積分スキーム（中央差分法）を採用している。変位場と損傷場を逐次的に解くための交互計算（staggered approach）が用いられる。各タイムステップにおいて、まず変位が更新され、続いて損傷が更新される。
4. 効率的な損傷ソルバー： アルゴリズムにおける主要な貢献は、損傷計算を最適化するための境界推定（先行研究に基づく）の使用である。プロセスは以下の通りである：
   • 純粋な局所最小化問題を解き、局所的な損傷場 ($d_{loc}$) を求める。
   • $d_{loc}$ をリプシッツ制約に基づいた上限および下限へと投影する。
   • これらの境界が異なる領域を特定する。
   • 非局所的（リプシッツ制約付き）な最小化問題を、境界が一致しない領域のみにおいて解く。境界が一致する領域では、局所解が採用される。これにより、全領域に対して非局所問題を解く場合に比べて計算コストが大幅に削減される。
5. 離散化： 空間離散化には、変位場に対して線形三角形要素を使用する。損傷場を保持・制約するために、別途「リップメッシュ」（要素重心から構成される）が使用され、これによりリプシッツ制約の数を最小限に抑えている。

主な貢献

• 2D動力学への拡張： 本論文は、1次元の動的破砕および準静的なケースで検証されていたリップフィールド・アプローチを、二次元の動的破壊問題へと成功裏に拡張した。
• 明示的動力学への統合： 高速衝突問題において計算効率が高い、明示的時間積分スキームを用いたリップフィールド・アプローチの具体的な実装方法を提示している。
• 境界による計算効率化： 境界推定を利用することで、高コストな非局所的最小化を小さなサブドメインに限定できることを示しており、リップフィールド・アプローチを動的シミュレーションにおいて計算競争力のあるものにしている。
• 変分的な区別： 本研究は、リップフィールド法が局所的なエネルギー汎関数を持ち、非局所的な制約を課すのに対し、フェーズフィールド法は非局所的なエネルギー汎関数を用いるという、両者の違いを明確にしている。

シミュレーション結果
提案されたモデルは、以下の二次元数値例を通じて検証された：

1. 単一エッジ切り欠き引張試験（SENT）： モデルは動的な亀裂進展および分岐を正常にシミュレートした。結果は、ポテンシャルおよび散逸エネルギーにおいてメッシュ収束を示した。亀裂の分岐現象は、アドホックな基準なしに捉えられ、亀裂先端速度は、分岐イベント中の過渡的なスパイクを除いて、概ねレイリー波速度以下に保たれた。結果は、文献に見られるフェーズフィールド法やTLSシミュレーションと同等であった。
2. 荷重強度の影響： 印加されるトラクションを変化させることで、モデルは複数の亀裂分岐を捉える能力を示し、荷重が高い場合には、よりノッチ先端に近い位置で分岐が発生することが示された。
3. カルトフ・ウィンクラー（Kalthoff-Winkler）試験： モデルは、動的なせん断荷重下での亀裂発生と進展をシミュレートした。
   • 低衝突速度（33 m/s）において、実験的な観察値（約70°）と一致する約62.4°のモードI亀裂進展角を予測した。
   • 高衝突速度（100 m/s）において、モデルは連続的な亀裂分岐と複数の亀裂を捉えた。
   • 限界事項： 著者らは、実験では高速域で延性破壊（約-10°のせん断破壊）への遷移が見られるものの、現在の弾性損傷モデルには塑性が欠けているため、この延性挙能を捉えることはできないと指摘している。

意義と主張
本論文は、リップフィールド・アプローチが、必要な長さスケールを導入しつつ、問題の数学的な妥当性を維持する堅牢な変分フレームワークを提供すると主張している。その主な意義は、境界ベースのソルバーを通じて達成される計算効率にあり、これにより非局所方程式をグローバルに解くことを回避している。著者らは、この手法がアドホックな基準なしに、亀裂の発生、分岐、および合体といった複雑な破壊現象を効果的に処理できると断言している。

論文は、リップフィールド・アプローチが脆性破壊と破砕には有効であるが、弾性損傷挙動に限定されていると控えめに結論付けている。高速歪速度下での延性破壊をシミュレートするには、弾塑性損傷モデルへの拡張が必要であるが、それは本研究の範囲外であるとしている。本手法は、特に計算効率が重要な明示的動力学において、フェーズフィールド法やTLS法の有力な代替案として提示されている。