斜拉桥塔柱水化热温度应力分析

2014-05-09赵亚龙彭敬磊阙育梅

交通科技 2014年3期

赵亚龙彭敬磊阙育梅

（中国公路工程咨询集团有限公司北京 100097）

大体积混凝土的温度裂缝控制问题是一项国际性的技术难题［1－2］，其研究也较为深入，但对于一些大尺度的薄壁混凝土结构则未引起足够重视，大尺度薄壁结构，其壁厚较小，散热较快，当热量从表面散发，在混凝土中产生温度梯度，且越靠近表面温度梯度越大，由此温度梯度导致的应力足以使薄壁表面产生裂缝［3－4］，在混凝土箱梁、空心桥墩、空心塔柱上都能看到这种由于温度引起的开裂现象。

1 水泥水化热与混凝土绝热升温

1.1 水泥水化生热

混凝土结构浇筑完成后，水泥的水化过程尚在持续中，释放大量的水化热使混凝土内部温度上升，通过与外界的热交换（边界上热量的不断散失、太阳或其他外部热源的辐射补充），其温度逐渐与周边环境的温度趋于平衡。这期间，温度变化引起结构的变形，若结构受到约束则会产生应力，假使温度变化引起的拉应力超出了该时刻（龄期）混凝土的抗裂强度，受约束结构件的混凝土必然破坏，即不可避免地出现温度裂缝。因此，温度裂缝产生的关键问题在于混凝土内部温度瞬态场的分布、边界条件的限制、温度变化量及其相应的温度变形之间的转化。

现浇混凝土构件置于现场环境中，除混凝土水化发展的水化热外，还要受到环境温度变化的影响。环境温度的来源主要是太阳辐射、施工后加热养护等。

水泥的水化热是影响混凝土温度应力的一个重要因素，实际上温度场计算中用的是混凝土的绝热温升。测定混凝土绝热温升有2种方法：①直接法，用绝热温升实验设备直接测定；②间接法，先测定水泥水化热、再根据水化热及混凝土的比热容、容重和水泥用量计算绝热温升。相比较而言，直接法较准确。

混凝土浇注后，水泥与水结合发生水化作用，同时也是混凝土的凝固过程。随着水化热的逐渐减少及热量的散发，混凝土的温度就会慢慢降低。桥塔施工中的混凝土体积大，其温度在浇筑后2～3 d呈上升趋势，以后温度逐渐下降，一般经过较长时间才能达到稳定温度。在这个过程中，水泥水化作用生成大量的热量，由于混凝土自身的热物理性能和外界不同情况下的气温状况，形成不同的温度分布，构成整个桥塔的温度场。

1.2 混凝土温度场的计算原理

1.2.1 水泥水化生热计算

水泥水化放出的热量与龄期关系密切，选取复合指数公式［5］：

式中：Q（τ）为龄期τ时积累的水化热，kJ／kg；Q0为τ→∞时的最终水化热，kJ／kg；τ为龄期，d；a，b为系数。最终水化热Q0和系数a，b见表1。

表1 水泥水化热常数

1.2.2 混凝土绝热升温

混凝土绝热温升函数：

式中：θ（τ）为τ时刻绝热温升；W 为水泥用量；c混凝土比热容；ρ为混泥土密度；F为混合料用量；Q（τ）为水泥水化热；k为折减系数，计算取k＝0.25。

1.2.3 混凝土的弹性模量

由于温度应力的数值与弹性模量成正比，而且混凝土浇注以后，水化热的散发、温度场的变化与混凝土弹性模量的变化是同步发展的，所以在混凝土应力计算中，混凝土弹性模量的数值以及它与龄期的关系是很重要的。

混凝土弹性模量是龄期的函数，对于常规混凝土一般采用下面3种表达式。

（1）指数修正式

式中：β与a均为常数。

（2）复合指数式

式中：E0，a ，b均为常数。

（3）双曲线式

式中：E0为混凝土最终弹性模量；q为常数，当τ＝q时，E（τ）＝E0／2。各计算公式中的常数见表2。

表2 弹性模量公式计算中的常数

2 实例分析

某斜拉桥，主塔顺桥向为单柱式，横桥向为钻石型结构，主塔塔柱采用空心矩形断面，下塔柱横向宽度自上而下由4 m变化至6.5 m，塔柱顺桥向宽度下横梁至塔根范围由6.5 m过渡到7 m，下塔柱为外侧向外倾斜度1／2.6、内侧向外倾斜度1／2.2的变截面矩形空心柱。

2.1 ANSYS有限元分析步骤［6－7］

ANSYS的热分析有稳态和瞬态之分。在工程上一般用瞬态热分析计算温度场，并将其作为热荷载对结构进行热应力分析，即先进行温度场分析再进行结构应力分析，分析采用热－应力耦合法分析，选用ANSYS中SOLID5热－应力耦合单元。①用APDL语言建立模型，共分为12 586个单元，16 974个节点；②定义混凝土的导热系数为2.37 kJ／（m·h·℃），比热容为0.92 kJ／（kg·℃）；③边界上存在空气和混凝土的热对流，属于热分析中的第三类边界条件，对流边界条件可以作为面荷载（具体输入参数为对流系数和空气温度）施加于实体模型的表面，来计算固体和流体间的热交换，对内外表面施加不同的的对流系数：sfa，area，conv，热对流系，温度；④对混凝土施加水化热，水泥水化生热采用式（1），则热生成，采用 APDL方式设定水泥的热生成率，将热生成率作为体荷载加于单元上来模拟水泥的化学反应；⑤考虑混凝土弹性模量随龄期的变化，计算弹模选用复合指数形式；⑥通过运用＊DD循环语句和生死单元来模拟施工过程，计算时间为7 d。有限元模型见图1。

图1 有限元模型

2.2 空心塔柱水化热温度场模拟

计算中考虑混凝土分层浇注、层厚、水化热温升变化、弹性模量变化、材料不同分区、浇注温度、对流边界条件等条件和因素，选取载荷步长为1 h。ANSYS有限元温度场计算结果见图2和图3。

图2 各时刻桥塔截面温度分布云图

图3 塔柱内、外表面及中间部位温度随时间变化曲线

由图3可见，各个时刻温度沿截面变化的趋势差不多一致。随着时间的增长，由于水泥水化放热，内外各点温度都处于上升阶段。内表面最高温度值发生在35 h左右，中间部位和外表面温度曲线的发展规律大致与之相同，但桥塔内表面的水化热温度最高值明显高于外表面的最高温度，温度上升段较陡，下降段平缓。假定混凝土入模时的温度均为20℃，根据计算混凝土最高升温约54℃，升温过程中温度变化速率为1.13℃／h，降温过程中温度变化速率为0.28℃／h。