分子动力学方法研究刃位错.pdf

分子动力学方法研究刃位错 1.LAMMPS 中晶体建模 2. 如何引入刃位错 3.通过加载使位错移动 LAMMPS 中晶体建模 units metal #定义单位系统 /doc/units.html dimension 3 #三个方向均为周期 性边界条件 boundary p p p atom_style atomic #定义颗粒类型 LAMMPS 中晶体建模 lattice bcc 2.8553 orient x 1 -1 1 orient y 2 1 -1 orient z 0 1 1 region box block 0 2 0 3 0 2 #定义晶格类型为体心立方，晶格常数为2.8553埃 #改变晶向，将晶格的（1 -1 1）方向指向盒子的x方向， （2 1 -1）方向指向盒子的y方向，（0 1 1）方向指向盒子的z 方向 #以上述晶格类型为基准单位，定义一个长方体区域，x方向范 围是0到2个基准,y方向上0到3个基准，z方向上0到2个基准 LAMMPS 中晶体建模 晶格取向： 1 -1 1 0 1 1 -1 0 0 -2 -1 1 x-y 平面 y-z 平面 LAMMPS 中晶体建模 create_box 1box #制造一个盒子，原子种类为1 create_atoms 1region box #在前面指定的区域内按前面 指定的晶体类型和晶向创造原子 replicate 30 6 1 #对样品进行复制，x,y,z方向 上各扩大30倍，6倍和1倍 LAMMPS 中晶体建模 region right block 55 INF INF INF INF INF #定义一个区域right，长方体，x方向范围是55到无限远， INF是无限远，单位仍是前面定义的晶格基准 delete_atoms region right #删除该区域的原子 region left block INF 5 INF INF INF INF delete_atoms region l #前面4行指令删除了盒子左右两侧的部分原子，制造出两 个自由表面 引入刃位错 region a block INF INF 17.5 INF INF INF region b block INF INF INF 0.3 INF INF region c block 25 37.5 INF INF INF INF units box group a region a group b region b group c region c group boundary union a b c #指定边界原子，先定义区域，然后把这个区域内的原子放在一个组 内，units box意味着长度单位不是默认的晶格基准，而是埃Angstrom 引入刃位错 *.cfg A = 1 Angstrom (basic length-scale) H0(1,1) = 296.731 A H0(1,2) = 0 A H0(1,3) = 0 A H0(2,1) = 0 A H0(2,2) = 83.9283 A H0(2,3) = 0 A H0(3,1) = 0 A