基于MATLAB语言的悬索桥施工过程颤振分析.pdfVIP

基于MATLAB 语言的悬索桥施工过程颤振分析 徐优 西南交通大学土木工程学院，成都（610031 ） E-mail ：stickxu@163.com 摘 要：MATLAB 语言由于其强大的矩阵运算功能，在桥梁结构有限元分析中具有广阔的 应用前景。本文首先基于大跨度悬索桥三维颤振的有限元分析方法，建立气动自激力作用下 的结构颤振运动方程，将颤振问题转化为求解复数矩阵的标准特征值问题，随后利用 MATLAB 语言编制相应程序。最后用该程序对某大跨度悬索桥的各个施工过程进行颤振分 析，据此评价该桥的颤振稳定性。 关键词：悬索桥，施工过程，颤振，MATLAB 1 引言 桥梁颤振是一种自激发散性振动，当桥梁结构由结构阻尼和因自激力引起的气动阻尼组 成的系统阻尼由正值趋向于负值时，振动系统无法消耗由于气流反馈作用吸收的能量，使得 【】 振幅逐步增大导致颤振 1 。当振幅大到一定程度时, 将导致结构的破坏。在大跨度悬索桥 的抗风设计中, 颤振是影响设计非常重要的因素之一。因此, 对悬索桥颤振问题进行研究具 有重大工程实际意义。 迄今为止, 国内外许多学者提出了一些行之有效的三维颤振分析方法。如Namini提出的 【】 【】 【】 2 3 4 PK-F法 ，谢霁明、项海帆提出了状态空间法 ，陈政清提出的M －S法 。由于计算机 技术的广泛运用使得三维颤振分析得到快速发展。本文将采用含有18个颤振导数的气动自激 力模型，并将其表示成复数的形式，然后将结构颤振运动方程转化成求解复数矩阵的标准特 征值问题，通过求解特征值进行颤振分析。 【～ 】 2 多模态参与颤振分析理论 5 6 采用Scanlan 提出的自激气动力模型，并将其表示成复数的形式，则作用在桥梁主梁单 L M 可以表示为如下的形式： 位长度上的升力 、阻力 、升力矩 D h p α 2 2 ⎡ * * * * * * ⎤ ωρ =+ + + α + + Lh B ⎢⎣(H 4 iH1 )h (H 3 iH 2 ) B (H 6 iH 5 )p ⎥⎦ 2 2 ⎡ * * * * * * ⎤ =+ + + + + D ω ρB ⎢(P iP )h (P iP ) p (P iP )Bα ⎥ p ⎣ 6 5 4 1 3 2 ⎦ （1） 2 2 ⎡ * * * * * * 2 ⎤ M B (A iA )Bh (A iA )Bp (A iA )B ω ρ ⎢ =+ +