斜拉桥动力特性分析

第２１卷第５期　四川理工学院学报（自然科学版）　Ｖｏ１．２１　Ｎｏ．５　Ｊ０ＵＲＮＡＬ　ＯＦ　ＳＩＣＨＵＡＮ　ＵＮＩＶＥＲＳＩＴＹ　ＯＦ　２００８年１０月　ＳＣＩＥＮＣＥ＆ＥＮＧＩＮＥＥＲＩＮＧ（ＮＡＴＵＲＡＬ　ＳＣＩＥＮＣＥ　ＥＤＩＴＩＯＮ）　０ｃｔ．２０ｏ８　文章编号：１６７３—１５４９（２００８）０５一Ｏ１０５—０３　斜拉桥动力特性分析　李文静　，张坤　（１．兰州交通大学土木工程学院，兰州７３００７０；２．甘肃科地工程咨询有限公司，兰州７３００２０）　摘要：以某斜拉桥为工程背景，基于ＡＮＳＹＳ建立了三维有限元分析模型，分析了斜拉桥自身的动力特性，　主要包括自由振动频率、振型特性。结果表明：大跨度斜拉桥自振周期长，频谱密集，自振特性表现出明显的三维性　和振型相互耦合的特点。　关键词：斜拉桥；大跨度；动力特性；有限元分析模型　中图分类号：Ｕ４４８　文献标识码：Ａ　引言　橡胶板抗震限位装置，边墩竖向支座采用耐久性好的抗　斜拉桥跨越能力强、造价经济、造型优美。近５Ｏ年来　震型球型钢支座，承载力８０００ｋＮ，横向设置２个。　得到迅速发展，并以其超强的跨越能力和出色的优点日　益向大跨度方向迈进。随着斜拉桥跨度的不断增大，：０２　００　４：￣ｏ　ｏ＿其　１９２ｏｏ　　，．结构刚度越来越柔，它在动力荷载（如风、地震和汽车荷　—载等）作用下的动力特性和结构性能倍受工程界关注，　很多问题也随之而来，其中桥跨结构的自振特性、抗震、　叠　－　－　抗风以及车辆荷载的冲击振动等动力学问题尤为突出。　图１斜拉桥立面布置图　斜拉桥的动力特性是反映其抗风、抗震性能的主要指标　２斜拉桥的动力特性分析　之一。　　．结构的动力响应取决于外部条件的干扰和结构自　１工程概述　身的动力特性，在分析结构在外荷载作用下，如地震响　某大桥为三跨全漂浮体系ｆ１９２ｍ＋４３４ｍ＋１９２ｍ），双　应、风振稳定等，必须先进行自振特性的分析。＼　●－　●　　塔双索面叠合梁式斜拉桥。全桥立面如图１所示。大桥　动力特性主要指固有频率、振型、阻尼等，是结构本　设计安全等级为一级，桥梁设计基准期为１００年，设计　身固有的属性。斜拉桥的动力特性受主梁截面性质、塔　车速为８０ｋｍ／ｈ，桥面全宽为３０ｍ，最大纵坡为１％，横向　与主梁刚度比、塔梁连接方式、结构的空间几何布置形　设置２％的横坡度。上部结构为两个钢箱梁和混凝土桥　式、辅助墩的数目与位置、塔梁连接方式等多种因素的　面板组成的叠合梁。桥面板厚０．２５ｍ，主梁为两分离高　影响。　２．８ｍ，宽２．３ｍ钢箱梁，钢箱梁之间，每４ｍ设置一道横隔　本文建立了该斜拉桥的动力计算模型如图２所示。　梁，横隔梁间设置两道小纵梁，端横梁处设置铸铁压重，　每侧重１６０４４ｋＮ。索塔为门形塔，钻孔桩高桩承台基础，　塔高１４１．３３１ｍ。塔柱设置三道横梁，塔柱均为空心薄壁　结构。下塔柱高２４．２４１ｍ，变截面单箱单室，下大上小，　中塔柱高５４．０９ｍ，上塔柱高５９ｍ，均为等截面箱形断面，　三道横梁亦为箱型截面，七、中横梁高４ｍ，下横梁高　６ｍ。全桥共１４０根斜拉索，采用半平行钢丝索。为了提　图２斜拉桥动力计算分析模型　高主桥的抗震性能，在桥塔处，纵向主箱外侧，设置横向　２．１模态分析数值方法　收稿１３期：２００８．０５．２９　作者简介：李文静（１９８０一），女，内蒙古赤峰人，硕士生，主要从事桥梁抗震方面的研究。　ｌ０６　四川理工学院学报（自然科学版）　２００８年１０月　根据有限元分析结构自由振动的动力平衡方程为：　５＋　６　（１）　式中：　一结构的质量矩阵　一结构的刚度矩阵　其解表示为：　６（ｆ）＝０ｓｉｎ（“　＋　）　（２）　将公式（２）代入公式（１）可得特征方程：　（　—ｏｊ２Ｍ）０＝０　（３）　根据Ｃｒａｍｅｒ法则，可得频率方程　ｌ　ｌＫ一　ｌｌ＝０　（４）　通过对方程（４）特征值求解，可得到各阶频率∞和　振型　。　数学上求解特征值问题的经典方法是雅可比法，其　优点是程序简单，数值稳定性好，但由于雅可比法一次　要求出全部的特征值和特征向量，这在矩阵的阶数很高　时是很不利的，将导致求解费用的增加。大跨度斜拉桥　的节点数和结构的总的自由度数很多，在工程实际中，　往往只需要了解其部分较低阶数的特征值及相应特征　向量即可。因此，在有限元分析中发展了一些适应上述　特点的效率较高的特征值求解方法，其中应用较广泛的　是矩阵逆迭代法和子空间迭代法。　子空间迭代法是目前求解大型矩阵特征值问题的　最常用的有效方法之一，已经广泛应用于大型结构动力　学的有限元分析中。子空间迭代法主要利用了质量和刚　度矩阵的带状特性，为精确地求得一个体系的Ｐ个振型　和频率，首先取ｑ个试探向量（ｇ　），位移可以表示为如　式（５）的多个形状的组合。　／Ｊ‘０　＝　‘０　Ｚ‘。　＝　‘。　（５）　其中初始广义矩阵仅仅是一个单位矩阵。对于Ｐ　个特征值和特征向量的广义特征值方程可写成　＝ｍ中Ａ　（６）　式中，矩阵　中的各列代表所要计算的特征向量，　Ａ是一个对角线元素为特征值的对角矩阵，即：　Ａ　Ａ＝　（７）　中＝｛　ｔ，　ｚ，…，　尸）　（８）　在方程（６）右边代入ｑ个试探向量，导得　‘　＝栅‘。　＝Ｗ‘。　（９）　求解式（９）得到未规格化的改进的形状，即　‘　＝　一　‘∞　（１０）　式（１０）改进的形状用于下一轮新的迭代之前，它们　必须用规格化和正交化修正。规格化使其在计算中数值　大小保持合理，正交化使得每一个向量收敛于不同振型　（而不是全部收敛于最低振型）。这些运算能按多种不同　的方式来实现，但是最方便而且两者立即就可以完成的　方法是进行一次Ｒｉｔｚ特征问题的分析。第一次循环的广　义坐标刚度和质量矩阵，计算如下：　：＝　¨’　¨　；沙¨　ｙｍ　’　（１１）　硝＝　“　”　（１２）　求解对应的特征问题　ｋ：２（１）：　２（１）　（１３）　得到广义坐标振型２¨　和频率　，。　重复整个迭代过程按公式（１０）求得未规格化的改　进振形　，然后求解对应的特征问题，求得规格化和　正交化的　∽）＿　等，迭代过程最后收敛于真实的　振型和频率即：　‘　）＿÷　２—　Ａ　当ｓ　∞　（１４）　２．２斜拉桥模态分析　本文采用ＡＮＳＹＳ程序中的Ｓｕｂｓｐａｃｅ（子空间）法对　斜拉桥进行模态分析，求解其动力特性。共计算了前３００　阶自振模态。前１０阶模态频率及振型特征见表１，前４　阶振型如图３至图６所示，其中第一阶振型为纵飘，侧　弯首次出现在第二阶振型，竖弯首次出现在第四阶振　型。图中振型图左上角为正等轴测视图，右上角为横桥　向视图，左下角为纵桥向视图，右下角为竖向视图。　表１前１０阶振型自振频率与特征　。　一　．　图３第一振型ｆ＝０．１０４３５Ｈｚ　第２１卷第５期　李文静等：斜拉桥动力特性分析　１０７　较小。此时，多维地震响应可以近似地简化为多个一维　地震响应的线性叠加。　３结束语　＼　～ｆ＝Ｏ．１９４８１Ｈｚ　＼　本文首先建立了斜拉桥三维有限元计算分析模型。　然后采用ＡＮＳＹＳ模态分析的子空间迭代法求解了斜拉　桥的动力特性。分析结果表明：　卿　嘲　（１）大跨度斜拉桥的柔度比较大，自振周期长，固有　频率低，具有密布的频谱，模态远较一般的结构密集。　（２）由于该斜拉桥采用全漂浮体系，而且跨度很大，　动力特性分析显示第一阶振型为纵飘，白振频率为　０．１０４３５Ｈｚ，符合飘浮体系斜拉桥的特征。　（３）主梁的对称竖向弯曲振型为纯竖向弯曲，主梁　跨中没有纵向位移，但主梁反对称竖向弯曲振动常与主　梁的纵向振动耦合在一起。　（４）双索面的采用和箱形截面的结合使用等均能提　高结构的抗扭刚度，从而提高了扭转频率。　（５）该桥的低阶振型在主方向间的耦合作用不明　显，前１０阶振型均表现为单一的振型形态。所以，当某　图５第三振型ｆ＝Ｏ．２６２６７Ｈｚ　一方向的地震波输入时只会引起结构沿本方向较大的　内力反应，而其他方向的反应则较小。此时，多维地震响　一　－＿＿　硎　＿　应可以近似地简化为多个一维地震响应的线性叠加。　参考文献：　【１］刘士林斜拉桥［ＭＩ．北京：人民交通出版社．２００２．　【２］林元培．斜拉桥【Ｍ］．北京：人民交通出版社．２００４．　【３】严国敏．现代斜拉桥［ＭＩ．成都：西南交通大学出版社．１９９６．　ｆ…一　一　一Ｉｆ＿一。。　＿　【４】范立础．桥梁抗震【Ｍ】．上海：同济大学出版社．１９９７．　【５】范立础，胡世德，叶爱君．大跨度桥梁抗震设计【Ｍ】．北京：人　民交通出版社．２００１．　图６第四振型ｆ＝Ｏ．２７２３５Ｈｚ　【６】范立础，王志强．桥梁减隔震设计［Ｍ］．北京．人民交通出版　从表１计算结果可以看出：该桥的低阶振型在主方　社．２００１．　向问的耦合作用不明显，前１０阶振型均表现为单一的　【７】周福霖．工程结构减震控制【Ｍ】．北京：地震出版社．１９９７．　振型形态。这就意味着，某一方向的地震波输入只会引　【８】尹晓明，朱益民．斜拉桥动力模型分析【Ｊ】．铁道标准设计，　起结构沿本方向较大的内力反应，而其他方向的反应则　２００３．（１０）：３８—３９．　Ｄｙｎａｍｉｃ　Ｃｈａｒａｃｔｅｒｉｓｔｉｃ　Ａｎａｌｙｓｉｓ　ｆｏｒ　Ｃａｂｌｅ－Ｓｔａｙｅｄ　Ｂｒｉｄｇｅ　ＬＩ　Ｗｅｎ－ｉｒｎｇ’，ＺＨＡＮＧ　Ｋｕｎ　（１．Ｓｃｈｏｏｌ　ｏｆ　Ｃｉｖｉｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，Ｌａｎｚｈｏｕ　Ｊｉａｏｔｏｎｇ　Ｕｎｉｖｅｒｓｉｔｙ，Ｌａｎｚｈｏｕ　７３００７０，Ｃｈｉｎａ；　２．Ｇａｎｓｕ　Ｋｅｄｉ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｃｏｎｓｕｌｔｉｎｇ　Ｃｏ．Ｌｔｄ．，Ｌａｎｚｈｏｕ　７３００２０，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ａ　ｃａｂｌｅ—ｓｔａｙｅｄ　ｂｒｉｄｇｅ　ｉｓ　ｔａｋｅｎ　ｆｏｒ　ｓｔｕｄｙ　ａｎｄ　ａ　３－Ｄ　ｆｉｎｉｔｅ　ｅｌｅｍｅｎｔ　ａｎａｌｙｓｉｓ　ｍｏｄｅ　ｗｉｔｈ　ＡＮＳＹＳ　ｉｓ　ｅｓｔａｂｌｉｓｈｅｄ，ｔｈｅ　ｓｐｅｃｉａｌ　ｄｙｎａｍｉｃａｌ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｔｈｉｓ　ｂｒｉｄｇｅ　ａｒｅ　ｉｎｖｅｓｔｉｇａｔｅｄ，ｉｎｃｌｕｄｉｎｇ　ｔｈｅ　ｎａｔｕｒａｌ　ｆｒｅｑｕｅｎｃｙ，ｔｈｅ　ｆｅａｔｕｒｅｓ　ｏｆ　ｎａｔｕｒａｌ　ｍｏｄｅｌｓ．Ｔｈｅ　ｒｅｓｕｌｔ　ｓｈｏｗｓ　ｔｈａｔ　ｔｈｅ　ｌａｒｇｅ　ｓｐａｎ　ｃａｂｌｅ－ｓｔａｙｅｄ　ｂｒｉｄｇｅ　ｉｓ　ｌｏｎｇ　ｎａｔｕｒａｌ　ｃｙｃｌｅ，ｓｐｅｃｔｒｕｍ－　ｉｎｔｅｎｓｉｖｅ，ｔｈｅ　ｖｉｂｒａｔｉｏｎ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｃｌｅａｒｌｙ　ｓｈｏｗ　ｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａｌ　ｍｏｄｅｓ　ａｎｄ　ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｍｕｔｕａｌ　ｃｏｕｐｌｉｎｇ．　Ｋｅｙ　ｗｏｒｄｓ：ｃａｂｌｅ－－ｓｔａｙｅｄ　ｂｒｉｄｇｅ；ｌｏｎｇ・・ｓｐａｎ；ｄｙｎａｍｉｃａｌ　ｃｈａｒａｃｔｅｒｉｓｔｉｃ；ｆｉｎｉｔｅ　ｅｌｅｍｅｎｔ　ａｎａｌｙｓｉｓ　ｍｏｄｅｌ