24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3038  |  回复: 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的回帖

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的回帖
查看全部 3 个回答

韩欢欢

金虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,请问你的问题解决了没啊?我遇到了类似的问题,从2000k降到300k,然后在300k恒温拉伸,但是在300k时开始拉伸时应变为0,应力不为0,是不是降温中产生了内应力,这个怎么解决,求教
2楼2013-08-12 13:51:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂求收留 +30 果然有我 2026-04-10 31/1550 2026-04-13 14:18 by 张zhihao
[考研] 本科生物信息学,总分362 求07 08调剂 +8 q小倩1210 2026-04-06 8/400 2026-04-13 11:58 by 414435937
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] 0854调剂 +10 长弓傲 2026-04-12 12/600 2026-04-13 08:49 by Sammy2
[考研] +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 283求调剂 086004考英二数二 +17 那个噜子 2026-04-10 18/900 2026-04-11 16:27 by 明月此时有
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 273求调剂 +51 麦小叮当 2026-04-06 58/2900 2026-04-10 15:54 by jiajinhpu
[考研] 一志愿华东师范生物学326分,求调剂 +8 刘墨墨 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 机械工程313分找工科调剂 +3 双一流本科机械 2026-04-08 3/150 2026-04-08 20:41 by 土木硕士招生
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
信息提示
请填处理意见