1. 纳米材料力学与分子动力学工程导论
纳米材料力学的核心困难，并不只是“尺寸变小了”，而是当特征尺寸进入 1–100 nm 区间后，材料的行为不再由连续介质假设单独支配。位错形核、表面应力、晶界滑移、界面键合、热涨落与自由表面效应会同时参与承载与失效过程。传统连续介质力学擅长描述宏观应力—应变场，但在纳米尺度上，弹性模量、屈服强度、断裂韧性都可能随尺寸、形状、表面状态和温度显著变化。分子动力学（Molecular Dynamics, MD）正是在这一尺度区间内，把材料离散为原子或粗粒化粒子，通过势函数描述相互作用，用牛顿方程推进轨迹，从而把“结构—能量—力学响应”直接联系起来。
本课程定位为“纳米材料力学 + 分子动力学工程”的交叉能力训练：既讲清纳米力学的物理问题，也讲清 MD 从建模、势函数选择、系综与边界条件设置，到力学量提取、误差评估与工程判读的完整链路。学习目标不是会点软件按钮，而是能够针对一个具体纳米力学问题，判断该用什么模型、需要哪些输入、输出哪些可验证的力学量，以及结果在什么条件下可信。
多尺度力学框架与 MD 的位置
纳米力学通常按特征尺度与主导机制分层理解。宏观连续介质层（>1 μm）以有限元、连续介质本构为主，关注应力集中、裂纹扩展与结构稳定性；介观层（10 nm–1 μm）常引入位错动力学、晶体塑性或相场方法，处理位错群、孪晶与相变；原子层（0.1–10 nm）则由 MD、蒙特卡洛或第一性原理承担，直接解析原子排布、键断裂与界面迁移。MD 的关键价值在于它是“力学量可计算”的原子级方法：给定原子坐标与势函数，即可计算势能、力、应力张量、温度与位移场，并进一步得到弹性常数、屈服强度、界面强度、摩擦系数与断裂行为。
核心问题：从“原子怎么动”到“材料怎么坏”
纳米力学要回答的核心问题可以归纳为四类。第一类是刚度与弹性常数：纳米线、纳米薄膜、纳米颗粒的杨氏模量是否与块体一致？表面原子配位不足会如何改变整体刚度？第二类是强度与塑性机制：当位错源稀缺时，纳米晶或纳米柱的屈服由表面位错形核、晶界滑移还是孪晶主导？第三类是界面与接触力学：异质界面、范德华界面、摩擦接触在原子尺度如何传递载荷？第四类是断裂与失效：裂纹尖端原子键断裂、空位聚集与应力腐蚀如何决定纳米结构的寿命？MD 的价值在于，它可以把上述问题转化为可计算的原子级力学量，而不是停留在定性推测。
MD 工程链路：输入、计算与输出
一个可复现的 MD 力学分析通常包含以下关键环节。首先是建模：确定晶体结构（FCC、BCC、HCP 或非晶）、晶向、尺寸、表面取向与缺陷类型。例如研究纳米铜柱压缩时，可建立截面 8 nm × 8 nm、高度 20 nm 的柱体，沿 [001] 方向加载，两端设置刚性压头或速度边界。其次是势函数选择：金属体系常用 EAM 或 MEAM，共价体系用 Tersoff 或 Stillinger-Weber，层状材料用 AIREBO 或 ReaxFF。势函数决定了弹性常数、层错能与界面能的准确性，必须与目标力学量匹配。第三是系综与边界条件：力学加载常用 NVT 或 NPT，温度 300 K，时间步长 1 fs，加载速率 10⁸–10¹⁰ s⁻¹。需要特别说明的是，MD 加载速率远高于实验，因此屈服强度通常偏高，需通过多速率外推或与实验趋势对比来判读。第四是力学量提取：通过维里应力公式计算原子级应力张量，对区域平均得到应力—应变曲线；通过共同邻域分析（CNA）或位错提取算法（DXA）识别位错、孪晶与相变；通过均方位移（MSD）判断扩散与熔化。最后是结果判读：弹性模量由应力—应变线性段斜率得到，屈服点由应力峰值或位错形核事件确定，界面强度由分离功—位移曲线峰值给出。
MD 在纳米力学中的角色边界
MD 不是万能的。它的时间尺度通常限制在纳秒到微秒，空间尺度限制在数百万原子以内，因此无法直接模拟宏观结构的长期疲劳或蠕变。它的势函数精度也决定了结果的可信度：EAM 对 FCC 金属的弹性常数和层错能较好，但对共价键断裂和化学反应描述不足。因此，工程实践中常采用“MD 提供本构参数与失效判据，连续介质模型承担宏观响应”的多尺度策略。理解 MD 的适用边界，比盲目追求大体系更重要。