帮忙看一下旋转力矩udf哪有问题？

chulx

一周前在fluent板块里发过求助，可能发错板块了，没有人回答，这里知道的朋友帮忙看看，Thanks

zone旋转中所需力矩UDF问题
想要用fluent求解风吹动叶片旋转，用的是sliding mesh，定义包含叶片的的内部圆柱区域进行旋转，zone的角速度通过UDF（define_zone_motion)，求解叶片上的力矩所得，叶片定义为旋转避免，速度相对于邻近单元为0。根据楼主的程序中力矩的计算公式，我写了自己udf，高手帮看看哪里有问题（绕Z轴旋转），计算中出现的问题：计算过程中区域不发生转动
#include "udf.h"
#define blade_id  27
  real Iz=6.25e-4;
  real W1= 0.0;
  real a0=0.0;
  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_CENTROID(x,f,t);      
       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;
}

[ 本帖最后由 chulx 于 2011-6-2 05:23 编辑 ]

gbd20022002

不好意思，我没见过这个宏DEFINE_ZONE_MOTION（），你可以说明来源和参数含义吗？

chulx

原帖由 gbd20022002 于 2011-5-25 15:18 发表

                               
登录/注册后可看大图

不好意思，我没见过这个宏DEFINE_ZONE_MOTION（），你可以说明来源和参数含义吗？

在fluent帮助中可以看到这个宏，下面是帮助中的解释；我以前看到过一篇文章，里面介绍的就是用ZONE_Motion的udf表示区域的运动，采用sliding mesh，这样就可以避免，内部叶轮旋转时的动网格问题了；

[ 本帖最后由 chulx 于 2011-6-1 07:43 编辑 ]

guidegege

a0是不是应该写在a1之后，dw之前？

hongfu2233

guidegege 发表于 2014-7-1 21:14
a0是不是应该写在a1之后，dw之前？

楼主解决没有？我也遇到类似问题了。。。