大家好，我是小馬老師。
前一篇文章介紹了單原子溫度的計算方法，有粉絲朋友問單原子應力計算及應力云圖計算方法。
單原子應力的計算要比單原子溫度計算方法簡單，lammps內嵌了原子應力計算命令：

compute stress/atom

這個命令可以計算某一原子組中單個原子的應力，如：

compute 1 mobile stress/atom NULL

需要注意的是，compute stress/atom計算出來的量帶有體積項，如果要轉換為pxx相同的單位，需要除以原子體積。
可使用compute voronoi/atom命令計算單原子體積。
下面以輸出單原子xyz三個方向應力為例，給出具體用法：

compute 1 all stress/atom NULL compute v all voronoi/atom variable stressx atom c_1[1]/c_v[1]/10000 variable stressy atom c_1[2]/c_v[1]/10000 variable stressz atom c_1[3]/c_v[1]/10000 dump 1 all custom 100 all.xyz id type x y z v_stressx v_stressy v_stressz

計算后的原子應力可通過dump輸出到軌跡文件中，導入到ovito即可繪制應力云圖。

更多lammps教程請關注公眾號：
lammps加油站

總結

以上是生活随笔為你收集整理的lammps教程：单原子应力计算及应力云图绘制方法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。