地铁活塞风作用下屏蔽门的气动特性分析

根据三维不可压缩Navier2Stokes方程和标准k2ε湍流方程,采用有限体积法模拟列车在区间隧道内运行过程,研究了屏蔽门的压力变化机理和压力时空变化特性,并分析了列车运行速度、阻塞比、活塞风井面积、风阀状态和区间隧道通风方案等诸多因素与屏蔽门压力之间的影响关系。

2009年　　6月

第30卷　第2期

郑州大学学报(工学版)

JournalofZhengzhouUniversity(EngineeringScience)

Jun1　2009Vol130　No12

文章编号:1671-6833(2009)02-0120-05

地铁活塞风作用下屏蔽门的气动特性分析

杨伟超,彭立敏,施成华,胡自林

1

1

1

2

(1中南大学土木建筑学院,湖南长沙410075;2广州市地下铁道设计研究院,广东广州510010)

ε湍流方程,采用有限体积法模拟列车在区间隧摘　要:根据三维不可压缩Navier2Stokes方程和标准k2道内运行过程,研究了屏蔽门的压力变化机理和压力时空变化特性,并分析了列车运行速度、阻塞比、活塞风井面积、结果表明:屏蔽主要影响因素;于降低屏蔽门压力..关键词:地铁;活塞风;屏蔽门;中图分类号:U450　引言

地铁屏蔽门(PlatformScreenDoor,PSD)作为一种安全设施,设置于站台边缘,将列车与站台候车室(厅)隔离开来,不仅能防止乘客有意或无意跌入轨道,也可以阻断区间隧道与站台候车区域之间的空气流通,降低区间隧道热负荷对车站候车环境的影响,节约车站环控系统的运营成本.自2003年广州地铁二号线引入屏蔽门系统以来,国内的上海、深圳、重庆等城市地铁线路相继采用.但在车站安装屏蔽门会出现一些新的问题,比如:引起区间内外空气的交换率下降,车厢内空气质量变差,CO2等有害气体浓度增大,列车的乘坐舒适性下降;另外由于屏蔽门对气流的阻隔作用,列车的活塞效应显著增强,列车的气压荷载也随之增加,加之通风系统的影响和地铁运行速度的不断提高,活塞风引起的交变压力荷载显著增大,设计时需要考虑活塞风压荷载作用下屏蔽门的结构

[1]

强度和使用寿命.但目前国内外对屏蔽门的研究主要集中在机械受力特性和环控系统方[2-3]面,对屏蔽门的压力特性研究相对较少,尚没有明确的活塞风压荷载标准,相关的设计也多根据工程经验估计,由于经验不同,各处理方法也存

[4]

在一定的差异.因此,为了能够对屏蔽门的结构强度和使用寿命等问题进行深入分析,有必要对地铁活塞风作用下屏蔽门的气动特性进行

收稿日期:2008-10-14;修订日期:2009-01-13

探讨.

笔者对列车由区间隧道进入车站过程的运行过程进行了模拟,分析了活塞风作用下屏蔽门的压力变化机理,并对影响屏蔽门压力变化的列车运行速度、阻塞比、活塞风井面积、风阀状态和区间隧道通风方案等诸多因素进行了分析.

为了论述方便,将文中使用的符号集中列出:

+-Vr为车速;压力波幅值ΔP=P+P,其中,P为正压力峰值,P为负压力峰值;β=Ar/At,其中,Ar,At分别为列车及隧道面积;Ad为活塞风

井面积.

+

-

1　计算模型

1.1　控制方程

根据以上假设,区间隧道内空气流动控制方程如下:

(1)+=J

9t9x9y9z

其中,

ρuρ

ρu

vU=ρ

2ρu+p-τxx

vu-τ,F=ρxy

,

ρw

ρ(e+V2/2ρwu-τxy

ρ(e+V2/2)+　　作者简介:杨伟超(1978-),男,河南许昌人,中南大学博士研究生,主要从事隧道通风及空气动力学研究.

根据三维不可压缩Navier2Stokes方程和标准k2ε湍流方程,采用有限体积法模拟列车在区间隧道内运行过程,研究了屏蔽门的压力变化机理和压力时空变化特性,并分析了列车运行速度、阻塞比、活塞风井面积、风阀状态和区间隧道通风方案等诸多因素与屏蔽门压力之间的影响关系。

　第2期杨伟超等　地铁活塞风作用下屏蔽门的气动特性分析　121

ρv

ρuv-τyx

v+p-τG=ρyy

2

ρw

ρuw-τzx

vw-τzy,H=ρ

2ρw+p-τzz

1.3　边界条件

,

ρwv-τyz

ρ(e+V2/2)+0

隧道外部的无穷远处大气按照黎曼不变量的

无反射边界定义.隧道侧壁及列车表面为静止墙边界,墙边界处无涡流,气体的法向速度un和法向压力梯度9p/9n为零.根据T.Cebeci粗糙管壁模型试验确定隧道壁面粗糙度.隧道及风井与外部大气的连接处按压力出口边界定义,计算区域内产生的压力波一部分透射出边界,一部分被反射.区间隧道通风口采用速度入口边界定义.列.1.4.考3种方案:方案1分析区间隧道阻塞比的影响,保持列车横断面不变,调整区间隧道断面,使区间隧道阻塞比β=0.35～0.65;方案2分析列车运行速度的影响,时速Vr=50～120km/h过站;方案3车站隧道风机正常运行,区间隧道风机送/排风,不同活塞风井风面积,前后相邻车站隧道风机送/排风.

[

5]

ρ(e+V2/2)+ρfx

fyJ=ρ

.

ρfz

ρ(ufx+vfy+wfz式中:ρ,V为空气密度和流速,u,v,w为V分量;e为单位内能;p为总压,fx,fy,fzτ量;M,N,O分别为pu-uxx-vxyxz-ττ-vw-uyy-wyz,pvzy[:湍流模型采用k(ρik)uεi(2)(μ=+G-ρeff+σ99xi9xi

2

(εutεεi)μ+=+c1G-c2(3)

9xi9xieffkk222

ε式中:Gk=2uj,j+(ui,j+uj,i);ut=ρcwk/;cu=

0.09;c1=1.44;c2=1.92;σk=1.0;σ.3;ujε=1

为流场中各点的速度分量;ρ为流体密度;p为压强;μ为动力黏性系数;μt为湍流黏性系数;k为

湍流动能;ε为湍流能量耗散率.

2　活塞风作用下屏蔽门压力荷载特性

2.1　活塞风作用下屏蔽门压力变化机理

1.2　计算模型及参数

活塞风作用下屏蔽门压力变化主要受列车的有压科特湍流和区间隧道内压力波两方面的影响.图2为列车即将进站时车体和站台屏蔽门附近的压力分布图.

忽略列车电弓、车体连接部位、转向架和铁轨等细部结构,列车及隧道壁平整.几何模型按某地铁线路实际条件设置,如图1所示.列车4动2拖编组,长度取120m,单车高3.8m,宽2.8m,5个车门,车体距轨面0.89m;屏蔽门尺寸为2.0m×2.2m,站台两端活动门尺寸为1.6m×2.2m;地铁站台段区间隧道为矩形,列车停车范围140m,两端各延伸11m,总长度164m,站台两端距屏蔽

2

门端门6.8m处设16m活塞风井.

图2　隧道内运行时车体表面压力分布

Fig.2　Airpressuredistributiononexternal

surfaceoftrain

列车在运行过程中,车头的压力较高,一般为正压;而车尾的压力较低,一般为负压,车身的压力则介于二者之间,并由车头和车尾缓慢下降,沿车体出现压力梯度.列车与隧道壁之间的空气在压力梯度、列车外表面切应力和隧道壁面切应力的联合作用下形成“有压科特Couette湍流”.在

图1　计算模型示意图

Fig.1　Schematicdiagramofcalculationmodel

列车的进站过程中,屏蔽门的压力不断升高,在车头经过的瞬时达到最大值,之后压力随着车头的

根据三维不可压缩Navier2Stokes方程和标准k2ε湍流方程,采用有限体积法模拟列车在区间隧道内运行过程,研究了屏蔽门的压力变化机理和压力时空变化特性,并分析了列车运行速度、阻塞比、活塞风井面积、风阀状态和区间隧道通风方案等诸多因素与屏蔽门压力之间的影响关系。

　122

郑州大学学报(工学版)2009年

-

离开而持续下降,在车尾经过的瞬时压力降到最低,达到常压以下,随着列车离去,屏蔽门的压力又缓慢回升至常压.一般地铁车站两端均设置有风井,列车过风井时会产生两组不同的压力[7]

波,如图3所示,压力波在车站区间隧道内传播,再次引起屏蔽门压力变化

.

驶的过程而持续降低,在车尾经过时降到最低P处(对应图4中的D点),期间伴随着区间隧道内压力波的交替小幅波动.随着列车离去,屏蔽门压力逐渐恢复至常压.2.3　屏蔽门气压的空间分布

+-图5为屏蔽门的正负压力峰值(P、P)和压力幅值ΔP沿站台的分布

.

图3　列车过风井生成的附加压力波

Fig.3　Pressurewaveduringtrainpassing流受列车头部挤压,波,转换成膨胀波,并以音速向进口端返回;当列车尾部经过站台进口端活塞风井时情况正好相反,形成向站台出口端传播的膨胀波,膨胀波到达站台出口端活塞风井处转换成压缩波,同样以音速向进口端返回.压缩波与膨胀波在活塞风井之间反复转换与反射,且同种类型的波叠加时,压力波幅值增加,不同类型的波叠加时,压力波幅值减小.2.2　屏蔽门压力的时程变化

图4为β=0.4,车速90km/h过站时屏蔽门轨道侧气压力和区间隧道内压力波传播的对照变化图

.

图5　不同位置屏蔽门气压分布

Fig.5　Airpressuredistributionalongthe

metroplatform

由图5可以看出,站台进口段屏蔽门的正压

力峰值较大,负压力峰值较小,而出口段正压力峰值较小,负压力峰值较大,而压力幅值除了靠近活塞风井处以外,其它位置压力变化不大.主要原因是车头空气在进口段受挤压,流速慢,因此进口段+

P较大,在出口段可沿活塞风井外溢,流速快,压力P较小;而车尾的情况则正好相反.

+

3　影响屏蔽门压力的诸因素分析

3.1　阻塞比对屏蔽门压力的影响

取列车速度Vr=70km/h,β=0.35～0.65,活塞风井面积取Ad/At=0.7,10、15和20门的压力幅值见表1.

表1　不同阻塞比条件下屏蔽门压力幅值

Tab11　Comparisonofblockageratio

kPa

#

###

β

0.

0.0.0.0.0.354045505560测点位置

100.6620.7420.7680.7910.8120.8790.1391.260

#

图4　屏蔽门压力与站台内压力波传播对照

Fig.4　TimehistoryofpressureonPSD

列车没有进入站台内时,屏蔽门基本保持恒定常压,当列车进入站台进口端活塞风井时(对应图4中的A点),隧道内空气被车头挤压,屏蔽门气压开始上升,随着车头向前推进而不断上升;

+

车头经过时,气压升到最高P(对应图4中的B

AvB150.6750.7570.7690.7920.8380.9020.1481.343

#

200.6730.7570.7660.7990.8340.9020.1491.340

由表1可以看出,当列车速度一定时,屏蔽门压力与隧道阻塞比呈线性变化,建立如下关系式P=Av+Bvβ,其中:Av近似看作列车在明线运行时车周气体的压力幅值,系数Bv用曲线拟合的方

点),屏蔽门气压随着车头离开而降低至常压(对

应图4中的C点);并随着车尾由远至近不断行

根据三维不可压缩Navier2Stokes方程和标准k2ε湍流方程,采用有限体积法模拟列车在区间隧道内运行过程,研究了屏蔽门的压力变化机理和压力时空变化特性,并分析了列车运行速度、阻塞比、活塞风井面积、风阀状态和区间隧道通风方案等诸多因素与屏蔽门压力之间的影响关系。

　第2期杨伟超等　地铁活塞风作用下屏蔽门的气动特性分析　123

法得到.Av、Bv拟合值见表1,单位kPa.

时速70km条件下屏蔽门压力变化公式应为:

β(4)P=0.139+1.3143.2　车速对屏蔽门压力的影响

由图6可以看出,增大活塞风井面积可有效降低屏蔽门的压力幅度,原因是随着活塞风井面积增大,区间隧道内空气沿着风井向隧道外溢出的速度加快,隧道内风压减小,列车的科特风效应减弱.当活塞风井断面积与隧道断面积之比在0～0.3之间时,随着活塞风井断面积的增大,测点处的降压效果较明显.当活塞风井断面积与隧道断面积之比超过0.5后,屏蔽门压力下降有减缓趋势.

当Ad/At=0时,塞风井,,因此,设置活塞风.为研究车速对屏蔽门压力的影响,取活塞风

井面积Ad/At=0.7,列车速度Vr=50～120km/h,β=0.45,各点压力幅值见表2.

表2　不同车速条件下屏蔽门压力幅值

Tab.2　ComparisonofTrainSpeed

kPa

0.4040.7660.9881.505.#

车速/(km h-1)

507080100120

0.4020.7680.9881.5032.139

#

0.4080.7690.9841.133

#

[8]

对屏蔽门压

由表2,屏蔽门压力,建立如下关系:P

2ρ=aVr.当Vr=0时,P=0,因此系数aβ+bβ应为

零;同样参数b通过拟合得出,阻塞比β=0.45时上式可写成:

2ρ(5)P=0.00016Vr考虑到实际运行时屏蔽门不可能完全密封,

实际压力值可能小于计算结果,可以认为列车速度在50～70km/h范围内运行时,在阻塞比β不超过0.45的条件下,屏蔽门的压力幅值在0.5kPa左右,当前的处理方法基本是合理的.但随着地铁运行速度的不断提高(广州部分线路已

力的影响,取活塞风井面积Ad/At=0.7,列车速度Vr=70km/h,β=0.45,计算如下2种工况:

工况①:车站隧道风机正常运行,活塞风井风阀开启,车站内区间隧道风机送风,前后相邻两个车站隧道风机排风,称隧道送风模式;

工况②:车站隧道风机正常运行,活塞风井风阀开启,车站内区间隧道风机排风,前后相邻两个车站隧道风机送风,称隧道排风模式.

测点压力见表3.

表3　不同的通风方式条件下屏蔽门的压力幅值

Tab.3　ComparisonofventilationmodekPa

通风方式压力

+P

测点位置

5

#

10

#

15

#

25

#

达90km/h),屏蔽门压力按车速的平方上升.在

列车面积不变的情况下,

新建地铁站台段的区间隧道部分的面积应适当扩大.

3.3　活塞风井面积对屏蔽门压力的影响

0.3970.7310.595

0.4800.7690.578

0.4620.7680.568

0.293-0.2110.7690.603-0.2040.8070.339-0.4650.804

无通风-P-0.234-0.288-0.276

ΔP

+P

隧道送风

为了研究活塞风井面积与屏蔽门压力之间的关系,取Ad/At=0.1～0.9,β=0.45,Vr=70km/h、80km/h、100km/h和120km/h,15门

#

-P-0.233-0.216-0.250

ΔP0.8280.7940.819+P

0.3860.838

0.3860.877

0.4460.862

-P-0.452-0.491-0.416

隧道排风

压力为准,屏蔽门压力幅值与竖井面积关系曲线

绘于图6(70km/h).

ΔP

由表3可以看出,与工况②和关闭区间隧道通风系统条件相比,工况①屏蔽门的活塞风压明显降低.原因是区间隧道风机与活塞风压出现叠加或抵消,在阻塞比、车速等条件一定的条件下,活塞风压变化是一定的,区间隧道送风时,隧道风机与活塞风压先叠加后抵消,致使屏蔽门正的峰值压力变大,负峰值压力变小;区间隧道排风时,

图6　竖井面积与屏蔽门压力关系

Fig.6　Relationshipbetweenpressureandshaftarea

隧道风机与活塞风压先抵消后叠加,致使屏蔽门正的峰值压力变小,负峰值压力变大.但总的压力

根据三维不可压缩Navier2Stokes方程和标准k2ε湍流方程,采用有限体积法模拟列车在区间隧道内运行过程,研究了屏蔽门的压力变化机理和压力时空变化特性,并分析了列车运行速度、阻塞比、活塞风井面积、风阀状态和区间隧道通风方案等诸多因素与屏蔽门压力之间的影响关系。

　124

郑州大学学报(工学版)2009年

幅值变化不大,但列车进站时车速是逐渐降低的,屏蔽门的压力变化主要表现为负压较大,因此结构强度和刚度的校核时应考虑区间隧道排风模式的影响.

(5)列车进站时车速逐渐减小,屏蔽门主要

受负压影响,区间隧道排风模式对屏蔽门压力属不利工况,结构强度和刚度的校核时应考虑区间隧道排风模式的影响.

4　结论

采用三维不可压缩非定常等熵N2S方程和标

ε湍流方程,分析了活塞风作用下屏蔽门的准k2

气动压力特性,并得到如下结论:

(1)屏蔽门压力变化主要由列车的有压科特湍流和区间隧道内压力波两方面因素引起,其中有压科特湍流是主要影响因素;屏蔽门的压力变化规律主要表现为先上升,在车头经过时达到最高,然后下降,在车尾经过时达到最低;

(2)沿着列车运行方向减小,,的屏蔽门外,;

(3)列车的压力变化幅度与阻塞比和列车速度的平方呈正比,70km/h以下的地铁运行速度,阻塞比β不超过0.45时,屏蔽门的压力幅值在0.5kPa左右,考虑到以后地铁运行速度的不断

参考文献:

[1]　中南大学土木建筑学院.广州市轨道交通三号线北

延段隧道内压力变化研究[R].长沙:中南大学土木建筑学院,2006.

[2]　陈海辉,胡小芳.[J].

华南理工大学然科学版,2005,33

(1)37-[3].[J].华

学报:自然科学版,2004,32(4):74

-77.

[4]　冯　炼,刘应清.地铁阻塞通风的数值模拟[J].中

国铁道科学,2002,23(3):120-123.

[5]　魏新利,张海红,王定标.旋风分离器流场的数值计

算方法研究[J].郑州大学学报:工学版,2005,26

(1):57-60.

[6]　CEBECIT,BRADSHAWP.Momentumtransferin

boundarylayers[M].NewYork:HemispherePublish2ingCorporation,1977.

[7]　RAGHUNATHANRS,KIMHD,SETOGUCHIT.

Aerodynamicsofhigh-speedrailwaytrain[J].Pro2gressinAerospaceSciences,2002,38:469-514.

提高,新建地铁的区间隧道部分面积应适当扩大;

(4)若站台区间有活塞风井时,打开站台两端的活塞风阀有利于降低屏蔽门压力,活塞风井的面积应控制在区间断面积的0.5倍以上;

AnalysisofAerodynamicCharacteristicsofPSDunder

thePistonWindInfluenceinSubway

YANGWei-chao,PENGLi-min,SHICheng-hua,HUZi-lin

searchInstitute,Guangzhou510010,China)

1

1

1

2

(1.DepartmentofCivilEngineering,CentralSouthUniversity,Changsha410075,China;2.GuangzhouMetroDesignandRe2

Abstract:Basedontheunsteadyincompressiblethree-dimensionalN-Sequationsandstandardk-

εturbulentmodel,thetravelprocesswassimulatedwiththefinitevolumemethod.Thepressurevariationmechanismanditsspatial-timechangecharacteristicofPSDwereanalyzed.Meanwhile,thispaperdiscussestherelationshipbetweenthepressureofPSDandsuchinfluencefactorswasastrainspeed,blockageratio,theareaandswitchconditionofpistonvent,aswellastheventilationschemeofthesubwaytunnel.Resultshowed:thepressuregradientalongwiththetrainandthepressurewavebetweenthepistonvents,comingintobeingduringthetraintravelinginsubwaytunnel,arethebasicfactorscausingthepressurefluctuationofPSD,especiallytheformer,whichiscalledtheplaneturbulentCouetteflow.ThescopeofPSDpressurefluctuationisproportionaltothetunnelblockageratioandthesquareoftrainspeed.Greaterareaofpistonventandswitchtheventvalve“ON”arefavorabletoreducethepressure.Itisnecessarytoconsiderthepressurevaria2tionininfluenceoftheexhaustventilationofthesubwaytunnel.

Keywords:subway;pistonwind;platformscreendoor(PSD);aerodynamiccharacteristic

本文来源：https://www.bwwdw.com/article/82t1.html