| 查看: 777 | 回复: 2 | |||
[交流]
计算出来的石墨烯弹性模量比理论值偏小 已有1人参与
|
| 大家好，我使用velocity来拉伸石墨烯，然后根据应力应变曲线计算出来的弹性模量只有150GPa左右，请问大家遇到过这种问题吗？ |
回复此楼» 猜你喜欢
退学或坚持读
已经有25人回复
-
遇见不省心的家人很难过
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有4人回复
-
免疫学博士有名额，速联系
已经有14人回复
-
面上基金申报没有其他的参与者成吗
已经有4人回复
-
多组分精馏求助
已经有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