| 查看: 6461 | 回复: 3 | ||
小讨厌123新虫 (小有名气)
|
[求助]
lammps中如何求每个原子的总能量 已有1人参与
|
| 用lammps模拟pt—ar系统，只想求一下ar的每个原子的平均能量，这个语句和输出文件语句怎么求？谢谢 |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
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
5
投票:
