| 查看: 2123 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
lammps软件无法dump出cfg/custom文件 已有1人参与
|
||
|
in文件如下： # FeNi units metal boundary p p p atom_style atomic #variable i loop 5 #variable x equal 300# read_data FeNi.lammps pair_style eam/alloy pair_coeff * * FeCuNi.eam.alloy Fe Ni timestep 0.002 thermo 100 neighbor 2.0 bin neigh_modify every 1 delay 0 check yes compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 3.0 velocity all create 2.5 17287 dist gaussian variable N equal step variable Ep equal pe variable Etotal equal etotal variable T equal temp #variable Press equal press #variable V equal vol variable K equal ke #velocity all create 2.5 825577 dist gaussian #timestep 0.002 thermo_style custom step temp etotal ke pe #dump 1 all xyz 1000 dump.xyz #dump 2 all custom 1000 dump.lammpstrj id type x y z c_3 c_4 #thermo 1000 #compute 3 all pe/atom #compute 4 all ke/atom #compute 5 all coord/atom 3.0 #dump 1 all custom 200 dump_*.lammpstrj id type xs ys zs c_3 c_4 c_5 #dump 1 all cfg 100 a*.cfg mass type xs ys zs #run 2000 dump 1 all cfg 100 layer*.cfg mass type xs ys zs run 2000 #fix 1 all nvt temp 2.5 2.5 0.2 #run 4000 #unfix 1 #fix 2 all npt temp 2.5 300 0.2 iso 0 0 10 #run 4000 #unfix 2 #fix extra all print 200 "${N} ${T} ${Etotal} ${Ep} ${K} " append data #fix 3 all nvt temp 300 300 0.2 #run 2000 #unfix 3 #unfix extra clear #next i #jump in.heating |
回复此楼» 猜你喜欢
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有13人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
4楼2019-02-15 15:02:21
2楼2016-07-07 12:30:32
江南小毛羊
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2380.6
- 散金: 346
- 红花: 3
- 帖子: 381
- 在线: 7024.2小时
- 虫号: 2054729
- 注册: 2012-10-11
- 专业: 凝聚态物性I:结构、力学和
3楼2016-07-10 01:13:56
28