24小时热门版块排行榜    

Znn3bq.jpeg
查看: 7788  |  回复: 2

cl1994861

金虫 (小有名气)

[求助] 关于fix ave/chunk command问题 已有1人参与

在运行lammps中examples 中的例子KAPPA中in.mp时出现  ”error :no value in fix ave/chunk command",我都没改它的代码,只是照着运行一下,这个例子的in文件的文本如下:
echo screen
variable        x equal 10
variable        y equal 10
variable        z equal 20

variable        rho equal 0.6
variable        t equal 1.35
variable        rc equal 2.5

units                lj
atom_style        atomic

lattice                fcc ${rho}
region                box block 0 $x 0 $y 0 $z
create_box        1 box
create_atoms        1 box
mass                1 1.0

velocity        all create $t 87287

pair_style        lj/cut ${rc}
pair_coeff        1 1 1.0 1.0

neighbor        0.3 bin
neigh_modify        delay 0 every 1

# 1st equilibration run

fix                1 all nvt temp $t $t 0.5
thermo                100
run                1000

velocity        all scale $t

unfix                1

# 2nd equilibration run

compute                ke all ke/atom
variable        temp atom c_ke/1.5

fix                1 all nve

compute layers all chunk/atom bin/1d z lower 1 units reduced
fix                2 all ave/chunk 10 100 1000 layers v_temp file profile.mp

fix                3 all thermal/conductivity 10 z 20

variable        tdiff equal f_2[11][3]-f_2[1][3]
thermo_style        custom step temp epair etotal f_3 v_tdiff

thermo                1000
run                20000

# thermal conductivity calculation
# reset fix thermal/conductivity to zero energy accumulation

fix                3 all thermal/conductivity 10 z 20

fix             ave all ave/time 1 1 1000 v_tdiff ave running
thermo_style        custom step temp epair etotal f_3 v_tdiff f_ave

run                20000

关于fix ave/chunk  command问题
M5B_V5XE5VZ9S0C`_7[8~U2.png
回复此楼

» 猜你喜欢

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

小小波bi

铁虫 (初入文坛)

楼主,compute chunk/atom 这个命令的问题解决了吗?
2楼2016-11-15 21:31:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2828512

铁虫 (初入文坛)

【答案】应助回帖

compute layers all chunk/atom bin/1d z lower 1 units reduced
这条命令中的1改成0.02,因为这条命令的意思是将all沿着z轴分成许多份,而你如果写的是1的话就是只有一份,应该输入小于1的值,例如输入0.02就是分成了1/0.02=50份
3楼2023-03-22 17:44:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cl1994861 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085410-273求调剂 +7 X1999 2026-04-10 7/350 2026-04-16 16:32 by noqvsozv
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +7 否极泰来2026 2026-04-15 9/450 2026-04-16 13:40 by liumingli7817
[考研] 290调剂生物0860 +38 哇哈哈,。 2026-04-11 44/2200 2026-04-16 09:52 by cuisz
[考研] 求调剂推荐 +8 小聂爱学习 2026-04-14 8/400 2026-04-16 07:22 by 学员JpLReM
[考研] 药学求调剂 +11 RussHu 2026-04-12 13/650 2026-04-15 19:07 by zhuwenxu
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 机械还有还有名额吗?太难了 +8 笑笑袁 2026-04-10 8/400 2026-04-14 08:44 by screening
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 346分,工科0854求调剂,专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 求调剂,一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
信息提示
请填处理意见