基于表面吸气的矩形截面涡激振动抑制及机理研究.pdf

第 卷第期 年月动 力 学 与 控 制 学 报J OUR NA LO FD YNAM I C SAN DC ON T R O LV o l N o J u n 文章编号:()D O I:/收到第稿,收到修改稿长沙理工大学大学生创新创业训练计划项目(),C h a n g s h aU n i v e r s i t yo fS c i e n c e&T e c h n o l o g yC o l l e g eS t u d e n t s I n n o v a t i v eE n t r e p r e n e u r i a lT r a i n i n gP l a nP r o g r a m()通信作者E m a i l:c e_h a n y a n c o m基于表面吸气的矩形截面涡激振动抑制及机理研究罗楚钰曾梦竹王响军韩艳(长沙理工大学 国际工学院,长沙 )(长沙理工大学 土木工程学院,长沙 )摘要随着桥梁跨度的增加,其风致振动问题逐渐凸显,发展简单、有效的风振控制措施对于保障桥梁结构的安全具有重要意义将桥梁断面简化为矩形,通过数值模拟研究了侧表面双气孔吸气对矩形截面的涡激振动的抑制作用研究表明:无量纲吸气流量是影响涡激振动抑制效果的决定性因素当无量纲吸气流量Qs 时,表面吸气对涡激振动的抑制效果达到最佳,矩形断面的振动幅值减小了 ,升力系数均方根值减小了 ,阻力系数平均值减小了 当Qs 时,矩形截面的振动幅值仍然被有效地抑制,但吸气流量过大,侧表面上的边界层厚度变小,边界层内的剪切效应变强,引起背风面尾流旋涡的增强,最终导致矩形截面升、阻力系数的增加关键词吸气控制,涡激振动,矩形断面,数值模拟中图分类号:U 文献标志码:AV I VS u p p r e s s i o na n dM e c h a n i s mA n a l y s i sB a s e do nt h eS u r f a c eA s p i r a t i o no fR e c t a n g u l a rS e c t i o nL u oC h u y uZ e n gM e n g z h uW a n gX i a n g j u nH a nY a n(I n t e r n a t i o n a lC o l l e g eo fE n g i n e e r i n g,C h a n g s h aU n i v e r s i t yo fS c i e n c e&T e c h n o l o g y,C h a n g s h a ,C h i n a)(S c h o o l o fC i v i lE n g i n e e r i n g,C h a n g s h aU n i v e r s i t yo fS c i e n c e&T e c h n o l o g y,C h a n g s h a ,C h i n a)A b s t r a c t W i t ht h ei n c r e a s eo fb r i d g es p a n,t h ep r o b l e mo fw i n d i n d u c e dv i b r a t i o nb e c o m e sm o r ea n dm o r ep r o m i n e n t T h ed e v e l o p m e n to f s i m p l ea n de f f e c t i v ew i n d i n d u c e dv i b r a t i o nc o n t r o lm e a s u r e s i so fg r e a t s i g n i f i c a n c e t oe n s u r e t h e s a f e t yo f b r i d g e s t r u c t u r e T h e s e c t i o no f t h eb r i d g e i s s i m p l i f i e da s a r e c t a n g l ea n d t h e s u p p r e s s i o no fV I Vo f t h e r e c t a n g u l a r s e c t i o nb ys u c t i o nw i t hd o u b l eh o l e s o n t h e s i d e s u r f a c e i ss t u d i e db yn u m e r i c a ls i m u l a t i o n T h er e s u l t ss h o wt h a tt h ei n s p i r a t o r yf l o wr a t ei st h ed e c i s i v ef a c t o ra f f e c t i n gt h es u p p r e s s i o ne f f e c to fV I V Wh e nt h ed i m e n s i o n l e s ss u c t i o nf l o wr a t eQs,t h es u r f a c es u c t i o nh a s t h eb e s ts u p p r e s s i o ne f f e c to nV I V T h ev i b r a t i o na m p l i t u d eo f t h er e c t a n g u l a rs e c t i o n i s r e d u c e db y ,t h er o o tm e a ns q u a r ev a l u eo f l i f t c o e f f i c i e n t i sr e d u c e db y ,a n dt h ea v e r a g ed r a gc o e f f i c i e n t i s r e d u c e db y Wh e nQs,t h ev i b r a t i o na m p l i t u d eo f t h e r e c t a n g u l a rs e c t i o n i ss t i l l e f f e c t i v e l ys u p p r e s s e d H o w e v e r,d u et ot h ee x c e s s i v e i n s p i r a t o r yf l o w,t h et h i c k n e s so f t h eb o u n d a r y l a y e ro n t h e s i d e s u r f a c eb e c o m e s s m a l l e r,a n d t h e s h e a r e f f e c t i n t h eb o u n d a r y l a y e rb e c o m e ss t r o n g e r,w h i c hl e a d s t o t h e e n h a n c e m e n t o f t h ew a k ev o r t e xo n t h e l e e w a r d s,a n d f i n a l l y l e a d s t ot h e i n c r e a s eo f t h e l i f t a n dd r a gc o e f f i c i e n to f t h er e c t a n g u l a rs e c t i o n K e yw o r d s s u c t i o nc o n t r o l,v o r t e x i n d u c e dv i b r a t i o n,r e c t a n g u l a rs e c t i o n,n u m e r i c a l s i m u l a t i o n第期罗楚钰等:基于表面吸气的矩形截面涡激振动抑制及机理研究引言钝体(如方柱、圆柱等)绕流常见于工程应用中,如风、水等流体作用下的输电线、斜拉索、高层建筑以及桥墩等当流体流过钝体时,会引发钝体尾流产生交替脱落的旋涡结构,当旋涡脱落频率与钝体固有频率满足一定关系时会引发大幅度振动,被称为“涡激振动”虽然涡激振动是一种限幅振动,但其常发生在低风速情况下,容易导致结构出现疲劳破坏,从而缩短结构的使用寿命因此,研究经济有效的控制方法来抑制钝体的涡激振动具有十分重要的工程应用价值在过去的几十年里,研究者提出了各种各样的方法抑制钝体构件的涡激振动,并减小构件所受的气动荷载根据是否需要外部额外提供能量,可以将控制措施分为两大类:主动控制方法和被动控制方法主动控制方法需要额外的能量输入,往往能够根据环境的变化自主调整控制参数,维持稳定的控制效果,但其控制装置通常较为复杂 W u等采用行波壁的方法成功地抑制了圆柱体尾流旋涡的交替脱落,使得圆柱体的平均阻力系数下降了 ,升力系数几乎降为零F e n g等 系统地研究了合成射流对圆柱体尾流旋涡交替脱落的抑制作用在圆柱体后驻点处布置射流,可以明显降低圆柱体的阻力,同时使得尾流旋涡的交替脱落完全被抑制C h e n等采用数值模拟方法研究了侧面稳定吸气对圆柱体表面 气 动 力 以 及 尾 流 旋 涡 交 替 脱 落 的 抑 制 作 用K o r k i s c h k o和M e n e g h i n i n在圆柱体两侧表面附近设置了旋转的小圆柱,为主圆柱的边界层注入动量,有效地延迟了边界层的分离,降低了阻力系数,同时也使得圆柱体的振动幅值减小了 B a o和T a o发现圆柱体尾流中沿流向做谐波振动的扁平水翼可以消除尾流的周期性脉动,同时阻力系数也减小了 最近,R e n等也在文章中指出圆柱体后方沿流向做谐波运动的分流板可以显著地抑制圆柱体的升、阻力系数以及横流向振动幅值与主动控制方法不同的是,被动控制方法的设备一般较为简单,且不需要额外的能量输入,但其控制效果会受外部环境的影响 A n d e r s o n和S z e w c z y k在圆柱体后驻点处设置刚性分离板来阻止尾流的横向流动,从而达到抑制尾流旋涡交替脱落的目的试验发现:分离板使得静止圆柱体表面动量边界层厚度减小了 当分离板的长度等于圆柱体直径时,顺流向阻力系数可以减小约 然而C u i等 通过研究发现:对于安装在弹性支座上的圆柱体而言,后驻点处的刚性分离板虽然可降低同步区的振动幅值,然而随着折算风速的继续增加,圆柱体的振动幅度却会线性增大,造成驰振的发生如果在后驻点处安装柔性分离板,不仅可以有效抑制圆柱体的涡激振动,还能避免驰振的出现L e e和K i m 在圆柱体表面等间距缠绕了三条螺旋线,圆柱体表面凸起的螺旋线引起近尾流区域内流动结构在展向上发生变化,从而减小了尾流旋涡结构在展向上的相关性,最终有效地抑制了圆柱体的横流向升力脉动Z h a n g等 采用大涡模拟方法研究了扰动波长为D的波浪矩形柱体在时尾流的流动特征,与普通矩形柱体相比,波浪矩形柱体的阻力系数平均值和升力系数脉动值分别减小了 和 N a w a z等 研究了倒角比()对方柱尾流非定常特性的影响,在R e 且时方柱的阻力系数减小了 H a mm a d等 研究了倒角尺寸的变化对长宽比为的矩形柱体的流动特性的影响倒角的增大可以在一定程度上降低阻力系数,抑制马蹄涡的形成在前人研究的基础上,本文将采用数值模拟方法研究侧表面双气孔稳定吸气对宽高比为的矩形柱体涡激振动的控制效果,并揭示其中潜在的物理机制数值模拟方法数值模型本文基于C F D(C o m p u t a t i o n a lF l u i dD y n a m i c s)商业计算软件AN S Y SF l u e n t 对宽高比为的矩形断面的涡激振动开展了二维数值模拟,并研究了侧表面双气孔吸气对横风向单自由度振动的矩形断面的涡激振动的抑制作用如图所示,矩形 断 面 的 宽 度(B)为 m,高 度(H)为 m,顺着流动方向,在矩形断面的两侧表面上各布置两个吸气孔根据无控情况下迎风面角点后方流动分离区的旋涡形成位置(见下文图),本文设置上(下)游吸气孔的中心离迎(背)风面距离等于矩形断面的高度H,使得上游气孔位于旋涡形动力学与控制学报 年第 卷成区域内,通过气孔吸气有效抑制大尺度旋涡的形成矩形断面可沿横流向振动,其质量(m)设置为 k g/m,固有频率(fnk/m,其中k为刚度系数)为 H z,阻尼比c/(m fn),其中c为阻尼系数为 图矩形断面示意图F i g S c h e m a t i co f t h er e c t a n g u l a r s e c t i o n计算域入口与矩形断面迎风面的距离为 B,出口与矩形断面背风面的距离为 B,计算域宽度为 B(阻塞比为)为解决动网格模拟中网格畸变、缠绕等问题,本文运用“刚性域动网格域静网格域”的方法对网格进行分块处理,如图所示在计算过程中,靠近矩形断面的网格跟随矩形断面一起运动,被称为刚性域,采用结构化网格划分刚性域外为动网格域,采用非结构化网格划分在结构振动过程中,动网格域内的网格会根据结构的振动幅度进行重构动网格域外为静网格域,采用结构化网格进行离散计算域入口边界条件设置为速度入口(v e l o c i t y i n l e t),出口为压力出口(p r e s s u r e o u t l e t),计算域两侧设置为对称边界条件(s y mm e t r y),矩 形 断 面 设 置 为 无 滑 移 壁 面(n o s l i pw a l l)图计算域网格划分示意图F i g S c h e m a t i co f t h eg r i db l o c k s i nt h ec o m p u t a t i o n a l d o m a i n本文采用k S S T方法求解雷诺平均方程组在每个时间步长内,程序先计算出整个计算域内的速度 场和压力场,然后通过UD F(U s e rD e f i n e dF u n c t i o n)提取矩形断面所受的横流向升力(Fy),接着采用N e wm a r k 方法求解矩形断面的运动方程(见式)m y c yk yFy()其中,y为矩形断面横流向位移,y和y 为矩形断面运动的速度与加速度,m、c和k分别是矩形断面竖向振动的等效质量、阻尼系数与刚度系数计算出下一时刻矩形断面的位置与速度后,程序对断面的位置与速度进行更新,同时刚性域的网格也跟随矩形断面同步运动,接着软件对动网格域内的网格进行更新,之后程序将进入下一个时间步的迭代计算,并如此循环往复工况设置在本文的数值模拟中,来流风速(U)选为 m/s,对应折算风速UrU/(fnH)(该风速下,矩形断面发生大幅度涡激振动,见下文节)文中吸气孔的位置保持不变,主要通过改变气孔的宽度(b)以及气孔的吸气速率(Us),从而找出该吸气方法中抑制矩形断面涡激振动的决定性控制参数根据吸气孔宽度的不同,本文将数值模拟工况分为三组,如表所示:表 数值模拟工况(其中b为气孔宽度,Us为气孔吸气速率)T a b l eN u m e r i c a l s i m u l a t i o nc o n d i t i o n s(bi sa i rh o l ew i d t h,Usi sa i rh o l es u c t i o nr a t e)工况b/HUs/US S S 数值模拟结果分析数值模拟有效性验证本文首先计算了无控矩形断面的风致振动幅值随来流折算风速的变化情况,并与刘志文等 的数值模拟结果以及M a t s u m o t o等 的风洞试验结果进行了对比,如图所示本文计算结果与参考文献中的结果十分吻合结合图可以看出:当尾流旋涡第期罗楚钰等:基于表面吸气的矩形截面涡激振动抑制及机理研究脱落频率(fv)小于断面振动固有频率(fn)时,fv随着折算风速(Ur)的增加而线性增大,此时矩形断面的振动幅值(A)较小当fv增加到与fn相接近时,fv基本稳定在fn左右,不再随Ur的增加而增大,即出现了涡激振动中著名的频率“锁定”现象 此时,断面的振动幅值显著增大当Ur继续增大并超过某一阈值后,fv又突然恢复到之前的线性增长趋势,此时矩形断面的振动幅值又迅速减小以上结果表明:本文数值计算准确地模拟出了矩形断面的涡激振动现象,并与前人的数值模拟结果以及风洞试验结果均十分吻合,从而验证了本文数值计算结果的准确性与可靠性图矩形断面振动幅值(A)随来流折算风速(Ur)的变化情况F i g T h ea m p l i t u d eo f ar e c t a n g u l a r s e c t i o n(A)v a r i e sw i t ht h er e d u c e dv e l o c i t i e s(Ur)图矩形断面旋涡脱落频率(fv)随Ur的变化情况F i g T h ev o r t e xs h e d d i n gf r e q u e n c y(fv)o f t h er e c t a n g u l a r s e c t i o n,v a r i e sw i t hUr振动幅值分析图给出了S、S以及S三组不同工况下矩形断面的振动幅值(A)随无量纲吸气速率(UsUs/U)的变化情况当无量纲吸气速率较小(如Us )时,矩形断面的振动幅值与无控情况相比反而略有增大,只有当Us超过某一阈值后,表面吸气才能对矩形断面的涡激振动产生抑制效果图表明:当气孔宽度(b)固定不变时,断面的A随着Us的增加而减小而当Us固定时,A随着b的增加而减小为了找出抑制矩形断面涡激振动幅值的决定性参数,图给出了矩形断面的振动幅值随无量纲吸气流量Qs(bUs)/(HU)的变化情况明显地,S 、S 以及S 三组工况下,图不同工况下,矩形断面振动幅值(A)随无量纲吸气速率(UsUs/U)的变化曲线其中A为无控矩形断面的振动幅值F i g T h ev i b r a t i o na m p l i t u d eo f t h e r e c t a n g u l a r s e c t i o n(A)v a r i e sw i t h t h e i n s p i r a t o r yr a t e(UsUs/U)u n d e rd i f f e r e n t c o n d i t i o n s Ai s t h ev i b r a t i o na m p l i t u d eo f t h eu n c o n t r o l l e dr e c t a n g u l a r s e c t i o n图不同工况下,矩形断面振动幅值(A)随无量纲吸气流量Qs(bUs)/(BU 的变化曲线F i g T h ev i b r a t i o na m p l i t u d eo f t h er e c t a n g u l a r s e c t i o n(A)v a r i e sw i t h i n s p i r a t o r yf l o wr a t eQs(bUs)/(BUu n d e rd i f f e r e n t c o n d i t i o n s动力学与控制学报 年第 卷矩形断面的A随Qs的变化曲线几乎完全重合,且在大部分区间内,矩形断面的随的增加而线性减小这反映出表面吸气对矩形断面涡激振动幅值A的控制效果由无量纲吸气流量Qs唯一决定,且在Qs ,范围内,A随着Qs的增加近似于线性减小气动力系数分析矩形断面侧表面上压力的周期性变化使得矩形断面受到横流向的升力(Fy)作用,从而引起矩形断面的横向振动而迎风面与背风面的压力差使得矩形断面受到顺流向的阻力(Fx)作用图给出了不同工况下,矩形 断面 的 升 力 系 数 均 方 根 值(c lr m s)与阻力系数平均值(c dm e a n)随Qs的变化情况与振幅A一样,各组工况下,矩形断面的c lr m s和c dm e a n随Qs的变化曲线完全重合结合图、图可以发现:当Qs小于 时,矩形断面的振动幅值、横流向升力系数均方根值以及顺流向阻力系数平均值均随着Qs的增加而减小Qs大于 后,矩形断面的涡激振动基本完全被消除然而,断面所受气动力系数随着的Qs增加而逐渐增大,其中阻力系数平均值甚至比无控情况下的阻力系数更大这反映出表面吸气控制方法对矩形断面的气动力的抑制效果并不是随着吸气流量的增加而单调增强在本文所计算工况中,Qs约等于 左右时,吸气控制方法获得最佳的控制效果,即断面的振动幅值与气动力系数均达到最小为了进一步揭示矩形断面的气动力系数随Qs变化的内在物理机制,图给出了S 情况下,Qs分别为、以及 三种工况(分别记为Qs、Qs和Qs)中矩形断面表面压力系数的分布情况,并对三种工况下矩形断面附近的旋涡结构的演化过程进行了追踪,如图所示在工况Qs(无吸气)情况下,来流在迎风面角点处发生流动分离,并导致角点后方分离区生成旋涡结构图(a)(h)显示:在一个振动周期内(以表面为例),上一个振动周 期 中 角 点 后 方 发 展 壮 大 的 旋 涡(记 为“V”)集聚了足够的能量,从而脱离角点的束缚开始向下游移动随之,在角点与V 之间的壁面附近会诱发出一对反向旋转的旋涡(分别记为V 和V,其中V 与V 的旋转方向相同),如图(b)所示随着V 继续向下游移动,V 逐渐衰弱直至消失;与之相反,V 因为不断地从上游来流中获取能量而逐渐增强,如图(d)所示之后,V 继续向下游运动,其强度逐渐衰减,最后在到达背风面角点附近时破碎成小尺度旋涡,并最终完全消散与此同时,V 持续从外层来流中吸收能量而发展壮大,最终发展为下一个振动周期的V 表面上旋涡的演变过程与表面相同,只是相位相差半个周期由于侧表面(表面与表面)上旋涡结构周期性地产生与消散,其上的压力脉动十分强烈因此,侧表面上的压力系数均方根值(c pr m s)是最大的,同时也使得迎风面上的c pr m s由中心向两侧迅速增大,而V 在下游角点附近的衰减与消散也使得侧表面尾部c pr m s迅速减小此外,尾流旋涡的交替脱落也造成了背风面上较大的压力脉动图不同工况下,矩形断面气动力系数随Qs的变化情况(a):升力系数均方根值(c lr m s);(b):阻力系数平均值(c dm e a s)F i g T h ea e r o d y n a m i cc o e f f i c i e n t sv a r yw i t hQsa td i f f e r e n t c a s e s(a):t h er o o t m e a ns q u a r eo f l i f t c o e f f i c i e n t(c lr m s);(b):t h em e a nd r a gc o e f f i c i e n t(c dm e a s)第期罗楚钰等:基于表面吸气的矩形截面涡激振动抑制及机理研究图S 工况组中三种不同吸气流量情况下矩形表面的压力系数分布(a):压力系数均方根值;(b):压力系数平均值其中工况Qs、Qs与Qs分别表示Qs/Q、和 插图为角点、的空间位置F i g P r e s s u r ec o e f f i c i e n td i s t r i b u t i o no nr e c t a n g u l a r s u r f a c eu n d e r t h r e ed i f f e r e n t i n s p i r a t o r yf l o wc o n d i t i o n s i nt h eS c o n d i t i o n(a):p r e s s u r ec o e f f i c i e n t r o o tm e a ns q u a r e;(b):a v e r a g ep r e s s u r ec o e f f i c i e n t H e r e,Qs、Qsa n dQsr e p r e s e n tQs/Q,a n d r e s p e c t i v e l y T h ei l l u s t r a t i o ns h o w s t h es p a t i a l p o s i t i o no f c o r n e rp o i n t,在工况Qs情况下,上游吸气孔的作用使得迎风面角点后方的旋涡结构无法发展壮大,因此稳定地依附在角点后方的侧表面上,不会脱离角点的束缚向下游移动因此,侧表面上压力的脉动被很好地抑制同时,下游吸气孔的吸气作用使得尾流变窄,背风面尾流中交替产生的旋涡也被明显削弱因此,如图所示:矩形断面侧表面上的压力系数的脉动几乎被完全抑制,背风面上的c pr m s也有大幅度的减小与无吸气工况(Qs)相比,工况Qs情况下,矩形断面的横流向升力系数均方根值减小了 ,而振动幅值下降了 在工况Qs情况下,随着吸气流量增加,上游吸气孔吸气对迎风面角点后方的旋涡结构的抑制增强与Qs情况相同,该旋涡结构无法摆脱角点的束缚,一直附着在侧表面上,使得矩形断面两侧的流动几乎对称分布由于吸气效应增强,大量外层高速流体被吸附到壁面附近,使得矩形截面侧表面的边界层厚度显著变小,这也导致了侧表面边界层内的剪切效应大幅增强侧表面边界层内被强化的剪切流到达背风面角点后发生分离,使得背风面尾流中旋涡结构交替地生成图(q)(x)中的涡量云图表明:Qs中背风面尾流中的旋涡强度远远大于Qs与Qs的尾流中的旋涡强度,因此,在工况Qs中,背风面的压力脉动比Qs与Qs中的都要强烈背风面尾流中压力的剧烈脉动向上游传递,从而使得侧表面后半部分的c pr m s显著增强因此,与工况Qs相比,工况Qs中矩形断面的c lr m s显著增大,如图(a)所示然而,c lr m s的增大却没有引起矩形断面振动幅值的明显增加(见图)这是由于吸气控制改变了矩形断面周围的流动形态,使得侧表面上压力的脉动频率发生了显著的变化,导致矩形断面的升力系数的波动频率远远大于其固有频率(见图),此时矩形断面无法发生涡激共振因此,虽然工况Qs中矩形断面的c lr m s达到了工况Qs中的 见图(a),然而其振动幅值仅为工况Qs中振动幅值的 综上所述,在工况Qs情况下,矩形断面侧表面的压力脉动导致了矩形断面的横风向振动,而侧表面的压力脉动是由于其上旋涡结构的交替产生、脱落与消散引起的在工况Qs情况下,上游吸气孔吸气使得迎风面角点后方的旋涡结构稳定地附着于侧表面上,从而极大地减小了侧表面的压力脉动同时,下游吸气孔吸气也在一定程度上抑制了背风面尾流的旋涡强度,使得背风面上的压力脉动大幅度减小因此,工况Qs中矩形断面的横风向振动被有效地抑制在工况Qs情况下,上游吸气孔吸气很好地抑制了迎风面角点后方旋涡结构的脱落(与工况Qs相同),使得侧表面前半部分的压力脉动大幅度减小然而,该工况下吸气流量过大,使得侧表面边界层厚度显著减小,边界层内的剪切效应增强,这使得背风面尾流中的旋涡强度明显增大,极大地增强背风面尾流的压力脉动尾流的压力脉动会向上游传递,从而导致侧表面后半部分的c pr m s显著增大 见图(a),因此矩形断面的c lr m s也随之变大 见图(a)动力学与控制学报 年第 卷第期罗楚钰等:基于表面吸气的矩形截面涡激振动抑制及机理研究图一个振动周期内矩形断面附近的流线与涡量演化过程其中第一列、第二列与第三列分别表示工况Qs、Qs和Qs从上至下分别表示t/T、/T、/T、/T、/T、/T、/T与T八个不同时刻F i g T h es t r e a m l i n ea n dv o r t i c i t ye v o l u t i o nn e a r t h er e c t a n g u l a r s e c t i o n i no n ev i b r a t i o np e r i o d T h e f i r s t t ot h i r dc o l u m n s r e s p e c t i v e l yr e p r e s e n t t h ew o r k i n gc o n d i t i o nQs,Qsa n dQs F r o mt o pt ob o t t o mr e s p e c t i v e l yr e p r e s e n t e i g h td i f f e r e n tm o m e n t s:t/T、/T、/T、/T、/T、/T、/Ta n dT图 Qs、Qs和Qs三种不同工况下矩形断面的升力系数(c l)时程曲线及频率谱其中fc l表示升力系数的脉动频率F i g T h e t i m eh i s t o r yc u r v e so f l i f t c o e f f i c i e n t s(c l)a n dt h e f r e q u e n c ys p e c t r u mo f r e c t a n g u l a r s e c t i o nu n d e r t h r e ed i f f e r e n tc o n d i t i o n so fQs,Qsa n dQs,fc lr e p r e s e n t s t h ep u l s a t i n gf r e q u e n c yo f l i f t c o e f f i c i e n t矩形断面的c dm e a n主要是由迎风面与背风面上的压差引起的如图(b)所示:在三种不同工况下,压力系数平均值(c pm e a n)在迎风面上的分布基本相同,而在背风面上却各不相同在工况Qs情况下,由于尾流旋涡的持续存在,消耗了尾流中的能量,从而使得背风面上的c pm e a n明显减小而在工况Qs情况下,由于下游吸气孔的吸气,使得尾流变窄,在一定程度上减小了背风面尾流旋涡的强度,从而使得背风面上的c pm e a n比工况Qs中的大,因此迎风面与背风面之间的压差变小,最终导致矩形断面c dm e a n的减小了 在工况Qs情况下,由于吸气孔吸气流量过大,致使外层的高速流体被吸入侧表面附近,使得壁面边界层厚度大幅度减小,而边界层内的剪切效应却显著增强当边界层内的流体在背风面角点处发生流动分离后,在尾流中形成的旋涡结构的强度明显高于工况Qs与Qs中的旋涡强度因此,工况Qs中背风面的c pm e a n最小,这造成了其c dm e a n比工况Qs与Qs中的大结论本文采用数值模拟研究了侧表面双气孔吸气对横风向单自由度振动的矩形断面的涡激振动的抑制作用通过改变吸气孔的宽度(b)与吸气速率(Us)发现,无量纲吸气流量(Qs)的大小决定了表面吸气方法对矩形断面涡激振动的控制效果当Qs 时,控制效果达到最佳,即矩形断面的振动幅值与气动力系数同时达到最小,此时矩形断面的升力系数均方根值比无控情况下的减小了 ,阻力系数平均值减小了 ,而振动幅值减小了 在无控情况(Qs)下,矩形断面的迎风面角点处发生流动分离,并在其后产生大尺度旋涡,该旋涡的周期性脱落造成了矩形断面侧表面上压力的强烈脉动,并最终导致矩形断面出现横风向涡激振动当Qs 时,由于上游吸气孔的吸气,迎风面角点处发生流动分离而产生的旋涡结构稳定地附着于侧表面上,从而使得矩形断面两侧的流动对称地向下游运动,因此极大地抑制甚至消除了矩形断面的横风向振动当Qs 时,上游吸气孔的吸气同样有效地抑制了迎风面角点后方旋涡结构的脱落,使得侧表面附近的流体平稳地向下游运动然而,由于吸动力学与控制学报 年第 卷气流量过大,使得外层高速流体进入壁面附近,导致侧表面上的边界层厚度显著减小,由此增强了边界层内的剪切效应当边界层内的流体运动到背风面角点处后发生流动分离,并在背风面尾流附近形成强有力的大尺度旋涡,而尾流旋涡的交替脱落引起了尾流压力的强烈脉动,并向上游传递使得侧表面后半部分压力系数的脉动明显增强,最终导致矩形断面的升力系数较Qs 时有显著增大当Qs 后,虽然矩形断面的升力系数逐渐增大,然而表面吸气致使矩形断面附近的旋涡结构发生了变化,从而使得升力系数的脉动频率远离了矩形断面的固有频率因此当Qs 后,矩形断面的振动幅值仍然维持在极低的水平本文仅研究了侧表面吸气对横风向单自由度振动的矩形断面的涡激振动的抑制作用,然而从本文的结果可以发现:矩形断面的侧表面上,压力的分布是极不均匀的,这将会引起矩形断面的扭转振动在未来的研究中,笔者将对横向振动与扭转振动耦合的双自由度振动的矩形断面的风致振动进行系统性研究参考文献WUCJ,W A N G L,WUJZ S u p p r e s s i o no f t h ev o nK r m nv o r t e xs t r e e tb e h i n dac i r c u l a rc y l i n d e rb yat r a v e l l i n g w a v eg e n e r a t e db yaf l e x i b l es u r f a c eJJ o u r n a l o fF l u i dM e c h a n i c s,:F E N G L H,WAN GJJ C i r c u l a rc y l i n d e rv o r t e x s y n c h r o n i z a t i o nc o n t r o l w i t has y n t h e t i cj e tp o s i t i o n e da tt h er e a rs t a g n a t i o np o i n tJ J o u r n a lo fF l u i dM e c h a n i c s,:F E N GLH,WANGJ J,P ANC P r o p e r o r t h o g o n a ld e c o m p o s i t i o na n a l y s i so fv o r t e xd y n a m i c so fac i r c u l a rc y l i n d e ru n d e rs y n t h e t i c j e tc o n t r o lJ P h y s i c so fF l u i d s,():F E NG L H,WANGJJ S y n t h e t i cj e tc o n t r o lo fs e p a r a t i o ni nt h ef l o wo v e rac i r c u l a rc y l i n d e rJE x p e r i m e n t s i nF l u i d s,():CHE N W L,C AOY,L IH,e ta l N u m e r i c a l i n v e s t i g a t i o no f s t e a d ys u c t i o nc o n t r o l o f f l o wa r o u n dac i r c u l a rc y l i n d e rJ J o u r n a lo fF l u i d sa n dS t r u c t u r e s,:KOR K I S CHKOI,ME N E GH I N IJR S u p p r e s s i o no fv o r t e x i n d u c e d v i b r a t i o n u s i n g m o v i n g s u r f a c eb o u n d a r y l a y e rc