一.引言.
365TRANBBS交通站365jt.com
365交通站365jt.com现代大跨度桥梁等工程结构的柔性特征已十分明显,对于这些结构考虑几何非线性的影响己必不可少。并且,计算机能力的大大提高也使得分析大型复杂结构的非线性问题成为可行。80年代国外对几何非线性问题的发展已相当完善[1,2],国内在这方面也做了不少的工作[4-6]
365交通站365jt.com
365交通站365jt.com在工程结构几何非线性分析中,按照参考构形的不同可分为TL(TotalLagranrian)法和UL(Updated Lagrangian)法[1]。后来,引入随转坐标系后又分别得出 CR(Co-rotational)-TL法和CR-LU法[2,3],在工程中UL(或CR-UL)法应用较多。以前的文献大都对结构的几何刚度矩阵进行了复杂而详细的推导。从文中的分析可以发现,结构几何刚度矩阵的精确与否并不实质性地影响迭代收敛的最终结果,求解几何非线性问题的关键在于如何由节点位移增量准确地计算出单元的内力增量,而这一点以前文献都没有提到过。因此,本文的重点放在论述单元内力增量的计算上。 
365交通站365jt.com
365交通站365jt.com工程上很早就开始使用拖动坐标系来求解大跨度桥梁结构的大挠度问题,本文则把它应用到单元内力增量的计算中。从实质上说,这里的拖动坐标系与上面提到的随转坐标系没有区别。因此,在理论方法上,目前文中的方法可以归类到CR-UL法。但由于本文重点不在于详细介绍这种方法的理论体系,所以论述中均不再使用该名词。本文的目的主要是通过简化复杂的几何非线性分析方法,推广该方法在实际工程中的应用。 
365交通站365jt.com
365交通站365jt.com二、非线性商限元求解过程
365交通站365jt.com
365交通站365jt.com对于工程结构的非线性问题,用有限元方法求解时的非线性平衡方程可写成以下的一般形式:
365交通站365jt.com
365交通站365jt.comFs(δ)-P0(δ)=0 (l)
365交通站365jt.com
365交通站365jt.com其中,为节点的位移向量;Fs(δ)为结构的等效节点抗力向量,它随节点位移及单元内力而变化;PO(δ)为外荷载作用的等效节点荷载向量,为方便起见,这里暂时假定它不随节点位移而变化。
365交通站365jt.com
365交通站365jt.com由于式(l)中的等效节点抗力一般无法用节点位移显式表示,故不可能直接对非线性平衡方程进行求解。但实际结构的整体切向刚度容易得到,所以通常应用Newton-Raphson迭代方法求解该问题。结构的整体切向刚度矩阵KT可表示如下
365交通站365jt.com
365交通站365jt.comdPO= KTdδ (2)
365交通站365jt.com
365交通站365jt.com式中,KT= KE十KG,其中KE为结构的整体弹性刚度矩阵,KG为几何刚度矩阵。
365交通站365jt.com
[1][2][3][4]下一页
[NextPage]
$False$