波包络法在水下辐射噪声的应用研究.pdf

波包络法在水下辐射噪声的应用研究 尹绪超 白兆宏程晓达 （哈尔滨工程大学船舶工程学院 黑龙江 哈尔滨150001 ） 摘要：本文针对无限大流场声学计算问题，介绍了波包络法用于外域声辐射问题的计算理论。基于ABAQUS 有限元软件，针对浸入无限流场的脉动球体模型，以解析解作为参照，对比分析了波包络法及无反射边界 条件方法计算精度，验证了波包络法在水下辐射噪声应用的可行性。 关键词： 无限元法 波包络法 计算声学 辐射噪声 1 引言 在计算声学领域，求解无界域声场计算效率问题在过去半世纪以来一直是一个艰巨的挑 战，期间学者们提出了很多解决方法，但其中边界元法和无限元法格外引人关注：边界元法 基于边界积分方程，将无界域问题转化弹性边界问题，这种转化优势是自动满足Sommerfeld 声辐射条件；此方法积分方程只包含一个未知结构与流场交界面上的单层势（或双层势）函 数，而非整个结构流体模型，因此降低了问题求解的维数，这是边界元方法的一个巨大优势。 然而边界元法存在振荡和奇异积分问题，需要额外程序处理，而且形成的线性方程矩阵为满 带宽阵，因此难以求解大规模的工程实际问题；无限元法基于偏微分方程采用有限人工边界 截断声场，人工边界和结构之间人采用有限单元填充，人工边界外引入一种几何上无限大的 “有限”单元模拟无限大流场，克服了近似人工边界条件带来的截断误差，同时又具备了限 元方法的优势，与有限元方法相比，无限元方法在计算精度及效率上都有明显提高。 无限元这一术语1977 年最初由Bettess 提出，他基本思想是保留有限元方法的型函数， 然后乘以一个衰减函数，以体现声波衰减特性和实现物理意义上合理声学反射边界；后来， 1983 年，Astley 根据Bettess 无限元思想引入了一个新的有限单元类型，在这一类单元中， 作者以形状函数的共轭作为 Garlerkin 加权残值法的权函数。这一改进带来了极大的简化， 消除了单元矩阵中的不确定积分，作者应用此方法用于管道声传播问题计算中取得了很好的 效果并将此新单元类型称之为波包络法。但是波包络法形状函数的构造理论基于几何上和物 理意义上都不存在的虚拟声源概念，没有严格的理论推导，其精度需要进一步得到验证。本 文基于波包络法针对具有解析解的二维脉动壳体模型，分析了波包络法应用于声辐射问题的 可行性。 2 控制方程 G 考虑一浸入在无限大流场R 的弹性体结构声辐射或散射问题， 表其内边界，利用半径 G 为R 的球形边界 R 截断无限大流场，若不考虑问题的时域特性，则在流场R 内声压P 需要 满足Helmholtz 方程： 2 2 在域R 内 （1） ? p k P = 0 k 2 式中 代表波数， 为拉普拉斯算子。针对散射或声辐射问题，问题的边界条件各不 Ñ 相同，可表示为：  P 散射 P 总声压 - P 入射波 声散射   辐射 总声压  P  P P 声辐射  （2 ）  散射+辐射 总声压 入射波  P P - P 两者兼有   Γ e R Ω 此时问题解仍不具有