力学学报  2018 , 50 (4): 766-775 https://doi.org/10.6052/0459-1879-18-057

Orginal Article

低雷诺数下串列三圆柱涡激振动中的弛振现象及其影响因素 <sup>^1)</sup>

陈威霖, 及春宁*, 许栋

天津大学水利工程仿真与安全国家重点实验室, 天津300072

GALLOPING IN VORTEX-INDUCED VIBRATION OF THREE TANDEM CYLINDERS AT LOW REYNOLDS NUMBERS AND ITS INFLUENCING FACTORS

Chen Weilin, Ji Chunning*, Xu Dong

State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China

中图分类号:  P751,TV131.2

文献标识码:  A

通讯作者:  *及春宁, 教授, 主要研究方向为流固耦合, 涡激振动, 泥沙运动. E-mail:cnji@tju.edu.cn*及春宁, 教授, 主要研究方向为流固耦合, 涡激振动, 泥沙运动. E-mail:cnji@tju.edu.cn

版权声明:  2018 力学学报期刊社 力学学报期刊社 所有

基金资助:  国家自然科学基金(51579175, 51779172), NSFC-广东联合基金 (第二期) 超级计算科学应用研究专项资助(U1501501)和国家超级计算广州中心支持项目.

展开

摘要

对间距比为1.2和雷诺数为100的串列三圆柱涡激振动进行数值模拟, 发现在某个折合流速之后, 三圆柱的响应均呈现为随着折合流速增大而增大的弛振现象, 平衡位置偏移、低频振动以及旋涡脱落与圆柱运动之间的时机三个因素共同决定了弛振现象的出现. 进一步的研究发现, 串列三圆柱的弛振现象仅出现在质量比不大于2.0和雷诺数不大于100的工况下. 当质量比较大时, 串列三圆柱的平衡位置固定不变, 且圆柱的振动不规律, 使得旋涡脱落与圆柱运动的时机处于变化之中. 当雷诺数较高时, 最上游圆柱的平衡位置在折合流速较大时回到初始位置, 不再参与对圆柱振动的调节, 使得圆柱的振动响应不再规律, 旋涡脱落与圆柱运动的时机也一直处于变化之中.

关键词: 涡激振动 ; 串列三圆柱 ; 弛振 ; 质量比 ; 雷诺数

Abstract

Vortex-induced vibrations of three tandem cylinders at the spacing ratio of 1.2 and the Reynolds number of 100 shows that when the reduced velocity is larger than a critical value, the galloping is observed where the amplitudes of three cylinders increase with the reduced velocity. Three factors, including shift of equilibrium position, low vibration frequency and timing between vortex shedding and motion of cylinders, determine the appearance of galloping. Further investigations show that the galloping occurs at a range with lower mass ratio no more than 2.0 and Reynolds number no more than 100. When the mass ratio is large, the equilibrium position is unchanged and the displacements are irregular, which leads to the variation of timing between vortex shedding and motion of cylinders. When Re is large, the shift of the most upstream cylinder is zero and no more accommodate cylinders’ vibration, which makes the vibration irregular and change of timing between vortex shedding and motion of cylinders.

Keywords: vortex-induced vibration ; three tandem cylinders ; galloping ; mass ratio ; Reynolds number

0

PDF (10472KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

陈威霖, 及春宁, 许栋. 低雷诺数下串列三圆柱涡激振动中的弛振现象及其影响因素 <sup>^1)</sup>[J]. 力学学报, 2018, 50(4): 766-775 https://doi.org/10.6052/0459-1879-18-057

Chen Weilin, Ji Chunning, Xu Dong. GALLOPING IN VORTEX-INDUCED VIBRATION OF THREE TANDEM CYLINDERS AT LOW REYNOLDS NUMBERS AND ITS INFLUENCING FACTORS[J]. Acta Mechanica Sinica, 2018, 50(4): 766-775 https://doi.org/10.6052/0459-1879-18-057

弛振现象一直是工程实际中一个值得关注的现象. 例如, 桥梁工程中的悬索桥, 由于其自身质量小、频率低、阻尼小, 在风的作用下极易发生大幅振动, 塔科马悬索桥的破坏就是一个典型的例子[1]. 此外, 工业中常见的热交换器管阵也容易发生弛振, 造成热交换器的疲劳破坏. 由于热交换管的直径较小, 且热媒多为密度较小的水蒸气, 其雷诺数较低. 因此, 对低雷诺数弛振现象的研究将具有非常重要的意义.

弛振现象最早发现于非圆形结构物的涡激振动中, 随着来流流速的增大, 结构的振幅持续增大. Den Hartog[2]基于对被冰包裹电缆的研究发现弛振一般伴随着低频的出现. Parkinson和Smith[3] 建立了一个准稳态的理论, 用来预测弛振下方柱的振幅. 该理论成功地预测了较大和较小阻尼下空气中方柱的弛振现象. Nemes 等[4]通过水槽实验研究发现, 当攻角为 α=0时, 振动的方柱表现为弛振, 角度设置如图1(a)所示. Zhao 等[5]发现, 三倍振动频率的升力成分对方柱弛振的出现有重要的影响. 目前的研究表明, 方柱的弛振现象仅在较高的雷诺数下出现. Barrero-Gril 等[6] 发现当Re < 159时, 振动的方柱不会出现弛振. Zhao 等[7] 对不同攻角( 0~45)下弹性支撑方柱的数值模拟发现, 当 Re=100时, 弛振没有出现. 此外, 对弹性支撑三角柱实验研究发现, 当雷诺数较高时, 振动的三角柱在某些攻角下表现为弛振现象[8,9]. Wang 等[10]对不同攻角( 0~60) 下振动三角柱的数值模拟研究表明, 当 α=60Re=100时, 攻角设置如图1(b)所示, 三角柱表现为弛振现象, 振幅随着折合流速的增大而增强.

图1   方柱和三角柱的攻角(α).

Fig.1   The attack angle of the square and triangular cylinder

由于圆柱本身绝对对称, 因此, 光滑单圆柱的涡激振动并未表现出弛振. Chang 等[11] 通过在光滑圆柱表面贴粗糙条的方式使其振动呈现出随流速增大而增强的弛振现象. 此外, 增加圆柱的个数和布置形式是另一种激发弛振现象的方式. Brika和Laneville[12]Re5×103~2.7×104 范围内串列双圆柱的研究发现, 当上游圆柱固定时, 下游圆柱会在一定的间距比范围内表现为弛振现象. Bokaian 和Geoola[13]Re7×102~2×103 范围内串列双圆柱的研究发现, 当下游圆柱静止不动并浸没于上游圆柱尾流中时, 振动的上游圆柱也可以表现为弛振. Assi 等[14,15,16]的实验研究表明, 下游圆柱的尾流弛振现象可以在L*处于2.0~5.6 之间时出现, 且顺流向振动与否对尾流弛振的出现没有明显的影响.

虽然振动的串列双圆柱可以在较高的雷诺数下表现出尾流弛振现象, 但目前的研究表明, 尾流弛振不会出现较低雷诺数下串列双圆柱的涡激振动中[17,18,19,20]. 此外, 对串列三圆柱涡激振动的研究也表明当间距比较大时, 尾流弛振现象也不会出现[21]. Werle[22]运用流动显示方法考察了间距比 L*为3.3下串列三圆柱的流场, 其中雷诺数为 Re=2×103. Igarashi和Suzuki[23] 通过风洞试验研究了间距比 L*1.0~4.0之间串列三圆柱绕流情况, 其中雷诺数为 1.09×104<Re<3.92×104. 此外, 采用烟流动显示方法考察了流动情况, 还测量了各圆柱的压力分布, 从而得到了各圆柱的阻力系数. Chen 等[24]L*处于1.2 ~5.0之间和 Re=100下振动串列三圆柱的数值模拟研究发现, 当 L*=1.2时, 串列三圆柱的振动均表现为弛振现象, 也就是说, 三圆柱的振幅均随着折合流速的增大而增大, 并一直持续到较大的折合流速. Chen 等[24]指出决定弛振现象的三个关键因素为平衡位置偏移、低频振动以及旋涡脱落与圆柱运动之间的时机. 平衡位置的偏移可以很好地调和圆柱之间的相互作用, 上游圆柱脱落的旋涡对下游圆柱的振动可以实现最大的促进作用. 由于三圆柱振动频率相同且均为低频, 也就是说旋涡从圆柱上脱落需要更多的时间, 因此, 可以对圆柱振动有更久的促进作用, 保证了圆柱具有较高的振幅. 由于串列三圆柱之间相位差的改变, 旋涡脱落与圆柱的振动达到了较好的时机, 当上游圆柱的旋涡运动到下游圆柱时, 恰好可以在下游圆柱的运动方向上产生低压区, 对圆柱的振动起到促进作用. 可以这么说, 三因素一起对串列三圆柱弛振的出现起到决定性作用, 因此, 三圆柱的弛振源自于圆柱之间强烈的稳定周期性相互作用, 这与高雷诺数下串列双圆柱的尾流弛振现象的机制明显不同[15]. 需要说明的是, 虽然Chen 等[24] 分析了尾流弛振现象的产生机理, 但是并未讨论一些关键因素对弛振现象的影响. 基于对该尾流弛振深入理解的需要, 本文重点研究质量比和雷诺数两个关键因素的影响, 找出该现象可能出现的范围.

1 数值方法

1.1 控制方程

流固耦合的数值模拟采用浸入边界法[25], 控制方程如下

ut=-(u)u-p+v2u+f(1)

u=0(2)

其中, u为速度, t为时间, p为压强, v为运动黏滞系数, 为梯度算子, f为附加体积力矢量, 代表流固耦合边界条件.

对以上控制方程采用二阶精度的Admas-Bashforth时间格式进行离散, 可得守恒形式如下

un+1=un+t32hn-12hn-1-32pn+12pn-1+

fn+12t

un+1=0

其中, h=-uu+vu+uT由对流项与扩散项组成, 上标 T为矩阵转置, 附加体积力表示为 fn+12t=D{Vn+1-I[un+t(32hn-12hn-1-32pn+12pn-1)]}, ID为插值函数, Vn+1为物面边界的速度, 上标 n+1, n+1/2, n, n-1为时间步.

针对传统浸入边界法施加边界条件精度不高的情况, Ji等[25]提出了基于嵌入式迭代的浸入边界法, 将浸入边界法嵌入到压强泊松方程的迭代求解中, 利用压强的中间解比初始值更接近真实值的特点, 迭代修正附加体积力, 在不显著增加计算耗时的前提下, 提高了整个算法的求解精度. 有关浸入边界法的细节, 请参考文献[25], 此处不再赘述.

对仅做横流向运动的刚性圆柱, 其运动方程可以用下述方程来描述

md2ydt2+cdydt+ky=Fy

其中, m为圆柱质量, c为结构阻尼, k为弹簧刚度系数, Fy为圆柱受到的横流向流体力. 方程采用标准的Newmark- β法求解.

1.2 问题描述

串列三圆柱涡激振动数值模拟的相关参数设置如下. 雷诺数为 Re=UD/v, 在 80~180之间取值, 质量比为 m*=m/mf, 在 2.0~10.0中取值, 其中 m为圆柱质量, mf为等体积流体质量. 相邻两圆柱的间距相等, L为圆柱中心的距离, D为圆柱直径, 取为 L*=L/D=1.2. 为使圆柱响应达到最大, 将阻尼比设为0. 此外, 为避免碰撞, 三圆柱均仅作横向振动. 计算域大小为 200D×50D, 如图2所示, 中间圆柱位于计算域的中心. 从上游到下游, 三圆柱依次编号为1, 2 和3. 为保证数值精度[26], 在圆柱周围采用加密网格, 加密区域大小为 12D×8D, 对应的网格尺度为 Δx/D=Δy/D=1/64.

图2   计算域和边界条件设置.

Fig.2   Computational domain and the boundary conditions

边界条件设置如下. 入口为Dirichlet型边界(u=U,v=0), 出口为Neumann边界(u/x=0,v/x=0), 上下为可滑移边界(u/y=0,v=0). 此外, 为保证数值收敛, 需满足CFL条件, 即UmaxΔt/Δx#x2264;0.5, 选定的时间步长为Δt=0.005.

2程序验证

表1给出了串列三圆柱绕流下各圆柱阻力均值与Harimi和Saghafian[30]结果的对比, 其中最大差别仅为6.7%, 从而验证本文数值方法和程序的正确性. 更多的验证算例包括单圆柱涡激振动、并列双圆柱绕流、串列双圆柱绕流和涡激振动等参见文献 [24,26-29].

表1   串列三圆柱绕流的阻力均值, 其中L*=2.0Re=100

Table 1   The mean drag coefficients of flow past three tandem circular cylinders at L*=2.0 and Re=100

C̅d1C̅d2C̅d3
Ref.[30]1.152-0.1070.127
present1.162-0.1140.119
difference/%0.96.16.7

新窗口打开

3 结果和讨论

3.1 弛振现象

图3所示, 当 L*=1.2, Re=100m*=2.0时, 串列三圆柱的最大振幅( Am*=Am/D, Am为最大振幅)均随折合流速 Ur的增大而增大, 在达到峰值之后会随着折合流速的增大而稍微减小. 但在某个折合流速后, 三圆柱的振幅又会随着折合流速的增大而增大, 并持续到非常大的折合流速, 因此称之为弛振现象. 需要说明的是圆柱1, 2 和3 的弛振分别开始于 Ur=17, Ur=13Ur=35. 此外, 随着折合流速的增大, 三圆柱振幅增大的幅度也并不相同, 如图3 所示, 圆柱2 振幅的增大幅度要稍大于其余两个圆柱的情况. 但需要指出的是, 与高雷诺数下串列双圆柱尾流弛振现象中振幅的增大幅度[14,15,16]相比, 层流中串列三圆柱的振幅增幅则要小了很多. 此外, 层流中串列三圆柱均呈现出弛振现象, 而在高雷诺数下串列双圆柱中, 仅下游或者仅上游圆柱呈现弛振现象[12,13]. 由图3 可知, 随着折合流速的增大, 三圆柱的振幅均可以达到较大的值, 其中圆柱2 和3的振幅接近 1.3D, 与高雷诺数下单圆柱涡激振动的最大振幅相当[30], 而圆柱1的振幅也达到了 0.8D, 比相同雷诺数下单圆柱涡激振动的最大值(0.57 D) 要大40%.

图3   当L*=1.2, Re=100m*=2.0时, 串列三圆柱最大振幅 (Am*)随折合流速变化情况.

Fig.3   The maximum amplitudes (Am*) of three cylinders versus Ur at L*=1.2, Re=100 and m*=2.0

Chen 等[24]对该弛振现象的产生给出详细的解释, 指出决定弛振现象的三个关键因素: 平衡位置偏移、低频振动以及旋涡脱落与圆柱运动之间的时机. 具体来看, 如图4 所示, 当上游圆柱向上运动时, 由于受到下游圆柱的推挤作用, 上游圆柱未分离的剪切层会堆积在其上侧(A~C), 因此会产生相应的低压区, 对其运动起到促进作用. 当下游两圆柱向下运动时(B~D), 由于此时中间圆柱的运动相位要领先于最下游圆柱, 使得上游圆柱分离的剪切层会依次重附着于下游两圆柱的下侧, 且出现的时刻恰好位于圆柱正加速向下运动的过程中[24], 因此对下游圆柱的振动起到了相应的激励作用. 由于中间圆柱本身也要分离一定强度的剪切层, 因此重附着于下游圆柱上的剪切层强度增大, 与下游圆柱更大的振幅相一致. 当下游两圆柱开始向上运动时(E~H), 上游圆柱分离的剪切层又恰好依次出现在下游两圆柱的上侧, 并对它们的运动再次产生促进作用. 在这个过程中, 串列三圆柱的平衡位置偏移、旋涡脱落与圆柱运动的时机都成为极其重要的因素. 由于圆柱的振动频率较低, 因此重附着的剪切层可以提供更长时间的促进作用. 由于Chen 等[24]并未提及质量比和雷诺数等对该现象的影响, 而这两个因素在涡激振动中又是极其重要的, 因此, 接下来本文将研究两者对弛振的影响, 进一步研究弛振现象可能出现的范围.

图4   一个周期内串列三圆柱涡激振动周围的涡量场情况, 其中m*=2.0, L*=1.2, Ur=40.0.

Fig.4   Vorticity fields around three tandem circular cylinders in one period, where m*=2.0, L*=1.2, Ur=40.0

3.2 质量比对弛振的影响

选取m*=5.0m*=10.0研究质量比对串列三圆柱流致振动中迟振现象的影响. 其中雷诺数为Re=100和间距比为L*=1.2, 保持不变. 此外, 采用的两个质量比下串列三圆柱的振动不再规律, 因此, 振幅的最大值不再适用于表示圆柱的振动响应, 本节采用无量纲的均方根值(Ar*=Ar/D, Ar为均方根振幅)来表示圆柱的振动响应.

图5所示, 当 m*=5.0时, 圆柱1和2的振幅随着折合流速的增大而先增后减, 在某个折合流速之后, 振幅随着折合流速的增大而持续增大, 这与质量比 m*=2.0时的情况相同. 然而圆柱3的振幅却没有再出现随折合流速增大而增大的现象, 仅在 Ur=35时出现稍微的增大. 也就是说, 当质量比 m*=5.0时, 仅圆柱1和2呈现出 “弛振” 现象. 当 m*=10.0时, 圆柱1 和3 的振幅呈现出随折合流速先增后减之后持续增大的现象, 圆柱2 的振幅则一直随着折合流速的增大而增大. 对比可以发现, 圆柱2 振幅的增大幅度要明显大于其余两个圆柱, 与质量比 m*=2.0时的情况相同. 因此, 在质量比 m*=10.0下, 串列三圆柱的振动均表现为明显的 “弛振” 现象. 需要指出的是, 实际上不论是质量比为5.0还是10.0下串列三圆柱涡激振动中的弛振现象仅出现在一定的折合流速范围内, 严格意义上来讲, 并不是弛振现象. 如表2所示, 当折合流速增加到80时, 质量比为5.0和10.0下三圆柱的振幅均出现明显下降, 几乎不再振动, 也就是说, 三圆柱振幅随折合流速增大而增大的弛振现象仅存在于一定的折合流速范围内. 不同的是, m*=2.0下串列三圆柱的振幅则一直随着折合流速的增大而增大. 因此, 可以这么认为弛振现象仅存在于 m*#x2264;2.0下串列三圆柱涡激振动中.

图5   不同质量比下串列三圆柱响应随折合流速变化情况.

Fig.5   Variations of the vibration of three cylinders with Ur at different mass ratio

表2   不同质量比和折合流速下串列三圆柱的振幅均方根

Table 2   The r.m.s. amplitudes of three cylinders at different mass ratios and reduced velocities

m*UrAr1*Ar2*Ar3*
*5.0401.051.150.81
800.020.030.06
*10.0400.721.271.00
800.000.010.02

新窗口打开

为进一步解释该现象不再存在于质量比为5.0和10.0下串列三圆柱涡激振动中的原因, 对该现象产生的决定性因素展开讨论. 首先, 当质量比较大时, 串列三圆柱的响应不再规律, 如图6所示( m*=5.0Ur=40). 也就是说三圆柱之间的相互作用是不稳定和非周期性的, 旋涡与圆柱运动之间的时机也是随时间发生变化的. 其次, 如图7所示,当 m*=2.0时, 串列三圆柱的平衡位置发生偏移, 不再固定在 y̅=0上; 相反, m*=5.0和10.0 下串列三圆柱的平衡位置一直在 y̅=0上, 并未发生变化, 因此也就无法通过调整圆柱之间的相对位置来调节圆柱之间的相互作用.

图6   当m*=5.0Ur=40.0时, 串列三圆柱振动的历时曲线.

Fig.6   Time histories of the displacements of three cylinders with m*=5.0 and Ur=40.0

图7   不同质量比下三圆柱平衡位置偏移随折合流速变化的情况.

Fig.7   Variations of the shift of the equilibrium position with Ur at different Reynolds numbers

串列三圆柱涡激振动中, 三圆柱的斯特劳哈尔数St总是相等[24], 因此, 本文仅给出圆柱1 的St.图8所示, 各质量比下, St在某个折合流速以后与自然频率相当, 圆柱响应进入锁定区间. 此后, 随折合流速的增大而逐渐减小. 当m*=10.0时, St 与自然频率最接近且频率最小. 因此, 可以这么说, 低频振动并不是串列三圆柱尾流弛振的唯一决定因素, 需要三个决定因素的共同作用. 此外, 还可以发现当St 以后, 单圆柱涡激振动的St要明显大于串列三圆柱的情况, 这是因为在较小间距比下的串列三圆涡激振动中, 上游圆柱分离的剪切层会附着于中间圆柱上并最终从下游圆柱脱落. 下游两圆柱的存在不仅使得上游圆柱剪切层的形成过程变慢, 而且脱落过程也比单圆柱涡激振动需要更长的时间. 因此, 串列三圆柱中St要明显小于单圆柱涡激振动的情况. 当St[24], 上游圆柱剪切层可以自由脱落形成旋涡, 此时串列三圆柱的St呈现出与单圆柱涡激振动相同的变化规律, 且Ur也与单圆柱涡激振动相近.

图8   不同雷诺数下串列三圆柱的m*随折合流速变化情况.

Fig.8   Variations of the Strouhal number of three cylinders with y̅=0 at different Reynolds numbers

综上, 虽然Re=80取5.0和10.0时串列三圆柱的振动在一定的折合流速范围内较大振幅, 但是由于此时圆柱的质量比较大, 惯性较大, 这使得上游圆柱分离的剪切层与下游圆柱运动相互作用对振动的调节降低, 串列三圆柱的平衡位置一直固定在初始位置(Re=180)上, 并未发生变. 相应地, 下游圆柱的振动不再规律, 旋涡脱落与圆柱运动的时机是一直变化, 不再是稳定的相互作用. 虽然质量比较大时, 串列三圆柱的振动频率较低, 由于脱涡与圆柱运动之间时机不稳定, 串列三圆柱涡激振动中的弛振现象并未出现在质量比较大的工况下.

3.3 雷诺数对弛振的影响

选取 m*=2.0L*=1.2研究雷诺数对串列三圆柱涡激振动中弛振现象的影响, 此时质量比为 Am*和间距比为 Re=80. 由于研究的雷诺数范围内, 串列三圆柱的振动规律, 因此用圆柱振动的最大幅值( Ur>4.0) 表示圆柱的振动响应.

图9(a)所示, 当 Re=100时, 串列三圆柱的振动在 Ur=40以后急剧增强并在达到峰值之后缓慢减弱, 之后三圆柱的振动在某个折合流速以后随着折合流速的增大而持续增强. 因此, 可以认为此时串列三圆柱涡激振动中出现了弛振现象. 与 Ur=80时的工况相比, 此时圆柱的振幅要更大一些, 随着折合流速增大而增大的幅度也要稍大一些. 表3 给出了折合流速 Re=80Ur=80时串列三圆柱最大振幅的对比结果, 可以发现, 雷诺数 Re=80下串列三圆柱的振幅增大可以持续到 Re=100时; 雷诺数 Ur时串列三圆柱的振幅要明显大于雷诺数 Re下的情况.

图9   不同雷诺数下串列三圆柱响应随折合流速变化的情况.

Fig.9   Variations of the vibration of three cylinders with m* at different Ur

表3   不同雷诺数和折合流速下串列三圆柱的最大振幅

Table 3   The maximum amplitudes of three cylinders at different Reynolds numbers and reduced velocities

Am1*Am2*Am3*Re=180Ur=9.0
*80400.871.331.35
800.981.481.50
*100400.801.251.27
800.841.281.28

新窗口打开

Ur=22.0时, 串列三圆柱的振幅从 Re以后开始增大, 均在 Re=180时达到最大值, 此后, 振幅随折合流速增大而减小. 如图9(b) 所示, 三圆柱的振幅均在 Re=80以后又开始增大并在 Re=100出现一个峰值. 当 Re=180以后三圆柱的振幅均几乎不再发生变化. 虽然三圆柱的振幅可以维持在较大的值上, 但是由于振幅不再随折合流速的增大而增大, 因此, 弛振现象不再出现.

为进一步解释尾流弛振现象出现在低雷诺数 ( Ur<30.0为80和100) 而没有出现在 Ur>30.0下串列三圆柱涡激振动中的原因, 对该现象出现决定性因素展开如下讨论.

图10所示, 各雷诺数下, 串列三圆柱的平衡位置均发生偏移. 随着雷诺数的增大, 平衡位置偏移出现更小的折合流速下. 对比 y̅=0Ur=40.0下平衡位置的偏移可以发现, 圆柱2和3 的平衡位置偏向同一个方向, 而圆柱1会偏向另一个方向. 当 Ur时, 圆柱1和2在 St时偏向同一个方向, 而圆柱3偏向另一个方向. 当 Re=180以后, 圆柱1 的平衡位置逐渐恢复到初始位置 ( Ur=4.0) 上. 当 Re时, 圆柱1 的平衡位置不再发生偏移, 而圆柱2 和3的平衡位置偏向了同一个方向, 且此时圆柱3 的平衡位置偏移要明显大于圆柱2的情况. 因此, 可以这么认为, 此时圆柱1不再具有调和圆柱2 和3 振动响应的作用了.

图10   不同雷诺数下三圆柱平衡位置偏移随折合流速变化的情况.

Fig.10   Variations of the shift of the equilibrium position with Ur=5.0 at different Reynolds numbers

图11所示, 各雷诺数下, 串列三圆柱的St非常接近. 当St时, 串列三圆柱在St时进入锁定区间, 而当Ur取80和100时, 三圆柱在Re=180时进入锁定区间, 之后, 随着折合流速的增大, Ur=80.0下降. 在锁定区间以后的折合流速范围内, 不同雷诺数下的Ur=80.0均非常接近.

图11   不同雷诺数下串列三圆柱的Re=180随折合流速变化情况.

Fig.11   Variations of the Strouhal number of three cylinders with Ur=80.0 at different Reynolds numbers

图12所示, 当Re=180Ur=80.0时, 串列三圆柱的振动不规律, 此时串列三圆柱之间的相互作用不稳定. 图13给出了Re=180时串列三圆柱后的尾流情况, 尾涡不规律, 旋涡脱落与圆柱运动的时机不固定, 弛振现象不再出现.

图12   当Ur=80.0Re=180时, 串列三圆柱响应的历时曲线.

Fig.12   Time histories of the displacements of three cylinders at Ur=80.0 and Re=100

图13   当m*=2.0L*=1.2时, 串列三圆柱后的涡量场.

Fig.13   Vorticity contour around three cylinders at m* and m*=2.0

4 结论

本文对串列三圆柱涡激振动中出现的弛振现象进行了深入研究, 分析了质量比和雷诺数对弛振现象的影响. 研究发现, 当Re#x2264;100, Re=80Re=100时串列三圆柱涡激振动中出现了弛振现象, 也就是说, 当折合流速大于某个值以后, 圆柱的振幅随着折合流速的增大而增大并持续到较大的折合流速. 决定串列三圆柱涡激振动中出现弛振现象的因素包括平衡位置的偏移、较低的振动频率和旋涡与圆柱振动的时机.

质量比对串列三圆柱中弛振现象的出现有显著影响. 当质量比较大 (5.0和10.0) 时, 圆柱的振动不再规律, 圆柱之间的相互作用也不再稳定, 使得旋涡与圆柱振动之间的时机处在不断的变化之中. 此外, 当质量比较大 (5.0和10.0) 时, 串列三圆柱的平衡位置不再发生偏移, 因此, 无法实现对圆柱振动的调节, 从而使圆柱之间的相互作用达到最优. 虽然在质量比 Re=180取5.0和10.0时串列三圆柱的振动频率仍较低, 甚至要低于 Re=180时的情况, 但由于其余两个决定性的因素无法满足, 使得弛振现象未出现在大质量比的串列三圆柱涡激振动中.

雷诺数同样是弛振的出现与否一个关键因素. 研究发现, 弛振现象仅出现在雷诺数较低( Re=80)的工况下. 当 y̅=0时, 串列三圆柱的平衡位置偏移出现在更大的折合流速下, 且平衡位置的偏移也要小于 3.0×104#x2264;Re#x2264;1.2×105时的情况. 当{Invalid MML}时, 串列三圆柱的振幅在{Invalid MML}以后便不再变化. 虽然当{Invalid MML}时, 串列三圆柱的平衡位置发生偏移, 且要早于{Invalid MML}和100的情况. 但当{Invalid MML}以后, 圆柱1 的平衡位置回到初始位置 ({Invalid MML}) , 因此, 不再参与对圆柱2 和3 振动的调节. 此外, 圆柱的振动和圆柱后的尾流变得不再规律, 旋涡脱落与圆柱之间的时机也就处于变化之中, 圆柱之间的耦合作用不再稳定.

综上可以判断, 串列三圆柱涡激振动中的弛振现象仅出现在质量比较小和雷诺数较小的参数空间内, 这与之前研究成果明显不同. 弛振现象的出现对结构物有非常显著的破坏作用, 这是因为此时结构物的振动会随着来流流速的增大而增大, 破坏作用也成倍增大. 之前研究结果表明弛振现象仅出现在非对称体, 比如某些攻角下的三角柱和方柱, 以及高雷诺数下串列双圆柱涡激振动中. 本文的研究结果表明, 弛振现象可以出现较低雷诺数下的低密度串列三圆柱涡激振动中. 此结论对实际工程结构物设计与防护具有明显的参考意义.

The authors have declared that no competing interests exist.


参考文献

[1] Steinman DB, Watson SR.

Bridges and Their Builders. Dover

, New York. 1957

[本文引用: 1]     

[2] Den Hartog JP.

Transmission line vibration due to sleet

.Transactions of the American Institute of Electrical Engineers, 1932, 51(4): 1074-1076

[本文引用: 1]     

[3] Parkinson GV, Smith JD.

The square prism as an aeroelastic nonlinear oscillator

.The Quarterly Journal of Mechanics and Applied Mathematics, 1964, 17(2): 225-239

[本文引用: 1]     

[4] Nemes A, Zhao J, Lo Jacono D, et al.

The interaction between flow-induced vibration mechanisms of a square cylinder with varying angles of attack

.Journal of Fluid Mechanics, 2012, 710: 102-130

[本文引用: 1]     

[5] Zhao J, Leontini JS, Lo Jacono D, et al.

Fluid-structure interaction of a square cylinder at different angles of attack

.Journal of Fluid Mechanics, 2014, 747: 688-721

[本文引用: 1]     

[6] Barrero-Gil A, Sanz-Andre´s A, Roura M.

Transverse galloping at low Reynolds numbers

.Journal of Fluids and Structures, 2009, 25: 1236-1242

[本文引用: 1]     

[7] Zhao M, Cheng L, Zhou T.

Numerical simulation of vortex-induced vibration of a square cylinder at a low Reynolds number

.Physics of Fluids, 2013, 25: 023603

[本文引用: 1]     

[8] Alonso G, Meseguer J.

A parametric study of the galloping stability of two-dimensional triangular cross-section bodies

.Journal of Wind Engineering and Industrial Aerodynamics, 2006, 94: 241-253

[本文引用: 1]     

[9] Alonso G, Sanz-Lobera A, Meseguer J.

Hysteresis phenomena in transverse galloping of triangular cross-section bodies

.Journal of Fluids and Structures, 2012, 33: 243-251

[本文引用: 1]     

[10] Wang H, Zhao D, Yang W, et al.

Numerical investigation on fl ow-induced vibration of a triangular cylinder at a low Reynolds number

.Fluid Dynamics Research, 2015, 47: 015501

[本文引用: 1]     

[11] Chang CC, Kumar RA, Bernitsas MM.

VIV and galloping of single circular cylinder with surface roughness at {Invalid MML}

.Ocean Engineering, 2011, 38(16): 1713-1732

[本文引用: 1]     

[12] Brika D, Laneville A.

The flow interaction between a stationary cylinder and a downstream flexible cylinder

.Journal of Fluids and Structures, 1999, 13: 579-606

[本文引用: 2]     

[13] Bokaian A, Geoola F.

Wake-induced galloping of two interfering circular cylinders

.Journal of Fluid Mechanics, 1984, 146: 383-415

[本文引用: 2]     

[14] Assi GRS, Meneghini JR, Aranha JAP, et al.

Experimental investigation of flow-induced vibration interference between two circular cylinders

.Journal of Fluids and Structures, 2006, 22: 819-827

[本文引用: 2]     

[15] Assi GRS, Bearman PW, Meneghini JR.

On the wake-induced vibration of tandem circular cylinders: The vortex interaction excitation mechanism

.Journal of Fluid Mechanics, 2010, 661: 365-401

[本文引用: 3]     

[16] Assi GRS.

Wake-induced vibration of tandem and staggered cylinders with two degrees of freedom

.Journal of Fluids and Structures, 2014, 50: 340-357

[本文引用: 2]     

[17] Papaioannou GV, Yue DKP, Triantafyllou MS, et al.

On the effect of spacing on the vortex-induced vibrations of two tandem cylinders

.Journal of Fluids and Structures, 2008, 24: 833-854

[本文引用: 1]     

[18] Carmo BS, Assi GRS, Meneghini JR.

Computational simulation of the flow-induced vibration of a circular cylinder subjected to wake interference

.Journal of Fluids and Structures, 2013, 41: 99-108

[本文引用: 1]     

[19] Borazjani I, Sotiropoulos F.

Vortex-induced vibrations of two cylinders in tandem arrangement in the proximity-wake interference region

.Journal of Fluid Mechanics, 2009, 621: 321-364

[本文引用: 1]     

[20] 及春宁, 陈威霖, 黄继露.

串列双圆柱流致振动的数值模拟及其耦合机制

. 力学学报, 2014, 46(6): 862-870

[本文引用: 1]     

(Ji Chunning, Chen Weilin, Huang Jilu, et al.

Numerical investigation on flow-induced vibration of two cylinders in tandem arrangements and its coupling mechanisms

.Chinese Journal of Theoretical and Applied Mechanics, 2014, 46(6): 862-870 (in Chinese))

[本文引用: 1]     

[21] Yu KR, Étienne S, Scolan Y-M, et al.

Flow-induced vibrations of in-line cylinder arrangements at low Reynolds numbers

.Journal of Fluids and Structures, 2016, 60: 37-61

[本文引用: 1]     

[22] Werle H.

Sur L’ecoulement Autour D’um Faisceau Tubulaire

.Revue Francaise de Mechanique, 1972, 41: 7-19

[本文引用: 1]     

[23] Igarashi T, Suzuki K.

Characteristics of the flow around three circular cylinders arranged in line

.Bulletin of the Japan Society of Mechanical Engineers, 1984, 27(233): 2397-2404

[本文引用: 1]     

[24] Chen W, Ji C, Williams J, et al.

Vortex-induced vibrations of three tandem cylinders in laminar cross-flow: Vibration response and galloping mechanism

.Journal of Fluids and Structures, 2018, 78: 215-238

[本文引用: 9]     

[25] Ji C, Munjiza A, Williams JJR.

A novel iterative direct-forcing immersed boundary method and its finite volume applications

.Journal of Computational Physics, 2012, 231: 1797-1821

[本文引用: 3]     

[26] Chen W, Ji C, Xu W, et al.

Response and wake patterns of two side-by-side elastically supported circular cylinders in uniform laminar cross-flow

.Journal of Fluids and Structures, 2015, 55: 218-236

[本文引用: 2]     

[27] Chen W, Ji C, Wang R, et al.

Flow-induced vibrations of two side-by-side circular cylinders: Asymmetric vibration, symmetry hysteresis and near-wake patterns

.Ocean Engineering, 2015, 110: 244-257

[28] 及春宁, 花阳, 许栋.

不同剪切率来流作用下柔性圆柱涡激振动数值模拟

. 力学学报, 2018, 50(1): 21-31

(Ji Chunning, Hua Yang, Xu Dong, et al.

Numerical simulation of vortex-induced vibration of a flexible cylinder exposed to shear flow at different shear rates

.Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(1): 21-31 (in Chinese))

[29] 陈威霖, 及春宁, 徐万海.

并列双圆柱流致振动的不对称振动和对称性迟滞研究

. 力学学报, 2015, 47(5): 731-739

[本文引用: 1]     

(Chen Weilin, Ji Chunning, Xu Wanhai.

Numerical investigation on the asymmetric vibration and symmetry hysteresis of flow-induced vibration of two side-by-side cylinders

.Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(5): 731-739 (in Chinese))

[本文引用: 1]     

[30] Harimi I, Saghafian M.

Numerical simulation of fluid flow and forced convection heat transfer from tandem circular cylinders using overset grid method

.Journal of Fluids and Structures, 2012, 28: 309-327

[本文引用: 2]     

[31] Williamson CHK, Govardhan R.

Vortex-induced vibrations

.Annual Review of Fluid Mechanics, 2004, 36: 413-455

/