用双互易边界元法进行声场特征值分析的研究

摘要

三维封闭空间(如船舶舱室、汽车驾驶室、工作间等)的噪声控制具有重要的实用价值.通过对声学系统的模态分析,可以得到其固有频率以及在此频率下声学量在声场空间上的分布,了解到声源及声场空间几何参数的变化对声学系统的相应的影响,为三维封闭空间的噪声控制提供了依据.双互易边界元法用于求解声场模态时形成的求解方程的系数矩阵与频率无关,并且可以很容易地形成广义特征值问题的求解方程,极大地方便了声场模态计算.本文对双互易边界元法计算声场模态问题进行了探讨.以Poisson方程为例对双互易边界元法的边界积分方程的推导以及数值实现进行了探讨.用一个二维算例验证了双互易边界元法的用于数值计算的可行性以及本文的编程思路,并初步讨论了近似函数和域内点对计算结果的影响.计算结果表明增加径向基函数的项数并不是改善结果的好办法,应该将更多的函数引入双互易边界元法中.在双互易边界元法中需要对由近似函数形成的矩阵求逆,这在计算大型的实际工程问题是很费时的.为此在声场模态的计算中引入不求矩阵逆的方法.分别用求矩阵逆与不求矩阵逆的方法,计算了两种边界条件下一个矩形二维声场的模态,并将结果与解析解和商业有限元分析软件ANSYS的结果作了比较.计算结果表明用求矩阵逆的双互易边界元法计算声场模态可以得到较好的结果,而用不求矩阵逆的双互易边界元法求解声场的特征值是可行的,但是不能将由此得到的特征向量直接用于描述特征频率下的声场分布.近似函数是影响双互易边界元法计算精度的主要因素,且各种函数在求解不同问题时的表现也是不同的.因此专门针对实际工程中的三维声场特征值问题对各种函数进行研究是非常必要的.本文选用多种近似函数对一个三维长方体声场模型在两种边界条件下的特征值做了计算.计算结果表明在相同的单元划分和相同的域内点的情况下,各种函数的结果是有差异的;对域内点的数量的敏感程度也是不同的.并以计算结果的相对误差为标准,讨论了所选各种函数的优劣,为在应用双互易边界元法计算三维声场特征值时如何选择近似函数指明了方向.

目录

文摘

英文文摘

独创性说明

1绪论

1.1本文的研究背景及意义

1.2国内外研究概况

1.3本文的主要内容

2计算结构声辐射的边界元法

2.1 引言

2.2 Helmholtz边界积分方程

2.3边界积分方程的离散

2.4奇异积分的处理

2.5特征频率问题的处理

3双互易边界元法

3.1引言

3.2双互易边界元法原理

3.2.1域内积分的产生

3.2.2域内积分转化为边界积分

3.2.3数值实现

3.2.4定义域内点

3.2.5近似函数

3.2.6需要注意的问题

3.3数值算例与分析

3.3.1数值算例

3.3.2结果分析

3.4小结

4基于边界元方法的声场特征值分析

4.1引言

4.2基于边界元理论的声场模态分析方法

4.2.1直接确定的方法(Direct Search Method or Determinant Search Method)

4.2.2计算域内积分的方法(Internal Cell Method)

4.2.3多互易边界元法(Multipli Reciprocity Mehtod)

4.2.4特殊积分方法(Particular Integral Method)

4.2.5级数展开法(Series Expansion Mothod)

4.3小结

5利用双互易边界元法进行声场模态分析

5.1引言

5.2基本原理

5.3不求矩阵逆的方法

5.4数值算例及结果分析

5.5小结

6采用双互易边界元法进行三维声场特征值分析的近似函数研究

6.1引言

6.2近似函数

6.2.1径向基函数(Radial Basis Function)

6.2.2TPS函数(Thin Plate Spline Function)

6.2.3 MQ函数(Multiquadric Function)

6.2.4全局函数(Global Function)

6.2.5混合函数(Hybrid Functions)

6.2.6用于无限域问题的函数

6.2.7计算三维声场特征值的近似函数的讨论

6.3数值算例及结果分析

6.4小结

7总结与展望

7.1本文工作总结

7.2本文工作中的一些问题

7.3关于相关问题的展望

参考文献

攻读硕士学位期间发表学术论文情况

致 谢

大连理工大学学位论文版权使用授权书