lammps教程：单原子应力计算及应力云图绘制方法

大家好，我是小马老师。
前一篇文章介绍了单原子温度的计算方法，有粉丝朋友问单原子应力计算及应力云图计算方法。
单原子应力的计算要比单原子温度计算方法简单，lammps内嵌了原子应力计算命令：

compute stress/atom

这个命令可以计算某一原子组中单个原子的应力，如：

compute 1 mobile stress/atom NULL

需要注意的是，compute stress/atom计算出来的量带有体积项，如果要转换为pxx相同的单位，需要除以原子体积。
可使用compute voronoi/atom命令计算单原子体积。
下面以输出单原子xyz三个方向应力为例，给出具体用法：

compute 1 all stress/atom