UDF程序文档格式.docx
- 文档编号:8383046
- 上传时间:2023-05-11
- 格式:DOCX
- 页数:3
- 大小:17.17KB
UDF程序文档格式.docx
《UDF程序文档格式.docx》由会员分享,可在线阅读,更多相关《UDF程序文档格式.docx(3页珍藏版)》请在冰点文库上搜索。
realvt[3];
realvn1[3];
realangle=0;
inti,idim=dim;
realNV_VEC(x);
#ifRP_2D
/*dimisalways2in2Dcompilation.Needspecialtreatmentfor2d
axisymmetricandswirlflows*/
if(rp_axi_swirl)
{
realR=sqrt(P_POS(p)[1]*P_POS(p)[1]+P_POS(p)[2]*P_POS(p)[2]);
if(R>
1.e-20)
{
idim=3;
normal[0]=f_normal[0];
normal[1]=(f_normal[1]*P_POS(p)[1])/R;
normal[2]=(f_normal[1]*P_POS(p)[2])/R;
}
else
for(i=0;
i<
idim;
i++)
normal[i]=f_normal[i];
}
else
#endif
for(i=0;
normal[i]=f_normal[i];
/*Computenormalvelocity.*/
for(i=0;
vn+=P_VEL(p)[i]*normal[i];
/*Computeparallelvelocity.*/
for(i=0;
i<
i++)
P_VEL(p)[i]
-=
vn*normal[i];
{
vt[i]=P_VEL(p)[i];
vt[i]=1/7*(5*vt[i]+P_DIAM(p)*angle);
angle=2*NV_MAG(vt)/P_DIAM(p);
Vtx=normal[1]*P_VEL(P)[2]-normal[2]*P_VEL(p)[1];
Vty=normal[0]*P_VEL(P)[2]-normal[2]*P_VEL(p)[0];
Vtz=normal[0]*P_VEL(P)[1]-normal[1]*P_VEL(p)[0];
Vt1=sqrt(Vtx*Vtx+Vty*Vty+Vtz*Vtz);
/*Computeanglevelocity.*/
If()
Vt2=1/7*(5*vt1+P_DIAM(p)*angle1);
Vn2=-en*vn;
angle2=2*ut2/P_DIAM(p);
/*Subtractoffnormalvelocity.*/
for(i=0;
P_VEL(p)[i]-=vn*normal[i];
/*Applytangentialcoefficientofrestitution.*/
P_VEL(p)[i]*=tan_coeff;
/*Addreflectednormalvelocity.*/
P_VEL(p)[i]-=nor_coeff*vn*normal[i];
/*StorenewvelocityinP_VEL0ofparticle*/
for(i=0;
P_VEL0(p)[i]=P_VEL(p)[i];
if(fabs(u))
if(p->
type==DPM_TYPE_INERT)
alpha=M_PI/2.-acos(MAX(-1.,MIN(1.,NV_DOT(normal,P_VEL(p))/
MAX(NV_MAG(P_VEL(p)),DPM_SMALL))));
if((NNULLP(t))&
&
(THREAD_TYPE(t)==THREAD_F_WALL))
F_CENTROID(x,f,t);
/*calculatethenormalcomponent,rescaleitsmagnitudeby
thecoefficientofrestitutionandsubtractthechange*/
/*Computenormalvelocity.*/
vn+=P_VEL(p)[i]*normal[i];
returnPATH_ACTIVE;
}
returnPATH_ABORT;
}
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- UDF 程序