Strength-degradation phase-field regularization of cohesive fracture: the antiplane case

ガラスの破片、粘土の塊、あるいは金属の梁がどのように壊れるかを予測しようとしていると想像してみてください。長い間、科学者たちはこの問題に対して、二つの全く異なる見方をしてきました。一方で、「脆性（ぜいせい）」の視点があります。これは、乾燥した小枝のように材料が突然パキッと折れるという考え方で、新しい亀裂の表面を作るためにどれだけのエネルギーが必要かのみが重要となります。他方には、「延性（えんせい）」の視点があります。これは、 taffy（タフィー）のように、壊れる前に材料がぐにゃりと曲がったり伸びたりするという考え方で、圧力がかかったときに材料がどれだけ降伏（変形）するかが鍵となります。

厄介なのは、現実世界の材料はしばなくの両方の性質を持っていることです。金属の破片は、最後にパキッと折れる前に、少し伸びたり（塑性）、形が変わったりすることがあります。数十年間、コンピュータモデルはこの混合状態を扱うのに苦労してきました。普及している「フェーズフィールド」モデル（亀裂がピクセル間をジャンプするのではなく、滑らかに成長するようにするデジタルシミュレーションの一種）には、大きな欠陥がありました。それは、材料の強度が、コンピュータ画面上で亀裂がどれくらい「ぼやけて（fuzzy）」見えるかに直接結びついてしまっていたことです。もし強い材料をシミュレートしたいなら、コンピュータ上の「ぼやけ」を極限まで小さくしなければならず、それにはスーパーコンピュータが必要でした。計算時間を節約するために「ぼやけ」を大きくすると、材料が魔法のように弱くなってしまうのです。これは、「強さ」「速さ」「現実的な亀裂の形状」のうち、二つを選ぶことはできても、三つすべてを同時に手に入れることはできないという、もどかしい「二択」のゲームでした。

この論文は、この結び目をようやく解き明かす、材料の破壊をモデル化する新しい方法を紹介しています。著者であるブレイズ・ブルダンとクラード・モーリニは、材料の強度が、その剛性や靭性（じんせい）と同じように、独立した別々の要素であるとするモデルを提案しました。彼らは、このアイデアを「アンチプレーンせん断」（トランプの束を横にスライドさせるようなもの）という簡略化された設定でテストし、この新しいモデルが非常にうまく機能することを見出しました。このモデルは、金属が伸び縮みする（塑性）ことから、きれいな突然の破断（脆性破壊）に至るまで、ゲームのルールを変えることなくシミュレートできます。最も重要な点は、コンピュータシミュレーションにおける「ぼやけ」が、数学を成立させるための数値的なツールに過ぎず、材料の強さを変えてしまう物理的な特性ではないことを示したことです。これにより、科学者たちは、グリッドサイズを調整しただけで結果が変わってしまう心配をすることなく、標準的なコンピュータで複雑な破壊挙動をシミュレートできるようになりました。

新しい「強度優先」のアプローチ

この論文の核心は、材料がどのように壊れるかをシミュレートするための新しい数学的なレシピです。著者らはこれを「強度劣化フェーズフィールド正則化（strength-degradation phase-field regularization）」と呼んでいます。これは非常に長い名前なので、簡単な比喩で説明しましょう。

材料を、手を繋いでいる人々の群衆だと想像してください。古いモデルでは、群衆がストレスを感じ始めると、人々の腕（剛性）がどんどん弱くなり、最終的に手を離してしまうというものでした。問題は、彼らが手を離すタイミングが、あなたがどれほど近くで見ているか（「正則化長」）に完全に依存していたことです。近くで見れば長く手を繋ぎ続け、遠くから見れば早く手を離してしまうのです。

この新しいモデルでは、人々は腕を弱めることはしません。その代わりに、彼らには「最大握力」があります。引力がその強さの下にある限り、彼らはしっかりと保持します。一度その限界に達すると、握力が低下し始め、最終的に手を離します。決定的なのは、この「最大握力」は人々の固定された特性であり、あなたの観察方法による副作用ではないということです。これにより、このモデルは、強くてあまり伸びない材料や、壊れる前にたくさん伸びる材料を、すべて同じ枠組みの中で扱うことができます。

破壊の三段階

著者らは、モデルがさまざまな条件下でどのように振る舞うかを調べるためにシミュレーションを行いました。彼らは、単一のモデルが、対象のサイズや材料の特性に応じて、自然に三つの異なる破壊段階を遷移できることを発見しました。

1. 「ぐにゃり」としたフェーズ（小規模降伏）： 小さな領域で亀裂が発生するとき、その先端付近の材料は単にパキッと折れるわけではありません。材料は伸び縮みして変形し、「プロセスゾーン」と呼ばれる降伏領域を作り出します。モデルはこの現象を完璧に捉え、古典的な塑性理論と一致する変形領域を示しています。
2. 「凝集（ぎょうしゅう）」フェーズ： 亀裂が開いていくにつれ、材料は即座に壊れるわけではありません。代わりに、ゆっくりと剥がれていく粘着テープのように、「凝集力」を持って踏みとどまります。隙間が広がるにつれて、引き離すために必要な力は減少していきます。これが「凝集亀裂」の挙動であり、現実の材料がどのように破壊するかを理解する上で不可欠な要素です。
3. 「脆性」フェーズ： 最終的に、隙間が十分に広くなると、材料は完全に諦めてパキッと折れます。亀裂は、最も単純な脆性破壊モデルと同様に、鋭くきれいな破断となります。

この研究の素晴らしさは、これら三つの段階を記述するために異なる数学的公式を切り替える必要がない点にあります。すべては一つの連続したプロセスです。モデルは、人間のプログラマーの選択ではなく、物理学に基づいて、ある状態から次の状態へと進化していくのです。

「サーフィン」実験と塑性ウェイク

論文の中で最も興味深い部分の一つは、著者らが「サーフィン」と呼ぶシミュレーションです。すでに亀裂が存在しており、それを一定の速度で材料の中をドラッグしていく様子を想像してください。まるで波に乗るサーファーのように。著者らはこのシミュレーションを二つの方法で行いました。一つは、材料の変形が可逆的（ゴムバンドのような）な場合、もう一つは、変形が不可逆的（形が変わったままの粘土のような）な場合です。

可逆的なケースでは、材料を壊すのに必要なエネルギーは理論の予測と正確に一致しました。しかし、不可逆的なケースでは、興味深いことが起こりました。亀裂が移動する際、それは「塑性ウェイク（航跡）」、つまり永久に変形した材料の跡を残しました。これはボートが通った後の航跡のようなものです。亀裂が通過するためには材料を押し潰さなければならないため、単純な脆性亀裂の理論が示唆するよりも多くのエネルギーを要しました。著者らはこの「実効的な靭性」を測定し、それが材料の基本となる靭性よりも約30%高いことを発見しました。これは、なぜ一部の材料が、基本的な公式が予測するよりも現実にはタフに見えるのかを説明しています。亀裂の先端における塑性変形が、追加の仕事をしているからです。

V字ノッチによる架け橋

モデルが複雑な形状に対しても機能することを証明するために、著者らは「V字ノッチ」（パックマンの角のような、材料への鋭い切り込み）のシミュレーションを行いました。彼らはこのノッチに力を加え、何が起こるかを観察しました。

荷重が増加するにつれ、材料は先ほどの三つの段階を経ることを発見しました。まず、先端に小さな塑性領域が形成され（「ぐにゃり」フェーズ）、次に凝集亀裂が成長し始め（「粘着テープ」フェーズ）、最後に脆性亀裂が突如として発生し、材料内を駆け抜けます。

ここで最も印象的な結果は、亀裂が発生するのに必要な力が「普遍的な法則」に従ったことです。ノッチが非常に鋭い（亀裂のような）場合でも、非常に鈍い（直線的なエッジのような）場合でも、モデルの予測は単一の曲線上に収束しました。この曲線は、材料の強度（壊すためにどれほど強く引く必要があるか）と、靭性（壊すためにどれほどのエネルギーが必要か）を滑らかに結びつけています。これは、モデルがどの公式を使うべきかを推測することなく、あらゆる形状の破壊を予測できることを示唆しています 있습니다。

なぜこれが重要なのか

著者らは、これらの結果が実際の材料を用いた物理的な実験ではなく、シミュレーションと数学的証明に基づいていることに注意を払っています。しかし、シミュレーションは堅牢であり、既知の理論的限界と完璧に一致しています。彼らは、強度の劣化方法を変更することで、柔軟性と正確さを兼ね備えたツールを作り出したことを示しました。

以前は、強くて脆い材料をシミュレートしたい場合、計算不可能（メッシュを細かくしすぎてコンピュータがクラッシュする）なモデルを使うか、間違った答えを出すモデルを使うかのどちらかしかありませんでした。この新しいアプローチでは、「正則化長」（コンピュータのぼやけ）を単なる数値的なツールとして扱うことができます。鋭い亀裂を得るために小さくすることも、計算力を節約するために大きくすることもでき、それでも材料の強度と靭性は全く同じまま保たれます。

要するに、この論文は破壊に関する統一された言語を提供しています。それは、塑性（押し潰し）、凝集破壊（剥がれ）、そして脆性破壊の世界を、一つの整合性のあるフレームワークへと統合したものです。私たちはこれらの異なる破壊形態のどれかを選ぶ必要はないのです。それらは、材料がどのように壊れるかという物語における、異なる章に過ぎないのです。エンジニアや科学者にとって、これは、橋、航空機、あるいは生物組織のような複雑な破壊を、これまで到達不可能だったレベルの詳細さと正確さでシミュレートできる可能性を意味しています。