《高度超静定斜拉桥的非线对性分析研究.doc》由会员分享，可在线阅读，更多相关《高度超静定斜拉桥的非线对性分析研究.doc（10页珍藏版）》请在三一文库上搜索。
1、高度超静定斜拉桥的非线性分析研究1.摘要一个拉索高度超静定的斜拉桥的非线性分析比较在研究中被实行。 包括桥的几何学和预应力分配的初始形状是使用双重迭代的方法决定的,也就是,一个平衡迭代和一个形状迭代。对于开始的形状分析，一个线性和一个非线性计算程序被建立。以前斜拉桥所有非线性被忽视，而且形状迭代是不考虑平衡而实行的。后来桥的所有非线性被考虑到，而且平衡和形状的重复都实行了。基于收敛于一点的起始形状由不同的程序决定,自振频率和震动模态也被详细地研究。数字的结果表明收敛于一点的起始形状能由二个环的重复方法快速地得到,合理的起始形状能由线性的计算程序决定，而且那样许多计算工作将被节省。在由线性的和非
2、线性计算程序决定的结果之间的几何学和预应力分配中只有很小的不同。然而，对于自振频率和震动模态的分析来说，基本的频率和震动模态将会有显著的不同，而且斜拉桥反应的非线性只出现在由非线性计算得到的初始形状的基础之上的模态中。2.序言在过去的三十年中斜拉桥分析和建筑中取得了飞速的进步。进步主要是由于计算机技术的领域发展,高强度的钢拉索,正交异性钢板和建筑技术产生的。既然第一座现代的斜拉桥1955年在瑞典被建造，他们的名声在全世界得到快速地增长。因为它的直立美学的外观，经济原因和便于直立，斜拉桥被认为是跨径范围从200m到大约1000m的最合适的建筑类型。 世界上现在最长的斜拉桥是日本的横跨岛海、连接本
3、州四国的多多罗桥。多多罗斜拉桥在1999年5月1日被开通，它有890m的一个中央跨径和1480m的总跨度。一座斜拉桥由三个主要的成分所组成，也就是主梁、索塔和斜拉索。主梁在沿纵向方向由拉索弹性支撑以使主梁能跨越一个更长的距离而不需要中间桥墩。主梁的永久荷载和车辆荷载通过拉索传递给索塔。很大的拉力存在于拉索中减小了索塔中大部分和梁的一部分压力。 斜拉桥的非线性的来源主要地包括拉索下垂,梁柱的偏压和大的偏转效应。因为在未施加活载前拉索中存在高度预应力，斜拉桥的初始形状和预应力由每条拉索决定。他们不能够被独立地看成是传统的钢或者是高强混凝土桥。因此开始的形状必须被在桥的分析之前正确的决定。只有基于正
4、确的起始形状才能得到一个正确的偏转和震动分析。这篇论文的目的要提供一个高度冗余的斜拉桥的非线性分析的比较，桥的开始形状将会由线性和非线性计算程序迭代来决定。基于开始的形状计算,桥的震动频率和模态被确定。3. 系统方程3.1 一般的系统方程当只有非线性在刚体中被考虑到，而且系统的衰减矩阵被认为是恒定的时候,在非线性动力学中结构的一个有限元模型才能从虚工作原则中得到，如下:Kjbj-Sjaj= Mq”+ Dq3.2 线性化系统方程为了要不断的解决更大的偏转问题,线性化系统的方程必需用到。通过泰勒的一般方程的扩展的最早的条目，对于一个小的时间（或荷载）间隔的线性化的方程便得到，如下： Mq”+Dq
5、+2Kq=p- up3.3 在静力学中的线性化系统方程在非线性静力学中，线性化系统方程变成：2Kq=p- up4非线性分析4.1. 起始形状分析 斜拉桥的初始形状提供了几何学的结构和桥在主梁和索塔的恒载、斜拉索的拉力作用下的预应力分配。作用的平衡条件，指定的边界条件和建筑的设计需求应该被满足。因为计算的形状,主梁和索塔的永久荷载必须被考虑，拉索的自重被疏忽，而且拉索下垂的非线性应包括在内。形状的计算通过使用二重迭代的方法运行,也就是,平衡重复和形状重复循环。这能用拉索中的任意小的张力开始。基于参考结构 (建筑设计形式)，没有歪斜和零的预应力在主梁和索塔中,斜拉桥平衡位置在恒载作用下是由迭代首先
6、确定的(平衡迭代)。虽然首先决定结构的是使平衡情况和边界情况得到满足，但是建筑的设计需求大体上没有得到实现。因为桥的跨径是很大的而没有预应力存在斜拉索中,相当大的偏转和非常大的弯矩可能在主梁和索塔中出现。那么另外的一个迭代有必要执行来减少偏转和使主梁的弯矩平滑并最后找出正确的初始形状。如此的一个迭代程序在这里命名为形状迭代。对于形状迭代,在先前步骤中确定的基本的轴线力将会被作为下个重复采取的初始基本力,这样一个新的平衡结构在恒载和这个初始力下再次被确定。在形状迭代的时候，一些控制点(主梁和拉索连接的点)将会被选择检验应力集中。在每次形状迭代过程中，主跨的控制点的垂直位移比率将会被检验。也就是,
7、 形状迭代将会重复直到应变可以达到所说的10-4。当应变达到的时候, 计算将会停止而斜拉桥的初始形状就找到了。数字的实验表明重复收敛于一点是没什么作用的，并且所有的三个非线性对最后的几何初始形状有比较少的影响。只有拉索下垂作用在确定初始形状分析中有显著作用，而偏压柱和大的偏转效应变则无关重要。 开始的分析能以二种不同的方式被实行: 一个线性和一个非线性计算程序。(1) 线性的计算程序:为了要找到桥的平衡结构,斜拉桥的所有非线性因素被疏忽，而只是线性的弹性拉索、梁单元、同等的线形的变形系数被使用。形状迭代是不考虑平衡迭代而实行的。合理的收敛于一点的起始形状被得到，而且许多计算的工作能被节省。(2
8、) 非线性计算程序:斜拉桥所有的非线性因素在整个的计算程序中被考虑。非线性拉索元素的下沉作用、主梁元素的稳定系数和非线性变形调整系数被应用。形状的迭代和平衡迭代都在非线性计算中实行。 牛顿-瑞普生方法在这里被用于平衡迭代。4.2 静态偏转分析基于确定的起始形状,斜拉桥在活载作用下的非线性静态偏转分析可通过模数或迭代运行。 荷载模数方法导致很大的数字错误是广为人知的。迭代方法比较适于非线性计算，而且需要的应变应能被达到。牛顿- 瑞普生的迭代程序将被使用。因为非线性分析较大或复杂的结构系统,一个完整的迭代程序(重复为一个单一全部荷载运行步骤)将会时常失败。一个模数-迭代程序高度地被推荐,荷载将会被
9、增加，而且重复将会在每个荷载步骤中实行。斜拉桥的静态偏转分析将会从使用线性或非线性计算程序决定的初始形状开始。斜拉桥静态的偏转分析的运算法则在第4.4.2节中被概述。4.3. 线性振动分析当一个结构系统是足够稳固而且外部的刺激不是太强烈,系统可能以一个确定的非线性的静态系数作一个小振幅振动，由振动引起的非线性静态系数的变化是很小的和可以忽略的。这种以一个非线性静态系数以一个小振幅的振动被称作线性化振动。线性化振动不同于线性振动，系统用很小的振幅以一个线性静态系数振动。非线性静态系数qa 能由非线性偏转分析决定。在决定qa之后，系统矩阵可能被建立有关于如此的一个非线性静态系数，线性化系统的等式如
10、下所示：MAq”+ DAq+ 2KAq=p(t)- TA上面的上标字母A 代表在非线性静态系数qa 被计算的数量。这个等式用恒定系数矩阵MA、DA、2KA表现第二的次序一组线性的一般差别的等式。这个等式能被模型的重叠方法，整体的变形方法或直接的整合方法解答。当减幅效应和荷载限制被忽略的时候,系统等式变成：MAq” + 2KAq=0 这个等式表现基于非线性静态系数qa 的不减幅的系统天然振动。天然振动的频率和模态可以从上面的等式运用程序，举例来说子空间重复方法来得到。对于斜拉桥,它的起始形状是非线性静态系数qa。 斜拉桥由于以开始的形状为基础小振幅振动的时候,天然的频率和模态能被找到来解决上述的
11、等式。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 求解增量系统的等
12、式以得到q。o 平衡迭代使用牛顿-瑞普生方法运行q。4. 变形迭代。5. 包括几何形状和基本力的初始的形状输出。6. 对于线性偏转分析，只有线性刚度单元和变形系数被采用且没有平衡迭代的实行。4.4.3. 振动分析1. 桥的几何的和物理的数据输入。2. 包括开始的几何和开始单元常受力的初始形状数据的输入。3. 建立以初始形状的自由振动的线性化系统等式。4. 运用子空间重复方法得到振动频率和模态,例如Rutishauser方法。5. 初始拉索受力估算：在王教授和林教授的最近研究中,小型的斜拉桥的通过任意小或任意大试验初始拉索应力来实现。如果不同的初始拉索受力试验评价被采用，在那里重复单调的迭代，而
13、最后的结论有相似的结果。然而对于大型斜拉桥，确定形状的计算变得更困难以达到一致。在非线性分析中，牛顿-瑞普生类型迭代计算能收敛到一点,只有当解决的被估计的价值是在真正的价值附近时才能实现。当斜拉桥分析的形状在王的论文中建议用任意小的初始拉索应力开始时，收敛到一点的困难可能会出现。因此,估计适当的试验开始的拉索应力来得到一致的结论对于形状确定分析变得重要起来。接下来,一些估计试验初始拉索应力的方法将会被讨论。5.1. 垂直荷载的平衡5.2. 零力矩的控制5.3. 零位移的控制5.4. 拉索等价系数比的概念5.5. 不对称的考虑如果估计的初始拉索应力是用上面介绍的方法对每条拉索独立地确定的，非对称
14、的斜拉桥索塔中可能会存在不平衡的水平受力。因为中间跨(主跨)和边跨的对称的拉索布置对索塔的水平分力的合力为零，也就是没有不平衡的水平力。而对于非对称的斜拉桥，拉索在中间跨和边跨分布是不对称的，拉索在索塔上产生的分力分别独立计算，很明显的索塔中的不平衡的水平力将会引起很大的弯矩和偏移。因此，非对称的斜拉桥的这个问题可以按如下解决。中间跨(主要部份)的拉索受力可以通过上面独立介绍的方法确定，其中上标m代表主跨，上标i代表第i条拉索。然后边跨上拉索受力通过与索塔连接的拉索的水平平衡方程确定，即 Tim cosi= Tis cosi, and Tis = Tim cosi/ cosi, i是指斜拉索与
15、主跨梁的夹角，i是指斜拉索与边跨梁的夹角。6. 例子在这项研究中，二座不同的类型小型斜拉桥从文学中取得,而且他们的起始形状将会用先前描述的形状确定方法使用线性和非线性程序来确定。最后，一座高度冗余的斜拉桥将会被研究。对于平衡迭代和形状迭代都采用应变10-4。重复的最大周期定为20。如果重复的循环数超过20则计算被认为是不收敛于一点的。接着的二个小型斜拉桥的开始形状在第6.1节和6.2首先采用任意小的试验拉索应力决定的。在这二个例子中收敛于一点是重复单调的。他们的收敛于一点的起始形状可以很容易地获得。由线性和非线性计算决定的开始形状之间只有很小的不同。收敛于一点的结论显示同样的结果，而且他们与试
16、验的拉索应力无关。7. 结论通过线性和非线性计算二重循环的建立而得到斜拉桥的初始形状。这个方法能达到建筑的设计形式有统一的预应力分配,而且使所有的平衡和边界情况满足。初始的形状确定是斜拉桥分析中最重要的工作。只有一个正确的起始形状，才能得到一个有意义的和正确的偏转及震动分析。基于研究的数字实验,一些结论概述如下:(1).对于小型斜拉桥初始形状的确定不会出现现在的困难，任意的初始试验拉索应力都能用来计算。然而对于大跨度的斜拉桥，循环的收敛于一点会产生很大问题。(2).在跨度的斜拉桥的形状确定的收敛于一点通常会产生困难，当拉索的试验应力通过垂直荷载平衡、零弯矩控制、零位移控制的方法给出时。(3).
17、 如果主跨的每条拉索的预应力符合到约Eeq的80%，通过两重循环的方法可以很快的找到收敛于一点的初始形状，而且边跨拉索的初始应力是由作用于索塔上的水平等式决定的。(4). 由线性程序和非线性的程序得到的初始形状的结果对于几何和预应力的分配的结果只有很小的不同。(5).使用线性的计算能提供一个合理的起始形状而且节省很多的计算工作, 所以在工程实践中它高度的被推荐。(6).在小型的斜拉桥中,由线性和非线性计算程序确定的初始形状为基础得到的自振频率只有很小的差别，而模态情形在两种情况中是一样的。(7). 基本频率和高度冗余刚性斜拉桥的显著不同在研究中被展示。只有基于非线性程序的初始形状确定的模态展示
18、非线性拉索下沉和主梁的偏转效应。举例来说，桥的第一和第三的模态被索塔的横向运动支配,而不是主梁。两个情形的基本频率的差别约为12%。只有通过非线性的计算而不是线性的计算确定的正确的初始形状才能得到斜拉桥的振动频率和模态的正确分析。Study on nonlinear analysis of a highly redundant cable-stayed bridge1AbstractA comparison on nonlinear analysis of a highly redundant cable-stayed bridge is performed in the study. The
19、 initial shapes including geometry and prestress distribution of the bridge are determined by using a two-loop iteration method, i.e., an equilibrium iteration loop and a shape iteration loop. For the initial shape analysis a linear and a nonlinear computation procedure are set up. In the former all
20、 nonlinearities of cable-stayed bridges are disregarded, and the shape iteration is carried out without considering equilibrium. In the latter all nonlinearities of the bridges are taken into consideration and both the equilibrium and the shape iteration are carried out. Based on the convergent init
21、ial shapes determined by the different procedures, the natural frequencies and vibration modes are then examined in details. Numerical results show that a convergent initial shape can be found rapidly by the two-loop iteration method, a reasonable initial shape can be determined by using the linear
22、computation procedure, and a lot of computation efforts can thus be saved. There are only small differences in geometry and prestress distribution between the results determined by linear and nonlinear computation procedures. However, for the analysis of natural frequency and vibration modes, signif
23、icant differences in the fundamental frequencies and vibration modes will occur, and the nonlinearities of the cable-stayed bridge response appear only in the modes determined on basis of the initial shape found by the nonlinear computation.2. IntroductionRapid progress in the analysis and construct
24、ion of cable-stayed bridges has been made over the last three decades. The progress is mainly due to developments in the fields of computer technology, high strength steel cables, orthotropic steel decks and construction technology. Since the first modern cable-stayed bridge was built in Sweden in 1
25、955, their popularity has rapidly been increasing all over the world. Because of its aesthetic appeal, economic grounds and ease of erection, the cable-stayed bridge is considered as the most suitable construction type for spans ranging from 200 to about 1000 m. The worlds longest cable-stayed bridg
26、e today is the Tatara bridge across the Seto Inland Sea, linking the main islands Honshu and Shikoku in Japan. The Tatara cable-stayed bridge was opened in 1 May, 1999 and has a center span of 890m and a total length of 1480m. A cable-stayed bridge consists of three principal components, namely gird
27、ers, towers and inclined cable stays. The girder is supported elastically at points along its length by inclined cable stays so that the girder can span a much longer distance without intermediate piers. The dead load and traffic load on the girders are transmitted to the towers by inclined cables.
28、High tensile forces exist in cable-stays which induce high compression forces in towers and part of girders. The sources of nonlinearity in cable-stayed bridges mainly include the cable sag, beam-column and large deflection effects. Since high pretension force exists in inclined cables before live l
29、oads are applied, the initial geometry and the prestress of cable-stayed bridges depend on each other. They cannot be specified independently as for conventional steel or reinforced concrete bridges. Therefore the initial shape has to be determined correctly prior to analyzing the bridge. Only based
30、 on the correct initial shape a correct deflection and vibration analysis can be achieved. The purpose of this paper is to present a comparison on the nonlinear analysis of a highly redundant stiff cable-stayed bridge, in which the initial shape of the bridge will be determined iteratively by using
31、both linear and nonlinear computation procedures. Based on the initial shapes evaluated, the vibration frequencies and modes of the bridge are examined.3. System equations3.1. General system equationWhen only nonlinearities in stiffness are taken into account, and the system mass and damping matrice
32、s are considered as constant, the general system equation of a finite element model of structures in nonlinear dynamics can be derived from the Lagranges virtual work principle and written as follows:Kjbj-Sjaj= Mq”+ Dq3.2. Linearized system equationIn order to incrementally solve the large deflectio
33、n problem, the linearized system equations has to be derived. By taking the first order terms of the Taylors expansion of the general system equation, the linearized equation for a small time (or load) interval is obtained as follows: Mq”+Dq +2Kq=p- up3.3. Linearized system equation in staticsIn non
34、linear statics, the linearized system equation becomes2Kq=p- up4. Nonlinear analysis4.1. Initial shape analysisThe initial shape of a cable-stayed bridge provides the geometric configuration as well as the prestress distribution of the bridge under action of dead loads of girders and towers and unde
35、r pretension force in inclined cable stays. The relations for the equilibrium conditions, the specified boundary conditions, and the requirements of architectural design should be satisfied. For shape finding computations, only the dead load of girders and towers is taken into account, and the dead
36、load of cables is neglected, but cable sag nonlinearity is included. The computation for shape finding is performed by using the two-loop iteration method, i.e., equilibrium iteration and shape iteration loop. This can start with an arbitrary small tension force in inclined cables. Based on a refere
37、nce configuration (the architectural designed form), having no deflection and zero prestress in girders and towers, the equilibrium position of the cable-stayed bridges under dead load is first determined iteratively (equilibrium iteration). Although this first determined configuration satisfies the
38、 equilibrium conditions and the boundary conditions, the requirements of architectural design are, in general, not fulfilled. Since the bridge span is large and no pretension forces exist in inclined cables, quite large deflections and very large bending moments may appear in the girders and towers.
39、 Another iteration then has to be carried out in order to reduce the deflection and to smooth the bending moments in the girder and finally to find the correct initial shape. Such an iteration procedure is named here the shape iteration. For shape iteration, the element axial forces determined in th
40、e previous step will be taken as initial element forces for the next iteration, and a new equilibrium configuration under the action of dead load and such initial forces will be determined again. During shape iteration, several control points (nodes intersected by the girder and the cable) will be c
41、hosen for checking the convergence tolerance. In each shape iteration the ratio of the vertical displacement at control points to the main span length will be checked, i.e., The shape iteration will be repeated until the convergence tolerance, say 10-4, is achieved. When the convergence tolerance is
42、 reached, the computation will stop and the initial shape of the cable-stayed bridges is found. Numerical experiments show that the iteration converges monotonously and that all three nonlinearities have less influence on the final geometry of the initial shape. Only the cable sag effect is signific
43、ant for cable forces determined in the initial shape analysis, and the beam-column and large deflection effects become insignificant.The initial analysis can be performed in two different ways: a linear and a nonlinear computation procedure.1. Linear computation procedure: To find the equilibrium co
44、nfiguration of the bridge, all nonlinearities of cable stayed bridges are neglected and only the linear elastic cable, beam-column elements and linear constant coordinate transformation coefficients are used. The shape iteration is carried out without considering the equilibrium iteration. A reasona
45、ble convergent initial shape is found, and a lot of computation efforts can be saved.2. Nonlinear computation procedure: All nonlinearities of cable-stayed bridges are taken into consideration during the whole computation process. The nonlinear cable element with sag effect and the beam-column eleme
46、nt including stability coefficients and nonlinear coordinate transformation coefficients are used. Both the shape iteration and the equilibrium iteration are carried out in the nonlinear computation. NewtonRaphson method is utilized here for equilibrium iteration.4.2. Static deflection analysisBased
47、 on the determined initial shape, the nonlinear static deflection analysis of cable-stayed bridges under live load can be performed incrementwise or iterationwise. It is well known that the load increment method leads to large numerical errors. The iteration method would be preferred for the nonline
48、ar computation and a desired convergence tolerance can be achieved. Newton Raphson iteration procedure is employed. For nonlinear analysis of large or complex structural systems, a fulliteration procedure (iteration performed for a single full load step) will often fail. An incrementiteration procedure is highly recomme