| 查看: 804 | 回复: 2 | |||
[交流]
计算出来的石墨烯弹性模量比理论值偏小 已有1人参与
|
| 大家好，我使用velocity来拉伸石墨烯，然后根据应力应变曲线计算出来的弹性模量只有150GPa左右，请问大家遇到过这种问题吗？ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
