基于分子动力学模拟的金属构件的弹-塑性分解方法

来源：计算数学杂志2020年第3期北京时间：

作者：崔俊芝, 余翌帆

单位：LSEC, ICMSEC, 中国科学院数学与系统科学研究院, 北京 100190;中国科学院大学数学科学学院, 北京 100049

摘要：针对金属多晶材料构件的分子动力学（MD）模拟，本文提出了一种新的弹-塑性分解方法.文章将MD运动轨迹分解为结构变形和热振动，给出了计算结构变形的方法和近似公式；针对金属多晶材料构件的当前构型，给出了基于FCC|BCC晶胞和四原子占位的四面体单元相组合的连续变形函数及计算变形梯度的算法；利用原子-连续关联模型，发展了计算当前构型应力场和弹性张量的算法.分析了当构件承受过大载荷后在材料内部所产生的微观缺陷，并将其分类标定为位错、层错、挛晶界、晶界和空位等；对于层错和挛晶界讨论了在弹性卸载过程中应满足的刚体运动约束方程；利用极小势能原理构造了基于当前构型的弹性卸载算法，进而给出了完整的基于MD模拟的计算弹-塑性应变的算法.最后，针对单晶铜纳米线拉伸的MD模拟，计算了弹-塑性应变场，验证了本文方法的合理性.本文提出的基于MD模拟的弹-塑性分解方法，为从微观到宏观的多尺度和多模型耦合计算提供了算法支撑.

关键词：分子动力学模拟,结构变形,弹-塑性分解,极小势能原理,弹性卸载算法

基金资助：<p>国家自然科学基金重点项目（51739007）资助.</p>