• <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
    紙杯變小魚
    在线观看免费视频日本深夜| 99re在线观看精品视频| 色老头精品视频在线观看| 啦啦啦中文免费视频观看日本| 最新在线观看一区二区三区| 十八禁人妻一区二区| a级毛片在线看网站| 国产伦理片在线播放av一区| 少妇粗大呻吟视频| 欧美精品亚洲一区二区| 精品国产国语对白av| 国产在线精品亚洲第一网站| 亚洲综合色网址| 日韩一卡2卡3卡4卡2021年| 亚洲五月婷婷丁香| 天堂中文最新版在线下载| 国产午夜精品久久久久久| 久久中文看片网| 国产xxxxx性猛交| 亚洲精品在线观看二区| 亚洲专区国产一区二区| 久久99热这里只频精品6学生| 国产成人啪精品午夜网站| 国产成人免费观看mmmm| 在线观看免费午夜福利视频| 久久精品国产99精品国产亚洲性色 | 真人做人爱边吃奶动态| 国产精品 欧美亚洲| 欧美黄色淫秽网站| 亚洲国产av新网站| 亚洲人成电影免费在线| 免费观看a级毛片全部| 中文字幕人妻丝袜一区二区| 18禁观看日本| 色婷婷久久久亚洲欧美| 久久久久久久久免费视频了| 国产成人一区二区三区免费视频网站| 天堂中文最新版在线下载| 欧美激情 高清一区二区三区| 自线自在国产av| 亚洲精品一卡2卡三卡4卡5卡| 黑人欧美特级aaaaaa片| 在线av久久热| 一边摸一边抽搐一进一小说 | 国产成人欧美| 99九九在线精品视频| 国产精品久久电影中文字幕 | 巨乳人妻的诱惑在线观看| 黄色视频在线播放观看不卡| 欧美精品av麻豆av| 日韩有码中文字幕| 国产精品久久久久久人妻精品电影 | 国产一区二区三区在线臀色熟女 | 国产av又大| 丰满饥渴人妻一区二区三| 一区二区三区国产精品乱码| 国产亚洲欧美精品永久| 成人手机av| 欧美激情高清一区二区三区| 人人妻人人爽人人添夜夜欢视频| 99久久人妻综合| 一二三四在线观看免费中文在| 黄色a级毛片大全视频| 日韩大片免费观看网站| 在线观看66精品国产| 亚洲 欧美一区二区三区| 国产精品亚洲av一区麻豆| 亚洲五月婷婷丁香| 天堂中文最新版在线下载| 啦啦啦 在线观看视频| 国产又爽黄色视频| 国产一区二区三区视频了| 欧美性长视频在线观看| 丰满迷人的少妇在线观看| 欧美日韩福利视频一区二区| 一级,二级,三级黄色视频| 老司机影院毛片| 久热这里只有精品99| 日韩欧美国产一区二区入口| 大香蕉久久网| 一区二区av电影网| 久久免费观看电影| a在线观看视频网站| 亚洲人成电影观看| 性高湖久久久久久久久免费观看| 美女午夜性视频免费| 日韩免费av在线播放| 一二三四在线观看免费中文在| 99久久99久久久精品蜜桃| 欧美大码av| 国产精品免费大片| 精品国产乱码久久久久久小说| 中亚洲国语对白在线视频| 精品国产一区二区三区久久久樱花| 一本一本久久a久久精品综合妖精| 精品久久久久久久毛片微露脸| 国产一区有黄有色的免费视频| 久久国产精品影院| 丁香六月欧美| 久久精品国产99精品国产亚洲性色 | 国产不卡av网站在线观看| 亚洲精品在线美女| 亚洲av成人一区二区三| 女人精品久久久久毛片| 最近最新中文字幕大全免费视频| 美女午夜性视频免费| 国产精品久久久久久精品电影小说| 亚洲国产欧美在线一区| 肉色欧美久久久久久久蜜桃| 热99久久久久精品小说推荐| 成人av一区二区三区在线看| av电影中文网址| 最新美女视频免费是黄的| 国产欧美日韩综合在线一区二区| 亚洲国产精品一区二区三区在线| 国产精品 国内视频| 久久性视频一级片| 欧美日韩一级在线毛片| 99精品久久久久人妻精品| 欧美日韩福利视频一区二区| 国产精品成人在线| 亚洲av国产av综合av卡| 最近最新免费中文字幕在线| 母亲3免费完整高清在线观看| 久久99一区二区三区| 国产伦人伦偷精品视频| 黄色视频,在线免费观看| 91成年电影在线观看| 亚洲中文av在线| 久久久水蜜桃国产精品网| 精品少妇黑人巨大在线播放| 美女高潮喷水抽搐中文字幕| 天堂动漫精品| 99在线人妻在线中文字幕 | 国产男靠女视频免费网站| 午夜免费成人在线视频| 欧美精品高潮呻吟av久久| 国产日韩欧美在线精品| avwww免费| 久久久欧美国产精品| 久热这里只有精品99| 9191精品国产免费久久| 高清欧美精品videossex| 亚洲自偷自拍图片 自拍| 夜夜骑夜夜射夜夜干| 下体分泌物呈黄色| 麻豆乱淫一区二区| 男女午夜视频在线观看| 精品少妇久久久久久888优播| 亚洲中文日韩欧美视频| 黄色成人免费大全| 精品国产超薄肉色丝袜足j| 国产深夜福利视频在线观看| avwww免费| 成年人免费黄色播放视频| 亚洲av美国av| 久久久欧美国产精品| 两个人免费观看高清视频| av欧美777| 国产不卡一卡二| 国产xxxxx性猛交| 啦啦啦免费观看视频1| 免费看a级黄色片| 亚洲伊人色综图| av天堂在线播放| 十八禁高潮呻吟视频| 十八禁高潮呻吟视频| 人人妻,人人澡人人爽秒播| www.自偷自拍.com| 天天躁夜夜躁狠狠躁躁| 午夜视频精品福利| kizo精华| 大香蕉久久网| 80岁老熟妇乱子伦牲交| 日韩欧美一区二区三区在线观看 | 精品少妇久久久久久888优播| 黄色a级毛片大全视频| 久久国产精品人妻蜜桃| 精品第一国产精品| 久久热在线av| 高潮久久久久久久久久久不卡| 我的亚洲天堂| 啦啦啦中文免费视频观看日本| 国产精品电影一区二区三区 | 五月开心婷婷网| 欧美黄色片欧美黄色片| 国产精品二区激情视频| 99在线人妻在线中文字幕 | 在线永久观看黄色视频| 国产一卡二卡三卡精品| 天堂动漫精品| 国产精品九九99| 欧美日韩中文字幕国产精品一区二区三区 | 在线av久久热| e午夜精品久久久久久久| 久久精品国产99精品国产亚洲性色 | 欧美日韩黄片免| 午夜视频精品福利| 亚洲伊人色综图| 亚洲国产毛片av蜜桃av| 日韩视频在线欧美| 另类亚洲欧美激情| 精品乱码久久久久久99久播| 露出奶头的视频| 久久中文字幕人妻熟女| 亚洲av日韩在线播放| 人成视频在线观看免费观看| 视频在线观看一区二区三区| 视频区图区小说| 麻豆国产av国片精品| 青草久久国产| 在线观看免费日韩欧美大片| 午夜福利在线观看吧| 最近最新中文字幕大全免费视频| 精品国产一区二区三区久久久樱花| av有码第一页| 国产一区二区激情短视频| 十八禁网站网址无遮挡| 国产精品久久久人人做人人爽| h视频一区二区三区| 精品视频人人做人人爽| 精品国产一区二区三区四区第35| 国产日韩欧美亚洲二区| 人人妻,人人澡人人爽秒播| 成人国产一区最新在线观看| 丰满人妻熟妇乱又伦精品不卡| 成人18禁在线播放| 高清在线国产一区| 两性午夜刺激爽爽歪歪视频在线观看 | 精品免费久久久久久久清纯 | 日本a在线网址| 欧美人与性动交α欧美精品济南到| 777米奇影视久久| 久久精品亚洲精品国产色婷小说| 国产亚洲av高清不卡| 亚洲精品在线观看二区| 十八禁网站免费在线| 久久精品国产亚洲av高清一级| 后天国语完整版免费观看| 成人国语在线视频| 黄片播放在线免费| 欧美日韩国产mv在线观看视频| 亚洲欧美日韩另类电影网站| 丝袜喷水一区| 成年人免费黄色播放视频| 日韩人妻精品一区2区三区| 三级毛片av免费| 黄色丝袜av网址大全| 一区二区三区乱码不卡18| 老汉色∧v一级毛片| 男女无遮挡免费网站观看| 91成年电影在线观看| 亚洲五月婷婷丁香| 免费久久久久久久精品成人欧美视频| 久久 成人 亚洲| 亚洲三区欧美一区| 美女午夜性视频免费| 视频区欧美日本亚洲| 国产精品成人在线| 法律面前人人平等表现在哪些方面| 成人永久免费在线观看视频 | 91成年电影在线观看| 久久九九热精品免费| 国产男女内射视频| 我要看黄色一级片免费的| 亚洲国产欧美日韩在线播放| 亚洲第一青青草原| 久久精品成人免费网站| 多毛熟女@视频| 免费日韩欧美在线观看| 少妇精品久久久久久久| 19禁男女啪啪无遮挡网站| 天堂中文最新版在线下载| 两性夫妻黄色片| 狂野欧美激情性xxxx| 亚洲欧美日韩高清在线视频 | 国产在线视频一区二区| 天堂中文最新版在线下载| 波多野结衣一区麻豆| 精品久久蜜臀av无| 亚洲欧美日韩高清在线视频 | 国产在线免费精品| 交换朋友夫妻互换小说| 女警被强在线播放| 2018国产大陆天天弄谢| 脱女人内裤的视频| 欧美精品亚洲一区二区| 天天躁夜夜躁狠狠躁躁| 一进一出好大好爽视频| 国产免费现黄频在线看| 一边摸一边抽搐一进一小说 | 久久性视频一级片| av又黄又爽大尺度在线免费看| 日本av手机在线免费观看| 亚洲午夜精品一区,二区,三区| 久久久久久久国产电影| 国产免费av片在线观看野外av| 人人妻人人澡人人爽人人夜夜| 韩国精品一区二区三区| 亚洲国产毛片av蜜桃av| 男女免费视频国产| 久久九九热精品免费| 精品一区二区三区四区五区乱码| 可以免费在线观看a视频的电影网站| 亚洲精品国产色婷婷电影| 9色porny在线观看| 男男h啪啪无遮挡| 激情在线观看视频在线高清 | 人人妻人人添人人爽欧美一区卜| 丰满少妇做爰视频| 久久久久久久国产电影| 十八禁网站网址无遮挡| cao死你这个sao货| 午夜激情久久久久久久| 热99国产精品久久久久久7| 一个人免费看片子| 99国产综合亚洲精品| 视频区图区小说| av国产精品久久久久影院| 国产无遮挡羞羞视频在线观看| 国产伦理片在线播放av一区| 这个男人来自地球电影免费观看| 国产精品久久电影中文字幕 | 青草久久国产| 久久ye,这里只有精品| 久久国产精品影院| 少妇被粗大的猛进出69影院| cao死你这个sao货| 人妻久久中文字幕网| 在线观看舔阴道视频| 亚洲少妇的诱惑av| 久久中文字幕人妻熟女| 亚洲熟女精品中文字幕| 黑人猛操日本美女一级片| netflix在线观看网站| 国产日韩欧美视频二区| 成年动漫av网址| 久久中文字幕一级| 久久这里只有精品19| 怎么达到女性高潮| 丝袜在线中文字幕| 欧美激情极品国产一区二区三区| 一级片免费观看大全| 99riav亚洲国产免费| 亚洲综合色网址| 99久久99久久久精品蜜桃| 女同久久另类99精品国产91| 国产精品久久久av美女十八| 欧美日韩中文字幕国产精品一区二区三区 | 丰满少妇做爰视频| 精品国产亚洲在线| 精品乱码久久久久久99久播| 国产精品 国内视频| 黄色丝袜av网址大全| 亚洲色图av天堂| 日韩视频在线欧美| 亚洲精品国产精品久久久不卡| 亚洲一码二码三码区别大吗| 国产极品粉嫩免费观看在线| www.熟女人妻精品国产| 亚洲欧美一区二区三区久久| 日本黄色视频三级网站网址 | 香蕉久久夜色| 在线亚洲精品国产二区图片欧美| 日韩精品免费视频一区二区三区| 国内毛片毛片毛片毛片毛片| 国产成人欧美在线观看 | 交换朋友夫妻互换小说| 极品人妻少妇av视频| 久久国产精品人妻蜜桃| 97人妻天天添夜夜摸| 人人妻,人人澡人人爽秒播| 99久久人妻综合| 精品人妻熟女毛片av久久网站| 亚洲精品自拍成人| 色婷婷av一区二区三区视频| 美女福利国产在线| 午夜精品国产一区二区电影| 亚洲人成电影免费在线| 日韩免费av在线播放| 黄色视频不卡| 免费观看av网站的网址| 欧美日韩视频精品一区| 国产亚洲午夜精品一区二区久久| 国产在视频线精品| 老司机靠b影院| 亚洲国产看品久久| 久久久国产精品麻豆| 麻豆av在线久日| 欧美日韩一级在线毛片| 国产片内射在线| 亚洲色图av天堂| 国产xxxxx性猛交| av网站免费在线观看视频| 午夜激情久久久久久久| 热99国产精品久久久久久7| 男女午夜视频在线观看| 成人免费观看视频高清| 国产极品粉嫩免费观看在线| 99精品在免费线老司机午夜| 丰满饥渴人妻一区二区三| 亚洲av美国av| 亚洲天堂av无毛| 国产精品电影一区二区三区 | 国产精品熟女久久久久浪| 在线观看免费视频网站a站| 国产黄色免费在线视频| xxxhd国产人妻xxx| 午夜免费成人在线视频| 欧美在线一区亚洲| 国产高清国产精品国产三级| 999久久久精品免费观看国产| 岛国毛片在线播放| 老司机午夜十八禁免费视频| www日本在线高清视频| a级片在线免费高清观看视频| 久9热在线精品视频| 欧美日韩国产mv在线观看视频| 日韩精品免费视频一区二区三区| 免费人妻精品一区二区三区视频| 日日夜夜操网爽| 丁香六月天网| 1024视频免费在线观看| 精品国内亚洲2022精品成人 | 91老司机精品| 男人舔女人的私密视频| 狠狠狠狠99中文字幕| 国产日韩一区二区三区精品不卡| 午夜精品国产一区二区电影| 狠狠婷婷综合久久久久久88av| 国产精品欧美亚洲77777| 老汉色∧v一级毛片| 日韩视频一区二区在线观看| 国产男女内射视频| 999精品在线视频| 中文字幕人妻丝袜制服| 99久久精品国产亚洲精品| 一级毛片精品| 日韩 欧美 亚洲 中文字幕| 制服诱惑二区| 老司机亚洲免费影院| 午夜福利,免费看| 亚洲国产欧美在线一区| 国产精品影院久久| 久久精品亚洲av国产电影网| 国产一区二区三区在线臀色熟女 | 国产精品久久久久成人av| 欧美在线黄色| 99riav亚洲国产免费| 中文字幕av电影在线播放| 丝袜喷水一区| 中亚洲国语对白在线视频| 久久青草综合色| 丁香六月天网| 18在线观看网站| 亚洲中文日韩欧美视频| 色在线成人网| 亚洲av国产av综合av卡| 丁香六月欧美| 少妇精品久久久久久久| 在线观看免费日韩欧美大片| 精品高清国产在线一区| 国产淫语在线视频| 后天国语完整版免费观看| 午夜福利影视在线免费观看| 欧美变态另类bdsm刘玥| 亚洲人成电影免费在线| www日本在线高清视频| www.999成人在线观看| 夜夜夜夜夜久久久久| 黑人巨大精品欧美一区二区mp4| 女性生殖器流出的白浆| 欧美另类亚洲清纯唯美| 12—13女人毛片做爰片一| 精品国产乱码久久久久久男人| 大片免费播放器 马上看| 国产91精品成人一区二区三区 | 国产精品久久久久久精品古装| 香蕉丝袜av| 成人18禁高潮啪啪吃奶动态图| 亚洲伊人色综图| 两性夫妻黄色片| 欧美精品亚洲一区二区| 欧美日本中文国产一区发布| 两人在一起打扑克的视频| 亚洲欧美激情在线| 99国产综合亚洲精品| av天堂久久9| 亚洲精品在线观看二区| 99riav亚洲国产免费| 欧美日韩国产mv在线观看视频| kizo精华| 精品国产乱子伦一区二区三区| 91老司机精品| 91成人精品电影| 汤姆久久久久久久影院中文字幕| 王馨瑶露胸无遮挡在线观看| 另类亚洲欧美激情| aaaaa片日本免费| 亚洲三区欧美一区| 久久人妻熟女aⅴ| av国产精品久久久久影院| 国产淫语在线视频| 人人妻人人澡人人看| 亚洲欧美精品综合一区二区三区| 国产av精品麻豆| 色视频在线一区二区三区| 久久国产精品男人的天堂亚洲| tocl精华| 日韩中文字幕欧美一区二区| 男女之事视频高清在线观看| 日日夜夜操网爽| 欧美变态另类bdsm刘玥| 桃花免费在线播放| 不卡一级毛片| 久久亚洲真实| 久久精品国产综合久久久| 亚洲,欧美精品.| av又黄又爽大尺度在线免费看| 午夜视频精品福利| 久久精品91无色码中文字幕| 精品久久久精品久久久| 国产在线视频一区二区| www.熟女人妻精品国产| 亚洲美女黄片视频| 国产欧美日韩综合在线一区二区| 又紧又爽又黄一区二区| 中文字幕另类日韩欧美亚洲嫩草| 精品少妇久久久久久888优播| 午夜福利乱码中文字幕| 国产区一区二久久| 老熟妇仑乱视频hdxx| 国产成人影院久久av| 757午夜福利合集在线观看| 国产成人精品在线电影| 国产在线精品亚洲第一网站| 一区二区日韩欧美中文字幕| 国产成人免费观看mmmm| 欧美日韩亚洲高清精品| 亚洲人成电影观看| 9热在线视频观看99| 亚洲人成77777在线视频| 人人妻人人添人人爽欧美一区卜| 日韩欧美一区视频在线观看| 国产成人免费观看mmmm| 纯流量卡能插随身wifi吗| 黄片大片在线免费观看| 午夜激情av网站| 一夜夜www| 亚洲午夜理论影院| 国产在线视频一区二区| 久久午夜综合久久蜜桃| 精品一品国产午夜福利视频| 两人在一起打扑克的视频| 老熟妇仑乱视频hdxx| 久久精品亚洲av国产电影网| 国产精品1区2区在线观看. | 免费在线观看影片大全网站| 青草久久国产| 999精品在线视频| 在线观看免费日韩欧美大片| 757午夜福利合集在线观看| 国产免费现黄频在线看| 两个人看的免费小视频| 国产不卡一卡二| 99热网站在线观看| 免费在线观看影片大全网站| 视频在线观看一区二区三区| 国产精品久久久久成人av| 女人久久www免费人成看片| 国产日韩一区二区三区精品不卡| 日韩视频在线欧美| 99久久精品国产亚洲精品| 女人被躁到高潮嗷嗷叫费观| 国产黄色免费在线视频| 色婷婷久久久亚洲欧美| 精品久久久久久久毛片微露脸| 欧美久久黑人一区二区| 精品一品国产午夜福利视频| 首页视频小说图片口味搜索| 12—13女人毛片做爰片一| 99riav亚洲国产免费| 下体分泌物呈黄色| 天堂中文最新版在线下载| 午夜91福利影院| 五月开心婷婷网| 黄色a级毛片大全视频| 久久久欧美国产精品| 高清av免费在线| 丁香六月天网| 国产麻豆69| 亚洲 国产 在线| 国产欧美日韩一区二区精品| 亚洲专区国产一区二区| 中文字幕人妻丝袜制服| 男男h啪啪无遮挡| 欧美日本中文国产一区发布| 黑人欧美特级aaaaaa片| 欧美变态另类bdsm刘玥| 免费久久久久久久精品成人欧美视频| 午夜福利视频在线观看免费| 99久久人妻综合| 欧美人与性动交α欧美精品济南到| 黄网站色视频无遮挡免费观看| 一级片免费观看大全| 天堂俺去俺来也www色官网| 久9热在线精品视频| 嫩草影视91久久| 视频区图区小说| 人人妻,人人澡人人爽秒播| 性少妇av在线| 亚洲欧美色中文字幕在线| 免费av中文字幕在线| 美女高潮喷水抽搐中文字幕| 国产日韩一区二区三区精品不卡|