网格基础概念：网格类型、单元形状与网格质量指标
在CAE仿真分析中，网格划分是连接几何模型与数值求解器之间的关键桥梁。网格的本质是将连续的空间离散化为有限数量的单元，在每个单元上近似求解偏微分方程。网格的类型、单元形状以及质量指标，直接决定了求解精度、计算效率乃至求解是否能够收敛。一个看似完美的几何模型，如果网格处理不当，轻则结果失真，重则计算发散。因此，深入理解网格基础概念，是每一位仿真工程师必须夯实的基本功。
核心原理：三类网格类型的本质差异
网格类型的选择是网格划分策略的第一步，它决定了后续所有工作的走向。从拓扑结构上，网格可分为结构化网格、非结构化网格和多面体网格三大类，每一类都有其独特的生成逻辑和适用场景。
结构化网格（Structured Mesh）最显著的特征是网格节点具有规则的拓扑连接关系，内部节点拥有相同的相邻单元数量。在二维中表现为四边形网格，三维中表现为六面体网格。其生成通常依赖映射法或扫掠法，要求几何模型具有较好的规则性。结构化网格的最大优势在于计算效率高、内存占用少，且网格方向与流动或应力梯度方向容易对齐，因此数值耗散小、精度高。然而，对于复杂几何（如带有大量倒角、圆孔、不规则曲面的零件），生成纯结构化网格极其困难，甚至不可能实现。
非结构化网格（Unstructured Mesh）则没有规则的拓扑连接限制，节点可以自由分布，单元形状以三角形（2D）和四面体（3D）为主。它最大的优点是几何适应性极强，几乎可以自动离散任意复杂的CAD模型，无需大量人工切分。现代前处理软件（如ANSYS Meshing、Fluent Meshing、HyperMesh）的自动网格划分功能，主要就是针对非结构化网格。其代价是单元数量通常较多，计算精度和效率略低于理想的结构化网格，尤其在边界层内需要加密处理。
多面体网格（Polyhedral Mesh）是近年来在CFD领域迅速普及的一种类型，其单元具有多个面（通常为10-15个面）。多面体网格通常由四面体网格转换而来，每个四面体单元的中心点与相邻单元共享面连接，形成多面体单元。这种网格的优势在于：相邻单元数量多，梯度计算更准确；对网格扭曲不敏感；在相同精度下，单元总数可比四面体减少约3-5倍，收敛速度更快。在Fluent等求解器中，多面体网格已成为处理复杂内流场和外流场的首选方案之一。
单元形状：从一维到三维的几何构成
单元形状是网格的基本组成元素，不同维度下的单元形状直接影响插值函数的阶次和求解精度。在三维实体网格中，最常见的四种单元形状为四面体、六面体、棱柱和金字塔。
| 单元类型 | 节点数（线性/二次） | 典型应用场景 | 优点 | 缺点 |
|---|---|---|---|---|
| 四面体（Tet） | 4 / 10 | 复杂几何自动网格、结构分析、CFD | 几何适应性强，自动生成 | 精度较低，单元数量多，计算量大 |
| 六面体（Hex） | 8 / 20 | 规则几何、扫掠网格、边界层 | 精度高，计算效率高，各向异性控制好 | 几何适应性差，需人工切分 |
| 棱柱（Wedge/Prism） | 6 / 15 | 边界层网格、薄壁结构 | 适合边界层粘性流动，与六面体兼容 | 几何适应性有限 |
| 金字塔（Pyramid） | 5 / 13 | 四面体与六面体之间的过渡 | 实现不同拓扑单元的连接 | 质量不易控制，尽量避免大量使用 |
在结构分析中，六面体单元（尤其是20节点六面体）由于具有二次位移插值，能够更精确地捕捉应力集中和弯曲变形。而在碰撞、冲击等大变形问题中，四面体单元（特别是10节点修正四面体）由于不易发生体积锁死，反而表现更优。在CFD中，棱柱层网格是模拟壁面边界层的标准做法，它能在近壁区生成高长宽比的单元，以解析粘性底层和过渡层。
判断方法：网格质量指标及其工程阈值
网格质量是衡量单元形状偏离理想状态程度的量化指标。低质量的网格会导致刚度矩阵病态、求解精度下降甚至发散。以下是工程中必须重点监控的几项核心指标。
- 偏斜度（Skewness）：定义为单元最优形状与当前形状的差异程度。对于三角形/四面体，理想值为0（等边三角形/正四面体），最大值为1（完全退化）。一般要求偏斜度小于0.85，对于边界层网格可放宽至0.9，但超过0.95的单元必须修复。
- 正交质量（Orthogonal Quality）：主要针对CFD网格，衡量单元各面法向矢量与从单元中心到相邻单元中心矢量之间的夹角。理想值为1（完全正交），最低为0。通常要求最小正交质量大于0.1，平均值大于0.5。
- 长宽比（Aspect Ratio）：单元最长边与最短边的比值。理想值为1。对于各向同性流动，长宽比应尽量接近1；但在边界层内，为了解析法向梯度，长宽比可达10:1甚至100:1，这是允许的。然而在结构分析中，过大的长宽比会显著降低精度，一般建议不超过5:1。
- 翘曲度（Warpage）：仅适用于四边形/六面体单元，衡量单元面偏离平面的程度。理想值为0，一般要求小于5度，最大不超过10度。
- 最小/最大内角：对于三角形单元，内角应介于20°~120°之间；对于四边形单元，内角应介于45°~135°之间。超出此范围会导致插值函数畸变。
案例拆解：某汽车水泵CFD网格划分策略
以某汽车冷却水泵的CFD仿真为例，该泵体几何包含螺旋形蜗壳、不规则进出水口和复杂的叶轮区域。初始方案采用全四面体网格，单元数量约1200万，在Fluent中计算时，残差在1e-3左右停滞，且叶轮出口处的压力脉动与实验值偏差达15%。
分析与调整过程：
- 问题定位：检查网格质量发现，蜗壳与叶轮交界处的四面体单元偏斜度高达0.92，正交质量仅为0.08，导致数值耗散严重。
- 策略调整：将叶轮区域划分为旋转域，采用扫掠方法生成六面体网格；蜗壳区域保持四面体，但在壁面处生成5层棱柱边界层网格（第一层高度0.02mm，增长比1.2）；在旋转域与静止域的交界面处，使用金字塔单元进行过渡。
- 结果对比：调整后网格总数降至850万，但最小正交质量提升至0.21，平均偏斜度降至0.24。计算收敛速度提升40%，叶轮出口压力脉动与实验值偏差缩小至4%以内。
复盘经验：网格划分的核心不是“单元数量少”，而是“关键区域质量高”。对于旋转机械，叶轮和蜗壳必须采用不同的网格策略，并通过金字塔单元实现拓扑过渡。同时，边界层网格的层数和第一层高度必须根据湍流模型（如k-ε或SST k-ω）的y+要求来确定，这直接决定了近壁面流动的解析精度。
常见误区与总结复盘
在网格划分的实践中，工程师常陷入以下误区，需要特别警惕。
- 误区一：网格越密越好。加密网格确实能提高精度，但存在“边际递减效应”。当网格密度达到一定程度后，精度提升微乎其微，而计算时间呈指数增长。正确的做法是进行网格无关性验证，找到精度与计算资源的平衡点。
- 误区二：忽略单元阶次的影响。线性单元（4节点四面体、8节点六面体）在弯曲问题中容易产生剪切锁死，而二次单元（10节点四面体、20节点六面体）精度更高但计算更慢。对于结构分析，推荐使用二次单元；对于CFD，通常使用线性单元配合多面体或棱柱层。
- 误区三：只检查质量指标，不检查单元法向。在CFD中，单元法向的一致性至关重要。法向不一致会导致通量计算错误。前处理软件中应检查并统一所有单元的法向指向计算域外部。
- 误区四：忽视网格密度分布。在应力集中区域（如圆角、缺口）或流动梯度大的区域（如边界层、尾迹区），必须局部加密。使用自适应网格细化（AMR）或基于曲率的网格尺寸控制，可以自动识别这些区域。
总结：网格划分是一项“权衡的艺术”。工程师需要根据物理问题特性、几何复杂度、计算资源三方面约束，综合选择网格类型、单元形状和质量目标。建议在项目启动初期，就制定网格划分策略文档，明确各区域的目标单元类型、尺寸和最小质量阈值，并在求解前执行严格的质量检查清单。只有将网格基础概念内化为工程直觉，才能在面对复杂模型时游刃有余。