在动网格中，对于那些既包含了运动也包含了变形的区域，可以通过UDF来指定区域中每一个节点的位置。这给了用户最大的自由度来指定网格的运动。在其他的动网格技术中（如重叠网格）则很难做到这一点。定义网格节点的位置是通过UDF宏DEFINE_GRID_MOTION来实现的。

步骤

可以通过以下步骤来指定既包含刚体运动又包含变形的区域：

http://i4.buimg.com/567571/426269c3a889376e.png

• 从Zone Names列表项中选择运动区域

• 选择Type为User-Defined

• 在Motion Attributes中，选择UDF函数。此时的UDF用的是DEFINE_GRID_MOTION宏。

宏的描述

DEFINE_GRID_MOTION ( name , d , dt , time , dtime)

参数：

name：自定义宏的名称，symbol name。

d：指向区域的指针，Domain *d。

dt：指向存储动网格属性的结构体指针，Dynamic_Thread *dt。

time：当前时间，real time。

dtime：时间步长。

宏中包含了5个参数：name，d，dt，time，dtime。其中name为用户自定义，而d，dt，time以及dtim均由Fluent传入。

示例文件

以下示例宏文件表达了一个简支梁的弯曲。

http://i4.buimg.com/567571/fb8d3ed2d9c4f804.png

#include "udf.h"



 DEFINE_GRID_MOTION(beam,domain,dt,time,dtime)

{

     Thread *tf = DT_THREAD(dt);

     face_t f;

     Node *v;

     real NV_VEC(omega), NV_VEC(axis), NV_VEC(dx);

     real NV_VEC(origin), NV_VEC(rvec);

     real sign;

     int n;



     SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));

     sign = -5.0 * sin (26.178 * time);

     Message ("time = %f, omega = %f\n", time, sign);

     NV_S(omega, =, 0.0);

     NV_D(axis, =, 0.0, 1.0, 0.0);

     NV_D(origin, =, 0.0, 0.0, 0.152);

     begin_f_loop(f,tf)

     {

        f_node_loop(f,tf,n)

        {

         v = F_NODE(f,tf,n);



         if (NODE_X(v) > 0.020 && NODE_POS_NEED_UPDATE (v))

         {



          NODE_POS_UPDATED(v);

          omega[1] = sign * pow (NODE_X(v)/0.230, 0.5);

          NV_VV(rvec, =, NODE_COORD(v), -, origin);

          NV_CROSS(dx, omega, rvec);

          NV_S(dx, *=, dtime);

          NV_V(NODE_COORD(v), +=, dx);

         }

      }

  }

  end_f_loop(f,tf);

}

下次会演示这个宏的效果。

更多关于CFD内容，可微信扫描下方二维码关注微信公众号。

https://ooo.0o0.ooo/2017/06/13/593fe2f2d0c06.jpg