• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    低雷諾數(shù)下串聯(lián)雙圓柱渦激振動機理的數(shù)值研究

    2014-09-05 06:43:14郭曉玲唐國強劉名名
    振動與沖擊 2014年4期
    關鍵詞:尾渦約化橫流

    郭曉玲, 唐國強,2, 劉名名, 呂 林,3, 滕 斌

    (1.大連理工大學 海岸和近海工程國家重點實驗室,大連 116024;2.韓國海洋科學技術研究所,大田 305-343;3.大連理工大學 深海工程研究中心,大連 116024)

    在海洋工程中,當水流經(jīng)過海洋平臺的樁柱、支撐結構、立管及海底管線等這些非流線型結構時,會在尾流區(qū)內產(chǎn)生交替的旋渦脫落,使作用在結構上的流體作用力表現(xiàn)出明顯的周期性特征,進而誘發(fā)渦激振動(Vortex-Induced Vibration,VIV)。當渦旋的脫落頻率與結構的固有頻率相接近時,將發(fā)生“鎖定”現(xiàn)象。此時,結構物會以較大的振幅振動,導致結構物的壽命急劇降低。在過去的幾十年中,均勻流中單個圓柱的渦激振動問題得到廣泛的關注[1-6]。人們通過對高雷諾數(shù)實驗研究發(fā)現(xiàn),在低質量-阻尼比(m*ξ)下,圓柱的位移響應可以分為三個階段[1]:初始分支(Initial branch)、上分支(Upper branch)及下分支(Lower branch),而在高質量-阻尼比(m*ξ)情況下,圓柱的位移響應只存在初始分支和下分支[2]。通過強迫振動的實驗研究發(fā)現(xiàn)[3]:對于僅發(fā)生橫流向受迫振動的圓柱結構,其尾渦脫落模式通常可劃分為2P、2S和P+S三種形式,從Khalak等[1]中可以看到,位移響應的初始分支對應的渦脫落模式為2S,下分支對應的為2P模式。

    對于海洋工程結構中近距離放置的多個圓柱的情況,由于圓柱之間存在明顯的水動力干涉作用,相應的渦激振動響應與單個圓柱的情況存在顯著的差別[7-11]。針對上游圓柱固定不動,下游圓柱發(fā)生渦激振動的實驗研究表明[7],由于受到上游圓柱尾跡區(qū)的影響,下游圓柱開始發(fā)生鎖定時所對應的約化速度大于單個圓柱的情況。并且,下游圓柱的鎖定帶寬與兩圓柱的間距有關。例如,Brika等[7]的結果表明,當兩圓柱間的相對圓心距L/D=10.0時,以約化速度表征的下游圓柱的鎖定區(qū)間約為單圓柱情況的2倍。Carmo等[8]分別計算了二維(Re = 150)和三維(Re = 300)條件下,彈性支撐圓柱在上游固定圓柱尾流干涉作用下的橫流向單自由度振動問題。計算條件為圓心間距分別為LX/D=1.5、3.0、5.0和8.0,質量比為m*=2.0,阻尼比為ξ=0.007。數(shù)值結果表明,在Re=150條件下,當發(fā)生共振時,下游圓柱的振幅較單圓柱情況增大近50%,而對于三維情況,在不同的圓心間距下,下游圓柱的振動幅值隨約化速度的變化趨勢與對應的二維情況類似。在以往的研究工作中,考慮到圓柱發(fā)生渦激振動時順流向位移的振幅比較小,因此,順流向的振動往往被忽略[8-9]。而Vandiver等[10]的研究工作表明,順流向和橫流向存在明顯的相互作用。陳文曲等[11]也對這一方面開展了數(shù)值分析工作。但是,人們對低雷諾數(shù)流動中,處于固定圓柱尾跡干涉作用下,同時具有順流向和橫流向雙自由度的圓柱渦激振動問題仍缺乏深刻的認識,特別是圓柱間距及質量比對圓柱受力和渦激振動響應的影響還需要更加細致的研究。

    以目前現(xiàn)有的計算能力,對大尺度(幾千米水深)范圍內的實際立管渦激振動開展全場三維數(shù)值分析還面臨很大的挑戰(zhàn),計算時間過長,不能滿足實際工程設計的需要。但有關低雷諾數(shù)下,多圓柱渦激振動問題的數(shù)值研究對于發(fā)展可靠的數(shù)值模型以及建立快捷的近似工程分析預報方法是有借鑒意義的。因此本文的工作主要在于加深對多圓柱渦激振動現(xiàn)象物理本質的認識,更側重于基礎理論方面的研究。本文采用三步有限元數(shù)值方法對二維Navier-Stokes方程進行了求解,并結合任意拉格朗日-歐拉(ALE)動網(wǎng)格方法,建立和開發(fā)了相應的數(shù)值模型和計算程序,以Re=150為代表,重點研究了兩個等直徑串聯(lián)圓柱的渦激振動問題,其中上游圓柱固定不動,下游圓柱在彈性支撐和阻尼作用下允許同時發(fā)生順流向和橫流向振動。通過數(shù)值計算,來研究圓柱間隙比、質量比以及約化速度等因素對下游圓柱渦激振動特性的影響。

    1 數(shù)值模型

    1.1 流動控制方程

    二維不可壓縮均勻粘性牛頓流體運動的基本控制方程為連續(xù)性方程和Navier-Stokes方程,在ALE參考坐標系下,可以表示成如下的無量綱形式:

    (1)

    (2)

    其中,xi表示笛卡爾坐標(二維情況下i= 1、2,分別對應x和y方向),ui為xi方向的流速分量,p為壓力,Re=UD/υ為雷諾數(shù),U為均勻來流流速,D為圓柱直徑,υ為流體的運動學粘性系數(shù),cj為對流速度,并有:

    (3)

    應用三步有限元方法[12-13],可得到如下離散形式的動量方程:

    (4)

    (5)

    (6)

    (7)

    對于時間推進,本文根據(jù)CFL(Courant-Friedrichs-Lewy)條件,采用如下的動態(tài)時間步長:

    式中,Se為網(wǎng)格的面積,ue為網(wǎng)格中心點的流速,min(*)表示在計算域內取最小值,Cs為安全系數(shù),文中取為Cs= 0.2。

    由于上述離散格式具有高階Tayler-Galerkin性質,同時速度與壓力通過投影(Projection)過程進行了解耦,因此式(4)~式(6)以及式(7)均可采用統(tǒng)一的標準Galerkin有限元方法進行空間離散,相關的離散方法可參見文獻[12~14]。

    當獲得流場和壓力場后,圓柱所受到的流體作用力可通過對圓柱表面的壓力及粘性剪切力進行表面積分求得,無量綱化的拖曳力系數(shù)CD和升力系數(shù)CL分別為:

    (8)

    (9)

    式中,p,u,v分別為圓柱表面上某點的壓力、x方向速度、y方向速度,θ為該點與圓心的連線與x軸正向之間的夾角。

    1.2 圓柱的運動方程

    (10)

    (11)

    當圓柱的位移確定后,本文通過微分網(wǎng)格變換方法進行網(wǎng)格更新[15],進而獲得新的網(wǎng)格坐標以及網(wǎng)格運動速度,用于下一時間步的流場求解。

    1.3 計算模型及邊界條件

    計算模型及邊界條件如圖1所示。兩等直徑圓柱串聯(lián)放置,無量綱的圓柱直徑為d=1。設坐標原點位于下游圓柱的初始圓心,兩圓柱圓心間距為LX,相應的無量綱參數(shù)為LX/D。計算域的順流向長度為55D,橫流向寬度為50D。在入口處指定無量綱的速度u=1,v=0;側壁采用了對稱邊界條件?u/?y=0,v= 0;出口處的速度邊界條件為?ui/?t+c?ui/?xi=0,其中c為局部平均流速;圓柱表面施加不可滑移邊界條件u=dx/dt,v=dy/dt。計算中,在出口處指定相對壓力p= 0,在其他邊界采用?p/?n=0的壓力邊界條件,n為指出流體域的外法向單位矢量。在初始時刻,流場中的速度及相對壓力分布均設為零(即初始速度場滿足連續(xù)方程)。

    圖1 計算域及邊界條件示意圖

    2 模型驗證

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

    表1 網(wǎng)格收斂性驗證結果

    2.2 圓柱橫流向渦激振動的數(shù)值驗證

    為進一步驗證本文數(shù)值模型的可靠性,同時考慮到在之前的工作中已對固定單圓柱繞流問題進行了詳細的驗證[13],以下僅對單個圓柱(Isolated)橫流向自激振動以及上游圓柱固定,下游圓柱(直徑同于上游圓柱)的橫流向自激振動問題分別進行計算。在計算的過程中,均采用了相同的雷諾數(shù)Re=150、質量比m*=2.0和阻尼比ξ=0.007。圖2給出了單圓柱橫流向最大振幅(Amax/D)隨約化速度(Ur)的變化及其與Carmo等[8]的對比結果。對于串聯(lián)雙圓柱的情況,考慮了LX/D= 1.5和3.0兩種條件下,下游圓柱的橫流向最大振幅對約化速度的依賴關系。從圖2中可以看出,本文的數(shù)值計算結果(present)與已有數(shù)值結果吻合良好,表明本文的數(shù)值模型具有良好的數(shù)值精度。

    3 計算結果及分析

    在Carmo等[8]的數(shù)值分析工作中,考慮了上游圓柱固定不動的情況下,下游圓柱只在橫流向發(fā)生渦激振動的問題。本文將對下游圓柱同時發(fā)生橫流向和順流向雙自由度渦激振動的問題進行研究。為便于與文獻[8]的結果進行對比分析,本文選擇與該文相同的雷諾數(shù)Re=150、阻尼比ξ=0.007以及圓心間距比LX/D=3.0、5.0和8.0,在約化速度Ur= 3.0 ~ 12.0的范圍內對下游圓柱的兩自由度渦激振動響應特性開展研究。同時,本文進一步以m*= 5.0、10.0、20.0為例,研究質量比對下游圓柱發(fā)生渦激振動時動力響應的影響。

    圖2 單圓柱及下游圓柱橫流向最大振幅Amax /D隨約化速度Ur的變化(與文獻[8]的對比,Re=150,m* = 2.0,ξ = 0.007)

    3.1 圓心間距比的影響

    根據(jù)固定串聯(lián)圓柱的不同渦脫落模態(tài)[16-17],本文選擇了三個代表性的圓心間距比,LX/D= 3.0、5.0和8.0,并與相同雷諾數(shù)Re=150下的單圓柱(Isolated)渦激振動情況進行了比較,分析了圓柱的位移響應、振動頻率、受力特性、尾渦脫落模式及運動軌跡的情況。

    3.1.1 下游圓柱的位移及振動頻率

    圖3(a)以質量比m*=10.0、阻尼比ξ=0.007為例,給出不同圓心間距下,下游圓柱橫流向位移均方根YRMS/D隨Ur的變化。從圖中可以看出,當Lx/D=3.0時,其橫流向位移響應曲線同單圓柱類似,只有一個峰,而Lx/D=5.0和8.0的YRMS/D曲線均出現(xiàn)了兩個峰值。在Lx/D=3.0時,下游圓柱的最大位移均方根比單圓柱的情況增大近50%,而對于Lx/D=5.0和8.0,最大位移均方根約等于單圓柱的位移。從圖3(a)中還可以看出,在Ur> 7.0的高約化速度范圍內,與單圓柱的位移響應相比,下游圓柱的橫流向位移依然維持在比較高的水平,并且Lx/D=5.0和8.0的位移均方根在幅值和變化趨勢上都很接近,二者均比Lx/D=3.0所對應的位移明顯偏大。

    圖3 不同圓心間距比LX / D下,下游圓柱渦激振動響應特性隨約化速度Ur的變化關系(Re = 150,m* = 10.0,ξ = 0.007)

    由圖3(b)中下游圓柱的順流向平均位移XM/D隨約化速度Ur的變化曲線可以看到,隨著Ur的增大,下游圓柱的順流向平均位移也逐漸增大。同時,由于受到上游固定圓柱遮蔽效應的影響,下游圓柱的順流向平均位移在各圓心間距下都小于單圓柱的情況。對于Lx/D=3.0的情況,在Ur≤5.5時,下游圓柱的XM/D出現(xiàn)負值,即下游圓柱向上游偏移,當發(fā)生鎖定時,XM/D明顯增大,在Ur>7.0后,XM/D基本隨Ur呈線性增加。對于Lx/D=5.0和8.0的情況,順流向平均位移與橫流向位移均方根的變化趨勢類似,在發(fā)生鎖定前各出現(xiàn)一個次峰,鎖定后,二者接近重合,并且隨著Ur的增大基本呈線性增加,在量值上與Lx/D=3.0的情況接近。

    圖4給出了不同圓心間距比下,下游圓柱相對振動頻率f/fn隨Ur的變化關系及其與單圓柱情況的比較。對于圓柱結構的渦激振動問題,當圓柱的振動頻率與固有頻率接近時(即f/fn趨近于1.0,以下以0.95≤f/fn≤1.05為界),一般可認為發(fā)生鎖定[18]。圖4表明,單圓柱的鎖定區(qū)間約為5.0 ≤Ur≤6.5,而對于串聯(lián)雙圓柱的情況,在Lx/D=3.0時,下游圓柱的鎖定區(qū)間為6.0≤Ur≤7.0,當Lx/D=5.0和8.0時,在本文考慮的約化速度范圍內,鎖定區(qū)間分別為 6.0≤Ur≤12.0和6.5≤Ur≤12.0。圖4關于鎖定區(qū)間的結果與圖3 (a)的位移響應呈現(xiàn)出良好的對應關系,即圓柱振動的鎖定區(qū)間對應著橫流向位移比較大的區(qū)間。從圖4中同時也可以發(fā)現(xiàn),隨著Lx/D的增大,開始發(fā)生鎖定的臨界約化速度也呈現(xiàn)出增大的趨勢,并且對于Lx/D=5.0和8.0的情況,鎖定區(qū)間明顯大于單圓柱及Lx/D=3.0的情況。值得注意的是,在Lx/D=5.0和8.0時,橫流向位移的第一個峰值出現(xiàn)在鎖定發(fā)生之前。

    圖4 不同圓心間距比LX / D下,下游圓柱橫流向相對振動頻率f/fn隨約化速度Ur的變化關系(Re = 150,m* = 10.0,ξ = 0.007)

    3.1.2 下游圓柱的受力特性

    圖5 不同圓心間距比LX / D下,下游圓柱受力特性隨約化速度Ur的變化關系(Re=150,m*=10.0,ξ= 0.007)

    3.1.3 尾渦脫落模式

    圖6分析了Lx/D及Ur對尾渦脫落模式的影響作用。對于Lx/D=3.0的串聯(lián)雙圓柱,數(shù)值結果表明,在Ur<6.0的范圍內,上游固定圓柱后方?jīng)]有渦旋脫落,圓柱系統(tǒng)尾流區(qū)內的旋渦以圖6(a1)和6(a2)的模式進行整體脫落;當Ur≥6.0,即開始發(fā)生鎖定后,上游圓柱后方開始出現(xiàn)渦脫落,如圖6(a3)所示;隨著Ur的進一步增大,下游圓柱不再發(fā)生鎖定時,其尾渦脫落模式(參見圖6(a4))相對鎖定區(qū)內的圖6(a3)未見明顯變化。當Lx/D=5.0時,在Ur<5.5的范圍內,典型尾渦脫落模式如圖6(b1)所示,上游固定圓柱后方發(fā)生交替的渦旋脫落,而下游圓柱的尾渦則整齊地排成兩排;當Ur=5.5時,尾渦模式如圖6(b2)所示,整體上與圖6(a1)及6(a2)的情況類似;當Ur≥6.0后,尾渦模式如圖6(b3)和6(b4)所示,上、下游圓柱尾流區(qū)內分別可見旋渦的脫落。從前文的圖3和圖5中看到,在Lx/D=5.0,Ur=5.5的條件下,其位移和受力會突然變小,這正是由于當Ur=5.5時,尾渦脫落模式的突然改變所致。對于Lx/D=8.0,由于圓柱間距比較大,隨著Ur的變化,上下游圓柱的尾流區(qū)內始終同時存在旋渦的脫落,參見圖6(c1)-(c4),上游圓柱脫落的渦旋經(jīng)充分發(fā)展后,與下游圓柱的尾渦發(fā)生更為明顯的干涉作用。

    3.1.4 下游圓柱的運動軌跡

    為說明下游圓柱在x和y兩個自由度方向運動的相互影響作用,圖7給出了不同圓心間距比和約化速度下的運動軌跡。其中,圖7(a1)、圖7(a2)、圖7(b1)和圖7(c1)中右下角分別給出了放大后的軌跡曲線。從圖7(a1)-(a4)中可以看到,在相同的圓心間距比Lx/D=3.0下,隨著Ur的增大,下游圓柱的運動軌跡始終是“8”字形的封閉曲線。而對于Lx/D=5.0和8.0,當Ur≤5.0時,下游圓柱的運動軌跡類似于Lx/D=3.0的情況,也為“8”字形。當Ur>5.0后,圓柱的軌跡則變?yōu)椴环忾]的且不規(guī)則的曲線,如圖7(b3)-(b4)及7(c3)-(c4)所示。這說明圓心間距比和約化速度對下游圓柱的運動有明顯的影響。當Lx/D比較小時,下游圓柱在兩自由度方向上運動的相關性并不隨約化速度有明顯的變化。而當Lx/D較大時,Ur> 5.0后,圓柱在x和y方向運動的相互影響作用則由簡單變得復雜。

    圖6 不同圓心間距比Lx /D及約化速度Ur下的尾渦脫落模式(Re = 150,m*=10.0,ξ=0.007)

    圖7 不同圓心間距比Lx / D及約化速度Ur下下游圓柱的運動軌跡(Re=150,m* = 10.0,ξ = 0.007)

    3.2 質量比對下游圓柱渦激振動的影響

    質量比是圓柱結構渦激振動響應的重要影響因素之一,但對于多圓柱問題,相對于前述的間隙比和約化速度,人們對質量比的影響作用還缺乏細致的研究和足夠的了解。在此,本文將以Lx/D=5.0為例,考慮三種不同的質量比m*= 5.0、10.0和20.0,對這一問題開展進一步的數(shù)值計算,分析質量比對圓柱的位移響應、振動頻率、鎖定區(qū)間、受力特性、尾渦脫落模式及運動軌跡的影響。

    3.2.1 下游圓柱的位移及振動頻率

    圖8首先給出了下游圓柱橫流向相對振動頻率f/fn隨約化速度Ur的變化情況。從圖中可以看出,在本文所考慮的參數(shù)范圍內,質量比對鎖定區(qū)間的影響作用不大,三種質量比下,鎖定區(qū)間所對應的約化速度范圍均為6.0≤Ur≤12.0。

    圖9給出了不同質量比下,下游圓柱的振動響應隨約化速度的變化情況。圖9(a)的結果表明,在鎖定區(qū)間內,下游圓柱的橫流向位移均方根總體上隨約化速度的增大而逐漸減小。同時,在相同Ur下,m*越小,YRMS/D越大。而在未發(fā)生鎖定的低約化速度范圍內,在m*=20.0的條件下,YRMS/D隨著Ur的增大持續(xù)增大,直到發(fā)生鎖定時達到極值;對于相對較小的質量比m*=5.0和10.0,二者的YRMS/D在發(fā)生鎖定前均出現(xiàn)跳躍。進一步的流場分析結果表明,m*=5.0,Ur=4.5和m*=10.0,Ur=5.5兩個跳躍點均由尾渦脫落模式的改變所致,這說明質量比的變化也會對渦激振動過程中的流動模式和位移響應產(chǎn)生影響作用。

    圖9(b)考察了下游圓柱順流向平均位移XM/D隨Ur的變化情況。從圖中可以看出,順流向平均位移總體上均隨著Ur的增大呈上升趨勢,質量比越大,圓柱向下游偏移的平均位移越小。對于m*=5.0和10.0,XM/D也分別在Ur=4.5和5.5處發(fā)生跳躍,這與圖9(a)的情況類似。

    值得說明的是,在m*=5.0和20.0的情況下,相對于m*=10.0,下游圓柱在x,y兩個方向的聯(lián)合位移形態(tài)沒有本質性的改變,與圖7(b1)-(b4)所給出的運動形式在總體上基本一致。這表明,在本文所考查的參數(shù)范圍內,質量比對下游圓柱縱向-橫向運動的相干性的影響作用有限。

    圖8 不同質量比m*下,下游圓柱橫流向相對振動頻率f / fn隨約化速度Ur的變化關系(Re = 150,Lx / D = 5.0,ξ = 0.007)

    圖9 不同質量比m*下,下游圓柱渦激振動響應特性隨約化速度Ur的變化關系(Re =150,Lx / D = 5.0,ξ = 0.007)

    3.2.2 下游圓柱的受力特性

    圖10 不同質量比m*下,下游圓柱受力特性隨約化速度Ur的變化關系(Re=150,Lx/D=5.0,ξ=0.007)

    4 結 論

    本文利用迎風有限元數(shù)值方法對不可壓縮粘性流體的Navier-Stokes方程進行了求解,并結合ALE動網(wǎng)格方法,對Re=150條件下,串聯(lián)雙圓柱的渦激振動問題開展了二維數(shù)值分析。其中,上游圓柱固定不動,下游圓柱在彈簧和阻尼約束下允許同時發(fā)生順流向和橫流向的運動。在對數(shù)值模型的可靠性進行了驗證的基礎上,重點研究了兩等直徑圓柱圓心間距、質量比以及約化速度對下游圓柱渦激振動響應的影響作用,包括位移響應、振動頻率特征、鎖定區(qū)間、流體作用力、尾渦模式以及運動軌跡等方面。數(shù)值結果表明,在本文的計算范圍內,隨著兩圓柱間相對圓心間距Lx/D的增大,下游圓柱開始發(fā)生鎖定時對應的約化速度增大,在較大的間距比下(Lx/D=5.0和8.0),鎖定區(qū)間較單圓柱及小間距情況(Lx/D=3.0)明顯變寬。在Lx/D=3.0的小間距下,下游圓柱的橫流向振動位移可達到單圓柱情況下的1.5倍,同時下游圓柱的順流平均位移出現(xiàn)負值,即下游圓柱的平衡位置向上游圓柱移動。數(shù)值結果也表明,隨著約化速度和間距比的變化,串聯(lián)圓柱的尾渦模式會發(fā)生顯著改變,進而深刻影響圓柱的受力和運動響應。本文還以Lx/D= 5.0為典型間距,研究了質量比對下游圓柱渦激振動響應和受力特性的影響。結果表明,質量比對鎖定區(qū)間的影響不明顯,但當鎖定發(fā)生時,在相同約化速度下,隨著質量比的增大,下游圓柱的橫流向運動位移響應會有所減小。進一步的流場分析表明,質量比的變化也會導致尾渦脫落模式的改變,進而影響到圓柱的受力和運動特性,這在進入鎖定區(qū)間前非常明顯。

    目前關于兩自由度渦激振動問題的研究,通常固定雷諾數(shù),考慮運動響應等對約化速度的依賴關系。但對實際問題而言,往往是圓柱結構形式和流體物性一定,來流流速發(fā)生改變,此時雷諾數(shù)和約化速度同時發(fā)生變化。因此,固定雷諾數(shù)而僅改變約化速度的研究方法不能滿足全面了解渦激振動特性的需要,這是在以后的研究工作中需要加以關注的問題。另一方面,隨著雷諾數(shù)的不斷增加,二維模型將不再適用,需要開展三維數(shù)值計算。本文的計算結果也表明,即便是在簡單的二維低雷諾數(shù)條件下,約化速度的微小變化也會引起渦激振動響應特性的劇烈改變,這就要求采用很小的約化速度增量來開展細致的數(shù)值分析工作,這些都會導致計算量的激增,對現(xiàn)有的計算條件提出了較大的挑戰(zhàn)。此外,在高雷諾數(shù)條件下,直接數(shù)值模擬(Direct Numerical Simulation)在理論上是可靠的,但計算量巨大。采用湍流模型雖可緩解計算量的壓力,但湍流模式本身的適用性問題至今仍是全人類共同面對科學難題之一。

    參 考 文 獻

    [1]Khalak A,Williamson C H K.Motions,forces and mode transitions in vortex-induced vibrations at low mass-damping [J].Journal of Fluids and Structures,1999,13(7-8):813-851.

    [2]Feng C C.The measurement of vortex induced effects in flow past stationary and oscillating circular and D-section cylinders [D].Vancouver:University of British Columbia,1968.

    [3]Williamson C H K,Roshko A.Vortex formation in the wake of an oscillating cylinder [J].Journal of Fluids and Structures,1988,2(4):355-381.

    [4]何長江,段忠東.二維圓柱渦激振動的數(shù)值模擬 [J].海洋工程,2008,26(1):57-63.

    HE Chang-jiang,DUAN Zhong-dong.Numerical simulation of vortex-induced vibration on 2D circular cylinders [J].The Ocean Engineering,2008,26(1):57-63.

    [5]方平治,顧 明.圓柱兩自由度渦激振動的數(shù)值模擬研究 [J].同濟大學學報(自然科學版),2008,36(3):295-298.

    FANG Ping-zhi,GU Ming.Numerical simulation for vortex-induced vibration of circular cylinder with two degree of freedoms [J].Journal of Tongji University (Nature Science),2008,36(3):295-298.

    [6]唐國強.立管渦激振動數(shù)值模擬方法及物理模型實驗 [D].大連:大連理工大學,2011.

    [7]Brika D,Laneville A.The flow interaction between a stationary cylinder and a downstream flexible cylinder [J].Journal of Fluids and Structures,1999,13(5):579-606.

    [8]Carmo B S,Sherwin S J,Bearman P W,et al.Flow-induced vibration of a circular cylinder subjected to wake interference at low Reynolds number [J].Journal of Fluids and Structures,2011,27(4):503-522.

    [9]Bao Y,Zhou D,Tu J H.Flow interference between a stationary cylinder and an elastically mounted cylinder arranged in proximity [J].Journal of Fluids and Structures.2011,27(8):1425-1446.

    [10]Vandiver J K,Jong J Y.The relationship between in-line and cross-flow vortex-induced vibration of cylinders [J].Journal of Fluids and structures,1987,1(4):381-399.

    [11]陳文曲,任安祿,李廣望.串聯(lián)雙圓柱繞流下游圓柱兩自由度渦致振動研究 [J].力學學報,2004,36(6):732-738.

    CHEN Wen-qu,REN An-lu,LI Guang-wang.The Numerical study of two-degree-of -freedom vortex-induced vibration of the downstream cylinder in tandem arrangement [J].Chinese Journal of Theoretical and Applied Mechanics,2004,36(6):732-738.

    [12]Jiang C B,Kawahara M.The analysis of unsteady incompressible flows by a three-step finite element method [J].International Journal for Numerical Methods in Fluids,1993,16(9):793-811.

    [13]Lu L,Qin J M,Teng B,et al.Numerical investigations of life suppression by feedback rotary oscillation of circular cylinder at low Reynolds number [J].Physics of Fluids,2011,23(3):033601.

    [14]Lu L,Cheng L, Teng B,et al.Numerical investigation of fluid resonance in two narrow gaps of three identical rectangular structures [J].Applied Ocean Research,2010,32(2):177-190.

    [15]呂 林.海洋工程中小尺度物體的相關水動力數(shù)值計算 [D].大連:大連理工大學,2006.

    [16]Zdravkovich M M.The effects of interference between circular cylinders in cross flow [J].Journal of Fluids and Structures,1987,1(2):239-261.

    [17]Carmo B S,Meneghini J R,Sherwin S J.Secondary instabilities in the flow around two circular cylinders in tandem [J].Journal of Fluid Mechanics,2010,644:395-431.

    [18]Leontini J S,Stewart B E,Thompson M C,et al.Predicting vortex-induced vibration from driven oscillation results [J].Applied Mathematical Modelling,2006,30(10):1096-1102.

    猜你喜歡
    尾渦約化橫流
    不同B-V頻率下的飛機尾渦數(shù)值模擬研究
    約化的(3+1)維Hirota方程的呼吸波解、lump解和半有理解
    高空巡航階段的飛機尾渦流場演化特性研究
    橫流熱源塔換熱性能研究
    煤氣與熱力(2021年3期)2021-06-09 06:16:20
    基于激光雷達回波的動態(tài)尾渦特征參數(shù)計算
    干擾板作用下飛機尾渦流場近地演變機理研究
    基于橫流風扇技術的直升機反扭驗證
    脊下橫流對PEMFC性能影響的數(shù)值分析
    電源技術(2015年11期)2015-08-22 08:50:26
    M-強對稱環(huán)
    (3+1)-維廣義Kadomtsev-Petviashvili方程的對稱約化與精確解
    免费黄色在线免费观看| 国产白丝娇喘喷水9色精品| 欧美精品一区二区大全| 国产亚洲欧美精品永久| 日韩免费高清中文字幕av| 国产精品蜜桃在线观看| 免费大片黄手机在线观看| 亚洲精品视频女| 亚洲一级一片aⅴ在线观看| 卡戴珊不雅视频在线播放| 另类精品久久| 日韩,欧美,国产一区二区三区| 丝袜人妻中文字幕| 亚洲av电影在线观看一区二区三区| 在线 av 中文字幕| 欧美av亚洲av综合av国产av | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美国产精品一级二级三级| 春色校园在线视频观看| 亚洲精品一区蜜桃| 国产av码专区亚洲av| 免费播放大片免费观看视频在线观看| 毛片一级片免费看久久久久| 自线自在国产av| 免费在线观看完整版高清| 老鸭窝网址在线观看| 免费不卡的大黄色大毛片视频在线观看| 毛片一级片免费看久久久久| 久久久国产欧美日韩av| 久久久久久久久久人人人人人人| 精品国产乱码久久久久久小说| 国产成人av激情在线播放| 在线免费观看不下载黄p国产| 日韩电影二区| 日韩av不卡免费在线播放| 国产成人精品婷婷| 欧美日韩av久久| 亚洲国产看品久久| 视频在线观看一区二区三区| 两个人看的免费小视频| 99热全是精品| 久久久久精品人妻al黑| 亚洲内射少妇av| 好男人视频免费观看在线| 18禁观看日本| 最近中文字幕2019免费版| 精品一区二区免费观看| 不卡av一区二区三区| 日本91视频免费播放| 老司机影院毛片| 天天操日日干夜夜撸| 波多野结衣一区麻豆| 日本猛色少妇xxxxx猛交久久| 又粗又硬又长又爽又黄的视频| 色视频在线一区二区三区| 男的添女的下面高潮视频| 纵有疾风起免费观看全集完整版| 极品人妻少妇av视频| 久久热在线av| 下体分泌物呈黄色| 高清在线视频一区二区三区| 九草在线视频观看| 视频在线观看一区二区三区| 国产精品亚洲av一区麻豆 | xxx大片免费视频| 国产免费视频播放在线视频| 久热久热在线精品观看| 午夜福利,免费看| 亚洲激情五月婷婷啪啪| 大码成人一级视频| 啦啦啦视频在线资源免费观看| 精品少妇一区二区三区视频日本电影 | 精品久久久久久电影网| 久久久久久久久免费视频了| 热re99久久国产66热| a级毛片黄视频| 97精品久久久久久久久久精品| 男人舔女人的私密视频| 午夜福利乱码中文字幕| 日韩,欧美,国产一区二区三区| 男女免费视频国产| 成人黄色视频免费在线看| 各种免费的搞黄视频| 美女xxoo啪啪120秒动态图| 欧美日韩国产mv在线观看视频| 99热全是精品| 免费黄频网站在线观看国产| 高清视频免费观看一区二区| 丝袜美腿诱惑在线| 久久99蜜桃精品久久| 久久影院123| 久久精品久久久久久久性| 精品少妇黑人巨大在线播放| 性少妇av在线| 精品久久久精品久久久| 少妇精品久久久久久久| 国产精品成人在线| 久久 成人 亚洲| 乱人伦中国视频| 中文字幕人妻丝袜制服| 97在线视频观看| 国产精品国产三级国产专区5o| 国产精品不卡视频一区二区| 精品一区二区免费观看| 久久久久网色| 欧美激情 高清一区二区三区| 午夜av观看不卡| 久久97久久精品| 久久久久久人妻| 黑丝袜美女国产一区| 中国三级夫妇交换| 国产乱来视频区| 中文字幕人妻丝袜一区二区 | 黄色毛片三级朝国网站| 久久婷婷青草| 婷婷色综合大香蕉| 在现免费观看毛片| 欧美人与善性xxx| 欧美bdsm另类| 天天躁日日躁夜夜躁夜夜| 国产一区有黄有色的免费视频| 久久精品国产亚洲av高清一级| 免费久久久久久久精品成人欧美视频| 婷婷成人精品国产| 国产男女超爽视频在线观看| 亚洲av福利一区| 少妇人妻精品综合一区二区| 人人妻人人澡人人爽人人夜夜| 99久久中文字幕三级久久日本| a级毛片黄视频| 熟女av电影| 最新中文字幕久久久久| 国产精品久久久久久精品电影小说| 国产不卡av网站在线观看| 国产精品无大码| 美女福利国产在线| 国产淫语在线视频| 欧美精品一区二区免费开放| 午夜91福利影院| 欧美97在线视频| 汤姆久久久久久久影院中文字幕| 国产免费一区二区三区四区乱码| 天天躁狠狠躁夜夜躁狠狠躁| 中文欧美无线码| 一二三四中文在线观看免费高清| 在线观看免费视频网站a站| 日产精品乱码卡一卡2卡三| 久久精品久久精品一区二区三区| tube8黄色片| 亚洲精品日韩在线中文字幕| av电影中文网址| 人妻 亚洲 视频| 中文欧美无线码| 99九九在线精品视频| 国产淫语在线视频| 精品国产一区二区久久| 少妇 在线观看| 亚洲国产日韩一区二区| 亚洲精品国产一区二区精华液| 国产精品国产三级国产专区5o| 精品酒店卫生间| 成年女人毛片免费观看观看9 | 观看美女的网站| 国产深夜福利视频在线观看| 男女免费视频国产| 在线观看免费高清a一片| 欧美精品高潮呻吟av久久| 91午夜精品亚洲一区二区三区| 天美传媒精品一区二区| 久久精品熟女亚洲av麻豆精品| 亚洲男人天堂网一区| 久久99蜜桃精品久久| 色播在线永久视频| 日韩电影二区| 亚洲伊人久久精品综合| 精品国产超薄肉色丝袜足j| 欧美激情 高清一区二区三区| 国产白丝娇喘喷水9色精品| 国产老妇伦熟女老妇高清| 欧美成人午夜精品| 欧美精品高潮呻吟av久久| 人人妻人人爽人人添夜夜欢视频| 成人毛片a级毛片在线播放| 桃花免费在线播放| 卡戴珊不雅视频在线播放| 国产精品蜜桃在线观看| 深夜精品福利| 91国产中文字幕| 国产亚洲最大av| 欧美日韩av久久| 男女无遮挡免费网站观看| 亚洲第一区二区三区不卡| 亚洲av男天堂| 男男h啪啪无遮挡| 纯流量卡能插随身wifi吗| 夫妻性生交免费视频一级片| 夫妻午夜视频| 波多野结衣av一区二区av| 国产极品天堂在线| 不卡av一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 啦啦啦在线免费观看视频4| 国产成人精品一,二区| av在线观看视频网站免费| 嫩草影院入口| 精品人妻熟女毛片av久久网站| 日韩中文字幕视频在线看片| 国产精品国产三级国产专区5o| 久久久久久久久久人人人人人人| 黄色怎么调成土黄色| 久久国产精品男人的天堂亚洲| 国产一区二区三区av在线| 春色校园在线视频观看| videosex国产| 久久人人爽av亚洲精品天堂| 久久精品aⅴ一区二区三区四区 | 纯流量卡能插随身wifi吗| 老司机影院毛片| 一区二区三区精品91| 91成人精品电影| 午夜福利,免费看| 大香蕉久久成人网| 亚洲av电影在线观看一区二区三区| 少妇 在线观看| 精品人妻一区二区三区麻豆| 看免费成人av毛片| 欧美精品av麻豆av| 在线观看一区二区三区激情| 最新中文字幕久久久久| 青春草国产在线视频| 99精国产麻豆久久婷婷| 高清不卡的av网站| 久久久欧美国产精品| 少妇人妻久久综合中文| 国产一级毛片在线| 超色免费av| 夜夜骑夜夜射夜夜干| 国产乱人偷精品视频| 视频在线观看一区二区三区| 国产老妇伦熟女老妇高清| 伦理电影大哥的女人| 男女高潮啪啪啪动态图| 丝袜喷水一区| 香蕉国产在线看| 国产老妇伦熟女老妇高清| 女性生殖器流出的白浆| 亚洲国产精品一区三区| 赤兔流量卡办理| 亚洲av.av天堂| 久久久精品区二区三区| 精品福利永久在线观看| 婷婷色综合www| 久久精品国产亚洲av涩爱| 午夜激情久久久久久久| 97在线人人人人妻| 亚洲av成人精品一二三区| 一级黄片播放器| 午夜福利在线观看免费完整高清在| 亚洲色图综合在线观看| 亚洲国产精品999| 精品久久蜜臀av无| 国产乱人偷精品视频| 亚洲av欧美aⅴ国产| 久久午夜综合久久蜜桃| 国产精品蜜桃在线观看| 亚洲四区av| 老熟女久久久| 麻豆av在线久日| 可以免费在线观看a视频的电影网站 | 母亲3免费完整高清在线观看 | 中文字幕亚洲精品专区| 男女下面插进去视频免费观看| 久久久久久久久免费视频了| 伊人久久国产一区二区| 纵有疾风起免费观看全集完整版| 免费观看性生交大片5| 国产欧美日韩一区二区三区在线| 寂寞人妻少妇视频99o| 只有这里有精品99| 热re99久久国产66热| 欧美变态另类bdsm刘玥| 99热国产这里只有精品6| a级毛片黄视频| 久久久国产精品麻豆| 各种免费的搞黄视频| 最近2019中文字幕mv第一页| 免费在线观看完整版高清| 国产在线免费精品| 亚洲国产av影院在线观看| 在现免费观看毛片| 欧美黄色片欧美黄色片| 一本—道久久a久久精品蜜桃钙片| 校园人妻丝袜中文字幕| 性色av一级| 成年人免费黄色播放视频| 香蕉国产在线看| 99久国产av精品国产电影| kizo精华| 成人国产麻豆网| 在线 av 中文字幕| 在线天堂中文资源库| 国产97色在线日韩免费| 国产免费一区二区三区四区乱码| 自线自在国产av| 日本欧美视频一区| 亚洲激情五月婷婷啪啪| www.熟女人妻精品国产| 日韩 亚洲 欧美在线| 国产成人精品福利久久| 一级毛片 在线播放| 一区二区三区激情视频| 肉色欧美久久久久久久蜜桃| 欧美精品av麻豆av| 在线观看一区二区三区激情| 国产精品一国产av| av女优亚洲男人天堂| 岛国毛片在线播放| 777米奇影视久久| 国产日韩欧美视频二区| 国产精品.久久久| 有码 亚洲区| 18禁动态无遮挡网站| 成人毛片60女人毛片免费| 国产一区二区 视频在线| 精品亚洲成a人片在线观看| 久久99精品国语久久久| 丰满少妇做爰视频| 久久毛片免费看一区二区三区| 国产一区二区在线观看av| 黄色视频在线播放观看不卡| 久久鲁丝午夜福利片| 国产男人的电影天堂91| 老汉色∧v一级毛片| 永久免费av网站大全| 看免费av毛片| 国产精品久久久av美女十八| 中文欧美无线码| 99久国产av精品国产电影| 十八禁网站网址无遮挡| 又粗又硬又长又爽又黄的视频| 老汉色∧v一级毛片| 亚洲欧美精品自产自拍| 成人漫画全彩无遮挡| 精品亚洲成a人片在线观看| 黄片无遮挡物在线观看| 夜夜骑夜夜射夜夜干| 久久av网站| 日韩大片免费观看网站| 黄网站色视频无遮挡免费观看| 国产男女超爽视频在线观看| 久久久久人妻精品一区果冻| 高清在线视频一区二区三区| 欧美日韩视频精品一区| 亚洲欧美清纯卡通| 久久国产精品大桥未久av| 丝袜美足系列| 黄色一级大片看看| 卡戴珊不雅视频在线播放| 精品酒店卫生间| 啦啦啦视频在线资源免费观看| 国产欧美亚洲国产| 色播在线永久视频| 你懂的网址亚洲精品在线观看| 好男人视频免费观看在线| 国产精品一区二区在线观看99| 啦啦啦在线免费观看视频4| 看免费成人av毛片| 好男人视频免费观看在线| 波野结衣二区三区在线| 午夜老司机福利剧场| 国产av码专区亚洲av| av电影中文网址| 肉色欧美久久久久久久蜜桃| 亚洲国产色片| 最新的欧美精品一区二区| 人妻少妇偷人精品九色| 欧美日韩综合久久久久久| 精品国产一区二区久久| 日本wwww免费看| 亚洲欧美清纯卡通| 欧美激情极品国产一区二区三区| 男女边吃奶边做爰视频| 久久久国产欧美日韩av| 两性夫妻黄色片| 人妻人人澡人人爽人人| 九草在线视频观看| 中文精品一卡2卡3卡4更新| 热re99久久精品国产66热6| 免费观看在线日韩| 啦啦啦在线免费观看视频4| 人妻 亚洲 视频| 99热网站在线观看| 国产精品久久久久久精品古装| 成人毛片a级毛片在线播放| 99九九在线精品视频| 久久久久精品人妻al黑| 亚洲精品国产色婷婷电影| 午夜影院在线不卡| 国产极品粉嫩免费观看在线| 亚洲五月色婷婷综合| 国产av码专区亚洲av| 最近中文字幕高清免费大全6| 精品视频人人做人人爽| 亚洲国产av新网站| 黄网站色视频无遮挡免费观看| 亚洲欧美精品自产自拍| 久久韩国三级中文字幕| 熟女av电影| 又粗又硬又长又爽又黄的视频| 韩国精品一区二区三区| 久久精品亚洲av国产电影网| 国产片特级美女逼逼视频| 美女高潮到喷水免费观看| 国产片内射在线| kizo精华| 免费观看无遮挡的男女| 你懂的网址亚洲精品在线观看| 五月天丁香电影| 天堂中文最新版在线下载| 啦啦啦在线免费观看视频4| 宅男免费午夜| 国产欧美亚洲国产| 男女啪啪激烈高潮av片| 国产1区2区3区精品| 久久婷婷青草| 菩萨蛮人人尽说江南好唐韦庄| 一级黄片播放器| 久久精品国产亚洲av高清一级| 午夜日韩欧美国产| 欧美+日韩+精品| √禁漫天堂资源中文www| 国产精品麻豆人妻色哟哟久久| 少妇被粗大猛烈的视频| 日韩一区二区三区影片| 成人毛片a级毛片在线播放| 亚洲精品国产av成人精品| 伦理电影大哥的女人| 亚洲情色 制服丝袜| 国产淫语在线视频| 各种免费的搞黄视频| 亚洲精品一区蜜桃| 最黄视频免费看| 天堂俺去俺来也www色官网| 日韩av不卡免费在线播放| 老司机影院成人| 成年美女黄网站色视频大全免费| 亚洲欧美一区二区三区久久| 99re6热这里在线精品视频| 伦理电影免费视频| 一区福利在线观看| 成年女人毛片免费观看观看9 | 免费观看在线日韩| 国产精品免费视频内射| 国产精品无大码| 国产欧美日韩一区二区三区在线| 亚洲精品乱久久久久久| kizo精华| 热re99久久国产66热| 国产成人精品久久久久久| 亚洲精品成人av观看孕妇| av免费观看日本| 男女午夜视频在线观看| av一本久久久久| 久久鲁丝午夜福利片| 一本久久精品| 亚洲激情五月婷婷啪啪| 精品99又大又爽又粗少妇毛片| 国产精品亚洲av一区麻豆 | 亚洲精品,欧美精品| 国产精品久久久久久精品古装| 精品一区二区三区四区五区乱码 | 午夜激情av网站| 母亲3免费完整高清在线观看 | 我的亚洲天堂| 亚洲成人一二三区av| 最新的欧美精品一区二区| 国产精品久久久久久精品电影小说| 亚洲经典国产精华液单| 青春草亚洲视频在线观看| 人妻 亚洲 视频| 亚洲av.av天堂| 精品99又大又爽又粗少妇毛片| 看非洲黑人一级黄片| 中文字幕制服av| 大香蕉久久成人网| 97精品久久久久久久久久精品| 成人毛片a级毛片在线播放| 午夜激情av网站| 老鸭窝网址在线观看| a 毛片基地| 亚洲激情五月婷婷啪啪| 午夜激情av网站| 母亲3免费完整高清在线观看 | 美女午夜性视频免费| 欧美日韩亚洲国产一区二区在线观看 | 老熟女久久久| 成人国产麻豆网| 日韩制服丝袜自拍偷拍| 看非洲黑人一级黄片| 三上悠亚av全集在线观看| 国产 一区精品| 亚洲精品在线美女| 香蕉国产在线看| 婷婷色av中文字幕| 国产精品一国产av| 飞空精品影院首页| 国产国语露脸激情在线看| 日韩制服丝袜自拍偷拍| 极品人妻少妇av视频| 国产亚洲av片在线观看秒播厂| av片东京热男人的天堂| 国产在线一区二区三区精| 视频区图区小说| 国产爽快片一区二区三区| 国产亚洲欧美精品永久| 老汉色av国产亚洲站长工具| 国产精品熟女久久久久浪| 国产精品一国产av| 国产又爽黄色视频| 日本爱情动作片www.在线观看| 你懂的网址亚洲精品在线观看| 国产乱来视频区| 美女国产视频在线观看| 热99久久久久精品小说推荐| 9191精品国产免费久久| 日韩av不卡免费在线播放| 九九爱精品视频在线观看| 一级片'在线观看视频| 精品少妇黑人巨大在线播放| 老司机亚洲免费影院| 日本色播在线视频| 国产淫语在线视频| 国产1区2区3区精品| 美女大奶头黄色视频| 国产爽快片一区二区三区| 色哟哟·www| 日韩熟女老妇一区二区性免费视频| 人妻少妇偷人精品九色| 精品国产一区二区久久| 热99久久久久精品小说推荐| 亚洲国产最新在线播放| 亚洲精品,欧美精品| 丝袜喷水一区| 亚洲精品久久午夜乱码| 在线看a的网站| 国产成人欧美| 免费久久久久久久精品成人欧美视频| 女性生殖器流出的白浆| 侵犯人妻中文字幕一二三四区| 国产一区二区在线观看av| 亚洲精华国产精华液的使用体验| 99九九在线精品视频| 欧美亚洲 丝袜 人妻 在线| 日日啪夜夜爽| 亚洲情色 制服丝袜| 亚洲第一区二区三区不卡| 香蕉精品网在线| 国产亚洲欧美精品永久| 国产白丝娇喘喷水9色精品| 少妇人妻精品综合一区二区| av片东京热男人的天堂| 国产精品二区激情视频| 国产精品熟女久久久久浪| 99香蕉大伊视频| h视频一区二区三区| 少妇猛男粗大的猛烈进出视频| 亚洲精品自拍成人| 极品少妇高潮喷水抽搐| 国产综合精华液| 永久免费av网站大全| www.自偷自拍.com| 中国国产av一级| 永久网站在线| 天天躁夜夜躁狠狠躁躁| 国产高清国产精品国产三级| 欧美激情高清一区二区三区 | 国产成人精品福利久久| 久久精品夜色国产| 日韩大片免费观看网站| 欧美日本中文国产一区发布| 美女国产高潮福利片在线看| 亚洲美女视频黄频| 国产伦理片在线播放av一区| 高清欧美精品videossex| 69精品国产乱码久久久| 少妇被粗大的猛进出69影院| 日韩免费高清中文字幕av| 观看av在线不卡| 少妇的逼水好多| 国产精品久久久久成人av| 亚洲美女黄色视频免费看| 日韩精品免费视频一区二区三区| 久久午夜福利片| 国产精品一国产av| 久久久精品免费免费高清| 亚洲综合精品二区| 男人操女人黄网站| 国产av一区二区精品久久| 亚洲欧美清纯卡通| 精品国产乱码久久久久久小说| 国产又色又爽无遮挡免| 婷婷色综合www| 免费在线观看黄色视频的| 1024香蕉在线观看| 女性被躁到高潮视频| 国产一区二区在线观看av| 99国产精品免费福利视频| 日韩 亚洲 欧美在线| 一本—道久久a久久精品蜜桃钙片| 九草在线视频观看| 国产精品国产三级国产专区5o| 又粗又硬又长又爽又黄的视频| 亚洲人成77777在线视频| 国产一级毛片在线| 国产视频首页在线观看|