|
|
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ 小木虫: 金币+0.5, 给个红包,谢谢回帖 云海孤狼: 金币+10 2013-08-20 10:30:31 xiegangmai: 金币+3, 谢谢参与! 2013-08-25 00:09:35
这个不算难,我帮你写了一个,有点仓促,你试一下,有问题再交流
#include “udf.h”
#include “math.h”
DEFINE_PROFILE(velocity,t,i)
{
real x[ND_ND];
real n=3; /*根据你自己的值更改*/
real h0=10;
real V0=3;
face_t f;
real h;
begin_f_loop(f,t)
{
F_CENTROID(x,f,t);
h = x[2]; /*此处假设是z方向*/
F_PROFILE(f,t,i) = V0*pow((h/h0),n);
}
end_f_loop(f,t)
} |
|