| 查看: 802 | 回复: 2 | |||
[交流]
计算出来的石墨烯弹性模量比理论值偏小 已有1人参与
|
| 大家好，我使用velocity来拉伸石墨烯，然后根据应力应变曲线计算出来的弹性模量只有150GPa左右，请问大家遇到过这种问题吗？ |
回复此楼» 猜你喜欢
085600材料与化工 求调剂
已经有10人回复
-
0703化学调剂，求各位老师收留
已经有6人回复
-
梁成伟老师课题组欢迎你的加入
已经有5人回复
-
中科大材料专硕319求调剂
已经有3人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有12人回复
-
伙伴们，祝我生日快乐吧
已经有21人回复
-
297一志愿上交085600求调剂
已经有5人回复
-
265求调剂
已经有4人回复
-
面上和青基一样限30页不合理
已经有7人回复
-
现在如何回避去年的某一个专家，不知道名字
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2021-04-29 11:25:04
|
你好！下面是我的脚本 # ------------------------ initialization ---------------------------- units metal boundary s s p atom_style full # ----------------- atom definition and force fields ----------------- read_data data.0429graphene pair_style airebo 3.0 pair_coeff * * ch.airebo c region upper block inf inf 77 inf inf inf units box region lower block inf inf inf -26 inf inf units box group upper region upper group lower region lower group boundary union lower upper group middle subtract all boundary variable yc equal ly variable l0 equal ${yc} variable strain equal "(ly - v_l0)/v_l0" variable strain1 equal "v_strain" variable stressz equal "-pyy/10000" # ----------------------------settings ------------------------------- timestep 0.001 thermo_style custom step time etotal pe ke temp lx ly lz press pxx pyy pzz density velocity middle create 298.15 53244 dist gaussian mom yes rot yes fix s1 boundary setforce 0 0 0 minimize 1.0e-10 1.0e-10 5000 10000 thermo 1 fix 222 all npt temp 298.15 298.15 1 z 0 0 0.1 dump 1 all atom 100 relax.g run 1000 unfix 222 undump 1 reset_timestep 0 velocity upper set 0.0 1 0.0 velocity middle ramp vy 0.0 1 y -26 77 sum yes fix 666 all nvt temp 298.15 298.15 0.001 thermo 1 thermo_modify lost ignore dump 888 all atom 100 tension.g run 150000 |
3楼2021-04-30 00:38:04
5