| 查看: 1297 | 回复: 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 |
回复此楼» 猜你喜欢
E口会评
已经有7人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有19人回复
-
生命口会评
已经有13人回复
-
明天E口面上会评
已经有9人回复
-
轴手性化合物的确定
已经有5人回复
-
MSER送审了还被拒稿
已经有6人回复
-
内心匮乏
已经有7人回复
-
祈祷青基必中
已经有16人回复
-
27届辽宁大学应届毕业生申博
已经有3人回复
-
关于如何从代码看上不上会
已经有23人回复
pxjdal
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 668.8
- 散金: 1
- 帖子: 91
- 在线: 100.9小时
- 虫号: 1017662
- 注册: 2010-05-13
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
2楼2018-11-22 10:49:41