文章快速检索     高级检索
  空气动力学学报  2021, Vol. 39 Issue (5): 151-161  DOI: 10.7638/kqdlxxb-2021.0101

引用本文  

王田天, 胡冲, 龚彦峰, 等. 扩大斜切式缓冲结构对时速400 km铁路隧道口微气压波缓解研究[J]. 空气动力学学报, 2021, 39(5): 151-161.
WANG T, HU C, GONG Y, et al. Mitigation of micro-pressure wave at 400 km/h railway tunnel exit by oblique enlarged tunnel-hood[J]. Acta Aerodynamica Sinica, 2021, 39(5): 151-161.

基金项目

国家自然科学基金青年&面上项目(51905547,52078199);国家数值风洞工程(NNW);湖南省研究生科研创新项目(CX20200132)

作者简介

王田天(1989-),男,副教授,研究方向:列车空气动力学,结构健康监测等. E-mail:wangtiantian@csu.edu.cn

文章历史

收稿日期:2021-08-03
修订日期:2021-09-05
优先出版时间:2021-10-25
扩大斜切式缓冲结构对时速400 km铁路隧道口微气压波缓解研究
王田天1,2 , 胡冲1 , 龚彦峰3,4 , 杨明智1 , 熊小慧1 , 踪敬良3,4 , 陆意斌1     
1. 中南大学 交通运输工程学院,轨道交通安全教育部重点实验室,长沙 410075;
2. 湖南大学 机械与运载工程学院,长沙 410082;
3. 中铁第四勘察设计院集团有限公司,武汉 430063;
4. 水下隧道技术国家地方联合工程研究中心,武汉 430063
摘要:针对时速400 km高铁隧道口的微气压波研究,建立了基于三维可压缩非定常N−S方程和RNG k-ε 双方程湍流模型的数值模拟方法,并通过动模型试验对计算方法进行了验证。文中分析了隧道长度、缓冲结构(带斜切的扩大段)长度和缓冲结构开孔数对微气压波的影响。研究结果表明:在一定隧道长度内,微气压波幅值随着隧道长度的增加逐渐增大;隧道口设置带斜切(30°)的等截面扩大段缓冲结构能有效缓解微气压波幅值,且随着缓冲结构的长度增加,微气压波幅值的缓解率先增大后减小,其中缓冲结构长88.56 m时缓解效果最佳,缓解率可达59.2%;该缓冲结构开孔后能进一步缓解微气压波幅值,缓冲结构开2个孔时的缓解效果最佳,缓解率可达70.9%。本文提出的缓冲结构能使得5 km及以下长度隧道的微气压波幅值达到国家规定标准,研究成果对时速400 km/h高速铁路隧道口缓冲结构的设计具有参考价值。
关键词时速400 km    高铁隧道    微气压波    缓冲结构    
Mitigation of micro-pressure wave at 400 km/h railway tunnel exit by oblique enlarged tunnel-hood
WANG Tiantian1,2 , HU Chong1 , GONG Yanfeng3,4 , YANG Mingzhi1 , XIONG Xiaohui1 , ZONG Jingliang3,4 , LU Yibin1     
1. Key Laboratory of Traffic Safety on Track, Ministry of Education, School of Traffic & Transportation Engineering, Central South University, Changsha 410075, China;
2. College of Mechanical and Vehicle Engineering, Hunan University, Changsha 410082, China;
3. China Railway Siyuan Survey and Design Group Co., Ltd, Wuhan 430063, China;
4. National & Local Joint Engineering Research Center of Underwater Tunnel Technology, Wuhan 430063, China
Abstract: Numerical simulations based on the three-dimensional, compressible, unsteady Navier-Stokes equations and the RNG k-ε two-equation turbulence model are used to study the micro-pressure generated by a high-speed train passing through a tunnel at the speed of 400 km/h, and the numerical results are further verified by the moving model test. The effects on the micro-pressure wave at the tunnel exit by the tunnel length, the tunnel-hood length (enlarged section with oblique tunnel portal) and the number of openings holes on the tunnel-hood, are analyzed. Results show that: when the tunnel length is in a certain range, the amplitude of the micro-pressure wave gradually increases with the increase of the tunnel length. The tunnel-hood with an enlarged section of equal cross-section with oblique tunnel portal (30°) at the tunnel entrance can effectively reduce the amplitude of the micro-pressure wave, the alleviative rate of the amplitude of the micro-pressure wave first increases then decreases with the increase of the tunnel-hood length, and the best mitigation effect is achieved when the length of the tunnel-hood is 88.56 m, with an alleviative rate up to 59.2%. The tunnel-hood can further alleviate the amplitude of the micro-pressure wave by opening holes, and the best mitigation effect is achieved when two holes are opened on the tunnel-hood, which results in an alleviative rate of 70.9%. The tunnel-hood proposed in this study can make the amplitude of the micro-pressure wave at the exit of a tunnel with a length less than 5 km meet the national standard. The present findings provide important guidance for the design of tunnel-hoods used in the 400 km/h high-speed railway.
Keywords: 400 km/h    high-speed railway tunnel    micro-pressure wave    tunnel-hood    
0 引 言

随着高速铁路的不断发展,中国高铁的最高运营速度已达350 km/h,但由于旅客出行需求和列车运行能力的矛盾依然突出,国家已开始研发设计时速400 km的高速列车。由于我国地形复杂,山地占国土面积的2/3,在高速铁路建设的过程中需修建大量铁路隧道,而车速的提高不可避免地会加剧车-隧耦合气动效应[1-2]。如图1所示,高速列车鼻尖进入隧道的瞬间,会在车前产生初始压缩波,该压缩波以声速沿着隧道长度方向传播并不断激化,到达出口时会以脉冲形式从隧道出口向外辐射,形成脉冲状的压力波,称之为隧道口微气压波或简称为微气压波[3-5]。微气压波幅值随着列车速度的增加而显著增大,当列车以时速小于250 km通过短隧道时,微气压波幅值基本与列车速度的三次方成正比[6];而随着车速和隧道长度的增加,微气压波幅值可与3倍以上的车速成正比[7]。过大的脉冲状压力不仅对隧道口附近的环境产生不利的影响,同时也会严重影响附近居民的正常生活[8-9]。因此,有效缓解时速400 km高速铁路隧道口微气压波对于铁路隧道的建设以及保障隧道口周围环境和居民安全具有重要的意义。


图 1 微气压波产生示意图 Fig.1 Schematic diagram of the generation of the micro-pressure wave

在近几十年间,国内外对缓解微气压波的方法进行了大量研究。降低微气压波的基本思路是在压缩波形成和传播阶段减小其压力梯度,现有的缓解方法主要分为两类,一是改变隧道结构(比如在隧道口加设缓冲结构或增大隧道断面面积或修建辅助坑道等)[6];二是改变列车头部形状(增加列车头部流线型长度和改变列车断面面积),从而达到减小初始压缩波最大压力梯度的效果[10]。Yamamoto[11]构建了一个称为“辐射立体角模型”(RSA模型)的预测模型,利用隧道出口周围的立体角来预测地形对微气压波的影响,后来Miyachi[12]扩展了山本的RSA模型,生成了微气压波的声学模型,该模型同时考虑了二阶声源和隧道口周围地形的声学效应。Nagai[13]和Sanetoshi[14]通过数值模拟和动模型试验研究了缓冲结构形式以及优化缓冲结构断面对压力梯度的影响,结果发现优化后能显著降低初始压缩波的压力梯度。贺旭洲等[15]利用数值模拟对斜切式洞门缓冲结构开口率进行了优化分析,并给出了不同车速下的最优开口率。Zhang等[16]使用1/20比例的动模型试验研究了倾斜的隧道口对微气压波的影响。Wang等[17]研究了改变隧道断面参数对压力波动的影响。此外,有学者对缓冲结构缓解微气压波的机理进行研究。周丹[18]认为等截面扩大段缓冲结构使列车前空气在进入隧道前提前被压缩,从而增大了压力上升到最大值的时间;且缓冲结构将压力上升分为了两个或三个阶段,从而压力梯度峰值得到减小。郑长青[19]和牛纪强[20]等认为开口式缓冲结构对初始压缩波幅值基本没有缓解效果,但是可以增加压力上升的时间,从而降低压力梯度达到缓解微气压波的效果。

近几年开始有学者对更高速度下列车隧道空气动力学问题进行研究。吴剑等[21]开展了时速 300 km至350 km高速铁路双线隧道微气压波激化作用及缓解措施研究。范胜利等[22]研究了一种带有排气管的扩大等截面式的新型缓冲结构,且给出了排气管的参数限制,认为排气管的长度应小于压缩波前长度的一半;单个排气管的断面积不宜超过缓冲结构断面面积的20%;并可以针对不同的车速调节排气管的长度或改变排气开孔的位置或数量以达到有效减缓微气压波的目的,但此新型缓冲结构是在车速360 km/h下进行的研究,其对更高车速下微气压波的缓解效果的普适性需进行深入研究,而本研究车速为400 km/h,且缓冲结构带有斜切和开孔,气动效应也更为复杂。魏雨生等[23]对高速列车以时速400 km通过隧道时的列车和隧道压力变化进行了数值模拟,发现车体表面和隧道壁面最大压力变化幅值比300 km时分别提升了90.4%和65.3%,但未对隧道口的微气压波进行研究。

综上,目前国内外针对高速列车以时速400 km通过隧道时的微气压波的缓解方法研究较少,存在较大空缺。故本文以时速400 km高速列车经过隧道引发的微气压波问题为导向,对隧道长度、缓冲结构长度和缓冲结构开孔个数等多个参数展开耦合研究,以期为缓解时速400 km高速铁路隧道口微气压波的缓冲结构设计提供参考。

1 数值计算方法 1.1 数值计算模型

高速列车模型如图2所示,由于微气压波主要与车头的形状相关,故选用2车编组(头车 + 尾车)的复兴号高速列车组(型号:复兴号400-AF)为研究对象,同时对转向架、受电弓、门窗等基本不影响计算准确性的结构进行简化,保留列车的典型气动外形。我国现有运营速度350 km/h的高速列车,主要在京沪、京张等线路上运行,其隧道选用断面积为100 m2的典型高铁双线隧道,两列车的线间距为5 m[24],如图3所示。本文选取此隧道断面进行研究。由于相关研究表明缓冲结构斜切角为30°、扩大段面积与隧道断面积比约为2∶1时,对微气压波有较好的缓解效果[16]。故本文选取了斜切角为30°的等截面扩大段(扩大段面积S = 200 m2)的缓冲结构为研究对象,分析其长度和开孔率对时速400 km高速铁路隧道口微气压波的缓解效果。如图4所示,缓冲结构斜切部分的长度为22.56 m,有研究表明列车以一定速度通过隧道,当等截面缓冲结构的长度超过4倍的隧道等效水力直径时,微气压波幅值不再随缓冲结构长度的增加而变化[18]。本文模型隧道的等效水力直径为11.6 m,且缓冲结构带有斜切和开孔,列车速度更高,所以本研究缓冲结构的总长度L 在本研究中设置4种,分别为68.56 m、78.56 m、88.56 m和98.56 m;每个开孔的面积为6 m × 4 m,两孔间间隔均为4 m,不同开孔数的工况设置为:1个开孔时为序号①开孔,2个开孔时为①②开孔,3个开孔时为①②③开孔,4个开孔时为①②③④开孔,5开孔时为①②③④⑤开孔。


图 2 两车编组高速列车模型 Fig.2 High-speed train model of a two-car marshalling


图 3 100 m2双线隧道断面 Fig.3 Cross-section of a 100 m2 double-track tunnel


图 4 研究缓冲结构长度和开孔个数影响时所采用的几何模型(S = 100 m2S = 200 m2 Fig.4 Geometric model used to study the influence of the tunnel-hood length and thenumber of opening holes (S = 100 m2S = 200 m2)
1.2 计算区域和边界条件

计算区域如图5所示,分为静止区域和滑移区域2个部分,静止区域又包含列车进口大气区域、隧道(长度L0因工况而定)和列车出口大气区域。考虑到流场和尾流扰动的充分发展,从隧道的两端分别向两侧各延伸450 m。计算域中列车进口大气区域的地面和与隧道接触的面定义为无滑移壁面,其他边界定义为压力出口;在列车出口大气区域中,地面和与隧道接触的面定义为无滑移壁面,其余边界均设置为压力远场边界;隧道和列车均设置为壁面边界条件[25];包含列车的滑移区域的两端分别定义为压力入口和压力出口边界,且在滑移区域和静止区域之间设置交换面以实现数据交换。另外,值得注意的是,由于本文中的流动为三维可压缩气体的不等熵流动,故所有壁面边界条件均为绝热壁面,所有边界初始温度均设置为300 K。


图 5 计算区域和边界条件(单位:m) Fig.5 Computational domain and boundary conditions (unit: m)
1.3 网格划分

图6所示,本研究采用的高速列车外形复杂,且在隧道口设置了带斜切的缓冲结构,因此选择混合网格对计算域进行离散。列车周围区域及斜切缓冲结构的部分区域离散为非结构化网格,其余部分离散为结构化网格,列车表面上的最小网格尺寸为0.05 m,隧道表面最小网格尺寸为0.15 m。为了对网格的无关性进行检验,按上述的方法调整网格参数,将下文2.1节中的动模型试验的列车和隧道划分为粗、中、细三种网格(其中中网格划分策略与本文选用的网格划分策略一致),网格单元总数分别为4.32×106、9.74×106、1.728×107。选取列车中部测点和隧道中部测点的压力峰-峰值与实验结果进行比较,三种网格的网格密度和压力峰-峰值的偏差结果列于表1。由表1可知,粗网格计算得到的列车表面和隧道壁面的压力峰-峰值较实验的偏差分别为−2.9%和−3.6%,而中等数量网格的列车和隧道的偏差分别为1.2%和−1.6%,细网格的偏差与中等网格基本一致,可以看出粗网格计算的结果与对比实验偏差较大。进一步细分网格,当网格数量超过9.74×106时,计算结果随网格数量增加变化较小,证明网格已经收敛,同时说明9.74×106的中等数量的网格可以同时满足计算精度和计算效率两方面的要求,所以本文选取中等网格划分策略进行研究。


图 6 计算网格分布 Fig.6 Computational grid distribution

表 1 网格密度及峰-峰值的偏差 Table 1 Grid density and the peak-to-peak value deviation
1.4 计算方法及求解设置

本研究选用的求解软件为基于有限体积法的ANSYS FLUENT求解器。高速列车通过隧道时,由于空气受到隧道和车体壁面的束缚,需考虑空气的压缩性[26-27]。因此,本文选用基于可压缩 N-S方程的RNG k-ε湍流模型,并采用滑移网格技术,对高速列车通过隧道时的复杂空气流场进行求解[28-29]。速度和压力的耦合方式选用SIMPLE 算法,基于格林高斯的单元法来控制梯度变化,时间离散采用二阶隐式格式,空间离散采用二阶迎风格式[10]。时间步长设置为0.0075 s,并设置50次内部迭代[16]

1.5 测点布置

因我国对微气压波的评判标准为:距隧道口20 m处的微气压波要小于50 Pa,距隧道口50 m处的微气压波要小于20 Pa[30]。故隧道口微气压波测点分布如图7所示,在距离隧道口20 m、50 m的隧道中线两侧分别对称布置2个测点,靠近列车一侧的测点编号为1和2,远离车一侧的测点编号为3和4,所有测点距离轨面高为1 m。


图 7 微气压波测点的布置 Fig.7 Layout of the micro-pressure wave measurement points
2 数值方法验证 2.1 动模型实验

本次动模型试验在中南大学轨道交通安全教育部重点试验室的动模型试验平台进行,该平台可用于单轨试验和双轨试验,主要由三部分组成,分别用于加速、测试和制动[31]。该平台拥有先进的列车穿越隧道和列车交会实验模拟研究系统,并获得CNAS国际认证(注册号CNAS L 10220)和CMA中国国家检测资质认证(证书号170021002479),能够对本文所采用的数值计算方法进行可靠性评价。

图8,动模型试验采用1∶20缩比模型,列车模型长度为2585 mm,高度为194.5 mm,宽度为163.3 mm;隧道模型对应的实际净空面积为100 m2,长164 m,线间距0.25 m。本次动模型试验的速度为400 km/h,在此平台上一共进行了15次弹射,其中发射速度不在目标速度1%以内的工况不被选用。为记录列车过隧道的产生的压力变化,如图9,在隧道壁面上布置5个测点分别记为S1-S5,其中S3为靠近列车一侧的测点;列车表面布置5个测点,记为T1-T5,如图10


图 8 动模型试验平台 Fig.8 Moving model test platform


图 9 隧道壁面测点(单位:m) Fig.9 Measurement points on the tunnel surface (unit: m)


图 10 列车表面测点(单位:m) Fig.10 Measurement points on the train surface (unit: m)
2.2 结果分析对比

按照上文的数值方法对试验模型进行网格划分和模拟计算。由于动模型试验采用的是1∶20的比例,所以试验时间要进行相应转换,同时由于雷诺数大于3.6×105,模型尺寸对列车周围的流场影响很小,因此试验测点的压力值可以代表全尺寸模型的压力值[32]图11(a)和图11(b)分别为隧道壁面测点S3和列车表面测点T1的压力时程曲线,发现与动模型试验相比,两测点的压力峰-峰值的偏差分别为−1.4% 和1.1%,证明本文的数值方法能够正确反映压力极值的变化规律。图11(c)为图11(a)中隧道壁面测点初始压缩波的压力梯度时程曲线,可以看出数值模拟和动模型试验的初始压缩波梯度曲线吻合较好,且两者压力梯度最大值的偏差仅为1.2%。由式(1)可知微气压波的压力幅值∆pw与初始压缩波的最大压力梯度 $\Bigg[\dfrac{\mathrm{d}p}{\mathrm{d}t}\Bigg]_{\rm{ex}}$ 成正比[1],所以可以得出数值模拟与动模型试验的微气压波幅值的偏差也为1.2%,这证明本文采用的数值方法具有较高的精度。


图 11 动模型实验和数值模拟所得压力时程曲线对比 Fig.11 Pressure time history comparison between the moving model test and the numerical simulation
$ \Delta {p_w} = \frac{{2S}}{{\varOmega ad}}{\left[ {\frac{{{\text{d}}p}}{{{\text{d}}t}}} \right]_{{\text{ex}}}} $ (1)
3 计算结果和讨论 3.1 无缓冲结构的隧道长度对微气压波的影响

截至2020年底,中国已投入运营的高速铁路隧道共3631座,总长约 6003 km,其中10 km以上的特长隧道只有87座[33],且赵勇等2017的调查结果显示中国高速铁路隧道长度在3 km以下的数量占高速铁路隧道总数的80%以上[34]。可见我国高速铁路隧道绝大多数都在3 km以内,5 km以上的隧道大部分位于群山峻岭,且对应的列车设计时速较小。因此本文选取1 km、2 km和5 km的隧道长度进行研究。

图12为时速400 km高速列车通过长度1 km隧道时,距隧道口20 m处同一截面上对称的两个测点(测点1和测点3)和距隧道口50 m处截面上对称的两个测点(测点2和测点4)的微气压波曲线图,可以发现同一截面上对称的两测点曲线基本吻合,这说明微气压波在同一截面上相同高度的强度相同。所以本文选取靠近列车一侧的测点1(距隧道口20 m)和测点2(距隧道口50 m)进行分析。


图 12 隧道口微气压波曲线图 Fig.12 Micro-pressure wave curves at the tunnel exit

图13(a、b)为高速列车以时速400 km通过长度(L0)为1 km、2 km、5 km的无缓冲结构隧道时,测点1和测点2的微气压波随时间变化的曲线图。可以发现隧道长度在5 km以内时,微气压波幅值随着隧道长度的增加而增大,这是因为初始压缩波在隧道内传播时,由于挤压和摩擦效应,空气密度和温度随之增加,引起声速的提高,使初始压缩波后部比前端传播的更快进而使其激化所导致。


图 13 不同隧道长度下微气压波曲线图 Fig.13 Micro-pressure wave curves for different tunnel lengths

表2表3分别为高速列车通过不同长度隧道时,测点1和测点2的微气压波幅值及相对1 km隧道的增长率。可以发现,在每个测点隧道长5 km的微气压波幅值及相对隧道长1 km的增长率都为最大,测点1和测点2微气压波幅值分别为165.7 Pa和68.3 Pa,微气压波幅值最大增长率分别为27.4%和28.4%,均超过国家相关标准。故下文开展斜切(30°)等截面扩大段式缓冲结构对微气压波的缓解研究,分析其长度和开孔数对微气压波的影响。

表 2 不同隧道长度下测点1微气压波幅值及增长率 Table 2 Amplitude and growth rate of the micro-pressure wave at measurement point 1 for different tunnel lengths

表 3 不同隧道长度下测点2微气压波幅值及相对增长率 Table 3 Amplitude and growth rate of the micro-pressure wave at measurement point 2 for different tunnel lengths
3.2 带斜切(30°)的缓冲结构长度对微气压波的影响

图14(a、b)为高速列车以时速400 km/h通过不同长度的缓冲结构隧道(隧道长1 km)时,测点1和测点2的微气压波随时间的变化曲线,其中缓冲结构长度L分别为68.56 m、78.56 m、88.56 m、98.56 m,且选用无缓冲结构的工况作为对比。可以发现与无缓冲结构相比,每个长度的缓冲结构对微气压波都有较好的缓解效果,最优缓冲结构长度为88.56 m。表4表5分别为测点1和测点2在不同长度缓冲结构下微气压波幅值及其相对无缓冲结构的缓解率,可见缓冲结构长88.56 m时,测点1和测点2的微气压波幅值的最大缓解率分别为59.2%和56.0%。图15为两个测点微气压波幅值随缓冲结构长度变化的曲线,可以看出微气压波幅值随着缓冲结构长度的增加先减小后又略微增大。


图 14 不同缓冲结构长度下隧道口微气压波曲线图 Fig.14 Micro-pressure wave curves at the tunnel exit for different tunnel-hood lengths

表 4 不同缓冲结构长度下测点1微气压波幅值及相对缓解率 Table 4 Amplitude and alleviative rate of the micro-pressure wave at measurement point 1 for different tunnel-hood lengths

表 5 不同缓冲结构长度下测点2微气压波幅值及相对缓解率 Table 5 Amplitude and alleviative rate of the micro-pressure wave at measurement point 2 for different tunnel-hood lengths


图 15 微气压波幅值随缓冲结构长度的变化 Fig.15 Variation of the amplitude of the micro-pressure wave with the tunnel-hood length

分析其原因,在周丹[18]的研究中可知:当缓冲结构长度在一定倍数的隧道等效水力直径之内时,微气压波幅值随缓冲结构长度的增加而减小;当其长度超过一定倍数的隧道等效水力直径后,微气压波基本不变,此长度定为临界长度Lh。其研究也发现了不同的车速对应不同的Lh,当缓冲结构长度超过各自速度对应的Lh时,速度越大,在该长度之后对应的微气压波增大的越明显,这与本文结果相符。在最优长度的缓冲结构下(88.56 m),测点1和测点2的微气压波幅值分别为53.1 Pa和23.4 Pa,仍不满足国家标准,故需进一步对微气压波进行缓解。

3.3 缓冲结构开孔个数对微气压波的影响

由上文知隧道长1 km时,不开口缓冲结构长为88.56 m对微气压波缓解效果最好,但仍不能达到标准要求,所以下文研究隧道长1 km时,斜切(30°)缓冲结构(L = 88.56 m)的开孔率对微气压波的影响,以寻求更好的缓解效果。

图16为高速列车以时速400 km通过不同开孔数的缓冲结构隧道时,测点1和测点2的微气压波随时间的变化曲线,可以发现相对于0开孔的缓冲结构,开1孔和开2孔可以使微气压波幅值减小,开3、4、5孔会使微气压波幅值增大,且微气压波曲线中的第一个峰值随着开孔数的增加向后移动。这是因为开孔的位置是沿着列车运行方向纵向排列的,随着列车的运行,开孔处发生泄气,延长了初始压缩波形成的时间,从而引起微气压波达到幅值的滞后。


图 16 不同开口数的缓冲结构下微气压波曲线图 Fig.16 Micro-pressure wave curves for tunnel-hoods with different numbers of opening holes

表6表7分别为测点1和测点2在不同开孔数的缓冲结构下微气压波幅值和相对无缓冲结构的缓解率。可见测点1和测点2具有相同的变化规律:微气压波幅值随开孔数的增加先减小后增大,在开孔数为2时微气压波幅值最小,这与Liu[35]等发现的“微气压波幅值随缓冲结构开孔率的增加先减小后增大”的规律基本相符,但其研究的最优开孔率为0.018,本文两个开孔对应的开孔率为0.0213,这主要是因为本研究的车速为400 km/h,缓冲结构带有斜切,且缓冲结构长度面积、隧道长度面积等都与文献[35]研究的不同,但整体变化规律是一致的。下面分析此变化规律的原因:孔①②③④⑤沿缓冲结构长度方向纵向排列,当开孔数较少时(2个开孔以内),列车进入缓冲结构被压缩的空气会在开孔处发生泄气,延长了初始压缩波最大值形成的时间,从而初始压缩波的最大压力梯度下降,且2个开孔时泄气量大于1个开孔,所以2开孔的初始压缩波最大压力梯度小于1开孔;当开孔数大于2时,因为①②开孔已经分散了大部分压缩气体,随着列车的运动,在开孔处会有进气和出气的复合作用,且在开孔处和扩大段与隧道交接处产生的复杂小波会相互叠加,在上述因素的共同作用下,相对于2开孔时,初始压缩波的最大压力梯度会增大。测点1和测点2微气压波幅值分别为37.9 Pa和18.1 Pa,相对无缓冲结构时的缓解率分别为70.9%和66.0%,都已满足标准要求。

表 6 不同开孔数的缓冲结构下测点1微气压波幅值及相对缓解率 Table 6 Amplitude and alleviative rate of the micro-pressure wave at measurement point 1 for tunnel-hoods with different number of opening holes

表 7 不同开孔数的缓冲结构下测点2微气压波幅值及相对缓解率 Table 7 Amplitude and alleviative rate of the micro-pressure wave at measurement point 2 for tunnel-hoods with different numbers of opening holes

图17(a)为采用斜切30°长88.56 m开2孔的等截面扩大段缓冲结构,隧道长分别为1 km、2 km、5 km时压缩波变化的曲线,图17(b)为对应压力梯度的变化。可以看出,压力幅值和压力梯度的幅值都受到了隧道在长度方向的激化,随隧道长度的增加而增加,这与3.1节中得到的微气压波受隧道长度的激化的规律一致。另外,Yamamoto[36]利用远场和低频的近似方法得到了微气压波幅值与压缩波的最大压力梯度成正比的关系,所以微气压波幅值与压力梯度幅值变化规律具有一致性。同时,本文旨在研究400 km时速下的隧道口微气压波的缓解措施,以使微气压波达到相关标准,所以直接监测隧道口的微气压波是可取的。


图 17 不同隧道长度下压缩波及压力梯度的变化 Fig.17 Variation of the compression wave and pressure gradient for different tunnel lengths

图18为隧道长1 km、2 km、5 km时,采用斜切30°长88.56 m开2孔的等截面扩大段缓冲结构,在测点1和测点2的微气压波幅值的曲线,表8是对应的微气压波幅值。可见采用上述缓冲结构可将长1 km、2 km、5 km隧道的微气压波幅值都降到标准以下。故本文提出的斜切角为30°、扩大段长度为88.56 m、开2孔的缓冲结构能使长为5 km及以下长度的时速400 km高速铁路隧道口微气压波缓解达到国家相关标准。


图 18 缓冲结构开2孔时微气压波幅值随隧道长度的变化 Fig.18 Variation of the amplitude of the micro-pressure wave with the tunnel length for the tunnel-hood with 2 opening holes

表 8 缓冲结构开2孔时不同隧道长度的微气压波幅值 Table 8 Amplitude of the micro-pressure wave for different tunnel lengths with 2 opening holes on the tunnel-hood
4 结 论

本文通过数值模拟和动模型实验,对时速400 km高速铁路隧道口微气压波进行分析,研究隧道长度、斜切角为30°的缓冲结构长度和开孔数对微气压波的影响。主要结论如下:

1)隧道长度在5 km以内时,微气压波幅值随着隧道长度的增加而增大,其中隧道长5 km与1 km相比,在距隧道口20 m和50 m处微气压波幅值分别增大了27.4%和28.4%;

2)微气压波幅值随着缓冲结构长度的增加先减小后又增大,随着缓冲结构开孔数的增加微气压波幅值先减小后增大,存在最优的缓冲结构长度和开孔数,最优长度为88.56 m,最优开孔数为2开孔;

3)最优长度与最优开孔数组合的缓冲结构,能使隧道长为1 km的微气压波幅值在距隧道口20 m和50 m处分别缓解70.9%和66.0%;并能使长为5 km及以下长度的时速400 km高速铁路隧道口微气压波缓解达到国家相关标准。

参考文献
[1]
田红旗. 列车空气动力学[M]. 北京: 中国铁道出版社, 2007: 273-281.
TIAN H Q. Train aerodynamics[M]. Beijing: China Railway Publishing House, 2007: 273-281. (in Chinese)
[2]
TIAN H Q. Review of research on high-speed railway aerodynamics in China[J]. Transportation Safety and Environment, 2019, 1(1): 1-21. DOI:10.1093/tse/tdz014
[3]
刘堂红, 田红旗, 金学松. 隧道空气动力学实车试验研究[J]. 空气动力学学报, 2008, 26(1): 42-46.
LIU T H, TIAN H Q, JIN X S. Experimental study of full-scale train on aerodynamics in tunnel[J]. Acta Aerodynamica Sinica, 2008, 26(1): 42-46. DOI:10.3969/j.issn.0258-1825.2008.01.008 (in Chinese)
[4]
MIYACHI T. Non-linear acoustic analysis of the pressure rise of the compression wave generated by a train entering a tunnel[J]. Journal of Sound and Vibration, 2019, 458: 365-375. DOI:10.1016/j.jsv.2019.06.033
[5]
何旭辉, 吉晓宇, 敬海泉, 等. 高速铁路全封闭声屏障列车压力波和微气压波数值模拟研究[J/OL]. 空气动力学学报: 1-14[2021-05-11].
HE X H, JI X Y, JING H Q, et al. Numerical simulation of the pressure waves and micro pressure waves of an enclosed sound barrier when a high-speed train passing through[J/OL]. Acta Aerodynamica Sinica: 1-14[2021-05-11].
[6]
马伟斌, 张千里, 程爱君, 等. 高速铁路隧道洞口微气压波影响因素与变化规律研究[J]. 铁道学报, 2013, 35(5): 97-102.
MA W B, ZHANG Q L, CHENG A J, et al. Study on influence factors and changing law of micro-pressure waves at tunnel portals of high-speed railways[J]. Journal of the China Railway Society, 2013, 35(5): 97-102. DOI:10.3969/j.issn.1001-8360.2013.05.016 (in Chinese)
[7]
赵有明, 马伟斌, 程爱君. 高速铁路隧道气动效应[M]. 北京: 中国铁道出版社, 2012: 152-157.
[8]
刘金通, 李人宪. 高速铁路隧道出口微气压波声学特性计算分析[J]. 声学与振动, 2018, 6(2): 62-69.
LIU J T, LI R X. Acoustic characteristics of micro-pressure wave at the high speed railway tunnel exit[J]. Open Journal of Acoustics and Vibration, 2018, 6(2): 62-69. DOI:10.12677/ojav.2018.62008 (in Chinese)
[9]
任杰. 更高速度下铁路隧道洞口微气压波特性研究[D]. 兰州: 兰州交通大学, 2018: 1-5.
REN J. Study on micro-pressure wave noise characteristics of railway tunnel portal at higher speed[D]. Lanzhou: Lanzhou Jiatong University, 2018: 1-5. (in Chinese)
[10]
ZHANG G, KIM T H, KIM D H, et al. Prediction of micro-pressure waves generated at the exit of a model train tunnel[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2018, 183: 127-139. DOI:10.1016/j.jweia.2018.10.015
[11]
YAMAMOTO A. Micro-pressure wave radiated from tunnel exit[C]//Preprint of the Spring Meeting of Physical Society of Japan, 1977: 137.
[12]
MIYACHI T. Acoustic model of micro-pressure wave emission from a high-speed train tunnel[J]. Journal of Sound and Vibration, 2017, 391: 127-152. DOI:10.1016/j.jsv.2016.09.031
[13]
NAGAI S, FUJITA H, IIDA M, et al. Study of micro-pressure wave radiated from tunnel exit with hood[C]//Proceedings of the 9th Annual Meeting of the Japan Society of Mechanical Engineers Kanto Branch, 2003: 425-426. doi:10.1299/jsmekanto.2003.9.425
[14]
SAITO S. Optimizing cross-sectional area of tunnel entrance hood for high speed rail[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2019, 184: 296-304. DOI:10.1016/j.jweia.2018.11.028
[15]
贺旭洲, 王英学, 付业凡, 等. 斜切式洞门缓冲结构开口率的优化分析[J]. 空气动力学学报, 2012, 30(2): 202-205.
HE X Z, WANG Y X, FU Y F, et al. Bamboo-truncating portal hood optimization of aperture ratio analysis[J]. Acta Aerodynamica Sinica, 2012, 30(2): 202-205. DOI:10.3969/j.issn.0258-1825.2012.02.011 (in Chinese)
[16]
ZHANG L, LIU H, STOLL N, et al. Influence of tunnel aerodynamic effects by slope of equal-transect ring oblique tunnel portal[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2017, 169: 106-116. DOI:10.1016/j.jweia.2017.07.011
[17]
WANG T T, WU F, YANG M Z, et al. Reduction of pressure transients of high-speed train passing through a tunnel by cross-section increase[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2018, 183: 235-242. DOI:10.1016/j.jweia.2018.11.001
[18]
周丹. 长大隧道、隧道群空气动力效应算法研究及应用[D]. 长沙: 中南大学, 2007.
ZHOU D. Research on the long tunnel and tunnel group's aerodynamic algorithm and its application[D]. Changsha: Central South University, 2007. (in Chinese)
[19]
郑长青. 高速铁路单线隧道断面优化及缓冲结构研究[D]. 成都: 西南交通大学, 2006.
ZHENG C Q. Research on optimization of the cross-section of the single-track high speed railway tunnel and tunnel hood[D]. Chengdu: Southwest Jiaotong University, 2006. (in Chinese)
[20]
牛纪强, 梁习锋, 周丹, 等. 明洞式隧道洞门开口率优化[J]. 哈尔滨工业大学学报, 2017, 49(3): 175-180, 188.
NIU J Q, LIANG X F, ZHOU D, et al. Optimization on the opening rate of the open-cut tunnel portal[J]. Journal of Harbin Institute of Technology, 2017, 49(3): 175-180, 188. DOI:10.11918/j.issn.0367-6234.2017.03.028 (in Chinese)
[21]
吴剑, 史宪明, 万晓燕. 时速300~350 km高速铁路双线隧道微气压波激化作用及缓解措施研究[J]. 土木工程学报, 2017, 50(S2): 209-214.
WU J, SHI X M, WAN X Y. Study on intensification and mitigation methods of micro pressure wave of double track tunnel in 300 to 350 km/h high speed railway[J]. China Civil Engineering Journal, 2017, 50(S2): 209-214. (in Chinese)
[22]
范胜利, 张慧玲, 朱勇, 等. 新型高速铁路隧道洞口微气压波减缓设施研究[J]. 四川建筑, 2020, 40(3): 73-75.
FAN S L, ZHANG H L, ZHU Y, et al. Research on a new type of facility for mitigating micro-pressure waves at the entrance of high-speed railway tunnels[J]. Sichuan Architecture, 2020, 40(3): 73-75. DOI:10.3969/j.issn.1007-8983.2020.03.028 (in Chinese)
[23]
魏雨生, 刘峰, 张翠平, 等. 高速列车以400 km/h通过隧道时气动效应的数值模拟与分析[J]. 中国科技论文, 2019, 14(6): 652-656, 704.
WEI Y S, LIU F, ZHANG C P, et al. Numerical simulation and analysis on aerodynamic effects of high-speed train passing through tunnels with 400 km/h[J]. China Science-paper, 2019, 14(6): 652-656, 704. DOI:10.3969/j.issn.2095-2783.2019.06.013 (in Chinese)
[24]
赵文成. 京沪高铁隧道洞门气动特性分析[J]. 路基工程, 2013(3): 114-116, 121.
ZHAO W C. Analysis on aerodynamic characteristics of tunnel portal at Beijing-Shanghai express railway[J]. Subgrade Engineering, 2013(3): 114-116, 121. DOI:10.3969/j.issn.1003-8825.2013.03.027 (in Chinese)
[25]
林鹏, 刘冬雪. 城际列车底部结构优化减小气动阻力研究[J]. 空气动力学学报, 2018, 36(5): 763-773.
LIN P, LIU D X. Reduction of aerodynamic drag on intercity railway through car bottom structure optimization[J]. Acta Aerodynamica Sinica, 2018, 36(5): 763-773. DOI:10.7638/kqdlxxb-2016.0176 (in Chinese)
[26]
RICCO P, BARON A, MOLTENI P. Nature of pressure waves induced by a high-speed train travelling through a tunnel[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2007, 95(8): 781-808. DOI:10.1016/j.jweia.2007.01.008
[27]
SUN Z K, WANG T T, WU F. Numerical investigation of influence of pantograph parameters and train length on aerodynamic drag of high-speed train[J]. Journal of Central South University, 2020, 27(4): 1334-1350. DOI:10.1007/s11771-020-4370-6
[28]
WANG T T, LEE C H, YANG M Z. Influence of enlarged section parameters on pressure transients of high-speed train passing through a tunnel[J]. Journal of Central South University, 2018, 25(11): 2831-2840. DOI:10.1007/s11771-018-3956-8
[29]
牛纪强. 高速列车通过隧道时产生的列车风研究[J]. 铁道科学与工程学报, 2015, 12(6): 1268-1276.
NIU J Q. Research on gusts caused by high-speed trains passing through tunnel[J]. Journal of Railway Science and Engineering, 2015, 12(6): 1268-1276. DOI:10.3969/j.issn.1672-7029.2015.06.002 (in Chinese)
[30]
国家铁路局. 中华人民共和国铁路运输行业标准: 高速铁路设计规范 TB 10621—2014[S]. 北京: 中国铁道出版社, 2015.
National Railway Administration. Railway & Train Standard of the People's Republic of China: code for design of high speed railway. TB 10621—2014[S]. Beijing: China Railway Publishing House, 2015. (in Chinese)
[31]
WANG J Y, WANG T T, YANG M Z, et al. Effect of localized high temperature on the aerodynamic performance of a high-speed train passing through a tunnel[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2021, 208: 104444. DOI:10.1016/j.jweia.2020.104444
[32]
LU Y B, WANG T T, YANG M Z, et al. The influence of reduced cross-section on pressure transients from high-speed trains intersecting in a tunnel[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2020, 201: 104161. DOI:10.1016/j.jweia.2020.104161
[33]
田四明, 王伟, 巩江峰. 中国铁路隧道发展与展望(含截至2020年底中国铁路隧道统计数据)[J]. 隧道建设(中英文), 2021, 41(2): 308-325.
TIAN S M, WANG W, GONG J F. Development and prospect of railway tunnels in China (including statistics of railway tunnels in China by the end of 2020)[J]. Tunnel Construction, 2021, 41(2): 308-325. (in Chinese)
[34]
赵勇, 田四明, 孙毅. 中国高速铁路隧道的发展及规划[J]. 隧道建设, 2017, 37(1): 11-17.
ZHAO Y, TIAN S M, SUN Y. Development and planning of high-speed railway tunnels in China[J]. Tunnel Construction, 2017, 37(1): 11-17. DOI:10.3973/j.issn.1672-741X.2017.01.002 (in Chinese)
[35]
LIU T H, TIAN H Q, LIANG X F. Design and optimization of tunnel hoods[J]. Tunnelling and Underground Space Technology, 2010, 25(3): 212-219. DOI:10.1016/j.tust.2009.12.001
[36]
YAMAMOTO A. Pressure variations, aerodynamic drag of train, and natural ventilation in Shinkansen type tunnel[J]. Quarterly Report of Railway Technical Research Institute (RTRI), 1974, 15(4): 207-214.