欢迎来到天天文库
浏览记录
ID:37594406
大小:775.14 KB
页数:7页
时间:2019-05-25
《基于虚拟激励法的结构随机振动声辐射分析》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、固体力学学报第32卷第6期Vol.32No.62011年12月CHINESEJOURNALOFSOLIDMECHANICSDecember2011*基于虚拟激励法的结构随机振动声辐射分析刘宝山1,22**1赵国忠李瑞勇(1中国石油大学(华东)储运与建筑工程学院工程力学系,青岛,266555)(2工业装备结构分析国家重点实验室大连理工大学工程力学系,大连,116024)摘要论文基于有限元法、边界元法和虚拟激励法,对随机激励下结构振动声辐射问题进行研究.提出了一种计算随机激励下结构振动声辐射问题的新方法,其中,有限元法用于计算结构谐
2、振响应,边界元法用于计算结构振动声辐射,虚拟激励法结合有限元和边界元计算随机激励下结构振动声辐射问题.数值算例表明,论文方法在计算精度上与传统方法等价,且更具高效性.关键词声辐射,虚拟激励法,有限元法,边界元法环境影响结构振动响应的问题.随机激励下结构振0引言动声辐射问题是非常复杂的问题,需要引进更为高[8-10]效的计算方法.虚拟激励法是处理随机振动问随机性是大量自然现象的重要规律之一,随机题非常有效的方法,它将平稳随机振动分析转化为激励作用下结构的振动是一种常见的结构振动形简谐振动分析,将非平稳随机振动分析转化为确定式,随机
3、振动的声辐射问题也在生活中处处可见.例性时间历程分析,从而使计算步骤大大简化,却仍保如,汽车行使时受到路面不平引起的振动,飞机飞行持理论的精确解.因此,本文基于虚拟激励法,结合中受到空气流引起的振动,以及轮船在水中受到波有限元法和边界元法对随机激励下结构振动声辐射浪冲击引起的振动,这些振动引起的声辐射问题均问题进行研究.属随机振动的声辐射问题.目前,随机激励下结构振动引起的声辐射问题1随机激励下结构振动的虚拟激励法[1]的研究还不是多见.刘钊和陈心昭提出将边界元方法和统计方法相结合的统计边界元法,并导出了若一个平稳随机激励作用于
4、一线性结构系统,求解随机振动物体辐射噪声声场的计算公式.王秀其方程为:[2,3]峰和陈心昭利用提出的统计边界点法对随机振M¨y+Cy+Ky=fx(t)(1)动结构声辐射的计算进行了研究,计算结果和理论其中M、K和C分别代表结构的质量矩阵、刚度矩解表明该方法计算随机振动结构声辐射问题在相当阵和阻尼矩阵;x(t)是平稳随机激励,其功率谱密宽的振动频率内,能给出良好的计算精度.Allen度为Sxx(ω);f是给定的矢量,代表随机激励分布[4,5]等采用传统的方法,利用有限元法和边界元法情况.根据虚拟激励法,令虚拟简谐激励x珟=研究了
5、随机激励下结构的声辐射问题,进一步提出iωt,并代入式(1):槡Sxx(ω)e了随机激励下结构振动声辐射灵敏度分析方法.纪[6]M¨y+Cy+Ky=fSiωt(2)纲等基于混合有限元法/边界元法,对浸没结构物槡xx(ω)e在随机作用力下的声辐射问题进行随机分析,并开其平稳解为:发了随机作用力下结构-声辐射系统响应的计算程y珘(t)=Y(ω)eiωt(3)[7]序.Chung等讨论了基于有限元法如何确定声学应用其前q阶模态:*国家自然科学基金(10772038,11072049),国家科技重大专项(2009ZX04001-021
6、-02)和中央高校基本科研业务费专项资金(12CX04071A,DUT10LK39)资助.2010-08-02收到第1稿,2010-12-07收到修改稿.**通讯作者.Tel:0411-84708769,E-mail:zhaogz@dlut.edu.cn.·582·固体力学学报2011年第32卷q程为:Y(ω)=∑γjHjj槡Sxx(ω)(4)j=1Gps=HRvxv(t)(15)nn22)-1,γTHj=(ωj-ω+2ijωωjj=jf(5)其中G和H分别为直接边界元法中系数矩阵,它们其中ωj、j、j和γj分别为第j阶的
7、结构圆频率、模均是不对称满秩复数矩阵[11];R是法向速度分布vn态、模态阻尼比和模态分量因子.根据虚拟激励法:矩阵;xv(t)自由度为N的法向速度随机矢量,其功nS(ω)=y珘*y珘T=Y(ω)*Y(ω)T(6)yy率谱密度Svv(ω)是已知的,而它的秩为r,和结构nn其中*代表复共轭,T代表转置.激励的相同,形成r个虚拟激励:将式(4)代入上式并展开得:v珓(t)=deiωt,(j=1,2,…,r)(16)nj槡vnjφvnjqq*TSyy(ω)={∑γjHjj}Sxx(ω){∑γkHkk}其简谐响应:j=1k=1p珟(
8、t)=P(ω)eiωt(17)sjsjqq*T=∑∑γjγkHjHkjkSxx(ω)(7)应用虚拟激励法,其表面声压响应的功率谱密j=1k=1度矩阵为:这就是传统的CQC公式,说明虚拟激励法采用式r(6)进行向量相乘的计算格式,与传统的需要进行矩S(ω)=P
此文档下载收益归作者所有