ProNas能量有限元方法在船舶中高频噪声预测的应用.docx

复杂结构的中高频噪声的控制一直以来都是各工业领域研究的重点与难点问题，尤其对于大型船舶其内环境相比其它工业产品更加独特口:结构形式纵横交错，舱室众多，噪声诱因复杂，声源品种繁多密集，噪声强度较大，噪声频域带宽且持续稳定，结构噪声与空气噪声相互转化。以上这些特点，就使得船舶噪声控制起来更加困难。2014年7月国际海事组织(IMO)签订生效的新的船上噪声等级规则要求居住区部分舱室声压级在旧规范的基础上降低5dB（八）,这就要求船舶工程设计人员需要采取更加有效的控制手段来降低船舶噪声。目前预测结构噪声及空气噪声的理论体系相对完善，并已将理论应用于大量实际工程中。按激励源频率及具体工业产品，可将振动和噪声的问题划分为低频、中频及高频问题。低频结构的响应具有随机和确定性，且为全局模态，工程中常用的数值方法为有限元(FEM)边界元(BEM)等。理论上，有限元方法可计算任意结构的振动场，但是为了能准确反映结构的振动特性，划分单元的网格长度要求要远低于结构中弯曲波的波长。这样随着频率的增加，结构弯曲波波长变小，结构网格需要进一步细化，对计算资源要求较高，且结构模态密集，重叠交错，响应表现出不确定性，具有了统计的概念，因此，有限元和边界元就不再适合解决中高频问题。统计能量分析(SEA)用于解决中高频问题，且模态越密集，其计算精度就越高。但统计能量分析不能保证子系统的空间变量信息的完整性2,难以精确预示子系统内能量密度分布且子系统的划分需要一定的经验，不易进行实际结构形态的设计与优化4,模态密度及耦合因子的准确与否直接影响结果的准确度5-6。所以为了更好的控制中高频噪声，就需要有更好的理论方法。能量有限元法7(EFEA)是一种预测中高频动响应的新方法，它是以波动理论为基础，将结构离散化，在单元之间建立能量密度关系式，从而求解得到所有节点的能量密度。能量有限元法以有限单元为对象，保证了模型的完整性，但是在结构突变处引入了大量重复节点，求解效率得不到平衡，能量有限元理论目前可以解决一些简单的有限元噪声问题，对于复杂问题的应用依然比较困难。ProNas能量有限元是在统计能量分析及能量有限元理论的基础上，以有限单元为研究对象，混合了SEA及EFEA理论，利用有限体积法及差分法推导出得类似于SEA的理论方程，联立求得每个有限单元的能量密度。PrONaS能量有限元法避免了在结构突变处引入重复节点，保证了求解精度和计算效率，是一种具有较高研究价值及应用价值的中高频噪声控制方法。本文着重介绍了Pr。NaS能量有限元理论原理且应用ProNaS商业软件对大型实际船舶的中高频噪声进行了仿真计算，求解得出船舶各舱室声压级。一、ProNas能量有限元基本理论ProNas能量有限元法可有效解决结构振动系统中高频噪声问题，它包含了能量有限元法(EFEA)及统计能量分析法(SEA)的优点，具有更大优势。ProNas能量有限元方法是在波动理论上建立的一种功率流方法，通过波动形式求解结构微元体运动方程，并将波动形式下的动能密度、势能密度及功率流进行周期内时间平均和局部空间平均，得到能量密度与功率流的关系，代入稳态下能量平衡方程中可以得到能量密度控制方程；采用一定的数值离散方法对控制方程进行离散，得到能量有限元方程；在结构存在耦合的地方考虑波的反射和折射，求出耦合处能量密度与功率流的关系；最后进行总系数矩阵的组集，得到总体的方程，从而求出能量密度。有限体积法的基础是将模型离散成若干个控制体，而这里的控制体可理解为有限元的“单元1.这样，ProNas能量有限元就是以单元为研究对象。弹性介质中的能量平衡方程网如下：式中，e为能量密度；为结构的应力张量；疔为位移向量；11in为输入功率;11diss为时间和空间的平均耗散功率。稳定状态下，子系统功率平衡，可以表示为：V+=7s声能量密度和耗散功率之间的关系如下：吗”=于加声能量密度和声强之间的关系如下：了=-Ve因此，可以得到能量密度方程：c（32e_式中，n为阻尼损耗因子；Cg为群速度，为圆频率。例如对于简单三段梁模型如图1所示。1 w2eQ卜加十.T一T图1简单三段梁示意图能量有限元方程如下：（%+叽）0_琮2必0jr£ir-限0,+/】+猿）。-脸0,/0-嬉。（+啕闯尻3一般性地，ProNaS能量有限元矩阵方程可以表述为:KE=11ii式中，K为总体能量密度系数矩阵；E为节点能量密度向量；11in为总体输入功率向量，求解该式就可以得到单元的能量密度。由上式可知，ProNaS能量有限元矩阵方程与SEA方程形式一致，为一次线性方程，易于求解。二、船舶模型建立1.船舶前处理模型本文以某客箱船为例。利用Patran前处理软件对某客箱船进行有限元建模，如图2所示，该模型共有单元114358个，节点94357个。该客箱船全长约200米，船宽约25米，采用的是双机双浆配置，主机为两台低速机。该船的前部为载货区域，后部为居住区域和机械处所。因对居住区域及机械处所噪声控制要求较高，同时为了节省计算时间，本案例截取局部模型，即仅对居住区域和机械处所进行研究。图2某客箱船有限元模型2.ProNas软件模型处理ProNas是拥有国内自主软件著作权的自主可控的中高频软件，是能量有限元分析(EFEA)和统计能量分析(SEA)领域的代表性解决方案，代表着振动噪声工程界新一代的前沿技术。本案例通过ProNas声学仿真软件分别对船舶的结构噪声与空气噪声进行了仿真计算。将NaStran格式模型导入ProNaS软件中，如图3所示。通过ProNaS软件声腔自动识别功能自动划分声腔子系统，如图4所示。软件自动耦合声腔子系统和其临近贴合的结构单元，这样能量即可在结构单元与声腔子系统之间进行交换。然后对各舱室结构单元进行声学材料设置，如甲板敷料、绝缘布置材料，设置门窗材料等，如图5所示。图3某客箱船ProNas有限元模型图4某客箱船各舱室声腔子系统图5某客箱船声学材料施加三、船舶模型激励输入船舶噪声9-10主要包括主机噪声、螺旋桨噪声及水动力噪声。其中主机噪声是船舶噪声中最强的噪声源，主机噪声主要分为空气辐射噪声、排气噪声和结构噪声，空气辐射噪声有燃烧噪声和机械往复运动碰撞噪声，成因复杂，一般实测得到。螺旋桨噪声也是船舶的主要噪声源,螺旋桨直接产生的噪声有空泡噪声,空泡噪声会引起船体结构噪声。水动力噪声主要是由于高速海流的不规则起伏作用于船体湿表面，激起船体的局部振动并向周围介质（空气与水）辐射的噪声。该客箱船模型根据实际情况所受激励载荷主要有24个，分别为主机、螺旋桨、发电机、空压机、空调、风机等。其中主机、发电机、螺旋桨与船体作用，船体会振动产生结构噪声。其它主要激励源通过空气传递到船舶舱室内，产生空气噪声。结构噪声源通过ProNas软件分别在各设备的有限单元表面以速度级及加速度级方式加载，如图6所示。空气噪声源通过PrONaS软件分别在各激励源所在声腔子系统内以声功率级方式加载，如图7所示。主机速度级激励如图8所示，螺旋桨加速度级激励如图9所示。图6某客箱船结构噪声载荷分布图图7某客箱船空气噪声载荷分布图主机SMCR40.020.0：=0.01IIM-1010010100Hz图8主机激励频谱图螺旋桨180160-WOooooooo1412108C&4C2加速度级10010100图9螺旋桨激励频谱四、船舶结构噪声及空气噪声分析结果与优化1 .船舶噪声分析结果在倍频程中心频率63-8000Hz下分别将该客箱船结构噪声模型与空气噪声模型提交ProNas软件求解器，经过ProNas软件计算，可得到各舱室声压级水平,图10、图11分别为该客箱船在中心频率为500HZ时结构噪声仿真得到的结构速度云图与舱室声压云图，图12与图13分别为该客箱船在中心频率为500Hz时的空气噪声仿真得到的结构速度云图和舱室声压云图。图10.结构速度云图（50OHZ.结构噪声）图11.舱室声压云图（500HZ-结构噪声）图12.结构速度云图（500HZ-空气噪声）图13.舱室声压云图（500HZ-空气噪声）该客箱船舱室声压分析结果见表1,部分舱室声压级不满足目标值。表1某客箱船舱室声压级（节选）96SA计权向历«dB（八）目除肉压IRdB（八）V63Hz125H2250HZSOoKZ100OHz2000H2400OHZ8000Hxroom-7138现77737S.S64.349.2387ss2669.2SS不8room-7137越977275.566.851.$38.8J26369.655TBroom713538.977.175wS663SUM751322.669.555三sroom-7134397727”66SU35.148317.3W.555room-7133M97775uS66351.1354s19694SSroom-7132417674.66S249.838.35Z723.868.S55room-811240.58889.882.665J41.349.32a483.855不禽aroom-S111408889.8S2.66Su843.3S4.3居883355room-81104087.989.782.665u741.75Q922983.755不禽Sroom*81094189.189.782.66M44457aas83<SSSroom-81M40.5899M93.1802S426313914SSiSroom-81334189.3928½17851359.735290.555不戏S2 .优化方案对于结构噪声超标的舱室，常见的优化方法为敷设阻尼。阻尼材料是将结构振动板的振动能量快速转化为热能，从而减弱金属板的弯曲振动，阻尼材料通过这种方式有效的控制金属板的结构辐射噪声。通过云图分析可知，本案例结构噪声超标，多为主机与螺旋桨结构噪声引起。由于该客箱船采用低速机，主机与船体钢板直接螺接，船体底部振动区域较大，且底部钢板厚度约为25mm-30mm,在实际应用中，阻尼层厚度一般为金属板厚度2-4倍，因此如果采用常规的在激励源处阻尼敷设方法，阻尼用量相对较大，成本较高。又由于船舶结构复杂，结构噪声和空气噪声相互转化，因此常规阻尼敷设方法对超标舱室噪声控制效果一般。针对以上特点，本案例通过Pr。NaS软件后处理界面可直观精确显示传递路径处能量分布，即在不达标舱室处直接敷设阻尼，采用这种敷设方法，阻尼用量及成本在可控范围内，且噪声控制效果显著。图14为常规敷设阻尼方法，图15为本案例阻尼敷设方法。图14常规敷设阻尼区域图15本案例敷设阻尼区域对于空气噪声不达标舱室，结合ProNaS软件后处理结果云图，本案例采用以下优化方案：第一、部分娱乐室及放映室等地面浮动地板均改为50mm岩棉+钢板+阻尼+钢板型，再在上侧加IOmm流平甲板敷料。第二、部分甲板厨房抽风机室的内部，腥、中、艄面增加4mm钢板做成双墙舱壁，形成密封夹层空间，两层钢壁都设阻尼涂料。第三、甲板房间合计26间处所的天花板需穿孔（穿孔面向室内），孔径按1.8mmo壁板不穿孔。第四、在部分甲板机舱棚两侧直接面向公共处所和房间的舱壁，在外侧增加4mm压筋板做成的舱壁，形成密封的夹层空间，夹层内设50mm隔音棉。第五、部分甲板均设9mm厚阻尼涂料处理。3 .实船仿真与测试结果对比经过以上方案优化并多轮仿真计算，优化以后各舱室噪声控制效果较好，将以上优化方案应用到实船设计中，并对该客箱船实船舱室进行了噪声测试，如表2所示，大多舱室噪声数值满足设计要求。4 2部分舱室仿真结果与测试结果（节选）航室港试铝累(dB（八）)仿真结果(dB（八）)仿真与流试差(dB)Room-811151S3.121Room-81125151.80.8Room-81334850.323Room-71335452.6-1.4Room-7134S454.70.7Room-71385357.54.5五、结论本文简要介绍了ProNas能量有限元方法，在此基础上应用ProNas软件建立了某客箱船声学模型，进行了整船的结构噪声与空气噪声仿真计算，对结构噪声及空气噪声不达标舱室进行仿真优化，并与实船测试结果进行对比，得出了以下结论：第一、ProNas能量有限元方法是一种用于预测和解决船舶及其它工业产品中高频噪声的可行性及有效性方法，具有很大的理论意义与现实价值。第二、采用ProNas能量有限元方法，可避免常规中高频方法子系统建模操作的复杂性与专业性，提高仿真计算精度。第三、ProNas能量有限元方法不需要判断子系统模态密度，不需要确定系统之间的耦合因子，既保证模型的完整性，又提高仿真的时效性。第四、ProNas软件为解决中高频噪声提供了新的方法，为改进船舶设计参数提供了科学依据，提高了船舶的噪声设计质量，缩短研发周期，节约开发成本。第五、ProNas软件仿真得到的船舶结果与实测结果的最大误差值的绝对值在5dB（八）内，满足模型的预测精度。附参考资料：能量有限元软件ProNaS综述如何有效解决中高频噪声问题目前仍是学术界和工程应用领域的难题之一。在当前解决中高频噪声的几种主要理论方法中，ProNas能量有限元方法作为一种全新的可行有效的中高频噪声控制理论，具有较强的理论和应用价值。ProNaS能量有限元方法克服了统计能量分析和能量有限元方法的不足之处,可用于求解强耦合、大阻尼等非保守系统，在降低工程应用人员的操作难度，缩短产品开发周期等方面都表现了极大的优势；并且，其核心算法，保证了仿真的精度与求解效率。在中高频噪声控制领域，ProNaS能量有限元方法很值得期待。基于ProNas能量有限元方法，安世亚太联合国际最先进的中高频专家资源共同开发了拥有国内自主软件著作权的中高频噪声仿真分析软件ProNas,助力解决中高频噪声控制难题。作为振动噪声工程界新一代前沿技术的代表，ProNas成功的破解了传统中高频方法面临的困境。ProNas能量有限元方法产生的背景当前，解决中高频噪声有几种主要理论方法：统计能量分析方法、能量有限元方法及ProNas能量有限元方法。统计能量分析是一种用于较宽频率范围内的随机噪声的统计方法。但统计能量分析的应用有大量前提假设，且统计能量分析不能保证子系统的空间变量信息的完整性,子系统的划分需要一定的经验，不易进行实际结构形态的设计与优化。在这样的背景下，能量有限元方法产生了。能量有限元方法以波动理论为基础，将结构离散，在单元之间建立能量密度的关系方程，求解得到结构上所有点的能量密度响应信息。能量有限元方法在结构突变处引入大量重复节点及能量密度，计算效率上得不到平衡。而且，目前的能量有限元方法在结构振动及声辐射问题上的应用还局限在简单的问题上。ProNas能量有限元方法应需而生。ProNas能量有限元方法在统计能量分析及能量有限元理论的基础上，将单元离散化，以有限单元为研究对象，并混合了SEA及EFEA理论，利用有限体积法及差分法推导出得类似于SEA的理论方程。ProNas能量有限元法避免了在结构不连续处引入重复节点，保证了求解精度和计算效率，是一种具有较高研究价值及工程应用价值的中高频噪声控制方法。ProNas能量有限元方法的原理ProNas能量有限元方法是在波动理论上建立的一种功率流方法，通过波动形式求解结构微元体运动方程，并将波动形式下的动能密度、势能密度及功率流进行周期内时间平均和局部空间平均，得到能量密度与功率流的关系，代入稳态下能量平衡方程中可以得到能量密度控制方程；采用一定的数值离散方法对控制方程进行离散，得到能量有限元方程；在结构存在耦合的地方考虑波的反射和折射，求出耦合处能量密度与功率流的关系；最后进行总系数矩阵的组集，得到总体的方程，从而求出能量密度。这里的有限体积法是将模型离散成若干个控制体，而控制体可理解为有限元的“单元”，这样，PrONaS能量有限元就是以单元为研究对象。弹性介质中的能量平衡方程如下：式中，为能量密度；为结构的应力张量；为位移向量；为输入功率；为时间和空间的平均耗散功率。稳定状态下，子系统功率平衡，可以表示为：声能量密度和耗散功率之间的关系如下:兀dM=7l8e声能量密度和声强之间的关系如下：I=-(cj)Vc因此，可以得到能量密度方程：式中，为阻尼损耗因子；为群速度，为圆频率。例如对于简单三段梁模型如图5所示。X图5.简单三段梁示意图ProNas能量有限元方程如下:-fJl0jr（小+埔21+哈）2-错QE2-能划（小+嘘）司与一般性地，PrONaS能量有限元矩阵方程可以表述为:KE=Q式中，K为总体能量密度系数矩阵；E为单元能量密度向量；为总体输入功率向量，求解该式就可以得到单元的能量密度。由上式可知，PrONaS能量有限元矩阵方程与SEA方程形式一致，为一次线性方程，易于求解。ProNas能量有限元方法的比较优势ProNas能量有限元方法可以预测整个系统在结构激励与声场激励下的中高频噪声，以有限单元为载体，具有以下先天性的优势：免去冗长子系统的建模过程，可直接利用强度或刚度有限元模型进行高频噪声计算；避"反复判断子系统模态密度，确保模型细节的完整性；免去了计算、测量和仿真结构阻抗的过程，可更精确的计算中高频结构噪声；没有大量前提假设，有更宽范的阻尼及耦合强度适用范围；降低了工程应用人员操作门槛，不需要工程人员具备专门知识；更高的求解效率及更精确的计算结果。以下为两个应用ProNas能量有限元理论计算的简单案例：案例一：分别采用SEA理论及ProNas能量有限元方法计算均质梁在简谐激励作用下的能量密度。由计算结果可知，SEA理论和ProNas能量有限元理论均可精确计算均质梁能量密度，但SEA理论对划分子系统的把握要求较高，而PrONaS能量有限元理论，对系统划分网格越细，求解精度越高。案例二：采用SEA理论与ProNas能量有限元理论去计算非连续梁强耦合工况的能量密度。图4非连续梁示意图-A。收的介瞬K恻f比葭:.“恍力&:一,方2扃*td万雄：ProNASr:”上MitattA<FV>：fttitttiA<A*VtftitV»图5SEA理论与ProNas能量有限元理论计算结果由上图计算结果可知，SEA理论在计算强耦合问题时，计算结果会有较大偏差如图5(b)所示，而ProNAS能量有限元理论计算结果与解析结果完全一致如图5（八）所示。因此，ProNAS能量有限元理论的产生为解决强耦合问题提供了方法，具有较强的理论优势。ProNas软件：融合EFEA-SEA的新一代前沿技术产品。以ProNas能量有限元理论为基础，安世亚太联合国际最先进的中高频专家资源共同开发了拥有国内自主软件著作权的自主可控中高频噪声仿真分析软件ProNasoProNas是目前国际上中高频噪声预测分析领域内最新技术的结晶，是能量有限元分析(EFEA)和统计能量分析(SEA)领域的代表性解决方案。混合EFEA-SEA技术和基于能量有限元法的工程开发与应用，能够破解传统能量有限元技术很难广泛和深入应用于实际工程项目的困局。下图所示为ProNas软件界面。图6ProNas软件界面基于能量有限元理论进行中高频计算，ProNas天然具有建模灵活、计算效率高的优势，其理论本身还具有宽泛的阻尼和耦合强度适用范围，和随之带来的简单易学的用户操作界面，非常适用于结构声学问题的可行性研究、灵敏度分析及优化设计等。ProNaS的应用领域包括：汽车、船舶、轨道车辆、航空航天、汽车零部件、起重机械、家用电器、声学材料、通用机械、环境保护、建筑声学设计等。参考文献：UJ陈小剑.船舶噪声控制技术IM1.上海：上海交通大学出版社，2012:12-107.12林志立.能量有限元法的模型及其在车身中的应用D.中国科学院大学博士学位论文，2018,宁波.引原凯，王建民，韩丽，等.能量有限元在振动与噪声预示中的研究进展U.强度与环境2015.42(3)：12-19.14|祝丹晖，解妙霞，孔祥杰，等.复杂机械结构中高频动响应能量有限元方法研究J.中国工程科学，2013,15(1):107-112.15庞剑，湛刚，何华.汽车噪声与振动-理论与应用M.北京：北京理工大学出版社，2005:54-62.16陈书明，王登峰，马一功，等.模态密度计算精度对车内噪声预测精度的影响J.吉林大学学报，2010,40(2):335-340.17孙丽萍，聂武.能量有限元法在船舶结构中的应用J.哈尔滨工业大学学报，2008,40(9):1491-1494.网BERNHARDRJ.ThefamilyofEFAequationsandtheirrelationshiptoSEACProceedingsInternationalConferenceNoiseandVibrationEnergyMethods.1.yon：s.n.2000：248-257.19翁长俭.船体振动学M.大连：大连海运学院出版社，1992101邱斌.高速船全频段舱室噪声预报与控制方法的研究D.武汉理工大学硕士学位论文，2010,武汉.