| 查看: 649 | 回复: 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 |
回复此楼» 猜你喜欢
各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。
已经有3人回复
-
化工学硕294分，求导师收留
已经有20人回复
-
327求调剂
已经有8人回复
-
297求调剂
已经有22人回复
-
085801电气专硕272求调剂
已经有11人回复
-
271求调剂
已经有33人回复
-
调剂
已经有27人回复
-
一志愿2110，化学学硕310分，本科重点双非求调剂
已经有20人回复
-
085600材料与化工，求调剂
已经有13人回复
-
本科郑州大学，一志愿华东师范大学282求调剂
已经有31人回复
10