| 查看: 1221 | 回复: 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 |
» 猜你喜欢
求调剂
已经有23人回复
301求调剂
已经有15人回复
304求调剂(085602,过四级,一志愿985)
已经有17人回复
302分求调剂 一志愿安徽大学085601
已经有12人回复
288环境专硕,求调材料方向
已经有23人回复
环境专硕调剂
已经有6人回复
22408 调剂材料
已经有6人回复
285求调剂
已经有12人回复
求调剂
已经有6人回复
085600材料与化工301分求调剂院校
已经有19人回复
2楼2024-01-26 08:48:41
3楼2024-01-27 09:35:09
luyaobao
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2457.9
- 散金: 1091
- 红花: 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














回复此楼