基于表面三角网格的四面体剖分-邓俊辉.PDF

基于表面三角网格的四面体剖分 杨义军（035456 ） 龚洁晖（035453 ） 王文珂（035452 ） yang-yj03@ 摘 要 四面体剖分可以表达模型的拓扑特性，简化模型的分析。对模型特别是工业模型进 行四面体剖分是有限元分析的一个重要步骤。本文详细描述了一种四面体剖分生成算法。首 先对 Gems™生成的模型用 Ansys™进行表面三角剖分；对 Ansys™处理过的表面三角网格 进行三角形的空球检测，对不满足空球特性的三角面片进行必要的边变换；在此基础上，根 据晶体网格(FCC) 加点规则在模型内部加入内部点；最后算法采用波前法和 Delaunay 准则 生成模型的四面体剖分。本文最后给出了部分剖分结果和中间结果，对算法中出现的一些问 题进行了分析。 关键词 晶体网格 八叉树 Delaunay 剖分 波前法 1 问题背景 图 1 （1） 正方体 （2 ）正方体一种四面体剖分 （3 ）另一种四面体剖分 模型的拓扑特性可以用四面体剖分很好的表示，所以四面体剖分是有限元分析等许多工 程应用的基础方法。一般模型的四面体不是唯一的，如图1，即使对一个正方体进行四面体 剖分，剖分结果也不唯一。一个模型在没有加入内部点的情况下是不是存在四面体剖分是 NP-Hard 问题。对于图2 中的多面体，没加入附加点的情况下不存在四面体剖分。但是每一 个多面体都能在加入内部点（最坏情况 O(n)）的条件下进行四面体剖分。一个好的四面体 剖分应满足以下要求： 图2 不可四面体剖分多面体 合法性： 四面体不能相交，在边界上的结点均应作为剖分的结点，不可丢弃。 相容性： 四面体必须落在待分区域内部，且所有四面体的并集等于待分区域。 逼近精确性：待分区域的顶点（包括特殊点）必须是单元的结点，待分区域的边界（包 括特殊边及面）被单元边界所逼近。 良好的形状：单元最佳形状是正多边形或正多面体，单元之间过渡应相对平稳。 自适应性： 四面体剖分的密度应符合待分区域的物理特性，在几何尖角处、应力温度 等变化大处较密集，其他部位较稀疏。 采用几何方法的四面体剖分技术主要有几何分解法、结点连元法。在产生结点的同时也 确定结点间拓扑关系的方法称为几何分解法，常用的有两种：递归法和迭代法。几何分解法 的关键技术在于新结点的生成方法，最大优点是在离散物体时考虑网格单元的形状和大小， 所生成的网格单元形状和分布均较好；最大缺点是自动化程度低，不利于复杂件网格生成。 结点连元法是先生成结点，然后连接结点构成单元。最常用的是DT 法和AFM法。DT法 的基本原理是利用 Voronoi 图和 Delaunay 三角剖分的对偶性质，连接相邻 Voronoi 多边形 的内核点构成 Delaunay 三角剖分。DT 法的最大优点是遵循 Delaunay 三角剖分“最小角最 大”和“空球”准则，同时满足全局和局部最优；但直接的 Delaunay 三角剖分只适用于凸 域，不适用于非凸域。AFM法的基本原理是假设区域的有向离散外边界集和前沿边界集已经 确定，按某种条件沿区域边界向区域内部削切四面体，同时不断更新前沿边界集，直到空集。 AFM法的关键技术有两个：一是区域的边界离散和内结点的合理生成。二是由当前面选择内 部点生成新的四面体的规则。AFM法最大优点是不仅在区域内部而且在区域边界所生成的网 格单元形状均优良，网格生成全自动。 本文实现了一种基于内部点添加和波前法的四面体剖分算法。算法利用了晶体网格的加 点规则和 Delaunay 准则，生成的四面体较好地保持了原来多面体的边界特性。算法的时间 复杂性为 O(n2)。 2 算法描述 本文首先对表面的三角形进行调整，使得每一个三角形的空球中不包含面上的其他点； 接着对面调整后的多面体根据晶体网格规则插入内部点；最后对整个多面体和内部点用波前 法和Delaunay 准则进行四面体剖分。 2.1 空间的八叉树划分 三角形调整时，为了加速三角形外接球包含点的检测以及内部点的添加，本文首先对多 面体进行了空间的八叉树划分。八叉树节点分成三类： 八叉树外部节点：完全位于多边体的外部，不包含多面体上的点，