WB+法在多域声学分析与应用.pdfVIP

第 4 卷第 3 期 动 力 学 与 控 制 学 报 Vol. 4 No. 3 2006 年 9 月 JOURNAL OF DYNAMICS AND CONTROL Sep. 2006 WB 法在多域声学分析中的应用 彭伟才 　何 　锃 (华中科技大学力学系 ,武汉 　430074) 摘要 　描述了振动声系统建模技术的基本概念. 根据域分解的连续性条件 ,讨论了界面的压力和速度连续 以及阻抗连续 ,应用加权余量法推导了两者的耦合模型. 并用 LMSSYSNOISERev 5. 5 进行了有限元数值模 拟 ,计算结果与有限元结果符合得较好. 通过比较两种连续性条件 ,发现前者更适合较小的计算模型而后者 更适合较大的计算模型. 最后对域分解提出了几个简单优化原则. 关键词 　声学 ,多域 ,域分解 ,Trefftz 法 用统计能量法 ,而中频段还没有比较适用的方法. 引言 WB 法收敛的一个充分非必要条件是所求声 ( ) ( ) 当今分析声和结构声 耦合 的数值方法中基 域为凸形 convex . 对于非凸形的声域要分解成许 ( ) 于单元技术的主要是有限元和边界元. 有限元法 多的子域 凸形 . 子域之间通过连续性条件进行耦 (FEM) [1 ] 将声域离散成大量单元 ,每一个单元上的 合. 本文讨论了这种方法在二维声问题中的应用. ( ) 动力学变量采用简单 、近似的形函数 多项式 来描 [2 ] 1 　基本理论 ( ) 叙. 边界元法 BEM 将声域的边界离散成大量单 元 ,每一个边界单元上的动力学变量采用简单 、近 1. 1 　定义 [7 ] 似的形函数来描叙 ,场变量由边界变量的近似值通 WB 方法的详细推导可参考文献 . 考虑一个 过边界积分方程得到. 2D 的凸形有界声域 Ω,如图 1 所示 ,其密度为 ρ和 由于每一个单元上的变量采用近似的形函数 声速为 c . 声域 Ω的载荷条件 :1) 在 rq 点受到源强 描叙 ,为了达到要求的精度 ,需要大量的单元. 随着 度幅值为 q 的点源 ;2) 给定的边界条件. 域 Ω内一 频率的增加 ,发散误差也增加 ,于是网格密度和迭 点 r ( x ,y ) 的压力 p 的非齐次 Helmholtz 方程 ( ) 2 代次数 为了达到要求的精度 也随着频率增加. 因 ( Δ ) ( ) ωδ( ) Ω ( )