|
|

楼主 |
发表于 2011-5-17 19:33:32
|
显示全部楼层
回复 2# czt1212 的帖子
谢谢您的提示!不过前面用同样的udf算过,虽然该动的movewall没动,但是没有出现上述错误提示。。不知道什么原因。。麻烦您帮我看看吧,udf如下:
#include "udf.h"
#include "unsteady.h"
#include "stdio.h"
#include "stdlib.h"
#include "memory.h"
DEFINE_CG_MOTION(movewall,dt,vel,omega,time,dtime)
{
Thread *t;
cell_t c;
real cx,cy;
real vx[3000],vy[3000];
int i;
int n=0;
if(!Data_Valid_P())
return;
t=DT_THREAD(dt);
cx=0;
cy=0;
begin_c_loop(c,t)
{
vx=C_U(c,t);
vy=C_U(c,t);
n++;
}
end_c_loop(c,t);
for(i=0;i<n;i++)
{
cx=cx+vx;
cy=cy+vy;
}
cx=cx/n;
cy=cy/n;
vel[0]=cx;
vel[1]=cy;
Message("\n n=%d, vx=%f, vy=%f\n",n,vel[0],vel[1]);
}
多谢啦! |
|