24小时热门版块排行榜    

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

2009云飞扬

铁杆木虫 (著名写手)

[求助] lammps中某个区域压强的求解

各位朋友,还是前面提到过的一个问题。我建立了一个立方体,中心挖去原子形成一个球孔,在球孔中填入了另一种原子A,想求解球孔区域A原子所形成的压强是多大,目前知道应该局部压强采用compute stress/atom 来计算(单位是能量——压强.体积),但是这样就需要求解球孔区域A原子所形成的体积有多大,计算体积目前是采用的compute  voronoi/atom 指令。但是在计算过程中总是觉得体积算的不对,相关的具体化指令如下,请知道的朋友帮忙指点指点,应该怎么计算才能正确得到。先谢谢啦。
compute   stress  voidA  stress/atom He_bubble
compute         stress1  voidA  reduce sum c_stress[1]
compute         stress2  voidA  reduce sum c_stress[2]
compute         stress3  voidA  reduce sum c_stress[3]
compute   volum   voidA  voronoi/atom only_group
variable  p equal -(${c_stress1}+${c_stress2}+${c_stress3})/(3.0*c_volum)
体积计算部分是不是有问题呢?希望知道的朋友能够详细说说,先谢谢啦
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2009云飞扬

铁杆木虫 (著名写手)

各位朋友体积计算哪儿 算出来的量值太大是怎么回事呢

发自小木虫Android客户端
4楼2019-01-26 11:22:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

2009云飞扬

铁杆木虫 (著名写手)

2楼2019-01-24 14:00:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2009云飞扬

铁杆木虫 (著名写手)

3楼2019-01-24 18:47:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ruibingqilin

铁虫 (小有名气)

你好,我刚接触lammps,按照例子在熟悉操作。从NIST官网上复制了势函数到text文件里,命名为Al99.eam.alloy,放在计算的文件夹里了,但运行后还是出现Cannot open EAM potential file Al99.eam.alloy 这个错误。请问应该怎么解决? 谢谢!
5楼2019-02-12 11:49:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见