fix deform算拉伸过程的应力应变关系和fix ave/chunk命令

小毛驴起飞

lammps新手想求助一下，本来是按照网上的案例想用fix deform算拉伸过程（x方向）的应力应变曲线，但是算出来pxx，pyy，pxx都是随应变增大而减小的趋势，完全不符合， 然后就想用fix ave/chunk 输出受力，in文件给出，initial.txt是有输出fx，fy，fz都是接近于0，但是stress.txt这里就没有输出受力的数据了，想请大家指导提点下是哪里出错了，感谢！！！
#NVT ensemble for relaxation
fix                1 all nvt temp 300 300 100
dump 1 all custom 500 nvt.xyz id type x y z vx vy vz
thermo_style        custom step temp press vol density etotal epair ebond eangle edihed pxx pyy pzz lx ly lz
thermo          500
timestep        0.5
reset_timestep         0
run                50000
unfix 1
undump        1
   

variable tmp equal "lx"
variable L0 equal ${tmp}
variable strain equal "(lx - v_L0)/v_L0"
variable stressx equal "pxx/10000"
dump 1 all custom 500 traj.lammpstrj id type x y z ix iy iz vx vy vz
fix                1 all npt temp 300 300 100 y 0 0 1000 z 0 0 1000 drag 2
fix                2 all deform 1 x erate 0.000005 units box remap x
fix     3 all print 500 "${strain} ${stressx}" file strain_stress.txt screen no
thermo_style        custom step temp v_strain v_stressx pxx pyy pzz lx ly lz etotal epair ebond eangle edihed
thermo          500
compute st all chunk/atom bin/1d x lower 40 units box
fix    11 all ave/chunk 1 1 5 st fx fy fz ave one file stress.txt
timestep   0.5
reset_timestep        0
run                20000
unfix 1
unfix 2
unfix 3
undump 1
unfix 11