|
![](static/image/common/ico_lz.png)
楼主 |
发表于 2017-6-19 08:05:23
|
显示全部楼层
# include "udf.h"
const int zone_ID=3;
DEFINE_DPM_BODY_FORCE(bforce,p,i)
{
int n;
real bforce,acc;
real A[ND_ND];
cell_t c;
Thread *tc,*tf;
Domain*d;
d=Get_Domain(1);
c=RP_CELL(&p->cCell);
tc=RP_THREAD(&p->cCell);
c_face_loop(c,tc,n)
{
tf=C_FACE_THREAD(c,tc,n);
if(zone_ID==THREAD_ID(tf))
{
if(i==1)
{
bforce=1.001e-4;
}
else
{
bforce=0;
}
}
else
{
bforce=0;
}
acc=bforce/P_MASS(p);
return acc;
}
}
这个是写的udf,编译出错。在编号为3的区域加载体积力。大神帮忙看看?
|
|