在当前的城市箱梁桥中,混凝土箱梁因受力合理、造价较低而占据较大的比例[1]。日照环境下,由于混凝土材料的导热性能较差,混凝土箱梁内部产生非线性的竖向温度梯度[2]。在这种温度效应下,箱梁局部会产生温度应力[3],而温度应力的累积效应会引起结构开裂[4-5],从而影响到桥梁结构的安全性[6]。
国内外学者对混凝土箱形梁的竖向温度梯度进行了大量的研究和讨论[6-8]。Lee等[7]在实验混凝土箱梁段中,记录并利用了整个一年周期的混凝土温度、太阳辐射、空气温度和风速数据,研究不同环境条件下的垂直和横向温度幅度。杨宏印等[8]指出桥梁结构温度变化主要依赖其所处的环境温度,主梁在纵桥方向上的不同位置具有一致的温度分布形式;卫俊岭等[9]用ANSYS软件建立了混凝土箱梁桥二维瞬态日照温度场模型,进一步模拟了混凝土箱梁桥的最大竖向温度梯度分布特征,分析了该温度分布模式对桥梁的作用效应。然而,研究成果大多是基于温度监测数据得到桥梁的竖向温度梯度,仅适用于桥址处所在的地区。现有桥梁温度场的仿真分析方法在热边界条件的处理上过于粗略[10],无法客观地考虑空气流动、辐射的多次反射、日照遮挡等因素。我国跨越多个气候区,但桥梁结构的温度梯度模式统一采用《公路桥涵通用规范》(JTG D60—2015)[11],该规范中的温度梯度模式是否适用于全国各地区的桥梁结构有待更深入地研究。
本文以某混凝土箱梁桥为研究对象,对箱梁典型截面的温度进行长期的监测和分析,以有限元软件联合使用的仿真分析方法建立其温度场时空模型,将模型计算值与现场监测值进行对比。在此基础上,研究极端天气条件下混凝土箱梁的最不利竖向温度梯度分布模式,分析各主要气象因素对混凝土箱梁温度梯度的影响。此外,文中的研究方法可为其他地区桥梁结构的温度场研究提供借鉴和参考。
1 工程背景
本文针对武汉市某立交匝道桥的第三联(20 m+20 m+20 m)混凝土连续曲线箱梁桥,桥梁跨径布置如图1所示。上部结构采用单箱单室截面,梁高1.50 m,顶板宽11.00 m,底板宽6.00 m,翼缘板均长2.50 m,翼缘板高度由0.15 m至0.40 m呈线性变化。桥面铺装由6 cm厚混凝土整平层和7 cm厚沥青混凝土组成。
<G:\武汉工程大学\2025\第2期\彭娅-1.tif>[2 000][2 000][2 000][温度传感器(29)][YC6][YC7][YC8][YC9][单位:cm]
图1 桥梁跨径布置图
Fig. 1 Bridge span layout
考虑到经济、效率,选择混凝土箱梁桥的中跨跨中截面合理布置温度传感器,共计29个,各温度监测截面尺寸及测点布置如图2所示,采用实时监测的方法,每10 min采集1次混凝土箱梁截面腹板、顶板及底板的温度数据。
2 混凝土箱梁桥温度场模型及验证
2.1 温度场时空模型建立
应用ANSYS、HYPERMESH和TAITHERM 3个有限元软件联合使用的仿真分析方法,建立了基于桥梁的基本参数的温度场时空模型。首先在ANSYS软件中建立有限元模型,然后导入HYPERMESH中对其进行网格划分、格式转换,最后再导入TAITHERM软件中,对桥梁结构的温度场进行仿真计算。温度场仿真计算流程如图3所示。
<G:\武汉工程大学\2025\第2期\彭娅-3.tif>[混凝土箱梁温度场][ANSYS
建立模型][HYPERMESH网格划分][TAITHERM温度场计算][否][温度监测值][温度场
仿真分析][是][仿真值与监测值是否吻合?]
图3 温度场仿真计算流程图
Fig. 3 Flow chart of temperature field simulation calculation
TAITHERM作为专业的热分析软件,可以综合模拟热传导、热对流和热辐射3种传热模式,具有高效的热建模和热求解能力。将含太阳辐射、气温、云量、风速等实时气象数据的天气文件作为热边界条件导入TAITHERM软件中,可准确地模拟环境因素对桥梁结构温度场的影响。在TAITHERM软件中建立的混凝土箱梁温度场时空模型如图4所示。
<G:\武汉工程大学\2025\第2期\彭娅-4.tif>[Model Size (mr):
X=33 191.5
Y = 56 896.9
Z= 1 290
Visible Size (mm):
X=33 191.5
Y = 56 896.9
Z =1 290
Visible Counts:
Parts=8
Elements=12 614
][X][Y][O][Z]
图4 混凝土箱梁温度场时空模型
Fig. 4 Spatio-temporal model of temperature field of
concrete box girder
该模型分为8个部分,总计12 614个单元,每个部分沿板厚方向均划分为6层,各部分可根据实际材料及厚度设定,材料的热物理参数如表1所示。在TAITHERM软件中将模型的地理位置、海拔等参数设置成武汉市的特征值,将WolframAlpha网站采集的气象数据作为天气文件导入软件中进行计算,得到混凝土箱梁桥各时刻的温度场分布。根据温度监测数据,1 h内桥梁温度变化不大,故计算模型的时间步长设为1 h。
2.2 温度场时空模型验证
为提高温度场时空模型的精度,基于温度监测数据对模型进行修正,进而得到混凝土箱梁温度场的精细化时空模型。将混凝土箱梁温度的计算值与监测值进行比较分析,从而验证温度场时空模型的准确性。根据桥梁实际的测点布置情况,在模型中提取对应单元的温度仿真计算数据,限于篇幅,综合考虑监测数据的完整性和天气状况,仅展示2020年12月典型测点(顶板T27、底板T11)的计算值和监测值的对比,如图5所示。
由图5可以看出,混凝土箱梁顶、底板测量点的温度计算值与监测值非常接近,变化趋势也很吻合。绝对误差均在3 ℃以内,满足实际工程的精度要求。因此,本文采用的温度场仿真计算方法与选取材料热物理参数合理准确。
3 混凝土箱梁竖向温度梯度分析
3.1 竖向温度梯度影响因素分析
选取混凝土箱梁中同一截面温度最低的数据作为基准温度,竖向温度梯度[12-13]由其余各点的温度与基准温度之差得到。通过分析各气象因素对混凝土箱梁竖向温度梯度的影响大小,可以更深入地研究日照下混凝土箱梁的竖向温度梯度。主要影响因素包括桥梁结构因素和气象因素[14]。然而,实际工程设计中一般以规范规定的温度梯度作为控制作用,在桥梁结构设计参数确定的情况下,其温度梯度已确定,太阳辐射、气温及风速等主要气象因素就成为影响桥梁结构温度场分布的变量。
为研究各主要气象因素对混凝土箱梁竖向温度梯度的影响,基于上述的温度场时空模型,选择2020年8月11日-20日天气较晴朗的10 d,采用控制变量法将这10天的气象数据依次控制,各参数按20%的幅度增减,分别计算太阳辐射、气温和风速对于混凝土箱梁竖向温度梯度的影响,计算结果见图6。
由图6可知,在各主要气象因素中,太阳辐射对日照作用下混凝土箱梁竖向温度梯度的影响最大,气温次之,风速的影响最小。结果表明:太阳辐射、气温与竖向温度梯度呈显著正相关,风速与竖向温度梯度呈显著负相关,气象因素对混凝土箱梁竖向温度梯度的影响主要集中在顶板和箱梁底部。
对上述计算结果进一步分析,得到各气象因素对混凝土箱梁竖向温度梯度最大值的影响程度,分析结果见图7。由图7可知,当太阳辐射和气温增加(减少)40%时,混凝土箱梁竖向温度梯度的最大值分别增加(减少)40%和15%;当风速增加(减少)40%时,混凝土箱梁竖向温度梯度的最大值减少(增加)10%。
3.2 最不利竖向温度梯度分布模式
为保证桥梁结构在设计基准期限内能正常运营,有必要分析极端天气条件对桥梁结构的影响[15]。对于混凝土箱梁桥在设计基准期限内的最不利竖向温度梯度,可以通过构建极端天气数据,结合建立的温度场时空模型进行预测计算。
构建极端天气数据,需考虑太阳辐射、气温及风速等气象参数的随机性[16]。根据《工业建筑供暖通风与空气调节设计规范》[17]确定室外空气代表参数的方法,确定武汉市(北纬30°、3级大气透明度等级)最不利太阳辐射的逐时代表值。以武汉市历年夏季的日最高气温为研究对象,通过概率统计分析方法得到最不利气温的逐时代表值[18]。因风速变化规律不明显、离散性较大,为计算方便,取武汉市历年夏季的日平均风速作为最不利风速的逐时代表值。
基于最不利太阳辐射、气温和风速等气象参数构建武汉市极端天气资料[19]。将数据导入TAITHERM软件,计算混凝土箱梁的温度场分布,得到极端天气条件下混凝土箱梁最不利的垂直温度梯度分布模式。箱梁截面沿梁高方向的最大温差为13.4 ℃,底板与腹板之间的温差为1.0 ℃。最不利竖向温度梯度分布拟合公式为:
[ty=13.4e-6y, 0≤y≤(H-0.2)ty=5.2y-6.76, (H-0.2)≤y≤H] (1)
式中:[ty]为箱梁截面梁高方向计算点处的温差;[y]为计算点至箱梁顶面的距离;[H]为箱梁截面的高度。
关于城市混凝土箱梁竖向温度梯度的取值,国内现行主要依据《公路桥涵通用设计规范》(JTG D60—2015)[10],其桥面板表面的温度t1规定见表2。该桥为沥青混凝土铺装,铺装层厚度为70 mm,根据线性内插法可以得到t1=17.6 ℃,t2=6.22 ℃。
表2 竖向日照正温差计算的温度基数
Tab. 2 The temperature base for calculating the positive
vertical solar radiation temperature difference
[结构类型 t1 / ℃ t2 / ℃ 水泥混凝土铺装 25 6.7 50 mm沥青混凝土铺装层 20 6.7 100 mm沥青混凝土铺装层 14 5.5 ]
根据预测的温度梯度值,将最不利的垂直温度梯度分布模式与公路桥规范规定的温度梯度分布进行对比,结果如图8所示。由图8可知,温度梯度预测值在混凝土箱梁顶板厚度范围内时均符合规范,在规范范围内,按设计规范取值是安全的。但腹板上部和底板的厚度范围内仍存在明显的垂直温度梯度,但规范中并无相关规定,特别是对混凝土箱梁温度效应的详细分析,应引起注意。因此,规范中的双折线温度梯度模式并不能完全反映所有桥梁的实际温度分布,在混凝土箱梁温度场及其温度效应分析中,应尽可能结合桥梁实际地理、气象等情况具体分析,以使结果更精确。
<G:\武汉工程大学\2025\第2期\彭娅-8.tif>[-0.3
0.0
0.3
0.6
0.9
1.2
1.5
1.8][距主梁顶部的距离 / m][-3 0 3 6 9 12 15 18 21
温度 / ℃][桥面板表面][预测温度梯度值
预测温度梯度分布
规范温度梯度分布]
图8 竖向温度梯度分布模式对比
Fig. 8 Comparison of vertical temperature gradient
distribution patterns
4 结 论
本文以某混凝土箱梁桥为工程背景,结合桥址处的天气资料,采用有限元软件联合使用的方法建立其温度场时空模型,基于温度场模型,对影响混凝土箱梁竖向温度梯度的主要气象因素进行了分析,研究了极端天气条件下混凝土箱梁的最不利竖向温度梯度分布模式,并将其与现行规范进行了对比。得到如下结论:
(1)混凝土箱梁桥温度的模型计算值与监测值吻合良好,说明本文所采用的有限元软件联合使用的仿真分析方法可较准确地计算混凝土箱梁温度场,该方法亦适用于其他地区和其他桥型。
(2)太阳辐射对日照作用下混凝土箱梁竖向温度梯度的影响较大,气温次之,风速的影响较小,且太阳辐射和气温与其竖向温度梯度呈正相关,风速与其竖向温度梯度呈负相关。
(3)基于极端天气数据和温度场时空模型,可计算得到极端天气条件下混凝土箱梁桥的最不利竖向温度梯度分布模式,其与规范中全国统一的双折线温度梯度模式存在一定的差异,规范未考虑腹板上部及底板处的竖向温度梯度,需引起重视。在混凝土箱梁温度场及其温度效应分析中,应尽可能结合桥梁实际地理、气象等情况具体分析,以使结果更精确。