24小时热门版块排行榜    

查看: 1076  |  回复: 0

hazheng

捐助贵宾 (著名写手)


[交流] lammps案例分析(1):石墨烯单轴拉伸之velocity方式

lammps模拟石墨烯拉伸过程有两个命令可选:deform和velocity,两个命令的原理不同。

deform是按照一定的速率拉伸box,在这个过程中,石墨烯的原子随着box的伸长而变化坐标,从而实现整体的拉伸。

velocity方式通常固定石墨烯的一端,给另一端一个固定的速度v,由这一端的原子带动其余原子运动。



本文采用velocity方式进行石墨烯的拉伸,下一篇文章将采用deform方式对石墨烯进行拉伸。

下面对石墨烯拉伸in文件进行详细分析。

#模型初始化,模拟单位为metal,原子方式为atomic,3维模型
units  metal
atom_style   atomic
dimension  3
#定义边界条件,x和y方向为周期性边界,y为拉伸方向,z方向为固定边界

boundary     p p f
#定义邻域原子定义,1个原子的邻域原子列表是以该原子为中心,以力截断半径+缓冲值(也就是下面设置的0.3)为半径的圆内原子列表

neior      0.3 bin
#设置模拟步长为1fs

timestep        0.001
#生成模拟box,box内含有3种原子,这三种原子都是c原子,为了后期线上着色方便,原子类型分为上中下三种

region   box block 0 175 0 175 -5 5 units box
create_box 3 box
#自定义石墨烯晶格参数

lattice  custom 2.4768 a1 1.0 0.0 0.0 &
a2 0.0 1.732 0.0 a3 0.0 0.0 1.3727 &
basis 0.0 0.33333 0.0 &
basis 0.0 0.66667 0.0 &
basis 0.5 0.16667 0.0 &
basis 0.5 0.83333 0.0
#生成石墨烯原子,类型为1

region  gp block 10 120 10 120 -0.5 0.5 units box
create_atoms 1 region gp
#设置原子质量

mass * 12
#定义区域,上下两部分固定区域和中间活动区域mobile

region upper block inf inf 117 inf inf inf  units box
region lower block inf inf inf 13 inf inf units box
group upper region upper
group lower region lower
group boundary union upper lower
group mobile subtract all boundary
#设定最上端原子类型为2,最下端原子类型为3,为后期着色方便

set group upper type 2
set group lower type 3
#设置势函数,airebo常用于石墨烯的模拟

pair_style  airebo 2.0
pair_coeff  * * ch.airebo c c c
#初始化温度

velocity mobile create 300.0 8877528
#固定上下两端原子

fix  1 boundary setforce 0.0 0.0 0.0
#以下四句代码,对模型在npt下进行弛豫,并将结果保存到文件中

thermo 100
fix  2 all npt temp 300.0 300.0 1 x 0 0 0.1 y 0 0 0.1
dump 1 all atom 100 gp_relax.lammpstrj
run 1000
#取消fix,dump设定,步数清零

unfix   2
undump 1
reset_timestep 0
#最上端设置y方向速度为1,沿y方向拉伸

velocity upper set 0.0 1.0 0.0
#按比例设置mobile部分的速度

velocity mobile ramp vy 0.0 1.0 y 8 52 sum yes
#在nvt下进行拉伸

fix  2 all nvt temp 300.0 300.0 0.01
#进行热力学输出,保存拉伸后原子坐标

thermo 1000
thermo_modify lost ignore
dump  1 all atom 1000 gp_tension.lammpstrj
run  10000
模拟最终结果如下图所示,本例比较简单,没有加入应力应变的统计,如要输出应力和应变,可用fix或者thermo语句输出pxx、pyy、pzz三个方向上的应力。

1.jpg



最后,欢迎大家通过微信公众号联系我们。

微信公众号:320科技工作室。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hazheng 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 最新消息:2024国自然下载文件名变了 +9 dream200606 2024-06-16 12/600 2024-06-17 00:34 by kobe0107
[硕博家园] 硕博巨婴,也许才刚刚开始 +30 SNaiL1995 2024-06-12 77/3850 2024-06-16 23:18 by 欢喜1998
[教师之家] 请问事业编制和年薪制冲突吗? +8 ZHONGWU_U 2024-06-14 8/400 2024-06-16 17:37 by spj860iuy
[找工作] 初始合伙人来啦!(生物试剂耗材标准品) +8 欢快的小科研人 2024-06-15 14/700 2024-06-16 16:00 by 欢快的小科研人
[基金申请] 面青地会评时间 +6 tanjydd 2024-06-15 6/300 2024-06-16 15:42 by 我是王小帅
[教师之家] 每次骚扰女学生的都是院系领导,而不是普通教师,小编们要注意措辞正确 +9 zju2000 2024-06-15 11/550 2024-06-16 14:49 by appleapple2
[基金申请] 博后基金,博管会会提前知道消息吗? +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[找工作] 成都产品质量检测研究院 200+3 鲸鱼663 2024-06-11 9/450 2024-06-16 10:08 by SNaiL1995
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
[论文投稿] 审稿人含糊拒稿,还需要回复吗?如何回复? 20+3 BruceChum 2024-06-15 17/850 2024-06-15 20:19 by arthas_007
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[基金申请] 有没有机械的前辈分享一下评上海优都是什么成果啊 +7 wulala800 2024-06-10 7/350 2024-06-15 09:33 by 晓目崇
[论文投稿] 投稿kbs被拒后,系统提供的推荐的其他期刊有用吗? +3 chenrui2015 2024-06-10 3/150 2024-06-14 10:30 by Tr.sjx1997
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了,能确定送审了吗? +5 潇洒怡惜 2024-06-12 10/500 2024-06-14 09:56 by 潇洒怡惜
[考博] 博导选择 +3 bing85977 2024-06-12 3/150 2024-06-13 15:34 by 我是邱尧
[有机交流] 原料反应完了,怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[硕博家园] 考博应该从哪里开始准备 +6 亮姐2019 2024-06-11 8/400 2024-06-13 12:53 by 帕秋莉GO
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
[基金申请] 博后特助这周出结果吗?往年都是啥时候啊? +13 jsqy 2024-06-12 17/850 2024-06-12 19:55 by Lynn212
信息提示
请填处理意见