Ni_Graphene复合材料纳米压痕模拟

#模型初始条件设置

units                metal

dimension            3

boundary            p p p

timestep             0.001

#生成模拟盒子

region                box block 0 40 0 40 0 40 units box

create_box            2 box

#生成Ni单晶

lattice               fcc 3.52

region                ni block 5 35 5 35 5 25 units box

create_atoms          1 region ni

#删除一层Ni原子

region                del block 0 25 0 25 18 20 units box

delete_atoms          region del compress yes

#自定义石墨烯晶格

lattice custom 2.4768 a1 1.0 0.0 0.0 a2 0.0 1.732 0.0 a3 0.0 0.0 1.3727 &

basis 0.0 0.33333 0.0 &

basis 0.0 0.66667 0.0 &

basis 0.5 0.16667 0.0 &

basis 0.5 0.83333 0.0

#定义石墨烯区域，生成石墨烯

region                 graphene block 5 35 5 35 18 21 units box

create_atoms        2 region graphene

#设置原子质量

mass                1 58.7

mass                2 12.0

#定义原子组，便于后期使用fix命令

region               lower block 0 35 0 35 5 8 units box

group                lower region lower

group                mobile subtract all lower

pair_style hybrid lj/cut 10 eam/alloy airebo 3.0

pair_coeff * * eam/alloy Ni99.eam.alloy Ni NULL

pair_coeff * * airebo CH.airebo NULL C

pair_coeff 1 2 lj/cut 0.023 2.852

#定义压头位置

variable   z equal "30 - 1*elaplong*dt"

velocity   mobile create 300 482748 temp new

fix       1 lower setforce 0.0 0.0 0.0

fix        2 all nvt temp 300 300 0.01

fix        3 all indent 10.0 sphere 20 20 v_z 5.0 units box

#输出轨迹

dump       1 all atom 1000 ni.lammpstrj

run        10000