非共形面声全息法研究水下弹性结构的振动和声辐射.pdfVIP

非共形庙礁景鎏研究壅至理壁绩垫的振动与妻擐孰 商德江 咕尔滨工挣大学 1．引言 水下结构振动与声辐射的研究对于潜艇及舰船之降噪有重要意义。水下结构与声藕合振动的 分析方法很多。常用解法有：有限元嘲洒数匹配展开法和T一矩阵浅以及有防乱边界元方法。其 中，有限元十边界元杰讫磕嘲中广泛使用的方法。其特点是：它以离数化的形式给出了结构表面振 动艟眵和声压的关系式利用数值计算技术进行求解。因此对结构黼的n侗形状和内部结构不规 则吐有较强的适应性。但是当计算漠型(如加肋双层壳潜艇)内部结构比较复蕊并且计算频率较 高时，有限元+边界元计算技术的前期建漠工作量很大计算时间较长。而利用非共形面声全息法可 较盘于j瑚端b此类问题尤其当我们只知雌自构外表面|《动分布及其辐射声场时，非共形面声全息法尤 为实用。 2．基本原理 所谓的非共形面声全息是指声测量面与重建场面是非共形的。本文采用的是利用边 界元技术的数字式的非共形面声全息。这种方法对水下结构表面形状有很强的适应性。其 基本原理是：首先，将结构的外表面用边界单元离散，利用表面Helmholtz积分公式，求出结 构外表面法向位移与表面声压之间的数值关系式： IP81=IABlCUnl L1、 其中：[Un]=结构表面法向位移：【Ps】=结构袭面复声压： [AB】=系数矩阵 之间的数值关系式： 【Pfl=(c1(Psl+lD3 其中：(Pf】=辐射场声压：It]，[D】=系数矩阵 将方程(1)代入方程(2)，可得： [Pf]=fFJfoil] (3) 其中：[F]=[C][AB]+[D] 最后，由方程(3)，利用奇异值分解和最小二乘法可求出以辐射场声压为已知量的结构表面 法向位移的矩阵关系式： 鼬n1=碑invtF)1 其中：[Piny(F)】=[F】的Moore～Penrose广义逆 方程(1)和方程(4)即为声全息基本关系式。因此，只要测量出结构模型辐射声场的近场复 声压就可计算出结构表面的振动位移和声压的分布，近而可计算出结构辐射的声功率和远 场指向性等各声学量。为便于工程上扫描测摄，可取全息测量面为位于结构两侧的双平行 平面。 ·117· 3．数值计算 为验证利用非共形面声全息法研究水下结构振动与声辐射特性的可行性，本文采用点 源替代的办法进行模拟。首先，假设结构模型外表面为一椭球面，将一点声源置于其中．由 点声源公式计算出椭球面上质点的法向位移，及椭球面两侧的双平面上的复声压。然后，利 用双平面上的复声压由声全息关系式(4)和(1)可计算出椭球面上的法向位移和声压，可将 其与由点声源公式直接计算的解析值相比较。即可验证非共形面声全息法的正确性。本文 模拟计算参数如下：椭球面长短轴之比为2(米)／I(米)，椭球表面分割为8个6节点三角形 单元和112个8节点四边形单元，共有480个节点；椭球两删全息面(8蠼米)到椭球中心的 距离为1．46米。相邻两点间距为0。53米(约为i／5波长)，共有512个测点：计算频率为 562Hz。计算结果误差见表(一)和表(二)。 衰(一】：声堂息幅值谩麓 量大使氆l 量小值∞ 平坶啦∞ 最面声压 9．62 Ot8 2．37 表面也穆 8．68 0．-ll 220 衰(：)：声全息相位误差