摘要
以典型流线型闭口箱梁为研究对象,利用节段模型风洞试验,获取原始断面与优化断面(栏杆扶手抑流板断面和检修轨道导流板断面)涡振响应,并分析其典型风速下断面周围分布气动力矩贡献,结合简化涡和数值模拟方法,推演断面周围流场演变特征,揭示了流线型箱梁扭转涡振及气动措施抑振机理,为主梁扭转涡振及抑振机理分析提供了一种思路。研究表明:原始断面存在明显扭转涡振现象,其振幅达0.112°,增设检修轨道导流板后振幅降低35.7%,增设抑流板后涡振现象消失。原始断面和导流板断面涡振时,上表面分布气动力矩对涡激力矩贡献值远大于下表面,二者均由上表面大尺度前缘分离涡主导,分离涡漂移时长约为2.5个断面振动周期,对应2阶扭转简化涡模态。增设导流板后,断面上表面分布气动力矩对涡激力矩的贡献显著减小,旋涡漂移模式与断面振动之间的相位关系发生改变,断面周围旋涡作用强度减小,故涡振振幅降低。增设栏杆扶手抑流板后,上表面分布气动力矩对涡激力矩贡献值显著减小且其波浪式分布消失,上表面前缘大尺度分离涡的形成得到抑制,故断面涡振现象消失。
涡振是大跨桥梁在常见风速下易出现的一种自限幅风致振动现象,由来流绕经主梁表面时产生的以某一固定时间间隔有规律脱落的旋涡引起。旋涡脱落、漂移及其非定常演化过程决定了断面表面气动力及其与结构运动之间的相位关系,实现风能向断面动能输入,进而导致桥梁涡振现象发生。涡振效应(振幅)又反过来决定旋涡脱落及漂移模式,如此循环往复。当气动力输能与结构阻尼比耗能达到平衡时,产生涡振极限环振动现
许多学者研究了主梁断面涡振性能,并针对性选用气动措施抑制涡振。XU
扭转涡振现象多见于风洞试验及数值模拟中,在实桥中较少观察到,因此关于断面扭转涡振效应机理的研究较少。SHIRAISHI
针对典型流线型闭口箱梁断面,采用风洞试验和数值模拟相结合的方法,阐述了增设气动措施对断面涡振效应的影响,结合断面表面分布气动力矩贡献和扭转涡振简化涡模态分析,揭示断面扭转涡振机理及气动措施的抑振机理。
研究对象为如

图1 主梁节段模型几何尺寸及测压孔布置(单位:mm)
Fig.1 Geometric sizes of the girder section model and pressure-measuring hole arrangement (Unit: mm)

图2 主梁附属设施布置
Fig.2 Layout of attachments on the girder

图3 气动抑振措施
Fig.3 Aerodynamic countermeasures for suppressing vortex‑induced vibration (VIV)

图4 原始断面及优化断面的扭转涡振响应对比
Fig.4 Comparison between torsional VIV responses of the original girder and the improved girders with aerodynamic countermeasures
综上所述,针对原始断面与优化断面,分别选取各自扭转涡振锁定区内位移的极值点作为典型工况分析该涡振锁定区。因此,后文选取原始断面在折减风速为0.65、导流板断面在折减风速为0.63和抑流板断面在折减风速为0.65下的三种工况进行分析。
分布气动力矩对涡激力矩的贡献同时取决于测点脉动风压大小、测点相对扭转中心的空间位置和分布气动力矩与涡激力矩的相关
(1) |
(2) |
式中,s为测点数量;为减去平均值后的第i个测点风压,测点处风压方向均垂直于断面表面,风压指向断面表面外法线方向为负压,风压指向断面表面内法线方向为正压;为测点每延米所占的权重面积,第测点处权重面积取其相邻两测点间面积的一半;为水平轴方向到壁面外法线方向沿顺时针方向的夹角,取值范围为0~2π;(,为扭转中心坐标;(,)为第i测压点坐标。涡激力矩示意图如

图5 涡激力矩示意图
Fig.5 Schematic diagram of the general vortex excited force (VEF)
箱梁表面分布气动力矩对涡激力矩的无量纲贡献值[
(3) |
式中,为来流风密度;为测点分布气动力矩标准差;为分布气动力矩与涡激力矩的相关系数,表示为:
(4) |
式中,和分别为涡激力矩和分布气动力矩的方差;为涡激力矩和分布气动力矩的协方差。
当为正时,表示第测点区域分布气动力矩对涡激力矩起增强作用;当为负时,表示第测点区域分布气动力矩对涡激力矩起减弱作用。

图6 分布气动力矩对涡激力矩的贡献值
Fig.6 Contribution values of distributed aerodynamic moment to general VEFs
注: X表示测压点的横坐标,即第i测压点坐标xi的集合
综上所述,原始断面与导流板断面上表面分布气动力矩对涡激力矩贡献远大于下表面,表明断面上表面分布气动力矩对涡振起主导作用。相较原始断面,增设导流板后,断面分布气动力矩对涡激力矩的贡献值减小,故扭转涡振振幅减小。增设抑流板后,上表面分布气动力矩对涡激力矩的贡献值急剧减小,其沿下游的波浪式分布消失,故涡振现象消失。因此,断面上表面分布气动力矩对涡激力矩贡献减小是涡振响应降低的关键原因。

图7 流线型闭口箱梁原始断面竖向及扭转涡振贡献值空间分布对比
Fig.7 Comparison of spatial distribution of contribution values of vertical and torsional VIV of streamlined close-box girder
为进一步探究气动措施抑制主梁断面扭转涡振机理,更直观观察增设气动措施前后断面流场变化,选用商用计算流体动力学软件Fluent对原始断面及两种优化断面进行数值模拟分析。计算域设置为长17.5B、宽7B的矩形区域,圆形区域整体为刚体区域,以外围结构化网格作为动网格区域,通过用户自定义函数(UDF)实现圆形区域网格整体的移动与更新。计算域上、下为滑移壁面,计算域右侧为速度入口边界,左侧为压力出口边界,如
(5) |
式中,为扭转涡振的振幅;为断面振动圆频率;为断面振动时刻;为断面扭转位移初相位。

图8 计算域整体及细部网格划分
Fig.8 Computing domain and local view of mesh
湍流模型选择SST k‑ω,空间离散采用2阶迎风格式,时间离散采用2阶隐式积分,以SIMPLE算法进行求解,计算总时间为60 s。原始断面和抑流板断面入口风速均为10.5 m/s,导流板断面入口风速取为10.3 m/s。
Q准则为一种应用较为广泛的旋涡识别方法,HUNT
(6) |
式中,和分别为流体微元和方向的速度矢量的导数项;和分别为流体微元和方向的速度矢量在和方向的交叉导数项。当时,流体微团旋转速率大于应变速率,该区域表现形式为旋涡;当时,流体微团旋转速率小于应变速率,流体单元的运动是非旋转的。
原始断面与导流板断面2.5个断面振动周期内典型时刻的Q值云图如

图9 不同断面振动周期时刻原始断面与导流板断面Q值云图
Fig.9 Cloud images of Q values of original girder and guide vane girder at different moments in the cycle of model vibration

图10 抑流板断面瞬时Q值云图
Fig.10 Cloud image of instant Q value of spoiler girder
对于原始断面下表面,上、下游检修轨道处分别存在一分离泡,上游分离泡尺度基本不变,下游分离泡在下游斜腹板处发展为大尺度尾涡。增设导流板后,上游检修轨道处分离泡尺度减小,下游斜腹板处尾涡尺度减小,增设导流板前后底板均无明显旋涡漂移。因此,可认为原始断面与导流板断面扭转涡振主要由上表面大尺度前缘分离涡主导。
原始断面上表面旋涡形成机理示意图如

图11 原始断面上表面旋涡形成机理示意图
Fig.11 Schematic diagram of formation mechanism of separated vortex on the surface of the original girder
综上所述,原始断面与导流板断面上表面存在大尺度旋涡漂移,且旋涡从产生到脱落经过的时长为2.5个断面振动周期,下表面无明显旋涡漂移,由于上表面区域分布气动力矩对涡激力矩贡献远大于下表面,可知上表面分布气动力矩特性及对应流场特征对扭转涡振起主导作用。增设导流板后上表面旋涡漂移模式与断面振动之间的相位关系改变,同时旋涡尺度减小,旋涡对断面作用强度减小,周围流场对断面的能量输入减少,涡振响应大幅减小。增设抑流板则抑制了上表面大尺度旋涡的产生,扭转涡振消失。
旋涡漂移可用表面压力与断面位移的相位差和涡激力矩与断面位移的相位差之间的差值Δφ表

图12 上表面压力与涡激力矩相位差对比
Fig.12 Comparison of phase difference between the pressure over upper surface and the general VEFs
结合3.2节流场特征可知,旋涡从上表面前缘分离并沿下游周期性漂移,且在漂移的过程中尺度逐渐增大,故涡振时原始断面和导流板断面与周围流场之间复杂的流固耦合现象可用简化涡方

图13 流线型箱梁简化涡示意图
Fig.13 Schematic diagram of simplified vortex model for streamlined closed-box girder
(7) |
式中,为旋涡漂移速度。
简化涡绕流模式中,旋涡漂移速度与断面振动周期之间关系为:
(8) |
式中,(n+0.5)表示旋涡从前缘产生漂移至尾缘脱落经过的断面振动周期数,n为扭转涡振简化涡模态阶数,当n=1时为1阶扭转涡振简化涡模态,当n=2时为2阶扭转涡振简化涡模态,以此类推。旋涡从前缘分离到尾缘脱落经过的时长为:
(9) |
旋涡漂移速度与来流风速之比为:
(10) |
通过简化涡方法得到潜在涡振锁定区的计算折减风速为:
(11) |
原始断面与导流板断面简化涡模态参数如
断面 | 涡振形式 | 简化涡模态阶数n | 旋涡漂移速度与来流风速比 | 计算折减风速 | 折减风速 | 误差/% |
---|---|---|---|---|---|---|
原始断面 | 扭转 | 2阶 | 0.59 | 0.61 | 0.65 | 6.2 |
竖向 | 1阶 | 0.51 | 1.76 | 1.99 | 11.6 | |
导流板断面 | 扭转 | 2阶 | 0.59 | 0.61 | 0.63 | 3.2 |
综上所述,原始断面与导流板断面扭转涡振由上表面前缘产生的分离涡主导,旋涡漂移时长为(2+0.5)个断面振动周期,对应2阶扭转简化涡模态;原始断面第3阶竖向涡振由上表面前缘产生的分离涡主导,旋涡漂移时长为1个断面振动周期,对应1阶竖向简化涡模态。增设抑流板抑制了上表面前缘分离涡的产生,破坏2阶简化涡模态,扭转涡振现象消失。增设导流板不改变上表面旋涡漂移对应的简化涡模态,仅改变旋涡漂移模式与断面振动之间的相位关系与前缘分离涡尺度,扭转涡振效应减弱。
针对典型流线型闭口箱梁断面开展大尺度节段模型测振测压风洞试验,分析典型风速下原始断面与优化断面分布气动力矩贡献,并与竖向涡振进行对比,结合数值模拟方法推演扭转简化涡模态,深入揭示了扭转涡振及检修轨道导流板和抑流板抑振机理。主要结论如下:
(1) 扭转涡振对于断面表面的附属设施十分敏感。增设检修轨道导流板可以降低扭转涡振响应,断面扭转涡振位移减少幅度达35.7%,增设抑流板后断面扭转涡振现象消失。
(2) 断面上表面分布气动力矩对涡振起主导作用。原始断面与导流板断面上表面分布气动力矩对涡激力矩贡献远大于下表面,增设导流板后,上表面分布气动力矩对涡激力矩的贡献值降低,故涡振响应减小。增设栏杆扶手抑流板后上表面分布气动力矩对涡激力矩的贡献值急剧减小,沿下游变化趋向平均分布,涡振现象消失。
(3) 扭转涡振由大尺度的前缘分离涡在断面上表面的周期性漂移主导。前缘分离涡漂移时长约为2.5个断面振动周期,对应2阶简化涡模态。增设导流板改变旋涡漂移模式与断面振动之间的相位关系及旋涡对断面作用强度,周围流场对断面输入能量减少,扭转涡振效应减弱;增设抑流板抑制了上表面前缘分离涡的产生,涡振现象得以消除。
简化涡方法可极大简化涡振时桥梁‑流体之间复杂的流固耦合关系,由贡献值等统计特征推演关键流场特征‑旋涡运动。然而,该方法难以准确表征主梁表面涡激气动力时空分布特征。未来将深入研究涡激气动力时空分布演变特征,提炼关键气动力时空分布模式,建立“气动力特性‑流场特征‑运动状态”全链条分析体系,深入揭示涡振效应物理机制。
参考文献
胡传新, 赵林, 龚玲珑,等. 基于旋涡漂移假设的流线型闭口箱梁竖向涡振机理研究[J]. 振动工程学报, 2024, 37(9): 1575‑1583. [百度学术]
HU Chuanxin, ZHAO Lin, GONG Linglong, et al. Vertical vortex-induced vibration mechanism of streamlined closed-box girder based on vortex drift hypothesis[J]. Journal of Vibration Engineering, 2024, 37(9): 1575‑1583. [百度学术]
OWEN J S, VANN A M, DAVIES J P, et al. The prototype testing of Kessock Bridge: response to vortex shedding[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1996, 60: 91-108. [百度学术]
赵林, 刘丛菊, 葛耀君. 桥梁结构涡激共振的敏感性[J]. 空气动力学学报, 2020, 38(4): 694-704. [百度学术]
ZHAO Lin, LIU Congju, GE Yaojun. Vortex-induced vibration sensitivity of bridge girder structures[J]. Acta Aerodynamica Sinica, 2020, 38(4): 694-704. [百度学术]
XU F Y, YING X Y, LI Y N, et al. Experimental explorations of the torsional vortex-induced vibrations of a bridge deck[J]. Journal of Bridge Engineering, 2016, 21(12): 4016093. [百度学术]
董佳慧, 黄林, 王骑, 等. ∏型钢-混叠合梁斜拉桥涡振性能及整流罩制振措施研究[J]. 振动与冲击, 2022, 41(16): 50-57. [百度学术]
DONG Jiahui, HUANG Lin, WANG Qi, et al. Vortex-induced vibration performance of a cable-stayed bridge with a ∏-shaped steel-concrete composite girder and aerodynamic countermeasure research[J]. Journal of Vibration and Shock, 2022, 41(16): 50-57. [百度学术]
朱思宇, 李永乐, 申俊昕, 等. 大攻角来流作用下扁平钢箱梁涡振性能风洞试验优化研究[J]. 土木工程学报, 2015, 48(2): 79-86. [百度学术]
ZHU Siyu, LI Yongle, SHEN Junxin, et al. Optimization of vortex-induced vibration of flat steel box girders at large attack angle by wind tunnel test[J]. China Civil Engineering Journal, 2015, 48(2): 79-86. [百度学术]
李春光, 张记, 樊永波, 等. 宽幅流线型钢箱梁涡振性能气动优化措施研究[J]. 桥梁建设, 2017, 47(1): 35-40. [百度学术]
LI Chunguang, ZHANG Ji, FAN Yongbo, et al. Study of aerodynamic optimization measures for vortex-induced vibration performance of wide streamlined steel box girder[J]. Bridge Construction, 2017, 47(1): 35-40. [百度学术]
孙延国, 廖海黎, 李明水. 基于节段模型试验的悬索桥涡振抑振措施[J]. 西南交通大学学报, 2012, 47(2): 218-223. [百度学术]
SUN Yanguo, LIAO Haili, LI Mingshui. Mitigation measures of vortex-induced vibration of suspension bridge based on section model test[J]. Journal of Southwest Jiaotong University, 2012, 47(2): 218-223. [百度学术]
宋锦忠, 林志兴, 徐建英. 桥梁抗风气动措施的研究及应用[J]. 同济大学学报(自然科学版), 2002, 30(5): 618-621. [百度学术]
SONG Jinzhong, LIN Zhixing, XU Jianying. Research and appliance of aerodynamic measures about wind-resistance of bridges[J]. Journal of Tongji University (Natural Science), 2002, 30(5): 618-621. [百度学术]
SHIRAISHI N, MATSUMOTO M. On classification of vortex-induced oscillation and its application for bridge structures[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1983, 14(1): 419-430. [百度学术]
刘圣源, 胡传新, 赵林, 等. 中央开槽箱梁断面扭转涡振全过程气动力演化特性[J]. 工程力学, 2020, 37(6): 196-205. [百度学术]
LIU Shengyuan, HU Chuanxin, ZHAO Lin, et al. Aerodynamic force evolution characteristics around the central-slotting box girder during the whole torsional vortex-induced vibration process[J]. Engineering Mechanics, 2020, 37(6): 196-205. [百度学术]
许福友, 林志兴, 李永宁, 等. 气动措施抑制桥梁涡振机理研究[J]. 振动与冲击, 2010, 29(1): 73-76. [百度学术]
XU Fuyou, LIN Zhixing, LI Yongning, et al. Vortex resonance depression mechanism of a bridge deck with aerodynamic measures[J]. Journal of Vibration and Shock, 2010, 29(1): 73-76. [百度学术]
李志国, 周强, 马存明, 等. 中央格栅抑制分离式双箱梁涡振的机理研究[J]. 桥梁建设, 2018, 48(1): 19-24. [百度学术]
LI Zhiguo, ZHOU Qiang, MA Cunming, et al. Mechanism of suppressing vortex-induced vibration of twin-box girder using central grids[J]. Bridge Construction, 2018, 48(1): 19-24. [百度学术]
HU C X, ZHAO L, GE Y J. A simplified vortex model for the mechanism of vortex-induced vibrations in a streamlined closed-box girder[J]. Wind and Structures, 2021, 32(4): 309-319. [百度学术]
HU C X, ZHAO L, GE Y J. Multiple-order vertical vortex-induced vibration mechanism of a typical streamlined closed-box girder[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2022,227:105066. [百度学术]
HU C X, ZHAO L, GE Y J. Mechanism of suppression of vortex-induced vibrations of a streamlined closed-box girder using additional small-scale components[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2019,189: 314-331. [百度学术]
HU C X, ZHAO L, GE Y J. Time-frequency evolutionary characteristics of aerodynamic forces around a streamlined closed-box girder during vortex-induced vibration[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2018,182: 330-343. [百度学术]
HUNT J C R, WRAY A A, MOIN P.Eddies, streams, and convergence zones in turbulent flows[J]. Center for Turbulence Research, 1988: 193-208. [百度学术]