(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211366341.1 (22)申请日 2022.11.03 (71)申请人 南京航空航天大 学 地址 210001 江苏省南京市秦淮区御道街 29号 (72)发明人 王逸斌　邓思强　覃宁　赵宁　 朱春玲　 (74)专利代理 机构 江苏圣典律师事务所 32 237 专利代理师 贺翔 (51)Int.Cl. G06T 17/20(2006.01) G06F 30/28(2020.01) G06F 113/08(2020.01) G06F 119/14(2020.01) (54)发明名称 一种基于Delaunay图的多块结构网格奇点 识别方法 (57)摘要 本发明公开了一种基于Delaunay图和 交叉 场的奇点识别方法， 该方法基于Delau nay图和交 叉场， 为生成高质量的四边形结构网格提供依 据， 只用到边界上的节点和添加的辅助点进行 Delaunay三角化， 在得到的Delaunay图中， 仅利 用这些节点上的交叉即可计算得到奇点的位置 和方向， 根据几何外形， 可利用奇点进行多块结 构网格划分； 该方法算法简单， 计算量小， 能够辅 助使用者划分高质量的多块结构网格 。 权利要求书2页 说明书6页 附图5页 CN 115471635 A 2022.12.13 CN 115471635 A 1.一种基于Delaunay图的多块结构网格奇点识别方法， 其特 征在于， 包括以下步骤： S1： 获取几何构型的边界节点； S2： 计算所述 边界节点上的交叉； S3： 在计算 域内添加辅助点， 并指定所述辅助点上的交叉； S4： 对所述 边界点和辅助点进行Delaunay三角化， 得到Delaunay图； S5： 循环Del aunay图中的所有单元， 判断其中是否存在奇点， 若存在奇点确定所述奇点 的类型； S6： 计算所有奇点的位置和方向， 并以此辅助结构网格分块。 2.根据权利要求1所述的多块结构网格奇点识别方法， 其特征在于， 所述步骤S1具体包 括： 将几何构型的曲线和直线离 散成节点形成边界节点。 3.根据权利要求1所述的多块结构网格奇点识别方法， 其特征在于， 所述步骤S2具体包 括： 交叉是相互垂直的四个向量的几何， 用一个角度 θ表示： 式中，cθ表示交叉， 表示交叉中的某个向量， (*)T表示矩阵的转置； 对于边界上的节点， 其与相邻两点组成的内角大小为 θc， 定义一个整数 ， 式中，round(x)函数表示对实数 x进行取整； 当 nc为偶数时， 交叉方向与内角平分线的方向 一致； 当nc为奇数时， 交叉 方向为角平分线方向逆时针旋转 45°。 4.根据权利要求1所述的多块结构网格奇点识别方法， 其特征在于， 在所述步骤S3中： 所述辅助点上 交叉方向应与周围交叉方向保持连续， 且相 邻两个交叉的连续要求其方向的 角度差的绝对值 不超过45度。 5.根据权利要求1所述的多块结构网格奇点识别方法， 其特征在于， 所述步骤S4进一步 包括： S41： 构造包 含所有散点的超级三角形， 储 存在三角形链 表中； S42： 依次插入散点， 在三角形链表中找出外接圆包含插入点的影响三角形， 删除影响 三角形的公共边， 连接插 入点与影响三角形的所有顶点； S43： 根据Delaunay准则优化得到的新 三角形， 并放入三角形链 表； S44： 重复步骤S42直到所有散点插 入完成， 并删除与超级三角形相关的三角形。 6.根据权利要求1所述的多块结构网格奇点识别方法， 其特征在于， 所述步骤S5进一步 包括： 循环所有Delaunay三角形， 其 三个节点上的交叉组成了一个局部交叉场； 利用交叉场的连续， 判断一个局部交叉场内是否存在奇点。 7.根据权利要求6所述的多块结构网格奇点识别方法， 其特征在于， 所述步骤S5进一步 包括： 在局部交叉场中， 三角形的三个顶点1、 2、 3以逆时针顺序排列， 对应的交叉分别为 θ1、θ2 和θ3， 则奇点的类型由下式计算，权　利　要　求　书 1/2 页 2 CN 115471635 A 2式中，k是奇点的类 型，∆θij=θj‑θi表示两个交叉之间的最小角度差， 其范围是[ ‑π/4， π/ 4]。 8.根据权利要求1所述的多块结构网格奇点识别方法， 其特征在于， 所述步骤S6进一步 包括： 对于包含奇点的三角形， 其顶点的坐标分别为( x1,y1)、 (x2,y2)、 (x3,y3)， 奇点的坐标 S (xs,ys)计算为 对于类型 k的奇点， 共有 k+4个方向； 奇点方 向的计算首先要确定一个初始方向 β0， 所述 奇点的其它方向则计算 为， 对于初始方向 β0， 首先确定所处三角形中的交叉之差最小的边， 所述边的中点求为 β (xβ,yβ)， 则β0即为矢量 的方向， 即 式中， (xs,ys)和(xβ,yβ)分别是点 S和β的坐标。权　利　要　求　书 2/2 页 3 CN 115471635 A 3