| 查看: 876 | 回复: 2 | |||
[交流]
计算出来的石墨烯弹性模量比理论值偏小 已有1人参与
|
| 大家好，我使用velocity来拉伸石墨烯，然后根据应力应变曲线计算出来的弹性模量只有150GPa左右，请问大家遇到过这种问题吗？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
40