薄壁空心高墩温度效应仿真分析

2013-12-18,

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

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

0 引言

随着我国桥梁建筑技术的不断提高，现代桥梁朝着薄壁、大跨、轻型的方向发展。而桥墩作为桥梁的重要结构，为了能更好地适应现代桥梁特点的需要，改变了以往粗、大、实的结构，转向轻型、薄壁、高强、注重造型的方向发展，其中薄壁空心高墩是一种较好的形式[1]。

此类桥墩裸露于大气中，必然要长期受到太阳辐射和气温变化的影响。由于外界环境的瞬变性以及混凝土结构本身的导热性能差，使得空心墩结构内外表面形成了比较大的温差，当由此温差产生的温度变形受到约束的时候，就会产生温度应力。这种由温差作用引起的温度变形和温度应力称为温差效应。现有研究结果证明[2-3]，在薄壁空心高墩结构中，日照温差效应产生的应力有时与恒载、活载产生的应力属于同一个数量级，甚至使桥墩结构发生开裂，导致结构损坏，影响桥梁的正常运行。因此为了预防工程事故的发生，很有必要对薄壁空心墩结构的温差效应进行研究。

目前ANSYS有限元计算软件已得到成熟发展，其计算结果具有较高的可信度和准确性。因此以晋济高速公路仙神河大桥150.07 m高的八边形薄壁空心高墩为研究对象，利用ANSYS软件强大的二次开发技术和热分析功能，编制了薄壁空心高墩结构日照温度场和温差效应分析的专用程序，并结合实例验证了二次开发成果的准确性。同时利用二次开发成果，对八边形薄壁空心高墩在日照温差荷载作用下的温差效应进行了仿真分析，并得出有益的结论。

1 ANSYS的二次开发

ANSYS具有强大的二次开发功能，用户可利用ANSYS自身提供的APDL、UPFs、UIDL等二次开发工具，根据各具体专业的需求，在ANSYS系统中开发出具有中文界面的、特定功能的专用模块[4]。本文通过利用二次开发工具UIDL语言来编制ANSYS界面的菜单部分和编写对话框，利用APDL语言来编写程序命令流，二者结合起来，在ANSYS软件中开发出一个关于薄壁空心高墩结构从建模到热分析再到温度应力计算的专用程序。

1.1 用UIDL语言编制用户图形界面

UIDL语言(User Interface Design Language)是ANSYS自身提供的用户界面设计语言，主要用来编写或更改ANSYS的图形界面，它可以在ANSYS中添加自定义的菜单，制定对话框，还可以将用APDL语言编辑的专业有限元程序放进ANSYS中，从而建立有专业针对性的有限元计算模块，实现ANSYS的二次开发[5-7]。修改和编写用户界面是通过用UIDL语言修改以.GRN为扩展名的控制文件来实现的。本文通过在控制文件UIMENU.GRN和UIFUNC2.GRN中添加编制的自定义菜单代码及对话框代码，实现了在ANSYS程序的Main Menu部分添加自定义菜单以及对话框的编制。对话框以热分析参数定义为例，其界面如图1所示。

1.2 用APDL语言编写命令流

APDL是ANSYS Parametric Design Language的缩写，即ANSYS参数化设计语言。它是类似FORTRAN的一种解释性语言，只能在ANSYS环境中运行，主要包括参数、矩阵操作、数组、流程控制、重复功能和缩写功能、宏文件以及用户子程序等程序设计语言以及1 000多条ANSYS命令。利用APDL语言编写了关于薄壁空心高墩从建模到网格划分，从热分析到结构分析的程序命令流，主要包括：

图1 热分析参数定义对话框

material.ans 热参数定义程序命令流

piermodel.ans 导入空心墩模型的程序命令流

duntaiyang.ans 计算太阳辐射的程序命令流

dunduiliu.ans 对流计算程序命令流

dunchangbo.ans 长波辐射计算程序命令流

dunjisuan.ans 综合计算热流密度、划分网格、热学分析程序命令流

yuntu.ans 获取截面温度云图程序命令流

dunjiemian.ans 获取截面温度数据命令流

Bridge.ans 导入截面数据程序命令流

Yinglijisuan.ans 温度应力计算程序命令流

houchuli.ans 结果后处理程序命令流

………………

1.3 二次开发成果

结合APDL语言和UIDL语言,按照上述方法进行二次开发的成果为“薄壁空心高墩温度应力分析系统”，程序界面如图2所示。

图2 二次开发程序界面

利用本文开发的“薄壁空心高墩温度应力分析系统”，只需要输入或改变必要的参数，便可以很方便、很直观地实现在任意时间和地点、任意季节以及任意截面形状等条件下的高墩结构在日照温差作用下的温度效应的计算。计算流程如图3所示。

图3 简要流程图

图4 仙神河大桥桥墩示意(单位：m)

1.4 实例验证

为了验证二次开发成果“薄壁空心高墩温度应力分析系统”的准确性和实用性，对一具有实测数据的薄壁空心高墩结构的温度场及温度效应进行了分析。表1和表2给出了距离墩顶15 m处截面向阳侧壁板一点在不同时刻的计算值与实测值以及向阳侧壁板沿壁厚方向分别距离外壁板0 m、0.20 m、0.40 m、0.60 m、0.80 m、1.00 m处共计6个点在13：00时刻的温度计算值与实测值。

表1 向阳侧一点不同时刻温度数据计算值与实测值℃时间计算值实测值0:0014.82914.03:0014.62714.16:0015.22414.69:0024.39123.612:0030.93831.213:0031.24431.914:0030.66931.615:0029.36330.316:0028.02628.918:0026.09626.421:0023.24222.124:0020.13918.6

表2 向阳侧壁板13:00沿壁厚方向温度数据计算值与实测值℃距离外壁板的距离/m计算值实测值0.031.24431.90.219.51320.70.415.38716.60.614.65815.40.814.74014.71.015.60815.6

从表1可以看出空心墩向阳侧壁板一点温度计算值和实测值24 h变化规律比较吻合，某一时刻的温度分布也大致相同，两者最大误差为1.539 ℃，符合计算规范要求，而且测点温度走势基本一致；由表2可以看出，沿着空心墩壁厚方向的计算温度分布与实测温度分布变化趋势基本一致。这说明八边形空心墩温度分布的仿真计算是准确的，本文的二次开发成果也是可信的。

2 基于ANSYS的薄壁空心高墩温差效应分析

2.1 桥梁概况

本文研究的仙神河大桥是一座较为典型的独塔预应力混凝土部分斜拉桥，地处山西省和河南省交界处，桥墩高150.07 m，采用的是正八边形薄壁空心高墩结构，为典型的薄壁空心高墩，桥墩布置如图4所示。主墩截面八边形内切圆直径为墩顶10.04 m，墩底16 m，自墩顶到墩底截面采用直线斜率变化，墩壁厚度均为1.2 m。该桥所在的地域属于晋城市山区河谷地区，昼夜温差变化较大，而且由于地形地貌等条件的限制，沿墩身各部分接受太阳辐射的时间也不一样，对于墩顶部分一天中接受太阳辐射时间要多一些，而墩底部分则受日辐射作用很少。因此温差对空心墩的影响比较大。

2.2 薄壁空心高墩温差效应ANSYS计算结果分析

为了研究日照下薄壁空心高墩结构的温差效应，根据相关气象资料提供的近几年山西晋城地区的天气情况，选择了夏季中气温相对较高，太阳辐射强烈，风速小的晴好天气进行分析。具体日期为2008年8月6日(晴)，最高气温33 ℃，最低气温16 ℃，风速3 m/s。经计算，西侧壁板下午17:00时刻温差最大，因此本文以下午17:00时刻的温度分布作为控制荷载，对沿墩身高度方向的各个截面，进行了的温差效应计算。限于篇幅，以距离墩顶2 m处截面为例，其在下午17:00时刻的应力云图和整体变形云图如图5和图6所示。

图5 距墩顶2 m处截面最不利时刻应力云图

图6 薄壁空心高墩结构整体变形位移云图

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

(1)对图5进行分析可知，沿着八边形空心墩截面周向方向沿墩壁外表面的应力梯度较小，而径向方向的应力梯度较大，并且背阳面和向阳面沿着径向的应力分布规律一致，应力大致呈对称分布。

(2)由图5(a)可知，截面的最大应力值在角隅处出现，且墩顶截面的最大等效应力值达到了5.32 MPa，最大拉应力达到了1.3 MPa。

(3)由图5(b)可以看出，空心墩截面沿着墩壁厚度由外表面向内，在0 ～0.4 m段内，墩壁主要受压应力，且应力值较大，在0.4 m处，应力由最大压应力变化到最大拉应力，变化梯度较大；而0.4～1.2 m厚度范围内壁板全部受拉，应力值相对较小，拉应力分布区域达0.8 m，占整个壁厚的2/3，但应力有所减小，变化也不大。

(4)由图6可以看出，结构的平动变形最大值和扭转变形最大值均出现在墩顶，结构整体平动变形最大值为1.253 4 cm，最大扭转变形值约为0.578′，平动变形较大，扭转变形很小；且由墩顶至墩中再到墩底，结构的变形逐渐变小，靠近墩底部，结构的变形值为0。