收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

复杂结构声学特性预测的快速多极子边界元法研究

崔晓兵  
【摘要】:边界元法是求解复杂结构声学问题的一种有效数值方法。由于传统边界元法(CBEM)形成满秩的系数矩阵,其计算量、存储量和计算时间为O(N2),限制了其在大尺度声学问题中的应用。相比之下,快速多极子边界元法(FMBEM)具有较高的计算效率,但是无法直接应用于带有薄壁结构以及含有吸声材料的复杂结构内部声场问题的计算。鉴于此,本文发展了两种适用于大尺度声学问题计算的快速算法:(1)适用于阻抗复合结构内部声场计算的子结构快速多极子边界元法(Sub-FMBEM),(2)适用于具有多薄壁元件的复杂抗性结构内部声场计算的快速多极子双重边界元法(FMDBEM)。所做的具体工作如下: 详细阐述了多级快速多极子边界元法(MLFMBEM)计算内部三维声场问题的理论基础和数值实施过程。研究表明,MLFMBEM的计算量和内存量均能达到O(NlnN)。应用三级FMBEM对简单结构消声器的传递损失进行了数值计算,通过和CBEM比较验证了FMBEM的计算精度及在高网格数时的计算效率。 将子结构技术应用到FMBEM中,形成了适用于阻抗复合结构内部声场计算的Sub-FMBEM,详细介绍了该方法的基本原理和总体矩阵向量积的计算过程,对不同吸声材料下单级和多级传递关系计算进行了考察。研究表明,展开式截断项数计算参数的选取在一定程度上影响了格林函数展开式的计算,而吸声材料的选取并非是影响其计算精度的主要因素。复波数的引入使展开式计算产生偏差,当复波数虚部与展开点间距的积k i i rLM大于13.5时,展开式值逐渐偏离理论值。因此,当k i i rLM值处于失真范围时,可将展开式中转移因子的值以0来近似,或者应用Sub-FMBEM将大尺寸结构模型划分为若干小尺寸结构分别计算。本文采用ILUT预条件处理技术和Bi-CGSTAB(l)迭代求解器进行Sub-FMBEM的迭代计算。在实际的编程计算中,未知量列向量的构建次序及离散节点编号顺序对迭代法的收敛速度有着重要影响,结合多子结构模型的声场计算给出了Sub-FMBEM整体系数矩阵的构建方法和原则。应用Sub-FMBEM对抗性和阻抗复合型消声器的声学性能进行了计算,通过与Sub-BEM计算结果以及实验测量结果的比较验证了该方法的正确性和计算精度。在计算效率方面,Sub-FMBEM在高网格数时有明显优势,而在低网格数时则略低于Sub-BEM。在迭代收敛特性方面,Sub-FMBEM的预处理时间和迭代计算时间均随着频率的增大而增加,而Sub-BEM随频率基本不变。网格数越多,Sub-FMBEM相对于Sub-BEM能体现高效性的频率范围也越大。无论是预处理阶段还是迭代阶段,Sub-BEM的计算时间随网格数的增长率远远大于Sub-FMBEM,从而验证了FMBEM在大尺度中高频声学问题计算中的优势。此外,吸声材料的引入不但在一定程度上降低了Sub-FMBEM的迭代次数,而且使之随着频率的变化更为稳定。因此,Sub-FMBEM更适合于阻性结构声场的计算。 由双重边界积分方程出发,推导出了适用于复杂抗性结构内部声场计算的边界积分控制方程,并结合快速多极子算法加速求解,创建了适用于含多薄壁结构的复杂抗性消声器内部声场计算的FMDBEM,给出了具体的求解过程以及超奇异积分的处理方法。FMDBEM的主要优点在于:无需针对薄壁结构划分子结构,并对薄壁边界和虚拟交界面重复离散,只需对所有薄壁结构离散一次,从而减少了网格数量;其缺点在于:不能直接应用于含有吸声材料的阻性消声器声场计算,而且要求抗性消声器中薄壁结构的厚度必须足够小到可以忽略不计。此外,由于控制方程中未知量系数计算存在超奇异积分,需要在计算中进行适当地弱奇异化处理。通过应用FMDBEM、Sub-FMBEM和Sub-BEM对复杂抗性消声器传递损失进行计算,验证了FMDBEM的计算精度和计算效率。在同样的网格尺寸下,由于FMDBEM减少了网格数量,从而在一定程度上节省了计算时间。 为进一步检验本文所发展的数值方法的正确性,使用阻抗管测量系统结合两负载法测量了直通穿孔管消声器和阻抗复合型消声器的传递损失。比较了实验测量结果与FMBEM的计算结果,它们在较宽频域内吻合良好,从而进一步验证了本文数值方法的正确性和在复杂结构内部声场计算中的适用性。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 卢习林;带横孔圆轴三维应力分析的边界元法[J];清华大学学报(自然科学版);1984年01期
2 伊增欣;边界元法应用的进展(二)[J];海洋预报;1986年03期
3 C.A.Brebbia ,杨伟 ,曾华;工程实践中的边界元法[J];力学进展;1986年03期
4 丁立,袁修干;人体温度场的边界元分析[J];北京航空航天大学学报;2000年06期
5 杨再朝;用边界元法进行区域重力测量外部改正的探讨[J];石油物探;1987年01期
6 朱浮声;间接边界元法的物理表述[J];东北大学学报(自然科学版);1987年04期
7 陈健华;;用边界元法计算光弹性应力分析中的主应力和[J];河北理工学院学报;1987年01期
8 潘安福;边界元法的边界效应问题[J];武汉大学学报(工学版);1988年06期
9 刘忠民;邓次平;;一种提高边界元法精度的有效方法——吻边边界元法[J];车辆与动力技术;1988年01期
10 叶建乔;;边界元法分析周边简支扁壳的一种简便计算[J];工程力学;1989年02期
11 孙淑苓,刘玉琦,周安宁,冯振兴,叶碧泉,沈成武;油液晃动及箱体耦联振动计算方法研究[J];航空学报;1990年05期
12 徐汉忠;;应用功的互等定理推导不连续位移基本解[J];河海大学学报(自然科学版);1990年04期
13 ;《物探化探计算技术》1990年1~4期总要目[J];物探化探计算技术;1990年04期
14 符力耘,牟永光;弹性波边界元法正演模拟[J];地球物理学报;1994年04期
15 殷少有,关志强;用边界元法求解吸气罩流场[J];湛江海洋大学学报;1997年01期
16 黄庆学,肖宏,申光宪;用边界元法分析轧制压力和摩擦力[J];计算力学学报;1998年02期
17 王文全,何宪刚;用积分方程离散——边界元法解含热源的稳态三维热应力问题[J];机械设计与制造;1998年02期
18 王小贞,臧跃龙;粘性液体小幅晃动常单元及线性元的分析比较[J];重庆建筑大学学报;2000年06期
19 蒲军平,姚振汉;二维接触问题边界元法的研究进展[J];天津理工学院学报;2000年01期
20 谭军安;边界元法在环境声学中的应用[J];噪声与振动控制;2002年03期
中国重要会议论文全文数据库 前10条
1 程长征;牛忠荣;胡宗军;叶建乔;;边界元法计算浅表面裂纹应力强度因子[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
2 冯正和;郝建斌;;用网络-边界元法解旋转对称结构的电磁波电路[A];1987年全国微波会议论文集(上)[C];1987年
3 王海涛;姚振汉;;3D快速多极边界元法的一些应用[A];北京力学会第11届学术年会论文摘要集[C];2005年
4 李世辉;;边界元法普及于隧道工程的途径[A];边界元法在岩石力学和工程中应用会议文集[C];1987年
5 王泳嘉;邢纪波;;离散元法及其与边界元法的耦合[A];边界元法在岩石力学和工程中应用会议文集[C];1987年
6 夏元友;徐长佑;;三维弹塑性边界元法的奇异积分处理[A];第二届全国结构工程学术会议论文集(上)[C];1993年
7 戴瑛;郑百林;贺鹏飞;;单纤维拔出实验模型的力学分析及脱粘判据研究[A];复合材料的现状与发展——第十一届全国复合材料学术会议论文集[C];2000年
8 姚振汉;王海涛;雷霆;王朋波;;快速多极边界元法研究进展[A];中国力学学会学术大会'2005论文摘要集(上)[C];2005年
9 张诗德;;广义传递矩阵和边界元组合法[A];第二届全国结构工程学术会议论文集(上)[C];1993年
10 张宏;;工程中渗流问题的边界元算法与分析[A];第五届全国结构工程学术会议论文集(第一卷)[C];1996年
中国博士学位论文全文数据库 前10条
1 郑昌军;三维声学敏感度分析的宽频快速多极边界元法研究[D];中国科学技术大学;2011年
2 王英俊;结构分析中的GPU并行快速多极边界元法研究[D];华中科技大学;2013年
3 王英俊;结构分析中的GPU并行快速多极边界元法研究[D];华中科技大学;2013年
4 崔晓兵;复杂结构声学特性预测的快速多极子边界元法研究[D];哈尔滨工程大学;2012年
5 刘士利;改进的边界元法及其在电场计算中的应用[D];华北电力大学(北京);2011年
6 张健飞;机群环境下的并行边界元法研究及其在水工结构分析中的应用[D];河海大学;2004年
7 刘永健;三维弹性体移动接触问题的边界元法研究[D];清华大学;2003年
8 胡宗军;边界元法中高阶单元奇异积分的一个新正则化算法及其应用研究[D];合肥工业大学;2012年
9 王海涛;快速多极边界元法研究及其在复合材料模拟中的应用[D];清华大学;2005年
10 刘清珺;Galerkin边界元法及其在结构软化分析中的应用研究[D];清华大学;1995年
中国硕士学位论文全文数据库 前10条
1 张晨利;边界元法中几乎奇异积分的正则化算法及其在弹性薄壁结构中的应用[D];合肥工业大学;2001年
2 谢昌林;机械结构振动噪声的仿真[D];西北工业大学;2004年
3 王锋;用双互易边界元法进行声场特征值分析的研究[D];大连理工大学;2005年
4 邓军;基于边界元法的变电站内工频电场计算方法研究[D];重庆大学;2010年
5 赵世伟;边界元法在中厚板筏式基础分析中的应用[D];兰州理工大学;2003年
6 程长征;几乎奇异积分正则化算法在多域及接触边界元法中的应用[D];合肥工业大学;2003年
7 刘朝霞;应用特解边界元法对热疗温度场的数值计算研究[D];北京工业大学;2000年
8 刘姜玲;合成绝缘子三维电场分区边界元法的研究[D];华北电力大学(北京);2003年
9 孟玮;边界元法在结构断裂分析中的应用[D];广西大学;2005年
10 王海涛;快速多极边界元法在二维弹性力学中的应用[D];清华大学;2002年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978