基于无限元2.5维方法预测轨道交通混凝土桥梁低频噪声.docVIP

基于无限元2.5维方法预测轨道交通混凝土桥梁低频噪声 　　摘要：提出了一种基于无限元的2.5维方法，可以快速预测桥梁结构噪声且不失精度。首先，建立沿桥梁纵向截面均匀的二维无限元模型，进行声传递向量计算。然后，联合桥梁三维模态分析结果，通过空间波数变换获得三维桥梁模型的模态声传递向量，并采用三维直接边界元的计算结果验证了所提出的2.5维方法的正确性。最后，进行车轨桥耦合振动计算，将获得的桥梁模态坐标与声模态传递向量结合起来预测桥梁振动辐射噪声。以上海某轨道交通为背景进行噪声现场测试，将现场实测结果与计算结果对比。分析表明，计算结果在近场与实测值吻合较好，但是由于忽略了相邻跨桥梁的影响，在远场数值计算结果小于实测值。 　　关键词：混凝土桥梁;低频噪声;2.5维;无限元法 　　中图分类号：U448.21;TB532 　　文献标志码：A文章编号：1004-4523（2015）06-0929-08 　　DOI：10.16385/j.cnki.issn.1004-4523.2015.06.010 　　引言 　　由于低频噪声对人的集中力、睡眠等有诸多负面影响[1，2]，交通系统产生的低频噪声问题已经引起了较多的关注[3，4]。由于低频噪声很难通过声屏障等措施得到有效的控制，越来越多的学者开始研究桥梁的低频噪声问题[5-11]。为了降低桥梁结构噪声，首先要提出有效的噪声预测方法。很多学者采用三维边界元法（Boundaryelementmethod，简称BEM）来进行结构噪声分析。Ghimire等[5]通过有限元谐响应分析与频域三维声学边界元结合分析了混凝土桥梁的噪声特性。李小珍[6]等通过结合时域车桥耦合振动和频域三维边界元来研究铁路桥梁噪声问题。由于车桥耦合振动是一个时变过程[7]，也有一些学者通过时域声学计算来研究低频桥梁噪声[8，9]，但是采用的振动和声学计算模型都较为简单;而对于较为复杂的模型，采用时域声学计算效率不高。李奇和韩江龙[10-12]提出了采用声模态传递向量的方法预测混凝土桥梁时域和频域的结构噪声，适合对不同车辆和车速下的桥梁辐射噪声进行参数分析。但是，用三维边界元计算较为耗时，以文献[12]中30m跨径的U形梁为例，计算200Hz以下的声模态传递向量需24h。对于更长的桥梁或更高的声学分析频率，三维边界元计算时间将更长。 　　2.5维方法适用于截面沿某一方向均匀的问题，通过对一系列二维声压方程的解进行傅里叶变换得到三维方程的解[13-16]。Duhamel[13]最早采用2.5维边界元方法对截面均匀的声屏障进行声学分析。Salonmons等[14]在研究交通噪声时，用2.5维边界元方法作为标准来验证射线模型的准确性。Hornikx和Forssén[15]将这一方法扩展到城市街道噪声的研究。上述研究主要基于点声源或者线声源的研究，并未直接将振动做为声源考虑。Nillson等[16]提出了一种波导有限元和2.5维边界元相结合的方法来预测钢轨的辐射噪声，钢轨的振动通过波导有限元获得，然后结合2.5维边界元计算得到其辐射声压。Nillson的模型中只考虑了钢轨在单个力作用下的振动辐射噪声，并未考虑车桥耦合效应。 　　由于传统边界元在计算外场噪声时存在特征频率不唯一等奇异问题，且边界元方程为满秩矩阵，计算量随着自由度的增加迅速增大，而无限元法的计算效率要比边界元高很多，因此适合求解较为复杂的结构[17]。本文将在文献[12]的基础上采用2.5维无限元方法来计算桥梁结构噪声。首先，采用2.5维无限元方法（Infiniteelementmethod，简称IFEM）计算了一座U形梁的声模态传递向量，并与三维边界元计算结果对比以验证其正确性。然后，建立车辆-钢轨-桥梁模型进行动力响应分析以获得桥梁的模态坐标。最后采用声模态传递向量方法获得桥梁结构辐射噪声，并与实测结果进行对比。 　　12.5维无限元理论 　　采用一个简单的模型来介绍2.5维无限元理论，三维模型截面沿z轴均匀，如图1所示，三维Helmholtz声波方程在频域内表达[18]式中为声压，k为波数。 　　最终，整个结构振动产生的任一断面z处的总声压可以通过对表面所有点的速度边界产生的声压叠加获得。注意到，当k2-k2z=0时，二维声波方程的解会出现无穷大或者无穷小的情况，此时公式（3）积分公式存在奇异点。为了避开奇异点，本文采用了Hornikx和Forssén[15]提出的数值方法进行求解。首先，将积分区间沿着kz轴分成几个区间，分别求解每个区间的二维解（kz）并对各区间积分，最后将各区间的积分结果相加本文取0.95。 　　2桥梁结构噪声 　　2.1声模态传递向量算法 　　三维桥梁模型的MATVs可以通过上述2.5维无限元方法获得，具体流程如下。首