剛性耦合三圓柱流致振動(dòng)特性和機(jī)制

劉旭菲， 陳威霖， 及春寧

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

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

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

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

1 數(shù)值方法

1.1 控制方程

黏性不可壓縮流動(dòng)的控制方程[18]如下

(1)

?·u=0

(2)

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

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

(3)

?·un+1=0

(4)

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

(5)

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

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

(6)

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

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

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

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

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

2 結(jié)果與討論

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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