24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2335  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 食品与营养(0955)271求调剂 +7 升格阿达 2026-04-12 7/350 2026-04-12 22:56 by yuyin1233
[考研] 295分求调剂 +13 ?要上岸? 2026-04-10 13/650 2026-04-12 15:37 by laoshidan
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成!
[考研] 299求调剂 +8 ZVVZ13 2026-04-08 8/400 2026-04-12 00:40 by 蓝云思雨
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 求调剂 +10 璃茉一定上岸 2026-04-10 10/500 2026-04-11 13:31 by 1005715100
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 一志愿中科院105500专业总分315求调剂 +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 一志愿哈工大,初试329,求环境科学与工程调剂! +11 余未辛 2026-04-06 11/550 2026-04-08 15:21 by screening
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
信息提示
请填处理意见