24小时热门版块排行榜    

CyRhmU.jpeg
查看: 735  |  回复: 2

汤圆喵

银虫 (小有名气)

[交流] 计算出来的石墨烯弹性模量比理论值偏小已有1人参与

大家好,我使用velocity来拉伸石墨烯,然后根据应力应变曲线计算出来的弹性模量只有150GPa左右,请问大家遇到过这种问题吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

随世流逝

新虫 (小有名气)

2楼2021-04-29 11:25:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

汤圆喵

银虫 (小有名气)

引用回帖:
2楼: originally posted by 随世流逝 at 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 汤圆喵 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见