24小时热门版块排行榜    

查看: 3606  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lcl1cxt

金虫 (小有名气)

[求助] lammps模拟中,dump输出错误 已有1人参与

错误 Dump custom compute does not calculate per-atom vector (../dump_custom.cpp:1171)怎么解决?
我用的compute s all stress/atom  myTemp
dump  2  all  cfg  200 dump.deform_*.cfg mass type xs ys zs id  c_s
手册中dump命令解释说可以输出才c_ID的 per-atom vector的啊,为什么输出不了呢? 大神们怎么输出每个原子的应力的呢?会的可以帮我解释一下吗?
回复此楼

» 收录本帖的淘帖专辑推荐

精华

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

努力
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lcl1cxt

金虫 (小有名气)

引用回帖:
2楼: Originally posted by xbb20121991 at 2015-08-06 12:57:04
compute s all stress/atom
计算出来结果有6个分量,分别对应xx,yy,zz,xy,xz,yz。你如果想dump输出,可能要写出c_s之类的。
不过你还要注意它的单位是stress*volume

谢谢你,我知道输出了c_s[1]是可以,但是怎么输出除以体积了的单个原子的应力?
努力
3楼2015-08-07 11:46:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

xbb20121991

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lcl1cxt: 金币+5, 有帮助 2015-08-07 11:17:50
compute s all stress/atom
计算出来结果有6个分量,分别对应xx,yy,zz,xy,xz,yz。你如果想dump输出,可能要写出c_s[1]之类的。
不过你还要注意它的单位是stress*volume
2楼2015-08-06 12:57:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbb20121991

银虫 (小有名气)

lammps手册写的
compute              peratom all stress/atom NULL
compute              p all reduce sum c_peratom[1] c_peratom[2] c_peratom[3]
variable     press equal -(c_p[1]+c_p[2]+c_p[3])/(3*vol)
thermo_style custom step temp etotal press v_press
4楼2015-08-07 13:00:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbb20121991

银虫 (小有名气)

不好意思,看错了你的问题。
手册上的是一个部分的应力

» 本帖已获得的红花(最新10朵)

5楼2015-08-07 13:02:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见