殷 彤， 及春寧， 宋立群， 袁德奎

(1. 天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300350; 2. 天津大學(xué) 機(jī)械工程學(xué)院，天津 300350)

當(dāng)水流經(jīng)過具有陡峭后緣的結(jié)構(gòu)時(shí)，會(huì)在結(jié)構(gòu)的背流側(cè)發(fā)生交替漩渦的生成與瀉放現(xiàn)象。與之相聯(lián)系，結(jié)構(gòu)受到橫流向(垂直于來(lái)流方向)周期性的流體力作用，激發(fā)結(jié)構(gòu)的渦激振動(dòng)(vortex-induced vibration，VIV)。關(guān)于均勻圓柱的渦激振動(dòng)問題，前人已經(jīng)取得了大量卓有成效的研究成果，對(duì)單個(gè)和多個(gè)均勻圓柱的振動(dòng)特性、尾流模式、干涉機(jī)理等進(jìn)行了全面而細(xì)致的詮釋，讀者可查閱綜述文獻(xiàn)[1-8]。然而實(shí)際工程中，圓柱狀結(jié)構(gòu)沿軸向的直徑往往不等，如海洋工程中常見的加裝浮力模塊的海洋立管，其外形是直徑不連續(xù)變化的階梯圓柱。

對(duì)于階梯圓柱的尾流和渦激振動(dòng)問題，當(dāng)前研究主要集中在：①直徑的不連續(xù)變化對(duì)階梯圓柱尾渦結(jié)構(gòu)的影響[9-13]；②直徑比、長(zhǎng)細(xì)比等因素對(duì)階梯圓柱渦激振動(dòng)的影響；③階梯圓柱不等直徑段之間的振動(dòng)激勵(lì)競(jìng)爭(zhēng)機(jī)制。

Dunn等結(jié)合頻譜分析與流場(chǎng)圖像，發(fā)現(xiàn)階梯圓柱后方沿展向形成了3個(gè)脫渦單元，即：小直徑段后方的S單元、大直徑段后方的N單元(靠近變截面處)和L單元(遠(yuǎn)離變截面處)。N單元的形成與變截面處的“下洗”運(yùn)動(dòng)有關(guān)，其脫渦頻率比S單元和L單元的脫渦頻率低。Morton等依據(jù)大直徑段的長(zhǎng)細(xì)比L/D(L和D分別為大直徑段的長(zhǎng)度和直徑)，進(jìn)一步對(duì)階梯圓柱的尾渦模式進(jìn)行劃分。當(dāng)L/D≥15時(shí)，階梯圓柱尾流中存在3個(gè)脫渦單元；當(dāng)8

對(duì)于彈性支撐的階梯圓柱，其渦激振動(dòng)特征與均勻圓柱的明顯不同。Ji等[14]發(fā)現(xiàn)，當(dāng)直徑比D/d=1.43時(shí)(D和d分別為階梯圓柱的大、小直徑)，階梯圓柱的振幅比均勻圓柱的減小了19%；當(dāng)直徑比D/d=2.0，大直徑段振動(dòng)鎖定時(shí)，小直徑段振動(dòng)可能處于非鎖定狀態(tài)，并且這種鎖定不同步的現(xiàn)象隨著直徑比的降低而逐漸消失。在不同直徑段振動(dòng)激勵(lì)競(jìng)爭(zhēng)機(jī)制研究方面，Lie等[15]發(fā)現(xiàn)：不同于均勻圓柱只受單個(gè)脫渦頻率激勵(lì)，加裝浮力塊的階梯立管同時(shí)受到兩個(gè)頻率的激勵(lì)作用。Vandiver等[16]進(jìn)一步指出：當(dāng)浮力塊覆蓋率超過50%時(shí)，立管振動(dòng)響應(yīng)由大直徑浮力塊的脫渦頻率主導(dǎo)。Rao等[17]提出了用以判斷激勵(lì)競(jìng)爭(zhēng)結(jié)果的預(yù)測(cè)模型，并認(rèn)為長(zhǎng)細(xì)比L/D是影響立管疲勞損傷率的重要因素。孫友義等[18-19]通過數(shù)值模擬手段，研究了不同浮力塊分布形式下的立管渦激振動(dòng)疲勞特性，對(duì)浮力塊配置方案進(jìn)行優(yōu)化。

綜上可知，已有針對(duì)階梯圓柱的研究成果多關(guān)注直徑比和長(zhǎng)細(xì)比對(duì)尾流模式的影響，而對(duì)階梯圓柱的渦激振動(dòng)特性研究較少。在海洋工程應(yīng)用背景下，相關(guān)學(xué)者開展了帶浮力塊細(xì)長(zhǎng)立管渦激振動(dòng)特性研究，但對(duì)振動(dòng)機(jī)理缺少深入的分析，對(duì)直徑比和長(zhǎng)細(xì)比的影響缺少系統(tǒng)的研究。為此，本文作者在低湍流度水槽中開展了均勻來(lái)流條件下階梯圓柱渦激振動(dòng)的試驗(yàn)研究，探究階梯圓柱渦激振動(dòng)的振動(dòng)特性和激勵(lì)機(jī)制，分析長(zhǎng)細(xì)比、直徑比對(duì)階梯圓柱渦激振動(dòng)響應(yīng)的影響。

1 試驗(yàn)裝置

1.1 階梯圓柱模型

階梯圓柱模型由均勻裸管與外套空心圓管組合而成，如圖1所示。均勻裸管總長(zhǎng)436.8 mm(其中淹沒段長(zhǎng)度l=272.8 mm)，外徑d=16 mm，質(zhì)量為71.5 g。裸管為兩端封閉的有機(jī)玻璃圓管，力學(xué)性能良好，在試驗(yàn)條件下可視作剛性圓柱。外接空心圓管由8200樹脂材料通過3D打印制成，材料密度為1.12 g/cm3。單個(gè)外接空心圓管長(zhǎng)度為34.1 mm，內(nèi)徑為16 mm，外徑D=20 mm,24 mm,32 mm，直徑誤差小于0.2 mm。通過改變外接空心圓管的個(gè)數(shù)，可組成不同長(zhǎng)細(xì)比的階梯圓柱。

1.2 試驗(yàn)設(shè)施與數(shù)據(jù)采集設(shè)備

圖2 試驗(yàn)裝置圖Fig.2 Sketch of experimental facilities

圖3 各測(cè)點(diǎn)位置的平均流速和脈動(dòng)流速均方根值Fig.3 Time averaged velocity and root-mean-square value of fluctuating velocity at different measuring points

試驗(yàn)設(shè)備安裝如圖2所示。振動(dòng)和測(cè)量系統(tǒng)支架由80 mm×80 mm鋁型材制成，力學(xué)性能良好，并與地面牢固連接。振動(dòng)和測(cè)量系統(tǒng)支架與試驗(yàn)水槽分離，避免水泵對(duì)振動(dòng)信號(hào)的干擾。為盡可能減小振動(dòng)系統(tǒng)的結(jié)構(gòu)阻尼，圓柱模型安裝在直線氣浮導(dǎo)軌(Newway S40-03075-038381)上。氣浮導(dǎo)軌由微孔介質(zhì)滑塊和鋁合金導(dǎo)軌構(gòu)成，滑塊與導(dǎo)軌之間形成氣膜，無(wú)直接接觸，滑動(dòng)摩擦小?；瑝K通過鋁合金連接件與圓柱模型剛性連接。直線氣浮導(dǎo)軌由空氣壓縮機(jī)供氣，氣壓波動(dòng)范圍0.55～0.58 MPa?；瑝K在導(dǎo)軌上的行程為305 mm，滑塊和連接件的質(zhì)量為764.8 g。通過自由衰減試驗(yàn)，測(cè)定振動(dòng)系統(tǒng)在空氣中的阻尼比平均值為5.08×10-3，標(biāo)準(zhǔn)差為2.0×10-4。用兩根線性彈簧連接滑塊與導(dǎo)軌兩端，兩根彈簧的彈性系數(shù)k=27.14 N/m。為減少圓柱末端三維繞流對(duì)振動(dòng)產(chǎn)生的干擾，遵循Fox等[20]的建議，在圓柱模型下方安裝圓矩形端板，并通過兩個(gè)流線型立板與水槽底固定。端板由一個(gè)半圓形薄板(直徑20 cm)和一個(gè)正方形薄板(邊長(zhǎng)20 cm)拼接而成，板厚2 mm，邊緣磨尖。端板上表面距離水槽底部5 cm，以消除水槽底部邊界層的影響。圓柱模型放置于半圓形薄板圓心處，與端板之間的間隙小于2 mm。

采用激光測(cè)距傳感器(Baumer OADM-2016480/S14F)及配套測(cè)量與控制系統(tǒng)(NI USB-6008)對(duì)圓柱位移進(jìn)行測(cè)量。傳感器量程為100～600 mm，分辨率為0.015 mm，線性誤差±0.05 mm。傳感器安裝于導(dǎo)軌一側(cè)，試驗(yàn)坐標(biāo)系原點(diǎn)位于傳感器光源處，x正向?yàn)轫樍飨颍瑈正向?yàn)闄M流向(激光射出方向)，z正向?yàn)樨Q直向上方向。試驗(yàn)中，圓柱僅沿橫流向(y向)振動(dòng)。每組試驗(yàn)以1 000 Hz頻率采樣120 s，并經(jīng)20 Hz低通濾波，得到圓柱振動(dòng)響應(yīng)。

1.3 試驗(yàn)系統(tǒng)可靠性驗(yàn)證

采用彈性支撐均勻圓柱單自由度渦激振動(dòng)對(duì)試驗(yàn)系統(tǒng)的可靠性進(jìn)行驗(yàn)證。驗(yàn)證試驗(yàn)中，圓柱模型為長(zhǎng)436.8 mm(其中淹沒段長(zhǎng)度l=272.8 mm)，直徑d=40 mm。試驗(yàn)設(shè)備與1.2節(jié)所述相同。試驗(yàn)采集了不同流速下的圓柱位移時(shí)程，得到了無(wú)量綱振幅A*隨折合速度Ur的變化趨勢(shì)(A*=A10/d，A10為前10%最大振幅；Ur=U/fn,wd，U為來(lái)流流速，fn,w為圓柱在水中的自然頻率)。將所得結(jié)果與Khalak等[21]、Zhao等[22]的試驗(yàn)結(jié)果進(jìn)行比較，無(wú)量綱參數(shù)見表1，其中：m*為質(zhì)量比；ζa為結(jié)構(gòu)阻尼比；m*ζa為質(zhì)量-阻尼聯(lián)合參數(shù)；A*為無(wú)量綱振幅。如圖4所示，本文結(jié)果呈現(xiàn)明顯的三分支特性，最大振幅在上端分支取得，最大振幅和鎖定區(qū)間范圍與Khalak等、 Zhao等的結(jié)果吻合良好，證實(shí)了試驗(yàn)系統(tǒng)具有良好的可靠性。

圖4 無(wú)量綱振幅A*隨折合流速Ur的變化Fig.4 Variation of the dimensionless amplitude A* with the reduced velocity Ur

表1 彈性支撐均勻圓柱渦激振動(dòng)試驗(yàn)無(wú)量綱參數(shù)Tab.1 Dimensionless parameters in the VIV test of a uniform elastically-supported cylinder

2 試驗(yàn)工況

采用控制變量法，以大直徑段覆蓋率R(R=L/l×100%，L為大直徑段長(zhǎng)度，l為圓柱淹沒長(zhǎng)度)和直徑比D/d為自變量，設(shè)計(jì)試驗(yàn)工況如下。

(1)覆蓋率：試驗(yàn)組1～6，R分別為0，12.5%，25%，50%，75%，100%。

(2)直徑比：試驗(yàn)組4，7，8，D/d分別為1.25，1.5，2。

依據(jù)工況要求選擇外接空心圓管的直徑和長(zhǎng)度，各試驗(yàn)組階梯圓柱，如圖5所示。試驗(yàn)參數(shù)如表2所示。其中：fn,a和fn,w分別為階梯圓柱在空氣和水中的自然頻率；m為滑塊和階梯圓柱的總質(zhì)量。由于空心圓管的密度很小，不同工況下階梯圓柱的質(zhì)量、阻尼比和固有頻率變化較小，其對(duì)階梯圓柱渦激振動(dòng)的影響可忽略不計(jì)。然而，由于不同工況下階梯圓柱的體積有較大變化，因此質(zhì)量比的變化較明顯(m*=6.5～15.2)，這不可避免會(huì)影響振動(dòng)響應(yīng)的幅值，但影響程度較低。Williamson等指出：隨著質(zhì)量比的增大，圓柱的振幅逐漸減小，鎖定區(qū)間收窄，并且振幅減小和區(qū)間收窄的幅度在低質(zhì)量比條件下明顯，而在中等質(zhì)量比條件(如本文的質(zhì)量比條件)下較不明顯。

圖5 各工況階梯圓柱示意圖(深色為外接空心圓管安裝位置)Fig.5 Stepped cylinders in each case (the deep rectangular indicates the installation location of the external pipe)

表2 彈性支撐階梯圓柱渦激振動(dòng)試驗(yàn)參數(shù)Tab.2 Parameters in the VIV test of a uniform elastically-supported cylinder

3 結(jié)果分析

本文對(duì)階梯圓柱渦激振動(dòng)響應(yīng)的無(wú)量綱化處理中統(tǒng)一采用等效直徑D*(D*=RD+(1-R)d)作為參考長(zhǎng)度，無(wú)量綱振幅定義為A*=A10/D*，無(wú)量綱頻率定義為f*=f/fn,w(f為階梯圓柱振動(dòng)頻率)，折合流速定義為Ur=U/fn,wD*(以大、小圓柱直徑為參考長(zhǎng)度，定義Ur,D=U/fn,wD，Ur,d=U/fn,wd)。對(duì)于每個(gè)工況，來(lái)流流速范圍為0.04～0.35 m/s，雷諾數(shù)ReD*范圍為450～6 000(ReD*=UD*/ν，ν為運(yùn)動(dòng)黏滯系數(shù))。

根據(jù)圓柱的動(dòng)力平衡方程

(1)

(2)

式中，ρ為水的密度。

3.1 覆蓋率對(duì)振動(dòng)響應(yīng)的影響

工況1～6中，大、小直徑段的直徑分別為D=24 mm和d=16 mm，覆蓋率R分別為0%，12.5%，25%，50%，75%，100%，其中0%和100%工況可分別視作直徑為16 mm與24 mm的均勻圓柱。圖6給出了各工況下階梯圓柱無(wú)量綱振幅與頻率隨折合流速的變化。

圖6 不同覆蓋率R條件下無(wú)量綱振幅A*和無(wú)量綱頻率f*隨折合流速Ur的變化Fig.6 Variation of the dimensionless amplitude A* and the dimensionless frequency f* with the reduced velocity Ur under different coverage rate R

根據(jù)已有研究結(jié)果，高m*ζa條件下，均勻圓柱的渦激振動(dòng)具有兩個(gè)響應(yīng)分支：初始分支(initial branch)和下端分支(lower branch)，無(wú)量綱最大振幅約為0.6，出現(xiàn)在初始分支上[23]。低m*ζa條件下，均勻圓柱的渦激振動(dòng)具有3個(gè)響應(yīng)分支：初始分支、上端分支(upper branch)和下端分支，無(wú)量綱振幅約為1.0，在上端分支取得。從本文的結(jié)果來(lái)看，R=0%的工況，質(zhì)量比較大，符合高m*ζa渦激振動(dòng)的特點(diǎn)；而R=100%的工況，質(zhì)量比較小，比較符合低m*ζa渦激振動(dòng)的特點(diǎn)。需要指出的是，與Khalak等的試驗(yàn)(m*ζa=1.08×10-2)相比，R=100%工況的m*ζa(m*ζa=3.71×10-2)仍較大，導(dǎo)致上端分支較窄，最大無(wú)量綱振幅略小。

對(duì)于階梯圓柱，當(dāng)覆蓋率較低時(shí)(R=12.5%)，振動(dòng)響應(yīng)與高m*ζa均勻圓柱的振動(dòng)響應(yīng)類似，但下端分支的振幅較大。隨著覆蓋率的增加(R=25%～50%)，在初始分支和下端分支之間出現(xiàn)過渡分支(transition branch)。在過渡分支上，階梯圓柱大、小直徑段的脫渦頻率鎖定在振動(dòng)頻率上，振動(dòng)頻率與固有頻率相近，無(wú)量綱振幅隨著折合流速的增大而增大，并在過渡分支末端取得最大振幅。初始分支和過渡分支具有不同的振幅增長(zhǎng)率，以R=25%工況較為明顯。對(duì)比R=25%和50%工況發(fā)現(xiàn)，R=50%時(shí)，初始分支較窄，過渡分支的振幅增長(zhǎng)更快，初始分支和過渡分支之間的界限更不明顯。此外，均處于鎖定狀態(tài)的過渡分支和下端分支的振動(dòng)頻率不同，下端分支的振動(dòng)頻率略大，這與在下端分支上大直徑段的脫渦頻率與振動(dòng)頻率鎖定有關(guān)。隨著覆蓋率進(jìn)一步增大，初始分支進(jìn)一步收窄，直至消失(R=75%)，過渡分支占據(jù)整個(gè)振幅上升段。

結(jié)構(gòu)振動(dòng)頻率在初始、過渡和下端分支(鎖定區(qū)間)上均呈現(xiàn)單頻特征，且與結(jié)構(gòu)的固有頻率相近，因而振動(dòng)響應(yīng)較為規(guī)律，振幅較大。在非鎖定區(qū)間，結(jié)構(gòu)的振動(dòng)頻率存在兩個(gè)主導(dǎo)頻率，其中的一個(gè)頻率與結(jié)構(gòu)的固有頻率相近，另外一個(gè)頻率與階梯圓柱的脫渦頻率有關(guān)，隨折合流速的增大而近似線性增大。這樣的結(jié)果與Khalak等在均勻圓柱渦激振動(dòng)中得到的結(jié)果一致。然而，需要說明的是，對(duì)于階梯圓柱工況(R=12.5%～75.0%)，與脫渦相關(guān)的無(wú)量綱振動(dòng)頻率f*(f*=fD*/U)低于均勻圓柱的無(wú)量綱脫渦頻率St(=fsD/U≈0.2)，這與階梯圓柱后方復(fù)雜的三維流場(chǎng)結(jié)構(gòu)引起的較低的脫渦頻率有關(guān)。當(dāng)R=0%和100%時(shí)(均勻圓柱工況)，與脫渦相關(guān)的無(wú)量綱振動(dòng)頻率f*與均勻圓柱的無(wú)量綱脫渦頻率St吻合良好。

進(jìn)一步對(duì)鎖定區(qū)間內(nèi)的無(wú)量綱振動(dòng)頻率隨折合流速的變化進(jìn)行分析，如圖7所示。當(dāng)R=0%和100%時(shí)(均勻圓柱工況)，振動(dòng)頻率具有一個(gè)鎖定區(qū)間，并且R=100%的振動(dòng)頻率高于R=0%的振動(dòng)頻率。這是因?yàn)?，?dāng)R=100%時(shí)，圓柱的質(zhì)量比較低，附加質(zhì)量效應(yīng)明顯，因此鎖定時(shí)振動(dòng)頻率較高。這與Khalak和陳威霖等[24-25]的結(jié)論一致。然而，當(dāng)R=12.5%～75.0%時(shí)(階梯圓柱工況)，振動(dòng)頻率出現(xiàn)兩個(gè)鎖定區(qū)間，隨折合流速的增大，依次為第一鎖定區(qū)間和第二鎖定區(qū)間。第二鎖定區(qū)間的振動(dòng)頻率高于第一鎖定區(qū)間的，兩者之間存在一個(gè)過渡區(qū)域。第一鎖定區(qū)間的出現(xiàn)伴隨著過渡分支的出現(xiàn)，而第二鎖定區(qū)間與下端分支對(duì)應(yīng)。本文稱此現(xiàn)象為階梯圓柱渦激振動(dòng)的“雙鎖定”(dual lock-in)現(xiàn)象。需要說明的是，R=12.5%工況的第一鎖定區(qū)間很短(Ur=5.8～6.3)，并未引起可觀察到的過渡分支(見圖6b)。

圖7 不同覆蓋率R條件下鎖定區(qū)間內(nèi)無(wú)量綱頻率 f*隨折合流速Ur的變化Fig.7 Variation of the dimensionless frequency f* with the reduced velocity Ur in the lock-in region under different coverage rates R

引起“雙鎖定”現(xiàn)象的原因與階梯圓柱小直徑段和大直徑段的脫渦頻率隨折合流速增大依次與圓柱的振動(dòng)頻率鎖定有關(guān)。眾所周知，在TrSL區(qū)間(Re=103～105)內(nèi)，圓柱的無(wú)量綱脫渦頻率St(St=fsD/U≈0.2)近乎不變，在相同的流速下，圓柱的脫渦頻率與圓柱的直徑呈反比，直徑越小，脫渦頻率越高，反之亦然。對(duì)于階梯圓柱，隨著來(lái)流速度的增加，小直徑段的泄渦頻率首先與圓柱的固有頻率接近，引起圓柱振幅增大，振動(dòng)響應(yīng)進(jìn)入初始分支。此時(shí)，由于大直徑段的脫渦頻率較低，與固有頻率相差較大，因此大直徑段的脫渦不引起明顯振動(dòng)。隨著流速的增大，振動(dòng)響應(yīng)進(jìn)入小直徑段的下端分支，此時(shí)小直徑段的泄渦鎖定在振動(dòng)頻率上。與此同時(shí)，大直徑段的脫渦頻率與圓柱的固有頻率接近，引發(fā)大直徑段的脫渦與振動(dòng)之間的鎖定。大小直徑段的鎖定同時(shí)發(fā)生，導(dǎo)致階梯圓柱的振幅進(jìn)一步增大。進(jìn)一步增大流速，小直徑段的脫渦與振動(dòng)脫離鎖定，而大直徑段的振動(dòng)進(jìn)入下端分支，脫渦與振動(dòng)保持鎖定。此時(shí)階梯圓柱的振幅開始減小，振動(dòng)進(jìn)入下端分支。再進(jìn)一步增大流速，大小直徑段的脫渦與振動(dòng)均不鎖定，階梯圓柱的振幅很小，振動(dòng)響應(yīng)進(jìn)入非鎖定區(qū)間。以上過程可結(jié)合圖8進(jìn)行理解。

圖8 階梯圓柱渦激振動(dòng)分支隨流速U的變化Fig.8 Variation of VIV branches of step cylinder with velocity U

除此以外，我們還發(fā)現(xiàn)：對(duì)于階梯圓柱，覆蓋率R越高，第一鎖定頻率越低，第二鎖定頻率越高。這可歸因于大直徑段脫渦和渦激振動(dòng)的影響;在第一鎖定區(qū)間內(nèi)，由于大直徑段的脫渦頻率較低，R越高時(shí)，大直徑段越長(zhǎng)，對(duì)階梯圓柱振動(dòng)頻率的拉低效應(yīng)越明顯;在第二鎖定區(qū)間，小直徑段處于非鎖定狀態(tài)，而大直徑段鎖定。大直徑段越長(zhǎng)，階梯圓柱的質(zhì)量比越低，附加質(zhì)量效應(yīng)越明顯，無(wú)量綱的鎖定頻率越高。此外，在第二鎖定區(qū)間內(nèi)，隨著覆蓋率的增加，階梯圓柱的振動(dòng)頻率逐漸由R=0%工況(小直徑均勻圓柱)向R=100%工況(大直徑均勻圓柱)過渡。

3.2 直徑比的影響

Vandiver等的研究表明，當(dāng)覆蓋率R較高時(shí)，大直徑段主導(dǎo)階梯圓柱的振動(dòng)響應(yīng)；反之，當(dāng)R較低時(shí)，小直徑段起主導(dǎo)作用。當(dāng)R=50%時(shí)，大直徑段和小直徑段存在復(fù)雜的競(jìng)爭(zhēng)關(guān)系。為了充分研究直徑比對(duì)階梯圓柱渦激振動(dòng)的影響，本節(jié)選取R=50%的工況進(jìn)行分析。圖9給出了覆蓋率R=50%不同直徑比下(D/d=1.25，1.5，2)無(wú)量綱振幅和頻率隨折合流速的變化。對(duì)比可知，直徑比對(duì)響應(yīng)分支的影響主要體現(xiàn)在初始分支與過渡分支上，直徑比D/d越大，初始分支越短，過渡分支越長(zhǎng)。

圖9 不同直徑比D/d條件下無(wú)量綱振幅A*和 無(wú)量綱頻率f*隨折合流速Ur的變化Fig.9 Variation of the dimensionless amplitude A* and the dimensionless frequency f* with the reduced velocity Ur at different diameter ratios D/d

圖10所示為鎖定區(qū)間內(nèi)不同直徑比階梯圓柱無(wú)量綱振動(dòng)頻率隨折合流速的變化。對(duì)比可知，不同直徑比工況下，階梯圓柱的振動(dòng)響應(yīng)均呈現(xiàn)“雙鎖定”特性，并且隨著直徑比的增大，第一鎖定頻率小幅降低，而第二鎖定頻率則明顯增加。其原因與覆蓋率對(duì)鎖定頻率的影響類似，此處不再贅述。

圖10 不同直徑比D/d條件下鎖定區(qū)間內(nèi) 無(wú)量綱頻率f*隨折合流速Ur的變化Fig.10 Variation of the dimensionless frequency f* with reduced velocity Ur in the lock-in region at different diameter ratios D/d

3.3 階梯圓柱渦激振動(dòng)的受力特性

根據(jù)式(1)和式(2)，可計(jì)算階梯圓柱在渦激振動(dòng)中所受升力的大小。本節(jié)以試驗(yàn)組4和6為例，說明升力系數(shù)均方根CL,rms、無(wú)量綱振幅A*和無(wú)量綱頻率f*隨折合流速的變化，如圖11所示。需要說明的是：試驗(yàn)組4的圓柱為D=24 mm，d=16 mm，R=50%的階梯圓柱，試驗(yàn)組6的圓柱為D=24 mm的均勻圓柱。

圖11 升力系數(shù)CL,rms、無(wú)量綱振幅A*和 無(wú)量綱頻率f*隨折合速度Ur的變化Fig.11 Variation of the lift coefficient CL,rms, the dimensionless amplitude A* and the dimensionless frequency f* with reduced velocity Ur

如圖11(a)所示，均勻圓柱的CL,rms在初始分支迅速增加，在上端分支快速減小，最大值在初始分支和上端分支間之間取得。在下端分支，CL,rms較小，隨著折合流速的增大而緩慢減小。上述與Khalak等的結(jié)論一致。

對(duì)于階梯圓柱來(lái)說，與振動(dòng)響應(yīng)相對(duì)應(yīng)，升力也呈現(xiàn)出類似的分支特性。如圖11(b)所示，階梯圓柱的CL,rms在初始分支快速增加；隨著過渡分支的出現(xiàn)，CL,rms出現(xiàn)增長(zhǎng)拐點(diǎn)(a點(diǎn))，增長(zhǎng)的斜率顯著降低。在過渡分支上，CL,rms先緩增后快減，最大值在Ur=5處取得(b點(diǎn))。在過渡分支和下端分支的交接處(c點(diǎn))，CL,rms出現(xiàn)明顯轉(zhuǎn)折。在下端分支上，CL,rms緩慢降低，值較小。對(duì)比可知，階梯圓柱所受升力隨折合流速的變化規(guī)律在初始分支和下端分支上與均勻圓柱的相似。雖然階梯圓柱的過渡分支和均勻圓柱的上端分支均處于初始分支和下端分支之間，但兩者具有明顯的差別。在過渡分支上，階梯圓柱的小直徑段處于下端分支，CL,rms隨折合流速快速降低，而大直徑段處于初始分支，CL,rms隨折合流速快速增加。兩種趨勢(shì)相互作用，造成了CL,rms在過渡分支上先增后減的變化趨勢(shì)。

4 結(jié) 論

本研究在天津大學(xué)流體力學(xué)試驗(yàn)室的低湍流度水槽中開展了均勻來(lái)流中彈性支撐階梯圓柱的渦激振動(dòng)試驗(yàn)，通過改變階梯圓柱的覆蓋率和直徑比，討論了其對(duì)振動(dòng)響應(yīng)、頻率和升力的影響，得到的研究結(jié)論總結(jié)如下：

(1)覆蓋率R和直徑比D/d對(duì)階梯圓柱渦激振動(dòng)特性影響顯著。當(dāng)R=25%～75%時(shí)，階梯圓柱的振動(dòng)響應(yīng)出現(xiàn) “過渡分支”。在過渡分支上，階梯圓柱的大、小直徑段均處于鎖定狀態(tài)，但小直徑段處于下端分支，而大直徑段處于初始和上端分支。在過渡分支上，階梯圓柱的振幅隨折合流速的增大而增大，階梯圓柱所受升力隨折合流速的增大先增后減。隨著覆蓋率R和直徑比D/d增加，過渡分支逐漸變寬，與此同時(shí)，初始分支逐漸收窄。

(2)發(fā)現(xiàn)了階梯圓柱渦激振動(dòng)“雙鎖定”現(xiàn)象。在鎖定區(qū)間內(nèi)，階梯圓柱的振動(dòng)頻率隨著折合流速的增加先后鎖定在兩個(gè)頻率上，且第一鎖定頻率略低于第二鎖定頻率。第一鎖定區(qū)間對(duì)應(yīng)過渡分支，而第二鎖定區(qū)間對(duì)應(yīng)下端分支。隨著覆蓋率R和直徑比D/d增加，第一鎖定頻率小幅降低，而第二鎖定頻率明顯增加。

本研究?jī)H對(duì)階梯圓柱渦激振動(dòng)的振動(dòng)響應(yīng)和受力特性進(jìn)行了研究，發(fā)現(xiàn)了區(qū)別于均勻圓柱渦激振動(dòng)的振動(dòng)特性，并給出了初步的機(jī)理解釋。然而，更深入地認(rèn)識(shí)階梯圓柱渦激振動(dòng)的物理規(guī)律需要更進(jìn)一步開展流動(dòng)顯示模型試驗(yàn)和數(shù)值模擬，對(duì)階梯圓柱的尾渦結(jié)構(gòu)以及尾流與柱體振動(dòng)之間的相互作用進(jìn)行研究。