动网格UDF案例.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,动网格,UDF,案例,问题运动,绿色的阀体作定速转动,黄色膜片作随时间做正,/,余旋的变形,动网格区域,阀体,DEFINE_CG_MOTION,膜片,DEFINE_GRID_MOTION,上部壁面,DEFINE_GEOM,DEFINE_CG_MOTION,阀体转动,#define omega 1.0,DEFINE_CG_MOTION(butterfly_flex_UDF,dt,cg_vel,cg_omega,time,dtime),cg_vel0=0.0;,cg_vel1=0.0;,cg_vel2=0.0;,cg_omega0=0.0;,cg_omega1=0.0;,cg_omega2=omega;,DEFINE_GRID_MOTION,DEFINE_GRID_MOTION(moving_arc,domain,dt,time,dtime),Thread*tf=DT_THREAD(dt);,face_t f;,Node*node_p;,real alpha,theta,x,ymag,yfull,y;,int n;,SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf);,alpha=omega*CURRENT_TIME;,theta=2.0*alpha+3.0*M_PI/2.0;,begin_f_loop(f,tf),f_node_loop(f,tf,n),node_p=F_NODE(f,tf,n);,if(NODE_POS_NEED_UPDATE(node_p),NODE_POS_UPDATED(node_p);,x =NODE_X(node_p);,ymag =sqrt(R*R-x*x)+0.03;,yfull=ymag-0.1;,y =-0.1+yfull*sin(theta);,NODE_Y(node_p)=y;,end_f_loop(f,tf);,DEFINE_GEOM,#define R 0.109,DEFINE_GEOM(plane,domain,dt,position),position1=R;,