| 查看: 6335 | 回复: 3 | ||
小讨厌123新虫 (小有名气)
|
[求助]
lammps中如何求每个原子的总能量 已有1人参与
|
| 用lammps模拟pt—ar系统，只想求一下ar的每个原子的平均能量，这个语句和输出文件语句怎么求？谢谢 |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
alundilong
木虫 (正式写手)
- 应助: 129 (高中生)
- 金币: 2019.1
- 红花: 18
- 帖子: 309
- 在线: 56小时
- 虫号: 595468
- 注册: 2008-09-08
- 性别: GG
- 专业: 传热传质学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
小讨厌123: 金币+10, ★有帮助 2015-11-10 12:32:31
感谢参与，应助指数 +1
小讨厌123: 金币+10, ★有帮助 2015-11-10 12:32:31
|
看一下 compute ke/atom 算出动能 http://lammps.sandia.gov/doc/compute_ke_atom.html compute pe/atom 算出势能 http://lammps.sandia.gov/doc/compute_pe_atom.html 用compute reduce，把两者加起来，就是总能量 http://lammps.sandia.gov/doc/compute_reduce.html 最后再计算出所有ar原子数(用variable中的count函数http://lammps.sandia.gov/doc/variable.html)，除一下就好了 |
2楼2015-11-10 07:00:29
小讨厌123
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 246
- 散金: 16
- 帖子: 55
- 在线: 39.6小时
- 虫号: 2922365
- 注册: 2014-01-09
- 性别: MM
- 专业: 工程热力学
|
谢谢你的回答 但是还有些地方不明白 1、计算了ke/atom和pe/atom不就是每个原子的能量吗？直接加起来他们不就是每个原子的能量吗， 为什么还要除ar的原子数呢 2、compute reduce这个命令没有看懂 compute 1 all reduce sum c_ke c_pe这样写对吗 3、我想输出的是平均能量随时间的变化，可是用fix ave/time有错误提示是fix ave/time compute does not calculate a scalar 我的程序是这样写的 compute ar_pe ar pe/atom #直接计算Ar中每个原子的动能和势能 compute ar_ke ar ke/atom compute ar_total ar reduce sum c_ar_ke c_ar_pe # 这样是不是将动能和势能的值加起来了 fix 1 all ave/time 1 1000 1000 c_ar_total file info.data #输出每一时刻的平均能量 可是这样写有错误 fix ave/time compute does not calculate a scalar 您能帮我改一下吗 看哪里有问题 ，该怎么改，我就是想要能量随时间的变化 谢谢 |
3楼2015-11-10 13:01:30
uy_naiT
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 66.1
- 散金: 79
- 帖子: 94
- 在线: 33.2小时
- 虫号: 6467863
- 注册: 2017-05-07
- 专业: 金属材料的合金相、相变及
4楼2018-03-21 16:12:02