摘要
为了提高结构健康监测的经济性,扩大其应用范围,提出了一种基于气象站大数据信息的大跨度桥梁结构温度场预测方法,从而实现轻量化和可持续性的桥梁温度场及温度变形计算.以某大跨钢拱桁架桥为研究对象,通过气象数据平台获取气象站大数据信息(包括天气、气温和风速信息),利用Elbadry等改进的凯尔别克辐射模型计算桥梁关键构件各表面受到的辐射强度,基于热交换传递理论计算结构热分析边界条件,结合有限元分析方法进行结构瞬态热分析,从而获得目标桥梁关键构件的时变温度场.温度场的计算结果与桥梁上实测温度在时间历程和结构竖向分布上均吻合良好,平均误差在3%以内.将计算得到的结构温度施加到桥梁结构整体有限元模型上,得到桥梁结构的温度变形,桥梁支座纵向位移的计算值与实测值吻合良好,平均误差在13 mm以内.验证了所提出的基于气象站大数据的桥梁结构温度场预测方法是可行的,利用该方法可进一步得到桥梁结构的温度变形.
环境温度场的变化将引起桥梁结构温度场的变化. 温度作用使桥梁结构发生热胀冷缩,复杂结构的连接和边界条件的约束导致结构变形-内力相互制约,温度作用或对桥梁结构产生不利影响. 日照辐射导致的温差效应以及异常温度影响耦合其他荷载和已有损伤等,导致桥梁结构严重变形、开裂甚至发生倒塌事
自20世纪60年代起,国内外学者针对桥梁结构温度效应进行了大量现场观测和理论研究,对桥梁结构温度场的认识从线性到非线
处于人工智能和大数据发展的时代,有效利用大数据成为各行业发展的新方向. 来自各地面气象站的气象大数据具有易获得、可追溯和全覆盖等特点. 近十年来,各气象数据平台的建设和运行极大地方便了社会公众查询和获取历史气象数据. 例如国家气象科学数据中心、NOAA全球气象数据平台(National Centers for Environmental Information )以及Weather Spark气象数据平台等,可以用于查阅并下载国内外的地面气象站和各民用机场气象站的历史气象数据. 其涉及的气象站点遍布全球,可以覆盖国内外绝大多数桥梁. 通过气象数据平台可获得的气象数据包括天气情况、气温、风速和部分降水量等. 较为完备的、可追溯的气象站大数据信息和较为成熟的结构温度场数值分析理论为桥梁结构温度场的快速计算和温度效应的评估预测提供了可能.
本研究通过充分挖掘和利用气象站大数据信息,结合热交换传递理论和有限元分析方法获得桥梁结构关键构件的温度场,将模拟计算得到的结构温度场输入桥梁结构整体的有限元模型中,进而计算得到桥梁结构的温度变形. 为实现结构健康监测轻量化的目标,本文提出了一种基于气象站大数据计算大跨度桥梁结构时变温度场的方法,该方法可进一步用于桥梁结构温度效应分析.
1 结构温度场计算理论基础
桥梁结构时刻以热对流和热辐射的方式与外部环境进行热量交换,如

图1 桥梁结构热量交换示意
Fig.1 Schematic of heat exchange in bridges
1.1 气象参数的确定
计算结构温度场所需的气象参数有气温、风速和辐射强度,气温和风速数据可从气象大数据平台获取,辐射强度则可通过计算得到.
太阳照射角度在一天中的变化造成结构表面的光照面和阴影面的变化,日照情况的差异可量化为辐射强度值.若已知时间、地点、结构材料特性及目标结构各表面法向角度,即可利用辐射计算模型实现结构各表面的辐射强度计算.桥梁结构关键构件受到的辐射如

图2 桥梁关键构件辐射计算
Fig.2 Radiation calculation of key component on bridge
当前常用的辐射计算模型有凯尔别克模型、Hottel模型和ASHRAE模型. ASHRAE模型和Hottel模型以美国地区的长期太阳辐射实测数据为依据,通过经验公式计算太阳辐射强度;凯尔别克模
(1) |
(2) |
(3) |
(4) |
(5) |
(6) |
式中:I为总辐射强度;、、、、分别为太阳直接辐射强度、太阳散射辐射强度、地面反射辐射强度、大气逆辐射和环境辐射强度;为太阳常数;为林克氏浑浊度系数;为大气相对气压;为太阳高度角;为地面反射系数;为所求面与水平面的夹角;、分别为大气逆辐射和环境辐射系数;为黑体辐射系数;用于换算摄氏度和绝对温度,取273.15 K;为气温.
式(1)~
1.2 综合换热边界条件
关于对流换热系数,不同学者提出了多种经验公式,本研究中采用应用广泛的凯尔别克公
(7) |
(8) |
式中:、分别为对流换热系数和辐射换热系数;v为风速;ε为结构材料的发射率.
对于只存在对流换热的结构,结构表面的热流密度可通过第三类换热边界条
(9) |
式中:n为结构与周围流体的换热边界s上所求点的外法线方向;为结构表面温度;和分别为结构周围流体的温度及换热系数.
基于
(10) |
(11) |
式中:h为综合换热系数;T为综合气温.
以上式(7)~
1.3 结构温度场计算方法
基于1.1节、1.2节,将基于气象大数据的桥梁结构温度场计算方法总结如下.
1)气象数据获取:从气象大数据平台获取目标桥梁附近气象站的气温、风速数据及天气情况.
2)辐射强度计算:将气温数据、桥址经纬度、时间、桥梁结构材料参数(包括辐射率和辐射吸收系数)输入辐射强度计算模型[式(1)~
3)换热系数计算:将所得气温和风速代入
4)换热边界条件计算:将换热系数以及气象参数代入
5)桥梁结构温度场计算:将结构材料特性(包括热导率和比热容)和换热边界条件输入建立的结构有限元模型中,计算结构温度场.
结构温度场计算方法如

图3 基于气象大数据的桥梁结构温度场计算方法
Fig.3 Temperature field calculation method of the bridge based on meteorological big data
2 工程算例
2.1 工程背景
本节以某高速铁路大桥为研究对象,该桥的相关信息来自文献[

(a) 桥梁几何外形

(b) 温度测点布置图

(c) 关键构件截面概况
图4 目标桥梁概况(单位:mm)
Fig.4 Overview of the target bridge (unit: mm)
2.2 桥址气象参数计算
本研究采用Weather Spark气象数据平台(网址:http://zh.weatherspark.com)提供的南京某机场的逐小时天气、气温和风速数据. Weather Spark气象数据平台基于为游客提供出行参考的初衷,提供全世界 145 449个地点详细的天气报告,气象数据来自各地机场或地面气象站. 在通过各气象数据平台能够获取到逐小时历史气象信息的气象站中,该机场气象站是距离目标桥梁最近的气象站. 该机场与目标桥梁的直线距离为16 km,假定两地气温相同. 值得指出的是,有研究表明两地风速可能存在一定差异,故风速数据仅供参考,不直接用于计算. 因气温是在时间和空间上的连续变量,而风速的时变性较强,且易受局部温差、下垫面条件、地形地势和构筑物等影响,目标桥梁所在地区为丘陵地形,地势起伏较大,故取该机场气象站测得气温以及平均风速对结构温度场进行计算. 本研究利用2015年3月5—16日的气象数据,如

图5 气温数据
Fig.5 Temperature data

图6 风速数据
Fig.6 Wind speed data
算例中的太阳-桥梁相对位置如

图7 3月份太阳-桥梁相对位置
Fig.7 The sun-bridge relative position in March

(a) 8时

(b) 12时

(c) 16时
图8 3月13日关键构件日照情况
Fig.8 Sun exposure of the key component on March 13
根据基于气象大数据的桥梁结构温度场计算方法第2步,利用辐射强度计算模型[式(1)~

(a) 太阳直接辐射强度

(b) 总辐射强度
图9 3月13日关键构件辐射强度分布情况
Fig.9 Distribution of radiation intensity of the key components on March 13
2.3 有限元模型建立
在ANSYS中选用PLANE55单元(

图10 PLANE55单元
Fig.10 Unit PLANE55

(a) 截面A

(b) 截面B

(c) 截面C

(d) 截面D
图11 关键构件有限元模型
Fig.11 Finite element model of key components
根据基于气象大数据的桥梁结构温度场计算方法第3和第4步,将计算得到的各气象参数依次代入
3 结构温度场模拟结果分析
3.1 结构温度场模拟结果
2015年3月5—16日气象站测量的风速在0~ 7 m/s内变化,平均风速为3.5 m/s. 由于缺少现场实际风速数据,计算时取气象站测得的平均风速.
通过基于气象大数据的桥梁结构温度场计算方法第5步,得到各温度测点的温度时程曲线如

(a) WA-S测点

(b) WB-N测点

(c) WB-S测点

(d) WC-N测点

(e) WC-S测点

(f) WD-N测点

(g) WD-S测点
图12 温度时程曲线
Fig.12 Temperature time history curves
由于3月13日的平均温度最高,对当天的结构温度场模拟结果准确性进行进一步验证. 将各计算截面上两个温度测点的平均温度记为该截面的平均温度,绘制桥梁结构的典型竖向温度梯度如

(a) 日出时刻到7:00
(b) 日温度梯度最大时刻到15:00

(c) 日落时刻到18:00
(d) 夜间到22:00
图13 3月13日桥梁典型竖向温度梯度
Fig.13 Typical vertical temperature gradient of the bridge on March 13 (set the center height of Section D to 0 m)
(将截面D中心高度定为0 m)
综上,本节从各测点的温度时程曲线和桥梁结构典型竖向温度梯度两个方面,对结构温度计算值在时间历程和结构竖向上分布的准确性进行了验证. 对比分析结果显示应用所提出的基于气象大数据的结构温度计算方法进行桥梁结构温度场计算是可行的.
3.2 温度场误差分析
对各温度测点的计算-实测值进行误差分析,如

(a) WA-S测点
(b) WB-N测点

(c) WB-S测点
(d) WC-N测点

(e) WC-S测点
(f) WD-N测点

(g) WD-S测点
图14 模拟温度值误差分析
Fig.14 Analog temperature value error analysis
其中WA-S和WC-S测点的误差基本控制在15%以内;WC-N和WD-N测点的误差基本在5%以内;截面B两个测点的温度计算值误差均较大. 截面B与算例中其他截面形式不同,计算其温度场时未考虑其自身对风的阻挡作用,或为截面B温度计算值误差过大的原因.
为探究各气象参数取值偏差对结构模拟温度值的影响,利用ANSYS计算目标桥梁上关键构件A在不同气象条件下的平均温度,对关键构件进行了气象参数灵敏度分析(

图15 气象参数灵敏度分析
Fig.15 Sensitivity analysis of meteorological parameters
辐射强度和风速的取值极大影响着结构温度场计算结果. 为了提高基于气象大数据的桥梁结构温度计算方法的计算精度,需要明确桥梁结构上各位置处的风速,同时需要进一步研究不同天气情况下的辐射强度计算方法.
4 桥梁结构温度变形计算
4.1 支座纵向位移计算
大跨桥梁结构的纵向位移主要由温度变化引起.在ANSYS中建立桥梁的有限元分析模型(

图16 桥梁结构有限元模型
Fig.16 Finite element model of the bridge
模型的边界条件设置为:中间支座处施加纵向和竖向约束,中心节点再多施加横向约束,其余支座处施加竖向约束.目标桥梁结构采用的主要材料为钢材,钢材型号包括Q420qE、Q370qE和Q345qD,其热膨胀系数均取为1.2×1
将计算所得截面B的平均温度施加于有限元模型中桥面板底面以及桥面板下纵横梁节点,对应桥梁结构上始终无法受到日照的部位. 将上部钢拱圈中截面A、B、C三者的平均温度施加于其他节点. 将整个结构的初始温度设置为桥梁支座位移传感器安装后投入使用时的综合气温. 桥梁各支座纵向位移的计算值(ANSYS)和实测值(SHM)如

(a) 1#支座

(b) 2#支座

(c) 3#支座

(d) 6#支座

(e) 7#支座
图17 桥梁支座纵向位移
Fig.17 Longitudinal displacement of bridge bearings
考虑到以传感器投入使用时刻为零点计算桥梁结构温度变形,其跨越的时长较长,或导致计算结果误差增加,故对3月13日一天中各支座的相对增量位移值进行了验证,如

(a) 1#支座
(b) 2#支座

(c) 3#支座
(d) 6#支座

(e) 7#支座
图18 桥梁支座纵向相对增量位移
Fig.18 Longitudinal relative incremental displacement of
bridge bearings
4.2 支座位移误差分析
对各支座纵向位移的计算-实测值进行误差计算,如
纵向位移误差 | 最大误差/mm | 平均误差/mm |
---|---|---|
1# | 16 | 0 |
2# | 22 | 4 |
3# | -21 | -13 |
6# | -17 | -4 |
7# | -29 | -13 |
除7#支座外,其余支座纵向位移的最大误差均在22 mm内,平均误差均在13 mm内. 该桥梁结构以4#支座为中心两边对称布置,其两侧支座的纵向位移应关于主跨固定支座呈反对称分布,但7#支座的实测位移值远小于其对称位置的1#支座. 或存在其他因素影响了该支座的纵向位移.
4.3 支座绝对位移计算
实际上,桥梁自合龙时刻起各支座开始发生位移,该桥梁健康监测系统测得的支座位移值是以传感器投入使用时刻为零点的相对位移值. 将整个结构的初始温度设置为该桥合龙时(2009年9月28日上午)的综合气温,计算得到桥梁各支座的绝对纵向位移如

图19 支座绝对纵向位移
Fig.19 Absolute longitudinal displacement of the bearings
实测1#支座相对位移的最大值为87 mm,计算得到1#支座的绝对位移最大值为140 mm,远大于实测值(

图20 1#支座绝对-相对位移
Fig.20 Absolute-relative displacement of 1# bearing
5 结 论
本文旨在探究直接利用气象大数据分析桥梁结构时变温度场以及温度变形的方法,以实现轻量化和可持续性的桥梁结构健康监测. 通过研究分析,得到以下主要结论:
1)提出了一种基于气象站大数据的桥梁结构温度场计算方法. 将温度场计算结果与桥梁上实测温度进行对比,两者在时间历程和结构竖向分布上均吻合良好,平均误差绝对值在3%以内.
2)提出了结构温度场数值模拟时风速的取值方法. 缺少桥址处实测风速数据时,可参考桥址附近气象站测得的风速,取平均风速计算. 计算夜间结构温度场时视情况可适当降低风速.
3)分析了各气象参数变化对结构温度的影响. 辐射强度对结构温度的影响大于气温,风速与结构温度成反比例函数关系. 风速低于3 m/s时,其对结构温度的影响甚至大于辐射强度.
4)基于计算温度场得到桥梁的温度变形. 各支座纵向位移计算值与实测值吻合良好,平均误差在13 mm内. 以桥梁合龙时刻为零点,计算得到各支座的绝对位移远大于实测的相对位移.
参考文献
ZHU S J,LEVINSON D,LIU H X,et al.The traffic and behavioral effects of the I-35W Mississippi River bridge collapse[J].Transportation Research Part A:Policy and Practice,2010,44(10):771-784. [百度学术]
KROMANIS R,KRIPAKARAN P,HARVEY B. Long-term structural health monitoring of the Cleddau bridge:evaluation of quasi-static temperature effects on bearing movements[J].Structure and Infrastructure Engineering,2016,12(10):1342-1355. [百度学术]
HUANG H B,YI T H,LI H N,et al.New representative temperature for performance alarming of bridge expansion joints through temperature-displacement relationship[J].Journal of Bridge Engineering,2018,23(7):04018043. [百度学术]
LIU H B,WANG X Q,JIAO Y B.Effect of temperature variation on modal frequency of reinforced concrete slab and beam in cold regions[J].Shock and Vibration,2016,2016:4792786. [百度学术]
闵志华,孙利民,淡丹辉.影响斜拉桥模态参数变化的环境因素分析[J].振动与冲击,2009,28(10):99-105. [百度学术]
MIN Z H,SUN L M,DAN D H. Effect analysis of environmental factors on structural modal parameters of a cable-stayed bridge[J].Journal of Vibration and Shock,2009,28(10):99-105.(in Chinese) [百度学术]
ZUK W. Thermal behaviour of composite bridges-insulated and uninsulated[J].Highway Research Record,1965,76(1):231-253. [百度学术]
PRIESTLEY M N.Design of concrete bridges for temperature gradients[J].ACI Journal Proceedings,1978,75(5):209-217. [百度学术]
ELBADRY M M,GHALI A.Temperature variations in concrete bridges[J].Journal of Structural Engineering,1983,109(10):2355-2374. [百度学术]
DILGER W H,GHALI A,CHAN M,et al.Temperature stresses in composite box girder bridges[J].Journal of Structural Engineering,1983,109(6):1460-1478. [百度学术]
葛耀君,翟东,张国泉.混凝土斜拉桥温度场的试验研究[J].中国公路学报,1996,9(2):76-83. [百度学术]
GE Y J,ZHAI D,ZHANG G Q. Temperature gradient test and research of RC cable-stayed bridges[J]. China Journal of Highway and Transport,1996,9(2):76-83.(in Chinese) [百度学术]
张元海, 李乔.斜支承箱梁的日照温度次内力研究[J].土木工程学报,2008,41(5):46-51. [百度学术]
ZHANG Y H,LI Q. Study on the secondary internal forces in skew box girders induced by sunshine temperature gradient[J].China Civil Engineering Journal,2008,41(5):46-51.(in Chinese) [百度学术]
叶见曙,贾琳,钱培舒.混凝土箱梁温度分布观测与研究[J].东南大学学报(自然科学版),2002, 32(5):788-793. [百度学术]
YE J S,JIA L,QIAN P S.Observation and research on temperature distribution in concrete box girders[J].Journal of Southeast University (Natural Science Edition),2002,32(5):788-793.(in Chinese) [百度学术]
王达, 张永健, 刘扬, 等.基于健康监测的钢桁加劲梁钢-混组合桥面系竖向温度梯度效应分析[J].中国公路学报,2015, 28(11): 29-36. [百度学术]
WANG D,ZHANG Y J,LIU Y,et al.Vertical temperature gradient effect analysis of steel-concrete composite deck system on steel truss stiffening girder with health monitoring[J]. China Journal of Highway and Transport,2015,28(11):29-36.(in Chinese) [百度学术]
丁幼亮, 王高新, 周广东, 等.基于现场监测数据的润扬大桥斜拉桥钢箱梁温度场全寿命模拟方法[J].土木工程学报,2013, 46(5): 129-136. [百度学术]
DING Y L,WANG G X,ZHOU G D,et al.Life-cycle simulation method of temperature field of steel box girder for Runyang cable-stayed bridge based on field monitoring data[J]. China Civil Engineering Journal,2013,46(5):129-136.(in Chinese) [百度学术]
XU Y L,CHEN B,NG C L,et al.Monitoring temperature effect on a long suspension bridge[J].Structural Control and Health Monitoring,2010,17(6):632-653. [百度学术]
刘扬,张海萍,邓扬,等.基于实测数据的悬索桥钢箱梁温度场特性研究[J].中国公路学报,2017,30(3):56-64. [百度学术]
LIU Y,ZHANG H P,DENG Y,et al.Temperature field characteristic research of steel box girder for suspension bridge based on measured data[J].China Journal of Highway and Transport,2017,30(3):56-64.(in Chinese) [百度学术]
LU Y,LI D J,WANG K,et al.Study on solar radiation and the extreme thermal effect on concrete box girder bridges[J].Applied Sciences,2021,11(14):6332. [百度学术]
赵人达,王永宝.日照作用下混凝土箱梁温度场边界条件研究[J].中国公路学报,2016,29(7):52-61. [百度学术]
ZHAO R D,WANG Y B.Studies on temperature field boundary conditions for concrete box-girder bridges under solar radiation[J].China Journal of Highway and Transport,2016,29(7): 52-61.(in Chinese) [百度学术]
ABID S R,TAYŞI N,ÖZAKÇA M.Experimental analysis of temperature gradients in concrete box-girders[J].Construction and Building Materials,2016,106:523-532. [百度学术]
LU H L,HAO J,ZHONG J W,et al.Analysis of sunshine temperature field of steel box girder based on monitoring data[J].Advances Civil Engineering,2020,2020:9572602. [百度学术]
ZHOU L R,LIANG C F,CHEN L,et al.Numerical simulation method of thermal analysis for bridges without using field measurements[J].Procedia Engineering,2017,210:240-245. [百度学术]
WANG Z W,ZHANG W M,ZHANG Y F,et al.Temperature prediction of flat steel box girders of long-span bridges utilizing in situ environmental parameters and machine learning[J].Journal of Bridge Engineering,2022,27(3):04022004. [百度学术]
FAN J S,LI B L,LIU C,et al.An efficient model for simulation of temperature field of steel-concrete composite beam bridges[J].Structures,2022,43:1868-1880. [百度学术]
凯尔别克. 太阳辐射对桥梁结构的影响[M]. 刘兴法, 译. 北京: 中国铁道出版社,1981:189-193. [百度学术]
KEHLBECK F.Effects of solar radiation on bridge structures[M]. Beijing: China Railway Press,1981:189-193.(in Chinese) [百度学术]
BRANCO F A,MENDES P A.Thermal actions for concrete bridge design[J].Journal of Structural Engineering,1993, 119(8): 2313-2331. [百度学术]
张兴中,黄文,刘庆国.传热学[M].北京:国防工业出版社,2011:55-62. [百度学术]
ZHANG X Z,HUANG W,LIU Q G.Heat transmission science[M].Beijing:National Defense Industry Press,2011:55-62.(in Chinese) [百度学术]
WANG G X,DING Y L,SONG Y S,et al.Detection and location of the degraded bearings based on monitoring the longitudinal expansion performance of the main girder of the Dashengguan Yangtze bridge[J].Journal of Performance of Constructed Facilities,2016,30(4):04015074. [百度学术]
WANG G X,DING Y L,SUN P,et al.Assessing static performance of the Dashengguan Yangtze bridge by monitoring the correlation between temperature field and its static strains[J].Mathematical Problems in Engineering,2015,2015:946907. [百度学术]
HUANG Q H,CROSETTO M,MONSERRAT O,et al.Displacement monitoring and modelling of a high-speed railway bridge using C-band Sentinel-1 data[J].ISPRS Journal of Photogrammetry and Remote Sensing,2017,128:204-211. [百度学术]
XIA C Y,ZHONG T Y.Numerical analysis of the Nanjing Dashengguan Yangtze River Bridge subjected to non-uniform seismic excitations[J].Journal of Mechanical Science and Technology,2011,25(5):1297-1306. [百度学术]