基于铁摩辛柯能量法的高墩大跨构桥稳定性分析

基于铁摩辛柯能量法的高墩大跨构桥稳定性分析

随着我国交通运输的快速发展，我国在千峰万谷、纵横变化地区修建的桥梁日益增多。这些峡谷地区的地形特点是高、宽、陡，所以修建的桥墩很高，一些超过100米。此外，这些桥梁通常由大型跨海桥和一条水道组成，桥梁结构的压力、结构强度和稳定性突出。高墩大跨径T构桥的高墩大多采用空心截面形式,为了跨越深谷,不得不设置高墩以满足线路的需要,为此,必须了解墩高与桥梁结构参数之间的关系。经典弹性稳定理论给出了等直杆临界力的精确公式,并通过长细比来界定等直杆的失效是稳定问题还是强度问题,但工程实际中的杆件多为变截面杆,没有对应的长细比计算公式可供参考,一般将稳定问题和强度问题一并考虑,按经验取值。目前,关于变截面高墩的文献大多利用有限元数值计算软件进行分析,给出解析解的结构截面形式较简单。在此,本文作者在经典稳定理论的基础上考虑强度问题,取施工阶段单壁墩T构桥为研究对象,推导出符合实际且简单的临界墩高和最大墩高理论计算公式,利用数学软件Maple编制程序,并给出界定大柔度压杆和小柔度压杆的坡度近似值。1截面几何特性分析T构桥有3种形式(图1),即跨中设铰的铰接T构桥、跨中带吊梁的T构桥和单跨T构桥,本文讨论施工阶段的T构静定体系(图2(a))。取变截面空心矩形桥墩进行分析(图2(b)与图3),按坡度变化给出截面几何特性沿墩高的表达式。主梁采用变截面箱形截面,顶板、腹板及底板皆按规范取值。为分析和计算方便,按设计图纸将主梁截面几何特性用数学软件Maple进行拟合,所有计算不计混凝土内钢筋。建立如图2(a)所示坐标系,桥墩与主梁的坐标原点均取在墩梁固接处截面形心处。图2(b)所示为变截面空心矩形桥墩墩顶截面形式,设墩顶截面横桥向(z方向)外壁宽a10,且与梁底等宽,内壁宽a20,顺桥向(y方向)外壁宽b10,内壁宽b20,墩身内外边缘均沿墩高按线性变化,横桥向墩身外壁坡度为k1,内壁坡度为m1,顺桥向墩身外壁坡度为k2,内壁坡度为m2,则墩身任意截面的截面尺寸为(图3):任意截面的横截面面积A和对y轴、z轴的惯性矩Iy和Iz分别为:图2(c)所示为主梁横截面形式和主梁坐标系,任意截面的横截面面积A1和对z1轴的惯性矩I1按设计图纸用数学软件Maple拟合为:2高码头结构分析2.1平衡状态的应变采用铁摩辛柯能量法进行分析,求结构的临界墩高,即结构失稳前所能达到的最大墩高。当结构体系处于稳定平衡状态,且给以微小扰动使其偏离原来的平衡位置时,体系的应变能将因此而增加ΔU,外力将增加作功ΔW。当结构体系从稳定平衡过渡到不稳定平衡而处在临界状态时,其能量关系为:目前,很多文献都认为成桥使用阶段的稳定性高于施工阶段的稳定性,因而在这里仅对施工时墩身自体施工阶段和最大悬臂施工阶段进行讨论。2.1.1构应变能及参数墩身结构如图2(a)所示,墩身自体施工阶段顺桥向失稳模态下的平面力学体系见图4(a),设墩身弯曲失稳的位移函数为:式(4)满足几何和力学边界条件,即结构应变能为:外力功为:其中:E为墩身弹性模量;γ为墩身混凝土容重。式(5)中的I在顺桥向失稳时取式(1)中的Iz,横桥向失稳时取式(1)中的Iy。将式(1),(4)~(6)一并代入式(3),可得l=f(c1,c2),因而,有:欲得到c1和c2的非零解,方程组(7)的系数行列式应等于0,即:其中:k11,k12,k21和k22为包含l和各已知参数的多项式,k12=k21。当给定截面参数以后,就可通过方程式(8)得到相应的临界墩高,临界墩高计算公式见文献。2.1.2主梁应力应变能上部结构施工至最大悬臂时最不安全,结构承受荷载有恒载、施工荷载、节段施工误差不平衡重力和风载。悬浇施工时,构件的顺桥向失稳变形情况包括挂篮正常工作和非正常工作(挂篮跌落)。在稳定公式推导时,以挂篮跌落为对象求临界墩高,结构如图2(a)所示,载荷如图4(b)所示,墩身弯曲失稳的位移函数与前面的相同,主梁位移函数设为:式(9)满足变形协调条件,即:在以下推导中,E1为主梁弹性模量,γ1为主梁混凝土容重,P为2倍挂篮自身重力(主梁一侧挂篮突然跌落,取动力放大系数2,即把失去的挂篮重力反向放大2倍,施加在原先挂篮位置,按静力问题计算),Mz为挂篮跌落时产生的不平衡弯矩。主梁应变能为主梁外力所作功为将式(4)~(6)和(9)~(11)代入能量法公式∆U+∆U1=∆W+∆W1,可得l=g(c1,c2),忽略方程组(12)的右边项,得到一个标准特征值问题,特征方程为其中:m11,m12,m21和m22为包含l和各已知参数的多项式,m12=m21。当给定截面参数以后,就可通过方程(13)得到相应的临界墩高。2.2强度分析强度分析时的2个阶段与稳定分析时的相同。2.2.1墩体施工阶段此阶段墩身最大应力满足下列方程:据此可计算出结构失效前所能达到的最大墩高。2.2.2有偏位时的不平衡弯矩此阶段需要考虑2种情况,即墩身顺桥向有偏位和横桥向有偏位。对应稳定分析时的顺桥向失稳和横桥向失稳,高墩偏位参照规范取l/3000,顺桥向有偏位时设偏位方向与挂篮跌落时产生的不平衡弯矩方向相同,最大应力满足下列方程:式中:Wz为墩身顺桥向抗弯截面模量;[σ]为许用应力,此时由于l/3000与l1相比非常小,可不计入Mz之中。由式(15)可计算出此阶段结构失效前所能达到的最大墩高。横桥向有偏位时设My为横桥向有偏位时产生的不平衡弯矩,最大应力满足下列方程:式中:Wy为墩身横桥向抗弯截面模量。由式(16)可计算出此阶段结构失效前所能达到的最大墩高。2.3各施工阶段墩身内、外壁坡度变化对临界墩高的影响工程实际中的桥墩多为空心变截面墩,壁厚沿墩高亦有变化,壁厚的变化可以通过改变内、外壁的坡度来实现。大部分桥墩的坡度变化沿墩高是唯一的,但有些桥墩的坡度变化沿墩高不唯一,如宜万线(湖北宜昌—四川万州)上的马水河桥,墩高较同类T构桥的高,其横桥向内、外壁坡度沿墩高分二段变化。本文假设坡度沿墩高唯一,并分别讨论自体施工阶段与悬臂施工阶段墩身内、外壁坡度的变化对临界墩高的影响。在施工过程中,失稳模态有横桥向和顺桥向2种,本文仅给出顺桥向失稳模态受力图(图4),强度计算中的受力图与稳定分析时的相同。内、外壁坡度变化时临界墩高的变化趋势见图6和图7,计算基本参数选取如下:E=3.25×1010Pa,E1=3.6×1010Pa,γ=γ1=26kN/m3,a10=7m,b10=10m,a20=3.8m,b20=6.8m,k1=0.005,k2=0.001,m1=0.002,m2=-0.0001,[σ]=18.4MPa。2.3.1外壁坡度较大时自体施工阶段墩高与坡度的关系见图6。从图6可见:自体施工阶段随横桥向外壁坡度k1的增加(图6(a)),横桥向失稳模态下的临界墩高的增长明显快于顺桥向失稳模态下的临界墩高以及强度条件下的最大墩高,说明横桥向外壁坡度的增加可以极大地提高横桥向的稳定性。因为增加k1可以使Iy迅速增长,而Iz和A只是略有增长,但临界墩高不能无限增长,这一点可以从强度条件下的最大墩高有上限看出。可见:针对理想压杆得到的稳定计算公式由于没有考虑初曲率或初偏心等缺陷,只能用来描述截面参数与临界墩高的关系,合理的墩高应该同时满足强度条件和稳定性要求。从图6(a)可见:外壁坡度较小时,墩的高度由稳定条件控制,外壁坡度较大时由强度要求控制。因为计算时采用的是工程实际中常用的空心墩,当外壁坡度较小时,壁厚沿墩高增加较小,自身重力增加较小,强度条件容易满足,同时,惯性矩较小稳定条件难于满足;而外壁坡度较大时,自身重力增长迅速,强度条件成为控制条件,而稳定性要求由于Iy的迅速增长不再成为控制条件。实际上外壁坡度较大时高墩已经不再是细长杆,对于经典稳定理论来说,非细长杆若误用了欧拉公式计算临界力,其结果比实际大(图5中虚线部分所示σcr与λ关系曲线),横截面上的正应力有可能超过比例极限,造成危险。横桥向失稳模态下临界墩高和考虑强度要求时最大墩高与k1的关系曲线在k1≈0.0126有交点,表明以此为界区分大柔度压杆和小柔度压杆,k1≥0.0126时桥墩成为小柔度压杆,因而可以解释为什么随横桥向外壁坡度的增长而临界墩高呈无限上升的态势,这与临界应力总图上误用欧拉公式计算小柔度压杆临界应力时曲线无限上升相符。随顺桥向外壁坡度k2的增加(图6(b)),可以得到与图6(a)基本类似的结论,只需将横桥向换成顺桥向,k1换成k2,Iy换成Iz即可。顺桥向失稳模态下临界墩高和考虑强度要求时最大墩高与k2的关系曲线在k2≈0.0112有交点,k2≥0.0112时桥墩成为小柔度压杆。自体施工阶段随横桥向内壁坡度m1的增加(图6(c)),横桥向失稳模态下的临界墩高先略有增加后缓慢减小,在m1≈0.0017处有极值;而顺桥向失稳模态下的临界墩高缓慢增加,强度条件下的最大墩高逐步减小。因为增加m1将使A和Iy及Iz均减小,A减小使得强度下降;但空心杆件的稳定性与Iy/A或Iz/A相关,A减小使得横截面积减小,自身重力减小,因而Iy及Iz的减小对高墩的稳定性影响不大。随顺桥向内壁坡度m2绝对值的增加(图6(d)),横桥向失稳模态下的临界墩高缓慢减小,顺桥向失稳模态下的临界墩高先略有增加后缓慢减小,在m2≈-0.0026取极值,强度条件下的最大墩高逐步增加趋于一个上限,因为m2绝对值的增加使得墩底截面趋于实心自身重力增加迅速,同时,Iy及Iz略增加。2.3.2各阶段验算强度条件的对比悬臂施工阶段墩高与坡度的关系见图7。从图7可见:悬臂施工阶段随横桥向外壁坡度k1的增加(图7(a))以及顺桥向外壁坡度k2的增加(图7(b)),可以得到与自体施工阶段基本类似的结论,只是对应同样的坡度,悬臂施工阶段的临界墩高均小于自体施工阶段的对应值,说明自体施工阶段的稳定性强于悬臂施工阶段的稳定性。悬臂施工阶段与自体施工阶段所不同的是:此阶段验算强度条件时需要考虑高墩的偏位,考虑横桥向偏位时最大墩高比考虑顺桥向偏位时的小。因为考虑横桥向偏位时高墩为双向弯曲与压缩的组合变形,考虑顺桥向偏位时高墩为单向弯曲与压缩的组合变形。横桥向失稳模态下临界墩高和考虑高墩偏位时最大墩高与k1的关系曲线在k1≈0.0131有交点,顺桥向失稳模态下临界墩高和考虑高墩偏位时最大墩高与k2的关系曲线在k2≈0.0117有交点。随着横桥向内壁坡度m1的增加,横桥向失稳模态下临界墩高和顺桥向失稳模态下临界墩高均缓慢减小,强度条件下不论横桥向偏位还是顺桥向偏位最大墩高都逐步减小(图7(c))。随顺桥向内壁坡度m2绝对值的增加,可以得到与自体施工阶段相类似的结果(图7(d))。2.4变截面墩理论墩与y的比较经典弹性稳定理论只能给出等截面墩在自体施工阶段的精确解(ql)cr=7.837EI/l2,当k1=k2=m1=m2时,变截面墩即退化为等截面墩。表1所示为等截面墩自体施工阶段在自身重力作用下的墩高。从表1可见:本文方法比精确解略大。这是因为假定的y不是真实的弹性曲线,相当于增加了约束使其实现,提高了临界荷载的值。将表1中数据与图6比较可以发现:在满足强度条件和稳定条件的前提下,变截面墩理论墩高比等截面墩大,但墩太高时,施工需满足竖直度等要求,若高墩的墩身由于施工原因出现了偏斜、弯曲等几何缺陷,将会使结构稳定性大大下降甚至产生整体失稳的严重后果。图6(a)和7(a)以及图6(b)和7(b)所示k1或k2界限值与临界应力总图上欧拉公式与经验公式界限点λp相对应,当k1或k2大于界限值时,临界墩高无限上升的态势与临界应力总图中误用欧拉公式计算小柔度压杆时,压杆横截面上应力无限增大相符,说明本文方法可以用于判别实际工程中的变截面压杆是大柔度压杆还是小柔度压杆,并藉此确定是否需要对工程中的变截面承压杆件进行稳定校核。2.5各级桥墩设计参数马水河桥是宜万线上1座铁路桥梁,位于四川境内建始县业州镇小溪口村,为(116+116)m预应力混凝土T构桥,其设计基本情况为:主梁箱梁横截面为单箱单室直腹板变截面,底板底缘按圆弧变化,圆曲线半径为631.01667m,中支点梁高12.20m,跨中梁高4.70m,顶宽10.70m,底宽6.00m,顶板厚0.34m,腹板宽0.45~0.90m,底板厚由0.32m变化至1.35m。桥墩采用空心矩形截面,墩高108.00m,墩顶横向宽7.00m,墩底横向宽28.00m,横向壁厚由墩顶1.40m变为距墩底7.00m处的7.50m,墩顶纵向宽10.00m,墩底纵向宽12.00m,墩纵向壁厚1.40m。按坡度唯一进行设计,k1=1/16.89=0.0592,k2=1/86.40=0.0116,m1=0.022,m2=-0.00127(其余参数与理论分析部分的相同)。k1>0.0126,k2>0.0112,故自体施工阶段不存在静力稳定问题;k1>0.0131,故悬臂施工阶段亦不存在静力稳定问题,这与文献中不计该桥静力稳定考虑动力稳定的结果相符。3正截面参数的变化(1)变截面空心高墩考虑稳定性要求时的临界墩高和满足强度条件时的最大墩高受截面几何性质影响很大,与等直墩相比,高墩截面参数特别是墩身内外壁坡度变化明显增加墩的高度,与内壁坡度相比,增加外壁坡度可以更有效地提高高墩的稳定性。(2)判别实际工程中的变截面压杆是大柔度压杆还是小柔度压杆的坡度界限值可以确定是否需要对工程中的变截面承压杆件进行稳定校核。该方法可推广至任意桥梁自体施工阶段的设计,亦可用