大家好,我是小马老师。

本文介绍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为旋转周期,单位为当前模拟体系的时间单位,周期越小,旋转速度越快。

lammps教程:工件旋转设置方法_微信

以金刚石球为例,旋转轴为z轴正方向,逆时针旋转,代码为:



fix 1 ball move rotate 0 0 0 0 0 1 10 units box

lammps教程:工件旋转设置方法_一对一_02

顺时针旋转:



fix 1 ball move rotate 0 0 0 0 0 -1 10 units box

lammps教程:工件旋转设置方法_一对一_03

旋转速度变慢(周期变长):



fix 1 ball move rotate 0 0 0 0 0 -1 50 units box

lammps教程:工件旋转设置方法_微信_04

需要注意的是同一组原子的fix move命令不能与系综同时设置,否则会造成模拟出错。

 



公众号案例代码以及学习交流,请加QQ群:754749935


点击关注微信公众号:lammps加油站