1. 引言
高速列车进入隧道时产生的压缩波在隧道内近似以当地声速沿列车运行方向向前传播,当压缩波传播至隧道出口时,由于隧道洞口的影响,一部分能量以膨胀波的形式在隧道内反向继续传播,另一部分能量快速向隧道出口外部释放形成微气压波。前人研究成果表明,微气压波属于次声波,通常情况下人耳不会感觉到,主要会引起隧道出口建筑的低频振动 [1] 。但在日本新干线和德国佛罗伦萨–罗马高速铁路线上,某些情况下隧道出口还会出现了人耳可听的音爆的现象。图1给出了隧道出口外微气压波产生过程的示意图。
目前,国内外研究人员针对铁路隧道出口微气压波现象的研究方法主要分为实测分析、理论推导和数值计算三类。在实测方面,日本学者在新干线铁路隧道中做了大量实验测试,得到了丰富的实验数据,实验数据表明隧道出口微气压波噪声主要能量集中在10 Hz以下,属次声声波 [2] ;德国学者M. Hieke [3] [4] 基于Euerwang隧道的测试结果给出了微气压波噪声的评价标准;王学英、高波等 [5] 搭建隧道空气动力学试验模型,分析了缓冲结构对微气压波减缓作用。理论推导方面,Yamamoto [6] 在线性声学理论的基础上,用低频近似和远场近似的方法得到了微气压波的理论计算公式,证明了微气压波幅值与隧道出口位置的压缩波波前梯度最大值成正比,与微气压波的传播距离成反比。数值计算方面,梅元贵 [7] 基于特征线法编写了隧道洞口微气压波的一维计算程序;杨志刚 [8] 用基于高阶谱差分的CAA技术对高速列车在隧道出口产生的微气压波进行二维数值模拟;史宪明 [9] 、张雷 [10] 基于三维有限体积法分别计算了隧道出口微气压波的不同影响因素。
上述研究主要分析了隧道出口微气压波的压力特性及其影响因素,但无法解释隧道外出现人耳可听的音爆现象的原因,更没有相关判定依据。因此,有必要对隧道出口外微气压波声学特性做进一步的分析。本文试图通过混合计算方法分析微气压波声学特性,即根据有限体积法求解隧道出口微气压波流场参数,并基于Lighthill声类比理论,运用有限元法求解隧道出口声场分布,以期对改善高速铁路隧道出口环境噪声提供参考依据。
![](//html.hanspub.org/file/6-2820057x9_hanspub.png)
Figure 1. The generation process of micro-pressure wave at tunnel exit
图1. 高速铁路隧道出口微气压波的产生过程示意图
2. 研究方法
计算微气压波声学特性时,采用基于Lighthill声类比理论的混合噪声计算方法,该方法将气动噪声计算分为流场参数计算和声场分布计算两部分。
2.1. 微气压波流场计算方法
计算微气压波流场参数时,利用动网格技术 [11] 实现列车相对于隧道的运动,并采用结合两方程
紊流模型的有限体积法 [12] 求解三维、粘性、可压缩N-S方程,得到隧道内压缩波传播至隧道出口时的流场参数。此时,为保留微气压波流场脉动参量,需要通过大涡模拟法求解微气压波流场进行模拟。
2.2. 微气压波声场计算方法
2.2.1. Lighthill声类比理论
Lighthill根据流体流动的连续性方程和动量方程,直接推导得到了将流体与声波动相联系的声学波动方程,并被命名为Lighthill方程。Ffowcs Williams和Hawkings仿照Lighthill方程的推导方法,考虑运动边界对噪声的影响,根据广义格林函数给出了物体表面存在非定向运动的声场FW-H方程 [13] :
(1)
FW-H方程右侧A项被视作四极子源项,B项视作偶极子源项,C项视作单极子源项。在计算隧道出口微气压波噪声时,由于微气压波是隧道内压缩波向出口外部空间快速释放产生的,并且当压缩波传播至隧道出口时,并不会引起气体明显的宏观流动,因此计算时可忽略偶极子源的影响;此外,单极子源与壁面的法向气流脉动速度相关,流场计算时将隧道出口壁面看作刚体,不会产生垂直于壁面方向的气流振动,因此计算时同样不考虑单极子源。
2.2.2. 有限元法求解声波方程
由于计算微气压波声学特性时主要关注的是声场分布状态,一般不考虑声波随时间的传递过程,可将简化后的FW-H方程通过分离变量法得到只包含空间变量的Helmholtz方程 [14] :
(2)
采用有限元法求解Helmholtz方程时,进行有限元格式离散后,得到数值形式的声压计算表达式 [15] :
(3)
式中,Qi表示输入的声源向量,Vni表示质点速度边界条件,Pi表示输入的声压边界条件,Fai表示声学激励,pi表示待求解的网格节点声压。
表示方程矩阵,属于稀疏矩阵。
2.2.3. 完美匹配层
能够在满足计算精度的前提下缩小计算域以减少计算工作量,在包含微气压波主要流场参数的空间边界处设置完美匹配层(Perfect Matched Layer,简称PML [16] )。当声波进入PML层后会被吸收,并以指数形式衰减。而远场空间的声辐射情况可通过基于Green函数的Kirchhoff积分计算得到 [17] 。
3. 微气压波声学计算模型的建立
进行隧道出口流场计算时,建立的列车模型参照日本E2系高速列车 [18] ,隧道阻塞比约为0.2。不考虑隧道出入口缓冲结构和山体的影响,对简化后得的列车-隧道模型进行划分网格,如图2所示。经
![](//html.hanspub.org/file/6-2820057x15_hanspub.png)
Figure 2. Partial grid of train-tunnel model at the entrance of tunnel
图2. 列车–隧道模型入口处局部网格
网格无关性验证后,流场计算网格总数量约为430万。
基于三维有限体积法求解N-S方程,得到列车以200 km/h的速度进入长度为500 m的隧道后,距离隧道出口中点不同距离测点处的微气压波波形,如图3所示。微气压波噪声计算分析范围为1.8~2.8 s。
根据隧道出口不同距离处微气压波幅值大小,为尽可能保留微气压波流场参数,设置隧道出口外半径为30 m,并向隧道内部延伸30 m的隧道出口附近空间作为微气压波声学计算域,如图4所示。
对于图5所示的微气压波噪声计算域边界,由于隧道出口实体边界处(隧道内壁面、地面和出口立面)的吸声属性没有实际的测试结果,且根据常用混凝土材料对次声声波吸声性能极低的特性,尝试性的将其定义为全反射边界条件。具体参数大小和其余边界条件类型如表1所示。
4. 隧道出口外部声场分析
当列车进入隧道时生成的压缩波传播至距离隧道出口35 m左右,采用大涡模拟法进行隧道出口微气压波流场瞬态计算。在得到稳定的隧道出口附近流场后,开始提取计算域内空间流场参数的脉动信号。为体现图3所示的隧道出口压力主要变化过程,在满足Nyquist采样定理的最高频率(取200 Hz)与时间步长关系条件下,定义采样时间步长为5 × 10−4 s。连续保存2000个时间步的时域脉动信号,采样总时长为1s。运用有限元法计算得到隧道出口中垂面和水平高度3.5 m处的横截面上微气压波噪声声压级分布情况,如图6所示。
隧道出口微气压波噪声声压分布云图中,次声频段内声压分布在隧道出口外呈现明显的球面形状,且噪声的次声频段集中了大部分能量,1 Hz的微气压波噪声最大声压级达177 dB。
读取隧道出口外隧道中线上高度为3.5 m距离洞口20 m的测点处不同频率对应声压级,如图7中所示。由于100 Hz以上频率噪声声压级低于20 dB,在计算微气压波噪声时,100 Hz以上的微气压波噪声极易被实际情况的环境噪声淹没,在后续计算分析时,主要分析100 Hz以下频率的微气压波噪声。
5. 隧道出口微气压波音爆成因分析
为分析高速铁路隧道出口音爆的发生机理,即判定隧道出口微气压波噪声能否可以被人耳听到,参考GB/T4963-2007《声学标准等响度级曲线》,我们引入人耳听阈曲线作为音爆的判定依据。
由于人耳对不同频率的声音的敏感程度不同,不同频率声音的声压级只有在图8中虚线所示的听阈曲线之上才能被人耳听到,并且同一频率的响度级(单位phon)越高,人耳对该频率声音的感觉越明显(图9)。
对比可知,当列车运行速度为200 km/h时,隧道出口附近20 Hz以上的微气压波噪声声压级位于人耳听阈曲线下方,即该速度条件下生成的微气压波无法被人耳听到。列车运行速度为300 km/h时,对应计算得到的声压级变化曲线会有部分进入到人耳听阈曲线以上的区域,但响度级低于10 phon。而当列车
![](//html.hanspub.org/file/6-2820057x16_hanspub.png)
Figure 3. Micro-pressure wave at the tunnel exit
图3. 隧道出口微气压波
![](//html.hanspub.org/file/6-2820057x17_hanspub.png)
Figure 4. Schematic of Micro-pressure wave calculate domain
图4. 隧道出口微气压波噪声计算域示意图
![](//html.hanspub.org/file/6-2820057x18_hanspub.png)
Figure 5. The acoustic boundaries of domain
图5. 声学计算域边界
![](Images/Table_Tmp.jpg)
Table 1. Boundary condition type and parameter setting of calculate domain
表1. 声学计算域边界条件类型和参数设定
![](//html.hanspub.org/file/6-2820057x19_hanspub.png)
Figure 6. Distribution of SPL of micro-pressure wave at tunnel exit at different frequencies (unit: dB)
图6. 不同频率下隧道出口微气压波噪声的声压级分布情况(单位:dB)
![](//html.hanspub.org/file/6-2820057x20_hanspub.png)
Figure 7. The SPL line of micro-pressure at tunnel exit
图7. 隧道出口外微气压波噪声的声压级曲线
![](//html.hanspub.org/file/6-2820057x22_hanspub.png)
Figure 9. The SPL line of micro-pressure at tunnel exit
图9. 隧道出口附近的微气压波噪声声压级变化曲线
运行速度为400 km/h时,20 Hz以上的部分微气压波频率噪声响度级将达30 phon,已经能被人耳听到。可见,当列车以高于300 km/h的车速进入长度为500 m、阻塞比为0.2时,隧道产生的压缩波传播至隧道出口后,产生的隧道出口外微气压波噪声将可能被人耳听到,即产生隧道出口音爆现象。
6. 结论
本文采用大涡模拟法计算隧道出口流场参数,基于Lighthill声类比理论运用有限元法计算隧道出口声场分布,并引入人耳听阈曲线分析了隧道出口外微气压波音爆现象的发生机理。得到结论如下:
1) 隧道出口微气压波噪声声源可类比为四极子源的发声;
2) 微气压波噪声频率主要集中在20 Hz以下,当频率高于100 Hz时微气压波噪声极易被环境噪声淹没;
3) 在长度为500 m、阻塞比为0.2的隧道中,当列车进入隧道的速度高于300 km/h时,引起的隧道出口外微气压波噪声将会出现人耳可听的音爆现象。
4) 计算分析是基于日本新干线E2车型建立的简化计算模型,若考虑缓冲结构,并且隧道长度和阻塞比的不同时,可能会得到不尽相同的结论。