网格质量 — CAE术语解说

好问题。比如在结构分析中使用扭曲的四边形单元。单元内的位移通过形状函数插值，但当单元极度扭曲时，形状函数的导数，即应变的计算就会出错。具体来说，雅可比矩阵的行列式接近零，导致应变能量被低估。这发生在应力集中处时，会显示出远低于实际的应力值。

考虑从自然坐标系(ξ, η)到实坐标系(x, y)的映射。位移在实坐标系中的导数为：

会出现，但机制略有不同。最重要的是「数值扩散」。例如，使用迎风差分法时，如果单元沿流向极度扭曲，人工粘性会被过度评估，导致涡和剪切层比实际衰减更快。另外，在求解压力-速度耦合的SIMPLE法等方法时，网格的非正交性过强会导致压力修正方程的系数矩阵条件数变差，可能发散或收敛极其缓慢。

许多求解器在组装单元刚度矩阵或质量矩阵前的前处理阶段检查这些指标并输出警告。但在计算本身中并不直接使用。不过，在查看Ansys Mechanical的迭代收敛历史时，有时会出现「由于条件数差，改用直接法而非迭代法求解」这样的消息。这说明单元形状差导致全局刚度矩阵条件数恶化，反复求解器（如PCG）无法收敛，不得不切换到计算成本高但数值稳定的直接法（稀疏求解器）。

在四面体单元的各积分点（通常是顶点或高斯点）计算雅可比矩阵的行列式，然后取最小值和最大值的比。理想的正四面体中所有点的值相同，比值为1。当单元扭曲时，某些点的|J|较大，其他点接近零。比值低于0.1时大多数求解器会警告，低于0.01时会出错。Abaqus/Standard的默认检查标准正是如此。

六面体由二维四边形面组成。翘曲是指该四边形面不在同一平面上的程度。比如，四边形的4个顶点扭曲成「蝴蝶结」而不是「菱形」。要评估这一点，将面分成两个三角形，计算其法向量的夹角。在Ansys ICEM CFD中，翘曲角推荐5度以下，超过15度则被判定为品质不良。歪斜度（Skewness）通常通过单元内最优内接圆和外接圆的中心偏离来定义。

设定优先级、分阶段检查很有效。第一阶段：使用软件的自动检查功能。Ansys Meshing中的「Mesh Metric」选项可以看到「Element Quality」的柱状图，确认是否存在0.1以下的单元。接下来，对几何学上关键的地方（肋的根部、小孔周围、接触面）进行重点的可视化检查。发现扭曲单元就局部重新划分网格。

不行，要根据分析类型和单元类型而改变。通用的结构静力分析（线性、弹性）中0.1通常可以接受，但非线性分析（塑性、接触）和特征值分析（座屈、振动）需要更严格的标准。例如，Abaqus Explicit手册警告说，为了确保与阻尼相关的稳定时间增分，应避免非常扭曲的单元。CFD的LES（大涡模拟）中，涡的分辨率直接受影响，所以斜度角最好控制在60度以下。

有两个实用的方法。一个是「网格敏感性分析」：将该区域的网格细化一级，检查结果（应力、热流等）的变化是否在例如2%以内。另一个是「结果的局部可视化」：观察品质差的单元周边的应力或温度梯度是否存在不自然的振动，分布是否光滑。如果该单元周围出现异常的「热点」或「轮廓不规则」现象，那该结果就不可信了。

差异很大。首先是指标定义的差异。Ansys多用「Element Quality」这个0到1的综合指标，它基于歪斜度、纵横比、翘曲等的独特计算公式。而Abaqus/CAE的「Verify Mesh」工具分别单独、非常细致地报告歪斜度（Skewness）和纵横比（Aspect Ratio）分布。COMSOL的「网格统计」显示「最小单元品质」「单元体积比」等，并能直接选择并可视化品质低的单元，这点很强。

Ansys Fluent Meshing的「Watertight Geometry」工作流中，在网格生成时通过「Size Function」设定急峻梯度，防止小间隙处网格过度细长。Abaqus/CAE的「Advancing Front」选项本身的网格生成算法更容易产生相对均匀的单元，但复杂几何中容易失败。COMSOL的「自由四面体网格」具有「细小单元大小分布」参数，可以调整大小单元的体积比，防止急剧过渡，从而抑制纵横比恶化。

主要有两点。第一，高级的手动编辑工具丰富。HyperMesh的「smooth」功能能在考虑相邻单元形状的前提下，优化单元节点位置，改善雅可比数。第二，具有行业和分析类型特化的模板和标准。Pointwise在航空航天CFD领域，能严格控制边界层网格的第一层高度和增长率，确保壁面剪应力计算精度。这些在通用CAE软件的标准网格划分器中很难做到的细致控制。

「Negative Jacobian」是单元被倒转（反向）了的状态。比如四边形单元的一个内角超过180度变成凹形。求解器在计算该单元的积分点雅可比时，行列式就会是负数。这是非常严重的问题，表示物理上不可能的状态。日志文件或错误消息通常会输出单元编号，所以可以将该单元可视化。在很多情况下，涉及大变形的非线性分析中，初始网格虽然良好，但变形过程中单元被挤压导致倒转，这就是「网格纠缠」。

有几个对策。第一，在Abaqus/Explicit中加入少量「人工应变能」来防止单元极度扭曲。第二，在Ansys Mechanical的非线性接触设置中打开「自动网格管理」，允许根据变形进行网格「重新划分」。第三，从一开始就用在变形方向有充分层数的结构化网格（有映射网格）。这些是橡胶密封件压缩分析和金属大变形加工模拟的必备技术。

首先看线性求解器的迭代次数和残差日志。在Ansys中，用/STATUS,CMD命令可以输出详细的求解器日志。如果出现「条件数差」之类的消息，或迭代次数卡在最大值不动，就要怀疑网格质量。接下来，简化模型做测试。比如将材料改为线性弹性，接触改为简单支座再运行。如果这样收敛了，那么问题很可能是网格质量，而非非线性。特别是塑性和超弹性等材料非线性会放大扭曲单元中的应力计算误差。

有三个方法。1. 重新划分网格：用Ansys SpaceClaim或Fusion 360等CAD软件把导入的曲面网格先转换成「设计曲面」，然后从中生成高品质网格。2. 平滑化：用HyperMesh的「smooth」功能或开源工具MeshFix这样的工具，调整节点位置来改善雅可比数。3. 过滤：设定品质阈值，自动删除品质差的单元，然后对产生的洞用补片法填补。但这个方法可能也会丧失几何特征（小边等），需要注意。