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

    魚游動渦結構PIV實驗研究

    2020-11-03 07:40:14王福君王洪平魏潤杰劉彥鵬
    實驗流體力學 2020年5期
    關鍵詞:尾跡魚尾旋渦

    王福君, 王洪平, 高 琪, 魏潤杰, 劉彥鵬

    (1.浙江大學航空航天學院 流體工程研究所, 杭州 310027; 2.中國科學院力學研究所 非線性力學國家重點實驗室, 北京 100190; 3.北京立方天地科技有限責任公司, 北京 100083; 4.北京航空航天大學航空科學與工程學院 流體力學研究所, 北京 100191)

    0 引 言

    魚類經過幾億年進化出高效、高機動的流線型身軀,能夠靈活應對復雜的水下環(huán)境。其特有的運動方式一直是流體力學的研究熱點,對水下航行器的推進有著重要指導意義和應用價值。因此,眾多學者通過理論分析、流動顯示、粒子圖像測速(PIV)實驗以及數值模擬等方法對該領域展開了研究[1]。

    現有對魚類的研究根據不同的游動速度將魚類游動模式大致分為了BCF(Body and/or Caudal Fin)和MPF(Median and/or Paired Fin) 2種[2]。前者依靠身體和尾鰭的擺動來提供動力,使得魚快速游動;后者依靠胸鰭、背鰭和腹鰭實現原地靜止游動。但是對單一的個體而言,魚的游動會同時表現出2種游動模式,以實現快速、高效的機動。根據Sfakiotakis等[2]的結論,BCF模式和MPF模式分別可細分出多種不同的模式,分別展示在圖1(a)和(b)中。BCF模式包含從波動(Undulatory)到擺動(Oscillatory)大致5種游動模式[2];MPF模式可分為兩大類,分別是鰭波動模式、鰭擺動模式[2]。中國船舶研究中心的張軍等[3]采用PIV技術對MPF模式下長鰭波動推進的流場流向結構進行實驗研究,給出了長鰭在不同截面的相位平均速度場。

    圖1 游動模式的演變圖[2]

    Smits等[1]對魚的波動和擺動進行了詳細的綜述。波動運動以長條體型魚的游動為主,典型代表為鰻鱺;擺動運動又可分為縱搖運動(Pitch motion)和垂蕩運動(Heave motion),以尾鰭發(fā)達且成月牙型的魚類為主,典型代表為金槍魚。事實上,自然環(huán)境中魚類的游動大多介于波動游動和擺動游動之間,而對于各類游動方式的研究大多僅限于數值模擬,實驗上大多也只采用形態(tài)學觀測或者平面PIV測量,少有定量的三維測量。Borazjani等[4-5]在不同雷諾數Re=ρVd/μ(ρ、V、d、μ分別表示流體密度、游動速度、特征長度、流體黏性)下對鲹科模式 (Carangiform)和鰻鱺科模式(Anguilliform)的游動進行數值模擬,發(fā)現魚游動時尾跡渦環(huán)呈現出單行和雙行特征。鰻鱺科魚在游動時由于尾部振幅較小,以致于脫落的渦環(huán)形成一條直線,稱為單行渦環(huán);而對于魚尾振幅較大的魚類游動,尾跡則形成雙行脫落的渦環(huán)。中國科學技術大學李龍等[6]對鲹科模式的魚類進行了魚尾模型實驗,發(fā)現魚尾轉動和平動的相位差會顯著影響魚尾的推力。藍鰓太陽魚的胸鰭運動模式是MPF模式中的刺鲀科模式,Lauder等[7]對藍鰓太陽魚的胸鰭進行PIV實驗研究,發(fā)現魚胸鰭流場的三維渦結構呈現出單個渦環(huán)的特征。

    魚類的游動模式與游動狀態(tài)密切相關。魚游動一般可簡化為2種游動狀態(tài):巡游狀態(tài)(Cruise)和猝發(fā)狀態(tài)(Burst)[8]。魚類巡游狀態(tài)方面的研究有很多,并且也取得了很多重要成果。Mwaffo等[8]通過平面PIV技術對斑馬魚在水中巡游時魚尾擺動頻率、游動速度、頭部移動方向進行了相關研究,發(fā)現頻率、速度、相位角在魚游動時相互影響。Mchenry等[9]則進行了成年斑馬魚的巡游時間、啟動速度、慣性等3個方面的實驗研究。模擬仿生方面,活魚巡游狀態(tài)一般定義為魚尾的正弦運動[10],因此在仿生魚尾的實驗時要控制仿生魚尾進行正弦擺動,從而可以初步研究機械魚的水動力學性能[11]。哈爾濱工程大學于凱等[12]進行了雙尾的仿生推進器的水動力學實驗,結果表明相對于單尾的仿生推進器,雙尾布局能夠提高推進器的穩(wěn)定性和靈活性。

    相對于巡游狀態(tài),魚猝發(fā)(Burst)狀態(tài)的實驗研究則成為國內外研究的焦點。魚的猝發(fā)運動狀態(tài)也稱為逃逸反應(Escape response)、C型啟動逃逸(C-start escape response)和C型轉彎[13-15]。實驗方面,Danos等[15]在不同的流體黏度情況下對魚類的逃逸反應進行了實驗研究,發(fā)現高黏度的流體對魚類逃逸反應的某一個階段有強烈的阻礙作用。Danos等[16]就斑馬魚轉彎時鰭的各種形態(tài)變化進行了實驗分析,揭示了斑馬魚的雙鰭、中鰭和骨架的形態(tài)變化情況是轉彎游動的重要因素。數值模擬方面,2009年Wang等[17]發(fā)現C型啟動不僅發(fā)生在逃逸反應時還發(fā)生在非逃逸反應時。藍鰓太陽魚逃逸反應的三維數值模擬結果清晰地顯示了尾流的渦結構[14],該結果可以說明魚類游動狀態(tài)的改變是尾流場變化的重要因素。C型游動模式的研究大多還停留在數值模擬和仿生學的研究上,定量的PIV實驗研究則很少有人進行。

    本文利用平面PIV技術和Tomo-PIV技術分別測量斑馬魚在不同游動狀態(tài)下的二維尾跡流場和曼龍魚巡游時的三維尾跡流場結構,并進行了分析。

    1 實驗設置及測量方案

    1.1 平面PIV實驗方案

    平面PIV實驗對象為成年黑色斑馬魚,體長BL=40 mm。實驗用水箱大小為400 mm×200 mm×200 mm(長×寬×高),水箱中用塑料板隔出寬50 mm的流道作為實驗測量區(qū)域。選用PH=7.2的純凈水,水溫在26±0.5 ℃(正常的熱帶魚生存溫度)。實驗前6 h魚不用進食。

    平面PIV實驗布置如圖2所示。水箱中均勻撒播20 μm的空心玻璃微珠作為示蹤粒子。為了照亮示蹤粒子,本實驗采用200 mJ的雙脈沖激光發(fā)生器(鐳寶Vlite-200)作為光源。其發(fā)出的激光波長為532 nm, 脈沖頻率為1 kHz,經擴束光路被擴展成為2 mm厚的片光源。擴束光路采用焦距50 mm凸透鏡和焦距200 mm凹透鏡的組合。圖像采集設備為一臺配置焦距200 mm尼康鏡頭的CMOS高速相機,分辨率為1024 pixel×1024 pixel。相機布置在水箱底部向上拍攝魚游流場。測量區(qū)域如圖2所示,平行于底壁且位于魚身中間。實驗時,當斑馬魚出現在拍攝區(qū)域,并且片光源照射在魚身中間截面時高速相機開始工作。實驗中對每種游動模式連續(xù)采集150~250幀粒子圖像,以便進行瞬時場分析。圖像放大率為0.03 mm/pixel, 對應測量區(qū)域大小為30 mm×30 mm。

    圖2 平面PIV實驗示意圖

    本實驗采用北京立方天地科技有限公司的MicroVec軟件進行PIV圖像分析。查詢窗口設置為32 pixel×32 pixel,步長8 pixel×8 pixel。最終速度場包含125×125個矢量。為研究流場結構與游動狀態(tài)的關系,需從粒子圖像中提取斑馬魚的二維形態(tài)特征。

    根據實驗圖像的灰度分布,在Matlab軟件里采用閾值分割法提取粒子圖像中斑馬魚的二維形態(tài)。首先對圖像進行中值濾波處理,剔除圖像上的示蹤粒子;然后采用閾值分割的方法提取平面斑馬魚形態(tài)特征。效果如圖3所示,其中(a)為原始圖像,(b)為提取的二維形態(tài)圖像。

    圖3 二維形態(tài)圖

    1.2 三維Tomo-PIV實驗方案

    Tomo-PIV實驗對象為成年曼龍魚,BL=60 mm,選擇曼龍魚作為三維實驗對象是因為曼龍魚體型大,尾跡流場三維結構明顯,形態(tài)紋理豐富,有利于進行三維形態(tài)重構;與之相比,斑馬魚體型較小,紋理較少,三維流動結構不明顯。實驗在環(huán)形水槽中進行,實驗段大小為400 mm×200 mm×100 mm(長×高×寬);在拍攝區(qū)域前端懸掛50 g魚食,選用PH=7.2的純凈水,水溫為26±0.5 ℃(正常的熱帶魚生存溫度);實驗前6 h魚不用進食。

    實驗時,環(huán)形水槽中均勻散播20~50 μm的空心玻璃珠作為示蹤粒子。本實驗同樣采用了200 mJ的固體激光器。與平面PIV實驗所不同的是,脈沖激光通過導光臂和擴束頭被擴展成厚度為60 mm的體光源。圖像采集設備為4臺配置45 mm 尼康移軸鏡頭的CCD相機,分辨率為2048 pixel×2048 pixel,采集頻率為10 Hz。整體實驗布置如圖4所示,4臺CCD相機呈一字型布置在觀測區(qū)域一側,調節(jié)相機的位置,測量魚游的尾跡流場。

    圖4 Tomo-PIV實驗示意圖

    Tomo-PIV的圖像分析采用MicroVec標準算法進行。三維重構體大小為60 mm×50 mm×50 mm,包含1090×909×909個體素。計算時采用三維體變形算法,迭代2次。第一次窗口大小為64 pixel×64 pixel×64 pixel,第二次窗口大小為48 pixel×48 pixel×48 pixel,最終步長為24 pixel×24 pixel×24 pixel。得到的速度場經過歸一化中值檢測、濾波等后處理,并進一步計算出旋渦強度等物理量以便于分析。為了同時提取魚的形態(tài)特征,根據雙目立體成像的原理[18]和標定的結果,從其中2個相機的粒子圖像中還原出魚運動的瞬時形態(tài)。圖5為三維重構實驗數據與重構結果圖,其中(a)、(b)為2個不同視角的相機圖像;在雙目立體成像原理和標定擬合的多項式公式下,通過圖像上曼龍魚形態(tài)上的獨有的特征點信息,重構出曼龍魚三維形態(tài),結果如圖5(c)所示。

    圖5 三維重構實驗數據與重構結果圖

    2 實驗結果及分析

    本節(jié)主要分析與斑馬魚不同游動狀態(tài)相對應的二維尾流場的渦結構,以及曼龍魚在巡游狀態(tài)下的三維尾渦結構。斑馬魚游動狀態(tài)分為啟動狀態(tài)、巡游狀態(tài)、巡游-C型逃逸過渡狀態(tài)、C型逃逸狀態(tài)。根據Williamso等[19]的結論,尾跡渦結構可以分為P模式、S模式、P+S模式。各模式對應的渦脫落模式如圖6所示。P模式為尾跡流場在半個周期內有1對運動方向相反的旋渦一同耗散消失;S模式為尾跡流場在半個或者1個周期內僅有1個旋渦出現,并且耗散消失;P+S模式為尾跡流場在1個周期內出現1對運動方向相反的旋渦,并且緊隨其后還會出現1個單一的旋渦。除了圓柱渦激振動實驗尾跡流場渦脫落模式外,也有學者在圓柱或機翼后安裝柔性平板,從而分析不同運動狀態(tài)下的渦結構[20-21],以上研究結果將被用于分析斑馬魚的尾跡流場。

    2.1 斑馬魚啟動狀態(tài)測量結果

    圖7為斑馬魚在水中從靜止狀態(tài)到游動狀態(tài)時刻的連續(xù)瞬態(tài)渦量場云圖,(a)~(h)分別為0、30、60、90、120、150、180和210 ms時刻的渦量云圖。從圖中可以看出:斑馬魚從靜止狀態(tài)到啟動狀態(tài),尾部輕微向下擺動,此時尾部后方出現了1個脫體渦1,同一時刻魚尾前部還存在1個前緣渦2(見圖7(c)),前緣渦在魚尾移動時逐漸增長。這2個渦配對表現出反卡門渦的特征,形成向后射流。當魚尾向下到達最大位移處后再向上移動時,脫體渦仍存在,前緣渦增長到最強并開始脫落,此時前緣渦2完全脫落,與之前的渦1運動方向相反(見圖7(e)),這2個渦結合在一起形成蘑菇狀的渦結構。在半個周期內流場中僅有1對渦脫落并且耗散,這就構成了渦脫落的P模式。這2個渦在形成以后逐漸耗散消失,之后魚尾在擺回至中心位置的半個周期內,仍會出現相同的渦變化特征。

    圖6 渦脫落模式圖

    圖7 啟動狀態(tài)渦量場云圖

    魚瞬態(tài)速度的計算以魚在圖像中質心處的速度為基準,如圖8所示。在啟動狀態(tài)下,總體而言魚游動速度表現為上升形式。剛開始啟動時的速度為0.305 BL/s,最大速度到達0.375 BL/s。魚尾的擺動在50和210 ms時為主要的運動狀態(tài),此時速度略微下降;正是此時魚尾的擺動為魚后面的運動提供了動量,使魚的速度快速增加。

    圖8 啟動狀態(tài)速度變化圖

    2.2 斑馬魚巡游狀態(tài)測量結果

    圖9為斑馬魚在水中巡游狀態(tài)時刻的連續(xù)瞬態(tài)渦量場云圖,(a)~(h)分別為0、10、20、30、40、50、60和70 ms時刻的渦量云圖。從圖中可以看出,斑馬魚在巡游狀態(tài)持續(xù)游動時,尾巴連續(xù)擺動,以魚水平中線為起點,魚尾向上擺動時尾后流場會出現1個前緣渦1(見圖9(c)),此渦對魚尾有推動作用,且當魚尾運動到最大位移處時該渦增強之后開始脫落;當魚尾開始向下擺動時,周圍會出現1個前緣渦2(見圖9(f)),而當魚尾恢復到接近水平位置時,前緣渦2增強后脫落,并與之前脫落的渦1形成P模式的渦脫落結構(見圖9(g))。兩脫落渦構成反卡門渦,并且都是由前緣渦生長而來,因此在魚尾運動的半個周期內均會出現相同的渦變化特征。

    圖9 巡游狀態(tài)渦量場云圖

    如圖10所示,魚在巡游過程中速度一直維持在0.650~1.125 BL/s。在60 ms處尾跡流場正好脫落成1對P模式的旋渦結構,此時游動速度增加到最大,說明魚尾前期的擺動為魚的游動提供了推力。在整個巡游期間,魚游動的平均速度為0.3625 BL/s(0.0145 m/s),魚尾擺動的平均振幅A約為5 mm。由于每個瞬時點間隔10 ms,因此魚尾擺動平均頻率f約為8.57 Hz,可知斯特勞哈爾數為2.96(St=lf/v,l為振幅,f為頻率,v為速度,St數為表征流動非定常性的相似準則)。

    圖10 巡游狀態(tài)速度變化圖

    2.3 斑馬魚巡游-C型逃逸過渡狀態(tài)測量結果

    圖11為斑馬魚在水中巡游-C型逃逸過渡狀態(tài)時刻的連續(xù)瞬態(tài)渦量場的云圖,(a)~(h)分別為0、20、40、60、80、100、120和140 ms時刻的渦量云圖。斑馬魚從巡游狀態(tài)轉向C型逃逸狀態(tài)的過程中,游動方向及游動速度未達到逃逸狀態(tài)時稱為過渡狀態(tài)。從圖11可以看出,斑馬魚在巡游-C型逃逸過渡狀態(tài)時,魚尾擺動幅度大、方向變化大。以魚尾開始擺動作為起始點,其向右上方擺動時周圍出現了脫體渦1(見圖11(b)),在魚尾逐漸到達最大位移處的過程中脫體渦1脫落,魚尾進而向左下方擺動,此時出現前緣渦2(見圖11(c)),以上兩渦構成了反卡門渦。在魚尾向左下方擺動的過程中前緣渦2脫落,而魚尾向左下方擺動到接近最大位移處時,又出現另一個前緣渦3(見圖11(d)),在魚繼續(xù)向右上方游動時前緣渦3脫落。整個周期內魚尾流場出現3個脫落渦,并且有1對渦同步耗散,另一個渦單獨耗散,因此表現出P+S模式的渦脫落結構。流場相對較為混亂,但總體的渦結構還是能看出來的。

    魚的運動角度是以魚在圖像中的質心位置至頭部方向的直線為基準,計算其與水平方向形成的夾角大小。從圖12可以看出,整個運動過程中魚的轉動角度為121°±1°,平均角速度為12.53 rad/s。在40 ms處魚的轉動角度最大,正好與圖12(c)處魚游動姿態(tài)成C形狀相對應。

    2.4 斑馬魚C型狀態(tài)測量結果

    圖13為斑馬魚在水中C型逃逸狀態(tài)時刻的連續(xù)瞬態(tài)渦量場云圖,(a)~(h)分別為0、40、80、100、120、160、200和240 ms時刻的渦量云圖。C型逃逸過程中魚形態(tài)從C型樣式變?yōu)橹本€樣式。從圖13可以看出,在C型逃逸狀態(tài)時,魚游動時轉動360°。以魚水平中線為起始點,魚尾在擺動時首先出現1個脫體渦1,由于此時魚尾速度為零,所以這個渦稱為啟動渦[11]。隨后由于魚身體的轉動和魚尾的擺動,魚尾附近的流場流動附著在魚尾的邊緣,流動方向與魚尾的運動方向大致相同,以至于引起一連串平面上看起來混亂的渦流動。由于整個周期內尾流場內部僅有1個旋渦脫落并且逐漸耗散消失,所以渦脫落的模式為S。如圖14所示,在整個過程中魚轉動角度為193°±1°,平均角速度為21.54 rad/s,在40和90 ms處各有1個峰值,說明在這兩處旋轉角變化最大,這也是魚形態(tài)呈C型的運動時刻。

    圖11 巡游-C型逃逸過渡狀態(tài)渦量場云圖

    圖12 巡游-C型逃逸過渡狀態(tài)角度變化圖

    圖13 C型逃逸狀態(tài)渦量場云圖

    圖14 C型逃逸狀態(tài)角度變化圖

    2.5 三維測試結果

    為便于確定三維尾跡與魚尾的相對位置,將三維流場和魚的形態(tài)疊加顯示,如圖15(a)所示。圖中可以清晰地看到曼龍魚巡游時尾流的三維結構,其中出現了經典的連續(xù)渦環(huán)鏈[22]。此時流動速度為30 mm/s。圖15(a)中的速度場切片云圖是根據渦環(huán)流向方向所畫的,此云圖表明通過每個渦環(huán)的射流不是水平的。流場內出現了相互連接的3個橢球形渦環(huán),與Flammang等[22]的實驗結果一致,所不同的是本實驗中魚尾擺動幅度較小、游動速度較低。

    巡游狀態(tài)下曼龍魚尾跡流場渦判據準則λci結構如圖15(b)所示。在游動期間內二維切面上有4個可見旋渦(V1、V2、V3、V4)、3個射流(J1、J2、J3)。圖15(c)為三維結構相互連接的渦環(huán)和射流的二維示意圖,根據渦環(huán)和射流的位置在圖中用虛線給出了魚運動的路徑,其中黑色的粗箭頭與圖15(b)中綠色的粗箭頭代表了射流的方向。

    為方便描述,將圖15(c)內的魚運動尾跡分為3個不同時刻T1、T2、T3。在T1時刻內魚尾擺動,產生射流J1和2個旋渦V1、V2,顯然渦環(huán)1(參考圖15(a))是非對稱的單獨渦環(huán),旋渦V1的強度略大于旋渦V2。射流J1的非均勻分布為魚質心的移動提供推力。

    在T2時刻旋渦V3出現,射流J2是由一對方向相反的旋渦V3、V4形成的。射流J2首先由Sakakibara等[23]發(fā)現,并命名為“側邊射流”。魚尾的擺動方向與射流的流動方向成一定的角度,從而獲得了角動量。從圖15(b)中可以發(fā)現旋渦V3轉動的方向與旋渦V2相同,與旋渦V1更接近,表明旋渦V3對射流J1的流動有推動作用。

    T3時刻僅能看到射流J3的存在,根據Ting等[24]的研究發(fā)現:魚進入巡游狀態(tài)時,魚尾產生射流J3的同時還伴有2個持續(xù)的旋渦。所以此時應該還有1個旋渦V5,可能是因為魚尾對三維粒子重構有影響,未能計算出來。從圖15(b)中可以發(fā)現魚尾與渦環(huán)3(參考圖15(a))等高,表明該渦環(huán)是尾鰭的尾緣渦。

    圖15 曼龍魚巡游狀態(tài)三維尾流場結構

    從魚尾的路徑可以發(fā)現所有的射流與尾部的方向均存在一定夾角,未出現射流方向與魚尾平行的情況,表明在巡游狀態(tài)下射流僅有部分動能為魚的游動提供動力。與二維結構相比,三維結構更能清晰地說明魚的尾跡流場特點。

    3 結 論

    本文以斑馬魚、曼龍魚為研究對象,針對游動的不同狀態(tài),利用平面PIV、Tomo-PIV技術分別測量了不同游動模式下的斑馬魚二維尾跡流場和曼龍魚巡游時的三維尾跡流場,根據瞬時渦量場的變化情況分析魚尾渦結構的變化規(guī)律。結果表明:

    (1) 斑馬魚在啟動狀態(tài)、巡游狀態(tài)下,半個周期內均出現2個旋渦,分別為起始時魚尾擺動時產生的渦、魚尾最大位移后脫落的渦。

    (2) 巡游-C型逃逸過渡狀態(tài)下整個周期內尾鰭流場出現3個旋渦,由于流場擾動導致最后一個渦先耗散。

    (3) 斑馬魚在C型逃逸狀態(tài)下,開始改變游動方向時,僅觀察到1個較為明顯的渦環(huán),此時魚尾外圍速度為零,產生的渦為啟動渦。

    (4) 根據渦脫落的不同特征,靜水中斑馬魚在啟動狀態(tài)、巡游狀態(tài)下渦脫落模式均為P,巡游-C型逃逸過渡狀態(tài)下渦脫落模式為P+S,而在C型逃逸狀態(tài)下,由于轉動方向過程中魚尾部流場較為混亂,渦脫落模式更為接近S。

    (5) 曼龍魚尾跡流場的三維結構較為清晰地說明,魚在巡游狀態(tài)下尾跡流動出現了經典三維連續(xù)環(huán)渦鏈結構。

    通過斑馬魚、曼龍魚在水中游動的PIV實驗與分析,對不同游動狀態(tài)下尾跡流場渦結構變化規(guī)律有了基本的認識。該研究為魚類推進的水動力學機理分析打下基礎,也為魚類游動的數值計算提供了驗證依據。下一步將進行更為深入的二維旋渦和三維渦環(huán)的定量研究。

    猜你喜歡
    尾跡魚尾旋渦
    一種基于Radon 變換和尾跡模型的尾跡檢測算法
    小心,旋渦來啦
    觀獅身魚尾像
    大班科學活動:神秘的旋渦
    旋渦笑臉
    山間湖
    揚子江(2019年1期)2019-03-08 02:52:34
    緊體魚尾裙魚尾展開處和圍度設置探討
    基于EEMD-Hilbert譜的渦街流量計尾跡振蕩特性
    蝸牛去拜訪蝴蝶的時候(外一首)
    少年文藝(2015年6期)2015-05-28 04:40:00
    紙杯變小魚
    黄色丝袜av网址大全| 两个人免费观看高清视频| 黄色女人牲交| 国产成人精品久久二区二区91| 欧美人与性动交α欧美精品济南到| 人妻久久中文字幕网| 韩国av一区二区三区四区| 麻豆一二三区av精品| 天天添夜夜摸| 国产亚洲精品久久久久久毛片| 成人18禁高潮啪啪吃奶动态图| 精品久久久久久久久久久久久 | 老鸭窝网址在线观看| 十八禁人妻一区二区| 啪啪无遮挡十八禁网站| 日本 欧美在线| 男人的好看免费观看在线视频 | 91成年电影在线观看| av有码第一页| 久久精品国产清高在天天线| 国产一区二区在线av高清观看| 亚洲中文av在线| 男女之事视频高清在线观看| 两性夫妻黄色片| 十八禁人妻一区二区| 成年女人毛片免费观看观看9| 国产亚洲精品av在线| 99久久久亚洲精品蜜臀av| 黄片小视频在线播放| 午夜福利在线观看吧| 欧美日韩亚洲国产一区二区在线观看| 老司机在亚洲福利影院| 黄色女人牲交| 俺也久久电影网| 国产精品永久免费网站| 久久香蕉精品热| 国产免费av片在线观看野外av| 亚洲aⅴ乱码一区二区在线播放 | 变态另类成人亚洲欧美熟女| 亚洲第一青青草原| 后天国语完整版免费观看| 最新美女视频免费是黄的| 欧美精品亚洲一区二区| 黄片小视频在线播放| 桃红色精品国产亚洲av| 高潮久久久久久久久久久不卡| 亚洲三区欧美一区| 久久久水蜜桃国产精品网| 成年人黄色毛片网站| 国产精品永久免费网站| 中文字幕高清在线视频| www.精华液| 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看黄色视频的| 18禁黄网站禁片午夜丰满| 亚洲精品av麻豆狂野| 欧美一区二区精品小视频在线| 一二三四社区在线视频社区8| 一区二区三区精品91| 观看免费一级毛片| 日本免费一区二区三区高清不卡| 成在线人永久免费视频| 欧美成人免费av一区二区三区| 老司机深夜福利视频在线观看| 国产三级在线视频| 日韩av在线大香蕉| 一本精品99久久精品77| 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出| 久久久精品国产亚洲av高清涩受| 老司机深夜福利视频在线观看| 国产精品影院久久| 欧美激情 高清一区二区三区| 90打野战视频偷拍视频| 两个人视频免费观看高清| 一本一本综合久久| 欧美成人性av电影在线观看| 欧美日韩乱码在线| 嫁个100分男人电影在线观看| 国产亚洲欧美精品永久| 国产伦一二天堂av在线观看| www.自偷自拍.com| 黄色女人牲交| 国产欧美日韩一区二区三| 国产亚洲欧美在线一区二区| 黄频高清免费视频| 亚洲五月天丁香| 欧美色欧美亚洲另类二区| 亚洲欧洲精品一区二区精品久久久| 看黄色毛片网站| 国产精品免费视频内射| 18禁国产床啪视频网站| 麻豆国产av国片精品| 啦啦啦韩国在线观看视频| 国产高清视频在线播放一区| 欧美日韩亚洲综合一区二区三区_| 亚洲免费av在线视频| 国产伦一二天堂av在线观看| 久久天堂一区二区三区四区| 99久久国产精品久久久| 又黄又爽又免费观看的视频| 香蕉av资源在线| 十八禁人妻一区二区| 欧美一级a爱片免费观看看 | 19禁男女啪啪无遮挡网站| 午夜激情福利司机影院| 超碰成人久久| 黄片大片在线免费观看| 老熟妇仑乱视频hdxx| 日韩欧美一区视频在线观看| 欧美性猛交黑人性爽| cao死你这个sao货| 国产av又大| 国产精品 欧美亚洲| 波多野结衣巨乳人妻| 亚洲熟女毛片儿| 亚洲国产欧美日韩在线播放| 一a级毛片在线观看| 久久久精品欧美日韩精品| 一本久久中文字幕| 国产精品美女特级片免费视频播放器 | 黄色视频,在线免费观看| 大型av网站在线播放| 中文字幕精品亚洲无线码一区 | 男人舔女人的私密视频| 在线永久观看黄色视频| 一进一出抽搐动态| 欧美 亚洲 国产 日韩一| 国内揄拍国产精品人妻在线 | www日本在线高清视频| 草草在线视频免费看| 高清毛片免费观看视频网站| 韩国av一区二区三区四区| 美女午夜性视频免费| 色精品久久人妻99蜜桃| 一级毛片高清免费大全| 巨乳人妻的诱惑在线观看| 久久天躁狠狠躁夜夜2o2o| 久久热在线av| 在线视频色国产色| 青草久久国产| 久久久久久久久中文| 亚洲国产精品999在线| 1024视频免费在线观看| 国产亚洲精品综合一区在线观看 | av视频在线观看入口| 伦理电影免费视频| 18禁观看日本| 欧美黑人巨大hd| 日韩大尺度精品在线看网址| 18禁黄网站禁片免费观看直播| 亚洲精华国产精华精| 啦啦啦观看免费观看视频高清| 亚洲全国av大片| 国产伦人伦偷精品视频| 欧美一级a爱片免费观看看 | 麻豆av在线久日| 人妻丰满熟妇av一区二区三区| 自线自在国产av| 大型av网站在线播放| 亚洲在线自拍视频| 国产亚洲精品第一综合不卡| 99riav亚洲国产免费| 国产激情偷乱视频一区二区| 国产欧美日韩一区二区三| 最近最新免费中文字幕在线| 日韩一卡2卡3卡4卡2021年| 在线播放国产精品三级| 狂野欧美激情性xxxx| 男男h啪啪无遮挡| 亚洲三区欧美一区| 男女床上黄色一级片免费看| 欧美一级毛片孕妇| 午夜福利高清视频| a级毛片在线看网站| 国产99久久九九免费精品| 国产野战对白在线观看| a在线观看视频网站| 麻豆一二三区av精品| 国产真人三级小视频在线观看| 狠狠狠狠99中文字幕| 91麻豆av在线| 午夜视频精品福利| 老鸭窝网址在线观看| 可以在线观看毛片的网站| 亚洲第一av免费看| 91国产中文字幕| 精品国产国语对白av| 草草在线视频免费看| 久久精品成人免费网站| 国产v大片淫在线免费观看| 最新美女视频免费是黄的| 国产伦一二天堂av在线观看| 女人高潮潮喷娇喘18禁视频| 91成年电影在线观看| 热99re8久久精品国产| 久久国产亚洲av麻豆专区| 成人国产一区最新在线观看| 又大又爽又粗| 色播在线永久视频| 男人舔奶头视频| 在线天堂中文资源库| 亚洲国产精品成人综合色| 在线av久久热| 看黄色毛片网站| 大型黄色视频在线免费观看| 欧美午夜高清在线| 午夜老司机福利片| 一二三四在线观看免费中文在| 国产三级黄色录像| 亚洲天堂国产精品一区在线| 国产一区在线观看成人免费| 两性夫妻黄色片| 两性午夜刺激爽爽歪歪视频在线观看 | 男女做爰动态图高潮gif福利片| 国产欧美日韩一区二区精品| 亚洲专区国产一区二区| 黄频高清免费视频| 亚洲av片天天在线观看| 国产精品二区激情视频| 久久伊人香网站| 亚洲全国av大片| 欧美黑人巨大hd| 亚洲中文日韩欧美视频| 18禁国产床啪视频网站| 熟女少妇亚洲综合色aaa.| 国产精品综合久久久久久久免费| www.www免费av| 搡老熟女国产l中国老女人| 精品国产美女av久久久久小说| 免费在线观看成人毛片| 国产一区二区三区视频了| 国产高清有码在线观看视频 | 天天一区二区日本电影三级| 午夜福利免费观看在线| 午夜久久久久精精品| 日本五十路高清| 黄片大片在线免费观看| 777久久人妻少妇嫩草av网站| 成人av一区二区三区在线看| 成人免费观看视频高清| 久久人妻av系列| 成人三级黄色视频| 亚洲人成伊人成综合网2020| 久久精品91蜜桃| 狠狠狠狠99中文字幕| 丁香六月欧美| 99精品欧美一区二区三区四区| 成年女人毛片免费观看观看9| 一区二区三区高清视频在线| 91国产中文字幕| 久久久久国产精品人妻aⅴ院| 日韩三级视频一区二区三区| 国产精品一区二区免费欧美| 久久精品国产亚洲av香蕉五月| 久久久精品国产亚洲av高清涩受| 少妇粗大呻吟视频| 亚洲成av人片免费观看| 白带黄色成豆腐渣| 日本免费a在线| 50天的宝宝边吃奶边哭怎么回事| 色综合亚洲欧美另类图片| 十分钟在线观看高清视频www| 国产一级毛片七仙女欲春2 | 国产精品亚洲一级av第二区| 色在线成人网| 在线观看免费视频日本深夜| 一区二区日韩欧美中文字幕| 又大又爽又粗| 久久久久久久久免费视频了| 女人高潮潮喷娇喘18禁视频| 亚洲精品久久国产高清桃花| 日韩大尺度精品在线看网址| 此物有八面人人有两片| 国产三级黄色录像| 亚洲人成伊人成综合网2020| www日本在线高清视频| 国内揄拍国产精品人妻在线 | 一个人免费在线观看的高清视频| 成年免费大片在线观看| 波多野结衣巨乳人妻| 国产又色又爽无遮挡免费看| 在线看三级毛片| 999久久久精品免费观看国产| 亚洲aⅴ乱码一区二区在线播放 | 男人舔女人下体高潮全视频| 老司机午夜十八禁免费视频| 大型黄色视频在线免费观看| 日韩有码中文字幕| 少妇粗大呻吟视频| 国产熟女午夜一区二区三区| 麻豆成人av在线观看| 欧美日韩亚洲综合一区二区三区_| 国产激情欧美一区二区| 99在线人妻在线中文字幕| 国产精品永久免费网站| 一区二区日韩欧美中文字幕| 亚洲男人的天堂狠狠| 久久久国产欧美日韩av| 十分钟在线观看高清视频www| 两个人免费观看高清视频| 亚洲精品色激情综合| 51午夜福利影视在线观看| 欧美色欧美亚洲另类二区| 12—13女人毛片做爰片一| 18禁裸乳无遮挡免费网站照片 | 午夜福利欧美成人| 日韩一卡2卡3卡4卡2021年| 悠悠久久av| a在线观看视频网站| 成年版毛片免费区| 大型av网站在线播放| 成人永久免费在线观看视频| 久久中文字幕一级| 欧美乱色亚洲激情| 日韩欧美免费精品| 国产成+人综合+亚洲专区| 久久久国产精品麻豆| 午夜免费鲁丝| 天堂动漫精品| 国内精品久久久久精免费| 亚洲激情在线av| 亚洲av电影不卡..在线观看| 啦啦啦免费观看视频1| 97超级碰碰碰精品色视频在线观看| 精品电影一区二区在线| 一进一出抽搐动态| 一边摸一边做爽爽视频免费| 国产一区二区三区在线臀色熟女| 亚洲国产精品久久男人天堂| 精品日产1卡2卡| 香蕉久久夜色| 国产精品一区二区三区四区久久 | 亚洲国产日韩欧美精品在线观看 | 成年免费大片在线观看| 麻豆成人av在线观看| 一个人免费在线观看的高清视频| www.精华液| 波多野结衣高清无吗| www.www免费av| 午夜两性在线视频| 在线观看舔阴道视频| 国产成人欧美| 一个人观看的视频www高清免费观看 | 伊人久久大香线蕉亚洲五| 又大又爽又粗| 夜夜看夜夜爽夜夜摸| 日日干狠狠操夜夜爽| 91九色精品人成在线观看| 亚洲精品粉嫩美女一区| 亚洲 欧美一区二区三区| 久久久久久久久久黄片| 91老司机精品| 国产精品一区二区三区四区久久 | 久久久久精品国产欧美久久久| 国产av一区在线观看免费| 国产精品自产拍在线观看55亚洲| 精品免费久久久久久久清纯| 99久久国产精品久久久| 国产亚洲精品一区二区www| 久久精品国产亚洲av香蕉五月| 久久午夜亚洲精品久久| 亚洲午夜理论影院| 亚洲国产看品久久| a级毛片在线看网站| 国产激情欧美一区二区| 一级a爱片免费观看的视频| 中亚洲国语对白在线视频| 丁香六月欧美| 天天躁狠狠躁夜夜躁狠狠躁| 久久亚洲真实| 高潮久久久久久久久久久不卡| 男人舔女人的私密视频| 搡老妇女老女人老熟妇| 高清在线国产一区| 国产爱豆传媒在线观看 | 啦啦啦 在线观看视频| 在线观看日韩欧美| 国产1区2区3区精品| 亚洲av成人不卡在线观看播放网| 日本a在线网址| 黄色a级毛片大全视频| 国产男靠女视频免费网站| 最近最新免费中文字幕在线| 成人精品一区二区免费| 啦啦啦免费观看视频1| 好男人在线观看高清免费视频 | 亚洲av日韩精品久久久久久密| 亚洲男人的天堂狠狠| 一区二区三区高清视频在线| 亚洲一码二码三码区别大吗| 午夜久久久在线观看| 国产精品亚洲美女久久久| 91成人精品电影| 美女高潮到喷水免费观看| 久久中文字幕一级| 男人的好看免费观看在线视频 | 亚洲精华国产精华精| 一卡2卡三卡四卡精品乱码亚洲| 色婷婷久久久亚洲欧美| 国产不卡一卡二| 1024视频免费在线观看| 久久热在线av| 一进一出抽搐gif免费好疼| 琪琪午夜伦伦电影理论片6080| 美女 人体艺术 gogo| 又黄又粗又硬又大视频| 一a级毛片在线观看| 中国美女看黄片| 日本在线视频免费播放| 少妇被粗大的猛进出69影院| 91大片在线观看| 男女床上黄色一级片免费看| 亚洲av中文字字幕乱码综合 | 国产乱人伦免费视频| 男女视频在线观看网站免费 | 精品国产美女av久久久久小说| 91大片在线观看| 精品国产美女av久久久久小说| 国产又黄又爽又无遮挡在线| 亚洲熟妇中文字幕五十中出| 色综合婷婷激情| 久久人妻福利社区极品人妻图片| 亚洲最大成人中文| 男女床上黄色一级片免费看| 淫妇啪啪啪对白视频| 亚洲熟妇中文字幕五十中出| 在线国产一区二区在线| 18禁美女被吸乳视频| 欧美国产精品va在线观看不卡| 18禁美女被吸乳视频| 精品福利观看| 禁无遮挡网站| 18禁观看日本| 少妇的丰满在线观看| 国产不卡一卡二| 久久久久精品国产欧美久久久| 亚洲av美国av| 精品久久久久久久人妻蜜臀av| 国产精品自产拍在线观看55亚洲| 久久久久国产精品人妻aⅴ院| 亚洲国产欧美一区二区综合| 亚洲中文字幕一区二区三区有码在线看 | 国产精品影院久久| 久久久久国产一级毛片高清牌| 国产高清激情床上av| 中文在线观看免费www的网站 | 欧美日韩福利视频一区二区| 2021天堂中文幕一二区在线观 | 波多野结衣高清作品| 看片在线看免费视频| 久久中文字幕人妻熟女| 白带黄色成豆腐渣| 桃红色精品国产亚洲av| 成人国产一区最新在线观看| 99久久综合精品五月天人人| 一进一出好大好爽视频| 一级毛片精品| 一夜夜www| 精品福利观看| 欧美色视频一区免费| 日韩精品青青久久久久久| 最新在线观看一区二区三区| 午夜福利一区二区在线看| 久久久久久久久免费视频了| 国产一区二区三区在线臀色熟女| 欧美激情久久久久久爽电影| 男人的好看免费观看在线视频 | 亚洲狠狠婷婷综合久久图片| 激情在线观看视频在线高清| 俺也久久电影网| 午夜福利高清视频| 国产v大片淫在线免费观看| 色综合亚洲欧美另类图片| 日韩大尺度精品在线看网址| 久久人妻福利社区极品人妻图片| 久热爱精品视频在线9| 制服诱惑二区| 欧美日本亚洲视频在线播放| 久久伊人香网站| 亚洲男人天堂网一区| 一级毛片高清免费大全| 99riav亚洲国产免费| 妹子高潮喷水视频| 桃色一区二区三区在线观看| 午夜福利一区二区在线看| 一级毛片女人18水好多| 极品教师在线免费播放| 国产精品1区2区在线观看.| 后天国语完整版免费观看| 国产av在哪里看| 色av中文字幕| 岛国视频午夜一区免费看| 看黄色毛片网站| 欧美zozozo另类| 精品少妇一区二区三区视频日本电影| 欧美人与性动交α欧美精品济南到| 男女午夜视频在线观看| 久久久久精品国产欧美久久久| 啦啦啦 在线观看视频| 欧美乱色亚洲激情| 俺也久久电影网| 波多野结衣高清无吗| 亚洲avbb在线观看| 亚洲精品美女久久av网站| 亚洲国产日韩欧美精品在线观看 | 18禁观看日本| 欧美日韩亚洲国产一区二区在线观看| 欧美黑人欧美精品刺激| 久久精品亚洲精品国产色婷小说| 俺也久久电影网| 国产精品久久久久久亚洲av鲁大| 熟女少妇亚洲综合色aaa.| 亚洲熟妇中文字幕五十中出| 免费在线观看视频国产中文字幕亚洲| 一边摸一边抽搐一进一小说| 波多野结衣高清无吗| www.精华液| 久久久水蜜桃国产精品网| 不卡av一区二区三区| 男人的好看免费观看在线视频 | 少妇裸体淫交视频免费看高清 | www.自偷自拍.com| 日本三级黄在线观看| 女生性感内裤真人,穿戴方法视频| 欧美性长视频在线观看| 成年人黄色毛片网站| 亚洲成人精品中文字幕电影| 久久人妻av系列| 听说在线观看完整版免费高清| 亚洲欧美精品综合一区二区三区| 99久久国产精品久久久| 婷婷六月久久综合丁香| 真人做人爱边吃奶动态| 亚洲国产欧美一区二区综合| 嫩草影视91久久| 久久伊人香网站| av有码第一页| 亚洲精品一卡2卡三卡4卡5卡| 国产麻豆成人av免费视频| 亚洲片人在线观看| 他把我摸到了高潮在线观看| 国产精品野战在线观看| 久久亚洲真实| 97碰自拍视频| 日本在线视频免费播放| 精品久久久久久久毛片微露脸| 亚洲欧美精品综合一区二区三区| 非洲黑人性xxxx精品又粗又长| av欧美777| 真人做人爱边吃奶动态| 久久精品91无色码中文字幕| 日韩精品青青久久久久久| 国产精品永久免费网站| 黑丝袜美女国产一区| 麻豆成人av在线观看| 欧美色欧美亚洲另类二区| 久久中文看片网| 国产熟女xx| 在线观看免费日韩欧美大片| 成人三级黄色视频| 免费在线观看日本一区| 性色av乱码一区二区三区2| 麻豆久久精品国产亚洲av| 特大巨黑吊av在线直播 | 女性被躁到高潮视频| 在线观看免费午夜福利视频| 黄色女人牲交| 亚洲一区二区三区不卡视频| 真人一进一出gif抽搐免费| 美女午夜性视频免费| 精品国内亚洲2022精品成人| 夜夜爽天天搞| 神马国产精品三级电影在线观看 | 精品一区二区三区视频在线观看免费| 成年女人毛片免费观看观看9| 国产av不卡久久| 欧美精品亚洲一区二区| 99久久精品国产亚洲精品| 亚洲男人天堂网一区| 757午夜福利合集在线观看| 一进一出抽搐动态| 久久精品国产清高在天天线| 亚洲自偷自拍图片 自拍| 国产欧美日韩一区二区精品| 国产高清激情床上av| 黑人操中国人逼视频| 欧美乱妇无乱码| avwww免费| netflix在线观看网站| 桃红色精品国产亚洲av| 国产又爽黄色视频| 一区二区三区国产精品乱码| a在线观看视频网站| 97人妻精品一区二区三区麻豆 | 国产精品,欧美在线| 亚洲专区字幕在线| 99精品久久久久人妻精品| 国产精品久久久人人做人人爽| 成人永久免费在线观看视频| 别揉我奶头~嗯~啊~动态视频| 久久99热这里只有精品18| 国产亚洲精品一区二区www| 久9热在线精品视频| 亚洲精品久久成人aⅴ小说| 9191精品国产免费久久| 亚洲色图 男人天堂 中文字幕| 日本 av在线| 国产亚洲精品一区二区www| 久久久国产成人免费| 一区二区日韩欧美中文字幕| 美国免费a级毛片| 在线观看日韩欧美| 国产一卡二卡三卡精品| 欧美日韩瑟瑟在线播放|