随着桥梁设计理论和建设技术革新,现代桥梁结构逐步向大跨径、高柔度、低阻尼
[1]方向发展,对风荷载变化也越发敏感.涡激振动
[2]是风在结构表面流经引发的限幅振动现象,具有强迫和自激性质, 虽不会对桥梁结构产生破坏,但会极大影响行车安全性与舒适性,并且增大桥梁产生裂纹或疲劳破坏风险
[3].为有效抑制涡激振动,保证桥梁正常使用, 提升桥梁韧性十分有必要.加设气动控制措施可优化桥梁断面气动外形,避免涡激振动,增强大跨径桥梁抵抗力.
π型断面作为典型开口钝体断面,按照主梁纵梁形式分为边箱式与边工字式,通过连接件将桥面和纵梁相连,两者共同受力,具有良好受力性能,同时π型主梁能充分利用混凝土抗压能力,与混凝土桥梁相比可减轻自重增大跨径,加上便于轻型化机械施工吊装,能极好地降低建设成本,被广泛应用于实际工程.但作为典型开口钝体断面,对风荷载作用极为敏感,其抗风稳定性能较差,易发生复杂涡激振动.目前,众多学者针对π型主梁断面涡激振动开展了许多研究.通过风洞试验,Kubo等
[4]研究π型主梁振动幅值和锁定区间与工字梁间距的关系;Lee等
[5]的研究也验证了π型叠合梁边主梁间距对断面涡振有显著贡献作用.通过风洞同步测振测力试验, 李加武等
[6]分析双层π型梁斜拉桥在不同梁间距和风攻角下涡振响应以及静三分力系数变化情况;潘辉
[7]利用Fluent 动网格技术模拟π型主梁竖弯双涡振锁定区间,从气动力、旋涡脱落演化角度来阐述双涡振特性;董锐等
[8]通过对三种π型结合主梁断面进行风洞试验,优化气动选型并总结出少量紊流对涡振周期性脱落有抑制作用;李春光等
[9]以双边主梁钢混叠合梁为工程背景,得出主梁断面上部构造对其涡振性能影响更显著;张国强
[10]对π型桥梁断面形式、栏杆布置、倾角高度等因素进行研究,提高断面整体稳定性;段青松等
[11]针对半开口和分离边箱开口断面分析风攻角、阻尼比等对其主梁涡振性能的影响;李锐
[12]通过选取三种宽高比断面,研究π型叠合梁涡振特性随断面宽高比的变化规律;李欢等
[13]通过调整隔流板和下稳定板等措施来控制施工阶段π型主梁涡振;张天翼等
[14-15]针对宽幅双箱叠合梁进行研究,指出在桥面板两侧加设竖向裙板抑振效果最佳;刘志文等
[16]以π型钢-混凝土结合梁为研究背景,得出与单一气动措施相比,采用组合气动措施优化效果最佳的结论;汪志雄等
[17]通过风洞试验分析中央稳定板、导流板等抑振措施对π型主梁斜拉桥弯扭耦合涡激共振的影响;黄林等
[18]研究导流板倾斜角度对π型叠合梁涡振性能的优化作用.
基于目前π型主梁研究现状
[19-21],导流板因控制措施经济、构造简单,且对π型主梁断面具有良好抑振效果
[22-23],故被广泛应用于实际工程.本文根据部分双边箱式π型断面桥梁的结构参数设计宽高比10∶1的双边箱式π型主梁断面刚体节段模型,进行节段模型同步测振测压试验,分析π型主梁模型涡振规律及表面压力分布变化情况,并且针对导流板措施抑振机理展开参数化研究.
1 节段模型风洞试验
1.1 试验参数
刚体节段模型同步测振测压试验在长安大学风洞试验室CA-01风洞进行.该风洞试验段宽3 m、高2.5 m、长15 m,流场紊流度Iu<0.3%,风洞试验风速在1.0~53 m/s范围连续可调.
图1为π型断面几何参数及测点布置示意图. 在模型1/2长度处布置82个测点,试验采用PSI电子压力扫描阀,采样频率为312.5 Hz,采样时长为45 s.试验时将电子压力扫描阀固定在模型两侧的二元端板上,固定好后计入模型质量并测得相应的阻尼和频率,之后再调整至目标频率.模型骨架使用铝合金以满足刚度要求,外衣采用木制框架,模型两侧安装二元端板保证端部气流两维特性.试验采用二自由度弹簧悬挂体系,见
图2,模型具体动力参数见
表1.
1.2 试验工况
节段模型试验在均匀流场中进行,研究π型断面在0°、±3°风攻角下的涡振特性,并设置多种导流板措施对断面涡振加以控制.导流板措施采用不同尺寸参数的水平导流板与倒L型导流板,其中,倒L型导流板由竖向裙板与水平导流板组成,设水平导流板长度为
b,竖向裙板高度为
d,边箱高度为
D,π型断面宽度为
B,如
图3所示.综合考虑气动优化效果、经济效益以及工程美观等因素,水平导流板选择5种尺寸参数进行试验,倒L型导流板措施参数采用水平悬挑长度 0.5
D、0.75
D、
D,加上两种尺寸参数0.33
D、0.5
D的竖向裙板,具体试验工况见
表2.
2 数值模拟
针对推荐气动措施方案,采用流体力学软件 FLUENT进行数值模拟研究,欧阳克俭等
[24]通过(PIV)粒子图像测速系统验证了在静态流场下桥梁断面数值模拟结果的正确性与合理性,因此采用数值模拟计算原始π型断面与加设气动措施后断面周围的流场特性.本文选取π型断面和加设倒L型导流板(
b=0.75
D,
d=0.5
D)断面作为代表进行数值模拟计算,研究原断面周围流场特性以及倒L型导流板对模型周围流场特性的影响.
依照与风洞试验模型1∶1比例关系进行几何建模,考虑模型在流场中沿顺风向阻塞率应小于5%,同时还原π型断面周围流场,设置计算区域尺寸15
B×10
B(
B为原设计断面模型宽度),采用分块划分网格,π型断面位于网格加密区,模拟风流经断面时周围 流场的变化,在模型顺风向下游设置尾流加密区,便于捕捉尾流涡脱后的变化趋势,而其他部分则设置 较为稀疏网格,减少计算量,提高计算效率,计算域如
图4所示.数值模拟借助FLUENT软件,采用SST
k-
ω湍流模型,设置湍流黏性比10%,湍流强度0.5%,采用SIMPLE算法求解,时间步长统一取0.000 5 s,计算20 000步,共计10 s.
图5为网格划分示意图,网格总数量在9万~10万个,加密区的网格厚度设置为 6.2×10
-5 m,两种断面
y+值均在2以下,近壁面厚度合理,网格质量良好.将CFD数值模拟π型断面在0°风攻角下气动升力系数时程图进行傅里叶变换,得到升力系数的卓越频率.对比风洞试验π型断面0°风攻角两个竖向涡振区间起振风速
v1=8.032 m/s,
v2=18.574 m/s,并通过计算将数值模拟与风洞试验所得的斯托罗哈数
St对比,计算结果所得误差在 10%以内,如
表3所示,表明数值模拟结果可较准确模拟计算π型断面周围气体绕流与旋涡脱落情况.
3 结果与分析
3.1 原始断面
图6为π型断面在0°、±3°风攻角下的涡振响应曲线,其中竖弯无量纲振幅是竖向振幅与主梁高之比(
A/
D),折减风速为风速与结构振动基频和主梁高度乘积之比[
U/(
fbD)].由
图6可知,π型断面在0°、±3°风攻角下均发生竖向与扭转涡激振动且响应幅值较大,同时存在两个竖弯与扭转涡振锁定区间,两者彼此耦合重叠.其中,竖弯涡振高风速锁定区间长度远大于低风速锁定区间,且高风速锁定区间最大竖向位移振幅约为低风速锁定区间的6倍.该π型断面在+3°风攻角下竖弯涡振振幅与锁定区间长度均最大,对π型断面最为不利.对于扭转涡振而言,高风速锁定区间扭转角度均大于低风速锁定区间,此时0°风攻角下扭转涡振角度最大,-3°风攻角下其涡振锁定区间最大.
3.2 气动控制措施
图7为π型断面在0°风攻角下采取5种水平导流板后涡振响应曲线.可以看出,措施1~5对竖弯涡振锁定区间长度、振幅峰值有不同程度的影响.随着水平导流板尺寸增加,高风速锁定区间振幅峰值呈逐次衰减趋势且锁定区间长度减少幅度均在50%以上,但对于低风速锁定区间影响较小且双锁定区间仍然存在,未能彻底抑制π型断面竖弯涡振.相较于竖弯涡振,水平导流板对于扭转涡振抑振效果十分显著,加设最小尺寸水平导流板后,扭转涡振扭转角度峰值迅速降低,随着水平导流板尺寸增加,低风速锁定区间优先得到抑制,高风速锁定区间长度减少,当导流板的尺寸大于边箱高度后,扭转涡振基本消失.因此,试验结果说明π型断面扭转涡振对水平导流板参数变化更加敏感且只加设水平导流板无法有效抑制π型断面涡振,故进一步考虑倒L型导流板措施.
图8为π型断面在0°风攻角下采取6种不同倒L型导流板后涡振响应曲线.可以看出,措施6~10 抑振效果明显,均有效抑制π型断面竖弯和扭转涡振,最大振幅显著减小,而措施11(
b=0.5
D,
d=0.33
D)虽降低锁定区间长度和涡振峰值,但仍发生竖向和扭转涡振,未达到理想抑振效果.
通过对比水平导流板与相同水平长度倒L 型导流板的试验结果,可知水平导流板、竖向裙板尺寸对主梁抑振效果影响显著,竖向裙板对于抑制π型断面涡激振动起决定性作用,竖弯及扭转涡振对竖向裙板长度变化敏感.当倒L型导流板水平长度b≥0.75D时,两种竖向裙板尺寸倒L型导流板均有效抑制π型断面涡振响应.当水平导流板b=0.5D时,两种措施涡振响应有明显区别,其中竖向裙板d=0.5D 时,低风速竖弯涡振响应仍然存在,但高风速涡振最大振幅显著减小.再对比措施9和措施11试验结果,当竖向裙板d=0.33D 时, 水平导流板长度b=0.75D倒L型导流板具有良好的抑振效果,而b=0.5D倒L型导流板无法抑制主梁断面涡振响应,说明水平导流板和竖向裙板尺寸取值均要求在一定范围内,两者组合抑制涡振最佳.
通过上述试验结果得出:倒L型导流板水平长度0.5D~D,竖向裙板长度0.33D~0.5D,在0°风攻角下对π型断面竖弯及扭转涡振均产生良好抑振效果.结合气动优化效果以及实际工程应用,选取设置措施8(b=0.75D, d=0.5D)π型断面进一步分析周围流场特性和表面气动力分布变化情况.
3.3 流场特性
CFD 数值模拟得到原π型断面以及加设导流板措施后断面周围流场瞬时涡量演化图,如
图9、
图10 所示.现选取其中一个旋涡脱落周期最长的涡量演化图,分析π型断面周围气流变化和旋涡在断面附近演变过程及规律,图中蓝色部分为顺时针旋涡,红色部分则为逆时针旋涡.由图可知,来流受到上游前缘边箱阻碍,上下表面气流发生分离流动,模型上表面形成“条状”涡量分布,旋涡在模型中央处发生附着,并且在经过模型前缘后涡量逐渐减小;在原断面下表面生成两个旋涡(W1、W2),旋涡 W1 一直附着在模型下表面,在运动过程中涡量逐步消减,旋涡 W2自
nT+
T/4时刻起,在下表面发生明显旋涡脱落现象,并且在断面后缘尾流区与上表面旋涡发生交替脱落.
加设气动措施后,倒L型导流板使得来流在边箱前缘部分分离程度明显减弱,模型下表面形成的旋涡W1 流速明显降低,与原断面同一位置相比涡量显著减小,旋涡W2 在模型后缘发生涡脱时涡量也减小,且在模型尾流区有效抑制上下表面旋涡的交替脱落.为验证上下表面压力差是激励π型主梁断面发生涡激振动的诱因,列出π型断面及设置导流板绕流场的瞬时压力云图.由
图11可知,在π型断面前缘,边箱迎风侧处于正压状态,上表面一直处在气流分离区内,呈负压状态,下表面旋涡周期性生成并运动,导致下表面该区域压力变化剧烈,进而形成周期性压力差.同时,断面前缘处的生成旋涡随来流运动至断面后缘发生脱落,导致断面后缘区域在一个周期内表现为正压负压交替出现,同样形成了巨大压力差.上下表面前后缘周期性压力差对于涡激振动具有促进作用.
再对比同一时刻(
nT+3
T/4)原断面与倒L型导流板断面的瞬时压力云图和流线图,如
图12所示.由于倒L型导流板对来流分离程度的削弱,旋涡的尺度显著降低,并且上表面前缘负压降幅60%左右,下表面后缘负压降幅40%左右,削弱了断面上下表面的压力差,进而有效抑制了π型断面的涡振.
3.4 表面局部风荷载变化特性
上述3.3节通过CFD数值模拟原断面及措施8(
b=0.75
D,
d=0.5
D)下π型断面流场特性,讨论分析该断面发生涡振的诱因,此节通过测压试验结果进一步分析总结断面周围流场变化后结构表面风压特性和风荷载变化情况.图
13~
15为设置不同导流板π型断面脉动风压系数分布,模型上表面脉动风压系数受模型前缘竖向裙板影响较大,与水平导流板尺寸无关,仅设置水平导流板的模型其表面脉动风压系数分布沿来流方向呈整体下降趋势.加上两种尺寸竖向裙板,倒L型导流板水平长度
b≥0.75
D时,模型上表面脉动风压系数分布沿来流风向保持一致,均有效削弱模型前缘脉动波峰,呈整体均匀分布.倒L型导流板水平长度
b=0.5
D时,两种措施下模型表面脉动风压系数分布有显著区别.
d=0.5
D时,有效地削减了模型表面脉动波峰,与措施6~9风压分布规律整体吻合;而
d=0.33
D时,对模型上表面前缘脉动风压削减程度较小,下表面脉动风压系数分布整体趋势与无竖向裙板时基本保持一致,未能有效降低π型断面表面脉动风压.
由于措施11对π型断面表面脉动风压影响较小,进一步分析该措施下模型表面局部气动力与整体涡激力相关性,见
图16.倒L型导流板水平宽度
b=0.5
D 时,两种措施下模型表面测点相关系数分布不同.
d=0.5
D 时,在倒L型导流板作用下,模型表面测点气动力相关性被破坏,相关系数趋近于零值;
d=0.33
D 时,模型表面相关系数分布与加设水平导流板时整体趋势基本一致,仅在模型下表面前缘下降,对模型中后部无影响.
结合3.3节原断面及加设措施后断面周围的流场特性,通过对比两者表面的流场特性和风压变化,分析可知:π型断面上下表面前缘旋涡的尺度、能量大,负压幅值大且空间分布广,同时,断面上表面前缘处、下表面后缘处的脉动风压幅值也较大,两者具有统一性,说明上表面相对稳定负压区形成较大的周期性压力差与下表面后缘周期性的旋涡脱落,是诱发 π 型断面涡振的主要原因.加设推荐的倒L型导流板措施后,上下表面的旋涡尺度、涡量明显降低,压力差减少,并且上下表面脉动压力峰值也显著下降,降幅达 200%~300%,其他区域的压力分布也呈下降趋势,局部气动力与整体涡激力的相关性被破坏,削弱了气动力对整体涡激力的贡献程度,抑振效果极佳.
通过对比数值模拟和测压试验的结果,进一步证明合适尺寸的倒 L 型导流板能有效削弱气流分离程度,抑制旋涡生成,降低模型上下表面压力差,从而抑制桥梁涡振,降低桥梁疲劳破坏风险.然而大尺寸倒L型导流板需要更大板宽及质量,对桥梁整体稳定性存在影响,需要后续进一步研究验证.
4 结 论
本文以边箱式π型断面为研究对象,通过刚体节段模型测振测压实验,研究π型断面涡振特性和不同参数导流板抑振性能,通过数值模拟分析断面流场特性,得到以下结论:
1)0º、±3º风攻角下,宽高比10∶1的π型断面存在显著涡激振动,有两个涡振锁定区间且高风速涡振锁定区间振幅较大,断面抗风性能较差,只加设水平导流板措施无法有效抑制涡振,而倒L型导流板措施能显著抑制涡振响应.
2)对于倒L型导流板措施,水平导流板与竖向裙板尺寸参数对抑制π型断面涡振有重要影响,尺寸参数过大或过小,其优化效果均不佳.倒L型导流板尺寸参数(b=0.75D,d=0.5D)抑振效果最好.
3)倒L型导流板措施可改善风流经断面前缘分离程度,优化气动外形,削弱旋涡涡脱能量,削减π型断面上下表面周期压力差,使得流场与风荷载分布合理.
国家重点研发计划资助项目(2021YFB2600600); 长安大学中央高校基本科研业务费专项资金资助项目(300102213101)