地铁活塞风相关计算.doc

第一章 活塞风的理论基础及风速计算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 气体流动的质量守恒方程（连续性方程）任何流动问题都必须满足质量守恒定律该定律可表述为：单位时间内流体微元体中质量的增加，等于同一时间间隔内流入该微元体的静质量按照这一定律，可以得到质量守恒方程（连续性方程）： (2-1) 上式中的第2，3，4项是质量流密度（单位时间通过单位面积的流体质量）的散度，可用矢量符号写出来： (2-2) 对于隧道中的气体流动连续性方程常按一维运动连续性微分方程考虑： (2-3) 式中 ——空气密度； A——隧道横断面积，m2；——气流平均速度，m / s；1.3.2 气体运动的微分方程及伯努力方程（不可压缩流体的能量方程）1、 恒定流的伯努利方程在隧道通风工程技术中，气流的密度变化可以忽略，可认为ρ为常量。

(2-4) υ——气流平均速度，m / s；P——气流静压强，Pa；ρ——气流密度，； g——重力加速度，m / s2；z——位置水头，m；a——过流断面风速不均匀所引起的动能修正系数；式（2-4)称为恒定流总流的伯努利方程式中各项为单位质量流体的平均能量，各项的单位均为据隧道通风模型试验和现场测试表明，在隧道横断面上各测点气流的流速分布比较均匀，测试断面各测点的流速与过流断面平均流速的比值比较接近，此时≈12、 非恒定流的伯努利当列车驶入隧道时，在列车的推动下，隧道内的空气发生流动严格来说流动速度是随着时间而发生变化的，这种流动称为非恒定流动 (2-5)式（2-5)称为非恒定流总流的伯努利方程式中 称为惯性水头，它表示单位质量流体的动能随时间的变化量，其中β为流速不均匀的修正系数对隧道中气流的≈1、β≈1，如果隧道的横断面积不变，即A1＝A2，则1＝2，，并且，式中l为沿流线由断面1-1至断面2-2的距离在这种情况下，式（2-5)变为 （2-6）在非恒定流计算列车通过隧道所引起的活塞风时，要应用式（2-5）。

1.4 鞍山道站活塞风速的计算及影响因素分析1.4.1 鞍山道站列车出站活塞风速鞍山道站相邻隧道均为设中柱的双跨矩形隧道，另外由于鞍山道站基本位于起点站和终点站之间，通常上下行列车同时到达鞍山道站，因此在鞍山道站相邻隧道不存在上下行列车相会的情况所以，鞍山道站列车出站活塞风速可以按照单线有竖井隧道的情况计算列车未经过竖井之前的风速[40]列车在隧道中运行时，竖井也成为气流通道，于是竖井两边的隧道段中的活塞风速不同设列车的速度为，列车的横断面积为A0，隧道横断面积为A，海拔高度为z，活塞风速为，列车与隧道壁之间的环状空间中气流的绝对速度（即相对于隧道壁的速度）为，在dt时间内，列车在隧道中移动所排开的空气体积为；而在列车前方，则有部分空气推至列车前方隧道，其体积为，另一部分空气通过列车与隧道壁之间的环状空间由列车前方流向列车后方，其体积为，如图2-3所示[41]图2-3鞍山道列车出站活塞风速简图鞍山道站相邻隧道全长为l＝1030m，竖井的位置距隧道进口为lz＝1000m，距隧道出口为ly＝30m由列车基本数据知，标准车长为19m，头车略长为19.19m，车厢两端的贯通连接通道长0.52m，若列车为6节编组，则列车长度m，若列车为4节编组，则列车长度l0=78.98m。

隧道的水力直径m当列车在竖井左边的隧道段中运行时，列车前方气流压力为正压，一部分气流由隧道出口排出列车后方的压力为负压，从而新鲜空气由隧道进口吸入设左边隧道段中气流流量为Qz，流速为υz ，右边隧道段中的气流流量为Qy ，流速为υy ，竖井向大气排出的气流流量为Qs ，流速为υs ，由连续性方程得即 （a）式中 A——隧道的横断面积，28m2；As——竖井的横断面积，9m2；为求得υs 与υy的关系，可先建立图2-3中的断面4-4与竖井出口断面5-5的能量方程 （b）式中 ——隧道气流分流入竖井的局部阻力系数，取0.5；——竖井中各局部阻力系数，取4.5；——竖井的水力直径，3m；——竖井沿程阻力系数，取0.3；——竖井的长度，20m；而 然后再写出断面4-4与隧道出口断面6-6的能量方程 （c）式中 ——气流流经竖井与隧道连接处的局部阻力系数，取0.03；而 当隧道外无自然风时，于是由式（b）和式（c）得 （d）以（d）式代入（a）式得或 （m / s） （e）式中 式（e）说明在有竖井的隧道中，竖井两边的隧道段中的活塞风速不相等。

为了求得活塞风速υz 与υy的值，可由图2-3中的断面1-1与断面2-2间的气流能量方程求得P2，再由断面3-3与断面6-6之间的气流能量方程求得P3，从而算得列车前端与列车尾端的压差以（e）式代入上式得 (f)式中 ——除环状空间外竖井左侧隧道段的阻力系数； ——竖井右侧隧道段的阻力系数； ——有竖井的折算阻力系数； 根据气流流动的连续性方程可写出于是 （m / s） （2-7）环状空间中的气流相对于列车的速度为 （m / s） （a’）式中 a——列车对隧道的阻塞比，现写出图2-3中的断面3-3与断面2-2之间气流相对于列车运动的伯努利方程：因此 （b’）式中 ——环状空间气流的沿程阻力系数； ——气流由列车前方的隧道段进入环状空间的进口局部阻力系数；——气流由环状空间进入列车后方隧道段的出口局部阻力系数；——列车的长度即环状空间的长度；d0——环状空间的水力直径考虑到我国单线隧道的横断面A和列车横断面A0都变化不大，环状空间的水力直径d0可认为是定值。

为计算方便起见，、可并到值中，于是可令 （c’）由现场试验得出[42]，当列车长度l0为150m至600m时，N=86×10-4m-1以（a’）式及（c’）式代入（b’）式得列车前方和后方的压力差 （g）式中 K——活塞作用系数[110] 由式(f)和式（g）得 因此,竖井左侧隧道段的活塞风速为 （m / s） （2-8） 因此当列车为6节编组时鞍山道站列车出站活塞风速为3.37m/s当列车为4节编组时鞍山道站列车出站活塞风速为3.31m/s1.4.2 鞍山道站列车进站活塞风速图2-4列车经过竖井后的（进站）活塞风简图应当按照列车经过竖井后在隧道中运行的情况计算列车出站引起的活塞风速列车经过竖井后在右边隧道段中运行，此时隧道中的气流情况发生变化，列车前方的气流由隧道出口排出，但列车后方的负压使得新鲜空气不仅经隧道进口吸入隧道内，而且也从竖井向隧道内吸入如图2-11所示[40]所以竖井两边隧道段的活塞风速和风量也相应的发生变化：左边隧道段的风速由变为，右边隧道段的风速由变为，竖井中的风速由变为并与方向相反。

现分别计算和如下：由连续性方程 （h）由图2-11中的断面5-5与断面4-4之间气流的伯努利方程 （i）由断面1-1与断面4-4之间气流的伯努利方程 （j）式中各符号的意义同2.2.3所述各式但 式中 ——竖井气流汇集流入隧道的局部阻力系数，取0.5； ——隧道气流与竖井气流汇流的局部阻力系数，取0.15当隧道外无自然风时，由式（h）、（i）及式（j）解得 （m / s） （k）式中 列车前端3-3断面与列车尾端2-2断面的压差以（k）式代入上式得。