双肢薄壁墩温度效应仿真分析

2014-03-22，

石家庄铁道大学学报(自然科学版) 2014年4期

，

(石家庄铁道大学土木工程学院，河北石家庄 050043)

0 引言

近几年，我国的交通事业取得了突飞猛进的成果，桥墩结构也不断地由原来的小跨径、实体墩、重型结构向轻型、薄壁、大跨径方向发展。随着对路线指标、造价和路线环境要求的日趋提高，高墩大跨连续刚构桥的应用已经屡见不鲜，目前双肢薄壁墩结构凭借着自身独具的优点，渐渐取代了单肢墩。它不仅能增加桥梁的刚度，允许桥梁纵向有较大的变位，而且两肢墩之间负弯矩为下凹的曲线，这样既减小了墩顶截面尺寸，又能充分发挥材料的受力性能，增加了桥梁美感。

桥梁结构长期暴露于自然环境中，受到太阳辐射、大气骤然降温和年温变化的影响，再加上混凝土自身导热性能差的特点[1]，使得结构物内外表面产生较大的温差，当由此温差产生的变形受到约束时，就会产生温度应力，从而使桥墩发生弯曲变形。根据理论实践证明[2]，此温度应力有时与荷载产生的应力属于同一个数量级，甚至引起桥墩结构的开裂，影响墩身结构的稳定性。在目前的规范中，对桥梁上部结构的温度规定很多，而对于下部结构很少[3]。由以下的计算也可以看出，在悬臂状态下的墩顶位移很大，当两端都受到约束时，即在铰支状态下的温度应力值也很大，这种温度应力达到一定值的时候，会使桥墩产生纵向裂缝。因此，对双肢薄壁空心墩结构进行温度效应的仿真分析具有必要性和紧迫性。

在前人研究的基础上，利用ANSYS的二次开发功能，编制了适合于计算双肢薄壁墩结构的有限元程序。以邢汾高速公路洺水大桥105.09 m高的双肢矩形薄壁空心墩为研究对象，利用编制的程序计算它的温度场和温度应力(悬臂状态)，并与实际数据进行对比，其数值相差不大，证明了二次开发成果的实用性。同时利用二次开发成果，对实际的双肢矩形薄壁空心墩进行温度场的模拟，得出了在悬臂状态下的墩顶位移比两端固定状态下的位移要大很多，因此有必要对日照作用下的墩顶位移进行施工控制。

1 双肢薄壁墩日照作用下的温度场计算

日照作用下的温度场需要同时考虑太阳辐射、桥墩各部分与外界之间的辐射换热(长波辐射)以及对流换热的影响，通过采用太阳辐射理论计算公式，并利用ANSYS的热分析和二次开发功能，编制了内嵌于ANSYS的集温度场计算和温度效应于一身的命令流[4-5]，使双肢薄壁墩结构的温度场仿真分析更简便。

1.1 太阳辐射

太阳辐射穿过地球大气层达到地面的辐射由于受大气分子、水蒸气、气溶胶等的散射和吸收而衰减，被大气散射的辐射有一部分到达地面，同时，到达地面的直接辐射和散射又被地面反射，因此，太阳辐射包括太阳直接辐射、大气散射和地面反射三部分[6]。下面对这三部分的理论计算公式进行说明。

1.1.1 太阳直接辐射

太阳直接辐射一年内有微小波动，可近似按下述公式计算[7]

(1)

式中,n表示自1月1日起算的日序数。

当太阳射线垂直于地面时的太阳直接辐射强度的近似计算公式为

Im=0.9mtuI

(2)

式中,m为经大气修正的大气光学质量；tu为大气浑浊度因子，在晴朗洁净的大气状态下取值1.8～3.3。

1.1.2 太阳散射

天空散射是从天穹的各个方向辐射到地面上的太阳辐射强度，主要是大气对太阳直接辐射的散射和反射形成的。它与结构物壁面的方位角、是否处于阴影状态无关，主要和大气的浑浊程度、壁面的倾角β和太阳高度角有关。水平面上的散射强度

Id=(0.271I-0.294Im)sinβs

(3)

式中,βs为太阳高度角。

任意壁面的散射强度

Iβ=Id(1+cosβ)/2

(4)

1.1.3 地面反射

太阳直接辐射和太阳散射到达地面时，一部分被地面反射，这种反射是一种短波辐射。对于结构物任意倾斜的接受面，反射辐射强度

If=υ(Im+Id)(1-cosβ)/2

(5)

式中,υ为地面反射系数，一般取0.2。

1.2 热辐射(长波辐射)

混凝土的表面通过电磁波的形式向外界发射能量，同时吸收来自大气及周围环境的辐射，从而形成了混凝土的辐射热交换过程。这种辐射热交换的热流密度

q=εaC0[(273+Ta)4-εd(273+Td)4]

(6)

式中,C0为Stefan-Boltzmann常数，其值为5.67×10-8(W/m2·K4)；εa表示混凝土表面辐射率，一般其值取0.88；εd表示大气的辐射率，一般取值εd=1-0.261exp[-7.776×10-4Td2]。

1.3 对流换热

对流换热是指由于固体表面与其接触的流体之间存在的温差而引起的热量交换。对流换热热流密度主要与接触面的性质、大小、流体速度、流动空间以及流体与接触面间的温度有关。根据牛顿公式，热流密度qr与接触面面积和温差成正比，即

qr=hf(Ta-Td)

(7)

式中,Ta表示固体表面的温度；Td表示周围流体的温度；hf为对流换热系数，与物体表面的光洁度和几何形状、风速等因素有关，通常由试验来确定。

2 内嵌于ANSYS的程序界面二次开发成果

在用ANSYS进行求解时，将太阳辐射、辐射换热和对流换热统一转换为热流密度，以TABLE表示，然后再施加到边界上。边界条件的气温、太阳辐射和风速都是时间的已知函数，那么热流密度总可以表示成时间和表面温度的函数。因此，可以给定fluxl在不同时间和表面温度下的值，定义一个随表面温度和时间变化的热流密度数表。然后用SFL命令对桥墩施加边界条件。因为定义fluxl时已经确定fluxl是随温度和时间变化的，因此在瞬态求解时，ANSYS就会根据温度和时间从fluxl表格中提取热流密度的值。

选用Beam189单元进行分析，通过热-结构耦合功能，将温度场单元PLANE77转化为结构场单元PLANE82，同时截面上的温度值通过材料号转换过来。将带有温度信息的截面定义为BEAM189单元的截面，从而可以计算桥墩结构的日照温度效应。

修改和编写用户界面是通过用UIDL语言修改以.GRN为扩展名的控制文件。TCL/TK是一套脚本语言，可以调用ANSYS数据库，以实现与ANSYS的数据交换，比UIDL功能更强大，但必须经UIDL和APDL调用。本文通过在控制文件UIMENU.GRN和UIFUNC2.GRN中修改编制的自定义菜单代码及对话框代码，用UIDL实现了在ANSYS程序的Main Menu菜单中修改及添加自定义的菜单[5,8-9]，用TCL/TK实现了对话框的编制。修改后的ANSYS界面如图1所示。

要实现从建模到分析再到求解的全过程，必须对Main Menu菜单中的各个新增加的菜单进行命令流的链接，即用APDL语言编写命令流。比如：如图1，“空心矩形截面平面模型导入”所链接的命令流就是“piermodel.ans”[10]。用TCL/TK编写的用于输入热分析参数数据的对话框如图2所示。

图1 二次开发程序界面图2 热分析参数对话框

3 双肢薄壁墩温度场实例验证及温差效应分析

3.1 桥梁概况

研究的洺水大桥处于温带大陆性季风气候，其四季分明、寒暑悬殊、春旱风大、夏热多雨、秋凉时短、冬寒少雪。年平均气温13.5 ℃，极端最低气温-22.4 ℃，极端最高气温41.8 ℃。年日照2 523 h，年平均风速1.4～1.9 m/s，最大风速21 m/s。现主要针对5,6号墩身结构进行研究：桥墩高105.09 m，采用双肢等截面矩形薄壁空心墩结构，两肢间的净距为5 m，壁厚0.7 m。为典型的双肢薄壁空心墩，桥墩布置如图3。桥墩周围的地貌属于剥蚀中低山地貌，地面标高一般为814.00～960.00 m，地形起伏较大，两侧山体较陡。沟底地形较平坦开阔，所以受太阳辐射较为均匀，可近似看做太阳射线平行的照射到墩身结构。

3.2 双肢矩形空心墩温度场实例验证

为了验证二次开发成果“双肢薄壁墩温度应力分析系统”的准确性和实用性，针对洺水桥这一双肢等截面矩形空心墩结构进行了温度应力的分析[11-12]。距离墩顶5 m处截面的有限元模型如图5所示。温度测定日期为2013年8月13号，当地10：00,12:00,15:00，17:00的温度分布云图如图6所示。表1和表2也分别列出了向阳侧一点不同时刻温度场以及沿壁厚0 m，0.2 m，0.4 m，0.6 m方向的温度分布。

图3 洺水1号大桥桥墩及温度计布置示意图 (单位：m) 图4 双肢墩横截面示意图图5 距墩顶10 m处截面有限元模型

图6 距墩顶10 m处截面不同时刻温度场

表1 向阳侧一点不同时刻温度计算值与实测值 ℃

表2 向阳壁板12:00沿壁厚方向温度数据计算值与实测值 ℃

由表1可看出，最高温差为下午17点，为1.002 ℃，符合规范要求。由表2也可以看出，沿壁厚方向温度分布的走势也基本一致，由此可以验证二次开发成果的准确性，为后面温度应力计算奠定了基础。

3.3 双肢薄壁墩温度效应计算结果分析

图7 双肢墩温度场模拟

为了研究双肢薄壁墩结构的温差效应，对洺水1号大桥进行了观测，研究发现，由于两肢墩的净间距为5 m，相互之间基本上没有遮挡，也就是说两肢墩之间互相没有影响，为了说明这一点，下面给出了双肢墩温度场模拟，如图7所示。

从图7中也可以看出，两个墩的温度场分布基本一致，因此可以单独研究一根墩，作为双肢墩变化的结果。选择在2013年8月13号对洺水桥进行分析，属于夏季中气温较高的一天，太阳辐射强烈，风速也小，桥墩处于最不利状态。具体情况为：最低气温18 ℃，最高气温32 ℃，风速2 m/s。由文献[4]可知，下午17:00时的温差最大，因此本文以17:00时距离墩顶5 m处截面分析，在悬臂状态下的温度应力云图和桥墩整体变形图如图8和图9。在铰支状态下的温度应力云图和桥墩整体变形图如图10和图11。

图8 距墩顶5 m处截面最不利时刻应力云图(悬臂状态)

图9 双肢薄壁墩结构整体变形位移云图(悬臂状态)

图10 距墩顶5 m处截面最不利时刻应力云图(铰支状态)

图11 双肢薄壁墩结构整体变形位移云图(铰支状态)

对上述温差效应计算结果进行分析，可以得出以下几点结论：

(1)由图8(a)和图10(a)可知，截面的最大应力在角隅处出现，且墩顶截面在悬臂状态下的最大应力为6.27 MPa，而在铰支状态下的最大应力达到了9.66 MPa。

(2)由图8可知，在悬臂状态下，截面的等效应力由向阳面向背阳面逐渐增大，且增大的梯度也较大，最大应力集中在角隅处，为6.27 MPa。径向应力背阳面主要以压应力为主，且沿墩壁厚度由外向内，压应力呈梯度变化-逐渐减小，在靠近墩壁内侧，压应力已变为拉应力，但此处的拉应力值较小。向阳面主要是拉应力为主，但拉应力值相对较小，且变化不大，最大拉应力集中在角隅处，为3.21 MPa。

(3)由图10可知，在铰支状态下，截面的等效应力主要以拉应力为主，且向阳面的拉应力几乎为零；背阳面的拉应力，沿墩壁厚度由内向外，拉应力值呈梯度逐渐增大，最大值集中在角隅处，为9.66 MPa。截面的径向应力主要以压应力为主，向阳面的压应力值很小，大部分区域还是以很小的拉应力为主；背阳面沿墩壁厚度由内向外，压应力值呈梯度逐渐增大，最大值集中在角隅处，值为9.66 MPa。

(4)由图9可知，在悬臂状态下，结构的平动变形最大值和扭转变形最大值均出现在墩顶，平动变形最大值为5.678 4 cm，最大扭转变形值约为3.623′，且由墩顶至墩中再到墩底，结构的变形逐渐减小，靠近墩底部，结构的变形值为0。

(5)由图11可以看出，在铰支状态下，结构的平动变形最大值和扭转变形最大值均出现在墩中附近，平动变形最大值为78.3×10-18cm,最大扭转变形值为(0.99×10-16)′，平动变形和扭转变形均较小，且由墩中向墩顶和墩底均逐渐减小，靠近墩顶部和墩底部，变形值为0。