无铺装层PK断面混凝土梁日照温度场分析
2016-05-17张玉平谢文昌李传习
张玉平,谢文昌 ,李传习
(长沙理工大学 土木与建筑学院,湖南 长沙 410004)
无铺装层PK断面混凝土梁日照温度场分析
张玉平,谢文昌 ,李传习
(长沙理工大学土木与建筑学院,湖南长沙410004)
摘要:针对我国桥涵设计规范未给出无铺装层PK断面混凝土梁日照梯度温度及国内外缺乏研究的现状,应用ANSYS求解了马鞍山长江公路大桥右汊斜拉桥无铺装层混凝土主梁温度场,并进行了现场实测,计算值与实测值吻合较好。对无铺装层混凝土梁日照温度场的参数分析表明,吸收率对无铺装层混凝土梁最大竖向温差影响较大,风速次之。计算得到了该桥主梁在施工过程中出现的最大竖向梯度温度及其分布的指数函数的梯度温度模式。该模式与现有各国规范的梯度温度模式比较表明,本文研究得到的无铺装层混凝土梁梯度温度分布趋势与各国规范已有相关规定的分布趋势基本吻合。但温度分布受地域影响,各国规范梯度温度取值存在差异,建议补充完善我国公路桥规中关于无铺装层混凝土梁梯度温度的规定。
关键词:桥梁工程;混凝土梁;PK断面;温度场;梯度温度
0引言
太阳辐射将使混凝土主梁截面产生非线性梯度温度。当这种梯度温度引起的变形被结构内、外约束阻碍时,会产生相当大的温差应力[1]。近年来,极端高温天气频繁出现,太阳辐射对桥梁结构的影响越来越受重视,这种影响定量分析的关键首先是合理确定其温度场。我国《公路桥涵设计通用规范(JTG D60—2004)》(简称《公路桥规》)[2]给出了桥梁结构梯度温度荷载的相应取值,但仅限于有铺装层的混凝土梁和带混凝土桥面板的钢结构,且是将这一荷载作为可变作用施加于成桥之后的结构。对于施工过程中无铺装层混凝土梁存在的梯度温度在设计计算中是否考虑、如何取值,并未规定。研究表明,桥面铺装对混凝土梁竖向日照温差有明显的削弱作用[3-4],施工过程中无铺装层混凝土梁的温度效应不容忽视。
近50年来,国内外学者基于热传导理论和实测数据,运用有限元方法对混凝土梁日照温度场进行了大量研究[5-11]。Soukhov[12]、肖建庄[13]和顾斌[14]先后通过对气象参数的概率分析,拟合得到桥址地区有铺装层混凝土梁出现的竖向梯度温度模式。文献[3]在对无铺装层混凝土梁实测数据的基础上提出了呈指数函数分布的梯度温度模式,但仅适用于桥梁所在地区。另外,国内外对混凝土梁温度场的研究主要是针对形状规则的箱形截面梁,对于平行设置的双幅箱梁桥,叶见曙[15]认为可按照单幅箱梁来考虑。文献[4]中发现双箱对腹板内侧温度没有大的影响,而对腹板外侧影响较大,这种差异主要与桥梁所在地理位置及桥梁走向有关。而PK断面(美国 Pasco-Kennewick 桥首次采用了双三角形边箱形式主梁,这种断面形式简称为 PK 断面)主梁的温度场与普通箱形截面梁分布是否相同,则值得研究。
本文以马鞍山长江公路大桥右汊斜拉桥为工程背景,借助有限元软件ANSYS对无铺装层PK断面混凝土梁的日照温度场进行计算,并与混凝土PK主梁的温度场实测数值进行了对比;对影响混凝土梁温度场分布的主要因素进行参数分析;采用指数函数的形式对其竖向梯度温度进行拟合,并与《公路桥规》、铁路桥规、英国规范等设计规范中有关混凝土梁梯度温度的规定进行比较研究。
1热传导基本理论及方程
处于自然环境中的混凝土梁日照温度场应按瞬态传热问题进行分析,在笛卡尔坐标系中,三维瞬态导热微分方程的一般形式为[16]
(1)
式中λ为导热系数;Φ为单位时间内单位体积中内热源产生的热量;ρ为混凝土的密度;c为比热容;t为时间。
对于混凝土桥梁结构,沿桥纵向温度均匀分布,故可简化为常物性、无内热源的二维瞬态导热微分方程
(2)
式中α=λ/ρc,为热扩散率(或称为导温系数),单位为m2/s。
求解混凝土梁的温度分布,实质上是在定解条件下求解微分方程(2)。对于瞬态热传导问题,定解条件包括初始条件和边界条件,初始条件即计算初始时刻混凝土梁的温度分布,可表示为
(3)
式中T0(x,y)为物体的初始温度分布函数。
边界条件则是反映影响物体热传导和温度状态的外部因素, 在日照作用下,混凝土桥梁与外界的热交换,主要有来自太阳辐射的热流,以及与周围环境之间的对流及辐射换热,其边界条件是由第二类和第三类边界条件组成的混合边界条件,可以表示为:
(4)
式中,Γ为物体边界;n为边界Γ的外法线方向;h为物体边界与周围流体的表面换热系数;q(t)为通过边界由外界流入物体内部的热流密度。
式(4)可以进一步改写成由太阳辐射热流密度qs、对流换热热流密度qc及辐射换热热流密度qr表示的传热边界条件:
(5)
qs一般按式(6)计算:
(6)
式中,α为表面热辐射吸收率(0<α<1);I为垂直于表面的太阳辐射强度。
qc根据Newton冷却定律表示为:
(7)
式中,hc为对流换热系数;Ta为大气温度(K);T为混凝土表面温度。
qr根据Stefen-Boltzmann定律表示为:
(8)
改写成对流换热的形式,即可把qr作为第三类边界条件来考虑:
(9)
(10)
式中,hr为辐射换热系数;ε为辐射率,对混凝土通常取0.88;σ为Stefen-Boltzmann常数,取值为5.677×10-8W/(m2·K4)。
将式(6)、(7)和(9)代入式(5),可将混凝土箱梁的箱外边界条件最终处理为第三类边界:
(11)
同理,箱内边界条件也可以处理为这种形式,由于箱内没有太阳辐射,相应的综合大气温度中不包含太阳辐射项。
2有限元模型的建立
采用有限元软件ANSYS模拟无铺装层PK断面混凝土梁的日照温度场,选用具有二维热传导能力的平面单元PLANE55模拟平面温度场[17],PLANE55单元具有4个节点,每个节点只有一个自由度,单元示意图见图1。
按马鞍山长江公路大桥右汊斜拉桥标准断面的实际尺寸建立模型,对其进行网格划分,整个模型由5 142个节点、4 335个单元组成,如图2所示。以6:00时刻的箱梁实测温度作为初始条件,将太阳辐射、对流及辐射换热按变换后的第三类边界条件施加。根据实测数据或参考相关文献,确定计算所需的混凝土材料物理参数值和边界参数值,见表1~表4[4,16,18-20]。利用Table数组定义荷载,计算时间为2013年9月1日6:00-9月2日6:00,全过程分为25个荷载步,每个荷载步分为3个子荷载步。部分计算结果见图3。
图1 PLANE55单元示意图Fig.1 Schematic diagram of PLANE55 element
图2 主梁断面模型图Fig.2 Model of girder section
参数密度ρ/(kg·m-3)比热容c/[J·(kg·K)-1]导热系数k/[W·(m·K)-1]吸收率α辐射率ε数值24729602.710.50.88
表2 混凝土梁各表面综合换热系数h[单位:W/(m2·K)]
注:综合换热系数按公式h=12.47+3.33v[19]计算得到,根据当天实测风速波动较小的特点,取各表面风速均值计算。
表3 混凝土梁各表面不同时刻太阳辐射强度I(单位:W/m2)
注:本文中太阳辐射强度是参照文献[16],根据桥梁所在地理位置和混凝土梁各表面朝向确定太阳常数、太阳时角和太阳入射角等参数,通过计算得到。
表4 桥址处大气温度Ta(单位:K)
注:箱梁内外温度是通过在相应位置设置测点,实测得到。
图3 实测温度和计算温度变化图Fig.3 Variations of measured and calculated temperatures
3温度场分析
图4 截面尺寸与测点布置图(单位:cm)Fig.4 Section dimensions and measuring point positions(unit:cm)
在马鞍山长江公路大桥右汊斜拉桥的施工过程中,选择马鞍山侧主梁(桥面尚未铺装)20#梁段中间断面进行温度场观测,测点布置如图4所示,全断面共布置了24个测点。测试元件为热电偶传感器。选择无云、微风条件下日照辐射较为直接的天气进行温度场周期性观测,选取一组实测数据(2013年9月1日6:00—9月2日6:00)和对应ANSYS 分析结果进行比较分析。图3表明,计算值与实测值吻合较好,ANSYS模型较好地模拟了混凝土梁的日照温度场及内部各点温度变化趋势。从图中可以看出,腹板温度与底板温度变化趋势相类似,腹板温度较底板要高,且腹、底板温度峰值时刻较顶板均延迟3 h左右。
根据《公路桥规》竖向梯度温度取值的规定,本文取PK断面混凝土梁各个腹板距顶面0.4 m范围内竖向温度梯度分布,见图5。图5所示的箱梁南侧边腹板、南侧中腹板、北侧中腹板和北侧边腹板的最大竖向温差分别为8.2,11.7,11.5 ℃和9.3 ℃,可见各腹板的最大竖向温差并不一致,中腹板最大竖向温差明显高于边腹板。
图5 竖向梯度温度分布图Fig.5 Distributions of vertical gradient temperatures
4混凝土梁日照温度场参数分析
无铺装层PK断面混凝土梁日照温度场的计算模型已得到实测结果的验证,可以用于无铺装层混凝土梁日照温度场的影响因素的参数分析。在无铺装层混凝土梁日照温度场的影响因素中,混凝土密度、比热容可以通过实测得到,导热系数可以通过实测数据反算[4],混凝土辐射率在各文献中的取值介于0.85~0.95之间,离散性均不大;太阳辐射强度可以通过桥址位置和表面朝向等参数计算得到;箱梁内外大气温度可以通过实测得到。受到表面粗糙度、颜色等因素影响,混凝土表面太阳辐射吸收率在各文献中的取值波动较大;影响混凝土表面对流换热系数的主要因素——风速存在一定的随机性。因此,下面研究吸收率和风速对无铺装层混凝土梁日照温度场的影响。应用控制变量法,按照本文取值±30% 的增幅改变二者的取值,对吸收率和风速进行参数分析。计算结果如图6所示。
图6 参数影响分析曲线(单位:%)Fig.6 Analysis curves of parameter impact(unit:%)
图6表明,无铺装层混凝土梁最大竖向温差变化与吸收率变化、风速变化基本呈线性关系,吸收率的影响较大,风速次之。吸收率决定混凝土对太阳辐射的吸收量,吸收率变化±30%会导致最大竖向温差变化±26%左右。风速则是影响无铺装层混凝土箱梁温度场的另一个重要参数,风速的大小影响着箱梁表面与大气对流换热的效率。
5竖向梯度温度的计算
研究温度场的主要目的是确定结构的竖向梯度温度,以便计算结构由于梯度温度产生的应力。对于无铺装层PK断面混凝土梁,可利用已得到验证的模型计算其在极端天气条件下的温度场,从而确定施工过程中存在的最不利梯度温度。根据桥址所在地区的气象数据,2013年当地极端天气出现在8月12日,最高气温和最低气温分别为41 ℃和30 ℃,微风。应用本文第2节中的方法计算太阳辐射强度、综合换热系数。由于每天的大气温度大致按正弦曲线规律变化,可根据最高气温和最低气温按下式计算各个时刻的大气温度[16]:
(12)
箱内温度变化较为平缓,可按实测数据中大气温度与箱内温度的差值关系进行换算;对于初始条件的确定,可先假设一初始温度,通过计算2个周期(48 h)得到趋于稳定的箱梁温度,以此时的箱梁温度作为初始条件,进行箱梁温度场计算。计算得到的箱梁中腹板竖向梯度温度可用指数函数17.1e-8.3y表示,从图7中可知,底板温度与腹板相差很小。对于PK断面,计算与实测结果均表明截面下缘基本无竖向温差,建议不考虑截面下缘梯度温度。
图7 竖向梯度温度拟合曲线Fig.7 Fitted curve of vertical gradient temperature
计算得到的竖向梯度温度模式与各国规范规定的梯度温度模式对比见表5及图8。对比结果表明,本文研究得到的无铺装层PK断面混凝土梁的梯度温度分布趋势与各国规范相关规定的趋势基本吻合,但具体数值有一定差别,本文取值与新西兰规范相差最大,温度最大值相差14.9 ℃。由于温度分布受地域影响,选取准确的温度荷载是保证结构温度计算准确的前提,且对结构而言无铺装层温度梯度应比有铺装层的温度梯度更加不利,因此有必要补充完善《公路桥规》中关于无铺装层混凝土梁梯度温度的规定。
表5 梯度温度模式对比
图8 梯度温度模式对比图Fig.8 Contrast of curves of gradient temperature modes
6结论
(1)通过建立ANSYS模型,计算了无铺装层PK断面混凝土梁日照温度场,同时采集了实桥温度场数据,计算值与实测值吻合较好。结果表明各腹板的最大竖向温差不一致,中腹板的最大竖向温差明显高于边腹板。
(2)参数分析结果表明,无铺装层混凝土梁最大竖向温差变化与吸收率变化、风速变化基本呈线性关系。吸收率对混凝土梁最大竖向温差影响较大,风速也会产生一定的影响。
(3)计算得到无铺装层混凝土梁竖向梯度温度的指数函数形式,其分布趋势与各国规范的无铺装层温度梯度分布趋势基本一致,但其数值与各国规范相关规定值相比存在一定的差别,建议补充完善《公路桥规》中关于无铺装层混凝土梁梯度温度的规定。
参考文献:
References:
[1]刘兴法.混凝土结构的温度应力分析[M].北京:人民交通出版社,1991:1-4.
LIU Xing-fa. Analysis of Thermal Stress for Concrete Structure[M]. Beijing: China Communications Press, 1991: 1-4.
[2]JTG D60—2004,公路桥涵设计通用规范[S].
JTG D60—2004,General Code for Design of Highway Bridges and Culverts[S].
[3]方志,汪剑.大跨预应力混凝土连续箱梁桥日照温差效应[J].中国公路学报,2007,20(1):62-67.
FANG Zhi, WANG Jian. Sun Light Thermal Difference Effect on Long-span PC Continuous Box Girder Bridge[J]. China Journal of Highway and Transport, 2007, 20(1): 62-67.
[4]聂玉东.寒区大跨径混凝土箱梁桥温度场及温度效应分析[D].哈尔滨:哈尔滨工业大学,2013.
NIE Yu-dong. Analysis of Temperature Field and Temperature Effect for Long Span Concrete Box Girder Bridges in Cold Regions[D]. Harbin: Harbin Institute of Technology, 2013.
[5]汪剑,方志.混凝土箱梁桥的温度场分析[J].湖南大学学报:自然科学版,2008,35(4):23-28.
WANG Jian, FANG Zhi. Temperature Variation of Concrete Box Girder Bridge [J]. Journal of Hunan University:Natural Science Edition, 2008, 35(4): 23-28.
[6]吴六政.混凝土箱梁桥温度场的模拟[J].公路交通科技,2011, 28(10):65-69.
WU Liu-zheng. Simulation of Temperature Field of Concrete Box-girder Bridge[J]. Journal of Highway and Transportation Research and Development,2011, 28(10): 65-69.
[7]ROBERTS-WOLLMAN C L, BREEN J E, CAWRSE J. Measurements of Thermal Gradients and Their Effects on Segmental Concrete Bridge[J]. Journal of Bridge Engineering, 2002, 8(3): 166-174.
[8]MONDAL P, DEWOLF J T. Development of Computer-Based System for the Temperature Monitoring of a Post Tensioned Segmental Concrete Box Girder Bridge[J]. Computer-Aided Civil and Infrastructure Engineering, 2007, 22(1): 65-77.
[9]YAN P. Research on Temperature Field and Temperature Stress of Large Span Concrete Cable-Stayed Bridge Girder[C]//Applied Mechanics and Materials. Yantai:[s.n.], 2012:2073-2080.
[10]LEE J H, KALKAN L. Analysis of thermal Environmental Effects on Precast, Prestressed Concrete Bridge Girders: Temperature Differentials and thermal Deformations[J]. Advances in Structural Engineering,2012,15(3):447-460.
[11]LARSSON O, THELANDERSSON S. Estimating Extreme Values of Thermal Gradients in Concrete Structures[J]. Materials and Structures,2011,44(8): 1491-1500.
[12]SOUKHOV D. Representative Values of Thermal Actions for Concrete Bridge[J]. Progress in Structural Engineering and Materials, 2000, 2(4): 495-501.
[13]肖建庄,宋志文,赵勇,等.基于气象参数的混凝土结构日照温度作用分析[J]. 土木工程学报,2010,43(4):30-36.
XIAO Jian-zhuang, SONG Zhi-wen, ZHAO Yong, et al. Analysis of Solar Temperature Action for Concrete Structure Based on Meteorological Parameters[J]. China Civil Engineering Journal, 2010, 43(4): 30-36.
[14]顾斌,陈志坚,陈欣迪.基于气象参数的混凝土箱梁日照温度场仿真分析[J].东南大学学报:自然科学版,2012,42(5):950-955.
GU Bin, CHEN Zhi-jian, CHEN Xin-di. Simulation Analysis for Solar Temperature Field of Concrete Box Girder Based on Meteorological Parameters[J]. Journal of Southeast University: Natural Science Edition,2012, 42(5): 950-955.
[15]叶见曙,贾琳,钱培舒.混凝土箱梁温度分布观测与研究[J].东南大学学报:自然科学版,2002,32(5):788-793.
YE Jian-shu, JIA Lin, QIAN Pei-shu. Observation and Research on Temperature Distribution in Concrete Box Girders[J]. Journal of Southeast University: Natural Science Edition, 2002, 32(5): 788-793.
[16]彭友松.混凝土桥梁结构日照温度效应理论及应用研究[D].成都:西南交通大学,2007.
PENG You-song. Study on Theory of Solar Radiation Thermal Effects on Concrete Bridges and Application[D]. Chengdu: Southwest Jiaotong University, 2007.
[17]美国ANSYS公司北京办事处.ANSYS热分析指南[M].北京:美国ANSYS公司北京办事处,2000.
Beijing Office of US ANSYS Inc. ANSYS Thermal Analysis Guide[M]. Beijing: Beijing Office of US ANSYS Inc., 2000.
[18]张建荣,徐向东,刘文燕.混凝土表面太阳辐射吸收系数测试研究[J].建筑科学, 2006,22(1):42-45.
ZHANG Jian-rong, XU Xiang-dong, LIU Wen-yan. A Test Study on the Solar Radiation Absorption Coefficient of Concrete Surface[J]. Building Science, 2006, 22(1): 42-45.
[19]刘照球.混凝土结构表面对流换热研究[D].上海:同济大学,2006.
LIU Zhao-qiu. Research for Convection Heat Transfer of Surface of Concrete Structure[D]. Shanghai: Tongji University, 2006.
[20]张玉平,杨宁,李传习.无铺装层钢箱梁日照温度场分析[J].工程力学, 2011, 28(6):156-162.
ZHANG Yu-ping, YANG Ning, LI Chuan-xi. Research of Temperature Field of Steel Box Girder without Pavement Caused by the Solar Radiations[J]. Engineering Mechanics, 2011, 28(6):156-162.
Analysis of Temperature Field of PK Section Concrete Girder without Pavement Caused by Solar Radiation
ZHANG Yu-ping,XIE Wen-chang,LI Chuan-xi
(School of Civil Engineering and Architecture, Changsha University of Science & Technology, Changsha Hunan 410004, China)
Abstract:Aiming at lack of solar gradient temperature of PK section concrete girder without pavement in Chinese highway bridge design specifications and few research has been done at home and abroad, the temperature filed of Maanshan Yangtze River expressway right branch cable-stayed bridge concrete girder without pavement is calculated by ANSYS software, and the result agrees well with the field measured values. The parameter analysis on the concrete girder without pavement shows that the absorptivity has significant influence on its maximum vertical temperature, while the wind speed is another important influencing factor. The vertical maximum gradient temperature of the girder in the worst environment and the mode of its exponential function mode in the course of the bridge construction are calculated. Compared with the modes of gradient temperature in the current specifications of different countries, it shows that the distribution trend of the gradient temperature obtained from our research basically tallies with those in existing regulations in the specifications of different countries. Owning to the effect of geography, the temperature values of gradient temperature are different. Therefore, it is suggested to add and improve the specification of gradient temperature for concrete girder without pavement in Chinese highway bridge regulations.
Key words:bridge engineering; concrete girder; PK section; temperature field; gradient temperature
中图分类号:U445.47
文献标识码:A
文章编号:1002-0268(2016)04-0059-07
doi:10.3969/j.issn.1002-0268.2016.04.010
作者简介:张玉平(1976-),男,辽宁建平人,副教授,工学博士.(zyp5032@163.com)
基金项目:国家重点研究发展计划(九七三计划)项目(2015CB057702);国家自然科学基金项目(51378080);长沙理工大学土木工程优势特色重点学科创新性项目(15ZDXK02)。
收稿日期:2015-05-04