| 查看: 1208 | 回复: 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 |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有8人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
pxjdal
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 668.8
- 散金: 1
- 帖子: 91
- 在线: 100.9小时
- 虫号: 1017662
- 注册: 2010-05-13
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
2楼2018-11-22 10:49:41