斜拉桥非线性分析综述.docVIP

斜拉桥非线性分析的研究现状 在结构分析领域，一直普遍存在着由结构几何变形、梁－柱效应和材料塑性引起的结构力－变形非线性关系。Turner等人，在1960年里发表的一篇论文里提出根据结构加载前已存在的应力建立刚度矩阵和在几何非线性分析中使用线性化方法与增量法的概念。1966年，Oden、Prazmieniecki提出了计算几何矩阵的方法。Turner、Oden等人中公式推导中，位移函数采用了简化的表达式，其分析计算实质上限制在大位移、小应变的范围内。从70年代初期，开始研究关于大位移、大应变大结构，将固体力学的分析方法引用到结构几何非线性有限元分析中，当位移与应变较大之后，通常用到修正坐标的方法，即所谓移动坐标法。大跨径斜拉桥是高次超静定结构，即使在正常荷载作用下，往往发生较大位移，结构几何形状发生显著的变化，整个结构由于有限变形而表现出明显的几何非线性行为。归纳起来，斜拉桥的几何非线性来自三个方面：斜拉索的垂度效应；主梁、索塔中轴力与弯矩相互作用而产生的梁－柱效应；大位移产生几何形状改变而引起的非线性效应。斜拉索作为柔性构件，在自重和轴力作用下，呈悬链线形状。其轴向刚度将随垂度的变化而变化，而斜拉索的垂度又取决于索中的拉力，因此斜拉索拉力与变形之间存在明显的非线性关系。对自锚固体系斜拉桥，斜拉索索力使主梁、桥塔等构件处于弯矩和轴力的组合作用下，桥塔和主梁变形过程中，由于横向挠度会使轴力产生附加弯矩，而弯矩又影响轴向刚度的大小，从而影响结构变形，由此产生所谓的梁－柱效应，使整个斜拉桥表现出非线性行为。大跨度斜拉桥的另一特点是由于柔性较大而产生较大的位移，此有限位移会使斜拉桥的几何形状产生较大的变化，从而使结构分析不能仅按未变形的初始几何形状进行，而应当随着位移的变化逐步修正结构的几何形状。此时，结构的几何刚度矩阵是几何变形的函数。因此平衡方程{F}=[K]{δ}=不再是线性关系，线弹性分析中的叠加原理也不再完全适用。斜拉索垂度效应产生的非线性效应随着索自重即水平投影长度增加而增加，随索中拉力减小而减小。1956年H.J.Ernst提出用直杆模拟斜拉索，用等效弹性模量来考虑斜拉索垂度的非线性效应，这就是今天斜拉桥分析中广泛采用的Ernst公式。 我国学者对斜拉桥的几何非线性也进行了广泛对理论分析与试验研究。1982年周上君用等效弹性模量考虑斜拉索垂度，考虑结构大位移效应，用小挠度全量平衡方程进行迭代计算。但他用小挠度平衡方程计算斜拉桥大位移效应，未考虑主梁与桥塔大梁－柱效应。1990年陈德伟引入稳定函数考虑梁单元的梁－柱效应，用 Ernst 公式考虑拉索的垂度效应，用拖动坐标系考虑大位移影响，求解斜拉桥的几何非线性问题。程国庆、潘家英等总结了斜拉桥几何非线性研究等现状，对各种斜拉桥几何非线性分析方法作了评述，指出：（1）等效弹性模量法用直杆单元模拟整根斜拉索，给斜拉桥的分析带来了很大的方便，但是当斜拉索两端节点位移相当大时，等效弹性模量法具有一定的近似性；（2）处理梁－柱效应可采用几何刚度矩阵和稳定函数法，其中稳定函数具有较高的精度，是处理梁柱非线性的经典方法。如将稳定函数法表达的单元刚度矩阵展开，取其一阶近似即可得到几何刚度矩阵；（3）目前已有的斜拉桥非线性计算理论可大致分为切线模量表达的增量求解法和割线模量表达的全量求解法；理论框架可以分为总体拉格朗日描述（T.L）和修正拉格朗日描述（U.L）。但是斜拉桥非线性有限位移理论在有限元格式的建立过程中大都作了不同程度的简化和近似，因此，所得到的计算模型也不尽相同。从现有的各种非线性计算方法存在的差异可以发现，大跨度斜拉桥的非线性计算理论还有待进一步深入探讨，这大致可以归纳为以下三个方面：1．斜拉桥各种非线性单元模式合理性及其精度的研究；2．斜拉桥几何非线性描述参考构型及非线性求解方法的研究；3．斜拉桥有限元离散方法、结构模型化方法对几何非线性分析结果的影响研究。 几何非线性分析理论垂度效应：斜拉索自重垂度引起的拉索拉力与变形之间的非线性关系；大变形效应：大位移产生的结构几何形状变化引起的几何非线性效应；梁-柱效应：由于斜拉索的拉力作用，主梁和索塔不仅承受弯矩而且还将承受巨大的轴向力，在主梁和索塔变形过程中，由于轴向力和弯矩相互影响，而产生所谓的梁-柱效应（P-Δ效应），使整个斜拉桥结构表现出几何非线性行为。 式中: 为考虑垂度影响的拉索换算弹性模量,kPa; 为拉索弹性模量,kPa; 为拉索换算容重,kN/m3; S为拉索长度,m; 为拉索与水平线的夹角; 为拉索应力,kPa。 柱效应和大位移效应。 在斜拉桥成桥阶段,斜拉索的初始拉力使主梁和桥塔存在较大应力,需要考虑单元初内力对单元刚度的影响。大位移对建立结构平衡方程有较大的影响,需要考虑大位移对单元刚度的影响。 1、完全的拉格朗