| 查看: 1397 | 回复: 0 | ||
[求助]
碰撞模型弛豫出现变形
|
|
建立一个球形弹丸撞击靶板的模型，弛豫时发生变形，请大家看看，有什么不妥？ log log.equil units metal dimension 3 boundary s p p atom_style atomic neighbor 2.0 bin neigh_modify every 5 delay 0 check yes lattice fcc 4.05 region box block 0 20 -20 20 -0.5 0.5 create_box 2 box create_atoms 1 box mass 1 27 mass 2 27 pair_style eam/fs pair_coeff * * Al_mm.eam.fs Al Al timestep 0.001 thermo 100 thermo_style custom step etotal pe ke temp vol press min_style cg minimize 0 1.0e-5 5000 20000 reset_timestep 0 velocity all create 300.0 47289 dist gaussian units box fix 1 all nvt temp 300.0 300.0 0.1 compute peatom all pe/atom compute keatom all ke/atom compute c all cna/atom 3.4571 dump 1 all custom 100 *.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom dump_modify 1 thresh c_c > 1 run 20000 unfix 1 region 1 sphere 2.0 0.0 0.0 2.0 side in group 1 region 1 set group 1 type 1 region 2 block 5 INF INF INF INF INF group 2 region 2 set group 2 type 2 group model union 1 2 group void subtract all model delete_atoms group void velocity model create 300.0 47289 dist gaussian units box fix 2 model nvt temp 300.0 300.0 0.1 dump 2 model custom 100 equil.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom dump 3 1 custom 100 equil1.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom dump 4 2 custom 100 equil2.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom dump 5 model custom 100 *.lammpstrj id type x y z vx vy vz c_c c_peatom c_keatom dump_modify 5 thresh c_c > 1 run 10000 unfix 2 undump 2 undump 3 undump 4 undump 5 write_restart restart.equil |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有17人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
越来越担心硕士毕业后没有竞争力了
已经有12人回复
-
闲聊
已经有6人回复
-
申博求助
已经有4人回复
点击这里搜索更多相关资源