24小时热门版块排行榜    

查看: 3609  |  回复: 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的回帖

2009云飞扬

铁杆木虫 (著名写手)

引用回帖:
8楼: Originally posted by lcl1cxt at 2015-08-11 13:34:25
谢谢分享,前两天有人给我介绍了一种:
compute s all stress/atom  myTemp
compute vorovol all  voronoi/atom
variable Press1 atom  "c_s/(10000.0* c_vorovol)"
你觉得是否可以参考一下 ...

楼主 你的这三指令,只是算到单个原子的压强吗?如果要算到整个组的原子形成的压强,你需要集合吗?
13楼2019-01-16 15:55:07
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖
信息提示
请填处理意见