百家争鸣：CAE前处理之网格划分评价指标控制要求（附直播教程）
导读：有限元网格划分是有限元分析前处理的重要环节，有限元单元类型的选择、单元尺寸大小控制、网格划分方式都对分析结果的准确性有很大影响。不同的分析人员网格划分的结果都不一样，这也直接导致应力分析的结果不唯一。目前对于网格划分的要求和评价标准，在压力容器行业没有明确规定，与其他行业相比，认识上也有较大的差异。本文从不同行业专家的认识、ANSYS软件的评价标准以及国内外相关规范的规定几个方面，做了对比分析。对于结构有限元应力分析，最常用的是四面体网格和六面体网格。在ANSYS中的四面体推荐使用10节点高阶的Solid187单元，该单元具有一下特点：如图1给出了四面体187单元的示意图，该单元为ANSYS单元库中的10节点高阶四面体单元，也是ANSYS Workbench环境中对复杂模型优先使用的单元。图2给出了四面体Solid285单元的示意图，该单元为ANSYS单元库中的4节点低阶四面体单元。在ANSYS中高阶六面体结构单元就是Solid186单元，如图3所示，六面体单元的特点：2）对于复杂模型，需要在前处理耗费较长时间，对复杂模型的适应性差；图4给出了六面体Solid185单元的示意图，该单元为ANSYS单元库中的8节点低阶六面体单元。宋博士在其新浪微博“宋博士的博客”2015年1月3日的文章《ANSYS WORKBENCH中划分网格的几种方法》中给出网格划分的观点如下：总体来说，对于空间物体而言，我们应当尽量使用六面体网格。当对象是一个简单的规则体时，使用扫掠网格划分是合适的；当对象是多个简单的规则体组成时，使用多域扫掠网格划分是合适的；接着尽量使用六面体主导的方式，它会在外层形成六面体网格，而在心部填充四面体网格。四面体网格是最后的选择。其中如果要忽略一些小细节，如倒角，小孔等，则使用patch independent算法;如果要要考虑一些小细节，则使用patch conforming算法。至于自动网格划分，是最傻瓜化的方式，一般对于初学者适用。仿真秀平台专栏作者张老师在“ANSYS空间”发表的文章《四面体单元与六面体单元的选择建议》中给大家的使用建议：（1）如果使用自动方法不能产生六面体单元的模型，建议使用四面体；（2）通过提升四面体单元的密度，来提升计算精度，这样处理可以大大减少模型前处理的时间，来实现仿真的时效性优势，关于计算精度对比，后面会通过实例进行对比。另外，在《ANSYS有限元网格划分的基本原则》一文中写道：单元质量评价一般可采用以下几个指标：（1）单元的边长比、面积比或体积比以正三角形、正四面体、正六面体为参考基准。理想单元的边长比为1，可接受单元的边长比的范围线性单元长宽比小于3，二次单元小于10。对于同形态的单元，线性单元对边长比的敏感性较高阶单元高，非线性比线性分析更敏感。（3）疏密过渡：网格的疏密主要表现为应力梯度方向和横向过渡情况，应力集中的情况应妥善处理，而对于分析影响较小的局部特征应分析其情况，如外圆角的影响比内圆角的影响小的多。（4）节点编号排布：节点编号对于求解过程中的总体刚度矩阵的元素分布、分析耗时、内存及空间有一定的影响。合理的节点、单元编号有助于利用刚度矩阵对称、带状分布、稀疏矩阵等方法提高求解效率，同时要注意消除重复的节点和单元。仿真秀平台专栏作者“水哥”在“ANSYS结构院”发布的文章《ANSYS网格质量评定指标介绍》一文中写道：不良的单元形状会导致不准确的结果，然而到目前为止，还没有一个比较通用的标准来判定单元形状的好坏。一种单元形状在一个分析中可能会带来不正确的结果，但在另外一个分析中又可能是完全能接受的，因而单元形状的好坏以及结果的准确性完全由用户根据经验或者相关行业规范进行判定和分析。所以平时我们在采用ANSYS进行有限元分析时，一个很重要的步骤便是查看网格质量，有时即使网格划分成功了，若质量太差，结果也不一定具有可信性，故而好的网格质量是保证分析结果精确的首要条件。（1）Aspect Radio（网格纵横比）：其值越接近1，说明网格质量越好。（2）Parallel Deviation（对边偏角差）：其值越接近0，说明网格质量越好。（3） Maximum Corner Angle（单元最大内角）：三角形，越接近60度越好；四边形，越接近90度越好。（4）Jacobian Ratio（雅克比比率）：其值越接近1，说明网格质量越好。（5）Wraping Factor（翘曲系数）：其值越接近0，说明网格质量越好。（6）Element Quality（网格质量系数）：其值越接近1，说明网格质量越好。（7）Skewness（偏斜系数）：其值越接近0，说明网格质量越好。（8）Orthogonal Quality（正交质量系数）：其值越接近1，说明网格质量越好。仿真秀专栏作者梁庆海老师在“承压设备分析设计”的文章《我们-网格篇》中写道：有时想要划分出六面体，却怎么也实现不了，一部分是六面体，一部分却只能是四面体，不敢给外人展示，好尴尬。报告里面如果给出漂漂的六面体网格，那是计算结果真实准确的体现，客户也会相信你的计算结果。网格质量的好坏是计算结果能否被接受的重要因素，如果网格质量很粗糙，那么这个计算很可能被否掉，直接拉回去毙了。我们做了很多整体分析设计，网格划分质量99%都是六面体，只有一个或两个体实在是无法划分出六面体，就像钻石的纯净度一样，只有很少很少的污点，瑕不掩玉。对于有些结构很难划分出六面体，常常会出现很奇葩的问题，比如体上的线都是断的，让划分网格变得很困难，自由划分也不行，多区域划分也不行，扫略更不行了，弄了一上午都划分不出来，哦，疯掉的感觉。陆明万、徐鸿在论文《分析设计中若干重要问题的讨论》中对压力容器分析设计网格划分提出一下建议：
第一、板壳结构采用实体单元进行模拟时，厚度方向应有4层线性单元或2层二次单元。
第二、1/4圆弧一般分为6~8个单元。
第三、四边形单元的长宽比一般在1:1到1:3之间。在板壳结构中，均匀的薄膜应力和弯曲应力区可以取1:5~7，最多可以放大到1:10。
第四、三角形单元的锐角一般不小于30°，最小到15°。四边形单元的最小夹角一般不小于45°。
第五、应力梯度大的地方网格要加密。在主要的递减（或递增）区间内至少3~4个单元。
第六，疏密网格的尺寸比一般不大于2。
第七，对尚无经验的算例，需要论证网格尺寸大小的合理性。做基本网格和尺度加密1倍的网格两个算例，若两者计算结果只差在3%~5%以下，则基本网格是可用的。在北京希格玛仿真技术有限公司发布的压力容器局部开发应力分析软件NSAS：选择高阶单元时，厚度方向2层，选择低阶单元时，厚度方向4层，均采用六面体单元。ANSYS Workbench中单元质量评价指标如下，其类别和含义即为前文中所描述。Workbench中增加了单元质量显示功能，如下图所示，能看到对应质量单元的数量及分布位置。GB/T33582-2017《机械产品结构有限元力学分析通用规则》对于7.4条网格划分要求如下：对于实体单元网格，在结构厚度方向上应确保三层以上。SYT7403-2018《油气输送管道应变设计规范》对于管道有限元应力分析单元尺寸的要求主要有以下两点：（1）当采用管单元时，直线段单元划分的长度不宜大于管径。弯管段单元划分长度不宜大于0.5倍管径且划分节点数不应小于3个；（2）当采用壳单元时，应根据管道各部分变形特点进行网格划分，包括单元密度、大小和形状。划分的单元边长不宜大于0.3倍管径。BS7910-2013《Guide to methods for assessing the acceptability of flaws in metallic structures》在对管道局部腐蚀减薄分析计算程序中，给出了建模、单元及网格的要求。G.9.4 网格设计相关内容及对应中文翻译如下：The element type should be adequate for the analysis being undertaken. Three dimensional continuum elements are preferred to shell elements. Second order, twenty noded, reduced integration hexahedral elements should be used. Severe element distortions, especially at the minimum ligament of the LTA, should be avoided. If shell elements are used, they should include variable wall thickness in the formulation and allow the section to be offset to represent internal or external material loss at the LTA. A minimum of eleven material integration points should be used through the thickness.单元类型应该适合正在进行的分析。三维连续单元优于壳单元，应使用二阶、二十个节点的缩减积分六面体单元。应该避免严重的单元变形，特别是在 LTA 的最小韧带处。如果使用壳单元，它们应在公式中包括可变壁厚，并允许截面偏移以表示 LTA 处的内部或外部材料损失。整个厚度至少应使用 11 个积分点。（5层高阶单元）DNV GL RP C208-2016《Determination of structural capacity by non-linear finite element analysis methods》对有限元网格要求做了比较详细规定，主要条款对应原文和中文翻译如下：（1）In general, structural parts welded together should be meshed using a continuous mesh.一般来说，焊接在一起的的结构件应该采用连续网格进行网格划分。（2）Connections and constraints such as bonded contact or kinematic coupling etc. should not be used for welded details in areas of interest unless the accuracy on stress and strain results is documented and accounted for when evaluating the results.除非验证了应力和应变结果的准确性，并在评估结果时对其进行了说明，否则相关区域的焊接细节不应使用粘结接触或运动耦合等连接。（3）Selection of element type and formulation is strongly problem dependent. Items to consider are:Recommended practice, DNVGL-RP-C208 – Edition September 2016 Page 18单元类型和公式的选择与分析问题密切相关。要考虑的项目是：推荐做法，挪威船级社DNVGL-RP-C208 – 2016 年 9 月版 第 18 页
（4）In general higher-order elements are preferred for accurate stress estimates; elements with simple shape functions (constant or linear) will require more elements to give the same stress accuracy as higher-order elements. Constant stress elements are not recommended used in areas of interest.一般来说，精确应力估计应该首选高阶单元；如果采用低阶单元，将需要更多数量的单元才能达到高阶单元相同的应力精度。不建议在感兴趣的区域使用低阶应力单元。（5）For large-displacement and large-rotation analyses, simple element formulations give a more robust numerical model and analysis than higher-order elements.对于大位移和大旋转分析，与高阶单元相比，低阶单元提供了更稳健的数值模型和分析。（6）Some types of elements are intended as transition elements in order to make the generation of the element mesh easier and are known to perform poorly. Typically 3-noded plates/shells and 4-noded tetrahedrons are often used as transition elements. These types of elements should if possible be avoided in areas of interest.某些类型的单元用作过渡单元，以便更容易生成网格，尽管已知其性能较差。通常，三节点板/壳和四节点四面体通常用作过渡单元。如有可能，应避免在感兴趣的区域使用这些类型的单元。（7）Proper continuity should be ensured between adjacent elements if elements of different orders are used in the same model.如果在同一模型中使用了不同阶的单元，则应确保相邻单元之间的适当连续性。（8）The element mesh should be sufficiently detailed to capture the relevant failure modes:— For ductility evaluations, preferably several elements should be present in the yield zone in order to have good strain estimates.
-对于韧性评估，最好在屈服区中存在几个单元，以便进行良好的应变估计。
— For stability evaluations, sufficient number of elements and degrees of freedom to capture relevant buckling modes, typically minimum 3 to 6 elements dependent upon element type per expected half wave should be used.
-对于稳定性评估，应使用足够数量的单元和自由度来捕获相关屈曲模式，通常应根据每个预期半波的单元类型至少使用3到6个单元。
The element aspect ratio should be according to requirements for the selected element formulation in the areas of interest. Typically, an aspect ratio close to unity is required in and nearby areas with large deformations.
单元纵横比应符合感兴趣区域内选定单元算法的要求。通常，在变形较大的区域及其附近需要接近1的纵横比。
Care is required in transitioning of mesh density. Abrupt transitioning introduces errors of a numerical nature.
在网格密度过渡时需要小心，差异太大会引入数值性质的错误。
Load distribution and load type also have an influence on the mesh density. Nodes at which loads are applied need to be correctly located, and in this situation can drive the mesh design, at least locally.
- 载荷分布和载荷类型也会对网格密度产生影响。施加荷载的节点需要正确定位，在这种情况下，可以至少在局部修改网格设计。
（9）Often it will be necessary to run mesh sensitivity studies in order to verify that the results from the analyses are sufficiently accurate.通常有必要进行网格敏感性研究，以验证分析结果是否足够准确。The analyst should make sure that the element mesh is adequate for representing all relevant failure modes. In the general case mesh refinement studies may be done by checking that convergence of the results are obtained e.g. by showing that the results are reasonably stable by rerunning the analysis with half the element size. See example in [8.2]. Note that geometric sharp corners represent singularities where convergence will never be obtained.分析员应确保单元网格足以表示所有相关故障模式。在一般情况下，可以通过检查结果是否收敛来进行网格细化研究，例如，通过使用一半的单元大小重新运行分析来证明结果是合理稳定的。参见[8.2]中的示例。请注意，几何尖角表示永远无法收敛的奇点。不同的行业，不同的失效模式，对网格质量的要求不一样，具体分析过程中，应兼顾行业规范规定及软件自身的质量评价指标。不同行业对网格要求的共同点，均要求做网格敏感性分析，规范中均没有明确规定偏差范围，建议参考陆明万、徐鸿在论文《分析设计中若干重要问题的讨论》中的要求，偏差在5%以内。11月1日-30日，仿真秀平台力邀压力容器行业资深行业问题解决专家，也是仿真秀平台的认证讲师，将为理工科学子和压力容器研发技术工程师带来实用性压力容器分析设计技术讲座。11月10日（周三）20时，笔者受仿真秀平台邀请在仿真秀官网和APP公开直播“压力容器应力分析前处理那些事”一起来聊聊有限元网格划分！期待您的光临！（扫描下方二维码可以报名，支持回看哦）
赵树炳，硕士研究生，原就职于中国石油天然气管道工程有限公司，现为北京应力分析科技有限公司总经理及技术负责人，主要从事应力分析、应力测试及ANSYS软件代理销售工作，作为技术专家多次进行ANSYS结构分析技术培训，具有丰富的软件使用和工程实践经验。先后为兰州兰石重型装备股份有限公司、山东济炼石化工程有限公司、中国大唐集团科学技术研究院、河南瑞昌石化工程有限公司、许昌远东传动轴股份有限公司、中国航天华创股份有限公司、株洲联诚集团控股股份有限公司、滕州市江鑫机械有限公司、胜利油田设计院、北京热力设计公司等单位进行培训服务。曾与中国石油和化学工业培训中心合作主讲《ANSYS结构分析及管道工程应用实例高级技术培训》，并主办《2018年压力容器、压力管道应力分析学术研讨会》。研究方向：压力容器、压力管道有限元应力分析及实验应力分析技术。作者简介：赵树炳，仿真秀专栏作者，硕士研究生，原就职于中国石油天然气管道工程有限公司，现为北京应力分析科技有限公司总经理及技术负责人。 获赞 10098粉丝 21559文章 3539课程 219
