(1.重庆大学 土木工程学院, 重庆 400045; 2. 山地城镇建设与新技术教育部重点实验室(重庆大学),重庆 400045)
(1.School of Civil Engineering, Chongqing 400045, China; 2.Key Laboratory of New Technology for Construction of Cities in Mountain Area(Chongqing University), Ministry of Education, Chongqing 400045, China)
DOI: 10.15986j.1006-7930.2018.04.004
为了研究裂缝对空心墩温度场和应力场的影响,以寒潮作用下的圆形空心桥墩为例,根据固体热传导微分方程和考虑温度应力的广义虎克定律,选择符合实际情况的边值条件,采用ANSYS软件计算了在寒潮作用下不同裂缝深度的混凝土桥墩温度场及其应力,计算结果说明:裂缝深度对混凝土圆形空心墩的温度场影响不大,空心墩外表面温度应力最大值出现在裂缝处,并且裂缝尖端的应力随着裂缝深度增加而增大,引起原有裂缝的继续发展,对混凝土桥墩产生不利影响.
In order to study the effect of crack on the temperature field and stress field of the hollow piers, this paper takes the circular hollow piers under the action of cold tide as an example according to the solid heat conduction differential equation and the generalized Hooke's law considering the temperature stress. In view of the differential equation of solid heat conduction and the generalized Hooke's law considering the temperature stress, the boundary conditions are chosen with the actual situation. The ANSYS simulation is performedfor calculating the temperature field and stress of concrete piers with cracks in different depths. Results show that the depth of the crack has little effect on the temperature field of the concrete circular hollow pier, and the stress at the crack tip increases with the increase of the crack depth, which causes the development of the original cracks and adversely affects the concrete piers.
对混凝土桥梁温度场及其效应的研究起源于20世纪50年代,美国的Zuk[1]推导出了结合梁顶、底面之间最大温差的近似方程,采用线性温度分布分析了结合梁的温度应力.意大利的Saetta[2]分析了混凝土结构在温度荷载作用下的温度应力是怎样计算的,提出了按二维平面应变有限单元法来分析箱型截面混凝土结构的温度应力.近年来研究人员对寒潮作用下的温度应力做了大量研究:河海大学的张子明等[3]研究了气温骤降引起大体积混凝土温度变化和温度应力的计算方法,给出气温骤降引起的温度变化规律及温度应力情况; 曾永平等[4]进行了铁路空心墩在寒潮作用下温度与应力场的实验研究,发现降温越快,对空心墩墩壁温度场影响深度越浅,桥墩外表面拉应力越大,越可能造成表面开裂; 长安大学的蒋国富[5]对洛河大桥的13#桥墩的温度场和温度应力进行了现场观测,推导出了适合陕北地区的混凝土薄壁空心桥墩的温度场表达式.东南大学的贾琳[6]指出我国公路桥涵设计规范所采用的温差值明显偏小.陈天地等[7]在高铁桥墩内预埋温度传感器,实测寒潮来临时墩壁温度场分布,总结出寒潮降温与墩壁内温度应力的对应关系.
以上工作都是在理想条件下研究混凝土结构的温度场,事实上混凝土结构广泛地存在着微小裂缝和缺陷,也有许多混凝土结构长期带裂缝工作,而大多数研究者并不太重视裂缝对混凝土结构温度场与应力场的影响.裂缝对温度荷载作用下墩身温度场及其应力的影响究竟如何?是否裂缝越深对温度荷载引起的混凝土墩身表面裂缝越不利?关于这些问题的研究资料目前还较少[8-13].又由于寒潮作用是铁路空心桥墩出现表面裂缝的主要原因之一[14],因此,本文以寒潮作用下的圆形空心桥墩为例,在空心墩外表面开设不同深度裂缝模拟实际开裂的桥墩,用有限元计算方法研究不同深度微裂缝对混凝土桥墩温度场及其应力的影响.
设圆形空心桥墩模型内径3 m,壁厚0.5 m,混凝土净保护层为3 cm.根据《铁路桥涵钢筋混凝土和预应力混凝土结构设计规范》中轴心受压构件的配筋构造要求,选取纵筋直径为20 mm,水平环状箍筋直径为12 mm,纵筋间距为20 cm,箍筋间距为15 cm; 连接内外纵筋的水平钢筋直径为10 mm,在竖向的间距为60 cm, 水平方向的间距约为50 cm, 呈梅花形布置.圆形空心桥墩模型钢筋示意图如图1所示,其中内侧纵筋和外侧纵筋间距均约为20 cm.
图1 圆形空心桥墩模型钢筋示意图
Fig.1 Schematic diagram of reinforcement in circle hollow pier model
圆形空心桥墩的热传导微分方程为[17]
式中,a=λ/(cpρ)称为固体材料的导温系数,单位为m2/s.若在寒潮作用下无内热源,则
边值条件是初始条件和边界条件的合称.寒潮作用下圆形空心桥墩的热传导过程是一个非稳态导热过程,求解该问题时必须说明初始条件.初始条件可以是各种各样的空间分布,即
T|τ=0=T0(x, y, z)(2)
寒潮作用下圆形空心桥墩壁内侧和外侧均与空气接触,其边界条件为第三类边界条件——假定经过混凝土表面的热流量与混凝土表面温度T(x, y, z)和气温Ta(x, y, z)之差成正比,即
式中:λ为材料的导热系数,W/(m·℃); β为对流系数,W/(m2·℃).
根据杜阿梅尔-诺伊曼的线性热应力理论,考虑温度应力的广义虎克定律,其形式如下[9].
εx-εx0=1/E[σx-μ(σy+σz)],
εy-εy0=1/E[σy-μ(σz+σx)],
εz-εz0=1/E[σz-μ(σx+σy)],
γxy=(τxy)/G, γyz=(τyz)/G, γzx=(τzx)/G(4)
式中,εx,εy和εz为物体任一点的实际应变值; εx0,εy0和εz0为该点不受约束自由热膨胀所产生的应变值.当物体在不均匀温度场作用下,如果暂时略去机械应力,则物体的实际应变可理解为某种平均温度意义上该点的应变值.所以在高温处ε-ε0<0,而受压应力; 在低温处ε-ε0>0,而受拉应力.
由于在日照作用下墩壁温度随高度的变化较小和寒潮作用下的空心桥墩内外空气温度几乎不随高度变化,所以可近似地认为相同结构形式和配筋的空心桥墩水平横截面上的温度及其应力分布不随高度变化,因此可将空间结构的温度及其应力问题简化为平面问题,采用ANSYS软件中的二维耦合场实体单元plane13进行温度及其应力的耦合计算.采用大型ANSYS有限元计算软件建模时,根据圆形空心桥墩结构的对称性和寒潮温度荷载的对称性,可选取1/4圆环建模,其有限元计算模型如图2所示.
图2 1/4圆形空心桥墩有限元模型
Fig.2 Quarter finite element model of circle hollow pier
图3 不同裂缝深度情况下的计算模型
Fig.3 Computational models with different crack depth
采用ANSYS有限元软件对上述圆形空心桥墩模型进行热力耦合瞬态计算.假设圆形空心桥墩模型的初始温度 26 ℃; 圆形空心桥墩壁内侧空气的对流系数为5,气温为26 ℃; 墩壁外侧空气的对流系数为12.5,气温为13 ℃.
为了节约篇幅,下面仅给出了寒潮来临后5 h和10 h四种圆形空心桥墩模型的温度场云图,见图4～5.从温度场云图可以看出,四种圆形空心桥墩模型的温度场随寒潮作用时间变化而变化,墩壁外侧温度变化较大而内侧温度变化很小.
图4 寒潮来临后5 h四种模型的温度场云图
Fig.4 Temperature cloud picture of four models 5 hours after cold wave
图5 寒潮来临后10 h四种模型的温度场云图
Fig.5 Temperature cloud picture of four models 10 hours after cold wave
寒潮作用下,圆形空心桥墩中产生的拉应力主要是沿水平圆周方向的.当该拉应力超过材料的强度极限时,桥墩将产生裂缝.因此,研究不同初始裂缝深度对混凝土圆形空心桥墩温度应力的影响时,温度应力主要是指环向应力.
对应上述温度云图,下面给出寒潮来临后5 h和10 h、四种圆形空心桥墩模型的温度应力云图,见图6～7.从应力云图可以看出,四种圆形空心桥墩模型的温度应力场随寒潮作用时间变化而变化,墩壁外侧环向应力均为拉应力,内侧环向应力均为压应力.
图6 寒潮来临后5 h四种模型的应力云图
Fig.6 Stress cloud picture of four models 5 hours after cold wave
图7 寒潮来临后10 h四种模型的应力云图
Fig.7 Stress cloud picture of four models 10 hours after cold wave
为了比较在寒潮作用下桥墩模型外侧温度应力随时间的变化情况,现将1-1截面上圆形空心桥墩壁外表面裂缝处温度应力随时间变化的有限元计算结果,绘制成曲线,如图8所示.由图8可知,寒潮作用初期三种桥墩模型外侧温度应力均随时间增加而快速增加; 寒潮作用5 h后外侧温度应力均随时间增加缓慢; 寒潮作用10 h后外侧温度应力达到最大值,然后再随时间增长而慢慢减小; 三种不同裂缝深度桥墩模型外侧温度应力有些差异:裂缝中心处应力以情况2、3、4的顺序增大,情况2最大应力为3.358 MPa,情况3最大应力为3.837 MPa,情况4最大应力为7.489 MPa.
图8 桥墩模型外侧应力随时间变化的对比曲线
Fig.8 Contrast curve of stress varied with time on the outside of pier models
本文采用大型有限元ANSYS软件对不同裂缝深度的混凝土空心桥墩和无裂缝混凝土空心桥墩进行了寒潮温度及其应力计算,并将计算结果进行了比较,得到以下结论:
①计算点到墩壁内表面越近,温度变化越小,不同裂缝模型之间温度差异也非常小;
②空心墩外表面温度最低,不同模型之间温度约有差异,情况1、2、3、4模型裂缝处温度依次有微小降低;
③混凝土寒潮温度应力由内到外增大,在圆形空心桥墩壁内侧应力最小,为压应力,外侧为拉应力,最大拉应力在裂缝中心处;
④空心墩壁裂缝处拉应力则因裂缝深度不同而异,例如在1-1截面上,墩壁外侧计算点拉应力以情况1、2、3和4模型的顺序增加,其中情况1模型壁外侧的拉应力最小,其值2.205 MPa大于C30混凝土轴心抗拉强度标准值2.01 MPa[10],情况4模型壁裂缝处的拉应力最大,其值7.489 MPa大于2.01 MPa.
该计算结果说明表面裂缝对空心墩温度场与应力的影响主要体现在裂缝尖端,不同深度的裂缝对空心墩温度场与应力影响不大,然而裂缝尖端的应力随着深度增加明显增大,并且远远超过相同深度其他位置的应力,对控制裂缝的继续发展非常不利.
