24小时热门版块排行榜    

查看: 3485  |  回复: 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努力

新虫 (初入文坛)

引用回帖:
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的回帖
查看全部 4 个回答

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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +3 负心者当诛 2026-03-07 3/150 2026-03-14 06:15 by JourneyLucky
[考研] 346分材料求调剂 +3 snow_反季节版 2026-03-07 3/150 2026-03-14 05:55 by JourneyLucky
[考研] 材料工程085601调剂求老师收留 +12 强木木木 2026-03-07 14/700 2026-03-14 05:37 by JourneyLucky
[考研] 一志愿天大化工(085600)调剂总分338 +6 蔡大美女 2026-03-09 6/300 2026-03-14 02:46 by JourneyLucky
[考研] 考研材料与化工,求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[文学芳草园] 伙伴们,祝我生日快乐吧 +12 myrtle 2026-03-10 17/850 2026-03-13 22:00 by angelyueyi
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 求0856调剂 +5 小力气珂珂 2026-03-08 5/250 2026-03-09 06:42 by houyaoxu
[考研] 求调剂,一志愿江南大学,食品科学与工程,总分,320 +3 yyyyyukino 2026-03-07 3/150 2026-03-08 23:07 by 清风月
[教师之家] 交大前校长王树国:现在最先进的科技并不在大学实验室,而是在企业研究院 +4 zju2000 2026-03-08 6/300 2026-03-08 19:15 by zju2000
信息提示
请填处理意见