摘要
本文强调了ANSYS有限元分析中网格划分与质量检查的重要性。高质量的网格是计算结果准确性的基础,需根据不同模型部位进行有针对性的设置和划分。同时,减少网格节点和单元数量能提升计算效率。文章介绍了ANSYS中八个主要网格质量评估参数,并通过实例展示了网格划分的优缺点和策略。这些内容有助于用户在有限元分析中合理设置网格,确保计算结果的准确性和效率。
正文
对于网格质量的判定,目前还没有一个比较通用的标准来判定好坏。一种单元形状在一个分析中可能会带来不正确的结果,但在另外一个分析中又可能是完全能接受的,因而单元形状的好坏以及结果的准确性一部分取决于用户根据肉眼观察、网格无关性验证、经验或者相关行业规范进行判定和分析。当然,ANSYS软件中也给出了评定网格质量的八个主要参数,分别是单元质量、纵横比、平行偏差、单元最大内角、雅克比、单元翘曲因子、偏斜度、正交质量。那么,这几个参数代表什么含义以及如何通过这几个参数对网格质量进行一个辅助判断呢?
上式中,C为常数,根据单元类型取值,C的取值在Ansys Help文件中有明确值。该系数为1时,表明网格质量最好,系数为0时,表明网格质量最差,在具体做分析时,可保证平均值不低于0.7,当然这是经验之谈,具体项目可具体分析。
Brick element 91 has anaspect ratio of 1.E+20, which exceeds the error limit of1000000。
这个就是提示单元的纵横比超限了,需改善网格,纵横比默认情况下达到20软件会发出警告,超过10^6时,则会直接弹出错误。
纵横比的计算分两种情况:三角形和四边形,对于实体单元,选取每一面的中线连接成三角形,然后在依次计算,如:
也即单元形状越不接近等边三角形或者正方形,则纵横比越大,最好的值为1,即等边三角形或正方形。
对于每一对相对的边,计算单位向量的点积,然后计算余弦为点积的角度(以度为单位)。平行偏差是这两个角中较大的一个。(在上图中,两个水平单位向量的点积为1,acos(1) = 0°。两个竖直向量的点积为0.342,acos(0.342) =70°。因此,该元件的平行偏差为70°。)
系统默认情况下:
(1)当无中间节点时,发出警告限值为70度,发出错误限值为150度;
(2)当由中间节点时,发出警告限值为100度,发出错误限值为170度;
故单元的形状越接近矩形,平行偏差角越小,单元形状越好,最好的值为0。
该参数主要是检查单元的最大内角值,参数很好理解,下面为参数在不同取值情况下的单元形状:
系统默认情况下:
(1)对于三角形,发出警告限值为165°,发出错误限值为179.9°;
(2)对于无中间节点四边形,发出警告限值为155°,发出错误限值为179.9°;
(3)对于有中间节点四边形,发出警告限值为165°,发出错误限值为179.9°。
故对四边形而言,90°最好;对三角形而言,60°最好。
雅可比是对给定单元形状与理想单元形状的一种度量。如果一个元素的雅可比质量不好,则该元素可能不能很好地从单元空间映射到实际空间,从而使这种单元的计算不可靠。这里可以简单的将雅克比系数理解为有限元模拟与实际情况的可靠度指标。雅克比的计算有基于MAPDL的,基于角节点的,基于高斯积分点的,比较繁杂,具体的方法可查询AnsysHelp文件。
下图为基于角节点的雅克比单元形状,雅克比越接近于1,则形状越规则,趋于等边三角形、正方形或矩形,网格质量则越好。
大家在网格划分时候也会经常遇到单元翘曲过度的情况,四边形壳单元和六面体单元、楔形体的四边形表面、金字塔单元的表面需要计算和检测翘曲因子,翘曲因子越高,表示底层单元公式不能很好地处理的条件,或者网格生成有缺陷。
下图分别为壳单元和实体单元在不同翘曲系数情况下的单元形状,由图可直观看出:无论对于壳单元或是实体单元,翘曲因子越小时,单元形状越规则,网格质量越好,翘曲因子为0时,网格质量则最好;值越大,则翘曲度越大。
系统默认情况下:
2、实体单元警告值为0.2,错误值为0.4。
该系数是衡量网格质量的主要指标之一。偏斜度决定了一个面或单元与理想(等边或等角)的接近程度。如下图左侧的等边三角形和矩形为理想状态,偏斜度为0,表明单元形状越规则。该值越小,表示越接近理想状态,网格质量越好;相反,该值越大,则表明网格质量越不好。
下图列出了偏斜度值得范围与单元网格质量的对应关系:
偏斜度为0时,表明单元形状是等边的,即理想状态;>0~0.25是优秀的;0.25~0.5是好的;0.5~0.75是一般可接受的;0.75~0.9是较差的;0.9~1,则是差的和完全不可接受的。故而,该偏斜度值越小越好,且不要超过0.75。
该值的计算方法采用基于相邻面的三个法向矢量计算,比较复杂,想进一步了解的可参考Ansys Help帮助文件。我们只需记住,正交质量的取值范围为0-1,其中:0为最差,1为最佳,网格质量最好。
Ansys软件中除过提供了上述八大网格质量判定参数外,新版的软件还增加了最小单元边长(Minimum ElementEdge Length)、最大单元边长(Minimum Element Edge Length)、最小三角形角度(Minimum Tri Angle)、最大三角形角度(Maximum Tri Angle)、最小四边形角度(Minimum Quad Angle)、最大四边形形角度(Maximum QuadAngle)等评判参数。具体可查看Ansys Help帮助文件。
本文部分内容来源于网络