zone旋转中所需力矩UDF问题

chulx

想要用fluent求解风吹动叶片旋转，用的是sliding mesh，定义包含叶片的的内部圆柱区域进行旋转，zone的角速度通过UDF（define_zone_motion)，求解叶片上的力矩所得，叶片定义为旋转避免，速度相对于邻近单元为0。根据楼主的程序中力矩的计算公式，我写了自己udf，高手帮看看哪里有问题（绕Z轴旋转），计算中出现的问题：计算过程中区域不发生转动#include "udf.h"
#define blade_id  27
static real Iz=6.25e-4;
static real W1= 0.0;
static real a0=0.0;
static real a1=0.0;
  real dw;
  real moment[3];
  real pd=0;
  real x0=0,y0=0;
  face_t f;
  Thread *t;
  Domain *domain;
  
DEFINE_ZONE_MOTION(fanmotion,omega,axis,origin,velocity,time,dtime)
{
  domain = Get_Domain(17);

  t = Lookup_Thread(domain,blade_id);

  begin_f_loop(f,t)
{      
       real x[ND_ND],A[ND_ND];
       F_AREA(A,f,t);
           
      
    pd=F_P(f,t);
   moment[2] +=(pd*(x[1]*A[0]-x[0]*A[1]));      a1 = moment[2]/Iz;
   dw = dtime * ((a1+a0)/2.0);
   a0=a1;
   W1 += dw;   
  
   *omega=W1;
         }
end_f_loop(f,t)

  
  N3V_D (velocity,=,0.0,0.0,0.0);
  N3V_S(origin,=,0.0);        /* default values, line could be omitted */
  N3V_D(axis,=,0.0,0.0,1.0);     /* default values, line could be omitted */
  return;
}