纳米材料的力学性能与其块体形式存在显著差异：当材料尺寸缩小至纳米量级时，表面原子比例急剧增加、位错源密度降低，材料表现出”越小越强”的尺寸效应。纳米线的屈服强度可达块体材料的数倍甚至一个数量级，这一特性对微机电系统、柔性电子与纳米增强复合材料具有重要意义。理解纳米线变形的微观机制（表面位错成核、孪生、无定形化）是设计高性能纳米材料的基础。
LAMMPS 分子动力学可对纳米线的拉伸变形进行原子级模拟，获得完整的应力-应变曲线并直接观察位错成核与滑移过程。本项目以铝（Al）纳米线为对象，用 LAMMPS（EAM 势）模拟单轴拉伸，分析应力-应变响应、屈服强度与直径的依赖关系。
金属铝采用嵌入原子方法（EAM）势描述，其总能量为对势与嵌入能之和：E = ΣᵢFᵢ(ρᵢ) + ½Σᵢⱼφᵢⱼ(rᵢⱼ)。EAM 势通过嵌入函数 F(ρ) 有效描述金属键的电子密度依赖，能准确复现铝的晶格常数、弹性常数、堆垛层错能与表面能，是金属纳米力学模拟的可靠选择。
纳米线沿轴向（[100] 或 [110]）施加恒应变率变形：每个时间步按设定应变率缩放轴向盒子长度并重新标定原子坐标，同时保持横向压力为零（自由表面）。原子受力由 EAM 势计算，通过积分获得轨迹。工程应力由维里应力或外载荷计算，工程应变为轴向长度变化与原长之比。典型的应变率取 1×10⁹ s⁻¹（MD 时间尺度的常用值，远高于实验但应变率效应可外推）。
弹性段应力随应变线性上升，斜率为杨氏模量；当应力达到屈服强度时，表面处成核位错并滑移，应力骤降或平台化，进入塑性流动。纳米线屈服强度显著高于块体，因为表面位错源稀缺、位错成核需要更高的临界应力。
图 9 的左面板给出 Al 纳米线的工程应力-应变曲线（蓝色曲线）。弹性段（应变 0–0.078）应力线性上升，斜率即杨氏模量约 68 GPa，与实验铝的弹性模量（约 70 GPa）接近；在应变 0.078 处应力达到屈服点 5.2 GPa（红色虚线标注），随后应力继续小幅上升至峰值 6.5 GPa（应变 0.12），之后应力下降并趋于塑性流动平台（约 4.5 GPa）。屈服强度 5.2 GPa 是块体铝（约 0.05–0.1 GPa）的数十倍，凸显纳米线”越小越强”的特征。弹性段结束后的应力跌落对应表面位错成核与第一批位错滑移，塑性平台的应力波动反映后续位错增殖与湮灭的动态平衡。
图 9 的右面板给出屈服强度随纳米线直径的变化（橙色曲线，直径 2–10 nm）。屈服强度从直径 2 nm 的约 7.8 GPa 单调下降至直径 10 nm 的约 4.6 GPa，呈现典型的”越小越强”尺寸效应。其物理机制为：直径越小，表面原子比例越高，位错源越稀缺，位错成核的临界应力越高；同时纳米线的表面积-体积比增大，表面应力对屈服贡献增强。该尺寸依赖关系与实验测得的金属纳米线强度-直径幂律行为定性一致，为纳米器件的强度设计提供了定量参考。
两面板结合，左面板给出纳米线完整的变形响应（弹性-屈服-塑性），右面板揭示强度对尺寸的依赖规律，共同完成 Al 纳米线拉伸力学的系统研究，为理解纳米尺度塑性、指导纳米结构强度设计提供原子级依据。
本项目用 LAMMPS 模拟了 Al 纳米线的单轴拉伸，获得杨氏模量约 68 GPa、屈服强度 5.2 GPa（直径 4 nm），屈服强度随直径减小从 4.6 升至 7.8 GPa。案例图以应力-应变曲线与尺寸效应曲线双面板呈现。
后续可拓展：(1) 用位错提取算法可视化位错成核与滑移，揭示塑性机制；(2) 研究温度与应变率对屈服强度的影响，建立本构模型；(3) 引入表面缺陷或孪晶界，评估其对强度与韧性的调控。
