lammps教程：FeNiCr合金多晶结构建模方法

大家好，我是小马老师。
本文介绍合金多晶模型的建模过程，使用的软件为Atomsk。
Atomsk是一个开源软件，可运行在windows、linux等多种操作系统上，本文使用windows版的Atomsk实现建模，其他操作系统的建模命令基本相同。
建模过程如下：
1. 生成单晶Fe晶胞

atomsk --create fcc 3.65 Fe Fe.xsf

2. 建立多晶节点文件polycrystal.txt
新建一个文本文件polycrystal.txt，文件内容为：

box 200 100 200
random	20

该文件设定盒子尺寸为200100200nm，盒子内随机生成20个晶粒。
3. 生成多晶文件final.lmp

atomsk --polycrystal Fe.xsf polycrystal.txt final.lmp -wrap

将Fe.xsf晶胞填充到20个晶粒空间内，最终多晶数据保存到文件final.lmp。
4. 修改final.lmp
使用文本编辑器打开final.lmp文件，做以下修改：
（1）原子类型由1种改为3种：

3 atom types

（2）添加Ni、Cr原子摩尔质量：

Masses
 
1   55.84500000             # Fe
2   58.69                   # Ni
3   51.96                   #Cr

5. 替换原子生成合金结构
编写in文件，在lammps中使用替换原子法，将部分Fe原子按照比例替换为Ni、Cr，得到合金多晶结构。
In文件代码如下：

units           metal
boundary        p p p
atom_style      atomic
timestep        0.001
neighbor        0.2 bin
read_data       final.lmp
set             type 1 type/ratio 2 0.33 8793
set             type 1 type/ratio 3 0.5 56332
write_data      Fe-Ni-Cr.data

---

更多lammps教程，关注微信公众号：
lammps加油站