A 2D discrete model with a phase-field regularisation applied to dynamic crack propagation

ガラスや硬化鋼のような脆い材料が壊れるとき、それは単に真っ二つに割れるのではありません。代わりに、亀裂が固体の中を猛スピードで駆け抜け、しばしば複雑な破断の網へと枝分かれし、物体を多くの破片へと粉砕します。これらの亀裂がどこを通るのか、そしてどれほどの速さで移動するのかを正確に予測することは、工学における最も困難な課題の一つです。もし、この破壊をシミュレートするために使用されるモデルに欠陥があれば、橋が実際には壊滅的に崩落するにもかかわらず、安全な方向に壊れると示唆してしまうかもしれません。数十年にわたり、科学者たちは材料の原子構造を模倣するコンピュータモデルの構築を試みてきました。それらは、小さな点同士がバネでつながれた格子として材料を扱うものです。これらのモデルは強力ではありますが、ある特有の欠陥を抱えています。それは、亀裂が進む経路が、材料自体の物理特性ではなく、コンピュータのグリッド（格子）がどのように回転しているかに依存してしまうという点です。それはまるで、本物のガラスに生じた亀裂が、エンジニアがコンピュータの画面を回転させたという理由だけで、突然進行方向を変えてしまうかのようです。

この問題を解決するために、スペインのネブリア大学の研究者であるマティアス・ブラウンは、動的破壊をシミュレートする新しい手法を開発しました。彼のアプローチは、従来の点の格子と「フェーズフィールド法」を組み合わせたものです。フェーズフィールド法とは、亀裂を点から点へと飛び跳ねる鋭い線としてではなく、材料の中に広がる滑らかで曖昧な損傷領域として扱う手法です。これら二つの考え方を融合させることで、ブラウンは、基礎となるコンピュータのグリッドの向きに関わらず、あらゆる方向に亀裂が進むことができるモデルを作り上げました。彼はこの新手法を、高速衝撃を受けた鋼鉄板、引き裂かれて破断する切り欠きのあるプレート、そして内部から爆発する厚い金属シリンダーという、性質の異なる3つのシナリオでテストしました。あらゆるケースにおいて、この新モデルは実世界の実験結果と一致し、コンピュータのグリッドによる人工的なバイアスを排除し、材料が真にどのように破壊されるかについて、より明確で信頼性の高い視点を提供しました。

ブラウンが取り組んだ問題の核心は、コンピュータが原子という目に見えない世界をどのように扱うかという点にあります。彼の以前の研究では、離散モデルを用いており、そこでは材料は線でつながれた三角形の格子（ドットのパターン）として表現されていました。材料にストレスがかかると、ドット間の接続が切れ、それが亀裂をシミュレートします。しかし、ドットが固定されたパターンで配置されているため、亀裂はその格子の線に従うことを強制されました。そのため、グリッドを回転させると、材料自体は変わっていないにもかかわらず、亀裂の経路が変わってしまうのです。これは重大な問題です。なぜなら、現実の亀裂は、コンピュータのメッシュの恣意的な角度ではなく、材料の特性によって決定される「抵抗の少ない経路」に従うべきだからです。

ブラウンの解決策は、「フェーズフィールド正則化」を導入することでした。特定のドット間の接続を壊すように強制する代わりに、彼は各ドットに対して、その地点がどの程度損傷しているかを表す値を割り当てました。この値は、材料が完全に無傷であることを示すゼロから、完全に破壊されていることを示す一までの範囲をとります。これにより、亀裂は値が徐々に増加していく領域となり、鋭い線ではなく滑らかな遷移を生み出します。これにより、亀裂はあらゆる方向に曲がったり枝分かれしたりすることが可能になり、実質的にグリッドの硬直した構造を無視できるようになりました。決定的なのは、ブラウンがこのシステムを、材料の2つの基本的な特性、すなわち「引張強度（亀裂を開始させるのに必要な力）」と「破壊エネルギー（亀裂を成長させ続けるために必要なエネルギー）」を用いて校正したことです。シミュレーションをこれらの実在する物理定数に結びつけることで、モデルの結果が数学ではなく、材料に依存することを保証しました。

この新手法が機能することを証明するために、ブラウンは3つの厳格なテストを行いました。最初のテストは、古典的な実験である「カルトホフ・ウィンクラー・テスト」です。これは、あらかじめ切り欠きを入れた鋼鉄板に、発射体が高速で衝突するものです。このシナリオでは、亀裂が切り欠きから特定の角度で飛び出します。従来のモデルでは、グリッドの回転に応じて誤った角度を予測することがよくありました。しかし、ブラウンのシミュレーションでは、亀裂は一貫して約70度の角度で進行し、実験の観察結果と完璧に一致しました。この結果は、グリッドがプレートに対して整列していても、あるいは15度や30度回転していても、同様に成立しました。亀裂の速度も安定しており、グリッドの向きに依存することなく、材料中の音速の約70パーセントで進行しました。

2番目のテストは、一定の力で引っ張られる単一エッジ切り欠きプレートを用いたもので、亀裂がゆっくりとどのように成長するかを研究するためのセットアップです。ここでの目的は、亀裂が切り欠きから反対側の端に向かって直線的に進むかどうかを確認することでした。従来の格子モデルでは、グリッドの配置に基づいて、亀裂がジグザグに進んだり逸れたりすることがよくありました。しかし、ブラウンのモデルは、あらゆる方向のグリッドに対して直線的な亀裂を生み出しました。プレートを破壊するために必要な力も、グリッドを回転させても2パーセント未満の変動にとどまり、一貫していました。これは、モデルがコンピュータメッシュの人工的な影響を受けることなく、材料の経路と強度を正確に予測できることを示しています。

最後にして最も複雑なテストは、内部爆発にさらされた厚いシリンダーの断片化です。このシナリオでは、単一の亀裂が切り欠きによって導かれるのではなく、多数の亀裂が自発的に発生し、外側に向かって猛烈に広がっていく必要があります。これは、グリッドの対称性が破片のパターンを決定してしまう可能性のある、混沌としたイベントです。ブラウンがシリンダーの軸を中心にグリッドを配置してシミュレーションを実行したとき、亀裂は対称なパターンを形成しました。しかし、グリンドを少しずらして中心から外れた状態にしたところ、亀裂は不規則に分布し、破片の数も変化しました。これは、モデルが単なるグリッドの対称性ではなく、物理的なセットアップに対して敏感であることを裏付けています。シミュレーションは12個から19個の破片を予測しましたが、これは他の高度な手法の範囲内ではありますが、やや低い数値でした。これはおそらく、モデルの滑らかな損傷領域が、近すぎる亀裂の形成を防いだためと考えられます。

これらのテストを通じて、研究者たちはシステムのエネルギーを注意深く追跡しました。彼らは、材料が壊れる際に失われるエネルギーと、新しい亀裂面を作るために必要なエネルギーを比較しました。いくつかのケースでは、損傷モデルによって計算されたエネルギーと、実際に切断された結合によって放出されたエネルギーとの間に、わずかな不一致が見られました。これは、ハイブリッドなアプローチにおける既知の限界であり、亀裂を駆動する力が、壊れた結合のエネルギーとは異なる方法で計算されることに起因します。しかし、機械的なエネルギーバランスは0.5パーセント以内の精度で正確であり、シミュレーションが物理的に妥当であることを裏付けました。これらの結果は、モデルが完璧ではないとしても、動的破壊の理解における重要な前進であることを示唆しています。

この研究の意義は、単なる優れたコンピュータグラフィックスの向上にとどまりません。グリッドの向きへの依存性を排除することで、エンジニアは、衝突テスト、構造的破壊、または材料試験のシミュレーションが、材料の真の物理現象を反映していると信頼できるようになります。このモデルは、原子の離散的な性質と、固体の連続的な挙動との間の溝を埋めることに成功しました。適切な数学的ツールを用いれば、コンピュータのグリッドに結果を左右されることなく、破壊という混沌とした激しいプロセスをシミュレートすることが可能であることを、この研究は示しています。結論として、このアプローチは、極限条件下での材料の破壊を予測するための堅牢な方法を提供し、より安全で弾力性のある構造物を設計するための明確な道筋を示すものです。