一种桥梁结构温度场监测方法技术

本发明专利技术公开了一种桥梁结构温度场监测方法，包括以下步骤：S1：利用实测气温和实时风速分别计算对流换热和辐射换热；S2：利用日序数和水平面上的太阳辐射计算太阳辐射；S3：根据对流换热、辐射换热和太阳辐射计算桥梁的热边界条件；S4：根据桥梁的热边界条件确定桥梁的热交换模型；S5：利用有限元数值模拟求解桥梁的热交换模型，得到桥梁的温度场，完成桥梁结构温度场的监测。本发明专利技术不需要过于依赖实测温度数据，只需要在桥址处安装包括温度传感器、风速传感器和总辐射表在内的气象站，初期使用时建议在桥梁上安装少数温度传感器，调整模型参数；成本低，安装使用方便；监测范围广，可以实现全桥温度监测。

【技术实现步骤摘要】
一种桥梁结构温度场监测方法
本专利技术属于桥梁结构
，具体涉及一种桥梁结构温度场监测方法。
技术介绍
在日照温差作用下，由于混凝土本身较差的导热性，外部温度迅速升高，内部温度基本不变，桥梁温度场呈现明显的非线性分布特征。对于混凝土高墩桥梁，空心高墩的非线性温差会引起桥梁变形，对桥梁的静力状态和动力特性产生不可忽视的影响。随着我国交通事业的发展，西南地区的路网建设不断完善，越来越多跨越山区峡谷的混凝土高墩桥梁开始建设运营。因此，专利技术一种该地区混凝土高墩桥梁的温度场监测方法显得尤为迫切。传统的桥梁温度场监测的方法要求在桥梁某些位置安装一定数量的温度传感器。这种方法存在三个缺点：第一，所需传感器和数据采集传输设备成本高，安装与调试都需要耗费大量人力物力；第二，只能获取桥梁有限位置的温度数据，无法把握全桥温度分布；第三，传感器寿命有限，表面传感器需要定期更换，预埋的传感器在几年之后则无法反馈数据。通过有限元模拟获取桥梁的温度分布，是桥梁日照温度场分析常用的方法。有限元模拟的难点在于准确计算以太阳辐射为主导的热边界条件，传统的计算方法是参考大气科学中的太阳辐射模型，如AHSRAE模型、Hottel模型和幂指数模型等。但是，理论计算模型中的系数需要依赖大量实测数据和相关的数据统计经验，如果用于模拟温度场，存在三点不足：第一，太阳辐射计算模型只适用于晴天，并不适用于所有天气情况，无法计算阴雨天的太阳辐射；第二，对于一座具体的桥梁，根据理论模型计算的辐射量，只能用于年尺度上的温度场模拟，具体到天尺度上，则无法考虑由于云层遮挡导致的地表太阳辐射变化，其精度尚有待提高；第三，随着时间的推移，我国工业化进程加快，以空气质量为代表的环境污染不可忽视，大气透明度发生了变化，三十多年前的经验系数用于现在，精度有所下降。因此，太阳辐射模型只可用于一般情况或极端情况下的温度场数值模拟，无法达到实时温度场监测的精度要求。如何准确计算桥址处太阳辐射，进一步确定桥梁各部分构件所处的热边界条件，是桥梁温度场精细化模拟的难点。
技术实现思路
本专利技术的目的是为了解决的问题，提出了一种桥梁结构温度场监测方法。本专利技术的技术方案是：一种桥梁结构温度场监测方法包括以下步骤：S1：利用实测气温和实时风速分别计算对流换热和辐射换热；S2：利用日序数和水平面上的太阳辐射计算太阳辐射；S3：根据对流换热、辐射换热和太阳辐射计算桥梁的热边界条件；S4：根据桥梁的热边界条件确定桥梁的热交换模型；S5：利用有限元数值模拟求解桥梁的热交换模型，得到桥梁的温度场，完成桥梁结构温度场的监测。进一步地，步骤S1中，对流换热qc的计算公式为：qc＝hc(Ta-T)其中，hc表示对流换热系数，Ta表示实测气温，T表示桥梁温度场，V表示实时风速；辐射换热qr的计算公式为：qr＝hr(Ta-T)-qrahr＝Cε[(T+273)2+(Ta+273)2]·[T+Ta+546]其中，hr表示辐射换热系数，Ta表示气温，T表示桥梁温度场，qra表示斜面天空辐射效应，C表示Stefan-Boltzmann常数，ε表示物体发射率。温度场数值模拟为瞬态分析，需要先假定整个桥梁的温度场T，作为桥梁温度状态的初始条件，以此来计算对流换热和辐射换热。一般采用日出前的环境温度作为桥梁温度场T的初始值，通过连续多天的计算，逐渐地消除初值选取带来的误差。进一步地，步骤S2包括以下子步骤：S21：根据日序数N计算太阳常数I0；S22：利用总辐射传感器测量水平面上的太阳辐射I；S23：根据太阳常数I0和水平面上的太阳辐射I计算大气清洁指数kT，根据水平面上的太阳辐射I计算地面发射辐射re×I，其中，re表示地面短波辐射反射率；S24：根据大气清洁指数kT和水平面上的太阳辐射I计算天空散射辐射IdH；S25：根据地面发射辐射re×I和天空散射辐射IdH计算投射到表面的太阳辐射量It；S26：根据投射到表面的太阳辐射量It计算太阳辐射qs。进一步地，步骤S21中，太阳常数I0的计算公式为：其中，N表示自每年1月1日起算的日序数，太阳常数I0的单位是W·m-2。进一步地，步骤S23中，大气清洁指数kT的计算公式为：其中，I表示水平面上的太阳辐射，I0表示太阳常数。进一步地，步骤S24中，天空散射辐射IdH的计算公式为：其中，kT表示大气清洁指数，I表示水平面上的太阳辐射。进一步地，步骤S25中，投射到表面的太阳辐射量It的计算公式为：其中，Ib表示太阳直接辐射；IdH表示天空散射辐射；I表示水平面上的太阳辐射，re表示地面短波辐射反射率，re×I表示地面反射辐射，φ表示太阳光线入射角，βn表示受照表面外法线与地平面的夹角。进一步地，步骤S26中，太阳辐射qs的计算公式为：qs＝αsIt其中，αs表示混凝土表面短波辐射吸收系数，It表示投射到表面的太阳辐射量。进一步地，步骤S3中，桥梁的热边界条件q的计算公式为：q＝qs+qc+qr其中，qs表示太阳辐射；qc表示对流换热；qr表示辐射换热。进一步地，步骤S4中，桥梁的热交换模型的计算公式为：其中，k表示混凝土材料热传导率，T表示桥梁温度场，x表示桥梁截面内坐标的横坐标，y表示桥梁截面内的纵坐标，nx表示边界表面横桥向单位法向量方向余弦，ny表示表示边界表面竖桥向单位法向量方向余弦，q表示桥梁的热边界条件。本专利技术的有益效果是：(1)本专利技术不需要过于依赖实测温度数据，只需要在桥址处安装包括温度传感器、风速传感器和总辐射表在内的气象站，初期使用时建议在桥梁上安装少数温度传感器，调整模型参数。(2)成本低，安装使用方便。(3)监测范围广，可以实现全桥温度监测。附图说明图1为桥梁结构温度场监测方法的流程图；图2为桥梁表面与外界环境的热交换的示意图；图3为太阳光线入射角的示意图；图4为受照表面外法线与地平面的夹角的示意图；图5为混凝土简支梁俯视图；图6为实施例的气象站示意图；图7为安装温度传感器的示意图；图8为主梁顶部的数据对比图；图9为主梁底部的数据对比图；图10为主梁朝南的数据对比图；图11为主梁朝北的数据对比图。具体实施方式下面结合附图对本专利技术的实施例作进一步的说明。如图1所示，本专利技术提供了一种桥梁结构温度场监测方法，包括以下步骤：S1：利用实测气温和实时风速分别计算对流换热和辐射换热；S2：利用日序数和水平面上的太阳辐射计算太阳辐射；S3：根据对流换热、辐射换热和太阳辐射计算桥梁的热边界条件；S4：根据桥梁的热边界条件确定桥梁的热交换模型本文档来自技高网...

【技术保护点】
1.一种桥梁结构温度场监测方法，其特征在于，包括以下步骤：/nS1：利用实测气温和实时风速分别计算对流换热和辐射换热；/nS2：利用日序数和水平面上的太阳辐射计算太阳辐射；/nS3：根据对流换热、辐射换热和太阳辐射计算桥梁的热边界条件；/nS4：根据桥梁的热边界条件确定桥梁的热交换模型；/nS5：利用有限元数值模拟求解桥梁的热交换模型，得到桥梁的温度场，完成桥梁结构温度场的监测。/n

【技术特征摘要】
1.一种桥梁结构温度场监测方法，其特征在于，包括以下步骤：
S1：利用实测气温和实时风速分别计算对流换热和辐射换热；
S2：利用日序数和水平面上的太阳辐射计算太阳辐射；
S3：根据对流换热、辐射换热和太阳辐射计算桥梁的热边界条件；
S4：根据桥梁的热边界条件确定桥梁的热交换模型；
S5：利用有限元数值模拟求解桥梁的热交换模型，得到桥梁的温度场，完成桥梁结构温度场的监测。


2.根据权利要求1所述的桥梁结构温度场监测方法，其特征在于，所述步骤S1中，对流换热qc的计算公式为：
qc＝hc(Ta-T)



其中，hc表示对流换热系数，Ta表示实测气温，T表示桥梁温度场，V表示实时风速；
辐射换热qr的计算公式为：
qr＝hr(Ta-T)-qra
hr＝Cε[(T+273)2+(Ta+273)2]·[T+Ta+546]
其中，hr表示辐射换热系数，Ta表示气温，T表示桥梁温度场，qra表示斜面天空辐射效应，C表示Stefan-Boltzmann常数，ε表示物体发射率。


3.根据权利要求1所述的桥梁结构温度场监测方法，其特征在于，所述步骤S2包括以下子步骤：
S21：根据日序数N计算太阳常数I0；
S22：利用总辐射传感器测量水平面上的太阳辐射I；
S23：根据太阳常数I0和水平面上的太阳辐射I计算大气清洁指数kT，根据水平面上的太阳辐射I计算地面发射辐射re×I，其中，re表示地面短波辐射反射率；
S24：根据大气清洁指数kT和水平面上的太阳辐射I计算天空散射辐射IdH；
S25：根据地面发射辐射re×I和天空散射辐射IdH计算投射到表面的太阳辐射量It；
S26：根据投射到表面的太阳辐射量It计算太阳辐射qs。


4.根据权利要求3所述的桥梁结构温度场监测方法，其特征在于，所述步骤S21中，太阳常数I0的计算公式为：
<...

【专利技术属性】
技术研发人员：陈钒，艾祖斌，侯圣均，王超，
申请(专利权)人：中电建路桥集团有限公司，
类型：发明
国别省市：北京;11