单晶Fe三向拉伸模拟

#-----------------模型设置---------

units        metal

atom_style      atomic

boundary      p p p

neighbor      0.3 bin

timestep      0.001

#-----------------建模------------

lattice        bcc 2.8552

region        box block 0 20 0 20 0 20

create_box      1 box

create_atoms    1 box

mass        1 56

#-----------------势能-------

pair_style      eam/fs

pair_coeff      * * Fe.eam.fs Fe

min_style      cg

minimize      1e-15 1e-15 5000 5000

#-----------------速度初始化-------

velocity      all create 300 876848 dist gaussian

#------------------npt弛豫---------

fix          1 all npt temp 300 300 0.5 iso 0 0 5

dump        1 all custom 1000 Fe_npt.xyz id type x y z

run          30000

unfix        1

undump       1

#-------------------应力应变设置----

variable       tmpx equal "lx"

variable           tmpy equal "ly"

variable           tmpz equal "lz"

variable        Lx0 equal ${tmpx}

variable        Ly0 equal ${tmpy}

variable         Lz0 equal ${tmpz}

variable        strainx equal "(lx - v_Lx0)/v_Lx0"

variable        strainy equal "(ly - v_Ly0)/v_Ly0"

variable        strainz equal "(lz - v_Lz0)/v_Lz0"

variable        stressx equal "-pxx/10000"

variable        stressy equal "-pyy/10000"

variable        stressz equal "-pzz/10000"

#-------------------------应力应变保存到文件---------------------------

fix         def3 all print 100 "${strainx} ${stressx} ${strainy} ${stressy} ${strainz} ${stressz}" &

          file stress-strain.dat screen no

#-------------------三向拉伸--------

fix          1 all nve

fix          2 all temp/rescale 10 300 300 10 1.0

fix          3 all deform 100 x erate 0.005 y erate 0.005 z erate 0.005 remap x units box

dump        1 all custom 100 Fe_tension.xyz id type x y z

run         50000