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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 石墨烯拉伸破坏
51/1返回列表
查看: 2992  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

charliecaili

铁虫 (初入文坛)

[求助] 石墨烯拉伸破坏 已有2人参与

固定石墨烯y轴上下边界一定的距离，沿Y方向拉伸，可是模型会在上下两端破坏，怎么会这样，程序哪里出来问题，求大神支招
## INITIALIZATION ######################
units metal
dimension 3
boundary p p f
atom_style atomic
newton on

##  ATOM DEFINITION  ######################
#read_data data.graphene

lattice  custom 2.4595 a1 1 0 0 a2 0 1.73203 0 a3 0 0 1  basis 0 0.666666666 0 basis 0.5 0.5 0 basis 0.5 0.16666666 0 &
basis 0 0 0

region box block -8 8 -8 8 -50 50
create_box         1 box
region mid block -8 8 -8 8 -0.5 0.5
create_atoms        1 region mid

mass        1  12.010000

##  SETTINGS  ######################
pair_style airebo 3.0
pair_coeff * * potentials/CH.airebo C

region lower block INF INF INF -6 INF INF units box
region upper block INF INF 6 INF INF INF units box
group upper region upper
group lower region lower
group boundary union upper lower
#group mobile subtract all boundary

# initial velocities
fix 1 all nvt temp 300 300 0.1
velocity boundary set 0.0 0.0 0.0 units box
fix 2 boundary setforce 0.0 0.0 0.0

min_style cg
minimize 1e-25 1e-25 5000 10000

unfix 1
unfix 2

timestep 0.001
fix 1 all npt temp 300 300 0.05 y 0 0 0.5
thermo 20

compute    myTemp all temp
compute 1 all stress/atom myTemp
compute 2 all pe/atom pair bond
compute 3 all reduce sum c_1[1] c_1[2] c_1[3]

shell mkdir 10.9dujjk
shell cd 10.9dujjk

thermo_style custom step temp pe ke etotal lx ly press pxx pyy c_3[1] c_3[2] c_3[3]
dump 1 all custom 1000 dump.chiyu.* id type x y z
##  RUNING-RELAXATION  ######################
run 5000

##  RUNNING-DEFORMATION  ######################
unfix 1
undump 1
fix 1 all nvt temp 300 300 0.05

fix 2 all ave/atom 1 1000 1000 c_1[1] c_1[2] c_1[3] c_2 fx fy fz
dump 1 all custom 1000 dump.new.* id type x y z vx vy vz c_1[1] c_1[2] c_1[3] &
c_2 f_2[1] f_2[2] f_2[3] f_2[4] f_2[5] f_2[6] f_2[7]

variable srate equal 1.0e10
variable srate1 equal "v_srate / 1.0e12"
fix 3 all deform 1 y erate ${srate1} units box remap x
###################
variable tmp equal "ly"
variable ly0 equal ${tmp}
print "Initial Length, ly0: ${ly0}"

variable strain equal "(ly- v_ly0)/v_ly0"

# Output strain and stress info to file
# for units metal, pressure is in [bars] = 100 [kPa] = 1/10000 [GPa]
# p2 are in GPa
variable p1 equal "v_strain"
variable p2 equal "-pyy/10000"
variable p3 equal "pe"

fix equill all print 500 "${p1} ${p2} ${p3}" file data.1.txt screen no
run 50000

###########
######################################
# SIMULATION DONE
print "All done"
回复此楼

» 收录本帖的淘帖专辑推荐

分子动力学模拟

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
我就是我
1楼 2015-10-06 21:25:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

charliecaili

铁虫 (初入文坛)
送红花一朵
引用回帖:
6楼: Originally posted by 岩土小辈 at 2015-10-07 18:05:45
我看了一下，你那个basis是不是弄错了？你进行拉伸是算什么？...

basis 是建模，固定Y轴两端的原子相当于是我们做拉伸实验的夹具位置
赞一下 回复此楼
我就是我
8楼2015-10-08 09:48:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

岩土小辈

新虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
你既然已经固定了y方向边界，为什么还在y方向拉伸？
赞一下 回复此楼
2楼2015-10-07 10:04:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

charliecaili

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by 岩土小辈 at 2015-10-07 10:04:45
你既然已经固定了y方向边界，为什么还在y方向拉伸？

我是想把原子固定，不让它在两端断裂，但是还是在Y轴的两端断裂
赞一下 回复此楼
我就是我
3楼2015-10-07 13:06:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

岩土小辈

新虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by charliecaili at 2015-10-07 13:06:06
我是想把原子固定，不让它在两端断裂，但是还是在Y轴的两端断裂...

关键都已经把两端的原子固定了，你还怎么进行拉伸？
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

charliecaili
4楼2015-10-07 17:04:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯，如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗？ 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后，还能找到好工作吗？ +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭