3. 裂纹扩展的物理图像与关键参量
在分子动力学（MD）模拟中，裂纹扩展并不是一个“原子键断裂”的简单事件，而是宏观连续介质力学参量与原子尺度离散行为之间相互映射的过程。理解这一映射，是设计模拟、解释轨迹、判断失效模式的前提。本节围绕四个核心概念展开：应力强度因子 K、能量释放率 G、裂纹尖端场，以及位错发射与解理的竞争，并说明它们如何在 MD 中建立联系。
一、宏观参量：K 与 G 的定义与关系
在线弹性断裂力学（LEFM）框架下，裂纹尖端附近的应力场具有 r-1/2 奇异性，其强度由应力强度因子 K 描述。对 I 型（张开型）裂纹：
σ_ij(r,θ) = K_I / √(2πr) · f_ij(θ) + 高阶项其中 r 为到裂纹尖端的距离，θ 为相对裂纹面的角度，f_ij(θ) 为角分布函数。当 K_I 达到材料断裂韧度 K_IC 时，裂纹起裂。能量释放率 G 则从能量角度描述裂纹扩展单位面积所释放的弹性能，对平面应变 I 型裂纹：
G_I = K_I² (1 - ν²) / E其中 E 为杨氏模量，ν 为泊松比。这一关系是宏观参量之间的桥梁，也是 MD 中把原子尺度能量变化“翻译”为宏观断裂参量的关键公式。
二、裂纹尖端场：从连续介质到原子尺度
LEFM 的 r-1/2 奇异场只在裂纹尖端外围（r 远大于原子间距、小于塑性区尺寸）成立。在 MD 中，裂纹尖端是一个由几十到几百个原子组成的离散区域，存在键断裂、原子重排、局部非晶化等过程。因此，模拟中需要区分三个尺度：
- 远场区：可用 LEFM 描述，K 或 G 是有效控制参量；
- 过程区：原子键断裂与位错发射发生，连续介质假设失效；
- 尖端核心区：几个原子间距内，应力集中、键角畸变、局部温度升高。
在 MD 中，通常通过原子级应力（Virial 应力）或位移场拟合来提取 K。例如，对 I 型裂纹，可在裂纹前方沿 y 方向取若干原子层，计算其 σ_yy，再按 σ_yy = K_I / √(2πr) 拟合得到 K_I。这一过程要求 r 落在 K 主导区内，否则拟合结果会偏离。
三、位错发射与解理竞争：微观机制的分岔
裂纹尖端在加载过程中面临两种竞争路径：
- 解理扩展：原子键直接断裂，裂纹沿特定晶面（如 BCC 的 {100}、FCC 的 {111}）向前推进；
- 位错发射：裂纹尖端发射位错，释放应力集中，使裂纹钝化，甚至停止扩展。
Rice 和 Thomson 提出的判据指出，位错发射的临界应力强度因子 K_IE 与解理临界值 K_IC 的相对大小决定竞争结果。若 K_IE < K_IC，则位错发射先发生，材料表现为韧性；反之则解理脆断。在 MD 中，这一竞争可通过观察裂纹尖端原子构型演化来判断：若尖端前方出现滑移面、原子层错或堆垛层错，则位错发射占优；若裂纹面直接分离、无位错形核，则解理占优。
四、宏观-微观参量映射：MD 中的实现路径
将宏观参量映射到 MD 模拟，通常遵循以下流程：
- 建立含裂纹的原子模型：选择晶向、裂纹面、裂纹长度 a，通常 a 需大于 5 nm 以减小边界效应；
- 施加远场加载：通过位移控制或力控制施加 I 型载荷，记录应力-位移曲线；
- 提取 K 或 G：用原子应力拟合 K，或用能量差计算 G；
- 判断竞争机制：观察裂纹尖端是否发射位错，计算 K_IE 与 K_IC 的相对大小；
- 反馈修正：若模拟结果与实验或理论不符，检查势函数、边界条件、加载速率等。
需要强调的是，MD 中的加载速率远高于实验（通常 108–1012 s-1），因此 K_IE 和 K_IC 的绝对值会偏高，但两者的相对关系仍可反映材料韧性趋势。此外，温度、势函数、裂纹取向都会显著影响竞争结果，需在模拟设计中明确控制变量。
五、模拟中的关键观察与判断依据
在 MD 运行过程中，判断裂纹扩展机制需要关注以下信号：
| 观察对象 | 解理占优 | 位错发射占优 |
|---|---|---|
| 裂纹尖端原子构型 | 裂纹面直接分离，无滑移 | 出现层错、孪晶或位错环 |
| 应力-位移曲线 | 脆性突降 | 屈服平台或锯齿状 |
| 局部温度 | 变化小 | 明显升高（塑性耗散） |
| K 提取值 | 接近 K_IC | 低于 K_IC，尖端钝化 |
这些判断依据需要结合势函数、温度、加载速率综合解读。例如，高加载速率会抑制位错发射，使材料表现得更脆；而高温则促进位错形核，增加韧性。因此，在 MD 模拟中，宏观参量映射不是一次性的计算，而是需要与微观机制观察反复对照的过程。