基于快速多极子边界元方法的大规模声学计算方法与应用研究

来源 :上海交通大学 | 被引量 : 11次 | 上传用户:qiaofei888
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
边界元方法是一种基于积分方程的数值方法,它可以将模型进行降维(三维变为二维,二维变为一维),且仅需要对边界进行离散。因积分方程能自动满足远场辐射条件,比较容易处理无限域问题的边界条件。作为一种半解析的数值方法,边界元方法在外部问题求解中比其它数值方法的精度更高,比较适用于声学分析计算。边界元方法虽然降低了模型的维数,但是生成的系数矩阵是满的、非对称的,甚至是奇异的,导致其系数矩阵存储量和求解计算量过大,限制了其在大规模声学计算中的应用。但目前声学分析和研究的问题逐渐向大规模、多物理场方向发展,如潜艇辐射和散射声场、大飞机起落噪声、生物组织声学特性及结构-声学优化等的模拟计算。这些模型的离散自由度将会在十万甚至百万的量级上,迫切地希望边界元方法能够快速地完成大规模结构声场的分析计算。快速多极子方法极大地降低了边界元方法的内存,提高了计算效率,为大规模声学问题的计算和分析提供了一种可能的选择。针对大规模声学计算对边界元方法的迫切需求,本文基于快速多极子算法,对自由空间、多联通域和半空间的声学问题,及声辐射模态理论进行了系统的分析,提出新的理论,发展相应的快速多极子边界元算法,以解决大规模声学计算中边界元方法内存占用量大和计算效率低的难题。并分别将这些新算法用于实验辐射声场的预测和大规模结构的声功率优化,充分显示所发展的快速多极子边界元方法在重大工程声学问题分析中的超强计算能力,及在以整体声学特性为目标的产品声学设计中的广阔应用前景。本文的主要研究内容如下:提出了一种解析的方法用来计算常数单元上奇异积分显式表达式中的两种线积分,并进行了误差分析。解决了常数单元离散下Burton-Miller方程中奇异积分的计算难题。对于高频快速多极子边界元方法,推导了常数单元和线性单元上的解析源点矩计算公式,提高了源点矩的计算效率和精度,从而进一步提高了高频快速多极子边界元方法的整体计算效率。将具有解析表达式的高频源点矩用于低频快速多极子边界元方法的上行传递计算中,在保留了低频快速多极子边界元方法内存使用量小的优点的同时,提高了源点矩的计算效率和精度。将所发展的低频和高频快速多极子边界元方法进行恰当的组合,改进了全频段快速多极子边界元方法。所发展的快速多极子边界元方法可以快速精确地对较宽频带内的大规模声学问题进行分析计算。将基于Burton-Miller方程的快速多极子边界元方法推广到多联通域声学分析。采用每个区域单独生成一个树状结构的方法对模型进行分组离散,解决了多联通域快速多极子边界元方法在算法实现上的困难。发展了一种基于边界块对角的预处理方法,改善了系数矩阵的条件数,减少了迭代求解的次数。所发展的快速多极子边界元方法大大提高了多孔吸声材料的设计、水下结构声学探测及生物组织中声学特性模拟等大规模多联通域声学问题的计算效率。发展了基于Burton-Miller方程的无限大阻抗平面上的快速多极子边界元方法。提出新的算法,将下行传递分成实栅格对实栅格和虚栅格两部分,推导了分段解析表达式用于M2L转换系数的精确、快速计算,从而巧妙地解决了Green函数中无穷积分项的快速多极子边界元方法的实现难题。将分段解析方法与多层树状结构相结合,进一步提高了M2L转换系数的计算效率,同时方便了临近栅格内单元间系数的直接计算。所发展的半空间快速多极子边界元方法拓宽了传统算法的应用范围,大大提高了半空间阻抗平面上快速多极子边界元方法的计算效率,从而可以对半空间阻抗平面上的大规模声学问题进行快速和精确的模拟计算。将快速多极子边界元方法和迭代特征值求解器相结合,发展了一种新的方法用于无限大障板上大规模结构的声辐射模态的快速计算。对于三维模型,基于等效源方法,Green函数的多极子展开式和边界积分方程,首次提出了映射声辐射模态的概念。理论上证明了Helmholtz方程的球函数基本解可以做为一组映射声辐射模态,并发展了一种较为简单的基于映射声辐射模态的辐射声功率计算方法。所发展的基于快速多极子边界元方法的映射声辐射模态,有效地解决了大规模辐射模型的声功率计算难题,并为基于辐射声功率的结构-声学优化提供了技术手段。基于快速多极子边界元方法对圆柱壳辐射声场进行了预测。提出新的基于Fourier级数和Lagrange函数的插值方法,更精确和方便地得到圆柱壳表面其它任意位置的表面加速度。然后使用所发展的快速多极子边界元方法,将插值得到的表面加速度作为输入,对声场进行分析计算。实测声场和预测声场的比较结果表明了所发展的快速多极子边界元方法及表面插值算法在圆柱壳结构辐射声场预测中的正确性。将所发展的快速多极子边界元方法,及基于此方法提出的映射声辐射模态理论,与优化算法相结合,以辐射声功率为目标函数,分别对压缩机壳体加筋位置和潜艇阻尼布局进行了优化。为进一步提高目标函数的计算速度,满足对潜艇模型的阻尼布局优化,将有限元商用软件和所开发的快速多极子程序在Matlab环境中集成,开发了并行算法。优化过程在可接受的时间内给出了较为理想的结果,充分显示了所发展的理论和算法在大型工程结构低噪声优化中的重要应用价值。开发了一款快速多极子边界元方法软件,为大规模声学分析提供一种强有力的计算工具。
其他文献
文章利用非条件分位数回归的分解方法对2011年城市外来劳动力群体的性别工资差距进行了分析。研究发现:(1)城市外来劳动力群体中,性别间受教育程度的差异已经变得不太显著,但
随着经济全球化脚步不断加快,对于我国港口航道工程进行信息化建设是非常重要的,不仅对时代进步有重大意义,更加强了我国港口航道事业的飞速发展。文章先分析了如今港口航道
随着世界全球化进程的加速,旅游业已然成为一个国际化产业。旅游文化资料的翻译目的是向外传播中国文化,吸引外国友人前来参观。本文以功能目的翻译理论为理论依据,以旅游胜
文章以2008—2010年民营上市公司为研究样本,实证分析了政治关联与过度投资的关系。研究发现:在其他同等条件下,相比没有政治关联的企业,有政治关联的企业投资支出更多,政治
随着近年国家对幼儿健康、幼儿教育的日趋重视,幼儿体育作为幼儿教育发展的重要组成部分,也成为我国教育领域密切关注的热点。但由于运动种类匮乏、运动场地稀缺等多因素制约,幼儿园体育活动的发展现状仍无法满足广大幼儿的运动需求,而其作为幼儿体育发展最为初始的基础环节,必须应加以重视。基于此,本文以炫彩球运动的研发及其可行性推广为研究对象,运用文献资料法、访谈法、问卷调查法、实验研究法、数理统计法等研究方法,
环境治理和能源开发成为了当代社会的两大主题,为了实现人与环境和谐发展,人们进行了多方位的探索和研究,随着不断地研究发现,半导体光催化技术对于有机污染物的分解和氢气的
针对蚁群算法运行参数选取问题,提出一种利用粒子群优化算法对蚁群算法的运行参数进行优化选择的方法。将蚁群算法的运行参数作为粒子群的位置信息,在算法迭代过程中使用粒子
无线传感器网络和无线射频识别技术通过IEEE802.15.4标准和ZigBee协议结合起来组成WSID网络,相对于现有的无线传感器网络与主动RFID,在通信距离、时钟同步、定位跟踪以及数据
本文以再读《生死场》为依据,放置三十年代"左翼"对它的评价,着重关注其女性身体体验、生存状况与命运的方面,以及战争在女性叙事中的矛盾性展示,重新思考其在当下的价值。
目的探讨不同运动疗法对膝骨性关节炎的治疗效果。方法 100例膝骨性关节炎患者随机分为治疗组(n=50例)和对照组(n=50例),两组均采用超短波电疗、关节活动度训练、有氧运动、A