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
返回小木虫查看更多
京公网安备 11010802022153号
出现断裂了。但是最大应力为180多GPa,感觉跟理论区别还是有点大的，而且断裂出现在左右固定端与移动的交界点。不知道应不应该是这样的
[图片]楼主应力是怎样转换的，我拉伸得到的force单位是kcal/mole*A怎样转换成牛顿�
，