1. 课程导入与工程问题界定

裂纹扩展是工程结构失效链条中最关键的一环。无论是航空发动机涡轮盘、核电压力容器接管，还是芯片封装中的焊点界面，绝大多数灾难性失效都不是因为材料整体强度不足，而是局部裂纹在循环载荷或腐蚀环境下萌生、稳定扩展，最终跨越临界尺寸引发失稳断裂。宏观断裂力学用应力强度因子 K、能量释放率 G 和 J 积分描述裂纹尖端的场强，能给出"是否扩展"和"扩展多快"的工程判据，但它把裂纹尖端处理为连续介质中的一个数学奇点，无法回答更底层的问题：裂纹究竟沿哪个原子面起裂？位错发射与解理断裂谁先发生？晶界偏析如何改变扩展路径？氢原子在尖端聚集后怎样降低键合强度？

分子动力学（MD）模拟正是在这一尺度缺口上发挥作用。它把材料离散为原子，用经验势函数或机器学习势描述原子间相互作用，通过牛顿方程积分得到每个原子随时间的位置与速度，从而在飞秒到纳秒、纳米到百纳米的时空尺度上直接观察裂纹尖端的原子级过程。需要明确的是，MD 不是万能工具，它有清晰的适用边界：时间尺度通常限制在纳秒量级，远低于工程疲劳的秒、小时量级；空间尺度限制在百万原子以内，难以直接模拟宏观构件；势函数的精度决定了结果的可信度，不同势对同一材料的断裂行为可能给出定性不同的结论。因此，MD 的定位是"机理揭示"而非"寿命预测"，它回答的是"为什么"和"如何"，而不是"还能用多久"。

本课程以一个贯穿全程的工程问题为牵引：在单晶 α-Fe 中，预置一条 I 型张开裂纹，在 300 K 温度下施加恒定应变率拉伸，观察裂纹尖端的原子尺度扩展行为，并分析温度、应变率与晶体取向对扩展路径和临界应力的影响。选择 α-Fe 是因为其势函数成熟（如 Embedded Atom Method 势）、晶体结构简单（BCC）、工程应用广泛（管线钢、压力容器钢的主要基体），适合作为入门体系。选择 I 型张开裂纹是因为其加载方式明确、边界条件易于实现，便于初学者建立"加载—响应—判据"的完整认知。

课程项目最终交付物包括四类：第一，一套可复现的 LAMMPS 输入脚本，涵盖建模、弛豫、加载、输出全流程；第二，一组关键模拟结果，包括裂纹尖端原子构型快照、应力—应变曲线、裂纹长度随时间的演化数据；第三，一份分析报告，用 OVITO 可视化结果，结合位错提取算法（DXA）识别尖端塑性变形机制，给出扩展路径与临界应力的判断依据；第四，一份参数记录表，完整记录势函数类型、时间步长、温度控制方式、应变率、边界条件等关键输入，确保他人可复现。这四类交付物对应"能跑、能看、能判、能复现"四项能力，也是本课程后续每一章逐步搭建的目标。

为了让项目目标可操作，需要把上述工程问题翻译成 MD 可执行的输入。具体来说，模拟盒子取 BCC 铁单晶，尺寸约 40a₀ × 40a₀ × 2a₀（a₀ ≈ 0.286 nm），总原子数约 12.8 万，三个方向均采用周期性边界条件以消除自由表面效应。裂纹通过在盒子中央删除一个原子层薄片实现，初始裂纹长度约 5 nm，裂纹面法向沿 y 方向，加载方向沿 y 方向（I 型张开）。原子间相互作用采用 EAM 势（如 Mendelev 2003 铁势），时间步长取 1 fs，先以 NPT 系综在 300 K 下弛豫 50 ps 使系统达到热力学平衡，再切换到 NVT 或 NVE 系综施加恒定应变率 10⁸ s⁻¹ 进行拉伸。每 0.5 ps 输出一次原子构型与应力张量，总模拟时间约 500 ps，对应约 5% 的工程应变。

判断裂纹是否扩展，不能只看原子构型"看起来裂了"。工程上需要三个可量化的观察量：第一，裂纹尖端位置随时间的演化，可通过识别裂纹面两侧原子配位数突变来确定；第二，系统应力—应变曲线，当应力达到峰值后出现明显下降，通常对应裂纹起裂或大规模塑性变形；第三，尖端区域位错密度或孪晶体积分数的变化，用 DXA 算法从原子构型中提取。若应力峰值对应的应变约为 8%~12%，且峰值后裂纹长度单调增加，同时尖端出现位错发射，则可判断为"塑性主导的延性扩展"；若应力峰值后裂纹长度突增而位错密度变化很小，则倾向于"解理脆性扩展"。这些判断依据将在后续章节中逐步展开，本章只需建立"输入—过程—输出—判据"的完整认知框架。

最后需要强调工程问题界定的边界意识。MD 模拟给出的临界应力、扩展路径和机制判断，是在特定势函数、特定温度、特定应变率下的结果，不能直接外推到工程构件的疲劳寿命。它的价值在于：当宏观实验观察到"某温度下断裂模式由脆转韧"时，MD 可以给出原子尺度的解释；当新材料或新工艺缺乏实验数据时，MD 可以快速筛选可能的机制方向。把 MD 放在"机理研究工具"的位置上，既不夸大也不贬低，才是工程上负责任的使用方式。本课程后续 29 章将围绕这一体系，从建模、势函数选择、弛豫、加载、输出、可视化、位错分析到参数敏感性，逐步构建完整的裂纹扩展 MD 模拟能力。