| 查看: 1220 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿211，0703化学305分求调剂
已经有9人回复
-
22408 调剂材料
已经有5人回复
-
273求调剂
已经有23人回复
-
一志愿华东理工085601材料工程303分求调剂
已经有5人回复
-
0703化学
已经有22人回复
-
华南理工0703化学，总分336求调剂
已经有11人回复
-
初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文
已经有10人回复
-
308求调剂
已经有10人回复
-
一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业
已经有10人回复
-
调剂
已经有11人回复
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
