应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 石墨烯单轴拉伸模拟
31/1返回列表
查看: 3080  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

应力-应变曲线
回复此楼

» 猜你喜欢

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

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

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

查看全部散金贴
1楼 2012-06-12 19:08:40
已阅   回复此楼   关注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    这些数字 代表什么？
谢谢回答
赞一下(1人) 回复此楼
3楼2017-05-10 09:19:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuzhihua1988 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] E0414, 我的本子有没有希望？ +7 布布和一二 2026-06-17 9/450 2026-06-21 19:29 by 布布和一二
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 19:18 by q1p9ixelzp
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:58 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:00 by q1p9ixelzp
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 4/200 2026-06-21 17:45 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:40 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:37 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 4/200 2026-06-21 17:25 by ky8v3skvft
[考博] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:20 by ky8v3skvft
[博后之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 6/300 2026-06-21 17:05 by ky8v3skvft
[公派出国] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:00 by ky8v3skvft
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +5 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:57 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 vlfdw50o0u 2026-06-20 5/250 2026-06-21 16:45 by ky8v3skvft
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:40 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 9/450 2026-06-21 16:37 by ky8v3skvft
[基金申请] 青A35岁以下通知答辩了吗 +3 暨阳一只柴 2026-06-17 3/150 2026-06-21 09:39 by kudofaye
[考博] 【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院 +3 robohaha 2026-06-15 4/200 2026-06-21 00:36 by lanbing1
[硕博家园] 一篇论文同时出现在两个期刊，一模一样，这算不算学术不端，请各位老师斧正。 +11 mahl005 2026-06-14 12/600 2026-06-18 17:01 by 770904134
[基金申请] 希望面上有个好结果 +7 碧水00 2026-06-16 7/350 2026-06-18 12:18 by wuke100666
[论文投稿] 三区计算机方向期刊推荐 +5 1457340941 2026-06-15 5/250 2026-06-17 13:16 by 会议编辑
信息提示
关闭
请填处理意见
关闭