学位论文 > 优秀研究生学位论文题录展示

高效可靠的三维约束Delaunay四面体有限元网格生成算法

作　者: 刘岩
导　师: 张洪武；关振群
学　校: 大连理工大学
专　业: 工程力学
关键词: 三维Delaunay三角化边界约束推进波前法(AFT) 有限元法
分类号: TP391.41
类　型: 硕士论文
年　份: 2010年
下　载: 192次
引　用: 0次
阅　读: 论文下载

内容摘要

有限元法是适应计算机使用而发展起来的一种有效的数值分析方法,它在工程实践中的作用已从分析校核扩展到优化设计,并且通过与计算机辅助设计(CAD)相结合成为计算机辅助工程(CAE)的重要组成部分。随着工程实践问题的复杂程度不断增加,集有限元法和误差分析方法于一体的自适应有限元法的应用不断普及。以曲面三角形网格描述的三维实体域为输入,生成四面体划分的三维实体网格是当今有限元网格生成和重要算法之一。可靠、高效的四面体生成算法是实现CAE的基础和前提。本文研究、改进和实现了非结构化自适应有限元网格自动生成的三维Delaunay四面体化方法。首先,本文在绪论中介绍了国内外有限元网格研究的现状和发展的方向。分析了几种当今流行的实体网格生成算法的特点。然后对工程应用最为普遍的Delaunay四面体化方法加以详细论述。在Delaunay四面体化方法的边界恢复、内部点生成和薄元消除这几个方面,本文介绍了其研究现状,并讨论已有算法的特点。第2章,介绍了本文Delaunay四面体化方法的算法框架。以下章节对框架中的三维Delaunay插点算法、约束边界恢复算法、边界一致性恢复算法以及内部点生成算法进行详细介绍。本文在已有工作基础上,对这些算法进行改进,提高Delaunay四面体化算法的稳定性。第3章,提出了一种面向大规模科学计算的三维Delaunay快速插点算法。这一算法提高了单机有限元建模的规模,可在PC计算机生成千万级有限元四面体网格。本文从提高点定位算法速度和提高建立新单元邻接关系算法效率两个方面实现三维Delaunay快速插点。针对三种在工程中最常使用的点集(空间任意点集、正交栅格点集和表面三角面片点集),提出相应的Delaunay快速插点算法,进而扩大算法在实际工程中的适用范围。第4章,介绍了一种有效的三维约束Delaunay三角化的边界恢复算法。本文算法以“任意精度”算法为基础,提出一种可靠的求解空间线段和三角形面片相交判断的算法。然后在这一算法之上,提出解决三维任意域的约束Delaunay三角化问题方法。该方法尽可能减少插入边界上的Steiner点,完成约束边界恢复,保持了实体边界的完整性,解决了经典Delaunay算法不能剖分凹域的问题,从而实现了复杂三维实体的网格剖分。第5章,提出一种能够满足多面体边界几何与拓扑约束的边界一致恢复算法,解决了任意多面体的边界一致四面体网格生成问题。在恢复多面体的几何约束时,边界上可能会引入Steiner点,这样就不满足拓扑约束。对此,本文采用动态规划方法将Steiner点从边界上消除；修复与其相关四面体单元的拓扑关系,以保持原多面体边界的拓扑完整性；并采用扩展的Laplacian光顺算法优化劣质单元。理论上本文算法能够保证完整地恢复任意多面体的边界。第6章,介绍了一种Delaunay四面体化方法与推进波前法(AFT)相结合的三维实体内部点生成算法。这种算法以约束Delaunay四面体化方法生成的模型初始划分为背景网格,并以模型的边界三角形为初始前沿,然后通过AFT生成内部结点,再将这些结点以约束Delaunay的方式插入模型体内。从而,有效的结合了AFT的布点优势和Delaunay整体网格质量较好的特点。最后总结全文,并展望了可以进一步开展的研究工作。

全文目录

摘要  4-6
Abstract  6-10
引言  10-11
1 绪论  11-19
  1.1 有限元法  11-12
  1.2 有限元网格生成算法  12-14
  1.3 三维Delaunay三角化方法研究现状  14-17
  1.4 有限元网格生成算法的度量标准  17-19
2 三维Delaunay三角化方法框架  19-21
3 三维Delaunay快速插点算法  21-33
  3.1 引言  21-22
  3.2 基本概念及算法流程  22-23
  3.3 快速查找BASE单元算法  23-26
    3.3.1 随机点集  23-24
    3.3.2 正交点集  24-25
    3.3.3 表面三角面片点集  25-26
  3.4 建立邻接关系的线性算法  26-27
  3.5 算法分析与实例  27-33
    3.5.1 时间复杂度分析  27-28
    3.5.2 程序效能分析  28-30
    3.5.3 "walk-through"算法平均路径长度  30
    3.5.4 算例  30-33
4 约束边界恢复算法  33-44
  4.1 引言  33-34
  4.2 算法概述  34
  4.3 约束Delaunay插点恢复边界算法  34-38
    4.3.1 约束Delaunay插点算法及其缺点  35-36
    4.3.2 改进的约束Delaunay插点约束边恢复算法  36-37
    4.3.3 改进的约束Delaunay插点约束面恢复算法  37-38
  4.4 基于任意精度插点的恢复边界算法  38-42
  4.5 数值算例  42-44
5 边界一致性恢复算法  44-57
  5.1 引言  44-45
  5.2 边界上Steiner点的消除算法  45-51
    5.2.1 Steiner点的消除算法  45-48
    5.2.2 判断线段的合法性  48
    5.2.3 任意多边形的最优三角形剖分  48-50
    5.2.4 启发式最优三角化准则  50-51
  5.3 边界的拓扑修复算法  51-54
    5.3.1 消除在约束面内的Steiner点并恢复网格拓扑  52
    5.3.2 消除在多个约束面共享的约束边上的Steiner点并恢复网格拓扑  52-54
    5.3.3 消除在仅有一个约束面的约束边上的Steiner点并恢复网格拓扑  54
  5.4 网格优化  54-55
  5.5 网格剖分实例  55-57
6 基于AFT的内部点生成算法  57-66
  6.1 引言  57-58
  6.2 模型尺寸场的构建  58-59
  6.3 前沿结点的生成算法  59-60
  6.4 前沿结点的插入算法  60-61
  6.5 前沿更新算法  61
  6.6 数值算例  61-66
结论  66-68
参考文献  68-73
攻读硕士学位期间发表学术论文情况  73-74
致谢  74-76

高效可靠的三维约束Delaunay四面体有限元网格生成算法

内容摘要

全文目录

相似论文