自动体网格生成共38页文档
ANSA自动生成四面体网格
wwwห้องสมุดไป่ตู้
Innovating through simulation
Logo
切换到HIDDEN 显示模式进行质量检查。
违反Skewness 准则的单元根 据右侧图标的颜色显示绿色。视图 中违反准则的单元数量显示在窗口 左侧的OFF 后面。
Innovating through simulation
Logo
接着弹出PROPERTIES 窗口，选择 一个Property赋给将要创建的体网格。因为 在列表中没有Property，点击NEW 创建一 个新的Property。弹出新的Solid Property 卡片。（ 注意到卡片的格式取决于当前的 DECK，在本例中是NASTRAN）。
Logo 4.1 ANSA中的四面体网格生成的步骤和方法
ANSA中的四面体网格生成的一般步骤：
读入零件的ANSA 文件。 为生成四面体网格，首先创建三角形面网格。 检查并提高面网格质量。 定义将要生成四面体网格的封闭区域。 生成四面体网格。 获取关于四面体网格和零件信息。
Innovating through simulation
Logo
在生成四面体网格之前首先定义VOLUME。
点击VOLUMEs>DEFINE功 能，在模式窗口中选择AUTO 选 项。
左键点击任何面以选中所有 的面形成封闭体。点击鼠标中键。
Innovating through simulation
Innovating through simulation
Logo
Innovating through simulation
体网格生成示例
体网格生成示例以生成f6外流场体网格为例讲述生成体网格的过程，整个过程包括几何导入、网格尺寸设置、表面网格生成、构建外包围盒、生成体网格、查看内部网格和边界条件设置。
1 导入几何在菜单栏中，选择“文件”> “添加文件…”> “添加几何文件…”，弹出选择几何文件对话框。
可以添加的几何文件类型有IGES、STEP、SA T、FLI和GM2(3)格式文件，其中前三种为标准数据交换格式，后两种为本软件自有几何文件格式。
在本示例中，我们导入IGES格式几何文件，文件导入后会弹出图1所示对话框用于设置全局容差。
图1 容差设置对话框如图2所示为导入到软件中的f6外形，图3所示为几何外形的三视图。
几何导入后，我们可以查看几何信息。
在菜单栏中，通过“窗口”> “拓扑信息”，可以查看几何外形包含的曲线数目、曲面数目和区域数目；通过“窗口”> “尺寸信息”，可以查看几何外形尺寸值。
通过工具栏上的“选择”下拉列表，我们可以使用鼠标左键点击选取几何元素，选中某一几何元素后，会在信息输出窗口显示选取的几何元素信息。
图2 f6外形在功能操作面板上，可以通过“快速设置视点”选择多种视图查看几何；也可以通过对象视图控制，设置可视的几何元素。
另外，点击“高级控制”按钮，在弹出“可视化控制”对话框中，我们可以控制特定类型和特定编号的几何元素的显示和隐藏。
在视图模式下，我们可以通过鼠标控制，实现对几何外形的旋转、放缩和平移等操作。
鼠标左键、中键和右键的功能分别为平移、旋转和放缩。
图3 几何外形三视图2 网格尺寸设置我们可以通过指定曲线离散段数、设置背景网格和网格源等方式实现对网格尺寸的控制。
三维项目中背景网格需要设置的唯一参数是“全局尺寸”，在“功能操作面板”的“网格”选项卡中的“密度控制”中进行设置，通过设置背景网格得到的网格尺寸是均匀的。
当需要控制网格局部尺寸时，我们可以采用网格源的方式，网格源包括点源、线源和面源。
网格生成技术
同时在物理平面的求解区域边界上规定ξ（x，y），η（x，y）的取值方法，于是就形成了物理平面上的第一类边界条件的Laplace问题.
②从计算平面上来看，如果从计算平面上的边值问题出发考虑，则情况就大为改观，因为在计算平面上可以永远取成一个规则区域。所谓计算平面上的边值问题，就是指在计算平面的矩形边界上规定x（ξ，η)，y（ξ,η)的取值方法，然后通过求解微分方程来确定计算区域内部各点的（x，y）值，即找出与计算平面求解区域内各点相应的物理平面上的坐标。实际上用椭圆型方程来生成网格时都是通过求解计算平面上的边值问题来进行的.为此需要把物理平面上的Laplace方程转换到计算平面上以ξ，η为自变量的方程.
物理平面求解区域内部的网格疏密程度要易于控制；
在贴体坐标的边界上，网格线最好与边界线正交或接近正交,以便于边界条件的离散化。
生成贴体坐标的过程可以看成是一种变换,即把物理平面上的不规则区域变换成计算平面上的规则区域，主要方法有微分方程法，代数生成法，保角变代以来发展起来的一种方法，基本思想是定义计算域坐标与物理域坐标之间的一组偏微分方程，通过求解这组方程将计算域的网格转化到物理域.其优点是通用性好，能处理任意复杂的几何形状，且生成的网格光滑均匀，还可以调整网格疏密，对不规则边界有良好的适应性,在边界附近可以保持网格的正交性而在区域内部整个网格都比较光顺；缺点是计算工作量大。该方法是目前应用最广的一种结构化网格的生成方法，主要有椭圆型方程法、双曲型方程法和抛物型方程法。
单元(Cell）是构成网格的基本元素.在结构网格中，常用的2D网格单元（图1。1)是四边形单元,3D网格单元(图1.2）是六面体单元。而在非结构网格中,常用的2D网格单元还有三角形单元，3D网格单元还有四面体单元和五面体单元，其中五面体单元还可分为棱锥型（楔形）和金字塔形单元等。
ICEM网格生成流程
ICEM网格生成流程预览说明:预览图片所展示的格式为文档的源格式展示,下载源文件没有水印,内容可编辑和复制Chapter 3二维非结构壳/面网格生成（2、3）1. 创建几何模型：Point --- Curve --- Surface --- Part --- Topology 2．定义网格参数2.1．定义全局网格参数2.1.1 定义网格全局尺寸：Scale factor、Max element2.1.2 定义全局壳网格参数：Mesh type、Mesh method2.2 定义Part网格尺寸3. 生成网格并导出3.1 生成网格，检查网格质量3.2 保存网格文件：Save mesh as…3.3 选择求解器：Output --- Select solver3.4 写入：Output --- Write inputChapter 4三维非结构自动体网格生成(自上而下)（2、3）1. 创建几何模型：Point --- Curve --- Surface --- Part --- Topology --- Body2．定义网格参数2.1．定义全局网格参数2.1.1 定义全局网格尺寸：Scale factor、Max element2.1.2 定义体网格全局参数：Mesh type、Mesh method2.1.3 定义棱柱网格全局参数：Grow Law、Initial height、Ratio、No.2.2 定义Part网格尺寸3. 生成网格并导出3.1 生成网格，检查网格质量3.2 保存网格文件：Save me sh as…*.uns3.3 选择求解器：Output --- Select solver3.4 写入：Output --- Write input三维非结构自动体网格生成(自下而上)（4）首先导入壳网格，在壳网格的基础上拉伸生成棱柱体网格，再填充棱柱体网格和远场边界之间的空隙。
（壳网格---棱柱体网格---体网格）。
Ansys 第七节 ICEM实例-棱柱体网格自动生成
9/9/05
ANSYS ICEMCFD V10
Inventory #002277 C4-13
平滑光顺 四面体/棱柱体 网格
平滑光顺 四面体/棱柱体 网格
– 首先只平滑 tets 和 tris（四面体/三角形） • 设置 PENTA_6 为 Freeze • 此时的要点是不要试图更改棱柱层结构
– 一旦tetra和tri单元足够的平滑, 则平滑所有单元 – 设置 PENTA_6 为 Smooth – 减低 Up to quality 值，以免过度扭曲棱柱体网格单元
9/9/05
ANSYS ICEMCFD V10
Inventory #002277 C4-14
切割棱柱体网格层
切割棱柱体网格层
– 如果需要许多棱柱体网格层, 这种方法更具鲁 棒性 – 通常更快 – 建立 “fat” layers，然后在编 辑网格时切割它们。
– Edit Mesh > Split Mesh > Split Prisms
ANSYS ICEMCFD V10
Inventory #002277 C4-10
高级选项 – Max Height Over Base
Max height over base – 限制棱柱体网格的纵横比 – 在棱柱体网格的纵横比超过指定值的区域棱柱层停止生长
– 棱柱层的数目在局部区域无法保证 – 在棱柱层边界网格融接为金字塔形网格
高级选项-Part Control
New volume part
– 指定新的part存放棱柱单元或者从已有的面或 体网格part中选择
Side part
– 存放侧面网格的part
Top part
– 存放最后一层棱柱顶部三角形面单元
自动网格生成法
自动网格生成法二维网格生成—Advancing Front方法从概念上来讲，Advancing front方法是最简洁的方法之一。
单位元素生成算法始于一个特殊边界条件所定义的“front”,此算法逐级地生成各个元素，同时“front”元素离散地前进，直至整个区域都被元素所覆盖。
网格生成过程包括三个主要步骤：1、在边界上生成节点，形成一个离散的区域边界。
2、在离散区域边界内生成元素（亦或节点）。
3、强化节点形状以提高网格图形清晰度。
在介绍这个方法之前我们先介绍以下有关于二维空间地几何表示。
一、二维网格的几何特征我们利用网格参数（一般是空间的函数）来表征网格的一些性质，诸如节点尺寸，节点形状和节点方向等等。
网格参数包括两个相互正交的单位矢量a1和a2表示的方向参数，和由两个相互正交代表节点形状的矢量的模值h1和h2。
前者表征网格节点伸展的方向，注意的是，只有在生成的是非各向同性的网格内，方向参数才有定义，否则方向矢量是常单位矢量，而尺寸参数有h1=h2，这样就定义了各向同性的平凡网格。
二、区域的几何表示边界曲线的表示：我们一般用组合参数样条线表示曲线边界单位，利用参数t，我们利用二维矢量函数表达出曲线边界：r t=x t,y t,0≤t≤1一般来讲，一条组合样条曲线至少是C1连续的，以保证边界曲线平滑和算法要求的数学连续性。
我们下面将要用厄米三阶样条线，当然还有许多就不一一举例了。
样条线的参数表达式如下：X t=H0t,H1t,G0t,G1t∗x0,x1,x,t0,x,t1T,0≤t≤1转置的前两项是曲线的两个端点，而后两项是它们对t求导现在端点处的值。
另外G和H分别是四个三阶厄米多项式：H0t=1−3t2+2t3 ; H1t=3t2−2t3G0t=t−2t2+t3 ; G1t=−t2+t3此时，参数表达式可以通过一个系数矩阵来描述：X t=1,t,t2,t3M x0,x1,x,t0,x,t1T,0≤t≤1其中M矩阵读者很容易写出，是一个4*4的方阵，而每一列是这些厄米多项式的系数排列而成。
第二章 网格生成
第二章网格生成本章要点●三种基本网格生成方法●程序的网格加工功能及使用方法●具体网格生成实例本章主要介绍网格生成的基本功能及使用方法。
先介绍MENTAT生成网格的三种基本方法：节点与单元的直接定义，或先定义结构的几何实体，再转换为单元和节点，以及对任意几何面、体自动生成单元的网格划分。
然后介绍MENTAT对已有网格进行加工、处理的功能及使用方法，最后通过具体实例介绍使用户更好地掌握MENTAT网格生成部分的菜单、命令使用方法。
与有限元分析相关的常用词ELEMENT （单元）由多个节点定义的用于分析的最基本区域。
NODE （节点）用于定义单元的点，具体位置由坐标值确定。
与几何实体相关的常用词POINT （点）描述曲线、曲面的控制点。
CURVE （曲线）线段、圆弧、样条等曲线的统称。
SURFACE （面）四边形面、球面、圆柱面等曲面的统称。
网格生成方法1（网格直接定义）网格生成主要是生成节点和单元，节点的位置由节点坐标决定，节点坐标可以由键盘输入，也可在屏幕上检取格栅点输入。
键盘输入方法简单，但使用起来却不方便，检取格栅点的方法更为人所常用，下面介绍一下格栅的使用方法。
格栅的显示格栅的显示分二步，首先必须检取格栅显示光钮，而且设置合适的格栅参数如间隔、大小等。
格栅的定义光钮在MESH GENERATION菜单的中部。
用户在MAIN菜单中检取MESH GENERATION，在MESH GENERATION菜单的中部的COORDINATE SYSTEM条目下有绿色的SET、RECTANGULAR、GRID光钮，如下图所示。
MESH GENERATIONNODES ADD REM EDIT SHOWELEMS ADD REM EDIT SHOWPTS ADD REM EDIT SHOWCRVS ADD REM EDIT SHOWSRFS ADD REM EDIT SHOWSOLIDS ADD REM SHOWBETWEEN NODE BETWEEN POINTSELEMENT CLASS QUAD(4)CURVE TYPE LINESURFACE TYPE QUADSOLID TYPE BLOCKCOORDINATE SYSTEMSET RECTANGULAR ∇GRID ○CLEAR MESH CLEAR GEOMATTACH AUTOMESHCHANGE CLASS CHECKCONVERT DUPLICATEEXPAND MOVERELAX RENUMBERREVOLVE SHELL EXPANDSOLIDS STRETCHSUBDIVIDE SWEEPSYMMETRY检取GRID，使之变为红色，就变为了光栅显示状态，在图形区将显示出一个田字型的格栅，隐含格栅大小为±1，点之间间隔为0.1，如下图所示。
网格生成技术概述
§9网格生成技术概述所谓网格划分就是把空间上连续的计算区域划分成许多子区域，并确定每个子区域中的节点。
网格划分的实质就是用一组有限个离散的点来代替原来连续的空间。
网格生成技术是计算传热学（NHT）和计算流体力学（CFD）的重要组成部分，在目前的CFD&NHT工作周期中，网格生成所需人力时间约占一个计算任务全部人力时间的60%左右，网格质量的好坏直接影响数值结果的精度，甚至影响数值计算的成败。
可见网格生成技术是CFD&NHT作为工程应用的有效工具需要解决的关键技术之一。
最初，因为主要从事理论研究，求解的方程通常是比较简单的模型方程。
对于二维问题，常在比较规则的区域内研究问题，此时针对具体的问题可用较简单的代数方法生成网格，并做简单的自适应，网格问题并不突出。
但是对于有实际应用价值背景的问题，如航空航天飞行中的高超声速流动、跨音速流动以及其它多介质、高温高压系统的计算流体力学问题。
这些问题所涉及的流场十分复杂，会出现各种形式的间断，必须采用非常密的网格才能对间断有较高的分辨，从而达到需要的计算精度。
事实上，计算流体力学的发展除了依赖于计算机和数值计算方法的发展以外，还在很大程度上依赖于网格技术的发展。
因此，近几十年来网格生成技术己受到越来越多的计算数学家、计算流体力学家的重视，并己经成为计算流体力学发展的一个重要分支。
1. 网格单元的分类单元（cell）是构成网格的基本元素。
在结构网格中，常用的2D网格单元是四边形单元，3D网格单元是六面体单元。
而在非结构网格中，常用的2D网格单元还有三角形单元，3D网格单元还有四面体单元和五面体单元，其中五面体单元还可分为棱锥形（或楔形）和金字塔形单元等。
图1和图2分别示出了常用的2D和3D网格单元。
图1 常用的二维网格单元图2 常用的三维网格单元2. 网格生成方法分类网格生成方法的分类表示于图3中。
（1）结构化网格自20世纪80年代开始，各国计算流体和工业界都十分重视网格生成技术的研究，首先发展了结构网格方法。
网格生成技术
双曲方程生成的内流网格
抛物方程生成网格
5 非结构网格
• 没有固定结构，节点编号命名无一定规则甚至随 意，节点的邻点个数也不是固定不变 • 二维和三维空间中最简单的形状是三角形和四面 体，任何平面或空间区域都可被三角形或四面体 填满。 • 除了每一单元及其节点的几何信息必须存储外， 相邻单元的编号也必须存储起来
椭圆型方程法 – 封闭边界情况
双曲型方程法 – 外部流动、内部流动 抛物型方程法 – 同上
A 保角变换法
• 根据复变函数中的保角变换理论，映射得到物理 域边界和计算域边界间的对应关系，进而利用边 界的对应关系生成内部节点。 • 可以保证物理平面上所生成的网格的正交性 • 仅适用于二维问题
B 代数法
网格自适应化的方法
• 网格细化法 (h型方法)：通过网格的进一步 细化来实现自适应目标 • 重新分布法 （r型方法）：指保持单元或节 点数不变而通过重新分布节点位置实现自 适应目标
通过边界点控制网格加密
网格局部加密
外掠后台阶流场计算的 网格自适应调整
控制容积凝聚法生成的粗密网格
8 区域扩充法
第9章 网格生成技术
为什么要研究网格生成技术？
工程上的流动与传热问题大多发生在复杂区域内
网格生成：计算流体和传热中十分重要的研究领域
数值计算的最终精度及效率，取决于：
生成的网格 采用的算法 高效的数值计算： 网格生成，求解算法 良好匹配
9.1 网格生成技术概述
正交曲线坐标系（14种）
保角变换法 边界规范化方法 双边界法
• 直角坐标系简单方便，不少研究者愿意在 直用阶梯形边界逼近真实边界
采用阶梯形边界逼近真实边界
• 在计算传热学发展的早期，曾广泛采用这种方法。 由于这种网格的构造简单，可以适用于任何形状 的物体，因而近年来又引起了许多研究者的兴趣， 特别在计算大规模问题时（如环境工程问题）经 常采用。并用局部加密的方法更好的逼近曲线边 界
多源多目标扫掠体的全六面体网格自动生成算法
多源多目标扫掠体的全六面体网格自动生成算法一、导论1.1 研究背景1.2 研究意义1.3 研究现状1.4 研究内容1.5 研究方法二、多源多目标扫掠体的建模方法2.1 扫掠体的形成2.1.1 曲线生成2.1.2 添加时间参数2.1.3 生成截面2.2 扫掠体的六面体网格化2.2.1 六面体网格生成2.2.2 自适应六面体剖分三、多源多目标扫描路径规划方法3.1 扫描路径规划基本原理3.2 多源扫描路径规划方法3.2.1 分支界定法3.2.2 遗传算法3.3 多目标扫描路径规划方法3.3.1 Pareto优化算法3.3.2 支配排序算法四、多源多目标扫掠体的自动六面体网格生成算法4.1 六面体网格生成流程4.2 自动六面体网格生成算法的实现4.2.1 六面体网格的构造4.2.2 六面体网格的优化五、多源多目标扫描体六面体网格自动生成算法的实现5.1 实验设置5.2 实验结果5.2.1 六面体网格自动生成时间5.2.2 六面体网格质量5.2.3 扫描路径规划效果六、总结与展望6.1 研究成果总结6.2 研究不足与展望6.3 研究的应用前景一、导论1.1 研究背景如今全六面体网格自动生成技术已被广泛应用于航空航天、汽车、电子、生物医学等工业领域，是建立虚拟样机的基础，有着广阔的市场前景。
而多源多目标扫描体的全六面体网格自动生成则是全六面体网格自动生成技术的重要扩展和拓展方向。
多源多目标扫描体是指由多个扫描源扫描得到的具有多个目标的三维物体，是典型的多目标优化问题。
全六面体网格自动生成技术的目标，则是要将三角网格模型转化为六面体网格模型，并兼顾六面体网格质量和自适应性能。
因此，对于多源多目标扫描体的全六面体网格自动生成技术的研究意义重大。
1.2 研究意义多源多目标扫描体的全六面体网格自动生成技术的研究，能够提高六面体网格自动生成技术的适用范围和实际应用水平，满足实际工程需求。
与此同时，该技术也可以为扫描源、物体形变和加工状况等提供更精细、更全面的分析与预测。
网格生成
对于连续介质系统，例如飞行器周围的气体，集中在障碍物上的压力，或者是回路中的电磁场，又或者是化学反应器中的浓缩物，都是可以用偏微分方程来进行描述的。
为了对这些系统进行模拟。
这些连续性方程需要基于一定数量的时间空间意义的点进行离散化，并且在这些点上对各种物理量进行计算。
通常的离散的办法又下列三种1 有限差分法2.有限体积法3.有限元法都是使用相邻的点来计算我们所需要的点。
所以在计算中网格点的概念被引出。
一般来说，有二种网格的类型，是通过连接点的方式的不同而被区分的。
结构网格是正交的处理点的连线，也就意味着每个点都具有相同数目的邻点，对于非结构网格来说，也就是不规则的连接，每个点周围的点的数目是不同的。
图一给出了每种网格的例子。
在一些情况中，也有着部分网格是结构的，部分网格是非结构的（例如在粘性流体中，边界一般都是使用结构网格，其他的部分是使用非结构网格。
在这章剩下的内容中，我们将针对不同的需求讲述不同的离散的办法。
第二章主要讲在简单域中使用普通的结构网格（代数和椭圆的方法）和在复杂区域中使用混合结构。
第三章主要讲产生非结构网格的方法，主要的针对产生三角形、四边形网格（奇怪四边形网格不是结构化网格吗？）的二种办法：advancing front 和delaunay triangulation。
第四章对适应性网格作了一个简单的介绍。
1.1离散和网格种类主要的离散的办法有，有限差分法，有限体积法（和有限差分法等同）和有限元法，为了说明这些方法，我们首先来考虑连续性方程。
其中ρ是密度，U是速度，S是源项，ρU表示各个方向上的质量流，有限差分公式模拟，是用下面的办法来达到对所需要的点的模拟的。
例如，对正交网格，矩形在横轴的长度是h.尽管不规则的网格可以运用有限差分法,但是规则的网格却是最适合使用它的,这样就产生了简单的差分格式和有效的算法(例如在向量机上(什么东东)),也可以在特殊的坐标系中对正交网格使用有限差分法(例如在球形极坐标中)在有限体积法中,物理空间被分成很多小的体积V,然后对每一个小的体积运用偏微分方程来积分.然后我们会用每个小体积中的每个所求量的平均值来代替我们要求的值,,用相邻体积中的变量的函数来表示流过每个体积的表面的流量.有限体积法来进行离散可以用于结构或者非结构网格,在非结构网格中,每个小面上的流通量依然可以用相邻的变量来进行很好的定义.有限元方法也是把空间分为很多小的体积, 相当于很多小的单元,然后在每个单元里,变量和流通量都用势函数来表示,计算的变量都势这些势函数中的系数.有限元的方法被普遍使用在非结构网格中,,使用结构网格也没有什么太特别的优势.1.2网格的性质.对很多种类的网格而言,有如下的特征我们非常关心的1.网格点的密度.很高的网格点的密度也就意味着很高的准确度,但是计算的时间也会长很多,这就导致了适应性网格方法的产生,这将在第四章中介绍.2.点分布的光滑程度.(我理解的是均匀的程度),对于网格的密度或者形状如果有比较大的变动,就会导致数值的发散,或者不发散但是大的误差,或者波动的传播的转向,这都会导致不稳定或者不准确3.网格形状,举个例子,流动中的边界层需要的网格是非常的扁平的(我理解的是非常的细小)并且垂直于流动的方向,若在有限元的方法使用的三角形单元,那么它的最大的角度就必须是受到限制的,并且应该严格的小于180来保障随着单元尺寸的减小这种办法的收敛性.对于简单的流动区域,对结构或者非结构的网格的选择主要是出于对离散办法的考虑, 但是,对于复杂的区域,(例如涡轮旋转机械),不规则的网格就比规则的网格的适用性更好,不规则的网格的生成(至少是用三角形或者四边形)可以充分的填满,并且速度迅速,对于规则的网格的生成则需要把区域划分为简单的可以自动划分网格的小块,这些小块的分解是半自动的(?),并且需要用户的思考努力.2.结构网格.这一章将就边界适应性网格和对他们展开PDEs 的离散方法展开讨论,然后,会处理对于简单区域中的网格生成方法(使用代数的或者差分方程的办法),并且解释用于复杂的区域的混合块的概念,2.1边界适应性网格结构网格是有着规则的特性的,网格的节点可以被标识,并且每个相邻的点都可以被计算而不是被寻找(例如,I,j这个点可以通过i+1.j 和i-1,j计算得到),在矩形区域的网格是非常容易生成的,(尽管在离散过程中对存在拐点的地方要特别地仔细),所以结构网格的生成的技术重点是放在了对带不规则的边界的区域进行的网格划分上面.例如,在血管中的流动,或者变形,金属在磨具中成型时候的流动和热传导,一般的网格都为了适应边界条件而生成,第一种办法是用一维坐标来构成边界(什么意思?)这样在边界层周围给出精确的解,并且可以使用快速的精确的求解方法,对于流体的流动,这些网格也可以支持简单的湍流模式的运用,而这些湍流模式一般是需要网格和边界保持重合的,另外的一种选择就是使用紧贴边界的矩形的网格(笛卡尔网格),在奇异的边界附近使用局部网格来完善,这样就可以减少在边界上的截取位数(?),并且可以使得网格单元方程的贴近边界,增加了算法的复杂程度,,笛卡尔网格产生的非常的快,并且可以用在欧拉空气动力学里面,但是对N-S方程的系统的适用性不好.最常见的产生贴体(?)网格的办法就是使用一个连续的网格(什么叫做连续的网格?)适应所有的边界,这样的作用是使得和矩形计算区域相邻的有着曲折边界的物体之间,可以比较好的连接.一般来说,在不产生倾斜的网格的基础上想要用正交计算区域中的分布办法来填充复杂的区域是非常的困难的.为了避免这个问题的产生,我们把主要的区域分解成为几个小块,然后对每个小块进行划分,,对各个小块的交接处,有连续性的要求,这就是混合块的概念了,把主要区域分解成为小块一般都需要用CAD技巧手动的进行实现,并且比较的慢.另外一种使用这种连续贴体多块网格的方法就是,在每个边界附近使用贴体网格,并且在内部区域采用简单的矩形网格，在连接处通过插值的办法实现，这叫做组合（嵌合网格）这样的网格对于多块网格而言比较容易产生,因为每个网格划分都是局部的,不需要和其他的网格进行匹配,独立的网格有着很高的质量(比较小的扭曲),尽管如此,差值连接是一个比较困难的事情,特别是多于二种以上的网格的搭接,并且增加了求解器的工作时间,搭接的网格不能在求解的时候太不同,否则在粘性流动计算的时候,就会在边界层附近的网格上出现问题.(?)但是组合网格方程的适合移动的边界,例如直升机的螺旋桨,或者多样的边界例如在流体种的颗粒,计算的时候,只有网格之间的差值会随着边界附近的网格的变动而改变,其他的大多数的部分都是固定不动的.组合网格的确是有一定的优点,,最近在差值收敛方面的一些工作扩充了它的使用的范围,尽管是这样,大块的结构网格的生成是基于混合块网格,我们接下来将把重点放在这上面.2.2在弯曲网格处难处理的解一旦网格生成了,接下来的问题就是离散并且求解网格,直接在网格上使用有限元或者有限体积公式来计算是可以的,但是这将会因为网格的不统一而减少FD/FV的截断位数(什么意思?),首选的办法是用问题的模型把方程转化成为可以计算的区域,因为物理的网格是根据矩形网格的转化而定义的(什么意思),这个过程是直接的,在计算区域内的方程就包括了坐标变化的内容,例如要推导的运动网格的时间(?),有限元法可以在计算网格上对转换后的方程进行求解,因为计算的网格通常都是规则的,高秩序的方法可以被使用,并且一旦FD/FV方程被转回到物理网格,截断的位数可以保存(总觉得很别扭什么保存乱七八糟的),(流动计算可能会向使用组合网格的方向发展,因为在每个网格之间的差值需要BOOK-KEEPING(?)尽管FD/FV方法的位数是可以保存的,但在解法的精确度方面曲线网格还是有一定的效果.1.网格的尺寸的变化导致数值的发散(或者不发散但是会不稳定),当网格尺寸小,求解的梯度大的时候,效果最差.(不是说网格越密越好吗)2.网格的非正交性,或者倾斜性都会增加截断误差系数,但是却不能改变位数,在相对于网格中心,小于45度角的倾斜都是可以接受的,但是一边的差分在必须保持正交性的边界上使用.这些对倾斜性的要求是多块划分思想应用的主要的约束.3.在多块系统中,各块在角落的汇合处一般使用非标准的连通来连接(例如对有限体积法离散这些点就需要在物理空间上进行离散)2.3在单块上的边界适应性网格对于简单的区域,一个简单的不过分倾斜的网格是可以被使用的,计算的取于是正交的网格或者立方的,或者至少有这正交的网格划分的边界,我们需要定义一种从计算空间转到物理空间的一种画法,这种办法正在被使用的有:1.代数网格生成法,网格是用物理边界的差值来计算的2.PDE方法3.变网格生成,一些性质(例如光滑)是最大或者最小(变网格生成在此不讨论)4.其他的办法例如投影画法2.4代数网格生成,差值法在网格生成中使用的代数的办法称为无限差值法TFI2.4.1TFI我们首先看二维的情况,取一个计算区域,基本单元矩形是[0,1]*[0,1],坐标是s和t,物理的计算区域的坐标是x y ,为了在物理计算区域内生成网格,我们可以在单元正方形中创造网格,并且将它移植到物理区域中,对这样的画法有二个要求1.必须是1*1的区间2.计算区域的边界必须可以转画到物理空间的边界上去,但是画的过程是任意的,尽管这样可能不能产生比较好的网格.TFI是画法的一种,它把物理坐标看作是计算坐标的一个函数,并且根据计算区域的边界上的值来进行差值计算,因为给出值的点的数量是不可数的(我理解的是一条边界点不可数) ,所以把这种差值的办法称为无限差值法.2维的差值是由2个一维的线性差值来组合成的(也叫做投影法?),首先定义混合函数φ0 φ1θ0 θ1 是用来表述差值的方向的,φ0 用来从S=0的边界向主要的区域差值,同理φ1是从边界s=1向主要的区域差值,相似的对于在t方向的θ0, θ1,对φ0 ( s) φ1( s )的要求是相似的对θ0 θ1在T方向,最简单的办法是进行线性差值,给出线性差值的表达式一维的无限差值可以用以下的构造他们的乘积为代表了X在计算区域四周的差值,二维的无限差值公式是它是对整个边界的差值,为了形成物理空间的网格,这个差值方法被用来画计算区域的规则网格的节点将TFI的方法推广,使得它的差值不仅仅限于边界,还可以借助几条坐标线,是可能的,这样就对控制网格的密度和保证矩形划分比较的有用,混合函数也会是立方体的线(这里翻译不确定),也可以匹配在边界上的正交的相应量,例如, dx/ ds ,这就可以容许在混合块连接的时候网格适度的倾斜.TFI也有一些问题,比如划分将传播边界处的奇异点到内部的区域,这将对流体的流动模拟产生不好的影响,更严重的问题是,如果我们的划分不是完全的正方形,溢出就会产生,这可以通过边界参量的重新确定来进行修正,或者在区域内添加约束线.2.5PDE网格生成在这些思想的背后都是想在根据计算空间坐标确定的物理空间坐标上来定义PDE,然后算出在计算区域上的网格的方程来对应出物理的空间的网格,PDE使用的主要的形式是椭圆的,(比较适合有着封闭边界的区域,对于无边界流动,在很远处的大边界常被引为假设)另外的一种是双曲线型和抛物线型方程.通常用于无边界的流动.2.5.1椭圆型网格的生成为了说明这种办法,我们来看,计算区域的坐标是ξ1 ξ2 相应的物理空间中的坐标是X1.X2,最简单的例子是拉普拉斯方程物理边界形成的等位面.,在二维情况下,极值原理保证了划分的是正方形,拉普拉斯算子是光滑的,所以边界处的网格的倾斜不会传入到内部.(什么意思)为了对当前的方程的形式求解,我们需要在物理区域中的网格,所以方程被转化了,所以X变成了不独立的变量,而ξ变成了独立的变量.使用下面的方程并且张量gij可以用来表示.这个类似于线性的方程可以在正交的iξ区间,用物理边界点的位置,对应于点在计算边界上的值来进行计算.任何标准的计算器都可以被使用,并且代数生成的网格可以用来做最开始的尝试以保证其收敛性.基于拉普拉斯方程的网格不是很容易弯曲,坐标线是等距离的,除了在凹凸不平的边界上,所以不容易控制边界上网格的倾斜(因为这是二次方程)为了控制网格的间距和倾斜度,我们增加了一个源项,(控制函数),叫做泊松方程.控制函数的作用在下图中做了说明.这个张量被分解位源项的形式表达如下变换方程就变成了通常在各个方向只需要做一维的拉伸,源项表达式是ii i ii P g 没有加,理论上,源项的表达式是根据物理空间中的坐标来定义的,但是实际上都是根据计算坐标来定义的,所以这就会导致网格点比固定的物理点会多了吸引或者排斥的问题.这使得直接使用控制函数比较的困难,通常用间接的办法来达到一系列的影响.1. 因为非常的光滑,在区域的内部坐标线是等距离分布的.为了让边界点的间距设置可以使影响传到内部,我们可以根据边界数据来增加源项以获得需要的间距,然后进行差值2. 边界上的点的间距是固定的,一维网格生成方程(拉普拉斯方程)是二次的,所以不能控制边界上的坐标线的倾斜度,一种可行的办法是使用双调和方程而放弃拉普拉斯方程,但是更好的办法是使用控制函数,在STEGER 和SORENSON 的办法中,控制函数是反复的调整直到椭圆方程可以解出,并且在边界上的网格线的倾斜度,边界上的网格点分布距离,还有边界上的第一条平行网格线距离,都可以被控制,这样就可以给出边界处的正交网格,也可以用来控制 用来修补的网格的倾斜度和间距分布的连续性(多块系统中的嵌入网格,稍后有一种更好,但是速度慢一些的匹配边界的方法)2.6三维执行在三维的情况下,我们的计算区域是立方体,并且网格生成过程有着顺序:线,面,体,首先我们选择计算网格,划分边界的使用一维泊松方程计算,然后面网格就可以产生了,有点类似于二维网格的生成,但是表面的弯曲处需要被考虑(面通常都会在CAD 系统中被定义位二维参量空间的函数,例如(不会翻 …….)最后,三维网格可以生成了,使用面网格作为数据来源.在椭圆网格的划分中有如下的一些困难1. 转化方程是非线型的,所以如果最开始给定一个好的初值,计算的速度会快很多,会很快的收敛2. 可能回收敛的很慢或者不收敛,对于比较大的控制函数,比较扭曲的角落,和高度拉伸的网格.3. 这种网格方法对于三维的网格划分而言太慢了,一般在三维问题中采用代数的技术.2.7其他的办法2.7.1双曲线的网格生成双曲线的网格生成办法,在二维中第一个方程强调了正交性,第二个方程控制了单元区域,这种方法在三维中也是可以运用的,他们有着这样的性质1.这样产生正交的网格2.不会自动的修复边界的中断3.可以应用于外部的区域(例如机翼绕流),网格坐标的解法器和边界的匹配4.对于凹边界可能会失败(因为计算进行时,网格线的下凹,单位体积会变成负数)5.速度比较快,是椭圆生成法的1倍2.7.2抛物线网格生成法抛物线方法,集合了双曲线生成法的速度,并且辅之于更好的光滑度,但是,似乎没有被广泛的应用2.8混合块理论上说,复杂的几何形状可以划分为一个正交的区域(这里可以使用于二维和三维)(可能不是简单的连接)但是这会导致不能接受的网格的单元的扭曲,实际上,物理区域被分解称为小片,每一个小片都有着相对简单的正交划分办法,这些小块被一些在它们的交互面上拥有连续性的网格联系在了一起,从无到所有(看起来象是单个的网格没有倾斜或者间距中断),所以网格生成的过程被分成了二部分,把物理的区域分解称为很多块,然后对每一块进行划分,分解的过程还不能做到完全的自动,需要和使用者的交互,(比如选择块的边界来匹配物体的边界)来生成好的网格,块的划分可以自动的进行,使用以上讨论的任意一种办法.在我们关注混合块之前,我们来看看一些可以用来适应目标的可以插入混合块使之称为一个块的简单的正交块.2.8.1C,O,H网格C.O.H网格,这样命名的原因是他们近似的形状,是正交块的简单的外延,计算的区域虽然一般都只是简单的矩形,但是计算区域的角落没有必要姚和物理空间的角落对应,反之亦然.这些经过转化后的点在差分公式计算的时候需要注意.计算的块可以被多样的连接.但是我们也可以进行剪切把它们变成简单的连接的块,用矩形来描述,这就产生了C,O,H网格(事实上H网格包括了这种多样连接的块)O型网格C型网格经常用于机翼这些有着圆滑的边界和锐利的尾部边界的形状,比H网格在这上面的应用要好.H网格是一个笛卡尔类型的网格,但是可以允许多样性的连接,也可以把内部的宽的片分成小的片]坐标的方向和形式在转化后都可以变化,这就意味着在物体的差分方程处理中需要做特殊的处理,这种复杂化可以通过使用一层halo点来避免,halo保持了变换前的连接点的复制,这就和在主要分解中为了平行处理而设置halo是一样的,但是在图十四中,边界是一条线/面的网格点,在有限体积格式中,数据点以网格单元为中心,在这个例子中halo和在区域分解中几乎一样.一旦halo设定了,相同的差分公式就可以在内部和切口上使用了.2.8.2混合块在合成块中使用COH,我们就可以生成符合的网格,混合块.第一步就是把物理区域分割称为很多四条边的小自区域,这个过程需要有经验的用户来进行区域的划分(例如有着急剧变化特征的边界)第二步是自动的,在每个子区域中生成网格,然后把它们聚合成一个整体,各个子区域之间的边界上的网格线的连续程度可以任意,甚至可以是网格点不在边界上.明显的,连续性较少,当处理子区域的时候需要做的差值计算就更多.这就就会在损失一些准确度,如果在每个块中使用代数的办法,网格线的连续性就可以得到保证,在产生线,面体的步骤中,每一步都使用上一次产生的边界点值,并且可以使用TFI来生成倾斜网格,这样保证了分布的连续性.如果PDE椭圆的办法被使用,方程就可以在整个区域上求解,使用halo把物理坐标信息从一个块转化到另外一个块,这就给出了网格坐标的连续性,并且区域边界也可以得到调整,解整个区域的画非常的耗时,所以另外一个办法就是解每一个小块,然后用STEGER和SORENSON型的控制函数来保证倾斜和分布的连续性.图十五显示了一个简单的混合块的例子,对于整边界适应混合块网格有如下的一些优点1.分块的时候需要用户的经验,是一个比较繁琐的过程,2.一个块上的几何形状的改变会导致其他的块的改变3.改变一个网点对块的贡献,例如,在变化剧烈的目标上增加一个点,将会对其他的块产生变化,如果保持网格点的连续性不变的话4.需要的网格点的连续性导致要在不增加相邻块的复杂程度上对一个块进行简化是非常困难的.(resolution怎么翻译)2.9小结结构网格的主要的产生的办法就在每个小块使用TFI或者椭圆生成法来组合成;混合块,结构网格的规则性决定了可以应用快速的解法,并行计算比较的方便(parallelization….),但是用户在组建混合块分解的时候需要很大的努力,目前复杂构造的趋势是使用非结构网格,它可以自动的生成,但是计算的要慢一些.。
基于栅格法的多体六面体网格自动生成
基于栅格法的多体六面体网格自动生成I. 引言A. 研究背景B. 研究意义C. 研究目的D. 研究内容II. 相关技术A. 栅格法概述B. 多体六面体网格生成技术C. 网格质量评估方法III. 栅格法模型拟合A. 模型建立B. 模型参数调整C. 拟合结果分析IV. 网格优化算法A. 网格初始生成B. 网格光滑处理C. 网格剖分处理V. 网格质量评估A. 基本质量指标B. 特定质量指标C. 质量评估方法比较VI. 结论与展望A. 结果回顾B. 存在问题C. 进一步研究建议D. 结论总结附：参考文献I. 引言在工程设计领域中，六面体网格作为一种最基本的网格类型，被广泛应用于有限元分析、流体力学模拟等领域。
随着计算机科学的发展，大规模的数字化模型制作和分析已经成为一种趋势。
因此自动化的六面体网格自动生成已成为当前探究的重点之一。
然而，在多体模型生成六面体网格的过程中，传统方法往往需要大量的人工操作，且时间耗费巨大，施工成本高。
因此，许多学者和工程师致力于研发一种高效、准确的自动化六面体网格生成方法，以降低施工成本和提高产品品质。
本文旨在介绍一种基于栅格法的多体六面体网格自动生成方法，并对其关键技术进行详细的讨论和探究。
本文以模型拟合、网格优化算法、网格质量评估方法等作为关键研究点，对方法进行了理论分析、实验验证和实际应用。
本文的组织结构如下：首先，介绍研究背景，包括六面体网格生成的意义和当前的研究现状。
其次，分析基于栅格法的多体六面体网格自动生成技术，包括栅格法概述、多体六面体网格生成技术、网格质量评估方法等。
接着，详细介绍栅格法模型拟合技术，包括模型建立、模型参数调整和拟合结果分析。
然后，阐述网格优化算法，包括网格初始生成、网格光滑处理和网格剖分处理。
最后，介绍网格质量评估，包括基本质量指标、特定质量指标和质量评估方法比较。
本文旨在为六面体网格生成领域的学者和工程师提供参考，让他们更好地了解基于栅格法的多体六面体网格自动生成方法，并为这一领域的进一步研究提供基础。
网格生成技术及应用PPT学习教案
引入与网格线相 交的边界点作为 附加的计算节点
凡是与直角坐标网格线倾斜 相交的边界，采用该网格的 对角线作为计算边界
无论网格划分的多 细，这些边界总是 充满锯齿形尖角
可改善模拟不规则边 界的光滑性，但易引 起计算数值不稳定性
第10页/共38页
实现了网格生成的自动化， 应用于有限分析法，计算了 具体问题，取得较好结果
微分方程法
微
椭圆型方程方法
应用最广
分
方
双曲型方程方法
程
法
抛物型方程方法
第21页/共38页
——微分方程法——
椭圆型方程方法
y
S3 S4 S1
S2 x
已知条件：
➢计算平面上ξ,η方向的节点总数和节点 位置；
➢物理平面计算区域边界上的节点设置 ，反映 出网格 疏密布 置；
S2 S4 S1
第22页/共38页
第23页/共38页
0,在S1上 1,在S 2上 ＝0 ( x, y),在S4上 ＝1 ( x, y),在S 3上
椭圆型方程方法
——微分方程法——
由于物理平面上的边界线都是曲线， 确定边 界条件 比较困 难，故 用ξ,η为 独立变 量， x,y为因变 量来建 立微分 方程， 推导过 程：
dx x d x d dy y d y d d xdx ydy d xdx ydy
第14页/共38页
代数法
———多面法
Z i 1 (r )
y
Zi (r)
x
s r
Vi (r, si )
ZN (r, sN )
在ZN，Z1两固定边界
之间生成辅助表面Z2…ZN-1
Zi1(r, si1 ) ，0<r<1，把相邻两表面上
在word中如何通过网格创建表格
在word中如何通过网格创建表格
在word中如何通过网格创建表格
表格由水平的行和垂直的列组成，行与列交叉形成的方框称为单元格。
那么在word中怎么通过网格创建表格?下面店铺就为大家详细介绍一下，来看看吧！
方法/步骤
1，启动word2007，单击“office”按钮，在弹出的对话框中，双击“另存为”选项。
2，在弹出的“另存为”对话框中，设置文档保存路径和名称。
3，单击“插入”选项卡下“表格”组中的“表格”下拉按钮。
4，这时则弹出一个下拉面板的表格。
5，在弹出的下拉面板的'网格中移动鼠标，选择4X4网格。
6，选择好网格后，单击鼠标左键，确定网格的选择，这时，则将4X4的表格插入到了文档中。
以上就是在word中怎么通过网格创建表格方法介绍，操作很简单的，大家学会了吗？希望不会的朋友看完本篇文章能知道怎么操作！
下载全文。