netgen對于網格的局部處理和細化

局部處理

這裡的局部處理指點，線，面的大小

nglib.h提供的接口

Ng_OCC_SetLocalMeshSize

meshclass.hpp提供的接口

RestrictLocalH

//使用occ需要在pro中增加
DEFINES +=OCCGEOMETRY \
MSVC_EXPRESS\

#include <occgeom.hpp>
//nType為id的類型
bool Widget::setMeshSizeH(int nType, int nId, double dh)
{
    Mesh* me = (Mesh*)m_pMesh;
    Ng_OCC_SetLocalMeshSize(pOccgeo,m_pMesh,&m_mparam);
    switch(nType)
    {
    case 0://面
    {
        //        Element2d& ele = me->SurfaceElement(nId);
        //        for(auto n  =0;n<ele.GetNP();n++)
        //            me->RestrictLocalH(RESTRICTH_POINT,n,dh);
        //擷取面，擷取面的中心點對與面的中心點進行大小設定
        const Element2d & sel = me->SurfaceElement(nId);
        Point3d p = Center(me->Point(sel.PNum(1)),
                           me->Point(sel.PNum(2)),
                           me->Point(sel.PNum(3)));
        me->RestrictLocalH (p, dh);
        break;
    }
    case 1://線
    {
    //擷取線設定線的Id
        const Segment & seg = me->LineSegment(nId);
        me->RestrictLocalHLine(me->Point(seg[0]),me->Point(seg[1]), dh);
        break;
    }
    case 3://點
    {
        me->RestrictLocalH((resthtype)nType,nId,dh);
        break;
    }
    default:
        break;
    }
    Ng_OCC_GenerateEdgeMesh(pOccgeo,m_pMesh,&m_mparam);
    Ng_OCC_GenerateSurfaceMesh(pOccgeo,m_pMesh,&m_mparam);
    Ng_GenerateVolumeMesh(m_pMesh,&m_mparam);
    int nNum = Ng_GetNE(m_pMesh);
    qDebug()<<"box size"<<nNum;
    return false;
}

           

局部點的控制效果

局部線的控制效果

局部面的控制效果

細化操作

Mesh* me = (Mesh*)m_pMesh;
    if(!me || !me->GetGeometry())
        return;
    me->GetGeometry()->GetRefinement().Refine(*me);