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

技术摘要：弹性网络的断裂阻力随缺陷密度的增加而增强，其规律类似于随机游走

问题陈述
无序弹簧网络是研究从陶瓷、聚合物网络到机械超材料等多种材料断裂现象的基础模型。断裂力学中的一个核心量是表观断裂能 $G_c$，它衡量了对预存主导裂纹扩展的阻力。虽然已确定无序性可通过局部不均匀性导致的裂纹钉扎来增加 $G_c$，但这种增韧效应与无序程度（特别是缺陷密度）之间的定量依赖关系仍知之甚少。先前的研究表明，在没有预存裂纹的情况下，缺陷通常会促进成核并降低断裂阻力。然而，在预裂纹系统中，这种关系更为复杂：在低缺陷密度下，增加无序性可能会增强对裂纹扩展的阻力，而在高密度下，材料则会弱化。在预存裂纹机制下，将微观缺陷统计与宏观断裂能联系起来的具体机制，需要进一步的理论澄清和定量说明。

方法论
作者采用二维三角弹簧网络在准静态加载下，研究随机移除键对裂纹扩展的影响。

• 模型系统：网络由具有平衡长度 $\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)$。他们应用了既定的平面裂纹钉扎理论，假设任意点的表观断裂能由裂纹路径上直至该点所遇到的最大局部断裂能决定。

主要结果

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）表现出指数尾部，而非纯高斯分布。这种指数尾部导致表观断裂能呈现 Gumbel 极值分布。因此，$G_c(a)$ 的中位数随裂纹扩展量 $a$ 对数增加（$G_c \propto \ln a$）。
5. 裂纹钉扎的主导作用：研究表明，在研究的机制中，由 $\Gamma_{loc}$ 局部波动引起的裂纹钉扎是主要的增韧机制。虽然裂纹路径可能会略微偏离直线（粗糙度），但平面裂纹钉扎模型能准确预测观察到的行为，表明与局部非均匀性诱导的钉扎相比，路径粗糙度的作用较小。

意义与主张
本文声称建立了微观结构无序与宏观断裂能之间的定量联系。通过证明局部断裂能的波动由于单个缺陷扰动的类随机游走叠加而按 $\sqrt{\nu}$ 标度，作者为预测无序材料的增韧提供了理论框架。识别出局部能量景观中的指数尾部，解释了断裂能随裂纹扩展的对数增加，这种行为与没有预存裂纹的系统截然不同。

作者指出，这些发现预计可扩展到其他晶格结构和连续介质系统，其中缺失键对应于微孔洞或微裂纹，前提是缺陷密度保持较低且断裂沿近乎直线的路径进行。他们明确指出了局限性，包括平面裂纹扩展的假设和键扰动的独立性，并承认这些假设在极高缺陷密度下或在发生裂纹前缘平均化的三维系统中可能会失效。这项工作旨在通过将微观缺陷统计与宏观韧性联系起来，为各种无序材料的断裂能定量预测铺平道路。