计算机对任意三维区域作正规网格剖分算法.pdf

第 ! 卷 第 # 期 湘潭师范学院学报（自然科学版） 8’9 ! -9# !$$% 年 % 月 !#$%’ ( )*%+,% -$.’ /%*01$2*,（3 -,#$’ 45*1%51 67*,*%） :$9 !$$% 计算机对任意三维区域作正规网格剖分算法! 黎自强， 杨奇为 （湘潭大学 信息工程学院，湖南 湘潭 ###$） 摘 要：详细介绍了一种三维复杂区域的正交差分网格的剖分方法。这种方法具有计算精度高、剖分速度快等特点。 关键词：三维区域；差分方法；偏微分方程式；网格剖分 中图分类号：’(%$#)* 文献标识码：+ 文章编号：#*,# - $!%（# !$$%）$# - $$#$ - $% # 问题提出 应用数值方法求解各类数学物理问题时，要求对求解区域进行有限剖分，然后在网格节点上将连续的问 题离散。然而给出网格节点的坐标及在节点上的信息是一种繁琐的工作。手工计算一般只适应二维区域， 三维区域则很麻烦，而且计算起来、填写数据很容易出错，如果要应用计算机辅助设计，选择最优结构参数时 需要对同一几何区域反复计算，则就需要一种快速而准确的网格剖分方法。因而网格剖分特别是任意三维 区域的网格剖分是工程计算一类急需解决的问题。 网格剖分有非正规网格剖分，如二维三角形网格和三维四面体网格以及二维区域的曲线网格，在这方面 的研究有文献［#］和文献［!］。另一种网格剖分是二维的矩形（正方形）网格和三维长方体（立方体）网格剖 分。关于这方面有文献［%］，其思想是从数据文件中读取几何实体的面方程、边的起末点、节点坐标及体界 盒、面界盒后存入相应的数组内，剖分从 方向开始，用 截几何实体，如果剖分第 层，则面方程是 . . / .012 3 . / （ ） ，先根据面界盒判断哪些面通过平面 （ ） ，然后求出这些面与平面 （ .012 4 5 - # 6. . / .012 4 5 - # 6. . / .012 4 5 - #）6. 的交线方程，这些线段必然构成一个或几个封闭的二维平面曲线。在 5 层的 7 - 8 平面上，同 . 方向分 割一样，按 方向用线 与二维封闭曲线相交，如果剖分到 层，则线方程为 （ ） ，首先判 8 8 / 8012 9 8 / 8012 4 9 - # 68 断封闭曲线段中哪些通过直线段 ，然后求出这些曲线段与直线段 （ ） 的交点，再使用排 8 / 89 8 / 8012 4 9 - # 68 序法判断并排除奇异点，最后给该点赋网格节点编号（1，9，5）。经过这样一些计算能将三维区域内部剖分出 计算网格。这种剖分方法： （#）需要反复解方程，工作量大； （!）当边界有一段是直线而又正好与网格线重合时，那么解为无穷多个，从而引起一些麻烦； （%）当区域为非凸区域时，网格线与边界的交点超过两个，这时判断哪些点在区域内也比较困难； （）如果区域边界给出的是离散点的数据时，这种方法无法直接进行剖分。 ! 剖分思想 设 是要剖分的任意三维区域，它的边界曲面用 （ ，）表示，网格的长、宽、高分别是 、 、 ，用一 :