[lammps案例]精确原子比高熵合金建模

在lammps中建立高熵合金模型是通过替换原子的方式实现的，其命令为set type/ratio, 其缺点为多次替换之后，各原子之间的比例将不准确。在此通过变脚本书写方式，建立准确原子比例的高熵合金。本文以FeNiCrCoCu高熵合金为例，如图1所示，其晶体结构为FCC，原子比例Fe:Ni:Cr:Co:Cu=25:15:20:25:15。

图1 FeNiCrCoAl模型

完整in文件如下：

# 模型参数设置units               metal
boundary            p p patom_style          atomic
variable lat equal 3.56
variable x_length equal 50
variable y_length equal 50
variable z_length equal 50
#计算总原子数variable N_atoms equal 4*v_x_length*v_y_length*v_z_length #在100 010 001 晶向中，晶胞原子数 FCC = 4 BCC = 2 variable N equal ${N_atoms} 
#计算2原子替换比例variable            2_ratio equal 0.15variable            2_replace equal ${2_ratio}#计算3原子替换比例variable            3_ratio equal 0.20variable            3_atoms equal v_N*v_3_ratio variable            1_atoms equal v_N-(v_N*v_2_ratio) variable            3_ratio1 equal v_3_atoms/v_1_atomsvariable            3_replace equal ${3_ratio1}#计算4原子替换比例variable            4_ratio equal 0.25variable            4_atoms equal v_N*v_4_ratiovariable            1_atoms equal v_N-(v_N*v_2_ratio)-(v_N*v_3_ratio)variable            4_ratio1 equal v_4_atoms/v_1_atomsvariable            4_replace equal ${4_ratio1}#计算5原子替换比例variable            5_ratio equal 0.15variable            5_atoms equal v_N*v_5_ratiovariable            1_atoms equal v_N-(v_N*v_2_ratio)-(v_N*v_3_ratio)-(v_N*v_4_ratio)variable            5_ratio1 equal v_5_atoms/v_1_atomsvariable            5_replace equal ${5_ratio1}#创建模型lattice             fcc ${lat} orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region              box block 0 ${x_length} 0 ${y_length} 0 ${z_length}create_box          5 boxcreate_atoms        1 region box set                 type 1 type/ratio 2 ${2_replace} 154648set                 type 1 type/ratio 3 ${3_replace} 121412set                 type 1 type/ratio 4 ${4_replace} 789546set                 type 1 type/ratio 5 ${5_replace} 564831mass                1 55.8mass                2 58.69mass                3 52mass                4 58.93mass                5 64write_data          hea.data

本期的分享到此结束，欢迎大家点赞、收藏和转发。如果觉得小编写的不错，不妨点个免费的赞。感谢大家！！！