lammps也玩凌波微步，原子组瞬间转移命令介绍

大家好，我是小马老师。
对于多部分组合的模型，如摩擦模型，在弛豫阶段，如果希望刀具和工件分开弛豫，或者在弛豫过程中，刀具和工件不会产生相互作用力，通用做法是把刀具设置到离工件相对较远的位置。
当刀具和工件处于力的截断半径之外时，不会产生相互作用力。
弛豫完成后，进入摩擦阶段，又希望刀具能够快速的移动到工件的表面。
这时候，最简单的方式是使用displace_atoms命令快速移动刀具。
displace_atoms命令语法规则为
displace_atoms group-ID move delx dely delz
displace_atoms设置move关键词可实现group组原子在x、y、z方向的坐标值按照增量“delx dely delz”变化，从而实现快速移动。
“delx dely delz”值为正表示沿坐标轴正方向移动，值为负表示沿坐标轴负方向移动。
下面以简化版的摩擦模型为例，演示displace_atoms具体使用方法。
模型底部为工件Cu，上部为金刚石刀具，初始模型为：

快速移动后，模型为：

全部in文件如下：

units metal
atom_style atomic
boundary p p p
#create box
region  box block 0 200 0 40 0 100 units box
create_box 2 box
#create Cu
region piece block 0 200 0 40 0 40 units box
lattice fcc 3.61 
create_atoms 1 region piece
#create ball
region ball sphere 100 20 85 15 units box
lattice       diamond  3.52 
create_atoms  2 region ball
group ball type 2
#mass define
mass    1 64
mass    2 12
#save data
write_data all1.data
#move ball
displace_atoms ball move 0 0 -30 units box
#save data
write_data all2.data

如果在使用lammps过程中遇到什么问题，可到公众号或者QQ群留言。
公众号案例代码以及学习交流，请加QQ群：754749935

扫描关注微信公众号：lammps加油站。