4. 原子尺度裂纹模型构建
在分子动力学（MD）中，裂纹扩展模拟的第一步不是“跑起来”，而是“把模型搭对”。一个可运行的原子尺度裂纹模型，本质上要同时回答四个问题：原子放在哪里（单晶还是多晶）、裂纹开在哪里（初始构型）、晶体朝哪个方向（取向）、自由表面怎么处理（表面能与边界）。这四者共同决定了后续加载时裂纹是沿解理面脆断、在晶界处偏折，还是发射位错钝化。
一、单晶与多晶：先确定“缺陷环境”
单晶模型是最基础的验证平台。以面心立方（FCC）金属为例，常用 LAMMPS 的 lattice fcc 3.615 建立约 40×40×4 个晶胞的薄板，原子数约 2.5 万，厚度方向采用周期性边界以近似平面应变。单晶的价值在于：裂纹路径唯一、应力场干净，便于与线弹性断裂力学（LEFM）的应力强度因子 KI 场对比，验证势函数和加载方式是否正确。
多晶模型则用于研究晶界作用。常用 Voronoi 镶嵌法：在模拟盒内随机撒 N 个种子点，每个种子对应一个晶粒，晶粒内原子按随机取向旋转后的晶格填充，落在不同 Voronoi 胞边界的原子被判定为晶界原子。例如 200×200×20 Å 的盒子放 12 个种子，平均晶粒尺寸约 50 Å，晶界体积分数可达 15% 以上。此时裂纹可能沿晶、穿晶或混合扩展，取决于晶界能与解理能之比。
二、裂纹初始构型：几何与原子删除
裂纹初始构型通常用“狭缝”表示：在板中央沿 x 方向删除一个长 2a、高约 2–3 个原子层厚的矩形区域，形成上下两个自由表面。以 100 Å×100 Å 的薄板为例，取 a=15 Å，裂纹尖端位于 x=±15 Å。删除原子后必须做能量最小化（如共轭梯度法），让尖端原子弛豫到平衡位置，否则初始构型中残留的原子重叠会在加载瞬间产生虚假应力波。
更真实的做法是构造“尖锐裂纹”：在尖端处按 1/r 奇异性布置原子，或直接采用椭圆孔洞再压缩。工程上常用后者：先删除椭圆区域，再沿 y 方向施加小应变使孔洞闭合为裂纹，这样尖端曲率半径可控制在 1–2 个晶格常数，接近理论解理裂纹。
三、晶体取向：决定解理面与滑移系
晶体取向直接决定裂纹沿哪个晶面扩展。对 FCC 金属，常见解理面为 {100}，滑移面为 {111}。若把 x 轴设为 [100]、y 轴设为 [010]、z 轴设为 [001]，则裂纹面为 (010)，扩展方向为 [100]，属于典型 I 型解理。若把 y 轴改为 [110]，裂纹面变为 (110)，此时滑移系与裂纹面夹角改变，尖端更易发射位错，表现为韧性行为。
取向设置必须在建模时一次性确定，并在整个模拟中保持一致。常用做法是先生成标准取向的单晶，再用旋转矩阵把晶格旋转到目标取向，最后按模拟盒坐标切割。注意：旋转后要重新检查周期性边界是否仍匹配，否则会在边界处产生人为失配。
四、表面能与自由表面处理
自由表面是裂纹模型的核心边界。表面能 γs 定义为创建单位面积自由表面所需的能量，FCC 金属典型值约 1–2 J/m²。在 MD 中，表面能通过“相同原子数、有表面与无表面”两块体系的能量差除以表面积得到。若势函数给出的 γs 明显偏离实验值，裂纹扩展的临界应力强度因子 KIC 也会失真。
自由表面处理包括三点：第一，裂纹上下表面不施加任何约束，允许原子自由弛豫；第二，外边界若为自由表面，需留足真空层（通常 ≥ 10 Å），避免原子与自身镜像相互作用；第三，若采用周期性边界，则裂纹面不能跨周期，否则裂纹会与自身镜像连通，失去物理意义。
五、一个可运行的建模流程示例
下面给出单晶 FCC 铝裂纹模型的 LAMMPS 关键片段，取向为 x[100]、y[010]、z[001]，裂纹沿 x 扩展，裂纹面为 (010)：
units metal
lattice fcc 4.05
region box block 0 100 0 100 0 20
create_box 1 box
create_atoms 1 box
# 删除裂纹区域：长30 Å，高6 Å
region crack block 35 65 47 53 0 20
delete_atoms region crack
# 自由表面：z方向非周期
boundary p p f
pair_style eam/alloy
pair_coeff * * Al99.eam.alloy Al
# 弛豫
min_style cg
minimize 1e-12 1e-12 10000 10000
运行后应检查：原子总数是否与预期一致、裂纹尖端是否出现明显原子重排、势能是否收敛到 -3.5 eV/atom 附近（铝的典型值）。若势能偏高，说明删除区域过大或弛豫不足。
完成上述步骤后，一个可运行的裂纹模型就建立起来了。此时不要急于加载，先做一次无载荷的 NVT 弛豫并输出原子构型，用 OVITO 观察裂纹尖端是否保持尖锐、晶界是否清晰、自由表面是否平整。只有确认初始状态干净，后续的应力-应变曲线和裂纹扩展路径才具有物理意义。