大跨度桥梁实用几何非线性分析(1).

大跨度桥梁实用几何非线性分析(1)

本文从简单实用的角度论述了空间杯系结构的几何非线性分析理论。文中分析了非线性 有限元方法的求解过程，特别强调决定几何非线性收敛结果的矢键问题，即由节点位移 增量计算单元的内力增量。通过引入随转坐标系，论述了平 面和空间梁单元小应变变形 时单元内力增量的计算问题O用本文方法可以分析大跨度桥梁结构的六位移大旋转问 题。并且用实桥算例进行了验证。

尖键词：大跨度桥梁几何非线性实用分析非线性有限元小应变理论江阴长江大桥 •引言・ 现代大跨度桥梁等工程结构的柔性特征已十分明显，对于这些结构考虑几何非线性的影 响己必不可少。并且，计算机能力的大大提高也使得分析大型复杂结构的非线性问题

成 为可行。80年代国外对几何非线性问题的发展已相当完善〔1，2〕，国内在这方面也 做了不少的工作〔4・6〕

在工程结构几何非线性分析中，按照参考构形的不同可分为TL (Total Lagranrian)法和 UL ( Updated Lagrangian)法［1 ］。后来，引入随转坐标系后又分别得出CR (Co- rotational)-TL 法和 CR- LU 法［2,3］，在工程中 UL(或

CR・UL法应用较多。以前的文献大都对结构的几何刚度矩阵进行了复杂而详细的推导。

从文中的分析可以发现，结构几何刚度矩阵的精确与否并不实质性地影响迭代收敛的最 终结果，求解几何非线性问题的尖键在于如何由节点位移增量准确地计算出单元的内力 增量，而这一点以前文献都没有提到过。因此，本文的重点放在论述单元内力增量的计 算上。

工程上很早就开始使用拖动坐标系来求解大跨度桥梁结构的大挠度问题，本文则把它应 用到单元内力增量的计算中o从实质上说，这里的拖动坐标系与上面提到的随转坐标系 没有区别。因此，在理论方法上，目前文中的方法可以归类到CR・UL法。但由于本文 重点不在于详细介绍这种方法的理论体系，所以论述中均不再使用该名词。本文的目的 主是通过简化复杂的几何非线性分析方法，推广该方法在实际工程中的应用。

二、非线性商限元求解过程对于工程结构的非线性问题，用有限元方法求解时的非线性 平衡方程可写成以下的一般形式：

Fs( S) -P0 (S) =0 (I)其中，为节点的位移向量；Fs( S )为结构的等效节点抗力向量，它 随

节点位移及单元内力而变化；PO(S)为外荷载作用的等效节点荷载向量，为方便起 见，这里暂时假定它不随节点位移而变化。 由于式(I)中的等效节点抗力一般无法用节点位移显式表示，故不可能直接对非线性平 衡方程进行求解。但实际结构的整体切向刚度容易得到，所以通常应用Newtom

Raphson迭代方法求解该问题。结构的整体切向刚度矩阵 如下

dP8 KTdS (2)式中，KT二KE十KG其中KE为结构的整体弹性刚度矩阵，KG为几何 刚

度矩阵。

用混合Newton- Raphson迭代方法求解结构非线性问题的基本过程如下：(1 )将等效节 点荷载PO分成n步，△PO=PO/n,计算并组集结构的整体切向刚度矩阵，进入加载步循 环；

(2) 求解节点位移增量；

(3) 计算各单元内力增量，修正单元內力；

(4) 更新节点坐标，计算节点不平衡力R ；

(5) 判断节点不平衡力R是否小于允许值，如满足条件，则进入下一个加载步；如不满足

条件，重新计算结构的整体切向刚度矩阵，用 R代替△ P0,回到第2 步；

(6 )全部加载步完成之后，结束。

从上述求解过程中可见，最为尖键的一步罡匸 是第3步，即由节点位移增量计算单 元的内力增量。也可以说是由这一步决定 了最终的收敛结果，以下将对此着重论 述。其实结构的整体切向刚度矩阵对结果并无实质性的影响，修正的NetwRa phson方 法正是利用这一点来节省迭代计算的时间o

以前的文献对空间梁单元几何刚度矩阵的推导方面论述较多，都建立在一些假定的基础 上，这里就不详细说明。考虑到结构的整体切向刚度矩阵精确与否并不改变最终结果， 仅影响迭代收敛的速度，并且不是越精确的整体切向刚度矩阵迭代收敛越快。

三、小应变时单元内力增百计算在一般情况下，工程结构的几何非线性都属于小应变大 位移(大平移、大转动)问题。对于这类问题，单元内力增量的计算比较简单。平面梁单 元是空间梁单元发展的基础，故这里先分析平面梁单元的情况。

平面梁单元在整体坐标系(OXY下从t到t十时刻的变形情况。定义随转坐标系 (oxy)的原点固定在单元的一端(i端)，x轴始终保持沿i Tj的直线方向。可见，在 随转坐

标系中平面梁单元的自由度减少为三个( uxO i 0 j)。

从随转坐标系中的三个自由度可以看出，它们反映的是单元的真实变形情况，与单元所 经历的刚性位移无矢。在用有限元方法求解非线性问题时，只将单元尺寸划分得适当 小，整

KT可表示 体坐标系下的小应变大位移问题在单元随转坐标系中就转化为小应变小位移问 题，这一点可从非线性连续介质力学给出证明。这样，随转坐标系下的受力变形情况就 可近似地接线性处理，单元内力增量的计算也就与线性情况一样，这里不再赘述。同时 也正说明了工程中常用拖动坐标法计算平面结构大变形问题的正确性o

四、算例分析结合以上论述，编制了相应的非线性有限元计算程序。为验证本文方法和 有限元程序，下面首先分析了 45度弯梁空间弯扭大位移问题。大跨度悬索桥在施工阶 段的几何非线性比较明显，因此，必须准确地考虑，否则计算结果可能不正确。作为实 桥算例，对江阴长江大桥在20%拼装率施工阶段的几何非线性问题进行了分析，并与 Ansys程序的计算结果相比较。

1・45度弯梁空间弯扭六位移分析 本例是ADINA中的45度弯梁大位移分析考题。该梁

位于X-y平面内，梁根固定，在自由端沿Z方向受一个集中荷载的作用，梁因此发生 空间弯扭大变形。 分析时将梁划分为8个单元，每步加载量为10.0 o分别用ADINA AnsyS和本文程序计 算了 60个加载步，各计算结果均基本上一致。梁自由端无量纲位置坐标在初始时刻， 加载30步与加载60步时的比较列于表1，可见三者相互较吻合。为了进行对比，都 没有考虑

剪切影响o为简洁起见，这里不指定专门的量纲单位。

2・江阴长江大桥非线性分析

江阴长江大桥是我国目前建成的最大跨度悬索桥。主跨跨度为 1385m主梁为 宽36.9m,M 3.0m的扁平状闭口钢箱梁。主缆相距32.5米，吊杆间距为1.6m，矢跨比 为I /10.5。桥塔为门式框架结构，南北桥塔高分别为 187m和

184m桥面波置为R二27710m的竖曲线。

摘本文从简单实用的角度论述了空间杯系结构的几何非线性分析理论。文中分析了非线 性有限元方法的求解过程

本篇论文是由3C0M文档频道的网友为您在网络上收集整理饼投稿至本站 的，论文 版权属原作者，请不用于商业用途或者抄袭，仅供参考学习之用，否者后果自负，如果 此文侵犯您的合法权益，请联系我们。

根据设计资料，建立了江阴长江大桥的计算模型。在成桥状态下，单根主缆的水平内力 约为23878t，单根吊杆的内力约为144t，考虑到悬索桥在施工时主缆与塔顶有相对位 移，计算模型中主缆与塔顶在顺桥向可自由移动，而其他方向均耦联o

悬索桥施工过程中分段安装主梁，小拼装率时各主梁段之间相互饺接。由于悬索桥在成 桥状态的位置和内力一般为已知，施工状态均从成桥状态通过拆除梁段的方法确定o江 阴长江大桥在成桥状态拆除两端梁段后，但未发生变形之前20%拼装率的初始状态。由 于该初始状态的节点位置和单元内力均为己知，用以上的非线性有限元程序可得出20% 拼装率变形后的施工平衡状态。在变形后的施工平衡状态下，跨中梁段随主缆发生了较 大的烧曲，主跨两端的主缆形状比成桥状态时变化较明显。跨中竖向向下的位移为 2.595m,塔顶处的主缆发生 向外0,669m的位移。为验证该程序的计算结果，对上述同样的工况用 AnsyS程 序进行了分析。计算结果中跨中向下的位移为2.581m塔顶处主缆向外的位移0.664m。 从本文方法的变形和内力结果与Ansys程序计算结果的比较发现，两者均较吻合，这就 验证了本文方法和非线性有限元程序的可靠性和有效性。

五、结语以上从简单实用的角度论述了空间杯系结构的几何非线性分析理论。通过对有 限元求解几何非线性问题过程的分析，特别强调了用选代方法求解杯系结构几何非线性 问题中的尖键问题，即由节点的位移增量计算单元内力增量的重性。在引人随转坐标系 之后，论述了小应变问题中单元内力增量的计算。从论述中可知，随转坐标系下的受力 变形情况可近似地接线性处理，单元内力增量的计算也与线性情况一样，同时也说明了 工程中常用拖动坐标法计算大跨度桥梁结构大变形问题的正确性。

对空间杆系结构用数值算例对本文方法进行了验证。为保证分析结果的正确 性，用多个 程序进行相互校核。

对江阴长江大桥在20%拼装率施工阶段的几何非线性问题进行了分析，分析结果与 An syS程序的计算结果吻合。从分析中可见，在小拼装率施工阶段，悬索桥跨中梁段随 主

缆发生了较大的挠曲，主跨两端的主缆形状比成桥状态时变化比较明显。