低雷诺数变截面细长柔性结构流固耦合能量交换特征分析.pdf

收稿日期:修回日期:录用日期:基金项目:海南省自然科学基金(M S )作者简介:邓秀兵(),高级工程师,主要从事道桥设计研究,电话(T e l):;E m a i l:d e n g x i u b i n g f o x m a i l c o m文章编号:()D O I:/j c n k i j s j t u 低雷诺数变截面细长柔性结构流固耦合能量交换特征分析邓秀兵,于曰旻,庞玺源(浙江省水利水电勘测设计院有限责任公司,杭州 ;海南大学 土木建筑工程学院,海口 ;上海交通大学 船舶海洋与建筑工程学院,上海 )摘要:波浪形变截面柱体结构具有良好的流动减阻特性,然而尚不清楚这种变截面柔性结构的流致振动特性基于高精度谱单元法建立了细长柔性结构的流固耦合力学模型和高精度算法,对在低速均匀流作用和驻波初始动力激励下的流致振动机理进行了数值分析,阐明了波浪形细长柔性结构的尾流特性、结构动力响应特性、能量传递规律、涡脱频率的展向变化特征,也对减阻、减振机理进行了初步探讨数值模拟结果表明,在合适截面扰动波高下波浪形变截面柔性结构能够大幅抑制涡激振动响应,并发现在波浪形细长结构两侧形成的特殊涡结构会稳定绕流剪切层,拉长涡旋形成区长度,从而降低尾流结构与结构的流固耦合效应,抑制涡激振动响应关键词:波浪形截面;流致振动;细长柔性结构中图分类号:O 文献标志码:AA n a l y s i so f F l u i d S t r u c t u r eC o u p l i n gE n e r g yT r a n s f e rC h a r a c t e r i s t i c sS l e n d e rS t r u c t u r ew i t hV a r i a b l eC r o s s S e c t i o na tL o wR e y n o l d sN u m b e rD ENGX i u b i n g,Y UY u e m i n,P ANGX i y u a n(Z h e j i a n gD e s i g nI n s t i t u t eo fW a t e rC o n s e r v a n c ya n dH y d r o E l e c t r i cP o w e rC o,L t d,H a n g z h o u ,C h i n a;S c h o o l o fC i v i l a n dA r c h i t e c t u r a lE n g i n e e r i n g,H a i n a nU n i v e r s i t y,H a i k o u ,C h i n a;S c h o o l o fN a v a lA r c h i t e c t u r e,O c e a na n dC i v i lE n g i n e e r i n g,S h a n g h a i J i a oT o n gU n i v e r s i t y,S h a n g h a i ,C h i n a)A b s t r a c t:T h ew a v yd e f o r m e dc r o s s s e c t i o nc y l i n d r i c a l s t r u c t u r eh a se x c e l l e n tp r o p e r t i e so fd r a gr e d u c t i o ni nf l u i df l o w,b u t t h ef l o w i n d u c e dv i b r a t i o nc h a r a c t e r i s t i c so f f l e x i b l es t r u c t u r ew i t hs u c hv a r i a b l ec r o s s s e c t i o na r es t i l lu n c l e a r I nt h i sp a p e r,b a s e do nt h eh i g h p e r f o r m a n c es p e c t r a le l e m e n tm e t h o d,af l u i d s t r u c t u r ec o u p l e dm e c h a n i s t i cm o d e l a n dan u m e r i c a l a l g o r i t h mf o r s l e n d e r s t r u c t u r e sa r ee s t a b l i s h e d T h ew a k ec h a r a c t e r i s t i c s,s t r u c t u r a ld y n a m i cr e s p o n s e s,e n e r g yt r a n s f e r s,a n ds p a n w i s ev a r i a t i o n so fv o r t e xs h e d d i n gf r e q u e n c i e sa r ed i s c u s s e d T h en u m e r i c a l s i m u l a t i o nr e s u l t ss h o wt h a t s l e n d e r s t r u c t u r ew i t ht h ew a v y d e f o r m e dc r o s s s e c t i o nc a ng r e a t l ys u p p r e s st h ev o r t e x i n d u c e dv i b r a t i o nr e s p o n s ea ta na p p r o p r i a t ec r o s s s e c t i o nd i s t u r b a n c ew a v eh e i g h t,a n dt h es p e c i a l v o r t e xs t r u c t u r e f o r m e do nb o t hs i d e so f t h ew a v y 第 卷 第 期 年 月上 海 交 通 大 学 学 报J OUR NA LO FS HANGHA I J I A OT ON GUN I V E R S I T YV o l N o N o v 第 期邓秀兵,等:低雷诺数变截面细长柔性结构流固耦合能量交换特征分析 s h a p e ds l e n d e rs t r u c t u r ec a ns t a b i l i z et h ef l o wa r o u n dt h es h e a rl a y e ra n de l o n g a t et h ev o r t e xf o r m a t i o nl e n g t h,t h e r e b yr e d u c i n gt h ef l u i d s t r u c t u r ec o u p l i n ge f f e c tb e t w e e nt h ew a k es t r u c t u r ea n dt h es l e n d e rs t r u c t u r e,a n ds u p p r e s s i n gt h ev o r t e x i n d u c e dv i b r a t i o nr e s p o n s e K e yw o r d s:w a v yc r o s s s e c t i o n;f l u i d i n d u c e dv i b r a t i o n;s l e n d e rs t r u c t u r e在土木工程领域,风、浪、流等环境动力荷载作用下的细长结构物受损甚至破坏屡见不鲜比如,悬索桥中的主缆索作为典型的细长柔性结构,在来风作用下将产生涡激振动,甚至发生疲劳损坏;又如,在主缆索施工过程中,缆索截面形式发生非对称性变化,当风速超过临界值后,将产生空气动力负阻尼,使得缆索振动逐渐增强,甚至超负荷振幅而破坏因此,科学评估和有效抑制细长结构的流致动力响应是风、浪、流敏感细长柔性结构设计必须考虑的内容然而,当前细长柔性结构的流固耦合效应方面研究不足,尤其精准模拟和分析细长结构流致振动响应的力学模型和算法方面仍较缺乏,且亟待深入理解其流固耦 合机理并 建立系 统 的 分 析 方 法 与手段波浪形变截面圆柱体结构具有良好的流动减阻效果近年来,国内外学者对该类型结构的绕流特性开展了较多研究,发现波浪形变截面能够延后剪切层的相互作用,从而有效减小阻力 A h m e d等采用实验方法研究了波浪形圆柱的边界层分离线和尾流结构,发现波浪形圆柱流动分离线呈现明显的三维特性,在分离节点附近形成流向涡,且发生边界层“上卷”现象,从而延迟或抑制剪切层中湍流的生成和发展 L a m等利用多种流动显示技术对波浪形圆柱的近尾迹进行了实验研究,给出了平均速度和波动速度分量沿流向、展向和横向的分布特征实验结果表明,波浪形柱涡旋平均形成长度比大于光滑圆柱湍流统计分析也表明,光滑圆柱尾迹中的涡街更为规则,而由于波浪形圆柱后面的涡具有较强的三维效应,使得由湍流掺混增强而尾迹表现出更强的非相干流动结构目前变截面结构绕流分析主要集中在静止或弹性支承的刚性波浪形结构的减阻和涡激振动特性方面然而在实际土木工程中,缆索等细长结构的刚度系数通常较小,因此这种结构的柔性作用对流动减阻作用以及涡激振动响应的影响需要进一步研究,包括对初始动力激励的敏感性方面也需要深入分析在文献 的基础上,针对这一问题,本研究采用基于高精度谱单元方法的流固耦合分析方法,建立了细长变截面柔性圆柱体结构的涡激振动力学模型,对其在低速均匀流作用和驻波初始动力激励下的流致振动机理进行了量化分析,获得了包括波浪形结构的尾流特性、结构动力响应特性、能量传递规律、涡脱频率展向变化特征,并对其在驻波扰动下的减阻、减振机理进行了深入探讨,对波浪形变截面细长结构的工程设计和应用提供参考 数值算法细长结构的流致振动是一种典型的流固耦合问题求解流固耦合的基本方法通常有两大类,分别为界面匹配法和非界面匹配法,包括经典的任意拉格朗日欧拉(A L E)方法和浸入边界法(I BM)本文采用另一种由D i m a s等提出的随体坐标系(B o d y f i t t e dC o o r d i n a t e s)方法,该方法在惯性坐标系中求解N S方程,再通过坐标转换至非惯性坐标系中,从而不需要动网格和浸泡边界近似N S方程和连续性方程在惯性坐标系(x,y,z)内可描述为u t(u)u f p R eu()u()式中:u(u,v,w),表示在惯性坐标系内个方向的流速;t为时间;p 表示动压力;f表示流体的密度;表示惯性坐标系内的梯度算子;雷诺数R eUD/,U为均匀来流速度大小,D代表波浪形截面柱的平均直径,为流体黏性系数由惯性坐标系至非惯性坐标系的坐标转换采用如下关系式:xx x(z,t),yy y(z,t)()式中:x(z,t)和y(z,t)分别表示结构在顺流向和横流向的位移相应地,速度项和压力项采用以下变换:uu xtw xzvv ytw yzww pp()因此,将式()和()代入N S方程式()和()上海交通大学学报第 卷可得:ut(u)u fpR euA(u,p,)()u()式中:为微分算子;附加项A(u,p,)是由坐标转换引起的附加加速度项,该值在个坐标方向的计算通过下式得到,AxdxdtR ez(uxzw)uzxzx ywxtzAydydtR ez(vyzw)vzyzx ywytzAzfxzpxfyzpyR e(wzwz)()在此定义:ddttuxvywzz zxzxyzyx yxy()此外,本文采用线性张力梁模型描述细长柔性结构的动 力学行为,该模型采用 小变形假设,可写成:ctktTzF(z,t)()式中:c,k和T分别表示结构的单位长度的质量、阻尼比和张力需要指出的是,张力T的大小将影响相速度c,c(T/c)/;(z,t)(x,y)表示结构在顺流向和横流向的位移;F(z,t)是作用于结构上的流体力,通过对压力和黏性力项沿结构固壁表面积分获得,F(z,t)p n(uuT)nds()n为指向结构外向的法线单位向量;s为结构表面微分此外,假设结构的动力响应满足沿展向的周期性条件,则 对 流 场 和 结 构 变 量 可 采 用 傅 里 叶 级 数(F o u r i e rE x p a n s i o n)表示,即u(x,y,z,t)v(x,y,z,t)w(x,y,z,t)p(x,y,z,t)M mum(x,y,t)vm(x,y,t)wm(x,y,t)pm(x,y,t)eim z()x(z,t)y(z,t)M mx,m(t)y,m(t)eim z()式中:/Lz表示展向的波数,Lz为结构长度;M为展开式中傅里叶模态个数;m为各阶模态将式()和()代入式()和()可获得解耦后的二维模态方程:umt(u)um fpmR e(x ym)umAm()um()式中:x,y,im相应地,对应于每阶模态m,结构的解耦常微分方程表述为cmc mmT mFm(t)()对结构运动的解耦方程式(),采用二阶N e wm a r k 方法进行求解对不定常流场模态方程式()和(),采用K a r n i a d a k i s等提出的高阶分步格式进行时间离散具体地,对每一时间迭代步对速度和压力进行解耦计算第步,考虑非线性项(u)u和A(u,p,),通过J阶强稳定化算法获得中间速度场um:umJ qqunqmtJ qq(u)uAnqm()式中:q和q均为与强稳定化积分(S t i f f l yS t a b l eI n t e g r a t i o n)有关的参数第步,将考虑压力梯度的作用修正速度场,并施加连续性约束条件和纽曼边界条件:umumtpm()第步,考虑黏性项更新下一步的速度场:un mumtR e(x ym)um()式中:为强稳定计算过程中向后差分系数(B a c k w a r d sD i f f e r e n t i a t i o nC o e f f i c i e n t)对上 述 时 间 离 散 后 的 模 态 方 程,根 据K a r n i a d a k i s等 将(x,y)平面的二维空间域离散为四边形 有 限 元 单 元,并 采 用 高 斯洛 巴 托勒 让 德(G L L)高阶正交拉格朗日多项式作为形函数进行空间离散 数值模型 物理模型考虑在均匀来流作用下的波浪形变截面圆柱第 期邓秀兵,等:低雷诺数变截面细长柔性结构流固耦合能量交换特征分析 (见图),其直径沿展向变化由下式确定:DzDAc o s z()式中:Dz为展向相应位置的圆截面直径由于波浪形柱 直 径 在 展 向 呈 余 弦 变 化,平 均 直 径 为D(Dm a xDm i n)/基于平均直径和均匀来流速度的雷诺数取值为 ,因此本文考虑的流动为层流流动定义波浪形柱直径最大截面所在的展向位置为几何节点,波浪柱直径最小截面所在的对应位置为几何鞍点A代表波 状表面的波 高,取值范围 为 D D;展向波长设置为 DA对应于光滑圆柱,也作为基准工况用于与波浪形柱工况进行对比分析参照N e wm a n等 的研究,可以通过规定圆柱的初始振幅和速度来确定其初始条件,本 研 究 中 初 始 扰 动 为 驻 波:y(z,t)ac o s(t)c o s(z/Lz),其中a为振幅,为振动频率,/(c Lz)图均匀流下截面直径沿展向余弦变化的波浪柱几何示意图F i g S c h e m a t i cv i e wo fu n i f o r mf l o wp a s tw a v yc y l i n d e rw i t ht h ed i a m e t e rv a r y i n ga l o n gt h es p a n w i s ed i r e c t i o n二维平面网格划分如图所示计算域几何形状为C形,由圆柱上游半径为 D的半圆弧和下游矩形组成,其流向长度为 D,在横流向长度设置为 D,圆柱放置在计算域横流向的中心位置在(x,y)平面内采用了 个基础网格;由于对每一基础单元再用阶正交多项式(P)对速度和压力场进行近似,所以每个基础网格再细分为个高阶网格此外,根据文献 ,沿展向采用 个傅里叶模态,满足本文雷诺数下的精度要求流体计算域的边界条件设置如下:入口边界为均匀来流条件,来流速度为(U,);在出口边界上采 用 纽 曼 边 界 条 件,令 出 口 边 界 上u/x图二维网格拓扑结构F i g Dm e s ht o p o l o g yv/xw/x以及p在上、下部边界上指定(u,v,w)(U,)的远场边界条件;在圆柱固壁采用无滑移和不可渗透边界条件,即流固耦合边界条件tv通过假设无限长波浪柱的展向周期性,在波浪柱的两端实施了只允许横向运动的自由边界条件 模型验证为了验证力学模型和数值模拟方法的适用性,首先获得了均匀来流作用下细长光滑圆柱涡激振动模拟结果,并与已有文献结果进行对比分析本模拟中质量比设置为m ,而其他模拟参数均与N e wm a n等 相同图所示为光滑圆柱在驻波和行波初始激励下的横向振幅(y)、阻力系数(Cd)和升力系数(Cl)模拟结果,与N e wm a n等 结果吻合良好尤其在驻波初始激励下本模拟结果与N e w m a n等 结果在展向节点(响应振幅最大的点)和反节点(响应振幅最小的点)位置坐标完全相同在行波初始激励下的行波响应幅值和行波速度等均与N e wm a n等 结果一致,故说明本文力学模型和模拟方法的可靠性接着正交多项式阶数对计算结果的影响进行验证分析通过变化多项式阶数P改变高阶网格的疏 上海交通大学学报第 卷图柔性光滑圆柱涡激振动拟结果F i g T e m p o r a l a n ds p a t i a l e v o l u t i o nc o n t o u ro f t r a n s v e r s ev i b r a t i o na n dh y d r o d y n a m i cc o e f f i c i e n to f s m o o t hc y l i n d e r密;在展向上通过改变傅里叶模态个数,改变展向网格的疏密由计算结果表明:P,M 时模拟结果足够收敛至精确解图y、Cd和Cl沿展向的均方根值分布F i g D i s t r i b u t i o no f r o o tm e a ns q u a r e(RM S)v a l u eo fy,Cd,a n dCl 结果与讨论 运动响应与水动力系数图所示为y、Cd和Cl沿细长结构展向的均方根值模拟结果阻力系数和升力系数分别定义为CdFd/(fUD)和ClFl/(fUD),其中Fd和Fl分别代表结构受到的顺流和横流向流体力分量随着扰动波高A的增加,结构的横向振幅逐渐减小,尤其是在A D时,结构振动几乎完全被抑制,与无表面扰动或光滑表面结构相比减小了 另一方面,当A D时横向振幅沿展向分布特性与光滑表面结构类似,并且幅值基本相同阻力系数均方根值,也随着扰动波高的增加而减小,并且在A D时,阻力系数相比于光滑表面结构明显降低;而当A D时其分布特征类似,但在展向各位置处的阻力系数均有小幅减小升力系数均方根计算结果显示,当A D时,升力系数相比于光滑表面结构有所升高,而A D时,升力系数会显著下降由此可知,变截面结构的表面扰动波高对振动幅值和水动力系数均有显著影响,并且在A D时,结构响应得到有效抑制,且水动力特性得到明显改善图为在A D,D时的结构横向振动响应时程模拟结果由图可知,两种工况下节点和反节点的展向位置与光滑截面结构相同,即节点位于z/D,反节点位于z/D,处,表明该工况下结构振动模式均与光滑表面结构相同阻力系数在A D时其随时间演化特征与光滑表面结构基本相同,均在节点处取得最大值,在反节点处取最小值然而A D工况下的阻力系数时空分布则与光滑表面结构完全不同,在阻力系数第 期邓秀兵,等:低雷诺数变截面细长柔性结构流固耦合能量交换特征分析 图结构响应与水动力系数计算结果F i g T e m p o r a l a n ds p a t i a l e v o l u t i o no fv i b r a t i o na n dh y d r o d y n a m i cc o e f f i c i e n t s分布中看不到明显的驻波效应在z/D,时阻力系数并未取得最大值,且最大阻力系数仅有光滑表面结构的,表明A D扰动波高具有良好的减阻效果进一步观察升力系数时程可知,A D工况下波形柱在展向z/D 范围内的升力系数较光滑表面结构有所下降,而在其他展向位置处升力系数有所升高,这与局部雷诺数的增大有关而当A D时升力系数时空演化特征发生显著改变,不仅最大升力系数绝对值明显降低,且在两个节点中间的拓扑结构发生显著变化根据波形柱在驻波初始激励下对横向振幅的抑制效果图展向涡量瞬时等值面(z )的透视图和横流向俯视图F i g P e r s p e c t i v ea n dc r o s s f l o wv i e wo f t h ei n s t a n t a n e o u s i s o s u r f a c eo f t h es p a n w i s ev o r t i c i t y(z )和对升、阻力特性的影响特征,大致确定出控制失效范围(A D)以及控制有效范围(DA D)尾流特性为了进一步了解扰动表面波高对结构振动响应的影响机理,需要对柔性波形柱和柔性光滑柱近尾迹三维涡结构进行对比分析图给出了不同扰动波高时展向涡量等值面(z )计算结果其中,图(a)(c)为柔性光滑柱展向涡量等值面透视图,图(d)(f)为俯视图,与N e wm a n等 计算结果吻合良好驻波初始激励下柔性光滑柱展向涡量中存在明显的交织结构这种交织涡结构的展向特征是,在z/D,截面处形成交错脱落的卡门涡街结构;而在z/D,截面处上、下表面同时脱落的完全对称型涡结构图中分别给出A D时模拟结果,显然在z/D,截面处的涡结构也呈现对称分布值得注意的是,在此工况下几何节点处的展向涡量幅值较大从而升力系数幅值也随之增大,而几何鞍点处涡量幅值较小从而升力系数幅值也随之减小,这与图所示结果一致图中分别给出A D时模拟结果,在该工况下 上海交通大学学报第 卷交织涡结构完全消失,弯曲涡管以交错方式从波浪形表面上、下侧分别脱离,并且在尾流中迅速消失更重要的是,与柔性光滑柱相比,此时波形柱表面剪切层的卷曲和相互作用明显较弱,使得涡形成区长度进一步增大 L i n等 用数值方法研究了亚临界雷诺数条件下具有相对较大展向波长的刚性波形柱周围的流动特征,也得到了类似的结果图种工况下的顺流向(x )与横流向(y )瞬时等值面F i g I n s t a n t a n e o u si s o s u r f a c eo ff l o wd i r e c t i o n(x )a n dc r o s sf l o wd i r e c t i o n(y )u n d e rt h r e ew o r k i n gc o n d i t i o n s图种工况下(x,y)平面内圆柱周围的瞬时等值面和相应的流向涡量示意图F i g S c h e m a t i cd i a g r a m o fi n s t a n t a n e o u si s o s u r f a c ea n dc o r r e s p o n d i n gs t r e a mw i s ev o r t i c i t ya r o u n dac y l i n d e r i nt h ep l a n e(x,y)u n d e r t h r e ew o r k i n gc o n d i t i o n s图所示为不同扰动波高条件下的涡量场顺流向(x )与横流向(y )分量瞬时等值面云图从该图可清晰分辨出顺流向涡结构和横流向涡结构在尾迹附近的演化特征显然,与柔性光滑柱和A D扰动波高工况相比,A D时结构附近涡结构完全不同且涡量强度相对较低,并且沿下游耗散较快 L i n等 指出,附加流向涡在波形柱表面几何鞍点和节点处分别诱导出上涌流和下涌流,使得在几何鞍点处产生宽尾流,而在几何节点处产生窄尾流这与本文模拟结果非常一致,表明波形结构的涡激振动抑制和减阻效应的力学机理有相同之处为了进一步研究流向涡对结构响应的影响,图中给出了种工况下(x,y)平面内流向涡量切片图;其中图中分别为柔性光滑柱和在A D,D时柔性波形柱的模拟结果该结果表明,柔性光滑柱两侧和A D波形柱两侧均形成一对同向旋转涡结构,而沿结构展向上、下侧分布的两对涡旋转方向刚好相反然而,在A D情况下,在同一展向位置处柔性波形柱两侧生成一对较强的反向旋转涡,在这对反向旋转涡的外侧还分布着一对与各自强涡旋转方向相反的弱涡 L i n等 指出,最优控制波形柱附加产生的反向流向涡对剪切层稳定性具有重要作用,能够防止剪切层与强涡结构的相互作用,使得涡旋形成长度增大在A D的情第 期邓秀兵,等:低雷诺数变截面细长柔性结构流固耦合能量交换特征分析 况下观察到的反向旋转涡旋也稳定了剪切层,拉长了涡旋形成区长度,与本文之前观察到的现象一致 能量转换特征与频谱特性结构的运动响应与流体与结构之间的能量转化特征密切相关,根据N e wm a n等 的研究,一个旋涡脱落周期的无量纲时均能量E(z)可定义为E(z)fUDTWldt()式中:代表无量纲脱落周期;Wl代表升力产生能量的功率,即Wly/t当E(z)为正值时,代表能量从流体转移至结构,E(z)为负值则转移方向相反 N e wm a n等 和Z h u等 研究指出对于驻波响应,E(z)取值与展向位置有关,而与无量纲时间无关同时,E(z)在一个脱落周期内沿整数波长的积分应该为,表明在一个旋涡脱落周期内,结构从流体中得到的能量和在流体中做功耗散的能量相同,在没有外界能量输入的前提下,整个流体和结构系统满足能量守恒因此,时均能量代表结构与流体进行能量交换的强度,其值越大则相应结构与流体之间的能量交换越频繁,反映到实际工程中,就会引起结构的疲劳,加速结构的老化图 光滑圆柱无量纲时均能量分布F i g S m o o t hc y l i n d e rp o w e r s p e c t r a l d e n s i t y图所示为所有工况下时均能量随展向分布计算结果由图可知,A D柔性波形柱的时均能量比柔性光滑柱显著增大,当扰动波面高度增至A D其时均能量相比于柔性光滑柱维持在同一水准而当A D时,E(z)始终在附近小幅振荡,说明在上述种工况下波形柱与流体之间的能量交换接近于,从能量角度验证了波形柱与流体之间的动力响应被显著抑制同时,在z/D 以及z/D 之间柔性光滑柱E(z)为正值,在其余位置E(z)为负值;而A D波浪形结构在两端E(z)为正值,而在跨中附近E(z)为负值,这与图中升力系数的时空演图细长柔性结构涡激振动时的无量纲时均能量分布F i g N o n d i m e n s i o n a lm e a nt i m ee n e r g yd i s t r i b u t i o no fs l e n d e rs t r u c t u r eu n d e r g o i n gv o r t e x i n d u c e dv i b r a t i o n化模拟结果相互一致完全消除初始瞬态效应后,对种典型工况下横向位移和升力系数时程进行了频域分析图 中分别为柔性光滑柱横向位移和升力系数功率谱密度(P S D)展向分布计算结果,可以据此确定 旋 涡 脱 落 频 率 为f D/U ,此 外,在P S D的展向分布云图中可观察到,谱峰幅值与圆柱展向波长密切相关 B o u r g u e t等 研究了长径比为 的细长柔性结构涡激振动问题,并也报道了类似现象由此可以进一步证实,结构的振动响应与波长相关此外,阻力频率(这里没有显示)是升力频率的两倍,说明尾迹是经典的 S模式图 中分别为A D柔性波浪柱横向位移和升力系数P S D展向分布模拟结果,可以得知其主频为f D/U ,略小于光滑柱时相应频率从图(b)可以看到,与光滑柱相比,A D柔性波形柱振动时程在跨中对应的P S D值明显减小,表明该位置附近升力的波动强度显著减小,与图显示的现象一致图 中分别为A D柔性波形柱横向位 上海交通大学学报第 卷图 A D波浪柱无量纲时均能量分布F i g W a v yc y l i n d e rp o w e r s p e c t r a l d e n s i t ya tA D图 A D波浪柱无量纲时均能量分布F i g W a v yc y l i n d e rp o w e r s p e c t r a l d e n s i t ya tA D移和升力系数P S D的展向分布计算结果对于A D柔性波形柱,除了主频f D/U 外,在f D/U 处也可识别出次峰频率,在图 用白色虚线表示这种低频次生频率,在P S D Cl中比在P S D y中更强次生频率产生的原因涉及流固耦合的非线性问题,其机理有待进一步探讨 结论采用基于高精度谱单元方法的流固耦合分析方法,建立了细长柔性波形柱体结构的涡激振动力学模型,对其在低速均匀流作用和驻波初始动力激励下的流致振动机理进行了量化分析结果显示,当扰动波高A D时,柔性波形柱横向振幅相比于光滑柱明显减小,升阻力系数显著降低,水动力特性得到明显改善,表明在合适波高下,柔性波形柱具有良好的流致振动抑制作用进一步比较扰动波高分别为A D,D的结构横向振幅和水动力系数的时空演变特征,大致确定出波形柱抑制涡激振动的控制失效范围(A D)以及控制有效范围(DA D)为了探究扰动表面波高对结构振动响应的影响机理,对扰动表面结构和无扰动表面结构近尾迹三维涡结构进行了对比分析结果表明,与柔性光滑柱相比,A D工况下交织涡结构完全消失,使得波形扰动截面剪切层的卷曲和相互作用明显较弱,涡形成区长度进一步增大对柔性光滑柱在A D,D时波形柱的流向涡结构瞬时等值面进行了模拟比较,结果表明在A D的情况下观察到的反向旋转涡旋也稳定了剪切层,拉长了涡旋形成区长度比较了所有工况下时均能量的展向分布计算结果,发现当A D时,E(z)始终在附近小幅振荡,从能量转移角度进一步验证柔性波形柱对涡激振动的抑制作用在完全消除初始瞬态效应后,对种典型工况下横向位移和升力系数时程进行了频域分析,确定了种工况下的漩涡脱落频率,对于A D波形柱工况下观察到的次生频率,其产生的原因涉及流固耦合的非线性问题,有待对其机理进行进一步探讨参考文献:陈伟民,付一钦,郭双喜海洋柔性结构涡激振动的第 期邓秀兵,等:低雷诺数变截面细长柔性结构流固耦合能量交换特征分析 流固耦合机理和响应J力学进展,():CHE N W e i m i n,F U Y i q i n,GUOS h u a n g x i R e v i e wo nf l u i d s o l i dc o u p l i n ga n dd y n a m i cr e s p o n s eo fv o r t e x i n d u c e dv i b r a t i o no fs l e n d e ro c e a nc y l i n d e r sJA d v a n c e s i nM e c h a n i c s,():AHME DA,KHAN MJ,B AY S MU CHMO R EBE x p e r i m e n t a li n v e s t i g a t i o n o f a t h r e e d i m e n s i o n a lb l u f f b o d yw a k eJA I A AJ o u r n a l,():L AM K,WAN GFH,S ORMC T h r e e d i m e n s i o n a ln a t u r eo f v o r t i c e s i nt h en e a rw a k eo f aw a v yc y l i n d e rJJ o u r n a lo fF l u i d sa n dS t r u c t u r e s,():Z HU HB,Z HOU D,B AO Y,e t a lW a k ec h a r a c t e r i s t i c so f s t a t i o n a r yc a t e n a r y r i s e r sw i t hd i f f e r e n t i n c o m i n gf l o wd i r e c t i o n sJ O c e a nE n g i n e e r i n g,:Z HU H B,WANG R,B AO Y,e t a lF l o wo v e ras y mm e t r i c a l l yc u r v e dc i r c u l a rc y l i n d e rw i t ht h ef r e es t r e a mp a r a l l e l t o t h ep l a n eo f c u r v a t u r ea t l o wR e y n o l d sn u m b e rJJ o u r n a lo fF l u i d sa n dS t r u c t u r e s,:WAN GR,Z HU H B,B AO Y,e t a lM o d i f i c a t i o no ft h r e e d i m e n s i o n a li n s t a b i l i t yi nt h ep l a n a rs h e a rf l o wa r o u n dt w oc i r c u l a rc y l i n d e r si nt a n d e mJP h y s i c so fF l u i d s,():WAN GR,B A OY,Z HOUD,e t a l F l o wi n s t a b i l i t i e s i n t h ew a k eo f a c i r c u l a r c y l i n d e rw i t hp a r a l l e l d u a l s p l i t t e rp l a t e sa t t a c h e dJ J o u r n a lo fF l u i d M e c h a n i c s,:D I MA SA A,T R I AN T A F Y L L OU GS B a r o c l i n i c b a r o t r o p i ci n s t a b i l i t i e so ft h eG u l fS t r e a me x t e n s i o nJ J o u r n a l o fP h y s i c a lO c e a n o g r a p h y,():KA R N I A D AK I SGE,I S R A E L IM,O R S Z A GSAH i g h o r d e rs p l i t t i n g m e t h o d sf o rt h ei n c o m p r e