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

    虹鱒(Oncorhynchus mykiss)魚體水動力學特性研究

    2022-07-07 02:53:04閆文彥黃六一李玉巖鄔騫力尤鑫星
    漁業(yè)現(xiàn)代化 2022年3期
    關鍵詞:虹鱒魚體流速

    閆文彥,黃六一,李玉巖,王 剛,鄔騫力,尤鑫星

    (中國海洋大學水產(chǎn)學院,青島 266003)

    魚類在游動過程中與流體產(chǎn)生非常復雜的相互作用[1]。魚類會借助側線等器官感覺周圍流體的變化,通過身體的自主擺動克服水流作用,以保持合適的魚體姿態(tài)和游泳速度[2-7]。在此過程中魚會消耗大量的能量,從而對其生存和生長產(chǎn)生重要影響[8-9]。研究魚體的水動力學特性,對了解魚類與水域環(huán)境的適應性、個體行為和集群特性、養(yǎng)殖技術和裝備開發(fā)以及仿生學具有重要的意義。

    目前國內(nèi)外學者對于魚體自身的肌肉運動、推進機理、游泳能力等方面開展較為豐富研究。Walters等[10]探究了鯖科魚類的身體形態(tài)和游泳性能;周仕杰等[11]對幾種幼魚曲線游泳能力進行了比較研究;Altringham等[12]研究了在模擬游泳條件下肌肉力量的體外輸出模式;Liao等[13]比較了虹鱒在不同直徑圓柱體后方的游泳狀態(tài)和在自由水流中的游泳狀態(tài),發(fā)現(xiàn)魚類可以從環(huán)境產(chǎn)生的漩渦中獲取能量,以維持其原有的位置;Liao等[14]從能量的角度解釋了采用卡門步態(tài)游泳的鱒魚的神經(jīng)肌肉控制;Lauder等[15-16]分析了魚類游泳的幾種波動推進的經(jīng)典模式;喬曉冬等[17]分析了魚類游動時反卡門渦街的形成和魚體推力的產(chǎn)生機理,結果表明,反卡門渦街的形成是魚體游動時產(chǎn)生推力的根本原因;周澍欣等[18]針對鯉魚的三維模型進行了數(shù)值模擬計算,對其水阻力情況、流場流速分布和壓力分布情況進行分析;Michele等[19]通過模擬肌肉功能來重現(xiàn)魚體游泳的主要特征,比較了不同游泳步態(tài)的性能,得到了頻率、波長、尾振幅等關鍵參數(shù)與所達到的游泳速度或產(chǎn)生的推力之間的關系,從而獲得其游泳性能;田曉潔等[20]針對仿金槍魚三維模型的流場受力進行對比研究,得到仿生機器魚可有效減少其水阻力,從而提高推進效率;王猛等[21]從亞科層面對中國鯉科魚類游泳能力進行了分析和總結;李明昊等[22]構建了仿生虎鯨直行和轉彎運動負載可靠性分析的狀態(tài)函數(shù),分析了仿生虎鯨的尾鰭設計變量對載荷的影響。之前大多數(shù)研究重點關注的是不同魚類的游泳模式與推進效率關系,為相關仿生學裝備設計提供服務。隨著野生漁業(yè)資源的衰退,人類通過養(yǎng)殖魚類,獲取高質(zhì)量蛋白質(zhì),已成為漁業(yè)的重要產(chǎn)業(yè)模式。養(yǎng)殖魚類自身的水動力特性與養(yǎng)殖水域的流速、能量消耗和養(yǎng)殖效率以及養(yǎng)殖技術和裝備設計密切相關,而目前針對魚體自身的水動力學特性的研究較少。

    虹鱒呈紡錘狀,體略側扁,因其肉質(zhì)蛋白含量高、膽固醇含量極低,具有很高的食用價值、保健價值及可期的市場前景[23-25],已成為中國水產(chǎn)養(yǎng)殖的重要品種之一,是黃海冷水團深遠海大型鋼結構網(wǎng)箱主要冷水魚養(yǎng)殖品種之一。

    本研究以虹鱒為研究對象,依據(jù)其魚體特征參數(shù)構建虹鱒三維模型,使用數(shù)值模擬方法探究虹鱒在不同流速和沖角下的水動力特性、流場流速分布以及魚體表面壓力分布,并通過水槽試驗驗證。研究結果可為虹鱒的個體和群體行為、養(yǎng)殖水域的選擇、養(yǎng)殖技術和裝備開發(fā)以及仿生學研究提供參考。

    1 材料與方法

    1.1 虹鱒模型

    預試驗發(fā)現(xiàn)小體長單個虹鱒模型在水流作用下的受力很小,儀器測量誤差對結果會產(chǎn)生較大影響。本研究選取大規(guī)格虹鱒(全長為84.22 cm,標準體長78.00 cm,體質(zhì)量6.80 kg),全長和體質(zhì)量參考《中國鮭鱒魚養(yǎng)殖》[26];體長、體高、頭長、體寬之間的關系參考李思忠等[27]對虹鱒的研究結果進行換算;除尾鰭外,其他部位魚鰭厚度較小,3D打印較難實現(xiàn),參照相關研究[20],忽略其他魚鰭。模型基本參數(shù)見表1。

    表1 虹鱒魚模型基本參數(shù)Tab.1 Basic parameters of rainbow trout model

    使用SolidWorks軟件構建1∶1的虹鱒魚體數(shù)值模型(圖1a),利用3D打印制作物理模型(圖1b)。3D打印材料選擇1.75 mm的聚乳酸(PLA),打印精度為層高0.1 mm,壁厚1.0 mm。

    圖1 虹鱒模型示意圖Fig.1 Schematic of rainbow trout model

    1.2 模型試驗

    1.2.1 儀器設備

    模型試驗在中國海洋大學循環(huán)動力水槽中進行,水槽觀察段尺寸為4.0 m × 1.0 m × 1.2 m;利用小威龍vectrino點式流速儀測定試驗流速,量程0~2 m/s,精度± 0.5%;穩(wěn)定流速范圍0.2~0.7 m/s;通過六分力儀傳感器測量試驗過程中魚體受到的水動力,量程0~50 kg,精度± 0.3%。

    1.2.2 試驗布設

    試驗布設如圖2所示,X方向沿水流方向,Y方向與水槽側壁垂直,Z軸垂直向上;通過兩根細螺紋桿(直徑d=10 mm)將虹鱒模型固定在六分力儀上,六分力儀處于觀察段中部,魚體中心線(吻端A尾鰭中心B的連線)距水面28 cm。魚體頭部的朝向-X方向。試驗水溫為20℃。魚體中心線所在的水平面為魚體中心平面,作為后續(xù)研究流場分布的參考面。

    圖2 試驗布設圖Fig.2 Experimental layout diagram

    1.3 數(shù)值模擬

    1.3.1 湍流模型

    對于復雜幾何的外部流動的問題,常采用基于雷諾平均N-S(RANS)方程組的k-ε模型。其控制方程為:

    (1)

    (2)

    式中:Ui表示i方向的濾波速度;uj表示j方向的濾波速度;P表示濾波壓力;ρ表示流體密度。

    本研究使用Realizablek-ε模型,可以在雷諾應力上保持與真實湍流的一致[28,29]。Realizablek-ε模型優(yōu)化了湍流粘度公式,從渦量擾動量均方根的精確輸運方程推導出ε方程。

    Realizablek-ε模型的湍動能及其耗散率輸運方程為:

    (3)

    (4)

    在上述方程中,Gk表示由于平均速度梯度引起的湍動能產(chǎn)生,Gb是由浮力產(chǎn)生的紊流動能;ym可壓速湍流脈動膨脹對總的耗散率的影響。C2和C1ε是常數(shù),σk,σε分別是湍動能及其耗散率的湍流普朗特數(shù),在 Fluent中,作為默認值常數(shù),C1ε=1.44,C2=1.9,σk=1.0,σε=1.2。

    1.3.2 計算域和邊界條件

    本研究使用的計算軟件為ANSYS 2021版中的Fluent模塊。計算域與試驗水槽觀察段尺度一致,設置為長度L=4.0 m、寬度W=1.2 m、高度H=1.0 m。虹鱒魚模型在計算域中的布設位置與模型試驗的布設相同。

    流體選為20℃下等溫不可壓縮水流,利用壓力-速度分離求解器(Semi-Implicit Method for Pressure Linked Equations-Consistent,SIMPLEC)開展穩(wěn)態(tài)計算。重力加速度為9.8 m/s2。計算域左側為速度入口邊界條件,右側為壓力出口邊界條件,其余為對稱壁面邊界條件(圖3),確保計算域側面邊界不會對魚體水動力分布產(chǎn)生無關干擾。N-S方程(式2)的對流項和擴散項在空間上的離散利用二階迎風格式。收斂標準設置為1×10-4,迭代計算步數(shù)為1 000步。

    圖3 計算域及邊界條件Fig.3 Calculation domain and boundary condition

    作為數(shù)值計算中的重要環(huán)節(jié),網(wǎng)格的質(zhì)量直接影響模擬結果的準確性。一般在數(shù)值仿真計算中采用的網(wǎng)格包括結構化網(wǎng)格和非結構化網(wǎng)格,非結構化網(wǎng)格具有網(wǎng)格生成速度快、質(zhì)量好、數(shù)據(jù)結構簡單等優(yōu)點。本研究采用的非結構化網(wǎng)格以六面體單元為主,利用Poly-Hexcore網(wǎng)格劃分方法(圖4),在魚體邊緣使用多面體網(wǎng)格與六面體網(wǎng)格實現(xiàn)共節(jié)點連接,以達到提升求解效率與精度的目的[30-33]??紤]到計算精度和計算資源的平衡,本研究劃分了5種網(wǎng)格進行網(wǎng)格無關性檢驗[34-36]。圖5表示0°沖角下0.4 m/s時魚體阻力系數(shù)隨著網(wǎng)格數(shù)的變化。由圖5可知,當網(wǎng)格數(shù)量達到314萬后,隨著網(wǎng)格數(shù)量的增加,阻力系數(shù)基本保持不變。綜合考慮目前計算機能力和計算精度要求,計算中保證網(wǎng)格數(shù)量大于314萬,根據(jù)不同的計算工況,實際使用的網(wǎng)格數(shù)量均處于400~500萬之間。

    圖4 網(wǎng)格示意圖Fig.4 Overview of the grids

    圖5 網(wǎng)格無關性驗證Fig.5 Mesh independence verification

    1.4 工況設置

    根據(jù)虹鱒的游泳速度、黃海冷水團常規(guī)流速以及水槽穩(wěn)定流速范圍,選取試驗流速為0.2~0.7 m/s,以0.1 m/s為間隔。沖角設定為來流方向與魚體中心線所在的縱向中剖面的水平夾角。將5條體長約19 cm的虹鱒放在水槽中試驗,從水槽頂部拍攝的圖像,獲取虹鱒在水流中投影線與來流方向的最大夾角為41.8°(圖6)。數(shù)值模擬時的沖角范圍據(jù)此設置為0°~45°范圍,以5°為間隔。

    圖6 虹鱒游動過程中魚體與水流的沖角Fig.6 The attack angle between the fish body andthe current during the swimming of therainbow trout

    在水槽中保持魚體沖角15°,測定不同來流速度下的魚體所受水阻力,與數(shù)值模擬結果進行對比驗證。在數(shù)值模擬驗證時,為方便對照,將水槽試驗固定模型的2根螺紋細桿一并設置在數(shù)值模型中,以減少誤差。在數(shù)值方法得到驗證基礎上,利用數(shù)值模擬方法開展所有工況的計算。

    1.5 數(shù)據(jù)處理方法

    魚體的阻力系數(shù)、側向力系數(shù)和雷諾數(shù)計算公式如下:

    (5)

    (6)

    (7)

    式中:Fx—水阻力,N;Fy—側向力,N;Cx—阻力系數(shù);Cy—側向力系數(shù);Re—雷諾數(shù);ρ—水介質(zhì)密度,取998.3 kg/m3;S—受流面積(魚體沿體長方向的中心截面),m2;U—入流速度,m/s;d—特征長度(體長),m;v—水的運動黏性系數(shù)(20℃水的運動黏性系數(shù)為1.006×10-6m2/s)。

    對監(jiān)測點速度無量綱處理,計算公式如下:

    (8)

    式中:K—速度比;Ui—速度場中任意點的速度;U—來流速度。

    2 結果分析

    2.1數(shù)值模擬與水槽試驗結果驗證

    圖7為15°沖角時虹鱒魚體水阻力實驗值與數(shù)值模擬結果的對比。由圖7可以看出,模型試驗和數(shù)值模擬得到的阻力曲線趨勢一致,均呈冪函數(shù)關系。在流速0.2~0.6 m/s時較為接近,流速0.4 m/s時相對誤差為2.56%,0.7 m/s時,相對誤差為14.02%。由此表明,本文采用的數(shù)值模擬方法具有較高的計算精度,能夠滿足后續(xù)的研究要求。

    圖7 15°沖角時虹鱒魚體水阻力實驗值 與數(shù)值模擬結果的對比Fig.7 Comparison of experimental and numerical resultson drag of rainbow trout at 15°

    2.2 魚體水動力特性

    2.2.1 水動力系數(shù)

    魚體的阻力系數(shù)、側向力系數(shù)與雷諾數(shù)關系如圖8和圖9所示。不同沖角條件下,阻力系數(shù)Cx和側向力系數(shù)Cy隨著雷諾數(shù)增大基本保持不變。表明雷諾數(shù)在1.6×105~5.8×105區(qū)間時,同一沖角條件下,水動力系數(shù)處于“自動模型區(qū)”。

    圖8 阻力系數(shù)Cx與雷諾數(shù)Re的變化關系Fig.8 Relationship between drag coefficient and Reynolds number

    圖9 側向力系數(shù)Cy與雷諾數(shù)Re的變化關系Fig.9 Relationship between side force coefficientand Reynolds number

    由于在試驗范圍內(nèi),在同一沖角條件下,阻力系數(shù)和側向力系數(shù)基本不變,采用阻力系數(shù)和側向力系數(shù)平均值得到圖10為阻力系數(shù)和側向力系數(shù)與沖角關系曲線。由圖可知,在0°~45°沖角范圍內(nèi),阻力系數(shù)Cx隨著沖角(α)增大呈現(xiàn)指數(shù)增長趨勢,其關系式為Cx=0.025 2 e0.073 1α(R2=0.994 6);側向力系數(shù)Cy隨著沖角(α)的增大而呈現(xiàn)線性增長,其關系式為Cy=0.015α(R2=0.977 8)。

    圖10 水動力系數(shù)與沖角的關系Fig.10 Relationship between hydrodynamic coefficientand inclinations

    2.2.2 魚體周圍速度場特性

    本研究取吻端與尾叉最深點連線所在的水平面為參考面(見圖2,魚體中心平面),探究該平面上的速度分布。將魚體周圍流速大于來流速度(無量綱流速大于1)的區(qū)域定義為高流速區(qū),低于來流流速(無量綱流速小于1)的區(qū)域定義為低流速區(qū)。圖11顯示了20°沖角時不同流速下魚體中心平面周圍無量綱流速分布。

    圖11 20°沖角時不同流速下魚體中心平面周圍流場圖Fig.11 Velocity fields around the fish center plane at different velocities at 20°

    由圖可知,當沖角固定時,不同流速下高流速區(qū)和低流速區(qū)分布趨勢基本一致,因此,下文以0.6 m/s為例說明由于沖角的變化而引起的速度場變化情況。

    不同沖角、流速0.6 m/s時,魚體中心平面周圍流場分布如圖12所示。由圖12可知,在0°沖角時,魚體左右兩側流速對稱分布,在魚體最寬處流速最大,形成兩側的高流速區(qū);頭部吻端和尾部附近流場流速較小,在頭部前端形成圓形的低流速區(qū),沿尾柄形成長條形低流速區(qū)。在沖角0°~15°時,迎流面高流速區(qū)出現(xiàn)于魚體最寬處到魚尾之間,說明沖角和魚體結構曲率對表面的流速分布具有協(xié)同影響。隨著沖角的逐漸增大,魚體背流面的高流速區(qū)逐漸向頭部前端移動,迎流面高流速區(qū)呈現(xiàn)逐漸縮小的趨勢,在0°沖角時,迎流面高流速區(qū)占魚體中心平面面積的2.30%,當沖角達到15°時,迎流面高流速區(qū)面積比例下降至1.14%。當沖角大于15°時,背流面高流速區(qū)形態(tài)在魚體頭部穩(wěn)定后,魚體尾鰭兩側開始逐漸形成影響范圍較小的高流速區(qū),隨著沖角的增大而增大,并沿著背流面逐漸從尾鰭尖端向軀干移動。這可能是由于邊界層在魚體頭部的背流面附近產(chǎn)生分離之后,魚體最寬處到魚尾之間曲率的變化導致了局部高流速區(qū)的產(chǎn)生。當沖角達30°時,迎流面高流速區(qū)消失。在沖角0°~45°時,迎流面低流速區(qū)隨著沖角增大而逐漸擴大,并由魚體前部吻端位置向后擴散,尾部迎流面低流速區(qū)隨著沖角增大而前移,最終低流速區(qū)覆蓋整個迎流面。周澍欣等[18]對鯉魚的三維模型進行了數(shù)值模擬探究,其流場結果與本研究一致,符合紡錘形魚體周圍流場的實際情況。

    圖12 不同沖角下魚體中心平面周圍流場分布圖Fig.12 Velocity fields around the fish center plane with varying inclinations

    2.2.3 魚體表面壓力分布

    流速0.6 m/s,不同沖角條件下,魚體表面壓力分布如圖13所示。0°沖角時,由于對稱于來流方向,魚體迎流面和背流面的壓力分布一致,負壓集中在體寬最大處。負壓是由于魚體的特定體型造成的結果,這與上節(jié)中流體流經(jīng)魚體體寬最大處時,兩側區(qū)域形成高流速區(qū)的結果相吻合,也與田曉潔等[20]對仿金槍魚的研究結果一致。

    圖13 不同沖角條件下魚體表面壓力分布云圖Fig.13 Surface pressure distribution nephogram offish with varying inclinations

    魚體迎流面正壓區(qū)在沖角小于20°時,主要集中在頭部和尾部,當沖角大于20°時,迎流面正壓區(qū)逐漸貫穿至整個軀干側面。隨著沖角的增大,迎流面上的負壓區(qū)域逐漸減少,負壓區(qū)逐漸向背、腹方向偏移。迎流面上的正壓區(qū)隨著沖角的增大逐漸增大,沖角0°時,正壓區(qū)占迎流面的23.90%,當沖角達到45°時,迎流面的71.66%均為正壓區(qū)。在沖角0°~20°時,魚體最寬處到魚尾之間由于其曲率的變化,壓力幅值逐漸減小,而上文中迎流面高流速區(qū)同樣出現(xiàn)于魚體最寬處到魚尾之間,這說明魚體表面的速度分布和壓力分布息息相關,符合伯努利原理。

    魚體背流面正壓區(qū)在沖角小于10°時,只存在于魚體頭部和尾柄部位,正壓區(qū)在魚體背流面的占比大于16.89%,當沖角大于10°時,頭部正壓區(qū)消失,尾柄部位正壓區(qū)隨著沖角的增大而減小,但在任意沖角條件下均為正壓,這也可以通過對照尾端速度的分布得到印證,當沖角達到45°時,背流面上正壓區(qū)僅有1.35%。魚體背流面負壓區(qū)在0°到15°沖角范圍內(nèi),隨著沖角的增加從最大體周處開始向頭部擴散,當沖角大于15°時,魚從魚體頭部開始,以魚體標準體長(L)為基準,提取魚體中心平面上的體表不同位置(0.1 L、0.25 L、0.5 L、0.75 L和0.9 L)的壓力數(shù)據(jù)(圖14和圖15),分析迎流面和背流面上魚體表面的壓力變化規(guī)律。由圖可知,隨著沖角的增大,迎流面表面壓力整體上均呈上升趨勢;從魚體位置角度分析,0.25 L之前(魚體頭部)以及0.9 L位置(尾柄),在所有沖角條件下,表面壓力整體上處于正壓力作用,其中,0.1 L處表面壓力均為最大,說明魚體頭部受到的表面壓力最大,最大壓力出現(xiàn)在沖角30°時,最大值約為194.22 Pa;從沖角的角度分析,沖角25°~45°,各個部位壓力均為正壓,除尾柄位置外,各位置壓力基本上隨著沖角的增加而增加,與圖13中的壓力分布云圖相一致。背流面表面壓力整體上隨著沖角的增加呈線性下降趨勢;沖角30°之后,各個部位表面壓力均為負壓,其中,頭部的壓力變化最大,由0°沖角的正壓102.31 Pa,下降至負壓-345.91 Pa。所有沖角條件下,0.9 L位置(尾柄)的壓力在0 Pa附近,變化不大,與圖13的壓力云圖相吻合。

    圖14 迎流面魚體不同部位的壓力變化Fig.14 The pressure in different parts of the fish bodychanges with varying inclinations in stoss side

    圖15 背流面魚體不同部位的壓力變化Fig.15 The pressure in different parts of the fish bodychanges with varying inclinations in slip slope

    3 結論

    本研究依據(jù)虹鱒魚體特征參數(shù),構建虹鱒魚三維模型,在數(shù)值模擬和水槽試驗驗證的基礎上,利用數(shù)值模擬方法探究了虹鱒個體在不同流況下的水動力特性、魚體中心平面周圍速度場分布以及魚體表面壓力分布情況。在試驗流速(0.2~0.7 m/s)范圍內(nèi),處于同一沖角條件下時,魚體水動力系數(shù)均處于自動模型區(qū)范圍內(nèi);在0°~45°沖角范圍內(nèi)阻力系數(shù)Cx隨著沖角增大從0.027增大至0.599,呈現(xiàn)指數(shù)增長趨勢,側向力系數(shù)Cy隨著沖角的增大而呈現(xiàn)線性增長,最大達到了0.71;由于虹鱒魚體本身的體型特征,隨著沖角的增大,速度場中高流速區(qū)和低流速區(qū)呈現(xiàn)一定的移動規(guī)律,隨著沖角的增大,魚體迎流面高流速區(qū)呈現(xiàn)逐漸后移并縮小的趨勢,背流面的高流速區(qū)逐漸向頭部前端移動,背流面高流速區(qū)形態(tài)在魚體頭部穩(wěn)定后,魚體尾鰭兩側開始逐漸形成影響范圍較小的高流速區(qū),隨著沖角的增大而增大,并沿著背流面逐漸從尾鰭尖端向軀干移動,說明沖角和魚體結構曲率對魚體周圍的流速分布具有協(xié)同影響;魚體結構與表面的壓力分布密切相關,對于魚體的迎流面而言,表面壓力在不同位置上整體呈上升趨勢,其中魚體頭部受壓最大(最大值可達194.22 Pa);對于背流面而言,背流面表面壓力整體上隨著沖角的增加呈線性下降趨勢,其中,頭部的壓力變化最大(102.31~-345.91 Pa)。以沖角的變化替代魚體的擺動,忽略魚鰭對魚體水動力學特性可能帶來的影響,此外魚體表面不同位置壓力分布和變化規(guī)律與魚體形狀密切相關,需要在今后開展更為深入的研究。

    猜你喜歡
    虹鱒魚體流速
    “流體壓強與流速的關系”知識鞏固
    『流體壓強與流速的關系』知識鞏固
    軸流泵內(nèi)魚體的運動行為與撞擊損傷分析
    淡水魚水平往復振動頭尾定向輸送方法
    山雨欲來風滿樓之流體壓強與流速
    大眾科學(2020年7期)2020-10-26 09:24:30
    淡水魚腹背定向裝置設計及試驗
    虹鱒人工繁育研究綜述
    愛虛張聲勢的水
    海水工廠化養(yǎng)殖虹鱒營養(yǎng)品質(zhì)分析與評價
    虹鱒魚的故事
    百科知識(2014年23期)2014-09-10 07:22:44
    国产三级黄色录像| av网站在线播放免费| 18禁裸乳无遮挡动漫免费视频| 久久精品国产99精品国产亚洲性色 | 国产精品影院久久| 成人免费观看视频高清| 久久亚洲真实| 亚洲avbb在线观看| 真人做人爱边吃奶动态| 麻豆成人av在线观看| 国产一区有黄有色的免费视频| 人人妻人人澡人人爽人人夜夜| 免费观看a级毛片全部| 久久久国产成人免费| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区精品视频观看| avwww免费| 久久人妻熟女aⅴ| 性少妇av在线| 一本大道久久a久久精品| 一二三四在线观看免费中文在| 婷婷成人精品国产| 国产一区二区激情短视频| 如日韩欧美国产精品一区二区三区| www.熟女人妻精品国产| 黄网站色视频无遮挡免费观看| 国产深夜福利视频在线观看| 叶爱在线成人免费视频播放| 夜夜夜夜夜久久久久| 一区二区三区国产精品乱码| 精品人妻1区二区| 午夜福利免费观看在线| 亚洲免费av在线视频| 亚洲成人手机| 91成人精品电影| 久久久国产精品麻豆| 日韩有码中文字幕| 黑人操中国人逼视频| 嫩草影视91久久| 亚洲精品国产精品久久久不卡| 欧美成狂野欧美在线观看| 亚洲精品国产色婷婷电影| 老司机在亚洲福利影院| 18禁国产床啪视频网站| 欧美精品一区二区免费开放| 国产精品电影一区二区三区 | 欧美丝袜亚洲另类 | 精品高清国产在线一区| 涩涩av久久男人的天堂| 满18在线观看网站| 亚洲精品国产一区二区精华液| 中文字幕精品免费在线观看视频| 天天影视国产精品| 1024视频免费在线观看| 50天的宝宝边吃奶边哭怎么回事| 中文字幕制服av| 国产一区二区激情短视频| av福利片在线| 精品久久久精品久久久| 午夜福利视频在线观看免费| 中文字幕人妻丝袜制服| 777米奇影视久久| 99国产精品一区二区三区| 天天影视国产精品| 亚洲国产欧美一区二区综合| 精品第一国产精品| 亚洲精品中文字幕在线视频| 免费高清在线观看日韩| 18禁裸乳无遮挡动漫免费视频| 又黄又粗又硬又大视频| 精品人妻在线不人妻| 母亲3免费完整高清在线观看| 国产成人啪精品午夜网站| 手机成人av网站| 久久亚洲精品不卡| 又大又爽又粗| 国精品久久久久久国模美| 一进一出抽搐gif免费好疼 | 久久久久久亚洲精品国产蜜桃av| 高清欧美精品videossex| 亚洲精品在线美女| 91精品三级在线观看| 免费观看a级毛片全部| 日日爽夜夜爽网站| 免费黄频网站在线观看国产| 母亲3免费完整高清在线观看| 成人18禁在线播放| 老熟妇仑乱视频hdxx| 久久香蕉国产精品| 18禁裸乳无遮挡动漫免费视频| 在线十欧美十亚洲十日本专区| 午夜福利欧美成人| 亚洲五月婷婷丁香| 久久久久精品国产欧美久久久| 一边摸一边做爽爽视频免费| 岛国毛片在线播放| 狂野欧美激情性xxxx| 久久久久国产精品人妻aⅴ院 | 人妻丰满熟妇av一区二区三区 | 日韩人妻精品一区2区三区| 高清在线国产一区| 18禁观看日本| 日韩视频一区二区在线观看| 法律面前人人平等表现在哪些方面| 女人爽到高潮嗷嗷叫在线视频| 在线视频色国产色| 日韩大码丰满熟妇| 中文字幕高清在线视频| 在线视频色国产色| 国产精华一区二区三区| 亚洲av成人一区二区三| av视频免费观看在线观看| 这个男人来自地球电影免费观看| 久久久久久久久免费视频了| av天堂在线播放| 亚洲人成77777在线视频| 成人国语在线视频| 国产亚洲精品久久久久久毛片 | 好看av亚洲va欧美ⅴa在| 国产一区有黄有色的免费视频| 久久亚洲真实| 久久国产亚洲av麻豆专区| 乱人伦中国视频| 天堂中文最新版在线下载| 亚洲avbb在线观看| 亚洲中文字幕日韩| 99国产综合亚洲精品| 99久久综合精品五月天人人| 妹子高潮喷水视频| 精品久久久久久久毛片微露脸| 久久久国产成人精品二区 | 国产av一区二区精品久久| 最新的欧美精品一区二区| 一级毛片女人18水好多| 搡老岳熟女国产| 日韩人妻精品一区2区三区| 亚洲在线自拍视频| 99久久人妻综合| 51午夜福利影视在线观看| 丰满人妻熟妇乱又伦精品不卡| 在线观看www视频免费| 日韩熟女老妇一区二区性免费视频| 首页视频小说图片口味搜索| av免费在线观看网站| 国产一区在线观看成人免费| 亚洲av日韩精品久久久久久密| 最新的欧美精品一区二区| 一级作爱视频免费观看| av天堂久久9| 午夜激情av网站| 国产精品久久久av美女十八| 捣出白浆h1v1| 视频在线观看一区二区三区| 精品午夜福利视频在线观看一区| 99国产精品99久久久久| 成人特级黄色片久久久久久久| 亚洲国产看品久久| 中文字幕另类日韩欧美亚洲嫩草| 女警被强在线播放| 怎么达到女性高潮| 亚洲欧美激情在线| 亚洲精品久久成人aⅴ小说| 狠狠婷婷综合久久久久久88av| 国产欧美日韩综合在线一区二区| 国产激情久久老熟女| 麻豆成人av在线观看| 久久香蕉国产精品| 黄色丝袜av网址大全| 欧美日韩亚洲高清精品| 亚洲七黄色美女视频| 中文字幕最新亚洲高清| 欧美av亚洲av综合av国产av| 午夜精品在线福利| 欧美精品啪啪一区二区三区| 国产午夜精品久久久久久| 国产97色在线日韩免费| 悠悠久久av| 国产亚洲精品久久久久5区| 中文字幕最新亚洲高清| 久久人人97超碰香蕉20202| 亚洲av第一区精品v没综合| xxx96com| 亚洲av成人不卡在线观看播放网| a在线观看视频网站| 看黄色毛片网站| 在线观看免费视频网站a站| 后天国语完整版免费观看| 咕卡用的链子| 欧美日本中文国产一区发布| 欧美中文综合在线视频| 日韩制服丝袜自拍偷拍| 嫁个100分男人电影在线观看| 丝袜美足系列| 精品国产美女av久久久久小说| 国产精品美女特级片免费视频播放器 | 很黄的视频免费| 九色亚洲精品在线播放| 国产精品.久久久| 大片电影免费在线观看免费| 亚洲av成人不卡在线观看播放网| 久久性视频一级片| 曰老女人黄片| 人人妻人人澡人人爽人人夜夜| 午夜成年电影在线免费观看| 免费黄频网站在线观看国产| 免费在线观看黄色视频的| 又黄又爽又免费观看的视频| 亚洲av电影在线进入| svipshipincom国产片| 伊人久久大香线蕉亚洲五| 久久亚洲真实| 欧美人与性动交α欧美软件| 99re6热这里在线精品视频| 首页视频小说图片口味搜索| 男人舔女人的私密视频| 黄片大片在线免费观看| 成年人午夜在线观看视频| 黄色视频不卡| 成人18禁高潮啪啪吃奶动态图| 人人妻,人人澡人人爽秒播| 在线观看午夜福利视频| 不卡av一区二区三区| 99精品欧美一区二区三区四区| 久久久国产一区二区| 婷婷成人精品国产| 天天操日日干夜夜撸| 久久青草综合色| 啦啦啦免费观看视频1| 欧美色视频一区免费| www日本在线高清视频| 黄色视频不卡| 黄片播放在线免费| 久久精品国产亚洲av香蕉五月 | avwww免费| 一进一出抽搐动态| 一二三四在线观看免费中文在| 人人妻人人澡人人爽人人夜夜| 成人影院久久| 无遮挡黄片免费观看| 久久性视频一级片| 又紧又爽又黄一区二区| 91成年电影在线观看| 久久ye,这里只有精品| 婷婷成人精品国产| 天天影视国产精品| 男女床上黄色一级片免费看| 人妻丰满熟妇av一区二区三区 | 久久久精品区二区三区| 国产免费av片在线观看野外av| 91av网站免费观看| 最新在线观看一区二区三区| 女性被躁到高潮视频| 久热这里只有精品99| 免费人成视频x8x8入口观看| www.精华液| 成人特级黄色片久久久久久久| 天堂√8在线中文| 一进一出抽搐动态| 又黄又爽又免费观看的视频| 国产精品1区2区在线观看. | 婷婷丁香在线五月| 国产1区2区3区精品| 国产精品香港三级国产av潘金莲| 超色免费av| 精品国产乱子伦一区二区三区| 久久婷婷成人综合色麻豆| 19禁男女啪啪无遮挡网站| 久久久国产一区二区| 中文字幕制服av| 一级片免费观看大全| 免费在线观看日本一区| 在线观看一区二区三区激情| 亚洲熟女精品中文字幕| 热99国产精品久久久久久7| 午夜福利一区二区在线看| а√天堂www在线а√下载 | 黄片小视频在线播放| 欧美日韩一级在线毛片| 精品国产国语对白av| 亚洲第一青青草原| 国产亚洲欧美98| 国产乱人伦免费视频| 国产av一区二区精品久久| 人人妻人人澡人人看| 少妇 在线观看| 女性生殖器流出的白浆| 老司机午夜十八禁免费视频| www日本在线高清视频| 精品一品国产午夜福利视频| 女警被强在线播放| 水蜜桃什么品种好| 午夜福利在线免费观看网站| 国产激情久久老熟女| 999精品在线视频| 18禁美女被吸乳视频| 99精品在免费线老司机午夜| 亚洲成人免费av在线播放| 亚洲aⅴ乱码一区二区在线播放 | 久久中文字幕一级| 91老司机精品| 啦啦啦免费观看视频1| 欧美日韩福利视频一区二区| 午夜91福利影院| 亚洲第一av免费看| 亚洲精华国产精华精| 精品国产国语对白av| 在线观看一区二区三区激情| 在线观看免费日韩欧美大片| 又黄又爽又免费观看的视频| 一本一本久久a久久精品综合妖精| av福利片在线| 女人爽到高潮嗷嗷叫在线视频| 国产xxxxx性猛交| 久久香蕉激情| 天天躁日日躁夜夜躁夜夜| 亚洲久久久国产精品| 女人被狂操c到高潮| 欧美日韩亚洲高清精品| 亚洲精华国产精华精| 咕卡用的链子| 天堂√8在线中文| √禁漫天堂资源中文www| av中文乱码字幕在线| 久久久久国内视频| 夜夜夜夜夜久久久久| 国产精品久久久久成人av| 中文字幕人妻丝袜一区二区| 日韩一卡2卡3卡4卡2021年| 欧美不卡视频在线免费观看 | 久久久久久久国产电影| 满18在线观看网站| 一级毛片精品| 国产一区有黄有色的免费视频| 欧美亚洲 丝袜 人妻 在线| 91在线观看av| 人妻 亚洲 视频| 国产精品秋霞免费鲁丝片| av免费在线观看网站| 国产精品久久久人人做人人爽| 99久久国产精品久久久| 欧美日韩av久久| 国产精品.久久久| 多毛熟女@视频| 99国产极品粉嫩在线观看| 麻豆av在线久日| 国产成人精品久久二区二区免费| 在线永久观看黄色视频| 国产成人影院久久av| 免费在线观看视频国产中文字幕亚洲| 男女免费视频国产| 国产乱人伦免费视频| 交换朋友夫妻互换小说| 国产精品免费视频内射| 欧美亚洲 丝袜 人妻 在线| 黑人巨大精品欧美一区二区mp4| 日韩人妻精品一区2区三区| 激情在线观看视频在线高清 | 精品国产亚洲在线| 国产av又大| 热re99久久精品国产66热6| 高清毛片免费观看视频网站 | 亚洲伊人色综图| 午夜老司机福利片| 精品人妻1区二区| 亚洲欧美激情综合另类| 久久国产乱子伦精品免费另类| 老司机午夜十八禁免费视频| 波多野结衣av一区二区av| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av第一区精品v没综合| 欧美 亚洲 国产 日韩一| 香蕉丝袜av| 看黄色毛片网站| 露出奶头的视频| netflix在线观看网站| 亚洲人成77777在线视频| netflix在线观看网站| 黄色视频,在线免费观看| 国产成人精品久久二区二区免费| 亚洲九九香蕉| 又黄又爽又免费观看的视频| 1024视频免费在线观看| 热re99久久精品国产66热6| 如日韩欧美国产精品一区二区三区| 国产99白浆流出| 欧美色视频一区免费| 亚洲欧美激情在线| 久久久国产成人免费| 中文字幕制服av| 99久久国产精品久久久| 亚洲成人免费av在线播放| 国产成人欧美在线观看 | 亚洲欧洲精品一区二区精品久久久| 国产精品久久久人人做人人爽| 中文字幕最新亚洲高清| 老汉色av国产亚洲站长工具| 国产成人免费无遮挡视频| 一个人免费在线观看的高清视频| 午夜日韩欧美国产| 国产精品av久久久久免费| 欧美人与性动交α欧美软件| av网站免费在线观看视频| 国产乱人伦免费视频| 1024视频免费在线观看| 婷婷成人精品国产| 亚洲人成伊人成综合网2020| 免费黄频网站在线观看国产| 午夜久久久在线观看| 日韩欧美一区视频在线观看| 制服人妻中文乱码| 久久香蕉国产精品| 亚洲专区字幕在线| 亚洲精华国产精华精| 午夜两性在线视频| 欧美日韩福利视频一区二区| 国产日韩欧美亚洲二区| 熟女少妇亚洲综合色aaa.| 一区二区日韩欧美中文字幕| 久久热在线av| 亚洲精品国产精品久久久不卡| 看黄色毛片网站| 日韩欧美免费精品| 欧美日韩精品网址| 久久午夜亚洲精品久久| 91字幕亚洲| 亚洲成a人片在线一区二区| 免费黄频网站在线观看国产| 免费av中文字幕在线| 国产精品免费一区二区三区在线 | 国产深夜福利视频在线观看| 亚洲一区高清亚洲精品| 在线观看免费视频网站a站| av天堂久久9| 国产一区二区三区综合在线观看| 丝袜美足系列| 大型av网站在线播放| 黑人巨大精品欧美一区二区mp4| 久久九九热精品免费| 亚洲国产欧美日韩在线播放| 国产成人免费观看mmmm| avwww免费| 最新美女视频免费是黄的| 成年人免费黄色播放视频| 91在线观看av| 人成视频在线观看免费观看| 在线观看www视频免费| 国产极品粉嫩免费观看在线| 在线视频色国产色| 欧美精品av麻豆av| 亚洲中文av在线| 日韩免费高清中文字幕av| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品国产精品久久久不卡| 亚洲午夜理论影院| 99国产精品一区二区蜜桃av | 久久这里只有精品19| 久久久国产欧美日韩av| 成人三级做爰电影| 黄片小视频在线播放| 久久 成人 亚洲| 免费观看人在逋| 黄色成人免费大全| 操出白浆在线播放| 久久国产精品男人的天堂亚洲| 正在播放国产对白刺激| 女人久久www免费人成看片| 国产在视频线精品| 又大又爽又粗| 日韩制服丝袜自拍偷拍| 欧美日韩亚洲国产一区二区在线观看 | 日韩视频一区二区在线观看| 成年版毛片免费区| 久久久国产成人免费| 91成人精品电影| 成人18禁在线播放| 美女高潮到喷水免费观看| 欧美人与性动交α欧美软件| 黄片大片在线免费观看| av天堂在线播放| 午夜久久久在线观看| 十八禁网站免费在线| 天天躁日日躁夜夜躁夜夜| av网站免费在线观看视频| 婷婷成人精品国产| 亚洲精品乱久久久久久| 欧美激情 高清一区二区三区| 一区二区三区国产精品乱码| 人成视频在线观看免费观看| 每晚都被弄得嗷嗷叫到高潮| 国产单亲对白刺激| 黄片小视频在线播放| 午夜日韩欧美国产| 亚洲五月天丁香| 老司机靠b影院| 国产麻豆69| 人妻 亚洲 视频| 日日爽夜夜爽网站| 男女下面插进去视频免费观看| 久久久久久免费高清国产稀缺| 国产麻豆69| 一区二区日韩欧美中文字幕| 国产1区2区3区精品| 久久99一区二区三区| 久久久精品区二区三区| 欧美丝袜亚洲另类 | 精品一区二区三卡| 精品少妇久久久久久888优播| 国产精品av久久久久免费| 亚洲五月色婷婷综合| 热99re8久久精品国产| 亚洲国产看品久久| 啦啦啦在线免费观看视频4| 免费高清在线观看日韩| 女人爽到高潮嗷嗷叫在线视频| 极品教师在线免费播放| 日本撒尿小便嘘嘘汇集6| 国产激情欧美一区二区| 亚洲视频免费观看视频| 亚洲精品国产一区二区精华液| 在线观看一区二区三区激情| 欧美一级毛片孕妇| 国产野战对白在线观看| 99热网站在线观看| 精品亚洲成a人片在线观看| 性少妇av在线| 极品少妇高潮喷水抽搐| 亚洲成a人片在线一区二区| 国产亚洲精品第一综合不卡| 国产区一区二久久| 在线视频色国产色| 久久久国产一区二区| 日韩欧美一区视频在线观看| 欧美国产精品一级二级三级| av视频免费观看在线观看| 欧美大码av| 老汉色av国产亚洲站长工具| 看片在线看免费视频| 午夜福利在线免费观看网站| 一本综合久久免费| 日本a在线网址| 老司机在亚洲福利影院| 麻豆av在线久日| 黄网站色视频无遮挡免费观看| 久久九九热精品免费| 欧美 日韩 精品 国产| 黄色片一级片一级黄色片| 精品视频人人做人人爽| 久久国产亚洲av麻豆专区| 电影成人av| 久久精品国产a三级三级三级| 国产高清videossex| 精品国产一区二区三区久久久樱花| 亚洲五月色婷婷综合| 国产区一区二久久| 91精品三级在线观看| 亚洲一区中文字幕在线| 午夜老司机福利片| 欧美老熟妇乱子伦牲交| 一本大道久久a久久精品| 成人黄色视频免费在线看| 午夜福利欧美成人| 男男h啪啪无遮挡| 中亚洲国语对白在线视频| 午夜免费成人在线视频| 成人18禁高潮啪啪吃奶动态图| 91av网站免费观看| 每晚都被弄得嗷嗷叫到高潮| 午夜成年电影在线免费观看| 成在线人永久免费视频| 亚洲中文字幕日韩| 精品乱码久久久久久99久播| netflix在线观看网站| 亚洲五月婷婷丁香| 国产又爽黄色视频| 看片在线看免费视频| 最新的欧美精品一区二区| 波多野结衣av一区二区av| 免费在线观看日本一区| 国产精品 国内视频| 国产日韩欧美亚洲二区| 一级毛片女人18水好多| 国产精品 国内视频| 日本五十路高清| 欧美国产精品一级二级三级| 91精品三级在线观看| 91在线观看av| 免费人成视频x8x8入口观看| 久久久精品区二区三区| 人人妻人人澡人人看| 99re在线观看精品视频| 制服人妻中文乱码| 日韩熟女老妇一区二区性免费视频| 免费在线观看视频国产中文字幕亚洲| 欧美精品人与动牲交sv欧美| 在线观看www视频免费| 老司机午夜十八禁免费视频| 欧美精品高潮呻吟av久久| 午夜久久久在线观看| 亚洲av成人av| 欧美 亚洲 国产 日韩一| 男人的好看免费观看在线视频 | 成人18禁高潮啪啪吃奶动态图| 久久久久久久午夜电影 | 中文亚洲av片在线观看爽 | 飞空精品影院首页| 日本vs欧美在线观看视频| 久久狼人影院| 精品一区二区三区av网在线观看| 久久久国产精品麻豆| 老熟妇乱子伦视频在线观看| 两性夫妻黄色片| 亚洲少妇的诱惑av| 免费av中文字幕在线| 在线天堂中文资源库| 欧美 亚洲 国产 日韩一|