| 查看: 1558 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者小毛驴起飞将赠送您 5 个金币 | ||
[求助]
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 |
回复此楼» 猜你喜欢
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
2楼2021-09-13 10:48:22