网格划分基础：网格类型与质量指标

各位同学，今天咱们聊聊网格划分的基础。说实话，网格划分这事儿，我做了十几年有限元分析，到现在也不敢说完全吃透了。但基础的东西，必须得讲清楚。

网格是什么？说白了，就是把连续的结构离散成一个个小单元。你想想看，我们没法直接求解一个连续体的偏微分方程，但我们可以把它切成小块，每个小块用简单的形函数去近似。这就是有限元法的核心思想。

网格类型：二维与三维的四种基本形态

先说说二维的。二维网格就两种：三角形和四边形。

三角形网格，我习惯叫它Tria。这东西适应性特别强。你碰到复杂的几何边界，比如圆孔、尖角、过渡区域，三角形网格能很好地贴合。我在做汽车发动机缸盖的网格时，那些冷却水道的复杂形状，全靠三角形网格来搞定。

但三角形也有缺点。它的刚度偏硬，说白了就是比实际结构要刚一些。为什么？因为常应变三角形单元只有三个节点，位移模式是线性的，没法描述弯曲变形中的高阶位移。所以，同样的网格密度，三角形算出来的位移往往偏小。

四边形网格，也就是Quad。这个我特别喜欢用。为什么？精度高啊。四边形单元有四个节点，位移模式是双线性的，能更好地模拟弯曲变形。我在做桥梁结构分析时，能用四边形的地方绝不用三角形。

但四边形有个毛病——适应性差。遇到复杂边界，你得花大量时间去切分几何，才能生成全四边形网格。有时候为了一个圆角，我得折腾半天。

三维的情况类似。四面体（Tet）相当于三维的三角形，六面体（Hex）相当于三维的四边形。

四面体网格生成快，自动化程度高。现在很多商业软件，你点一下按钮，四面体网格就出来了。但精度嘛，说实话，不太行。我记得有一次做涡轮叶片的模态分析，用四面体网格算出来的前五阶频率，跟实验值差了将近15%。后来换成六面体，误差控制在3%以内。

六面体网格，这才是精度之王。但代价是什么？建模时间。一个复杂的结构件，生成六面体网格可能需要几天甚至几周。值不值得？看你分析的目的。如果是做趋势分析，四面体够了。如果是做精确的应力评估，尤其是疲劳寿命预测，我建议你咬咬牙，上六面体。

网格质量指标：三个你必须掌握的数字

网格画完了，怎么判断好不好？不能光靠眼睛看。我给大家讲三个最常用的质量指标。

1. 长宽比（Aspect Ratio）

这个最简单。对于四边形或六面体，长宽比就是最长边和最短边的比值。理想值是1，也就是正方形或正立方体。

长宽比太大意味着什么？单元被拉得很长。这种单元在长边方向的刚度会被严重低估。我在做薄壁结构时经常遇到这个问题——壁厚方向只有一层单元，长宽比能到100以上。结果算出来的局部应力完全不对。

一般来说，长宽比控制在5以内是比较安全的。对于重点区域，我建议控制在3以内。

2. 偏斜度（Skewness）

偏斜度衡量的是单元偏离理想形状的程度。对于三角形，理想形状是等边三角形；对于四边形，理想形状是正方形。

偏斜度的计算方式有很多种，最常用的是基于单元内角的。比如四边形，理想内角是90度。如果某个内角变成了120度，那偏斜度就大了。

偏斜度大的单元会带来什么问题？数值刚度矩阵的条件数会变差，求解器迭代次数增加，甚至不收敛。我曾经遇到过一个模型，怎么算都不收敛，折腾了两天。最后发现是某个区域的四面体偏斜度达到了0.95。重新划分后，一次就收敛了。

经验值：偏斜度小于0.5是优秀，0.5-0.8是可接受，大于0.8就要小心了。

3. 正交性（Orthogonality）

正交性衡量的是单元边与边之间的垂直程度。对于四边形，理想情况下相邻边夹角是90度。正交性越接近1越好。

正交性差的单元，会导致虚假的剪切锁定。什么意思？就是单元在弯曲时，本该产生的剪切变形被抑制了，结果算出来的位移偏小，应力偏大。这在薄板弯曲问题中特别明显。

我建议正交性不要低于0.7。如果低于0.5，这个单元基本就是废的。

质量指标	理想值	可接受范围	警告值
长宽比	1	1-5	>10
偏斜度	0	0-0.5	>0.8
正交性	1	0.7-1	<0.5

知识体系总览

下面这张图，我把本章的核心逻辑梳理了一下。你可以把它当作一个快速参考。

嗯，这张图把今天的内容串起来了。左边是网格类型，右边是质量指标。你选择什么类型的网格，就要用对应的质量指标去检查。两者缺一不可。

最后说一句，网格划分没有绝对的对错。同一个结构，不同的人会画出不同的网格。关键是你要理解每种网格的优缺点，以及质量指标背后的物理意义。这样你才能做出合理的判断。

好了，今天就到这儿。下一节我们讲网格密度怎么控制，以及收敛性分析怎么做。到时候我会拿一个实际案例来演示。

---

公众号:蓝海资料掘金营，微信deep3321