大跨度钢管混凝土拱桥非线性整体稳定性研究.docVIP

大跨度钢管混凝土拱桥非线性整体稳定性研究 　　摘要：本文从非线性有限元的基本原理出发，以有限元软件Ansys为平台，采用计算模型详尽地对某桥施工阶段及运营阶段的线弹性及弹塑性稳定问题进行分析和计算，得到了各种工况下第一类稳定安全系数、失稳模态及第二类稳定安全系数，并比较了多种非线性因素下结构的稳定安全系数及各类非线性问题对稳定性的影响。 　　关键词：钢管混凝土拱桥;整体稳定;非线性有限元;初始缺陷;安全系数 　　1 引言 　　对于钢管混凝土拱桥结构稳定分析，目前都是采用有限元等数值计算方法。从本构关系来考虑稳定问题可以分为弹性稳定与弹塑性稳定，而弹性稳定又从是否考虑几何非线性、初始缺陷等因素又分为线弹性稳定和非线性弹性稳定。线弹性屈曲分析假设结构失稳状态为弹性小变形，结构内力与外荷载成比例关系，结构的稳定分析就转化为求特征值问题，求得最小特征值即为失稳临界荷载。线性屈曲分析计算简便、概念清晰，但其理论基础是分支点稳定理论，只适用于理想结构。由于施工环节会存在不可避免的施工误差，最后成形的拱轴线与设计的理想轴线会有偏离;此外，拱肋在结构自重及外荷载作用下，将产生较大的变形，稳定计算必须计入初始缺陷及大位移的影响，所以基于极值点失稳为理论基础的计入双重非线性的弹塑性稳定问题越来越来重要[1]。 　　本文对拱肋采用统一理论模型进行模拟，通过求解结构从加载开始到失稳全过程的结构响应，得到全过程荷载位移曲线[2]，从而探讨几何非线性和材料非线性对整体稳定性影响。 　　2 稳定理论及基于ANSYS的应用 　　2.1第一类稳定分析 　　第一类稳定可归结为如下特征方程： 　　求解时，先对结构施加一个参考荷载，求出对应的几何刚度矩阵，然后代入（1）式，求解广义特征值，解出最小特征值，即可得出临界荷载。且令为第一类稳定问题的稳定安全系数[3]。在基于ANSYS进行线弹性屈曲分析中有以下几点需要注意： 　　（1）线弹性屈曲稳定分析前要先进行线弹性静力分析，在此过程中必须要打开预应力效应开关，因为这样才能计入参考荷载所对应的几何刚度矩阵。其对应的相应命令为：Pstres。 　　（2）第一类稳定问题在数学方法上可以化解成矩阵特征值的问题。对于求高阶矩阵特征值，主要采用子空间迭代法（Subspace Method）和兰索斯分块法（Block Lanczos）[4]。 　　（3）特征值对所有的荷载都作相应的缩放。如果某些荷载是常数，例如结构的自重，而其它荷载是可变的。则在分析中必须要确保从常数荷载得到的刚度在特征值求解时不被缩放。为了达到这一目的通常调整可变荷载，直到特征值为1.0或者接近1.0。 　　（4）通常我们需要观察结构的屈曲模态形状，所以无论采用何种方法提取特征值都必须对解作展开。 　　2.2第二类稳定分析 　　极值点失稳问题的本质在于求解结构在受荷载全过程中荷载-位移曲线。通常我们认为求荷载-位移曲线可用采用两种方法：荷载增量法与位移增量法[3]。本文仅探讨荷载增量法求荷载-位移曲线在ANSYS软件中具体实现需注意的问题[4]： 　　（1） 线弹性屈曲分析得到的是屈曲荷载的上限，所以通常采用此值作为非线性屈曲分析给定的荷载，有时也将此值提高10%～20%。 　　（2） 由于无法跟踪到理想结构屈曲后的路径。通常采用一阶屈曲模态作为初始缺陷，其对应的命令为：Upgeom[4]。 　　（3）一个非收敛的解，并不一定意味结构达到其最大载荷，它也可能是由于数值不稳定引起的，通常采用弧长法扩展到后屈曲的范围，以此来穿过那些发生了“阶跃”响应的区域。 　　3 某大跨度拱桥非线性稳定计算与分析 　　3.1 工程概况 　　某钢管混凝土拱桥为一孔净跨径为280m的下承式钢管混凝土系杆拱，主拱轴线为悬链线，拱轴系数m=1.543，净矢跨比为1/5。拱肋为等截面钢管混凝土桁架结构，由四根100012mm的钢管内包混凝土组成，高5.5m，宽2.4m，拱脚为外包混凝土部分高6.5m，宽3m的实心截面。两拱肋间设11道空间格构横撑，确保拱肋横向稳定。系杆采用40束钢绞线，其中OVM15-19锚固的体外索36根，OVM15-13锚固的体外索4根。 　　3.2 材料本构关系的选取 　　（1）钢材应力―应变关系模型 　　钢材应力―应变曲线采用理想弹塑性模型 　　式中，为钢材的应力;为钢材的弹性模量;为钢材的应变;为钢材的屈服应力;为钢材的屈服应变。 　　（2） 钢管混凝土的应力―应变关系模型[5] 　　根据统一理论模型并假定弹塑性阶段应力―应变曲线为直线，由此得到约束混凝土的轴向应力―应变三段折线模型。 　　式中，为钢管混凝土组合材料的轴向应力比例极限;为钢管混凝土组合材料的弹性变形比例极限;为钢管混凝土组合材