24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者小毛驴起飞将赠送您 5 个金币

小毛驴起飞

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 18
帖子: 7
在线: 4.8小时
虫号: 26825730
注册: 2021-07-19
专业: 聚合物共混与复合材料

[求助] lammps新手想求助一下，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

回复此楼