地铁活塞风相关计算

第一章 活塞风的理论基础及风速计算1.1 活塞风的基本概念当列车在隧道中运行时，隧道中的空气被列车带动而顺着列车运行前进的方向流动，这一现象称为列车的活塞作用，所形成的气流称为活塞气流。

列车在空旷的地面上运行时，列车前面的空气可毫无阻挡地被排挤到列车的两侧和上方，然后绕流到列车的后面。

列车在隧道中运行时，由于隧道壁所构成空间的限制，列车所推挤的空气不能全部绕流到列车后方，必然有部分空气会被列车向前推动，排出到隧道出口之外；而列车尾端后方存在着负压涡旋区域，因此也必然会有相应空气经开口被引入到隧道中，由此形成活塞风。

如图2－1。

地铁活塞风的大小与列车在隧道内的阻塞比、列车的行驶速度、列车行驶时的空气阻力、空气与隧道壁面间的摩擦力等因素有关。

隧道壁隧道壁图2-1 活塞风成因示意图1.2 活塞风模型的简化由于地铁隧道中活塞风的影响因素较多且活塞风速的计算复杂，在对计算结果误差影响较小的情况下，本文的计算中对活塞风的简化如下：（1）根据流体力学的基本原理，当气流速度小于音速时，流体密度的变化很小，流体的压缩性可以忽略不计（在标准状况下，如果气流速度不超过60m/s ，则不考虑压缩性所引起的相对误差不大于1％[37]）。

地铁车辆最大行驶速度一般不超过35m/s （126公里/每小时），产生的活塞风速远小于音速，因此在本论文中，如无特殊说明，所进行分析的地铁隧道活塞风气流均认为是不可压缩流体。

（2）根据管内流动的基本性质，当流体的雷诺数Re﹤2000时，管内流动称为层流，粘性力起主要作用，空气横断面上的流速梯度明显。

而当流体的雷诺数Re﹥2000时，管内流动逐渐转化为紊流。

在靠近壁面的一个薄层内，流动仍活塞风基本达到稳定流状态,活塞风压稳定不变,与列车走行位置无关[38~39]相对而言,地铁隧道长度远大于列车长度,故在本论文中,地铁隧道活塞风可按恒定流计算。

（4）为简化计算模型,本文按一个区间内仅有一列车行驶考虑,且只考虑计算区段前后各两座活塞风井的作用,忽略相邻其他(前端及后端) 区段及列车的影响。

1.3活塞风空气动力学基本理论方程空气的流动要受到物理守恒定律的支配，其理论基础是空气动力学原理，即空气流动过程中的质量守恒、能量守恒和动量守恒定律。

空气流动过程中的质量守恒、能量守恒和动量传递的定律是隧道通风的理论基础，流体在隧道中的运动应遵循空气动力学的基本方程。

1.3.1气体流动的质量守恒方程（连续性方程）任何流动问题都必须满足质量守恒定律。

该定律可表述为：单位时间内流体微元体中质量的增加，等于同一时间间隔内流入该微元体的静质量。

按照这一定律，可以得到质量守恒方程（连续性方程）：()()()()0A Au Av Aw t x y zρρρρ∂∂∂∂+++=∂∂∂∂ (2-1) 上式中的第2，3，4项是质量流密度（单位时间通过单位面积的流体质量）的散度，可用矢量符号写出来：()div()0A AU tρρ∂+=∂ (2-2) 对于隧道中的气体流动连续性方程常按一维运动连续性微分方程考虑：()d i v ()0A Av tρρ∂+=∂ (2-3) 式中 ρ——空气密度；A ——隧道横断面积，m 2；v ——气流平均速度，m / s ；1.3.2 气体运动的微分方程及伯努力方程（不可压缩流体的能量方程）1、恒定流的伯努利方程在隧道通风工程技术中，气流的密度变化可以忽略，可认为ρ为常量。

22'112212a a 22f p pgz gz h υυρρ++=+++ (2-4) υ——气流平均速度，m / s ； P ——气流静压强，Pa ；ρ——气流密度，3/m Kg ；g ——重力加速度，m / s 2； z ——位置水头，m ；a ——过流断面风速不均匀所引起的动能修正系数；式（2-4)称为恒定流总流的伯努利方程。

式中各项为单位质量流体的平均能量，各项的单位均为22/m s 。

据隧道通风模型试验和现场测试表明，在隧道横断面上各测点气流的流速分布比较均匀，测试断面各测点的流速与过流断面平均流速的比值比较接近，此时α≈1。

2、非恒定流的伯努利当列车驶入隧道时，在列车的推动下，隧道内的空气发生流动。

严格来说流动速度是随着时间而发生变化的，这种流动称为非恒定流动。

ds th gz p gz p s f ⎰∂∂++++=++υβρυαρυα2222112122(2-5) 式（2-5)称为非恒定流总流的伯努利方程。

式中 ds ts ⎰∂∂υβ称为惯性水头，它表示单位质量流体的动能随时间的变化量，其中β为流速不均匀的修正系数。

对隧道中气流的α≈1、β≈1，如果隧道的横断面积不变，即A 1＝A 2，则υ1＝υ2，dt d t υυ=∂∂，并且s d ds lt dtυυ∂=∂⎰，式中l 为沿流线由断面1-1至断面2-2的距离。

在这种情况下，式（2-5)变为dtdulh gz p gz p f +++=+2211ρρ（2-6） 在非恒定流计算列车通过隧道所引起的活塞风时，要应用式（2-5）。

1.4 鞍山道站活塞风速的计算及影响因素分析 1.4.1 鞍山道站列车出站活塞风速鞍山道站相邻隧道均为设中柱的双跨矩形隧道，另外由于鞍山道站基本位于起点站和终点站之间，通常上下行列车同时到达鞍山道站，因此在鞍山道站相邻隧道不存在上下行列车相会的情况。

所以，鞍山道站列车出站活塞风速可以按照单线有竖井隧道的情况计算列车未经过竖井之前的风速[40]。

列车在隧道中运行时，竖井也成为气流通道，于是竖井两边的隧道段中的活塞风速不同。

设列车的速度为0υ，列车的横断面积为A 0，隧道横断面积为A ，海拔高度为z ，活塞风速为z υ，列车与隧道壁之间的环状空间中气流的绝对速度（即相对于隧道壁的速度）为0ω，在dt 时间内，列车在隧道中移动所排开的空气体积为00A dt υ；而在列车前方，则有部分空气推至列车前方隧道，其体积为z A dt υ，另一部分空气通过列车与隧道壁之间的环状空间由列车前方流向列车后方，其体积为()0-A A dt ω，如图2-3所示[41]。

压力线图2-3鞍山道列车出站活塞风速简图鞍山道站相邻隧道全长为l ＝1030m ，竖井的位置距隧道进口为l z ＝1000m ，距隧道出口为l y ＝30m 。

由列车基本数据知，标准车长为19m ，头车略长为19.19m ，车厢两端的贯通连接通道长0.52m ，若列车为6节编组，则列车长度019.1921940.525116.98l =⨯+⨯+⨯=m ，若列车为4节编组，则列车长度l 0=78.98m 。

隧道的水力直径()7444 5.09742A d P ⨯==⨯=+⨯m 。

当列车在竖井左边的隧道段中运行时，列车前方气流压力为正压，一部分气流由隧道出口排出。

列车后方的压力为负压，从而新鲜空气由隧道进口吸入。

设左边隧道段中气流流量为Q z ，流速为υz ，右边隧道段中的气流流量为Q y ，流速为υy ，竖井向大气排出的气流流量为Q s ，流速为υs ，由连续性方程得y s z Q Q Q +=即 A A A y s s z υυυ+= （a ） 式中 A ——隧道的横断面积，28m 2；A s ——竖井的横断面积，9m 2；为求得υs 与υy 的关系，可先建立图2-3中的断面4-4与竖井出口断面5-5的能量方程222554445455222s s s s s s s p l p k p gz gz gz d υυυξλξρρρ⎛⎫++=++++=++ ⎪⎝⎭（b ） 式中 45ξ——隧道气流分流入竖井的局部阻力系数，取0.5；s ξ——竖井中各局部阻力系数，取4.5；s d ——竖井的水力直径，3m ；s λ——竖井沿程阻力系数，取0.3；s l ——竖井的长度，20m ；而 4520.+0.3+4.5=7.03s s ss s l k 05d ξλξ=++=⨯ 然后再写出断面4-4与隧道出口断面6-6的能量方程2226644464661222y y y y l k p p p gz gz gz d υυυξλρρρ⎛⎫++=++++=++⎪⎝⎭（c ） 式中 46ξ——气流流经竖井与隧道连接处的局部阻力系数，取0.03； 而 46301=0.03+0.31=2.805.09y y l k dξλ=++⨯+ 当隧道外无自然风时，5656p p gz gz ρρ+=+于是由式（b ）和式（c ）得 sy y s k k υυ= （d ）以（d ）式代入（a ）式得y s y s y s ysy z A k k A A k k A A υυυυ⎪⎪⎭⎫ ⎝⎛+=+= 或z y bυυ==（m / s ） （e ）式中11.b 120=+== 式（e ）说明在有竖井的隧道中，竖井两边的隧道段中的活塞风速不相等。

为了求得活塞风速υz 与υy 的值，可由图2-3中的断面1-1与断面2-2间的气流能量方程求得P 2，再由断面3-3与断面6-6之间的气流能量方程求得P 3，从而算得列车前端与列车尾端的压差220343246232204623222310.5()220.51()22()22y y z z y y z z yz zyl l l l p p g z z d d l l l g z z d d k k g z z ρυρυξλλρρυρυλξλρρυρυρ⎛⎫--⎛⎫-=+++++- ⎪ ⎪⎝⎭⎝⎭⎛⎫-⎛⎫=+++++- ⎪ ⎪⎝⎭⎝⎭=++- 以（e ）式代入上式得222232231()222z zzy tsp p k k g z z b ρυρυρυρξ⎛⎫-=++-= ⎪⎝⎭(f ) 式中 z k ——除环状空间外竖井左侧隧道段的阻力系数； 01000-116.980.50.5+0.3.5.09z z l l k 5254d λ-=+=⨯=y k ——竖井右侧隧道段的阻力系数；461y y l k dξλ=++ts ξ——有竖井的折算阻力系数；2211.+2.8=54.481.2ts z y k k 5254b ξ⎛⎫⎛⎫=+=⨯ ⎪ ⎪⎝⎭⎝⎭根据气流流动的连续性方程可写出000()z A dt A dt A A dt υυω=+-于是 00zA A A A υυω-=- （m / s ） （2-7）环状空间中的气流相对于列车的速度s υ为000000011z z zs A A A A A a A υυυυυυυωυυ---=+=+==--⎛⎫- ⎪⎝⎭（m / s ） （a’）式中 a ——列车对隧道的阻塞比，03 3.80.4174A a A ⨯===⨯。

现写出图2-3中的断面3-3与断面2-2之间气流相对于列车运动的伯努利方程：()()22202102022d 2z z ss p l p υυυυυξλξρρ--⎛⎫+=++++ ⎪⎝⎭因此 203210202s l p p d ρυξλξ⎛⎫-=++ ⎪⎝⎭ （b’） 式中 0λ——环状空间气流的沿程阻力系数；1ξ——气流由列车前方的隧道段进入环状空间的进口局部阻力系数；2ξ——气流由环状空间进入列车后方隧道段的出口局部阻力系数；0l ——列车的长度即环状空间的长度；d 0——环状空间的水力直径。