24小时热门版块排行榜    

查看: 2395  |  回复: 1

yaohonghao

金虫 (小有名气)

[求助] Invalid attribute in dump custom command问题求助 已有1人参与

本人毕业设计做界面力学性能模拟,使用候怀宇和张新平老师编写的《材料科学与工程中的计算机应用》学习LAMMPS。
练习例5-3 Cu晶体拉伸过程的分子动力学模拟时,出现了Invalid attribute问题。查看了manual中dump命令没有发现问题。
如果把dump删了程序就出现Invalid Syntax in variable formula错误
求各位大神抽空帮个忙看看~~~
# ------------------------ INITIALIZATION ----------------------------
units                metal
atom_style        atomic
boundary         p p p

# ----------------------- ATOM DEFINITION ----------------------------
variable         a equal 3.615
lattice         fcc $a
region                 mybox block 0 10 0 6 0 6
create_box         1 mybox
create_atoms        1 box

# ------------------------ FORCE FIELDS ------------------------------
pair_style         eam
pair_coeff        1 1 Cu_u3.eam

# ------------------------- SETTINGS ---------------------------------

timestep         0.001
thermo                 200
thermo_style        custom step temp pe vol lx ly lz press
velocity         all create 300 72437 dist gaussian
fix                 1 all npt temp 300 300 0.5 iso 0 0 0.8 drag 1.0
run                 3000
unfix                 1
variable         L equal lx
variable        L0 equal $L
reset_timestep 0
fix                 npt all npt temp 300 300 0.5 y 0 0 0.8 z 0 0 0.8 drag 1.0

fix                 def  all deform 1 x erate 0.005 units box remap x

variable        ss equal (lx-${L0}/)/${L0}
variable         p1 equal v_ss
variable         p11 equal -pxx
variable         p22 equal -pyy
variable         p33 equal -pzz

compute                 p all pe/atom
compute         cs all centro/atom fcc
compute                s all stress/atom NULL

thermo                 1000
thermo_style         custom step temp v_ss v_p11 v_p22 v_p33 pe etotal lx press

dump                 1 all cfg 200 cu * .cfg mass type xs ys zs id c_p c_cs c_s[1] c_s[2] c_s[3]
dump_modify         1 element Cu

fix                 pp all print 100 "${p1} ${p11} ${p22} ${p33}" file ex3_ss.txt screen no
run                  200
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

alundilong

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
每个case运行都会产生日志文件log.lammps
打开log.lammps看看
看最后一行说啥
2楼2016-03-20 22:14:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yaohonghao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +3 gouxfjh 2026-08-28 3/150 2026-08-29 01:01 by jurkat.1640
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见