詹皇 发表于 2017-3-5 14:17:46

关于圆柱简谐运动

最近在做关于管道振动的仿真，管道里流动的是液固两相流，udf程序如下：
#include "udf.h"

/**计算*/
DEFINE_GRID_MOTION(pipe,domain,dt,time,dtime)
{
    Thread *tf = DT_THREAD(dt);
    face_t f;
    Node *v;
   real scale;
    real sign;
        real diam=0;
      real x=0;
    real z=0;
    int n;
    /* set deforming flag on adjacent cell zone */
       

    SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));
   sign = 0.004 * sin (10 * time);
   
    begin_f_loop(f,tf)
    {
       f_node_loop(f,tf,n)
       {
                  
      v = F_NODE(f,tf,n);

               if(NODE_Y(v)<=0.15)diam=0.03;
               elsediam=-(NODE_Y(v)-0.2)*0.6;
                  
                scale=(diam+sign)/diam;/**直线振动变化*/

      
      if (NODE_POS_NEED_UPDATE (v))
      {x=NODE_X(v);
               z=NODE_Z(v);
         
         NODE_POS_UPDATED(v);
               NODE_X(v)=diam*scale*(x/pow((x*x+z*z),0.5));
               NODE_Z(v)=diam*scale*(z/pow((x*x+z*z),0.5));
               /*Y轴固定*/
         }
   }
}
end_f_loop(f,tf);
}
请问大神程序有没有问题，dynamic mesh里该如何设置？为什么我的仿真结果里固体相在出口端呈阶梯式减少，结果与不震动情况下相差太大

望山跑 发表于 2017-3-6 13:12:03

你的圆柱是柔性圆柱么？

打俩鸡蛋 发表于 2017-8-27 18:14:57

这个模型及代码我见过，是收缩运动，但我没算
页: [1]
查看完整版本: 关于圆柱简谐运动