地铁活塞风相关计算

第一章活塞风的理论基础及风速计算1.1活塞风的基本概念当列车在隧道中运行时，隧道中的空气被列车带动而顺着列车运行前进的方向流动，这一现象称为列车的活塞作用，所形成的气流称为活塞气流。列车在空旷的地面上运行时，列车前面的空气可毫无阻挡地被排挤到列车的两侧和上方，然后绕流到列车的后面。列车在隧道中运行时，由于隧道壁所构成空间的限制，列车所推挤的空气不能全部绕流到列车后方，必然有部分空气会被列车向前推动，排出到隧道出口之外；而列车尾端后方存在着负压涡旋区域，因此也必然会有相应空气经开口被引入到隧道中，由此形成活塞风。如图2—1。地铁活塞风的大小与列车在隧道内的阻塞比、列车的行驶速度、列车行驶时的空气阻力、空气与隧道壁面间的摩擦力等因素有关。隧道壁活塞风——活塞风隧道壁活塞风——活塞风图2-1活塞风成因示意图1.2活塞风模型的简化由于地铁隧道中活塞风的影响因素较多且活塞风速的计算复杂，在对计算结果误差影响较小的情况下，本文的计算中对活塞风的简化如下：（1）根据流体力学的基本原理，当气流速度小于音速时，流体密度的变化很小，流体的压缩性可以忽略不计（在标准状况下，如果气流速度不超过60m/s，则不考虑压缩性所引起的相对误差不大于1%[37]）。地铁车辆最大行驶速度一般不超过35m/s（126公里/每小时），产生的活塞风速远小于音速，因此在本论文中，如无特殊说明，所进行分析的地铁隧道活塞风气流均认为是不可压缩流体。根据管内流动的基本性质，当流体的雷诺数Re<2000时，管内流动称为层流，粘性力起主要作用，空气横断面上的流速梯度明显。而当流体的雷诺数Re>2000时，管内流动逐渐转化为紊流。在靠近壁面的一个薄层内，流动仍保持层流状态，称为层流底层。层流底层之外则是紊流区。地铁活塞风的雷诺数Re>2000，属于紊流流动。而在紊流状态下，整个空气横断面上的速度分布比较均匀(如图2-2所示)，测试时可以用壁面附近的测点速度值代表隧道断面的活塞风速。因此本论文中，将地铁活塞风近似看成是沿隧道方向的一维流动。管内层流管内紊流图2-2管内流动示意图列车进入入隧道后的一段时间内，活塞风压的压源是随列车而移动的,活塞风速随时间增大。但当隧道足够长时，一段时间后活塞风速便趋于一稳定值,活塞风基本达到稳定流状态，活塞风压稳定不变，与列车走行位置无关［38〜39相对而言,地铁隧道长度远大于列车长度，故在本论文中,地铁隧道活塞风可按恒定流计算。为简化计算模型，本文按一个区间内仅有一列车行驶考虑,且只考虑计算区段前后各两座活塞风井的作用，忽略相邻其他(前端及后端)区段及列车的影响。1.3活塞风空气动力学基本理论方程空气的流动要受到物理守恒定律的支配，其理论基础是空气动力学原理，即空气流动过程中的质量守恒、能量守恒和动量守恒定律。空气流动过程中的质量守恒、能量守恒和动量传递的定律是隧道通风的理论基础，流体在隧道中的运动应遵循空气动力学的基本方程。1.3.1气体流动的质量守恒方程(连续性方程)任何流动问题都必须满足质量守恒定律。该定律可表述为：单位时间内流体微元体中质量的增加，等于同一时间间隔内流入该微元体的静质量。按照这一定律，可以得到质量守恒方程（连续性方程）:（2-1）8（pA）5（pAu）合（pAv）合（pAw）+++=0dtdxdydz（2-1）上式中的第2,3,4项是质量流密度（单位时间通过单位面积的流体质量）的散度，可用矢量符号写出来:8（pA）（2-2）一+div（pAU）=0所（2-2）对于隧道中的气体流动连续性方程常按一维运动连续性微分方程考虑:（2-3）8（pA）（2-3）+divpAv=）0dt式中p——空气密度；A——隧道横断面积，m2；v气流平均速度，m/s；1.3.2气体运动的微分方程及伯努力方程（不可压缩流体的能量方程）1、恒定流的伯努利方程在隧道通风工程技术中，气流的密度变化可以忽略，可认为p为常量。au2pau2p—2^+—i+gz=~2^+t+gz+h'（2-4）u气流平均速度，m/s；P一气流静压强，Pa；p——气流密度，Kg,m3；g重力加速度，m/s2；z位置水头，m；a——过流断面风速不均匀所引起的动能修正系数；式（2-4）称为恒定流总流的伯努利方程。式中各项为单位质量流体的平均能量，各项的单位均为m2/s2。据隧道通风模型试验和现场测试表明，在隧道横断面上各测点气流的流速分布比较均匀，测试断面各测点的流速与过流断面平均流速的比值比较接近，此时a^1。2、非恒定流的伯努利当列车驶入隧道时，在列车的推动下，隧道内的空气发生流动。严格来说流动速度是随着时间而发生变化的，这种流动称为非恒定流动。a^2+乙+gz=a^2+%+gz+h+gj竺（2-5）2p12p2fsdt式（2-5）称为非恒定流总流的伯努利方程。式中pf^^ds称为惯性水头，它表示sdt单位质量流体的动能随时间的变化量，其中乃为流速不均匀的修正系数。对隧道中气流的a21、&el,如果隧道的横断面积不变，即*=A2,则口1=口2，火=竺，并且fdlds=l史，式中l为沿流线由断面1-1至断面2-2的距离。dtdtsdtdt在这种情况下，式（2-5）变为P+gz=%+gz+h+1四（2-6）p1p2fdt在非恒定流计算列车通过隧道所引起的活塞风时，要应用式（2-5）。1.4鞍山道站活塞风速的计算及影响因素分析1.4.1鞍山道站列车出站活塞风速鞍山道站相邻隧道均为设中柱的双跨矩形隧道，另外由于鞍山道站基本位于起点站和终点站之间，通常上下行列车同时到达鞍山道站，因此在鞍山道站相邻隧道不存在上下行列车相会的情况。所以，鞍山道站列车出站活塞风速可以按照单线有竖井隧道的情况计算列车未经过竖井之前的风速［4叽列车在隧道中运行时，竖井也成为气流通道，于是竖井两边的隧道段中的活塞风速不同。设列车的速度为u°，列车的横断面积为A0,隧道横断面积为A，海拔高度为z，活塞风速为uz，列车与隧道壁之间的环状空间中气流的绝对速度（即相对于隧道壁的速度）为气，在dt时间内，列车在隧道中移动所排开的空气体积为A°u°dt；而在列车前方，则有部分空气推至列车前方隧道，其体积为Audt，另一部分空气通过列车与隧道壁之间的环状空间由列车前方流向列车后方，其体积为（A-A°）3dt，如图2-3所示［41］。图2-3鞍山道列车出站活塞风速简图鞍山道站相邻隧道全长为/=1030m,竖井的位置距隧道进口为L=1000m,距隧道出口为ly=30m。由列车基本数据知，标准车长为19m,头车略长为19.19m,车厢两端的贯通连接通道长0.52m,若列车为6节编组，则列车长度l=19.19x2+19x4+0.52x5=116.98m,若列车为4节编组，则列车长度0l0=78.98m。隧道的水力直径d=4-=4x『7x4=5.09m。P(7+4)x2当列车在竖井左边的隧道段中运行时，列车前方气流压力为正压，一部分气流由隧道出口排出。列车后方的压力为负压，从而新鲜空气由隧道进口吸入。设左边隧道段中气流流量为q,流速为％，右边隧道段中的气流流量为q,流速为％，竖井向大气排出的气流流量为'qs，流速为气，由连续性方程得Q广Q+Qy即uA=uA+uA(a)式中A——隧道的横断面积，28m2；As——竖井的横断面积，9m2；为求得气与％的关系，可先建立图2-3中的断面4-4与竖井出口断面5-5的能量方程''pu2p二l如pko2乌+gz+寸=f+gz+&+人丁+&言=f+gz+w

P42P5"45sdsJ2P52s式中&45——隧道气流分流入竖井的局部阻力系数，取0.5；&:一竖井中各局部阻力系数，取4.5；

d,——竖井的水力直径，3m；人:——竖井沿程阻力系数，取0.3；[—竖井的长度，20m；而k=&+人L+&=0.5+0.3x20+4.5=7.0s45sds3然后再写出断面4-4与隧道出口断面6-6的能量方程（C）pLn\=-6+g<+&46+x方+1ku（C）式中&46——气流流经竖井与隧道连接处的局部阻力系数，取0.03；而k=&+人L+1=0.03+0.3x^+1=2.80y46d5.09当隧道外无自然风时，§+g‘5=§+gZ6于是由式0）和式（于是由式0）和式（C）得=u'ky商s（d）以（d）式代入（1）式得UzA+U/tr-—\kUzA+U/tr-—\kAislkk式中b=1+%!s（m/s）（e）式（e）说明在有竖井的隧道中，竖井两边的隧道段中的活塞风速不相等。为了求得活塞风速％与七的值，可由图2-3中的断面1-1与断面2-2间的气流能量方程求得P2,再由断面3-3与断面6-6之间的气流能量方程求得P3,从而算得列车前端与列车尾端的压差

P3-P2PU2y2r.i\=kE46+方+1JL.i\+&46，+1kUJPU2+Pg(z_—z3)PU2土+Pg(zP3-P2PU2y2r.i\=kE46+方+1JL.i\+&46，+1kUJPU2+Pg(z_—z3)PU2土+Pg(z2—Z3)式中kz除环状空间外竖井左侧隧道段的阻力系数；l—l1000-116.98式中kz5.09k=0.5+人zdo=0.5+0.3x=52.54ky——竖井右侧隧道段的阻力系数；有竖井的折算阻力系数;=52.54+2.8x2=54.485.09=52.54+2.8x2=54.48根据气流流动的连续性方程可写出Audt=Audt+(A一A於dt于是(m于是(m/s)（2-7）环状空间中的气流相对于列车的速度七为u=w+uAu—Au=w+uAu—Au0z+uA—A00u—u1-k1J(m/s)（a’）..A式中a——列车对隧道的阻塞比，a=普=3x38=0.41。7x4现写出图2-3中的断面3-3与断面2-2之间气流相对于列车运动的伯努利方程：p(o-O)2O+—0z—P2p(o-o=―2+―0P2(.l，一、+&1+*0才+&2VdJ0O2―s-2TOCp(o-O)2O+—0z—P2p(o-o=―2+―0P2(.l，一、+&1+*0才+&2VdJ0O2―s-2'07式中人0——环状空间气流的沿程阻力系数；{——气流由列车前方的隧道段进入环状空间的进口局部阻力系数；&2——气流由环状空间进入列车后方隧道段的出口局部阻力系数；10一列车的长度即环状空间的长度；d0——环状空间的水力直径。考虑到我国单线隧道的横断面A和列车横断面A0都变化不大，环状空间的水力直径d0可认为是定值。为计算方便起见，§、&2可并到人0值中，于是可令（c’）由现场试验得出［42］，当列车长度10为150m至600m时，N=86x10-4m-1。以（。’）式及（c’）式代入（b^）式得列车前方和后方的压力差（g）N1P(o-o)2Kp(o-o（g）"EG=—L式中K——活塞作用系数m。］。K=8；x10j1086x10-4xK=8；x10j10(1-0.41》=2.89由式（/）和式（g）得Kp（O0-o»=§2O因此，竖井左侧隧道段的活塞风速为18O=3.37(m/s)（2-8）18O=3.37(m/s)（2-8）因此当列车为6节编组时鞍山道站列车出站活塞风速为3.37m/s。当列车为4节编组时鞍山道站列车出站活塞风速为3.31m/s。1.4.2鞍山道站列车进站活塞风速图2-4列车经过竖井后的（进站）活塞风简图图2-4列车经过竖井后的（进站）活塞风简图应当按照列车经过竖井后在隧道中运行的情况计算列车出站引起的活塞风速。列车经过竖井后在右边隧道段中运行，此时隧道中的气流情况发生变化，列车前方的气流由隧道出口排出，但列车后方的负压使得新鲜空气不仅经隧道进口吸入隧道内，而且也从竖井向隧道内吸入如图2-11所示［40。所以竖井两边隧道段的活塞风速和风量也相应的发生变化：左边隧道段的风速由u［变为u；，右边隧道段的风速由uy变为u；,竖井中的风速由u、变为u.并与u、方向相反。现分别计算u'和u'如下：y由连续性方程u'A=uA+u'A（h）由图2-11中的断面5-5与断面4-4之间气流的伯努利方程=P+gz+u2+上*

p42s2（i）pu2,(u')2

T+gz+艾+k—-z—（j=P+gz+u2+上*

p42s2（i）pu2,(u')2

T+gz+艾+k—-z—（j）p+gz=p+gz+二+(0.5+人L+&)也^2p1P42d142式中各符号的意义同2.2.3所述各式。但k'=&+人匕+&'=0.5+0.3x20+4.5=7.0S54Sds3

k'=0.5+人七+&=0.5+0.3X116^+0.15=7.54zd145.09式中七一一竖井气流汇集流入隧道的局部阻力系数，取0.5；§4一一隧道气流与竖井气流汇流的局部阻力系数，取0.15。当隧道外无自然风时，p+gq=p+gz5由式（h）、（i）及式（/）解得式中ub=1+式中ub=1+^As(m/s)（k）列车前端3-3断面与列车尾端2-2断面的压差/芬O(u')2l—lP(u')2,、P—P=(0.5+J+&)L+(V^-^+1)—+Pg(z—z)32d142d223以（k）式代入上式得式中k'=Xy1、,式中k'=Xy1、,'zb'yP(u')2,、,P(u')2—+Pg(z—z)=&'—223ts2（I）500-116.98509+1=23.57g'=k，(一)2+k'=7.54xtszb'y(1¥——+23.57=27.83由相对于列车的能量方程得（m）KP(u—u'（m）—+Pg(J7由（l）式及（m）式得右边隧道段的活塞风速u'=——y1+gs1u'=——y1+gs1+18,=4.38(m/s)27.832.