| 查看: 1898 | 回复: 1 | ||
[求助]
关于lammps里examples的问题（shear） 已有1人参与
|
|
各位大神，希望能指点一下。 在equilibrate之前，lower与upper的温度并没有设定，只是设定了mobile的温度，请问那其余两组原子的初速度怎样？还有就是原子的质量好像没有给啊，不用给定吗？ 以下是shear的程序： # 3d metal shear simulation units metal boundary s s p atom_style atomic lattice fcc 3.52 region box block 0 16.0 0 10.0 0 2.828427 create_box 3 box lattice fcc 3.52 orient x 1 0 0 orient y 0 1 1 orient z 0 -1 1 & origin 0.5 0 0 create_atoms 1 box pair_style eam pair_coeff * * Ni_u3.eam neighbor 0.3 bin neigh_modify delay 5 region lower block INF INF INF 0.9 INF INF region upper block INF INF 6.1 INF INF INF group lower region lower group upper region upper group boundary union lower upper group mobile subtract all boundary set group lower type 2 set group upper type 3 # void #region void cylinder z 8 5 2.5 INF INF #delete_atoms region void # temp controllers compute new3d mobile temp compute new2d mobile temp/partial 0 1 1 # equilibrate velocity mobile create 300.0 5812775 temp new3d fix 1 all nve fix 2 boundary setforce 0.0 0.0 0.0 fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0 fix_modify 3 temp new3d thermo 25 thermo_modify temp new3d timestep 0.001 run 100 # shear velocity upper set 1.0 0 0 velocity mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes unfix 3 fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0 fix_modify 3 temp new2d #dump 1 all atom 100 dump.shear #dump 2 all image 100 image.*.jpg type type & # axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0 #dump_modify 2 pad 4 #dump 3 all movie 100 movie.mpg type type & # axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0 #dump_modify 3 pad 4 thermo 100 thermo_modify temp new2d reset_timestep 0 run 3000 |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有8人回复
-
关于filecode
已经有4人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
mgqqlwq
铜虫 (著名写手)
- 应助: 12 (小学生)
- 金币: 12831.1
- 散金: 108
- 红花: 14
- 帖子: 1282
- 在线: 207.6小时
- 虫号: 2436148
- 注册: 2013-04-25
- 专业: 凝聚态物性I:结构、力学和
2楼2015-12-10 18:59:12
100