《基于Rhino的有限元网格离散技术研究_毕业论文.doc》由会员分享，可在线阅读，更多相关《基于Rhino的有限元网格离散技术研究_毕业论文.doc（61页珍藏版）》请在三一文库上搜索。
1、基于Rhino的有限元网格离散技术研究 上海工程技术大学毕业设计论文上海工程技术大学毕业设计论文题 目：基于Rhino的有限元网格离散技术研究学 院： 电子电气工程学院 专 业： 自动化 58目 录摘要 1ABSTRACT 20 引言 41 有限元法 51.1 有限元法的基本原理 51.1.1 有限元法基本思想 51.1.2 有限元法几个基本概念 51.2 有限元求解的基本步骤 61.2.1 区域离散 61.2.2 插值函数的选择 81.2.3 方程组公式的建立 81.2.4 方程组的求解 111.3 有限元法的应用 132 Rhino建模技术及网格剖分 142.1 Rhino软件简介 142
2、.2 Rhino实体建模技术 162.2.1基本几何体创建 162.2.2 实体工具 182.3 网格剖分 212.3.1 有限元网格剖分的基本原则 212.3.2 有限元网格生成的方法 232.3.3 网格质量的度量准则 272.3.4 Rhino网格工具 282.4 模型及四面体网格剖分实例 332.4.1 矩形腔体 332.4.2 圆柱腔体 352.4.3 阶梯模型 363 有限元计算 393.1 矩形腔体计算 393.2 圆柱腔体计算 404 结束语 415 参考文献 426 译文 447 原文说明 55摘 要本文首先对有限元法的基本原理进行简要的叙述，指出了四面体网格在有限元分析中的
3、使用及四面体网格剖分的难点和现状，使用了Rhino软件建立三维几何模型并对模型进行四面体网格自动剖分，提取网格数据，再转换成有限元计算所需数据，并进行有限元计算。本文的主要内容有以下几个方面：（1）对有限元法基本原理的学习,其中重点偏向三维有限元分析。以及有限元网格剖分的研究，其中重点是对四面体网格自动剖分进行了研究和叙述，并提出了现今四面体网格单元质量的度量准则和几种常用的算法。（2）在Rhino建模功能的基础上，建立简单矩形腔体模型和其他复杂模型，本利用Rhino的网格工具对其进行四面体网格自动剖分，并提取网格数据。在Fortran语言环境下将提取网格数据转换成有限元计算所需的网格数据。（
4、3）最后就是在运用前面所提到有限元分析的基本求解步骤和网格剖分的基础理论条件下，利用我们提取转换的有限元计算所需的网格数据对矩形腔体模型进行相关的有限元计算。关键词：有限元分析，网格自动剖分，Rhino建模，四面体网格Study of Finite Element Mesh Discretization Technique in RhinoABSTRACTFirstly,basic principles of the finite element method provides a brief description in this paper. It was then pointed out
5、 the advantage of tetrahedral mesh on analyze with the finite element method and the difficulty of the automatic generation technique of tetrahedral mesh. a three-dimensional geometric model and the model tetrahedral mesh automatic generation in Rhino, and extraction grid data, converted to the requ
6、ired data finite element calculation and finite element calculations. The main contents of the following aspects:(1) Become skilled at the basic principles of the finite element method , Which tend to focus on three-dimensional finite element analysis .And the finite element mesh research, which foc
7、uses on tetrahedral mesh automatic generation were studied and described. And put forward today tetrahedral mesh element quality measurement criteria, and several commonly used algorithms.(2) Based on modeling capabilities in Rhino, the establishment of a simple rectangular cavity and other complex
8、models, and use the tools for their Rhino tetrahedral grid automatic mesh grid and grid data extraction. In Fortran language environment will extract data into a finite element mesh to calculate the required grid data.(3) The final step is in the use of the previously mentioned finite element analys
9、is of the basic steps and solve the basic theory of meshing conditions, using our transformation to extract the desired grid finite element data model for the rectangular cavity associated with FEM .Key words: finite element, automatic mesh subdivision, Rhino modeling, tetrahedral mesh基于Rhino的有限元网格离
10、散技术研究0 引言得益于计算机技术的高速发展，有限元法已迅速从工程结构强度分析计算扩展到几乎所有的科学技术领域，成为一种应用广泛并且实用高效的数值分析方法。如今在工程分析法中，有限元分析法成为一种非常有效的方法，能够解决多种复杂工程问题。有限元法应用已经由弹性力学平面问题扩展到空间问题，由静力平衡问题扩展到稳定问题，动力问题和波动问题。从固体力学发展到流体力学、传热学、电磁学等连续介质力学领域。随着计算机辅助前处理技术和计算机辅助后处理技术的发展，以及有限元法与计算机的结合，有限元法必将会获得更广泛的应用。然而在进行有限元分析之前，必须对将要分析的几何模型进行网格离散，为保证计算的准确性，网格
11、剖分需满足原来的模型以及有限元分析的要求，所以有限元网格的合理剖分在有限元法中占据着重要的位置。但是和有限元高度发展的计算功能不相符的是，有限元与处理功能的网格自动剖分很弱。因此，对有限元网格生成技术的分析和研究，对发展有限元技术及实用化有重要的作用。本文利用Rhino的建模功能建立复杂几何模型，再对其进行有限元四面体网格剖分，并提取网格数据，并把网格数据转换成高质量的网格数据文件，再用在有限元计算。1 有限元法1.1 有限元法的基本原理1.1.1 有限元法基本思想我们实际要处理的对象都是连续体，在传统设计思维和方法中，是通过一些理想化的假定后，建立一组偏微分方程及其相应的边界条件，从而求出在
12、连续体上任一点上未知量的值。因为点是无限多的，存在无限自由度的问题，很难直接求解这种偏微分方程用来解决实际工程问题，因此需要采用近似方法来处理。其中最主要的是离散化方法，把问题归结为只求有限个离散点的数值，把无限自由度问题变成有限个自由度。把一个连续体分割成有限个单元，即把一个复杂的结构看成由有限个通过结点相连的单元组成的整体，先进行单元分析，然后再把这些单元组合起来代表原来的结构，以得到复杂问题的近似数值解。这种方法称为有限元法。有限元法是一种以计算机为手段，通过离散化将研究对象变换成一个与原始结构近似的数学模型，再经过一系列规范化的步骤以求解相应参数的数值计算方法。所谓离散化就是将一个连续
13、体分割成若干个通过结点相连的单元，这样一个有无限个自由度的结构就变换成一个具有有限个自由度的近似结构。该过程还包括对单元和结点进行编码以及局部坐标系和整体坐标系的确定。1.1.2 有限元法几个基本概念（1）单元：将求解工程结构看成是由许多小的、彼此用点连接的基本构件，称为单元。在有限元法中，单元用一组结点间相互作用的数值和刚度系数矩阵来描述。（2）结点：单元与单元之间的连接点，称为结点。在有限元法中，结点就是空间中的坐标位置，它具有物理特性，且存在相互物理作用。（3）有限元模型：有限元模型真实系统理想化的数学抽象。有一些形状简单的单元组成，单元之间通过结点连接。每个单元的特性是通过一些线性方程
14、来描述的。作为一个整体，所有单元的组合就形成了整体结构的数学模型。（4）有限元分析：利用数学近似的方法对几何模型进行模拟。并利用简单而又相互作用的单元，用有限数量的未知量去逼近无限未知量的真实系统。1.2 有限元求解的基本步骤1.2.1 区域离散在任何有限元分析中，区域离散（设置区域为）是第一步，也是最重要的一步，因为区域离散的方式将影响计算内存需求、计算时间和数值结果的精度。在这一步骤中，全域被分成许多小区域，用（e=1,2,3,M)表示，这里M表示子域总数。这些子域通常被称为单元。对于实际上是直线或曲线的一维区域单元通常是短直线，如图1.1（a)。对于二维区域，单元通常是小三角形或矩形，如
15、图1.1（b）。对于三维区域，可剖分成四面体、三棱柱或矩形块，如图1.1（c），其中四面体是最简单、最适合离散任意体积区域的单元。在图1.2中，给出了表示三维区域用四面体网格有限元离散化的一个立方体。(a) 直线（一维）(b) 三角形和矩形（二维）(c) 四面体、三棱柱和矩形块（三维）图1.1 基本有限元图1.2 立方体剖分效果在多数有限元解中，问题是用与单元有关的结点上的为指数函数表达的。例如，线性三角形单元有三个结点，线性四面体单元有四个结点。再次有必要描述一下这些结点。一个结点的完整描述应包括它的坐标值、局部编码和全局编码。结点的局部编码表示它在单元中的位置，而全局编码表示它在整个系统中
16、的位置。因为区域离散过程完全可以与其它步骤分开，所以，通常将它当作一项预处理工作。1.2.2 插值函数的选择有限元分析的第二步是选择能近似表达一个单元中未知解的插值函数。通常，插值函数可选择为一阶（线性）、二阶（二次）、或高阶多项式。尽管高阶多项式的精度较高，但通常得到的公式也比较复杂。因此，简单且基本的线性插值仍被广泛采用。一旦选定了多项式的阶数，我们就能导出一个单元未知解的表达式。以单元为例，可得到下列形式： (1.1)式中，是单元中的结点数；是单元结点的；是插值函数，通常也称为展开函数或基函数，的最高阶被称为单元的阶。例如，若是线性函数，则单元是线性单元。函数的重要特征是：它们只有在单元
17、内才不为零，而在单元外均为零。1.2.3 方程组公式的建立第三步（也是有限元分析中的主要步骤）是导出方程组的公式。里兹变分和伽辽金方法均可用于这种目的。求解公式前我们先说明边值问题，典型的边值问题可用控制微分方程和边界条件来定义。微分方程可表示为： （1.2）式中，是微分方程，激励或加强函数，是未知量。（1） 里兹方法的求解公式让我们考虑（1.2)式中的问题。为了简单起见，假设问题是实数值得。可得泛函数表示为： （1.3）式中是组成全域的单元数，而 （1.4）将（1.1）式代入（1.4）式，得到 （1.5）写成矩阵形式 （1.6）其中是矩阵，是的列向量，它们的元素为 （1.7）和 （1.8）因
18、为是自伴算符，所以，单元矩阵是对称的。将（1.6）式代入（1.3）式，可得到 （1.9）通过求和运算，并采用全局结点编码，上式可写成 （1.10）式中是对称矩阵；N是未知量或结点总数；是的未知向量，其中元素是未知展开系数；是的已知向量。应用驻点条件，令对的偏导数为零： （1.11）即可的到方程组。因为是对称的，因此，（1.11）式变成 （1.12）或写成矩阵形式： （1.13）（2） 伽辽金方法的求解公式上述方程组也可以用过伽辽金方法导出，对于（1.2）式，第e个单元的残数加权为 （1.14）将（1.1）式代入（1.14）式后得到 （1.15）也可写成矩阵形式 （1.16）这里，矩阵元素和分别
19、于（1.7）式和（1.8）式的相同。因为这里的算符不必是自伴的，所以，单元矩阵也不必是对称的。既然与一结点有关的展开函数遍及所有和该结点直接相连的单元，那么，与结点有关的残数加权是对所有直接和结点相连的单元求和。所以，利用局部和全局编码的关系，我们可以扩展（1.16）式，然后将它对每一单元求和，得到 （1.17）式中，。（1.17）式中求和号内的所有向量和矩阵都按前面描述的方式展开了。令（1.17）式等于零，可得方程组 （1.18）在方程组（1.18）式能求得定解之前，必须应用所需的边界条件。有两类边界条件经常出现：一是狄利克雷边界条件，另一类是齐次罗曼边界条件。第一类边界是必要边界条件，它必
20、须显式地加强加在计算中；第二类边界条件通常在求解过程中隐含地自动满足。这是这种原因，第二类边界条件通常称为自然边界条件。可以看出，在这一步骤中，实际上有三个步骤。首先，应用两种方法中的任一种，写出单元方程组（1.6）式或（1.16）式。其次，将单元方程组对所有单元求和，得到方程组，这个过程叫做组合。最后。我们应用边界条件来得到方程组的最终形式。1.2.4 方程组的求解方程组的求解是有限元分析的最后一步。最终的方程组是下列两种形式之一： （1.19）或者 （1.20）方程（1.19）式是确定型的，它是从非齐次微分方程或非其次边界条件或从它们两者兼有的的问题中导出的。在电磁学中，确定性方程组通常与
21、散射、辐射以及其它存在源或激励的确定性问题有关。而方程（1.20）式是本征值型的，它是从齐次控制微分方程和齐次边界条件导出的。在电磁学中，本征方程组通常与诸如波导中波传输和腔体中谐振等无源问题有关，在这种情形下，已知向量为零，矩阵可写成的形式，这里表示未知的本征值。一旦解出的方程组，就能计算出所需的参数。例如电感、电容、输入阻抗和散射或辐射图等，并能用曲线、图形或彩色图片等形式表示结果，这些形式更有意义并易于解释。这个最后步骤也称为后处理过程，完全可以同其它步骤分开处理。对电磁学以及其他工程领域中的许多问题，有限元公式的最终结果是一组线性代数方程组如（1.19）式那样，有许多的算法可解（1.1
22、9）式。然而并不是所有的解法都适合用于求解有限元方程组，这里我们介绍一下常用的三角分解法和共轭梯度法。（1） 三角分解法三角分解法是将原正方矩阵分解成一个上三角形矩阵或是排列的上三角形矩阵和一个下三角形矩阵，这样的分解法又称为LU分解法。它的用途主要在简化一个大矩阵的行列式值的计算过程，求反矩阵，和求解联立方程组。不过要注意这种分解法所得到的上下三角形矩阵并非唯一，还可找到数个不同的一对上下三角形矩阵，此两三角形矩阵相乘也会得到原矩阵。（2） 共轭梯度法共轭梯度法的基本思想是把共轭性与最速下降方法相结合，利用已知点处的梯度构造一组共轭方向，并沿这组方向进行搜素，求出目标函数的极小点。根据共轭方
23、向基本性质，这种方法具有二次终止性。共轭梯度法不仅是解决大型线性方程组最有用的方法之一，也是解大型非线性最优化最有效的算法之一。在各种优化算法中，共轭梯度法是非常重要的一种。其优点是所需存储量小，具有步收敛性，稳定性高，而且不需要任何外来参数。正是基于共轭梯度法的优点，本文所涉及的有限元算例都是运用共轭梯度法求解。1.3 有限元法的应用目前，有限元法计算在许多工程计算中都有应用，例如流体力学、电磁学、工程结构等。我们这里就简单地介绍一下，有限元法在电磁电气方面的应用，如电机的电磁分布、电磁力、变形、转子运动、动态变化过程以及电子电力装置相结合等情况下的分析和特性预测及电机参数的计算等。有近些年
24、比较新的内容。如传动调节器的计算等。还有变压器及其他电力系统元件，如高压线圈、输电线电缆、接地系统、输电线配电线的外部磁场、电晕等的分析。也有一些研究项目，例如核磁共振成像系统磁体、磁流体发电、电池物体发射以及电磁对人体影响等计算分析这些应用各具特色，它们的计算分析方法也不尽相同。总的来说，可以根据这些应用中电磁现象的变化快慢程度分为三大类型，即静磁场、涡流场和高频场等。三维有限元静磁场计算目前还是比较成熟的。三维有限元涡流场的计算目前的目标主要是如何来减少计算量同时提高计算精度。高频场的有限元计算有其特殊性。出现了相应的一些新的方法。对于高频场，位移电流不能忽视。其有限元法计算采用矢量磁位，
25、仍具有介质交界处连续的特点，因此介质交界条件自动满足，不用考虑。对于高频场的处理有两种方法：一种是类似于涡流场情况，这是求解获得的是对应某个频率下的场量分布；另一种专门用于高频场的分析，即为特征值和特征矢量的计算。通过对特征值的计算分析，既能获得对于某个频率下场分布的完整的解，而且也能考察系统的其它一些特点如共振特性等。它特别适合于谐振腔、天线等装置设备的电磁有限元计算。2 Rhino建模技术及网格剖分2.1 Rhino软件简介Rhino，中文名称犀牛，是一款超强的三维建模工具，大小才几十兆，硬件要求也很低。不过不要小瞧它，它包含了所有的NURBS建模功能，用它建模感觉非常流畅，所以大家经常用
26、它来建模，然后导出高精度模型给其他三维软件使用。强大的曲面建造和优质的模型结构，使得它能够轻松自如的创建非常精确的工业产品、建筑、首饰、家具的三维模型。从设计稿、手绘到实际产品，或是只是一个简单的构思，Rhino所提供的曲面工具可以精确地制作所有用来作为彩现、动画、工程图、分析评估以及生产用的模型。Rhino可以在Windows系统中建立、编辑、分析和转换NURBS曲线、曲面和实体。不受复杂度、阶数级以及尺寸的限制。Rhino也支援多边形网格和点云特色包括：不受约束的自由造形3D建模工具以往您只能在二十至五十倍价格的同类型软件中找到这些工具。让您可以建立任何可以想象的造形。精确完全符合设计、快
27、速成形、工程、分析和制造从飞机到珠宝所需的精确度。兼容性兼容于其它设计、制图、CAM、工程、分析、着色、动画以及插画软件。读取和修复难以处理的IGES档案。容易使用非常容易学习使用，让您可以专注于设计与想象而不必分心于软件的操作上。高效率不需要特别的硬设备，即使在一般的笔记型计算机上也可以执行。经济实惠普通的硬设备，容易上手，价格相当于一般的Windows软件，并且不需额外的维护费用。目前广泛的应用与韩国、台湾等工业设计发达的地区。最后我们给出两个用Rhino创造的3D作品，如图2.1所示的汽车模型和图2.2所示的变形金刚。图2.1 汽车模型图2.2 变形金刚模型2.2 Rhino实体建模技术
28、2.2.1基本几何体创建Rhino中的基本几何体包括：立方体、球体、圆柱体、椭球体、圆管、圆锥、抛物面椎体、圆环等。如图2.3所示。图2.3 Rhino中的基本几何体关于这些基本几何体命令的使用，操作起来很简单，因此只简单叙述下，见图2.4所示。图2.4 实体命令以上的命令中，简单地介绍一下操作拉伸曲线出实体命令建立一个三棱体实体，具体过程，先用Rhino的直线工具在主体操作窗上，画一个如图2.5所示的三角形，再点击拉伸曲线出实体命令之后，选中三角形，按完成键，拖动三角形到合适的高度即可。建立的三棱体如图2.6所示。图2.5 三角形图2.6 三棱体2.2.2 实体工具实体工具集中主要是有三类工
29、具：布尔运算工具、抽面工具、倒角工具。实体工具集功能及主要使用方法如图2.7所示。图2.7 实体工具命令关于实体工具的运用，我们用实际几个的例子来说明，实体布尔运算、自动建立实体和实体倒角等实体工具命令的使用。（1） 布尔运算实例图2.8 立方体和球体图2.9 实体布尔运算并集后的效果（2） 实体倒角图2.10 实体倒角前图2.11 实体倒角后的效果2.3 网格剖分网格剖分作为建立有限元模型的一个重要环节，要求考虑的问题多，需要的工作量大，不同的网格剖分方式会对计算规模、计算结果和计算精度产生很大的影响。对网格剖分基本原则是否遵守、网格生成算法的不同选择和网格剖分单元形状等都会影响有限元计算的
30、结果。故而，对有限元网格剖分的研究十分必要2.3.1 有限元网格剖分的基本原则（1） 网格数量有限元分析计算的精确度和规模的大小与网格的数目有关。网格数目变大，虽然会提高计算的精确度，但相应会扩大计算规模，因此应充分考虑这两个因素的影响来确定网格的多少。首先要确保计算的精度。对于结构形状相对简单的实体模型选择较多网格；相反，为了不降低计算的精确度而又要保证网格的数量不要过多，一般用分布式计算法或其他的方法降低模型的计算规模。（2） 网格密度网格疏密是指网格在结构的不同位置大小不相同。疏密不同的网格剖分基于以下原则：在计算分析的数据变化较大的地方，采用较密集的网格，来更好的反应数据变化的规律；相
31、反，在计算分析的数据变化较小的地方，采用相对稀疏的网格，来减小模型的规模。基于以上原则，实体结构的网格剖分密度便表现出疏密不同。（3）单元阶次有限元分析计算的精确度与单元阶次有关，单元分线性、二次和三次三种形式,后两种形式的单元称为高阶单元，高阶单元的曲线或曲面边界可以近似的逼近实体模型的曲线或曲面边界,因此对于结构形状不规则或很复杂的实体模型可采用高阶单元,以便提高计算的精确度由于高阶单元有很多结点,网格数目相同时,由其构成的模型的规模会很大,所以在选择时要合考虑计一算的精确度和规模，结点数量相同时,选择高阶次的单元会更好。（4） 网格单元形状网格单元形状的好坏对计算精度有着很大的影响,在网
32、格剖分时应保证合理的单元形状,即使只有一个单元形状很差或畸形(畸形网格通常有单元结点交叉编号、结点重合、内角大于180和单元面积为零等类型)时,也可能给计算结果带来很大的误差, 甚至使得计算无法进行下去。（5）网格布局当实体模型的结构形状对称时,为使其表现出这种对称性,生成的网格最好也保持对称，动态分析时,实体的结构形状对称,若生成的网格不保持对称,那生成的网格质量也不具有对称性，利用集中质量矩阵计算质量分布时,采用图2.12中的对称网格,结果是对称的；采用图2.13中所示的不对称网格,结果是不对称的,这将不利于动力计算。图2.12 网格对称图2.13 网格不对称（6）结点和单元编号网格结点和
33、网格单元的编号会影响结构总刚度矩阵的带宽与波前数,从而影响计算时间和存储空间的大小，合理的结点编号有利于提高计算速度和计算精度,这对于大型的复杂结构更为重要。2.3.2 有限元网格生成的方法（1） 映射法映射法的基本思想是：通过适当的映射函数将待剖分物理域映射到参数空间中形成规则参数域；对规则参数域进行网格剖分；将参数域的网格反向映射回物理空间，从而得到物理域的有限元网格。其主要特征是采用了“调配函数”概念。产生的网格整齐划一，非常规则，但同时这也是该法的缺点，尤其是当待剖分区域不规则时，如瓶颈形状，得到的网格形状是畸变的。另外，映射法是非全自动方法，必须通过人工交互方式，将剖分对象先剖分成具
34、有简单拓扑关系的子域。但映射法处理曲面问题很有效。映射法的优点是：算法简单、速度快、单元质量好、密度可控制，它既可生成结构化网格又可生成非结构化网格，既可生成四边形单元网格又可生成六面体单元网格，可用于曲面网格生成，可与形状优化算法集成等。映射法一般可直接处理单连通域问题，但对于复杂多连通域问题需要首先用手工或自动方法将待剖分域分解成几何形状规则的可映射子区域，然后在每个子区域内应用映射法。然而在实践中仍有几个难点需要克服：如何自动地将复杂的不可映射的待剖分域分解成简单的可映射的子区域；如何满足某些物理问题中对网格疏密过渡的要求；如何满足子区域之间的网格相容性要求。（2）几何分解法在产生结点的
35、同时，也确定了结点之间的连接关系的网格剖分方法称为几何分解法。这类算法基于问题域的拓扑几何描述，通过从域中逐个移去单元而生成有限元网格。它较多地考虑了待分域的几何特征，确保生成质量较好的网格单元，通常有两种方法：递归法和迭代法。几何分解法可实现从实体几何描述到初始网格生成之间的自动转换，并允许网格密度变化，但只能通过边界点的分布来控制网格规模，网格质量不高，且很难实现局部自适应加密。几何分解法的最大优势是离散时考虑了网格的形状和大小，因此，所生成的网格单元形状和分布比较好。但是，这种方法自动化程度比较低，也不利于复杂件的网格生成。（3） 基于栅格法基于栅格法的基本剖分流程如下：首先用一组不相交
36、的尺寸相同或不同的栅格覆盖在目标区域上面，保留完全或部分落在目标区域之内的栅格，删除完全落在目标区域之外的栅格；然后对与物体边界相交的栅格进行调整、剪裁、再分解等操作，使其更准确地逼近目标区域；最后对内部栅格和边界栅格（特别是后者）进行栅格级的网格剖分，进而得到整个目标区域的有限元网格。这种方法预先产生网格模板，然后将要进行网格化的物体加到其上，并在实体内部尽可能多地填充规则的长方体或正方体网格，在实体的边界上根据实体边界的具体特征更改网格的形状和相互连接关系，使边界上的单元尽可能无限地逼近物体的边界形状。这种方法能实现网格生成的自动化，网格的生成速度也非常快，能够生成的单元类型很多，剖分简单
37、，效率较高。其最大缺点是物体边界单元的质量较差；另一个缺点是所生成的单元尺寸相近，网格密度很难得到控制。栅格法首先用交互方式将物体剖分为形状简单的子区域，每个子区域分别用定形的网格模板作为规整的部分，再采取适当的措施，使得相邻子区域在结合面共享公共的结点，并且网格相容。（4） 结点连元法结点连元法一般分为两步：在物体的边界和有效区域内按照网格密度的要求均匀布点；根据一定的准则将这些结点连接成三角形或四面体网格。最早提出的结点生成法是随机布点法，该方法不能保证布点的均匀性，同时为了确定可用结点，要进行点距检查和多次试探,算法效率低。近年出现的长方形网格直接布点法和距水平扫描线直接布点法等方法在一
38、定程度上提高了结点连元法的计算效率，但Delaunay三角剖分方法(简称DT法)的出现才真正使得结点连元法成为目前主流的全自动网格剖分方法之一。DT方法仅适用于凸区域，对于非凸区域则不能保其边界的完整性,因此对非凸区域必须引入一个恢复边界的步骤。但恢复了边界完整性的三角剖分并不能严格满足Delaunay准则，因此称为约束DT。DT有两个重要特性：最大-最小角特性和空外接圆特性。DT的最大-最小角特性使它自动地避免了生成小内角的长薄单元,因此特别适用于有限元网格生成。空外接圆特性，是指在DT中的每个三角形单元或四面体单元的外接圆(二维)或外接球(三维)都不包含其它结点。（5） 网格模板法网格模板
39、法生成单元主要分两步（以三维实体为例）：第一，将待剖分网格的实体用适当大小的立方体（树根）完全包容，按照“一化八”的原则递归离散，通过对每个八分块分类，形成和八分块的并集，称为网格模板模型。第二，根据得到的网格模板模型，再进行网格剖分处理。网格模板法的优点是网格生成完全自动，网格剖分速度快，非常适用于自适应网格生成；主要缺点是边界单元形状难于完全保证。另外，网格模板法对物体的方向较敏感。（6） 推进波前法近年来，推进波前法也已经成为目前最流行的通用的全自动网格生成方法之一。其基本原理是：设区域的有向离散外边界集和边界前沿点集已经确定，按某种条件沿区域边界向区域内部扣除三角形（四面体）直到区域为
40、空集。的关键技术有两个：区域的边界离散与内部结点合理生成，扣除三角形条件。而扣除三角形条件有多种：最短距离条件、最大角条件、最大形状质量条件、最小外接圆条件等。AFT可以全自动地在平面或曲面上生成网格，用户可控制生成单元的几何分类：四边形或三角形，或者四边形和三角形混合网格。这种自动计算网格的方法是一次生成一个单元，从区域的边界向内部逐渐生成全域网格。由它生成的网格同样有着很好的几何尺寸和形状且疏密过渡平滑。当网格疏密过渡较剧烈时，它也同样能够生成高质量的网格。AFT方法的特点之一是能够在生成结点的同时生成单元，这样就可以在生成结点时对结点的位置加以控制，从而控制单元形状、尺寸以达到质量控制、
41、局部加密及网格过渡的要求。大量文献提出了各种不同的结点生成方法及单元生成方法。AFT方法的特点之二是在生成新单元时需要进行大量的相交判断、包含判断，以及为了保证单元的质量而进行的距离判断。相交判断包括线段之间的相交判断,线段与三角形面片之间的相交判断；包含判断主要指单元是否包含前沿结点的判断；距离判断包括线段与线段的距离，线段与前沿结点的距离以及线段与三角形面片的距离。2.3.3 网格质量的度量准则为了保证有限元分析计算的精确度。我们期望生成的网格的质量尽量高对。网格质量的度量，是网格生成和网格优化过程中的重要环节，网格质量的好坏将直接影响有限元分析的计算结果，因此网格剖分后生成的网格单元的质
42、量是一个好的网格生成算法必须考虑的方面。三维空间中，网格剖分的单元主要有两种：四面体单元和六面体单元，四面体网格简单灵活，对复杂的边界有较强的适应能力，故成为三维有限元网格剖分中常用的单元，网格单元质量的度量标准应满足以下原则：单元的旋转、反射、平衡、均匀绽放都不应该改变其度量值；当四面体的体积趋近于零时，其度量值也近似于零；当且仅当四面体是正面体时，其度量值取得最大值。下面介绍几种常用的度量准则：(1) 半径比 （2.1）其中，和分别为四面体内切圆和外接圆的半径，其表达式 （2.2） （2.3）其中，、是四个面的面积，、分别为对边边长之积，是四面体的体积。(2) 系数Q （2.4）其中=18
43、32.8208，与为连接顶点和的边的长度，的取值可以保证正四面体的质量度量值取最大值。fff（3） 系数 （2.5）（4） 最小立体角 （2.6）其中由下面的公式确定 （2.7）其中是由顶点1，2，3，4，组成的四面体体积，为连接顶点和的边的长度，、是由指标轮换得到。2.3.4 Rhino网格工具Rhino也和许多CAD软件一样，自身都带有网格剖分的功能，而且网格剖分都分自动的剖分和手动剖分，本文所涉及的网格剖分都是网格自动剖分，所用的模型都是实体模型。Rhino中网格剖分工具的使用是比较简单的，Rhino中网格工具的按键如图2.14所示。在Rhino中网格剖分都是自动完成的，因此我们只需要设置好网格的参数就行了。网格参数的设置可分为简易设置和高级手动设置，简易设置界面，如图2.15所示，它的作用一目了然，调节滑竿来调节多边形的网络密度。网格面较少时模型粗糙，不精确，但面少，文件和占用资源小。图2.14 Rhino网