Chinese Journal of Ship Research
基于几何声学的船舶舱室声学设计方法
冯爱景1，魏强2，3，张大海 2，3
1 201108上海船用柴油机研究所，上海2 430064中国舰船研究设计中心，湖北 武汉3 430064中国舰船研究设计中心 船舶振动噪声重点实验室，湖北 武汉 摘 要：［目的］为了在船舶舱室的众多噪声传递路径中选取最佳噪声控制位置及控制措施，［方法］基于几何声学理论中的声线跟踪法，考虑舱壁声透射的作用，提出声线搜索法。模拟船舶多舱室声场的分布，计算舱室声压。通过搜索目标舱室的供能声线，计算不同位置舱壁对目标舱室噪声的声灵敏度，根据灵敏度计算结果，设计船舶舱室降噪方案，优化舱室中高频噪声。［结果］利用该方法优化典型舱室噪声，噪声降低了7.3 dB。
［结论］通过与统计能量法的对比分析，验证该方法可行，可指导船舶舱室降噪精细化设计。关键词：舱室降噪；几何声学；声线跟踪法；声线搜索法；优化设计中图分类号：U661.44 文献标志码：A DOI：10.3969/j.issn.1673-3185.2017.04.008
Acoustic design method of ship's cabin based on geometrical acoustics FENG Aijing1，WEI Qiang2，3，ZHANG Dahai2，3 1 Shanghai Marine Diesel Engine Research Institute，Shanghai 201108，China 2 China Ship Development and Design Center，Wuhan 430064，China 3 National Key Laboratory on Ship Vibration and Noise，China Ship Development and Design Center， Wuhan 430064，China
Abstract：In light of the question of how to select the best noise control position and measures in the large noise transmission path of the cabins of a ship, based on the acoustic ray-tracing method in the theory of geometrical acoustics, and by considering the effect of the sound transmission of the bulkhead, this paper proposes the sound line search method. It is used to calculate the sound pressure of a ship's cabin, allowing the sound field distribution of multiple compartments to be simulated. The paper proposes a sound ray-searching method in which the acoustic sensitivity of different positions of the bulkhead to the noise of the target cabin is calculated by searching for the sound ray passing the target cabin. According to this, a cabin noise reduction plan can be designed to optimize medium and high frequency cabin noise. With this method, the noise of a typical cabin can be optimized and reduced by 7.3 dB. Through comparative analysis with the statistical energy method, it is proven that the method is feasible and can guide the refined design of noise reduction in ships' cabins. Key words：noise reduction of ship's cabin；geometrical acoustics；ray-tracing method；ray-searching method；optimal design
0引言
随着船舶行业科技的发展，对船舶性能的要求越来越高，例如，为保障人员身心健康以及设备正常工作，国内外对船舶舱室噪声的要求日益严格，船舶舱室噪声的控制优化成为船舶设计不得不考虑3的问题。船舶舱室噪声控制的 个主要措施包括： 1）舱室布局。将高噪声要求的舱室布置在远离声源处或进行隔离。2）低噪声设备选用。根据总体指标及结构，选择符合标准的主机设备。3 ）声源控制。对声源，如主机、空调等做减振、降噪、消声处理。对于噪声水平还未达标的舱室，就需进一步开展噪声控制优化设计，目前主要的降噪措施是在舱室内敷设吸声和隔声材料。然而船舶具有多舱室、多噪声源、多传递路径等特点，且舱室空间以及舱壁承重存在限制，盲目敷设吸、隔声材料并不一定能取得良好的降噪效果。针对这一问题，众多学者开展了关于舱室声学优化的研究。Koo［1］利用边界元方法分析了结构形状变量对声参数Le Bot 等［2 ］引入热力学中熵的概念，的灵敏度。基于损耗因子贡献法，提出了熵赋权图法，可计算统计能量法系统中的主要能量传递路径。高处等［3］在熵赋权图法的基础上，采用前K条最短路径算法计算了统计能量系统中的前K条主要能量传递路径。解决舱室噪声优化问题的关键是求得舱壁各参数对目标舱室的声灵敏度，其主要方法有：直接求导法、摄动法、有限差分法和矩阵特征值反问题求解法［4-6］。由于船舶舱室声学模型复杂，目前还不能从有限元法和统计能量法参数矩阵入手来直接定量计算参数灵敏度。鉴于此，本文拟引入声场空间计算的另一主要方法——几何声学，基于声线跟踪法，提出跨空间、多舱室声场计算方法以及声学优化方法——声线搜索法。通过实例计算，与专业声学软件的对比分析，证明该方法可有效优化中、高频噪声，能指导船舶舱室降噪精细化设计。
1 基于声线跟踪法的多舱室声学模拟
几何声学是用射线的观点研究声学问题，被广泛应用于建筑行业。主要方法包括：声线跟踪法、虚声源法和混合法，除此之外，还有声束跟踪法、声辐射度法等［7］。
1.1 声线跟踪法
声线跟踪法是假设声音沿直线传播，不考虑 其波动性，通过计算声线在传播过程中能量的变拟［8］。由声源化及路径来对目标区域进行声场模处发射多条携带相同能量的声线，声线在空间中射［9］，散射以直线传播，遇到壁面，发生反射和散一定的概率随机发生。与壁面碰撞时，声线能量1-减至数字原声线的 α 倍（ α为壁面吸声系数）10 。［ ］此后，声线携带剩余声能沿着新方向继续传播。当声线能量与初始声线数能量的比值小于设定的Energy Discontinuity Percent，EDP）［11］能量限值（时，认为该声线反射与透射充分，声线剩余能量可忽略，停止跟踪。通过遍历跟踪所有声线传播路径，叠加所有声线能量来计算空间声强。该方法是采用单个频率逐次计算，若某个频率下的声源功率为 I ，则该频率下的初始声线 In 可表示为1 In = ns(x 1 y z) （ ）式中：n为声线数量；s(x y z)为初始声线由声源发射时的初始传播路径。以球状声源为例，生成由球心均匀向外发散1的声线，如图 所示。声线采用等立体角均匀划分的方法生成，采用球坐标系，声线角度公式如式（2）~式（8 ）所示［ 9，12 ］。式中： θ 为声线与xoy 平面的夹角；φ为声线在xoy平面上的投影同 y 坐标式（2）~式（8）中轴之间的夹角。 各参数的含义详见文献［9，12］。
1.2 声线跟踪法多空间计算
在建筑行业，墙壁隔声量较大，跨越空间的噪声源较少，可单独对每个房间做声学优化设计。声线跟踪法在声场计算时只考虑声线反射、散射和吸收损失，计算一个封闭空间内的声场分布情况。而在船舶行业，舱室之间以舱壁隔断，壁板较薄，声透射量较大，噪声源通常影响多个舱室，所以采用声线跟踪法的基本原理，在声线传播计算时，考虑声线透射的情况，其透射计算遵循薄板隔声原理。假设当第n条声线传播至壁板m时，发生透2射、吸声和反射（散射）作用，其路线如图 所示。透射声线沿原方向传播，能量减少 倍（ 为壁板m的透射系数）。同时产生一条新的声线 Inm ，其能量为原声线的 (1 - 倍。其中 为壁板m
的吸声系数，吸声系数随声线频率的不同而变化。 由于船舶舱壁一般为薄板，其透射符合“质量定律”及“吻合效应”。当 sin θ= c 0/c 时，产生“吻b合效应”，其吻合频率为［13］
基于上述计算原理，即可实现跨越舱壁、多空间的声场模拟计算。
２ 声线搜索法
用声线跟踪法描述声能传播的方式简便，易于计算，故根据其计算原理，提出了舱室噪声的优化算法——声线搜索法。在声线跟踪法的计算原理中，声能传播以声线传播模拟，声线传播路线即是声能传递路径，经过目标舱室的声线数量越多，
携带能量越大，舱室噪声便越大。反之，则越少。在这部分对目标舱室有贡献的声线的传播过程中，若减少这部分声线携带的能量，就可有效降低舱室噪声值，从而达到降噪的目的。声线搜索法就是搜索这部分“有效”声线经过的位置并记录，计算舱壁不同位置对这部分声线的影响，计算舱壁对目标舱室噪声的声学灵敏度，并根据灵敏度大小制定舱室隔声降噪最优方案。该方法的流程如下： 1）遍历跟踪由声源处发射的全部声线，记录通过目标舱室的声线编号、声线能量变化以及传递路径，计算目标舱室声强I。2）分析声线经过的舱壁位置，并按舱壁位置分组。同一条声线可经过不同的舱室，因此，同一条声线可同时分布在多个组别。3）非相干叠加计算每组声线到达目标舱室时的声强之和，第n组声强之和为 In = Ik 。4 ）计算舱壁对目标舱室声压影响的灵敏度。第n组声学灵敏度为
(1I- 10 lg ¶p p0 -p I - In 1/IL) （15） = = IL·¶S IL·S IL·S式中：p0 为没有插入损失的目标舱室的声压；p为有插入损失的目标舱室的声压；IL为舱壁插入损失或引入降噪措施的隔声量；S 为单位面积。5）根据灵敏度给出降噪方案，计算降噪效果。
3 声场计算及噪声优化
3.1 声线跟踪法在多空间声场计算中的准确性及收敛性分析
声线跟踪法是声线搜索法的前期计算，是声线搜索法的理论基础，其准确性是声线搜索法这一噪声优化算法有效的保证。声线跟踪法不考虑声的波动性，适用于中、高频噪声，因此采用统计能量法进行对比分析。算例采用船舶中典型的舱1 2室结构，由相邻的声源舱 和受声舱 构成。声腔长3m 3.5 m，高3m 5mm ，宽 ，舱壁采用 厚钢板。声1号声腔（0.1，0.2，0.5）处，为 105.6 dB源位于 白噪31.5~声（计算时换算成声功率），计算频率为8 000 Hz 2 3所，目标舱室为 号声腔。结构如图=344 m/s，钢示。空气密度 ρ0 =1.29 kg/m3，声速 c0 =7 800 =0.3，阻尼损耗密度为 ρ kg/m3，泊松比 μ st =0.3，弹性模量 E=2×1011 N/m2。各频率的因子 η 1.2隔声量采用 节中方法计算得到，吸声系数在全0.1。此算例中，A，B频段取 为一块钢板，不做区分。 31.5~8 000 Hz采用声线跟踪法数值仿真计算2 4，0.5，1.2），频段内 号声腔的声压级，具体为（ （4，3，1.2），（4.5，1.7，1.5），（5.5，0.5，1.8），（5.5，3， 1.8）这 5 7 200，点的声压级均值，初始声线数为EDP 0.001，使用统计能量法软件进行计算比对为 4分析，计算结果如图 所示。 4可见，2由图 种算法的变化趋势一致，声压2 500 Hz级随频率的增加而减小，在约 处出现了2.85 dB，相对误差吸收峰值。其声压级平均误差6%，其中误差在低频处和2 000~2 600 Hz处差距较大。在低频阶段，噪声波长大于舱室尺寸，声线法计算结果误差较大，统计能量法也仅适用于中、2高频噪声，因此这 种方法均不适用于低频段，不能准确计算舱室声压级。此结构钢板的临界频率=2 458.2 Hz，发生“吻合效应”，板透声较大，为 fcr 2 500 Hz 2 000~因此在 处出现了声压级峰值。在2 600 Hz 2处，这 种方法计算差距较大，其原因是在临界频率附近，板隔声量计算误差较大。若板隔声量以实际实验数据计算，则会增加其计算准确率。声线跟踪法是以带能量的声线模拟声波传递，初始声线数越大，其声线分布密度越大，每条声线所携带的能量越小，越符合实际情况。然而受计算量的限制，初始声线数不能无限大，为保证结果的有效性，选取合适的初始声线数，本文选取
500，630 ，800，1 000，1 600，3 150，5 000 8 000 Hz和8 1 350，1 800，这 个频率，分别计算初始声线数为2 592，4 050，7 200，11 250时目标舱室的噪声值， 5其结果如图 所示。 根据收敛性计算结果，不同声线数的计算结果在定值附近波动，频率越低波动越大，频率越高越趋于稳定。根据收敛计算结果，可以认为初始7 200声线数取 时，声线反射和透射充分，计算结果有效。
3.2 声线搜索法优化计算
2 1m 3个相连声腔的长宽高均为 ，结构如图所示。B 10 dB），A 0.7 mm处为门（隔声量取 为 厚y=0.5钢板，两者分界线位于 处，采用声线搜索法500 Hz A，B计算频率为 时 两处对目标舱室降噪1的灵敏度，结果如表 所示。 1 B由表 的计算结果可知，舱壁 的声学灵敏4.209 2×10-2，大 A度为 于舱壁 处的声学灵敏度， B根据声线搜索法理论，在舱壁 处引入降噪措施， 2如添加隔声材料，对舱室 的降噪效果最佳。A，B 10 dB，分别在 处引入降噪措施，插入损失2 2舱室 处的声压计算结果如表 所示。2 B 10 dB根据表 的计算结果，在舱壁 处添加2 7.283 dB，而的插入损失，舱室 的噪声下降 在舱A 10 dB 2室 处添加 的插入损失，舱室 的声压下降0.191 1 dB。在舱壁B仅为 处引入降噪措施远好A 1、表 2于 处。两者（表 ）结果相符，采用声线搜索法可以计算舱室对目标舱室的降噪灵敏度，方
4结语
考虑透射情况的声线跟踪法可用于计算多空间声场，适用于中、高频情况，该方法可计算声场中某点或某区域的声压级情况，弥补了统计能量法只能计算子系统平均声压级的不足。与有限元法相比，其在中、高频的计算更加准确。依据声线跟踪法的计算过程，提出的声线搜索法可直接、快速计算各个位置的舱壁对目标舱室噪声的灵敏度，计算结果表明，声线搜索法有效，可用于指导舱室中、高频段噪声优化设计。该方法与改变结构参数、迭代计算模型来求得灵敏度相比，用时大幅减少。同时，该方法还可以计算多声源时多舱壁对目标舱室的灵敏度，在复杂结构中有着良好的应用前景。
参考文献：
［1］ KOO B U. Shape design sensitivity analysis of acoustic problems using a boundary element method［J］. Com⁃ puters & Structures，1997，65（5）：713-719. 2 LE BOT A，CARCATERRA A，MAZUYER D. Statisti⁃ ［ ］ cal vibroacoustics and entropy concept［J］. Entropy， 2010，12（12）：2418-2435. 3 高处，杨德庆. ［］船舶舱室噪声传递路径分析的声振
熵赋权图法［J］. 上海交通大学学报，2014，48（4）： 469-474. GAO C，YANG D Q. Vibroacoustical entropy weighted graph method for sound transmission path analysis of ship cabin noise［J］. Journal of Shanghai Jiaotong Uni⁃ versity，2014，48（4）：469-474（in Chinese）. 4 薛晓理. 薄板结构的声功率及其灵敏度研究［D］. ［］ 南昌：南昌航空大学，2014. XUE X L. Analysis of the sound power and its sensitivi⁃ ty of thin plate structure［D］. Nanchang：Nanchang Hangkong University，2014（in Chinese）. 5］ 彭伟才，何锃. ［ 基于统计能量模型的辐射声功率灵
敏度分析［J］. 中国舰船研究，2009，4（6）：11-14. PENG W C，HE Z. Radiated sound power sensitivity
54 analysis based on statistical energy analysis model［J］. Chinese Journal of Ship Research， 2009， 4（6）： 11-14（in Chinese）. ［6］ BOROBIA A，MORO J，SOTO R. A unified viewon compensation criteria in the real nonnegative inverse ei⁃ genvalue problem［J］. Linear Algebra and its Applica⁃ tions，2008，428（11/12）：2574-2584. 7］ . 40 年（1968～ ［ 曾向阳 室内声场计算机模拟发展2008）［J］. 声电技术，2008（增刊1）：12-17，23. ZENG X Y. Forty years of development of room acous⁃ tic computer simulation［J］. Audio Engineering，2008 （Supp 1）：12-17，23（in Chinese）. 8］ 蔡铭，陈志斌，邹竞芳. ［ 一种结合室内空间剖分的声线跟踪法［J］. 应用声学，2013，32（1）：34-38. CAI M，CHEN Z B，ZOU J F. A ray tracing method based on indoor space partitioning［J］. Applied Acous⁃ tics，2013，32（1）：34-38（in Chinese）. 9］ 蒋忠进，邱小军. ［ 基于声线跟踪法的稳态声场空间算［J］. 2007，分布计 南京大学学报（自然科学版）， 43（1）：56-65. JIANG Z J，QIU X J. Spatial distribution imaging of steady-state sound field with the ray-tracing method ［J］. Journal of Nanjing University（Natural Sciences）， 2007，43（1）：56-65（in Chinese）. 10］ 曾向阳，王海涛，杜博凯. ［ 基于声粒子分布积分的J］. 2016，35无网格声场计算方法［ 应用声学， （1）：84-89. ZENG X Y，WANG H T，DU B K. A meshless meth⁃ od based on the particle-distribution integration for simulationsof sound fields in small enclosures［J］. 40 （上接第 页） ［9 . 船舶舱室声学优化设计理论与方法研究［J］. ］ 高处上海：上海交通大学，2015. GAO C. Optimization theory and methods of ship cab⁃ ins acoustics desigh［D］. Shanghai：Shanghai Jiao Tong University，2015（in Chinese）. 10］ 高处，杨德庆. ［ 船舶舱室噪声传递路径分析的声振J］. 2014，48熵赋权图法［ 上海交通大学学报，
12 Journal of Applied Acoustics，2016，35（1）：84-89 （in Chinese）. ［11］ YANG L N，SHIELD B M. The prediction of speech intelligibility in underground stations of rectangular cross section［J］. Journal of the Acoustics Society of America，2001，109（1）：266-273. 12］ 李久会，袁泉，吴中立，等. ［ 室内声场中点声源的
J］.算法与模拟［ 辽宁工业大学学报（自然科学版），2009，29（3）：208-210. LIJH ，YUAN Q， WUZL ，et al. Algorithm and sim⁃ ulation of point source in acoustic fields［J］. Journal of Liaoning University of Technology（Natural Science Edition），2009，29（3）：208-210（in Chinese）. 13］ 黄玉盈. 结构振动分析基础［M］. ［ 武汉：华中理工大学出版社，1988. ［14］ SEWELL E C. Transmission of reverberant sound through a single-leaf partition surrounded by an infi⁃ nite rigid baffle［J］. Journal of Sound and Vibration， 1970，12（1）：21-32. 15］ 陆红艳. ［ 板结构的振动声辐射及其隔声性能研究［D］.武汉：武汉理工大学，2003. LU H Y. The vibration & acoustic radiation and noise muffling characteristics of plate structures［D］. Wu⁃ han：Wuhan University of Technology，2003（in Chi⁃ nese）. ［16］ WARNOCKAC C，FASOLD W. Sound insulation： airborne and impac［t M］//CROCKER M J. Encyclope⁃ dia of Acoustics. New York：Wiley-Interscience Pub⁃ lication，1997：1129-1161. （4）：469-474. GAO C，YANG D Q. Vibroacoustical entropy weight⁃ ed graph method for sound transmission path analysis of ship cabin noise［J］. Journal of Shanghai Jiaotong University，2014，48（4）：469-474（in Chinese）. ［11］ WALL M. GAlib—a C ++ library of genetic algorithm components［EB/OL］.（1996-08-01）［2017-03-01］. http：//lancet.mit.edu/ga/.