桥墩二类稳定分析报告

1、桥墩二类稳定分析报告2009.12.11结构稳定分析理论按照结构在逐渐加载过程中平衡形式是否发生质变的观点，将结构的失稳区 分为第一类稳定问题和第二类稳定问题，如图 I-图2所示。前者表示在加载过 程中，构件的平衡状态将出现分枝现象，使原有的平衡状态失去稳定性而转向新 的平衡状态；而后者在加载过程中平衡形式并不发生质变。 在第一类稳定同题中， 当加载至Per时，表示平衡的分枝即将出现，称 Per为压屈荷载。在第二类稳定 问题中，当加载到Per时，表示构件的承载能力即将降低，称为压溃荷载。工程 上通常把两者统称为 失稳的临界荷载。工程问题中研究结构稳定问题的目的， 在 于寻求相应的临界荷载及其临

2、界状态，防止不稳定平衡状态的发生，从而确保结 构安全。结构的第一类稳定问题，在数学上归结为广义特征值问题,Kd Kg、J式中，Kd为弹性刚度矩阵，Kg为几何刚度矩阵，它只与结构的轴力有关；为结构的位移增量；为载荷稳定系数。式（1）数学上表现为广义特征问题，应 用各种迭代方法，如逆矢量迭代法、子空间迭代法等都可以很方便地求解。结构的几何非线性是指大位移 问题。在载荷作用下，当位移大到足以使得结 构的几何形状发生显著的改变时，必须按已变形的位置建立平衡方程。弹性大变 形问题与变形的历史不相关，可以采用全量方法研究，也就是直接求在已知载荷 的约束下的总变形和应力，结构全量形式的非线性平衡方程可表示为

3、：K（、）、-工 1 =0（2）式中：K（、J为非线性刚度矩阵，为节点位移，F为等效节点载荷。非线性稳定分析的基本方法是逐步地施加一个恒定的载荷增量直到解变得 开始发散为止。由于特征值屈曲载荷是预期的线性屈曲载荷的上限，故可以作为非线性屈曲的给定载荷，在渐进加载达到此载荷前，非线性求解应该发散。设为特征值屈曲分析时的稳定系数，在这里取非线性屈曲分析的给定载荷系数。 将 给定载荷分为n个载荷增量，0 = 0 ： 1t r n = (3)在每一个载荷步内对非线性方程(2)进行线性化，可得增量形式的平衡方程：Kt_： -F =0(4)心二 K。KiF =( t - t)F式中：Kt为切线刚度矩阵，为

4、等效节点载荷增量矩阵，为节点位移 增量，Ko为小位移线性刚度矩阵，Ki为大位移矩阵，K二为初应力矩阵，F为 等效节点载荷矩阵。在每一个载荷步内，为了改进求解的精度，可以应用牛顿法 进行迭代。结构的极限承载力在开始发散的载荷和在此前一级已收敛的载荷之 间。如载荷增量步分得较细，可以偏于安全地认为是前一级载荷， 而避免更加复 杂的计算。目前二类稳定分析都是采用有限元等数值计算方法。从本构关系来考虑稳定 问题可以分为弹性稳定与弹塑性稳定， 而弹性稳定又从是否考虑几何非线性、 初 始缺陷等因素又分为线弹性稳定和非线性弹性稳定。由于施工环节会存在不可避免的施工误差，桥墩轴线可能出现偏差，高墩在水平荷载作

5、用下将产生较大的变 形等等，稳定计算必须计人初始缺陷及大位移的影响， 所以基于极值点失稳为理 论基础的计人双重非线性的弹塑性稳定问题对于工程是必须的。2桥墩ANSYS分析2.1单元模拟为在ANSYS软件中，采用beam189梁单元模拟桥墩和横梁，beam189为3-D二 次有限应变梁单元，可以很好的模拟桥墩力学行为。2.2材料本构材料非线性本构关系是进行二类稳定计算的重要依据，尤其是极限阶段混凝 土的强化过程。研究中只考虑了混凝土抗压部分，选用美国学者Hongn estad提出的本构曲线。具体表达式为：55- - 0；-；0式中，6为棱柱体极限抗压应力，p为对应于二0的应变，取0.002， ；

6、u通常取 0.0038。2.3稳定的评价（1）一类稳定实质上是特征值失稳，计算公式见下式，一般不小于45=豆_45PT（2）第二类稳定实质上是极值点失稳。用极限状态法设计桥梁时，稳定与最终的极限承载力是统一的。因为二类稳定分析的最终状态也是桥梁达到塑性变 形而破坏，因此，桥梁结构的二类稳定安全系数与强度安全系数也是一致的。按照公路钢筋混凝土及预应力混凝土桥涵设计规范JTGD622004），取荷载系数为1.2；设计强度对应的混凝土安全系数为1.25,结构工作条件系数为0.95,则要求钢筋混凝土结构的整体安全系数应不小于1.58,计算见下式。由于目前规范为规定结构二类稳定系数，根据上述分析，二类稳

7、定系数可以参考取1.5 8。3分析项目1.20 1.25Z 0.95= 1.58该桥引桥上部结构为40米T梁，先简直后连续，下部结构墩型有矩形实腹 墩、矩形空心墩和双柱式墩，墩身基础为桩基础。王桥为 85+160+85m预应力混 凝土连续刚构桥，主桥桥墩为双柱式空心薄壁墩，最高墩高为83.423m，基础为桩基础。首先对结构进行特征值分析，得到的最小特征值作为第一类稳定安全系数， 并将该值作为后继非线性分析的参考荷载因子； 第一阶屈曲变形要作为初始缺陷 加入原结构；同时考虑结构的材料非线性、几何非线性及初始缺陷。定义材料的 本构模型，打开大变形开关，将一阶特征值屈曲变形作为结构的初始缺陷加入原

8、结构，采用逐步加载的方式求解结构的极限荷载。进行后处理分析，得到结构的二类稳定安全系数。3.1实心矩形桥墩根据设计单位提供的一跨40m梁恒+活荷载为14435 kN。将该荷载作于与桥 墩墩顶，进行稳定分析。(1)计算墩高40米，截面尺寸6.0X2.0;分析结果：桥墩一类稳定为10.4, 二类稳定为6.8，失稳模态为横桥向侧倾失稳图4有限元模型图5求解历程曲线NODAL SOLUTIONSTEP=1SUB =12TIME=6.88UZ (AVG)RSYS=0DMX =.303629SMN =-.159E-03SMX =.302476DEC 7 200911:37:04.067093.134345

9、.201597.26885.033467.100719.167971.235223.302476图6墩身位移图1ELEMENT SOLUTIONSTEP=1SUB =12TIME=6.88EPTOX (NOAVG)RSYS=0DMX =.303629SMN =-.773E-03SMX =.135E-03DEC 7 200911:34:31-.773E-03-.672E-03III卩MX-.571E-03-.369E-03-.168E-03.341E-04-.470E-03-.268E-03-.668E-04.135E-03图7墩身应变图8.2,(2)计算墩高50米，截面尺寸6.0X2.2，分析

10、结果：桥墩一类稳定为二类稳定为6.3，失稳模态为横桥向侧倾失稳墩顶位移（m）图8求解历程曲线1NODAL SOLUTIONSTEP=1SUB =14TIME=6.343UZ (AVG)RSYS=0DMX =.363807SMN =-.151E-03SMX =.362739DEC 7 200911:48:30-.151E-03.04017.080491.120812MX.161133.241776.201455.322418.282097.362739图9墩身位移图ELEMENT SOLUTIONSTEP=1SUB =14TIME=6.343EPTOX (NOAVG)RSYS=0DMX =.36

11、3807SMN =-.661E-03SMX =.112E-03DEC 7 200911:49:16z JMjMX-.661E-03-.575E-03-.489E-03-.403E-03-.318E-03-.146E-03-.232E-03.257E-04-.602E-04.112E-03图10墩身应变图(3) 计算墩高63米，截面尺寸6.0X2.5。分析结果：桥墩一类稳定为6.8, 二类稳定为5.6，失稳模态为横桥向侧倾失稳。图11求解历程曲线1NODAL SOLUTIONSTEP=1SUB =15TIME=5.632UZ (AVG)RSYS=0DMX =.299997SMN =-.138E-

12、03SMX =.299059ANDEC 7 200913:45:41-.138E-03.0331061ELEMENT SOLUTIONSTEP=1SUB =15TIME=5.632EPTOX (NOAVG)RSYS=0DMX =.299997SMN =-.451E-03SMX =.970E-05MX.06635.132838.199326.099594.166082.265814.23257.299059图12墩身位移图AIMDEC 7 200913:46:15-.451E-03-.349E-03-.246E-03-.144E-03-.400E-03-.298E-03-.195E-03-.41

13、5E-04 -.927E-04.970E-05图13墩身应变图(4) 结论根据稳定评价标准，该类桥墩满足稳定要求3.2空心矩形桥墩(1)计算墩高60米，截面尺寸6.0X3.0m，壁厚0.55m。分析结果：桥墩一1DEC 7200914:00:43图14有限元模型类稳定elEmen12.1，二类稳定为7.9，失稳模态为横桥向侧倾失稳1ELEMENTSDEC 7 200914:02:08图15局部有限元模型DEC 7 2009MX14:10:150.000.100.200.300.400.50图16求解历程曲线iNODAL SOLUTIONSTEP=1SUB =12TIME=7.975UZ (AV

14、G)RSYS=0DMX =.403955SMN =-.402E-03SMX =.401005-.402E-03.044199.0888.133401.178001.267203.222602.311804.401005ELEMENT SOLUTIONSTEP=1SUB =12TIME=7.975EPTOX (NOAVG)RSYS=0DMX =.403955SMN =-.943E-03SMX =-.128E-03zHMXDEC 7 200914:10:40-.943E-03-.762E-03-.580E-03-.399E-03-.218E-03-.852E-03-.671E-03-.490E-0

15、3-.309E-03-.128E-03图18墩身应变图(2)计算墩高70米，截面尺寸6.0X3.0m，壁厚0.55m。分析结果：桥墩一 类稳定为8.6，二类稳定为6.6，失稳模态为横桥向侧倾失稳。图19求解历程曲线MXDEC 7 200914:18:20NODAL SOLUTIONSTEP=1SUB =14TIME=6.606UZ (AVG)RSYS=0DMX =.515908SMN =-.339E-03SMX =.513484.5134841ELEMENT SOLUTIONSTEP=1SUB =14TIME=6.606EPTOX (NOAVG)RSYS=0DMX =.515908SMN =-

16、.835E-03SMX =-.665E-04ANDEC 7 200914:19:09ZMX.835E-03-.665E-03-.494E-03-.323E-03-.152E-03-.750E-03-.579E-03-.408E-03-.237E-03-.665E-04-.339E-03.113844.228027.34221.456393.056753.170936.285119.399301图20墩身位移图图21墩身应变图(3)计算墩高80米，截面尺寸6.0X3.0m，壁厚0.55m。分析结果：桥墩一 类稳定为5.6，二类稳定为4.0。图22求解历程曲线1NODAL SOLUTIONSTEP

17、=1SUB =13 TIME=4.064UZ (AVG)RSYS=0DMX =.436872SMN =-.211E-03SMX =.435352一ANDEC 7 2009-.211E-03.096581.048185MX14:23:30.386956.193372.290164.144977.241768.33856.435352ELEMENT SOLUTIONSTEP=1SUB =13TIME=4.064EPTOX (NOAVG)RSYS=0DMX =.436872SMN =-.509E-03SMX =-.521E-04-.509E-03-.407E-03-.306E-03-.204E-03

18、-.103E-03-.458E-03-.357E-03-.255E-03-.154E-03-.521E-04图24墩身应变图(4)结论根据稳定评价标准，该类桥墩满足稳定要求。3.3双柱式桥墩计算墩高40米，单墩截面尺寸直径2.0m圆截面。分析结果：桥墩一类稳定 为4.29, 二类稳定为3，失稳模态为横桥向侧倾失稳。分析显示，该类桥墩满足 稳定要求。ELEMENTSNOV 29 200911:08:20图25有限元模型数系定稳0.50.000.100.200.300.40墩顶位移（m1NODAL SOLUTIONSTEP=1SUB =13TIME=3.068UZ (AVG)RSYS=0DMX =

19、.488927SMN =-.148E-03SMX =.488129ANDEC 72009 14:35:32-.148E-03.108358.054105.162611.216864.32537.271117.433876.379623.488129图27墩顶位移图1ELEMENT SOLUTIONSTEP=1SUB =13TIME=3.068EPTOX (NOAVG)RSYS=0DMX =.488927SMN =-.001016SMX =.425E-03MXANDEC 7 200914:33:55.265E-03.425E-03-.001016-.696E-03-.376E-03-.552E-

20、04-.856E-03-.536E-03-.215E-03.105E-03图28墩身应变图3.4主桥桥墩(1)计算墩高84米，单墩截面尺寸直径6.0X3.0截面，壁厚0.7m。该桥采用悬臂施工法，稳定最不利阶段为边跨合拢前一阶段，计算的该阶段梁恒载重64655 kN。将该荷载施加在桥墩上，进行稳定分析。分析结果：桥墩一类稳定为11.1,二类稳定为6.53,失稳模态为横桥向侧倾失稳。分析显示，该类桥墩满足NOV 29 200921:45:04图29有限元模型1ANELEMENTSDEC 7200914:39:21图30局部有限元模型数系定稳0 0.000.200.400.600.80 1.00墩顶位移（m图31求解历程曲线1NODAL SOLUTIONSTEP=1SUB =31TIME=6.525UX (AVG)RSYS=0DMX =.835145SMN =-.001203SMX =.826879DEC 7 200916:14:30-.001203.090806.182815.274824.366833.458842.550851.64286.73487.826879ELEMENT SOLUTIONSTEP=1SUB =31TIME=6.525EPTOX (NOAVG)RSYS=0DMX =.835145SMN =-.001985SMX =.198E-04DE