一种改进的Delaunay三角形化剖分方法.pdfVIP

2001年第25卷 石油大学学报(自然科学版) vol25No2 第2期 oftl】e of Apr．2001 JoumaIUnivershyPend锄，Chi呻 文章编号：1000—5870(2001)02．0100．06 一种改进的Delaunay三角形化剖分方法 徐明海1， 张俨彬2， 陶文铨3， 摘要：提出了一种基于B。w弦wa幻m算法的平面区域Ⅸ岫三角化剖分的改进方法。它结合了前沿推进法的 内部结点生成技术和Dda衄y联点网格生成技术，使得每插入一点所破坏的单元尽可能地少。采用适当的数据结构， 使IH删y搜索过程限于局部，算法大为简化，易于编程，浮点计算量少，同时也避免了使用函数递归调用。采用在基 网格上定义网格步长的办法控制阿格的疏密，使网格疏密易于控制。几个算例表明，该算法是行之有效的。 关键词：非结构网格；三角形单元；Dd虬m町剖分；数值模拟 121 中围分类号：TK 文献标识码：A 法。 引 言 1 Dda吼ay网格生成算法原理 非结构化网格是相对于结构化网格而言的，它具 有分布灵活、适应复杂区域剖分及局部任意加密的优 Delaunay原理要求每个三角形的外接圆内不包 点。非结构网格需要较多的前期处理工作，网格形成 含其他节点，这样的三角形系统具有最小角最大的 一般不能用手工操作，须采用计算机自动剖分。三角 特点。根据这一原理，可以构造出各种生成二维区 形网格自动剖分以DdⅢw三角化和前沿推进法两 域网格的算法。 类方法为基础。为了提高自动剖分程度、减小剖分的 Sharn。s和H)eyL川提出了分治法的思想，L＆ 手工工作量和运行机时，近20年内发展了许多剖分 和scllachter【“1对其作了改进和完善，把某一点集 技巧。P锄一“给出了基于两点的前沿推进网格生 voronoi图的三角网生长法【l“，其基本步骤是以任 成的基本方法，Eyrnardl2J给出了基于三点的前沿推进 意点为起始点找到距其最近的点构成一个基本边， 法。Joe【3，41和w∞出商1I【51给出了Delalmy剖分的各 种改进方法。两类方法相比较，前沿推进法所用到的 然后找第三点形成Dblaullay三角形，依次处理所有 数据结构较为简单，但浮点运算工作量较大，容易因 边。其找点过程按照Delaunay原则进行，故生成 硬件舍人误差引起失效。De】“my剖分的浮点运算 主要集中在搜索插人新点所破坏三角形的过程中，相 的基本思想，脚yer等人【”191对其作了完善， 对而言，其浮点运算较少，但一般需要全域枚举或采 用递归技术局部寻找被破坏的三角形。枚举寻找工 使之能够剖分大型网格系统。逐点插入算法的基本 作量极大，网格节点多时尤其如此。通过适当的数据 思想是在Delaumy三角形网格系统中插入一点P， 结构实现局部搜索，可使计算工作量大为减少，但须 找到外接圆包含点P的三角形集合s(破坏的三角 采用允许递归调用的程序语言编程，否则缩程很复 形集合)，去掉集合S中的三角形，构成插点空腔 杂。但递归调用使其程序运行内存占用量大为增加， H，连接点P和H的顶点而形成新的三角形网格系 J。Deku-nay算法的另一个缺点 执行效率大为降低【7 统。具体步骤如图1所示。 是判断一点在圆内还是圆外，是基于浮点运算的结 DelauIlay三角形化算法的浮点运算主要集中在 果，有时会因舍入误差导致判断失误。本文结合文献