高度超静定斜拉桥的非线性分析研究.doc

此文档收集于网络，如有侵权，请联系网站删除高度超静定斜拉桥的非线性分析研究1.摘要一个拉索高度超静定的斜拉桥的非线性分析比较在研究中被实行。 包括桥的几何学和预应力分配的初始形状是使用双重迭代的方法决定的,也就是,一个平衡迭代和一个形状迭代。对于开始的形状分析，一个线性和一个非线性计算程序被建立。以前斜拉桥所有非线性被忽视，而且形状迭代是不考虑平衡而实行的。后来桥的所有非线性被考虑到，而且平衡和形状的重复都实行了。基于收敛于一点的起始形状由不同的程序决定,自振频率和震动模态也被详细地研究。数字的结果表明收敛于一点的起始形状能由二个环的重复方法快速地得到,合理的起始形状能由线性的计算程序决定，而且那样许多计算工作将被节省。在由线性的和非线性计算程序决定的结果之间的几何学和预应力分配中只有很小的不同。然而，对于自振频率和震动模态的分析来说，基本的频率和震动模态将会有显著的不同，而且斜拉桥反应的非线性只出现在由非线性计算得到的初始形状的基础之上的模态中。2.序言在过去的三十年中斜拉桥分析和建筑中取得了飞速的进步。进步主要是由于计算机技术的领域发展,高强度的钢拉索,正交异性钢板和建筑技术产生的。既然第一座现代的斜拉桥1955年在瑞典被建造，他们的名声在全世界得到快速地增长。因为它的直立美学的外观，经济原因和便于直立，斜拉桥被认为是跨径范围从200m到大约1000m的最合适的建筑类型。 世界上现在最长的斜拉桥是日本的横跨岛海、连接本州四国的多多罗桥。多多罗斜拉桥在1999年5月1日被开通，它有890m的一个中央跨径和1480m的总跨度。一座斜拉桥由三个主要的成分所组成，也就是主梁、索塔和斜拉索。主梁在沿纵向方向由拉索弹性支撑以使主梁能跨越一个更长的距离而不需要中间桥墩。主梁的永久荷载和车辆荷载通过拉索传递给索塔。很大的拉力存在于拉索中减小了索塔中大部分和梁的一部分压力。 斜拉桥的非线性的来源主要地包括拉索下垂,梁柱的偏压和大的偏转效应。因为在未施加活载前拉索中存在高度预应力，斜拉桥的初始形状和预应力由每条拉索决定。他们不能够被独立地看成是传统的钢或者是高强混凝土桥。因此开始的形状必须被在桥的分析之前正确的决定。只有基于正确的起始形状才能得到一个正确的偏转和震动分析。这篇论文的目的要提供一个高度冗余的斜拉桥的非线性分析的比较，桥的开始形状将会由线性和非线性计算程序迭代来决定。基于开始的形状计算,桥的震动频率和模态被确定。3. 系统方程3.1 一般的系统方程当只有非线性在刚体中被考虑到，而且系统的衰减矩阵被认为是恒定的时候,在非线性动力学中结构的一个有限元模型才能从虚工作原则中得到，如下:Kjbj-Sjaj= Mq”+ Dq3.2 线性化系统方程为了要不断的解决更大的偏转问题,线性化系统的方程必需用到。通过泰勒的一般方程的扩展的最早的条目，对于一个小的时间（或荷载）间隔的线性化的方程便得到，如下： Mq”+Dq +2Kq=p- up3.3 在静力学中的线性化系统方程在非线性静力学中，线性化系统方程变成：2Kq=p- up4非线性分析4.1. 起始形状分析 斜拉桥的初始形状提供了几何学的结构和桥在主梁和索塔的恒载、斜拉索的拉力作用下的预应力分配。作用的平衡条件，指定的边界条件和建筑的设计需求应该被满足。因为计算的形状,主梁和索塔的永久荷载必须被考虑，拉索的自重被疏忽，而且拉索下垂的非线性应包括在内。形状的计算通过使用二重迭代的方法运行,也就是,平衡重复和形状重复循环。这能用拉索中的任意小的张力开始。基于参考结构 (建筑设计形式)，没有歪斜和零的预应力在主梁和索塔中,斜拉桥平衡位置在恒载作用下是由迭代首先确定的(平衡迭代)。虽然首先决定结构的是使平衡情况和边界情况得到满足，但是建筑的设计需求大体上没有得到实现。因为桥的跨径是很大的而没有预应力存在斜拉索中,相当大的偏转和非常大的弯矩可能在主梁和索塔中出现。那么另外的一个迭代有必要执行来减少偏转和使主梁的弯矩平滑并最后找出正确的初始形状。如此的一个迭代程序在这里命名为形状迭代。对于形状迭代,在先前步骤中确定的基本的轴线力将会被作为下个重复采取的初始基本力,这样一个新的平衡结构在恒载和这个初始力下再次被确定。在形状迭代的时候，一些控制点(主梁和拉索连接的点)将会被选择检验应力集中。在每次形状迭代过程中，主跨的控制点的垂直位移比率将会被检验。也就是,形状迭代将会重复直到应变可以达到所说的10-4。当应变达到的时候, 计算将会停止而斜拉桥的初始形状就找到了。数字的实验表明重复收敛于一点是没什么作用的，并且所有的三个非线性对最后的几何初始形状有比较少的影响。只有拉索下垂作用在确定初始形状分析中有显著作用，而偏压柱和大的偏转效应变则无关重要。 开始的分析能以二种不同的方式被实行: 一个线性和一个非线性计算程序。(1) 线性的计算程序:为了要找到桥的平衡结构,斜拉桥的所有非线性因素被疏忽，而只是线性的弹性拉索、梁单元、同等的线形的变形系数被使用。形状迭代是不考虑平衡迭代而实行的。合理的收敛于一点的起始形状被得到，而且许多计算的工作能被节省。(2) 非线性计算程序:斜拉桥所有的非线性因素在整个的计算程序中被考虑。非线性拉索元素的下沉作用、主梁元素的稳定系数和非线性变形调整系数被应用。形状的迭代和平衡迭代都在非线性计算中实行。 牛顿-瑞普生方法在这里被用于平衡迭代。4.2 静态偏转分析基于确定的起始形状,斜拉桥在活载作用下的非线性静态偏转分析可通过模数或迭代运行。 荷载模数方法导致很大的数字错误是广为人知的。迭代方法比较适于非线性计算，而且需要的应变应能被达到。牛顿- 瑞普生的迭代程序将被使用。因为非线性分析较大或复杂的结构系统,一个完整的迭代程序(重复为一个单一全部荷载运行步骤)将会时常失败。一个模数-迭代程序高度地被推荐,荷载将会被增加，而且重复将会在每个荷载步骤中实行。斜拉桥的静态偏转分析将会从使用线性或非线性计算程序决定的初始形状开始。斜拉桥静态的偏转分析的运算法则在第4.4.2节中被概述。4.3. 线性振动分析当一个结构系统是足够稳固而且外部的刺激不是太强烈,系统可能以一个确定的非线性的静态系数作一个小振幅振动，由振动引起的非线性静态系数的变化是很小的和可以忽略的。这种以一个非线性静态系数以一个小振幅的振动被称作线性化振动。线性化振动不同于线性振动，系统用很小的振幅以一个线性静态系数振动。非线性静态系数qa 能由非线性偏转分析决定。在决定qa之后，系统矩阵可能被建立有关于如此的一个非线性静态系数，线性化系统的等式如下所示：MAq”+ DAq+ 2KAq=p(t)- TA上面的上标字母A 代表在非线性静态系数qa 被计算的数量。这个等式用恒定系数矩阵MA、DA、2KA表现第二的次序一组线性的一般差别的等式。这个等式能被模型的重叠方法，整体的变形方法或直接的整合方法解答。当减幅效应和荷载限制被忽略的时候,系统等式变成：MAq” + 2KAq=0 这个等式表现基于非线性静态系数qa 的不减幅的系统天然振动。天然振动的频率和模态可以从上面的等式运用程序，举例来说子空间重复方法来得到。对于斜拉桥,它的起始形状是非线性静态系数qa。 斜拉桥由于以开始的形状为基础小振幅振动的时候,天然的频率和模态能被找到来解决上述的等式。4.4. 斜拉桥计算运算法则分析斜拉桥的形状的确定计算、静态偏转分析和振动分析的计算法则简短的概述如下：4.4.1 起始形状分析1. 桥的几何和实际的数据输入。2. 主梁和索塔的恒载的输入而且适当地估计了起始拉索中的受力。3. 确定平衡位置(i)线性的程序o 运用线性的拉索和主梁刚性单元。o 运用线性的恒定变形调整系数aj。o 建立线性系统刚度矩阵K通过排列元素的刚度矩阵。o 求解线性等式得到q。(平衡位置)o 没有平衡的迭代被实行。(ii)非线性程序o 非线性下垂效应的拉索和主梁单元被使用。o 非线性变形调整系数aj; aj，被使用。o 建立接触的系统刚度矩阵2K。o 求解增量系统的等式以得到q。o 平衡迭代使用牛顿-瑞普生方法运行q。4. 变形迭代。5. 包括几何形状和基本力的初始的形状输出。6. 对于线性偏转分析，只有线性刚度单元和变形系数被采用且没有平衡迭代的实行。4.4.3. 振动分析1. 桥的几何的和物理的数据输入。2. 包括开始的几何和开始单元常受力的初始形状数据的输入。3. 建立以初始形状的自由振动的线性化系统等式。4. 运用子空间重复方法得到振动频率和模态,例如Rutishauser方法。5. 初始拉索受力估算：在王教授和林教授的最近研究中,小型的斜拉桥的通过任意小或任意大试验初始拉索应力来实现。如果不同的初始拉索受力试验评价被采用，在那里重复单调的迭代，而最后的结论有相似的结果。然而对于大型斜拉桥，确定形状的计算变得更困难以达到一致。在非线性分析中，牛顿-瑞普生类型迭代计算能收敛到一点,只有当解决的被估计的价值是在真正的价值附近时才能实现。当斜拉桥分析的形状在王的论文中建议用任意小的初始拉索应力开始时，收敛到一点的困难可能会出现。因此,估计适当的试验开始的拉索应力来得到一致的结论对于形状确定分析变得重要起来。接下来,一些估计试验初始拉索应力的方法将会被讨论。5.1. 垂直荷载的平衡5.2. 零力矩的控制5.3. 零位移的控制5.4. 拉索等价系数比的概念5.5. 不对称的考虑如果估计的初始拉索应力是用上面介绍的方法对每条拉索独立地确定的，非对称的斜拉桥索塔中可能会存在不平衡的水平受力。因为中间跨(主跨)和边跨的对称的拉索布置对索塔的水平分力的合力为零，也就是没有不平衡的水平力。而对于非对称的斜拉桥，拉索在中间跨和边跨分布是不对称的，拉索在索塔上产生的分力分别独立计算，很明显的索塔中的不平衡的水平力将会引起很大的弯矩和偏移。因此，非对称的斜拉桥的这个问题可以按如下解决。中间跨(主要部份)的拉索受力可以通过上面独立介绍的方法确定，其中上标m代表主跨，上标i代表第i条拉索。然后边跨上拉索受力通过与索塔连接的拉索的水平平衡方程确定，即 Tim cosi= Tis cosi, and Tis = Tim cosi/ cosi, i是指斜拉索与主跨梁的夹角，i是指斜拉索与边跨梁的夹角。6. 例子在这项研究中，二座不同的类型小型斜拉桥从文学中取得,而且他们的起始形状将会用先前描述的形状确定方法使用线性和非线性程序来确定。最后，一座高度冗余的斜拉桥将会被研究。对于平衡迭代和形状迭代都采用应变10-4。重复的最大周期定为20。如果重复的循环数超过20则计算被认为是不收敛于一点的。接着的二个小型斜拉桥的开始形状在第6.1节和6.2首先采用任意小的试验拉索应力决定的。在这二个例子中收敛于一点是重复单调的。他们的收敛于一点的起始形状可以很容易地获得。由线性和非线性计算决定的开始形状之间只有很小的不同。收敛于一点的结论显示同样的结果，而且他们与试验的拉索应力无关。7. 结论通过线性和非线性计算二重循环的建立而得到斜拉桥的初始形状。这个方法能达到建筑的设计形式有统一的预应力分配,而且使所有的平衡和边界情况满足。初始的形状确定是斜拉桥分析中最重要的工作。只有一个正确的起始形状，才能得到一个有意义的和正确的偏转及震动分析。基于研究的数字实验,一些结论概述如下:(1).对于小型斜拉桥初始形状的确定不会出现现在的困难，任意的初始试验拉索应力都能用来计算。然而对于大跨度的斜拉桥，循环的收敛于一点会产生很大问题。(2).在跨度的斜拉桥的形状确定的收敛于一点通常会产生困难，当拉索的试验应力通过垂直荷载平衡、零弯矩控制、零位移控制的方法给出时。(3). 如果主跨的每条拉索的预应力符合到约Eeq的80%，通过两重循环的方法可以很快的找到收敛于一点的初始形状，而且边跨拉索的初始应力是由作用于索塔上的水平等式决定的。(4). 由线性程序和非线性的程序得到的初始形状的结果对于几何和预应力的分配的结果只有很小的不同。(5).使用线性的计算能提供一个合理的起始形状而且节省很多的计算工