• <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国产精品久久久久影院| 久久精品久久久久久噜噜老黄| 黄色配什么色好看| 国产乱人偷精品视频| 欧美3d第一页| 久久久久久久精品精品| 少妇人妻久久综合中文| 国产精品国产三级专区第一集| 久久久久久久久久人人人人人人| 中文字幕人妻熟人妻熟丝袜美| 精品一区二区三区视频在线| 夫妻午夜视频| 91午夜精品亚洲一区二区三区| 国产欧美亚洲国产| 一级爰片在线观看| 老司机影院成人| 听说在线观看完整版免费高清| 91狼人影院| 久久久久久久久久久丰满| 乱系列少妇在线播放| 我要看日韩黄色一级片| 亚洲国产欧美在线一区| freevideosex欧美| 天堂中文最新版在线下载 | 精品国产一区二区三区久久久樱花 | 91久久精品国产一区二区三区| 丝袜喷水一区| 久久午夜福利片| 26uuu在线亚洲综合色| av在线app专区| 久久久久国产精品人妻一区二区| 成人无遮挡网站| 69av精品久久久久久| 我要看日韩黄色一级片| 狂野欧美激情性xxxx在线观看| 五月玫瑰六月丁香| 久久久久久久久大av| 91久久精品电影网| 国产成人freesex在线| 欧美97在线视频| videos熟女内射| 九九爱精品视频在线观看| 久久久久网色| videossex国产| 色播亚洲综合网| 免费观看性生交大片5| 视频中文字幕在线观看| 激情五月婷婷亚洲| 亚洲自偷自拍三级| 视频区图区小说| 视频区图区小说| 久久久久精品久久久久真实原创| 在线观看av片永久免费下载| 最近中文字幕2019免费版| 五月开心婷婷网| 不卡视频在线观看欧美| 国产伦在线观看视频一区| 成年女人在线观看亚洲视频 | 男女下面进入的视频免费午夜| 日韩欧美精品v在线| 欧美激情久久久久久爽电影| 又黄又爽又刺激的免费视频.| 少妇 在线观看| 国产一区亚洲一区在线观看| 免费电影在线观看免费观看| 久久99蜜桃精品久久| 天天躁夜夜躁狠狠久久av| 成人免费观看视频高清| 日韩一区二区视频免费看| av黄色大香蕉| 男人狂女人下面高潮的视频| 国产淫片久久久久久久久| 午夜亚洲福利在线播放| 男女国产视频网站| 一级片'在线观看视频| av免费在线看不卡| 久久99热这里只有精品18| 精品人妻视频免费看| 国产黄色免费在线视频| av天堂中文字幕网| 国产又色又爽无遮挡免| 在线观看一区二区三区激情| 亚洲在线观看片| 亚洲天堂国产精品一区在线| 免费黄网站久久成人精品| 亚洲图色成人| 久久精品国产鲁丝片午夜精品| 日本黄色片子视频| 午夜激情福利司机影院| 国产亚洲av嫩草精品影院| 亚洲精品影视一区二区三区av| 一个人观看的视频www高清免费观看| 一级二级三级毛片免费看| 97超视频在线观看视频| 亚洲精品色激情综合| 2021天堂中文幕一二区在线观| 久久99蜜桃精品久久| 免费观看a级毛片全部| 街头女战士在线观看网站| 午夜福利在线在线| 免费高清在线观看视频在线观看| 国产精品福利在线免费观看| 国产午夜精品久久久久久一区二区三区| 日产精品乱码卡一卡2卡三| 在线播放无遮挡| 欧美变态另类bdsm刘玥| 久久亚洲国产成人精品v| 少妇的逼水好多| 人妻 亚洲 视频| 国产精品人妻久久久久久| 亚洲精华国产精华液的使用体验| 久久久欧美国产精品| 少妇的逼水好多| 欧美老熟妇乱子伦牲交| 黄色配什么色好看| 三级国产精品欧美在线观看| 在线看a的网站| 精品熟女少妇av免费看| 麻豆久久精品国产亚洲av| 成人毛片a级毛片在线播放| 人妻夜夜爽99麻豆av| 麻豆久久精品国产亚洲av| 国产永久视频网站| 蜜臀久久99精品久久宅男| 国产精品偷伦视频观看了| 欧美日韩在线观看h| av播播在线观看一区| 91aial.com中文字幕在线观看| 97人妻精品一区二区三区麻豆| 亚洲欧美精品专区久久| 日韩伦理黄色片| 久久99热这里只有精品18| 亚洲综合色惰| 一本色道久久久久久精品综合| 日日撸夜夜添| 亚洲自拍偷在线| 亚洲欧美精品专区久久| 亚洲av成人精品一二三区| 日韩三级伦理在线观看| 男女边吃奶边做爰视频| 国产亚洲午夜精品一区二区久久 | 国产在视频线精品| 国产在视频线精品| 在现免费观看毛片| 大香蕉97超碰在线| 亚洲国产色片| 一级黄片播放器| 五月天丁香电影| 亚洲欧美中文字幕日韩二区| 午夜老司机福利剧场| 一区二区三区乱码不卡18| 免费大片18禁| 亚洲精品影视一区二区三区av| 一区二区三区免费毛片| 国精品久久久久久国模美| 一级毛片电影观看| 国产精品三级大全| 国产一区二区亚洲精品在线观看| 一级毛片 在线播放| 高清av免费在线| 国产男人的电影天堂91| 国产欧美另类精品又又久久亚洲欧美| 久久99热这里只有精品18| 夫妻午夜视频| 午夜亚洲福利在线播放| 一区二区三区乱码不卡18| 51国产日韩欧美| 男人狂女人下面高潮的视频| 免费黄频网站在线观看国产| videos熟女内射| 中文资源天堂在线| 免费看光身美女| 国产毛片a区久久久久| 国产午夜福利久久久久久| 国产爱豆传媒在线观看| 免费看a级黄色片| 日本午夜av视频| 色吧在线观看| 成人特级av手机在线观看| 嫩草影院新地址| 一区二区三区乱码不卡18| 国产精品一区www在线观看| 免费在线观看成人毛片| 午夜激情福利司机影院| 99热全是精品| 超碰97精品在线观看| 亚洲久久久久久中文字幕| 久久久久久久大尺度免费视频| 欧美高清性xxxxhd video| 亚洲人成网站高清观看| 美女被艹到高潮喷水动态| 在线播放无遮挡| 国产精品99久久久久久久久| 熟女av电影| 波多野结衣巨乳人妻| 91精品一卡2卡3卡4卡| 欧美日本视频| 99精国产麻豆久久婷婷| 久久热精品热| 少妇人妻 视频| 一本色道久久久久久精品综合| 欧美精品国产亚洲| av国产精品久久久久影院| 亚洲va在线va天堂va国产| 中文乱码字字幕精品一区二区三区| 欧美日韩国产mv在线观看视频 | 91aial.com中文字幕在线观看| eeuss影院久久| 国产成年人精品一区二区| 欧美性猛交╳xxx乱大交人| 18禁裸乳无遮挡免费网站照片| 51国产日韩欧美| 亚洲av在线观看美女高潮| 精品久久久久久久人妻蜜臀av| 国产精品一区二区三区四区免费观看| 国产午夜福利久久久久久| 在线观看一区二区三区| 毛片一级片免费看久久久久| 亚洲精品第二区| 狠狠精品人妻久久久久久综合| 男人添女人高潮全过程视频| 欧美激情国产日韩精品一区| 精品国产一区二区三区久久久樱花 | 一级毛片电影观看| 日日摸夜夜添夜夜爱| 99久国产av精品国产电影| 又爽又黄无遮挡网站| 成人美女网站在线观看视频| 搞女人的毛片| 亚洲电影在线观看av| 亚洲精品成人久久久久久| 久久久久久久午夜电影| 亚洲高清免费不卡视频| 亚洲欧美日韩另类电影网站 | 国产精品三级大全| 在线观看免费高清a一片| 熟女人妻精品中文字幕| 特大巨黑吊av在线直播| 精品视频人人做人人爽| 亚洲欧美日韩东京热| 国产成人精品久久久久久| 亚洲自拍偷在线| 国产精品人妻久久久影院| 一区二区三区四区激情视频| 男人爽女人下面视频在线观看| 97在线人人人人妻| 中文资源天堂在线| 国产女主播在线喷水免费视频网站| 日本午夜av视频| 97在线视频观看| 亚洲国产最新在线播放| 建设人人有责人人尽责人人享有的 | 白带黄色成豆腐渣| av在线蜜桃| 亚洲精品成人久久久久久| 丰满乱子伦码专区| 亚洲国产精品成人综合色| 人妻少妇偷人精品九色| 精品亚洲乱码少妇综合久久| 国产免费又黄又爽又色| 国产毛片a区久久久久| 午夜福利在线在线| 国产成人精品婷婷| 少妇裸体淫交视频免费看高清| 水蜜桃什么品种好| 国产成人a区在线观看| 亚洲国产日韩一区二区| 男人舔奶头视频| 国产国拍精品亚洲av在线观看| 亚洲aⅴ乱码一区二区在线播放| 下体分泌物呈黄色| 中文在线观看免费www的网站| 十八禁网站网址无遮挡 | 乱系列少妇在线播放| av黄色大香蕉| 一个人看视频在线观看www免费| 欧美精品人与动牲交sv欧美| 国产成人91sexporn| 男男h啪啪无遮挡| 91精品国产九色| 十八禁网站网址无遮挡 | 亚洲精品日本国产第一区| 伊人久久国产一区二区| 国产毛片a区久久久久| 国产午夜精品久久久久久一区二区三区| 亚洲精品国产成人久久av| 精品久久久精品久久久| 中文字幕免费在线视频6| 中国国产av一级| 久久精品综合一区二区三区| 国语对白做爰xxxⅹ性视频网站| 久热久热在线精品观看| 少妇的逼水好多| 午夜福利高清视频| 国产一区亚洲一区在线观看| 精品久久久精品久久久| 嘟嘟电影网在线观看| 亚洲性久久影院| av在线观看视频网站免费| 舔av片在线| 亚洲欧美中文字幕日韩二区| 久久久a久久爽久久v久久| 99精国产麻豆久久婷婷| av免费观看日本| 免费大片18禁| 最近的中文字幕免费完整| 联通29元200g的流量卡| 欧美日韩综合久久久久久| 在线观看一区二区三区| 91在线精品国自产拍蜜月| 性色avwww在线观看| 久久国内精品自在自线图片| 国产乱人偷精品视频| 欧美最新免费一区二区三区| 97人妻精品一区二区三区麻豆| 亚洲av电影在线观看一区二区三区 | 亚洲欧美一区二区三区国产| 97在线人人人人妻| 人人妻人人澡人人爽人人夜夜| 亚洲国产av新网站| 一区二区三区免费毛片| 久久精品国产鲁丝片午夜精品| 欧美日本视频| 成人高潮视频无遮挡免费网站| 免费av观看视频| 亚洲色图av天堂| 久久精品国产自在天天线| 26uuu在线亚洲综合色| 夜夜看夜夜爽夜夜摸| 精品久久久噜噜| 婷婷色av中文字幕| 极品教师在线视频| 中文字幕久久专区| 国产老妇女一区| 国产精品一及| 成人无遮挡网站| 丝袜喷水一区| 看十八女毛片水多多多| 在线观看一区二区三区| 女人十人毛片免费观看3o分钟| 色视频www国产| 麻豆成人av视频| 街头女战士在线观看网站| 午夜免费鲁丝| 热re99久久精品国产66热6| 欧美高清成人免费视频www| 婷婷色麻豆天堂久久| 在线观看美女被高潮喷水网站| 超碰97精品在线观看| 久久久久久久久久成人| 亚洲欧美日韩无卡精品| 校园人妻丝袜中文字幕| 性色av一级| 搞女人的毛片| 视频区图区小说| 亚洲欧美清纯卡通| 久久精品熟女亚洲av麻豆精品| 久久热精品热| 国产精品爽爽va在线观看网站| 麻豆乱淫一区二区| 久久久精品免费免费高清| 汤姆久久久久久久影院中文字幕| 99九九线精品视频在线观看视频| 香蕉精品网在线| 女人久久www免费人成看片| 搡女人真爽免费视频火全软件| 性插视频无遮挡在线免费观看| 高清毛片免费看| 国产男女超爽视频在线观看| 一级毛片 在线播放| 美女cb高潮喷水在线观看| 午夜精品一区二区三区免费看| 边亲边吃奶的免费视频| 亚洲精品国产色婷婷电影| 久久精品久久久久久久性| 亚洲在久久综合| 国产伦精品一区二区三区视频9| 春色校园在线视频观看| 免费人成在线观看视频色| 国产极品天堂在线| 午夜免费鲁丝| 国产乱人偷精品视频| 国内精品美女久久久久久| 综合色av麻豆| 超碰97精品在线观看| 一个人看的www免费观看视频| 亚洲国产av新网站| 亚洲精品国产av成人精品| 亚洲av成人精品一区久久| 国产成人精品久久久久久| 国产精品三级大全| 麻豆乱淫一区二区| 亚洲精品一二三| 美女主播在线视频| 国产精品久久久久久精品电影小说 | 欧美xxxx黑人xx丫x性爽| 老女人水多毛片| 三级国产精品欧美在线观看| 国产高清国产精品国产三级 | 亚洲人与动物交配视频| 欧美激情国产日韩精品一区| 亚洲怡红院男人天堂| 大香蕉97超碰在线| 少妇的逼水好多| av免费观看日本| 欧美日韩国产mv在线观看视频 | 成年av动漫网址| 天堂中文最新版在线下载 | 精品久久久精品久久久| 国产成人a区在线观看| 男人舔奶头视频| 一级毛片久久久久久久久女| 如何舔出高潮| 夫妻性生交免费视频一级片| 777米奇影视久久| 国产美女午夜福利| 欧美人与善性xxx| 日本午夜av视频| av卡一久久| 国产 一区 欧美 日韩| 91久久精品电影网| 青青草视频在线视频观看| 免费黄网站久久成人精品| 国产精品久久久久久久电影| 精品久久久久久久久av| 国产v大片淫在线免费观看| 国产白丝娇喘喷水9色精品| 精品久久久久久电影网| 黄色视频在线播放观看不卡| 国产精品不卡视频一区二区| 五月伊人婷婷丁香| 国模一区二区三区四区视频| 青春草视频在线免费观看| 一区二区三区精品91| 美女xxoo啪啪120秒动态图| 亚洲国产精品国产精品| 人妻夜夜爽99麻豆av| 校园人妻丝袜中文字幕| 欧美极品一区二区三区四区| 亚洲三级黄色毛片| 中文天堂在线官网| 成人漫画全彩无遮挡| 香蕉精品网在线| 日韩av在线免费看完整版不卡| 久久久久网色| 伦精品一区二区三区| 最新中文字幕久久久久| 国产在视频线精品| 成人一区二区视频在线观看| 国产精品成人在线| 在线观看免费高清a一片| 亚洲天堂国产精品一区在线| 交换朋友夫妻互换小说| 色视频在线一区二区三区| 尤物成人国产欧美一区二区三区| 国产永久视频网站| 青青草视频在线视频观看| 国产色爽女视频免费观看| 久久精品国产自在天天线| 熟女人妻精品中文字幕| 一个人看视频在线观看www免费| 精品一区二区免费观看| 国产日韩欧美在线精品| 国产日韩欧美在线精品| 2021天堂中文幕一二区在线观| 91久久精品国产一区二区三区| 亚洲精品乱码久久久v下载方式| 欧美日韩一区二区视频在线观看视频在线 | 国产免费一级a男人的天堂| 国产综合精华液| 国产伦精品一区二区三区四那| 久久久久精品久久久久真实原创| 色婷婷久久久亚洲欧美| 搞女人的毛片| 国产国拍精品亚洲av在线观看| 热re99久久精品国产66热6| av天堂中文字幕网| 中国三级夫妇交换| 成人鲁丝片一二三区免费| 成年女人看的毛片在线观看| 国产综合精华液| 99久久精品一区二区三区| 26uuu在线亚洲综合色| 一级黄片播放器| 亚洲成人av在线免费| 超碰av人人做人人爽久久| 高清在线视频一区二区三区| 日韩欧美精品v在线| 国产片特级美女逼逼视频| 人人妻人人澡人人爽人人夜夜| 一级a做视频免费观看| 日本三级黄在线观看| 波野结衣二区三区在线| 噜噜噜噜噜久久久久久91| av线在线观看网站| 国产亚洲5aaaaa淫片| 久久久色成人| 国产毛片a区久久久久| 丰满乱子伦码专区| 国产一区有黄有色的免费视频| 亚洲欧美日韩无卡精品| 亚洲国产色片| 亚洲色图av天堂| 欧美日韩视频高清一区二区三区二| 国产淫语在线视频| av在线观看视频网站免费| 免费观看的影片在线观看| 国产成人freesex在线| 激情五月婷婷亚洲| 中文在线观看免费www的网站| 久久99蜜桃精品久久| 男人和女人高潮做爰伦理| 一个人看的www免费观看视频| 视频中文字幕在线观看| 男人添女人高潮全过程视频| 又大又黄又爽视频免费| 国产人妻一区二区三区在| 国产高清不卡午夜福利| 综合色丁香网| 一级毛片aaaaaa免费看小| 精品一区二区免费观看| 亚洲精品456在线播放app| 亚洲精品乱码久久久久久按摩| 2021少妇久久久久久久久久久| 欧美最新免费一区二区三区| 黄色怎么调成土黄色| 国产成人午夜福利电影在线观看| 亚洲一级一片aⅴ在线观看| 国产久久久一区二区三区| 日本与韩国留学比较| 亚洲成人久久爱视频| 亚洲电影在线观看av| 亚洲aⅴ乱码一区二区在线播放| 麻豆成人av视频| 亚洲最大成人av| 亚洲三级黄色毛片| 日韩大片免费观看网站| 亚洲人成网站在线播| 国产黄频视频在线观看| 秋霞在线观看毛片| 日韩av在线免费看完整版不卡| 男人和女人高潮做爰伦理| 久久久久久久久久人人人人人人| 国精品久久久久久国模美| 国产毛片在线视频| 舔av片在线| 成人欧美大片| 亚洲av福利一区| av在线老鸭窝| 日日撸夜夜添| 尤物成人国产欧美一区二区三区| 3wmmmm亚洲av在线观看| 少妇熟女欧美另类| 国内少妇人妻偷人精品xxx网站| 免费观看的影片在线观看| 三级男女做爰猛烈吃奶摸视频| 欧美zozozo另类| 伊人久久精品亚洲午夜| 18禁在线播放成人免费| 精品一区二区三卡| 亚洲自偷自拍三级| 最近最新中文字幕大全电影3| 免费av毛片视频| 亚洲精品中文字幕在线视频 | 亚洲av福利一区| 国产精品麻豆人妻色哟哟久久| 亚洲国产精品专区欧美| 97在线人人人人妻| 亚洲美女视频黄频| 国产91av在线免费观看| 亚洲精品,欧美精品| 国产精品一区www在线观看| 亚洲av日韩在线播放| 国产亚洲av嫩草精品影院| 精品熟女少妇av免费看| 国模一区二区三区四区视频| 老司机影院毛片| 18禁在线无遮挡免费观看视频| 伦理电影大哥的女人| 啦啦啦中文免费视频观看日本| 久久久欧美国产精品| 日日撸夜夜添| 91精品一卡2卡3卡4卡| 成人特级av手机在线观看| 免费看不卡的av| 国产亚洲最大av| 看黄色毛片网站| 国国产精品蜜臀av免费| 日韩欧美一区视频在线观看 | 中文字幕久久专区| xxx大片免费视频| 在线观看人妻少妇| 亚洲av男天堂| tube8黄色片| 97超视频在线观看视频| 国产成人a区在线观看| 亚州av有码| 真实男女啪啪啪动态图| xxx大片免费视频| 亚洲av成人精品一区久久| 国产成人免费观看mmmm| 超碰97精品在线观看| 免费观看a级毛片全部| 极品教师在线视频| 九九久久精品国产亚洲av麻豆| 免费黄网站久久成人精品| 26uuu在线亚洲综合色| 国产91av在线免费观看| 成人国产麻豆网| 国产乱人偷精品视频| 日韩大片免费观看网站| 尤物成人国产欧美一区二区三区| 国产精品av视频在线免费观看| 九九爱精品视频在线观看|