文章快速检索    
  同济大学学报(自然科学版)  2019, Vol. 47 Issue (6): 755-763.  DOI: 10.11908/j.issn.0253-374x.2019.06.003
0

引用本文  

孙建渊, 谢津宝. 基于等效龄期的钢管拱内混凝土硬化过程热应力[J]. 同济大学学报(自然科学版), 2019, 47(6): 755-763. DOI: 10.11908/j.issn.0253-374x.2019.06.003.
SUN Jianyuan, XIE Jinbao. Thermal Stress of Concrete-filled Steel Tube Arch During Hardening Process Based on Equivalent Age Method[J]. Journal of Tongji University (Natural Science), 2019, 47(6): 755-763. DOI: 10.11908/j.issn.0253-374x.2019.06.003

基金项目

国家自然科学基金(51778468)

第一作者

孙建渊(1966—),男,副教授,工学博士,主要研究方向为混凝土桥梁、钢-混凝土组合结构、桥梁计算机辅助设计.E-mail:sunjy@tongji.edu.cn

通信作者

谢津宝(1995—),男,硕士生,主要研究方向为混凝土桥梁、桥梁计算机辅助设计.E-mail:xiejinbao@tongji.edu.cn

文章历史

收稿日期:2018-09-17
基于等效龄期的钢管拱内混凝土硬化过程热应力
孙建渊 , 谢津宝     
同济大学 土木工程学院,上海 200092
摘要:为揭示大直径钢管拱内混凝土硬化过程中力学性能增长的温度依赖性因素对组合结构热力作用效应的影响机理,采用有限元程序模拟钢管拱内混凝土的水化热传导过程,并与实测温度场数据进行对比,随后基于等效龄期法考虑其对管内混凝土弹性模量增长的影响,在此基础上结合热弹性力学理论得到了硬化过程中组合结构热应力的变化规律,并与未考虑温度依赖性影响的计算结果进行比较分析.结果表明,水化热温度场加快了管内混凝土硬化过程中弹性模量的增长速度,进而导致混凝土温度应力明显增大,截面径向、环向以及纵向温度应力增幅分别可达1.3倍、1.3倍和1.4倍,但对钢管应力的影响可忽略不计.因此,在分析大直径钢管拱内混凝土硬化过程中的热力作用效应时,必须考虑水化热温度场对管内混凝土弹性模量增长的影响.
关键词钢管混凝土拱桥    硬化过程    水化热    温度依赖性    等效龄期法    热力作用    
Thermal Stress of Concrete-filled Steel Tube Arch During Hardening Process Based on Equivalent Age Method
SUN Jianyuan , XIE Jinbao     
College of Civil Engineering, Tongji University, Shanghai 200092, China
Abstract: In order to reveal the influence mechanism of temperature dependence of development of mechanical properties of concrete in large-diameter steel tube arches on the thermodynamic effect of composite structures during hardening process, a finite element model was established to simulate the hydration heat transfer process of the concrete in the steel tube arch, which was compared with the measured temperature field data. Subsequently, the effect of hydration heat on the development of the modulus of elasticity of the concrete was further considered based on the equivalent age method, by which the variation of thermal stress of composite structure during hardening process were revealed combining with thermal elasticity mechanics theory, and compared with the calculation results without considering the influence of the temperature dependence. The results show that the hydration heat temperature field accelerates the growth rate of the modulus of elasticity of the concrete in the tube during hardening, which leads to a significant increase in the thermal stress of the concrete during hardening process, and the increases in the radial, circumferential and longitudinal directions can reach 1.3 times, 1.4 times and 1.3 times, respectively, but the effect on the stress of the steel tube is negligible. Hence, in the analysis of the thermodynamic effect during the hardening of concrete in large-diameter steel tube arches, the effect of hydration heat temperature field on the development of the modulus of elasticity of the concrete in the steel tube must be considered.
Key words: concrete-filled steel tube arch bridge    hardening process    hydration heat    temperature dependence    equivalent age method    thermodynamic effect    

我国的钢管混凝土(concrete-filled steel tube, CFST)拱桥在近二十多年中得到了大量的应用,截止2015年1月,已建和在建的钢管混凝土拱桥(跨径在50 m及以上)已达413座[1-2].我国仍处于大规模修建钢管混凝土拱桥的时期,对其需求量还很大,跨径不断增加,采用的钢管管径随跨径而增大,目前钢管混凝土拱桥跨径已经超过500 m,可望在不久的将来达到700 m.

随着钢管混凝土拱跨径不断增加,同时钢管拱肋的直径不断加大,最大钢管拱肋管径已达1.5 m[3],大直径钢管灌注的混凝土具有大体积混凝土水化反应特性[4].处于密封条件下的钢管混凝土水化热引起较大的温度变化,不仅严重影响管内混凝土力学性能的增长规律,而且使得钢管与混凝土硬化过程产生复杂相互作用,导致与其相对应的截面应力以及成桥后的受力状态也受到影响[5].并且由于管内大体积砼水化热造成的混凝土温度裂缝内包在钢管内,人们难以发现或检测钢管混凝土的开裂情况,而这些温度裂缝将对大量钢管混凝土拱的承载力、耐久性等都会留下永久的严重质量隐患.

关于钢管混凝土拱桥水化热问题虽然展开了一定的研究,但研究数量整体偏少,并且现行研究主要偏向于水化热温度场研究[3, 6-8],鲜有对水化热作用下的结构温度应力进行研究,研究尚不深入.而且,由于大直径钢管混凝土管内泵送的混凝土是高强度微膨胀混凝土,掺加了微膨胀剂、粉煤灰等,材料组成较为复杂,其水化放热及其力学性能增长规律异于普通混凝土,而已有研究多集中在小直径钢管混凝土拱桥上[8-10],不能真实地反映大直径钢管混凝土拱桥硬化过程中的热力作用效应.特别是研究硬化过程温度应力时,受水化热温度场不均匀分布的影响,管内混凝土各部位的水化反应速率以及水化反应程度不同,进而导致同一龄期不同位置处混凝土力学性能(比如弹性模量)的发展存在显著差异,但已有研究未能综合考虑水化热温度场分布对管内混凝土力学性能增长的影响[10-17].

以某铁路上承式大直径钢管混凝土拱桥为工程背景,以热传导理论结合有限元法建立了大直径钢管拱内微膨胀高强度混凝土水化过程中的热传导分析模型,并通过引入等效龄期法进一步考虑水化热温度场对管内混凝土弹性模量增长的影响,在此基础上结合热弹性力学理论对钢管拱内混凝土硬化过程中的结构温度应力进行了热力作用数值模拟与比较分析.

1 钢管砼拱水化热温度场计算原理 1.1 水化热计算模型

大直径钢管混凝土硬化的过程是一个复杂的放热过程.假定浇筑后混凝土为均质、各向同性,即水泥砂浆、骨料等呈现均匀分布状态,其水化热的释放与空间位置无关,只是时间的函数,因此,应该计算出各龄期的累积水化热量,以便求导得到水化放热速率,并将其作为内热源代入后续瞬态温度场的计算中.而水泥累积水化热依赖于龄期,常见的计算模型有指数式(式(1))、双曲线式(式(2))和复合指数式(式(3))水化热表达式三种[18].

$ Q(t)=Q_{0}\left(1-\mathrm{e}^{-m t}\right) $ (1)
$ Q(t)=\frac{Q_{0} t}{n+t} $ (2)
$ Q(t) = {Q_0}\left( {1 - {{\rm{e}}^{ - p{t^q}}}} \right) $ (3)

式(1)~(3)中:Q(t)为在龄期t时的累计水泥水化热;Q0t→∞时的最终水泥水化热;t为龄期;m为常数; n为常数,它是水化热达到一半时的龄期,d; pq为常数.

钢管内灌注的复合胶凝材料混凝土水化热最好由实验测定,在缺乏直接测定的资料时, 建议按照式(4)进行估算[18-19].

$ Q(\tau)=Q(t)(W+k F) $ (4)

式中:Q(τ)为混凝土水化热;W为按配合比计算的水泥用量;F为按配合比计算的矿物掺合料用量;k为折减系数,对粉煤灰可取0.25[18].

1.2 水化热温度场模拟

假定没有日照作用、环境温度均匀分布,同时材料热工性能参数是各向同性的,由于有内热源的圆钢管混凝土的几何形状、热力学边界条件都对称于中心轴,因此,圆钢管混凝土水化热瞬态温度场问题可简化为轴对称热传导问题,采用柱坐标系来描述,以Oz轴为轴向对称轴.由于钢管混凝土温度沿环向θ分布不变,即T/θ=0,则得到管内混凝土的轴对称热传导方程如下:

$ \frac{{\partial T}}{{\partial \tau }} = {a_{\rm{c}}}\left( {\frac{{{\partial ^2}T}}{{\partial {r^2}}} + \frac{1}{r}\frac{{\partial T}}{{\partial r}} + \frac{{{\partial ^2}T}}{{\partial {z^2}}}} \right) + \frac{{\partial {\theta _c}}}{{\partial \tau }} $ (5)

式中:r为半径;τ为时间;T为温度;θc是混凝土绝热温升,℃;ac=kc/ccρc,为混凝土的导温系数;其中,kc为混凝土的热传导系数;cc为混凝土的比热容;ρc为密度.

$ \frac{\partial \theta_{\mathrm{c}}}{\partial \tau}=\frac{W q_{\mathrm{v}}}{c_{\mathrm{c}} \rho_{\mathrm{c}}} $ (6)

式中:W为单位体积混凝土的水泥用量;qv为单位时间单位质量水泥释放的水化热.

钢管壁外表面与内表面之间存在热传导,则钢管内热流密度qs表达式为

$ q_{\mathrm{s}}=\frac{k_{\mathrm{s}}}{t_{\mathrm{s}}}\left(T_{\mathrm{so}}-T_{\mathrm{si}}\right) $ (7)

式中:ks为钢材导热系数;ts为钢管壁厚;Tso为钢管壁外表面温度;Tsi为钢管壁内表面温度.

为了确定所求温度场的唯一解,热方程还需要满足下列初始条件和边界条件.前者为

$ \left.T\right|_{t=0}=T_{0}(r, \theta, z) $ (8)

边界条件给出导热物体在边界上的温度或换热情况,反映影响物体热传导和温度状态的外部因素.钢管混凝土水化热温度场分析主要存在下列三种边界条件(图 1).图中, qce为管内边缘混凝土热流密度;qsi为钢管壁内表面热流密度;Tce为管内边缘混凝土温度;Tf为大气环境温度;β为钢管壁外表面在大气环境中的放热系数; qso为钢管壁外表面热流密度.

图 1 钢管混凝土水化热温度场边界条件 Fig.1 Boundary conditions of hydration heat temperature field of CFST

边界条件C1,钢管壁内表面与管内边缘混凝土之间假设接触良好,在接触面上温度和热流密度都是连续的,用公式表示为

$ \left. \begin{array}{l} {T_{{\rm{si}}}} = {T_{{\rm{ce}}}}\\ {q_{{\rm{si}}}} = {q_{{\rm{ce}}}} \end{array} \right\} $ (9)

边界条件C2,钢管壁外表面与大气环境进行热对流,则钢管壁外表面热流密度qso表达式为

$ {q_{{\rm{so}}}} = \beta \left( {{T_{\rm{f}}} - {T_{{\rm{so}}}}} \right) $ (10)

式中:Tso为钢管壁外表面温度;β与风速v有关,可用式(11)计算[7]

$ \left\{ \begin{array}{l} 粗糙表面:\beta = 23.9 + 14.50v\\ 光滑表面:\beta = 21.8 + 13.53v \end{array} \right. $ (11)

边界条件C3,钢管混凝土拱肋的拱脚与拱座之间传热较少,即:

$ {\left. {\frac{{\partial T}}{{\partial n}}} \right|_\mathit{\Gamma }} = 0 $ (12)

式中:n为表面外法线方向.

热传导某一瞬间的钢管壁内热流量为稳定状态,则钢管内热流密度qs与钢管壁内、外表面热流密度qsiqso相等,联立式(7)、(9)、(10)可得式(13):

$ {q_{{\rm{so}}}} = {q_{\rm{s}}} = {q_{{\rm{si}}}} = {q_{{\rm{ce}}}} = \frac{1}{{\frac{1}{\beta } + \frac{{{t_{\rm{s}}}}}{{{k_{\rm{s}}}}}}}\left( {{T_{\rm{f}}} - {T_{{\rm{ce}}}}} \right) $ (13)

上述问题可以采用热传导有限元法来求解,并在时间离散上采用向后差分法,可以得到[20-21]

$ \left[\boldsymbol{H}+\frac{1}{\Delta t_{n}} \boldsymbol{R}\right] \boldsymbol{T}_{n+1}-\frac{1}{\Delta t_{n}} \boldsymbol{R} \boldsymbol{T}_{n}+\boldsymbol{F}_{n+1}=0 $ (14)
$ H_{i j}=\sum\limits_{e}\left(\boldsymbol{h}_{i j}^{\mathrm{e}}+\boldsymbol{g}_{i j}^{\mathrm{e}}\right) $ (15)
$ R_{i j}=\sum\limits_{e} \boldsymbol{r}_{i j}^{\mathrm{e}} $ (16)
$ F_{i}=\sum\limits_{e}\left(-f_{i}^{\mathrm{e}} \frac{\partial \theta}{\partial t}-p_{i}^{\mathrm{e}} T_{f}\right) $ (17)

其中:

$ \left. \begin{gathered} h_{ij}^{\text{e}} = \iint\limits_{\Delta R} {\left( {\frac{{\partial {N_i}}}{{\partial r}}\frac{{\partial {N_j}}}{{\partial r}} + \frac{{\partial {N_i}}}{{\partial z}}\frac{{\partial {N_j}}}{{\partial z}}} \right)r{\text{d}}r{\text{d}}z} \hfill \\ r_{ij}^{\text{e}} = \frac{1}{a}\iint\limits_{\Delta R} {{N_i}{N_j}r{\text{d}}r{\text{d}}z} \hfill \\ f_i^{\text{e}} = \frac{1}{a}\iint\limits_{\Delta R} {{N_i}r{\text{d}}r{\text{d}}z} \hfill \\ p_i^{\text{e}} = \int\limits_\mathit{\Gamma } {\frac{\beta }{{{k_{\text{s}}}}}{N_i}r{\text{d}}s} \hfill \\ g_{ij}^{\text{e}} = \int\limits_\mathit{\Gamma } {\frac{\beta }{{{k_{\text{s}}}}}{N_i}{N_j}r{\text{d}}s} \hfill \\ \end{gathered} \right\} $ (18)

式(14)~(18)中:HRF分别为整体热传导矩阵、整体热容矩阵和整体温度荷载向量;T为待求解的节点温度向量;Δtntn时刻的时间步长;hije为单元热传导矩阵;gije为热交换边界条件对单元热传导矩阵的修正;rije为单元热容矩阵;$f_i^{\rm{e}}\frac{{\partial \theta }}{{\partial t}}$为水化热源引起的单元温度荷载;pieTa为热交换引起的单元温度荷载;a为单元导温系数, 并区分混凝土单元和钢管单元;N为形函数矩阵;ΔR为单元求解域;Г为钢管外表面所在单元的边界.当Tn为已知、Tn+1为未知时,式(14)是关于Tn的线性方程组.由于初始温度是已知的,基于方程式(14)求解即可得到各节点在tn+1时刻的温度Tn+1.

2 钢管砼拱温度应力场计算原理 2.1 基于等效龄期的混凝土弹性模量增长模型

研究管内混凝土硬化过程温度应力时,受水化热温度场不均匀分布的影响,管内混凝土各部位的水化反应速率以及水化反应程度不同,进而导致同一龄期不同位置处混凝土弹性模量的发展存在显著差异,因此,需要采用等效龄期计算方法进一步综合考虑水化热温度分布对管内混凝土弹性模量整体增长的影响.通过等效龄期法换算,可以用参考温度下混凝土弹性模量的发展历程表征任意温度历程下混凝土弹性模量发展历程[22].

等效龄期的概念建立在混凝土成熟度理论基础上,Nurse[23]和Saul[24]认为在其他养护条件相同的情况下(比如湿度等),不同的温度历程达到相同的成熟度时,同一混凝土的力学性能相近.用Nurse-Saul方程表征的成熟度理论的积分形式为

$ M = \int_0^t {\left( {T\left( \tau \right) - {T_{{\text{ref}}}}} \right){\text{d}}\tau } $ (19)

式中:M为成熟度;τ为养护龄期;Tref为参考温度;T(τ)为龄期τ时的温度.

与成熟度的概念一致,等效龄期是指在不同温度历程时混凝土达到相同成熟度而需要在参考温度下养护的时间[25],即:

$ t_{\mathrm{e}}=\int_{0}^{t} \frac{K(T(\tau))}{K\left(T_{\mathrm{ref}}\right)} \mathrm{d} \tau $ (20)

式中:K(Tref)为参考温度下的反应速率系数;K(T(τ))为龄期为τ、温度为T时的反应速率系数.

而混凝土的水化是一个放热反应过程,可以用化学反应中的Arrhenius方程来表征温度对水泥化学反应速率的影响,即:

$ K(T) = A{{\text{e}}^{\frac{{ - {E_{\text{a}}}}}{{\mathit{\boldsymbol{R}}(\mathit{\boldsymbol{T}} + 273)}}}} $ (21)

式中:R为理想气体常数8.314,J·mol-1·K-1A为前因子;Ea为水泥的化学活化能,J·mol-1.

根据Arrhenius方程,Hansen和Pedersen等[26]提出了新的成熟度理论的表征方法,以标准养护温度20 ℃(293 K)为参考温度,即等效龄期teq

$ {t_{{\text{eq}}}} = \int_0^t {\exp \left[ {\frac{{{E_{\text{a}}}}}{R}\left( {\frac{1}{{{T_{{\text{ref}}}} + 273}} - \frac{1}{{T\left( \tau \right) + 273}}} \right)} \right]{\text{d}}\tau } $ (22)

在有限元计算的过程中, 采用时间步累计叠加的形式计算相应的等效龄期,即把龄期等分为n段,并假设每段时间间隔内的温度为常数,则等效龄期的计算简化为

$ {t_{{\text{eq}}}} = \sum\limits_{i = 1}^n {\exp \left( {\frac{{{E_{\text{a}}}}}{R}\left( {\frac{1}{{293}} - \frac{1}{{T\left( \tau \right) + 273}}} \right)} \right)\left( {{t_i} - {t_{i - 1}}} \right)} $ (23)

水化热温度场各节点每时刻的温度均不相同,因此各节点均要建立这样的等效龄期公式.根据离散化的等效龄期定义式(23), 在每个时间段内分别累积计算各节点的等效龄期te, 并代入式(24)[27]中求出该节点此计算步内的弹性模量值.

$ {E_{\text{c}}}\left( {{t_{\text{e}}}} \right) = \left( {1 - \exp \left( { - at_{\text{e}}^b} \right)} \right) \cdot {\left( {{E_{\text{c}}}} \right)_{\max }} $ (24)

式中:ab为常数;(Ec)max表示混凝土最终能达到的极限弹性模量值;Ec(te)表示等效龄期te时的混凝土弹性模量.

因此,对于配合比一定的同类型混凝土,在其他养护条件相同的情况下,不管温度与时间如何组合,只要在参考温度下的等效龄期相等,其对应的弹性模量也是相等的[25].

2.2 水化热温度应力场有限元分析

混凝土终凝之后完全丧失可塑性,强度和弹性模量开始增长,符合弹性变形规律,并假定混凝土是均质、各向同性的.因此,钢管混凝土水化热引起的结构温度应力可作为拟静态热应力问题求解,时间只是作为确定温度场的参数,任意时刻的温度应力都可看成是一个纯空间域问题[21].

物体热膨胀只产生线应变,剪切应变为零.这种由于热应变εT可以看做是物体的初应变,对于钢管混凝土的轴对称热应力问题,εT的表达式是:

$ {\mathit{\boldsymbol{\varepsilon }}_T} = {\left[ {{\varepsilon _r}{\varepsilon _z}{\varepsilon _\theta }{\gamma _{rz}}} \right]^{\rm{T}}} = \alpha \left( {\phi - {\phi _0}} \right){\left[ {\begin{array}{*{20}{c}} 1&1&1&0 \end{array}} \right]^{\rm{T}}} $ (25)

式中:α为单元材料的线膨胀系数;ϕ0为结构的初始温度场;ϕ为结构的瞬态温度场,可由温度场分析得到的单元结点温度ϕi通过插值求得,即:

$ \phi = \sum\limits_{i = 1}^{{n_e}} {{N_i}{\phi _i}} = N{\phi ^{\rm{e}}} $ (26)

在物体中存在初应变的情况下,应力应变关系可表示为[21]

$ \mathit{\boldsymbol{\sigma }} = {\mathit{\boldsymbol{D}}_{{\rm{eq}}}}\left( {\varepsilon - {\varepsilon _{\mathit{}{T}}}} \right) $ (27)

其中:

$ {\mathit{\boldsymbol{D}}_{{\rm{eq}}}} = {\lambda _{{\rm{eq}}}}\left[ {\begin{array}{*{20}{c}} {\frac{{1 - \mu }}{\mu }}&1&1&0\\ 1&{\frac{{1 - \mu }}{\mu }}&1&0\\ 1&1&{\frac{{1 - \mu }}{\mu }}&0\\ 0&0&0&{\frac{{1 - 2\mu }}{{2\mu }}} \end{array}} \right] $ (28)
$ {\lambda _{{\rm{eq}}}} = \frac{{{E_{{\rm{eq}}}}\mu }}{{\left( {1 + \mu } \right)\left( {1 - 2\mu } \right)}} $ (29)

式(28)~(29)中:μ为单元材料的泊松比;Eeq为引入等效龄期方法计算后的弹性模量,其中,对于混凝土单元弹性模量采用Ec(te),而钢管单元的弹性模量Es为常数;λeq为引入等效龄期方法计算后的Lame系数;Deq为同时考虑龄期和温度的单元弹性矩阵.

根据最小位能原理,只考虑温度应力问题,即假设结构不受外力和体力作用,并对泛函求驻值可得[21]

$ {\rm{ \mathsf{ δ} }}\mathit{\Pi }\left( u \right) = \int_\mathit{\Omega } {\delta {\mathit{\boldsymbol{\varepsilon }}^{\rm{T}}}\mathit{\boldsymbol{\sigma }}{\rm{d}}\mathit{\Omega }} = \delta {\mathit{\boldsymbol{a}}^{\rm{T}}}\int_\mathit{\Omega } {{\mathit{\boldsymbol{B}}^{\rm{T}}}{\mathit{\boldsymbol{D}}_{{\rm{eq}}}}\left( {\mathit{\boldsymbol{Ba}} - {\mathit{\boldsymbol{\varepsilon }}_T}} \right){\rm{d}}\mathit{\Omega }} $ (30)

式中:B为单元应变—位移关系矩阵;a为整体节点位移列向量.

再考虑到虚位移δaT是任意的,所以:

$ \int_\mathit{\Omega } {{\mathit{\boldsymbol{B}}^{\rm{T}}}{\mathit{\boldsymbol{D}}_{{\rm{eq}}}}\left( {\mathit{\boldsymbol{Ba}} - {\mathit{\boldsymbol{\varepsilon }}_T}} \right){\rm{d}}\mathit{\Omega }} = 0 $ (31)

将式(31)按单元逐一积分再叠加, 即得到温度应力问题的有限元基本方程如下:

$ \boldsymbol{K a}=\boldsymbol{P}_{T} $ (32)

式中:K为整体刚度矩阵;PT为温度变化引起的荷载列向量,且有:

$ \left. \begin{array}{l} {\mathit{\boldsymbol{K}}_{\rm{T}}} = \sum\limits_e {\int_{{\mathit{\Omega }_e}} {{\mathit{\boldsymbol{B}}^{\rm{T}}}\mathit{\boldsymbol{DB}}{\rm{d}}\mathit{\Omega }} } \\ {\mathit{\boldsymbol{P}}_{\rm{T}}} = \sum\limits_e {\int_{{\mathit{\Omega }_e}} {{\mathit{\boldsymbol{B}}^{\rm{T}}}\mathit{\boldsymbol{D}}{\mathit{\boldsymbol{\varepsilon }}_T}{\rm{d}}\mathit{\Omega }} } \end{array} \right\} $ (33)
3 工程案例 3.1 工程背景

以某铁路上承式钢管混凝土拱桥为工程背景[3],桥梁主跨360 m,矢跨比1/6,拱轴线为悬链线,拱轴系数为2.5.主拱肋采用钢管混凝土结构,截面形式为四肢桁式(图 2),2片桁架拱肋向中间倾斜8°,通过相贯线焊接形式的K形钢管连接.桁架拱肋由上下两个横置哑铃通过H型弦杆连接而成,哑铃截面单根钢管外径为1.5 m,钢管壁厚30~35 mm.拱肋钢结构采用Q370qD,拱肋内部填充C50微膨胀混凝土.首先进行1号管的浇注,管内混凝土达到90%设计强度后浇注另一根钢管混凝土.

图 2 钢管混凝土拱肋截面(单位:mm) Fig.2 Cross section of concrete-filled steel tube arch (unit: mm)
3.2 有限元模型

钢管混凝土拱肋水化热热力耦合分析采用有限元计算程序进行.由于物体的几何形状、热力学参数和边界条件等与θ角无关(即对称于轴线),并假定没有日照作用、环境温度均匀分布,同时材料热工性能参数是各向同性的,则可以推论温度分布与θ角亦无关.这就只需要取出物体的子午面进行轴对称单元建模,其有限元轴对称和三维扩展模型如图 3所示.图中,a、b、c为测点.热力耦合单元采用具有温度自由度和结构位移自由度的轴对称单元.混凝土单元与钢管单元共用节点、完全连接以满足界面热传导连续的要求,温度场分析时沿两端截面取为绝热边界,应力场分析时钢管混凝土两端截面固结.

图 3 水化热分析有限元模型 Fig.3 Finite element model of hydration thermal analysis

针对拱肋水化热分析,主要有以下假设:管内混凝土与钢管壁接触良好,不存在脱空;大气温度均匀分布,不考虑日照等因素引起的非均匀温度场[3];钢管混凝土拱桥选用的混凝土一般添加缓凝剂,且以泵送的方式浇注,浇注完成时间很短,可以忽略不计,因此假设在一根拱肋全部灌注完成后,混凝土才开始水化放热.

3.3 初始条件和材料参数

初始条件为浇注时刻混凝土和钢管的内部温度分布,考虑钢管相对混凝土不仅传热较快而且体积小很多,因此钢管温度会较快与管内混凝土一致,可以认为混凝土和钢管的初始温度都为定值,即t=0,T=T0(起始时刻为浇注完管内混凝土开始放热的时刻,T0取为混凝土浇注温度20 ℃).

实际工程中钢管混凝土拱肋水化热温度变化在20 ℃~80 ℃范围,温度对钢管混凝土两种材料热工性能影响可忽略不计[28].考虑到材料是各向同性的,有限元模拟分析时钢管混凝土材料热工性能参数取值见表 1[3, 29].

下载CSV 表 1 钢管混凝土材料热工性能参数 Tab.1 Thermal performance parameters of concrete-filled steel tube materials

采用双曲线式水化放热模型.根据实验室测定,水泥3 d水化热为364.7 kJ·kg-1,通过C50微膨胀混凝土配合比(表 2)拟合计算[3],即可得到管内混凝土水化放热曲线.钢管外表面与大气之间的换热系数按风速1.0 m·s-1选取9.814 W·m-2·℃-1 [18].

下载CSV 表 2 C50微膨胀混凝土配合比 Tab.2 Mix proportion of C50 micro-expansion concrete

水化反应的活化能Ea在0 ℃~100 ℃之间大致为常数,计算中Ea/R取为恒定值2 700 K[28, 30].对于管内灌注的C50微膨胀高强度混凝土而言,其弹性模量成长历程不同于普通的硅酸盐混凝土,参考温度20℃下弹性模量增长拟合公式[11-12, 14]

$ {E_{\rm{c}}}\left( {{t_{\rm{e}}}} \right) = 43.2\left( {1 - \exp \left( { - 0.806\;5t_{\rm{e}}^{0.609\;2}} \right)} \right) $ (34)
3.4 结果分析

为考察有限元结果计算的正确性,选取测点a,b,c的温度计算值与文献中实测值进行对比,如图 4所示.1.2节有限元模型计算的钢管混凝土拱肋截面温度场,与文献[3]中的实测温度场基本相符,最大误差不超过5 ℃,说明本计算模型可正确模拟钢管混凝土水化热温度场.

图 4 温度场数值计算结果与实测结果的对比 Fig.4 Comparison of the calculated and measured temperatures of typical points

钢管混凝土拱截面各点温度在浇筑完成后10 d内的变化曲线,如图 5所示.在水化热影响下,前期混凝土水化短时间内急剧放热,拱肋截面温度迅速升高,管内混凝土中心温度在混凝土灌注后26 h左右达到峰值70.5 ℃,且在截面上存在内高外低的温度差:越靠近截面边缘,节点温度峰值越低,最外侧节点(即钢管外表面)温度峰值为41.7 ℃,截面最大温差接近30 ℃.后期混凝土水化放热减少,拱肋截面温度相对缓慢下降,最终趋于环境温度20 ℃左右,内侧节点受环境温度影响小,外侧节点随环境温度变化最为明显.

图 5 水化热截面温度场分布 Fig.5 Cross section temperature field distribution caused by hydration heat

图 6是根据水化热温度场计算得到的混凝土截面上各点等效龄期随浇筑后实际龄期的变化规律.计算等效龄期选取的参考温度与标准养护温度一致,均为20 ℃,从图 6曲线上得出的混凝土截面各点等效龄期的成熟度可以认为与标准养护下该点达到该龄期时的成熟度一致.图 6表明,受水化热温度场影响,早龄期管内混凝土成熟度呈现不均匀增长的趋势,当实际龄期相同时,管内混凝土b、c点的成熟度均滞后于a点,且这一现象随龄期增长更显著,管内混凝土内外侧的等效龄期最大可相差4 d.

图 6 实际龄期与等效龄期关系 Fig.6 Relationship between actual age and equivalent age

图 7对比了采用两种不同计算方式所得管内混凝土弹性模量Ec增长的差异,其中曲线a,b,c分别代表截面上各点(测点a, b, c)基于等效龄期法考虑水化热温度场影响的计算结果,曲线Ec, ref代表未考虑该因素影响的计算结果.结果表明,相比参考温度下管内混凝土弹性模量的增长情况,考虑水化热温度场影响的弹性模量增长速度显然更快,水化放热开始后两者弹性模量增长差值迅速增大,在26 h时最大达到11 GPa,之后逐渐减小,在10 d左右趋于一致.并且,由于受到水化热内高外低的温度场分布的影响,混凝土弹性模量沿截面径向增长不同,混凝土中心弹性模量明显增长加快,靠近边侧增长变慢.

图 7 管内混凝土各点弹性模量发展规律 Fig.7 Developments of elastic modulus of the concrete in the steel tube at typical points

图 8a8b8c分别给出了采用两种不同计算方式所得管内混凝土径向温度应力σrc、纵向温度应力σzc以及环向温度应力σθc在前10 d内的变化规律,其中曲线a,b,c分别代表截面上各点(测点a, b, c)采用等效龄期法(即考虑了水化热温度场对其弹性模量影响)的计算结果,曲线arefbrefcref分别代表截面上各点未采用该方法的计算结果.

图 8 管内混凝土各点弹性模量发展规律 Fig.8 Temperature stress variations of the concrete in the steel tube at typical points

图 8a可知,考虑水化热温度场对弹性模量增长影响的条件下,随着前期水泥水化放热,升温膨胀的混凝土受到钢管的套箍作用产生径向温度压应力,并且快速增大,在38 h左右核心混凝土a点应力峰值达到4.4 MPa,滞后水化热温度峰值12 h,后期逐渐降低,到10 d左右消失.从整个截面应力场来看, 在钢管混凝土水化过程中, 由于管内混凝土的温度高于钢管温度, 致使整个截面呈现出径向压应力, 在钢管和混凝土的界面上不存在拉应力, 不会因出现脱粘现象而影响该钢管混凝土后期的正常使用性能.并且,混凝土径向温度压应力与水化热温度场分布趋势一致,均由核心混凝土中心向外衰减.从径向温度应力的对比来看,采用等效龄期法的计算结果明显更大,且越靠近截面中心差异越大,核心混凝土a点的径向温度压应力约为不采用时的1.3倍.

图 8b表明,与径向温度压应力的变化规律一致,管内混凝土截面纵向温度压应力呈现先增大到应力峰值后逐渐减小,并最终消失的趋势,但混凝土截面纵向温度压应力整体上均高于径向温度压应力,核心混凝土a点应力峰值可达18.8 MPa.并且,在混凝土截面上存在较大的纵向温度压应力梯度,其最大差值可达10.5 MPa.后期混凝土水化放热结束后,混凝土存在约0.5 MPa的残余温度拉应力.对比纵向温度应力,考虑温度对弹性模量增长影响的计算结果明显更大,且越靠近截面中心,两者差异越大,核心混凝土a点计算结果差异可达5.8 MPa.

图 8c可知,与径向、纵向温度压应力不同,水化热温度作用下混凝土整个截面同时存在环向温度压应力和拉应力,靠近核心处混凝土产生环向温度压应力,压应力峰值约为4.4 MPa,由于混凝土与钢管壁按共节点考虑,靠近钢管壁处混凝土产生一定的环向温度拉应力.从环向温度应力的对比来看,采用等效龄期法的计算结果明显比不采用时更大,核心混凝土a点的环向温度压应力约增大1.3倍.

图 9对比了采用两种不同计算方式所得钢管径向温度应力σrs、纵向温度应力σzs以及环向温度应力σθs在前10 d内的变化规律,其中曲线σrsσzsσsθ代表采用等效龄期法的计算结果,曲线σr, refsσz, refsσθ, refs代表未采用该方法的计算结果.

图 9 钢管温度应力时程曲线 Fig.9 Temperature stress variations of the steel tube

图 9表明,钢管由于在水化热温度作用下未受到径向约束,因此钢管径向温度应力大小趋近于零,有限元计算结果为0,与实际情况吻合.水化热释放过程中,钢管环向温度拉应力峰值为23.4 MPa,纵向温度压应力峰值为48.1 MPa.后期受变化的环境温度影响,钢管温度拉应力可达5 MPa.对比钢管三个方向的温度应力,采用等效龄期法的计算结果比不采用时略微偏大,但影响很小,可忽略不计.

4 结语

(1) 大直径钢管拱内混凝土硬化过程中的力学性能增长有较强的温度依赖性,而水化放热过程最高温度可达70.5 ℃,截面最大温差接近30 ℃,严重影响管内混凝土硬化过程中弹性模量的增长规律,导致增长速度加快,并沿截面径向增长不均匀,弹性模量增长差异与标准养护条件下相比可达11 GPa,但约10 d两者逐渐趋于一致.

(2) 钢管混凝土水化热温度场及受其影响的力学性能增长共同决定了硬化过程中组合结构热应力的变化规律.其中,在水化热释放过程中,整个混凝土截面呈现出径向温度压应力,靠近钢管壁处混凝土因受压而不会脱粘;混凝土截面纵向温度压应力最大,在核心混凝土中心处其峰值可达18.8 MPa;管内混凝土硬化过程中截面各点温度压应力与其至核心混凝土中心的距离成负相关,同等条件下,距离越远则温度压应力越小,其中混凝土截面纵向温度压应力最大梯度可达10.5 MPa.

(3) 相比未考虑温度依赖性影响的传统方法,采用等效龄期法考虑温度依赖性影响的计算结果发生显著的改变,其对混凝土径向温度应力、环向温度应力以及纵向温度应力的影响一致,发生显著的改变,均使硬化过程中混凝土应力明显增大,峰值应力增幅可达1.3倍、1.3倍和1.4倍,但对钢管应力的影响可忽略不计.因此,在分析大直径钢管拱内混凝土硬化过程中的热力作用效应时,有必要综合考虑水化热温度场对管内混凝土弹性模量增长的影响.

参考文献
[1]
陈宝春, 韦建刚, 周俊, 等. 我国钢管混凝土拱桥应用现状与展望[J]. 土木工程学报, 2017, 50(6): 50
CHEN Baochun, WEI Jiangang, ZHOU Jun, et al. Application status and prospect of CFST arch bridge in China[J]. Civil Engineering Journal, 2017, 50(6): 50
[2]
ZHENG J L, WANG J J. Concrete-filled steel tube arch bridges in China[J]. Engineering, 2018, 4(1): 143
[3]
高卫卫. 大跨度钢管混凝土拱桥水化热的试验与数值分析[J]. 铁道建筑, 2016(8): 35
GAO Weiwei. Experimental and numerical analysis on hydration heat of long-span concrete-filled steel tube arch bridge[J]. Railway Construction, 2016(8): 35 DOI:10.3969/j.issn.1003-1995.2016.08.08
[4]
中华人民共和国住房和城乡建设部.大体积混凝土施工规范: GB 50496-2009[S].北京: 中国计划出版社, 2009.
Ministry of Housing and Urban-Rural Development of the People's Republic of China. Mass concrete construction specification: GB 50496-2009[S]. Beijing: China Planning Press, 2009.
[5]
陈宝春. 钢管混凝土拱桥设计与施工[M]. 北京: 人民交通出版社, 1999
CHEN Baochun. Design and construction of concrete-filled steel tube arch bridge[M]. Beijing: China Communications Press, 1999
[6]
韩林海, 杨有福, 李永进, 等. 钢管高性能混凝土的水化热和收缩性能研究[J]. 土木工程学报, 2006, 39(3): 1
HAN Linhai, YANG Youfu, LI Yongjin, et al. Hydration heat and shrinkage properties of high performance concrete-filled steel tube[J]. China Civil Engineering Journal, 2006, 39(3): 1 DOI:10.3321/j.issn:1000-131X.2006.03.001
[7]
冯斌.钢管混凝土中核心混凝土的水化热、收缩与徐变计算模型研究[D].福州: 福州大学, 2004.
FENG Bin. Research on the calculation models of hydration heat, shinkage and creep of core concrete in the steel tube[D]. Fuzhou: Fuzhou University, 2004.
[8]
刘振宇.钢管混凝土拱肋截面温度场研究[D].福州: 福州大学, 2006.
LIU Zhenyu. Research on the temperature field of CFST arch rib section[D]. Fuzhou: Fuzhou University, 2006.
[9]
林春姣, 郑皆连, 黄海东. 钢管混凝土拱计算合龙温度试验研究[J]. 广西大学学报(自然科学版), 2010, 35(4): 601
LIN Chunjiao, ZHENG Jielian, HUANG Haidong. Experimental study on the calculated closure temperature of concrete-filled steel tube arch[J]. Journal of Guangxi University (Natural Science Edition), 2010, 35(4): 601 DOI:10.3969/j.issn.1001-7445.2010.04.014
[10]
张治成. 钢管混凝土拱桥混凝土灌注阶段的受力仿真分析[J]. 工程力学, 2007, 24(2): 146
ZHANG Zhicheng. Force analysis of CFST arch bridge during concrete casting[J]. Engineering Mechanics, 2007, 24(2): 146 DOI:10.3969/j.issn.1000-4750.2007.02.025
[11]
孙国富.大跨度钢管混凝土拱桥日照温度效应理论及应用研究[D].济南: 山东大学, 2010.
SUN Guofu. Theory and application of solar radiation effect of long-span CFST arch bridge[D]. Jinan: Shandong University, 2010.
[12]
刘恒.施工阶段钢管拱非线性温度影响研究[D].西安: 长安大学, 2012.
LIU Heng. Research on the effect of nonlinear temperature of steel tube arch during construction period[D]. Xi'an: Chang'an University, 2012.
[13]
宣纪明, 向华伟, 芦可琪. 钢管混凝土拱桥拱肋水化热温度场和温度应力分析[J]. 桥梁建设, 2010(3): 29
XUAN Jiming, XIANG Huawei, LU Keqi. Analysis on hydration heat temperature field and temperature stress of CFST arch bridge[J]. Bridge Construction, 2010(3): 29
[14]
王江龙.哑铃型钢管混凝土拱桥温度场及温度效应研究[D].西安: 长安大学, 2015.
WANG Jianglong. Research on the temperature field and thermal effects of dumbbell-shape CFST arch bridge[D]. Xi'an: Chang'an University, 2015.
[15]
YANG B, HUANG J H, LIN C J, et al. Experimental study on temperature fields of dumbbell-shape section of CFST arch rib and its effects[J]. Advanced Materials Research, 2011, 163
[16]
王友彪.大跨度钢管混凝土温度场及温度效应分析[D].成都: 西南交通大学, 2013.
WANG Youbiao. Research on the temperature field and thermal effects of long-span CFST arch bridge[D]. Chengdu: Southwest Jiaotong University, 2013.
[17]
荆旭.钢管混凝土拱管内混凝土水化热及灌注方案研究[D].重庆: 重庆交通大学, 2012.
JING Xu. Research on the hydration heat and casting scheme of CFST arch[D]. Chongqing: Chongqing Jiaotong University, 2012.
[18]
朱伯芳. 大体积混凝土温度应力与温度控制[M]. 第2版. 北京: 中国水利水电出版社, 2012
ZHU Bofang. Thermal stresss and temperature control in mass concrete[M]. 2nd ed. Beijing: China WatersPower Press, 2012
[19]
张庆欢.粉煤灰在复合胶凝材料水化过程中的作用机理[D].北京: 清华大学, 2006.
ZHANG Qinghuan. Mechanism of action of fly ash in hydration process of composite cementitious materials[D]. Beijing: Tsinghua University, 2006.
[20]
孔祥谦. 有限单元法在传热学中的应用[M]. 第3版. 北京: 科学出版社, 1998
KONG Xiangqian. Application of finite element method in heat transfer[M]. 3rd ed. Beijing: Science Press, 1998
[21]
王勖成. 有限单元法[M]. 北京: 清华大学出版社, 2003
WANG Xucheng. Finite element method[M]. Beijing: Tsinghua University Press, 2003
[22]
张增起, 石梦晓, 王强, 等. 等效龄期法在大体积混凝土性能预测中的准确性[J]. 清华大学学报(自然科学版), 2016, 56(8): 806
ZHANG Zengqi, SHI Mengxiao, WANG Qiang, et al. The accuracy of equivalent age method in predicting the performance of mass concrete[J]. Journal of Tsinghua University (Natural Science Edition), 2016, 56(8): 806
[23]
NURSE R W. Steam curing of concrete[J]. Magazine of Concrete Research, 1949, 1(2): 79 DOI:10.1680/macr.1949.1.2.79
[24]
SAUL A G A. Principles underlying the steam curing of concrete at atmospheric pressure[J]. Magazine of Concrete Research, 1951, 2(6): 127 DOI:10.1680/macr.1951.2.6.127
[25]
高原, 张君, 罗孙一鸣. 基于水泥水化度的混凝土早期力学性能发展预测[J]. 工程力学, 2013, 30(10): 133
GAO Yuan, ZHANG Jun, LUO Sunyiming. Prediction of development of early mechanical properties of concrete based on cement hydration[J]. Engineering Mechanics, 2013, 30(10): 133
[26]
HANSEN P F, PEDERSEN E J. Maturity computer for controlled curing and hardening of concrete[J]. Nordisk Betong, 1977, 41(19): 21
[27]
马跃峰.基于水化度的混凝土温度与应力研究[D].南京: 河海大学, 2006.
MA Yuefeng. Research on temperature and stress of concrete based on hydration degree[D]. Nanjing: Hohai University, 2006.
[28]
Neville A M. Properties of concrete[M]. New York: John Wiley&Sons Inc, 1996
[29]
LI Fangyuan, SHEN Yin. Full-scale test of the hydration heat and the curing method of the wet joints of a precast segmental pier of a bridge[J]. Revue Francaise de Génie Civil, 2017, 21(3): 1
[30]
王甲春, 阎培渝. 基于等效龄期的粉煤灰混凝土抗压强度计算模型[J]. 中山大学学报(自然科学版), 2014, 53(4): 83
WANG Jiachun, YAN Peiyu. Computational model of compressive strength of fly ash concrete based on equivalent age[J]. Journal of Sun Yat-sen University (Natural Science Edition), 2014, 53(4): 83