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

compute stress/atom

這個(gè)命令可以計(jì)算某一原子組中單個(gè)原子的應(yīng)力，如：

compute 1 mobile stress/atom NULL

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

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

計(jì)算后的原子應(yīng)力可通過dump輸出到軌跡文件中，導(dǎo)入到ovito即可繪制應(yīng)力云圖。

更多l(xiāng)ammps教程請(qǐng)關(guān)注公眾號(hào)：
lammps加油站

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。