| 查看: 518 | 回复: 1 | |||
| 【悬赏金币】回答本帖问题，作者dwademvp将赠送您 20 个金币 | |||
[求助]
求大佬帮我看看这个udf
|
|||
|
计算壁面条件的，这个udf是仿照帮助文档写的，能加载，但碰撞的粒子都abort了，无反弹的粒子，想知道哪里出问题了。 #include "udf.h" DEFINE_DPM_BC(my_dpm_bc,p,t,f,f_nomal,dim) { real vn=0.; real alpha; real vc; real vn1; real nor_coeff=0.8; real tan_coeff=0.3; real nomal[3]; real diam; real NV_VEC(x); diam=P_DIAM(p); int i,idim=dim; for(i=0;i<idim;i++) nomal=f_nomal; for(i=0;i<idim;i++) vn+=p->state.V*nomal; vn1=fabs(vn); vc=pow((2*0.00000018/diam/45.3/45.3*pow((45.3+pow(vn1,0.718)),2)),1.43); if(vn1<vc) { return PATH_ABORT; } else if(p->type==DPM_TYPE_INERT) { alpha=M_PI/2.-acos(MAX(-1.,MIN(1.,NV_DOT(nomal,p->state.V)/ MAX(NV_MAG(p->state.V),DPM_SMALL)))); if((NNULLP(t))&&(THREAD_TYPE(t)==THREAD_F_WALL)) F_CENTROID(x,f,t); for(i=0;i<idim;i++) p->state.V-=vn*nomal; for(i=0;i<idim;i++) p->state.V*=tan_coeff; for(i=0;i<idim;i++) p->state.V-=nor_coeff*vn*nomal; for(i=0;i<idim;i++) p->state0.V= p->state.V; return PATH_ACTIVE; } } |
回复此楼» 猜你喜欢
一志愿哈尔滨工业大学0856材料与化工，前三科206，总分283，求调剂
已经有7人回复
-
A区一本交叉课题组，低分调剂，招收机械电子信息通信等交叉方向
已经有48人回复
-
288求调剂085600材料与化工
已经有7人回复
-
278求调剂
已经有5人回复
-
物理学求调剂
已经有3人回复
-
0703求调剂
已经有3人回复
-
322分 085600求调剂，有互联网+国金及主持省级大创经历
已经有3人回复
-
材料专硕323求调剂
已经有4人回复
-
招调剂
已经有15人回复
-
调剂的同学，走过路过，不要错过
已经有8人回复
2楼2022-06-09 15:04:43
