| 查看: 3478 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者lupei努力将赠送您 5 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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 |
回复此楼» 猜你喜欢
288求调剂
已经有7人回复
-
310 070300化学求调剂
已经有5人回复
-
A区一本交叉课题组，低分调剂，招收机械电子信息通信等交叉方向
已经有67人回复
-
320材料与化工，求调剂
已经有13人回复
-
新疆大学地质与矿业工程学院招生
已经有27人回复
-
找博导
已经有3人回复
-
一志愿南大化学339分求调剂，四六级已过，有比赛，有文章
已经有7人回复
-
0703求调剂
已经有11人回复
-
327分求调剂086
已经有6人回复
-
复试调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
yuedaixz
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 154
- 帖子: 123
- 在线: 15.4小时
- 虫号: 19796653
- 注册: 2019-11-18
- 性别: GG
- 专业: 其他无机非金属材料
3楼2020-06-28 21:01:03
2楼2020-06-24 15:44:08
4楼2020-07-10 15:56:16
100