大家好，我是小马老师。本文介绍lammps模拟冲击的一个案例。冲击是指一个物体以较快的速度冲向另一个物体，这个过程会有温度的升高，因此，在lammps冲击模拟中，冲击过程一般采用nve系综。冲击的模拟相对比较简单，在弛豫完成后，设置基体不动，冲击物体设置一个较高的速度冲向基体即可。
下面给出Al块的冲击案例代码以供参考，建模在lammps in文件中进行，不需要读取data文件。
#模拟基本参数
units metal
boundary s p p
atom_style atomic
timestep 0.001
neighbor 2 bin
neigh_modify every 1 delay 0
#建立盒子
lattice fcc 4.05
region box block 0 20 0 10 0 10
create_box 2 box
#生成左右两部分Al模型
region left block INF 3 INF INF INF INF
create_atoms 2 region left
region right block 8 INF INF INF INF INF
create_atoms 1 region right
#分组
group left type 2
group right type 1
#边界分组
region boundary block INF INF INF INF 19 INF
group boundary region boundary
#势函数设置
pair_style eam/fs
pair_coeff * * Al1.eam.fs Al Al
#能量最小化
minimize 1e-10 1e-10 10000 10000
#温度初始化
velocity all create 300 98788
#固定边界层
velocity boundary set 0 0 0
fix 01 boundary setforce 0 0 0
#弛豫10ps
dump 1 all atom 1000 relax.xyz
fix 1 all nvt temp 300 300 0.1
run 10000
unfix 1
undump 1
#冲击
dump 1 all atom 100 shock.xyz
fix 1 all nve
velocity left set 10 0 0 units box
run 2500
如需模拟其它材料的冲击，仅需更换结构模型和势函数即可
END
公众号案例代码以及学习交流，请加QQ群：754749935
关注微信公众号：lammps加油站
