高阻尼板系统混合射线追踪-统计能量分析（Ray-Tracing-SEA）振动声学预测方法的数值与实验研究

《Advances in Mechanical Engineering》：Hybrid Ray-Tracing-SEA method for vibroacoustic prediction in heavily damped plate systems: A numerical and experimental study

　　

高阻尼系统的精确模拟在振动声学工程中仍面临重大挑战，这源于其能量分布不均匀和高衰减特性，违反了经典统计能量分析（Statistical Energy Analysis, SEA）的基本假设。为解决这一问题，研究人员提出了一种混合射线追踪-统计能量分析（Ray-Tracing-SEA）方法，在SEA框架内整合了基于射线的精细化能量追踪技术。研究人员针对一个三铝板系统在无阻尼和重阻尼两种条件下开展了数值与实验研究。通过功率输入法（Power Input Method, PIM）测量了各板的总损耗因子（Total Loss Factor, TLF），并记录了板间能量级差（Energy Level Difference, ELD）。通过将经典SEA与混合方法的数值预测结果同实验数据进行对比，评估了预测精度。研究人员开展灵敏度分析，量化了混合方法在射线数量、反射次数以及输入损耗因子扰动条件下的收敛性与稳健性，并通过弯曲模态计数（bending mode counts）和模态重叠因子（modal overlap factors）对低频适用界限进行了定量考察。结果表明，当2号板处于重阻尼状态时，混合方法在500 Hz至5000 Hz频率范围内的频带平均能量级差误差保持在约2.5 dB以内，而经典SEA的误差高达约7 dB。该混合方法为采用高阻尼处理的结构分析提供了一种更为可靠的工具。

结构传声振动控制对工程结构的完整性与运行可靠性至关重要。未受控的振动会导致材料疲劳、服役寿命缩短以及交通运输系统的运行故障。无源阻尼处理，包括自由层阻尼（free-layer damping）、约束层阻尼（constrained-layer damping）以及调谐质量阻尼等，仍是工程中最广泛采用的减振降噪措施，其内损耗因子（internal loss factor）可高达0.21至0.53，具备显著的宽频带能量耗散能力。然而，当板结构损耗因子超过0.2时，强烈波动衰减导致振动能量集中于激励点与连接处附近，形成显著的空间能量梯度。这种非均匀能量分布直接违反了经典统计能量分析（Statistical Energy Analysis, SEA）所依赖的能量均分（equipartition）与漫射声场（diffuse field）假设，导致经典SEA预测精度随阻尼增大而显著恶化。尽管混合射线追踪-统计能量分析（Ray-Tracing-SEA）方法通过在SEA框架内引入基于射线追踪（ray tracing）的能量追踪以放松均分假设，在轻阻尼结构中已展现潜力，但现有研究仍存在三方面空白：其一，该公式在重阻尼区（总损耗因子Total Loss Factor, TLF大于0.2）的验证仍停留在纯数值阶段，缺乏严格实验支撑；其二，方法关键数值参数（射线数量、反射次数）及输入损耗因子测量不确定度对预测稳健性的影响从未被量化；其三，文献中区分轻、重阻尼的0.2经验阈值缺乏与频率、子系统几何关联的定量物理依据。这些不足限制了该方法在高阻尼处理工程结构中的可信应用。

针对上述问题，研究人员围绕一个由单块铝板切割弯折而成的三铝板系统（three-plate aluminium system）开展联合实验与数值研究，分别测试无阻尼与重阻尼工况。研究人员通过功率输入法（Power Input Method, PIM）实验测定各板总损耗因子，并测量外侧轻阻尼板间的能量级差（Energy Level Difference, ELD）。主要研究结论如下：首次在重阻尼条件下（2号板TLF在1250 Hz以上达0.39–0.47）完成了混合Ray-Tracing-SEA方法的实验验证；定量建立了射线数量、反射次数及输入损耗因子扰动的灵敏度分析；基于弯曲模态计数（bending mode counts）与模态重叠因子（modal overlap factors）确定了低频适用下限；将0.2经验损耗因子阈值与约13 dB的单向弯曲波衰减定量关联，证明该阈值是随频率变化的折中值而非绝对物理常数。研究证实，混合方法为含高阻尼处理的结构提供了更可靠的振动声学分析工具。

为开展上述研究，研究人员综合运用四类关键技术方法：其一，数值建模方面，建立经典SEA模型与混合Ray-Tracing-SEA模型；后者将重阻尼板视为耦合单元，采用路径追踪（path tracing）在板边界布置均匀点源并发散射线，依据弯曲波能量衰减规律与角部连接处角度相关透射系数递归更新每条射线的反射与透射能量，进而在界面上按能量加权平均得到平均透射系数，计算等效耦合损耗因子（equivalent CLF）后代入SEA矩阵求解。其二，实验测量方面，利用PIM通过力锤激励与加速度响应测量获取各板TLF，依据输入能量与子系统总振动能量之比计算损耗因子。其三，ELD实验测量方面，采用单点激振器施加白噪声，在外板各布设20个空间测点，通过空间均方加速度之比计算板间能量级差。其四，系统灵敏度分析方面，在保持其他参数不变条件下，定量改变射线数量、最大反射次数及输入TLF扰动量，评估预测ELD的收敛性与稳健性，并基于模态密度与模态重叠因子确定低频有效性判据。

研究结果

总损耗因子（Total Loss Factor）
研究人员通过PIM对三铝板系统各板进行测量。在无阻尼预实验中，3号板实测TLF与SEA预测值在250 Hz以上吻合良好，证实了PIM测量流程的可靠性。施加约束层阻尼（constrained layer damping）与颗粒阻尼（granular damping）联合处理后，2号板TLF在大部分频段超过0.2，被明确归类为重阻尼系统；1、3号板TLF在两种工况下最大偏差低于8%，表明该阻尼处理未显著改变外板固有耗散特性。实测TLF被直接用作后续数值模型的输入参数。

能量级差（Energy Level Difference）
在无阻尼条件下，经典SEA与混合方法的预测值在250 Hz以上与实验值均吻合良好，因为轻阻尼时等效CLF与经典CLF差异甚微。当2号板处于重阻尼状态时，500 Hz至5000 Hz范围内，混合方法的频带平均ELD误差保持在约2.5 dB以内（平均绝对偏差约1 dB），而经典SEA误差高达约7 dB（平均绝对偏差约3 dB）。具体而言，经典SEA在500–1000 Hz表现为对实测ELD低估约5 dB，而在1250–6300 Hz转为高估，最大偏差出现在2 kHz附近；混合方法在两种频段均准确跟踪实测值。在315 Hz和400 Hz过渡频段，因1、3号板模态重叠因子仍低于1，两种模型均高估实测ELD 4–8 dB；315 Hz以下两模型均失效，此归因于弯曲模态数不足而非方法本身缺陷。5000 Hz以上，由于2号板TLF回落至0.2以下且激励带宽限制影响，预测ELD逐渐回归无阻尼状态。

灵敏度分析（Sensitivity Analysis）
在无阻尼构型下，混合方法预测与经典SEA最大差异为0.6 dB，验证了算法实现正确性。针对平均透射系数τ₂₃的计算表明，射线总数从50×90增至200×360时，其变化小于0.01%，结果对射线离散化极不敏感；而低于50×90的粗网格因混叠矩形板周期性轨道族，偏差可达约10%。反射次数的影响仅在无阻尼时显著：5000 Hz无阻尼条件下，5次反射高估τ₂₃达1.4%，20次反射时偏差降至0.1%以下；阻尼板条件下5次反射即收敛。损耗因子扰动分析显示，2号板TLF改变±5%和±10%时，预测ELD变化分别小于0.20 dB与0.41 dB；三板TLF同时存在10%测量不确定度时，预测ELD最大变化为0.8 dB，即使20%不确定度也不超过1.6 dB。结果表明混合方法对数值参数与输入测量误差均具有足够稳健性。

讨论与结论

低频限制（below 315 Hz）
基于弯曲模态密度公式n(f) = S√3/(hc_L)计算，每三分之一倍频程模态数N = 0.018f。满足Hopkins准则N ≥ 5的起始频率为277 Hz，故315 Hz为首个有效频段。315 Hz处弯曲波长为24 cm，沿0.9 m板边不足四个波长，射线概念仅边际适用。2号板模态重叠因子M从100 Hz的1.5增至2 kHz的73，呈波动主导特性；而1、3号板模态重叠因子在1 kHz以下保持低于1，其频带平均能量仍保留与具体实现相关的模态成分，导致315 Hz以上仍存在不可消除的频带间散布。315 Hz以下三板模态密度均不足，SEA统计基础完全丧失，两种方法均失效。

经典SEA误差的符号反转机制
该现象可通过两种误差机制定量解释。在1250 Hz以上，2号板TLF达0.39–0.47，平均透射系数τ₂₃升至0.46–0.47，比漫射场假设值1/3高出38%–41%，经典SEA严重低估板间耦合，导致ELD高估2–7 dB。在315–1000 Hz范围内，1、3号板模态重叠因子仍低于1，每频带仅含6–18个弯曲模态，离散模态主导能量分配，经典SEA高估有效耦合，低估ELD达5 dB。约1 kHz处的误差交叉点对应外板模态重叠因子穿过1以及2号板TLF超过0.35，标志着从模态主导误差向阻尼主导误差的过渡。

0.2损耗因子阈值的定量阐释
对0.9 m板长，弯曲波单向跨越2号板的衰减因子为exp(-ηk_BL/2)。当η = 0.2时，100 Hz、500 Hz与1 kHz处衰减分别为5.8 dB、12.9 dB与18.2 dB，在中频段超过13 dB（95%能量被耗散）。定义子系统内空间能量变化超过一个数量级的临界阈值η_crit为250 Hz时0.29、500 Hz时0.20、1 kHz时0.14，故0.2实质是中频折中值，此时经典SEA漫射场误差已成为主导不确定性。该阈值与13 dB准则是工程惯例而非物理常数；若采用10 dB准则，η_crit将降低约23%。射线追踪结果独立验证：η = 0.2时，τ₂₃在200 Hz–2 kHz内偏离漫射场值22%–34%，导致ELD误差达1–3 dB。

研究结论（Conclusions）
本研究表明，联合约束层阻尼与颗粒阻尼处理可显著拓宽板结构有效阻尼频带。在重阻尼工况下，混合Ray-Tracing-SEA方法将500–5000 Hz范围内的频带平均ELD预测误差控制在约2.5 dB以内，显著优于经典SEA最高约7 dB的误差；该精度提升源于其对阻尼板内非均匀能量分布与局部衰减特性的解析能力。约315 Hz以下，因模态密度不足及SEA统计假设失效，两种方法精度均显著下降。整体而言，混合方法为高阻尼结构振动声学建模提供了更为可靠的工具，有助于指导减振降噪设计。需强调的是，当前验证仅在子系统能量级差层面完成；板内空间能量衰减这一射线追踪修正所针对的核心物理机制尚未被直接测量。后续研究应通过沿阻尼板两连接点间直线扫描获得的频带平均加速度衰减率，以及子系统内最大与空间平均能量密度之比等可量化指标，实现对板内空间能量分布的直接验证。