维普资讯 http://www.cqvip.com 第25卷第2期 空气动力学学报 Vo1.25.No.2 2007年o6月 文章编号:0258.1825(2007)02—0169—07 ACTA AERoDYNAMICA SINICA Jun.,2007 斜拉桥拉索风雨激振理论模型和机理研究 顾 明,李寿英,杜晓庆 (同济大学土木工程防灾国家重点实验室,上海200092) 摘要:斜拉桥拉索风雨激振的机理研究是国际风工程和桥梁工程领域的著名难题。本文设计、制作了可方便调 节拉索模型倾角和风向角的试验装置以及用于测压试验的带人工雨线的拉索模型,并在风洞中进行了细致的试 验。试验得到了具有典型倾角的拉索在不同风向角下,拉索和水线模型上的平均风压和脉动风压系数,以及气动 力系数。在此基础上,建立了拉索风雨激振新的理论模型,计算分析了节段三维拉索的风雨激振响应及其机制。 关键词:斜拉桥拉索;风雨激振;气动力系数;风洞试验;理论模型;机制 中图分类号:TU312 文献标识码:A 0 引 言 1998年大跨桥梁空气动力学国际学术会议上将 斜拉桥拉索风雨激振的机理研究确定为较长一段时 间内国际风工程和桥梁工程领域的4个重点研究问 拉索模型的测压试验研究,得到了带人工雨线三维拉 索模型上的压力系数和气动力系数(包括拉索和水线 上的气动力)。在此基础上.建立了拉索风雨激振的 理论模型,计算了节段三维拉索的风雨激振响应,并 分析了机制。 题之一。斜拉桥拉索在风雨共同作用下常发生剧烈 振动,国内外对此已有很多报道【1-5 J。在一定的雨量 和8m/s~18m/s左右风速下,有些拉索的振幅可达1 1试验装置及试验工况 试验在同济大学TJ一3大气边界层风洞的均匀流 风场中进行。该风洞试验段尺寸为宽15m、高2m、长 14m。风速范围从0.2m/s~17.6m/s连续可调。用皮 米以上,为其直径的十倍左右。拉索的大幅振动造成 了多座斜拉桥拉索的损坏,数座大桥不得不更换部分 甚至全部拉索,造成了重大的经济损失。 十多年来,研究人员认识到雨线(主要是上雨线) 是产生拉索风雨激振的根本原因,并在风洞中进行了 一托管和微压计测量风速。风压测量及数据处理系统 为美国Seanivalve扫描阀公司的量程为±254mm和± 508mm水柱的DSM3000电子式压力扫描阀系统、Pc 机、以及自编的信号采集及数据处理软件。 为了测得作用在拉索表面上雨线上的气动力,需 要在上雨线表面布置足够多的测压点。但从实际观 测和风洞试验结果可知[11-13】,拉索上雨线的尺寸很 小。因此,本试验采用放大拉索模型直径因而也放大 人工雨线的方法。设计人工雨线的形状和尺寸见图 1。拉索模型采用有机玻璃材料。拉索模型直径为 350mm,约为实际拉索直径的2~3倍。雨线尺寸也 放大2~3倍。模型全长3.5m。为了保证模型和实 际拉索的 数相同(约为7.04 X 10 ~2.35 X 105). 试验风速设定在3m/s~10m/s,约为实际发生拉索风 雨激振风速的1/2~1/3。 系列带人工雨线拉索的气动力试验研究【∞】,为解 释其机理以及开展初步理论分析【6,8-10】提供了依据。 但这些测力试验所采用的拉索模型均是二维模型(即 拉索的倾角和风向角均为0。),因而理论模型中采用 的气动力也是二维的[ .8. ]。实际拉索发生风雨激 振时,由于风向角不为零,拉索与来流之间的相对位 置是三维空间关系,二维拉索模型获得的气动力不能 准确反映三维拉索的气动力特性,二维理论模型也不 能完全反应实际拉索风雨激振的特征。此外,过去的 研究中,由于没有测量水线气动力,所以在理论模型 中假设了水线的运动规律,而没有建立合理的水线运 动方程。 针对以上问题,本文首先进行了带人工雨线三维 ・收稿日期:2006-03.18;修订日期:2006-06—22. 基金项目:国家自然科学基金创新研究群体基金资助. 作者简介:顾明(1957一),男,江苏兴化人。教授,长江学者。从事结构抗风研究.Email:minggu@mail.tongii.edu.en 维普资讯 http://www.cqvip.com 维普资讯 http://www.cqvip.com 第2期 顾 明等:斜拉桥拉索风雨激振理论模型和机理研究 l7l 拟合出了拉索上的气动力公式 ’引,用于建立理论模 型。限于篇幅,气动力公式这里从略。 一平均风压系数一脉动风压系数 /,, \ / ~平均风压系数一均方根风压系数 j E 《 (c)日t=46 ̄,%=246。 图5带人工雨线拉索模型表面风压系数的分布 及典型测点的风压功率谱 Fig.5 Wind pressure coefficients and PSD of wind pressure at typical point of model (3)水线上的气动力系数 对上雨线表面测得的风压系数沿雨线表面积分, 得到上雨线所受的气动力。图7给出了风向角为35。 时, 数为1.2e5、倾角30。的拉索的上雨线的气动力 系数Cd和Cf随上雨线位置不同而变化的情况。由 图可见,水线上的气动力系数随位置也有较大变化。 根据试验结果,同样拟合出了水线上的气动力公 式【l 4l。 图6拉索的气动力系数 Fig.6 CD and CL of cable 图7上水线上的气动力系数 Fig.7 Aerdoyrmnficfo ̄'e coeficients of rivulet 3斜拉索风雨激振理论模型 根据现场实际观测和风洞试验结果,为建立斜拉 索风雨激振理论模型,提出如下基本假定: (1)水线在运动过程中外形和大小保持不变;(2) 由于拉索涡激共振风速远离风雨激振风速,故准定常 假设适用;(3)水线与拉索之间的阻尼力包括库仑阻 尼力和线性阻尼力;(4)拉索仅发生面内振动,不考虑 面外振动;(5)忽略下水线的影响,仅考虑上水线的影 响。如无特殊说明,本文水线均表示上水线。斜拉索 风雨激振力学模型示意图见图8。 坐标系oxy, 轴水平向右,Y轴竖直向下。风速 为 ,拉索倾角a和风向角 (见图3),来流与拉索 的相对风速为 , ,与水平方向的夹角为 。 拉索刚度为 ;阻尼系数C ;单位长度质量为 。在Y轴向上,拉索受到的荷载有回复力KvY、阻 尼力c y、惯性力M 、气动力Fy、重力 。此外,忽 略水线运动惯性力和水线对拉索的作用力。Y轴方 向上拉索的运动微分方程如下: .. . Y+2 +∞2∥=一 F (1) 』rJ 式中,位移Y为拉索离开静平衡位置的距离; 为拉 维普资讯 http://www.cqvip.com 172 空气动力学学报 第25卷 索模型的阻尼比;∞,为拉索模型的圆频率; 的正 向为Y轴的负方向: Fy= 1 作为计算对象。拉索直径为120mm;单位长度拉索质 量为6kg/m;拉索固有频率为1.0Hz;拉索阻尼比为 0.1%;单位长度水线质量为0.01kg/m;水线与拉索 之间的线性阻尼系数C,设为0.0008;单位长度水线 IDu。2D(cLc。s +CD¥in ) (2) 其中 为 为拉索的平均气动力;』D为空气密度; 上的库仑阻尼力F0取为单位长度水线重力的20%, D为拉索直径;气动力系数采用本文前述试验结果。 ‘ , 磊 图8斜j立累风雨激振力学模型 Fig.8 Mechanical model for rain.wind induced vibration of stay cables 水线上在拉索表面形成并做周向振荡时,在拉索 表面的切向方向上,作用在水线上的荷载有:水线与 拉索表面之间的阻尼力F0+c 、水线自身运动产 生的惯性力mR 0、拉索运动产生的附加惯性力m ycos0、水线重力rng COSGcosO、气动力 。其中 为 拉索截面半径。则在拉索表面切向方向上水线的运 动微分方程为: mR 0+F0+c 0= +mycos0一mg COSGcosO (3) 水线上的气动力 由下式确定: I = IDu02 (Cdsin( + )+Czcos( + )) (4) 为水线的圆弧半径;Cd、Cz为水线上的平均气动力 系数,已在前文给出。特别要说明的是,本文提出水 线与拉索表面之间的阻尼力为库仑阻尼力 和线 性阻尼c 之和,这和以往仅考虑线性阻尼是不同 的,也更符合实际。 4拉索风雨激振响应和机理分析 选取和文献[12]试验模型具有相同参数的拉索 即0.02N/m;水线初始位置0=1.0tad及初速度为 零;拉索倾角为30。;风向角为35 ̄。 首先分析拉索和水线的响应随风速的变化。计 算结果表明,当风速小于7.3m/s时,水线不能在拉索 表面形成,拉索振幅很小。风速在7.6m/s~9.0m/s 范围内时,水线的振动比较规则,水线平衡位置都在 0=46.0。左右,且随着来流风速的增大。在这个风速 范围内,拉索的振幅都较大。当风速从9.0m/s增大 时,水线平衡位置突然从48.4。增大到59.0。,拉索振 幅也突然变得很小。图9是两种典型风速(7.2m/s 和8.0m/s)时拉索和水线的振动时程。 60 50 40 30 20 20 0 20 40 60 80 100120140l60 l0 20 0 20 40 60 80 100 120140160 rime/s Tim ̄/s (a)Uo=7.2m/s ■一(b)uo=8.0m/s 图9拉索(左)和水线(右)的振动时程 Fig・9 Time histories of cable(1ift)and rivulet(r;ght) 图10给出了拉索振幅与风速的关系,图中同时 给出了文献[12]通过人工降雨试验得到的拉索的振 幅结果。由图可见,拉索在7.3m/s~9m/s之间有较 大振幅,反映了风雨激振“限速”的特点。本文选用的 基本参数与文献[12]进行人工降雨试验的参数基本 相同,计算得到的起振风速范围与试验值基本相同, 振幅也具有很好的可比性。 维普资讯 http://www.cqvip.com 第2期 顾 明等:斜拉桥拉索风雨激振理论模型和机理研究 l73 图11是典型风速(8.0m/s)时拉索和水线振动响 O O 2 2 ● ● 5 O 5 O O 应的功率谱分析结果。由图可见,拉索振动频率主要 以固有频率为主(1Hz),而水线振动频率除拉索固有 频率外,还包括能量较低的倍频成分。 E \ 一 《 《 口 暑 =A E 《 6 7 8 9 lO l1 l2 Uo/rigs 图lO拉索振动幅度与风速的关系 Fig.10 Vibration amplitudes of cable vs.wind velocity O.O3 U —O O2 }0.01 《 O 0 0.5 1.0 1.5 2.O 2.5 3.0 3.5 4 0 Fmquency/Hz (a)拉索 0 j .兰 一 U : 一△ E 《 Frequency/Hz (b) 水线 图11 拉索和水线振动响应的功率谱( :8.Ore/s) Fig.1 1 PSD of displacement of cable and rivulet 根据计算结果,我们对拉索风雨激振的过程作出 如下解释:水线初始位置如果在“不稳定区域”(即水 线在震荡中将运动到“危险区域”的区域)以外时,拉 索和水线不会发生不稳定振动;而当水线的初始位置 位于“不稳定区域”内时,外激励引起拉索的微小振动 将使水线沿拉索周向振荡并加剧,当水线在某一瞬时 运动到“危险区域”(即驰振不稳定区)时,气动力对拉 索做正功且频率相同,导致拉索和水线的振动幅值的 突然增加,也即拉索风雨激振发生。作者在文献[8] 中分析二维拉索风雨激振特性是得到相同的结论。 以上讨论可能解释拉索风雨激振的机制。 5 结 语 通过对带人工雨线的拉索模型的风洞试验,得到 了具有典型倾角的拉索在不同风向角下拉索和水线 模型上的气动力系数。在此基础上,建立了拉索风雨 激振新的理论模型,计算分析了节段三维拉索的风雨 激振响应和机制。基于本文理论模型的计算结果和 人工降雨条件下三维拉索刚性节段模型的振动试验 结果吻合。研究表明,拉索气动力的突变区域对应于 拉索风雨激振的大幅振动区域,这类似于经典驰振。 但由于水线的运动,拉索风雨激振的不稳定区范围远 大于驰振不稳定区。 参 考 文 献: [1]HIKAMI Y,SHIRAISHI N.Raln-wind induced vibrations of cables in cable stayed bridges[J].Journa/of E, ̄si— neering andIndustrialAerodynamics,1988,29:409—418. 【2 J MAIN J A,JONES N P.FuU—scale measurements of stay ca— ble vibration【J]. Engineering into the 21 st Century, 1999,Balkema,Rotterdalu,963.970. [3]顾明,刘慈军,罗国强,项海帆,林志兴.斜拉桥拉索的 风(雨)激振及控制[J].上海力学,1998,19(4):281— 288. [4]陈政清,柳成荫,倪一清,王修勇,伏晓宁.洞庭湖大桥 拉索风雨激振中的风场参数[J].铁道科学与工程学 报,2004,1(1):52.57. 【5 J PERSOON A J,NOORLANDER K.Full-scale measurements on the Erasmus Bridge after rain/wind induced cable vibra— tions【J].阢,ld Engineering into the 21 st Cetnury,1999, Balkema,Rotterdam,1019—1026. [6] YAMAGUCHI H.Analyticla study on growth mechanism of rain vibration fo cables[J].Journal foWind Engineering and Industrial Aerodynamics,1990,33:73—80. 【7] GU MING,LIU CIJUN,XU YOULIN,XIANG HAIFAN. Response characteristics of wind excited cables with artiifcial irvulet【J].Applied Math.and Mech.,2002,23(1O): 1176.1187. 【8] GU MING,LU QIANG.Theoreitcal analysis of wind—rian in— duced ivbration fo cables fo cbale—stayed bridges[J].Journa/ fo W/nd Eng/ener/ ̄,2001,89:125—128. 维普资讯 http://www.cqvip.com l74 [9] 空气动力学学报 第25卷 XU YOULIN,WANG UNYAN.Analytical study of wind・ wind-induced vibration of cables in cable—stayed bridges and rain-induced cable vibration:SDOF model[J].Journal of Wind Engieerin ̄27—40. its miitgation[J].Journal f oAerodynamics.2005,93:79—95. Engineering and Idunstrial and Industrial Aerodynamics,2003,91: n—wind induced vibration of cables[J]. [13] FLAMAND O.RaiW.Simple model of rain—wind-in— [1O] WILDE K,WITKOWS ̄Journal of Wind Engineering and Idunstrial Aerdynoamics, 1995,57:353—362. duced vibrations of stayed cables[J].Jouma/of Wind Engi- eerni ̄andIndustrial Aerodynamics,2003,91:873・891. MATSUMOTO M,et a1.Response characterlsdes of rain—wind [14] 顾明,杜晓庆.不同风向角下斜拉桥拉索模型测压试 验研究[J].振动与冲击,2005,24(6):5—9. [15] 杜晓庆,顾明.斜拉桥拉索上雨线的气动力特性的试 验研究[J].同济大学学报,2005,33(5):585.589. [16] 顾明,杜晓庆.带固定人工雨线的斜拉桥拉索模型测 压试验研究[J].空气动力学学报,2005,(4):419—425. induced vibration of stay—cable of cable-stayed bridges【J]. ourJnal of Wind D咖1995,57:323—333. 而 and Idunstrial Aerodynamics, [12] GU MING,DU XIAOQIN.Experimental investigation of rain一 Testing study on wind pressure distributions of stayed cables with a rLxed artiicialf rivulet GU Ming,LI Shou—ying,DU Xiao-qing (State Keylaboratory,or Disaster Reduction in Civil Engineering, ,j University,Shanghai 200092,China) Abstract:Rain・・wind induced vibration of cables of cable・・stayed bridges is a problem of great concern in communities of wind engineering and bridge engineering.In this paper.a cable model’vith an artiifcial rivulet is tested with the help of a newly designed test set-up on which the cable model Can be easily adjusted to required inclined angle and wind angle. Mean and fluctuating pressures and aerodynamic force coefficients of the cable model and the artiiciafl rivulet under the con- ditions ftoypical cable inclined nglaes and wind angles arc obtained.On the basis ofthe test I℃sults.a new theoretical mod- e1 for describing rain- rind induced vibration of stayed cables is established.with which the rain—wind induced responses a typical cable model are computed and the mechanism is analyzed. Key words:cable of cable-stayed bridge;rain-wind induced vibration;aerodynamic force coefficient;wind tunnel test;theoretical model;mechanism