宏DEFINE_GRID_MOTION用来移动任意边界和流体区域内的网格节点。它提供了对节点和网格最大限度的操作，可以将刚体运动、变形和相对运动等结合起来。但是使用此UDF时，每一个时间步都必须执行。

DEFINE_GRID_MOTION限制：

1.  给定节点的运动是基于绝对坐标的，不可基于相对坐标的。

2.  给定节点的新位置是基于该节点上一个时间步的位置，而不是该节点的初始（t=0）位置

3.  宏DEFINE_GRID_MOTION中没有用户定义的节点位置内存记忆

4.  当使用宏DEFINE_GRID_MOTION时，节点移动的区域的连接关系不能改变。

 

摘自《ANSYS FLUENT技术基础与工程应用：流动传热与环境污染控制领域》，源代码有误，对原UDF进行了修正，添加注释

圆管在垂直流动方向上的振动方程如下：

式中，A表示振动圆柱的振幅，ω表示振动角频率，f表示振动频率。作为示例，在本例中，设A=4mm，f=25Hz

#include"udf.h"

#include"math.h"

#define omeg 50*3.14159265358979323846

#define A 0.004    //振幅

 

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

{

    Thread* tf=DT_THREAD(dt);     //DT_THREAD 就是获取网格运动所在的 thread 指针

    face_t f;

    Node* v;    //指向网格节点的指针

    int n;

    

    SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));

    

    Message("\nprevious time=%.5f\t current time=%.5f\n",PREVIOUS_TIME,CURRENT_TIME);    

    

    Message("time=%.5f\t dtime=%.5f\t n=%d",time,dtime,N_TIME);

    

    begin_f_loop(f,tf)

    {

        

        f_node_loop(f,tf,n)

        {

            

            v=F_NODE(f,tf,n);

            

            if(NODE_POS_NEED_UPDATE(v))        //在UDF手册当中没有查询到NODE_POS_NEED_UPDATED和NODE_POS_UPDATED宏

            {

                

                NODE_POS_UPDATED(v);

      NV_D(NODE_COORD(v),=,NODE_X(v),NODE_Y(v)+A*sin(omeg*(N_TIME+1)*CURRENT_TIMESTEP)-A*sin(omeg*N_TIME*CURRENT_TIMESTEP),NODE_Z(v));    

            }

        }

    }

    end_f_loop(f,tf);

}

 二维算例（点击图片查看动态图）：

三维算例（点击图片查看动态图）