Peridynamic modeling of the crack velocity dependence via an incubation time fracture criterion

本研究は、ホマライト100に関するラヴィ・チャンダールとナウスの実験をモデル化するために、孵化時間破壊基準を備えたペリダイナミックアプローチを採用し、一定のき裂速度におけるモードI応力拡大係数の変動およびより高い速度におけるマイクロ分枝の発生が、動的破壊におけるき裂速度依存性の本質に関する新たな洞察をもたらすことを明らかにした。

要約

脆いプラスチック、例えばホマライト100のシートを想像してください。その上を亀裂が走っている様子を思い浮かべてみましょう。物理学の昔の時代、科学者たちは、亀裂の移動速度が分かれば、それを推進する「応力」（圧力）の大きさを正確に計算できると考えていました。まるで、車が時速60マイルで走行しているなら、エンジンが正確に100馬力を発生しているに違いない、と考えるようなものです。シンプルですね？

しかし、1980年代の実験はこれが真実ではないことを示しました。ある時、亀裂は全く同じ速度で進んでいても、それを押し進める圧力は劇的に異なっていたのです。まるで、2台の車がどちらも時速60マイルで走っているのに、一方は小さなエンジンを持ち、他方はロケットブースターを搭載しているかのようでした。科学者たちは困惑しました：なぜ同じ速度なのに「押し力」が異なるのか？

この論文は、著者が新しい種類のコンピュータシミュレーションを用いてこの謎を解く探偵物語です。

探偵の道具：ペリダイナミクス

亀裂のコンピュータモデルの多くは、ドミノの連鎖のようなものです。1つのドミノが倒れると、次のドミノを押します。しかし、もしドミノが欠けていれば（亀裂があれば）、連鎖は途切れ、数学は行き詰まります。

著者が用いたのは「ペリダイナミクス」と呼ばれる手法です。これを連鎖ではなく、蜂の群れとして考えてみてください。すべての蜂は、中央に隙間があっても、一定の距離内にある他のすべての蜂と会話できます。もし1匹の蜂が飛び去れば（亀裂が形成されれば）、他の蜂はその蜂との会話をやめますが、群れの残りの部分は問題なく動き続けます。これにより、コンピュータは混乱することなく破壊や亀裂の発生を処理できます。

秘密の材料：「潜伏時間」

この論文における真の画期的な点は、亀裂が実際に破壊される「タイミング」をどのように決定したかです。

従来の方法では、圧力が十分に高くなると、材料は瞬時に破壊されました。しかし、著者は「潜伏時間基準」と呼ばれる規則を用いました。

乾いた小枝を折ろうとしている様子を想像してください。引っ張って瞬時に折れるわけではありません。引っ張り、繊維が伸びて弱まるまで一瞬の間、そこに保持し、そして初めて折れます。その一瞬の間が「潜伏時間」です。

著者は、コンピュータ上の蜂の群れに、直前の数マイクロ秒間の圧力を記憶させるようにプログラムしました。材料が破壊されるのは、その短い「潜伏」期間における圧力の「平均」が十分に高い場合に限られます。これにより、材料が破壊を「決意」するにはわずかな時間が必要であるという事実が考慮されました。

彼らが発見したもの

彼らは、実際の実験と同様に、プラスチックの板を引き裂くシミュレーションを行いました。彼らが発見したのは以下の通りです。

1. 速度と圧力の謎： 実際の実験と同様に、彼らのコンピュータは、同じ亀裂速度であっても、圧力（応力拡大係数）は単一の数値ではないことを示しました。それはある範囲にありました。時には低く、時には高く、変動していました。
2. 「マイクロ分枝」効果： 亀裂がゆっくりと移動しているときは、まっすぐ進みました。しかし、速度が速まり（時速400メートルを超えると）、振動し始めました。それは、木の枝が小枝に分かれるように、微小な側面亀裂を sprout（芽生え）させ始めました。
   • 比喩： 走っているランナーを想像してください。一定のジョギングでは、まっすぐ走ります。しかし、最高速度でスプリントすると、バランスを保つためにわずかに揺れ、ジグザグに動き始めます。
   • 結果： これらの小さな「揺れ」（マイクロ分枝）が、圧力測定値を激しく上下させる原因となりました。これが、与えられた速度に対して圧力が一意ではない理由を説明しました。亀裂は、走っている間に物理的にその形状をわずかに変化させていたのです。

結論

この論文は、同じ亀裂速度に対して異なる圧力値が見られる理由は、亀裂が滑らかで完璧な直線ではないからだと結論付けています。それは、変動する混沌とした生き物なのです。

• 低速時： 亀裂は安定しており、圧力も比較的安定しています。
• 高速時： 亀裂は「マイクロ分枝」（微小な側面亀裂の芽生え）を始めます。この混沌が圧力を跳ね回りさせ、実験で見られた散らばりを生み出します。

この「蜂の群れ」（ペリダイナミクス）と「待機期間」（潜伏時間）を組み合わせることで、著者たちは、現実の実験が数十年にわたって示してきた、亀裂速度と圧力の間の messy（ごちゃごちゃした）、非一意な関係を成功裡に再現しました。彼らは、データ内の「ノイズ」が誤りではなく、高速で移動する亀裂の振る舞いにおける真の物理的特徴であることを証明しました。

技術概要

技術サマリー：孵化時間破壊基準を介したき裂速度依存性のペリダイナミックモデル化

問題提起
本研究は、動的破壊力学における中心的な未解決課題、すなわち、瞬間的なモードI応力拡大係数（SIF、$K_I$）がき裂伝播速度（$\dot{a}$）に依存する問題に取り組む。古典的な線形弾性破壊力学（LEFM）およびフレンドとロサキスに基づく動的基準は、しばしば$K_I$と$\dot{a}$の間の一意の関係を示唆する。しかし、ラヴィ・チャンダとナウスによる脆性ホマライト -100 ポリマー板の実験データは、単一のき裂速度が SIF 値の範囲に対応する非一意の依存性を示している。さらに、より高い速度ではマイクロ分枝が発生し、SIF 測定値に大きなばらつきが生じる。課題は、経験的に調整された速度依存材料特性に頼ることなく、この非一意性と関連する変動を数値的に再現することにある。

手法
著者らは、ラヴィ・チャンダとナウスの実験条件を再現するため、ホマライト -100 における動的き裂伝播をシミュレートするために、非局所的なペリダイナミック枠組みを採用する。手法の中核的な革新は、ペリダイナミック定式化に孵化時間破壊基準（ITFC）を統合することである。

• ペリダイナミック枠組み：シミュレーションは、オープンソースソフトウェア Peridigm で実装された通常の状態ベースのペリダイナミックモデルを利用する。材料は、球状のホライズン（$\delta$）内で相互作用する点に離散化され、空間微分演算子なしにき裂のような不連続を自然に扱うことを可能にする。
• 孵化時間破壊基準（ITFC）：瞬時破壊を仮定する古典的基準とは異なり、ITFC は、破壊には微細き裂などの準備過程に有限時間（$\tau$）が必要であると仮定する。この基準は、孵化期間$\tau$および特徴長さ$d$における時間平均応力を評価する。
  • 特徴長さ$d$は、臨界静的 SIF（$K_{Ic}$）との整合性を確保するため、準静的アイリウ基準から導出される。
  • ペリダイナミック実装は、遠隔応力（RS）破壊基準を適応させる。結合の破損は、孵化時間$\tau$（式 11）にわたって法線方向結合応力を積分することで決定される。
  • 計算複雑性を管理し、完全なマクロ分枝ではなくマイクロ分枝に焦点を当てるため、この基準は主き裂軸に平行（$0^\circ$）および垂直（$90^\circ$）な平面でのみ評価される。
• シミュレーション設定：モデルは、50 mm の予ノッチを有する 300 mm $\times$ 500 mm の試料をシミュレートする。動的負荷は、25 $\mu$s の立ち上がり時間を有する台形体積力パルスによって付与される。動的 SIF は、動的 J 積分のための移動長方形輪郭を用いて計算され、積分経路が主要なき裂先端を囲みつつ、発達中のマイクロ分枝を除外することが保証される。

主要な結果
数値シミュレーションは、$K_I$-$\dot{a}$関係に関する実験的に観察された現象を成功裡に再現する。

1. 非一意な$K_I$-$\dot{a}$依存性：結果は、ほぼ一定のき裂伝播速度に対して瞬間的 SIF が顕著な変動を示すことを確認する。その結果、モデルは、実験で観察された非一意性を反映して、与えられた速度に対する SIF 値の範囲を予測する。
2. 速度の安定化と過渡現象：実験ではき裂速度がほぼ瞬時に安定化することが示されるが、シミュレーションは定常状態に達する前に約 30–40 $\mu$s 持続する初期加速段階を示す。この過渡段階中、SIF の履歴は顕著な変動を示す。
3. マイクロ分枝とばらつき：より高い負荷率（400 m/s を超える速度をもたらす）において、シミュレーションはマイクロ分枝の発生を示す。この現象は、SIF 振動の振幅の著しい増加と相関する。高速度における SIF 値の数値的ばらつき（0.53–1.25 MPa$\sqrt{m}$）は、実験範囲（0.45–1.5 MPa$\sqrt{m}$）と合理的に一致する。
4. 機構の検証：本研究は、観測された SIF のばらつきが測定誤差のアーティファクトではなく、特に加速段階とマイクロ分枝の発生に起因するき裂先端の動力学によって物理的に正当化されることを実証する。

意義と主張
本論文は、ITFC のペリダイナミック実装が、広範な負荷率にわたる動的破壊靭性の評価およびき裂伝播のモデル化にとって価値あるアプローチを提供すると主張する。主な意義は、任意の速度依存材料パラメータを導入することなく、脆性破壊の速度感応性および$K_I$-$\dot{a}$関係の非一意性をモデルが捉える能力にある。

著者らは、モデルが実験的傾向を定性的に捉える一方で、速度安定化のタイミングおよび SIF 変動の正確な大きさに関して定量的な不一致が残ると控えめに指摘する。彼らは、非一意な依存性の成功した予測を、応力の時間的履歴を考慮する孵化時間の導入、および動的破壊に関連する特異場および不連続を自然に処理するペリダイナミックの非局所性に帰因する。本研究は、き裂加速段階、SIF 変動、およびマイクロ分枝の発生との相互作用が、実験的に観察された$K_I$-$\dot{a}$挙動の支配的機構であると結論づける。