3. 网格质量指标详解

非结构网格的质量评价不能只看“有没有负体积”，而要从多个几何维度同时刻画单元形状偏离理想形态的程度。偏斜度、长宽比、雅可比、正交性和扭曲度是工程中最常用的五类指标，它们分别从角度、尺度、映射有效性、面法向一致性和三维畸变角度描述单元质量。理解每个指标的计算对象与物理含义，才能在看云图时判断“这个单元到底坏在哪”，并选择正确的修复手段。

一、偏斜度：角度偏离理想形态的程度

偏斜度（Skewness）衡量单元内角或面角偏离理想值的程度。对三角形单元，理想内角为 60°；对四边形和六面体，理想面角为 90°。常用定义之一为：

Skewness = max( (θ_max − θ_e) / (180° − θ_e), (θ_e − θ_min) / θ_e )

其中 θ_e 为理想角（三角形 60°，四边形/六面体 90°），θ_max、θ_min 为单元实际最大、最小角。该式把偏斜度归一化到 0～1：0 表示理想等角单元，接近 1 表示出现极尖或极扁的角。工程上通常要求三角形偏斜度小于 0.8，六面体小于 0.85，超过 0.95 的单元往往导致求解器矩阵条件数急剧上升。偏斜度只反映角度，不反映尺度，因此一个又长又窄但角度正常的单元可能偏斜度不高，却长宽比很差。

二、长宽比：尺度各向异性的度量

长宽比（Aspect Ratio）描述单元在不同方向上的尺度差异。对四边形或六面体，常取最长边与最短边之比，或最长特征尺寸与最短特征尺寸之比；对三角形，可用外接圆半径与内切圆半径之比（等边三角形为 2）。长宽比接近 1 表示各向同性，数值越大表示单元越“扁长”。

长宽比并非越小越好。在边界层内，为捕捉法向梯度，刻意使用长宽比 100～1000 的扁单元是合理的；但在各向同性流动区域出现长宽比大于 20 的单元，则会造成数值扩散方向失真。因此判断长宽比要结合所在区域：边界层内看增长率和法向尺寸，远场看是否接近各向同性。

三、雅可比：映射有效性的局部判据

有限元与有限体积方法把标准单元（如单位三角形、标准六面体）通过形函数映射到实际单元。雅可比矩阵 J 描述这一映射的局部伸缩与旋转，其行列式 det(J) 的符号和大小直接决定单元是否有效：det(J) > 0 表示映射保向、单元未翻转；det(J) = 0 表示单元退化（压扁成面或线）；det(J) < 0 表示单元局部翻转，是必须清除的致命错误。

工程中常报告“最小雅可比”或“归一化雅可比”：将单元内各积分点的 det(J) 除以该单元最大 det(J)，得到 0～1 的比值。六面体最小雅可比低于 0.3 通常视为质量较差，低于 0.1 或为负则必须修复。雅可比对节点位置非常敏感，一个节点被误拖到相邻单元内部，往往首先表现为局部雅可比为负。

四、正交性：面法向与连线方向的夹角

正交性（Orthogonality）主要用于评价相邻单元公共面两侧的几何关系。对内部面，取面法向 n 与两侧单元中心连线 d 的夹角 α，定义正交性角为 90° − α，或直接报告夹角偏离 90° 的度数。理想正交时面法向与中心连线平行，α = 0°；偏离越大，离散格式中非正交修正项越重要。

在非结构网格中，正交性差会引入较大的非正交修正，影响梯度计算和压力方程收敛。工程上常要求最大非正交角小于 60°～70°，超过 75° 时需检查是否由尖锐几何、面网格过渡过急或节点投影不当造成。正交性是“面—对”级别的指标，修复时通常通过局部光滑节点或调整面网格分布改善，而不是简单加密。

五、扭曲度：三维单元的翘曲与畸变

扭曲度（Warpage / Twist）针对四边形面和六面体单元，描述面内节点不共面或单元棱边相对扭转的程度。对四边形面，可将四个节点拟合为平面，取各节点到该平面的最大距离与特征边长之比作为翘曲度；对六面体，还可比较相对面的扭转角。扭曲度大意味着面不再近似平面，体积计算和面通量计算误差增大。

扭曲度与偏斜度、长宽比并不冗余：一个面可能四边等长、角度接近 90°，但四个角点不共面，此时偏斜度和长宽比都不报警，扭曲度却明显偏大。因此在六面体主导的网格中，扭曲度是必须单独检查的指标。

六、指标之间的关系与判读顺序

五个指标并非彼此独立。偏斜度关注角度，长宽比关注尺度，雅可比关注映射是否有效，正交性关注面与中心连线的方向关系，扭曲度关注三维面外畸变。实际检查时建议按“先致命、后一般”的顺序：先查最小雅可比是否为负或接近零，排除翻转和退化；再看偏斜度和正交性是否超过求解器建议阈值；最后结合区域判断长宽比和扭曲度是否合理。只盯一个指标容易漏判，例如边界层内长宽比很大但雅可比和正交性正常，属于有意设计而非缺陷。

七、一个可复现的判读示例

假设对某汽车外流场模型做质量检查，得到如下统计（示例数据，用于说明判读方法）：三角形偏斜度最大值 0.93，六面体最小归一化雅可比 0.18，最大非正交角 71°，边界层第一层长宽比 320，过渡区出现长宽比 45 的单元，扭曲度最大值 0.42。判读时：偏斜度 0.93 已接近上限，需定位尖角单元；雅可比 0.18 偏低但未翻转，属于可修复范围；非正交角 71° 接近 70° 阈值，应检查尖锐特征附近；边界层长宽比 320 合理，但过渡区 45 的长宽比偏大，说明尺寸场过渡过快；扭曲度 0.42 提示部分六面体面翘曲，需拆分或改用四面体填充。对应处理是：对偏斜度大的面做局部重划，对雅可比低的节点做光滑回退，对过渡区放宽尺寸增长率，对翘曲六面体区域改用金字塔或四面体过渡。

掌握这五类指标的计算对象和物理意义后，下一步才能针对性地选择修复策略：偏斜度和正交性多通过面网格重划与节点光滑改善，长宽比通过尺寸场和增长率调整，雅可比通过节点回退与局部重建恢复，扭曲度则往往需要改变单元类型或拆分过渡。指标不是打分工具，而是定位问题的坐标。