基于复数矢径的波叠加法解声辐射问题.docVIP

Ξ基于复数矢径的波叠加法解声辐射问题宇1 ,2黄玉盈2向( 1 广西工学院汽车工程系 ,柳州 ,545006 Ξ 基于复数矢径的波叠加法解声辐射问题 宇1 ,2 黄玉盈2 向 ( 1 广西工学院汽车工程系 ,柳州 ,545006) ( 2 华中科学大学土木工程与力学学院 , 武汉 , 430074) 摘 要 利用波叠加法与结构动力分析中的相似性 ,提出了一种在波叠加法中克服解非唯一的通用方法 ,即 在虚拟源强系统中加入一定的虚拟阻尼从而能获得全波数域内的唯一解 ,并以此为基础提出了一种新的加入虚拟 阻尼的方法 ———复数矢径波叠加法. 文中给出了脉动球和摆动球两个数值算例 ,计算结果表明 :本文方法不仅能有 效解决数值求解过程中解非唯一的问题 ,且计算时间只与标准波叠加法相当 ,计算精度却比同类方法高. 关键词 虚拟边界积分方程 ,波叠加法 ,复数矢径波叠加法 ,声辐射 ,虚拟阻尼 0 引言 进行线性组合形成新的边界积分方程 ,解决了特征 频率处解的非唯一性问题. 但由于法向求导增加了 积分核的奇异性 ,数值积分处理十分繁杂 . 1989 年 , Koopmann 等 人 根 据 Burgess 和 自上 世 纪 60 年 代 以 来 , 利 用 边 界 积 分 方 程 (B IE) 分析任意形状空穴的振动声辐射和散射问题 一直是国内外学者致力研究的一个课题. 边界积分 方 程 由 于 具 有 降 低 求 解 空 间 维 数 、自 动 满 足 Helmholtz 方程以及 Sommerfeld 辐 射 条 件 等 优 点 , 因 而在分析无穷域的声辐射和散射问题上 ,该方法比 有 限 元 素 法 ( FEM) 更 为 有 效 . 但 是 , 利 用 常 规 的 Helmholtz 边界积分方程法求解振动空穴的外声场问 题 ,存在两个数值计算上的困难 : (1) 奇异积分处理 ; (2) 在 Dirichlet 内问题的特征频率处 ( 本征值) ,解的 唯一性不能保证1 ,2 ,尤其是在中 、高频段 ,特征频率 分布密集 , 导致其 求 解 结 果 严 重 失 真 , 甚 至 完 全 错 误 . 为解决这一问题 ,近几十年来人们提出了多种改 进方 法 , 其 中 最 具 代 表 性 的 有 Schenck3 提 出 的 CHIEF 法 和 由 Burton 、Miller 提 出 的 Burton2Miller 5 Mahajerin 求解弹性力 学 问 题 的 叠 加 法 思 想 , 提 出 了求解声辐射和散射问题的一 种 间 接 积 分 方 程 法 6 ———波叠加法 ,该方法是将一系列虚拟源强配置 在空穴内部的一虚拟区域内 ,然后根据振动空穴表 面给定的法向振速 ,采用配点法或最小二乘法计算 出虚拟源强的强度值 ,进而得出虚拟源强系统的模 拟外声场. 在该方法中由于所选择的虚拟源强是满 足 Helmholtz 方程和 Sommerfeld 辐射条件的波函数 , 由解的唯一性定理可知 ,该虚拟源强系统所产生的 模拟外 声 场 即 为 原 振 动 空 穴 的 辐 射 或 散 射 声 场 . Benthien 和 Schenck 等人指出7 : 由于在波叠加法中 所选择的虚拟源强配置区域与 实 际 振 动 空 穴 不 重 合 ,在求解中不存在奇异积分处理 ,因此与其他边界 积分方程方法相比较该方法具有较高的计算精度和 效率. 在实际应用中 ,为计算方便 ,一般将源强配置 区域 选 为 空 穴 内 部 的 一 条 封 闭 曲 面 ( 虚 拟 边 界) , J eans 、Mathews8 和 Wilton 等 人9 指 出 : 正 是 由 于 采 用一条封闭曲面作为虚拟源强的配置区域 ,导致其 在虚拟曲面上相应内问题的特征频率处存在非唯一 解 ,他们 提 出 采 用 复 数 形 式 的 Burton2Miller 型 组 合 层势法克服在相应内域问题的特征频率处解的非唯 一性. 但计算时间比采用标准单层势法或双层势法 增加 50 %左右7 . 法4 . CHIEF 法是以常规 Helmholtz 边界积分方程为 基础 ,通过在内域中选取若干个点 ( CHIEF 点) 而形 成内点积分方程作为补充方程 ,与边界积分方程联 立形成超定方程组 ,在最小二乘意义下求得问题的 唯一解. 但当所选择的内点与 Dirichlet 内问题的驻 波节点重合时 ,该方法失效. 尤其是在高频段 ,驻波 结点分布密集 ,欲选取内点使其远离驻波节点是很 困难的 ,因此从理论上讲 CHIEF 法并没有完全解决 计算中解唯一性的问题 . Burton2Miller 法是