|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
在学习fluent 动网格中遇到问题了
有谁知道 udf错哪里了??????
#include "udf.h"
#include "unsteady.h"
#include "memory.h"
#include "stdlib.h"
#include "stdio.h"
#include "metric.h"
#include "mem.h"
#include "global.h"
DEFINE_CG_MOTION(movewall,dt,vel,omega,time,dtime)
{ Thread *t; /*************10**********/
cell_t c;
real cx,cy,x,y;
//CArray<real,real&> vx,vy;
int i;
int m,n=0; /*************15**********/
t=DT_THREAD(dt);
//real xx[ND_ND];
//real x,y;
x=y=0; /*************20**********/
Node *node;
//int m;
//if(!Data_Valid_P())
//return; /*************25**********/
cx=0.;
cy=0.1;
//real zhongxin=0;
begin_c_loop(c,t)
{ /*************30**********/
//y=x+y;
// x.Add(C_U(c,t));
// y.Add(C_U(c,t));
cx=cx+C_U(c,t);
cy=cy+C_V(c,t); /*************35**********/
n++;
c_node_loop(c,t,m)
{
y=y+NODE_Y(node);
// x=x+NODE_X(node); /*************40**********/
//x=C_NODE(c,t,m);
//C_CENTROID(xx,c,t);
Message("\n m= %d , x= %f , y= %f , \n",m,x,y);
//Message("\n m= %d , x= %f , y= %f , zhongxin= %f \n",m,x,y,zhongxin);
} /*************45**********/
}
end_c_loop(c,t);
cx=cx/n;
cy=cy/n; /*************50**********/
cx=0;
cy=0;
vel[0]=cx;
vel[1]=cy;
Message("\n n=%d, vx=%f, vy=%f, vy=%f\n",n,vel[0],vel[1]);/*************55**********/
}
哪里有问题,为什么不能编译,我想计算单元网格的中心位置
谢谢 |
|