A Combined Microbeam and Phase-Field Approach to Identify the Toughness and Ultimate Strength of Amorphous Silica

技术摘要：结合微束与相场法确定非晶二氧化硅的韧性与极限强度

问题陈述
硅酸盐玻璃的机械可靠性从根本上受限于其对表面缺陷的高度敏感性，这些缺陷作为应力集中源，会在材料达到其理论固有强度之前触发断裂。尽管已有大量研究表征了玻璃的脆性断裂响应，但很少有研究能够同时量化断裂表面能（临界能量释放率，$G_c$）和固有抗拉强度（$\sigma_c$）。传统的断裂力学通常依赖于无限尖锐裂纹的假设，这会导致奇异性问题，从而无法预测无预置裂纹完整样品的失效。相反，基于强度的准则无法预测存在裂纹时的失效。此外，大多数微观力学测试方法仅能得出应力强度因子（$K_{IC}$），使得固有强度和材料的特征长度尺度（$\ell_c$）在很大程度上难以获取。因此，需要一个一致性的框架来弥合脆性断裂表征与强度控制机制之间的鸿沟。

方法论
作者提出了一种结合聚焦离子束（FIB）微纳加工与相场断裂建模的新颖方法，用于评估非晶二氧化硅（$SiO_2$）的脆性性质。

• 实验设计： 利用 FIB 铣削技术从二氧化硅玻璃板中制备了两种互补的微梁几何结构：

  1. 单缺口微梁： 具有尖锐缺口，用以模拟格里菲斯（Griffith）脆性机制，在此机制下，断裂由裂尖处的应力强度主导。
  2. 骨形微梁： 一种新设计的几何结构，具有无缺口的窄颈部及钝化奇异点。该设计探测的是强度控制机制，即失效受材料固有抗拉强度而非预置裂纹严重程度的影响。
     两种几何结构的微梁均在空气中接受了三点弯曲测试，测试使用纳米压痕仪进行。

• 数值框架： 实验结果通过耦合相场公式的有限元分析（FEA）进行处理。该方法利用连续损伤变量（$d$）使裂尖附近的应力场正则化，从而无需预先设定显式的裂纹路径。该模型纳入了两个关键材料参数：临界能量释放率（$G_c$）和固有材料长度尺度（$\ell_c$）。

  • 尖锐缺口梁允许独立于 $\ell_c$ 对 $G_c$ 进行标定。
  • 骨形梁由于具有钝化奇异点，对 $\ell_c$ 具有敏感性，从而实现了对其的测定。
  • 与两种几何结构断裂力一致的参数空间交集，得出了 $G_c$ 和 $\ell_c$ 的唯一值。
  • 随后，利用均匀相场解从这些参数中推导出极限抗拉强度（$\sigma_c$）。

主要结果
这种结合实验与数值的方法得出了在空气中测试的非晶二氧化硅的定量属性：

• 临界能量释放率 ($G_c$): $5.1 \text{ J/m}^2$。
• 应力强度因子 ($K_{IC}$): $0.61 \text{ MPa}\cdot\text{m}^{1/2}$（由 $G_c$ 推导得出）。
• 固有长度尺度 ($\ell_c$): $9.1 \text{ nm}$。
• 极限抗拉强度 ($\sigma_c$): $6.8 \text{ GPa}$。

这些数值与先前文献报道的范围一致。值得注意的是，本研究报告了第一个实验得出的非晶二氧化硅正则化长度尺度（$\ell_c$）的值。结果还证实，骨形梁的断裂力大约是缺口梁的两倍，验证了从断裂主导向强度主导机制的转变。

意义与主张
本文声称，这种结合的方法论成功地弥合了传统韧性表征与强度控制机制之间的鸿沟，而后者通常是仅依赖 $K_{IC}$ 的技术所无法触及的。

• 同步评估： 不同于仅能得出 $K_{IC}$ 的传统微观力学方法，该方法通过利用互补的几何结构和不预设显式裂纹路径的相场框架，同时提供了 $G_c$ 和 $\sigma_c$。
• 物理洞察： 研究表明，推导出的长度尺度 $\ell_c$ 不仅仅是一个数值正则化参数，而是与玻璃的原子结构物理相关。具体而言，$\ell_c$ 的值与材料网络内弥散损伤（键断裂和刚度降低）的空间范围相关，而这受到非晶结构中环尺寸分布的影响。
• 参数化理解： 参数化研究表明，缺口几何形状（深度和尖端半径）决定了从强度控制向格里菲斯脆性机制的转变。这为设计旨在针对特定失效模式的微梁实验提供了实践指导。
• 敏感性分析： 作者强调，尺寸测量误差（特别是梁高度和缺口深度）会显著影响确定的 $G_c$ 和刚度的准确性，强调了精确几何表征的必要性。

作者得出结论，该框架为评估包括其他非晶系统在内的脆性固体的韧性和强度提供了一条通用的途径，将宏观断裂行为与内在材料长度尺度及原子结构联系起来。