• <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国产av综合av卡| 日韩在线高清观看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩人妻高清精品专区| 成人免费观看视频高清| 肉色欧美久久久久久久蜜桃| 亚洲精品aⅴ在线观看| 日韩一本色道免费dvd| 久久久久久久久大av| 国产男人的电影天堂91| 亚洲内射少妇av| 久久影院123| 一二三四中文在线观看免费高清| 日韩欧美精品免费久久| 2021少妇久久久久久久久久久| 国产 精品1| 久久精品国产鲁丝片午夜精品| 欧美日本中文国产一区发布| 精品久久久久久久久亚洲| 青春草亚洲视频在线观看| 国产美女午夜福利| 青春草视频在线免费观看| 久久久久久久久久人人人人人人| 精品一区二区免费观看| 插阴视频在线观看视频| 99热这里只有精品一区| 成人免费观看视频高清| 久久精品久久久久久噜噜老黄| 国产视频首页在线观看| 99久久综合免费| 久久99热6这里只有精品| 国产av一区二区精品久久| 高清午夜精品一区二区三区| 中国三级夫妇交换| 99久久精品国产国产毛片| 十八禁高潮呻吟视频 | 一区二区三区免费毛片| 免费黄频网站在线观看国产| 欧美区成人在线视频| 亚洲美女视频黄频| 好男人视频免费观看在线| 纯流量卡能插随身wifi吗| 伊人久久精品亚洲午夜| 国产一区二区在线观看日韩| 亚洲婷婷狠狠爱综合网| 国产色婷婷99| 免费高清在线观看视频在线观看| 国产亚洲最大av| 狂野欧美白嫩少妇大欣赏| 久久这里有精品视频免费| 成人毛片60女人毛片免费| 国产欧美另类精品又又久久亚洲欧美| videossex国产| 精品亚洲乱码少妇综合久久| 免费观看av网站的网址| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品无大码| 欧美性感艳星| 免费看日本二区| 欧美人与善性xxx| 91aial.com中文字幕在线观看| 国产日韩一区二区三区精品不卡 | 欧美成人午夜免费资源| 国产精品久久久久久久电影| 亚洲欧洲日产国产| 97精品久久久久久久久久精品| 一级毛片aaaaaa免费看小| 国产永久视频网站| 七月丁香在线播放| 建设人人有责人人尽责人人享有的| 黑丝袜美女国产一区| 精品酒店卫生间| 欧美精品一区二区大全| 街头女战士在线观看网站| 97超视频在线观看视频| 亚洲国产精品一区三区| 夜夜骑夜夜射夜夜干| 免费看日本二区| 国产av一区二区精品久久| 人妻夜夜爽99麻豆av| 各种免费的搞黄视频| 青青草视频在线视频观看| 一边摸一边抽搐一进一出视频| 日本五十路高清| 国产av一区二区精品久久| 最近中文字幕2019免费版| 欧美亚洲 丝袜 人妻 在线| 久久女婷五月综合色啪小说| 亚洲精品自拍成人| 国产97色在线日韩免费| 国产一区二区 视频在线| 一本大道久久a久久精品| 天堂俺去俺来也www色官网| 亚洲少妇的诱惑av| 亚洲,欧美精品.| 免费观看av网站的网址| 别揉我奶头~嗯~啊~动态视频 | 中文字幕制服av| 国产成人系列免费观看| 大香蕉久久成人网| 午夜福利视频精品| 啦啦啦免费观看视频1| av网站免费在线观看视频| 亚洲欧洲精品一区二区精品久久久| 黄网站色视频无遮挡免费观看| 十八禁网站免费在线| 免费在线观看黄色视频的| 精品高清国产在线一区| 精品一品国产午夜福利视频| 亚洲天堂av无毛| 精品高清国产在线一区| 精品视频人人做人人爽| 国产精品成人在线| 午夜91福利影院| av在线老鸭窝| 免费人妻精品一区二区三区视频| 9色porny在线观看| 欧美日韩av久久| 久久久久久久国产电影| 欧美日韩中文字幕国产精品一区二区三区 | 欧美精品高潮呻吟av久久| 肉色欧美久久久久久久蜜桃| 三级毛片av免费| 男女国产视频网站| 国产成人精品久久二区二区免费| 下体分泌物呈黄色| 91成人精品电影| 久热这里只有精品99| 国产精品.久久久| 老鸭窝网址在线观看| 在线观看免费午夜福利视频| 亚洲国产中文字幕在线视频| 免费一级毛片在线播放高清视频 | 精品久久蜜臀av无| 久久影院123| 亚洲自偷自拍图片 自拍| 一区二区三区激情视频| 亚洲成人免费av在线播放| 免费在线观看完整版高清| 性色av乱码一区二区三区2| 欧美人与性动交α欧美精品济南到| 老司机福利观看| 一级毛片精品| av免费在线观看网站| 亚洲伊人色综图| 久久亚洲国产成人精品v| 性色av乱码一区二区三区2| 国产av又大| 久久久国产一区二区| 日日夜夜操网爽| 国产精品免费视频内射| 精品人妻一区二区三区麻豆| 捣出白浆h1v1| 交换朋友夫妻互换小说| 99久久精品国产亚洲精品| 久久久久久久久免费视频了| 老司机影院成人| 国产亚洲精品第一综合不卡| 久久亚洲国产成人精品v| 精品免费久久久久久久清纯 | 午夜影院在线不卡| 亚洲精品粉嫩美女一区| 久久免费观看电影| 桃花免费在线播放| 淫妇啪啪啪对白视频 | av欧美777| 91麻豆精品激情在线观看国产 | 一区二区三区四区激情视频| 国产成人影院久久av| 国产区一区二久久| 国产一区二区 视频在线| 不卡av一区二区三区| 男女之事视频高清在线观看| 99久久国产精品久久久| 精品少妇内射三级| 99久久人妻综合| 亚洲九九香蕉| 黄片大片在线免费观看| 日本猛色少妇xxxxx猛交久久| 热99国产精品久久久久久7| av在线老鸭窝| 国产有黄有色有爽视频| 极品人妻少妇av视频| 在线观看免费高清a一片| 亚洲精品第二区| 国产av一区二区精品久久| av网站免费在线观看视频| 法律面前人人平等表现在哪些方面 | av在线播放精品| 一区二区三区乱码不卡18| 精品一区二区三区av网在线观看 | 亚洲av电影在线观看一区二区三区| 又紧又爽又黄一区二区| 国产不卡av网站在线观看| 亚洲人成电影观看| 中国国产av一级| 一级黄色大片毛片| 99精品欧美一区二区三区四区| 亚洲黑人精品在线| 十分钟在线观看高清视频www| 国产精品欧美亚洲77777| 国产老妇伦熟女老妇高清| 一边摸一边做爽爽视频免费| 国产成人欧美| 国产免费av片在线观看野外av| 性高湖久久久久久久久免费观看| 天天添夜夜摸| 亚洲一码二码三码区别大吗| 超色免费av| 高清av免费在线| 国产免费视频播放在线视频| 日韩免费高清中文字幕av| 老司机在亚洲福利影院| 国产无遮挡羞羞视频在线观看| 999久久久国产精品视频| 一区二区日韩欧美中文字幕| 视频区欧美日本亚洲| 国产亚洲精品久久久久5区| 丝袜人妻中文字幕| 99精品久久久久人妻精品| 熟女少妇亚洲综合色aaa.| 99久久综合免费| 美女国产高潮福利片在线看| 性色av一级| 美女中出高潮动态图| 女性生殖器流出的白浆| 欧美日韩亚洲综合一区二区三区_| 国产福利在线免费观看视频| 午夜精品久久久久久毛片777| 18禁裸乳无遮挡动漫免费视频| 2018国产大陆天天弄谢| 精品第一国产精品| 国产成+人综合+亚洲专区| 女人高潮潮喷娇喘18禁视频| 亚洲av美国av| 性高湖久久久久久久久免费观看| 国产色视频综合| 亚洲伊人久久精品综合| 精品一区在线观看国产| 国精品久久久久久国模美| 久久人人爽人人片av| 999久久久国产精品视频| 亚洲第一av免费看| 日本黄色日本黄色录像| 中文字幕另类日韩欧美亚洲嫩草| 亚洲一区二区三区欧美精品| 久久香蕉激情| 精品久久蜜臀av无| 9色porny在线观看| 日韩免费高清中文字幕av| 黑人巨大精品欧美一区二区mp4| 午夜免费成人在线视频| 日韩电影二区| 男女边摸边吃奶| 操美女的视频在线观看| 中文字幕最新亚洲高清| 麻豆乱淫一区二区| 亚洲精品一区蜜桃| 午夜免费观看性视频| 国产一区二区三区在线臀色熟女 | 国产欧美日韩一区二区精品| 国产日韩一区二区三区精品不卡| 99久久99久久久精品蜜桃| 韩国高清视频一区二区三区| 一个人免费看片子| 欧美激情极品国产一区二区三区| 黄频高清免费视频| 成人黄色视频免费在线看| 精品第一国产精品| 99国产精品一区二区三区| 老司机靠b影院| 青春草视频在线免费观看| 大香蕉久久网| 免费观看a级毛片全部| 国产精品免费大片| 精品高清国产在线一区| 午夜两性在线视频| 午夜福利视频在线观看免费| 十八禁网站网址无遮挡| 麻豆av在线久日| 欧美精品av麻豆av| 亚洲 国产 在线| 久久中文字幕一级| 日本wwww免费看| 午夜福利视频精品| 美女高潮喷水抽搐中文字幕| 久久久国产精品麻豆| 亚洲精华国产精华精| 最新在线观看一区二区三区| 在线观看舔阴道视频| 91成人精品电影| 亚洲久久久国产精品| 欧美在线一区亚洲| 国产精品亚洲av一区麻豆| 亚洲专区字幕在线| 精品人妻1区二区| 亚洲久久久国产精品| 亚洲一区二区三区欧美精品| 日本猛色少妇xxxxx猛交久久| 国产成人精品久久二区二区免费| 久久久久久免费高清国产稀缺| 欧美午夜高清在线| 成年人黄色毛片网站| 国产成人啪精品午夜网站| av电影中文网址| 99国产极品粉嫩在线观看| 午夜福利,免费看| 欧美在线一区亚洲| 国产日韩一区二区三区精品不卡| 亚洲va日本ⅴa欧美va伊人久久 | 国产欧美日韩精品亚洲av| 美女大奶头黄色视频| 亚洲国产av影院在线观看| 黑人欧美特级aaaaaa片| 亚洲国产精品成人久久小说| 久久人人97超碰香蕉20202| 91麻豆精品激情在线观看国产 | 麻豆国产av国片精品| 亚洲欧洲精品一区二区精品久久久| 欧美激情久久久久久爽电影 | 看免费av毛片| 丝瓜视频免费看黄片| 欧美黑人欧美精品刺激| 精品久久久久久电影网| 亚洲 欧美一区二区三区| 又紧又爽又黄一区二区| 欧美97在线视频| 丝袜人妻中文字幕| 亚洲avbb在线观看| 欧美日韩av久久| 久久精品国产亚洲av香蕉五月 | 9热在线视频观看99| 欧美黑人精品巨大| 99国产精品一区二区蜜桃av | 国产精品久久久av美女十八| 色精品久久人妻99蜜桃| 国产精品 欧美亚洲| 99国产精品一区二区三区| 热99久久久久精品小说推荐| 91国产中文字幕| 人妻人人澡人人爽人人| 午夜免费成人在线视频| 精品人妻1区二区| 久久精品人人爽人人爽视色| 黄色片一级片一级黄色片| 我要看黄色一级片免费的| 久久精品国产亚洲av香蕉五月 | 俄罗斯特黄特色一大片| 两性夫妻黄色片| 久久久久久久久免费视频了| 亚洲欧美日韩另类电影网站| 免费在线观看完整版高清| 这个男人来自地球电影免费观看| 黄色a级毛片大全视频| 亚洲熟女毛片儿| 在线永久观看黄色视频| 日韩大码丰满熟妇| 午夜成年电影在线免费观看| 精品国产超薄肉色丝袜足j| 五月开心婷婷网| 一进一出抽搐动态| 久久国产亚洲av麻豆专区| www.自偷自拍.com| 日韩熟女老妇一区二区性免费视频| 黄色片一级片一级黄色片| 岛国在线观看网站| 两个人免费观看高清视频| 国产精品久久久久久人妻精品电影 | 精品免费久久久久久久清纯 | 在线观看免费视频网站a站| 精品免费久久久久久久清纯 | 考比视频在线观看| 中文欧美无线码| 青青草视频在线视频观看| 超碰97精品在线观看| 亚洲男人天堂网一区| 嫁个100分男人电影在线观看| 欧美大码av| 性色av一级| 捣出白浆h1v1| 操出白浆在线播放| 国产成人免费无遮挡视频| 人人妻,人人澡人人爽秒播| 69av精品久久久久久 | 欧美日韩亚洲国产一区二区在线观看 | 在线十欧美十亚洲十日本专区| www.999成人在线观看| 热99re8久久精品国产| 国产视频一区二区在线看| 日韩欧美免费精品| 国产精品一区二区在线不卡| 亚洲国产毛片av蜜桃av| 香蕉国产在线看| 超碰成人久久| 国产伦理片在线播放av一区| 搡老岳熟女国产| 狠狠狠狠99中文字幕| 国产成人av激情在线播放| 丝瓜视频免费看黄片| 免费女性裸体啪啪无遮挡网站| 老司机午夜福利在线观看视频 | 国产亚洲午夜精品一区二区久久| 日韩一区二区三区影片| 国产一区二区三区综合在线观看| 夜夜夜夜夜久久久久| 久久精品亚洲熟妇少妇任你| 日韩视频在线欧美| 欧美日韩成人在线一区二区| 国产成人a∨麻豆精品| 精品亚洲成a人片在线观看| 母亲3免费完整高清在线观看| 色视频在线一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲九九香蕉| 我的亚洲天堂| 99久久人妻综合| 男女下面插进去视频免费观看| 精品国内亚洲2022精品成人 | 精品国产乱码久久久久久男人| 国产人伦9x9x在线观看| 亚洲专区中文字幕在线| 美女脱内裤让男人舔精品视频| 看免费av毛片| 男女无遮挡免费网站观看| 两性夫妻黄色片| 精品国产国语对白av| 国产伦理片在线播放av一区| 操出白浆在线播放| 岛国在线观看网站| 大香蕉久久网| 在线av久久热| 国产精品国产三级国产专区5o| 国产在线免费精品| 老司机福利观看| 极品少妇高潮喷水抽搐| 天天操日日干夜夜撸| 久久国产精品人妻蜜桃| 波多野结衣av一区二区av| 大香蕉久久成人网| 丝袜美足系列| 久热这里只有精品99| 最黄视频免费看| 不卡av一区二区三区| 成人免费观看视频高清| 大码成人一级视频| 亚洲国产欧美日韩在线播放| 黄频高清免费视频| 亚洲精品久久午夜乱码| 国产日韩欧美亚洲二区| 欧美在线一区亚洲| 一级毛片精品| 青春草视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 下体分泌物呈黄色| 精品国产乱码久久久久久男人| 50天的宝宝边吃奶边哭怎么回事| 午夜福利视频在线观看免费| 欧美亚洲 丝袜 人妻 在线| 搡老熟女国产l中国老女人| 国产精品久久久人人做人人爽| 精品国产一区二区久久| 美女国产高潮福利片在线看| 亚洲av日韩在线播放| 久久热在线av| 国产精品麻豆人妻色哟哟久久| 操出白浆在线播放| 亚洲国产精品成人久久小说| www.999成人在线观看| 欧美一级毛片孕妇| 婷婷色av中文字幕| 国产伦人伦偷精品视频| 国产99久久九九免费精品| 国产日韩欧美在线精品| 夜夜夜夜夜久久久久| 国产精品久久久久久精品古装| 久久人妻熟女aⅴ| 人人妻人人添人人爽欧美一区卜| 天天影视国产精品| 精品国产乱码久久久久久小说| 亚洲熟女毛片儿| 最近最新免费中文字幕在线| 美女福利国产在线| 人妻一区二区av| 日本猛色少妇xxxxx猛交久久| 纯流量卡能插随身wifi吗| 中亚洲国语对白在线视频| 欧美激情久久久久久爽电影 | 一二三四在线观看免费中文在| 国产一区二区激情短视频 | 精品一品国产午夜福利视频| 久久这里只有精品19| 天天躁夜夜躁狠狠躁躁| 久久久久国内视频| 老熟妇乱子伦视频在线观看 | 99香蕉大伊视频| 97精品久久久久久久久久精品| cao死你这个sao货| 91精品伊人久久大香线蕉| 亚洲久久久国产精品| 亚洲欧洲日产国产| 日本91视频免费播放| 日本av免费视频播放| 丰满饥渴人妻一区二区三| 久久精品国产亚洲av高清一级| 99九九在线精品视频| 人妻一区二区av| 午夜成年电影在线免费观看| videos熟女内射| 亚洲一卡2卡3卡4卡5卡精品中文| 下体分泌物呈黄色| 在线看a的网站| 久久精品国产综合久久久| 中文字幕人妻丝袜一区二区| 亚洲精品国产色婷婷电影| 日本91视频免费播放| 操美女的视频在线观看| 中文字幕制服av| 在线观看免费日韩欧美大片| 欧美精品啪啪一区二区三区 | 国产麻豆69| 成人国产一区最新在线观看| 国产欧美日韩综合在线一区二区| 成年av动漫网址| 热99re8久久精品国产| 嫩草影视91久久| 久久久久久亚洲精品国产蜜桃av| av在线老鸭窝| 久久这里只有精品19| 成年av动漫网址| 欧美黑人精品巨大| 色综合欧美亚洲国产小说| 午夜免费成人在线视频| 国产精品秋霞免费鲁丝片| 亚洲国产欧美日韩在线播放| 桃红色精品国产亚洲av| 久久国产精品影院| 男女下面插进去视频免费观看| 欧美日韩国产mv在线观看视频| 欧美日本中文国产一区发布| 成人18禁高潮啪啪吃奶动态图| 男女无遮挡免费网站观看| 久久久久久亚洲精品国产蜜桃av| 久久久久久久国产电影| 成年人黄色毛片网站| 黑人巨大精品欧美一区二区mp4| 在线看a的网站| 夜夜骑夜夜射夜夜干| 午夜福利免费观看在线| 人人妻人人添人人爽欧美一区卜| 日韩欧美一区二区三区在线观看 | 亚洲精品一卡2卡三卡4卡5卡 | 一级片'在线观看视频| 国产在线一区二区三区精| 色婷婷av一区二区三区视频| 亚洲全国av大片| 麻豆国产av国片精品| 十八禁高潮呻吟视频| 国产精品秋霞免费鲁丝片| 午夜精品国产一区二区电影| tocl精华| 亚洲精品国产区一区二| 18禁黄网站禁片午夜丰满| 每晚都被弄得嗷嗷叫到高潮| 青草久久国产| 亚洲成人免费电影在线观看| 久久精品成人免费网站| 又紧又爽又黄一区二区| 亚洲国产欧美一区二区综合| 成年人午夜在线观看视频| 精品福利永久在线观看| 在线十欧美十亚洲十日本专区| 免费观看av网站的网址| 国产老妇伦熟女老妇高清| 日本a在线网址| 黄色视频在线播放观看不卡| 精品免费久久久久久久清纯 | 伊人亚洲综合成人网| 丝袜脚勾引网站| 男女边摸边吃奶| 亚洲专区国产一区二区| 两性夫妻黄色片| 999久久久国产精品视频| 日韩有码中文字幕| 97在线人人人人妻| 夜夜夜夜夜久久久久| 亚洲欧美色中文字幕在线| 欧美一级毛片孕妇| 80岁老熟妇乱子伦牲交| 欧美少妇被猛烈插入视频| 在线永久观看黄色视频| 正在播放国产对白刺激| 一本一本久久a久久精品综合妖精| 亚洲欧洲精品一区二区精品久久久| 精品乱码久久久久久99久播| 天堂8中文在线网| 十八禁网站网址无遮挡| 男男h啪啪无遮挡| 黑丝袜美女国产一区| 美国免费a级毛片| 欧美在线一区亚洲| 九色亚洲精品在线播放| 国产一区二区三区av在线| 婷婷色av中文字幕| 在线观看免费午夜福利视频| 精品人妻1区二区| 中文字幕色久视频| av欧美777| 精品少妇久久久久久888优播| 国产伦理片在线播放av一区| 午夜免费鲁丝| 国产一区二区三区在线臀色熟女 | 成年人午夜在线观看视频| 永久免费av网站大全| 国产视频一区二区在线看|