运行速度对高速列车火灾热释放速率的影响

2021-12-02 06:39刘方吉毕海权高慧翔曾清珣
制冷与空调 2021年5期
关键词:峰值数值速率

刘方吉 毕海权 高慧翔 曾清珣 王 菁

(1.西南交通大学机械工程学院 成都 610031;2.中铁二院工程集团有限责任公司 成都 610031)

0 引言

近年来,随着我国经济的持续快速发展,高速列车的安全运营越来越受到人们的重视。由于高速列车内部空间封闭、人员密集,一旦发生火灾,极易造成群死群伤事故,产生严重后果。因此,世界各国对列车防火安全都提出了很高的要求。高速列车着火后,应急救援原则是尽可能将高速列车行驶至安全地点再组织乘客疏散逃生。高速列车继续行驶时产生的空气流场必然会对着火车厢内火灾燃烧过程产生影响,一方面,空气流动可以为客室内的燃烧提供更多的新鲜空气;另一方面,高速行驶的列车产生的空气流场也可能会带走较多热量,增加火灾燃烧区域的热量损失,进而降低可燃物的燃烧速率。因此,列车着火后继续行驶,既有风险,又具有一定的合理性,其关键在于列车的运行速度对着火列车安全性的影响。

热释放速率是描述列车火灾过程最重要的参数,它体现了火灾中能量释放的多少,是决定列车火灾危险性的重要因素。火灾热释放速率是可燃物燃烧时单位时间内释放出的热量[1],大部分可燃物燃烧时,热释放速率随时间发生变化[2]。在火灾现场中,它能够表征火灾的大小,预测火焰蔓延速度[3]、热量及烟气产生量,提供可用于安全疏散的时间以及扑灭火灾可釆取的最有效方法等信息,同时,火灾热释放速率也是比较火灾行为以及评价火源对附近其他可燃物影响最可靠的数据。在火灾数值模拟中,火场的能见度[4]、温度场、氧气浓度、烟气层高度、火羽流的卷吸特性[5]等参数都是以火灾热释放速率为基础进行计算的。在消防工程设计中确定火灾热释放速率的变化规律可以定量评估分析建筑物的火灾危险性。在消防性能化设计[6]中的关键步骤是火灾场景的设定以及定量分析火灾发生过程,这些都是以热释放速率随时间变化为基础而展开的,热释放速率设定越接近实际情况,计算得到的结果越准确可靠。因此,在我国铁路事业高速发展的情况下,进行高速列车热释放速率[7]研究具有非常重要的现实意义。

本文在以某型二等座车厢为研究对象,采用美国国家标准与技术研究室发布的FDS 火灾动态模拟软件,对高速列车在不同时速30km/h~80km/h运行时的热释放速率大小进行研究。

1 数值计算模型

1.1 数值模型

如图1 和图2 所示分别为某型二等座车厢客室断面和平面布局。车厢尺寸为24.175m×3.257m×3.890m,座椅采用2+3 布置,定员85 人,两侧墙上设置6 个紧急逃生窗,如图3 车辆侧视图所示。

图1 客室断面图Fig.1 Section map of the passenger compartment

图2 客室平面布局Fig.2 Layout of the passenger compartment

图3 车辆侧视图Fig.3 Side view of the vehicle

高速列车内部构成极其复杂,拥有大量零部件。其可燃烧的部分主要包括木地板、座椅、侧墙板、顶板、间壁、电子线路等。本文主要建立了车厢内部空间的物理模型。最终,运用FDS 建立的高速列车仿真模型如图4 所示。

图4 列车客室火灾数值计算模型Fig.4 Numerical calculation model of the train compartment fire

1.2 材料参数

锥形量热仪(cone calorimeter)是20 世纪80年代早期美国国家标准与技术研究院开发的一种基于耗氧原理的实验室尺度材料燃烧行为测试仪器,如图5 所示。锥形量热仪具有精度较高、参数测定值受外界因素影响小、与大型实验结果相关性好等优点,被广泛应用于火灾科学领域的研究。锥形量热仪可以检测到燃烧样品的多个性能参数,如点燃时间(TTI)、热释放速率(HRR)、热释放速率峰值(PHRR)、热释放速率峰值对应时间(T)、总释热量(THR)、质量损失速率(MLR)、比消光面积(SEA)、有效燃烧热(EHC)以及相关燃烧气体的毒性和腐蚀性等[8]。因此,本文采用符合ISO5660 国际标准的锥形量热仪作为材料燃烧特性测试的标准实验方法。

图5 锥形量热仪结构及外形Fig.5 Cone calorimeter schematic structure and actual view

所用实验样品面积尺寸均为10cm×10cm,厚度为车体的真实厚度。实验过程中除了样品顶面露置在空气中外,其余各面均用铝箔包覆。实验前,将样件放置在23℃、相对湿度50%环境下调节48小时至质量恒定(质量变化小于0.1g)。如图6 所示,在实验进行时,对样品照射加热,使其燃烧。

图6 样品材料实验测试Fig.6 Experimental test of sample materials

高速列车客室内装材料比较复杂,同一部件不同区域材料组成也不尽相同,例如,侧墙采用的防寒材在不同部位可能选用不同材料。为简化计算模型,本文主要考虑座椅、侧墙板、顶板、地板及间壁等占比较大、对客室火灾蔓延有重要作用的可燃材料。锥形量热仪测试得到主要可燃材料单位面积热释放速率曲线如图7 所示。

图7 主要可燃材料单位面积热释放速率曲线Fig.7 Heat release rate curve of main combustible materials per unit area

1.3 网格尺寸与计算域

本文采用的是FDS 湍流模型中的大涡模拟,网格尺寸应小到足以包含大涡运动的湍流尺度。但是,受计算机资源限制,过小的网格尺寸会导致计算时间增长,过大的网格尺寸不能保证计算结果的精确性。在进行火灾数值模拟之前,需要对网格尺寸进行评估分析。

在火灾模拟中,需要精确的模拟出火灾场中的温度分布、蔓延情况及热释放速率变化等关键参数,因此,模拟的网格需要足够精细。FDS 软件用户指南给出了火灾特征直径的估计值D*:

式中,D*为火源特征直径,m;Q为总热释放速率,kW;ρ∞为环境空气密度,此处取1.204kg/m³;Cp为环境空气比热,此处取1.005,kJ/kg·K;T∞为环境空气温度,K;g为重力加速度,m/s2。

由式(1)可知,影响火灾特征直径D*的关键参数为火灾热释放速率Q。火灾特征直径D*随火灾热释放速率增大而增大。

许多学者对网格尺寸的选取进行了研究,如胡隆华[9]在隧道火灾研究中心将网格尺寸取在0.1~0.12D*之间。YZ Li[10]在利用FDS 对瑞典SP全尺寸地铁列车火灾实验进行数值计算研究中将网格尺寸取为0.1m,计算值与实验结果比较相符。参考以上研究,本文决定将网格尺寸设为0.08m。

高速列车燃烧过程中,车厢内热烟气与车厢周围新鲜空气相互交换,故计算区域也是影响列车火灾数值计算的一个重要因素。为保证计算结果的准确性,需将开口处的计算区域向外延展适当距离。王升[11]通过对不同计算区域扩展量时我国高速列车火灾热释放速率的计算表明,当列车外计算区域扩展至1.5m 后,继续扩展计算区域对火灾热释放速率的影响较小。本文参考其研究结果来确定高速列车客室火灾计算区域的设置。如图8 所示为本文高速列车火灾数值计算模型计算区域设置示意图。

图8 计算区域设置Fig.8 Calculation locale

利用FDS 中的supply 边界即进口速度边界条件来模拟高速列车运行时从列车周围扰流的空气速度场。高速列车运行过程中车门均关闭,为考虑空气流场与车厢内燃烧的相互作用,考虑实际情况中可能出现乘客击碎列车逃生窗的情形,将列车火源附近的逃生窗全部打开。

2 数值计算结果及分析

如图9 所示为计算得到的各运行速度下列车客室火灾热释放速率曲线。从图中可以看出,0~200s 时各运行速度下的热释放速率增长趋势几乎一致,200s 之后各运行速度下列车热释放速率曲线可以明显地分为两类:30km/h、40km/h、50km/h时火灾热释放速率曲线较为接近,60km/h、70km/h、80km/h 时的火灾热释放速率曲线也较为接近。以80km/h 时的火灾热释放速率曲线为例,火灾热释放速率曲线明显可以分为两个阶段:初期增长阶段和快速增长阶段。0~200s 即为初期增长阶段,火灾热释放速率在这段时间内增长速度较慢,热释放速率始终保持在较低水平,200s 时热释放速率仅仅为2.96MW。200s 后火灾逐渐进入快速增长阶段,热释放速率急剧上升,最终在590s 左右达到峰值28MW。同时热释放速率曲线表明列车客室火灾在快速增长阶段热释放速率近似呈线性增长,增长速度极快。

图9 不同运行速度下列车火灾热释放速率图Fig.9 Heat release rate diagram of train fire under different running speed

如图10 所示为不同运行速度下列车火灾热释放速率峰值以及相邻两运行速度间的热释放速率峰值增加率曲线。从图中可以看出,30km/h~80km/h运行速度下的列车火灾热释放速率峰值分别为9.56MW、11.2MW、12.9MW、23.3MW、25.8MW和28MW。以60km/h 为界线,运行速度大于等于这一数值时的列车火灾热释放速率峰值超过了低于这一数值时的两倍。同时,热释放速率峰值随着列车运行速度的提高而增大。从图10 中还可以看出,运行速度在30km/h~50km/h 时,火灾热释放速率峰值的增加率较低,即这三种不同运行速度下热释放速率的峰值较为接近,同样的规律在60km/h~80km/h 这个运行速度区间内也有体现。但是,当运行速度由50km/h 提高到60km/h 时,火灾热释放速率峰值急剧增加。这和图9 中热释放速率曲线明显分为两类的结论相符合。

图10 速度-HRR 峰值/HRR 峰值增加率Fig.10 Sspeed-HRR peak/HRR peak increase rate

将热释放速率达到峰值的20%~80%这一阶段视为火灾快速增长阶段,分别得到各运行速度下列车在火灾快速增长阶段的热释放速率增长速率如图11 所示。从图中可以看出,尽管70km/h 时的火灾增长速率比60km/h 时要低,但列车火灾增长速率总体上随着列车运行速度的提高而增大。

图11 速度-火灾增长速率Fig.11 Speed-Fire growth rate

3 结论

通过对高速列车在30km/h~80km/h 时速下运行时的热释放速率进行研究,可以得出以下结论:

(1)列车在30km/h~80km/h 时速下运行时的火灾热释放速率峰值分别为9.56MW、11.2MW、12.9MW、23.3MW、25.8MW 和28MW。

(2)列车火灾热释放速率峰值随着列车运行速度的提高而增大。

(3)运行速度在30km/h~50km/h 时,列车火灾热释放速率的峰值较为接近,同样的规律在60km/h~80km/h 这个运行速度区间内也有体现,当运行速度由50km/h 提高到60km/h 时,火灾热释放速率峰值急剧增加。

(4)尽管70km/h 时的火灾增长速率比60km/h时要低,但列车火灾增长速率总体上随着列车运行速度的提高而增大。

猜你喜欢
峰值数值速率
犊牛生长发育对成年奶牛高峰奶产量和峰值日的影响
体积占比不同的组合式石蜡相变传热数值模拟
数值大小比较“招招鲜”
舰船测风传感器安装位置数值仿真
铝合金加筋板焊接温度场和残余应力数值模拟
锚杆锚固质量等级快速评级方法研究
盘点高考化学反应速率与化学平衡三大考点
化学反应速率与化学平衡考点分析
通过提高心理速率改善记忆