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

    表面 PIV 在潛航體興波伴流場測量中的應(yīng)用

    2016-11-15 14:26:05周文進(jìn)蔣小勤王建中方頻捷海軍工程大學(xué)理學(xué)院湖北武漢430033
    艦船科學(xué)技術(shù) 2016年6期
    關(guān)鍵詞:興波內(nèi)波水面

    周文進(jìn),蔣小勤,王建中,方頻捷(海軍工程大學(xué) 理學(xué)院,湖北 武漢 430033)

    表面 PIV 在潛航體興波伴流場測量中的應(yīng)用

    周文進(jìn),蔣小勤,王建中,方頻捷(海軍工程大學(xué) 理學(xué)院,湖北武漢 430033)

    介紹潛航體興波伴流場的表面 PIV 測量方法,對(duì)表面 PIV 的標(biāo)定,示蹤粒子的選擇,背景干擾流場的消除方法以及表面 PIV 測量精度等問題開展分析討論。實(shí)驗(yàn)表明:運(yùn)用表面 PIV 技術(shù)可靈敏地探測到潛航體微弱的 V 型興波伴流場,該興波伴流場可以用潛航體的‘源-匯'效應(yīng)解釋。本文采用普通 CCD 數(shù)碼攝像機(jī)以及普通光源攝取表面粒子圖像,與傳統(tǒng)的激光 PIV 設(shè)備相比,不僅經(jīng)濟(jì)而且使用方便,測量結(jié)果精確可靠,運(yùn)用表面 PIV 技術(shù)測量表面微弱流場的精度可達(dá) 0.01 pixel。

    潛航體;興波伴流;表面PIV

    0 引 言

    潛艇的非聲學(xué)探測已成為一個(gè)熱門課題,其中基于潛艇水動(dòng)力學(xué)尾跡特征的探測是其重要的分支,相應(yīng)的探測手段孕育而生,如海洋內(nèi)波影響海面毛細(xì)波的粗糙度,從而導(dǎo)致微波散射的變化,據(jù)此形成了運(yùn)用合成孔徑雷達(dá)(SAR)從空中捕捉內(nèi)波的探測技術(shù)。為了準(zhǔn)確分析評(píng)估潛航體激發(fā)的內(nèi)波對(duì)海洋表面的調(diào)制作用,在實(shí)驗(yàn)室中開展?jié)摵襟w水面效應(yīng)的定性或定量的研究,是一個(gè)不可或缺的重要環(huán)節(jié)[1-3]。

    20 世紀(jì) 90 年代末,馬暉揚(yáng)等[4]在水-鹽水、水-柴油 2 種分層的水槽中進(jìn)行拖曳球的實(shí)驗(yàn),展示了密度突躍層和水面的變形。常煜等[5]基于 RANS 方程和VOF 方法的非定常粘性數(shù)值方法,模擬了均勻流體和2 層流體中運(yùn)動(dòng)潛艇激發(fā)內(nèi)波的水面效應(yīng)。張效慈[6]利用 Tuck 方法,數(shù)值計(jì)算了幾種潛艇運(yùn)動(dòng)時(shí)潛艇尾跡傳播到海面時(shí)的映波幅值。在 PIV 實(shí)驗(yàn)方面,S.I.Voropayev等[7]在雷諾數(shù) Re>104條件下測量表面射流和拖曳尾流的表面旋渦,給出了旋渦大小的渦量判斷條件;秦朝峰等[8]利用 PIV 技術(shù),在分層流體自由面下 5 mm 處,探測并計(jì)算了流場的分布規(guī)律、散度和旋度等;Rottman 和 Broutman[9]在水中分層處撒聚合物粒子,用激光片光照明,利用 PIV 技術(shù)對(duì)內(nèi)波場進(jìn)行研究;姚志崇[10]采用撒聚合物粒子的方法,但使用的是普通白熾燈光源,從垂向方向觀測,對(duì)小球的尾流激發(fā)內(nèi)波在流體內(nèi)部的演化規(guī)律進(jìn)行研究。從實(shí)驗(yàn)上來說,大多數(shù)學(xué)者利用 PIV 技術(shù)很好地得到了流體內(nèi)部流場的演化規(guī)律,但并沒有直接測量水面的流場;陳祥瑞[10]利用線陣 CCD 相機(jī)直接拍攝拖曳潛航體興波,反射光的強(qiáng)度對(duì)水面變形非常敏感,通過線陣 CCD 相機(jī)對(duì)水面掃描獲取的是水面各處反射光的強(qiáng)度信息,定性地得到了水面興波紋理圖像及其隨速度的演化規(guī)律,但未能給出水面興波波高、表面流場等定量測量結(jié)果[11]。

    無論是在密度均勻流體還是分層流體中,潛航體的水面效應(yīng)都屬于微弱流場[12-13],線性分層流內(nèi)波理論預(yù)言潛航體激發(fā)內(nèi)波運(yùn)動(dòng)會(huì)在水面上產(chǎn)生輻聚輻散流,該輻聚輻散流與表面微尺度波的相互作用是微波雷達(dá)(SAR)探測內(nèi)波的物理基礎(chǔ)[14-15]。從實(shí)驗(yàn)研究的角度目前至少有 2 個(gè)問題待解決:1)在均勻流與分層流中潛航體的表面興波伴流場有無明顯的差異;2)如果存在明顯差異,就需要定量的測量,給出伴流場的強(qiáng)度與分層流、潛航體特征參數(shù)、航速、航深等的關(guān)系。運(yùn)用表面粒子圖像測速(簡稱表面 PIV 或SPIV)技術(shù)測量潛航體興波、以及分層流內(nèi)波的表面效應(yīng)的研究尚處于起步階段,許多問題如示蹤粒子如何選擇、如何撒布、如何從復(fù)雜背景流場中提取微弱的興波伴流信號(hào)等尚待研究。

    本文以均勻流體中潛航體興波伴流場的定量測量為例,介紹 SPIV 在定量測量表面微弱流場、表面微幅波中的應(yīng)用。

    1 實(shí)驗(yàn)設(shè)備及模型

    1.1實(shí)驗(yàn)環(huán)境及實(shí)驗(yàn)水槽

    實(shí)驗(yàn)在 50 m2左右的實(shí)驗(yàn)室中進(jìn)行,實(shí)驗(yàn)前盡量將門窗關(guān)上形成密閉空間,室溫在 14 ℃~22 ℃,室內(nèi)空氣流動(dòng)較緩慢。

    潛航體興波伴流場 SPIV 測量實(shí)驗(yàn)在一長方體型鋼化玻璃水槽中進(jìn)行,水槽置于實(shí)驗(yàn)室中間。水槽長2 400 mm、寬 800 mm、高 700 mm,水深 320~395 mm。實(shí)驗(yàn)時(shí)在水槽上方架設(shè)一面傾斜 45°的鏡子,通過鏡面反射可以從水平方向拍攝獲取垂向觀測圖像。

    1.2模型及拖曳系統(tǒng)

    拖曳模型頭部為半橢圓形,尾部為細(xì)長流線錐形。模型長度為 Lm=200 mm,最大直徑 Dm=30 mm,控制塔高度為 Th=13 mm,長度為 Tc=41 mm,控制塔距模型前端為 Td=62 mm,該模型的縮比為 1∶350。實(shí)驗(yàn)時(shí),模型距水槽底面高度 L3=295 mm,通過增減水量來調(diào)節(jié)模型距水面的距離,若無特殊說明,后文中模型距水面(或底面)距離均指拖曳線距水面(或底面)距離。

    模型采用微型減速電機(jī)驅(qū)動(dòng),額定電壓為12 V。控制電壓為 8~18 V 直流可調(diào)電源,控制精度為 0.01 V,供電壓和拖曳速度成線性關(guān)系。電壓與拖曳速度的方程如下:

    其中 U 為拖曳速度,mm/s;V 為控制電壓,V。

    為增加拖曳過程中模型的穩(wěn)定性,在模型控制塔上設(shè)置一根帶圓孔的引導(dǎo)桿,并且在拖曳線上方在按設(shè)一根與之平行且穿過導(dǎo)引桿圓孔的導(dǎo)引線,導(dǎo)引線一端固定在水槽壁上,另一端懸掛裝 1 kg 重物使之繃直以穩(wěn)定模型。拖曳方式分采用有控制塔拖曳方式(T)和無控制塔拖曳方式(NT)2 種。

    1.3表面 PIV 測量系統(tǒng)

    拍攝系統(tǒng)為一架水平放置的 CCD 相機(jī),利用45°反射鏡鏡像對(duì)水面進(jìn)行拍攝。相機(jī)鏡頭距鏡面反射點(diǎn) 1 970 mm,鏡面反射點(diǎn)距水面 770 mm。實(shí)驗(yàn)采用的CCD 相機(jī)是 Canon PowerShot SD400 IS,分辨率為1 920×1 080,攝像速度為 24 幀/s。將拍攝范圍對(duì)焦在水表面,采用連續(xù)攝像的方式來記錄水體表面粒子流動(dòng)信息,最后截取視頻圖像進(jìn)行 PIV 計(jì)算。示蹤粒子為聚酰胺(俗稱尼龍)白色粉末,粒徑 150 sm 左右。示蹤粒子由于水面張力的作用,在水面成懸浮狀態(tài)。觀測表明:撒布粒子一段時(shí)間之后,逐漸出現(xiàn)團(tuán)聚現(xiàn)象,粒子團(tuán)的直徑大多在 0.5~1 mm 之間。光源系統(tǒng)采用 2 臺(tái)幻燈投影儀,分別從左右 2 個(gè)方向以寬度約100 mm 的片光方式照亮水面上的示蹤粒子。在保證拍攝視場內(nèi)光照度均勻的前提下,盡可能減小背景雜散光的影響。

    2 表面 PIV 測量技術(shù)

    2.1測量原理

    懸浮在水面的示蹤粒子會(huì)隨表面興波伴微弱流場而產(chǎn)生運(yùn)動(dòng),SPIV 測量技術(shù)通過圖像處理方法,測量表面示蹤粒子的運(yùn)動(dòng)速度來間接得到水面興波流場分布。

    SPIV 技術(shù)和傳統(tǒng) PIV 技術(shù)共同點(diǎn)是都采用粒子圖像互相關(guān)算法計(jì)算速度場,計(jì)算軟件和處理程序具有通用性。本文采用的程序是 J.Kristian Sveen 所編寫的MatPIV v1.61,在 Matlab 平臺(tái)上運(yùn)行。不同點(diǎn)是 SPIV觀測的是飄浮在水面上的粒子運(yùn)動(dòng),傳統(tǒng) PIV 觀測的是懸浮在水中粒子的運(yùn)動(dòng)。分層流 PIV 測量的一個(gè)難點(diǎn)是難以篩選密度和待測剖面密度相等的粒子,導(dǎo)致粒子難以穩(wěn)定懸浮在理想的待測位置,其優(yōu)勢是通過調(diào)節(jié)片光源的位置,能測量流體中任意截面的速度場。SPIV 的優(yōu)點(diǎn)在于粒子受表面張力作用能穩(wěn)定漂浮在水面,測量表面速度場的精度較高,缺點(diǎn)是受表面張力、有機(jī)油污等的影響,粒子的運(yùn)動(dòng)跟隨性變差。此外,本實(shí)驗(yàn)采用的 SPIV 技術(shù)使用的是普通 CCD 相機(jī),而非價(jià)格昂貴的專業(yè) PIV 相機(jī),所使用的光源也是普通白熾燈光源,而非大功率激光光源,實(shí)驗(yàn)成本較低,實(shí)驗(yàn)設(shè)計(jì)更加靈活,簡單易行。

    2.2背景流場及消除方法

    為保證各組實(shí)驗(yàn)的獨(dú)立性,每組實(shí)驗(yàn)間隔時(shí)間不少于 30 min,然而水槽中的流體不可能處于絕對(duì)平靜,產(chǎn)生干擾有以下幾個(gè)原因:

    1)水體上下層熱對(duì)流。水面溫度隨室溫變化較大,水體內(nèi)部溫度變化較小,從而產(chǎn)生冷熱部分的相互對(duì)流,且屬于整個(gè)水槽的大范圍對(duì)流,特別是在新鋪水后熱對(duì)流較為明顯。水體穩(wěn)定之后通常水面與室溫相差 5 ℃ 以內(nèi),水面與水體內(nèi)波溫差 2℃以內(nèi);

    2)風(fēng)壓擾動(dòng)。風(fēng)壓擾動(dòng)為空氣流動(dòng)引起水面流動(dòng),包含室內(nèi)冷熱空氣對(duì)流和人為擾動(dòng);

    3)模型啟動(dòng)振蕩波干擾。由于拖曳模型從靜止?fàn)顟B(tài)突然加速,勢必會(huì)產(chǎn)生擾動(dòng)波;

    4)由于實(shí)驗(yàn)室設(shè)在 5 樓,高樓的晃動(dòng)和外界低頻振動(dòng)對(duì)水槽的影響。

    在實(shí)驗(yàn)時(shí),各種干擾因素通?;旌显谝黄穑芴綔y的背景流場強(qiáng)度可達(dá) 0.01~0.3 mm/s,運(yùn)動(dòng)潛航體興波伴流場強(qiáng)度與之同等數(shù)量級(jí),。

    假設(shè)干擾因素引起的背景流場與興波伴流場具有線性疊加性,因此可以嘗試直接將測量流場與背景流場相減。反過來說,如果減背景流場能夠有效得出實(shí)際流場的特征,則也可以說明各流場之間的關(guān)系具有線性疊加性。實(shí)驗(yàn)表明背景流場形態(tài)主要有大尺度漩渦、低頻振蕩流和表面紋波波動(dòng)。若背景由多種流場疊加,單純減掉之前某一時(shí)刻的背景不合理,甚至?xí)饘?shí)驗(yàn)后期處理的人為誤差,因此通過對(duì)減后流場是否符合特征規(guī)律來判斷背景選取是否合適。圖1是減背景前后的對(duì)比。其背景場主要為風(fēng)壓擾動(dòng)的旋渦,背景場與測量場間隔 2 s,因此背景平移量也很小,沒有明顯的波動(dòng)現(xiàn)象,且 2 幅背景相減后速度場均方根降低一個(gè)數(shù)量級(jí),所以這類情況采用減單幅背景場的方法十分有效。后文展示計(jì)算結(jié)果均為減背景后的速度場,不再加以說明。

    2.3測量精度及標(biāo)定誤差評(píng)估

    2.3.1測量范圍和精度評(píng)估

    SPIV 將圖像分割成小塊像素查詢區(qū)域,對(duì)每一塊查詢區(qū)域?qū)ふ亦徲蛑谢ハ嚓P(guān)率最高的查詢區(qū),由此確定移動(dòng)距離。計(jì)算結(jié)果表明:通過圖像處理能夠獲得的互相關(guān)平移距離精度為 0.01 pixel。若取標(biāo)定系數(shù)為0.2 mm/pixel,因此有效精度為 0.002 mm??赏扑銣y量范圍為 380 mm×220 mm,標(biāo)定系數(shù)隨具體實(shí)驗(yàn)條件來確定。

    圖1 拖曳深度 46 mm,拖曳速度為120 mm/s,拖曳方式NT,速度場空間分辨率為 6.4 mm×6.4 mm,標(biāo)定系數(shù)為 0.19 mm/pixelFig.1 Depth=46 mm,towed speed=120 mm/s,towed model is NT,(1-a)background velocity field,(1-b)substraction between two background fields,(1-c)measured velocity field,(1-d)actual velocity field after background substracted,velocity spatial resolution is 6.4 mm×6.4 mm,comap is 0.19 mm/pixel

    速度場空間分辨率為一速度矢量代表實(shí)際面積,與查詢區(qū)大小、重疊率以及標(biāo)定系數(shù)密切相關(guān)。表1列舉標(biāo)定系數(shù)為 0.2 mm/pixel 的幾種查詢區(qū)和重疊率對(duì)應(yīng)的速度場空間分辨率。實(shí)驗(yàn)時(shí)能清晰明辨的特征粒子有限,采用較小查詢區(qū)容易導(dǎo)致計(jì)算誤差,因此采用大查詢區(qū)和高重疊率來保證計(jì)算準(zhǔn)確性和提高分辨率。

    表1 標(biāo)定系數(shù) 0.2 mm/pixel 的幾種查詢區(qū)和重疊率對(duì)應(yīng)速度場空間分辨率Tab.1 Velocity spatial resolution under several interrogation regions and overlap with comap=0.2 mm/pixel

    2.3.2標(biāo)定誤差評(píng)估

    嚴(yán)格來說,模型與水面有一定距離,如圖2所示,水對(duì)光線的折射和相機(jī)成像視角會(huì)造成標(biāo)定誤差。圖中 Lm為實(shí)際運(yùn)動(dòng)距離,Lm''為觀測運(yùn)動(dòng)距離,光線從水中進(jìn)入到空氣中的入射角為β,折射角為α,且β<α,設(shè)標(biāo)定系數(shù)誤差為 δ,經(jīng)測量,鏡頭距水面l=2 740 mm,模型深度 h=25~100 mm,水的折射率約為 n=1.3。假設(shè)光線從水進(jìn)入空氣不發(fā)生折射,用一個(gè)不等式來說明標(biāo)定誤差 δ 的范圍:

    計(jì)算誤差最大為 3.6%,所以在本實(shí)驗(yàn)中予以忽略,更多的是要考慮溫度、風(fēng)壓等環(huán)境因素帶來的影響。

    3 測量結(jié)果分析

    圖2 標(biāo)定誤差示意圖Fig.2 Error of calibration

    本實(shí)驗(yàn)采用 SPIV 技術(shù)成功的測量了潛航體激發(fā)的微弱興波伴流場。本文使用基于潛深 h 的傅氏數(shù)和雷諾數(shù)作為無綱量相相似參數(shù),其中 20 ℃淡水中 v 為 10-2cm2/s,U 為航速,g為重力加速度。因此傅氏數(shù)范圍在 0.086<Fr<0.390,2 550<Re<5 790。在此實(shí)驗(yàn)條件下,潛艇所產(chǎn)生的水面興波波幅十分微弱,肉眼幾乎看不見水面的隆起或變形。但由潛艇排水效應(yīng)所產(chǎn)生的水面流場變化經(jīng)過SPIV 計(jì)算卻是十分明顯的速度場特征。

    圖3 所示為 2 種拖曳方式(T 和 NT),深度 50 mm條件下的速度場矢量圖。由于拍攝范圍有限,quiver畫矢量的長度具有相對(duì)性,在不同的 2 幅圖中,矢量長度并不代表值的大小,同一幅圖中矢量長度越長值越大。相機(jī)由手動(dòng)調(diào)節(jié)縮放,標(biāo)定系數(shù)略有差異。為了顯示清晰采用較低速度場分辨率。

    圖3 深度 50 mm,2 種拖曳模式測量結(jié)果展示Fig.3 Depth is 50mm,measured result of two towed models

    由于潛艇控制塔高 17 mm,因此控制塔距水面33 mm,在此深度下,控制塔對(duì)水面流場的影響十分明顯。從圖3(a)可看出,在控制塔上方速度場明顯加劇。圖中用黑色橢圓大致標(biāo)示出模型所在位置,模型的頭部表現(xiàn)為‘源',而尾部表現(xiàn)為‘匯'的特征,模型的周圍形成 V 型流場結(jié)構(gòu)。較于 NT 模式,控制塔和模型主體共同產(chǎn)生體積效應(yīng),導(dǎo)致控制塔上方流場強(qiáng)度更加強(qiáng)烈。

    下面以拖曳方式 T(深度 45 mm,拖曳速度120 mm/s,有控制塔)為例,給出速度場的偽彩色分布圖,規(guī)定從左往右為 x 軸速度場正方向,從上往下為 y 軸速度場正方向,如圖4所示。

    圖4 h=45 mm,v=120 mm/s,F(xiàn)r=0.179,Re=3600,拖曳方式為 T,速度場空間分辨率為 3.2 mm×3.2 mm。Fig.4 h=45 mm,v=120 mm/s,F(xiàn)r=0.179,Re=3 600,towed model is T,velocity spatial resolution is 3.2 mm*3.2 mm.4-a is vector of velocity field,4-b is pcolor image of united velocity field,4-c is pcolor image of velocity_u field if xaxis direction,4-d is pcolor image of velocity_v field if yaxis direction

    由圖4 可知:1)由于運(yùn)動(dòng)模型體積排水效應(yīng),頭尾處成一對(duì)源匯流場,且上下分布對(duì)稱;2)速度場 v偽彩圖為蝴蝶狀,模型后方流場成 V 型,是 Kelvin波;3)速度場 u 偽彩圖中矩形框內(nèi)疑似橫波波系,肉眼可見 1~2 個(gè)橫波,速度場 v 偽彩圖中矩形框內(nèi)疑似側(cè)波波系,肉眼可以見兩列側(cè)波。

    提取圖5模型軸線 x 方向和尾跡(圖4(d)所示方框內(nèi))y 方向的速度,如圖5所示。圖5(a)中箭頭所指的位置推測存在 4 個(gè)波峰,對(duì)應(yīng) 4 個(gè)橫波波系,圖5(b)中箭頭所指的位置出推測存在 2 個(gè)波峰,對(duì)應(yīng)一個(gè)散波波系,要采用濾波手段去掉強(qiáng)源匯效應(yīng),才能進(jìn)一步分析。

    陳祥瑞等[11]利用線陣 CCD 拍攝拖曳細(xì)長回轉(zhuǎn)體的水面興波,在 Fr=0.28,Re=14 000 條件下,線陣 CCD幾乎探測不到表面興波信息,而 SPIV 技術(shù)在 Fr=0.18,Re=3 600 條件下能清晰獲取水面興波伴流場(見圖4)。

    圖5 直線上速度分量變化曲線Fig.5 Change of component of u,v along line in velocity field

    4 結(jié) 語

    本實(shí)驗(yàn)主要得出以下結(jié)論:

    1)SPIV 技術(shù)使用普通 CCD 相機(jī)、普通白熾燈光源以及聚合物粒子粉末,在低弗洛德數(shù)和低雷諾數(shù)條件下,能有效的測量出模型運(yùn)動(dòng)興波伴流場的表面流場,可以取代昂貴的專業(yè) PIV 測量系統(tǒng),也較線陣 CCD的測量能靈敏的興波信息。采用 64×64 或 32×32 查詢區(qū)實(shí)際計(jì)算中發(fā)現(xiàn) SPIV 計(jì)算精度能達(dá)到 0.01 pixel。速度場空間分辨率與標(biāo)定系數(shù)、查詢區(qū)大小和重疊率密切相關(guān),通常采用標(biāo)定系數(shù)為 0.2 mm/pixel,由此可計(jì)算速度場空間分辨率為 3.2 mm×3.2 mm,實(shí)際位移計(jì)算精度可達(dá)到 0.002 mm,在此速度場空間分辨率下展示出的流場細(xì)節(jié)構(gòu)能夠滿足分析流場特征的需要;

    2)背景干擾流場的來源有大范圍熱對(duì)流、風(fēng)壓擾動(dòng)、啟動(dòng)擾動(dòng)波,主要擾動(dòng)來源為風(fēng)壓擾動(dòng)。在水面的表現(xiàn)形式為大尺度漩渦、低頻振蕩流、表面紋波波動(dòng)。假設(shè)干擾流場具有線性疊加性,根據(jù)實(shí)際情況采用減單幅背景的方法能有效減去旋渦背景;

    3)列舉了 T 和 NT 兩種模式,F(xiàn)r=0.154 和 0.279,Re=3 240 和 5 850 條件下的表面興波速度場計(jì)算結(jié)果,定量觀測到了表面 V 型波速度場,以及潛航體運(yùn)動(dòng)產(chǎn)生的‘源匯'效應(yīng)。對(duì)速度場分量 u 和 v 偽彩圖中可以看到速度場分量 v 為蝴蝶狀,可探測 2 列側(cè)波,速度場分量 u 中可探測到至少 1 列橫波。

    今后還將在具體分析水面速度場特征隨航速、深度的演化規(guī)律以及水面張力的變化對(duì)速度場的影響規(guī)律。

    [1]張軍,張效慈,趙峰,等.源于水動(dòng)力學(xué)的潛艇尾跡非聲探測技術(shù)研究之進(jìn)展[J].船舶力學(xué),2003,7(2):121-128.

    [2]魏崗,戴世強(qiáng).分層流體中運(yùn)動(dòng)源生成的內(nèi)波研究進(jìn)展[J].力學(xué)進(jìn)展,2006,36(1): 111-124.

    [3]CHEN Y Q,F(xiàn)ENG J,ZHU M H.Detection methods of submerged mobile using SAR images[C]//Proceedings of 2005 IEEE International Geoscience and Remote Sensing Symposium.Seoul: IEEE,2005,3: 1717-1720.

    [4]馬暉揚(yáng),麻柏坤,張人杰.分層流體中運(yùn)動(dòng)物體與自由面相互作用的實(shí)驗(yàn)研究[J].流體力學(xué)實(shí)驗(yàn)與測量,2000,14(3): 7-11,17.

    [5]洪方文,常煜.均勻流中潛艇水下運(yùn)動(dòng)表面尾跡的數(shù)值模擬[J].船舶力學(xué),2005,9(4): 9-17.

    [6]張效慈.潛艇內(nèi)波尾跡物理場在海面映波量值的確定[J].船舶力學(xué),2005,9(4): 25-30.

    [7]VOROPAYEV S I,F(xiàn)ERNANDO H J S,SMIRNOV S A,et al.On surface signatures generated by submerged momentum sources[J].Physics of Fluids,2007,19(7): 076603.

    [8]秦朝峰,陳旭,于佳.躍層中運(yùn)動(dòng)小球產(chǎn)生的內(nèi)波對(duì)表層流場的影響[J].中國海洋大學(xué)學(xué)報(bào),2009,39(S): 303-310.

    [9]ROTTMAN J W,BROUTMN D,SPEDDING G R,et al.Internal wave generation by a horizontally moving sphere at low Froude number[C]//Proceedings of the 25th Symposium on Naval Hydrodynamics,St.John's.Canada: Newfoundland and Labrador,2004.

    [10]姚志崇,趙峰,梁川,等.分層流體中拖曳球體尾流及輻射內(nèi)波試驗(yàn)研究[J].船舶力學(xué),2014,18(11):1275-1283.

    [11]陳祥瑞.基于分層流水面CCD影像的運(yùn)動(dòng)潛體尾跡結(jié)構(gòu)特征實(shí)驗(yàn)及其應(yīng)用研究[D].南京: 南京理工大學(xué),2014.

    [12]馬暉揚(yáng),麻柏坤,張人杰.分層流體中物體運(yùn)動(dòng)尾跡的理論和實(shí)驗(yàn)研究[J].中國科學(xué)技術(shù)大學(xué)學(xué)報(bào),2000,30(6): 677-684.

    [13]魏崗,樂嘉春,戴世強(qiáng).有限深度兩層流體系統(tǒng)中運(yùn)動(dòng)點(diǎn)源生成的內(nèi)波及其與自由面的相互作用[J].應(yīng)用數(shù)學(xué)和力學(xué),2003,24(9): 906-918.

    [14]HUDIMAC A A.Ship waves in a stratified ocean[J].Journal of Fluid Mechanics,1961,11(2): 229-243.

    [15]GASPAROVIC R F,APEL J R,KASISCHKE E S.An overview of the SAR internal wave signature experiment[J].Journal of Geophysical Research: Oceans (1978-2012),1988,93(C10):12304-12316.

    The application of surface PIV in the measurement of wave making of underwater moving body

    ZHOU Wen-jin,JIANG Xiao-qin,WANG Jian-zhong,F(xiàn)ANG Pin-jie
    (Department of Physics,Naval University of Engineering,Wuhan 430033,China)

    This paper is mainly to introduce an experimental method for measuring wave making of underwater moving body by PIV.Apart from of analysing the principle of Surface-PIV,defining coordinate,and how to eliminate interference from background,the precision of Surface-PIV is also analysed in detail.Experiment show that 'V' type wave and sourcesink effect are successfully detected.Comparing with the traditional PIV,Surface-PIV is not only economical and practicable under the laboratory conditions,but also accurate and reliable.Effective accuracy of Surface-PIV can reach 0.01pixel.

    underwater moving body;wave making;surface PIV

    P714.3;O353.2

    A

    1672-7619(2016)06-0075-06

    10.3404/j.issn.1672-7619.2016.06.015

    2015-09-15;

    2015-10-13

    基金支助:海軍工程大學(xué)自然科學(xué)基金資助項(xiàng)目(HJGSK2014G122)

    周文進(jìn)(1990-),男,碩士研究生,從事分層海水內(nèi)波現(xiàn)象研究。

    猜你喜歡
    興波內(nèi)波水面
    孤立內(nèi)波對(duì)過渡海域聲場干涉結(jié)構(gòu)的影響分析
    船舶興波運(yùn)動(dòng)的非線性現(xiàn)象算法研究與仿真
    基于Friendship的高速客船特征參數(shù)靈敏度分析
    內(nèi)波與死水,連潛艇都怕的海浪
    基于MODIS 遙感影像的安達(dá)曼海內(nèi)波特征參數(shù)分布及生成周期研究
    水黽是怎樣浮在水面的
    創(chuàng)造足以亂真的水面反光
    爭奪水面光伏
    能源(2016年3期)2016-12-01 05:11:02
    一塊水面
    滇池(2014年5期)2014-05-29 07:30:25
    SAR極化特征圖像與σ0圖像的海洋內(nèi)波可視性對(duì)比
    99国产极品粉嫩在线观看| 男人舔女人下体高潮全视频| 97碰自拍视频| 一级a爱片免费观看的视频| 禁无遮挡网站| 给我免费播放毛片高清在线观看| 香蕉国产在线看| 日韩一卡2卡3卡4卡2021年| 久久久精品国产亚洲av高清涩受| 精品欧美国产一区二区三| 午夜精品久久久久久毛片777| 非洲黑人性xxxx精品又粗又长| 成人国语在线视频| 可以在线观看的亚洲视频| 久久久国产成人免费| 人妻久久中文字幕网| 成人18禁在线播放| 制服人妻中文乱码| 看片在线看免费视频| 欧美日韩瑟瑟在线播放| 超碰成人久久| 国产亚洲精品久久久久久毛片| 动漫黄色视频在线观看| 男人舔奶头视频| 他把我摸到了高潮在线观看| 久久久水蜜桃国产精品网| 午夜亚洲福利在线播放| 国产黄色小视频在线观看| 日本免费a在线| 性色av乱码一区二区三区2| 精品免费久久久久久久清纯| 欧美在线一区亚洲| 1024视频免费在线观看| 国产精品久久久久久人妻精品电影| 一级片免费观看大全| 黄片大片在线免费观看| 高清毛片免费观看视频网站| 别揉我奶头~嗯~啊~动态视频| 女警被强在线播放| e午夜精品久久久久久久| 成人欧美大片| 自线自在国产av| 长腿黑丝高跟| 日本一区二区免费在线视频| 久久精品国产综合久久久| 91字幕亚洲| 日韩三级视频一区二区三区| 中文字幕精品亚洲无线码一区 | 国产精品1区2区在线观看.| 国产精品久久久av美女十八| 免费观看人在逋| 老司机午夜十八禁免费视频| 一本综合久久免费| 女人高潮潮喷娇喘18禁视频| 国产精品久久久av美女十八| 18禁国产床啪视频网站| 88av欧美| 视频区欧美日本亚洲| 日本免费a在线| 18美女黄网站色大片免费观看| 国产精华一区二区三区| 亚洲成av人片免费观看| 国产真实乱freesex| 久久久久久亚洲精品国产蜜桃av| 亚洲av电影不卡..在线观看| 天堂影院成人在线观看| 成人一区二区视频在线观看| 亚洲一区二区三区不卡视频| 777久久人妻少妇嫩草av网站| 国产一级毛片七仙女欲春2 | 午夜免费激情av| 黄色片一级片一级黄色片| 人妻久久中文字幕网| 色尼玛亚洲综合影院| 1024视频免费在线观看| 亚洲成av人片免费观看| 欧美色欧美亚洲另类二区| 午夜久久久在线观看| 成人18禁在线播放| 真人做人爱边吃奶动态| 两性午夜刺激爽爽歪歪视频在线观看 | 母亲3免费完整高清在线观看| 免费看日本二区| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜a级毛片| 日韩av在线大香蕉| 熟女少妇亚洲综合色aaa.| 亚洲午夜精品一区,二区,三区| 国产精品香港三级国产av潘金莲| 中亚洲国语对白在线视频| 国产单亲对白刺激| 成人18禁高潮啪啪吃奶动态图| 午夜日韩欧美国产| 中文字幕av电影在线播放| 在线观看舔阴道视频| 欧美性长视频在线观看| 女同久久另类99精品国产91| 亚洲国产欧美日韩在线播放| 国产男靠女视频免费网站| 国产高清激情床上av| 亚洲av成人av| 国产国语露脸激情在线看| 久久 成人 亚洲| 亚洲电影在线观看av| 一进一出好大好爽视频| 日韩大码丰满熟妇| 中文字幕久久专区| 婷婷精品国产亚洲av在线| 亚洲狠狠婷婷综合久久图片| 国产视频内射| 亚洲一区中文字幕在线| 欧美黑人精品巨大| 后天国语完整版免费观看| 脱女人内裤的视频| 午夜影院日韩av| 免费观看人在逋| 欧美成人免费av一区二区三区| 手机成人av网站| 久久婷婷人人爽人人干人人爱| 色精品久久人妻99蜜桃| 啦啦啦观看免费观看视频高清| 老熟妇仑乱视频hdxx| 久久久久久大精品| 一卡2卡三卡四卡精品乱码亚洲| 婷婷亚洲欧美| 88av欧美| 精品无人区乱码1区二区| 女人爽到高潮嗷嗷叫在线视频| 2021天堂中文幕一二区在线观 | 亚洲va日本ⅴa欧美va伊人久久| 国产成人系列免费观看| 一本久久中文字幕| 精品久久久久久久人妻蜜臀av| 欧美+亚洲+日韩+国产| 嫩草影院精品99| 成人国语在线视频| 亚洲人成网站在线播放欧美日韩| 亚洲熟妇熟女久久| 久久久国产成人精品二区| 欧美日韩瑟瑟在线播放| 亚洲熟妇中文字幕五十中出| e午夜精品久久久久久久| 亚洲最大成人中文| 亚洲 国产 在线| av在线天堂中文字幕| 老司机午夜十八禁免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成人国产一区在线观看| 18禁黄网站禁片午夜丰满| 国产精品亚洲美女久久久| 国产精品香港三级国产av潘金莲| 青草久久国产| 国产亚洲精品av在线| 97人妻精品一区二区三区麻豆 | 国产精品亚洲美女久久久| 国产真人三级小视频在线观看| 久久久久国内视频| 热re99久久国产66热| 国产野战对白在线观看| 久久久久国内视频| 欧美日本视频| 一进一出好大好爽视频| 夜夜夜夜夜久久久久| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 亚洲中文日韩欧美视频| 啦啦啦免费观看视频1| 丝袜人妻中文字幕| 国产亚洲精品一区二区www| av在线播放免费不卡| 久久久久九九精品影院| 一区二区三区国产精品乱码| 欧美黑人精品巨大| 午夜激情福利司机影院| 欧美精品啪啪一区二区三区| 亚洲精品美女久久av网站| 欧美色欧美亚洲另类二区| 亚洲国产精品合色在线| 女生性感内裤真人,穿戴方法视频| 欧美黄色淫秽网站| 精品一区二区三区视频在线观看免费| 一夜夜www| 欧美大码av| 成人国产一区最新在线观看| a级毛片在线看网站| 精品国产乱子伦一区二区三区| 日本五十路高清| 欧美乱妇无乱码| 久久久久久久久中文| 丁香欧美五月| 精品熟女少妇八av免费久了| 777久久人妻少妇嫩草av网站| 99在线人妻在线中文字幕| 日本撒尿小便嘘嘘汇集6| 婷婷亚洲欧美| 亚洲人成77777在线视频| 色综合欧美亚洲国产小说| 熟女电影av网| 淫秽高清视频在线观看| 久久午夜综合久久蜜桃| 亚洲国产毛片av蜜桃av| 亚洲成av人片免费观看| 老熟妇乱子伦视频在线观看| 成人亚洲精品一区在线观看| 一本精品99久久精品77| 亚洲第一av免费看| 女警被强在线播放| 黄网站色视频无遮挡免费观看| 亚洲狠狠婷婷综合久久图片| 黄色a级毛片大全视频| 亚洲精品在线观看二区| 久久午夜亚洲精品久久| 又黄又爽又免费观看的视频| 黄色视频不卡| 美女扒开内裤让男人捅视频| 久久香蕉国产精品| 欧美大码av| 国产三级在线视频| 久久精品成人免费网站| 熟妇人妻久久中文字幕3abv| 天堂√8在线中文| 欧美丝袜亚洲另类 | 无限看片的www在线观看| 免费av毛片视频| 三级毛片av免费| 婷婷精品国产亚洲av| 琪琪午夜伦伦电影理论片6080| 亚洲avbb在线观看| 日韩欧美国产一区二区入口| 青草久久国产| 国产精品一区二区精品视频观看| 叶爱在线成人免费视频播放| 怎么达到女性高潮| 国产蜜桃级精品一区二区三区| 国产成人av教育| 欧美精品亚洲一区二区| 黄色视频,在线免费观看| 一本精品99久久精品77| 久久狼人影院| 在线看三级毛片| 国产亚洲精品久久久久久毛片| 极品教师在线免费播放| 免费看日本二区| 少妇熟女aⅴ在线视频| www日本黄色视频网| 18禁美女被吸乳视频| www.精华液| 成人18禁高潮啪啪吃奶动态图| 日韩大尺度精品在线看网址| 亚洲成人免费电影在线观看| 欧美午夜高清在线| 在线观看午夜福利视频| 制服人妻中文乱码| 亚洲 欧美 日韩 在线 免费| 1024香蕉在线观看| 桃红色精品国产亚洲av| 国产人伦9x9x在线观看| 一级a爱视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆 | 女生性感内裤真人,穿戴方法视频| 亚洲五月婷婷丁香| 免费在线观看日本一区| 婷婷六月久久综合丁香| 亚洲成国产人片在线观看| 精品一区二区三区视频在线观看免费| 国产成人精品无人区| 欧美最黄视频在线播放免费| 欧美激情 高清一区二区三区| 久久草成人影院| 在线观看www视频免费| 亚洲一码二码三码区别大吗| 91九色精品人成在线观看| 波多野结衣高清作品| 精品日产1卡2卡| 日韩免费av在线播放| 在线观看午夜福利视频| 午夜激情av网站| 精品久久久久久久毛片微露脸| 天堂影院成人在线观看| 黄色女人牲交| 亚洲成国产人片在线观看| 日韩免费av在线播放| 亚洲精品久久国产高清桃花| av电影中文网址| 亚洲熟妇中文字幕五十中出| 中文在线观看免费www的网站 | 久9热在线精品视频| 美女国产高潮福利片在线看| 欧美性猛交╳xxx乱大交人| 黄色a级毛片大全视频| 亚洲精品久久成人aⅴ小说| 美女 人体艺术 gogo| 国产亚洲精品久久久久久毛片| 狠狠狠狠99中文字幕| 在线国产一区二区在线| 变态另类丝袜制服| 看免费av毛片| 美女高潮喷水抽搐中文字幕| 国产人伦9x9x在线观看| 久久99热这里只有精品18| 国产区一区二久久| 自线自在国产av| 香蕉国产在线看| 老熟妇乱子伦视频在线观看| 满18在线观看网站| 又黄又粗又硬又大视频| 久久久精品欧美日韩精品| 夜夜躁狠狠躁天天躁| 国产av一区二区精品久久| 最新在线观看一区二区三区| av中文乱码字幕在线| 国产熟女午夜一区二区三区| 亚洲精品中文字幕在线视频| 久久九九热精品免费| 国产精品1区2区在线观看.| 搡老妇女老女人老熟妇| 中文字幕精品免费在线观看视频| 国产精品日韩av在线免费观看| 亚洲 欧美 日韩 在线 免费| 久久久久国内视频| 在线观看免费视频日本深夜| 欧美黑人巨大hd| 亚洲自偷自拍图片 自拍| 国产精品日韩av在线免费观看| 夜夜躁狠狠躁天天躁| 欧美午夜高清在线| 亚洲欧洲精品一区二区精品久久久| 啦啦啦免费观看视频1| 欧美色欧美亚洲另类二区| 久久久久久人人人人人| 可以在线观看的亚洲视频| 久久久久久免费高清国产稀缺| 久久天躁狠狠躁夜夜2o2o| 老司机深夜福利视频在线观看| 国产一区二区三区在线臀色熟女| 一级毛片高清免费大全| 天堂√8在线中文| 亚洲人成网站高清观看| 美女高潮喷水抽搐中文字幕| 国产精品一区二区免费欧美| 宅男免费午夜| 久久国产精品人妻蜜桃| 亚洲,欧美精品.| 午夜成年电影在线免费观看| 啦啦啦免费观看视频1| 欧美黄色淫秽网站| 丁香欧美五月| 午夜福利成人在线免费观看| 久久人人精品亚洲av| www.www免费av| netflix在线观看网站| 亚洲精品久久成人aⅴ小说| 久久精品成人免费网站| 色播在线永久视频| 日本 av在线| 夜夜爽天天搞| 精品高清国产在线一区| 久久精品国产亚洲av香蕉五月| 99精品久久久久人妻精品| 免费无遮挡裸体视频| 国产黄a三级三级三级人| 久久精品成人免费网站| 国产又黄又爽又无遮挡在线| 最好的美女福利视频网| 大型av网站在线播放| 国产精品永久免费网站| 美女国产高潮福利片在线看| 免费av毛片视频| 国产av一区在线观看免费| 少妇熟女aⅴ在线视频| 久久中文看片网| x7x7x7水蜜桃| 国产成+人综合+亚洲专区| 国产又色又爽无遮挡免费看| 午夜福利视频1000在线观看| 99久久99久久久精品蜜桃| 老司机靠b影院| 国产一卡二卡三卡精品| 1024手机看黄色片| 少妇熟女aⅴ在线视频| 人人妻人人澡欧美一区二区| 日本免费a在线| 两个人看的免费小视频| 午夜激情av网站| 黑人欧美特级aaaaaa片| 妹子高潮喷水视频| 国产精品九九99| www.精华液| 国产高清视频在线播放一区| 国产精品精品国产色婷婷| 亚洲欧美精品综合久久99| 亚洲五月色婷婷综合| 黄色视频不卡| 欧美黄色片欧美黄色片| 2021天堂中文幕一二区在线观 | 亚洲熟女毛片儿| 99国产综合亚洲精品| 欧美久久黑人一区二区| 国产精品综合久久久久久久免费| 免费搜索国产男女视频| 欧美乱色亚洲激情| 亚洲熟妇中文字幕五十中出| 成熟少妇高潮喷水视频| 精品人妻1区二区| 在线观看日韩欧美| 一边摸一边做爽爽视频免费| 99热只有精品国产| 又大又爽又粗| 国产1区2区3区精品| 国产三级在线视频| 欧美乱妇无乱码| 久久热在线av| 亚洲欧美精品综合一区二区三区| 19禁男女啪啪无遮挡网站| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产欧美日韩在线播放| 69av精品久久久久久| 成人一区二区视频在线观看| 国产av一区二区精品久久| 操出白浆在线播放| 一级毛片女人18水好多| 欧美在线黄色| 免费在线观看成人毛片| 99re在线观看精品视频| 在线播放国产精品三级| 亚洲成av片中文字幕在线观看| 亚洲 欧美一区二区三区| 亚洲精品国产精品久久久不卡| 久久国产亚洲av麻豆专区| 久久久精品欧美日韩精品| 亚洲自偷自拍图片 自拍| 免费人成视频x8x8入口观看| 日本一本二区三区精品| 国产一区二区在线av高清观看| 精品久久久久久久人妻蜜臀av| av在线天堂中文字幕| www日本在线高清视频| 老鸭窝网址在线观看| 又紧又爽又黄一区二区| tocl精华| 国产精品久久久久久亚洲av鲁大| 欧美丝袜亚洲另类 | 国产精品免费视频内射| 禁无遮挡网站| 男人舔奶头视频| 国产单亲对白刺激| 国产欧美日韩一区二区三| 亚洲欧美日韩高清在线视频| 成在线人永久免费视频| 曰老女人黄片| 国产色视频综合| 午夜视频精品福利| 国产精品亚洲av一区麻豆| 国产不卡一卡二| 久久久久精品国产欧美久久久| 在线十欧美十亚洲十日本专区| 久久香蕉精品热| 大香蕉久久成人网| 999久久久国产精品视频| 最新在线观看一区二区三区| 男女做爰动态图高潮gif福利片| 欧美大码av| 国产午夜精品久久久久久| 亚洲国产欧美网| 男女之事视频高清在线观看| 亚洲五月色婷婷综合| 久久久国产成人免费| 国产蜜桃级精品一区二区三区| 正在播放国产对白刺激| 日韩三级视频一区二区三区| 一本综合久久免费| 视频在线观看一区二区三区| 精品电影一区二区在线| 狂野欧美激情性xxxx| 国产av一区二区精品久久| 日韩大尺度精品在线看网址| 欧美激情高清一区二区三区| 久久国产精品影院| 午夜精品久久久久久毛片777| 免费看a级黄色片| 午夜久久久在线观看| 中文字幕人妻丝袜一区二区| netflix在线观看网站| 免费人成视频x8x8入口观看| 麻豆久久精品国产亚洲av| 亚洲人成伊人成综合网2020| 亚洲片人在线观看| 午夜免费鲁丝| 欧美激情高清一区二区三区| 女同久久另类99精品国产91| 美女 人体艺术 gogo| 午夜福利欧美成人| 日韩大尺度精品在线看网址| 欧美日韩乱码在线| 黄网站色视频无遮挡免费观看| 无人区码免费观看不卡| 国内揄拍国产精品人妻在线 | 国内少妇人妻偷人精品xxx网站 | 成年人黄色毛片网站| 丝袜在线中文字幕| 欧美中文日本在线观看视频| 18禁国产床啪视频网站| 波多野结衣巨乳人妻| 十分钟在线观看高清视频www| 在线观看一区二区三区| 怎么达到女性高潮| 欧美激情高清一区二区三区| 国产视频一区二区在线看| 日韩有码中文字幕| 性欧美人与动物交配| 亚洲美女黄片视频| www国产在线视频色| 亚洲激情在线av| 看片在线看免费视频| 最近最新中文字幕大全免费视频| 麻豆成人av在线观看| 禁无遮挡网站| 999精品在线视频| 久9热在线精品视频| 成人av一区二区三区在线看| 日韩国内少妇激情av| 欧美在线黄色| 一级作爱视频免费观看| 亚洲av五月六月丁香网| 黄色视频不卡| 色播在线永久视频| 成人免费观看视频高清| 国产亚洲精品一区二区www| 一夜夜www| 欧美日韩精品网址| 国产精品久久视频播放| 一区二区三区激情视频| 欧美乱妇无乱码| 国内少妇人妻偷人精品xxx网站 | 国产精品久久久久久亚洲av鲁大| 久久久精品国产亚洲av高清涩受| 久久欧美精品欧美久久欧美| 亚洲成av人片免费观看| 国产爱豆传媒在线观看 | 亚洲欧美精品综合久久99| 午夜a级毛片| 男女那种视频在线观看| 亚洲av美国av| 国产99久久九九免费精品| 日韩三级视频一区二区三区| 国产亚洲精品综合一区在线观看 | 久久九九热精品免费| 一区二区三区高清视频在线| 国产成人av教育| 欧美日韩福利视频一区二区| 色在线成人网| 色综合欧美亚洲国产小说| 啦啦啦免费观看视频1| 在线播放国产精品三级| 长腿黑丝高跟| 黑人操中国人逼视频| 欧美性猛交黑人性爽| 天天一区二区日本电影三级| 黑人巨大精品欧美一区二区mp4| 午夜精品在线福利| 12—13女人毛片做爰片一| 在线观看免费午夜福利视频| 国产精品久久久久久亚洲av鲁大| 老司机午夜十八禁免费视频| 国产野战对白在线观看| 一卡2卡三卡四卡精品乱码亚洲| 免费av毛片视频| 精品第一国产精品| 啦啦啦 在线观看视频| 久久精品国产综合久久久| 最近最新中文字幕大全免费视频| 亚洲自偷自拍图片 自拍| 夜夜爽天天搞| 又紧又爽又黄一区二区| 18美女黄网站色大片免费观看| 精品不卡国产一区二区三区| 久久久国产成人免费| www.自偷自拍.com| 亚洲一区高清亚洲精品| 久久亚洲精品不卡| 日本一本二区三区精品| 老汉色∧v一级毛片| 亚洲成人久久爱视频| 午夜激情福利司机影院| 午夜视频精品福利| 国产成人av教育| 欧美av亚洲av综合av国产av| 精品欧美一区二区三区在线| 日韩 欧美 亚洲 中文字幕| 成年人黄色毛片网站| 国产一级毛片七仙女欲春2 | 一区二区三区激情视频| 1024手机看黄色片| 日本五十路高清| 欧美丝袜亚洲另类 | 满18在线观看网站| 一二三四在线观看免费中文在| 亚洲成av人片免费观看| 午夜福利在线在线| 国产v大片淫在线免费观看| 精品国产美女av久久久久小说| 黄片播放在线免费| 亚洲狠狠婷婷综合久久图片| 在线十欧美十亚洲十日本专区| 亚洲aⅴ乱码一区二区在线播放 | 一边摸一边做爽爽视频免费| 国产成人精品无人区| 日本成人三级电影网站| 18美女黄网站色大片免费观看| 欧美日韩黄片免| 久久精品成人免费网站| 国产高清视频在线播放一区| 国产97色在线日韩免费| 两个人免费观看高清视频| 日本一本二区三区精品|