悬索桥静风扭转发散有限元精细化分析.docVIP

悬索桥静风扭转发散有限元精细化分析 　　摘要：基于大跨度悬索桥刚度退化及静风扭转发散的机理，选取了两个评价标准分别用以判断均匀流场和紊流场中悬索桥的刚度退化情况.在均匀流场中，由于主缆的变形形态与演变规律都很简单，因此可以选取主缆中点的竖向位移作为评价刚度退化及扭转发散的标准.当主缆中点向上的竖向位移恰好使其应力松弛时，结构将会出现扭转发散现象，这一竖向位移称为临界竖向位移.然而，在紊流场中，紊流中的脉动成分往往会引起结构显著的多模态耦合叠加的复杂响应.此时，前述的标准不再适用.为此，本文提出了一种基于识别时域范围内主缆长度的方法，即在时域范围内，当任意一条主缆的长度的最小值达到或十分接近无应力长度时，主缆将软化进而引起间歇式扭转发散.静、动力有限元分析表明，上述方法可以较好地解释大跨度悬索桥在不同流场中的扭转发散现象. 　　关键词：悬索桥；扭转发散；刚度退化；评价标准；有限元分析 　　中图分类号：U448.25 文献标识码：A 　　对于大跨度悬索桥而言，静风扭转失稳可能先于动力失稳出现.因此，大跨度悬索桥必须具备良好的静风稳定性能.扭转发散概念来源于航空工程学中机翼的静风扭转发散现象，即机翼的扭转变形随着风速的逐渐增加而缓慢增加，当风速达到某一值时，扭转变形急剧跃增至一个很大的程度而导致结构发生破坏[1-2].Simiu与Scanlan在大跨度桥梁抗静风研究中首次引用了这一概念，并提出了一种简单估算主梁扭转发散临界风速的线性方法[3]，但是它是基于线性假定而提出的，忽略了所有的非线性效应，这显然不能准确地评估桥梁的静风稳定性能. 　　为了考虑非线性因素的影响，Boonyapinyo等提出了一种静力有限位移法来分析大跨度斜拉桥的耦合屈曲问题[4-5].不久，Chen等采用此方法对大跨度悬索桥的静风失稳问题进行了研究[6-7].在这些研究中，均采用了静力有限元方法，并考虑了结构的几何非线性、材料非线性及风荷载非线性的影响.需要提出的是，静力有限元方法可以较好地揭示桥面扭转发散的机理与发展的全过程，但是无法考虑紊流引起的随机响应对扭转发散的影响.众所周知，大气边界层附近的流场都是紊流.因此，有必要采用动力有限元方法来探究紊流对大跨度悬索桥静风扭转发散的影响. 　　已有研究表明，随着作用在结构上向上荷载的逐渐增减，主缆系统将会出现软化现象，这表明向上的运动对结构的刚度会产生不利的影响[8].为了讨论紊流的影响，张志田等采用了动力有限元方法研究了一座大跨度悬索桥的静风稳定性[9].结果表明，紊流引起的随机动力响应明显地降低了桥梁的静风稳定性能.但是大跨度悬索桥的静风扭转发散机理那时依然不清楚.后来.张志田等基于风洞试验成果并建立了广义梁索分析模型，较好地揭示了大跨度悬索桥扭转发散的本质原因是主缆系统扭转刚度的退化[10-11].对于大跨度悬索桥而言，结构体系的扭转刚度主要来源于由两根主缆组成的主缆系统而不是桥面[12].因此，当主缆系统的扭转刚度退化后，桥面极易在风荷载的作用下发生扭转破坏. 　　本文基于大跨度悬索桥刚度退化及静风扭转发散的机理，采用静、动力有限元分析方法分别研究了均匀流场及紊流场中矮寨大桥的静风稳定性能，并提出了这两种流场各自适用的刚度退化的评价标准.在均匀流场中，由于主缆的变形形态与发展规律明确且易于描述，因此可以选取主缆中点的临界竖向位移作为评价标准.当任意一条主缆向上的竖向位移达到临界值时，主缆将发生软化并致使主缆系统丧失了绝大部分的扭转刚度，这一过程即为刚度退化过程.刚度退化后，桥梁结构只残余主梁提供的扭转刚度.然而这一部分刚度很容易被风速引起的气动扭转负刚度抵消，最终导致桥面发生扭转破坏.本文采用静力有限元的方法求解了临界竖向位移值，并评估其他一些响应对该临界值的影响.然而，在紊流场中，脉动风将会引起结构显著的随机动力响应，主缆的响应也变得十分复杂.此时，前述的标准将不再适用.为此，本文提出了一种基于识别时域范围内主缆长度的方法，即在时域范围内，当任意一条主缆的长度的最小值达到或十分接近无应力长度时，主缆将发生软化而使得桥面出现间歇式扭转发散现象.本文将采用动力有限元分析方法并结合这一标准来解释紊流场中大跨度悬索桥扭转发散特性. 　　3有限元分析 　　3.1工程概况 　　以矮寨大桥为有限元算例.矮寨大桥为一座山区峡谷大跨度悬索桥，采用钢桁架梁方案.桥梁主跨1 176 m，加劲梁全长1 000.5 m；主缆矢跨比1∶9.6，两条主缆形心距27 m；梁高7.5 m，桁架节间距7.25 m.全桥布局如图1所示.表1列出了矮寨大桥的主要计算参数.图2为全桥ANSYS有限元模型，采用的是单主梁模型，共有471个节点和687个单元.图3给出了加劲梁断面的静力三分力系数.计算时主梁上同时施加阻力、升