一种全四边形网格生成方法——改进模板法.pdfVIP

第19卷第1期 　 计 算 力 学 学 报 　 V o l. 19,N o. 1 20 0 2 年2 月 Ch inese Journa l of Com puta tiona l M echan ics February 2002 文章编号:(2002) 一种全四边形网格生成方法——改进模板法 1 2 1 2 李　华 , 　李笑牛 , 　程耿东 , 　吴　杰 ( ( ) ) 1 大连理工大学 工程力学系, 大连 116024; 2 吉林大学 南岭校区 力学系, 长春 130025 摘　要: 首先对全四边形单元网格自动剖分算法中的模板法进行了探讨, 并提出了相应的改进方法。在此基础上 提出了一种新的全四边形单元网格自动生成方法。该方法允许在两个方向上存在网格疏密过渡, 并可以根据单 元的密度要求自动计算亲单元每条边上的结点数, 有效地对局部实施加密处理。 关键词: 网格自动剖分; 全四边形单元; 模板法 中图分类号: O 242. 21　　　文献标识码:A 生成方法——模板法, 在一定程度上解决了第二个 1　引　言 问题。本文在此基础上, 简化了亲单元上模板的分 随着有限元技术的应用范围不断扩大, 越来越 类方法, 并利用单元密度要求来自动计算亲单元的 多并不熟悉有限元技术的人需要使用这门技术。有 每条边的单元分割数, 可以有效地在两个方向上实 限元方法在应用中的一个主要困难是建立有限元 现网格疏密过渡、对局部点和区域进行局部网格加 模型, 即生成有限元网格。建立计算模型需要耗费 密处理, 而对没有网格疏密过渡要求的区域可以生 很大的工作量, 因此如何快速地生成高质量的有限 成和映射法一致的结构化网格。 元网格是一个重要的研究课题。 [ 13 ] 2　算法描述 在目前存在的网格生成算法中 , 三角形单 元网格的生成技术比较成熟。如D elaunay 三角化 2 1　单个正方形亲单元剖分方法 方法、推进网阵法等。相比之下, 虽然全四边形单元 设正方形亲单元每条边上节点数分别为N 1、 的生成方法也得到了一定的发展, 并出现了一些直 N 2、N 3、N 4, 代表的边分别为上、左、右、下, 如图 ( ) 接生成全四边形单元网格的算法, 如四分法、铺路 1 a 所示。每条边的分割数为节点数减 1。在存在 法、中轴转换法等, 但这些方法中需要大量基于角 结点过渡问题时, 为了保证最后生成的单元是全四 度和长度计算的搜索, 效率低、速度慢, 并且容易产 边形单元, 边界上的节点数之和应为偶数。即N 1、 生过渡单元和形状不好的边界单元。 N 2、N 3 和N 4 应满足条件 映射法以其算法简单、速度快、单元质量好、网