lammps模拟喷丸过程

大家好，我是小马老师。
本文介绍一个lammmps模拟喷丸过程的案例。
lammps模拟喷丸过程的方法比较多，比如速度冲击法、fix deposit沉积法等。
原理基本类似，给球体一个特定的速度撞击基体。
本案例采用velocity速度冲击法，金刚石小球高速撞击Si基体，撞击后在Si表面留下一个撞击坑。
模拟动画如下图：

案例代码比较简单，抛砖引玉给大家提供一个思路，省略了弛豫过程以及撞击后的数据采集过程，有兴趣的可以添加撞击之后基体的应力分布、温度分布、球体的能量损失等等。
in文件代码：

units metal
atom_style atomic
boundary p p f
neighbor 2 bin
neigh_modify every 1 delay 0
timestep 0.001

lattice diamond 5.43
region box block 0 15 0 15