24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者lupei努力将赠送您 5 个金币

lupei努力

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 27
帖子: 8
在线: 5.7小时
虫号: 20030296
注册: 2019-12-05
专业: 微/纳机械系统

[求助] 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

回复此楼