地铁活塞风相关计算

第一章 活塞风的理论基础及风速计算 1.1 活塞风的基本概念 当列车在隧道中运行时，隧道中的空气被列车带动而顺着列车运行前进的方向流动，这一现象称为列车的活塞作用，所形成的气流称为活塞气流。 列车在空旷的地面上运行时，列车前面的空气可毫无阻挡地被排挤到列车的两侧和上方，然后绕流到列车的后面。列车在隧道中运行时，由于隧道壁所构成空间的限制，列车所推挤的空气不能全部绕流到列车后方，必然有部分空气会被列车向前推动，排出到隧道出口之外；而列车尾端后方存在着负压涡旋区域，因此也必然会有相应空气经开口被引入到隧道中，由此形成活塞风。如图2－1。地铁活塞风的大小与列车在隧道内的阻塞比、列车的行驶速度、列车行驶时的空气阻力、空气与隧道壁面间的摩擦力等因素有关。

隧道壁隧道壁活塞风活塞风列车

图2-1 活塞风成因示意图 1.2 活塞风模型的简化 由于地铁隧道中活塞风的影响因素较多且活塞风速的计算复杂，在对计算结果误差影响较小的情况下，本文的计算中对活塞风的简化如下： （1）根据流体力学的基本原理，当气流速度小于音速时，流体密度的变化很小，流体的压缩性可以忽略不计（在标准状况下，如果气流速度不超过60m/s，则不考虑压缩性所引起的相对误差不大于1％[37]）。地铁车辆最大行驶速度一般不超过35m/s（126公里/每小时），产生的活塞风速远小于音速，因此在本论文中，如无特殊说明，所进行分析的地铁隧道活塞风气流均认为是不可压缩流体。 （2）根据管内流动的基本性质，当流体的雷诺数Re﹤2000时，管内流动称为层流，粘性力起主要作用，空气横断面上的流速梯度明显。而当流体的雷诺数Re﹥2000时，管内流动逐渐转化为紊流。在靠近壁面的一个薄层内，流动仍保持层流状态，称为层流底层。层流底层之外则是紊流区。地铁活塞风的雷诺数Re﹥2000，属于紊流流动。而在紊流状态下，整个空气横断面上的速度分布比较均匀（如图2-2所示），测试时可以用壁面附近的测点速度值代表隧道断面的活塞风速。因此本论文中，将地铁活塞风近似看成是沿隧道方向的一维流动。

管内层流管内紊流 图2-2 管内流动示意图 （3）列车进入入隧道后的一段时间内,活塞风压的压源是随列车而移动的,活塞风速随时间增大。但当隧道足够长时,一段时间后活塞风速便趋于一稳定值,活塞风基本达到稳定流状态,活塞风压稳定不变,与列车走行位置无关[38~39]相对而言,地铁隧道长度远大于列车长度,故在本论文中,地铁隧道活塞风可按恒定流计算。 （4）为简化计算模型,本文按一个区间内仅有一列车行驶考虑,且只考虑计算区段前后各两座活塞风井的作用,忽略相邻其他(前端及后端) 区段及列车的影响。

1.3 活塞风空气动力学基本理论方程 空气的流动要受到物理守恒定律的支配，其理论基础是空气动力学原理，即空气流动过程中的质量守恒、能量守恒和动量守恒定律。空气流动过程中的质量守恒、能量守恒和动量传递的定律是隧道通风的理论基础，流体在隧道中的运动应遵循空气动力学的基本方程。

1.3.1 气体流动的质量守恒方程（连续性方程） 任何流动问题都必须满足质量守恒定律。该定律可表述为：单位时间内流体微元体中质量的增加，等于同一时间间隔内流入该微元体的静质量。按照这一定律，可以得到质量守恒方程（连续性方程）： 0AAuAvAwtxyz (2-1)

上式中的第2，3，4项是质量流密度（单位时间通过单位面积的流体质量）的散度，可用矢量符号写出来： div()0AAUt (2-2)

对于隧道中的气体流动连续性方程常按一维运动连续性微分方程考虑： div()0AAvt (2-3)

式中 ——空气密度；

A——隧道横断面积，m2； v——气流平均速度，m / s；

1.3.2 气体运动的微分方程及伯努力方程（不可压缩流体的能量方程）

1、恒定流的伯努利 方程在隧道通风工程技术中，气流的密度变化可以忽略，可认为ρ为常量。 22'1122

12

aa22fppgzgzh (2-4)

υ——气流平均速度，m / s；

P——气流静压强，Pa； ρ——气流密度，3/mKg；

g——重力加速度，m / s2； z——位置水头，m； a——过流断面风速不均匀所引起的动能修正系数； 式（2-4)称为恒定流总流的伯努利方程。式中各项为单位质量流体的平均能量，各项的单位均为22/ms。据隧道通风模型试验和现场测试表明，在隧道横断面上各测点气流的流速分布比较均匀，测试断面各测点的流速与过流断面平均流速的比值比较接近，此时≈1。 2、非恒定流的伯努利 当列车驶入隧道时，在列车的推动下，隧道内的空气发生流动。严格来说流动速度是随着时间而发生变化的，这种流动称为非恒定流动。 dsthgzpgzpsf22221121

22 (2-5)

式（2-5)称为非恒定流总流的伯努利方程。式中 dsts称为惯性水头，它表示单位质量流体的动能随时间的变化量，其中β为流速不均匀的修正系数。对隧道中气流的≈1、β≈1，如果隧道的横断面积不变，即A1＝A2，则1＝2，

dtdt，并且sddsltdt，式中l为沿流线由断面1-1至断面2-2的距离。

在这种情况下，式（2-5)变为

dtdulhgzpgzpf2211

 （2-6）

在非恒定流计算列车通过隧道所引起的活塞风时，要应用式（2-5）。

1.4 鞍山道站活塞风速的计算及影响因素分析 1.4.1 鞍山道站列车出站活塞风速 鞍山道站相邻隧道均为设中柱的双跨矩形隧道，另外由于鞍山道站基本位于起点站和终点站之间，通常上下行列车同时到达鞍山道站，因此在鞍山道站相邻隧道不存在上下行列车相会的情况。所以，鞍山道站列车出站活塞风速可以按照单线有竖井隧道的情况计算列车未经过竖井之前的风速[40]。列车在隧道中运行时，竖井也成为气流通道，于是竖井两边的隧道段中的活塞风速不同。 设列车的速度为0，列车的横断面积为A0，隧道横断面积为A，海拔高度

为z，活塞风速为z，列车与隧道壁之间的环状空间中气流的绝对速度（即相对于隧道壁的速度）为0，在dt时间内，列车在隧道中移动所排开的空气体积为

00Adt；而在列车前方，则有部分空气推至列车前方隧道，其体积为zAdt，另

一部分空气通过列车与隧道壁之间的环状空间由列车前方流向列车后方，其体积为0-AAdt，如图2-3所示[41]。 ννν列车

压力线 图2-3鞍山道列车出站活塞风速简图 鞍山道站相邻隧道全长为l＝1030m，竖井的位置距隧道进口为lz＝1000m，距隧道出口为ly＝30m。由列车基本数据知，标准车长为19m，头车略长为19.19m，车厢两端的贯通连接通道长0.52m，若列车为6节编组，则列车长度

019.1921940.525116.98lm，若列车为4节编组，则列车长度

l0=78.98m。 隧道的水力直径74445.09742AdPm。

当列车在竖井左边的隧道段中运行时，列车前方气流压力为正压，一部分气流由隧道出口排出。列车后方的压力为负压，从而新鲜空气由隧道进口吸入。设左边隧道段中气流流量为Qz，流速为υz ，右边隧道段中的气流流量为Qy ，流速为υy ，竖井向大气排出的气流流量为Qs ，流速为υs ，由连续性方程得

yszQQQ 即 AAAyssz （a） 式中 A——隧道的横断面积，28m2； As——竖井的横断面积，9m2；

为求得υs 与υy的关系，可先建立图2-3中的断面4-4与竖井出口断面5-5的能量方程

222

554445455222ssssss

s

plpkp

gzgzgzd



（b）

式中 45——隧道气流分流入竖井的局部阻力系数，取0.5； s——竖井中各局部阻力系数，取4.5； sd——竖井的水力直径，3m； s——竖井沿程阻力系数，取0.3；

sl——竖井的长度，20m；

而 4520.+0.3+4.5=7.03ssssslk05d 然后再写出断面4-4与隧道出口断面6-6的能量方程 222

6644

464661222yyyylkpppgzgzgzd





（c）

式中 46——气流流经竖井与隧道连接处的局部阻力系数，取0.03； 而 46301=0.03+0.31=2.805.09yylkd 当隧道外无自然风时，5656ppgzgz 于是由式（b）和式（c）得 syyskk （d） 以（d）式代入（a）式得 ysysysysyzAkkAAkkAA

或 1xzyyssbkAAk （m / s） （e）

式中 92.8011.287.1ysskAb120Ak 式（e）说明在有竖井的隧道中，竖井两边的隧道段中的活塞风速不相等。 为了求得活塞风速υz 与υy的值，可由图2-3中的断面1-1与断面2-2间的气流能量方程求得P2，再由断面3-3与断面6-6之间的气流能量方程求得P3，从而算得列车前端与列车尾端的压差