第7卷第6期 2010年12月 现代交通技术 Modern Transportation Technology VO1.6 NO.7 Dec.2010 大跨度桥梁涡激振动及其半经验模型介绍 朱亚飞 ,鲜 荣2,刘运志, (1.江苏省交通科学研究院长大桥梁健康检测与诊断技术交通行业重点实验室,江苏南京211112; 2.广东公路建设有限公司,广东广州510600;3.中国三峡新能源公司,北京100053) 摘要:大跨度桥梁在低风速下可能会频繁发生较大振幅的涡激共振而导致结构破坏,因而涡激振动已成为大 跨度桥梁抗风设计关键问题之一。该文列举了国外大跨度桥梁典型涡激振动实例,并讨论了具有代表性的一些 半经验模型,最后针对大跨度桥梁涡激振动提出了基于风洞试验和半经验模型的研究思路。 关键词:大跨度桥梁;涡激振动;半经验模型;风洞试验 中图分类号:U441.3 文献标识码:A 文章编号:1672—9889(2010)05~0039—06 Vortex-induced Vibration of Long Span Bridges and Introduction of Semi-emperical Model Zhu YafeiI,Xian Rong2,Hu Yunzhi (1.Key Laboratory of Long-span Bridge Health Detection and Diagnostic Technology of Communications,Jiangsu Transportation Research Institute,Nanjing 21 I 1 12,China;2.Guangdong Road Construction Limited Company, Guangzhou 5 1 0600,China;3.China Three Gorges New Enery gCorp.,Beijing 1 00053,China) Abstract:The bridge could be forced into vibration with large amplitude in low wind speed because of vortex shedding, and is easy to be damaged.Therefore,vortex—induced vibration(VIV)has been one of the most important questions in wind resistance design of long span bridges.This paper introduces some examples of VIV of long span bridges in abroad, and then discusses a few of representative semi-experienced models of VIV.Lastly,one way based Off wind tunnel test and semi・experienced model is given here to research VIV。of long span bridges in the future. Key words:long span bridge;vortex—induced vibration;semi—emperical model;wind tunnel test 2l世纪桥梁将向更大跨度发展,结构更轻柔、 阻尼更小、自振频率更低且分布密集,低风速下即 将频繁发生较大振幅的涡激共振而影响结构正常 运营甚至导致结构破坏,因而涡激振动已成为大跨 度桥梁抗风设计关键问题之一。涡激振动带有自激 性质,但振动的结构反过来会对涡脱形成某种反馈 振动研究提出了基于半经验模型和风洞试验的研 究思路,为大跨度、超大跨度桥梁涡激振动研究给 出了一种研究模式。 一定条件下.一个固定的钝体会脱落出交替的 旋涡,其主频率Ⅳ可以由斯托劳哈(Strouha1)关系式 得出: :作用。使得涡振振幅受到。因此涡激共振是一 种带有自激性质的风致限幅振动。尽管涡激振动不 象颤振、弛振一样是发散的毁灭性的振动,但由于 s 式中:.s 为斯托劳哈数;D为物体的横风向尺寸; 是低风速下容易发生的振动。且振幅之大足以影响 行车安全,因而在施工或者成桥阶段避免涡激共振 为浸没物体均匀流动的平均速度。当弹性支承物体 在尾流中脱落旋涡作用下,将被周期性地驱动,但 是这种驱动力只引起很小的响应,除非交替变化压 力的斯托劳哈频率接近物体的横流向固有频率。当 或其振幅在可接受的范围之内具有十分重要 的意义[ 。本文首先介绍国外典型大跨度桥梁涡激 振动实例。然后介绍了涡激振动研究方法并重点论 述了半经验模型的应用,最后针对大跨度桥梁涡激 接近这一频率时,物体的机械频率控制了旋涡脱落 现象。甚至当风速的变化使名义斯托劳哈频率已经 作者简介:朱亚飞(1979一),男,江苏东海人,工程师,主要从事桥梁荷载试验、施工监控及桥梁检测工作 第6期 朱亚飞,等:大跨度桥梁涡激振动及其半经验模型介绍 。41・ 头录制到全过程。这些现象导致该桥结构动力性 能、气动弹性,以及桥梁的完善性、适用性和总体安 全等均引起公众极大关注。 图5巴西Rio—Niter6i桥 1.3 丹麦Great Belt East桥涡激振动 丹麦Great Belt East桥也出现了明显的涡激振 动现象,图6所示为该桥涡激共振时拍得的照片。 图中车辆是静止在桥面上的,左右两张照片是在同 一位置不同时间所拍摄。从车辆位置变化可知其振 幅之大^10]。此外,Great Belt East桥引桥在修建之处 亦发生了明显的涡激振动,Larsen对其作了相关研 究[I1] 图6丹麦GreatBeltEast桥涡振 此外,英国Kossock斜拉桥出现了涡激振动现 象;对意大利Messina桥的分析同样表明主梁涡激 振动必须引起重视。综上所述,大跨度桥梁主梁涡 激振动发生频率高、危害大,在抗风设计时必须加 以考虑。 2涡激振动半经验模型 涡激振动研究可以大致归纳为解析方法、数值 方法、基于风洞试验的半经验模型方法。 至今没有提出一种完全成功的解析方法能从 基本流动原理出发,表现在旋涡脱落作用下一个弹 性钝体响应动作的全过程。那么理论上要解决复杂 外形钝体涡振更是十分困难,而简化模型也会对结 果产生一定影响,所以对于复杂钝体外形得到涡振 解析解难以实现…。虽然近年来随着计算机技术的 进步数值方法已经取得了很大进展,但在处理雷诺 数达10 量级的3维复杂钝体桥梁涡激振动问题上 仍然不能实现。 已经证实非常有效的办法是根据风洞试验结果 建立半经验模型,通过合理地选择参数,使模型的特 性与真实情况相吻合。涡振半经验模型大致可分为3 种形式:尾流一振子模型(Wake.Oscillator Models)、单 自由度模型(Single degree—of-freedom models)和力分 解模型(Force.decomposition models)。本文主要介绍 尾流一振子模型和力分解模型。 2.1尾流一振子模型(Wake Oscillator Models) Bishop和Hassan率先采用Van der Pol型振子 描述柱体涡脱气动力时程[ ],Hartlen和Currie建立 了最具代表性的涡振模型( ],利用升力van der Pol 型弱非线性振子方程与柱体运动速度线性耦合。升 力垂直于来流方向和柱体轴向。柱体受线性弹簧和 线性阻尼约束,模型采用1对耦合无量纲微分方程。 cLn 0+2 r02 t +2 t r+上=口 ccL ) +to ̄c£=bx ll f,)、 n,0 J 式中:对无量纲时间求导I"=00 ̄t,簖为柱体无量纲位 移,C 为升力系数,O)0为Strouhal频率与柱体自然 频率之比000=3 ̄ , 为材料阻尼系数,参数口为无 量纲常数。对于待定参数0、 和b,只需要选择其中 2个参数与试验数据拟合,因为a、 之间相关关系 为c, ̄(4oa3y)I/2 c伪为固定柱体c 振幅。值得提及 的是(2)式等号左边第2项表示c 的增加,第3项 则其增长范围,这些项具有十分重要作用,因 为涡激振动大振幅特性为伴随一个明显而有限的 升力系数增长。在适当选择参数后,Hartlen和 Currie模型可以从量化描述许多从试验结果中得到 的特性,例如当涡脱频率与柱体自然频率接近时。 会出现较大共振响应振幅。振动频率在共振区域近 似保持常量且与柱体自然频率接近。 Skop和Griifn改进发展了Hartlen和Currie模 型[H],改进了之前模型参数与系统物理参数无关的 现象。Griffin等比较了半经验模型结果与弹性支撑 柱体风洞试验结果[15].在雷诺数范围350~1000之 间对振幅和频率进行测量,所有的测量均在涡脱频 率和柱体振动频率间同步现象中进行。 1wan和Blevins得出了流场振子方程[16].该模 型基于一个隐式流场变量z表示了问题的流体动 力效应。Landl在他的2方程钝体涡激振动模型中 添加了5阶非线性启动阻尼项[I7],作者认为在流 场振子方程中加入阻尼项 c 可更好地描述 滞后效应。 ・42・ 现代交通技术 2010盎 Facchinetti等对2维涡振流场振子动力学模型 作了很好的回顾[ j。检验了结构作用力与流场振子 问3种耦合项。研究发现位移和速度耦合不能描述强 迫振动柱体涡脱试验中出现的升力相位,位移耦合不 能描述锁定区域时的升力增大现象,亦不能描述低 Skop—Griifn参数s 下涡激振动的所有显著特I生,而 加速度耦合可以成功定量模拟涡激振动特I生。 1wan将1wan—Blevins模型扩展应用于描述不 同边界条件张拉索和柱体梁的最大涡振振幅E .推 导了结构第n阶振型下的最大响应振幅为 ~ 一 、/4/3 『i【 !ae二 2+ ’ cv2 s J1 , / 有趣的是上式中与结构相关的仅有无量纲振 型参数 ,质量比IX =4m/p'rrD 和结构模态阻尼比 , 为 之和,ai为无量纲模态参数。 Skop和Griifn从他们自己的流场一振子模型人 手,以几乎相同的形式将其拓展应用与弹性柱体 。 他们得到第i阶纯模态下的最大振幅如下 Y/,maDx( )一A一(一sG、/7 , ) ( )I , 式中: ( )为归一化模态, 一为模态响应参数5。 他 的函数。总体来说,当0.5<S。<2时该理论结果 与弹性结构试验结果较为一致。 尾流振子模型具有以下共同特性:(1)振子具 有自激和自限特性。振子自然频率与自由流速成比 例并且满足Strouhal关系,柱体运动与振子耦合; (2)柱体与振子耦合表明柱体运动将强烈影响升力, 同时升力亦将反作用于柱体;(3)模型假定柱体绕 流为2自由度特性(即流场完全相关),因而模型将 不能大响应振幅;(4)这类模型通常没有具体分 析流场,仅用于模拟试验结果,所以这类模型通常 与现象逻辑学相关;(5)建立模型的学者希望得到柱 体振动方程和流体振子方程的形式。并利用二 者共同描述流固耦合系统响应。 2.2力一分解模型(Force—decompostion Models) Sarpkaya提出了力一分解模型[ .弹性支撑刚性 柱体升力分解为与柱体位移相关的流场惯性力,和 与柱体速度相关的流场阻尼力。升力系数 it"r /D\ , CL=Cml,fi2昔( j sin(wt)一 2 — 2 嘉 ( )( )c0s( (5) 式中:c 为惯性系数, 为阻力系数, 为柱体横向 振动周期,Um=2rrA/D,V,=VT/D为折算速度,A为柱 体运动最大振幅, 为来流速度。式(5)可以和弹性 支撑、线性阻尼或周期强迫振动柱体运动方程结 合,柱体运动方程为: 一 . 1 f ,+2 ,+ :JD (C.e sin(f2t))一— ,c cos(O r)(6) J 丌 式中:Xr ̄x[D, 为柱体振动频率与其自然频率之比 ,P,为流体密度与柱体密度之比 。,r=tont。如 果同步状态下 和 可由试验识别,那么式(6) 即可求解,因而柱体同步状态响应亦可求得。 Griifn和Koopman[ 、Griifn[231将流体力分解 为激励项和反馈项,反馈项中包含所有运动相关 力,该模型弹性支撑刚性柱体无量纲方程为: y+2w +6 , =,z ( —CR) (7) 式中: 为升力系数, 为反馈系数,∞ 为Strouhal 频率, 为结构粘性阻尼, 为质量参数为Skop— Griifn参数的倒数。流体动力反馈力与柱体速度反 相且与流速相关,升力与柱体速度同相,向振动柱 体传递能量。 Wang等E ]提出一均匀流下弹性支撑刚性柱体 涡激振动非线性流体力模型。柱体横向振动和顺向 振动通过2自由度结构模型分析,柱体方程由 Euler—Bernoulli梁理论和基于模态分析的模态振动 效应得到,非耦合的运动方程无量纲形式为 一 . 、 X (z,t)+2L (£J,mX z,£)+∞, X ( , )= ( ,t)/2Mr i .. . Y z, )+2 (EJ加Y z,t)+tO ̄o2Y ( ,t) ( ,t)/2Mr Jf (8) 式中: 为柱体轴向, 为结构第n阶模态阻尼, 为结构质量比, 棚为静止柱体自然频率 ( ,f)为 第n阶模态对应横向流体力系数。对于弹性支撑刚 性柱体( =1), ( )=c。( )一C (£)l, ( ),fx (f)=cD( )Y ( )一C (T),其中CD( )为阻力系数,eL(£)为升力系数。 对于静止刚性柱体。横向和顺向流体力系数与升力 和阻力系数相对应。采用2自由度模型是因为竖 向振动对横向振动有一定影响,高次谐波项表示 流固耦合非线性用流体力非线性项表达。模型中 流体力分量由自由振动柱体振幅和频率数据获 得,采用自回归移动平均(ARMA)技术对结构振动 时程谱分析。 3大跨度桥梁涡激振动半经验模型应用 上述半经验模型大多集中于1维竖向振动问 题描述.极少考虑跨向气动力变化及振幅变化情 况,即很少考虑涡激力沿跨向相关性问题。一些研 第6期 朱亚飞,等:大跨度桥梁涡激振动及其半经验模型介绍 ‘43・ 究也仅用1维涡振气动力线性组合的方式粗略地 考虑振动沿跨向变化问题[ ],且其作者本身也承 认仅仅是为重大桥梁结构初步设计提供一些参考 依据。因此目前尚欠缺适用于大跨度桥梁的涡激气 动力相关性函数。 Wilkinson试验研究了方柱涡激气动力相关性. 试验采用刚性节段模型风洞试验,利用不同振幅下测 压获取相关性函数并拟合了涡激力相关性半经验公 式[273。Ehsan结合Scanlan非线性模型和Wilkinson 时采用了耦合方法,主要有时间步法(TMT, Time.Marching Technique)、直接数值模拟法(DNS, Direct Numerical Simulation) 和涡元法(VIC。 Vortex—In.Cell Method)。大多数数值方法通常受限 于雷诺数,但是大涡模拟法(LES,Large Eddy Simulation)已被用于处理雷诺数高达2.4x10 的强 迫振动问题,和中等雷诺数达8 000的自激振动问 题。从数值角度讲,模拟大高宽比复杂外形钝体沿 跨向流致振动仍然受到,故在近期内尚不能实 相关f生函数推导了主梁沿跨向涡激振动振幅公式[28]。 Wilkinson涡激力相关性函数由方柱节段模型 试验获得,若要描述大跨度桥梁主梁沿跨向涡激振 动.首先需要研究主梁断面的涡振特性.由2维涡 振分析扩展到沿跨向涡振分析。通过风洞模型试验 研究模型沿跨向涡振气动力相关性。最终在2维涡 振半经验模型的基础上加人沿跨向涡激气动力相 关性函数,最终形成大跨度桥梁主梁沿跨向涡振分 析体系。 基于半经验模型和风洞试验的大跨度桥梁涡 激振动研究思路为:(1)利用节段模型风洞试验测试 涡激振动响应和涡激气动力;(2)基于几个具有代表 性的半经验模型如经验线性模型、Scanlan经验非线 性模型、升力一振子模型以及Larsen广义非线性模 型,识别主梁断面涡激气动力参数;(3)分析适用于 描述大跨度桥梁主梁涡激振动的半经验模型;(4) 设计适当的主梁气动弹性模型,测试沿跨向涡激 力,并最终拟合得到涡激力相关性函数;(5)最后综 合半经验模型以及涡激力相关性函数,建立大跨度 桥梁沿跨向涡激振动描述体系。 4国内外对涡激振动研究现状及发展趋势 目前国内外对涡激振动的研究可以大致归纳 为解析方法、数值方法、基于风洞试验的半经验模 型方法。 首先,至今没有提出一种完全成功的解析方法 能从基本流动原理出发,表现在旋涡脱落作用下一 个弹性钝体响应动作的全过程。那么理论上要解决 复杂外形钝体涡振更是困难,而简化模型亦会对结 果产生一定影响,所以对于复杂钝体外形得到涡振 解析解难以实现。 其次,数值方法为求解完全耦合钝体涡激振动 问题的另一种途径,数值模拟流致振动必需考虑4 个基本问题:(1)流场模拟;(2)结构振动模拟;(3) 流固耦合作用模拟:(4)数据分析。最近在处理涡振 现主梁沿跨向涡激振动数值模拟。此外.数值方法 虽然近年来随着计算机技术的进步已经取得了很 大进展。但在处理雷诺数达10 量级的3维复杂钝 体桥梁涡激振动问题上仍然不能实现。 最后。已经证实非常有效的办法是文中介绍的 根据风洞试验结果建立半经验模型.通过合理地选 择参数,使模型的特性与真实情况相吻合。上文中 已经详细对各种主要适用的涡激振动模型进行了 介绍,不再赘述。 对于涡激振动研究发展趋势,数值方法应该是 未来研究的主要途径。但这必须基于计算硬件速度 的提高、对气动弹性试验结果的不断丰富以及理论 研究发展不断完善,最终将理论计算结果、数值计 算结果、风洞试验结果以及实际桥梁涡激振动现象 相对照统一,建立起一套完善有效的钝体涡激振动 分析体系。 5结语 由于在大跨度主梁锁定状态下主要发生横风 向振动,大多数情况下可只考虑横风向气动力及其 相关横风振动,称之为单自由度模型。但顺风向也 有风的作用,因此结合顺风向、横风向涡振研究也 十分必要。相应的涡振气动力称为2自由度模型.这 方面的研究相对较少,在为数不多的相关研究中也 是将顺风向、横风向结构运动分离开来进行计算。忽 略其中的耦合现象。而对于实际大跨度桥梁钝体结 构或细长钝体结构而言.涡振现象属于沿跨向范畴 内的3维问题。沿结构跨向方向涡激力并非完全相 关,因此应该将涡振气动力扩展至沿跨向框架下研 究。国外在柱体沿跨向涡脱现象进行了一系列试验。 发现了许多规律性的结果,其中最重要最普遍统一 的结论是:旋涡脱落在沿轴向是不完全相关的,其规 律十分复杂。沿跨向框架下对大跨度主梁涡激振动 研究既十分必要,又存在着相当大的难度。本文从一 些典型大跨度桥梁涡激振动现象出发,介绍了具有 现代交通技木 2010盎 代表l生的涡激振动半经验模型,并为将来大跨度桥 梁涡激振动研究给出了一种研究思路。 参考文献 for vortex..induced osci.. [16]W D 1wan,R D Blevins. A modelnal of Applied Mechanics, 1lation of structures[J] Jour1974,41(3):581—586. l17 IR Land1.A mathematical model orf vortex.excited vibrations [1]陈政清.桥梁风工程[M].北京:人民交通出版社.2005. [2]埃米尔・希缪,罗伯特・H・斯坎伦(刘尚培、项海帆、谢霁 明译),风对结构的作用一风工程导论(第二版)[M].上 海:同济大学出版社.1992. ofbluffbodies[J].Journal ofSound and Vibration,1975,42 (2):219—234. [18]M L Facchinetti,E.de Langre,F Biolley.Coupling of stn—r cture and wake oscillators in vortex—induced vibrations[J3. Journal of Fluids and Structures,2004,19(2):123—140. [3]Strouhal V C.On a particular way of tone generation(in German).Wiedmann s Annalen der Physik und Chemic(new [19]W D 1wan.The vo ̄ex—induced oscillation of elastic stnc.r seires)5,1878:216—251. [4]Zdravkovieh M M.Modiifcation of vortex shedding in the synchronization range[J]_ASME Journal of Fluids En— gineering,1982,104:513-517. [5]Roshko A.On the development ofturbulent wakes from vortex streets[R].National Advisory Committee for Aeronautics, NACA Tech Report 1 191.1954. [6]Berger E.The determination of hydrodynamic parameters of a KdrmOn vortex street from hot.wire measurements at low Reynolds numbers(in Gemran). Zeitsehriftftir Flug— wissenschaften 12,1964:41—59. [7]Fujino Y.Wind—induced vibration and control of Tran—Tokyo Bay Crossing Bridge[J].J Stru.Eng,2002,8:1012—1025. [S I R C Battista et a1.Global analysis of the structural behaviour of the central spans of Rio・Niter6i bridge[R].PONTE SA Contract Report,1993,3. [9]J Upstone,D Reily.Construction ofthe navigation spans ofthe Rio—Niter6i bridge[J].Proceedings of the Institution of Civil Engineers,1979,66(1):227-246. [10]A Larsen,S Esdahl,J E Andresen,et a1.Storebaeh suspension bridge—vortex shedding excitation and mitigation by guide vanes[J].J.Wind Eng.Ind.Aerodyn,2000,(88):283-296. 『1 1]Gtinter Schewe,Allan I.,arsen.Reynolds number effects in the flow around a bluff bridge deck CEOSS section. [12]R E D.Bishop,A Y Hassan.The litf and drag forces on a circular cylinder in a lfowing lfuid[J].Proceedings of the Royal Society Series A.1963.277:32—50. f13]R T Hartlen.I G Currie.Lift.oscillator model of vortex. induced vibration[J]. Journal of the Engineering Mechanics,1970,96(5):577-591. 『14]R A Skop.O M Grimn. A model for the vortex.excited resonant response of bluff cylinders[J].Journal of Sound and Vibration,1973,27(2):225-233. [1510 M Griifn,R A Skop,G H Koopman.The vortex—excited resonant vibrationsofcircularcylinders[J].Journal ofSound andVibration,1973,31(2):235—249. tures[J].Journal of Engineering for Industyr,1975,97: 1378—1382. [20JR A Skop,O M Griifn.On a theory for the vortex—excited oscillations of lfexible cylindrical structures[J].Journal of Sound andVibration,1975,41(3):263—274. [21]T Sarpkaya.Fluid forces on oscillating cylinders[J].Journal of Waterway,Port,Coastal and Ocean Engineering,1978, 104:275—291. [22]O M Griifn,G H Koopman.The vortex—excited 1ift and reaction forces on resonantly vibrating cylinders[J].Journal ofSound andVibration,1977,54(3):435—448. [23]O M Grifin.Vortex・excited cross-lfow vibrations of a single cylindrical tube[J].Journal of Pressure Vessel Technology, 1980,102:158-166. [24]X Q Wang,R M C So,K T Chan.A non—linear lfuid force model ofr vortex-induced vibration of an elastic cylinder[J]. Journal ofSound and Vibration,2003,260(2):287—305. [25]Sangsan Lee,Jae Seok Lee,Jong Dae Kim.Prediction of vortex—induced wind loading on long-span bridges[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1997,67—68:267—268. [26]T Yoshimuraa,Y Mizutaa,F Yamadab,et a1.Prediction of vortex・-induced oscillations of a bridge girder with span・・wise varying geometyr[J].Journal of Wind Engineeirng and Industrial Aerodynamics,200 ,89:1717—1728. [27]Wilkinson,R H.Fluctuating pressures on an oscillating square prism.Part II.Spanwise correlation and loading[J]. Aero.Quarterly,1981,32(2):111-125 『28 1Fazl Ehsan,R H Scanlan.Vortex。Induced Vibrations of Flexible Bridges[J]_Journal of Engineering Mechanics, 1990,1 16:1392-1410. [29]Allan Larsen.A generalized model ofr assessment of vortex— induced vibrations offlexible stnrctures[J].Jouma1 ofWind Engineering and Industrial Aerodynamics,1995(57):281— 294. (收稿日期:2010—05—12)