24小时热门版块排行榜    

查看: 3408  |  回复: 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的回帖

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

lupei努力

新虫 (初入文坛)

出现断裂了。但是最大应力为180多GPa,感觉跟理论区别还是有点大的,而且断裂出现在左右固定端与移动的交界点。不知道应不应该是这样的
2楼2020-06-24 15:44:08
已阅   回复此楼   关注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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见