| 查看: 1384 | 回复: 4 | ||
| 【悬赏金币】回答本帖问题，作者白马踏平川将赠送您 50 个金币 | ||
[求助]
水合物单向拉伸的时候屈服强度大了十倍是怎么回事？
|
||
|
各位老师，做水合物拉伸时候屈服强度大了10倍，单位是GPa，数值是6-8左右，这是怎么回事那？希望各位老师不吝赐教 以下是in文件 units real atom_style full dimension 3 boundary p p p pair_style lj/cut/tip4p/cut 1 2 1 1 0.1546 10.0 8.5 pair_modify mix geometric pair_coeff 1 1 0.1852 3.1589 pair_coeff 2 2 0.0 0.0 pair_coeff 3 3 0.05589 2.757 pair_coeff 4 4 0.15998 3.033 bond_style harmonic bond_coeff 1 450 0.9572 bond_coeff 2 570 1.149 angle_style harmonic angle_coeff 1 55 104.52 angle_coeff 2 295.41109 108 #PROSET timestep ${dt} neigh_modify every 1 delay 0 check yes #minimize min_style sd minimize 1.0e-20 1.0e-20 100000 100000 fix mom all momentum 10 linear 1 1 1 angular fix 1 all box/relax iso 0.0 velocity all create 270 $R mom yes rot yes dist gaussian velocity all scale 270 run 0 unfix 1 fix shake all shake 0.0001 20 0 b 1 a 1 #NVT thermo_style custom time atoms step temp press vol lx ly lz enthalpy pe ke etotal density evdwl ecoul epair ebond pxx pyy pzz pxy pxz pyz thermo 200 fix NVT all nvt temp 270 270 100 dump 1 all xyz 200 nvt.xyz dump_modify 1 element OW HW C O run ${runnvt} unfix NVT undump 1 reset_timestep 0 #npt弛豫 fix 1 all npt temp 270 270 100 x 98.69 98.69 1000 y 98.69 98.69 1000 z 0 0 1000 thermo_style custom time atoms step temp press vol lx ly lz enthalpy pe ke etotal density evdwl ecoul epair ebond pxx pyy pzz pxy pxz pyz thermo 200 dump 1 all xyz 200 npt.xyz dump_modify 1 element OW HW C O run 1000000 undump 1 unfix 1 reset_timestep 0 #stress-strain variable tmpx equal "lx" variable tmpy equal "ly" variable tmpz equal "lz" variable Lx0 equal ${tmpx} variable Ly0 equal ${tmpy} variable Lz0 equal ${tmpz} variable strainz equal "(lz - v_Lz0)/v_Lz0" variable stressz equal "-pzz/10000" fix 2 all npt temp 270 270 100 x 98.69 98.69 1000 y 98.69 98.69 1000 fix 3 all deform 100 z erate 0.00000002 remap x units box thermo_style custom time atoms step temp press vol lx ly lz density pxx pyy pzz pxy pxz pyz thermo 20000 dump 1 all xyz 20000 tension.xyz dump_modify 1 element OW HW C O #保存应力应变数据 fix def3 all print 20000 "${strainz} ${stressz}" screen no file stress-strain_${name}.dat run 10000000 |
回复此楼» 猜你喜欢
时间戳又变了
已经有15人回复
-
科研孤儿太难了
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
今天基金会出结果吗？20260819
已经有15人回复
2楼2024-01-26 08:48:41
3楼2024-01-27 09:35:09
luyaobao
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2457.9
- 散金: 1876
- 红花: 23
- 沙发: 8
- 帖子: 1533
- 在线: 150.7小时
- 虫号: 7847640
- 注册: 2018-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2024-02-18 01:34:31
5楼2024-02-24 11:52:29
