| 查看: 2595 | 回复: 3 | |||
| 【悬赏金币】回答本帖问题，作者nulixiaoda将赠送您 15 个金币 | |||
[求助]
lammps拉伸时变形特别大了还没有拉断 已有1人参与
|
|||
|
我在用lammps进行水合物的拉伸，但拉伸时三个尺寸都发成变形，而且水合物变形特别大了还是不会断裂。导出来的应力应变曲线也是不规律。请问各位大佬知不知道这是什么原因呢，小弟在这跪求了，谢谢各位！下面是我的in文件。 #水合物拉伸 #初始化 echo screen units real dimension 3 boundary p p p atom_style full #势函数 pair_style lj/cut/coul/long 14 kspace_style ewald/disp 1.0e-4 bond_style harmonic angle_style harmonic dihedral_style harmonic pair_modify mix arithmetic #定义文件 read_data "shw.data" #公共近邻 neighbor 0.3 bin neigh_modify every 2 delay 10 check yes page 100000 velocity all create 273 12345 mom yes rot no #定义变量 计算值 #应力 variable stress equal "-pxx/10000" #应变 variable l_x equal lx variable lx0 equal ${l_x} variable strain equal (lx-v_lx0)/v_lx0 #驰豫 timestep 0.001 fix 1 all npt temp 243 243 0.1 z 50 50 0.1 y 50 50 0.1 drag 0.1 thermo 1000 thermo_style custom step lx ly lz press pxx pyy pzz pe temp etotal thermo_modify lost ignore run 20000 #拉伸 reset_timestep 0 unfix 1 fix 1 all npt temp 243 243 0.1 y 50 50 0.1 z 50 50 0.1 drag 0.1 fix 2 all deform 10 x erate 0.001 units box fix 3 all print 100 "${strain} ${stress}" file strain_stress.txt screen no thermo 1000 thermo_style custom step lx ly lz press pxx pyy pzz pe temp etotal thermo_modify lost ignore dump traj all custom 500 shw.xyz id type x y z dump_modify traj sort id timestep 0.001 run 200000 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
2楼2022-07-04 17:19:20
|
您好我想请教您几个问题，以下是我的in文件 #npt弛豫 units real timestep 1 #npt弛豫 fix NPT all npt temp 270 270 100 x 98 98 1000 y 98 98 1000 z 0 0 1000 thermo_style custom time atoms step temp press vol lx ly lz enthalpy pe ke etotal density pxx pyy pzz pxy pxz pyz thermo ${thermo} dump 1 all xyz ${thermo} npt_${name}.xyz dump_modify 1 element O H C fix pt all print 200 "${istep} ${etol} ${han} ${T0} ${V} ${PE} ${P0}" screen no file npt_${name}.dat run ${runmd} undump 1 unfix pt unfix NPT reset_timestep 0 #stress-strain variable tmpz equal "lz" variable Lz0 equal ${tmpz} variable strainz equal "(lz - v_Lz0)/v_Lz0" variable stressz equal "-pzz/10000" #z向拉伸 fix 2 all npt temp 270 270 100 x 98 98 1000 y 98 98 1000 fix 3 all deform 100 z erate 0.00000001 remap x units box thermo_style custom time atoms step temp press vol lx ly lz density pxx pyy pzz pxy pxz pyz thermo ${thermo} dump 1 all xyz ${thermo} tension_${name}.xyz dump_modify 1 element O H C #保存应力应变数据 fix def3 all print 100 "${strainx} ${stressx} ${strainy} ${stressy} ${strainz} ${stressz}" screen no file stress-strain_${name}.dat run ${runten} 1.您所说的x方向拉伸，在NPT时，使x方向的压强调节为零，是在弛豫阶段设定为0吗？我是z方向拉伸，这么设置对吗？但是我这么设置应力依旧不是从0开始的，我检查了log文件发现我的模拟过程控不住压，压力波动巨大，是否是这个原因那？应该怎么解决？ 2.我的应力应变文件的屈服强度很大，有8GPa，这明显是有问题的，这是什么原因那？也是因为控制不住压力导致的吗？ 非常希望能得到您的指点，感谢。 |
3楼2024-01-11 17:06:33
4楼2024-01-11 17:13:36
