3 网格生成技术30.docx
- 文档编号:23479189
- 上传时间:2023-05-17
- 格式:DOCX
- 页数:23
- 大小:194.12KB
3 网格生成技术30.docx
《3 网格生成技术30.docx》由会员分享,可在线阅读,更多相关《3 网格生成技术30.docx(23页珍藏版)》请在冰豆网上搜索。
3网格生成技术30
3网格生成技术(通用汽车徐工为辅,约30页,周老师为主)
3.1概述
对流动与传热问题进行数值计算的第一步是生成网格,即要对空间上连续的计算区域进行剖分,把它划分成许多个子区域,并确定每个区域中的节点。
由于工程上所遇到的流动与传热问题大多发生在复杂区域内,因而不规则区域内网格的生成是计算流体力学与计算传热学中一个十分重要的研究领域。
实际上,流动与传热问题数值计算结果的最终的精度及计算过程的效率,主要取决于所生成的网格与所采用的算法。
现有的各种生成网格的方法在一定的条件下都有其优越性及弱点,各种求解流场的算法也各有其适应范围。
一个成功而高效的数值计算,只有在网格的生成及求解流场的算法这两者之间有良好的匹配时才能实现。
自从1971年Thompson等人提出生成适体坐标的方法以来,网格生成技术在计算流体力学及传热学中的作用日益被研究者所认识到。
从1986年召开第一届国际计算流体力学网格生成会议以后,该会议每隔2一3年召开次,一直延续至今,就可以看出网格主成技术在计算流体与计算传热学中的地位及这一方面研究的活跃程度。
文献中现有的生成复杂计算区域中网格的方法大致可以按图2-1所示方式来分类。
从总体上来说,流动与传热问题数值计算中采用的网格可以大致分为结构化网格与非结构化网格两大类。
一般数值计算中正交与非正交曲线坐标系中生成的网格都是结构化网格.其特点是每一节点与其邻点之间的联结关系固定不变且隐含在所生成的网格中,因而我们不必专门设置数据去确认节点与邻点之间的这种联系。
生成适体坐标的方法原则上都是些特定的变换,即把物理空间上的一些不规则区域变换成为计算空间上的规则区域。
3.2结构化网格
在结构化网格中,每一个节点及控制容积的几何信息必须加以存储,但该节点的邻点关系则是可以依据网格编号的规律而自动得出的,因而不必专门存储这一类信息,这是结构化网格的一大优点。
但是,当计算区域比较复杂时,即使应用网格生成技术也难以妥善地处理所求解的不规则区域,这时可以采用组合网格,又叫块结构化网格。
在这种方法中,把整个求解区域分为若干个小块,每一块中均采用结构化网格,块与块之间.可以是并接的,即两块之间用一条共公边联接。
也可以是部分重叠的。
这种网格生成方法既有结构化网格的优点,同时又不要求一条网格线贯穿在整个计算区域中,给处理不规则区域带来不少方便,目前应用很广。
这种网格生成中的关键是两块之问的信息传递。
结构化网格的生成源于美国LaurenceLivemore国家实验室的Winslow和Crowley在20世纪60年代后期的工作,以及前苏联科学家Godunov和Prokopov在同一时期的工作。
结构化网格是指物理空间中的三簇网格线可以映射为一个逻辑空间中三簇网格线,单元和节点的相对位置在映射后保持不变。
其特征为同一簇中的曲线不能相交,不同簇中的两曲线仅相交一次,物理平面上
的节点与计算平面上的节点一一对应。
物理空间(x,y,z)中单元或节点的标志可以用三个独立的逻辑坐标(i,j,k)表示,i=1,2,…,nx,j=1,2,…,ny,k=1,2,…nz。
其中,nx,ny,nz分别是物理空间中沿三簇网格线的节点数,计算沿i,j,k三个方向顺序进行。
当计算区域形
状复杂时,通常要将计算区域扩充到包围复杂形体之外,用阶梯形逼近复杂形体的不规则边界。
这种“区域扩充”使得实际形体边界之外还包围了大量的虚单元,因而计算效率较低。
KIVA-II就采用了这种网格。
3.3非结构化网格
非结构化网格由干对不规则区域的特别适应性而自20世纪80年代以来得到迅速的发展,在这种网格中单元与节点的编号无固定规则可遵循,因而除了每一单元及其节点的几何信息必须存储外,与该单元相邻的那些单元的编号等也必须作为联接关系的信息存储起来,使非结构化网格的存储信息量比较大。
在有限容积法中引入非结构化网格后,使有限容积法与有限元法间的差别缩小了,在某些情形下两者是等价的,因而文献称它们是一对“好朋友”。
“非结构化直角坐标法”是20世纪90年代发展起来的处理不规则计算区域的方法,
它利用直角坐标系中的网格生成的最大优点—概念上的简便,通过局部加密网格及边界上的一些特殊处理来适应不规则的区域。
将结构化与非结构化网格的各自特点结合起来,对某些及几何形状采用混合网格的做法也是近年来网格生成技木中的一种发展方向。
非结构化网格是相对于结构化网格而言的,它不存在上面结构化网格的限制。
它的物理空间和逻辑空间不再一一对应,节点和单元随意分布,控制容积可以任意选取(例如,在二维情况下,
控制容积可以取为三角形、任意四边形等,而在三维情况下,控制容积可以取为四面体或六面体等)。
由于它不必保证网格的正交性与没有邻近控制容积总数的限制,所以,其特征表现为网格划分的自由度较大,对不规则边界也有可能较好地描述。
多种非结构化网格生成方法在其生成过程中都要采用一定的准则进行优化判定,因而生成网格的质量较高,且比结构网格更易加密,生成时间也较少,这是因为一旦在边界上指定网格的分布,在两个边界之间可以自动生成网格,无须分块或用户干预。
非结构化网格对复杂形状的计算区域具有良好适应性这一突出特点,它使得在结构分析领域占据统治地位的有限元法中得到了广泛的应用。
正是由于上述优点,
近年来对非结构网格生成技术给予了高度重视,有了很大的发展。
但也由于非结沟网格的数据结构不规则,控制方程的离散化及计算程序都较为复杂,从而计算时间会大大增加,计算精
度不够高,内存占用量较大,而且在粘流分析中要求贴近物面的网格更密集,在边界层中要求使用拉伸比很大的网格单元,而通常的Delauney法和推进阵面法按其本性是难于完成此任务的,因而,非结构化网格在粘流计算及某些其他领域中仍是具有挑战性的课题。
这些都意味着非结构化网格可能在计算流体力学领域中还会有进一步的发展。
分块结构化网格是介于上述两种网格之间的一种新的网格生成技术,兼有二者的优点,其思想始见于1982年Robbert和Lee的论文;1985年Weatherill和Forsey的工作真正掀起了分块结构化网格的研究热潮,之后,国内外都对其进行了不少的研究,得到了进一步的发展,可以说开创了CFD工程应用的新局面。
这种方法按照研究对象将计算区域分成若干块,对单块区域采用结构化网格,通过贴体坐标系统用Poisson方程来生成。
这样通过拼接而集成的整体网格,具有灵活多样的特性,它既能适应复杂形体的需要,又大大地减少了虚单元的含量,提高了计算效率。
这里采用了分块结构化网格。
分块结构化网格是CFD工程实践中应用最广泛的网格类型,也是发展最为成熟的一种网格生成技术。
如KIVA-3通过采用分块结构化网格消除了保留大量冻结单元的必要性,克服了低效性,并且减少了复杂形状的存储单元和计算时间。
采用间接寻址,利用系统的结构化矢量,允许单元完全任意粘接。
每个单元用一个而不是三个下标标志,相邻单元的标志存放在数组中,数组规定了与计算单元有公共面的六个相邻单元的下标。
这样整个网格被伪单元包围,从而简化了进/出流边界条件的使用。
采用分块结构化网格,相应的流场求解算法比较成熟,执行效率高,特别适用于粘流计算;可以处理复杂的几何外形,降低了对计算机内存的要求,方便区域分裂并行计算过程的实施,适于采用并行计算的流动模拟;在不同的块中可以采用不同的求解方程,以提高效率。
当然,分块结构化网格生成过程相对较繁琐,对于外形十分复杂的实体其自动分块技术还有待于进一步研究:
分块结构化网格生成技术途径如下:
(1)分块方案设计。
确定网络拓扑,做出分块草图,进行子块编号,设置曲线和曲面上拟分布的网格节点数和网格节点分布等总体参数。
块的划分是分块结构化网格生成的第一步,也是最关键的一步。
一个好的分块方案是生成高品质网格的前提。
分块方案包括各个子块采用何种网格拓扑、子块与子块之问将如何粘接、每个子块网格分布多少网格点、网格点如何分布等内容。
分块方案的设计可以说是一种“艺术”,有较强的技巧性,需要综合考虑流场物理特性、几何特性、计算程序对网格的要求以及计算机资源等因素。
(2)单块网格生成。
单块网格生成过程一般包括曲线网络生成、曲面网格生成和体网格生成等三个步骤。
曲线网格生成就是生成各块网格的边界曲线网格。
曲面网格生成是通过曲线网格生成各网格块的边界曲面网络。
当曲面为平面时,可根据曲线网格采用合适的二维网格生成方法生成曲面网格;对特殊曲面,如平移面、旋转面等,可采用特殊方法(比如代数法)生成;当曲面为自由曲面时,曲面网格很难生成,此时可利用计算机进行几何建模(详见2.4),将几何模型数据导出,并读入网格生成系统,获取表面网格点坐标。
或通过三维坐标测量仪获取表面网格点坐标,最后进行网格点重排得到所需求的网格曲面。
体网格生成是根据每个子块已求得的六个表而网格,通过插值生成初始网格,由椭圆方程控制网格点分布和网格的光滑,通过Sorenson椭圆方程控制网格的边界正交,以此生成三维网格。
(3)块合并。
将各子块合并生成整体网格(详见4-2),进而进行网格质量分析。
3.4常用网格生成软件
计算机辅助网格生成应在几何建模平台和网格生成软件支撑环境中完成。
目前国内外已开发出许多商业几何建模平台和网格生成器,几何建模平台如Pro/E,UG,IDEAS,CATIA.等(详见2.2),网格生成器有EAGLEVIEW、GridGen、Geme++、GRAPEVINE、CAEGrid、ICEM/CFD等。
ICEMCFD技术
ICEMCFD作为CFD的前处理系统,采用独特的网格生成方法,可以根据初始三维几何形体控制模型直接生成网格,无需另外建立单独的分析模型,且几何形体控制模型可以用任何商业CAD软件创建的三维CAD数据、STL(stereolithography)数据、或点数据表示,其CAD界面允许用户在自带的CAD系统中操作,用户能够选择面、曲线和点,并将这些元素组合成簇,为网格生成进行网格大小分配。
选取的集合形体可以直接转换为ICEMCFD的几何形体数据库,在这个数据库中用户可以进入网格模块。
这些界面可用于I一DEAS,Pro/E,UnigraPhics和CATIA,该界面提供的环境将CAD系统中的参数化几何形体创建工具与ICEMCFD软件中的计算网格生成、后处理和网格优化工具相连接。
ICEMCFD通过标准的数据格式(IGES、VDA/FS、DF)由自带的CAD软件导入CAD数据,或在ICEMCFDDDN中建立和修改几何形体控制模型。
ICEMCFD是一个功能强大的网格生成器,它提供了精确获取几何形状、自动进行网格划分、后处理及网格优化等功能,利用其中高度自动化、独特的网格生成工具箱可以生成高质量的结构化网格,非结构化网格和混合网格,具有多种网格输出接口,如Nastran、Ansys、FluentKiva-3(v))等,输出网格可以直接应用于有限元计算和CFD分析,其所花时间只是其他网格生成器的一部分,并且在设计过程中网格特性可以完全进行控制。
不同类型的计算网格分别通过不同的网格生成模块生成,这些网格生成模块包括Hexa、Tetra、Prism、Quad、AutoHexa、Global、P-CubeMuCad等。
其中Hexa模块是一个基于对象的半自动三维曲面和体网格生成器,可以生成多块结构化六面体网格和多块非结构化六面体网格。
ICEMCFD集成了一组输出过滤器,利用这些输出过滤器可以自动生成约100种CFD或CAE求解器的输入文件。
这样,IcEMCFD在几何形体与计算分析之间提供了一个直接的接口。
ICEMCFD网格生成器与CAD等有关求解器之间的关系如图所示。
在用Hexa生成网格过程中涉及到两种类型的实体:
块拓扑(Blocktopology)和几何形体(Geometry)。
块拓扑模型直接在CAD几何形体上生成。
在交互式地创建了等价几何形体的块拓扑
模型后,块拓扑可以通过边(edge)、面(face)和块(block)的分裂进一步划分块。
块的顶点可以单个或成组移动到相关的曲线或CAD曲面上,一些特定的块、边、线可以与重要的CAD曲线相
关联以捕捉网格中重要的几何特征。
Heax在生成网格后,首先生成块拓扑中各个块的网络,然后将块的六个面网格投射到最新的CAD几何形体表面。
如果投射面不是期望的表面,则可采用与投射面关联的方式,将特定的块拓扑面与指定的儿何形体表面关联起来,生成正确的期望网格。
用ICEMCFDHeax生成网格时一般包括以下步骤:
(1)族的创建。
几何形体文件导人,生成tin文件;对几何形体进行必要的修改,完成族(family)的创建。
(2)块拓扑模型的创建。
修改几何形体并创建族之后,再通过块分裂、合并、边/面修改及顶点移动交互地创建块拓扑模型,不管几何形体如何复杂,Heax最初都是生成一个规则的六面体块(长方体块),再根据几何形体形状,通过块的分裂和合并划分子块。
(3)网格生成。
边界线网格划分参数的分配,不同投影方式网格的生成。
(4)网格的质量检查与改善。
网格质量与具体问题的几何特性、流动特性以及流场求解算法有关,最终要由计算结果来评判。
误差分析及经验表明,CFD对计算网格有一些基本要求,如光滑性、正交性、网格单元的正则性以及在流动变化剧烈区域网格分布的疏密性等。
对于复杂几何外形所生成的网格,这些要求往往不能同时完全满足。
例如,给定边界网格点分布,采用Laplace方程生成的网格是最光滑的,但是最光滑的网格不一定满足物面边界正交性条件,其网格点分布也很有可能不能捕捉流动特征,因此,最光滑的网格不一定是最好的网格;常用网格质量参数包括扭角(skewangle)、纵横比(aspectratio)、Laglacian以及弧长(arclength)等。
考虑到网格光滑性对于截断误差的影响以及椭圆算子的光滑性和网格边界正交性对有限差分法和有限体积法的密切相关性,一般认为网格点分布足以反映流动特征的网格质量较好,如椭圆方程生成的边界正交网格。
(5)确定边界条件,网格文件输出。
这里,可采用IGE格式实现Pro/E与ICEMCFD两个不同平台间的数据传递,即由Pro/E导出三维模型的IGES文件,再由ICEMCFD的IGES后处理器将IGES文件转换成tin文件,因为tin文件包含了模型的曲线、曲面信息,由此可创建必要的曲线和曲面,删除多余的曲面。
对修改后的几何形体的曲线、曲面和由这些曲面围成的封闭空间可创建族(family)。
这里的族是划分网格时指定的关联曲面和边界条件。
10.1ICEMCFD介绍
ICEMCFD是ANSYS公司的产品,主要应用于计算流体力学、结构分析等数值模拟领域,为相应的工程实际问题提供了方便的前处理、求解与后处理器。
ICEMCFD的前处理模块是极为优秀的专业网格划分程序;求解器模块是解决某些专业CFD问题的计算机数值求解程序;后处理模块是专业的CFD/FEA求解结果的可视化程序。
为集中解决问题,本文主要针对ICEMCFD的前处理器进行详细介绍,其他模块将在以后的学习和应用中逐步介绍。
ICEMCFD前处理器主要包括CAD几何建模处理、网格生成处理、网格优化处理以及网格输出处理4大模块,每一个模块_又根据不同需要分成独立的几个模块。
这些模块之间结合紧密、使用方便,并配有大量的教程可供参考,因此ICEMCFD前处理器具有系统性强、建模方便、界面友好、模块众多、网格划分思路清晰、运算速度快、接口众多、学习方便等其他网格划分软件无法比拟的优点。
很多人认为ICEMCFD难学,主要原因在于没有掌握基本的几何与拓扑知识和相关的ICEMCFD处理思路,本文将在前面几节介绍这些基础知识以便十初学者快速入门。
ICEMCFD的CAD几何处理模块为用户提供了方便的CAD数据处理方案,其中包括标准CAD数据转换器、直接CAD界面以及非结构化CAD数据转换器。
ICEMCFD的网格划分模块为用户提供了多种网格划分方案,其中包括各种类型的二维、三维网格划分模块,每一模块都具有十分强大的功能。
ICEMCFD的网格优化模块可以对网格生成模块创建的各种类型网格进行质量检查与优化等操作。
ICEMCFD的网格输出模块支持一百多种流行的CFD/CAE求解器,还提供了包括CGNS格式在内的标准格式,这也是其他任何网格生成软件都无法比拟的优点。
10.2ICEMCFD的处理思路
从整体上看,ICEMCFD主要用于数值模拟的前处理环节,求解与后处理模块的应用不如前处理模块广泛。
为实现网格划分,ICEMCFD首先要接收CAD软件创建的模型,因此ICEMCFD的处理思路是和CAD软件直接相关的。
导入CAD几何模型之后,就可以使用ICEMCFD进行前处理的最重要一步了—网格划分与修改。
网格划分与修改之后,再通过转换器将网格发送到相应的CAE/CFD求解器进行问题求解。
得到求解结果之后就可以通过后处理程序进行显示和报告。
这个流程构成了ICEMCFD的处理思路,由此思路展开就构成了本章的主要内容。
如图10-1所示为ICEMCFD处理思路的框图,一般说来,使用ICEMCFD进行处理都遵循这个思路进行。
由图可见,本文主要将ICEMCFD及相关模块分成5大部分:
CAD处理模块
ICEMCFD主模块、网格处理模块、网格优化模块以及求解模块。
其中,ICEMCFD主模块还可以细分为网格生成模块、ICEM到求解器的转换模块以及后处理模块等。
在具体操作上可以将ICEMCFD的处理思路总结如下:
(1)采用相应的几何建模工具(CAD软件)进行建模,并在建模过程中为相应的曲面、曲线和点命名。
(2)将所需要的网格尺寸信息与相应的几何元素的名字联系起来,这样可以方便地进行参数化网格划分。
只要几何元素命名没有改变,相应尺寸的任何改变都可以直接反映在网格的变化中。
(3)采用相应的网格生成模块来生成所需要的各种计算网格,同时实现网格的可视化。
(4)检查网格的质量并对网格进行修改与优化处理。
(5)将边界条件信息与相应的几何图形的名字关联起来,实现各种类型求解器的边界条
件定义。
(6)通过输出转换器将生成的网格文件转换为各种求解器可以识别的文件类型。
并采用流动求解器进行后处理。
(7)采用ICEMCFD后处理模块对流动求解器输出的数据进行处理,或者读入到网格优化工具实现自适应网格。
3.5GAMBIT的使用
3.5.1GAMBIT的图形用户界面
3.5.2GAMBIT绘制几何图形
3.5.3GAMBIT绘制网格
3.6FLUENT的使用
3.6.1FLUENT的安装与运行
FLUENT的安装
(1)首先在C盘上创建目录C:
\flexlm,然后将软件公司提供的license.dat文件拷贝到
这个目录下。
(2)将FLUENT的安装光盘放入光盘驱动器,然后双击setup.exe文件启动安装程序,
设置好安装目录、工作组名称等选项后,安装程序开始向硬盘拷贝文件。
(3)文件拷贝结束后,在弹出的工作组窗口中双击运行setenvironment,设置好环境
变量,再将license(软件使用许可证)服务器的地址设置完毕,重新启动计算机,就可以启动FLUENT了。
这里需要说明的是:
FLUENT与我们以往使用的软件不同之处是其安装过程中不需要使用安装密码,而是在每次启动软件之前,通过互联网在license服务器上进行软件合法性验证,只有通过验证才能使用。
启动FLUENT
以WINDOWSNT为内核的操作系统包括WINDOWS2000和WINDOWSXP,其启动
方式有两种:
(1)从WINDOWS的开始菜单中进行启动,即顺序点击:
开始一>程序一>FluentInc一>FLUENT6.1
就可以启动FLUENT。
(2)从DOS终端窗口启动,即在命令行中:
1)键入”`fluent2d"启动二维单精度计算。
2)键入”fluent3d"}启动三维单精度计算。
3)键入“`fluent2ddp"启动二维双精度计算。
4)键入“`fluent3ddp"启动三维双精度计算。
如果想启动并行计算模式,可以在上述4个命令后而加上一tx参数,其中x是并行计算
的CPU数量,例如键入“fluent3d-t3”意思是在三个处理器上运行三维计算。
3.6.2FLUENT的用户界面
FLUENT同时采用了图形用户界而和文字用户界两种界而形式进行操作,下面分别进行介绍。
1.3.1图形用户界面
FLUENT采用WINDOWS风格的图形用户界而,其界而类型主要分为下列四种:
(1)主控窗口。
(2)控制参数而板。
(3)消息窗口。
(4)图形显示窗口。
下面分别进行介绍。
主控窗口
启动FLUENT后首先进入的窗口就是主控窗口,见图1-20主控窗口分为菜单栏和文字信息窗口两个部分。
用户与FLUENT软件之间的互动可以用两种方式进行,一种是在下拉菜单中直接选择相应的菜单命令进行操作,或者在弹出的参数控制而板上进行参数选择操作,同时系统通过消息窗口将系统信息直接反映给用户;另一种是在文字信息窗口中键入命令进行操作,系统信急也以文字信息的方式反映在文字信息窗口中,就象在DOS环境中一样。
与微软公司的WORD等软件一样,当窗口中信息量超过一个页面的容量时,页面内容将自动上移,同时窗口右端的滚动条被激活,用户可以通过移动滚动条上的滑块查看滚动出页面的信息。
在命令执行过程中,可以按组合键Ctrl+C终止命令的执行。
菜单命令的执行除了可以用鼠标左键点选执行的方式外,还可以通过热键方式进行。
比如点击打开File菜单的动作,既可以通过将光标移动到File上,再按左键点击打开的方式完成,也可以通过热键方式完成。
所谓热键方式就是用键盘操作代替鼠标操作的方式,
比如上述打开File菜中一的动作可以用组合键Alt+F来完成,而不必移动鼠标。
请注意在可
以使用热键打开的菜单项中带下划线的字母就是打开这个菜单的“热键”,比如File菜单中字母下而带有下划线,即表示可以通过按组合键Alt+F的方式打开。
同样,在Grid菜单
项中,字母G下面带有下划线,表示可以通过按Alt+F组合键的方式打开,其他菜中一与此
相同。
在文字信息窗口中可以对显示的信息、命令进行剪切、拷贝、粘贴、删除等编辑操作,
方法与微软的OFFICE系列软件中‘样,即:
(1)首先将光标移动到要进行编辑区域的角。
(2)按下鼠标左键并拖动到编辑区域的另角,此时编辑区域被反白显示。
控制参数面板
控制参数面板类似于一个表格,是FLUENT软件中用于输入控制参数的主要形式之一
图1-3是用于网格显示操作的一个控制参数面板。
这个面板中包含了一些常见的可视化操作组件,如选择按钮、单选按钮、多选单和文本输入框,以及可以打开其他面板的按钮,如图中的Display(显示)按钮和Colors…(
色)按钮,以及全选按钮和不选按钮。
选择按钮和单选按钮被点击选中后,表示在计算中
将采用相应的参数设置,再次点击则可以取消该项设置。
多选单用于可以进行多重选择的场合。
具体参数的设置需要人工决定时,则可以使用文本输入框输入相应的参数。
点击全选按钮可以将多选框中所有选项选中,点击不选按钮则可以取消所有选择。
除了图1-3中所示的组件外,FLUENT中还使用下拉菜单(图1-4a)和单选菜单(图
1-4b)进行单一选择操作。
FLUENT除了使用直接的文字输入进行参数设置,还使用整数
增减框(图1-5a)和滑尺(图1-5b)进行数字的输入操作。
如同其他WINDOWS风格的软
件一样,FLUENT中也使用“标签”组件,以便在同一个面板中选择不同种类或群组的输
入操作(图1-6)。
整数增减框通过其右
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 网格生成技术30 网格 生成 技术 30