网格划分.docx
- 文档编号:5144654
- 上传时间:2022-12-13
- 格式:DOCX
- 页数:83
- 大小:784.49KB
网格划分.docx
《网格划分.docx》由会员分享,可在线阅读,更多相关《网格划分.docx(83页珍藏版)》请在冰豆网上搜索。
网格划分
第一章模型生成概述
1.1什么是模型生成?
有限元分析的最终目的是要还原一个实际工程系统的数学行为特征,换句话说分析必须是针对一个物理原型准确的数学模型。
广义上讲,模型包括所有的节点、单元、材料属性、实常数、边界条件,以及其它用来表现这个物理系统的特征。
在ANSYS术语中,模型生成一般狭义地指用节点和单元表示空间体域及实际系统连接的生成过程。
因此,在这里讨论的模型生成指模型的节点和单元的几何造型。
ANSYS程序为用户提供了下列生成模型的方法:
·在ANSYS中创建实体模型
·利用直接生成方法
·输入在计算机辅助设计(CAD)系统创建的模型。
1.2ANSYS中建模的典型步骤
通常的建模过程应该遵循以下要点:
·开始确定分析方案。
在开始进入ANSYS之前,首先确定分析目标,决定模型采取什么样的基本形式,选择合适的单元类型,并考虑如何能建立适当的网格密度。
·进入前处理(PREP7)开始建立模型。
多数情况下,将利用实体建模创建模型。
·建立工作平面。
·利用几何元素和布尔运算操作生成基本的几何形状。
·激活适当的坐标系。
·用自底向上方法生成其它实体,即先定义关键点,然后再生成线、面和体。
·用布尔运算或编号控制将各个独立的实体模型域适当的连接在一起。
·生成单元属性表(单元类型、实常数、材料属性和单元坐标系)。
·设置单元属性指针。
·设置网格划分控制以建立想要的网格密度,这个步骤并不总是必要的,因为进入了ANSYS程序有缺省单元尺寸设置存在(参见§7)。
(若需要程序自动细化网格,此时应退出前处理(PREP7),激活自适应网格划分。
)
·通过对实体模型划分网格来生成节点和单元。
·在生成节点和单元之后,再定义面对面的接触单元,自由度耦合及约束方程等。
·把模型数据存为Jobname.DB
·退出前处理。
1.2.1实体建模和直接生成的比较。
可以用两种方法来生成模型:
实体建模和直接生成。
对于实体建模,需要描述模型的几何边界,建立对单元大小及形状的控制,然后令ANSYS程序自动生成所有的节点和单元。
与之对比,用直接生成方法,在定义ANSYS实体模型之前,必须确定每个节点的位置,及每个单元的大小、形状和连接。
尽管有些数据自动生成是可能的,直接生成方法基本上是依次传递的,这种方法要求在建立有限元网格时记录所有的节点号。
这种详细的记录对于大模型来说是乏味的,并很可能出错。
实体建模一般比直接生成方法更加有效和通用,是一般建模的首选方法。
尽管实体建模有诸多优点,有时会碰到直接生成更方便的情形。
用户可以在直接生成与实体建模间方便地来回转换,对模型的不同部分采取适当的不同建模技术。
有关实体建模与直接生成的讨论可分别参考§5和§9。
为便于用户在给定条件下判别使用哪种方法更合适,现将这两种方法的相对优缺点总结如下。
1.2.1.1实体建模
实体建模的优点:
·对于庞大或复杂的模型,特别是对三维实体模型更合适。
·相对而言需处理的数据少一些。
·容许对节点和单元不能进行的几何操作(如拖拉和旋转)。
·支持使用面和体体素(如多边形面和圆柱体)及布尔运算(相交、相减等)以顺序建立模型。
·便于使用ANSYS程序的优化设计功能。
·是自适应网格划分所需的。
·为便于施加载荷之后进行局部网格细化所要求的。
(实体模型加载也需要如此)。
·便于几何上的改进。
·便于改变单元类型,不受分析模型的限制。
实体建模的缺点:
·有时需要大量的CPU处理时间
·对小型、简单的模型有时很繁琐,比直接生成需要更多的数据。
·在某些条件下可能会失败(程序不能生成有限元网格)
1.2.1.2直接生成
直接生成的优点
·对小型或简单模型的生成较方便。
·使用户对几何形状及每个节点和单元的编号有完全的控制。
直接生成的缺点
·除最简单的模型外往往比较耗时,大量需要处理的数据可能令人难以忍受。
·不能用于自适应网格划分。
·使优化设计变得不方便。
·改进网格划分十分困难(诸如面网格细化,SmartSizing等工具均不能使用)。
·可能是十分乏味的,需要用户留意网格划分的每一个细节;更容易出错。
1.3从CAD系统中输入实体模型。
代替在ANSYS中建模,可在用户擅长的CAD系统中建模,存成IGES文件格式或其他ANSYS接口产品之一的文件格式,并把它输入ANSYS中进行分析。
利用CAD软件包建模有如下优点:
·避免了重复对现有CAD模型的劳动而生成待分析的实体模型。
·工程师可利用熟悉的工具去建模。
但是,从CAD系统中输入的模型如果不适于网格划分则需要大量的修补工作。
关于从IGES文件输入实体模型的更多信息,参见§6。
对于从其他类型文件输入实体模型,参见《ANSYSConnectionUsersGuide》。
第二章规划分析方案
2.1规划的重要性
当开始建模时,用户将(有意地或无意地)作许多决定以确定如何来对物理系统进行数值模拟;分析的目标是什么?
模型是全部或仅是物理系统的部分?
模型将包含多少细节?
选用什么样的单元?
有限元网格用多大的密度?
总之,你将对要回答的问题的计算费用(CPU时间等)及结果的精度进行平衡考虑。
你在规划阶段作出的这些决定将大体上控制你分析的成功与否。
2.2确定分析目标
确定分析目标的工作与ANSYS程序的功能无关,完全取决于用户的知识、经验及职业技能,只有用户才能确定自己的分析目标,开始时建立的目标将影响用户生成模型时的其它选择。
2.3选择模型类型(二维、三维等)
有限元模型可分为二维和三维两种。
可以由点单元、线单元、面单元或实体单元组成,当然,也可以将不同类型的单元混合使用(注意要保证自由度的相容性)。
例如,带筋的薄壳结构可用三维壳单元离散蒙皮,用三维梁单元来离散蒙皮下的筋。
对模型的尺寸和单元类型的选择也就决定生成模型的方法。
线模型代表二维和三维梁或管结构,及三维轴对称壳结构的二维模型。
实体建模通常不便于生成线模型,而通常由直接生成方法创建。
二维实体模型在薄平板结构(平面应力),等截面的“无限长”结构(平面应变)或轴对称实体结构。
尽管许多二维分析模型用直接生成方法并不困难,但通常用实体建模更容易。
三维壳模型用于描述三维空间中的薄壁结构,尽管某些三维壳模型用直接生成方法创建并不困难,但用实体建模方法通常会更容易。
三维实体分析模型用于描述三维空间中截面积不等,也不是轴对称的厚结构。
用直接生成的方法建立三维实体模型较复杂,实体建模会使其变得容易些。
2.4线性和高次单元的选择
ANSYS程序的单元库包括两种基本类型的面和体单元:
线性单元(有或无特殊形状的)和二次单元。
这些基本单元类型如图2-1所示,下面来探讨这两种基本类型单元的选择。
图2-1面和体类型。
(a)线性等参元
(b)特殊形状的线性等参元
(c)二次单元
2.4.1线性单元(无中间节点)
对结构分析,带有附加形函数的角点单元会在合理的计算时间内得到准确的结果。
当使用这些单元时,要注意防止在关键区域的退化形式。
即避免在结果梯度很大或其它关注的区域使用二维三角形单元和楔形或四面体形的三维线单元。
还应避免使用过于扭曲的线性单元,对于非线性结构分析,如果使用线性单元细致地而不是用二次单元相对粗糙的进行网格划分,那么将以很少的花费获得很好的精度。
图2─2网格的比较
(a)线性单元和(b)二次单元的例子如图2-2。
当对弯曲壳体建模时,必须选用弯曲的(二次的)或平面(线性)的壳单元,每种选择都有其优缺点,对于多数的实际情况,主要问题利用平面单元以很少的计算时间,即可获得很高精度的结果。
但是,必须保证使用足够多的平面单元来创建曲面。
明显地,单元越小,准确性越好。
推荐三维平面壳单元延伸不要超过15度的弧,圆锥壳(轴对称线)单元应限制在10度的弧以内(或离Y轴5度)。
对多数非结构分析(热、电磁等),线性单元几乎与高次单元有同样好的结果,而且求解费用较低。
退化单元(三角形和四面体)通常在非结构分析中产生准确结果。
2.4.2二次单元(带中间节点)
对于用退化的单元形式进行的结构分析(即二维三角形单元和楔形或三维四面体单元),二次单元通常会以比线性单元的求解费用更低且产生良好的结果。
可是,为正确地使用这些单元,需要注意它们的特殊的性质:
·对于分布载荷和面压力不象线性单元按一般意义上分配到单元节点上(见图2-3),单元的中间节点对反力也表现出相同的非直观的解释。
·三维带中间节点的热流单元在承受对流载荷时按固定模式分配热流,在中间节点沿一个方向流动而在角点又沿另外方向的流动。
·对于结构单元,中点节点的温度如果在两相邻角点温度范围之外则要重新定义为这两角点的平均温度。
·由于中间节点的质量也大于角节点的质点,所以通常将中间节点选为主自由度(对于减缩自由度分析)。
图2─3节点分配的平衡
(a)二维单元
(b)三维单元
(c)三维三角形单元
·由于质量分配不均匀,在动力分析中感兴趣的波传波技术不推荐使用带中间节点的单元。
·不要在有中间节点(CONTAC12,COMBIN40,CONTAC48,CONTAC49,andCONTAC52)的边定义节点为基础的接触单元,也不要将间隙单元与带中间节点的边连接。
类似地,对热问题,不要应用辐射连接或非线性对流表面到带有中间节点的边。
节点为基础的接触要同有中间节点的表面接触,中间节点应该去掉。
对面对面接触单元不用担心(TARGE169,TARGE170,CONTA171,CONTA172,CONTA173,andCONTA174)。
划分实体模型时提供了一些方法忽略一些中间节点。
·当约束一个单元(或表面)的边缘自由度,包括中间节点在内的边缘上所有的节点都要约束。
图2—4在间隙或接触表面处避免中间节点
·单元的角点只能与单元的角点相连,而不能与相邻单元的中间节点相连。
相邻的单元应该有相连(或共同的)中间节点
图2─5避免单元间中间节点与角点相连
·对于有中间节点的单元,通常希望每一个这样的中间节点在相应角点之间连线的中点位置,可是,有时却希望出现在其它地方:
─节点沿着弯曲的几何边界通常可产生更准确的分析结果─所有的ANSYS网格划分器缺省地将它们放在那里。
─有的内边界也不得不弯曲以防止单元倒置或过于扭曲,ANSYS网格划分器有时也产生这种弯曲。
─用带有故意将中间节点偏离中心四分之一点可以模拟裂纹尖端的奇异性,利用ANSYS的KSCON命令可以产生这种特殊的面网格。
(MainMenu>Preprocessor>-Meshing-SizeCntrls>-ConcentratKPs-Create)
·中间节点位置可由下面描述的单元形状测试进行检查(对于控制单元形状的检查信息可参见本手册的§7章)。
──除三节点三角形和四节点四面体外的所有实体和壳单元都要进行实三维空间与单元本身的自然坐标空间一致性映射的测试。
雅可比比值过大表明单元过于扭曲,可能是由中间节点的位置设置不当引起的。
关于雅可比比值测试的细节,参见《ANSYS,Inc.TheoryReference》中的单元形状测试部分。
·如果不给中间节点指定位置,程序会自动按线性笛卡尔坐标插值将中间节点放在两角点的中间,按此法放置的节点的节点坐标系旋转角度也是按线性插值得到。
·在相连单元的公共边应有相同的节点数,当混合单元类型时有必要从一个单元去除中间节点。
例如,图2-6中的8节点单元与一个4节点单元相连时应把8节点单元的N节点去掉(或在生成这个单元时给它一个0节点号)。
图2─6避免单元相交时中点节点不匹
注意:
程序在下列情况会自动地将线性和二次单元共同一侧的中间节点去掉:
一个面(或体)用线性单元划分网格〔AMESH、VMESH、FVMESH〕,然后相邻面(或体)用二次单元划分网格。
如果网格划分的次序颠倒了,中间点节则不能去掉(先分二次单元,后分线性单元)。
·去掉了中间节点意味着边缘仍保持直的,相应地导致刚度增加,建议只在过渡区域使用去掉中间节点的单元,而不在增加了形函数的简化线单元处使用。
如果需要,那么在产生单元之后可利用下列命令增加或去掉中间节点:
命令:
EMID
GUI:
MainMenu>Preprocessor>Move/Modify>AddMidNodes
MainMenu>Preprocessor>Move/Modify>RemoveMidNd
命令:
EMODIF
GUI:
MainMenu>Preprocessor>Move/Modify>ModifyNodes
·二次单元并不比线性单元的积分点多。
因此,在非线性分析中优先使用线性单元。
·诸如平面PLANE82和SHELL93的高阶四边形单元的一种网格可由于零变形能而产生奇异。
·对后处程序只用截面的角点和隐藏线显示,类似地,节点应力结果的输出和后处理只能对角节点进行。
·在图形显示时,曲边形的中间节点单元显示为直线段(除非使用Powergaphics)模型因此看起来比实际的要粗糙些。
2.5不同单元连接的限制
在连接有不同自由度(DOFs)的单元时必须小心,因为在界面处可能会发生不协调的情况,当单元彼此不协调时,求解时会在不同单元之间传递不适当的力或力矩。
为保证协调,两个单元必须有相同的自由度。
例如,它们必须有相同数目和类型的位移自由度及相同数目和类型的旋转自由度,而且,自由度必须是耦合的,即它们必须连续地穿过界面处单元的边界。
考查三个应用了不协调单元的例子:
·单元有不同个数的自由度是不协调的,SHELL63和BEAM4单元每个节点有三个平动和三个转动自由度。
Solid45单元每个节点有三个平动自由度,但缺少转动自由度。
如果Solid45单元与SHELL63或BEAM4单元相连,相应平动自由度的节点力会传到实体块单元上。
但是,相应SHELL63和BEAM4单元的转动自由度的节点位移则不会传递给SOLID45单元。
·单元有相同个数的自由度也不总是协调的。
BEAM3(二维弹性梁)单元和SHELL41(薄膜壳)单元每个节点都各有三个自由度。
可是,壳单元有三个平动自由度(UX、UY和UZ),而梁单元只有两个平动自由度(UX和UY),因为只有UZ的结果能反映壳单元的刚度,而且壳单元没有转动自由度(ROTZ)而梁单元有。
与梁单元旋转自由度对应的节点位移将不会传递给壳单元,界面表现为梁象是被钉住了。
·三维梁单元与三维壳单元每个节点都有6个自由度,可是,壳单元的ROTZ自由度是与平面内旋转刚度相联系的,这是一个虚构的刚度;即它不是数学计算的真实刚度,因此壳的ROTZ自由度不是真实的自由度。
(例外是当SHELL43或SHELL63单元(两者都有KEYOPT(3)=2)的Allman旋转刚度被激活时),因此对三维梁单元仅有一个节点与三维壳单元相连导致梁单元的旋转自由度与壳单元的ROTZ自由度对应是不协调的,不应该将梁单元与壳单元按此方式连接起来。
类似的不协调也可在不同数目和(或)类型自由度的单元之间存在。
这些问题并不会使分析无效,但至少应该注意到两个不同类型单元交界的条件。
2.6找到利用对称性的办法
许多物体都有某种对称,重复对称(如一个长管道上平均分布的散热片,反对称(如塑料容器的模子),或轴对称(如灯炮)。
当一个物体在所有方面都是对称的(几何、载荷、约束和材料属性),可以利用用述事实减少模型的大小和范围。
图2-7对称的例子
(A)重复对称(散热等)(B)镜像对称(模铸容器)(C)轴对称(灯炮)
2.6.1轴对称结构的说明
任何表现在为对一中心轴几何对称的结构(如旋转壳或实体)是轴对称。
例如直管、圆锥、圆盘、圆盖等。
三维轴对称结构的模型可用二维形式等效。
可以想象得到二维轴对称分析较相应的三维分析更精确。
由定义可知,一个完全轴对称模型只能施加轴对称载荷。
但在很多情况下,轴对称结构实际受的并不是轴对称载荷,必须利用一种特殊的轴对称谐波单元,来建立轴对称结构二维模型的非轴对称载荷。
详见《ANSYSElementsReference》中的§2.10一节。
2.6.1.1对轴对称结构的一些特殊要求
对轴对称结构模型的特殊要求有:
·对称轴必须与总体笛卡尔坐标的Y轴重合。
·不允许出现负的X方向节点坐标。
·总体笛卡尔坐标Y方向代表轴向,总体笛卡尔坐标X方向代表径向,总体笛卡尔坐标Z方向代表周向。
·模型应当用适当的单元类型组合在一起:
——对轴对称模型,可用二维实体且KEYOPT(3)=1,和(或)轴对称壳体。
而且,可用不同的连接,接触,组合及表面单元并可将轴对称实体和壳单元容在一个模型里。
(除非是轴对称实体或壳单元否则程序不会识别其它的单元)如果《ANSYSElementsReference》没有对一个特殊单元类型讨论其在轴对称结构中的应用,那么不要使用此种单元类型进行轴对称分析。
——对轴对称谐波模型,只能使用轴对称谐波单元。
·SHELL51和SHELL61单元不应位于总体Y轴上。
·对包含二维实体单元的模型剪切影响是重要的,在厚度方向上至少要使用二个单元。
2.6.1.2进一步的须知和限制。
如果结构沿对称轴包含有孔,不要忘记在Y轴和二维轴对称模型间留适当的距离(见图X方向的偏移表示一个轴对称孔。
)对轴对称载荷的讨论参见《ANSYSBasicAnalysisGuide》的§2。
图2-8X方向的偏移表示一个轴对称孔
2.7决定包含多少细节
一些小的细节对分析来说不重要,不必在模型中体现,因为它只会使你的模型过于复杂。
可是对有些结构,小的细节如倒角或孔可能是最大应力位置之所在,可能非常重要,取决于用户的分析目的。
必须对结构的预期行为有足够的理解以对模型应包含多少细节做出适当的决定。
有些情况下,仅有一点微不足道的细节破坏了结构的对称。
那么,可以忽略这些细节(或相反的将它们视为对称的),以利于用更小的对称模型,必须权衡模型简化带来的好处与精度降低的代价来确定是否对一个非(拟)对称结构故意忽略其非对称细节。
2.8确定合适的网格密度
有限元分析中经常碰到的问题是单元网格应划分得如何细致才能获得合理的好结果?
不幸的是,还没有人能给出确定的答案;你必须自己解决这个问题,关于这个问题的解决可求助于以下一些技术:
·利用自适应网格划分产生可满足能量误差估计准则的网格(此技术只适用于线性结构静力或稳态热问题,对什么样的误差水平是可接受的判断依据于你的分析要求)。
自适应网格划分需要实体建模。
·与先前独立得出的实验分析结果或已知解析解进行对比。
对已知和算得结果偏差过大的地方进行网格细化。
(对所有的由四面体组成的面或体网格可用NREFINE、EREFINE、KREFINE、LREFINE和AREFINE命令(MainMenu>Proprocessor>-Meshing-ModifyMesh>-RefineAt-entitytype)进行局部网格细化)
·执行一个你认为是合理的网格划分的初始分析,再在危险区域利用两倍多的网格重新分析并比较两者的结果。
如果这两者给出的结果几乎相同,则网格是足够的。
如果产生了显著不同的结果,应该继续细化你网格直到随后的划分获得了近似相等的结果。
·如果细化网格测试显示只有模型的一部分需要更细的网格,可以对模型使用子模型以放大危险区域。
网格划分密度很重要,如果网格过于粗糙,那么结果可能包含严重的错误,如果网格过于细致,将花费过多的计算时间,浪费计算机资源,而且模型可能过大以致于不能在你的计算机系统上运行,为避免这类问题的出现,在生成模型前应当考虑网格密度问题。
第三章坐标系
3.1坐标系的类型
ANSYS程序提供了多种坐标系供用户选取。
·总体和局部坐标系用来定位几何形状参数(节点、关键点等)的空间位置。
·显示坐标系。
用于几何形状参数的列表和显示。
·节点坐标系。
定义每个节点的自由度方向和节点结果数据的方向。
·单元坐标系。
确定材料特性主轴和单元结果数据的方向。
·结果坐标系。
用来列表、显示或在通用后处理(POST1)操作中将节点或单元结果转换到一个特定的坐标系中。
工作平面与本章的坐标系分开讨论,以在建模中确定几何体素,参见§4中关于工作平面的详细信息。
3.2总体和局部坐标系
总体和局部坐标系用来定位几何体。
缺省地,当定义一个节点或关键点时,其坐标系为总体笛卡尔坐标系。
可是对有些模型,定义为不是总体笛卡尔坐标系的另外坐标系可能更方便。
ANSYS程序允许用任意预定义的三种(总体)坐标系的任意一种来输入几何数据,或在任何用户定义的(局部)坐标系中进行此项工作。
3.2.1总体坐标系
总体坐标系统被认为是一个绝对的参考系。
ANSYS程序提供了前面定义的三种总体坐标系:
笛卡尔坐标、柱坐标和球坐标系。
所有这三种系统都是右手系。
且由定义可知它们有共同的原点。
它们由其坐标系号来识别:
0是笛卡尔坐标,1是柱坐标,2是球坐标(见图总体坐标系)
图3-1总体坐标系
·(a)笛卡尔坐标系(X,Y,Z)0(C.S.0)
·(b)柱坐标系(R,θ,Zcomponents)1(C.S.1)
·(c)球坐标系(R,θ,φcomponents)2(C.S.2)
·(d)柱坐标系(R,θ,Ycomponents)5(C.S.5)
3.2.2局部坐标系
在许多情况下,有必要建立自己的坐标系。
其原点与总体坐标系的原点偏移一定的距离,或其方位不同于先前定义的总体坐标系(如图3-2所示用局部、节点或工作平面坐标系旋转定义的一个坐标系的例子)。
用户可定义局部坐标系,按以下方式创建:
图3-2欧拉旋转角
·按总体笛卡尔坐标定义局部坐标系。
命令:
LOCAL
GUI:
Utility Menu>WorkPlane>LocalCoordinateSystems>CreateLocalCS>AtSpecifiedLoc
·通过已有节点定义局部坐标系。
命令:
CS
GUI:
Utility Menu>WorkPlane>LocalCoordinateSystems>CreateLocalCS>By3Nodes
·通过已有关键点定义局部坐标系。
命令:
CSKP
GUI:
Utility Menu>WorkPlane>LocalCoordinateSystems>CreateLocalCS>By3Keypoints
·在当前定义的工作平面的原点为中心定义局部坐标系。
命令:
CSWPLA
GUI:
Utility Menu>WorkPlane>LocalCoordinateSystems>CreateLocalCS>AtWPOrigin
·通过激活的坐标系由CLOCAL命令定义局部坐标系。
(见§3.2.3一节)(没有与CLOCAL相应的GUI)
当用户定义了一个局部坐标系后,它就会被激活。
当创建了局部坐标系后,分配给它一个坐标系号(必须是11或更大),可以在ANSYS进程中的任何阶段建立(或删除)局部坐标系。
若要删除一个局部坐标,利用下列方法:
命令:
CSDELE
GUI:
UtilityMenu>WorkPlane>LocalCoordinateSystems>DeleteLocalCS
要查看所有的总体和局部坐标系使用下列方法:
命令:
CSLIST
GUI:
UtilityMenu>List>Other>LocalCoordSys
与三个预定义的总体坐标系类似,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 网格 划分