The fracture resistance of elastic networks increases with the density of defects like a random walk

技術的サマリー：弾性ネットワークの破壊抵抗は、ランダムウォークのように欠陥密度の増加に伴って増大する

問題提起
無秩序なバネネットワークは、セラミックスやポリマーネットワークから機械的メタマテリアルに至るまで、多様な材料における破壊を研究するための基礎的なモデルとして機能する。破壊力学における中心的な物理量は、既存の支配的き裂の進展に対する抵抗を測定する見かけの破壊エネルギー $G_c$ である。無秩序性が局所的な不均一性によるき裂の停止を通じて $G_c$ を増大させることは確立されているが、この強化が無秩序の度合い（具体的には欠陥密度）にどのように定量的に依存するかは、依然として十分に理解されていない。既存のき裂がない場合、欠陥は一般的に核生成を促進し、破壊抵抗を低下させることが以前の研究で示されている。しかし、予き裂を有する系では、関係はより複雑である：低欠陥密度では、無秩序の増加がき裂進展に対する抵抗を高める可能性があるが、高密度では材料が弱化する。予き裂の領域において、微視的欠陥統計と巨視的破壊エネルギーを結びつける特定のメカニズムについては、さらなる理論的および定量的な明確化が必要である。

手法
著者らは、ランダムに配置された欠損結合がき裂進展に及ぼす影響を調査するために、準静的負荷下にある 2 次元三角形バネネットワークを採用している。

• モデル系： ネットワークは、平衡長 $\ell_r$ と剛性 $k$ を持つバネで構成される。バネはひずみが閾値 $\varepsilon_{max} = 1$ を超えると急激に破断する。
• 欠陥の導入： 完全な格子から、結合の割合 $\nu$ がランダムに除去される。
• シミュレーションプロトコル： 系は、一様な垂直変位増分を通じて負荷される。ネットワークは、FIRE 最小化アルゴリズムを用いて静的平衡状態に緩和される。ひずみ閾値を超えるバネは、最も伸びたものから順次除去され、き裂進展をシミュレートするための特定のプロトコルに従って除去される。
• 主要指標： 見かけの破壊エネルギー $G_c(a)$ は、初期長さ $L_c$ から $L_c + a$ へき裂を進展させるために必要な臨界弾性エネルギー解放率として定義される。本研究は、特定のき裂進展 $a_f$ における開始エネルギー $G_c(0)$ と破壊エネルギー $G_c(a_f)$ に焦点を当てている。
• 理論的枠組み： 著者らは、欠損結合によって引き起こされる弾性不均質性を、等価な局所破壊エネルギー地形 $\Gamma_{loc}(a)$ onto マッピングする。彼らは、任意の点における見かけの破壊エネルギーが、その点までのき裂経路で遭遇した最大局所破壊エネルギーによって決定されると仮定し、確立された平面き裂停止の理論を適用する。

主要な結果

1. 欠陥による強化： 欠損結合の濃度 $\nu$ を固定すると、見かけの破壊エネルギー $G_c(a)$ はき裂進展 $a$ とともに増大する。最初の結合が破断した瞬間に破壊する完全なネットワークとは異なり、無秩序なネットワークは、不安定な成長と停止が交互に現れる間欠的なき裂進展を示す。
2. 欠陥密度とのスケーリング： 等価な局所破壊エネルギー $\Gamma_{loc}$ の変動の標準偏差は、欠損結合割合の平方根 $\sqrt{\nu}$ に比例して増大する。その結果、平均破壊エネルギー $\langle G_c(a_f) \rangle$ も $\sqrt{\nu}$ に比例してスケーリングする。
3. ランダムウォーク機構： この $\sqrt{\nu}$ スケーリングは、個々の欠損結合による摂動を独立した確率変数としてモデル化することで説明される。局所破壊エネルギーへの総摂動は、これらの個々の効果の重ね合わせであり、ランダムウォークに類似している。総摂動の分散は欠損結合の数に比例して線形に増大し、標準偏差が $\sqrt{\nu}$ に比例するようになる。
4. 指数尾部と対数的成長： 局所破壊エネルギー $\Gamma_{loc}$ の確率密度関数 (PDF) は、純粋なガウス分布ではなく、指数尾部を示す。この指数尾部は、見かけの破壊エネルギーに対してギューベル極値分布をもたらす。その結果、$G_c(a)$ の中央値はき裂進展 $a$ に対して対数的に増大する（$G_c \propto \ln a$）。
5. き裂停止の支配性： 本研究は、$\Gamma_{loc}$ の局所変動によるき裂停止が、研究対象の領域における支配的な強化機構であることを実証している。き裂経路は直線からわずかに逸脱する（粗さ）可能性があるが、平面き裂停止モデルは観測された挙動を正確に予測しており、経路の粗さは局所不均質性によって引き起こされる停止に比べて二次的な役割しか果たしていないことが示唆される。

意義と主張
本論文は、微視的構造の無秩序性と巨視的破壊エネルギーとの間の定量的な関連性を確立することを主張している。個々の欠陥からの摂動のランダムウォーク的な重ね合わせにより、局所破壊エネルギーの変動が $\sqrt{\nu}$ としてスケーリングすることを示すことで、著者らは無秩序材料の強化を予測するための理論的枠組みを提供している。局所エネルギー地形における指数尾部の同定は、予き裂を有しない系とは異なる挙動である、き裂進展に伴う破壊エネルギーの対数的増加を説明する。

著者らは、これらの知見は、欠損結合が微細空隙や微細き裂に対応し、かつ欠陥密度が低く、破壊がほぼ直線的な経路で進行する限り、他の格子構造や連続体系にも拡張されると予想されると指摘している。彼らは、平面き裂進展の仮定と結合摂動の独立性を含む制限を明示的に述べており、これらの仮定は、非常に高い欠陥密度や、き裂前面の平均化が生じる 3 次元系では破綻する可能性があることを認めている。この研究は、微視的欠陥統計と巨視的靭性を結びつけることで、多様な無秩序材料における破壊エネルギーの定量的予測への道を開くことを目指している。