大家好,我是小马老师。
本文介绍一个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
其中:
- Vx,Vy,Vz 为摩擦球运动速度
- Px,Py,