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

    自由面對潛艇尾流場流動(dòng)特性影響研究1)

    2021-03-24 06:13:32郭春雨徐鵬韓陽王超郭欣雨鄶云飛
    力學(xué)學(xué)報(bào) 2021年1期
    關(guān)鍵詞:盤面尾流航速

    郭春雨 徐鵬 韓陽,2) 王超 郭欣雨 鄶云飛

    ?(哈爾濱工程大學(xué)船舶工程學(xué)院,哈爾濱 150001)

    ?(水下測控技術(shù)重點(diǎn)實(shí)驗(yàn)室,大連 116013)

    引言

    潛艇的隱身性能不僅是保障完成其作戰(zhàn)使命的基本性能,更有助于潛艇在海戰(zhàn)中更好地發(fā)揮其強(qiáng)大的突襲能力.隨著如今作戰(zhàn)任務(wù)多樣性的提升,潛艇不可避免地會(huì)涉及到近自由面航行,艇體的分離渦與自由面相互作用,致使渦結(jié)構(gòu)會(huì)在原有結(jié)構(gòu)基礎(chǔ)上產(chǎn)生位置偏移,進(jìn)一步加劇流場的復(fù)雜程度,這對潛艇的隱身性能來說是一項(xiàng)更大的挑戰(zhàn).由于戰(zhàn)略戰(zhàn)術(shù)的需求,有關(guān)潛艇近水面航行的研究逐漸受到重視.近年來隨著水動(dòng)力學(xué)研究的深入,越來越多的科研人員通過流場流動(dòng)細(xì)節(jié)來尋找影響宏觀力及噪聲的本質(zhì)原因.對于精細(xì)流場的獲取主要有兩種途徑:試驗(yàn)測試與數(shù)值模擬.其中數(shù)值模擬由于成本較低且能獲取全場流動(dòng)信息而受到廣泛應(yīng)用,但其可信程度稍低,往往需要試驗(yàn)數(shù)據(jù)加以驗(yàn)證.隨著流場測試技術(shù)的不斷發(fā)展,粒子圖像測速(particle image velocimetry,PIV)因其無接觸、全局、瞬態(tài)測量等優(yōu)點(diǎn),迅速成為研究者的利器[1-4].

    國內(nèi)外公開發(fā)表的運(yùn)用PIV 技術(shù)測量潛艇繞流場的文章很少.Ashok 等[5-6]對SUBOFF 標(biāo)模進(jìn)行了系列試驗(yàn)并進(jìn)行了分析.Fureby 等[7]對直航和10?斜航角狀態(tài)下的潛艇流場進(jìn)行測量,并與數(shù)值模擬結(jié)果進(jìn)行對比,兩者吻合良好.Lee 等[8]通過SPIV 技術(shù)測量了潛艇10?斜航角狀態(tài)下指揮臺圍殼分離渦情況.Wang 等[9]研究了自由面的存在對螺旋槳尾流的影響,并與CFD 計(jì)算結(jié)果相互驗(yàn)證.在國內(nèi),張軍等[10-11]利用隨車式PIV 系統(tǒng)測量了不同主體與附體連接形式的水下航行體尾部流場.張國平等[12]應(yīng)用隨車式2D-3C 水下PIV 測量系統(tǒng)進(jìn)行了SUBOFF潛艇流場測試,然而并未公開相關(guān)測試結(jié)果.薛慶雨等[13-15]在空泡水筒中測量了SUBOFF 模型直航、10?攻角、尾附體10?舵角等三種工況下槳盤面流場湍流強(qiáng)度和湍流積分長度.海軍工程大學(xué)張文照等[16]對潛艇模型的尾流速度場進(jìn)行測量,得到潛艇模型的尾流渦量場分布,為研究潛艇水動(dòng)力噪聲與尾流渦結(jié)構(gòu)之間的關(guān)系提供了參考.

    由于流場測量困難,鮮有考慮自由面對潛艇尾流場流動(dòng)特性影響的試驗(yàn)研究,大都采用數(shù)值模擬的方法探究艇體與自由面相互作用的影響.Dawson[17]是詳細(xì)開展?jié)撏Ы杂擅婧叫袛?shù)值模擬較早的研究者之一.Shariati 等[18-19]應(yīng)用RANS 方法對SUBOFF 模型多種工況、航速、潛深做了詳細(xì)研究.近期,Amiri等[20-22]對自由面如何影響潛艇水動(dòng)力性能進(jìn)行了系列研究.2019 年Carrica 等[23]考慮靜水環(huán)境及波浪環(huán)境對Joubert BB2 潛艇近自由面自航的水動(dòng)力性能進(jìn)行了研究.Tian 等[24]研究了不同波浪作用以及不同潛深對水下航行器水動(dòng)力的影響.在國內(nèi),張楠等[25-28]開展了大量關(guān)于潛艇近水面航行數(shù)值模擬研究,奠定了國內(nèi)自由面對潛艇尾流影響數(shù)值研究的基礎(chǔ).隨后,逐漸受到各科研高校的關(guān)注,涌現(xiàn)出大批有關(guān)自由面對潛艇尾流影響的論文[29-34].

    采用精細(xì)流場測量技術(shù)對潛艇繞流場的試驗(yàn)研究很少,有關(guān)自由面對潛艇尾流場流動(dòng)特性影響的研究更是少之又少.本文依托哈爾濱工程大學(xué)船模水池實(shí)驗(yàn)室,采用先進(jìn)的水下三維粒子圖像測速(stereoscopic particle image velocimetry,SPIV) 系統(tǒng)測量潛艇尾流場微觀結(jié)構(gòu),輔以數(shù)值模擬技術(shù)觀測近自由面航行時(shí)表面波系,以此來研究SUBOFF 潛艇無限潛深與近水面航行時(shí)尾流場特性,通過試驗(yàn)與數(shù)值模擬結(jié)果相互驗(yàn)證、補(bǔ)充,得到潛艇尾流場流動(dòng)特征以及隨航速與潛深變化的規(guī)律,以期為未來潛艇設(shè)計(jì)提供依據(jù).

    1 試驗(yàn)方案

    1.1 試驗(yàn)設(shè)備與試驗(yàn)?zāi)P瓦x取

    1.1.1 試驗(yàn)設(shè)備

    本文依托哈爾濱工程大學(xué)船模拖曳水池開展SUBOFF 模型尾流場SPIV 試驗(yàn),基本實(shí)驗(yàn)設(shè)備如下.

    (1)拖曳水池:長度108 m;寬度7 m;水深3.5 m.

    (2) 拖車:車速V≤6.5 m/s;穩(wěn)速范圍:0.1~6.5 m/s;精度:0.1%;通過自主設(shè)計(jì)并制作的固定裝置將試驗(yàn)?zāi)P凸潭ㄔ谕宪嚿?固定方式采用艇尾固定.

    (3)水下粒子圖像測速系統(tǒng):2D-3C SPIV 系統(tǒng).

    1.1.2 模型選取

    本試驗(yàn)選用的模型是美國泰勒水池(David Taylor Model Basin,DTMB) 設(shè)計(jì)的SUBOFF 模型,該模型具有大量的試驗(yàn)結(jié)果進(jìn)行驗(yàn)證,因此比較適合做相關(guān)研究.由于水池深度及PIV 測試能力的限制,模型長度選為1 m,模型如圖1 所示,其主要參數(shù)如表1所示.

    圖1 SUBOFF 模型圖Fig.1 SUBOFF model

    試驗(yàn)所用模型由聚乳酸(PLA) 材料經(jīng)3D 打印技術(shù)制作而成.聚乳酸密度為1.20~1.30 kg/L,彈性模量為3000~4000 MPa,拉伸強(qiáng)度為40~60 MPa.由于打印機(jī)工作區(qū)域限制,整個(gè)艇體分三段打印,最后拼接而成.首先將艇體表面打磨光順并噴涂防水涂層,然后對艇體表面進(jìn)行亞光黑處理,最后對測量截面附近進(jìn)行噴涂吸光材料,避免試驗(yàn)過程中激光強(qiáng)反光影響.由于本次試驗(yàn)采用的潛艇模型縮尺比較小,故在艇體首部L/20 處加裝激流裝置以加快艇體表面的層流轉(zhuǎn)淚.經(jīng)3D 打印的實(shí)物模型以及噴涂效果如圖2 所示.

    圖2 SUBOFF 實(shí)物模型Fig.2 SUBOFF solid model

    1.2 試驗(yàn)工況的確定

    為分析自由面對潛艇尾流場的影響,進(jìn)行三個(gè)潛深、多個(gè)航速下的潛艇尾流場測量.由文獻(xiàn)[17]可知,潛深在3.3 倍最大直徑D以上時(shí)可以忽略自由面的影響;同時(shí)由另一文獻(xiàn)[35]可知,除了自由面影響外,池底壁面也存在影響,要完全避免限制水域的干擾需要水下航行體距池底的距離不宜小于1/2 倍艇長.因此本次試驗(yàn)選取4D潛深作為無限潛深工況,此時(shí)可以忽略自由面和池底的影響.同時(shí),在潛艇靠近自由面時(shí)存在自由液面反光干擾,為減小反光現(xiàn)象,選取1.5D潛深作為近自由面工況,選取2.5D潛深作為中間工況.試驗(yàn)航速選取為Fr=0.3~0.7,傅汝德數(shù)Fr定義如下

    式中,U為拖曳航速,L為潛艇的特征長度.

    測量截面取為槳盤面處,模型參數(shù)及試驗(yàn)工況如表1 所示.

    表1 模型參數(shù)及試驗(yàn)工況Table 1 Model parameters and test conditions

    1.3 數(shù)值模型與網(wǎng)格劃分

    采用基于有限體積法的STAR-CCM+商業(yè)軟件進(jìn)行潛艇流場數(shù)值模擬,應(yīng)用隱式非定常時(shí)間推進(jìn)格式進(jìn)行非穩(wěn)態(tài)求解,湍流模型選取分離渦模型(DES),時(shí)間步長的選取則充分考慮CFL~1 的條件,綜合網(wǎng)格空間尺度及時(shí)間尺度,最終選為0.001 s.由于涉及到近自由面航行,自由表面采用兩相體積VOF 方法建模.控制方程采用基于壓力的耦合求解,其中對流項(xiàng)采用二階迎風(fēng)格式進(jìn)行空間離散,耗散項(xiàng)采用二階中心差分格式進(jìn)行離散.計(jì)算時(shí)長為定為2 s.

    計(jì)算域如圖3 所示,具體設(shè)定如下:速度入口邊界距艇首2 倍艇長L,壓力出口距艇尾4 倍艇長L,其余邊界設(shè)為壁面邊界條件,距艇中心10 倍艇徑D.

    圖3 計(jì)算域示意圖Fig.3 The calculation domain schematic

    網(wǎng)格劃分采用軟件中自帶的切割體網(wǎng)格形式,在艇體周圍進(jìn)行三層局部加密,在主附體交接部分及分離渦發(fā)展區(qū)域(指揮臺圍殼后端及十字舵后端)進(jìn)行局部加密,對潛艇近水面航行時(shí)的自由表面和開爾文波形等物理現(xiàn)象也進(jìn)行加密.邊界層網(wǎng)格劃分則采用棱柱層網(wǎng)格形式,并采用Ally+壁面處理,同時(shí)根據(jù)不同航速進(jìn)行棱柱層厚度調(diào)整,保證y+值處于1 左右范圍內(nèi).根據(jù)不同航速下的興波程度有針對性地加密自由表面及開爾文波形,最終將網(wǎng)格數(shù)控制在1000 萬左右,艇體周圍各加密區(qū)細(xì)節(jié)如圖4 所示.

    2 潛艇尾流場測量試驗(yàn)

    2.1 水下SPIV 測量系統(tǒng)

    圖4 艇體周圍加密區(qū)示意圖Fig.4 The mesh refinement area

    本試驗(yàn)所采用的流場測量設(shè)備為DANTEC 公司為哈爾濱工程大學(xué)船模拖曳水池定制的隨車式水下SPIV 系統(tǒng),布置于拖車右側(cè),SPIV 系統(tǒng)規(guī)格參數(shù)如下:CCD 相機(jī)分辨率為2048×2048 像素;雙脈沖Nd-YAG 激光器,最大脈沖頻率為15 Hz;最大脈沖激光能量為1200 mJ;常規(guī)工作能量為200 mJ;激光束持續(xù)時(shí)長為4 ns;激光波長為532 nm;脈沖激光片厚度為0.6 mm;最大測量范圍為400 mm×400 mm;示蹤粒子為聚酰胺顆粒.

    該套系統(tǒng)包括一個(gè)浸沒水中的雷體、兩個(gè)翼型支架、脈沖式激光器、同步器、圖像采集裝置、粒子播撒裝置等.雷體包含相機(jī)及眾多反射鏡和柱面鏡結(jié)構(gòu),其內(nèi)部詳細(xì)光學(xué)結(jié)構(gòu)布置如圖5 所示.雷體內(nèi)的兩臺相機(jī)為Dantec 公司的FlowSense 4M MKII 12位精度相機(jī)組,其使用高性能的逐行掃描行間CCD芯片,相機(jī)的分辨率為2048×2048 像素并且具有相同數(shù)目的存儲(chǔ)單元.兩臺相機(jī)的鏡頭參數(shù)各不相同,靠近前端的Camera1 鏡頭為50 mm f/1.8 的佳能鏡頭,靠近后端的Camera2 鏡頭為85 mm f/1.8 的佳能鏡頭.支柱1 中包含激光導(dǎo)光臂,支柱2 中包含相機(jī)線纜支柱結(jié)構(gòu).

    圖5 SPIV 系統(tǒng)中各元器件詳細(xì)布置圖[36]Fig.5 Detailed arrangement of the probes in the underwater SPIV system[36]

    圖6 SPIV 系統(tǒng)實(shí)物圖Fig.6 Stereoscopic particle image velocimetry system

    2.2 流場數(shù)據(jù)采集、測量與分析

    采用平均直徑為50μm 聚酰胺顆粒(PSP-50)作為示蹤粒子,密度為1.03 g/cm3,具有良好的跟隨性.由于拖曳水池尺度較大,如何在拖曳水池進(jìn)行大范圍粒子均勻播撒是一項(xiàng)技術(shù)挑戰(zhàn),本試驗(yàn)應(yīng)用一套自主設(shè)計(jì)的粒子播撒裝置進(jìn)行示蹤粒子播撒,播撒裝置如圖7 所示.裝置上方裝有攪拌電機(jī)的水箱可以存儲(chǔ)大量的水并均勻攪拌粒子,打開水箱下方的閥門后,粒子沿下方不同管道播撒到水池中,下方每條管道都開有若干小孔,其潛深可以在水下50 cm 范圍內(nèi)自由調(diào)節(jié),完全符合潛艇的潛深要求,在不同深度下均能保證潛艇模型周圍存在適量粒子.每航次試驗(yàn)結(jié)束后進(jìn)行示蹤粒子補(bǔ)充,通過相機(jī)觀測窗口進(jìn)行粒子濃度判別,保證每次測量過程中平均粒子濃度約為150 g/m3即每個(gè)查問區(qū)域具有15~30 個(gè)粒子,這樣能夠增大互相關(guān)分析的精確度.

    在試驗(yàn)過程中,脈沖激光照亮待測區(qū)域示蹤粒子,相機(jī)雙幀拍攝,采集頻率調(diào)至最大采樣頻率7.5 Hz,根據(jù)不同航速,調(diào)整激光(相機(jī))在相鄰兩幀之間的時(shí)間間隔,以保證在此間隔內(nèi)粒子最大位移不超過查詢區(qū)域的1/4,每次測量采集的圖片數(shù)控制在150~200 對之間以滿足數(shù)據(jù)收斂性的要求.試驗(yàn)過程如圖8 所示.

    圖7 粒子播撒裝置Fig.7 The particle seeding device

    圖8 工作中的SPIV 系統(tǒng)Fig.8 The SIV system in operation

    采用Dantec Studio 6.2 軟件對試驗(yàn)圖片進(jìn)行批量處理與后續(xù)分析.自適應(yīng)互相關(guān)算法獲得流場矢量信息,查詢區(qū)域大小設(shè)置為64×64 像素,32×32 像素,16×16 像素,重疊率為50%.鄰域驗(yàn)證設(shè)置為5×5 矩陣,通過比較矩陣內(nèi)的向量并啟用替換,可以有效過濾掉偽向量并用鄰域向量的中值進(jìn)行插值.通過實(shí)時(shí)監(jiān)測迭代步數(shù)對結(jié)果收斂性的影響,迭代次數(shù)選擇為10,檢測最少有效粒子數(shù)設(shè)置為10,最終得到的相鄰速度矢量間距為2.45 mm.

    吳鐵成[36]應(yīng)用該SPIV 系統(tǒng)在拖曳水池中進(jìn)行了大量船舶艉流場試驗(yàn),對該系統(tǒng)不確定度進(jìn)行分析,并對KCS 標(biāo)模多工況下伴流場進(jìn)行測量,與公開數(shù)據(jù)進(jìn)行對比,驗(yàn)證了該系統(tǒng)的準(zhǔn)確性.

    3 模型試驗(yàn)準(zhǔn)確性驗(yàn)證

    應(yīng)用泰勒水池的SUBOFF 實(shí)驗(yàn)數(shù)據(jù)作為參照,對本次模型試驗(yàn)及數(shù)值計(jì)算進(jìn)行驗(yàn)證.由于測試水池以及PIV 測試能力的限制,模型縮尺比略有不同,兩者主尺度參數(shù)見表2.測量航速由Fr相似換算得到,由于此時(shí)浸深較大,自由面影響較小,使用Fr相似來確定本試驗(yàn)航速不夠嚴(yán)謹(jǐn).受試驗(yàn)條件中模型長度及拖曳航速的限制,在今后的試驗(yàn)研究中應(yīng)盡量滿足雷諾數(shù)相似條件,使用更大的模型及航速.

    表2 模型差異Table 2 Model differences

    提取0.6R圓周范圍內(nèi)若干點(diǎn)上的潛艇槳盤面軸向伴流速度與標(biāo)準(zhǔn)實(shí)驗(yàn)值進(jìn)行對比.取點(diǎn)原則為周向間隔5?取一個(gè)點(diǎn),在0?,90?,180?的±10?鄰域內(nèi)角度間隔加密為2?.其中,周向角0?規(guī)定指向艇體正上方(十二點(diǎn)鐘位置),從艇艉前視,順時(shí)針旋轉(zhuǎn)為正方向.最終,軸向速度使用u/U作無量綱化表示,其中U為拖曳航速.軸向速度對比結(jié)果如圖9 所示.槳盤面上方速度值徑向分布如圖10 所示.

    圖9 槳盤面0.6R 處無量綱軸向速度值對比圖Fig.9 Comparison diagram of dimensionless axial velocity value at 0.6R of propeller disk

    圖10 槳盤面上方速度徑向分布圖Fig.10 Radial velocity distribution above the propeller disk

    PIV 試驗(yàn)得到的軸向速度值普遍小于泰勒水池實(shí)驗(yàn)值,將其歸結(jié)于模型大小引起的尺度效應(yīng).本試驗(yàn)采用的潛艇模型尺度較小,因此相同F(xiàn)r下的Re較小,導(dǎo)致邊界層稍厚,測量得到的速度要小一些.以泰勒水池的實(shí)驗(yàn)值作為參照計(jì)算本試驗(yàn)測量的誤差值,最大誤差在5%以內(nèi).對比發(fā)現(xiàn)CFD 計(jì)算的軸向速度值偏小,其原因?yàn)榇蟛糠志W(wǎng)格位于自由液面區(qū)域,艇體周圍網(wǎng)格只能達(dá)到四五百萬量級,另外模型太小,模型本身產(chǎn)生的流場特征也較小,經(jīng)過一些耗散實(shí)際測量到的量就會(huì)小很多,而數(shù)值模擬條件比較理想化,模擬出與試驗(yàn)環(huán)境相同的真實(shí)物理環(huán)境還存在差異,即CFD 結(jié)果過高估計(jì)了馬蹄渦的強(qiáng)度,使得馬蹄渦的影響區(qū)域增大,造成了槳盤面處偏低的軸向速度.通過以上流場對比結(jié)果來看,本次試驗(yàn)和數(shù)值計(jì)算結(jié)果與泰勒水池標(biāo)準(zhǔn)實(shí)驗(yàn)結(jié)果吻合良好,因此本試驗(yàn)結(jié)果能夠應(yīng)用于后續(xù)的研究.

    4 試驗(yàn)結(jié)果分析

    4.1 不同潛深下螺旋槳盤面速度場分析

    圖11 給出Fr=0.3~0.7 范圍內(nèi)不同潛深工況下槳盤面軸向速度場,左為1.5D潛深結(jié)果,中間為2.5D潛深結(jié)果,右為4D潛深結(jié)果.其中空白區(qū)域?qū)?yīng)槳軸遮擋處的無數(shù)據(jù)區(qū)域.

    從圖11 中可以看出,本次PIV 試驗(yàn)大致描繪出了潛艇槳盤面處流場的幾項(xiàng)重要特征.首先,速度值呈現(xiàn)出沿徑向連續(xù)增大的圓圈式形狀;其次,4 個(gè)十字舵角度(0?,90?,180?,270?) 最外側(cè)的速度等值線向外凸出,速度等值線0.8~0.95 范圍表現(xiàn)出了十字舵后端由于速度削減效應(yīng)產(chǎn)生的低速區(qū)的存在;另外,內(nèi)側(cè)的速度等值線朝向槳軸中心凹陷,對應(yīng)速度等值線范圍為0.4~0.8,這是由于主附體交接部位馬蹄渦作用造成動(dòng)量交換,使得中心位置動(dòng)量增大,兩側(cè)動(dòng)量降低,從而形成近似“V”型的速度等值線結(jié)構(gòu).在徑向上表現(xiàn)為馬蹄渦將高動(dòng)量的流體卷攜至徑向內(nèi)部,低動(dòng)量的流體卷攜至徑向外部.大潛深工況下,隨著Fr增大,流場結(jié)構(gòu)特征基本一致.

    在1.5D潛深工況下,隨著Fr增大,潛艇尾流場依舊存在上述幾項(xiàng)重要特征.進(jìn)一步比較Fr對軸向速度的影響,可以發(fā)現(xiàn)隨著Fr增大潛艇上方外側(cè)速度等值線有不斷向下的趨勢,最外層0.95 等值線輪廓位于Y/R=1 處.與此同時(shí),上方等值線整體趨于扁平化,由指揮臺圍殼脫落的馬蹄渦而形成的速度等值線凸起越來越不明顯.這是由于潛艇近水面航行時(shí),不同F(xiàn)r下興波高度不同,自由面與艇體伴流相互作用.

    圖11 不同潛深工況下槳盤面軸向速度云圖(左側(cè)為1.5D 潛深工況,中間為2.5D 潛深工況,右側(cè)為4D 潛深工況)Fig.11 Axial velocity of propeller disk under different submergence conditions(1.5D submergence depths on the left,2.5D submergence depths in the middle and 4D submergence depths on the right

    圖11 不同潛深工況下槳盤面軸向速度云圖(左側(cè)為1.5D 潛深工況,中間為2.5D 潛深工況,右側(cè)為4D 潛深工況)(續(xù))Fig.11 Axial velocity of propeller disk under different submergence conditions(1.5D submergence depths on the left,2.5D submergence depths in the middle and 4D submergence depths on the right(continued)

    不同潛深工況下潛艇尾流場流動(dòng)特征基本類似.對比4D潛深、2.5D潛深與1.5D潛深槳盤面軸向速度分布,可以發(fā)現(xiàn)在低Fr時(shí),1.5D潛深的上部分速度分布呈現(xiàn)較為飽和的半圓形分布,4D潛深的上部分速度有壓縮現(xiàn)象,2.5D潛深介于1.5D潛深和4D潛深之間.隨著Fr增大,自由面主要影響流場的上半?yún)^(qū)域,三種潛深工況下的速度均有壓縮趨勢,而1.5D潛深的壓縮現(xiàn)象更為明顯,2.5D潛深次之.在試驗(yàn)過程中,觀察到高速工況對應(yīng)的興波越發(fā)明顯,速度等值線體現(xiàn)出的特征應(yīng)與興波密切相關(guān),自由面產(chǎn)生了一種擠壓作用使得上半?yún)^(qū)域流場結(jié)構(gòu)整體下移.由于2.5D潛深工況整體特性介于1.5D潛深和4D潛深之間,因此下文僅對1.5D潛深和4D潛深工況展開分析.

    4.2 不同航速及潛深下潛艇尾部脈動(dòng)速度場分析

    由一系列的瞬時(shí)速度云圖分析得到脈動(dòng)速度,并進(jìn)行無量綱化處理,得到uRMS/U脈動(dòng)速度云圖.4D潛深與1.5D潛深工況不同航速下脈動(dòng)速度云圖分別如圖12 和圖13 所示.

    從圖12 和圖13 中可以看出,脈動(dòng)速度云圖與軸向速度云圖較為相似,為一圈圈向外延伸的等值線結(jié)構(gòu),不同點(diǎn)在于向外的脈動(dòng)速度值不斷降低,在外流場達(dá)到0.02~0.03 的量級,同時(shí)在槳盤面內(nèi)半徑區(qū)域以及尾舵梢部附近即Y/R=1 出現(xiàn)較大的脈動(dòng)速度,這對應(yīng)于馬蹄渦及梢渦的位置.

    對比不同潛深下的脈動(dòng)速度值,可以看出1.5D潛深下的脈動(dòng)速度云圖受自由面影響較大,在局部上體現(xiàn)出脈動(dòng)速度極大值點(diǎn),同時(shí)隨著Fr增大,槳盤面上方速度等值線明顯下移,表明受自由面的影響,較大的脈動(dòng)速度結(jié)構(gòu)整體下移,體現(xiàn)出與軸向速度等值線相似的規(guī)律.

    圖12 1.5D 潛深不同航速下脈動(dòng)速度云圖Fig.12 Fluctuating velocity at different speeds of 1.5D depth

    圖13 4D 潛深不同航速下脈動(dòng)速度云圖Fig.13 Fluctuating velocity at different speeds of 4D depth

    進(jìn)一步提取0.6R范圍內(nèi)的脈動(dòng)速度值,如圖14所示,由于槳軸反光區(qū)域遮擋,60?~120?范圍數(shù)據(jù)缺失.從圖中可以看出不同潛深工況下槳盤面脈動(dòng)速度隨Fr增大所體現(xiàn)出的趨勢是相同的.首先,在每個(gè)90?范圍內(nèi),脈動(dòng)速度呈現(xiàn)出先下降后上升的趨勢,即對應(yīng)每個(gè)尾舵后方都會(huì)出現(xiàn)較大的速度脈動(dòng),具有一定周期性.另外上半部分(0?~90?范圍)的脈動(dòng)速度值略大于下半部分(90?~180?范圍)的脈動(dòng)速度值,這就表明槳軸上方的流場會(huì)受到指揮臺尾流的影響,從而產(chǎn)生更大的脈動(dòng)速度.相比于4D潛深工況,1.5D潛深時(shí)局部上出現(xiàn)脈動(dòng)速度極大值點(diǎn),且不同F(xiàn)r間的脈動(dòng)速度差值也有所增大,因此可知自由面的引入增大了流動(dòng)的脈動(dòng)特性.

    4.3 對于PIV 測試結(jié)果的CFD 補(bǔ)充與擴(kuò)展

    圖14 槳盤面0.6R 處脈動(dòng)速度對比Fig.14 Fluctuating velocity at 0.6R of propeller disk

    由于PIV 測試難以實(shí)現(xiàn)對于潛艇興波的精確探測,僅能夠?qū)Υ笾虏ㄐ芜M(jìn)行觀察以完成粗略的對比,而CFD 計(jì)算可以精確地模擬出波浪的細(xì)節(jié)信息.本小節(jié)基于CFD 模擬結(jié)果,對潛艇近水面航行所興起的波浪進(jìn)行量化分析,并給出中縱切面的波浪高度,從而在興波的角度對潛艇繞流場隨航速與潛深的變化規(guī)律進(jìn)行解釋.

    通過數(shù)值模擬發(fā)現(xiàn),在4D潛深時(shí),興波高度隨航速增大不斷提高,但不同航速對應(yīng)的波浪峰值差距僅為毫米量級,完全可以忽略自由表面的影響,因此下文僅對1.5D潛深工況進(jìn)行分析.

    從圖15 可以看出,不同航速的潛艇興波具有以下規(guī)律:首先潛艇首部后方存在一個(gè)波峰,且航速越高,波幅越大.由于波長隨Fr增大而增大,在航速增大的情況下,興波波峰的起始位置會(huì)向艇后方移動(dòng).同時(shí)在潛艇尾肩部存在一波谷,隨著航速變大,波谷位置同樣向后移.因此在槳盤面位置總體表現(xiàn)為隨航速增大而降低的液面高度.

    圖15 1.5D 潛深不同速度下潛艇中縱剖面波浪高度對比Fig.15 Comparison of wave height of longitudinal section in centre plane of submarine at different speeds of 1.5D

    給出不同潛深及不同航速下的潛艇中縱切面速度云圖,如圖16 所示,左為1.5D潛深工況,右為4D潛深工況.

    圖16 不同航速下潛艇中縱剖面速度云圖(左為1.5D,右為4D)Fig.16 Velocity nephogram of longitudinal section in centre plane of submarine at different speeds(1.5D submergence depths on the left and 4D submergence depths on the right)

    通過對比不同潛深工況下速度云圖發(fā)現(xiàn),自由液面存在對流場的最大影響體現(xiàn)在上半?yún)^(qū)域.同一航速下,近水面航行時(shí),艇體與自由面間流場速度明顯增大,特別在波谷區(qū)域(對應(yīng)槳盤面區(qū)域),流場速度明顯提升.這可以由文丘里效應(yīng)合理地進(jìn)行解釋,即受限流動(dòng)在經(jīng)過縮小的過流斷面時(shí),流體流速會(huì)增大.由前文可知,隨著Fr提升,槳盤面處的液面高度越來越低,這就導(dǎo)致了槳盤面位置出現(xiàn)更大的流體速度,而流速的增大勢必會(huì)導(dǎo)致伴流場結(jié)構(gòu)的變化.

    5 結(jié)論

    本文應(yīng)用SPIV 系統(tǒng)對潛艇尾部流場進(jìn)行測試,同時(shí)輔以數(shù)值模擬對試驗(yàn)結(jié)果進(jìn)行補(bǔ)充與擴(kuò)展.通過對不同潛深工況不同航速的分析得到以下結(jié)論:

    (1)在1.5D潛深工況下,隨著Fr增大,由于自由面產(chǎn)生的擠壓作用使得槳盤面上方速度等值線有不斷向下的趨勢,上方等值線整體趨于扁平化,由指揮臺圍殼脫落的馬蹄渦形成的速度等值線凸起越來越不明顯.

    (2)每個(gè)尾舵后方都會(huì)出現(xiàn)較大的脈動(dòng)速度,且具有一定周期性,在周期內(nèi),整體脈動(dòng)速度呈現(xiàn)先下降后上升的趨勢.受指揮臺尾流影響,槳盤面上方流場體現(xiàn)出更大的脈動(dòng)速度.自由面會(huì)增大流體的脈動(dòng)特性,在1.5D潛深工況的脈動(dòng)速度云圖中存在局部脈動(dòng)速度極大值點(diǎn),同時(shí)隨著Fr增大,槳盤面上方0.04 以上的速度等值線明顯下移,表明受自由面的影響,較大的脈動(dòng)速度結(jié)構(gòu)整體下移.

    (3) 由于波長隨Fr增大而變大,隨著航速增大,興波波峰的起始位置向艇后方移動(dòng).同時(shí)在潛艇尾肩部存在一波谷,隨著航速增大大,波谷位置同樣向后移.因此在槳盤面處總體表現(xiàn)為自由面高度隨航速增大而降低.

    (4)自由面存在時(shí),艇體與自由面間流場速度明顯增大,特別是波谷區(qū)域,也恰好對應(yīng)槳盤面區(qū)域,流場速度明顯提升.隨著Fr增大,槳盤面處的液面高度越來越低,這就導(dǎo)致了槳盤面位置出現(xiàn)更大的流體速度,即造成了槳盤面伴流場擠壓現(xiàn)象的發(fā)生.

    猜你喜歡
    盤面尾流航速
    VLCC在波浪中的航速優(yōu)化與能效優(yōu)化分析
    提升全回轉(zhuǎn)港作拖輪航速的有效途徑
    水上消防(2022年1期)2022-06-16 08:06:56
    淺談制動(dòng)盤面端跳動(dòng)測量
    低速水面目標(biāo)航速精度分析及精確解算
    飛機(jī)尾流的散射特性與探測技術(shù)綜述
    數(shù)獨(dú)競猜
    數(shù)獨(dú)競猜
    基于CFD的波浪滑翔機(jī)航速預(yù)測
    錐形流量計(jì)尾流流場分析
    水面艦船風(fēng)尾流效應(yīng)減弱的模擬研究
    国产成人91sexporn| 少妇被粗大猛烈的视频| 成人免费观看视频高清| 在线观看免费视频网站a站| 女的被弄到高潮叫床怎么办| 国产亚洲精品久久久com| 黄色毛片三级朝国网站| 久久久久人妻精品一区果冻| 啦啦啦在线观看免费高清www| 日本vs欧美在线观看视频| 日韩av在线免费看完整版不卡| 99国产综合亚洲精品| 久久毛片免费看一区二区三区| 一个人免费看片子| 伊人久久国产一区二区| 久久久精品区二区三区| 熟女人妻精品中文字幕| 免费观看无遮挡的男女| 久久久久久久精品精品| 极品人妻少妇av视频| 男人操女人黄网站| 亚洲精品视频女| 夫妻性生交免费视频一级片| 深夜精品福利| 80岁老熟妇乱子伦牲交| 97在线人人人人妻| 99热全是精品| 亚洲国产精品成人久久小说| 午夜老司机福利剧场| 久久99热这里只频精品6学生| 国产亚洲最大av| 亚洲一级一片aⅴ在线观看| 久久精品国产自在天天线| 久久99精品国语久久久| 人成视频在线观看免费观看| 欧美日韩综合久久久久久| 国产成人精品福利久久| 色5月婷婷丁香| 久久精品国产鲁丝片午夜精品| 亚洲欧美日韩另类电影网站| 日韩免费高清中文字幕av| 精品卡一卡二卡四卡免费| 久久 成人 亚洲| 欧美日韩视频高清一区二区三区二| 国产成人精品婷婷| 亚洲av福利一区| 日韩中文字幕视频在线看片| 久久久久久伊人网av| 视频区图区小说| 婷婷色综合大香蕉| 欧美bdsm另类| 五月伊人婷婷丁香| 国产黄频视频在线观看| 亚洲av中文av极速乱| 久久久久久久久久久久大奶| 在线天堂最新版资源| 精品酒店卫生间| 亚洲中文av在线| 少妇高潮的动态图| www.av在线官网国产| 久久精品国产亚洲av涩爱| 免费在线观看完整版高清| 80岁老熟妇乱子伦牲交| 青春草视频在线免费观看| 日本欧美视频一区| 91成人精品电影| 久久久久久久国产电影| 三级国产精品片| 国产乱来视频区| 三级国产精品片| 久久精品久久久久久久性| 韩国av在线不卡| www日本在线高清视频| 99九九在线精品视频| 91在线精品国自产拍蜜月| 成人漫画全彩无遮挡| 男女高潮啪啪啪动态图| 草草在线视频免费看| 日本与韩国留学比较| 2022亚洲国产成人精品| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成色77777| 欧美激情极品国产一区二区三区 | 中国三级夫妇交换| 欧美精品一区二区大全| 90打野战视频偷拍视频| 国产精品.久久久| 99国产精品免费福利视频| 在线精品无人区一区二区三| 国产成人一区二区在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人av激情在线播放| 成人18禁高潮啪啪吃奶动态图| 国产精品欧美亚洲77777| 九色成人免费人妻av| 精品久久国产蜜桃| 多毛熟女@视频| 久久影院123| 搡老乐熟女国产| av天堂久久9| 丰满迷人的少妇在线观看| 久久人人爽av亚洲精品天堂| 天堂8中文在线网| av在线播放精品| 亚洲色图综合在线观看| 美女福利国产在线| 日韩欧美精品免费久久| 亚洲四区av| 中文乱码字字幕精品一区二区三区| 亚洲精品第二区| 超碰97精品在线观看| 在线亚洲精品国产二区图片欧美| 99久久人妻综合| 免费少妇av软件| 亚洲人成网站在线观看播放| 91精品三级在线观看| 黄色配什么色好看| 国产成人精品福利久久| 黑丝袜美女国产一区| 国产精品久久久久久久电影| 满18在线观看网站| 日韩av在线免费看完整版不卡| 大香蕉久久成人网| 哪个播放器可以免费观看大片| 视频区图区小说| 成人国产av品久久久| 亚洲欧美一区二区三区黑人 | 最近中文字幕2019免费版| 亚洲人成77777在线视频| 18+在线观看网站| kizo精华| 国产午夜精品一二区理论片| 两个人免费观看高清视频| 少妇被粗大的猛进出69影院 | 高清黄色对白视频在线免费看| 亚洲欧美精品自产自拍| 精品熟女少妇av免费看| 男的添女的下面高潮视频| 最近最新中文字幕大全免费视频 | 国产免费现黄频在线看| 国产精品久久久久久久电影| 国产永久视频网站| 欧美日韩成人在线一区二区| 成人国产麻豆网| 18禁裸乳无遮挡动漫免费视频| 少妇的逼好多水| 欧美日韩视频高清一区二区三区二| 大片免费播放器 马上看| 性色avwww在线观看| 国产精品一区www在线观看| 啦啦啦啦在线视频资源| a级毛色黄片| 免费看av在线观看网站| 丝瓜视频免费看黄片| 成人综合一区亚洲| 久久毛片免费看一区二区三区| 高清视频免费观看一区二区| 99re6热这里在线精品视频| 欧美人与善性xxx| 色94色欧美一区二区| 欧美日韩一区二区视频在线观看视频在线| 亚洲一码二码三码区别大吗| 最新的欧美精品一区二区| 成人黄色视频免费在线看| 国产成人91sexporn| 大码成人一级视频| 精品一区在线观看国产| 啦啦啦中文免费视频观看日本| 免费观看无遮挡的男女| 大香蕉97超碰在线| a级毛片黄视频| 青青草视频在线视频观看| 丰满迷人的少妇在线观看| 各种免费的搞黄视频| 中文精品一卡2卡3卡4更新| 男女国产视频网站| 国产1区2区3区精品| 好男人视频免费观看在线| 一本久久精品| 国产 精品1| 99热全是精品| 亚洲精品国产av蜜桃| 丁香六月天网| 考比视频在线观看| 岛国毛片在线播放| 亚洲第一区二区三区不卡| 午夜福利网站1000一区二区三区| 99香蕉大伊视频| 亚洲一级一片aⅴ在线观看| 欧美国产精品一级二级三级| 青青草视频在线视频观看| 国产精品久久久久久av不卡| 成年人午夜在线观看视频| 如何舔出高潮| 青春草视频在线免费观看| 亚洲精品美女久久久久99蜜臀 | av有码第一页| 母亲3免费完整高清在线观看 | a级片在线免费高清观看视频| 国产精品.久久久| 国产免费一级a男人的天堂| 丝袜人妻中文字幕| 日韩精品有码人妻一区| 国产精品三级大全| 大片免费播放器 马上看| 少妇熟女欧美另类| 亚洲少妇的诱惑av| 亚洲精华国产精华液的使用体验| 波野结衣二区三区在线| 国产黄色免费在线视频| 亚洲国产av影院在线观看| 欧美国产精品一级二级三级| 99九九在线精品视频| 啦啦啦啦在线视频资源| 一级爰片在线观看| 99国产精品免费福利视频| 亚洲一码二码三码区别大吗| 免费观看a级毛片全部| 黑人高潮一二区| 视频中文字幕在线观看| 有码 亚洲区| 9色porny在线观看| 一级片'在线观看视频| 黄片播放在线免费| 最后的刺客免费高清国语| 99热6这里只有精品| 成年av动漫网址| 日韩欧美一区视频在线观看| 捣出白浆h1v1| 国产熟女欧美一区二区| 国产精品成人在线| 中国三级夫妇交换| 在线观看一区二区三区激情| 欧美激情 高清一区二区三区| 免费高清在线观看日韩| 成人毛片60女人毛片免费| 中文字幕人妻熟女乱码| 国产男女内射视频| 人人妻人人添人人爽欧美一区卜| 日韩欧美精品免费久久| 亚洲欧美日韩另类电影网站| av在线app专区| 免费在线观看完整版高清| 日本免费在线观看一区| 只有这里有精品99| 亚洲成国产人片在线观看| 青青草视频在线视频观看| 成人无遮挡网站| 黄色 视频免费看| 青春草视频在线免费观看| 不卡视频在线观看欧美| 中文字幕人妻丝袜制服| 国产精品国产三级专区第一集| 色吧在线观看| av国产精品久久久久影院| 欧美国产精品va在线观看不卡| 国产一级毛片在线| 80岁老熟妇乱子伦牲交| 熟女电影av网| 国产 一区精品| 夫妻午夜视频| 国产精品熟女久久久久浪| 日日啪夜夜爽| 看免费成人av毛片| 91aial.com中文字幕在线观看| 大片电影免费在线观看免费| 亚洲欧美日韩卡通动漫| 侵犯人妻中文字幕一二三四区| 精品亚洲乱码少妇综合久久| 蜜桃在线观看..| 伊人亚洲综合成人网| 国产极品粉嫩免费观看在线| 欧美少妇被猛烈插入视频| av电影中文网址| 性色av一级| 国产午夜精品一二区理论片| 欧美xxⅹ黑人| 人人妻人人澡人人看| 激情五月婷婷亚洲| 午夜福利影视在线免费观看| 熟妇人妻不卡中文字幕| 欧美激情 高清一区二区三区| 在线天堂最新版资源| 亚洲熟女精品中文字幕| 大片电影免费在线观看免费| 国产成人a∨麻豆精品| 免费看av在线观看网站| 日韩电影二区| 99国产综合亚洲精品| 亚洲熟女精品中文字幕| 在线观看国产h片| 夜夜骑夜夜射夜夜干| 内地一区二区视频在线| 另类亚洲欧美激情| 国产av精品麻豆| 成人黄色视频免费在线看| 久久午夜综合久久蜜桃| 国产精品一区www在线观看| 乱人伦中国视频| 精品少妇黑人巨大在线播放| 午夜日本视频在线| 爱豆传媒免费全集在线观看| 久久午夜福利片| 国产欧美日韩一区二区三区在线| av视频免费观看在线观看| 天堂中文最新版在线下载| 最近最新中文字幕免费大全7| 卡戴珊不雅视频在线播放| 国产亚洲一区二区精品| 久久久精品免费免费高清| 永久网站在线| av片东京热男人的天堂| 卡戴珊不雅视频在线播放| 另类精品久久| 亚洲国产精品999| 精品国产露脸久久av麻豆| 亚洲av.av天堂| 99九九在线精品视频| av黄色大香蕉| 99精国产麻豆久久婷婷| 亚洲国产av新网站| 91精品伊人久久大香线蕉| 亚洲美女搞黄在线观看| 成人亚洲精品一区在线观看| 看非洲黑人一级黄片| 人人澡人人妻人| 久久99热6这里只有精品| 久久av网站| 免费高清在线观看视频在线观看| 成人影院久久| 久久久久国产精品人妻一区二区| 亚洲丝袜综合中文字幕| 午夜福利网站1000一区二区三区| 女人精品久久久久毛片| 一区二区三区乱码不卡18| videos熟女内射| av.在线天堂| 精品久久蜜臀av无| 中文字幕精品免费在线观看视频 | av播播在线观看一区| 天堂8中文在线网| 国产精品久久久久久精品古装| 国产探花极品一区二区| 观看av在线不卡| 91在线精品国自产拍蜜月| 亚洲精品国产av蜜桃| 成人免费观看视频高清| 日韩不卡一区二区三区视频在线| 色视频在线一区二区三区| 亚洲av综合色区一区| 成年人免费黄色播放视频| 91aial.com中文字幕在线观看| 搡老乐熟女国产| 午夜激情av网站| 成人综合一区亚洲| 国产女主播在线喷水免费视频网站| 99久久中文字幕三级久久日本| 亚洲国产最新在线播放| 韩国精品一区二区三区 | 少妇的逼好多水| av女优亚洲男人天堂| 秋霞伦理黄片| 伊人亚洲综合成人网| 中文精品一卡2卡3卡4更新| 成年人午夜在线观看视频| 日本欧美视频一区| 午夜福利在线观看免费完整高清在| 国产一区亚洲一区在线观看| 婷婷色综合www| av在线老鸭窝| 午夜福利在线观看免费完整高清在| 亚洲国产欧美在线一区| 女性被躁到高潮视频| 久热这里只有精品99| 日韩人妻精品一区2区三区| 久久影院123| 欧美成人精品欧美一级黄| 九草在线视频观看| 久久这里只有精品19| 国语对白做爰xxxⅹ性视频网站| 亚洲国产精品国产精品| 老司机影院毛片| 中国三级夫妇交换| 热99国产精品久久久久久7| 欧美97在线视频| 午夜福利网站1000一区二区三区| 国产激情久久老熟女| 国产精品嫩草影院av在线观看| 男人添女人高潮全过程视频| 午夜福利网站1000一区二区三区| 18禁裸乳无遮挡动漫免费视频| 亚洲四区av| 成年动漫av网址| 视频在线观看一区二区三区| 亚洲国产精品一区三区| 亚洲精品中文字幕在线视频| 成人免费观看视频高清| 亚洲av男天堂| 亚洲欧美中文字幕日韩二区| 高清视频免费观看一区二区| 婷婷色综合大香蕉| 一区二区三区乱码不卡18| 一级爰片在线观看| 亚洲精品,欧美精品| 国产成人免费无遮挡视频| 国产色爽女视频免费观看| 天天躁夜夜躁狠狠久久av| 青青草视频在线视频观看| 高清毛片免费看| 欧美xxⅹ黑人| 国产亚洲最大av| 国产爽快片一区二区三区| 香蕉丝袜av| 男女边吃奶边做爰视频| 丰满饥渴人妻一区二区三| 黄色一级大片看看| 桃花免费在线播放| 51国产日韩欧美| 久久鲁丝午夜福利片| 一本色道久久久久久精品综合| a 毛片基地| 日韩,欧美,国产一区二区三区| 少妇人妻 视频| 国产片特级美女逼逼视频| 久久久久国产网址| 大码成人一级视频| 婷婷色av中文字幕| 永久网站在线| 国产精品女同一区二区软件| 国产在线一区二区三区精| 亚洲精品自拍成人| 少妇猛男粗大的猛烈进出视频| 精品一品国产午夜福利视频| 成人黄色视频免费在线看| av不卡在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产淫语在线视频| 99香蕉大伊视频| 亚洲国产精品成人久久小说| 欧美日本中文国产一区发布| 国产激情久久老熟女| 久久久久久人妻| 亚洲国产精品一区二区三区在线| 久久人人爽av亚洲精品天堂| 国产又爽黄色视频| av福利片在线| 久久精品夜色国产| 大陆偷拍与自拍| 两个人看的免费小视频| 免费黄频网站在线观看国产| 亚洲精品视频女| 99久久精品国产国产毛片| 亚洲在久久综合| 在线观看三级黄色| 99久久人妻综合| 91国产中文字幕| 欧美国产精品va在线观看不卡| 亚洲精品美女久久av网站| 国产一区亚洲一区在线观看| 久久久久网色| 18禁观看日本| 欧美成人午夜精品| 亚洲国产日韩一区二区| 成人国语在线视频| 性色avwww在线观看| 18禁动态无遮挡网站| 日韩成人伦理影院| 丝袜人妻中文字幕| 一级片'在线观看视频| 激情五月婷婷亚洲| 国内精品宾馆在线| 亚洲少妇的诱惑av| 亚洲第一av免费看| 18+在线观看网站| 亚洲精品自拍成人| 少妇的逼好多水| 免费黄色在线免费观看| 老司机影院毛片| 日韩大片免费观看网站| 国产1区2区3区精品| 免费看不卡的av| 亚洲精品,欧美精品| 一本大道久久a久久精品| 亚洲三级黄色毛片| 又大又黄又爽视频免费| 午夜老司机福利剧场| 另类精品久久| 成年美女黄网站色视频大全免费| 日韩大片免费观看网站| 看免费av毛片| 精品亚洲乱码少妇综合久久| 天天躁夜夜躁狠狠久久av| 精品久久国产蜜桃| 午夜福利影视在线免费观看| 国产在线一区二区三区精| 久久久久久久久久久免费av| 久久这里只有精品19| 黑人猛操日本美女一级片| 女性被躁到高潮视频| 巨乳人妻的诱惑在线观看| 肉色欧美久久久久久久蜜桃| 日本vs欧美在线观看视频| 王馨瑶露胸无遮挡在线观看| 人妻少妇偷人精品九色| 久久女婷五月综合色啪小说| 色婷婷av一区二区三区视频| 精品久久久精品久久久| a级毛色黄片| 欧美人与善性xxx| 国产极品天堂在线| 狂野欧美激情性xxxx在线观看| 亚洲av综合色区一区| 18禁国产床啪视频网站| av又黄又爽大尺度在线免费看| 中文字幕精品免费在线观看视频 | 国产成人午夜福利电影在线观看| 国产成人精品无人区| 亚洲精品乱久久久久久| 夜夜骑夜夜射夜夜干| 在线观看一区二区三区激情| 国产欧美日韩综合在线一区二区| 少妇猛男粗大的猛烈进出视频| 久久久久久久国产电影| 久久久久精品人妻al黑| 日韩精品有码人妻一区| 精品少妇内射三级| 国产在线一区二区三区精| 亚洲精品视频女| 成人漫画全彩无遮挡| 啦啦啦啦在线视频资源| 国产探花极品一区二区| 精品一区二区三区四区五区乱码 | 激情五月婷婷亚洲| 最近的中文字幕免费完整| www.av在线官网国产| 国产精品三级大全| 999精品在线视频| 熟女电影av网| 大片免费播放器 马上看| 亚洲第一av免费看| 国产淫语在线视频| 夫妻性生交免费视频一级片| av免费在线看不卡| 色婷婷久久久亚洲欧美| 国产熟女欧美一区二区| 国产精品免费大片| av黄色大香蕉| 久久免费观看电影| 亚洲av免费高清在线观看| 欧美日韩亚洲高清精品| 22中文网久久字幕| 两个人看的免费小视频| 亚洲国产成人一精品久久久| 国产精品成人在线| 一本久久精品| 亚洲人成77777在线视频| 又粗又硬又长又爽又黄的视频| 亚洲精品一区蜜桃| 成人毛片60女人毛片免费| 国产高清国产精品国产三级| 日韩一区二区视频免费看| 91国产中文字幕| 视频中文字幕在线观看| 亚洲,欧美精品.| 一区二区av电影网| 狠狠精品人妻久久久久久综合| 国产免费福利视频在线观看| 久久精品国产综合久久久 | 国产成人午夜福利电影在线观看| 亚洲av福利一区| 亚洲国产欧美日韩在线播放| 免费人妻精品一区二区三区视频| 精品少妇内射三级| 另类亚洲欧美激情| 夫妻午夜视频| 国产69精品久久久久777片| 中文字幕人妻丝袜制服| 亚洲综合精品二区| 午夜免费男女啪啪视频观看| 五月伊人婷婷丁香| 2021少妇久久久久久久久久久| 国产精品久久久av美女十八| 亚洲av男天堂| 热re99久久精品国产66热6| 黄片无遮挡物在线观看| 亚洲综合精品二区| 国产精品久久久久久久电影| 久久久久久久精品精品| 欧美+日韩+精品| 爱豆传媒免费全集在线观看| 最近2019中文字幕mv第一页| 亚洲在久久综合| 久久久国产一区二区| 国产日韩欧美在线精品| 国产欧美日韩综合在线一区二区| 欧美 日韩 精品 国产| 男女免费视频国产| 18禁在线无遮挡免费观看视频| 女的被弄到高潮叫床怎么办| 18+在线观看网站| 七月丁香在线播放| av不卡在线播放| 99国产精品免费福利视频| 最近中文字幕高清免费大全6| 日韩视频在线欧美| 久久久久视频综合| 久久精品久久久久久久性| av播播在线观看一区| 国产精品熟女久久久久浪| 在线看a的网站| 成人无遮挡网站| 一级毛片黄色毛片免费观看视频| 精品一区在线观看国产| 国产精品人妻久久久久久| 亚洲精品乱码久久久久久按摩|