3. 声腔几何建模与网格划分策略
声学仿真的精度与效率,很大程度上在网格划分阶段就已经被决定。声腔几何建模与网格划分并不是把CAD模型“丢进软件自动划一下”那么简单,它涉及几何清理、网格类型选择、单元尺寸控制、质量检查与修复等一系列相互制约的决策。本节围绕从CAD到声学网格的转换链路,说明网格尺寸与最高分析频率的定量关系,比较六面体与四面体网格的适用场景,并给出网格质量检查与修复的实用方法。
一、从CAD到声学网格:转换链路与几何清理
声学网格的源头通常是结构CAD模型,但结构模型面向制造,包含大量对声学分析无意义的细节:倒角、圆角、螺纹孔、安装凸台、微小缝隙、非流形边等。这些细节若直接进入声学网格,会造成单元数量激增、局部尺寸畸变,甚至导致求解器无法生成体网格。因此,从CAD到声学网格的第一步不是划网格,而是几何清理与声腔抽取。
典型处理流程是:先识别声腔边界——对于封闭空腔(如车内乘员舱、音箱箱体内部),需要从结构模型中抽取内表面并封闭开口;对于半开放声场(如管道、消声器),则需明确入口、出口和壁面。随后进行几何简化:抑制小于目标网格尺寸的特征,缝合自由边,修复面与面之间的间隙和重叠,确保形成水密(watertight)的封闭体积。若几何不水密,体网格生成会失败或产生穿透单元,声学求解时表现为质量矩阵异常、模态频率漂移。
一个实用的判断标准是:几何清理完成后,应能用“封闭体积检查”确认声腔是一个连通的单连通域(或按需保留多个独立声腔),且不存在零厚度面、重复面和自由边。对于复杂装配,可先抽取声腔包络面,再在声学前处理软件中重建为封闭实体。
二、网格尺寸与最高分析频率的关系
声学网格尺寸不是随意设定的,它由最高分析频率和声速共同决定。声波在介质中传播时,每个波长内需要足够多的单元来分辨波峰与波谷。工程上普遍采用的经验规则是:每个波长至少包含6个单元,对精度要求较高时取10个单元。由此得到单元尺寸上限:
h_max = c / (N * f_max)
其中:
h_max —— 最大单元尺寸(m)
c —— 介质声速(m/s),空气常温约343 m/s
f_max —— 最高分析频率(Hz)
N —— 每波长单元数,一般取6,高精度取10
以空气介质、c=343 m/s、N=6为例:若最高分析频率为2000 Hz,则 h_max = 343/(6×2000) ≈ 0.0286 m,即约28.6 mm;若最高频率提高到5000 Hz,则 h_max ≈ 11.4 mm。可以看到,频率提高一倍,网格尺寸需缩小一半,单元数量约增加8倍(三维情况)。这就是高频声学仿真计算量急剧上升的根本原因。
实际划分时,还应考虑局部波长变化:若声腔内部填充不同介质(如空气与多孔材料),声速不同,需按各自声速分别控制网格尺寸。此外,网格尺寸应尽量均匀过渡,相邻单元尺寸比不宜超过1.5,避免因尺寸突变引入数值反射和虚假模态。
三、六面体与四面体网格的选择
声学体网格主要有六面体和四面体两类,选择依据是几何复杂度、计算精度要求和可用工时。
六面体网格的优势在于:单元形状规则,数值耗散小,相同节点数下精度更高;在规则声腔(如长方体管道、简单箱体)中可用较粗网格获得准确模态。其缺点是自动生成困难,对复杂几何需要大量手工分块和映射,耗时高。对于以中低频模态分析为主、几何较规则的声腔,六面体是首选。
四面体网格的优势是适应任意复杂几何,自动划分成功率高,前处理速度快。其代价是相同精度下需要更多单元,且存在数值刚度过硬的问题,高频段误差更明显。对于复杂铸件内部声腔、带大量曲面和分支的消声器,四面体往往是唯一可行选择。若采用四面体,建议使用二阶单元(10节点四面体)来补偿精度损失,但计算量会显著增加。
实际工程中常采用混合策略:主体规则区域用六面体,复杂过渡区用四面体或棱柱体,交界处设置金字塔单元过渡。无论哪种类型,都应保证声腔边界层附近网格与壁面贴合,避免出现阶梯状边界。
四、网格质量检查与修复
网格生成后必须进行质量检查,否则求解器可能报错或给出不可信结果。声学网格重点关注以下指标:
- 单元尺寸:检查最大单元是否超过 h_max,局部高频区域是否加密。
- 雅可比(Jacobian):衡量单元扭曲程度,一般要求大于0.3,低于0.1的单元应修复或删除。
- 长宽比(Aspect Ratio):六面体建议小于5,四面体建议小于3,过大长宽比会降低精度。
- 翘曲度(Warpage):面单元偏离平面的程度,过大时影响声学边界积分精度。
- 自由边与重复单元:自由边意味着声腔未封闭,重复单元会导致质量矩阵重复累加。
修复策略按优先级排列:先处理自由边和重复单元,再修复雅可比过低的单元,最后优化长宽比。常用手段包括局部重划、节点合并、单元删除后补洞、光顺(smoothing)和拓扑优化。修复后应重新检查,直到所有指标进入可接受范围。一个实用的验收标准是:自由边数量为0,雅可比最小值大于0.2,长宽比最大值小于8,且声腔体积与CAD理论体积偏差小于1%。
综上,从CAD到声学网格的转换是一条“几何清理—尺寸计算—类型选择—质量修复”的完整链路。网格尺寸由最高分析频率通过 h_max = c/(N·f_max) 定量确定;六面体适合规则声腔且精度高,四面体适合复杂几何且前处理快;质量检查以自由边、雅可比、长宽比为核心指标,修复后需复检。掌握这条链路,才能在后续模态与声场仿真中获得可信的频率响应和声压分布结果。