基于stl格式文件的全四边形网格生成方法

《基于stl格式文件的全四边形网格生成方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、基于STL格式文件的全四边形网格生成方法*国家973计划(2004CB719402)和教育部跨世纪优秀人才计划资助项目。.*联系作者：陈涛Email:danielchen2005@gmail.com陈涛+,高晖,李光耀(汽车车身先进设计制造国家重点实验室湖南大学长沙410082)摘要:提出一种以STL格式文件所描述的离散几何模型为基础，使用改进的铺路法自动生成全四边形网格的方法。重建STL文件的拓扑结构数据，而后进行模型的内外边界搜索及初始化布点。算法依次向模型内部加入新的节点以生成新的四边形网格单元，直至把模型内部全部覆盖。原始几何模型中的特征线被提取出来，并在网

2、格生成阶段将其作为内部孔洞处理，减小了生成网格模型所导致的离散误差，在特征所在区域生成质量较佳的网格。使用一种联合Laplacian方法与小种群遗传算法（μGA）的网格光顺方法，可以有效地纠正反转单元、退化单元等形态质量很差的单元。多个算例验证了本文提出方法的有效性。关键词:STL；网格剖分；特征提取；网格光顺；小种群遗传算法0.引言随着汽车碰撞有限元仿真分析在国际上各大汽车公司的广泛应用，对有限元网格模型的要求也越来越高。一方面，CAD模型中大量的细节需要在网格模型中得以保留，以便更加真实的模拟汽车碰撞过程；另一方面，为了提高碰撞仿真计算精度，保证求解的稳定性，要

3、求网格模型中的单元具有更高的质量。目前在主流CAD（ComputerAidedDesign计算机辅助设计）软件使用NURBS（非均匀有理b样条）作为几何建模内核，商业化网格生成软件需要通过接口软件从软件中获取模型进行网格剖分。但通过接口软件进行数据交换时常常会产生数据丢失或数据错误，严重影响了网格剖分的进行。近年来，面片格式的几何模型表述方式得到了广泛应用，逐渐成为NURBS表述的一种替代方式，如在快速原型制造领域使用的STL（Stereolithography）文件格式，它使用三角形来表述几何模型，格式简单且不易出错，因此更多的研究者把倾向于使用面片模型作为网格剖

4、分的输入模型[1]。在汽车碰撞和薄板冲压等强非线性问题的有限元仿真计算中，优先使用四边形单元，因为相比于三角形单元，四边形单元具有更高的求解精度和计算效率。非结构化四边形网格的生成方法分为直接法和间接法两大类[2]。直接法包括影射法、铺路法、前沿推进法、四叉树法等；间接法则首先将待划分的区域划分为三角形单元，接着用不同的算法将三角形单元转换成四边形单元。Blacker[2]在1991年提出铺路法(Paving)属于直接法，具有适用性强，生成网格模型的质量高等优点，现已被广泛运用到商业化有限元分析软件的前处理程序中，如Ansys、HyperMesh等。提出一种以STL

5、格式几何模型的全四边形网格生成算法。克服了使用NURBS几何模型时所遇到的繁琐的模型修补问题。针对三角面片的特点，改进了原有的铺路算法，并能够更有效的控制网格单元尺寸的疏密过渡。原始几何模型中的重要几何特征在网格模型中均能够较好的保留。本文提出联合使用Laplacian方法与小种群遗传算法（μGA）的网格光顺方法，在合理的计算时间内，大幅提高了模型的单元质量。1．STL文件的拓扑重建和特征提取STL格式的数据文件中包含的三角面片，并不能满足有限元计算要求。从图1可以看到，模型包含有大量的狭长的三角形，无法满足碰撞仿真计算的精度要求。因此在本文所进行的研究中，三角面片

6、模型仅作为生成四边形网格所需的背景网格。Béchet[2]使用STL格式的数据文件生成了满足有限元计算要求的三角形网格，而本文的研究关注于全四边形网格的生成。图1某车型行李箱盖面片模型1．1模型的拓扑重建与边界初始化STL格式的数据文件中提供了每个三角形的三个顶点的坐标和三角形的法向矢量。从模型的拓扑层次结构上看，STL格式中只有顶点—面片两个拓扑层次，二者之间也仅仅存在从属关系（顶点从属于某个单元）。进行STL模型拓扑重构的目的是要建立顶点—边—面片的拓扑层次结构，并提取拓扑实体之间的从属、相邻关系。依次从文件中读取三角面片，记录它的三个顶点的坐标和法向，并搜寻已

7、存储的节点中是否存在重叠顶点，如果有则进行合并。由此建立顶点与顶点、面片与面片之间的相邻关系。读入完成之后，根据面片之间的相邻关系创建边，建立顶点—边—面片的完整拓扑层次结构关系。模型的边界包括内边界(指模型内部的孔洞)和外边界。边界由仅从属于一个三角形的边组成的。边界搜索分三步进行：1．首先搜索出模型中所有的仅从属于一个三角形中的边；2．将它们按照首尾相接的顺序排列，构成一条或多条边链，每一条链代表一个边界；1．最后要正确区别内外边界。假定模型的法向是Z轴正向，则外部边界链是逆时针方向的，而内部孔洞是顺时针方向的。铺路算法需要首先对所有的边界布上偶数个点，并且