[lammps第十讲] CuNi合金熔化模拟L

公众号：LAMMPS交流站javascript:void(0);

大家好，本文是小编根据自己所学的知识，因能力有限，将介绍CuNi合金的lammps熔化过程，本教程适合于新手学习，希望能给读者提供具体工程意义还读者结合实际同时希望专业前辈提出不足。

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

01

基本思路

首先建立CuNi合金模型

在npt系综下1550K保温

在npt系综下1550K升温至2550K

在npt系综下2550K降温至2550K

02

完整代码

完整in文件如下：

units           metalboundary        p p patom_style      atomicread_data  melting.in.datlattice         fcc 3.65 region          box block -30 30  -5 5 -5 5#create_box      1 box#create_atoms  mass            1  64.55mass            2  58.69 pair_style      eam/alloypair_coeff      * * CuNi.lammps.eam Cu Niregion          1 block INF -0.01  -6 6  -6 6group           left region 1region          2 block -0.01 INF  -6 6 -6 6group           right region 2timestep        0.002velocity        all create 1550.0 4928459 dist gaussianfix             1 all npt temp 1550 1550 0.1 iso 0 0 0.2dump            melt all atom 100 dump.dat#dump            melt all custom 1000 dump.atom.slow tag type xu yu zu vx vy vz  sxx syy szzrestart         50000  restart.*thermo          100run             4000#thermo_style    custom step temp etotal pe press vol pxx pyy pzz pyz pxz pxyunfix   1fix             1 right npt temp 1550 2500 0.1 x 0 0   0.2run             2000unfix   1fix             1 right npt temp 2500 1550 0.1 x 0 0   0.2run             44000unfix           1fix 1     all nph iso 0.0 0.0 0.2  run             100000

03

运行结果

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

本文中所用到的data文件及势文件未给出，请关注公众号获取