| 查看: 1253 | 回复: 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 |
回复此楼» 猜你喜欢
085600材料与化工329分求调剂
已经有5人回复
-
复试调剂
已经有18人回复
-
310求调剂
已经有16人回复
-
274求调剂
已经有7人回复
-
284求调剂
已经有8人回复
-
245求调剂
已经有7人回复
-
有爆料，一个青年教师卖房得400万，然后换了一个四青帽子
已经有11人回复
-
290求调剂
已经有15人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有16人回复
-
一志愿西交机械专硕求调剂
已经有9人回复
pxjdal
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 668.8
- 散金: 1
- 帖子: 91
- 在线: 100.9小时
- 虫号: 1017662
- 注册: 2010-05-13
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
2楼2018-11-22 10:49:41
