自适应网格生成

CFD网格及其生成方法概述

CFD网格及其生成方法概述作者：王福军网格是CFD模型的几何表达形式，也是模拟与分析的载体。网格质量对CFD计算精度和计算效率有重要影响。对于复杂的CFD问题，网格生成极为耗时，且极易出错，生成网格所需时间常常大于实际CFD计算的时间。因此，有必要对网格生成方式给以足够的关注。 1 网格类型网格（grid）分为结构网格和非结构网格两大类。结构网格即网格中节点排列有序、邻点间的关系明确，如图1所示。对一于复杂的儿何区域，结构网格是分块构造的，这就形成了块结构网格（block-structured grids）。图2是块结构网格实例。图1 结构网格实例图2 块结构网格实例与结构网格不同，在非结构网格（unstructured grid）中，节点的位置无法用一个固定的法则予以有序地命名。图3是非结构网格示例。这种网格虽然生成过程比较复杂，但却有着极好的适应性，尤其对具有复杂边界的流场计算问题特别有效。非结构网格一般通过专门的

程序或软件来生成。图3 非结构网格实例 2 网格单元的分类单元（cell）是构成网格的基本元素。在结构网格中，常用的ZD网格单元是四边形单元，3D网格单元是六面体单元。而在非结构网格中，常用的2D网格单元还有三角形单元，3D 网格单元还有四面体单元和五面体单元，其中五面体单元还可分为棱锥形（或楔形）和金字塔形单元等。图4和图5分别示出了常用的2D和3D网格单元。图4 常用的2D网格单元图5 常用的3D网格单元

3 单连域与多连域网格网格区域（cell zone）分为单连域和多连域两类。所谓单连域是指求解区域边界线内不包含有非求解区域的情形。单连域内的任何封闭曲线都能连续地收缩至点而不越过其边界。如果在求解区域内包含有非求解区域，则称该求解区域为多连域。所有的绕流流动，都属于典型的多连域问题，如机翼的绕流，水轮机或水泵内单个叶片或一组叶片的绕流等。图2及图3均是多连域的例子。对于绕流问题的多连域内的网格，有O型和C型两种。O型网格像一个变形的圆，一圈一圈地包围着翼型，最外层网格线上可以取来流的条件，如图6所示。C型网格则像一个变形的C字，围在翼型的外面，如图7所示。这两种网格部属于结构网格。图6 O型网格图7 C型网格 4 生成网格的过程

ANSYS之三维自动体网格-弯管流动

ANSYS之三维自动体网格生成—弯管流动实例问题描述：弯管流动是工程实际中常遇到的问题，弯管如图所示，管内流体为水，入口速度为10m/s，出口为自由流出口，弯管具体几何尺寸如图所示。一、创建几何模型（1）设定工作目录 File→ChangeWorking Directory,选择文件存储目录（2）创建Point Step1：创建Point0，选择Geometry，如图所示，输入X=0，Y=0，Z=0，创建Point0

Step2：创建Point1，选择Geometry，如图所示，输入DX=0，DY=0，DZ=500选择Point0为基准创建Point1 Step3：采用相似的方法创建其余各点，其中创建不同点时的基准点和偏移量如表所示

（3）创建Curve，如图所示，选择Geometry标签栏，创建Curve Step1：创建Curve1，右击模型树Model→Points，选择Blank Points，再选择Point1，中键确定，此时Point1被隐藏。将鼠标置于主窗口Z轴坐标附近，当出现+Z时，单击中键，主窗口将显示Z轴视图。如图所示，勾选Radius复选框，输入半径为600,。单击选择圆心Point0，然后在其附近任意点两点，得到Curve1。 Step2：创建Curve2，如图所示，依次选择Point0和Point1，中键确定。

Step3：创建Curve3，在Method下拉列表中选择From 3Points，单击依次选择Point1，Point2，Point3，中键确定。 Step4：创建Curve4，与Step2的方法相同，连接Point3和Point4创建Curve。（4）创建Surface，选择Geometry标签，单击进入创建Surface的操作。 Step1：创建弯管的侧面。如图所示，单击Curve2为Driving curve，单击选择Curve1为Driven curve，中键确定，如下图所示。 Step2：采用与Step1相同的方法，依次生成剩余圆柱的侧面。

ansys教程之自适应网格划分

ansys教程之自适应网格划分 [摘要]：ANSYS程序提供了近似的技术自动估计特定分析类型中因为网格划分带来的误差。（误差估计在ANSYS Basic Analysis Procedures Guide第五章中讨论。）通过这种误差估计，程序可以确定网格是否足够细。如果不够的话，程序将自动细化网格以减少误差。这一自动估计网格划分误差并细化网格的过程就叫做自适应网格划分，然后通过一系列的求解过程使得误差低于用户指定的数值（或直到用户指定的最大求解次数）。自适应网格划分的先决条件 ANSYS软件中包含一个预先写好的宏，ADAPT.MAC，完成自适应网格划分的功能。用户的模型在使用这个宏之前必须满足一些特定的条件。（在一些情况下，不满足要求的模型也可以用修正的过程完成自适应网格划分，下面还要讨论。）这些要求包括：标准的ADAPT过程只适用于单次求解的线性静力结构分析和线性稳态热分析。模型最好应该使用一种材料类型，因为误差计算是根据平均结点应力进行的，在不同材料过渡位置往往不能进行计算。而且单元的能量误差是受材料弹性模量影响的。因此，在两个相邻单元应力连续的情况下，其能量误差也可能由于材料特性不同而不一样。在模型中同样应该避免壳厚突变，这也可能造成在应力平均是发生问题。模型必须使用支持误差计算的单元类型。（见表3－1）模型必须是可以划分网格的：即模型中不能有引起网格划分出错的部分。表3-1 自适应网格划分可用单元 2-D Structural Solids PLANE2 2-D 6-Node Triangular Solid PLANE25 Axisymmetric Harmonic Solid PLANE42 2-D 4-Node Isoparametric Solid PLANE82 2-D 8-Node Solid PLANE83 Axisymmetric Harmonic 8-Node Solid

网格生成及修正技巧

网格生成及修正技巧 1引言网格是CFD 模型的几何表达形式，也是模拟与分析的载体。网格质量对CFD 计算精度和计算效率有着重要的影响。对于复杂的CFD 问题，网格的生成极为耗时，并且极易出错，生成网格所需的时间常常大于实际CFD 计算的时间。因此，有必要对网格生成以及修正方法进行足够的研究。考虑到目前的CFD 计算多是通过专用的网格生成软件来划分所需要的网格，因此，本文就如何利用专用前处理软件GAMBIT 来介绍网格的生成和修正技巧。 2 网格类型网格主要有两种：结构网格和非结构网格[1] [2]在结构网格中，常用的2D 网格单元是四边形单元，3D 网格单元是六面体单元。而在非结构网格中，常用的2D 网格单元还有三角形单元，3D 网格单元还有四面体单元和五面体单元，其中五面体单元还分为棱锥形(或楔形)和金字塔形单元等。结构网格的最大特点在于网格中节点排列有序，邻点间关系明确，结构简单，构造方便，与计算机语言自然匹配，容易计算，网格生成速度快，质量好，数据结构简单等优点；缺点是适用的范围比较窄，只适用于形状规则的图形，对复杂几何形状的适应能力差。非结构网格舍去了网格节点的结构性限制，易于控制网格单元的大小、形状及节点位置，灵活性好，对复杂外形的适应能力强——流场变化比较大的地方，可以进行局部网格加密。但其无规则性也导致了在模拟计算中存储空间增大，寻址时间增长，计算效率低于结构化网格，计算时间长等缺点。 [1]。（a ）三角形（b ）四边形图1 常用的2D 网格单元（a ）四面体（b ）六面体（c ）五面体（凌锥）（d ）五面体（金字塔）图2 常用的3D 网格单元 3 单连域与多连域网格网格区域分为单连域和多连域两类。所谓单连域是指求解区域边界线内不包含有非求解

有限元设计软件生成网格的PAVING算法

有限元设计软件生成网格的PAVING算法一、简介使用有限元软件分析计算几何体的物理性质，其计算的过程可以划分为几个大的模块，输入几何体区域→为该区域生成一个网格→对生成的网格施加一个干扰→从受到干扰的网格开发分析数据→确定几何体的物理行为。分析计算流程图如图1所示。图1 有限元分析的模块在有限元分析的前处理模块中，网格生成时很重要的一个步骤。生成网格的质量会影响后处理计算结果的精度。当前，行业内流行多种网格生成的算法，各有各自的特点，该部分内容在本文国内外研究现状一节中已经详细阐述。其中paving算法健壮性良好，计算速度快，而且生成的网格质量好。本节主要阐述采用C++语言实现paving算法的实现过程。如图2所示，采用paving算法生成网格的算法流程图，从输入边界数据到最后输出划分好的网格，其中主要有生成新行，平滑处理，缝合处理，边界相交处理等几个子模块。

图2 Flow chart of paving algorithm 为了清晰理解上述paving算法的流程，以图3所示为例，图a当中为输入的原始外边界数据，围成待划分网格的区域。选择边界上的一行节点为基础，添加生成一行新的浮动节点，生成顺序为沿着外边界按逆时针方向进行。对新生成的浮动节点进行平滑处理，使节点围成的单元的internal angle以及aspect ratio变得更为合理，单元更趋近于规则四边形。对剩余的待划分网格区域进行缝合，检查单元是否相交，对相交的单元进行处理，对单元进行调整，直到整个区域生成高质量的网格为止。图3 paving算法铺筑单元示意图依据图2所示流程图，生成相应的伪代码： Do Row choise While add row is not complete Add row portion Smooth row portion Seam boundary If intersection occurs then Connect overlaps Seam boundary End if Row adjustment

自适应网格

ALE adaptive mesh单元： AC1D2, AC1D3, AC2D3, AC2D4, AC2D6, AC2D8, AC3D4, AC3D6, AC3D8, AC3D10, AC3D15, AC3D20, ACAX3, ACAX4, ACAX6, ACAX8 CPS4, CPS4T, CPS3 CPE4, CPE4H, CPE4T, CPE4HT, CPE4P, CPE4PH, CPE3, CPE3H CAX4, CAX4H, CAX4T, CAX4HT, CAX4P, CAX4PH, CAX3, CAX3H C3D8, C3D8R, C3D8H, C3D8RH, C3D8T, C3D8HT, C3D8RT, C3D8RHT, C3D8P, C3D8PH, C3D8RP, C3D8RPH 从列表来看，ALE自适应网格不适用于壳（S4、S4R等），另外对于实体单元也不适用于四面体（C3D4）。问题1: The requested number of domains cannot be created due to restrictions in domain decomposition. 措施：job---Editjob---Parallelization---Number of domains: 设为1 问题2：ALE算法和CEL算法有什么区别？措施：①CEL只能用于explicit，AEL在implicit（声畴、冲蚀、磨损）和explicit都能用； ②ALE方法最初出现于数值模拟流体动力学问题的有限差分方法中。这种方法兼具 Lagrange方法和Euler方法二者的特长，即首先在结构边界运动的处理上它引进了 Larange方法的特点，因此能够有效的跟踪物质结构边界的运动；其次在内部网格的划分上，它吸收了Euler的长处，即是使内部网格单元独立于物质实体而存在，但它又不完全和Euler网格相同，网格可根据定义的参数在求解过程中适当调整位置，使得网格不致出现严重的畸变。 CEL是欧拉-拉格朗日耦合，用于固体液体之间的耦合。说法1：ALE是arbitary lagrange euler 算法 CEL couple lagrange euler 流固耦合的设置应该不属于算法的范畴问题3：为什么odb转换输出坐标系后，只有S11等应力分量改变，而像Mises Equivalent 等都不变呢？措施：①看变量情况：S11指的是沿一方向的力，改了坐标系，值也会变。但是如果是CPRESS，接触压力的话，它是指垂直于接触面上的力，与你的坐标没有关系，这样，你改了坐标系，自己它不会变。其它的变量如mises等同理。

网格生成技术

I 目录 1 概述 (1) 2 结构网格 (3) 2.1 贴体坐标法 (3) 2.2 块结构化网格 (11) 3 非结构网格 (16) 3.1 概述 (16) 3.2 阵面推进法 (16) 3.3 Delaunay三角划分 (19) 3.4 四叉树（2D）/八叉树（3D）方法 (21) 3.5 阵面推进法和Delaunay三角划分结合算法 (22) 4 其他网格生成技术 (23) 4.1 自适应网格 (23) 4.2 混合网格 (25) 4.3 动网格 (26) 4.4 曲面网格 (27) 4.5 重叠网格 (28) 5 网格生成软件 (29) 5.3 Gambit (29) 5.2 ICEM CFD (30) 5.1 TrueGrid (32) 5.2 Gridgen (34)

1 概述计算流体力学作为计算机科学、流体力学、偏微分方程数学理论、计算几何、数值分析等学科的交叉融合，它的发展除依赖于这些学科的发展外，更直接表现于对网格生成技术、数值计算方法发展的依赖。在计算流体力学中，按照一定规律分布于流场中的离散点的集合叫网格（Grid），分布这些网格节点的过程叫网格生成（Grid Generation）。网格生成是连接几何模型和数值算法的纽带，几何模型只有被划分成一定标准的网格才能对其进行数值求解，所以网格生成对CFD至关重要，直接关系到CFD计算问题的成败。一般而言，网格划分越密，得到的结果就越精确，但耗时也越多。1974年Thompson等提出采用求解椭圆型方程方法生成贴体网格，在网格生成技术的发展中起到了先河作用。随后Steger等又提出采用求解双曲型方程方法生成贴体网格。但直到20世纪80年代中期，相比于计算格式和方法的飞跃发展，网格生成技术未能与之保持同步。从这个时期开始，各国计算流体和工业界都十分重视网格生成技术的研究。上个世纪90年代以来迅速发展的非结构网格和自适应笛卡尔网格等方法，使复杂外形的网格生成技术呈现出了更加繁荣发展的局面。现在网格生成技术已经发展成为CFD的一个重要分支，它也是计算流体动力学近20年来一个取得较大进展的领域。也正是网格生成技术的迅速发展，才实现了流场解的高质量，使工业界能够将CFD的研究成果——求解Euler/NS方程方法应用于型号设计中。随着CFD在实际工程设计中的深入应用，所面临的几何外形和流场变得越来越复杂，网格生成作为整个计算分析过程中的首要部分，也变得越来越困难，它所需的人力时间已达到一个计算任务全部人力时间的60%左右。在网格生成这一“瓶颈”没有消除之前，快速地对新外形进行流体力学分析，和对新模型的实验结果进行比较分析还无法实现。尽管现在已有一些比较先进的网格生成软件，如ICEM CFD、Gridgen、Gambit等，但是对一个复杂的新外形要生成一套比较合适的网格，需要的时间还是比较长，而对于设计新外形的工程人员来说，一两天是他们可以接受的对新外形进行一次分析的最大周期。要将CFD从专业的研究团体中脱离出来，并且能让工程设计人员应用到实际的设计中去，就必须首先解决网格生成的自动化和即时性问题，R.Consner等人在他们的一篇文章中，详细地讨论了这些方面的问题，并提出：CFD研究人员的关键问题是“你能把整个设计周期缩短多少天？”。而缩短设计周期的主要途径就是缩短网格生成时间和流场计算时间。因此，生成复杂外形网格的

体网格生成示例

体网格生成示例以生成f6外流场体网格为例讲述生成体网格的过程，整个过程包括几何导入、网格尺寸设置、表面网格生成、构建外包围盒、生成体网格、查看内部网格和边界条件设置。 1 导入几何在菜单栏中，选择“文件”> “添加文件…”> “添加几何文件…”，弹出选择几何文件对话框。可以添加的几何文件类型有IGES、STEP、SA T、FLI和GM2(3)格式文件，其中前三种为标准数据交换格式，后两种为本软件自有几何文件格式。在本示例中，我们导入IGES格式几何文件，文件导入后会弹出图1所示对话框用于设置全局容差。图1 容差设置对话框如图2所示为导入到软件中的f6外形，图3所示为几何外形的三视图。几何导入后，我们可以查看几何信息。在菜单栏中，通过“窗口”> “拓扑信息”，可以查看几何外形包含的曲线数目、曲面数目和区域数目；通过“窗口”> “尺寸信息”，可以查看几何外形尺寸值。通过工具栏上的“选择”下拉列表，我们可以使用鼠标左键点击选取几何元

素，选中某一几何元素后，会在信息输出窗口显示选取的几何元素信息。图2 f6外形在功能操作面板上，可以通过“快速设置视点”选择多种视图查看几何；也可以通过对象视图控制，设置可视的几何元素。另外，点击“高级控制”按钮，在弹出“可视化控制”对话框中，我们可以控制特定类型和特定编号的几何元素的显示和隐藏。在视图模式下，我们可以通过鼠标控制，实现对几何外形的旋转、放缩和平移等操作。鼠标左键、中键和右键的功能分别为平移、旋转和放缩。

图3 几何外形三视图 2 网格尺寸设置我们可以通过指定曲线离散段数、设置背景网格和网格源等方式实现对网格尺寸的控制。三维项目中背景网格需要设置的唯一参数是“全局尺寸”，在“功能操作面板”的“网格”选项卡中的“密度控制”中进行设置，通过设置背景网格得到的网格尺寸是均匀的。当需要控制网格局部尺寸时，我们可以采用网格源的方式，网格源包括点源、线源和面源。二维点源是两个同心圆，三维点源是两个同心球，线源和面源分别有两个点源和三个点源组成。设置网格源需要设置的参数有：内径、外径和网格尺寸值（如图4所示），此外，在创建网格源过程中，点源需要选择1个几何点，线源需要选择2个几何点，面源要选择3个几何点。图4 网格源设置

自动网格生成法

自动网格生成法二维网格生成—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?2t3 G0t=t?2t2+t3 ; G1t=?t2+t3 此时，参数表达式可以通过一个系数矩阵来描述： X t=1,t,t2,t3M x0,x1,x,t0,x,t1T,0≤t≤1 其中M矩阵读者很容易写出，是一个4*4的方阵，而每一列是这些厄米多项式的系数排列而成。我们把这个表示称之为样本表示。每个边界都包含n个这样的数据点： x i,i=1,2,3,……,n 利用内插法可以构造出如下形式的关系式： X u=H0t x u i?1+H1t x u i+Δi G0t x,t u i?1+Δi G1t x,t u i 其中Δi是单位区间的长度。同时参数t也变为离散的取值是单位区间从原点到任意点所有的个数。如果参数的离散取值正好是i，那么u的表达式将简化为：

基于映射法的六面体网格生成算法

基于映射法的六面体网格生成算法王东风,翟建军,陈文亮 (南京航空航天大学机电学院,江苏南京　210016) 摘要:六面体网格划分技术是三维有限元仿真软件处理的关键环节之一,等参映射法既可适应特殊的区域边界形状,又可控制所生成单元的形状和密度。对基于等参映射法的六面体网格划分原理进行了深入研究,并在此研究基础上对等参映射法的计算过程进行了细致的分析,利用VC++开发了该算法的相应程序,最后给出了2个等参映射法具体的应用实例,计算结果表明该程序的计算精度已经达到了工程要求。关键词:等参映射法;六面体网格;有限元中图分类号:TP391 文献标识码:A 文章编号:1672-1616(2009)05-0025-03 在有限元仿真过程中,单元类型的选择对整个有限元仿真的计算效率、自动化程度、计算精度都将产生重要影响。六面体单元由于变形特性好、计算精度高等优点在三维有限元仿真领域中得到了广泛应用[1]。映射法是三维网格划分中最早使用的方法,和扫略法、基于栅格法等其他方法相比,该方法生成网格速度快、生成的网格单元质量好、网格密度可控制[2～4]。映射法对复杂实体生成三维有限元网格有两大难点,一是子区域划分问题,二是子区域之间网格相容性问题。Price与Armstrong等提出中面法,将三维复杂区域分解成可映射子区域[5～7],但是该算法存在一些问题,特别是几何适应能力问题。李华和程耿东提出了三维组合式模板,一定条件下解决了子区域之间的网格相容性问题[2]。还有学者提出了Embedded Voronoi Graph[8]和BLOBs[9],对复杂实体利用映射法划分六面体网格。映射法在众多有限元分析软件中占有重要地位,美国Altair公司Hyper-Mesh软件中的Solid Mesh Panel就是利用映射法生成六面体网格。本文对基于等参映射法的六面体网格划分技术进行了详细研究。通过形函数映射技术将物理域映射到参数空间域,对规则参数域进行网格剖分,将参数域的网格反向映射回物理空间,从而得到物理空间六面体网格。利用VC++实现了映射过程,在输入边界信息和划分信息后,即得到了六面体网格的节点信息和单元信息。 1　映射法生成四边形网格和六面体网格本文主要讨论的是怎样在一个子区域中划分六面体网格。这里的子区域指的是具有6个面12条边,每条边的特征点已知的区域。求子区域六面体网格节点的步骤:(1)利用积累弦长参数化法对每条边进行参数化。(2)利用拉格朗日插值公式求边界函数。(3)利用边界函数,由双线性混合孔斯曲面片公式求曲面节点坐标。 (4)由孔斯线性混合插值公式求子区域节点坐标。在计算的过程中要用到2种坐标系,即笛卡尔坐标系和自然坐标系。笛卡尔坐标系用x,y,z表示,自然坐标系用一组不超过1的无量纲参数r,s, t表示,边界点分别对应自然坐标等于1或0的点。如图1所示。图1　自然坐标和笛卡尔坐标之间的变换收稿日期:2008-08-08 作者简介:王东风(1979-),男,河南商丘人,南京航空航天大学硕士研究生,主要研究方向为CAD/CAM/CAE。

ANSYS命令流学习笔记12-自适应网格及其在WB中运用的对比

!ANSYS命令流学习笔记12-自适应网格及其在WB中运用的对比 !学习重点： !1、网格收敛的重要性由于应力集中（区别于应力奇异）的存在，在结构不连续处存在应力较大，而且随着网格质量数量的增加，应力值趋于收敛，据说收敛与否的应力差值可以很大，所以说重要细节结构的网格收敛十分重要。 !2、WorkBench中网格收敛的实现 WorkBench中在solution选项中设置网格循环次数，关键点网格优化系数，在求解结果选项下插入convergence，定义deformation或者stress的收敛系数。计算前后的网格对比虽然六面体网格变成四面体网格进行细分，但是初始网格划分的尺寸，对结果仍然有一定影响。而且优化的方式和APDL中也有一定差异，此例与APDL结果相同，是有一定运气成分。此处如果网格继续细化，肯定是fix处的尖角处应力奇异。（所以如何指定优化区域也是个问题） !3、APDL中网格收敛的实现 (1) 建模，注意不要划分网格，而且3D模型只能使用4面体单元网格； (2) 加载边界条件，由于没有网格，边界条件只好由面或者线确定； (3) 启动ADAPT宏命令，指定能量收敛误差，最大循环次数，网格优化系数；看起来很厉害的样子，但是使用方法和命令一样，只是输入命令框的不提示有此命令存在；. (4) 后处理查看结果。 !4、对网格收敛的一些疑问 (1) 宏命令调用：通过help文件查询到ADAPT的命令含义，但是不懂宏的内容，权且当做命令处理。 (2) ADAPT收敛误差：help中说是结构能量误差（SEPC），如果在热分析是热能量误差（TEPC），SEPC等效于应变能量误差（strain energy error ）。由于APDL和WorkBench收敛准则的不同，收敛结果无法对应，不明所以。 (3) 网格划分方式：实体单元只能采用非结构化网格形式，WB和APDL都是如此，WB 即使划分了结构和单元也会无效；但是APDL和WB指定网格初始尺寸有意义？；虽说非结

ANSYS自适应网格划分

ANSYS自适应网格划分 (1) 何为网格自适应划分？ ANS YS程序提供了近似的技术自动估计特定分析类型中因为网格划分带来的误差。（误差估计在ANSYS Basic Analysis Procedures Gui第五章中讨论。）通过这种误差估计，程序可以确定网格是否足够细。如果不够的话，程序将自动细化网格以减少误差。这一自动估计网格划分误差并细化网格的过程就叫做自适应网格划分，然后通过一系列的求解过程使得误差低于用户指定的数值（或直到用户指定的最大求解次数）。自适应网格划分的先决条件 ANSYS软件中包含一个预先写好的宏，ADAPT.MAC完成自适应网格划分的功能。用户的模型在使用这个宏之前必须满足一些特定的条件。（在一些情况下，不满足要求的模型也可以用修正的过程完成自适应网格划分，下面还要讨论。）这些要求包括：标准的ADAPT过程只适用于单次求解的线性静力结构分析和线性稳态热分析。模型最好应该使用一种材料类型，因为误差计算是根据平均结点应力进行的，在不同材料过渡位置往往不能进行计算。而且单元的能量误差是受材料弹性模量影响的。因此，在两个相邻单元应力连续的情况下，其能量误差也可能由于材料特性不同而不一样。在模型中同样应该避免壳厚突变，这也可能造成在应力平均是发生问题。模型必须使用支持误差计算的单元类型。模型必须是可以划分网格的：即模型中不能有引起网格划分出错的部分。自适应网格划分可用单元 2-D Structural Solids PLANE2 2-D 6-Node Triangular Solid PLANE25 Axisymmetric Harmonic Solid

有限元网格剖分方法概述

有限元网格剖分方法概述在采用有限元法进行结构分析时，首先必须对结构进行离散，形成有限元网格，并给出与此网格相应的各种信息，如单元信息、节点坐标、材料信息、约束信息和荷载信息等等，是一项十分复杂、艰巨的工作。如果采用人工方法离散对象和处理计算结果，势必费力、费时且极易出错，尤其当分析模型复杂时，采用人工方法甚至很难进行，这将严重影响高级有限元分析程序的推广和使用。因此，开展自动离散对象及结果的计算机可视化显示的研究是一项重要而紧迫的任务。有限元网格生成技术发展到现在, 已经出现了大量的不同实现方法，列举如下：映射法映射法是一种半自动网格生成方法，根据映射函数的不同，主要可分为超限映射和等参映射。因前一种映射在几何逼近精度上比后一种高，故被广泛采用。映射法的基本思想是：在简单区域内采用某种映射函数构造简单区域的边界点和内点，并按某种规则连接结点构成网格单元。也就是根据形体边界的参数方程，利用映射函数，把参数空间内单元正方形或单元三角形(对于三维问题是单元立方体或单元四面体)的网格映射到欧氏空间，从而生成实际的网格。这种方法的主要步骤是，首先人为地把分析域分成一个个简单可映射的子域，每个子域为三角形或四边形，然后根据网格密度的需要，定义每个子域边界上的节点数，再根据这些信息，利用映射函数划分网格。这种网格控制机理有以下几个缺点：（1）它不是完全面向几何特征的，很难完成自动化，尤其是对于3D区域。（2）它是通过低维点来生成高维单元。例如，在2D问题中，先定义映射边界上的点数，然后形成平面单元。这对于单元的定位，尤其是对于远离映射边界的单元的定位，是十分困难的，使得对局部的控制能力下降。（3）各映射块之间的网格密度相互影响程度很大。也就是说，改变某一映射块的网格密度，其它各映射块的网格都要做相应的调整。其优点是：由于概念明确，方法简单，单元性能较好，对规则均一的区域，适用性很强，因此得到了较大的发展，并在一些商用软件如ANSYS等得到应用。 2 。拓扑分解法拓扑分解法较其它方法发展较晚, 它首先是由Wordenwaber提出来的。该方法假设最后网格顶点全部由目标边界顶点组成, 那么可以用一种三角化算法将目标用尽量少的三角形完全分割覆盖。这些三角形主要是由目标的拓扑结构决定, 这样目标的复杂拓扑结构被分解成简单的三角形拓扑结构。该方法生成的网格一般相当粗糙, 必须与其它方法相结合, 通过网格加密等过程, 才能生成合适的网格。该方法后来被发展为普遍使用的目标初始三角化算法, 用来实现从实体表述到初始三角化表述的自动化转换。单一的拓扑分解法因只依赖于几何体的拓扑结构使网格剖分不理想，有时甚至很差。 3．连接节点法这类方法一般包括二步：区域内布点及其三角化。早期的方法通常是先在区域内布点, 然后再将它们联成三角形或四面体, 在三角化过程中, 对所生成的单元形状难于控制。随着Delaunay三角化(简称为DT ) 方法的出现, 该类方法已成为目前三大最流行的全自动网格生成方法之一。 DT法的基本原理：任意给定N个平面点Pi(i=1,2,…,N)构成的点集为S，称满足下列条件的点集Vi为Voronoi多边形。其中，Vi满足下列条件： Vi ={ X：|X- Pi|(|X- Pj|,X(R2,i(j,j=1,2,…,N }Vi为凸多边形，称{ Vi}mi=1为Dirichlet Tesselation

Fluent的自适应网格问题

加密网格的话有两种参考标准一种是y+值，一种是y*值，一般来说，要加密网格主要是为了是y+值满足需求，具体的情况看楼主你的需要... 根据y+值来加密网格的步骤如下：运行fluent，导入cas and dat 文件后，点击adapt——Yplus/Ystar..。，之后出现选择界面，一般情况可以保持默认界面，当然也可以根据自己的需求选择选项，一般type项选择Yplus，然后点击compute，在min及max项会出现你的选择壁面的Y+值，在其下方，有minallowed 和maxallowed，输入你所需要的Y+值范围，点击Mark按钮，会标记出不符合要求的部分，然后点击adapt，就可以了，这部分区域的网格会加密，以适应你的要求 Y*的步骤也是这样的但是前提是要知道你的计算的y+值范围，而这个值一般是估计值，且跟计算有关的，是个不确定量，所以一般只作参考用希望能帮到你......另外，希望给加分啊，呵呵追问我点完adpat，Yplus/Ystar这个是灰的，不能点。。回答额，你计算了吗或者说你导入的是cas & dat 文件吗如果不是，你都没有一个y+值的范围，怎么可能让软件给你加密网格...（这是基本条件）追问当然计算了，我保存完再导入cas& dat也不行回答那你试试计算完，直接点adapt试试.....还真没遇到过你说的情况追问 adapt都能点只是里面的Yplus/Ystar不能点，是灰色的 fluent里的常见问题（一） (2011-02-26 09:44:43) 1什么叫松弛因子松弛因子对计算结果有什么样的影响它对计算的收敛情况又有什么样的影响 1、亚松驰（Under Relaxation）：所谓亚松驰就是将本层次计算结果与上一层次结果的差值作适当缩减，以避免由于差值过大而引起非线性迭代过程的发散。用通用变量来写出时，为松驰因子（Relaxation Factors）。《数值传热学-214》 2、FLUENT中的亚松驰：由于FLUENT所解方程组的非线性，我们有必要控制的变化。一般用亚松驰方法来实现控制，该方法在每一部迭代中减少了的变化量。亚松驰最简单的形式为：单元内变量等于原来的值加上亚松驰因子a 与变化的积, 分离解算器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方程，包括耦合解算器所解的非耦合方程（湍流和其他标量）都会有一个相关的亚松驰因子。在FLUENT中，所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题，但是对于一些特殊的非线性问题（如：某些湍流或者高Rayleigh数自然对流问题），在计算开始时要慎重减小亚松驰因子。使用默认的亚松驰因子开始计算是很好的习惯。

有限元网格划分和收敛性

一、基本有限元网格概念 1.单元概述几何体划分网格之前需要确定单元类型。单元类型的选择应该根据分析类型、形状特征、计算数据特点、精度要求和计算的硬件条件等因素综合考虑。为适应特殊的分析对象和边界条件，一些问题需要采用多种单元进行组合建模。 2?单元分类选择单元首先需要明确单元的类型，在结构中主要有以下一些单元类型：平面应力单元、平面应变单元、轴对称实体单元、空间实体单元、板单元、壳单元、轴对称壳单元、杆单元、梁单元、弹簧单元、间隙单元、质量单元、摩擦单元、刚体单元和约束单元等。根据不同的分类方法，上述单元可以分成以下不同的形式。 3. 按照维度进行单元分类根据单元的维数特征，单元可以分为一维单元、二维单元和三维单元。一维单元的网格为一条直线或者曲线。直线表示由两个节点确定的线性单元。曲线代表由两个以上的节点确定的高次单元，或者由具有确定形状的线性单元。杆单元、梁单元和轴对称壳单元属于一维单元，如图 1?图 3所示。二维单元的网格是一个平面或者曲面，它没有厚度方向的尺寸。这类单元包括平面单元、轴对称实体单元、板单元、壳单元和复合材料壳单元等，如图 4所示。二维单元的形状通常具有三角形和四边形两种，在使用自动网格剖分时，这类单元要求的几何形状是表面模型图1捋果詰柯与一维杆单无犠型（直豉） &2桁舉第构石一隼杆早死撲型（曲线） B3毀姑构与一纯梁单元除世（直疑和呦疚〕

或者实体模型的边界面。采用薄壳单元通常具有相当好的计算效率。

洞丨伍金哉钩和潯壳社电三维单元的网格具有空间三个方向的尺寸，其形状具有四面体、五面体和六面体，这类单元包括空间实体单元和厚壳单元，如图5所示。在自动网格划分时，它要求的是几何模型是实体模型（厚壳单元是曲面也可以）。图5三址乙勺久和父侬草无 4. 按照插值函数进行单元分类根据单元插值函数多项式的最高阶数多少，单元可以分为线性单元、二次单元、三次单元和更高次的单元。线性单元具有线性形式的插值函数，其网格通常只具有角节点而无边节点，网格边界为直线或者平面。这类单元的优点是节点数量少，在精度要求不高或者结果数据梯度不太大的情况下，采用线性单元可以得到较小的模型规模。但是由于单元位移函数是线性的，单元着应力突变，如图6所示。 S6錢41吕节点点单无fu节庖实体羊元

ICEM网格生成流程

Chapter 3 二维非结构壳/面网格生成（2、3） 1. 创建几何模型：Point --- Curve --- Surface --- Part --- Topology 2．定义网格参数 2.1．定义全局网格参数 2.1.1 定义网格全局尺寸：Scale factor、Max element 2.1.2 定义全局壳网格参数：Mesh type、Mesh method 2.2 定义Part网格尺寸 3. 生成网格并导出 3.1 生成网格，检查网格质量 3.2 保存网格文件：Save mesh as… 3.3 选择求解器：Output --- Select solver 3.4 写入：Output --- Write input Chapter 4 三维非结构自动体网格生成(自上而下)（2、3） 1. 创建几何模型：Point --- Curve --- Surface --- Part --- Topology --- Body 2．定义网格参数 2.1．定义全局网格参数 2.1.1 定义全局网格尺寸：Scale factor、Max element 2.1.2 定义体网格全局参数：Mesh type、Mesh method

2.1.3 定义棱柱网格全局参数：Grow Law、Initial height、Ratio、 No. 2.2 定义Part网格尺寸 3. 生成网格并导出 3.1 生成网格，检查网格质量 3.2 保存网格文件：Save mesh as…*.uns 3.3 选择求解器：Output --- Select solver 3.4 写入：Output --- Write input 三维非结构自动体网格生成(自下而上)（4）首先导入壳网格，在壳网格的基础上拉伸生成棱柱体网格，再填充棱柱体网格和远场边界之间的空隙。（壳网格---棱柱体网格---体网格）。 1.创建（导入）几何模型 2.创建生成（导入）壳网格 3.生成棱柱体网格 3.1定义棱柱网格参数： Growth law、Initial height、Ratio、No.、New volume part（表征体网格的材料，相当于自上而下中的body） 3.2指定生成棱柱边界层的Surface（定义Part网格尺寸） 3.3生成棱柱体网格：Mesh --- Compute mesh --- Prism mesh 4.生成棱柱网格与远场边界之间的体网格

流体机械CFD中的网格生成方法进展

文章编号: 1005 0329(2010)04 0032 06 技术进展流体机械CFD中的网格生成方法进展刘厚林,董亮,王勇,王凯,路明臻 (江苏大学,江苏镇江 212013) 摘要: 网格生成技术是流体机械内部流动数值模拟中的关键技术之一,直接影响数值计算的收敛性,决定着数值计算结果最终的精度及计算过程的效率;本文在分析大量文献的基础上,首先,对流体机械CFD中的网格生成方法即结构化网格、非结构化网格、混合网格进行了比较全面的总结,系统地分析这些网格划分方法的机理、特点及其适用范围;其次,对特殊的网格生成技术,如曲面网格生成技术、动网格技术、重叠网格生成技术、自适应网格技术进行了阐述;再次,指出了良好的网格生成方法应具备的特点;最后提出了网格生成技术的发展趋势。关键词: 流体机械;网格生成;计算流体动力学;动网格;自适应网格中图分类号: TH311 文献标识码: A do:i10.3969/.j i ssn.1005-0329.2010.04.008 Overvie w onM esh Generati o n M et hods i n CF D of F lui d M achinery L IU H ou-lin,DONG L iang,W ANG Y ong,W ANG K a,i LU M i ng-zhen (Jiangsu U n i v ers it y,Zhenji ang212013,Ch i na) Abstrac t: M esh genera ti on techno logy i s one of the cr iti ca l technology f o r fl u i d m ach i nery fl ow nume rica l s i m u l at-i on,and d-i rectly i nfl uence t he astr i ngency o f nume rical si m u l a ti on,wh ich has an i m portan t e ffect on the nu m er ica l s i m u l a tion results,fi na l precision and the effi c i ency o f compu tati onal process.O n the bas i s o f analyzi ng a great dea l litera t ures,firstl y,m esh genera ti on m ethods and t heory of fluid m ach i nery are comprehens i ve l y su mm ar i zed such as structured mesh,unstructured mesh,hybrid gr i d and respecti ve re lati ve m erits and the pr i nciple,charac teristcs and scopes of t hese m ethods we re sy stema ti ca lly ana l ysed.Second-ly,Spec i a lm esh generation m ethod w ere su mm ar i zed,such as surface m eshi ng,m ov ing gr i d,adapti ve gr i d and especiall y i ntro-duced the pr i nci p le and app licati on areao f adapti ve g ri d.T h irdly,the character i sti c o f m esh g enerati on m e t hod w ere pion ted out. F i na lly,t he trends of mesh generati on are presen ted,and the tre m endous d ifference i s analyzed i n mesh au t om atic gene ra tion at a-broad and the necessary o f exp l o iti ng CFD soft w are and resea rchi ng the m esh auto m atic gene ration techn i que i n our country are put forwa rd. K ey word s: fl uids m achi nery;m esh g enerati on;co m puta ti ona l fl u i d dyna m ics;mov i ng gr i d;adaptive gr i d 1 前言计算流体动力学(CFD)中,按一定规律分布于流场中的离散点的集合叫网格,产生这些节点的过程叫网格生成。网格生成是连接几何模型和数值算法的纽带,几何模型只有被划分成一定标准的网格时才能对其进行数值求解,一般而言,网格划分越密,得到的结果就越精确,但耗时也越多。数值计算结果的精度及效率主要取决于网格及划分时所采用的算法[1],它和控制方程的求解是数值模拟中最重要的两个环节。网格生成技术已经发展成为流体机械CFD的一个重要分支。现有的网格生成方法主要分为结构化网格、非结构化网格和混合网格三大类。收稿日期: 2009 11 04 基金项目: 国家杰出青年基金(50825902);国家 863 计划(2006AA05Z250)