第二章船舶辐射噪声与计算方法.pptVIP

如果从能量传递的角度来研究就可以不必分清能量是哪阶模态传递的， 即可以从统计的角度分析密集模态平均的振动能量传递水平， 也就避免了CMA的局限性， 这就是统计能量分析( Statistical Energy Analysis，SEA)的基本思想。 潜艇及船舶结构辐射噪声计算方法 1．有限元法及边界元法原理 将有限元法应用到声振耦合问题时则需要在结构的外部紧贴结构表面“构造”一个流体结构，对结构和流体都要进行FEM网格离散，建立FEM方程进行分析。 有限元法及边界元法 1．有限元法及边界元法原理 假定流体是理想的声学介质，且满足如下条件： (1)假定流体是可压的，但只允许压力与平均压力相比有较小的变化，流体是各向同性、均匀的。 (2)声波动过程是绝热的。 (3)假定流体为非流动并且无黏的（即黏性不引起耗散作用）。 (4)假定流体平均密度和平均压力不变，计算中求解的压力是偏离平均压力的相对压力而不是绝对压力。 有限元法及边界元法 用FEM对模型进行分析是通过对模型进行离散，同时求解波动方程以及结构运动方程来实现的。在流体介质中声波应满足波动方程 分别对结构和流体进行单元离散最后可得到矩阵形式的流体波动方程： 有限元法及边界元法 式中： 为流体的质量矩阵； 为流体的刚度矩阵； 为流体面上的耦合矩阵。 有限元法及边界元法 在流体与结构的交界面上，声压对结构同时也产生一个面力的作用，将其变换到节点上，于是结构方程可以写成下列形式： 用统一的矩阵形式表示有 有限元法及边界元法 在无限远处要满足Sommerfeld辐射条件： 利用波动方程，可以转化为单频声场的Helmholtz积分方程： 有限元法及边界元法 在流固耦合面S上，存在边界条件： 式中：E，S，I分别表示位于结构外部、表面上和内部； Q为结构表面上的点，p为空间中的点； G(P，Q)为自由空间的格林函数， ； 为边界表面的法向速度； n为结构表面外法向单位矢量； 为波数。 有限元法及边界元法 利用BEM法进行离散，即可得到BEM求解方程： 式中：B和C为系数矩阵；p为节点声压矢量； 为节点法向速度矢量。 得到表面压力后，结构面上i节点处单频声场的声强In由下式给出： 作面积积分即可得到结构辐射声功率： 有限元法及边界元法 ANSYS与VirtualLab声学分析 ANSYS软件处理多种物理场耦合的功能很强，并且单元类型丰富，可满足一般工程问题的需要。其缺点是声场的后置处理功能较弱。 而Virtual Lab软件是有关声学和振动的专用软件，它既有有限元技术，又含有边界元技术，可计算一般复杂弹性结构的水下耦合振动问题。其优点是对声场的后置处理功能很强，但它没有网格划分的功能。 根据ANSYS和Virtual Lab各自的特点，将两者结合起来使用。即首先在ANSYS中对结构进行有限元建模，计算模型在激励力作用下的动力响应，再利用结构表面的位移作为边界元的边界条件，在Virtual Lab软件中计算各种声学物理量。 有限元法及边界元法数值实现方法 1)单元类型的确定 在ANSYS建模中采用的三维声学单元类型有Fluid30和Fluid130。 Fluid130是一个表面单元，用于模拟有限截断流体域的边界，每个节点只有一个声压自由度。 当声波传到边界上时能被Fluid130单元“吸收”，只有少量反射回流体域，这样就能较精确地模拟出无限水深的实际情况。 有限元法及边界元法 1)单元类型的确定 ANSYS中已确定吸收单元远离结构或振动源区域以外大约0.2λ放置能产生准确的结果， 其中：λ=c/f为压力波的主波长， c为流体中的声速， f为压力波的主颇。对于直径为D的圆柱壳体，封闭边界的半径则至少应为 D/2+0.2λ。 有限元法及边界元法 2)网格精度的确定 一般来说，网格越是密集计算结果越精确，随着网格密度的增加，计算结果会收敛于精确解，但同时也会带来计算量的急剧增加。 一般来说，为了描述振动波在结构中的传播，一个波长内至少应该有5个节点，即4个单元。 有限元法及边界元法 3.计算结果分析 试验中将模型置于水下，在激励力的作用下，模型最终会产生一个稳态振动，并向周围辐射噪声。 周围的声场反过来由会影响结构的振动，它们之间的交互作用最终会形成一个稳定的声振耦合系统。 有限元法及边界元法 试验模型 有限元法及边界元法 该模型在轴向中心处沿周向角度分别为0°、45°、90°、135°、180°的距离耐压壳