关于difine_cg_motion求高手指点

geneexz2005

这是一个在原点让转子转动udf编译如果沿着x方向移动10mm,在中心不是(0.0)的情况下,在udf里面需要怎么修改?

#include "udf.h"
#define r 0.006625 /*Eccentric distance m*/
#define k 57.2727 /*rad to degree*/
#define w 125.6571429 /*angular velocity in rad/sec (3600rpm)*/

/* for grid motion */
/* rotor is profile name in lab */
DEFINE_CG_MOTION(rotor, dt, vel, omega,time, dtime)
{
real beta;

NV_S(vel, =, 0.0);
NV_S(omega, =, 0.0);

omega[2] = -125.6571429; /*angular velocity for rotor in rad/sec*/
beta = 376.9714286; /*angular velocity of eccentricshaft in rad/sec*/

vel[1] = beta*r*sin(beta*time);
vel[0] = -beta*r*cos(beta*time);
vel[2]= 0;

}