基于Rhino的有限元网格离散技术研究_毕业论文.doc
《基于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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 Rhino 有限元 网格 离散 技术研究 毕业论文
链接地址:https://www.31doc.com/p-3921927.html