24小时热门版块排行榜    

查看: 3476  |  回复: 3
【悬赏金币】回答本帖问题,作者lupei努力将赠送您 5 个金币

lupei努力

新虫 (初入文坛)

[求助] lammps使用velocity方法进行石墨烯拉伸

最近使用lammps做石墨烯之字形方向的拉伸模拟,按照网上各种资源设计了自己的in文件,一直出现拉不断的情况,边界原子键角已经很平了,可是就是不断裂。 应力已经到了170gpa,应变到了0.3还没断。我看很多文献之字形的理论抗拉强度是120,GPA,130GPA等等。不知道什么地方出现问题。下面是我的in文件
boundary  s s s                             
units       metal                             
atom_style  full                             
read_data  perfect.data
region          left block INF 5 INF INF INF INF
region          right block 93 INF INF INF INF INF
group           left region left
group           right region right
group           boundary union left right
group           mobile subtract all boundary

# potentials
pair_style airebo 1.92 1 1
pair_coeff * * CH.airebo C
neighbor        2.0 nsq
neigh_modify every 1 delay 0 check yes
# define computes
compute    myTemp all temp
compute    disp all displace/atom
compute    pot all pe/atom
compute    ke all ke/atom
compute    s1 all stress/atom NULL
compute    press  all pressure thermo_temp
compute    lx  all reduce sum c_s1[1]
variable  tmp  equal "lx"
variable  L0   equal  ${tmp}
variable  strain  equal  "(lx-v_L0)/v_L0"               
thermo_style custom step temp  c_myTemp etotal ke pe lz ly lx vol press pxx pyy pzz v_tmp v_L0 v_strain
thermo 100
thermo_modify    lost  ignore                                                     
timestep    0.001   
min_style sd
minimize 1.0e-6 1.0e-6 1000 1000            
#relax  for equilibrium
velocity all  create 300 4928459 mom yes rot no
fix 1 all  nve
fix 2 all  langevin 300.0 300.0 0.08 48279
run  40000                                                                                                   
dump                5 all custom 1000 atom.lammpstrj id type x y z vx vy vz xs ys zs c_pot  c_ke  c_disp[1] c_disp[2]
c_disp[3] c_s1[1] c_s1[2] c_s1[3] c_s1[4] c_s1[5] c_s1[6]
#set the deformation process
velocity        boundary  set 0.0 0.0 0.0
fix             3 boundary setforce 0.0 0.0 0.0
velocity        right set 0.0995 0.0 0.0  sum yes
velocity        mobile ramp vx 0.0 0.0995 x 5 93 sum yes                                   
run    2000000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lupei努力

新虫 (初入文坛)

出现断裂了。但是最大应力为180多GPa,感觉跟理论区别还是有点大的,而且断裂出现在左右固定端与移动的交界点。不知道应不应该是这样的
2楼2020-06-24 15:44:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuedaixz

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by lupei努力 at 2020-06-24 15:44:08
出现断裂了。但是最大应力为180多GPa,感觉跟理论区别还是有点大的,而且断裂出现在左右固定端与移动的交界点。不知道应不应该是这样的

[图片]楼主应力是怎样转换的,我拉伸得到的force单位是kcal/mole*A怎样转换成牛顿啊
3楼2020-06-28 21:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lupei努力

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by yuedaixz at 2020-06-28 21:01:03
楼主应力是怎样转换的,我拉伸得到的force单位是kcal/mole*A怎样转换成牛顿啊...

不好意思,我是直接采集应力的。

发自小木虫Android客户端
4楼2020-07-10 15:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lupei努力 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +3 上班不着吉 2026-03-09 3/150 2026-03-09 16:52 by 紫曦紫棋
[考研] 0832食品科学与工程293调剂 +3 东东不東 2026-03-07 3/150 2026-03-09 15:58 by ls_13456
[考研] 化工284求调剂 +8 小黑想变白 2026-03-07 8/400 2026-03-09 14:35 by 30660438
[考博] 求材料读博院校 +8 yanglei131 2026-03-08 8/400 2026-03-09 10:54 by 馥安馥安
[考博] 申博 +5 添菜了哈 2026-03-04 8/400 2026-03-09 10:52 by jamesqi
[考研] 一志愿中科大材料与化工。353分能调剂到哪? +3 MayUxw1 2026-03-08 3/150 2026-03-09 10:13 by EBSD
[考研] 0856求调剂 +3 squirtle11 2026-03-07 3/150 2026-03-09 09:54 by @飒飒飒飒
[考研] 290 材料与化工求调剂 +7 Nebulala 2026-03-08 7/350 2026-03-09 05:21 by houyaoxu
[考研] 招调剂 +20 帆船哥 2026-03-04 23/1150 2026-03-08 16:25 by taycat
[考研] 308求调剂 +7 倘若起风了呢 2026-03-05 9/450 2026-03-06 22:59 by 倘若起风了呢
[考研] 085600,一志愿郑州大学,280分求调剂 +7 Wuqi725 2026-03-05 7/350 2026-03-06 08:12 by 杨杨杨紫
[考研] 328求调剂 +6 vuzhdkfjkx 2026-03-04 7/350 2026-03-05 23:28 by Leeding1356
[考研] 材料085601一志愿哈工大317 +4 压迫感行 2026-03-04 4/200 2026-03-05 20:08 by 黑衣馒头人
[考研] 材料328求调剂 +10 一个萝卜02 2026-03-03 10/500 2026-03-05 09:31 by oxidpl
[考研] 304分材料专硕求调剂 +11 qiuzhigril 2026-03-03 14/700 2026-03-05 09:30 by vfdbghn
[考研] 0856材料专硕274能调剂去哪里? +3 22735 2026-03-04 4/200 2026-03-05 09:06 by 斩魂滴兔子!
[考研] 22408-273求调剂-担任3个项目的负责人-1篇国际期刊论文(一作)1篇核心期刊论文在投。 +3 没想好取什么名 2026-03-03 3/150 2026-03-05 08:38 by Destiny688
[考研] 085601 材料305分求助 +4 泡泡邮件 2026-03-03 6/300 2026-03-04 13:00 by 无际的草原
[考研] 331求调剂 +3 zzZ&zZ 2026-03-03 3/150 2026-03-04 10:37 by aaa显卡批发王总
[考研] 化工335求调剂 +5 摸摸猫猫头 2026-03-02 5/250 2026-03-03 09:10 by houyaoxu
信息提示
请填处理意见