计算化学公社

标题: fix deform算拉伸过程的应力应变关系和fix ave/chunk命令 [打印本页]

作者
Author: 小毛驴起飞 时间: 2021-9-7 21:25
标题: 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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)