| 查看: 1246 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿武理材料工程348求调剂
已经有6人回复
-
08工科 320总分 求调剂
已经有11人回复
-
一志愿华中农业071010，总分320求调剂
已经有6人回复
-
323求调剂
已经有5人回复
-
284求调剂
已经有5人回复
-
317求调剂
已经有16人回复
-
287求调剂
已经有9人回复
-
求调剂
已经有6人回复
-
308求调剂
已经有3人回复
-
环境学硕288求调剂
已经有6人回复
pxjdal
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 668.8
- 散金: 1
- 帖子: 91
- 在线: 100.9小时
- 虫号: 1017662
- 注册: 2010-05-13
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
2楼2018-11-22 10:49:41
