大家好,我是小马老师。
本文介绍lammps模拟中工件旋转的设置方法。
在扭转、抛光等lammps模拟中,需要设置工件的旋转动作。
实现旋转的命令为fix move rotate,命令格式为:
fix ID group-ID move rotate Px Py Pz Rx Ry Rz period
其中:
group-ID:工件原子组
Px,Py,Pz:旋转轴上某点的坐标
Rx,Ry,Rz:旋转轴矢量
period:旋转周期
旋转设置的原理是确定旋转轴、旋转方向以及旋转周期。
Px,Py,Pz为旋转轴上的一个点的坐标,比如球心或者圆柱中心轴的任意一个点坐标。
Rx,Ry,Rz确定旋转轴矢量,如指向z轴正方向可设为“0 0 1”,旋转方向符合右手法则。
period为旋转周期,单位为当前模拟体系的时间单位,周期越小,旋转速度越快。
以金刚石球为例,旋转轴为z轴正方向,逆时针旋转,代码为:
fix 1 ball move rotate 0 0 0 0 0 1 10 units box
顺时针旋转:
fix 1 ball move rotate 0 0 0 0 0 -1 10 units box
旋转速度变慢(周期变长):
fix 1 ball move rotate 0 0 0 0 0 -1 50 units box
需要注意的是同一组原子的fix move命令不能与系综同时设置,否则会造成模拟出错。
公众号案例代码以及学习交流,请加QQ群:754749935
点击关注微信公众号:lammps加油站