24小时热门版块排行榜    

查看: 2328  |  回复: 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
[考研] 298求调剂 +6 上岸6666@ 2026-03-20 6/300 2026-03-22 20:21 by edmund7
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 305分求调剂(食品工程) +4 Sxy112 2026-03-21 6/300 2026-03-22 15:26 by 无懈可击111
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
信息提示
请填处理意见