APP下载

近塔桥面火灾对三塔悬索桥结构性能影响研究

2017-06-19马如进崔传杰马明雷

关键词:中塔火源温度场

马如进,崔传杰,马明雷

(同济大学桥梁工程系, 上海 200092)



近塔桥面火灾对三塔悬索桥结构性能影响研究

以某三塔悬索桥中间钢塔为例,研究了火灾对钢桥塔和桥梁结构性能的影响.首先,通过火灾场景数值模拟,得出桥面近塔区域不同类型火灾场景下的火焰温度分布规律,以及钢结构桥塔的温度分布特性.其次,通过非线性分析获得了三塔悬索桥在不同的近塔火灾场景下结构静力性能的变化.结果表明,在大型车辆火灾作用下,钢中塔有超过80 m2的区域温度超过800 ℃,最高温度达到1 000 ℃以上,中塔性能受到显著影响.钢中塔的应力、变形发生较为明显的变化:钢塔产生4.6 mm的竖向残余变形和53.8 mm的侧向残余变形,火灾区域应力折减达70 MPa,而主缆、主梁的应力和变形变化较小.

结构分析;三塔悬索桥;中间钢塔;火灾;管理养护

桥梁火灾是桥梁结构在其寿命周期内常见的和严重的灾害之一.根据纽约交通局2008年对1 746起桥梁事故的调查结果,其中有52座桥梁因火灾事故遭到破坏,而由于地震导致破坏的仅有19座[1].近年来桥梁火灾事故也呈现增加趋势,对桥梁的运营安全带来了很大的挑战.2007年4月,连接美国旧金山和奥克兰两大城市的高速公路上发生油罐车火灾,造成旧金山地区一座钢梁桥熔化倒塌;2014年10月在建的赤石特大桥塔内火灾引起9根斜拉索断裂导致桥面一侧下沉.可见针对桥梁抗火的研究应引起设计者与运营管理者的充分重视.

自1957年Blinov等[2]率先开展火灾实验研究以来,国内外很多学者对建筑火灾进行了相关方面的理论研究,且建筑火灾及抗火设计规范日益成熟.而相比于建筑领域,桥梁火灾事故的理论研究却是刚刚起步.近十年来针对桥梁火灾进行的研究有:Mendes等[3]通过二维平面模型模拟了Vasco da Gama大桥船舶火灾事故的发生过程;Dotreppe等[4]使用计算机程序SAFIR对比利时Vivegnis大桥火灾事故做了数值模拟;Kodur等[5]使用计算机程序SAFIR针对钢-混凝土组合结构桥梁的加劲梁防火性能做了案例分析;IGNACIO等[6]采用LUSAS软件研究了由美国联邦公路管理局设计的12.20 m简支梁桥的抗火性能.长安大学李利军则尝试进行了公路火灾温度场的模拟,并对大跨缆索承重桥梁进行了火灾分析.但这些研究多着眼于主梁和缆索系统火灾场景下的分析,而对缆索承重桥梁关键承力构件的桥塔,还缺乏深入研究.

本文在总结既有桥梁火灾研究工作的基础上,以泰州长江大桥为工程实例,对桥梁的常见火灾场景进行了基于FDS(Fire Dynamics Simulator)的数值模拟,在此基础上求解了中间钢桥塔受到火灾影响时的温度场分布,研究了各类火灾场景对全桥结构静力性能的影响.泰州长江大桥为三塔特大跨悬索桥,中塔的力学性能对全桥结构至关重要,本文研究结果对了解钢塔结构在火灾作用下的性能变化以及对全桥受力行为的影响有较好的参考意义.

1 桥梁火灾数值模拟

1.1 火源模型和火灾场景

桥梁火灾多由车辆事故引起,属于非稳态模型,相应火灾发展过程大致要经历3个阶段:初始增长阶段、稳定燃烧阶段和减弱阶段.目前国际上常用的非稳态火源模型有3种[7].本次研究中采用Heskestad所提出的t2火源模型(图1)[8].设火灾持续时间为tmax,热释放速率在td时刻达到最大值Qd,由于实际分析中往往只需关心初始增长和稳定燃烧两个阶段,该火源模型可表述为:

(1)

式中:Q为热释放速率,kW;α为增长系数,kW/s2;t为时间,s.由于桥梁火灾的相关试验研究和统计资料仍处在起步阶段,而桥梁火灾与隧道火灾在火源类型和火灾场景等诸多方面具有相似之处,本文借鉴了隧道火灾在火源模型和火灾场景等方面的研究成果.根据车辆的类型与燃烧生成的热量分类,将火灾分为小汽车火灾、客车火灾、载重卡车和油罐车火灾4类[9].Cafaro等人[10]通过大量隧道火灾试验发现小汽车火灾热释放速率为4~8 MW,增长系数可取为0.010~0.012 kW/s2;客车火灾热释放速率为15~30 MW,增长系数可取0.10~0.15 kW/s2;载重卡车火灾热释放速率为15~150 MW,增长系数取0.50 kW/s2.而油罐车火灾由于火源庞大,相关试验研究较少,但一般认为在5~10 min内即可达到最高热释放速率100~300 MW[11].对于火灾持续时间,法国CETU(隧道研究中心)[12]对小汽车、客车、载重卡车、油罐车的建议取值分别为2 700 s,5 400 s,6 000 s和6 000 s,Haack在一项隧道火灾的专题研究中则建议对油罐车火灾取7 200 s的燃烧时间[13].对于火源尺寸,Heselden[14]研究认为小汽车、客车、载重卡车、油罐车的尺寸可分别取为1.5 m×4 m,2 m×6 m~2.5 m×10 m,4 m×6 m,4 m×12 m;世界道路协会常设协会(PIARC)[15]规定小汽车、客车和油罐车的火焰尺寸分别为1.5 m×4 m,2 m×6 m,2.5 m×10 m.结合上述研究成果,4类火灾的参数取值如表1所示,其中场景A,B,C和D分别对应一辆普通小汽车、一辆公共汽车、一辆大型载重卡车以及一辆油罐车的燃烧场景.

图1 非稳态火源火灾场景的发展过程

火灾场景增长系数/(kW·s-2)持续时间/s热释放速率/MW火源尺寸/(m×m)A0.011270051.5×4B0.1235400302×6C0.56000804×6D172002004×12

1.2 火源模型数值模拟

本研究对火源模型和火灾场景的模拟采用基于CFD理论的FDS软件进行.FDS基于大涡模拟理论,以火灾中的流体运动为主要模拟对象,通过满足火灾场景下的质量守恒、动量守恒、能量守恒和组分守恒,结合相关状态方程和压力方程,可较好地用于开放环境下的桥梁火灾模拟[16].FDS在火灾动力模拟中的可靠性已经得到许多证实[16-17].将火灾源设置在泰州大桥钢中塔附近,在模拟过程中假定周边边界氧气供应充足且周围不存在绝热障碍物.通过4类火灾场景的模拟,可以得到焰心不同高度处的温度分布,图2给出了场景A和D的火焰中心温度分布.图3给出了场景D的火灾场景的模拟效果.

从分析可知,场景A火焰根部温度约为700 ℃,随着火焰高度增加,温度迅速降低,5 m处火焰温度已经降低到200 ℃左右.场景B,C,D火焰最高温度相仿,均为1 000~1 100 ℃左右.但随着火焰高度增加,其温度降低幅度存在较大差异,5 m高度时场景B已低于600 ℃,场景C,D为800 ℃左右.场景D与场景C各高度处火焰温度均相差不大,但场景D火焰尺寸(4×12 m)远大于场景C(4×6 m),可以预见火灾场景D对结构可能造成的破坏最大.

t/s

t/s

图3 火灾场景D模拟展示

1.3 FDS数值模拟可靠性验证

目前国际上广泛采用的火灾升温曲线有ISO834曲线[18]、ASEM-EI19曲线[19]、HC曲线[20]等.本文将数值模拟的最高温度与3种升温曲线的计算结果进行了比较,图4给出了场景A与场景D的比较结果.通过比较可以发现,数值模拟整体上具有较高的吻合度.但由于ISO834曲线和ASEM-EI19曲线主要是针对建筑火灾制定的,相比较处在开放环境下、氧气供应充足的桥梁火灾,在大型火灾场景下的温度计算结果偏低;HC曲线主要应用在石油化工领域,计算曲线考虑了火灾初期的爆炸和热冲击效应,对场景A这类小型火灾,计算结果明显偏高.总体看来模拟结果可以较合理地反映类似桥面特征的开放式空间在火灾作用下的温度场特征.

t/s

t/s

2 火灾下中塔温度场分布

对处于开放环境下的桥梁火灾而言,近塔桥面火灾主要通过热对流和热辐射的形式对桥塔产生影响[16].热对流是由于流体的运动而传递的热量,对流可分自然对流和强迫对流两种,桥梁火灾对流形式为自然对流,是由于温度不均匀而引起的.其传热满足[21]:

q″=hΔT

(2)

式中:q″为单位面积的热流量,W/m2;h为热对流传热系数,W/(m2·K);ΔT为温差,K.

热辐射是由于物体内部微观粒子的热运动(或者说由于物体自身的温度)而使物体向外发射辐射能的现象.桥梁火灾处在开放环境中,烟气容易扩散,受热辐射影响较大,其净热量传递可以用斯蒂芬-波尔兹曼方程计算:

(3)

式中:Q为热流率;ε为吸收率(黑度);σ为斯蒂芬-波尔兹曼常数,W/(m2·K4);A1为辐射面1的面积,m2;F12为从辐射面1到辐射面2的形状系数;T1,T2分别为辐射面1,2的表面温度,K.

由式(2),式(3)可以看出,热对流与热辐射的热量传递均与当前时刻的表面温度有关,因此火源-结构传热模型的建立是一个复杂的时程分析,目前国内外多采用数值模拟的手段,尚无针对火灾-结构传热模型的简化或经验公式研究.本文借助Ansys有限元分析软件的热力学模块,选取受火灾影响的钢中塔节段以shell131单元进行温度场模拟.分析中假设火灾作用于桥塔顺桥向中心处,火焰外侧边缘距离桥塔0.5 m.钢结构在高温下的导热系数和比热按照欧洲规范EC3[22]及EC4[23]的规定取值,如式(4),式(5)所示:

(4)

Cs=

(5)

式中:λs,Cs和Ts分别表示钢材的导热系数、比热和温度.

通过分析得到4类火灾场景下在中塔部分节段最大温度场分布如图5所示.

图5 4类火灾场景下中塔温度场分布对比图

从图5可以看出,4类火灾场景下中塔的温度场分布规律基本与火源模型的模拟结果相一致.相比之下,场景A,B的温度场受影响程度远小于场景C,D.场景C,D在桥塔中心处沿高度方向的分布规律基本一致,但场景D在顺桥向的影响范围远大于场景C,场景D的影响范围会波及到桥塔侧表面,中塔有超过80 m2范围内温度超过800 ℃.

3 火灾作用前后结构静力分析

3.1 高温下钢结构性能

目前,欧洲EC3[22],EC4[23]、美国AISC(American Institute of Steel Construction)[24]、英国BS5950等[25]国家规范都对高温作用下钢结构的材料特性作出了规定.其中EC3和EC4的方法计算简单、适用性强,本研究对钢结构材料特性的计算以此规范为参考进行.

EC3和EC4中定义的屈服强度和弹性模量随温度变化结果如图6所示.

温度/℃

由图6可以看出,当温度达到800 ℃时,钢结构的弹性模量与屈服强度几乎只剩10%.结合图5所示的温度场分析结果,可以预见车辆火灾发生时,特别是场景C,D这类大型火灾,会对结构产生较强烈的影响.

3.2 结构静力性能分析

针对三塔两跨悬索桥,钢中塔是结构的关键.在考虑火灾影响下,传统的梁式模型不能考虑其局部影响效应.为此,本文将钢中塔采用板壳单元模拟,如图7所示.考虑到火灾发生时,一般桥梁管理中严控车辆通行,因此本文重点比较了火灾影响时恒载下结构行为的变化.

图7 泰州长江大桥有限元模型图

图8和表2给出了4类火灾场景下,中塔的变形趋势以及各关键位置的位移增量数值,其中va,ha代表远离火灾一侧的塔顶竖向和侧向位移,vb,hb代表火灾发生一侧的塔顶竖向和侧向位移.

图8 火灾下中塔变形趋势图

位移增量/mm场景A场景B场景C场景Dva-0.1-1.9-3.7-4.6vb0.11.73.64.5ha0.321.748.052.9hb0.422.849.053.8

以火灾场景D为例,将中塔位移按照150倍比例放大,单独观察桥塔受火面附近的变形和应力变化,如图9所示.

可以看到,受火灾影响区域由于材料性能受损,其应力值相比周边单元大幅降低,受火灾影响区域的变形也比未受火灾的一侧桥塔显著得多,在火灾作用下,受火面横桥向变形明显呈现出外凸的趋势,而中塔整体上呈现向远离火灾一侧倾斜的变形趋势.火灾场景D作用下,中塔最大竖向位移增量4.6 mm,

最大横桥向位移增量53.8 mm.温度最高区域应力降低约70 MPa,中塔静力行为受到较为明显的影响.

图9 场景D中塔火灾区域应力云图及变形图

分别提取不同火灾形式下加劲梁竖向挠度增量图,如图10和图11所示.图12给出了主缆不同位置处的应力增量变化图.

主梁顺桥向坐标(m,0代表中塔中心处)

主梁顺桥向坐标(m,0代表中塔中心处)

由分析结果可以看出,火灾作用对主梁竖向挠度以及主缆应力的影响很小,主梁最大竖向位移增量仅为5 mm,主缆最大应力增量0.16 MPa.结合图9可以看出,即使在最严重的场景D作用下,火灾主要造成直接受到辐射和对流影响的部分区段产生应力重分布,而对对流辐射区域外的区段应力分布影响很小.本桥主跨1 080 m,垂跨比1/9,仅考虑恒载作用下,由塔顶位移计算结果可知主缆最大竖向和横桥向变形率分别为1/26 087和1/20 074.因此仅考虑恒载作用下,火灾对主梁竖向挠度和主缆应力影响很小是合乎逻辑的.

图12 不同火灾场景下主缆应力增量图

4 考虑中塔火灾影响的管养策略

根据以上分析结果,针对中塔位置所发生的车辆火灾,本文给出以下管养策略:

1)从t2火源模型可以看出,火灾从发生到达到稳定燃烧要经过一定的时间(场景A,B,C和D达到稳定燃烧所需的时间分别为11 min,8 min,7 min和7 min),因此,车辆火灾发生后保证消防救援的及时性极为重要,为此必须设立专门的消防应急通道,建议消防救援到达时间不得超过15 min.

2)火灾场景A(普通小汽车火灾)对桥梁的影响可基本忽略不计,因此无需采取专门针对性的后续加固措施;火灾场景B,C和D均会对中塔造成较为明显的损伤,必须采取合理的加固措施.中塔火灾对主缆和加劲梁的影响很小,无需专门处理.

5 结 论

本文结合对桥梁火灾的已有认识,分析确定了4类桥梁火灾场景的基本参数,通过FDS得到了具有通用价值的火灾数值模拟结果,并与现有火灾模型曲线进行对比,验证了数值模拟的可靠性.在此基础上,以泰州长江大桥为工程实例,对中间钢塔火灾作用下的全桥静力性能进行了系统分析,研究得到的主要结论如下:

1)近塔桥面火灾作用对中塔影响十分明显,火灾作用区域由于材料性能受损,应力值相比正常状态下大幅降低;火灾作用下,受火区域横桥向变形明显呈现出外凸的趋势,中塔整体上呈现向远离火灾一侧倾斜的变形趋势,最大侧移超过50 mm.

2)近塔桥面火灾作用对加劲梁竖向挠度以及主缆应力的影响很小,基本可忽略不计.

3)结合火灾作用下中塔温度场分布以及全桥静力行为的分析结果,提出了具有针对性的防火管养策略,以供有关部门参考借鉴.

[1] MARIA G, IGNACIO P Z, VENKATESH K,etal. Fire hazard in bridges: review, assessment and repair strategies[J]. Engineering Structures, 2012 (35):89-98.

[2] BLINOV V I, KHUDYAKOV G. On certain regularities of diffusion burning of liquids[J]. Doklady Akademii Nauk SSSR, 1957(113) :241-244.

[3] MENDES P A, VALENTE J C, BRANCO F A. Simulation of ship fire under Vasco da Gama Bridge[J],ACI Struct J, 2000, 97 (2): 285-290.

[4] DOTREPPE J C, MAJKUT S, FRANSSEN J M. Failure of a tied-arch bridge submitted to a severe localized fire[C]//IABSE Symposium. Lisbon: Structures and Extreme Events, 2006:15-22.

[5] KODUR V R, GU L, GARLOCK M E. Review and assessment of fire hazard in bridges[C]//Journal of the Transportation Research Board No. 2172. Washington DC: Transportation Research Board of the National Academies,2010: 23-29.

[6] IGNACIO P Z, GARLOCK M E M. A 3D numerical analysis of a typical steel highway overpass bridge under a hydrocarbon fire[C]//Proceedings of the Fifth International Conference for Structures in Fire. East Lansing(USA) :2010:11-18.

[7] 程远平,陈亮,张孟君.火灾过程中火源热释放速率模型及其实验测试方法[J].火灾科学,2002, 11(2):70-74.

CHENG Yuanping, CHEN Liang, ZHANG Mengjun. The models and experimental testing method of heat release rate of fuel during the development of fire[J]. Fire Safety Science, 2002, 11(2):70-74. (In Chinese)

[8] LI G Q, WANG P J. Advanced analysis and design for fire safety of steel structures[M]. Hangzhou:Zhejiang University Press, 2012:15-17.

[9] HAACK A. Fire protection in traffic tunnels: general aspects and results of the eureka project [C] //Tunnellingand Underground Space Technology.Amsterdam: Elsevier Science Ltd, 1998: 377-381.

[10]CAFARO E, BERTOLA V. Fires in tunnels: experiments and modelling[J]. The Open Thermodynamics Journal, 2010(4) :156-166.

[11]倪照鹏,陈海云.国内外隧道防火技术现状及发展趋势[J]. 交通世界,2003,3(2):28-31.

NI Zhaopeng, CHEN Haiyun. Present situation and development trend for tunnel fireproof technigue at home and abroad[J]. Transpo World,2003,3(2):28-31.(In Chinese)

[12]LACROIX D. New French recommendations for fire ventilation in road tunnels[C]//Conf. on Aerodynamics and Ventilaion of Road Tunnels. Aosta(Italy):Mechanical Engineering Publications,1997:103-124.

[13]HAACK A.Thematic network on fires in tunnels-Technical report part1-design fire scenarios[M]. Brussels(Belgium) :WTCB,2001:61-72.

[14]HESELDEN A. Studies of fire and smoke behaviour relevant to tunnels[M]. East Acton :Building Research Establishment, Fire Research Station, 1978:66-78.

[15]LACROIX D. The new piarc report on fire and smoke control in road tunnel[C]//Third International Conference on Safety in Road and Rail Tunnels. Nice: Comité Technique 5 Tunnels Routiers, 1999:185-197.

[16]马明雷,马如进,陈艾荣. 桥面火灾条件下斜拉桥拉索及全桥结构安全性能研究[J]. 华南理工大学学报:自然科学版,2014,42(10):117-124.

MA Minglei, MA Rujin, CHEN Airong. Safety of cables and full structure of a cable-stayed bridge exposed to fires on deck[J]. Journal of South China University of Technology:Natural Science Edition, 2014, 42(10): 117-124.(In Chinese)

[17]李念平,肖书博,关军. 风管中聚集物燃烧时极限温度的分析研究[J]. 湖南大学学报:自然科学版,2009,36(6): 18-21.

LI Nianping, XIAO Shubo, GUAN Jun. Analysis and study of the limit temperature of the burning sediment in ventilation ducts[J]. Journal of Hunan University:Natural Sciences, 2009, 36(6): 18-21.(In Chinese)

[18]ISO 834-11:2014 Fire-resiestance tests elements of building construction, amendent 1, amendment 2[S]. Geneva: ISO, 2014:11-23.

[19]BABRAUSKAS V, ASCE manuals and reports on engineering practice no. 78: structural fire protection[C]// LIE T T. ASCE Manuals and Reports on Engineering Practice. New York: American Society of Civil Engineers, 1992: 96-122.

[20]OLST D V, BOSCH R V D. Behaviour of electrical cables during fire in tunnels[C]// SAVEUR J. Proceedings of ITA World Tunnel Congress 2003. Amsterdam: ITA, 2003: 991-992.

[21]汪剑,方志. 混凝土箱梁桥的温度场分析[J]. 湖南大学学报:自然科学版,2008,35(4): 23-28.

WANG Jian, FANG Zhi. Temperature variation of concrete box girder bridge[J]. Journal of Hunan University:Natural Sciences, 2008, 35(4): 23-28.(In Chinese)

[22]ENV 1993-1-2 Eurocode 3, Design of steel structures, part 1.2: structural fire design[S]. Brussels: European Committee for Standardization, 2005:18-23.

[23]ENV 1994-1-2 Eurocode 4, Design of composite and concrete structures, part 1.2: structural fire design[S]. Brussels: European Committee for Standardization, 2005:29-32.

[24]AISC 360-05 Specification for steel building [S]. Chicago: AISC-American Institute of Steel Construction Inc,2005:238-241.

[25]BS 5950-8:2003 The structural use of steelwork in buildings, part 8: code of practice for fire resistant design [S]. British: British Standards Institution (BSI), 2003:3-13.

Structural Performance Analysis of a Three-pylons Suspension Bridge Considering Fire Accident in the Vicinity of Middle Steel Pylon

MA Rujin†, CUI Chuanjie, MA Minglei

(Department of Bridge Engineering, Tongji University, Shanghai 200092, China )

The yield strength of the steel component gradually decreases under high temperature. As a result, the structural performance of the steel component on the bridges will be greatly affected by the on-deck fire disaster. Taking the steel pylon of a three-pylon suspension bridge as an example, this paper mainly studied the influence of the fire disaster on the structural performance of the steel pylon and the entire bridge. First of all, the temperature distribution inside the fire and steel pylon was obtained by the numerical simulation of different on-deck fire scenarios near the middle steel pylon. Secondly, the change of the structural static performance of the bridge under different fire scenarios was obtained by a nonlinear analysis. The results showed that under the effect of a serious vehicle fire, larger than 80 m2area of the middle steel pylon was under the temperature up to 800 ℃, and the maximum temperature would exceed 1 000 ℃, which can greatly affect the structural performance of the steel pylon. Obvious changes of stress and formation of the middle steel pylon were discorverd, which means 4.6 mm vertical residual deformation, 53.8 mm lateral residual deformation and 70 MPa stress reduction. However, the stress and deformation of the main cable and girder only showed minor changes.

structural analysis; three-pylons suspension bridge; middle steel pylon; fire accident; management and maintenance

1674-2974(2017)05-0088-08

10.16339/j.cnki.hdxbzkb.2017.05.011

2016-03-19 基金项目:“十二五”交通运输部重大科技专项(2011318494160),Major Science and Technology Projects of Ministry of Transport(2011318494160) 作者简介:马如进(1978-),男,江苏南通人,同济大学 副教授,博士 †通讯联系人,E-mail:rjma@tongji.edu.cn

马如进†,崔传杰,马明雷

(同济大学桥梁工程系, 上海 200092)

U445.7

A

猜你喜欢

中塔火源温度场
双火源隧道火灾数值模拟
不同火源位置情况下的内天井结构建筑
中外汉语教师教学情况的比较与思考
铝合金加筋板焊接温度场和残余应力数值模拟
火源位置对轻型门式刚架竖向位移的影响
吉林省主要森林火灾火源的时间变化特征
2219铝合金激光电弧复合焊接及其温度场的模拟
卡拉苏口岸:中塔贸易风生水起
MJS工法与冻结法结合加固区温度场研究
中塔企业成功签约