| 查看: 702 | 回复: 0 | ||
[求助]
简单Ni单晶的拉伸，老是出现原子丢失，可有人答疑解惑？
|
|
# in.ni units metal boundary p p p atom_style atomic dimension 3 lattice fcc 3.52 orient y 0 1 1 orient z 0 -1 1 region box block 0 24 0 13 0 2 create_box 1 box create_atoms 1 box region left block INF 0 INF INF INF INF region right block 23 INF INF INF INF INF group left region left group right region right group boundary union left right group mobile subtract all boundary region void cylinder z 42.24 32.4815 7.465 INF INF units box delete_atoms region void pair_style meam pair_coeff * * library.meam Ni ni.meam Ni neighbor 2 bin neigh_modify every 5 check yes timestep 0.001 velocity all create 300 345678 fix 1 all npt temp 300 300 0.1 y 0 0 5 z 0 0 5 thermo_style custom step temp etotal pe ke press pxx pyy pzz thermo 1000 run 10000 fix 1 all npt temp 300 300 0.1 y 0 0 5 z 0 0 5 fix tensile all deform 1 x erate 0.01 units box # define strain variable tmp equal "lx" variable L0 equal ${tmp} variable strain equal "(lx - v_L0)/v_L0" variable p1 equal "-pxx/10000" variable p2 equal "-pyy/10000" variable p3 equal "-pzz/10000" variable st equal "step" variable p0 equal "v_strain*100" fix def1 all print 1000 "${st} ${p0} ${p1} ${p2} ${p3}" file x-Ni.txt screen no dump 1 all custom 1000 dump.dat id x y z thermo 1000 run 41000 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
