为了正常的体验网站,请在浏览器设置里面开启Javascript功能!

雷诺数效应对斜拉索风致振动的影响

2012-05-22 6页 pdf 7MB 21阅读

用户头像

is_384009

暂无简介

举报
雷诺数效应对斜拉索风致振动的影响 http://qks.cqu.edu.cn 第33卷第6期 土 木 建 筑 与 环 境 工 程 Vol.33No.6 2011年12月 JournalofCivil,Architectural&EnvironmentalEngineering Dec.2011 雷诺数效应对斜拉索风致振动的影响 刘庆宽,王 毅,郑云飞,马文勇 (石家庄铁道大学 风工程研究中心,石家庄050043) 收稿日期:20110405 基金项目:基金项目:国家自然科学基金(50878135);河北省自然科学基金(E2008000442);...
雷诺数效应对斜拉索风致振动的影响
http://qks.cqu.edu.cn 第33卷第6期 土 木 建 筑 与 环 境 工 程 Vol.33No.6 2011年12月 JournalofCivil,Architectural&EnvironmentalEngineering Dec.2011 雷诺数效应对斜拉索风致振动的影响 刘庆宽,王 毅,郑云飞,马文勇 (石家庄铁道大学 风工程研究中心,石家庄050043) 收稿日期:20110405 基金项目:基金项目:国家自然科学基金(50878135);河北省自然科学基金(E2008000442);河北省科技支撑 (09215626D);教育部新世纪优秀人才支持计划(NECT100130) 作者简介:刘庆宽(1971),男,教授,博士,主要从事桥梁与结构的风荷载、风致振动与控制研究,(Email)lqk@stdu.edu.cn。 摘 要:为了从雷诺数效应的角度斜拉索发生大幅振动的机理,利用不同粗糙度表面无人工水 线的斜拉索模型和光滑表面贴有人工水线的斜拉索模型,进行了测力和测振风洞试验,得到了斜拉 索模型气动力系数、驰振表达式的数值和气动稳定性的结果。通过分析相关的关系表明:在临界雷 诺数区域,力系数及流场的变化特性可能导致斜拉索发生大幅振动,这可能是干索驰振的机理;斜 拉索表面粗糙度能改变雷诺数效应的强弱和临界雷诺数的区域;水线的存在能改变雷诺数效应,一 些水线位置在特定的雷诺数范围内,由于力系数的特殊变化规律导致d犆F/dα为负,斜拉索发生符 合驰振判据的振动;在高雷诺数区域,振动还可能与力系数急剧下降和流场的不稳定有关。 关键词:斜拉索;风致振动;雷诺数效应;水线;阻力系数;升力系数 中图分类号:TU317.1  文献标志码:A  文章编号:16744764(2011)06010606 犚犲狔狀狅犾犱狊犖狌犿犫犲狉犈犳犳犲犮狋狅狀犠犻狀犱犻狀犱狌犮犲犱犞犻犫狉犪狋犻狅狀狅犳犛狋犪狔犮犪犫犾犲狊 犔犐犝犙犻狀犵犽狌犪狀,犠犃犖犌犢犻,犣犎犈犖犌犢狌狀犳犲犻,犕犃犠犲狀狔狅狀犵 (WindEngineeringResearchCenter,ShijiazhuangTieDaoUniversity,Shijiazhuang050043,P.R.China) 犃犫狊狋狉犪犮狋:InordertoanalyzethelargeamplitudevibrationmechanismofcablesbasedonReynoldsnumbereffect, windforcesmeasurementandfreevibrationtestswerecarriedoutoncablemodelswith/withoutartificialwater rivulet.Inthetests,thecablemodelswithoutartificialwaterrivuletwereindifferentsurfaceroughnessandthose withartificialwaterrivuletattachedwereinsmoothsurface,andartificialwaterrivuletwasattachedindifferent positionsrespectively.Windforcecoefficients,DenHartoggallopingcoefficientsandfreevibrationamplitudeswere obtained.ItisshownthatincriticalReynoldsnumberrange,thereisthepossibilitythatthespecialwindforce coefficientcharacteristicsandflowpatterninducevibration,andthisisthepossiblemechanismfordrycable galloping.TheroughnessofcablesurfacecanchangetheintensityofcriticalReynoldsnumbereffectandtherange ofcriticalReynoldsnumber.TheexistenceofwaterrivuletcanchangethecriticalReynoldsnumbereffect.For somewaterrivuletpositions,incertainReynoldsnumberrange,犱犆犉/犱α<0,vibrationhappens,whichagreeswith theDenHartoggallopingcriterion.InhighReynoldsnumberrange,besidesDenHartoggalloping,vibrationsare relevanttosignificantdecreaseofwindforcecoefficientsandunstableflowpattern. 犓犲狔狑狅狉犱狊:cables;windinducedvibration;Reynoldsnumbereffect;waterrivulet;dragforcecoefficient; liftforcecoefficient   由于斜拉桥斜拉索的长细比大、阻尼低等特点, 导致在大风或风雨联合作用下经常发生大幅振动, 对桥梁结构的安全造成危害。与传统的振动类型相 比,风雨激振和干索驰振(DryCableGalloping)由 http://qks.cqu.edu.cn 于振幅大、破坏严重,是值得深入研究的问。关于 风雨激振,自20世纪80年代被发现以来,各国学者 通过现场观测与风洞试验[14]、理论分析[57]及CFD 等手段已经进行了较广泛的研究,针对其振动机理 提出了水线驰振理论、高风速涡致振动理论、轴向流 理论、弯扭两自由度耦合理论、水线摆动理论、卡门 涡被抑制导致大幅振动等理论,建议采用的气动抑 振措施有缠绕螺旋线和索的表面设置凹坑[8]等。 在斜拉索风雨振的现场观测[910]、足尺模型观 测[11],风洞试验[1213]研究中,研究人员发现在没有 降雨、或者降雨量很小不足以形成水线、或降雨已经 停止的状态下,斜拉索也可能发生大幅振动,尤其是 日本Sunbridge桥的斜拉索发生大幅振动[10]以来, 这种振动现象及其机理引起了高度重视。目前 Cheng 等[14] 利 用 Den Hartog 驰 振 理 论、 Matsumoto[15]利用卡门涡的脱落抑制理论等对干索 驰振的机理进行了解释,但是无论其发生的机理、还 是与风雨激振的内在关系,以及抑振措施等,都需要 进一步的深入研究。 该文从雷诺数效应入手,通过测力和测振风洞试 验,以不同雷诺数下的阻力系数、升力系数和振幅为 参数,研究了雷诺数效应是否导致振动、不同粗糙度 表面模型的雷诺数效应及与振动的关系、水线影响雷 诺数效应及通过影响雷诺数效应导致振动的机理。 1 试验介绍 为了研究斜拉索的雷诺数效应和气动稳定性的 关系,共进行了2大类的风洞试验,其一是两端固定 刚性模型的测力试验,其二是两端弹簧支撑刚性模 型的测振试验,试验对象分别是表面粗糙度不同、没 有人工水线的斜拉索模型和表面光滑的贴有人工水 线的斜拉索模型。试验在石家庄铁道大学风工程研 究中心的双试验段回/直流大气边界层风洞内进行, 其低速试验段转盘中心宽4.4m,高3.0m,长 24.0m,最大风速大于30.0m/s,背景湍流度犐≤ 0.4%;高速试验段宽2.2m,高2m,长5.0m,最大 风速大于80.0m/s,背景湍流度犐≤0.2%。风洞结 构如图1所示[16]。本试验在高速段内进行。 图1 风洞平面图   为了实现不同雷诺数下的雷诺数效应,使用了 4个刚性斜拉索模型,具体参数如表1所示。 表1 模型参数 模型编号 直径犇/mm 质量犿/(kg·m-1)长度犔/m 表面状态 M1 150 7.008 1.910 光滑无水线 M2 150 7.083 1.910 小粗糙无水线 M3 151 7.354 1.910 大粗糙无水线 M4 150 7.125 1.910 光滑有水线 模型的材质为有机玻璃,原型为表面光滑的圆 柱,两端设置端板,由中间贯穿的刚性圆管支撑在风 洞两侧的支架上。M1模型为光滑表面,M2模型是 将光滑圆管用 P24号砂纸均匀打磨而成,M3模型 是将光滑圆管用表面粗糙的壁纸包裹而成,M4模 型是在 M1模型的基础上分别在不同的位置粘贴人 工水线而成。对于人工水线,利用有机塑料加工成 圆弧外形,粘贴在斜拉索的表面。斜拉索表面和水 线的形状如图2所示(为了对比表面的粗糙度,拍照 时模型表面放置了最大直径约7mm铅笔)。水线 的位置用θ表示,是从前驻点到水线中心转过的圆 心角。 测力模型为两端固定支撑,端部安装美国 ATI 公司生产DELTA系列六分力高频天平,测试采样 时间为60s。 测振模型两端分别用4根弹簧支撑,弹簧刚度 的选取以系统的振动频率与实际相同直径斜拉索的 振动频率一致为原则。为了便于起振,模型系统的 犛犮数比实际斜拉索的稍小。振动过程中记录瞬态 位移。 701第6期 刘庆宽,等:雷诺数效应对斜拉索风致振动的影响 http://qks.cqu.edu.cn 图2 模型表面状态和水线尺寸 测力和测振模型的空间位置为:竖直倾斜角 α=0°,水平倾斜角β=0°,即斜拉索模型在水平面 内,与来流风向垂直。 试验中的控制风速由安装在试验段入口的传感 器与控制台组成的稳风速控制系统控制,模型处的来 流风速澳大利亚TurbulentFlowInstrumentation公 司生产的4孔眼镜蛇探头(4holeCobraProbe)测 试,安装位置为模型中心上游1.05m、下方0.47m 处,采样频率2000Hz,测力的采样时间60s,测振 的记录时间为整个振动观察的时间范围。 因为不同粗糙度模型对应的临界雷诺数不同, 所以4个模型的试验雷诺数范围不同,并且为了准 确反应力系数的变化情况,针对各个模型的临界雷 诺数区域分别加密了测试工况。 水线位置从10°开始,以2.5°为步长增加到70°。 为了从机理上进行研究,采用的水线位置比实际可 能形成水线的位置范围要大。3个无水线模型对应 的试验雷诺数范围、步长如表2所示,有水线模型的 试验工况如表3所示。 表2 无水线模型(犕1—犕3)测力和测振试验工况 M1 犚犲范围/ 万 犚犲步长/ 万 M2 犚犲范围/ 万 犚犲步长/ 万 M3 犚犲范围/ 万 犚犲步长/ 万 12~30 2.0 12~26 2.0 12~22 2.0 30~35 1.0 26~30 1.0 22~27 1.0 35~42 0.5 30~38 0.5 27~34 0.5 42~51 1.0 38~47 1.0 34~46 1.0 表3 有水线模型(犕4)测力和测振试验工况 测力试验 水线位置/ (°) 犚犲范围/ 万 犚犲步长/ 万 测振试验 水线位置/ (°) 犚犲范围/ 万 犚犲步长/ 万 10~70 8~15 1.0 10~50 8~15 1.0 17~47 2.0 17~47 2.0 52.5~70 8~21 1.0 23~47 2.0 2 无水线模型雷诺数效应对气动稳定 性的影响 21 无水线模型的雷诺数效应 使用3个不同粗糙度模型测得的平均阻力系数 和平均升力系数随雷诺数的变化曲线如图3所示。 由图可知,在临界雷诺数区域,模型的平均阻力系数 下降,平均升力出现。平均升力开始出现时的雷诺 数,基本对应平均阻力系数开始下降时的雷诺数;平 均阻力系数大约下降到整个下降幅度一半的时候, 平均升力系数取得最大值,之后随着雷诺数的增大 平均升力系数开始减小,当平均阻力系数下降到最 小值的时候,平均升力系数基本恢复到零值。 图3 平均阻力系数和平均升力系数曲线   随着模型表面粗糙度的增加,平均阻力系数的 下降幅度减小,平均升力系数的最大值减小,即雷诺 数效应减弱。同时,随着粗糙度的增加,临界雷诺数 区域整体向低雷诺数方向移动了一定数值,即在较 小的雷诺数数值时就进入了临界雷诺数的状态。 22 无水线模型的气动稳定性 无水线模型的测振结果如图4所示。对照图3 可以发现,各个模型在亚临界雷诺数区域,振动的振 幅都很小,可以认为是稳定的;雷诺数到达临界区域 时,振幅显著增大,光滑斜拉索模型的最大振幅达到 801 土 木 建 筑 与 环 境 工 程                 第33卷 http://qks.cqu.edu.cn 了12.3cm,0.82犇(其犇为斜拉索的直径),超过临 界雷诺数区域之后,振幅又下降至很小,基本可以认 为是恢复到了稳定状态。最大升力系数出现时的雷 诺数,对应最大振幅。随着模型表面粗糙度的增加, 在临界区域发生振动的振幅减小。结合图3综合分 析可以推知,在临界雷诺数区域,由于阻力系数的减 小和升力的出现、流场从卡门涡周期脱落状态变为 不规则的状态等因素,导致了振动的发生,但是振动 发生具体的机理尚待流场分析等进一步的研究。 图4 振幅与雷诺数曲线 3 有水线模型雷诺数效应对气动稳定 性的影响 31 水线对雷诺数效应的影响 对粘贴有人工水线的斜拉索模型的阻力系数、 升力系数、自由振动振幅分别进行了测试,利用阻力 系数和升力系数,计算如(1)所示的表达式的 值,即DenHartog驰振准则。如果值为负,则可能 发生驰振。 d犆F dα= d犆L dα+犆D (1) 分析表明:随着雷诺数从小到大的变化,水线位 置不同,阻力系数、升力系数的变化也不同,相应的 d犆F/dα的值和自由振动的状态也不同。选取水线 位置θ=15°、25°、55°3个有代表性的工况进行分 析,其阻力系数和升力系数分别如图5—7所示。 由图5可知,水线位置θ=15°时,阻力系数从 犚犲=17万左右开始随着犚犲的增长呈现阶段性下降 趋势,其中25万至37万之间基本保持不变,从37 万开始急剧下降,一直到犚犲=43万左右降到最低; 与此对应,升力系数从犚犲=17万左右开始上升,在 犚犲=25万和37万之间保持较大值,之后急速下降, 到犚犲=43万左右基本降到最低。 从犚犲=17万到43万之间阻力系数下降、出现 较大升力的现象,与光滑模型的临界雷诺数区域的 特征一致,可以判断这个区域为临界雷诺数区域。 对照图3(a)的结果可知,2个工况临界雷诺数区域 结束时的雷诺数(43万)基本一致,而水线的存在, 大大提前了临界雷诺数区域开始时的雷诺数数值 (无水线时35万,15°水线时17万)。 图5 水线位置15°的阻力系数和升力系数 图6 水线位置25°的的阻力系数和升力系数 比较25°和15°时的阻力和升力系数(图6、图5) 可知,25°水线时临界雷诺数区域的升力系数整体比 比15°时的小(25°水线时升力系数最大值在1.3~ 1.1之间,15°的在1.0~0.5之间),相同阶段的阻力 系数比15°时的大,临界雷诺数效应开始时的雷诺数 数值比15°时的小,力系数降到最低时的雷诺数(41 万)比15°时(43万)的小。 901第6期 刘庆宽,等:雷诺数效应对斜拉索风致振动的影响 http://qks.cqu.edu.cn 图7 水线位置55°的的阻力系数和升力系数 比较55°水线时的情况可知,在犚犲=10~35万 范围内,升力系数基本为零值,阻力系数保持在1.5 左右,从图上已经看不到升力系数上升阻力系数分 阶段下降的情况。从35万开始,2个力系数急剧下 降,39万左右基本降到最低。 综上,可以总结出力系数随水线升高的变化特 征为:力系数急剧下降基本出现在犚犲=35万,随着 水线位置的升高,力系数降到最低时的雷诺数数值 逐渐减小,临界雷诺数区域对应的升力系数逐渐减 小,阻力系数逐渐增大,临界雷诺数效应发生的范围 扩大。即水线位置不同,阻力系数、升力系数的大 小、变化规律和雷诺数效应不同。 32 水线对气动稳定性的影响 粘贴有人工水线的斜拉索模型各个雷诺数下的 d犆F/dα值和测振结果如图8—10所示。 随着水线位置的变化,阻力系数和升力系数的 变化有可能导致d犆F/dα<0,发生驰振。 由15°水线的d犆F/dα和振动图(图7)可知,犚犲 在27万到28万之间,d犆F/dα由正变负,在该雷诺 数附近模型开始发生大幅振动,d犆F/dα为负的区域 与发生大幅振动的区域基本一致。 水线位置为25°时,虽然从 犚犲=15万之后 d犆F/dα一直为负值,但是大幅振动却分别发生在2 个雷诺数区域,一个是犚犲=15~25万之间,另一个 是犚犲=37~42万之间。前一个区域的发生机理,可 能是单纯由于水线位置上下变动时阻力系数、升力 系数的数值不同,因而d犆F/dα<0引起,后一个区域 发生振动的机理,结合无水线时力系数突然下降时 发生大幅振动的结果(图3(a)和图4(a)),推断除了 d犆F/dα<0之外,是否可能与力系数的急剧下降和 流场的不稳定有关,该机理尚需要进一步研究才能 明确。 水线位置为55°时发生振动的状况同25°时类 似,也是发生在2个区域,只不过在力系数急剧下降 时,发生大幅振动的区域变窄。2个区域振动的机 理分析同25°工况。 图8 水线位置15°的犱犆犉/犱α和振动振幅 图9 水线位置25°的犱犆犉/犱α和振动振幅 图10 水线位置55°的犱犆犉/犱α和振动振幅 4 结论 利用不同粗糙程度无人工水线的斜拉索模型和 光滑表面贴有人工水线的斜拉索模型,通过测力和 测振风洞试验,得到了斜拉索模型气动力、气动稳定 性的结果。通过分析相关的关系,得到了以下结论: 1)在临界雷诺数区域,力系数及周围流场的变 化特性可能导致斜拉索发生大幅振动,这可能是干 索驰振的机理。 011 土 木 建 筑 与 环 境 工 程                 第33卷 http://qks.cqu.edu.cn 2)随着斜拉索表面粗糙度的增加,临界雷诺数 区域整体向小的方向移动,雷诺数效应减弱,在临界 雷诺数区域发生振动的振幅减小。 3)水线的存在能改变雷诺数效应,一些水线位 置在特定的雷诺数范围内,由于力系数的特殊变化 规律导致d犆F/dα为负,斜拉索发生了符合驰振判据 的振动。 4)在高雷诺数区域,斜拉索发生振动的机理除 d犆F/dα为负之外,还可能与力系数急剧下降及流场的 不稳定有关,该部分机理尚需进一步的研究进行明确。 参考文献: [1]MACDONALDJHG,LAROSEGL.Twodegreeof freedom inclined cable gallopingPart 1: General formulationandsolutionforperfectlytunedsystem[J]. Journal of Wind Engineering and Industrial Aerodynamics,2008,96(3):291307. [2]MACDONALDJHG,LAROSEGL.Twodegreeof freedominclinedcablegallopingPart2:Analysisand preventionforarbitraryfrequencyratio[J].Journalof WindEngineeringandIndustrialAerodynamics,2008, 96(3):308326. [3]刘庆宽.斜拉桥斜拉索风雨振时索表面水线摆动作用及规 律的试验研究[J].土木工程学报,2007,40(7):6267.   LIUQINGKUAN.Experimentalstudyonmovement ofwaterrivuletoncablesurfaceinrainwindinduced vibrationofstaycables[J].ChinaCivilEngineering Journal,2007,40(7):6267. [4]李永乐,卢伟,陶齐宇,等.斜拉桥拉索风雨振动特性风 洞试验研究[J].试验流体力学,2007,21(4):3640.   LIYONGLE,LUWEI,TAOQIYU,etal.Studyon rainwindinducedvibrationofcablesincablestayed bridgesby windtunneltest[J].Experimentsand MeasurinFluidMechanics,2007,21(4):3640. [5]李寿英,顾明,陈政清.运动水线三维连续弹性拉索风 雨激振理论模型[J].湖南大学学报:自然科学版,2009, 36(2):17.   LISHOUYING,GU MING,CHENZHENGQING. Ananalyticalmodelofrainwindinducedvibrationof threedimensionalcontinuousstaycable with actual movingrivulet[J].Journalof Hunan University: NaturnalScience,2009,36(2):17. [6]张琪昌,李伟义,王炜.斜拉索风雨振的动力学行为研 究[J].振动与冲击,2010,29(4):173176.   ZHANG QICHANG,LI WEIYI, WANG WEI. Nonlinear dynamic behavior of rainwindinduced vibrationofastaycable[J].JournalofVibrationand Shock,2010,29(4):173176. [7]李暾,陈政清,李寿英.连续弹性拉索风雨激振理论模 型研究[J].振动工程学报,2010,23(4):380388.   LITUN,CHEN ZHENGQING,LISHOUYING. Analyticalstudyofrainwindinducedcablevibrations basedoncontinuouselasticcablemodel[J].Journalof VibrationEngineering,2010,23(4):380388. [8]刘健新,李哲.气动措施对斜拉索风荷载及结构响应的 影响[J].建筑科学与工程学报,2010,27(3):8993.   LIUJIANXIN,LIZHE.Influenceofaerodynamic measureonwindloadandstructuralresponseofstayed cable [J]. Journal of Architecture and Civil Engineering,2010,27(3):8993. [9]VIRLOGEUX M.Cable vibrationsin cablestayed bridges [C ]//Bridge Aerodynamics. Balkema: Rotterdam,1998:213233. [10]MATSUMOTO M.Therole of waterrivuleton inclinedcableaerodynamics[C]//Proceedingsofthe6th AsiaPacificConferenceon WindEngineering.Seoul, Korea,2005:6377. [11]MATSUMOTO M,SHIRATO H,YAGIT,etal.Field observationofthefullscalewindinducedcablevibration [J].Journal of Wind Engineering and Industrial Aerodynamics,2003,91(1/2):1326. [12]MATSUMOTOM,YAMAGISHIM,AOKIJ,etal. Various mechanism ofinclined cable aerodynamics [C]//Proceedingof9thInternationalConferenceon WindEngineering,NewDelhi,India,1995:759770. [13]CHENGS,LAROSE G L,SAVAGE M G,etal. Experimentalstudyonthewindinducedvibrationofa dryinclinedcablePartI:Phenomena[J].Journalof WindEngineeringandIndustrialAerodynamics,2008, 96(12):22312253. [14]CHENGS,IRWINPA,TANAKAH.Experimental studyonthewindinducedvibrationofadryinclined cablePartII:Proposed mechanisms[J].Journalof WindEngineeringandIndustrialAerodynamics,2008, 96(12):22542272. [15]MATUMOTO M,YAGIT,HATSUDA H,etal. Drygalloping characteristicsandits mechanism of inclined/yawedcables[J].JournalofWindEngineering andIndustrialAerodynamics,2010,98(6/7):317 327. [16]刘庆宽.多功能大气边界层风洞的与建设[J].实 验流体力学,2011,25(3):6670.   LIUQINGKUAN.Aerodynamicandstructuredesigh of multifunction boundarylayer wind tunnel[J]. JournalofExperimentsinFluidMechanics,2011,25 (3):6670. (编辑 胡英奎) 111第6期 刘庆宽,等:雷诺数效应对斜拉索风致振动的影响
/
本文档为【雷诺数效应对斜拉索风致振动的影响】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。 本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。 网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。

历史搜索

    清空历史搜索