3. 壳单元网格划分基础

在 HyperMesh 的 NVH 前处理流程中,壳单元网格是车身、内饰支架、薄壁覆盖件等结构最主要的离散形式。壳单元用二维面几何承载三维薄壁结构的弯曲、膜与面内剪切行为,因此网格划分质量直接决定模态频率、声腔耦合与传递函数的可信度。本节围绕四个相互关联的问题展开:选什么单元、划多大尺寸、用什么方式划、划完怎么判合格。

一、单元类型选择:先定阶次与积分,再定拓扑

HyperMesh 中壳单元主要来自 OptiStruct / Nastran 模板下的 CTRIA3、CQUAD4(一阶)与 CTRIA6、CQUAD8(二阶)。NVH 分析通常关注 20~2000 Hz 的中低频段,推荐以 CQUAD4 一阶四边形为主力单元:它计算成本低、刚度矩阵稳定,配合足够的网格密度即可准确捕捉弯曲波长。二阶单元(CQUAD8)在相同节点数下精度更高,但会显著增加自由度与求解时间,一般只在薄壁高频局部或需要精确捕捉曲率时局部使用。

三角形单元(CTRIA3)在纯弯曲工况下存在“剪切锁定”与刚度过硬问题,NVH 中应尽量避免大面积使用,仅用于过渡区或复杂几何填充,且单个三角形区域不宜超过总单元数的 5%~10%。此外,壳单元必须赋予正确的单元属性:PSHELL 卡片中的厚度 T、材料 MID、弯曲惯性矩与剪切修正系数。厚度方向若用中面抽取,需保证偏置(offset)设置正确,否则弯曲刚度会出现系统性偏差。

判断依据:若结构厚度 t 与最小特征尺寸 L 满足 t/L < 1/10,且以弯曲变形为主,优先壳单元;若 t/L 接近 1/5 或存在明显三维应力状态,应考虑实体单元或壳-实体耦合。

二、网格尺寸规划:从分析频率反推单元长度

网格尺寸不是凭经验随意取 5 mm 或 10 mm,而应由目标分析频率与材料波速反推。弯曲波波长公式为 λ = c / f,其中 c 为弯曲波速,近似 c ≈ √(E t² / (12ρ(1-ν²))) 的平方根形式,工程上更常用经验准则:每个弯曲波长至少布置 6~8 个单元,即单元长度 L ≤ λ/6。

以钢制车身板为例,E = 210 GPa,ρ = 7.85×10⁻⁹ t/mm³,ν = 0.3,板厚 t = 1.0 mm。在 200 Hz 时,弯曲波长约为 120~150 mm,按 λ/6 计算,单元长度应不大于 20~25 mm;若分析上限提高到 1000 Hz,波长缩短至约 50~60 mm,单元长度需控制在 8~10 mm。因此,NVH 白车身模型常用 5~10 mm 网格,局部接头与安装点加密至 3~5 mm。

尺寸规划还需兼顾过渡与连接:不同尺寸区域之间应设置 2~3 层过渡单元,避免相邻单元尺寸比超过 2:1;焊点、粘胶、螺栓连接处应保证连接点落在单元节点上,必要时局部细化。网格过密会带来求解规模爆炸,过疏则丢失高阶模态,建议先按目标频率上限估算,再通过一次模态试算验证。

壳单元网格划分决策与质量判读关系图 单元类型选择 CQUAD4 一阶四边形(主力) CQUAD8 二阶(局部高精度) CTRIA3 三角形(≤10%) PSHELL:厚度 / 材料 / 偏置 t/L < 1/10 → 壳单元 网格尺寸规划 L ≤ λ / 6(每波长 6~8 单元) 200 Hz → L ≈ 20~25 mm 1000 Hz → L ≈ 8~10 mm 过渡区尺寸比 ≤ 2:1 连接点落在节点上 划分方式 交互式:几何清理 + 手动 自动:Automesh 面网格 混合:先自动后局部修 关键区域手动控制 规则区域批量自动 网格质量检查指标(NVH 常用阈值) 翘曲度 Warpage ≤ 5°(理想) ≤ 10°(可接受) 长宽比 Aspect ≤ 5(理想) ≤ 10(可接受) 雅可比 Jacobian ≥ 0.7(理想) ≥ 0.5(可接受) 内角 Angle 四边形 45°~135° 三角形 30°~120° 偏斜 Skew ≤ 40°(理想) ≤ 60°(可接受)

三、交互式与自动网格划分:分工与配合

HyperMesh 提供两条主要路径。交互式划分以 Geom → Edge Edit、Quick Edit 清理几何为前提,再用 2D → Automesh 的 interactive 模式逐面控制:可指定单元尺寸、单元类型(quad / tria / mixed)、偏置方向与密度线。它适合接头、圆角、加强筋等几何复杂且对网格走向敏感的区域,能保证单元沿主应力方向排列。

自动网格划分则通过 Automesh 的 automatic 模式或 Batch Mesher 批量处理规则面,按全局尺寸与质量准则一次性生成。实际工程中推荐“先自动、后交互”:先用自动网格快速铺满大面积平板与规则曲面,再对关键连接区、过渡区和质量超标单元进行交互式修补。这样既保证效率,又保留关键区域的可控性。

操作要点:在 Automesh 面板中,element size 输入目标尺寸,mesh style 选择 size and bias 或 edge deviation;勾选 interactive 后可逐边调整节点数。生成后立即用 Check Elements 查看失败单元,不要等到装配完成再回头修网格。

四、网格质量检查指标与判读

HyperMesh 的 Check Elements(快捷键 F10)提供一整套质量指标。NVH 分析中应重点关注以下五项,并建议在 Settings → Quality Criteria 中按项目要求设定阈值:

  • 翘曲度 Warpage:四边形四个节点不共面的程度。翘曲过大导致弯曲刚度虚假增大,模态频率偏高。理想 ≤ 5°,可接受 ≤ 10°。
  • 长宽比 Aspect Ratio:单元最长边与最短边之比。过大使刚度矩阵病态,影响低频模态精度。理想 ≤ 5,可接受 ≤ 10。
  • 雅可比 Jacobian:单元映射的畸变程度,反映积分点质量。低于 0.5 时结果不可信,理想 ≥ 0.7。
  • 内角 Angle:四边形内角宜在 45°~135°,三角形宜在 30°~120°,避免出现尖角或钝角单元。
  • 偏斜 Skew:单元偏离理想形状的角度,理想 ≤ 40°,可接受 ≤ 60°。

判读原则:先看失败数量,再看失败位置。若失败单元集中在非关注区域且数量占比小于 1%,可暂时保留;若出现在焊点、悬置安装点、声腔耦合面等关键路径上,必须局部重划。检查完成后用 Quality Index 面板按颜色分级显示,红色单元优先处理,黄色单元评估影响后决定是否修。

常见误区:只追求“全部绿色”而过度细化,导致模型规模失控。NVH 网格质量应服务于目标频率与关注区域,非关键区域允许适度放宽,但连接点与载荷路径附近必须严格达标。

五、可执行检查清单(示例)

检查项目标值判读动作
单元类型占比四边形 ≥ 90%三角形超限区域局部重划
目标频率对应单元长度L ≤ λ/6按 200 Hz / 1000 Hz 分别核算
翘曲度≤ 10°超标单元用 Edit Element 调整节点
长宽比≤ 10长条单元拆分或重新布种
雅可比≥ 0.5畸变单元局部重划
连接点节点对齐100%焊点 / 螺栓处强制节点重合

完成上述四步后,壳单元网格才具备进入装配、连接与求解阶段的基本条件。下一阶段将在此基础上处理焊点、粘胶与螺栓连接,并建立声腔与结构耦合关系。