函数Compute_Force_And_Moment

2017-01-16 by:CAE仿真在线来源:互联网

FLUENT中可以调用函数 Compute_Force_And_Moment()进行力和力矩的计算。该函数的声明位于f_wall.h头文件内。没有任何对该函数的官方什么。

该函数的描述:

Compute_Force_And_Moment (d, t, x_cg, f_body, m_body, TRUE)

fluent软件中一个内部函数,主要是用以计算三维物体的受力。其中t指向所要求解的物体边界,x_cg是物体边界的中心位置,f_body是物体所受水作用力,m_body是物体所受外力对于中心位置的力矩。

函数说明:

1. 这个宏不需要执行循环操作;

2. domain是计算区域,tf1是面索,也就是你要计算受力的面,x_cg是面的形心位置,f_glob与m_glob分别指在惯性系下的力与力矩分量;

3. 该宏既可以计算3D受力,也可用于2D受力计算;

4. TRUE与FALSE是逻辑类型的量,TRUE就是非零的整数即可,当然也可以直接写1,FALSE当然就是0了;如果是1则调用该宏,0的话则不调用;

5. 可以在并行版中应用,不过得注意将UDF改写为并行版本的。

以上资料来源:http://www.efluid.com.cn/bbs/showtopic-4306.aspx

范例:

#include "udf.h"
#include "mem.h"
#include "dynamesh_tools.h"

static real v_x = 0.0,v_y=0.0,omega_z=0.0;
DEFINE_CG_MOTION(piston,dt,vel,omega,time,dtime)
{
face_t f;
cell_t c;
real f_glob[ND_ND],m_glob[ND_ND],x_cg[ND_ND],dv_x,dv_y,domega_z;
Domain *domain= Get_Domain (1);
Thread *tf1 = Lookup_Thread (domain, 3);
int i;

NV_S(vel, =, 0.0);
NV_S(omega, =, 0.0);
if (!Data_Valid_P())
return;

for(i=0;i<=ND_ND;i++)
{
f_glob[i]=0;
m_glob[i]=0;
}

for(i=0;i<ND_ND;i++)
x_cg[i]=DT_CG(dt)[i];

if(time<=0.001)
{
vel[1]=34.5;
Message ("ntime=%f ,x=%.10lf ,y=%.10lf , force_x=%.1f ,force_y=%.1f ,moment_x=%f,moment_y=%f,moment_z=%f ,Vx=%f ,Vy=%f ,omega_z=%fn", time,x_cg[0],x_cg[1],f_glob[0],f_glob[1],m_glob[0],m_glob[1],m_glob[2],v_x,v_y,omega_z);
return;
}

Compute_Force_And_Moment (domain, tf1, x_cg, f_glob, m_glob, TRUE);

dv_x = dtime * f_glob[0] / 39400;
v_x += dv_x;

dv_y = dtime * f_glob[1] / 39400;
v_y += dv_y;

domega_z=dtime*m_glob[2]/390000;
omega_z+=domega_z;

Message ("ntime=%f ,x=%.10lf ,y=%.10lf , force_x=%.1f ,force_y=%.1f ,moment_x=%f,moment_y=%f,moment_z=%f ,Vx=%f ,Vy=%f ,omega_z=%fn", time,x_cg[0],x_cg[1],f_glob[0],f_glob[1],m_glob[0],m_glob[1],m_glob[2],v_x,v_y,omega_z);