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);