24小时热门版块排行榜    

查看: 2329  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂,英语六级551,有科研经历 +6 生物工程调剂 2026-03-17 10/500 2026-03-22 20:22 by edmund7
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 08工科 320总分 求调剂 +11 梨花珞晚风 2026-03-17 11/550 2026-03-22 17:42 by luoyongfeng
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +3 vv迷 2026-03-22 3/150 2026-03-22 15:18 by 杨杨杨紫
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
信息提示
请填处理意见