毕业论文(设计)-《大跨度桥梁实用几何非线性分析》.docVIP

PAGE w 大跨度桥梁实用几何非线性分析 摘要： 本文从简单实用的角度论述了空间杯系结构的几何非线性 分析 理论 。文中分析了非线性有限元 方法 的求解过程，特别强调决定几何非线性收敛结果的关键 问题 ，即由节点位移增量 计算 单元的内力增量。通过引入随转坐标系，论述了平面和空间梁单元小应变变形时单元内力增量的计算问题。用本文方法可以分析大跨度桥梁结构的六位移大旋转问题。并且用实桥算例进行了验证。 关键词： 大跨度桥梁 几何非线性 实用分析 非线性有限元 小应变理论 江阴 长江大桥 一 . 引言． 现代 大跨度桥梁等工程结构的柔性特征已十分明显，对于这些结构考虑几何非线性的 影响 己必不可少。并且，计算机能力的大大提高也使得分析大型复杂结构的非线性问题成为可行。 80 年代国外对几何非线性问题的 发展 已相当完善［ 1 ， 2 ］，国内在这方面也做了不少的工作［ 4 － 6 ］ 在工程结构几何非线性分析中，按照 参考 构形的不同可分为 TL （ Total Lagranrian) 法和 UL （ Updated Lagrangian) 法 [1] 。后来，引入随转坐标系后又分别得出 CR （ Co-rotational)-TL 法和 CR － LU 法 [2,3] ，在工程中 UL （或 CR-UL ）法 应用 较多。以前的 文献 大都对结构的几何刚度矩阵进行了复杂而详细的推导。从文中的分析可以发现，结构几何刚度矩阵的精确与否并不实质性地影响迭代收敛的最终结果，求解几何非线性问题的关键在于如何由节点位移增量准确地计算出单元的内力增量，而这一点以前文献都没有提到过。因此，本文的重点放在论述单元内力增量的计算上。 工程上很早就开始使用拖动坐标系来求解大跨度桥梁结构的大挠度问题，本文则把它应用到单元内力增量的计算中。从实质上说，这里的拖动坐标系与上面提到的随转坐标系没有区别。因此，在理论方法上， 目前 文中的方法可以归类到 CR － UL 法。但由于本文重点不在于详细介绍这种方法的理论体系，所以论述中均不再使用该名词。本文的目的主要是通过简化复杂的几何非线性分析方法，推广该方法在实际工程中的应用。 二、非线性商限元求解过程 对于工程结构的非线性问题，用有限元方法求解时的非线性平衡方程可写成以下的一般形式： Fs( δ） -P0 （δ） =0 （ l ） 其中，为节点的位移向量； Fs( δ ) 为结构的等效节点抗力向量，它随节点位移及单元内力而变化； PO （δ）为外荷载作用的等效节点荷载向量，为方便起见，这里暂时假定它不随节点位移而变化。 由于式（ l ）中的等效节点抗力一般无法用节点位移显式表示，故不可能直接对非线性平衡方程进行求解。但实际结构的整体切向刚度容易得到，所以通常应用 Newton-Raphson 迭代方法求解该问题。结构的整体切向刚度矩阵 KT 可表示如下 dPO ＝ KTd δ （ 2 ） 式中， KT ＝ KE 十 KG ，其中 KE 为结构的整体弹性刚度矩阵， KG 为几何刚度矩阵。 用混合 Newton － Raphson 迭代方法求解结构非线性问题的基本过程如下： (1 ）将等效节点荷载 PO 分成 n 步，Δ P0 ＝ PO ／ n ，计算并组集结构的整体切向刚度矩阵，进入加载步循环； (2) 求解节点位移增量； (3) 计算各单元内力增量，修正单元内力； (4) 更新节点坐标，计算节点不平衡力 R ； (5) 判断节点不平衡力 R 是否小于允许值，如满足条件，则进入下一个加载步；如不满足条件，重新计算结构的整体切向刚度矩阵，用 R 代替Δ P0 ，回到第 2 步； (6 ）全部加载步完成之后，结束。 从上述求解过程中可见，最为关键的一步是第 3 步，即由节点位移增量计算单元的内力增量。也可以说是由这一步决定了最终的收敛结果，以下将对此着重论述。其实结构的整体切向刚度矩阵对结果并无实质性的影响，修正的 NetwRaphson 方法正是利用这一点来节省迭代计算的时间。 以前的文献对空间梁单元几何刚度矩阵的推导方面论述较多，都建立在一些假定的基础上，这里就不详细说明。考虑到结构的整体切向刚度矩阵精确与否并不改变最终结果，仅影响迭代收敛的速度，并且不是越精确的整体切向刚度矩阵迭代收敛越快。 三、小应变时单