FLUENT動網格宏——DEFINE_CG_MOTION宏

2016-08-29 by:CAE仿真在線來源:互聯網

譯自《FLUENT UDF手冊》,2.6.1節

1、描述

通過使用DEFINE_CG_MOTION宏來指定特定動態區域的運動。該宏能夠指定每一時間步的線速度與角速度。ANSYS FLUENT利用這些速度更新動態區域的節點位置。該宏只能用于編譯型。

2、語法說明

DEFINE_CG_MOTION(name,dt,vel,omega,time,dtime)

(1)name:宏名,由用戶指定,需要符合UDF中宏命名規范。

(2)dt:類型為Dynamic_Thread *dt,是一個包含了用戶指定的動網格屬性結構的指針。

(3)vel:數據類型為real vel[],速度向量,vel[0]表示x方向線速度,vel[1]表示y方向,vel[2]為Z方向。

(4)omega:real omega[],角速度向量。索引為0~2,分別表示x,y,z方向角速度。

(5)time:類型為real,表示當前時間。

(6)dtime:類型為real,表示時間步

該宏的返回類型為void

此宏包含有六個參數,其中name是用戶提供,其他的參數如dt,vel,omega,time,dtime均有ansys fluent求解器傳遞。線速度與角速度能返回至ansys fluent。

3、實例

以下的例子中,線速度是通過物體上的x方向的力平衡達到的。表達形式為:

此處v為速度,F為外力,m為質量。使用顯示歐拉格式表達t時刻速度為:

程序如下:

#include “udf.h”

static real v_prev=0.0;

DEFINE_CG_MOTION(piston,dt,vel,omega,time,dtime)

{

Thread *t;

face_t f;

real NV_VEC(A);

real force,dv;

NV_S(vel,=,0);

NV_S(omega,=,0);

if(!Data_Valid_P())

{

     return ;
    }

t=DT_THREAD(dt);

force = 0.0;

begin_f_loop(f,t)

{

    F_AREA(A,f,t);

    force += F_P(f,t)* NV_MAG(A);

}

end_f_loop(f,t)

dv = dtime * force /50.0;

v_prev += dv;

Message(“time =%f, x_vel=%f, force=%f\n”,time,v_prev,force);

vel[0] = v_prev;