大家帮帮忙吧，动网格。

41906957

我想得到类似于红色图的翼型上表面正弦运动的动网格，但是在我计算中总是算着算着就出现负体积了，下面是我的cas文件，和udf文件，哪位高手给我看一下吧。谢谢.
#include  "udf.h"
#define pi 3.1415926
DEFINE_GRID_MOTION(beam,domain,dt,time,dtime)
{
    Thread  *tf  =  DT_THREAD(dt);
    face_t  f;
    Node  *v;
    real  NV_VEC(dx);
    int  n;

    SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));
    begin_f_loop(f,tf)
        {
    f_node_loop(f,tf,n)
        {
    v  =  F_NODE(f,tf,n);
/*  只让翼型上表面的一部分边界运动  */
    if  (NODE_X(v)  >  0.20  && NODE_X(v)  <=  0.50 && NODE_Y(v)  >  0.0 && NODE_POS_NEED_UPDATE  (v))
        {

        NODE_POS_UPDATED(v);
        dx[0]=0;
                dx[1]=0.01*sin(pi*NODE_X(v)/0.1);
                dx[2]=0;
        NV_V(NODE_COORD(v),  +=,  dx);
        }
        }
        }
    end_f_loop(f,tf);
}

[ 本帖最后由 41906957 于 2012-1-9 19:35 编辑 ]

zhanghejn

负体积是网格质量造成的，逐步光滑网格很重要。

shmily0527

网格问题