ProNas能量有限元法在船舶中高频噪声预测的应用
本文，结合具体应用，介绍ProNa能量有限元法在船舶中高频噪声预测中的应用。以ProNas能量有限元理论为基础，建立了船舶的ProNas能量有限元计算模型，采用安世亚太大型商用软件ProNas对复杂激励在船舶各舱室产生的中高频结构噪声及空气噪声进行仿真计算，得到船舶各舱室声压级，并利用ProNas软件后处理功能确定激励源及传递路径处的能量分布云图，对不满足噪声目标的舱室进行了声学优化，最终解决了大型实际船舶工程的中高频噪声预测与控制问题。
复杂结构的中高频噪声的控制一直以来都是各工业领域研究的重点与难点问题，尤其对于大型船舶其内环境相比其它工业产品更加独特：结构形式纵横交错，舱室众多，噪声诱因复杂，声源品种繁多密集，噪声强度较大，噪声频域带宽且持续稳定，结构噪声与空气噪声相互转化。以上这些特点，就使得船舶噪声控制起来更加困难。
2014年7月国际海事组织（IMO）签订生效的新的《船上噪声等级规则》要求居住区部分舱室声压级在旧规范的基础上降低5dB(A),这就要求船舶工程设计人员需要采取更加有效的控制手段来降低船舶噪声。
传统的以有限元（FEA）、边界元（BEA）、统计能量分析（SEA）等算法为基础而发展起来的商用软件工具，在计算效益上存在不足和瓶颈，很难满足来自噪声振动工程界及学术科研的越来越复杂、精细及多学科综合解析优化的工程设计和技术发展需求。
能量有限元-统计能量混合模块（ProNAS/EFEA-SEA）涉及的基本变量是平均的能量或能量密度，既可以快速建模预测及优化振噪特性，也可以直接使用现有普通有限元网格进行中高频振动的分析和模拟，从而大大节省工程设计人员的建模时间，使工程人员在设计初期能够有效地进行工程预测和优化。
一、船舶模型建立
1、船舶前处理模型
本文以某客箱船为例。利用Patran前处理软件对某客箱船进行有限元建模，如图1所示，该模型共有单元114358个，节点94357个。该客箱船全长约200米，船宽约25米，采用的是双机双浆配置，主机为两台低速机。该船的前部为载货区域，后部为居住区域和机械处所。因对居住区域及机械处所噪声控制要求较高，同时为了节省计算时间，本案例截取局部模型，即仅对居住区域和机械处所进行研究。
图1.某客箱船有限元模型
2、ProNas软件模型处理
通过ProNas声学仿真软件分别对船舶的结构噪声与空气噪声进行了仿真计算。将Nastran格式模型导入ProNas软件中，如图2所示。通过ProNas软件声腔自动识别功能自动划分声腔子系统，如图3所示。软件自动耦合声腔子系统和其临近贴合的结构单元，这样能量即可在结构单元与声腔子系统之间进行交换。然后对各舱室结构单元进行声学材料设置，如甲板敷料、绝缘布置材料，设置门窗材料等，如图4所示。
图2.某客箱船ProNas有限元模型
图3.某客箱船各舱室声腔子系统
图4.某客箱船声学材料施加
二、船舶模型激励输入
船舶噪声主要包括主机噪声、螺旋桨噪声及水动力噪声。其中主机噪声是船舶噪声中最强的噪声源，主机噪声主要分为空气辐射噪声、排气噪声和结构噪声，空气辐射噪声有燃烧噪声和机械往复运动碰撞噪声，成因复杂，一般实测得到。
螺旋桨噪声也是船舶的主要噪声源，螺旋桨直接产生的噪声有空泡噪声，空泡噪声会引起船体结构噪声。水动力噪声主要是由于高速海流的不规则起伏作用于船体湿表面，激起船体的局部振动并向周围介质(空气与水)辐射的噪声。
该客箱船模型根据实际情况所受激励载荷主要有24个，分别为主机、螺旋桨、发电机、空压机、空调、风机等。其中主机、发电机、螺旋桨与船体作用，船体会振动产生结构噪声。
其它主要激励源通过空气传递到船舶舱室内，产生空气噪声。结构噪声源通过ProNas软件分别在各设备的有限单元表面以速度级及加速度级方式加载，如图5所示。空气噪声源通过ProNas软件分别在各激励源所在声腔子系统内以声功率级方式加载，如图6所示。主机速度级激励如图7所示，螺旋桨加速度级激励如图8所示。
三、船舶结构噪声及空气噪声分析结果与优化
1、船舶噪声分析结果
在倍频程中心频率63-8000Hz下分别将该客箱船结构噪声模型与空气噪声模型提交ProNas软件求解器，经过ProNas软件计算，可得到各舱室声压级水平，图9、图10分别为该客箱船在中心频率为500Hz时结构噪声仿真得到的结构速度云图与舱室声压云图，图11与图12分别为该客箱船在中心频率为500Hz时的空气噪声仿真得到的结构速度云图和舱室声压云图。
该客箱船舱室声压分析结果见表1，部分舱室声压级不满足目标值。
2、优化方案
对于结构噪声超标的舱室，常见的优化方法为敷设阻尼。阻尼材料是将结构振动板的振动能量快速转化为热能，从而减弱金属板的弯曲振动，阻尼材料通过这种方式有效的控制金属板的结构辐射噪声。
通过云图分析可知，本案例结构噪声超标，多为主机与螺旋桨结构噪声引起。由于该客箱船采用低速机，主机与船体钢板直接螺接，船体底部振动区域较大，且底部钢板厚度约为25mm-30mm，在实际应用中，阻尼层厚度一般为金属板厚度2-4倍，因此如果采用常规的在激励源处阻尼敷设方法，阻尼用量相对较大，成本较高。
又由于船舶结构复杂，结构噪声和空气噪声相互转化，因此常规阻尼敷设方法对超标舱室噪声控制效果一般。针对以上特点，本案例通过ProNas软件后处理界面可直观精确显示传递路径处能量分布，即在不达标舱室处直接敷设阻尼，采用这种敷设方法，阻尼用量及成本在可控范围内，且噪声控制效果显著。图13为常规敷设阻尼方法，图14为本案例阻尼敷设方法。
对于空气噪声不达标舱室，结合ProNas软件后处理结果云图，本案例采用以下优化方案：
部分娱乐室及放映室等地面浮动地板均改为50mm岩棉+钢板+阻尼+钢板型，再在上侧加10mm流平甲板敷料。
部分甲板厨房抽风机室的��部，艉、中、艏面增加4mm钢板做成双墙舱壁，形成密封夹层空间，两层钢壁都设阻尼涂料。
甲板房间合计26间处所的天花板需穿孔（穿孔面向室内），孔径按1.8mm。壁板**孔。
在部分甲板机舱棚两侧直接面向公共处所和房间的舱壁，在外侧增加4mm压筋板做成的舱壁，形成密封的夹层空间，夹层内设50mm隔音棉。
部分甲板均设9mm厚阻尼涂料处理。
3、实船仿真与测试结果对比
经过以上方案优化并多轮仿真计算，优化以后各舱室噪声控制效果较好，将以上优化方案应用到实船设计中，并对该客箱船实船舱室进行了噪声测试，如表2所示，大多舱室噪声数值满足设计要求。
四、结论
本文基于能量有限元方法的ProNas软件建立了某客箱船声学模型，进行了整船的结构噪声与空气噪声仿真计算，对结构噪声及空气噪声不达标舱室进行仿真优化，并与实船测试结果进行对比，得出了以下结论：
ProNas能量有限元方法是一种用于预测和解决船舶及其它工业产品中高频噪声的可行性及有效性方法，具有很大的理论意义与现实价值。
采用ProNas能量有限元方法，可避免常规中高频方法子系统建模操作的复杂性与专业性，提高仿真计算精度。
ProNas能量有限元方法不需要判断子系统模态密度，不需要确定系统之间的耦合因子，既保证模型的完整性，又提高仿真的时效性。
ProNas软件为解决中高频噪声提供了新的方法，为改进船舶设计参数提供了科学依据，提高了船舶的噪声设计质量，缩短研发周期，节约开发成本。
ProNas软件仿真得到的船舶结果与实测结果的最大误差值的绝对值在5dB(A)内，满足模型的预测精度。
