面向工程与科学计算的表面网格处理方法的若干问题研究

面向工程与科学计算的表面网格处理方法的若干问题研究

论文摘要

网格生成是诸如有限单元法(Finite Element Method,FEM)等各类数值方法的前处理过程,是计算机辅助工程(Computer Aided Engineering,CAE)的共性支撑技术之一,所用网格的质量和生成速度跟数值方法的精度、效率甚至成败都密切相关。作为几何网格的一个子分支,表面网格(包括曲面网格和平面网格)有其特殊性和重要性,不但它的生成跟几何造型直接相关,其自身也是体网格生成程序的输入和各类边界条件的直接加载对象,因而它对最终得到的体网格质量和数值计算结果都有很大影响。在对工程中所用的很多薄壳结构进行壳分析时也要求有高质量的表面网格来支持。此外表面网格在其他很多领域中也有广泛应用。为满足工程与科学计算中对高质量曲面网格的需求,本文系统研究了曲面网格生成、重生成和优化等相关热点和难点问题,提出或改进了一系列相应算法;随着大规模工程与科学计算需求的日益迫切,并行网格生成已成为一个新的研究热点,我们还实现了一个通用的并行平面网格生成框架。具体地说,本文从以下几个方面展开论述:第1章是绪论,简述了本文的研究背景和意义、研究内容及基本框架。第2章是研究综述,回顾了有限元网格生成方法的国内外研究现状,其中跟本文密切相关的几个部分是论述重点,包括网格单元的质量控制与网格过渡、并行网格生成、曲面网格生成和重生成以及网格优化。第3章给出了一些平面全四边形有限元网格的拓扑优化策略。由于单纯的光滑化效果受到网格拓扑结构的限制,在进行网格优化时常常需要将拓扑优化和光滑化结合起来使用以得到更高质量的网格。基于有限元网格的局部拓扑结构,这部分给出的拓扑优化策略被组织成“型—操作”的形式,其中“型”是指一类满足一定约束条件的局部区域网格,而“操作”则是指与特定型相对应的拓扑变换,它能优化局部网格中某些节点的度值,从而优化该局部网格质量,最终实现整体网格的优化。第4章给出了一个通用的并行平面网格生成框架,它包括如下模块:序列化或并行几何分解器、子域图管理模块、并行平面网格生成器以及可选的网格重划分模块。这个并行框架具有可扩展、稳定和高效等基本特性,同时还具有其他一些良好性质:它可对任意序列化及并行几何分解算法、序列化平面平格生成算法实现完全的代码复用;通过引入子域图的概念,并结合静态或动态图划分策略,可使得并行生成的分布式网格具有很好的划分质量,从而减少或消除传统方法中所需要的网格重划分代价,进而提高整个并行模拟过程的效率。存该框架基础上,实现了一个高效、可扩展的并行平面Delaunay网格生成器PDMG-2D,它能利用中等规模并行资源在几分钟内生成包含上亿三角形单元的平面网格。第5章给出了一个复杂组合参数曲面网格生成框架。从系统层面来看,它有自己的面向网格生成的几何模型及其拓扑结构定义,并有相应的几何建模功能,且能和商业CAD文件进行数据交换,并具备常用的几何验证和修复功能。从算法层面来看,它包含两种应用不同参数平面网格生成技巧的间接表面网格生成方法,分别是已有的基于变换矩阵的前沿推进法(Advancing Front Technique,AFT)和新提出的基于黎曼度量的Delaunay方法,它们所生成的表面网格质量都较高。本章详细介绍了框架和算法的各个方面,包括几何模型及其拓扑结构定义;物理空间和参数空间上的黎曼度量和变换矩阵及其相互关系,它们可用来控制参数平面上的网格生成;边界裁剪曲线离散;边界裁剪曲线离散节点在参数平面上的投影计算;参数平面上边界节点方向调整;参数平面网格生成以及对最终得到的整体曲面网格进行方向调整。这些不同环节之间只要满足一定的数据交换标准便可进行模块化封装。第6章给出了一个改进的曲面网格拉普拉斯光滑化算法。拉普拉斯光滑化算法是最基本也最流行的一种方法,它实现简单,效率很高,但也存在一些致命缺陷,如可能产生畸形单元或单元自交现象,会使曲面收缩及曲面几何特征消失。本章针对这些问题一一提出了相应的改进方法,它检测曲面的几何特征并在后面的光滑化过程中进行特别处理以保证它们不被丢失;它考虑了一些的新的几何因素,并根据微分几何中的极小曲面原理,通过求解一个带约束的优化问题来防止无效单元的产生;通过投影算法保证光滑化后的网格节点仍落在原始网格上,因此曲面也不会被收缩。第7章给出了一个保特征的曲面网格重生成方法,它可以根据用户指定的密度通过拓扑变换对网格的不同区域分别进行粗化细化,这些拓扑变换直接在离散曲面上操作而没有一个连续曲面支持,因而实现较为简单,最后通过边交换操作和本文第6章中给出的改进的拉普拉斯光滑化算法对网格质量进行优化,得到最终高质量的计算网格。本文最后的第8章总结全文,并展望了可以进一步开展的研究工作。

论文目录

  • 摘要
  • Abstract
  • 第1章 绪论
  • 1.1 研究背景和意义
  • 1.1.1 几何网格和数值方法
  • 1.1.2 结构化网格和非结构化网格
  • 1.1.3 网格质量控制的基本要求
  • 1.1.4 曲面网格生成、重生成和优化
  • 1.1.5 并行网格生成
  • 1.2 研究内容和主要贡献
  • 1.3 基金资助
  • 第2章 研究综述
  • 2.1 引言
  • 2.2 有限元网格生成算法研究现状
  • 2.2.1 网格单元的质量控制与网格过渡
  • 2.2.2 通用的结构化有限元网格生成方法
  • 2.2.3 通用的非结构化有限元网格生成方法
  • 2.3 并行网格生成的基本概念
  • 2.3.1 并行网格生成算法的评价标准和分类
  • 2.3.2 数据并行和任务并行
  • 2.3.3 数据分解
  • 2.3.4 数据映射和负载平衡策略
  • 2.4 曲面网格生成和重生成
  • 2.4.1 面向网格生成的几何模型获取
  • 2.4.2 有限元网格和曲面网格数据结构介绍
  • 2.4.3 曲面网格生成
  • 2.4.4 曲面网格重生成
  • 2.5 有限元网格优化方法
  • 第3章 平面全四边形有限元网格的拓扑优化策略
  • 3.1 引言
  • 3.2 方法概述
  • 3.3 内部单元优化
  • 3.3.1 减点操作
  • 3.3.2 增点操作
  • 3.3.3 对角线交换操作
  • 3.3.4 组合操作
  • 3.4 边界处理
  • 3.4.1 减点操作
  • 3.4.2 增点操作
  • 3.5 算法实现
  • 3.6 实验结果
  • 3.7 结论
  • 第4章 一种通用的并行平面网格生成框架
  • 4.1 引言
  • 4.2 总体框架
  • 4.3 并行几何分解
  • 4.4 基于子域图分解的并行网格生成
  • 4.5 实验结果分析
  • 4.5.1 并行框架的可扩展性和效率
  • 4.5.2 并行框架的稳定性
  • 4.5.3 分布式网格的划分质量
  • 4.5.4 数值模拟例子
  • 4.6 小结
  • 第5章 复杂组合参数曲面网格生成的框架与算法实现
  • 5.1 引言
  • 5.2 面向网格生成的几何模型描述
  • 5.2.1 参数曲线和Ferguson参数曲线
  • 5.2.2 参数曲面和Ferguson参数曲面
  • 5.2.3 几何模型的拓扑信息描述
  • 5.3 网格单元质量指定方式
  • 5.3.1 黎曼度量
  • 5.3.2 变换矩阵
  • 5.3.3 参数曲面的黎曼度量
  • 5.4 组合参数曲面网格生成总体框架和算法流程
  • 5.5 三维曲线离散
  • 5.6 计算曲面区域边界离散节点存参数平面上的投影
  • 5.7 调整参数平面上边界节点的方向
  • 5.8 用调整后的AFT方法在参数平面上生成网格
  • 5.9 用黎曼度量下的Delaunay方法生成参数平面网格
  • 5.9.1 黎曼度量下的2维Delaunay网格生成方法总体流程
  • 5.9.2 内部节点产生
  • 5.9.3 黎曼度量下的Delaunay准则和内部节点插入
  • 5.9.4 网格优化
  • 5.10 对最终得到的整体曲面网格进行方向调整
  • 5.11 数值算例
  • 5.11.1 AFT方法生成的网格例子
  • 5.11.2 Delaunay方法生成的网格例子
  • 5.12 小结
  • 第6章 改进的曲面网格拉普拉斯光滑化方法
  • 6.1 引言
  • 6.2 光滑化方法总体流程
  • 6.3 节点分类算法
  • 6.4 节点移动算法
  • 6.4.1 计算初始偏移向量
  • 6.4.2 调整偏移向量
  • 6.4.3 将新节点投影到原始曲面网格
  • 6.5 实验结果分析
  • 6.6 小结
  • 第7章 保特征的曲面网格重生成
  • 7.1 引言
  • 7.2 几何特征检测
  • 7.3 曲面网格重生成
  • 7.4 网格优化
  • 7.5 实验结果
  • 7.6 结论
  • 第8章 总结与展望
  • 8.1 总结
  • 8.2 展望
  • 参考文献
  • 攻读博士学位期间主要的研究成果
  • 作者简历
  • 致谢
  • 相关论文文献

    标签:;  ;  ;  ;  ;  ;  ;  ;  ;  

    面向工程与科学计算的表面网格处理方法的若干问题研究
    下载Doc文档

    猜你喜欢