1. 结构声耦合问题导论
当一台发动机在台架上运转时,机体表面的振动会推动周围空气形成声波;反过来,这些声波又会以压力的形式作用回结构表面,使结构振动发生改变。这种结构与声场之间双向、同时发生的相互作用,就是结构声耦合。它不是"结构先振动、声场后响应"的单向因果链,而是一个闭合的反馈回路:结构壁面法向速度决定声场质点速度边界条件,声场压力又作为分布载荷叠加到结构动力学方程右端。理解这一点,是后续所有建模与预测工作的起点。
从数学上看,未耦合的结构方程可写作 M·ü + C·u̇ + K·u = F,未耦合的声腔方程可写作 Ma·p̈ + Ca·ṗ + Ka·p = q。耦合项出现在两处:结构方程右端增加 −S·p(S 为耦合面积矩阵,压力以法向载荷形式作用),声腔方程右端增加 ρ0·ST·ü(结构加速度作为声源项)。两个方程必须联立求解,这正是耦合问题的本质。
耦合的强弱可以用一个无量纲量来初步判断:结构模态频率 fn 与声腔模态频率 fa 的接近程度。当两者接近时,能量在结构与声场之间高效交换,耦合效应显著;当两者远离时,可近似解耦处理。工程中常用 模态耦合系数 来量化这种匹配程度,它取决于结构模态振型与声腔模态振型在耦合界面上的空间重叠积分。
典型工程场景各有其耦合特征。汽车 NVH 中,车身地板、顶棚等薄壁板件与乘员舱大空腔构成耦合系统,关注频段通常在 20–200 Hz,此时声腔模态密集,板件模态与声腔模态容易发生共振匹配,产生令人不适的低频轰鸣。飞机舱内噪声则面临螺旋桨或发动机的强谐波激励,蒙皮与框段组成的加筋结构在轻质化设计下刚度较低,与客舱声腔耦合后会在特定谐波频率上放大噪声。船舶舱室的问题更复杂:船体板既与舱室空气耦合,又与水介质耦合,形成"水—结构—空气"三重耦合路径。家电噪声虽然尺度小,但壳体薄、空腔小,耦合频率往往落在人耳敏感的 200–2000 Hz 范围。
噪声预测的核心挑战可以归纳为五个方面。第一,宽频带:从几十赫兹的低频模态响应到几千赫兹的声学短波长,单一方法难以覆盖全频段。第二,多尺度:结构细节(筋条、螺栓连接)与声学大空间尺度差异可达三个数量级。第三,强耦合:当结构与声场模态频率接近时,必须联立求解,不能分别计算后叠加。第四,边界条件不确定:实际安装状态、密封条件、材料阻尼与仿真假设存在偏差。第五,计算规模大:三维声腔有限元/边界元模型自由度可达数百万,直接求解代价高昂。
从预测方法的角度看,低频段(约 20–200 Hz)通常采用结构有限元 + 声学有限元的耦合求解,直接组装耦合矩阵并求解特征值或频响;中频段(约 200–1000 Hz)可采用混合 FE-SEA 方法,将确定性结构与统计声学处理结合;高频段(1000 Hz 以上)则更多依赖统计能量分析(SEA),以能量为变量描述子系统间的功率流。方法的选择本质上取决于耦合强度、模态密度和计算资源三者的平衡。
以一个简化的汽车乘员舱模型为例:假设舱室尺寸约 2.5 m × 1.5 m × 1.2 m,空气密度 1.21 kg/m³,声速 343 m/s,则第一阶声腔模态频率约为 f = c/(2L) = 343/(2×2.5) ≈ 68.6 Hz。若车身地板的一阶弯曲模态恰好落在 65–75 Hz 附近,两者将发生强耦合,在 68 Hz 左右出现明显的声压峰值。此时若忽略耦合、仅做单向激励计算,预测的声压级可能偏低 3–6 dB,导致降噪方案设计偏差。
理解结构声耦合的物理本质,意味着在建模之初就要回答三个问题:结构模态与声腔模态是否接近?耦合界面的面积与形状如何影响能量交换?关注频段内哪种方法最合适?这三个问题的答案,将贯穿后续所有章节的建模与预测流程。