fluent14计算3D嵌入滑移网格遇到问题，求助

参考沙哥的一个动网格例子， http://blog.sina.com.cn/s/blog_6d84e3950101a5rr.html
采用define zone motion计算三层流流体相对旋转的流动，瞬态计算
外层不动；中间层旋转速度500rad/s,旋转中心为原点指向z轴；内层在x轴方向偏心4mm，跟着中间层移动，并且自转速度1000rad/s；
内层采用的udf如下，
#include "udf.h"
DEFINE_ZONE_MOTION(rotor,omega,axis,origin,velocity,time,dtime)
{
real theta0,theta1,omegaC,omegaR,radR,ox1,oy1;
theta0=0;
omegaR=1000.0;/*内层转速**/
omegaC=500.0; /*中间层转速
theta1=omegaC*time;内层旋转中心角度位置
radR=4/1000;
ox1=radR*cos(theta0+theta1);内层旋转中心x坐标
oy1=radR*sin(theta0+theta1);内层旋转中心y坐标
origin[0]=ox1;
origin[1]=oy1;
origin[2]=0.0;
axis[0]=ox1;
axis[1]=oy1;
axis[2]=10;
*omega=omegaR;
}
初始网格见图1，过几个时间步之后的网格见图2，已经发生了错位，
调试了半天，觉得udf中定义内层旋转中心的语句无效，即origin和axis语句无效
实际内层绕坐标原点（0,0,1） 以omegaR角速度 旋转，但这不是我想要的，请高手看看出什么问题，谢谢