楊沖霄，袁昱超，薛鴻祥，唐文勇

(上海交通大學(xué) 海洋工程國家重點實驗室，上海 200240)

隨著我國海洋油氣資源開發(fā)逐漸向深海邁進(jìn)，海洋工程結(jié)構(gòu)物的安全性能越來越得到重視。海洋立管系統(tǒng)作為連接水面浮式裝置和海底設(shè)備的導(dǎo)管，在復(fù)雜的載荷作用下，極易出現(xiàn)碰撞、振動及疲勞等破壞形式，是深海工程裝備中最薄弱的環(huán)節(jié)[1]。在來流作用下，立管的交替泄渦現(xiàn)象產(chǎn)生振蕩的流體力，引起結(jié)構(gòu)的渦激振動[2]。早期國外學(xué)者對該問題開展了大量的試驗研究[3-5]。其中Gopalkrishnan[3]對Re為10 000下的圓柱受迫振動進(jìn)行了模型試驗，并首次構(gòu)建了受迫振動流體力系數(shù)庫。在數(shù)值方法中，隨著計算機(jī)技術(shù)的發(fā)展，計算流體力學(xué)方法(computational fluid dynamics， 簡稱CFD)應(yīng)運(yùn)而生，為了提高計算效率，將立管簡化為二維圓柱是較為常見的研究手段。Muhamad和Krish[6]使用RANS k-ω湍流模型對Re=10 000來流下的圓柱進(jìn)行單自由度受迫振動數(shù)值模擬，結(jié)果表明漩渦尾跡對圓柱振動頻率較為敏感，當(dāng)振動頻率接近Strouhal頻率時，會發(fā)生漩渦的相位變換。王亞非[7]數(shù)值模擬了雙自由度彈性支撐圓柱體的渦激振動問題，結(jié)果表明振蕩流擴(kuò)大了自激振動的鎖定范圍，并且使得順流向振動幅值大大增加，甚至超過了橫向所能激發(fā)的最大幅值。付博文等[8]基于切片理論，通過使用徑向基函數(shù)法作為OpenFOAM中的動網(wǎng)格策略，模擬了長細(xì)比為1 000的柔性立管在橫流向和順流向的振動，數(shù)值模擬重現(xiàn)了高階主控模態(tài)及主控模態(tài)的頻繁變換等大長細(xì)比柔性立管的渦激振動特性。王凱鵬[9]通過對均勻來流和剪切來流兩種來流形式下的靜止圓柱、橫向受迫振動圓柱繞流及圓柱渦激振動問題展開數(shù)值模擬，首次系統(tǒng)地分析了來流形式對圓柱繞流和圓柱渦激振動問題的影響，并指出來流形式是很重要的一個影響因素。總體來說目前的研究主要集中在定常流。

實際生產(chǎn)過程中，海洋浮式結(jié)構(gòu)物在風(fēng)、浪、流的聯(lián)合作用下會帶動立管在水中往復(fù)運(yùn)動，此時立管的遭遇流場可等效為振蕩流場[10]。相比定常流，振蕩流場引起的圓柱繞流問題更為復(fù)雜，需要更加深入的研究，目前已有學(xué)者獲得了相關(guān)成果。Zhao等[11]模擬了振蕩流和均勻流共同作用下圓柱的渦激振動，探討了流動比例a對圓柱的響應(yīng)，并發(fā)現(xiàn)了在a=0.8，Vr=7時，漩渦在一個振蕩流周期中經(jīng)歷2S、2P和2T三種泄渦模式。鄧躍[12]對低雷諾數(shù)時單自由度彈性支撐圓柱在均勻流和振蕩流共同作用下的受迫和自激振動進(jìn)行了數(shù)值模擬，結(jié)果發(fā)現(xiàn)當(dāng)有振蕩流參與時，結(jié)構(gòu)的鎖定范圍和振動幅值等都有明顯的變化。鄧迪等[13]采用OpenFOAM開源軟件對在靜水中做周期性振蕩運(yùn)動的二維剛性圓柱渦激振動進(jìn)行數(shù)值模擬，發(fā)現(xiàn)圓柱的橫向振動加劇了升力系數(shù)的變化，使得泄渦方向和圓柱表面漩渦分離點的位置發(fā)生了明顯變化。

目前對于振蕩流的研究主要集中在低雷諾數(shù)或小KC數(shù)工況，而已有認(rèn)知表明，不同雷諾數(shù)及KC數(shù)對圓柱繞流特性影響顯著[14-16]。文中重點研究高雷諾數(shù)條件并兼顧大KC數(shù)工況。首先，對均勻流下雷諾數(shù)10 000時的圓柱受迫振動進(jìn)行了數(shù)值模擬，通過與試驗結(jié)果進(jìn)行對比，驗證了基于CFD方法復(fù)現(xiàn)受迫振動試驗并研究流體力系數(shù)的可行性。之后開展振蕩流下圓柱受迫振動數(shù)值仿真，分析不同KC數(shù)下圓柱的水動力特性及漩渦形態(tài)，并歸納出大KC數(shù)和小KC數(shù)下升力系數(shù)和漩渦發(fā)放特點。

1 數(shù)值分析方法

1.1 控制方程

考慮不可壓縮流場，采用雷諾時均方法(Reynolds Averaged Navier-Stokes， 簡稱RANS)，結(jié)合剪切應(yīng)力輸運(yùn)模型(k-ω SST)求解N-S方程，控制方程：

(1)

(2)

k-ω SST湍流模型的輸運(yùn)方程：

(3)

(4)

式中：μt為渦黏性；Sij為平均速度應(yīng)變率張量；τtij表示雷諾應(yīng)力的渦黏性模型；σk、β*、σω、σω2均為經(jīng)驗參數(shù)；Pω代表交錯擴(kuò)散項。F1為混合函數(shù)，在近壁處采用Wilcox k-ω模型，邊界層邊緣和自由剪切層采用k-ε模型，中間通過F1函數(shù)實現(xiàn)過渡。

1.2 水動力系數(shù)識別方法

在Gopalkrishnan[3]的試驗中，圓柱進(jìn)行受迫振動，位移函數(shù)：

y(t)=Asin(2πf0t)

(5)

式中：A為振動幅值；f0為振動頻率；t為運(yùn)動時間。

由此得到圓柱運(yùn)動的速度函數(shù)為：

v(t)=2πf0Acos(2πf0t)

(6)

當(dāng)圓柱以式(5)進(jìn)行振動時，垂直來流方向的升力可表示為:

L=L0sin(2πf0t+φ0)+Lssin(2πfst+φs)

(7)

式中：L0、Ls為升力幅值；φ0、φs為相位角；fs代表Strouhal頻率，其中下標(biāo)為0表示與圓柱振動頻率相關(guān)，下標(biāo)為s表示與Strouhal頻率相關(guān)。

當(dāng)圓柱處于鎖定狀態(tài)時，Strouhal頻率消失，并且圓柱并未以頻率fs在振動，即該成分不參與到流體與結(jié)構(gòu)的能量傳遞中，故在實際分析時，升力得到簡化：

L=L0sin(2πf0t+φ0)

(8)

升力幅值及相位角可通過對時歷曲線進(jìn)行傅里葉擬合得到：

(9)

計算可得參數(shù)：

(10)

(11)

(12)

升力幅值和相位角可表示為：

(13)

(14)

升力系數(shù)由下式計算得到：

(15)

式中：l為圓柱高度，在二維計算中取單位高度；D為圓柱直徑；U為遭遇流速。

升力系數(shù)還可進(jìn)一步分解得到激勵力系數(shù)和慣性力系數(shù)。升力系數(shù)中和圓柱運(yùn)動速度同相位的部分定義為激勵力系數(shù)：

CL_V0=CL0sinφ0

(16)

激勵力系數(shù)為正代表能量從流體輸入到圓柱結(jié)構(gòu)中，此時易發(fā)生圓柱的渦激振動，激勵力系數(shù)為負(fù)代表能量從結(jié)構(gòu)輸出到流體中。

升力系數(shù)中和圓柱運(yùn)動加速度同相位的部分定義為慣性力系數(shù)：

CL_A0=CL0(-cosφ0)

(17)

無因次振幅定義為：

A*=A/D

(18)

無因次頻率定義為：

f*=f0U/D

(19)

振蕩流定義如下：

u(t)=2πfBsin(2πft)

(20)

式中：B為振蕩流幅值；f為振蕩流頻率。

在均勻流中，采用雷諾數(shù)定義流體特征：

Re=UD/ν

(21)

振蕩流中，引入柯萊根—卡彭特數(shù)(Keulegan-Carpenter number，KC)表征來流：

KC=Umax/fD

(22)

Umax=2πfB

(23)

式中：Umax為振蕩流最大速度。

1.3 計算網(wǎng)格和邊界條件

計算域如圖1所示。X軸平行來流方向，Y軸垂直來流方向。圓柱直徑D=0.025 4 m，上下邊距圓柱中心8D，水平方向50D。由于圓柱周圍流場變化較為劇烈，在圓柱中心3.5D范圍內(nèi)進(jìn)行加密，如圖2所示。左端為速度入口，均勻流流速0.4 m/s；上下兩端為對稱邊界；右端為自由出流；圓柱采用無滑移固壁條件；壓力—速度耦合采用SIMPLE算法；時間項采用二階隱式積分方法；對流項采用二階迎風(fēng)離散格式。

圖1 計算域Fig. 1 Computational domain

圖2 局部加密Fig. 2 Partial encryption

2 均勻流工況試驗驗證

2.1 網(wǎng)格收斂性驗證

首先對定常流速0.4 m/s、Re=10 000下的靜止圓柱進(jìn)行了數(shù)值模擬，給出了普通網(wǎng)格、加密網(wǎng)格及其他學(xué)者計算結(jié)果的對比，如表1所示。

表1 兩套網(wǎng)格參數(shù)Tab. 1 Two sets of grid parameters

根據(jù)對比可知，兩套網(wǎng)格的計算結(jié)果均與試驗數(shù)據(jù)對應(yīng)較好，在綜合考慮計算精度和計算成本的前提下，選擇普通網(wǎng)格完成振蕩流流動計算。

2.2 流體力系數(shù)對比

計算了均勻流振幅比為0.3時各無因次頻率下的流體力系數(shù)，并與Gopalkrishnan[3]的試驗數(shù)據(jù)進(jìn)行了對比，升力系數(shù)、激勵力系數(shù)、慣性力系數(shù)及相位角結(jié)果如圖3所示。

圖3 流體力系數(shù)對比Fig. 3 Comparison of fluid force coefficient

由圖3可知，隨著無因次頻率的增加，升力系數(shù)幅值呈現(xiàn)逐漸增加的趨勢，文中模擬結(jié)果在趨勢和數(shù)值上和試驗結(jié)果吻合較好。激勵力系數(shù)變化較為復(fù)雜，在無因次頻率0.17處出現(xiàn)正峰值，在0.21處出現(xiàn)負(fù)峰值，模擬結(jié)果基本可復(fù)現(xiàn)這一現(xiàn)象，但數(shù)值上略有差距。慣性力系數(shù)變化與無因次頻率呈現(xiàn)負(fù)相關(guān)關(guān)系，文中CFD模擬結(jié)果較好的吻合了試驗數(shù)據(jù)。相位角在低頻率處為負(fù)值，在頻率0.15附近發(fā)生相位突變現(xiàn)象，之后隨著無因次頻率增加逐漸趨于0，CFD模擬結(jié)果也基本符合這些特征。但在低頻率處，文中模擬結(jié)果與試驗存在一定誤差，這可能是因為Gopalkrishnan所進(jìn)行的試驗結(jié)果包含了復(fù)雜的三維效應(yīng)，故二維模型暫未全面反映三維試驗結(jié)果。文獻(xiàn)[18]針對二維數(shù)值模擬與試驗誤差也提出了相似的分析。

綜上所述，文中所采用的網(wǎng)格劃分和CFD設(shè)置方法復(fù)現(xiàn)了圓柱受迫振動的試驗結(jié)果，在水動力系數(shù)模擬上具有較高的準(zhǔn)確性，可運(yùn)用到振蕩流中開展研究。

2.3 均勻流漩渦形態(tài)分析

均勻流工況漩渦發(fā)放形態(tài)如圖4所示，選取一個脫落周期的漩渦發(fā)放結(jié)果。圖4(a)時，圓柱向上方運(yùn)動，尾渦在圓柱下流向產(chǎn)生，并逐漸向后延伸。當(dāng)圓柱向下運(yùn)動到圖4(b)位置時，漩渦從圓柱的下尾渦末端脫落，此時上尾渦的末端也分離出了即將脫落的漩渦。圖4(c)時，漩渦從上尾渦末端分離。脫落周期持續(xù)0.37 s，與圓柱此時的振動周期是吻合的?？梢钥吹剑鰷u在一個周期內(nèi)呈現(xiàn)上下交替脫落現(xiàn)象，為典型的2S脫落模式。

圖4 均勻流下圓柱漩渦脫落(A*=0.5，f*=0.17)Fig. 4 Vortex shedding under uniform flow(A*=0.5，f*=0.17)

3 振蕩流結(jié)果與分析

3.1 工況設(shè)置

為了與前文均勻流形成對比，振蕩流流速幅值定為0.4 m/s，振蕩流可表示為：

v=0.4×sin(2πft)

(24)

基于第2節(jié)的研究基礎(chǔ)，通過改變振蕩流振蕩周期，研究不同KC數(shù)下的二維圓柱受迫振動特性。其中，振蕩流工況設(shè)置情況如表2所示。為研究圓柱振動幅度和振動頻率對計算結(jié)果的影響，同一振幅比下選擇8組無因次頻率進(jìn)行計算，無因次振幅比設(shè)置6組，共計算240個振蕩流工況。

表2 振蕩流工況設(shè)置Tab. 2 Setting of oscillation flow condition

3.2 流體力系數(shù)分析

提取0.50D和1.00D兩組振幅比的計算結(jié)果進(jìn)行分析，升力系數(shù)幅值變化情況見圖5(a)。在低振動頻率時，升力系數(shù)呈現(xiàn)較小的值，隨著振動頻率增加，升力系數(shù)幅值也逐漸增加。在該工況下，不同KC數(shù)對升力系數(shù)幅值影響不大。激勵力系數(shù)變化情況見圖5(b)。激勵力系數(shù)變化較為復(fù)雜，且對KC數(shù)的敏感性較高。此時不同KC數(shù)下的激勵力系數(shù)隨無因次頻率變化基本呈現(xiàn)先增加后減小的趨勢，峰值區(qū)出現(xiàn)在0.18附近，相比均勻流的0.17略有變化，符合文獻(xiàn)[12]的研究結(jié)論。在振動頻率0.15～0.20之間，KC數(shù)對激勵力系數(shù)影響較大，該區(qū)間也屬于均勻流試驗測得的鎖定發(fā)生的關(guān)鍵區(qū)域。由圖5(b)可知，小KC數(shù)下，激勵力系數(shù)變化較大。除在峰值區(qū)，KC=31.5時激勵力系數(shù)基本保持在較低水平。這可能是因為此工況下流速及流向變化較快，漩渦尚未脫落便遭遇反向流速，無法形成穩(wěn)定的漩渦脫落周期。KC數(shù)增大后，激勵力系數(shù)也逐漸增加，但隨著KC數(shù)變大，流態(tài)逐漸趨于定常流，此時激勵力系數(shù)值將保持相近。

圖5 低振幅比(0.50D)流體力系數(shù)對比Fig. 5 Contrast of fluid force coefficients at low amplitude ratio (0.50D)

圖6給出了振動幅值1.00D時，升力系數(shù)幅值和激勵力系數(shù)隨KC數(shù)的變化情況。隨著振動幅值的增加，KC數(shù)對流體力系數(shù)的影響效應(yīng)并不顯著，這與小振幅工況時是不同的。圖6(b)表明，激勵力系數(shù)隨著無因次頻率的增加逐漸減小。大KC數(shù)工況相對較為穩(wěn)定，小KC數(shù)工況KC=31.5時振蕩流周期變化較快，漩渦脫落相對不穩(wěn)定，規(guī)律性較弱，故在低無因次頻率下呈現(xiàn)出了下降趨勢。但整體來看，各工況下激勵力系數(shù)均小于0，此時流體對結(jié)構(gòu)振動起阻尼作用，能量由結(jié)構(gòu)傳向流體。激勵力系數(shù)由升力系數(shù)幅值和相位角計算得到，由Gopalkrishnan[3]均勻流試驗可知，低振幅比與高振幅比時相位角變化具有明顯不同。振蕩流工況下高低振幅比時相位角變化趨勢也不同，1.00D時相位角均保持為負(fù)值，故激勵力系數(shù)變化趨勢與0.50D存在較大差別。

圖6 高振幅比(1.00D)流體力系數(shù)對比Fig. 6 Contrast of fluid force coefficients at high amplitude ratio (1.00D)

綜上可知，在小振動幅值時，當(dāng)圓柱的無因次頻率處于渦激振動鎖定區(qū)間內(nèi)，KC數(shù)對激勵力系數(shù)影響較大，且大、小KC數(shù)下呈現(xiàn)不同的規(guī)律性。在大振動幅值時，KC數(shù)對該雷諾數(shù)下流體力系數(shù)的影響將逐漸減小，激勵力系數(shù)保持負(fù)值。

3.3 漩渦形態(tài)分析

基于振幅比0.50D計算結(jié)果，分析KC數(shù)對漩渦發(fā)放形態(tài)的影響。前文分析可得，振動頻率在0.15～0.20之間流體力系數(shù)受KC數(shù)變化較為敏感，因此選擇頻率0.17作為分析重點。將5組KC數(shù)分為兩類，其中，小KC數(shù)工況為組1和組2(KC=31.5和KC=63.0)，大KC數(shù)工況為組3、組4和組5(KC=126.0、KC=252.0和KC=503.9)。由于大小兩類工況下不同KC數(shù)的漩渦形態(tài)具有相似性，故大KC數(shù)和小KC數(shù)各選擇一組典型形態(tài)進(jìn)行分析。

3.3.1 大KC數(shù)工況

圖7給出KC=503.9、無因次頻率0.17、一個振蕩流周期內(nèi)漩渦脫落發(fā)展?fàn)顟B(tài)。隨著流速的逐漸增大，尾渦在圓柱壁面上產(chǎn)生，并逐漸向圓柱右側(cè)延伸。當(dāng)?shù)竭_(dá)圖7(a)時，有漩渦從尾渦后方脫落。圖7(b)對應(yīng)流速達(dá)到峰值，在由(a)到(b)的過程中，速度進(jìn)一步增加，圓柱后方的尾渦也逐漸伸長。此階段內(nèi)伸長的尾渦與之前脫落的漩渦連在一起，并未出現(xiàn)顯著漩渦脫離現(xiàn)象。經(jīng)過圖7(b)后，流速逐漸下降，此時圓柱后方的漩渦開始上下交替脫落，呈現(xiàn)出典型的2S泄渦模式，如圖7(c)所示。圓柱運(yùn)動到圖7(d)時，此時圓柱后方依舊有漩渦脫落，但已不再是2S模式，并且漩渦的大小和強(qiáng)度均開始下降。

圖7 大KC數(shù)下漩渦脫落(A*=0.5，f*=0.17，KC=503.9)Fig. 7 Vortex shedding at large KC number (A*=0.5，f*=0.17，KC=503.9)

當(dāng)流速繼續(xù)下降時，脫落的漩渦也逐漸在圓柱后方消散。圖7(e)顯示此時流速很小，圓柱壁面已無法產(chǎn)生尾渦。之后流速反向并逐漸增加。漩渦開始在圓柱的左側(cè)產(chǎn)生，但此時尾渦較短，漩渦在離圓柱較近的位置脫落，如圖7(f)所示。流速進(jìn)一步增加，尾渦逐漸延伸，但無法觀察到脫落的漩渦。這一階段類似圖7(a)至7(b)。圖7(g)時流速剛過最大值，漩渦開始從尾渦上脫落，并呈現(xiàn)出2S泄渦模式，與流速正向不同的是，這一階段2S模式的持續(xù)時間要更小。圖7(h)后，流速逐漸下降，圓柱的尾渦長度慢慢變短，強(qiáng)度逐漸減弱。

3.3.2 小KC數(shù)工況

圖8為KC數(shù)為31.5時的尾渦演化。圖8(a)時，流速剛經(jīng)過最大值，尾渦在圓柱后側(cè)產(chǎn)生延長。隨著流速下降，從圖8(b)可以發(fā)現(xiàn)有單個漩渦從圓柱壁面上脫落，但由于沒有長時間的單側(cè)流向，難以觀察到漩渦一個接一個在圓柱的下流向消散的現(xiàn)象。圖8(c)時，流速反向，脫落的漩渦被反向流速帶回到圓柱附近，導(dǎo)致圓柱周圍漩渦分布復(fù)雜化，這也是小KC數(shù)時的典型狀態(tài)。之后流速增加，尾渦在圓柱的左側(cè)產(chǎn)生。

圖8 小KC數(shù)下漩渦脫落(A*=0.5，f*=0.17，KC=31.5)Fig. 8 Vortex shedding at low KC number (A*=0.5，f*=0.17，KC=31.5)

在f*=0.17下，振蕩流隨著KC數(shù)的增加，漩渦脫落的模式各有不同。小KC數(shù)下流體速度變化較快，漩渦從圓柱壁面脫落后難以保持穩(wěn)定泄渦，當(dāng)流速反向后，脫落的漩渦又被沖向圓柱，使得圓柱周圍壓力變化更加復(fù)雜。在大KC數(shù)時，漩渦以典型的2S模式從圓柱壁面脫落，且隨著速度下降，漩渦強(qiáng)度逐漸下降，脫落的漩渦也可在流速反向前消散完成。

3.4 升力系數(shù)時歷分析

為更好地分析漩渦產(chǎn)生發(fā)展過程和升力系數(shù)之間的關(guān)系，提取了A*=0.5，f*=0.17對應(yīng)的升力系數(shù)時歷曲線，并對照漩渦發(fā)放過程，探討二者之間的內(nèi)在聯(lián)系。

圖9展示了大KC數(shù)時的升力系數(shù)時歷曲線。3種大KC數(shù)工況均觀察到了振幅調(diào)制現(xiàn)象。圖9(c)給出了KC=503.9時升力系數(shù)變化情況。該曲線可以看到兩個升力系數(shù)增大區(qū)，即振蕩流速度幅值附近(72 s和88 s)。當(dāng)運(yùn)動到64 s時，振蕩流速度由0開始逐漸增大，升力系數(shù)也逐漸增加。當(dāng)圓柱運(yùn)動至72 s時，進(jìn)入到2S發(fā)放模式，升力系數(shù)進(jìn)入第一個峰值區(qū)。80 s時，流速降低為0，可以觀察到此時升力系數(shù)要小于峰值區(qū)。流速反向后，速度開始增加，再次進(jìn)入到2S發(fā)放模式，此時升力系數(shù)同樣達(dá)到峰值區(qū)。

圖9 大KC數(shù)升力系數(shù)曲線Fig. 9 Amplitude curve of lift coefficient at large KC number

KC=126.0時，峰值區(qū)出現(xiàn)在18 s和22 s附近，并且可以看到第二個峰值區(qū)的范圍要略大于第一個，且峰值更高。KC=252.0時，峰值區(qū)出現(xiàn)在36 s和44 s，該工況下振幅調(diào)制現(xiàn)象較為對稱。

圖10給出了小KC數(shù)及均勻流工況下升力系數(shù)時歷曲線。當(dāng)KC數(shù)較小時，流體速度變化較快，在一個振蕩流周期內(nèi)無法觀察到明顯的兩次振幅調(diào)制現(xiàn)象，但振幅仍隨時間呈現(xiàn)波動規(guī)律。圖10(c)給出了均勻流下10個漩渦脫落周期內(nèi)升力系數(shù)變化情況，升力系數(shù)基本不隨時間發(fā)生變化。

圖10 小KC數(shù)及均勻流升力系數(shù)曲線Fig. 10 Amplitude curve of lift coefficient at low KC number

4 結(jié) 語

基于CFD方法對振蕩流下的振動圓柱進(jìn)行了數(shù)值模擬，分析了不同KC數(shù)工況圓柱的水動力特性及泄渦形態(tài)，得到以下結(jié)論：

1) 采用CFD方法復(fù)現(xiàn)均勻流圓柱受迫振動試驗的研究手段是可行的，對比流體力系數(shù)可知，CFD模擬結(jié)果在趨勢和數(shù)值上均與試驗結(jié)果吻合較好。

2) 振蕩流下，在低振幅比、低振動頻率時，振動圓柱的激勵力系數(shù)會隨著KC數(shù)改變發(fā)生顯著變化。小KC數(shù)時，流速變化較快，激勵力系數(shù)受KC數(shù)影響較大，隨著KC數(shù)逐漸增加，流態(tài)趨向均勻流，激勵力系數(shù)數(shù)值差異減小。隨著圓柱振幅比增加，不同KC數(shù)下的流體力系數(shù)基本保持一致。

3) 大KC數(shù)工況下由于流體在單向上可保持較長時間，漩渦有足夠時間發(fā)展和脫落，可觀察到明顯的2S泄渦模式；小KC數(shù)工況時，流體變化加快，脫落的漩渦還未完全消散便被反向流速沖回圓柱，漩渦發(fā)放規(guī)律性較弱。

4) 大KC數(shù)工況下的升力系數(shù)時歷曲線可在同一振蕩流周期中觀察到兩次明顯的振幅調(diào)制現(xiàn)象，分別對應(yīng)正向、反向兩個流速階段；小KC數(shù)工況時，流速頻繁改變大小和方向，升力系數(shù)時歷不像均勻流下穩(wěn)定，但變化不具備明顯規(guī)律性。