CN113094935A - 一种用于“v”形坡隧道内顶板最高温度预测的方法及应用 - Google Patents
一种用于“v”形坡隧道内顶板最高温度预测的方法及应用 Download PDFInfo
- Publication number
- CN113094935A CN113094935A CN202110507149.9A CN202110507149A CN113094935A CN 113094935 A CN113094935 A CN 113094935A CN 202110507149 A CN202110507149 A CN 202110507149A CN 113094935 A CN113094935 A CN 113094935A
- Authority
- CN
- China
- Prior art keywords
- tunnel
- fire source
- slope
- fire
- highest temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Fire Alarms (AREA)
Abstract
本发明涉及一种用于“V”形坡隧道内顶板最高温度预测的方法及应用。所述最高温度预测方法包括:根据隧道坡度确定火灾场景模拟基础设定参数;改变火源位置与变坡点之间的距离,通过火灾动力学模拟软件进行数值模拟计算,获得不同条件下隧道内顶板最高温度值;根据数值模拟实验所得出的最高温度值,确定所述“V”形坡隧道内顶板最高温度变化相关参数及模拟基础参数对隧道温度变化的影响因子;根据所述相关参数及影响因子,构建“V”形坡隧道内顶板最高温度预测模型。本发明适用于不同火源功率、不同火源位置下,不同坡度的V形坡隧道火灾发生时隧道纵向中心线的顶板最高温度预测。
Description
技术领域
本发明涉及隧道火灾防排烟技术领域,具体涉及一种用于“V”形坡隧道内顶板最高温度预测的方法及应用。
背景技术
由于隧道的狭长封闭结构导致隧道内温度快速上升。据事故数据统计表明,烟气是建筑火灾中对人类生命安全威胁最大的因素,而隧道拱顶下方烟气层的温度会影响到钢筋混凝土的力学性能,当隧道温度达到400℃以上,混凝土抗压强度显著下降,混凝土衬砌结构发生开裂、脱落等现象,不仅阻碍救援行动,还会破坏隧道整体结构,带来更严重的事故后果。
现有研究表明,V形隧道内火灾热释放速率、火源位置及纵向通风等因素对隧道火灾最大烟气温度均有影响。
然而对于隧道内交通事故引发的火灾事故,前人的研究多集中在水平隧道或单坡隧道,对双坡耦合隧道的研究是有限的。同时,由于火源位置的变化会导致双坡耦合隧道内产生“烟气竞争”效应,导致隧道内的温度分布规律发生变化。且与单坡隧道相比V形隧道内火源位置变化对隧道火灾烟气运动状况有一定影响,前人提出的最高温度计算公式并不完全适用于V字形隧道。因而,双坡耦合作用下的“V”形坡隧道火灾最高温度预测模型还需进行进一步的研究,“V”形坡隧道火灾最高温度的准确预测能为预防事故损害提供科学依据。
发明内容
针对现有技术存在的不足,本发明的目的是提供一种用于“V”形坡隧道内顶板最高温度预测的方法,旨在搭建双坡耦合作用下的“V”形坡隧道火灾最高温度预测模型。
为实现上述目的,本发明通过技术方案来实现:一种用于“V”形坡隧道内顶板最高温度预测的方法,包括以下具体步骤:
步骤1,“V”形坡隧道内顶板最高温度影响因素的量纲分析:
V形坡隧道内发生火灾时,影响隧道中心线上最大温升ΔTmax的影响因素包括,恒定影响因素:隧道内纵向风速Vt、隧道高度H、火源所在坡度隧道长度L、环境温度T0、环境空气密度ρ0、空气比定压热容Cp和重力加速度g,偶发影响因素:火源功率Q、火源与V形坡变坡点的距离l、隧道左侧坡度θ1及右侧坡度θ2,V形坡隧道左右两侧坡度相同,即|θ1|=|θ2|=θ*,则关系式如下:
ΔTmax=f(Vt,H,L,Q,T0,ρ0,Cp,g,l,θ*) (6)
从式(6)的11个量中选取Vt、H、L、To、ρ0为基本物理量,根据π定理,可以变为:
f(π1,π2,π3,π4,π5,π6)=0 (7)
通过量纲分析求解得到以下6个无量纲项:
式(9)可以化为:
ΔTmax/T0=f(Q*,l/L,θ*) (11)
式(11)包括无量纲影响因素:无量纲隧道内最大温升ΔTmax/T0、无量纲火源功率Q*、无量纲火源位置l/L、隧道坡度θ*,隧道内最高温度Tmax=ΔTmax+T0;
步骤2,确定模拟火灾场景的基础设定参数,基础设定参数是与经步骤1的量纲分析得到的无量纲项相关的偶发影响因素:隧道坡度θ*、火源位置l和火源功率Q;
步骤3,以步骤2选取的偶发影响因素作为变量使用FDS软件设置多种隧道模型的燃烧工况,并给FDS软件设置网格尺寸;
所述隧道模型是“V”形隧道,“V”形隧道为两端隧道对称布置并且自然通风的隧道,采用混凝土材质,隧道模型还包括参数:隧道模型总长L,隧道横断面宽度、高度H,燃烧工况是不同偶发影响因素组合后的多种工况,所述隧道坡度θ*表示隧道的水平方向变化距离与其垂直方向变化距离的比值,隧道坡度θ*的范围设定为5%-15%,所述火源位置l表示隧道中心线上稳态火源与变坡点之间的距离,火源位置l与“V”型坡单侧隧道长度L的比值l/L为无量纲参数,所述火源功率Q的范围设定为5-10MW;
步骤4,在步骤3设置的各种燃烧工况中,选取不同火源功率和各隧道坡度θ*的组合,每种组合中改变火源与变坡点之间的距离l,通过火灾动力学模拟软件FDS进行各燃烧工况的隧道中心线上温度T的数据模拟计算,提取出各燃烧模型中燃烧过程为准稳态时隧道中心线上最大温度Tmax,将各燃烧工况中得到的隧道中心线上最大温度Tmax转换为步骤1得到的无量纲隧道内最大温升ΔTmax *,
步骤5,根据步骤4获得的拟合系数A和B的数据集,根据所述基础设定参数与拟合系数A和B的关系曲线,得到拟合系数A与隧道坡度θ*和火源功率Q的关系,和拟合系数B与隧道坡度θ*的关系:
A=f(Q*,θ*) (15)
B=f(θ*) (16);
由已知的水平隧道中隧道顶板最高温度与火源功率关系,拟合各种隧道坡度θ*下,拟合系数A与火源功率Q之间的相关性,并使用幂函数公式进行数据拟合,得到多条拟合曲线,每条拟合曲线均由下式表达
A=a·Q*b (17);
将幂指数b取为不同隧道坡度θ*时各拟合曲线幂指数的平均值b’,修正式(17),由式(15)和式(17)的修正式可得A·Q*-b‘=f(θ*) (19);
按式(19)拟合各隧道坡度θ*与修正后的A值之间的关系曲线,得到A·Q*-b’=d·θ*c(20),d、c为常数;
通过将无量纲火源功率Q*项从等式左侧移到右侧,可以将式(20)写为:
A=d·θ*cQ*-b′ (21);
根据式(16),拟合拟合系数B与隧道坡度θ*的关系曲线,得到
B=h+g·θ*(22),h、g为常数;
将式(21)、式(22)代入式(14),整合后得出:
优选地,步骤3和步骤4之间还包括:
式中,D*为火源特征直径;Q为火源功率;T0为环境温度;ρ0为空气密度,取ρ0=1.2kg/m3;cp为空气的定压热容cp=1.02kJ/(kg·K);g为重力加速度,g=9.81m/s2;
步骤3.2,采用步骤3.1计算得到的不同网格尺寸分别作为FDS软件的最小计算单元,利用FDS软件在隧道模型中,以一种隧道坡度θ*、火源功率Q、火源位置l作为计算条件,计算得到不同网格尺寸对应的火源位置l处的竖向温度分布曲线,选取重合度高的多条竖向温度分布曲线中的网格尺寸的最大值作为待验证网格尺寸;
步骤3.3,采用待验证网格尺寸作为FDS软件第最小计算单元,利用FDS软件在隧道模型中对步骤3选取的一种燃烧工况进行模拟,得到顶板最高温度模拟数据,使用现有的坡度隧道火灾中温度顶板最高温度预测模型对相同燃烧工况进行模拟得到预测值,将预测值与模拟数据进行对比,若两种数据吻合,则确定待验证网格尺寸作为FDS软件的最小计算单元。
进一步地,所述网格尺寸为0.25m。
进一步地,所述隧道模型总长600m,隧道横断面宽12m、高5m。
进一步地,所述“V”形坡隧道内顶板最高温度预测模型为:
进一步地,上述方法得到的“V”形坡隧道内顶板最高温度预测模型应用于坡度范围为5%~15%,及火源功率范围为5~10MW的“V”形坡隧道内顶板最高温度值预测。
有益效果
本发明是一种用于“V”形坡隧道内顶板最高温度预测的方法,综合考虑了双坡耦合隧道的特殊性,确定了火源功率Q、隧道坡度及火源位置对隧道内烟气温度的影响,适用于多种不同坡度、不同火源功率及不同火灾位置下的“V”形坡隧道内顶板最高温度预测。通过火灾模拟软件仿真模拟实验发现,最高温度数值与模型预测值吻合度较高,验证了此预测模型的科学性与准确性,补充了前人隧道火灾研究中未考虑的“V”形坡隧道情况,完善了坡度隧道内顶板最高温度的预测公式。
附图说明
表1为本发明中的数值模拟工况汇总表;
图1为本发明用于“V”形坡隧道内顶板最高温度预测的方法的流程图;
图2为本发明中的“V”形坡隧道火灾数值模拟示意图;
图3为本发明中5MW火源功率下隧道火源位置变化时顶板最大无量纲温升曲线图;
图4为本发明中7.5MW火源功率下隧道火源位置变化时顶板最大无量纲温升曲线图;
图5为本发明中10MW火源功率下隧道火源位置变化时顶板最大无量纲温升曲线图;
图6为本发明中隧道坡度与系数A的关系图;
图7为本发明中隧道坡度与系数B的关系图;
图8为本发明中不同隧道坡度下A与Q*的变化关系曲线图;
图9为本发明中隧道坡度与修正后的A值之间的关系曲线图;
图10为本发明中FDS模拟值与计算值对比图;
图11为不同网格尺寸下距火源25m处的竖向温度分布曲线;
图12为坡度隧道火灾时预测模型预测结果与燃烧模型模拟实验结果的对比图。
附图中,各标号代表的不见列表如下:
1、隧道,2、火源。
具体实施方式
为使本发明实现的原理和特征易于了解,下面将结合具体实施方式进行举例,进一步阐述本发明。
如图1至图10所示,本发明提供一种用于“V”形坡隧道内顶板最高温度预测的方法,包括以下具体步骤:
步骤1,“V”形坡隧道内顶板最高温度影响因素的量纲分析:
V形坡隧道内发生火灾时,影响隧道中心线上最大温升ΔTmax的影响因素包括,恒定影响因素:隧道内纵向风速Vt、隧道高度H、火源所在坡度隧道长度L、环境温度T0、环境空气密度ρ0、空气比定压热容Cp和重力加速度g,偶发影响因素:火源功率Q、火源与V形坡变坡点的距离l、隧道左侧坡度θ1及右侧坡度θ2,由于本研究基于V形坡隧道,左右两侧坡度相同,即|θ1|=|θ2|=θ*,则可列关系式如下:
ΔTmax=f(Vt,H,L,Q,T0,ρ0,Cp,g,l,θ*) (6)
式(6)的11个量中,选取Vt、H、L、T0、ρ0为基本物理量,根据π定理,可以变为:
f(π1,π2,π3,π4,π5,π6)=0 (7)
其中,
通过量纲分析求解得到以下6个无量纲项:
其中隧道内风速Vt为自然通风引起的风流,式(6)可以变为:
式(9)可以化为:
ΔTmax/T0=f(Q*,l/L,θ*) (11)
式(11)给出了无量纲影响因素:无量纲隧道内最大温升ΔTmax/T0与无量纲火源功率Q*、无量纲火源位置l/L以及隧道坡度θ*相关,因此只需要确定ΔTmax/T0与每一个参数之间的关系,即可得出V形坡隧道内发生火灾时,隧道中心线上最大温升的无量纲计算公式,而隧道内最高温度Tmax=ΔTmax+T0。
步骤2,确定模拟火灾场景的基础设定参数,基础设定参数是与经步骤1的量纲分析得到的无量纲项相关的偶发影响因素:如表1中的:隧道坡度θ*、火源位置l和火源功率Q,表1中的坡度为左右两侧隧道的坡度。
步骤3,以步骤2选取的偶发影响因素作为变量使用FDS软件设置多种隧道模型的燃烧工况,并给FDS软件设置网格尺寸。
隧道模型依据武汉黄鹤楼隧道建立,采用模拟尺寸与实体隧道尺寸比例为1:1,“V”形坡隧道1变坡点两侧隧道对称倾斜且长度相同,隧道模型总长600m,隧道横断面宽12m、高5m,即设置为“CONCRETE”属性。所述隧道两端口开启方式均设置为“OPEN”,即为自然通风状态,所述隧道坡度θ表示隧道的水平方向变化距离与其垂直方向变化距离的比值,隧道坡度θ的范围设定为5%-15%,如图2所述火源2位置l表示隧道中心线上,火源与变坡点之间的距离设置为0~250m,烟流可分为火源的上游和下游部分,火源位置l与“V”型坡单侧隧道长度L的比值l/L为无量纲参数,在城市道路隧道中,客车或厢式货车是最常见的交通工具。一辆小汽车火灾的放热速率约为5MW,2-3辆小汽车碰撞或1辆客车的放热速率为8-15MW,因此,所述火源功率Q的范围设定为5-10MW;设置火源为稳态火,在模拟过程中放热速率是恒定的。隧道模型内环境温度等物理量在运行200s后达到准稳态。因此提取280~300s时的平均值进行代表准稳态的平均值。
选取3种火源功率、5种隧道坡度、9种火源位置。共设置145组燃烧工况,如表1所示。
表1为燃烧工况基础设定参数汇总表
使用FDS软件进行燃烧模拟前以火源特征直径D*为参照进行网格划分,网格为燃烧模型的最小计算单元,网格尺寸δx设置为1/16D*~1/4D*,
式中,D*为火源特征直径(m);Q为火源功率(kW);T0为环境温度(K);ρ0为空气密度(kg/m3),取ρ0=1.2kg/m3;cp为空气的定压热容(kJ/(kg·K),取cp=1.02kJ/(kg·K);g为重力加速度(m/s2),取g=9.81m/s2。
选取以下燃烧工况:上游坡度为-5%,下游坡度为5%,火源功率为Q=5MW,根据式(12)可得Q=5MW时的合适网格尺寸δx为0.07m~0.29m,选取δx=0.1m、0.125m、0.2m、0.25m、0.3m这5种网格尺寸作为最小计算单元,使用FDS在隧道模型中对选取的燃烧工况进行模拟,对得到的下距火源25m处的竖向温度进行比较,竖向温度分布曲线如图11所示。由图11可知当网格尺寸逐渐减小时,温度曲线逐渐重合;网格尺寸为0.25m时,其温度曲线与0.1、0.125m、0.2m的曲线几乎一致,考虑到计算的准确度以及计算耗时,火源上下游30m范围内的网格尺寸选为0.25m。
为进一步验证选取网格尺寸后得到的燃烧模型的准确性,将燃烧模型得到的模拟数据与Yi等人开发的坡度隧道火灾中温度顶板最高温度预测模型的预测值进行对比。
ks=1-0.065α,式中,ks是坡度修正系数,α是隧道的斜率。
从式(4)可以得到火源位于坡度隧道时顶板最高温度预测值,图12为预测模型预测结果与燃烧模型模拟实验结果的对比图。结果表明,Yi等经验公式预测的顶板最高温度与FDS模拟实验结果吻合较好。
步骤4,在步骤3设置的各种燃烧工况中,选取不同火源功率和各隧道坡度θ*的组合,每种组合中改变火源与变坡点之间的距离l,通过火灾动力学模拟软件FDS进行各燃烧工况的隧道中心线上温度T的数据模拟计算,提取出各燃烧模型中燃烧过程为准稳态时隧道中心线上最大温度Tmax,将各燃烧工况中得到的隧道中心线上最大温度Tmax转换为步骤1得到的无量纲隧道内最大温升ΔTmax *,获得不同火源功率条件下各隧道坡度θ*对应的无量纲隧道内最大温升ΔTmax *与l/L的关系曲线,如图3~5,根据该关系曲线的形态,可知该曲线对应指数函数模型,得到拟合方程:其中A和B是拟合系数;
步骤5,根据步骤4获得的拟合系数A和B的数据集,如表2,表2中不同的A和B拟合值,其中R是模拟数据和拟合曲线之间的皮尔逊相关系数。由表2可知,A、B的所有值均为正,这表明最大无量纲温升随着无量纲火源位置参数的增大而增大。此外,所有的相关系数R均大于0.95,拟合曲线的相关性较好;
表2变量A、B拟合值
A和B值与隧道坡度的关系如图6和图7所示。A值在一定的放热速率下随隧道坡度增加而减小。B值在一定的放热率下随隧道坡度的增加而增加,并在一定的斜率下、不同的放热率之间相对恒定。总之,从图6和图7中可以看出,系数A的变化与隧道坡度和放热率有关,而系数B是仅与隧道坡度有关的函数。因此,这两个系数可以表示为:
A=f(Q*,θ*) (15)
B=f(θ*) (16)
由前人的水平隧道中隧道顶板最高温度与火源功率关系的研究可知,隧道顶板最高温度与火源功率有关,因此,首先拟合各种隧道坡度下,系数A与放热率之间的相关性,并使用幂函数公式对数据进行拟合,得到拟合曲线图8。考虑到系数A的无量纲特性,使用式(17)将参数A与无量纲放热率数据相关联,得到:
A=a·Q*b (17)
如图8,幂指数b的参数值相差较小,令式(17)中b值取不同坡度时的平均值0.77,式(17)修正为:
A=a·Q*0.77 (18)
由式(15)、式(18)可以得出:
A·Q*-0.77=f(θ*) (19)
按式(19)隧道坡度与修正后的A值之间的关系得到曲线图9;
将图9与图6进行比较可以看出,通过无量纲放热率修正后的A值从完全不同的值变得更为相似,且图9中的拟合曲线可以表示为:
A·Q*-0.77=0.47·θ*-0.7 (20)
通过将无量纲火源热释放速率Q*项从等式左侧移到右侧,可以将式(20)写为:
A=0.47·θ*-0.7Q*0.77 (21)
经计算得出式(21)与表2中A值的相关系数为0.952,拟合相关性较好。
根据式(16),B的值仅与隧道坡度有关。
拟合图7中的数据点得出:
B=1.2+7.9·θ* (22)
将式(21)、式(22)代入式(14),整合后得出:
式(23)表明,式中包括相关参数:T0为环境温度;ρ0为环境空气密度;Cp为空气比定压热容;g为重力加速度;H为隧道高度。
无量纲最大温升参数与无量纲放热速率的0.8次方成正比,且随无量纲火源位置的增大而增大,与隧道坡度的关系是非线性和非单调的。由于没有对V形隧道顶板最大温升进行任何实体试验,因此使用等式(23)对无量纲最大温升进行计算,并将其与FDS模拟温升值进行比较,比较结果如图10所示。从图10中可以看出,数据点分散在等效线附近,式(23)与FDS模拟值基本一致,且拟合得出其相关系数R2均在0.95以上,式(23)在预测非变坡点处火灾时的隧道内最大温升效果较好。
结果表明式(23)该公式可用于5%~15%坡度范围及5~10MW火源功率范围内的“V”形坡隧道内顶板最高温度值预测。
以上所述仅为本发明的较佳实施案例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、同等替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种用于“V”形坡隧道内顶板最高温度预测的方法,其特征在于,包括以下具体步骤:
步骤1,“V”形坡隧道内顶板最高温度影响因素的量纲分析:
V形坡隧道内发生火灾时,影响隧道中心线上最大温升ΔTmax的影响因素包括,恒定影响因素:隧道内纵向风速Vt、隧道高度H、火源所在坡度隧道长度L、环境温度T0、环境空气密度ρ0、空气比定压热容Cp和重力加速度g,偶发影响因素:火源功率Q、火源与V形坡变坡点的距离l、隧道左侧坡度θ1及右侧坡度θ2,V形坡隧道左右两侧坡度相同,即|θ1|=|θ2|=θ*,则关系式如下:
ΔTmax=f(Vt,H,L,Q,T0,ρ0,Cp,g,l,θ*) (6)
从式(6)的11个量中选取Vt、H、L、T0、ρ0为基本物理量,根据π定理,可以变为:
f(π1,π2,π3,π4,π5,π6)=0 (7)
通过量纲分析求解得到以下6个无量纲项:
式(9)可以化为:
ΔTmax/T0=f(Q*,l/L,θ*) (11)
式(11)包括无量纲影响因素:无量纲隧道内最大温升ΔTmax/T0、无量纲火源功率Q*、无量纲火源位置l/L、隧道坡度θ*,隧道内最高温度Tmax=ΔTmax+T0;
步骤2,确定模拟火灾场景的基础设定参数,基础设定参数是与经步骤1的量纲分析得到的无量纲项相关的偶发影响因素:隧道坡度θ*、火源位置l和火源功率Q;
步骤3,以步骤2选取的偶发影响因素作为变量使用FDS软件设置多种隧道模型的燃烧工况,并给FDS软件设置网格尺寸;
所述隧道模型是“V”形隧道,“V”形隧道为两端隧道对称布置并且自然通风的隧道,采用混凝土材质,隧道模型还包括参数:隧道模型总长L,隧道横断面宽度、高度H,燃烧工况是不同偶发影响因素组合后的多种工况,所述隧道坡度θ*表示隧道的水平方向变化距离与其垂直方向变化距离的比值,隧道坡度θ*的范围设定为5%-15%,所述火源位置l表示隧道中心线上稳态火源与变坡点之间的距离,火源位置l与“V”型坡单侧隧道长度L的比值l/L为无量纲参数,所述火源功率Q的范围设定为5-10MW;
步骤4,在步骤3设置的各种燃烧工况中,选取不同火源功率和各隧道坡度θ*的组合,每种组合中改变火源与变坡点之间的距离l,通过火灾动力学模拟软件FDS进行各燃烧工况的隧道中心线上温度T的数据模拟计算,提取出各燃烧模型中燃烧过程为准稳态时隧道中心线上最大温度Tmax,将各燃烧工况中得到的隧道中心线上最大温度Tmax转换为步骤1得到的无量纲隧道内最大温升ΔTmax *,获得不同火源功率条件下各隧道坡度θ对应的无量纲隧道内最大温升ΔTmax *与l/L的关系曲线,根据该关系曲线的形态,可知该曲线对应指数函数模型,得到拟合方程:其中A和B是拟合系数;
步骤5,根据步骤4获得的拟合系数A和B的数据集,根据所述基础设定参数与拟合系数A和B的关系曲线,得到拟合系数A与隧道坡度θ*和火源功率Q的关系,和拟合系数B与隧道坡度θ*的关系:
A=f(Q*,θ*) (15)
B=f(θ*) (16);
由已知的水平隧道中隧道顶板最高温度与火源功率关系,拟合各种隧道坡度θ*下,拟合系数A与火源功率Q之间的相关性,并使用幂函数公式进行数据拟合,得到多条拟合曲线,每条拟合曲线均由下式表达
A=a·Q*b (17);
将幂指数b取为不同隧道坡度θ*时各拟合曲线幂指数的平均值b’,修正式(17),由式(15)和式(17)的修正式可得A·Q*-b′=f(θ*) (19);
按式(19)拟合各隧道坡度θ*与修正后的A值之间的关系曲线,得到A·Q*-b′=d·θ*c(20),d、c为常数;
通过将无量纲火源功率Q*项从等式左侧移到右侧,可以将式(20)写为:
A=d·θ*cQ*-b′ (21);
根据式(16),拟合拟合系数B与隧道坡度θ*的关系曲线,得到
B=h+g·θ* (22),h、g为常数;
将式(21)、式(22)代入式(14),整合后得出“V”形坡隧道内顶板最高温度预测模型:
2.根据权利要求1所述的用于“V”形坡隧道内顶板最高温度预测的方法,其特征在于,步骤3和步骤4之间还包括:
式中,D*为火源特征直径;Q为火源功率;T0为环境温度;ρ0为空气密度,取ρ0=1.2kg/m3;cp为空气的定压热容cp=1.02kJ/(kg·K);g为重力加速度,g=9.81m/s2;
步骤3.2,采用步骤3.1计算得到的不同网格尺寸分别作为FDS软件的最小计算单元,利用FDS软件在隧道模型中,以一种隧道坡度θ*、火源功率Q、火源位置l作为计算条件,计算得到不同网格尺寸对应的火源位置l处的竖向温度分布曲线,选取重合度高的多条竖向温度分布曲线中的网格尺寸的最大值作为待验证网格尺寸;
步骤3.3,采用待验证网格尺寸作为FDS软件第最小计算单元,利用FDS软件在隧道模型中对步骤3选取的一种燃烧工况进行模拟,得到顶板最高温度模拟数据,使用现有的坡度隧道火灾中温度顶板最高温度预测模型对相同燃烧工况进行模拟得到预测值,将预测值与模拟数据进行对比,若两种数据吻合,则确定待验证网格尺寸作为FDS软件的最小计算单元。
3.根据权利要求2所述的用于“V”形坡隧道内顶板最高温度预测的方法,其特征在于,所述网格尺寸为0.25m。
4.根据权利要求3所述的用于“V”形坡隧道内顶板最高温度预测的方法,其特征在于,所述隧道模型总长600m,隧道横断面宽12m、高5m。
6.如权利要求1~5的方法得到的“V”形坡隧道内顶板最高温度预测模型应用于坡度范围为5%~15%,及火源功率范围为5~10MW的“V”形坡隧道内顶板最高温度值预测。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110013389 | 2021-01-06 | ||
CN2021100133893 | 2021-01-06 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113094935A true CN113094935A (zh) | 2021-07-09 |
CN113094935B CN113094935B (zh) | 2022-04-08 |
Family
ID=76664878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110507149.9A Active CN113094935B (zh) | 2021-01-06 | 2021-05-10 | 一种用于“v”形坡隧道内顶板最高温度预测的方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113094935B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113653535A (zh) * | 2021-09-16 | 2021-11-16 | 招商局重庆交通科研设计院有限公司 | 一种防火板保护下沉管隧道壁厚方向火灾温度场估计方法 |
CN113669740A (zh) * | 2021-08-23 | 2021-11-19 | 清华大学 | 一种预测燃烧室内局部最高温度的方法 |
CN113779804A (zh) * | 2021-09-16 | 2021-12-10 | 招商局重庆交通科研设计院有限公司 | 一种沉管防火板安装空腔高度设计方法 |
CN114492248A (zh) * | 2022-02-25 | 2022-05-13 | 西南交通大学 | 一种蛇型曲线隧道火灾纵向温度的确定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2405438A1 (en) * | 2010-07-07 | 2012-01-11 | Crocus Technology S.A. | Method for writing in a MRAM-based memory device with reduced power consumption |
CN110887583A (zh) * | 2019-11-07 | 2020-03-17 | 武汉科技大学 | 一种用于水下“v”形坡隧道内温度分布预测的方法 |
CN111027176A (zh) * | 2019-11-07 | 2020-04-17 | 武汉科技大学 | 一种隧道侧向重点排烟火灾烟气逆流长度的计算方法 |
CN111477085A (zh) * | 2020-04-21 | 2020-07-31 | 合肥工业大学 | 一种分岔隧道群的火灾燃烧和污染物窜流测量系统 |
-
2021
- 2021-05-10 CN CN202110507149.9A patent/CN113094935B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2405438A1 (en) * | 2010-07-07 | 2012-01-11 | Crocus Technology S.A. | Method for writing in a MRAM-based memory device with reduced power consumption |
CN110887583A (zh) * | 2019-11-07 | 2020-03-17 | 武汉科技大学 | 一种用于水下“v”形坡隧道内温度分布预测的方法 |
CN111027176A (zh) * | 2019-11-07 | 2020-04-17 | 武汉科技大学 | 一种隧道侧向重点排烟火灾烟气逆流长度的计算方法 |
CN111477085A (zh) * | 2020-04-21 | 2020-07-31 | 合肥工业大学 | 一种分岔隧道群的火灾燃烧和污染物窜流测量系统 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113669740A (zh) * | 2021-08-23 | 2021-11-19 | 清华大学 | 一种预测燃烧室内局部最高温度的方法 |
CN113669740B (zh) * | 2021-08-23 | 2022-08-12 | 清华大学 | 一种预测燃烧室内局部最高温度的方法 |
CN113653535A (zh) * | 2021-09-16 | 2021-11-16 | 招商局重庆交通科研设计院有限公司 | 一种防火板保护下沉管隧道壁厚方向火灾温度场估计方法 |
CN113779804A (zh) * | 2021-09-16 | 2021-12-10 | 招商局重庆交通科研设计院有限公司 | 一种沉管防火板安装空腔高度设计方法 |
CN114492248A (zh) * | 2022-02-25 | 2022-05-13 | 西南交通大学 | 一种蛇型曲线隧道火灾纵向温度的确定方法 |
CN114492248B (zh) * | 2022-02-25 | 2024-04-05 | 西南交通大学 | 一种蛇型曲线隧道火灾纵向温度的确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113094935B (zh) | 2022-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113094935B (zh) | 一种用于“v”形坡隧道内顶板最高温度预测的方法及应用 | |
Zhao et al. | An experimental study on smoke back-layering and critical velocity in tunnel fires with canyon cross wind | |
Bari et al. | Simulation of airflow and pollution levels caused by severe traffic jam in a road tunnel | |
Yang et al. | Multiple patterns of heat and mass flow induced by the competition of forced longitudinal ventilation and stack effect in sloping tunnels | |
CN111075493B (zh) | 一种高原隧道火灾烟气控制风速设计方法 | |
Zhao et al. | Re-direction of smoke flow in inclined tunnel fires | |
CN113849998B (zh) | 一种基于bim的隧道火灾预警处理方法及系统 | |
Ko et al. | Determination of the applicable exhaust airflow rate through a ventilation shaft in the case of road tunnel fires | |
Xu et al. | Study on smoke exhaust performance in tunnel fires based on heat and smoke exhaust efficiency under the lateral centralized mode | |
Zhao et al. | Numerical investigation on fire-induced indoor and outdoor air pollutant dispersion in an idealized urban street canyon | |
Kong et al. | A numerical study on smoke behaviors in inclined tunnel fires under natural ventilation | |
Ma et al. | Characteristics of fires in coal mine roadways and comparative analysis of control effectiveness between longitudinal ventilation and cross-section sealing | |
Zhu et al. | Rapid prediction of fire-induced heat flux on the liners of horseshoe and circular tunnels with longitudinal ventilation at critical velocity | |
Zhao et al. | Experimental study on the characteristics of fire smoke longitudinal crossflow in adjacent tunnels | |
Yang et al. | Numerical study on the smoke flow characterization and phenomenon of plug-holing under lateral smoke exhaust in tunnel fire | |
Xu et al. | Investigating the influence of outdoor temperature variations on fire-induced smoke behavior in an atrium-type underground metro station using hybrid ventilation systems | |
Wu et al. | Spread behavior of spill fire under lateral wind and temperature distribution of PE sandwich panels on facade | |
Berkoe et al. | CFD Modeling of the Fjardaal smelter potroom ventilation | |
Lo et al. | A numerical study of external smoke spread in designated refuge floor | |
Zhao et al. | Effects of ambient pressure on fire-induced buoyancy driven plume dispersion and re-entrainment behavior in a street canyon | |
Benkoussas et al. | Numerical Security Assessment in Case of Fire in Underground Transport Spaces. | |
Yu et al. | Study on the combination effect of tunnel slope and longitudinal fire location on the asymmetric flow fields in a naturally ventilated tunnel | |
Chow et al. | Assessment of alternative ventilation schemes at public transport interchange | |
Fan et al. | Prediction of back-layering length in subway tunnel with on-fire train running | |
Zhu et al. | Numerical simulation & scene experiment of metro tunnel with crossover interval under fire conditions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |