[lammps教程] 合金拉升模拟

原创 YJ学长 LAMMPS交流站javascript:void(0); 大家好，本文是小编根据自己所学的知识，因能力有限，将介绍一种高熵合金的lammps简单拉伸过程，本教程只适合于新手学习，希望能给读者提供具体工程意义还读者结合实际同时希望专业前辈提出不足。

    如下图为本文建立的合金试件模型。

01

基本思路

本文采用deform命令对高熵合金进行一个拉伸操作

02

完整代码

完整in文件如下：

#初始模拟系统设置units metal #单位boundary p p p #边界条件atom_style atomic #原子类型timestep 0.001 #时间步长#构建模型lattice fcc 3.56 #设置晶格参数region box block 0 100 0 20 0 20 #生产长方形区域create_box 5 boxcreate_atoms 1 box#将type为1的原子转换为其他原子set type 1 type/ratio 2 0.2 87393 #将1原子转换为2原子set type 1 type/ratio 3 0.5 87393 #将1原子转换为3原子set type 1 type/ratio 4 0.5 87393 #将1原子转换为4原子set type 3 type/ratio 5 0.5 87393 #将1原子转换为5原子#设置原子类型mass 1 27 #Al原子mass 2 56 #Fe原子mass 3 28 #Si原子mass 4 64 #Cu原子mass 5 24 #Mg原子#拉伸设置fix 1 all deform 2000 x erate 0.05 remap x  #设置拉伸方式dump 1 all atom 500 lasheng.xyz  #输出文件run 100000

03

运行结果

    接着，在cmd窗口运行in文件，将生成的xyz文件拖到OVITO软件中，并在OVITO软件中进行渲染，就产生了如下结果：

接着

在加上使用输出的应力应变数据（本文中in文件未给出）在origin中绘制应力-应变曲线，如下图所示：

更多案例尽在公众号：LAMMPS交流站