1 基本概念
1.1 网格的定义
网格是将连续几何区域离散化后形成的单元集合，通常由节点、边、面和体单元组成。它用于把复杂的连续问题转化为可计算的离散问题，是数值分析中最基础的表示方式之一。
在不同应用中，网格的形式有所差异。二维问题中常见三角形与四边形网格，三维问题中则常见四面体、六面体、棱柱体与多面体网格。网格既可以严格按规则排列，也可以根据几何特征自由分布。
1.2 网格质量的含义
网格质量是指网格在几何外形、尺寸分布、单元排列以及数值性质等方面满足计算要求的程度。它不仅包括单元是否形状规整，还包括网格是否连续平滑、是否贴合边界以及是否存在明显畸变。
一般而言，质量较高的网格能够更稳定地表达几何特征，并使离散方程具有更好的数值表现。相反，质量较差的网格容易引入较大误差，甚至导致求解失败。
1.3 网格质量在工程分析中的作用
在工程分析中，网格质量直接影响仿真的可靠性。高质量网格有助于提高有限元、流体力学和热传导等计算中的精度，使求解器更容易收敛。
同时，网格质量也是模型可信度的重要基础。若网格在局部存在扭曲、翻转或尺寸突变，即便材料参数和边界条件设定正确，计算结果仍可能偏离真实情况。因此，网格质量常被视为前处理阶段的核心指标。
1.4 网格与离散化误差的关系
离散化误差来源于将连续问题近似为离散模型的过程，而网格质量会显著影响这种误差的大小和分布。网格越能准确表示几何边界，并在空间上合理分配单元，误差通常越容易控制。
若网格过于粗糙，可能无法捕捉局部梯度变化；若网格质量差，则误差可能因单元畸变而放大。因而，网格质量并不等同于网格密度，但二者共同决定离散计算的有效性。
2 质量评价指标
2.1 几何类指标
几何类指标主要衡量单元外观是否接近理想形状，以及相邻单元之间是否保持合理关系。这类指标直观且易于理解，常作为质量评估的基础。
2.1.1 单元形状
单元形状用于描述网格单元与标准单元之间的接近程度。理想情况下，单元应尽量规则，如三角形接近等边三角形、四边形接近正方形、四面体接近正四面体。
当单元被拉长、压扁或出现尖角时，形状指标会明显下降。形状不佳通常意味着该单元更容易带来数值误差集中。
2.1.2 角度分布
角度分布反映网格单元内部角与理想角度之间的偏离程度。较均匀的角度分布通常有利于插值和梯度计算。
如果单元内角过小或过大，容易造成局部数值不稳定，尤其在高梯度区域更为明显。因此，角度范围常被作为判断网格优劣的重要参考。
2.1.3 长宽比
长宽比用于衡量单元在不同方向上的尺度差异。对于规则区域，长宽比过大通常意味着单元过度拉伸。
适度的各向异性在某些边界层问题中是有益的，但若长宽比过高且与物理变化方向不匹配，则会降低计算精度，并影响求解收敛。
2.2 数值类指标
数值类指标从代数和映射关系的角度衡量网格是否适合进行数值计算，常用于更深入地判断网格有效性。
2.2.1 雅可比行列式
雅可比行列式反映参考单元到实际单元的映射是否保持局部体积或面积的正向关系。其值为正通常表示单元未翻转，且映射有效。
当雅可比行列式接近零时，说明单元可能发生塌陷，数值计算会变得非常不稳定；若其为负，则往往表示单元已经翻转，属于严重质量问题。
2.2.2 扭曲度
扭曲度描述单元相对于理想几何形态的偏离程度，尤其适用于评估边、面或体单元的空间变形。该指标越大，说明单元越不规则。
扭曲会影响局部插值关系，使计算中变量的变化不再平滑。对于曲面贴合网格，扭曲度尤为重要。
2.2.3 正交性
正交性主要衡量单元边、面或坐标方向之间的垂直程度。正交性好的网格通常更便于离散导数和通量项的计算。
在流体和传热问题中，正交性不足常导致数值扩散增大，甚至引发迭代振荡。因此，许多网格生成过程都会尽量保持较高的正交水平。
2.3 尺寸类指标
尺寸类指标用于描述单元大小及其空间分布是否合理，通常与局部精细化和整体计算成本密切相关。
2.3.1 单元边长
单元边长是最直接的尺度指标，用于表示网格单元的局部大小。边长较小的区域通常能够更精细地捕捉变化。
但单元过小会增加自由度和计算量，因此边长需要在精度与效率之间取得平衡。
2.3.2 尺寸变化率
尺寸变化率指相邻单元尺寸之间的变化速度。若变化过快，容易在网格中形成“突变带”，影响数值传递的连续性。
平缓过渡的尺寸分布更有利于保持误差传播稳定，也能降低局部不适定现象。
2.3.3 局部加密程度
局部加密程度反映某些关键区域内网格细化的强弱。对于边界层、尖角、接触区或高梯度区域，适当加密通常是必要的。
不过，加密应与整体尺寸过渡协调，否则可能引入新的质量问题，使局部改进抵消于全局不平衡。
3 网格类型与质量差异
3.1 结构化网格
结构化网格具有明确的行列或层状拓扑关系，节点连接方式规则，常用于规则几何或可映射区域。
3.1.1 规则性与优势
结构化网格的主要优势在于编号清晰、数据存储紧凑、计算效率较高。其规则性有助于提高求解器的局部访问效率，也方便进行方向性分析。
在简单几何上，这类网格常能获得较好的正交性和尺寸连续性，因此在工程上应用广泛。
3.1.2 典型质量问题
结构化网格的主要问题通常出现在几何复杂区域，例如网格线挤压、局部扭曲或单元拉伸。若映射不当，还可能出现单元极端变形。
此外，结构化网格对边界形状的适应性有限，一旦几何变化剧烈，质量控制难度会明显上升。
3.2 非结构化网格
非结构化网格的单元连接关系更自由，适合复杂几何建模，生成灵活且便于局部加密。
3.2.1 三角形网格
三角形网格常用于二维不规则区域，具有良好的填充能力和适配性。对于复杂边界，它比规则四边形网格更容易生成。
其质量问题主要体现在细长三角形、尖锐角以及局部尺寸不均等方面，这些问题会影响插值稳定性。
3.2.2 四面体网格
四面体网格是三维复杂几何中最常见的形式之一，生成相对方便，适合快速建模和自动划分。
不过，四面体网格在高质量控制上较为敏感，容易出现长条形单元、局部塌陷或正交性不足等情况。
3.2.3 多面体网格
多面体网格由多边形或多面体单元组成，通常具备较好的数值稳健性和较少的单元数量需求。其单元连接关系较为灵活，也便于处理复杂区域。
这类网格在控制质量时重点关注面数分布、单元体积平衡与边界贴合情况，以保证计算性能。
3.3 混合网格
混合网格结合了多种单元类型，常用于兼顾几何适应性与数值效率的场景。
3.3.1 适用场景
混合网格适合复杂工程模型，尤其在局部区域需要高精度、而整体又希望保持较少单元数时更为常见。常见组合包括四边形与三角形、六面体与四面体等。
它能针对不同区域采用不同网格策略，从而提高整体建模灵活性。
3.3.2 质量控制重点
混合网格的质量控制重点在于过渡区处理。不同类型单元之间需要保持连接顺畅，避免出现尺度断裂或拓扑冲突。
若过渡设计不合理，可能造成局部应力集中式的数值异常，降低整体计算可靠性。
4 常见质量问题
4.1 单元畸变
单元畸变是指网格单元偏离理想形状，表现为拉伸、压扁、弯曲或局部扭斜。它会削弱单元的插值性能，并使误差在局部积累。
畸变通常由复杂几何、边界拖拽或生成算法不足引起，是最常见的网格问题之一。
4.2 单元翻转
单元翻转是较严重的质量缺陷，指单元的朝向发生逆转，导致局部映射失效。此时数值计算往往无法正常进行。
这类问题通常对应负雅可比行列式，常被视为必须立即修正的错误。
4.3 网格过粗
网格过粗意味着单元尺寸过大，无法捕捉局部几何细节或物理变化。它会导致分辨率不足，误差增加。
在边界层、尖角、接触区或高梯度场中，过粗网格尤其容易造成结果失真。
4.4 网格过密
网格过密会显著增加计算成本，包括内存占用、运算时间和后处理负担。若细化区域并无明显物理必要，过密网格属于资源浪费。
此外，过密并不必然带来更好结果，如果局部质量差，单纯增加单元数量也难以改善数值表现。
4.5 边界贴合不足
边界贴合不足是指网格未能准确表达几何边界，尤其在曲面、圆角或细小特征处更为明显。此问题会直接影响边界条件的施加精度。
当边界被过度简化时，几何误差会转化为数值误差，进而影响整体分析结果。
4.6 局部突变
局部突变表现为单元尺寸、方向或形状在短距离内快速变化。它会破坏网格的平滑性，使误差传播不连续。
此类问题常出现在局部加密与全局网格衔接不当的区域，需要通过渐变过渡加以控制。
5 网格生成方法
5.1 基于几何划分的方法
基于几何划分的方法先对几何区域进行分区，再在每个子区域内分别生成网格。此类方法对结构清晰的模型较为适用。
它的优点是边界控制较强，能较好地保证关键区域的形状质量，但对前期几何处理要求较高。
5.2 基于生长的方法
基于生长的方法从边界或已有单元出发，逐步向外扩展生成网格。该方法常用于边界层或具有方向性的区域。
其特点是能够控制单元排列方向，并在局部形成较规则的层状结构。不过，如果生长策略不当，也容易出现单元挤压或交叠。
5.3 基于映射的方法
基于映射的方法通过将规则参数域映射到实际几何域来生成网格。它适用于拓扑较规整、可建立对应关系的模型。
这种方法通常能获得较高的一致性和质量，但对几何形状和边界条件有较强依赖。
5.4 基于重构的方法
基于重构的方法通常先从离散点、面片或已有粗网格出发，再重新组织单元结构。它适合处理复杂边界或需要后期修复的模型。
该方法在几何适应性上较强，常用于修补局部缺陷或重新优化网格分布。
6 网格优化与修正
6.1 平滑处理
平滑处理旨在改善单元位置分布，减少局部尖锐变形，使网格更均匀、更连续。
6.1.1 拉普拉斯平滑
拉普拉斯平滑通过将节点调整到邻域节点的平均位置，达到平衡单元形状的目的。它实现简单，应用广泛。
但该方法可能过度收缩或导致边界偏移，因此通常需要配合约束条件使用。
6.1.2 质量驱动平滑
质量驱动平滑以提升特定质量指标为目标，结合单元形状、角度、体积等因素进行优化。相比简单平均，它更关注数值效果。
这类方法计算更复杂，但往往能在保持边界形状的同时获得更好的局部质量。
6.2 重划分与重构
当原网格局部问题较严重时，单纯平滑往往不够，需要重新划分或重构网格。该过程可针对缺陷区域进行整体替换。
重划分常用于局部失效区域的修复，而重构则更适合从头建立符合要求的单元布局。
6.3 局部调整
局部调整是对特定区域进行针对性修改，以尽量少地影响整体网格结构。它通常用于解决小范围质量异常。
6.3.1 节点移动
节点移动通过改变节点位置来改善单元形状和尺寸分布。若移动得当，可以显著减轻畸变和突变问题。
但节点移动需要遵守边界约束和几何限制，否则可能引发新的缺陷。
6.3.2 单元合并与拆分
单元合并与拆分通过调整网格拓扑来改变局部密度与形状。合并适用于消除过密区域，拆分则适用于局部加密。
这种方法灵活性较高，常与其他优化手段结合使用。
6.4 自适应改网格
自适应改网格根据误差估计或解变量变化自动调整网格分布，是提高计算效率与精度的重要手段。它能够将资源集中于关键区域。
自适应策略通常分为加密、粗化和重构三类，适用于多次迭代求解过程，尤其在非均匀场中效果明显。
7 质量评估流程
7.1 前处理检查
前处理检查主要包括几何完整性、边界闭合性、重复节点、悬挂边以及单元连接关系等基础验证。此步骤用于在计算前排除明显错误。
若几何本身存在问题，即使网格划分正确，也可能无法得到可靠结果。
7.2 指标计算
指标计算阶段会对单元形状、尺寸、角度、雅可比值等参数进行统计分析。通常会输出最大值、最小值、平均值及分布情况。
通过这些数据，可以判断网格是否满足不同求解器的要求。
7.3 阈值判定
阈值判定是将计算结果与既定标准进行比较，识别是否存在超限单元。不同软件、不同物理问题可能采用不同阈值。
这一环节有助于快速定位风险区域，并决定是否需要进一步修正。
7.4 可视化诊断
可视化诊断借助颜色映射、剖切图和局部放大图，直观展示质量分布。它能帮助工程人员迅速识别问题集中区域。
相较于纯数值表格，可视化更便于发现空间结构上的异常模式。
7.5 质量报告生成
质量报告通常汇总网格统计信息、异常单元位置、主要指标结果及建议处理方式。它是网格审查和协同沟通的重要文档。
在工程项目中，质量报告也常作为仿真归档的一部分，便于后续复核。
8 典型应用场景
8.1 有限元分析
有限元分析高度依赖网格质量，尤其在应力集中、接触和材料非线性问题中更为明显。良好的网格有助于提高位移与应力计算的稳定性。
若单元质量较差，可能出现局部应力异常放大，从而影响设计判断。
8.2 计算流体力学
在计算流体力学中，网格质量关系到速度、压力与湍流量的离散精度。边界层附近的网格组织尤其关键。
正交性、尺寸渐变和边界贴合程度往往决定了流动模拟的可信度。
8.3 热分析
热分析中，网格质量影响温度场和热流分布的计算效果。对于热源集中或导热路径复杂的结构，局部加密尤为重要。
若网格过粗或变化突兀，可能掩盖真实的温度梯度。
8.4 结构优化
在结构优化中，网格不仅用于计算响应，也影响设计变量更新的准确性。优化过程中，几何会不断变化，因此网格质量需要持续维护。
若网格随形变而恶化，优化结果可能失真，甚至中断迭代。
8.5 多物理场耦合
多物理场耦合问题往往涉及多个变量与多个计算域，对网格一致性要求更高。不同物理场之间的交界面通常是质量控制重点。
良好的网格有助于保持耦合信息传递稳定，减少界面误差。
9 软件与工具
9.1 网格生成软件
网格生成软件用于创建、编辑和优化离散网格，支持不同几何和单元类型。常见功能包括自动划分、局部加密、边界层生成和缺陷修复。
这类软件通常与前处理流程紧密结合，是仿真建模的重要工具。
9.2 网格检查工具
网格检查工具主要用于发现翻转单元、负体积、过大扭曲和尺寸异常等问题。它们通常能输出统计报表和问题定位结果。
在复杂模型中，这类工具有助于快速完成质量筛查。
9.3 后处理与可视化工具
后处理与可视化工具可展示网格分布、质量分级及局部异常区域。通过图形化方式，用户能够更直观地理解网格状态。
这些工具也常与结果场分析联动，用于比较网格质量与计算结果之间的关系。
9.4 自动化脚本与批处理
自动化脚本与批处理可实现网格生成、检查、修正和报告输出的流程化操作。它们适合大规模模型或重复任务。
借助脚本化处理，可以减少人工操作误差，提高工程效率和一致性。
10 相关概念
10.1 几何建模
几何建模是构建待分析对象形状信息的过程，通常是网格生成的前置步骤。几何模型的完整性和简洁性会直接影响后续划分质量。
若几何中存在缝隙、重叠或小特征，网格处理往往会更加复杂。
10.2 数值稳定性
数值稳定性是指离散求解过程中误差是否会被控制在可接受范围内。网格质量是影响稳定性的关键因素之一。
较差的网格可能导致迭代发散、振荡或结果不连续。
10.3 收敛性分析
收敛性分析用于判断随着网格加密，计算结果是否逐渐趋于稳定。它常被用来验证网格是否足够细致。
若不同网格下的结果差异明显，说明当前离散程度仍不足。
10.4 误差估计与控制
误差估计与控制旨在识别误差来源并采取措施限制其影响。网格自适应、局部加密与质量优化都属于常见手段。
在工程仿真中，这一概念与网格质量密切相关，共同服务于更可靠的数值预测。