| 查看: 693 | 回复: 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 |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
20
投票: