FLUENT動(dòng)網(wǎng)格宏——DEFINE_CG_MOTION宏

2016-08-29 by:CAE仿真在線來(lái)源:互聯(lián)網(wǎng)

譯自《FLUENT UDF手冊(cè)》,2.6.1節(jié)

1、描述

通過(guò)使用DEFINE_CG_MOTION宏來(lái)指定特定動(dòng)態(tài)區(qū)域的運(yùn)動(dòng)。該宏能夠指定每一時(shí)間步的線速度與角速度。ANSYS FLUENT利用這些速度更新動(dòng)態(tài)區(qū)域的節(jié)點(diǎn)位置。該宏只能用于編譯型。

2、語(yǔ)法說(shuō)明

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

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

(2)dt:類型為Dynamic_Thread *dt,是一個(gè)包含了用戶指定的動(dòng)網(wǎng)格屬性結(jié)構(gòu)的指針。

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

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

(5)time:類型為real,表示當(dāng)前時(shí)間。

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

該宏的返回類型為void

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

3、實(shí)例

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

此處v為速度,F為外力,m為質(zhì)量。使用顯示歐拉格式表達(dá)t時(shí)刻速度為:

程序如下:

#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;