| 查看: 1455 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者feng_le将赠送您 10 个金币 | ||
[求助]
lammps建模，不移动，模拟求助
|
||
|
我建了一个简单模型，上面是金刚石，中间是铜，下面是镍，在跑完驰豫之后我让金刚石向下移动，但是没有按我设想的带着铜移动，而是在接触到铜之后直接给铜压碎了，望指教，是哪里没有编写正确。图片是我建模的，和我期望的样子 # melt exercise C-Cu-Ni echo screen units metal dimension 3 boundary p p f atom_style atomic # create atoms region box block 0 80 0 50 -100 100 units box create_box 3 box lattice fcc 3.524 region s1 block 0 80 0 50 -11 0 units box create_atoms 1 region s1 units box lattice fcc 3.615 region s2 sphere 40 25 22.5 10 units box create_atoms 2 region s2 units box lattice fcc 3.567 region s3 block 0 80 0 50 45 56 units box create_atoms 3 region s3 units box mass 1 58.69 mass 2 64 mass 3 12 pair_style hybrid/overlay lj/cut 1.0 eam tersoff eam/alloy morse 2.5 pair_coeff 1 3 morse 0.1 2.2 2.4 pair_coeff 2 3 lj/cut 0.0342 2.903 pair_coeff 1 1 eam Ni_u3.eam pair_coeff 2 2 eam Cu_u3.eam pair_coeff * * tersoff SiC.tersoff NULL NULL C pair_coeff * * eam/alloy CuNi.eam.alloy Ni Cu NULL group Ni region s1 neighbor 2.0 bin neigh_modify every 2 delay 10 check yes exclude group Ni Ni group Cu region s2 group C region s3 neighbor 2.0 bin neigh_modify every 2 delay 10 check yes exclude group C C velocity all create 300 5812775 dist gaussian compute 1 Cu temp #/com qu diao zhi xin fix 1 Cu nve fix 2 Cu temp/berendsen 300.0 300.0 0.1 fix 3 Ni setforce 0 0 0 #nve/noforce fix 4 C setforce 0 0 0 #nve/noforce timestep 0.001 run 3000 write_data chiyu.data write_restart chiyu.restart unfix 1 unfix 2 unfix 4 fix 1 Cu nvt temp 300 300 0.1 fix 5 C move linear 0 0 -1 run 300000 write_data dong.data write_restart dong.restart thermo 100 thermo_style custom step temp c_1 ke pe press dump 1 all custom 100 test.lammpstrj id type x y z  QQ图片20190511084501.jpg  QQ图片20190511084508.png  QQ图片20190511084535.png@oxox6085@smutao |
回复此楼» 猜你喜欢
基金开始函评了吗？
已经有10人回复
-
召集青年编委
已经有9人回复
-
博士论文被抄袭
已经有39人回复
-
期末给学生划重点都是什么话术啊
已经有18人回复
-
with editor 两个月了，什么原因？
已经有8人回复
-
普通院校药学硕士，做合成的，感觉找不到工作
已经有11人回复
-
求博导
已经有10人回复
-
Nature一直在编辑手里，考虑好几天了，是悬了吗
已经有10人回复
-
博士白读了
已经有41人回复
-
一直找不到审稿人
已经有6人回复
2楼2019-05-11 08:49:06
3楼2019-07-15 13:56:30
4楼2020-05-27 19:04:17
5