| 查看: 814 | 回复: 2 | |||
[交流]
计算出来的石墨烯弹性模量比理论值偏小 已有1人参与
|
| 大家好，我使用velocity来拉伸石墨烯，然后根据应力应变曲线计算出来的弹性模量只有150GPa左右，请问大家遇到过这种问题吗？ |
回复此楼» 猜你喜欢
277求调剂 数一104分
已经有7人回复
-
288环境专硕,求调材料方向
已经有16人回复
-
277工科求调剂
已经有8人回复
-
266求调剂
已经有19人回复
-
环境专硕调剂
已经有3人回复
-
081200-11408-276学硕求调剂
已经有5人回复
-
材料调剂
已经有13人回复
-
工科 22408 267求推荐
已经有4人回复
-
305求调剂
已经有3人回复
-
270分求调剂
已经有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
5