3. 网格划分策略总览：全局网格与局部网格、网格密度控制原则、计算精度与资源消耗的平衡

网格划分是CAE仿真分析中承上启下的关键环节——它既决定了物理模型能否被数值方法准确离散，也直接决定了求解器的计算规模与收敛行为。许多工程师在完成几何清理后，往往凭经验或默认设置直接生成网格，结果要么因网格过密导致计算资源浪费，要么因关键区域网格过疏而错过重要物理特征。本章从全局网格与局部网格的职责划分入手，系统阐述网格密度控制的核心原则，并给出计算精度与资源消耗之间的工程平衡方法。

3.1 全局网格与局部网格的职责划分

3.1.1 核心原理：两层网格体系的互补逻辑

全局网格（Global Mesh）定义了整个计算域的基础离散尺度，它决定了模型的总单元数量级和整体求解规模。全局网格尺寸通常由最小几何特征尺寸、边界层厚度或波传播波长等因素决定。例如，在结构分析中，若最小螺栓孔直径为6mm，全局网格尺寸通常取2~3mm以保证孔周至少有2~3层单元；在CFD分析中，全局网格尺寸需满足湍流积分尺度的分辨率要求。

局部网格（Local Mesh）则是在全局网格基础上，针对特定区域进行加密或粗化处理。其核心价值在于：在不显著增加全局单元数量的前提下，提高关键区域的求解精度。局部加密区域通常包括：应力集中区（圆角、缺口、焊缝）、几何突变处（台阶、薄壁）、物理量梯度大的区域（温度边界层、激波面、接触界面）以及需要精确捕捉的边界层区域。

3.1.2 判断方法：如何确定全局与局部的边界

工程实践中，建议按以下流程确定网格策略：

1. 几何特征扫描：统计模型中的最小孔、最小圆角、最薄壁厚、最大长宽比等特征尺寸，取最小特征尺寸的1/3~1/5作为全局网格基准尺寸。
2. 物理场预估：根据载荷类型和边界条件，预估应力集中系数、温度梯度范围或流速变化区域，将这些区域标记为局部加密候选区。
3. 全局网格试算：先以较粗的全局网格（约为基准尺寸的2倍）进行一次快速试算，观察整体物理场分布是否合理，识别出未预判到的高梯度区域。
4. 局部加密迭代：在试算结果基础上，对高梯度区域施加局部网格控制，逐步加密直至关键结果收敛。

3.1.3 案例拆解：某支架结构的网格策略制定

某汽车底盘支架，材料为铝合金，承受交变疲劳载荷。初始几何包含：4个直径8mm的安装孔、2处R3mm的过渡圆角、1处厚度仅2.5mm的薄壁区域。若直接采用全局2mm网格，单元总数约180万，单次求解耗时45分钟；若全局4mm+局部加密，单元总数降至65万，求解耗时12分钟。

方案B与方案A的最大应力偏差仅0.17%，但计算效率提升近4倍。方案C虽然最快，但应力结果偏低15.6%，无法满足工程精度要求。此案例验证了“全局控制规模、局部保证精度”策略的有效性。

3.1.4 常见误区与总结

• 误区一：认为全局网格越细越好——实际上当网格细化到一定程度后，精度提升趋于平缓，而计算代价呈指数增长。
• 误区二：局部加密区域设置过大——过度加密会浪费计算资源，加密区应控制在模型总体积的10%~20%以内。
• 误区三：忽略全局网格与局部网格之间的过渡——相邻区域网格尺寸突变会导致数值误差和收敛困难。

3.2 网格密度控制原则

3.2.1 核心原理：物理梯度与几何曲率双驱动

网格密度的本质是对物理场离散误差的空间分布控制。有限元/有限体积方法的误差估计理论表明，局部离散误差与单元尺寸的p次方成正比（p为插值阶数），且与物理量的高阶导数（梯度变化率）直接相关。因此，网格密度控制应遵循“梯度大则密、梯度小则疏”的基本原则。

具体而言，需要加密的区域包括：

• 应力/应变集中区：缺口、圆角、裂纹尖端等，应力梯度可达平均值的5~10倍。
• 热边界层/流体边界层：近壁面速度与温度梯度极大，需布置边界层网格。
• 接触/冲击区域：接触压力分布高度不均匀，需要细化接触面网格。
• 几何曲率大的表面：曲率半径小于5倍网格尺寸时，需加密以保证几何保真度。

3.2.2 操作流程：网格密度控制的五步法

1. 梯度预判：根据载荷路径和约束位置，在CAD模型上标注预期的高梯度区域。
2. 尺寸设定：为每个加密区设定目标单元尺寸，通常为全局尺寸的1/2~1/5。
3. 过渡控制：设置相邻区域网格尺寸过渡比不超过1.5~2.0，避免尺寸突变。
4. 边界层布置：对壁面区域设置边界层网格，第一层厚度满足y+要求（CFD）或接触精度要求（结构）。
5. 质量检查：检查偏斜度（Skewness）<0.85、长宽比（Aspect Ratio）<10、最小角度>18°等指标。

3.2.3 案例拆解：电子散热器热应力分析

某功率模块散热器，铝制翅片结构，翅片厚度1.2mm，间距3mm，底部与IGBT接触。热-结构耦合分析中，翅片根部热梯度极大，而翅片顶部温度梯度较小。初始方案采用全局1mm网格，单元数420万，求解耗时3.5小时。优化后采用全局2mm网格，翅片根部加密至0.5mm，接触面设置0.3mm边界层，单元数降至180万，求解耗时55分钟。关键结果对比：翅片根部最大热应力从初始方案的187MPa变为优化后的192MPa（偏差2.7%），但计算效率提升近4倍。更重要的是，优化方案捕捉到了翅片根部更真实的应力分布形态，为疲劳寿命评估提供了更可靠的输入。

3.2.4 常见误区与总结复盘

• 误区一：只关注单元尺寸而忽略单元质量——一个偏斜度0.95的坏单元可能比10个正常单元带来的误差更大。
• 误区二：加密区域边界处理不当——加密区与粗网格区之间应设置至少3~5层过渡单元。
• 误区三：忽略网格方向性——对于各向异性物理场（如薄壁结构的厚度方向），应采用各向异性网格而非均匀细化。

3.3 计算精度与资源消耗的平衡

3.3.1 核心原理：网格无关性验证与收敛判据

计算精度与资源消耗的平衡，本质上是寻找“满足工程精度要求的最小计算代价”。工程上通常采用网格无关性验证（Mesh Convergence Study）来确定最优网格密度。其核心思想是：逐步细化网格，观察关键结果（如最大应力、最大变形、流量、温度等）的变化趋势，当结果变化小于设定阈值（通常为1%~3%）时，认为网格已达到无关性要求。

网格无关性验证的数学基础是Richardson外推法：若已知三套不同密度网格的计算结果，可估算精确解的渐近值，从而判断当前网格的误差水平。实际操作中，通常采用“粗-中-细”三套网格进行对比分析。

3.3.2 判断方法：网格无关性验证的标准流程

3.3.3 案例拆解：某压力容器接管区的网格无关性验证

某压力容器，设计压力15MPa，接管与筒体连接处存在结构不连续。分别采用三套网格方案进行弹性应力分析：

• 粗网格：全局20mm，接管区局部8mm，单元数8.2万，最大等效应力312MPa。
• 中网格：全局15mm，接管区局部5mm，单元数23.5万，最大等效应力356MPa。
• 细网格：全局10mm，接管区局部3mm，单元数67.8万，最大等效应力361MPa。

中网格与细网格的应力偏差仅1.4%，满足3%的收敛判据。因此选择中网格方案作为最终计算网格。该方案相比细网格节省了65%的计算时间，同时保证了工程精度。进一步分析发现，粗网格偏差达12.3%，主要原因是接管区应力集中梯度未被充分捕捉。

3.3.4 常见误区与总结复盘

• 误区一：只做一次网格无关性验证就认为所有工况都适用——不同载荷工况下，高梯度区域位置可能不同，需分别验证。
• 误区二：过度追求“完美收敛”——工程上3%~5%的偏差通常已足够，过度细化只会浪费资源。
• 误区三：忽略硬件资源约束——即使网格无关性验证表明需要更细网格，若超出可用内存或计算时限，需通过子模型技术或降阶模型来平衡。

3.4 综合策略：从全局到局部的完整决策框架

3.4.1 核心原理：分层决策与迭代优化

综合网格策略应遵循“先全局、后局部、再验证”的三层递进逻辑。第一层确定全局网格尺寸，控制整体规模；第二层施加局部加密，聚焦关键区域；第三层通过网格无关性验证确认精度达标。整个过程需要2~3轮迭代，每轮迭代都应记录关键结果和计算耗时，形成决策依据。

3.4.2 操作流程：网格策略制定的完整清单

1. 几何预处理：清理小特征（倒角、小孔）、简化非承载细节，保留关键几何特征。
2. 全局尺寸设定：基于最小特征尺寸和物理场预估，设定初始全局网格尺寸。
3. 局部加密区域标记：在CAD模型中创建命名选择（Named Selections），标记加密区域。
4. 网格生成与质量检查：生成网格后，检查单元质量指标，修复不合格区域。
5. 试算与结果评估：运行初步求解，观察物理场分布，识别未预判的高梯度区。
6. 网格调整与再计算：根据试算结果调整局部加密设置，重新求解。
7. 网格无关性验证：按3.3.2节流程进行验证，确认结果收敛。
8. 最终确认与文档记录：记录网格参数、求解设置、结果数据，形成可追溯的分析报告。

3.4.3 案例拆解：新能源汽车电池包结构分析

某新能源汽车电池包下壳体，尺寸约1800mm×1200mm×80mm，材料为压铸铝合金。分析目标为验证振动载荷下的结构完整性。初始方案采用全局5mm网格，单元数高达520万，单次模态分析耗时6小时，且内存占用接近工作站上限（64GB）。优化策略如下：

• 全局网格粗化：将全局尺寸调整为8mm，单元数降至210万。
• 局部加密：对螺栓连接区域（直径10mm）加密至2mm，对加强筋根部圆角加密至1.5mm。
• 子模型技术：对电池包与车身连接的关键支架，采用子模型技术进一步细化至0.5mm。

优化后单元总数降至180万，求解时间缩短至1.5小时。关键结果对比：第一阶固有频率从初始方案的23.8Hz变为优化后的23.5Hz（偏差1.3%），支架根部最大应力从初始方案的89MPa变为优化后的94MPa（偏差5.6%）。虽然应力偏差略大，但通过子模型验证确认94MPa更接近真实值，且满足安全系数要求。该案例表明，通过合理的全局-局部配合，可以在保证精度的前提下将计算效率提升4倍。

3.4.4 总结与经验提炼

网格划分策略的最终目标不是追求“最细的网格”，而是追求“最合适的网格”——在满足工程精度要求的前提下，以最小的计算代价获得可靠的仿真结果。这需要工程师对物理问题有深入理解，对网格技术有系统掌握，并通过持续的案例积累形成自己的判断直觉。本章所阐述的全局与局部配合、密度控制原则、精度与资源平衡方法，为建立这种判断力提供了系统化的方法论基础。