桥墩二类稳定分析报告
2009.12.1
1.结构稳定分析理论
按照结构在逐渐加载过程中平衡形式是否发生质变的观点，将结构的失稳区分为第一类稳定问题和第二类稳定问题，如图l －图2所示。前者表示在加载过程中，构件的平衡状态将出现分枝现象，使原有的平衡状态失去稳定性而转向新的平衡状态；而后者在加载过程中平衡形式并不发生质变。在第一类稳定同题中，当加载至P cr 时，表示平衡的分枝即将出现，称P cr 为压屈荷载。在第二类稳定问题中，当加载到P cr 时，表示构件的承载能力即将降低，称为压溃荷载。工程上通常把两者统称为失稳的临界荷载。工程问题中研究结构稳定问题的目的，在于寻求相应的临界荷载及其临界状态，防止不稳定平衡状态的发生，从而确保结构安全。
结构的第一类稳定问题，在数学上归结为广义特征值问题，
()0d g K K λδ+= (1)
式中，d K 为弹性刚度矩阵，g K 为几何刚度矩阵，它只与结构的轴力有关；
δ为结构的位移增量；λ为载荷稳定系数。式(1)数学上表现为广义特征问题，应
用各种迭代方法，如逆矢量迭代法、子空间迭代法等都可以很方便地求解。
结构的几何非线性是指大位移问题。在载荷作用下，当位移大到足以使得结构的几何形状发生显著的改变时，必须按已变形的位置建立平衡方程。弹性大变形问题与变形的历史不相关，可以采用全量方法研究，也就是直接求在已知载荷的约束下的总变形和应力，结构全量形式的非线性平衡方程可表示为：
{}()0K F δδ−= (2)
式中：()K δ为非线性刚度矩阵，δ为节点位移，F 为等效节点载荷。
非线性稳定分析的基本方法是逐步地施加一个恒定的载荷增量直到解变得开始发散为止。由于特征值
屈曲载荷是预期的线性屈曲载荷的上限，故可以作为非线性屈曲的给定载荷，在渐进加载达到此载荷前，非线性求解应该发散。设λ为特征值屈曲分析时的稳定系数，在这里取非线性屈曲分析的给定载荷系数。将给定载荷分为n 个载荷增量，
010t n λλλλλ=<⋅⋅⋅⋅<⋅⋅⋅⋅<= (3)
在每一个载荷步内对非线性方程(2)进行线性化，可得增量形式的平衡方程：
0t K F δΔ−Δ= (4) 0t l K K K K σ=++ (5) 1()t t F F λλ−Δ=− (6)
式中：t K 为切线刚度矩阵，F Δ为等效节点载荷增量矩阵，δΔ为节点位移增量，0K 为小位移线性刚度矩阵，l K 为大位移矩阵，K σ为初应力矩阵，F 为等效节点载荷矩阵。在每一个载荷步内，为了改进求解的精度，可以应用牛顿法进行迭代。结构的极限承载力在开始发散的载荷和在此前一级已收敛的载荷之间。如载荷增量步分得较细，可以偏于安全地认为是前一级载荷，而避免更加复杂的计算。
目前二类稳定分析都是采用有限元等数值计算方法。从本构关系来考虑稳定问题可以分为弹性稳定与弹塑性稳定，而弹性稳定又从是否考虑几何非线性、初始缺陷等因素又分为线弹性稳定和非线性弹性稳定。由于施工环节会存在不可避免的施工误差，桥墩轴线可能出现偏差，高墩在水平荷载作用下将产生较大的变形等等，稳定计算必须计人初始缺陷及大位移的影响，所以基于极值点失稳为理论基础的计人双重非线性的弹塑性稳定问题对于工程是必须的。
2.桥墩ANSYS 分析
2.1 单元模拟为
在ANSYS 软件中，采用beam189梁单元模拟桥墩和横梁，beam189为3-D 二次有限应变梁单元，可以很好的模拟桥墩力学行为。
2.2 材料本构
材料非线性本构关系是进行二类稳定计算的重要依据，尤其是极限阶段混凝土的强化过程。
研究中只考虑了混凝土抗压部分，选用美国学者Hongnestad 提出的本构曲线。具体表达式为：
2000000210.15u εεσεεσεεσεε⎧⎡⎤
⎛⎞⎛⎞⎪⎢⎥−⎜⎟⎜⎟⎪⎢⎥⎝⎠⎝⎠⎪⎣⎦=⎨⎡⎤⎪⎛⎞−−⎢⎥⎪⎜⎟−⎝⎠⎪⎣⎦⎩
0u εεεεε≤<≤ 式中，0σ为棱柱体极限抗压应力， 0ε为对应于0σ的应变，取0.002，u ε通常取0.0038。
图3 混凝土本构模型
2.3 稳定的评价
（1）一类稳定实质上是特征值失稳，计算公式见下式，一般不小于4~5。
（2）第二类稳定实质上是极值点失稳。用极限状态法设计桥梁时，稳定与最终的极限承载力是统一的。因为二类稳定分析的最终状态也是桥梁达到塑性变形而破坏，因此，桥梁结构的二类稳定安全系
数与强度安全系数也是一致的。按照《公路钢筋混凝土及预应力混凝土桥涵设计规范》(JTGD62一2004)，取荷载系数为1.2；设计强度对应的混凝土安全系数为1.25，结构工作条件系数为0.95，
4~5cr
T
P P λ=
≥
则要求钢筋混凝土结构的整体安全系数应不小于1.58，计算见下式。由于目前规范为规定结构二类稳定系数，根据上述分析，二类稳定系数可以参考取1.58。
3.分析项目
该桥引桥上部结构为40米T 梁，先简直后连续，下部结构墩型有矩形实腹墩、矩形空心墩和双柱式墩，墩身基础为桩基础。主桥为85+160+85m 预应力混凝土连续刚构桥，主桥桥墩为双柱式空心薄壁墩，最高墩高为83.423m ，基础为桩基础。
首先对结构进行特征值分析，得到的最小特征值作为第一类稳定安全系数，并将该值作为后继非线性分析的参考荷载因子；第一阶屈曲变形要作为初始缺陷加入原结构；同时考虑结构的材料非线性、几何非线性及初始缺陷。定义材料的本构模型，打开大变形开关，将一阶特征值屈曲变形作为结构的初始缺陷加入原结构，采用逐步加载的方式求解结构的极限荷载。进行后处理分析，得到结构的二类稳定安全系数。
3.1 实心矩形桥墩
根据设计单位提供的一跨40m 梁恒+活荷载为14435 kN 。将该荷载作于与桥墩墩顶，进行稳定分析。
（1）计算墩高40米，截面尺寸6.0X2.0；分析结果：桥墩一类稳定为10.4，二类稳定为6.8，失稳模态为横桥向侧倾失稳。
1.20 1.25
1.58
0.95
λ×≥
=
图4 有限元模型
图5 求解历程曲线
