高熵合金FeNiCrCoAl磨削，金刚石球旋转+直线同时运动

大家好，我是小马老师。

本文介绍一个lammps模拟高熵合金摩擦的一个案例代码。

传统的摩擦过程中，磨头只能存在一种运动，在本案例中，磨头同时存在两种运动：旋转+直线运动。

磨头沿y直线摩擦的时候，自身也发生旋转，效果图如下图所示：

模拟文件分两部分，in.relax是高熵合金建模和弛豫代码，in.friction是摩擦代码。

旋转和直线运动用到的命令为fix move transrot，这是lammps 2020新版本才有的一个命令，所以，运行下面的代码要更新到最新的lammps版本。

fix move transrot的命令格式为

fix id group_id move transrottransrot Vx Vy Vz Px Py Pz Rx Ry Rz period
1

其中：

• Vx,Vy,Vz 为摩擦球运动速度
• Px,Py,Pz 为旋转轴的初始位置
• Rx,Ry,Rz 为旋转轴的方向
• period 为旋转周期

在本例中，先计算摩擦球的重心坐标，

variable    ball_x equal xcm(<