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

    附加整流裝置的圓柱體渦激振動數(shù)值研究

    2022-01-27 14:14:58姚宗鍇韓曉雙
    振動與沖擊 2022年1期
    關(guān)鍵詞:整流罩方根值渦激

    周 波, 王 杰, 姚宗鍇, 劉 輝, 韓曉雙

    (1. 大連理工大學(xué) 工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室 運載工程與力學(xué)學(xué)部船舶工程學(xué)院, 遼寧 大連 116024;2. 港口航道泥沙工程交通行業(yè)重點實驗室, 南京 210024; 3. 大連海事大學(xué) 輪機(jī)工程學(xué)院, 遼寧 大連 116026)

    在一定來流作用下,圓柱體兩側(cè)交替瀉放的漩渦引起的渦激振動(vortex-induced vibration, VIV)現(xiàn)象是導(dǎo)致結(jié)構(gòu)疲勞損傷的主要原因。為減輕VIV對結(jié)構(gòu)造成損害,延長使用壽命,國內(nèi)外學(xué)者在相關(guān)領(lǐng)域進(jìn)行了大量研究[1-3]。常見的VIV抑制裝置可分為兩大類:主動控制和被動控制。主動控制利用計算機(jī)技術(shù)對流場進(jìn)行實時監(jiān)控,通過引入外部擾動來調(diào)節(jié)漩渦脫落。被動控制則是通過改變結(jié)構(gòu)截面形狀或附加額外的裝置來改變結(jié)構(gòu)周圍流場,達(dá)到抑制漩渦的形成和發(fā)展[4]。由于被動控制操作簡單,效率高且成本低,在工程領(lǐng)域得到廣泛應(yīng)用,常見的被動控制裝置如圖1所示[5]。

    (a) 整流罩

    (b) 分隔板圖1 被動控制裝置Fig.1 Passive control device

    譚波等[6]計算了亞臨界雷諾數(shù)范圍內(nèi)附加分隔板的靜止圓柱體繞流。結(jié)果表明,附加不同尺寸分隔板時圓柱體抑制效果存在顯著差異。王嘉松等[7]設(shè)計了一款可旋轉(zhuǎn)整流罩,克服了傳統(tǒng)渦激振動抑制裝置對于來流方向的限制。盛磊祥等[8]利用LES計算了Re=2.5×105時,附加螺旋列板的隔水管渦激振動。結(jié)果表明安裝螺旋列板能夠顯著減小立管所受的渦激升力,但是由于螺旋列板后方壓力驟降,平均阻力顯著增加。水慶象等[9]研究帶橫隔板圓柱體繞流時發(fā)現(xiàn),橫隔板能有效的抑制尾渦的脫落,提高圓柱體背面的壓力,減小圓柱體上下表面的壓力差。Sahu等[10]研究了附加剛性分隔板的圓柱體流激振動時發(fā)現(xiàn):在渦激振動區(qū)間內(nèi),圓柱體的響應(yīng)幅值和鎖定區(qū)間隨著分隔板長度的增加而增加。與此同時,在渦激振動區(qū)間和幅值激增區(qū)間內(nèi),隨著分隔板的增加,圓柱體的振動頻率減小。

    現(xiàn)有研究表明,附加整流裝置能夠有效改善靜止圓柱體尾渦脫落和受力情況。但實際工程情況下,圓柱體并非固定不動。本文主要研究了附加整流裝置的圓柱體橫向單自由度渦激振動。重點分析了3個代表性頂角θ以及分隔板長度Ls/D時,圓柱體渦激振動特性。整流罩頂角θ和分隔板長度Ls具體定義如圖2所示。

    圖2 整流罩頂角與分隔板長度定義Fig.2 Definition of the top angle of fairing and the length of splitter

    1 數(shù)值計算

    1.1 控制方程

    二維不可壓縮流體流控制方程通常包括連續(xù)性方程和動量方程,具體表達(dá)式為

    ?·u=0

    (1)

    ?u/?t+(u·?)u=-1/ρ?p+ν?2u

    (2)

    式中:u為來流速度;p為壓力;ρ為流體密度;ν為流體運動黏性系數(shù)。本文采用FVM對計算區(qū)域進(jìn)行離散;采用SIMPLE算法求解壓力-速度的耦合;壓力項采用二階格式進(jìn)行離散;對流向采用二階迎風(fēng)格式;擴(kuò)散項采用中心差分格式;時間項采用一階隱式離散。各參數(shù)的殘差設(shè)置為1×10-5;時間步長為Δt=0.005 s。

    圓柱體橫向單自由度渦激振動運動方程為

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    K3)/6

    (9)

    1.2 流域網(wǎng)格劃分

    計算流域如圖3所示,圓柱體的直徑為D,坐標(biāo)原點位于圓柱中心,x軸正方向指向下游,y軸垂直于來流方向;計算流域長度為42D,流場寬度24D,圓柱中心距入流口和出流口分別為12D和30D,圓柱的中心距上下邊界均為12D,阻塞率為4.5%。研究表明,阻塞率小于5%時,流場寬度對圓柱體響應(yīng)的影響就可以忽略了。為了更好地捕捉圓柱體周圍及尾部流動特征,將流域劃分為3個部分:隨動區(qū)域,變形區(qū)域和靜止區(qū)域。圓柱體周圍4D范圍內(nèi)的網(wǎng)格隨圓柱體一起運動,隨動區(qū)域內(nèi)的網(wǎng)格不隨圓柱體的運動發(fā)生變化;變形區(qū)域內(nèi)的網(wǎng)格在每個時間步內(nèi),根據(jù)動網(wǎng)格守恒進(jìn)行更新;靜止區(qū)域內(nèi)的網(wǎng)格在計算過程中保持不變。

    圖3 計算流域劃分示意圖Fig.3 Schematic of the computational zone division

    1.3 邊界條件設(shè)置

    (1) 入口邊界條件:來流速度設(shè)為常數(shù),即u=ux,v=0;(2) 出口邊界條件:為充分發(fā)展流出邊界條件,?u/?x=0,?v/?x=0;(3) 上下邊界條件:為對稱邊界條件,?u/?y=0,v=0;(4) 無滑移壁面邊界條件:壁面上的速度為零,即u=0,v=0。

    2 網(wǎng)格獨立性驗證和數(shù)值計算結(jié)果

    2.1 網(wǎng)格獨立性驗證

    網(wǎng)格質(zhì)量和數(shù)量會很大程度上影響計算精度和效率。首先計算了Re=200時靜止圓柱體繞流,計算結(jié)果如表1所示。最大相對誤差出現(xiàn)在升力幅值為0.9%,兼顧計算效率和精度,選擇Mesh2作為計算網(wǎng)格。此外,本文計算結(jié)果與文獻(xiàn)[11-12]結(jié)果吻合度較高,說明本文計算網(wǎng)格的可靠性。

    表1 雷諾數(shù)Re=200時,不同網(wǎng)格圓柱體繞流計算結(jié)果

    2.2 附加整流裝置圓柱體靜止繞流

    計算了Re=150時,附加不同幾何參數(shù)整流裝置的圓柱體靜止繞流,計算結(jié)果如表2所示。

    表2 雷諾數(shù)Re=150時,附加整流裝置的圓柱體計算結(jié)果

    對于附加整流裝置的圓柱體來說,其尾部凸起能夠有效的抑制剪切層相互作用,從而減小圓柱體的受力。對于附加整流罩的圓柱體來說,隨著θ增加,圓柱體受力增大。隨著θ增加,整流罩尾部抑制剪切層相互作用的凸起越來越小。當(dāng)θ=90°時,其幾何外形和單圓柱體十分接近,圓柱體的受力、壓力分布、尾渦脫落等和單圓柱體相近。圓柱體受力的最小值出現(xiàn)在θ=45°時,其升力均方根值減小55%,阻力均值減小17%。對于附加分隔板的圓柱體來說,隨著Ls增加,圓柱體的受力減小。隨著Ls的增加,圓柱體尾部瀉渦中心偏離主圓柱體,對圓柱體的影響減弱。圓柱體受力的最小值出現(xiàn)在Ls/D=0.5時,其升力均方根值減小27%,阻力均值減小13%。此外,加裝整流裝置的圓柱體St較單圓柱體減小,表明圓柱體的瀉渦頻率減小,漩渦脫落周期變大,說明加裝整流裝置能夠有效地抑制圓柱體尾渦的脫落。

    圖4為附加整流裝置的靜止圓柱體繞流壓力云圖。為了方便對比,將單圓柱體計算結(jié)果加入。附加整流罩的靜止圓柱繞流壓力云圖,整流罩頂角θ依次為:45°,60°和90°;分隔板長度Ls/D依次為0.25,0.35和0.5。

    圖4 附加整流裝置的固定圓柱體周圍壓力云圖Fig.4 Pressure field in proximity of a fixed cylinder with additional rectifiers

    從圖4壓力云圖中可以看出,隨著θ減小,整流罩上下兩側(cè)壓力分布更加均勻,由單圓柱體的低壓集中在一側(cè)擴(kuò)散到兩側(cè),導(dǎo)致圓柱體所受的升力均方根值減??;隨著θ增加,圓柱體尾部的低壓中心向上游方向移動,即負(fù)壓中心越來越靠近圓柱體后部,圓柱體前后兩端的壓力差增大,從而圓柱體所受的阻力增大。

    對于附加分隔板的圓柱體來說,隨著Ls增加,圓柱體上側(cè)的低壓區(qū)域越大,上下兩側(cè)的壓力差減小,導(dǎo)致圓柱體所受的渦激升力減??;隨著Ls增加,圓柱體尾部的低壓中心發(fā)生后移,前后兩端的壓力差減小,導(dǎo)致圓柱體所受的阻力減小。

    2.3 附加整流罩的圓柱體橫向單自由度渦激振動

    首先計算了Re=150時,彈性支撐的單圓柱體橫向單自由度渦激振動。質(zhì)量比m*=4m/ρfπD2=2;結(jié)構(gòu)阻尼系數(shù)ζ=0.007。整個計算過程中來流速度保持不變,采用改變fn來改變約化速度Ur=U/fnD。

    圖5 附加不同頂角整流罩的圓柱體動振動幅值(左)和無因次化振動頻率(右)變化曲線Fig.5 Variation of the amplitude (left) and dimensionless vibration frequency (right) of a cylinder with fairing at different top angles

    附加整流罩的圓柱體振動頻率隨著θ增加而增加,但始終小于單圓柱體的振動頻率,說明附加整流罩的圓柱體在一定程度上能夠減輕結(jié)構(gòu)的振動頻率。當(dāng)θ=45°和60°時,圓柱體的振動頻率隨著Ur變化平緩。當(dāng)θ=90°時,隨著Ur增加鎖定到圓柱體固有頻率。當(dāng)Ur≥10時,脫離鎖定區(qū)間,重新鎖定到瀉渦頻率上。

    圖6為附加不同頂角整流罩的圓柱體升力均方根值(左)和阻力均值(右)變化曲線。對于單圓柱來說,隨著Ur增加,升力系數(shù)的均方根值在Ur=4達(dá)到最大值后急劇下降。在Ur=6時,升力均方根值達(dá)到最小值。之后,隨著Ur增加,升力均方根值緩慢增加并趨于平緩。當(dāng)θ=45°時,升力均方根值隨著Ur增加急劇增加,在Ur=6時達(dá)到最大值,隨后緩慢減小。當(dāng)θ=60°時,隨著Ur增加,升力均方根值顯著增加,在Ur=5達(dá)到最大值,隨后逐漸減小。當(dāng)θ=90°時,升力均方根值變化趨勢和單圓柱體類似。從圖7,Ur=4壓力云圖中可以看出,θ=45°時圓柱體上下兩側(cè)壓力較單圓柱體分布更加均勻,所受的渦激升力小于單圓柱體。當(dāng)θ=60°和90°時,圓柱體下側(cè)的低壓逐漸減小,低壓主要集中在圓柱體的上側(cè),且隨著θ增加,低壓值逐漸增大,而整流罩下側(cè)壓力變化不大,使得圓柱體上下兩側(cè)的壓力差逐漸增大。從圖7,Ur=8壓力云圖中可以看出,隨著θ增加,圓柱體上側(cè)低壓值減小,右下側(cè)低壓值增大,上下兩側(cè)的壓力差迅速減小,導(dǎo)致圓柱體所受的渦激升力減小。

    圖6 附加不同頂角整流罩的圓柱體升力均方根值(左)和阻力均值(右)變化曲線Fig.6 Variation of the r.m.s. lift coefficient (left) and mean drag coefficient (right) of a cylinder with fairing at different top angles

    圖7 約化速度Ur=4(上)和8(下)時,附加整流罩的圓柱體橫向單自由度渦激振動壓力云圖Fig.7 Pressure field in proximity of a cylinder undergoing transverse-only VIV with a fairing at Ur=4 (up) and 8 (down)

    對于附加整流罩圓柱平均阻力來說,在Ur較小時,帶整流罩的圓柱體所受的阻力均值較單圓柱體減?。恢械萓r時,阻力均值較單圓柱體增大;Ur較大時,阻力均值和單圓柱體相近。同時,隨著θ增加,圓柱體受力的最大值向更小的Ur偏移。從圖7,Ur=4壓力云圖可以看出,當(dāng)θ=45°時,盡管其前端高壓駐點壓力值大于θ=60°和90°,但是其尾部壓力值遠(yuǎn)大于θ=60°和90°,故使得圓柱體前后端壓力差小于θ=60°和90°。當(dāng)θ=60°和90°時,前端高壓值相近。θ=60°時,尾部壓力比θ=90°大,故θ=60°時圓柱體所受的平均阻力小于θ=90°。從圖7,Ur=8壓力云圖可以看出,隨著θ增加,圓柱體前端高壓駐點壓力和尾部壓力越來越大,前端增加的較尾部大,使得壓力差越來越大,平均阻力越來越大。

    圖8為約化速度Ur=4(上)和8(下)時,附加整流罩的圓柱體橫向單自由度渦激振動瞬時渦量云圖。當(dāng)Ur=4時,單圓柱體表現(xiàn)為緊湊的雙列線渦。當(dāng)θ=45°時,其尾渦結(jié)構(gòu)表現(xiàn)為2S模式。與單圓柱體不同,漩渦從整流罩的尾部脫落,漩渦脫落遠(yuǎn)離主圓柱體,主圓柱體遠(yuǎn)離漩渦的影響。由于振幅較小,整流罩對于瀉渦的影響比較弱。當(dāng)圓柱體幅值為正時,整流罩的下表面形成一個逆時針的次渦,并逐漸擴(kuò)散到整個整流罩尾部;當(dāng)圓柱體幅值為負(fù)時,整流罩的上表面形成一個順時針的次渦,并逐漸擴(kuò)散到整個整流罩尾部。次渦的存在使得圓柱體上下表面的壓力差減小,導(dǎo)致圓柱體所受的渦激升力減小。當(dāng)θ=60°時,圓柱體的橫向振幅略有增加,整流罩尾部對漩渦的脫落產(chǎn)生明顯的影響。在整流罩的作用下,漩渦被拉長,相鄰的漩渦中心距離變長。當(dāng)θ=90°時,其幾何外形和單圓柱體相近,故其尾渦結(jié)構(gòu)和單圓柱體相近表現(xiàn)為緊湊的雙列線渦且一直延續(xù)到遠(yuǎn)場。

    當(dāng)Ur=8時,對于單圓柱體來說,振幅迅速減小,第一個漩渦中心較約化速度低時向下游偏移,其尾渦表現(xiàn)為2S模式。當(dāng)θ=45°時,在整流罩的作用下,圓柱體的尾渦被拉長,圓柱體兩側(cè)的漩渦能量較Ur小時增強(qiáng),使得圓柱體的振幅增大。當(dāng)θ=60°時,整流罩對于尾流的影響比θ=45°小,尾渦被拉長,每個周期內(nèi)脫落兩個旋向相反的漩渦形成一個漩渦對。當(dāng)圓柱體向上運動時,漩渦對偏向圓柱體下側(cè);當(dāng)圓柱體向下運動時,漩渦對偏向圓柱體上側(cè)。當(dāng)θ=90°時,圓柱體后的尾渦被拉長,瀉渦模式表現(xiàn)為2S模式。

    圖8 約化速度Ur=4(上)和8(下)時,附加整流罩的圓柱體橫向單自由度渦激振動瞬時渦量云圖Fig.8 Instantaneous vorticity of a cylinder undergoing transverse-only VIV with a fairing at Ur=4 (up) and 8 (down)

    2.4 附加分隔板的圓柱體橫向單自由度渦激振動

    附加不同長度分隔板的圓柱體振動振動幅值(左)和無因次化振動頻率(右)變化曲線如圖9所示。當(dāng)Ls/D=0.25時,圓柱體橫向振動幅值隨著Ur增加先增大后減小,在Ur=7時取得最大值0.857,比單圓柱體的大約50%。當(dāng)Ls/D=0.35時,圓柱體橫向振動幅值變化趨勢和Ls/D=0.25相近,在Ur=10時取得最大值。當(dāng)Ls/D=0.5時,圓柱體的橫向振動位移隨著約化速度的增加而增加??偟膩碚f,附加分隔板的圓柱體振幅在約化速度較小時,其橫向振幅小于單圓柱體;在約化速度較大時,圓柱體的橫向振幅遠(yuǎn)大于單圓柱體,且隨著Ls增加,圓柱體的振幅越來越大,最大橫向位移向更高的約化速度偏移。附加分隔板的圓柱體振動頻率整體上隨著Ur增加緩慢增加,且隨著Ls增加,振動頻率越來越小,說明附加分隔板的圓柱體在一定程度上能夠減輕結(jié)構(gòu)的振動頻率。

    圖10為附加不同長度分隔板的圓柱體渦激振動升力均方根值(左)和阻力均值(右)變化曲線。附加分隔板的圓柱體升力均方根值總的來說比單圓柱體大。當(dāng)Ur<5時,隨著Ls增加,圓柱體所受的升力均方根值越來越?。划?dāng)Ur>5時,隨著Ls增加,圓柱體所受的升力均方根值越來越大。從圖11,Ur=4的壓力云圖可以看出,圓柱體運動到最大振幅時,低壓中心出現(xiàn)在圓柱體的上側(cè),且隨著Ls增加低壓絕對值減小,而圓柱體的下側(cè)壓力值幾乎相等,圓柱體上下兩側(cè)的壓力差減小,導(dǎo)致圓柱體所受的渦激升力減小。從圖11,Ur=8的壓力云圖可以看出,對于單圓柱體,上下兩側(cè)的壓力差小于附加分隔板的圓柱體。附加分隔板的圓柱體低壓中心主要集中在圓柱體的上側(cè)和右下側(cè),隨著Ls增加,右下側(cè)的壓力值增大,上下兩側(cè)的壓力差增大,導(dǎo)致圓柱體所受的渦激升力增大。

    圖9 附加不同長度分隔板的圓柱體振動振動幅值(左)和無因次化振動頻率(右)變化曲線Fig.9 Variation of the amplitude (left) and non-dimensional vibration frequency (right) of a cylinder with splitter in different lengths

    圖10 附加不同長度分隔板的圓柱體渦激振動升力均方根值(上)和阻力均值(下)變化曲線Fig.10 Variation of the r.m.s. lift coefficient (up) and mean drag coefficient (down) of a cylinder with splitter in different lengths

    對于平均阻力來說,當(dāng)約化速度較小時,阻力均值遠(yuǎn)小于單圓柱體;中等約化速度時,阻力均值較單圓柱體增大;約化速度較大時,阻力均值和單圓柱體相近。從圖11,Ur=4的壓力云圖可以看出,對于單圓柱體,其后部低壓區(qū)域使得前后壓力差遠(yuǎn)大于附加分隔板的圓柱體。附加分隔板的圓柱體運動到最大幅值時,圓柱體前端高壓駐點由圓柱體正前方偏移到圓柱體的左下方,圓柱體尾部的壓力值隨著Ls增加逐漸增加,圓柱體前后壓力差越來越小,導(dǎo)致圓柱體所受的平均阻力越來越小。從圖11,Ur=8的壓力云圖可以看出,對于單圓柱體,其前端和尾部壓力值均大于附加分隔板的圓柱體,但是前端高壓幅度小于后端低壓大的值,反而使得圓柱體前后壓力差小于附加分隔板的圓柱體。附加分隔板的圓柱體高壓中心回到正前方附近,圓柱體前后兩端的壓力差幾乎相近,故圓柱體所受的阻力均值很接近。

    圖11 約化速度Ur=4(上)和8(下)時,附加分隔板的圓柱體橫向單自由度渦激振動壓力云圖Fig.11 Pressure field in proximity of a cylinder undergoing transverse-only VIV with a splitter at Ur=4 (up) and 8 (down)

    從圖12,Ur=4瞬時渦量云圖中可以看出,圓柱體幅值為正時,在分隔板的下表面形成一個逆時針的次渦,并逐漸擴(kuò)張到整個圓柱體后半部分;當(dāng)圓柱體幅值為負(fù)時,在分隔板的上表面形成一個順時針的次渦,并逐漸擴(kuò)張到整個圓柱體的后半部分。隨著Ls增加,分隔板末端的漩渦越來越強(qiáng),使得圓柱體上下兩側(cè)的壓力差越小,圓柱體所受的渦激升力越來越小,圓柱體振動幅值越來越小。圓柱體的瀉渦模式表現(xiàn)為2S模式,且隨著Ls增加,圓柱體后的尾渦中心距離圓柱體中心越遠(yuǎn),相鄰的兩個漩渦之間的距離增加,漩渦的脫落周期增大。從圖12,Ur=8瞬時渦量云圖中可以看出,Ls/D=0.25時,近圓柱體的瀉渦模式表現(xiàn)為2Po模式。隨著時間的推移向下游移動,能量較弱的逆時針漩渦消散,形成P+S模式。當(dāng)Ls/D=0.35時,圓柱體由位移最大處向下運動時,兩個順時針漩渦和一個逆時針漩渦脫落;圓柱體向上移動時,在圓柱體的尾部脫落一個逆時針的渦帶貫穿整個渦道。當(dāng)Ls/D=0.5時,圓柱體的瀉放表現(xiàn)為多個周期,每個周期內(nèi)的漩渦脫落較分隔板長度較短時變的復(fù)雜,尾渦結(jié)構(gòu)變的混亂。

    圖12 約化速度Ur=4(上)和8(下)時,附加分隔板的圓柱體橫向單自由度渦激振動瞬時渦量云圖Fig.12 Instantaneous vorticity of a cylinder undergoing transverse-only VIV with a splitter at Ur=4 (up) and 8 (down)

    3 結(jié) 論

    本文通過數(shù)值模擬方法對附加整流裝置圓柱的流場特性進(jìn)行了研究,得出如下結(jié)論:

    (1) 附加整流裝置的圓柱體橫向單自由度渦激振動特性比附加整流裝置的固定圓柱體復(fù)雜的多,使用靜止圓柱體繞流結(jié)果來表征整流裝置好壞有些欠妥。

    (2) 附加整流裝置的圓柱體振動頻率較單圓柱體顯著減小,說明附加整流裝置能夠減輕結(jié)構(gòu)的振動頻率。

    (4) 整流裝置尾部形成的次渦對結(jié)構(gòu)受力影響較大。當(dāng)Ur較低時,次渦的存在使得圓柱體上下表面的壓力差減小,導(dǎo)致圓柱體所受的渦激升力減小。當(dāng)Ur較大時,在整流裝置和次渦共同作用下,圓柱體的所受的渦激升力增強(qiáng)。此外,圓柱體尾部低壓中心發(fā)生明顯后移,圓柱體前后兩端壓力差減小,最終導(dǎo)致圓柱體所受的阻力均值減小。

    猜你喜歡
    整流罩方根值渦激
    支撐結(jié)構(gòu)對整流罩分離特性的影響分析①
    不同間距比下串聯(lián)圓柱渦激振動數(shù)值模擬研究
    磁流變彈性體減振單元動力學(xué)分析
    整流罩高空開傘完整回收——邁出整流罩落區(qū)控制與精確回收重要一步
    渦激振動發(fā)電裝置及其關(guān)鍵技術(shù)
    客車平順性仿真及優(yōu)化
    法爾肯-9火箭整流罩回收方案
    盤球立管結(jié)構(gòu)抑制渦激振動的數(shù)值分析方法研究
    電子制作(2018年14期)2018-08-21 01:38:42
    柔性圓管在渦激振動下的模態(tài)響應(yīng)分析
    Studyof velocityfluctuations in the plenum of a 3/4open jet automotive wind tunnel
    我要搜黄色片| 免费一级毛片在线播放高清视频| 偷拍熟女少妇极品色| 美女高潮喷水抽搐中文字幕| 三级男女做爰猛烈吃奶摸视频| 好看av亚洲va欧美ⅴa在| 丰满人妻熟妇乱又伦精品不卡| 在线观看日韩欧美| 男女之事视频高清在线观看| 国产精品香港三级国产av潘金莲| 亚洲 国产 在线| 日本黄大片高清| 国产黄片美女视频| 99久久99久久久精品蜜桃| 内地一区二区视频在线| 免费在线观看成人毛片| 精品欧美国产一区二区三| 欧美一区二区亚洲| 亚洲精品粉嫩美女一区| 久久精品国产自在天天线| 黄色成人免费大全| 国产高清三级在线| av片东京热男人的天堂| 99riav亚洲国产免费| 老鸭窝网址在线观看| 嫩草影院精品99| 国产精品一区二区三区四区久久| 男人舔女人下体高潮全视频| 无遮挡黄片免费观看| 一区二区三区激情视频| 在线国产一区二区在线| 国产精品 欧美亚洲| а√天堂www在线а√下载| 国产精品自产拍在线观看55亚洲| 国产午夜福利久久久久久| 人妻久久中文字幕网| 一本久久中文字幕| 色综合亚洲欧美另类图片| 婷婷精品国产亚洲av| 精品久久久久久久人妻蜜臀av| 久久久久久久亚洲中文字幕 | 嫩草影视91久久| av福利片在线观看| 日韩欧美在线乱码| 丰满人妻熟妇乱又伦精品不卡| 成人鲁丝片一二三区免费| 国产精品嫩草影院av在线观看 | 日本一二三区视频观看| 国产伦人伦偷精品视频| 国产一区二区三区在线臀色熟女| 草草在线视频免费看| 久久精品人妻少妇| 亚洲国产高清在线一区二区三| 亚洲av二区三区四区| 好看av亚洲va欧美ⅴa在| 亚洲天堂国产精品一区在线| 婷婷丁香在线五月| 99riav亚洲国产免费| 国产黄色小视频在线观看| 久久久国产精品麻豆| 日本 av在线| 97碰自拍视频| 麻豆成人av在线观看| 亚洲成av人片在线播放无| 日本在线视频免费播放| 国产精品美女特级片免费视频播放器| 免费大片18禁| 12—13女人毛片做爰片一| 一级黄色大片毛片| 熟女人妻精品中文字幕| 欧美性猛交╳xxx乱大交人| 天堂√8在线中文| 国产精品日韩av在线免费观看| 淫秽高清视频在线观看| 九九久久精品国产亚洲av麻豆| 亚洲国产精品sss在线观看| 国产精品综合久久久久久久免费| 亚洲色图av天堂| 亚洲午夜理论影院| 可以在线观看毛片的网站| 中文在线观看免费www的网站| 成人精品一区二区免费| 少妇高潮的动态图| 欧美精品啪啪一区二区三区| 99久久99久久久精品蜜桃| 日本在线视频免费播放| 国内精品久久久久久久电影| 禁无遮挡网站| 男插女下体视频免费在线播放| 制服丝袜大香蕉在线| 免费高清视频大片| 国产精品av视频在线免费观看| 成人三级黄色视频| 欧美又色又爽又黄视频| 我的老师免费观看完整版| 又黄又爽又免费观看的视频| 欧美日韩瑟瑟在线播放| 麻豆一二三区av精品| 黑人欧美特级aaaaaa片| 免费在线观看成人毛片| 欧美色视频一区免费| 日韩亚洲欧美综合| 1024手机看黄色片| 男插女下体视频免费在线播放| 亚洲五月婷婷丁香| 久久精品91无色码中文字幕| 搡老妇女老女人老熟妇| 女生性感内裤真人,穿戴方法视频| 一进一出好大好爽视频| 美女大奶头视频| 村上凉子中文字幕在线| 男人舔奶头视频| 成人特级黄色片久久久久久久| 成人永久免费在线观看视频| 精品无人区乱码1区二区| 成人午夜高清在线视频| 天天添夜夜摸| 女人十人毛片免费观看3o分钟| 亚洲七黄色美女视频| 成人午夜高清在线视频| 日日干狠狠操夜夜爽| 国产老妇女一区| 午夜日韩欧美国产| 午夜精品一区二区三区免费看| 99热这里只有是精品50| 高潮久久久久久久久久久不卡| 黄色丝袜av网址大全| 国产色爽女视频免费观看| 最新美女视频免费是黄的| 高清在线国产一区| 精品国内亚洲2022精品成人| 亚洲在线观看片| 天堂av国产一区二区熟女人妻| 最新美女视频免费是黄的| 日韩欧美在线二视频| 看免费av毛片| 怎么达到女性高潮| 中文字幕久久专区| 成人一区二区视频在线观看| 免费观看精品视频网站| 丁香六月欧美| 久久精品国产综合久久久| 麻豆一二三区av精品| 看片在线看免费视频| 脱女人内裤的视频| 午夜精品一区二区三区免费看| 淫妇啪啪啪对白视频| 一级毛片女人18水好多| 一区二区三区国产精品乱码| 国产精品电影一区二区三区| 在线国产一区二区在线| 老司机深夜福利视频在线观看| 国内精品一区二区在线观看| 久久香蕉精品热| 国内精品美女久久久久久| 国产三级中文精品| 丁香欧美五月| 久久精品国产清高在天天线| 亚洲片人在线观看| 亚洲片人在线观看| 亚洲,欧美精品.| 99在线人妻在线中文字幕| 国产欧美日韩精品亚洲av| 女生性感内裤真人,穿戴方法视频| 综合色av麻豆| 亚洲国产精品成人综合色| 中国美女看黄片| 日韩有码中文字幕| 夜夜夜夜夜久久久久| 日本一本二区三区精品| av专区在线播放| 国产免费男女视频| 国产极品精品免费视频能看的| 亚洲成人免费电影在线观看| 欧美日韩综合久久久久久 | 中文亚洲av片在线观看爽| 麻豆成人av在线观看| 免费在线观看影片大全网站| 中文字幕av成人在线电影| 成人三级黄色视频| 国产精品av视频在线免费观看| 精品国产美女av久久久久小说| 波多野结衣高清作品| 亚洲七黄色美女视频| 欧美xxxx黑人xx丫x性爽| 日韩欧美免费精品| 午夜福利在线观看吧| 一本一本综合久久| 无人区码免费观看不卡| 女人被狂操c到高潮| svipshipincom国产片| 欧美一级a爱片免费观看看| 欧美最黄视频在线播放免费| 国产精品野战在线观看| 欧美精品啪啪一区二区三区| 亚洲av五月六月丁香网| 国产免费一级a男人的天堂| 国产午夜精品久久久久久一区二区三区 | 一个人看的www免费观看视频| 国产一级毛片七仙女欲春2| 久久欧美精品欧美久久欧美| 中文字幕人成人乱码亚洲影| 久久久久久国产a免费观看| 免费大片18禁| 青草久久国产| 热99re8久久精品国产| 露出奶头的视频| 丁香欧美五月| 国产私拍福利视频在线观看| 亚洲国产欧美人成| 69av精品久久久久久| 国产亚洲av嫩草精品影院| 亚洲人成电影免费在线| 首页视频小说图片口味搜索| 丰满人妻熟妇乱又伦精品不卡| 18+在线观看网站| 变态另类成人亚洲欧美熟女| 日本黄色视频三级网站网址| 动漫黄色视频在线观看| 亚洲 国产 在线| 亚洲av美国av| 日日干狠狠操夜夜爽| 亚洲人与动物交配视频| 午夜福利成人在线免费观看| 免费人成在线观看视频色| 亚洲av二区三区四区| 又粗又爽又猛毛片免费看| 少妇丰满av| 成人特级av手机在线观看| 中亚洲国语对白在线视频| 久99久视频精品免费| 亚洲第一欧美日韩一区二区三区| 看片在线看免费视频| 好看av亚洲va欧美ⅴa在| 国产精品98久久久久久宅男小说| 国产高清视频在线播放一区| 亚洲精品日韩av片在线观看 | 舔av片在线| 性色avwww在线观看| 日本撒尿小便嘘嘘汇集6| 国内久久婷婷六月综合欲色啪| 亚洲无线观看免费| 悠悠久久av| 岛国视频午夜一区免费看| 亚洲av熟女| 色视频www国产| 真实男女啪啪啪动态图| 日韩国内少妇激情av| 国产高潮美女av| 国产不卡一卡二| 国产精品日韩av在线免费观看| 天天躁日日操中文字幕| 别揉我奶头~嗯~啊~动态视频| 欧美日韩亚洲国产一区二区在线观看| 成人精品一区二区免费| 夜夜看夜夜爽夜夜摸| 色在线成人网| 久久中文看片网| 俄罗斯特黄特色一大片| 国产精品久久视频播放| 精品一区二区三区视频在线 | 日本黄大片高清| 最好的美女福利视频网| 五月玫瑰六月丁香| 日韩欧美免费精品| 国产欧美日韩一区二区精品| 亚洲不卡免费看| 欧美绝顶高潮抽搐喷水| 波多野结衣高清无吗| 欧美午夜高清在线| 香蕉丝袜av| 国产美女午夜福利| 制服人妻中文乱码| 欧美日韩瑟瑟在线播放| 中文字幕av成人在线电影| 久久精品国产亚洲av香蕉五月| 我要搜黄色片| 久久久国产成人免费| 亚洲av不卡在线观看| 少妇的逼水好多| 老司机在亚洲福利影院| 国产野战对白在线观看| 女生性感内裤真人,穿戴方法视频| 嫁个100分男人电影在线观看| 欧美3d第一页| 亚洲欧美日韩高清在线视频| 久久精品亚洲精品国产色婷小说| 三级国产精品欧美在线观看| 亚洲国产精品999在线| 波多野结衣高清作品| 午夜福利视频1000在线观看| 首页视频小说图片口味搜索| 一级毛片女人18水好多| 99久久九九国产精品国产免费| 精品电影一区二区在线| 国产免费av片在线观看野外av| 亚洲av电影在线进入| 青草久久国产| 欧美成人免费av一区二区三区| 精品国内亚洲2022精品成人| 全区人妻精品视频| 国产午夜精品久久久久久一区二区三区 | 国产黄色小视频在线观看| 亚洲国产色片| 日韩欧美在线二视频| 国产亚洲精品一区二区www| 国产一级毛片七仙女欲春2| 男女做爰动态图高潮gif福利片| 亚洲无线在线观看| 禁无遮挡网站| 国产精品久久电影中文字幕| netflix在线观看网站| 国产极品精品免费视频能看的| 午夜a级毛片| 午夜亚洲福利在线播放| 欧美丝袜亚洲另类 | 18禁黄网站禁片免费观看直播| 午夜福利免费观看在线| 在线观看日韩欧美| 天堂网av新在线| 欧美日韩综合久久久久久 | 97超视频在线观看视频| 99国产精品一区二区蜜桃av| 成年免费大片在线观看| 日日夜夜操网爽| 日本五十路高清| 88av欧美| 性色avwww在线观看| 国产v大片淫在线免费观看| 国产av麻豆久久久久久久| 三级男女做爰猛烈吃奶摸视频| 亚洲av成人不卡在线观看播放网| АⅤ资源中文在线天堂| 老鸭窝网址在线观看| 蜜桃亚洲精品一区二区三区| 在线播放无遮挡| 我的老师免费观看完整版| 国产av不卡久久| 18禁黄网站禁片午夜丰满| 免费在线观看日本一区| 国语自产精品视频在线第100页| 国产精品av视频在线免费观看| 欧美性感艳星| 国产单亲对白刺激| 国产精品国产高清国产av| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品456在线播放app | 亚洲av免费在线观看| 日韩高清综合在线| 高清在线国产一区| 熟女电影av网| www.色视频.com| 国产成人影院久久av| 久久人妻av系列| 国产99白浆流出| 免费高清视频大片| 亚洲欧美日韩高清专用| 欧美一级a爱片免费观看看| 夜夜看夜夜爽夜夜摸| 日韩欧美一区二区三区在线观看| 老鸭窝网址在线观看| 亚洲人成网站在线播| 亚洲最大成人中文| 亚洲精品日韩av片在线观看 | 中国美女看黄片| 久久国产精品人妻蜜桃| 999久久久精品免费观看国产| 五月玫瑰六月丁香| 一进一出抽搐动态| 亚洲一区二区三区不卡视频| 18禁国产床啪视频网站| 最近在线观看免费完整版| 久久人妻av系列| 亚洲中文字幕一区二区三区有码在线看| 国内揄拍国产精品人妻在线| 老鸭窝网址在线观看| 欧洲精品卡2卡3卡4卡5卡区| 手机成人av网站| 日韩欧美一区二区三区在线观看| 国产成人影院久久av| 国产在线精品亚洲第一网站| 91字幕亚洲| 一夜夜www| 级片在线观看| 老司机在亚洲福利影院| 宅男免费午夜| 欧美色视频一区免费| av专区在线播放| 亚洲精品在线观看二区| 亚洲七黄色美女视频| 亚洲精品乱码久久久v下载方式 | 国产精品久久久人人做人人爽| 国产淫片久久久久久久久 | tocl精华| 国产探花在线观看一区二区| 波多野结衣高清无吗| 在线观看美女被高潮喷水网站 | 床上黄色一级片| 亚洲成a人片在线一区二区| 国产精品98久久久久久宅男小说| 国产精品免费一区二区三区在线| 女警被强在线播放| 又黄又爽又免费观看的视频| 欧美黄色片欧美黄色片| 非洲黑人性xxxx精品又粗又长| 丁香六月欧美| www日本黄色视频网| 久久久久九九精品影院| 免费在线观看日本一区| 嫩草影视91久久| 亚洲av一区综合| 精品久久久久久久久久久久久| 夜夜看夜夜爽夜夜摸| 变态另类丝袜制服| 97超级碰碰碰精品色视频在线观看| 久久久久免费精品人妻一区二区| 国产高清激情床上av| 久久午夜亚洲精品久久| 麻豆一二三区av精品| 精华霜和精华液先用哪个| 男人的好看免费观看在线视频| 9191精品国产免费久久| 久久久久久久久大av| 国产日本99.免费观看| 亚洲精品一卡2卡三卡4卡5卡| www.999成人在线观看| 免费在线观看成人毛片| 成年版毛片免费区| 一二三四社区在线视频社区8| 一级黄片播放器| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av美国av| 老司机午夜福利在线观看视频| 麻豆成人av在线观看| 19禁男女啪啪无遮挡网站| 一卡2卡三卡四卡精品乱码亚洲| 精品一区二区三区人妻视频| 欧美一区二区亚洲| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 青草久久国产| 国产亚洲精品av在线| 丰满人妻一区二区三区视频av | 小蜜桃在线观看免费完整版高清| 男女下面进入的视频免费午夜| 亚洲 国产 在线| 亚洲中文日韩欧美视频| 国产精品亚洲美女久久久| 全区人妻精品视频| 久久久久亚洲av毛片大全| 亚洲,欧美精品.| 亚洲国产精品成人综合色| 久久伊人香网站| 日韩欧美 国产精品| 欧美性感艳星| 最近最新中文字幕大全电影3| 精品人妻1区二区| 91久久精品国产一区二区成人 | 18+在线观看网站| 精品国产三级普通话版| 男人和女人高潮做爰伦理| 搡老熟女国产l中国老女人| 久久这里只有精品中国| 一本一本综合久久| 欧美日本视频| 丰满的人妻完整版| 51国产日韩欧美| 成年版毛片免费区| 亚洲国产欧美人成| 99久久综合精品五月天人人| 熟女电影av网| or卡值多少钱| 亚洲男人的天堂狠狠| 1000部很黄的大片| 久久6这里有精品| 天堂网av新在线| 午夜福利欧美成人| 欧美成狂野欧美在线观看| 欧美日韩精品网址| 国产一区二区在线av高清观看| 免费在线观看成人毛片| 网址你懂的国产日韩在线| 国产一区二区激情短视频| 亚洲欧美精品综合久久99| 日韩中文字幕欧美一区二区| 欧美一区二区亚洲| www.熟女人妻精品国产| 婷婷精品国产亚洲av在线| 嫩草影院入口| 亚洲av中文字字幕乱码综合| 热99在线观看视频| 可以在线观看的亚洲视频| 亚洲成av人片在线播放无| 久久99热这里只有精品18| 国产精品亚洲av一区麻豆| 午夜两性在线视频| 日日干狠狠操夜夜爽| 99精品在免费线老司机午夜| 色噜噜av男人的天堂激情| 日韩亚洲欧美综合| 特级一级黄色大片| 在线观看一区二区三区| 国产成人aa在线观看| 国产精品亚洲一级av第二区| 黄色视频,在线免费观看| 观看免费一级毛片| 婷婷精品国产亚洲av在线| 日本一二三区视频观看| 99国产精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 欧美成人免费av一区二区三区| 久久久久久人人人人人| 亚洲黑人精品在线| 午夜免费男女啪啪视频观看 | 麻豆成人av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美国产一区二区入口| 成人无遮挡网站| 国产伦在线观看视频一区| 国产精品一区二区三区四区久久| 成年免费大片在线观看| 成年女人看的毛片在线观看| 69av精品久久久久久| 亚洲国产欧美网| 亚洲,欧美精品.| 欧美日韩亚洲国产一区二区在线观看| 男人舔奶头视频| 精品欧美国产一区二区三| 欧美日韩精品网址| 综合色av麻豆| xxx96com| 在线观看一区二区三区| 国产日本99.免费观看| 成年女人看的毛片在线观看| av视频在线观看入口| 小蜜桃在线观看免费完整版高清| 国产免费av片在线观看野外av| 在线免费观看不下载黄p国产 | 日本免费a在线| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 久久久国产成人精品二区| 成人av一区二区三区在线看| 日韩欧美精品免费久久 | 国内精品久久久久精免费| 色噜噜av男人的天堂激情| 免费在线观看亚洲国产| 日本黄大片高清| 亚洲黑人精品在线| 亚洲国产欧美网| 黄片小视频在线播放| 午夜精品久久久久久毛片777| 久久精品国产综合久久久| 婷婷六月久久综合丁香| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区三区四区免费观看 | 免费av不卡在线播放| 久久精品影院6| 久久久久性生活片| 无限看片的www在线观看| 午夜福利欧美成人| 一区二区三区免费毛片| 最新美女视频免费是黄的| 成人av一区二区三区在线看| 中文字幕熟女人妻在线| 色av中文字幕| 中文亚洲av片在线观看爽| 日韩高清综合在线| 老司机深夜福利视频在线观看| 亚洲av成人av| 性欧美人与动物交配| 国产视频一区二区在线看| 亚洲五月婷婷丁香| 精品一区二区三区视频在线观看免费| 国产精品久久视频播放| 免费无遮挡裸体视频| 亚洲黑人精品在线| 成年版毛片免费区| 69av精品久久久久久| 天天一区二区日本电影三级| 国内毛片毛片毛片毛片毛片| 99久久99久久久精品蜜桃| 久久精品综合一区二区三区| 天堂网av新在线| 亚洲精品色激情综合| 丰满人妻熟妇乱又伦精品不卡| 国产不卡一卡二| 亚洲人成网站在线播| 国产不卡一卡二| 成熟少妇高潮喷水视频| 国产不卡一卡二| 国产成人啪精品午夜网站| 国产伦精品一区二区三区视频9 | 真实男女啪啪啪动态图| 精品国产美女av久久久久小说| 亚洲国产中文字幕在线视频| av欧美777| 日本在线视频免费播放| 免费电影在线观看免费观看| 国产综合懂色| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 日韩欧美三级三区| 男女那种视频在线观看| 国产91精品成人一区二区三区| 99久久成人亚洲精品观看| 欧美午夜高清在线| 最新美女视频免费是黄的| 亚洲av成人av| 亚洲av中文字字幕乱码综合| 床上黄色一级片| 亚洲中文日韩欧美视频| 亚洲av电影在线进入| 成人精品一区二区免费|