劉旭菲， 陳威霖， 及春寧

(1. 浙江水利水電學(xué)院 水利與環(huán)境工程學(xué)院，杭州 310018;2. 天津大學(xué) 水利工程仿真與安全國家重點實驗室，天津 300072)

多圓柱流致振動常見于海洋工程、橋梁以及核工程等，如海洋立管、橋梁纜索和熱交換管等，對該問題的研究具有重要的科學(xué)和工程意義。作為經(jīng)典的流固耦合問題，圓柱流致振動問題在過去的幾十年中得到了廣泛的關(guān)注。Sarpkaya等[1-2]對單圓柱渦激振動的研究成果進(jìn)行了詳細(xì)的綜述。由于復(fù)雜的相互作用，多圓柱流致振動的振動特性和誘發(fā)機(jī)制十分復(fù)雜，明顯不同于單圓柱的情況。Brika等[3-4]發(fā)現(xiàn)在串列布置中，在上游圓柱尾流的影響下，下游圓柱的振動受到顯著的促進(jìn)作用，出現(xiàn)了弛振的現(xiàn)象。陳威霖等[5-6]詳細(xì)討論了小間距比下串列雙圓柱渦激振動中的振動響應(yīng)、流體力及尾流模式等。Huera-Huarte等[7-8]發(fā)現(xiàn)當(dāng)間距比不大時，并列雙圓柱之間存在強(qiáng)烈的相互作用。Chen等[9]發(fā)現(xiàn)在并列布置下，對稱布置的雙圓柱會出現(xiàn)不對稱振動現(xiàn)象。Behara等[10-11]發(fā)現(xiàn)在三角排列三圓柱中，下游兩圓柱在上游圓柱尾流的作用下出現(xiàn)了較大振幅的振動，且出現(xiàn)了顯著的多頻現(xiàn)象。Chen等[12-13]發(fā)現(xiàn)在三角形布置下，隨著間距比(P/D，P為圓柱中心之間的距離，D為圓柱的直徑)的變化，振動響應(yīng)可以劃分為5個不同的區(qū)域，每個區(qū)域內(nèi)圓柱的振動情況與P/D有密切關(guān)系。Ma等[14-15]通過拖曳水池試驗研究了三角排列細(xì)長柔性圓柱的振動模態(tài)隨折合流速變化的情況，對三圓柱之間的尾流/臨近干涉和順流向與橫流向振動之間的相互作用進(jìn)行了細(xì)致地分析。

除了上述單圓柱和獨立支撐多圓柱流致振動問題，剛性耦合地多圓柱振動也常見于實際工程，如海洋工程中采用剛性管卡連接的立管束以及核工程中熱交換管等。與獨立支撐多圓體相比，在剛性耦合下多圓柱振動同步，其振動特性和振動機(jī)理與獨立支撐多圓柱明顯不同。由于各圓柱在耦合系統(tǒng)中的位置不同，它們對系統(tǒng)振動的貢獻(xiàn)是重要的研究方向。本文側(cè)重于揭示剛性耦合三圓柱的振動機(jī)制以及各圓柱對振動系統(tǒng)流體力的貢獻(xiàn)，有助于促進(jìn)對剛性耦合系統(tǒng)的理解。在剛性耦合多柱體的研究上，Munir等[16]對剛性耦合并列雙圓柱的流致振動進(jìn)行了數(shù)值模擬。雷諾數(shù)為Re=1 000，質(zhì)量比為m*=2，間距比為P/D=1.5～4.0，折合流速為Ur=1～20。研究發(fā)現(xiàn)，當(dāng)P/D=1.5～2.0時，系統(tǒng)振動分為兩個分支，鎖定分支內(nèi)振動和升力同頻，而非鎖定分支內(nèi)，振動頻率為升力頻率的一半。當(dāng)P/D=4.0 時，尾流之間的相互作用變?nèi)?，振動鎖定區(qū)域與單圓柱渦激振動相當(dāng)。Han等[17]研究了不同來流攻角下等邊三角形排列剛性耦合三圓柱的流致振動。其中，間距比固定為P/D=4.0，來流攻角α=0°～60°，折合流速Ur=1～30，雷諾數(shù)Re=1 000～3 000。研究發(fā)現(xiàn)，來流攻角和折合流速對渦激振動和脫渦模式都有顯著的影響。當(dāng)α=45°時，振動系統(tǒng)鎖定范圍最廣，而當(dāng)α=60°時，鎖定范圍最窄，在α=0°工況下出現(xiàn)弛振現(xiàn)象(Ur>12)。

盡管前人對剛性耦合多圓柱體的流致振動開展了一些研究，但由于此類問題的影響因素眾多，已有研究僅能窺見一斑?；诖耍疚尼槍嶋H工程中出現(xiàn)的等邊三角形排列剛性耦合三圓柱系統(tǒng)的流致振動進(jìn)行數(shù)值模擬研究，以期對其振動、受力特性和激勵機(jī)制進(jìn)行分析和探討。

1 數(shù)值方法

1.1 控制方程

黏性不可壓縮流動的控制方程[18]如下

(1)

?·u=0

(2)

式中：u為流速；t為時間；p為壓強(qiáng)；ρ為流體密度；ν為運(yùn)動黏滯系數(shù)； ?為梯度算子；f為附加體積力矢量，代表流固耦合邊界條件。

時間離散后的控制方程守恒形式如下

(3)

?·un+1=0

(4)

式中：δt為時間步長； 組合項h=?·[-uu+ν(?u+?uT)]由對流項與擴(kuò)散項組成， 上標(biāo)T為矩陣轉(zhuǎn)置，附加體積力表示為

(5)

式中：Ih和Dh為引入的插值函數(shù)；Vn+1為物面邊界點速度； 上標(biāo)n+1，n+1/2，n，n-1為時間步。

僅做橫流向運(yùn)動的剛性耦合圓柱體，其運(yùn)動可以用式(6)來描述

(6)

1.2 參數(shù)設(shè)置

剛性耦合三圓柱流致振動的模擬中，計算域流向X和橫向Y的寬度分別為Lx和H。Lx=Lu+Ld，三圓柱圓心間距為P，圓心的中點距上游邊界Lu=100D，距下游邊界Ld=100D，H=200D，D為圓柱的直徑，如圖1所示。此外，為保證模擬結(jié)果的準(zhǔn)確性，在圓柱周圍設(shè)置一個大小為8D×8D的加密區(qū)域，區(qū)域內(nèi)流向和橫向網(wǎng)格大小均為1/64D。加密區(qū)域外的網(wǎng)格等比增大。計算域邊界條件設(shè)置如下：入口為Dirichlet型邊界條件；流向速度u=U∞；橫向速度v=0；其中，U∞為無窮遠(yuǎn)處的流速。出口為Neumann邊界條件，上下邊界為自由可滑移邊界。本文計算域和網(wǎng)格的設(shè)置與Chen等的研究相同，其參數(shù)敏感性已經(jīng)在Chen等的研究中證明。

圖1 正三角形排列剛性耦合三圓柱流致振動的計算域和邊界條件示意圖Fig.1 Sketch of computational domain and boundary conditions for flow-induced vibration of three rigidly coupled cylinders in equilateral-triangular arrangements

剛性耦合三圓柱按等邊三角形布置，上游為兩個并排圓柱，下游為一個圓柱，三圓柱的間距比為P/D=1.0～4.0，質(zhì)量比為m*=2，雷諾數(shù)為Re=100和折合流速為Ur=3～30。為了保證圓柱能以更大振幅振動，阻尼比設(shè)為零。耦合三圓柱僅在橫向自由振動。

需要說明的是，本文數(shù)值方法的模擬精度已經(jīng)在獨立支撐的三圓柱渦激振動研究中得到了充分驗證，讀者可自行查閱，此處不再贅述。

2 結(jié)果與討論

2.1 振幅與響應(yīng)頻率

不同間距比下剛性耦合三圓柱的無量綱振幅A/D和頻率比f*隨折合流速變化情況，如圖2所示。由圖2(a)可知：不同間距比下，剛性耦合三圓柱的振動呈現(xiàn)出明顯不同的變化趨勢。當(dāng)P/D=1.0時，振動響應(yīng)表現(xiàn)為典型的弛振響應(yīng)[19-20]，振幅隨折合流速的增加而單調(diào)增加，而振動頻率比f*(f*=fv/fn，fv為振動頻率，fn為系統(tǒng)固有頻率)則一直小于1.0。當(dāng)P/D=1.0時，位移歷時曲線的情況，如圖3(a)～圖3(e)所示。圖3中橫軸為無量綱時間t*=tU∞/D，縱軸為無量綱位移y/D?？傮w上看，圓柱的位移呈現(xiàn)正弦變化，比較規(guī)律，這與相似雷諾數(shù)下三角柱的馳振響應(yīng)類似。當(dāng)Ur=25時，圓柱的振動出現(xiàn)微弱的拍振特征。所謂拍振是指振動響應(yīng)中包含兩個或以上振動分量，并且各個分量的頻率非倍數(shù)關(guān)系，造成了振幅的大小隨時間周期性變化。

圖2 不同間距比下正三角形排列剛性耦合三圓柱振幅和頻率隨折合流速的變化Fig.2 Variations of the vibration amplitude and frequency of three rigidly coupled cylinders in equilateral-triangular arrangements with the reduced velocity under different spacing ratios

當(dāng)P/D=1.6時，振幅明顯小于P/D=1.0時的工況。由圖2(a)可知：當(dāng)Ur≤7.5時，振幅幾乎為零。之后，隨著折合流速的增加，振幅持續(xù)增加，并在很寬的折合流速區(qū)間內(nèi)維持穩(wěn)定的振幅，最大振幅為A/D=0.86。當(dāng)Ur>22以后，振幅開始緩慢下降。由圖2(b)可知：當(dāng)Ur≤7.5時，振動頻率隨折合流速線性增加，對應(yīng)靜止工況下的泄渦頻率斯特勞哈爾數(shù)St。當(dāng)Ur>7.5時，振動頻率偏離St，而沿著一個更小的斜率近似線性增加。在鎖定區(qū)間7.5

當(dāng)P/D=2.5時，以Ur=5.5為界，振動可以分為兩個鎖定區(qū)間。兩個區(qū)間內(nèi)圓柱的振動均呈現(xiàn)為先增后減的趨勢，最大振幅相當(dāng)，但是后一個區(qū)間的寬度明顯大于前一個。由圖2(a)可知：在第一個區(qū)間(Ur=3.5～4.5)內(nèi)，圓柱振動響應(yīng)隨折合流速變化較為劇烈，上升段和下降段均較陡，最大振幅(A/D=0.50)出現(xiàn)在Ur=4處。由圖2(b)可知：在第二個區(qū)間(Ur=6～15)內(nèi)，振幅的上升段和下降段均較緩，最大振幅(A/D=0.55)在Ur=10處取得。與P/D=1.0，P/D=1.6工況相比，振幅明顯更小一些。振動頻率在兩個區(qū)間內(nèi)也呈現(xiàn)不同的變化趨勢。在第一個區(qū)間內(nèi)，振動頻率維持在f*=0.8附近，振動處于鎖定狀態(tài)。在第二個區(qū)間內(nèi)，振動頻率維持在f*=0.8～1.2內(nèi)，振幅也較大一些。在兩個鎖定區(qū)間之間(Ur=5.5)和之外(Ur=3和Ur>15)，振幅較小，振動頻率與靜止工況下的St接近，隨著折合流速線性增加，振動處于非鎖定狀態(tài)。當(dāng)P/D=2.5時圓柱位移的歷時曲線，如圖3(k)～圖3(o)所示。在兩個鎖定區(qū)間內(nèi)部，圓柱的振動均較為規(guī)律，而在區(qū)間交界處，振動響應(yīng)呈現(xiàn)拍振特征?？偟膩砜?，當(dāng)P/D=2.5時，剛性耦合三圓柱的振動響應(yīng)呈現(xiàn)出雙鎖定區(qū)間的特征，這與單圓柱的單鎖定區(qū)間形成了鮮明的對比。此外，兩鎖定區(qū)間的尾渦與圓柱之間的相互作用有明顯差異，產(chǎn)生了截然不同的激勵作用。

當(dāng)P/D=4.0時，振幅隨折合流速的變化情況與P/D=2.5工況的類似，可劃分為振幅較大的第一鎖定區(qū)間(Ur=4.0～8.5)和振幅較小的第二鎖定區(qū)間(Ur=10～12)。由圖2(a)可知：在第一區(qū)間內(nèi)，隨著折合流速的增加，振幅迅速增加，并在Ur=5.5時達(dá)到最大值(A/D=0.63)。之后，隨著折合流速的增加，振幅緩慢下降，直到第一個區(qū)間結(jié)束；在第二區(qū)間內(nèi)，振幅的變化明顯小于第一區(qū)間的情況，最大振幅僅為A/D=0.25。由圖2(b)可知：雖然第一區(qū)間的振幅明顯大于第二區(qū)間，兩個區(qū)間的振動頻率均在f*=0.8～1.2內(nèi)。兩個區(qū)間之間，振動頻率出現(xiàn)了明顯的跳躍和不連續(xù)現(xiàn)象。當(dāng)P/D=4.0時，圓柱位移的歷時曲線，如圖3(p)～圖3(t)所示。在第一個鎖定區(qū)間內(nèi)，振動響應(yīng)均為規(guī)律的正弦曲線；而在第二個鎖定區(qū)間內(nèi)，振動響應(yīng)表現(xiàn)為規(guī)律的非對稱振動，呈現(xiàn)出明顯的多頻特征，這與圓柱尾流之間復(fù)雜的相互作用有關(guān)。

圖3 不同折合流速下正三角形排列剛性耦合三圓柱的位移歷時曲線Fig.3 Displacement histories of three rigidly coupled cylinders in equilateral-triangular arrangements at different reduced velocities

綜合以上的結(jié)果可知，不同間距比下三角形排列剛性耦合三圓柱流致振動響應(yīng)可大致分為兩類：一類是間距比較小情況下的馳振響應(yīng)(不穩(wěn)定性激勵)，圓柱振動幅值隨折合流速近似線性增大，振動頻率小于體系的固有頻率，這與直邊迎流的三角柱的振動響應(yīng)類似，換句話說，由于圓柱之間的間隙流較弱，3個剛性耦合的圓柱可近似視為一個擴(kuò)展了的三角柱；另一類為間距比適中和較大情況下的渦激振動響應(yīng)(共振激勵)，但不同間距比下會產(chǎn)生單鎖定區(qū)間和雙鎖定區(qū)間的振動響應(yīng)。由于圓柱之間距離較小，P/D=1.6工況下僅有一個鎖定區(qū)間，而且對比鎖定區(qū)間內(nèi)的頻率特征可以發(fā)現(xiàn)，P/D=1.6工況的單鎖定區(qū)間是與P/D=2.5，P/D=4.0雙鎖定的第一鎖定區(qū)間相對應(yīng)。

2.2 流體力系數(shù)

圖4 不同間距比下三角形排列剛性耦合三圓柱的流體力系數(shù)隨折合流速的變化Fig.4 Variations of the hydrodynamic force coefficients of three rigidly coupled cylinders in equilateral-triangular arrangements with the reduced velocity under different spacing ratios

不同間距比下脈動阻力C′D,t隨折合流速變化的情況，見圖4(b)。P/D=1.0工況的脈動阻力明顯大于其余間距比的情況。與阻力均值相似，脈動阻力在Ur>14以后幾乎不再變化，維持在C′D,t=2.2附近。當(dāng)P/D=1.6時，脈動阻力隨折合流速的增加先增加后減小，呈現(xiàn)單峰特征。當(dāng)P/D=2.5，P/D=4.0時，脈動阻力隨折合流速的變化趨勢與阻力均值相似，分別呈現(xiàn)雙峰和單峰特征。

不同間距比下脈動升力C′L,t隨折合流速變化的情況，見圖4(c)。從總體上看，除個別折合流速，P/D=1.0工況的脈動升力要比其余工況的更大一些，隨著折合流速的增加，脈動升力先增后減，最大值(C′L,t=3.21)出現(xiàn)在Ur=7處。而當(dāng)P/D=1.6時，脈動升力在Ur=8.5時達(dá)到最大值(C′L,t=1.19)，之后開始緩慢下降。當(dāng)P/D=2.5，P/D=4.0時，脈動升力均表現(xiàn)出雙峰特征，盡管P/D=4.0工況的第二個峰值(位于Ur=11處)很小且平緩。

2.3 尾流模式和激勵機(jī)制

為深入分析各間距比下圓柱振動與流場之間的耦合作用，本節(jié)結(jié)合尾流模式給出詳細(xì)的機(jī)理解釋。當(dāng)P/D=1.0時，三圓柱緊密靠在一起，圓柱之間不存在間隙流，如圖5所示。圖5中，x/D和y/D分別為無量綱的流向和橫向坐標(biāo)。剪切層在剛性耦合三圓柱的上、下兩側(cè)形成旋渦，并交替泄放。由圖5(a)可知： 當(dāng)Ur=5時，旋渦從圓柱的兩側(cè)脫落，形成典型的2S泄渦模式。隨著圓柱振幅的增加，單個振動周期內(nèi)脫落的旋渦個數(shù)明顯增多，形成了多渦模式，見圖5(c)～圖5(e)。需要說明的是，由于泄渦長度較長，圖5(c)～圖5(e)所示范圍未能顯示一個完整振動周期內(nèi)的旋渦泄放。例如，圖5(e)所示僅為半個振動周期內(nèi)的旋渦泄放過程。馳振頻率與泄渦頻率相差很大，因此馳振不是由交替泄渦激發(fā)的，而是由三圓柱運(yùn)動時上、下表面的時均壓差引起和維持的。

當(dāng)P/D=1.6時，圓柱間隙間形成剪切層。由于受到下游圓柱的束縛，上游圓柱的間隙剪切層僅發(fā)生小幅擺動，并不能脫落旋渦。上游圓柱的外側(cè)剪切層可脫落旋渦，并與下游圓柱的剪切層融合，形成擴(kuò)展的2S泄渦模式。對比P/D=1.0的工況可以發(fā)現(xiàn)，間隙流的出現(xiàn)使得P/D=1.6工況下旋渦更弱，且距離圓柱更遠(yuǎn)。此工況下，振動頻率和升力頻率相等，圓柱的振動由交替泄渦(擴(kuò)展2S模式)驅(qū)動。

當(dāng)P/D=2.5時，間隙流的作用更加顯著，上游圓柱柱間剪切層大幅擺動，與下游圓柱發(fā)生復(fù)雜相互作用，一些情況下能脫落旋渦。不同間隙比條件下，柱間剪切層和下游圓柱的相互作用存在兩種截然不同的模式。當(dāng)Ur=4時(第一鎖定區(qū)間)，上游圓柱間隙剪切層交替脫落旋渦，并“抨擊”在下游圓柱上，隨后發(fā)生分裂，形成兩個較小的旋渦，如圖6(a)所示。此時，圓柱受到的脈動升力和阻力均較大(見圖4(b)～圖4(c))，振幅較大。上、下游圓柱的泄渦發(fā)生復(fù)雜的融合后，在下游形成較寬的平行渦街。但是，當(dāng)折合流速更大時(Ur=6，Ur=10，Ur=12，Ur=20，第二鎖定區(qū)間，如圖6(b)～圖6(e)所示)，上游圓柱的剪切層不與下游圓柱發(fā)生直接作用，而是與下游圓柱的剪切層發(fā)生融合，形成擴(kuò)展的交替渦街?？梢?，P/D=2.5工況下存在兩種振動激勵模式：一種是低折合流速條件下的剪切層重附著激勵模式；另一種是中、高折合流速條件下的交替泄渦(擴(kuò)展2S模式)激勵模式。

圖5 不同折合流速下正三角形排列剛性耦合三圓柱的尾渦模式(P/D=1.0，P/D=1.6)Fig.5 Wake patterns of three rigidly coupled cylinders in equilateral-triangular arrangements at different reduced velocities (P/D=1.0, P/D=1.6)

當(dāng)P/D=4.0時，旋渦可以自由地從3個圓柱上脫落。當(dāng)Ur=5時(第一鎖定區(qū)間)，從上游圓柱內(nèi)側(cè)脫落的旋渦交替地重附著于下游圓柱上，如圖6(f)所示。此時，圓柱的脈動升力較大，振動由剪切層重附著驅(qū)動。旋渦經(jīng)過復(fù)雜的融合后，在下游形成較寬的平行渦街。當(dāng)Ur=8時，泄渦不規(guī)律，如圖6(g)所示。當(dāng)Ur=11時(第二鎖定區(qū)間)，上下游圓柱各自交替脫落旋渦，但上游圓柱的泄渦近似同相位，如圖6(h)所示。從上游圓柱脫落的旋渦不直接與下游圓柱作用，而是與下游圓柱的旋渦共存，形成雜亂無章的尾渦模式。此時，振動由單個圓柱的交替泄渦驅(qū)動。進(jìn)一步增大折合流速(Ur=15，Ur=25)，兩個上游圓柱的泄渦近似反相位(如圖6(i)～圖6(j)所示)，上游圓柱的升力相互抵消，圓柱振幅很小。

圖6 不同折合流速下正三角形排列剛性耦合三圓柱的尾渦模式(P/D=2.5,P/D=4.0)Fig.6 Wake patterns of three rigidly coupled cylinders in equilateral-triangular arrangements at different reduced velocities (P/D=2.5, P/D=4.0)

3 結(jié) 論

本文對等邊三角形排列剛性耦合三圓柱的渦激振動問題進(jìn)行了數(shù)值模擬研究，其中上游為兩個并排的圓柱，而下游為一個圓柱。對不同間距比下振幅、頻率和各圓柱的流體力貢獻(xiàn)進(jìn)行了詳細(xì)討論，最后結(jié)合尾渦模式，對圓柱與流場之間的相互作用進(jìn)行了分析。

研究發(fā)現(xiàn)三角形排列剛性耦合三圓柱的振動模式分可以為兩類，分別為P/D=1.0的馳振模式和P/D=1.6～4.0的渦振模式。當(dāng)P/D=1.0時，無間隙流動，三圓柱形成一個擴(kuò)展的整體，系統(tǒng)振動表現(xiàn)為直邊迎流的三角柱相似的弛振特征。當(dāng)系統(tǒng)振動較小時，從三圓柱表面交替脫落的旋渦形成2S模式，隨著振幅的增加，單個振動周期內(nèi)會脫落更多的旋渦。當(dāng)P/D=1.6時，間隙流較弱，間隙內(nèi)剪切層小幅擺動，但不會形成旋渦，三圓柱后方形成2S模式尾渦，并驅(qū)動三圓柱發(fā)生渦激振動。振動僅出現(xiàn)單個較寬的鎖定區(qū)間。當(dāng)P/D=2.5，P/D=4.0時，圓柱的振動呈現(xiàn)為兩個鎖定區(qū)間，分別對應(yīng)不同的振動激勵機(jī)制。第一個鎖定區(qū)間內(nèi)，上游圓柱的間隙剪切層脫落旋渦，并撞擊到下游圓柱上，引起下游圓柱較大的脈動升力，從而激發(fā)圓柱振動，即剪切層重附著激勵模式；第二個鎖定區(qū)間內(nèi)，上游圓柱的剪切層不與下游圓柱發(fā)生直接作用，而是與下游圓柱的剪切層發(fā)生融合，形成交替渦街，激發(fā)圓柱振動，即交替泄渦驅(qū)動模式。

需要說明的是由于本文僅研究了層流條件下剛性耦合三圓柱的振動特性和機(jī)理，然而在實際工程中流動為紊流，雷諾數(shù)多在105～107量級，因此，研究成果推廣到實際工程具有一定的局限性。然而，盡管高雷諾數(shù)條件下剛性耦合三圓柱系統(tǒng)具有更大振幅，各圓柱尾流之間的相互作用更加強(qiáng)烈，但是一些現(xiàn)象以及機(jī)理是相通的。例如，本文的研究發(fā)現(xiàn)在間距比不大時，剛性耦合三圓柱表現(xiàn)為弛振現(xiàn)象，而在Han等的高雷諾數(shù)類似工況下同樣也出現(xiàn)了弛振現(xiàn)象，其機(jī)理也與本文提到的剪切層的重附著有關(guān)。綜上，本文研究成果對于深入理解實際工程的高雷諾數(shù)剛性耦合三圓柱渦激振動的現(xiàn)象和機(jī)制有一定幫助。