lammps案例：单晶Cu压缩案例讲解

大家好，我是小马老师。

最近几天好几位粉丝朋友咨询如何用lammps模拟压缩过程。

因此，本文分享一个lammps模拟单晶铜压缩案例代码。

压缩过程其实很简单，就是反向的拉伸过程，把拉伸方向反过来就可以了。

lammps模拟拉伸有deform拉伸法、velocity拉伸法、addforce拉伸法，同理，这些方法也可以应用到压缩模拟。

在压缩过程中，可使用nvt或者npt系综，压缩应力-应变输出方法参考之前的文章：

​​lammps教程：应力-应变数据的计算与输出​​

单晶Cu的压缩采用deform拉伸法，沿着z轴进行压缩，压缩速率为0.01/ps。

在代码中，压缩速率设置为-0.01实现压缩（设置为0.01为拉伸）。

拉伸过程中使用npt系综，x和y方向压力为0，周期性边界条件。

压缩结果如下图所示（Ovito-CNA）：

案例in文件代码如下，代码已经注释，供初学者参考。这个案例不需要data文件，力场文件也是lammps自带：

END

公众号案例代码以及学习交流，请加QQ群：754749935

关注微信公众号：lammps加油站