消声器声学性能预测的子结构快速多极子边界元法.pdfVIP

第27卷 第5期 计 算 物 理 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Vo1．27．NO．5 2010年 9月 Sep．，2010 文章编号 ：1001—246X(2010)05．0711—06 消声器声学性能预测的子结构快速多极子边界元法 崔晓兵 ， 季振林 (哈尔滨23程大学 ，动力与能源工程学院，黑龙江 哈尔滨 150001) 摘 要 ：鉴于快速多极子边界元法 的应用主要局限于单区域声学 问题计算 ，发展基于子结构技术 的快速多极子 边界元法以计算 多区域声场问题 ，介绍基本原理 、具体实施过程 以及优缺 点．以带有插进 口管 的膨胀腔消声器 为例 ，应用子结构快速多极子边界元法和传统边界元法计算其传递损失，通过与实验测量结果的比较 ，验证方 法 的有效性和计算精度．研究表 明，快速多极子边界元法与传统边界元法相 比，节点数越多，其在节省计算时 间，减少计算量等方面的优势越 明显． 关键词 ：子结构快速多极子边界元法 ；多区域声场 ；消声器 ；传递损失 中图分类号 ：TB535 文献标识码 ：A 0 引言 边界元法 以其能降低问题维数 ，计算精度高等优点得 以广泛应用．然而 ，传统边界元法 (CBEM)形成的 是非对称 的满秩矩阵，若采用间接法求解 ，其计算量和计算时间均为 O(N)，限制了其在处理中高频声学和 大区域声场问题中的应用 ，成为边界元法应用的一个瓶颈．为了解决大尺度问题的计算 ，快速多极子法 “ (FastMultipoleMethod—FMM)近年来得 以发展 ，它是 由Rokhlin提出的一种快速算法 ，目前已被广泛地应用到 电磁场 、流体力学和应用物理等诸多领域 ，被美国应用数学物理学会评为20世纪计算程序前 10名． 近年来 ，为了提高边界元法在节点数 目较多时的适用性 ，将 FMM 引入边界元法，形成了快速多极子边 界元法 (FMBEM)．通过矩阵预处理技术 ，FMBEM 的总计算 时间和全部计算量与 NlnN成正 比，极大地减少 了计算机的处理时间和计算量 ，成为计算大尺度 问题 的有力工具 J．消声器声学性能计算与分析是 FMBEM 的一个典型应用领域 ，而将 FMBEM应用于实 际消声器的声学性能计算还有一些具体问题需要处 理：①FMBEM 目前主要应用于单区域声场问题的计算，而消声器内部通常含有薄壁结构，形成奇异边界，限 制了FMBEM的直接应用 ；②对于细长结构的管道和消声器系统，直接应用 FMBEM会引起多极子展开式过 大的累积误差 ，降低计算精度 ；③对于复杂结构的消声器 ，若直接对整体结构表面划分网格 ，会使得节点数过 于庞大 ，即使进行系数矩阵的预处理 ，其迭代收敛速度仍不理想 ，使其高效性受到极大影响． 本文研究的 目的就是寻找一种能够有效解决上述 问题 的实际方法．鉴于消声器的结构特点 ，研究基于子 结构技术的快速多极子边界元法 以计算复杂结构 内部 的声场和消声器 的声学性能．以带有插管的膨胀腔消 声器为例计算其传递损失 ，进而讨论 FMBEM和 CBEM 的计算精度和效率． 1 快速多极子边界元法 快速多极子算法的原理为 ：首先把多个单元归结到一个区域 ，其次把 区域 内单元的影响都集聚到区域的 中心，然后计算各区域之间的相互影响关系 ，最后计算得到的中心点的值再返 回各单元．应用快速多极子算 法于边界元法中，加速其边界节点之间影响系数的计算速度 ，减少计算量 ，便形成了FMBEM． 1．1 Helmholtz方 程基 本解 的多极 子展 开式 基本解的多极子展开式是将 FMM应用于边界元法的前提 ，根据 Gegenbauer附加定理 [93，Helmhohz方程 收稿 日期 ：2009—09—09；修 回日期 ：2009—12—23 基金项 目：国家 自然科学基金 “大尺度 声学问题 计算 的快速多极子边界元法及其应用研 究”资助项 目 作者简介 ：崔晓兵 (1983一)，男 ，河北 景县 ，博士生 ，从事振动与噪声控制研究 ，哈尔滨工程 大学动力 与能源工程学 院