快速多极边界元与有限元耦合分析结构声学敏感度的算法研究

来源 :中国科学技术大学 | 被引量 : 0次 | 上传用户:Z_PEPSI
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
浸没在流体中的弹性结构的振动声辐射与散射是一项非常重要的研究课题。比如潜艇或船舶的水下辐射噪声影响了其隐蔽性,限制了海军装备及民用水声设备的使用性能。与在空气中不同,水下结构的机械阻抗并没有远远大于水的声阻抗而只是其数倍或数十倍。因此进行像潜艇或船舶这样大型复杂结构的振动或辐射、散射声场分析,必须考虑流体声场与结构振动的相互耦合作用。在可压缩流体介质中结构振动产生辐射声压,同时声场又反作用于结构,引起结构的附加振动。单独求解结构动力学方程或流体运动方程无法得到正确结果,必须结合结构与流体动力学方程,通过求解这一声振耦合方程组,方可得到结构振动响应和流体声场的正确结果。基于解析求解方法的声振耦合分析为声振耦合问题的机理研究奠定了坚实的理论基础。然而解析方法难以应用到复杂实际结构的声振耦合分析,因此发展合适的数值方法计算这一问题是十分重要的。有限元法在结构振动响应分析方面得到了广泛的应用,而边界元法在无限域声场分析方面有独特的分析优势。因此通过结合有限元与边界元法形成的耦合有限元/边界元法(FEM/BEM)进行结构和流体耦合场分析是十分有效的。本文主要进行水下结构振动辐射声场或散射声场的快速算法研究,开发FEM/快速多极非连续边界元法算法进行水下声振耦合问题的求解,进而建立结构声学感度分析算法,为基于梯度的结构优化分析提供必要的理论和算法基础。论文的主要内容和创新点包括:(1)基于快速多极边界元法的二维无奇异声学及其敏感度分析的算法建立。本文使用Burton-Miller法克服解的非唯一性问题,针对该方法产生的各阶奇异积分,采用常量单元离散声学边界,同时使用Cauchy主值积分与Hadamard有限部分积分法直接计算奇异积分项,推导出一组二维无奇异声学及其敏感度边界积分方程。针对传统边界元法形成的是非对称满系数矩阵,求解与存储这个满阵需要花费大量的计算时间和内存的问题。本文使用宽频快速多极算法(FMM)加速系统方程的求解,在低频处使用低频算法,在高频处使用高频算法,保证了快速算法在整个频段的稳定性与有效性,明显提高了边界元法在声学及其敏感度分析上的计算效率。(2)基于快速多极非连续边界元法的三维声学及其敏感度分析的算法研究。对于非连续边界元,插值节点放置在单元内部,其连续性有保证,避免了角点问题的处理,是一种广受重视的实用单元。本文针对不同类型的非连续边界元,推导了三维无奇异声学及其敏感度表达式,并通过算例对比连续元与非连续元的计算精度,给出了最优单元类型。最后引入快速多极算法,形成快速多极非连续边界元算法,进行声学及其敏感度计算,大幅提高了计算效率和降低了内存占有量,使得边界元法在大规模实际问题的敏感度分析上的应用变得可能。(3)基于FEM/非连续FMBEM耦合算法的结构声学分析算法研究。推导出了用于计算耦合声场分布、适于引入快速算法的耦合边界元方程,并进而推导出了水下结构振动辐射声功率计算表达式。通过数值算例考察FEM/非连续BEM耦合算法的计算精度与节点位置参数的关系,得到一组用于计算声振耦合问题的优化节点位置参数值,并对比耦合连续单元与耦合非连续单元的计算效率,得到计算效率最高的耦合单元类型。本文通过结合FEM和非连续FMBEM求解声振耦合问题,为复杂结构声振耦合方程的求解提供了有效的数值分析工具。(4)基于FEM/非连续FMBEM耦合算法的结构声学单与多设计变量的敏感度分析算法研究。本文推导出基于直接微分法的结构振动辐射或散射声场敏感度表达式,设计变量可以选为流体与结构材料参数、结构形状尺寸参数,例如流体密度、结构密度、结构泊松比、杨氏模量、壳厚度、结构形状尺寸参数。针对不同的设计变量,本文推导出相应的结构声学感度表达式。对于单一设计变量的感度分析,使用直接微分法能有效得到高精度的计算结果;对于多设计变量的感度分析,推导得到基于伴随变量法的结构振动辐射声功率感度表达式。引入不含对设计变量的感度项的伴随方程,在进行多设计变量的感度分析时,该伴随方程只需求解一次,即可用于不同设计变量的感度计算,因此大幅提高了本文FEM/非连续FMBEM耦合算法的计算效率。因此本文的研究可以明显提高结构辐射和散射声场的计算效率,给出了计算声振耦合问题的最优耦合有限元与边界元单元类型,并为更大规模更接近实际的数值分析提供有效分析算法。并通过结构声学感度分析进行结构形状优化设计,得到具有最低辐射声功率的优化结构外形,以大幅提高水下结构的声隐声性能。
其他文献
在现代化、信息化高速发展的今天,办公空间在人们的工作和生活中,扮演者越来越重要的角色。人们办公环境也随之发生了重大变化,单调的办公环境已经不再适应人们的审美能力和
什么是“丹书铁券”“铁券”.外形如筒瓦状的铁制品.民间俗称“免死牌”。其制法为熔铁铸瓦.状如蒸笼,断而为二,形如板瓦.是封建帝王颁发给功臣、重臣的一种带有奖赏和盟约性质的凭
近年来中国轨道交通技术发展迅速,尤其是地铁、重载以及高铁的建设与运营技术更是占据了世界至高点。在铁路运营以及技术装备输出的过程中,安全成为讨论最多,最为重要的方面
案情:石某欲向江某以每套500元的价格购买银行卡,即俗称的“四件套”(一张银行卡、一个U盾、持卡人身份证照片、一张绑定银行卡的手机卡),江某答应,并到农业银行某支行办理一张储蓄
报纸
目的:探讨丙戊酸钠治疗癫痫、持续状态的疗效。方法将2010年1月~2011年12月收治的60例癫痫持续状态急诊患者随机分为治疗组和对照组.每组各30例。其中,治疗组给予丙戊酸钠治疗,对
当前,世界各国都通过出台财政监督相关法律制度来保障和促进财政监督活动正常运行。借鉴国际上的先进做法对于加强我国财政监督制度建设,推动财政监督工作科学发展具有重要意
兖州矿业(集团)公司机电设备制造厂通过对聚晶金刚石复合片(PDC)钻头失效形式及相关工艺过程研究,总结出了PDC钻头加工使用过程应该注意的问题,对提高PDC钻头加工水平及延长使用寿
如何对企业家的高报酬现象做出合理解释,比赛理论在委托代理理论的框架下,从比赛竞争的角度指出,高报酬作为一种显示机制,能激励企业家及其竞争者的努力水平,降低挑选企业家的信息
综述了土壤有机碳分组方面的研究结果,主要讨论了在不同利用类型和不同经营方式条件下,土壤有机碳及其不同组分分解、蓄积、含量、存量、分配比例等的动态变化,为制定合理的经营
针对K中心点算法的初始聚类中心可能过于临近、代表性不足、稳定性差等问题,提出一种改进的K中心点算法。将样本集间的平均距离与样本间的平均距离的比值作为样本的密度参数,