24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3041  |  回复: 2

wuzhihua1988

新虫 (初入文坛)


[交流] 石墨烯单轴拉伸模拟

这是用airebo势函数,300K条件下的拉伸试验,下面附带有IN文件以及应力-应变曲线图。个人觉得问题有:1.驰豫过程中,温度是从300K开始升到两千多,然后快到驰豫结束时才又降为300K,不知道这个现象是否正常?;2 应力-应变曲线起伏波动比较大,没有一般文献的那么标准,不知道是否哪一快设置没有好?;3。拉伸断裂用VMD观察时,没有缓慢断裂而是比较突然,从应力应变曲线中可以相应看出。且断裂后,Y方向的长度是非常明显的变化,不知道是不是因为300K这个温度,使原子断裂后受热的扰动影响作用而导致。4.应变为零时,应力没法从零开始,大概是20GPA左右,不知道怎么可以能够改进(PS:平衡到最后,应力值就是保持在20几GPA)
希望各位能够不吝赐教。
dimension 3
> boundary p s s
>
> units metal
> atom_style atomic
> newton on
>
> neighbor 2.0 nsq
> neigh_modify every 10 check yes
>
> read_data data.graphene
>
> # potentials
> pair_style airebo
2.0 0 0
> pair_coeff * * CH.airebo C
>
> # groups
> region Ledge block INF 3 INF INF INF INF units box
> region Redge block 69 INF INF INF INF INF units box
>
> group left region Ledge
> group right region Redge
> group mobile subtract all left right
>
> thermo_style custom step lx ly lz press pxx pyy pzz pe temp
> thermo 1000
>
> # initial velocities
> velocity        left zero linear
> velocity        right zero linear
> velocity mobile create 300.0 5481314 dist gaussian mom yes rot yes loop
> local units box
>
> # nve equilibrate
> dump 1 all atom 100000 grain.lammpstrj
>
> fix 0 right store/force
> fix 1 left setforce 0.0 NULL 0.0
> fix 2 right setforce 0.0 NULL 0.0
> fix 3 mobile nvt temp 3000.0 300.0 0.1
>
> timestep 0.001
>
> min_style cg
> minimize
1.0e-15 1.0e-5 1000 1000
>
> run 1000
> unfix 3
> undump 1
>
> # fracture
> #fix           5 all temp/berendsen 300.0 300.0 0.1
> #fix           3 all nve
> fix 3 all nvt temp 300.0 300.0 0.1
> run_style verlet
>
> dump 1 all atom 10000 grain.lammpstrj
> dump_modify 1 append yes first yes sort id
>
> label loading
> variable flag loop 150
>

fix 6 all deform 1000 x erate 0.004 units box
> run 1000
> unfix         6
>
> next flag
> jump in.graphene loading

应力-应变曲线
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

韩欢欢

金虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,请问你的问题解决了没啊?我遇到了类似的问题,从2000k降到300k,然后在300k恒温拉伸,但是在300k时开始拉伸时应变为0,应力不为0,是不是降温中产生了内应力,这个怎么解决,求教
2楼2013-08-12 13:51:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

honest199174

铜虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主问题解决了吗?
我想用lammps 模拟石墨烯拉伸断裂 用airebo势函数,想把robo项设置为2  LJ项设置为6.8,请问该怎么设置。谢谢。  另外pair_style      airebo 2.0  0 0 和pair_style airebo 3.0 1 0    这些数字 代表什么?
谢谢回答
3楼2017-05-10 09:19:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuzhihua1988 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 273求调剂 +4 白居不易. 2026-04-09 6/300 2026-04-14 19:54 by 逆水乘风
[考研] 求助调剂,跨调 +18 X十甫寸Y 2026-04-11 19/950 2026-04-14 19:26 by Art1977
[考研] 08工学 309分求调剂 +12 Yin DY 2026-04-08 12/600 2026-04-14 17:49 by lhj2009
[考研] 297工科调剂? +13 河南农业大学-能 2026-04-13 13/650 2026-04-14 16:46 by Art1977
[考研] 296求调剂 +8 汪!?! 2026-04-09 9/450 2026-04-14 15:15 by 逆水乘风
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 314求调剂 +24 wakeluofu 2026-04-09 25/1250 2026-04-13 08:58 by lhj2009
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
信息提示
请填处理意见