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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 石墨烯单轴拉伸模拟
31/1返回列表
查看: 3037  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工349分求调剂 +9 李木子啊哈哈 2026-04-12 10/500 2026-04-13 08:45 by Sammy2
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 22专硕求调剂 +6 haoyun上岸 2026-04-11 8/400 2026-04-11 23:21 by labixiaoqiao
[教师之家] 请问地理、遥感方面，可以做哪些横向项目啊，纵向完不成考核啊 +3 锦衣卫寒战 2026-04-07 5/250 2026-04-11 20:51 by 豫椒
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 一志愿北理工298英一数二已上岸，感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[考研] 机械还有还有名额吗？太难了 +6 笑笑袁 2026-04-10 6/300 2026-04-10 11:54 by 高维春
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 材料化工总分334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 20:19 by maddjdld
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
[考研] 338求调剂 +5 小猪红色 678 2026-04-06 6/300 2026-04-07 21:18 by 乔哒哒哒
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
信息提示
关闭
请填处理意见
关闭