| 查看: 1360 | 回复: 1 | |||
[求助]
用velocity命令进行拉伸时，为什么只有左侧两个原子跑出界外，其他原子都不动？
|
|
units metal dimension 3 boundary s s s atom_style atomic variable latparam equal 5.42 # ----------------------- ATOM DEFINITION ---------------------------- lattice diamond ${latparam} region box block 0 6 0 6 0 15 #units box create_box 1 box create_atoms 1 box mass 1 28.0855 # ------------------------ FORCE FIELDS ------------------------------ pair_style tersoff pair_coeff * * Si.tersoff Si # ------------------------- SETTINGS --------------------------------- ###################################### # EQUILIBRATION reset_timestep 0 ##################################### region 1 block INF INF INF INF INF 3 group lower region 1 region 2 block INF INF INF INF 12 INF group upper region 2 group boundary union lower upper group mobile subtract all boundary compute new mobile temp ##################################### fix 1 lower setforce NULL NULL 0.0 #固定下边的原子 fix 3 upper setforce NULL NULL 0.0 fix 2 all nvt temp 1.0 1.0 1 velocity lower set 0.0 0.0 0.0 velocity upper set 0.0 0.0 0.3 velocity mobile ramp vz 0.0 0.3 z 3 17 sum yes thermo 200 timestep 0.003 thermo_modify temp new dump 1 all atom 500 dump.crack run 300000 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有15人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有13人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有20人回复
pxjdal
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 668.8
- 散金: 1
- 帖子: 91
- 在线: 100.9小时
- 虫号: 1017662
- 注册: 2010-05-13
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
2楼2018-11-22 10:49:41
20