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

    螺旋槳梢渦及梢渦空泡數(shù)值模擬

    2019-05-14 01:36:40劉芳遠傅慧萍
    船舶力學 2019年4期
    關鍵詞:空泡空化螺旋槳

    劉芳遠,傅慧萍,李 杰

    (1.上海交通大學 船舶海洋與建筑工程學院,上海200240;2.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240)

    0 引 言

    螺旋槳空泡是一種汽化空泡,即水因降壓到飽和蒸汽壓力導致汽化,水汽通過界面,進入氣核并使之膨脹。螺旋槳槳葉上出現(xiàn)的空泡有渦空泡、泡狀空泡、片空泡和云霧空泡4 種類型。螺旋槳空泡會對螺旋槳的水動力性能和噪聲性能產(chǎn)生重要影響。準確的數(shù)值模擬方法可以對螺旋槳性能進行預報,并為螺旋槳的設計提供理論依據(jù)。

    近廿年來,國內(nèi)外對空化流提出的數(shù)值模擬方法有很多,大致可以分為勢流方法和Navier-Stokes方程求解及歐拉—拉格朗日混合方法三種。盡管作為一個水動力學熱點問題已經(jīng)被研究多年,但囿于目前的數(shù)值模擬準確性,仍然有眾多學者在空化模型及與空化相匹配的湍流模型,以及數(shù)值計算方法等方面開展相關的研究。尤其是對于螺旋槳梢渦及梢渦空泡的模擬,目前仍然是一個難點。德國波茲坦水池(SVA)曾于2011 年在漢堡主辦了第二屆船舶推進器國際研討會(SMP’11),該研討會采用盲算的方式進行比較分析,其中發(fā)布的第3 個空化案例(Case 2.3.3)在物理模型試驗中出現(xiàn)了顯而易見的梢渦空泡,數(shù)十家與會單位竟無一能算出。

    國內(nèi),韓寶玉等(2011)[1]針對梢渦空化問題,應用雷諾平均的N-S 方程并結合改進的VOF(volume fraction)模型,對橢圓水翼梢渦空化特性進行了研究,采用的湍流模型為代數(shù)雷諾應力模型。為了準確模擬梢渦渦核內(nèi)系統(tǒng)旋轉和流線曲率影響,將旋轉和曲率修正方法應用于湍流模型中。但是對水翼的水動力沒有進行定量的校核,只是給出了所謂消失空泡數(shù)。Shi 等(2015)[2]應用分離渦模型(DES)對螺旋槳尾流中梢渦流場分布進行了數(shù)值研究,但還沒有算過空化。Ji 等(2013)[3]采用大渦模擬(LES)和基于質(zhì)量輸運方程的空化模型預測了一個3 維扭曲水翼的非定??栈?。辛公正 (2014)[4]采用數(shù)值RANS 方法模擬了水翼和螺旋槳的梢渦流動,并不求解空泡,而是尋找最小壓力系數(shù)與起始空泡數(shù)之間的關系。Liu 等(2016)[5]基于OpenFOAM 平臺,比較了線性與非線性湍流模型對螺旋槳梢渦空化的數(shù)值模擬能力,發(fā)現(xiàn)非線性k-ε 模型能夠更加準確地捕捉到梢渦空化現(xiàn)象。但沒有同時對螺旋槳的推力系數(shù)和扭矩系數(shù)進行校核。將螺旋槳空泡水動力性能計算做得比較好的有馮學梅等(2012)[6],他們基于FLUENT 軟件,采用多塊結構化網(wǎng)格對業(yè)內(nèi)選作標槳的E779A 槳和PPTC 槳進行了均勻來流下的空泡數(shù)值模擬。

    目前主流的空化計算方法是結合空化模型,求解兩相流的RANS 方程。歐拉—拉格朗日混合方法主要用于研究空化初生[7]。Hsiao 等(2005)[8]以水翼為研究對象,采用特定區(qū)域發(fā)放氣核的研究方法,考慮氣核尺寸分布進行梢渦空化初生噪聲尺度效應研究。真正與螺旋槳相關是在2008 年第27 屆海軍水動力學會議上,他們發(fā)表了敞水螺旋槳的梢渦空化初生尺度效應研究[9]。最近的研究成果是關于螺旋槳及氣體擴散作用對氣核分布的影響[10]。國內(nèi)熊鷹等(2013)[11]也采用這種方法對螺旋槳空化初生進行尺度效應研究,但均不直接模擬空泡。

    本文基于螺旋槳空泡數(shù)值模擬的現(xiàn)狀,以PPTC 槳為研究對象,通過一種新穎的網(wǎng)格劃分方式,求解均勻來流下的螺旋槳空泡及其水動力。目標是同時求解出與試驗相符的空泡形狀,以及螺旋槳的推力和扭矩系數(shù)。

    1 控制方程

    采用多相流模型中的均質(zhì)混合流模型,此時氣液兩相流被當作包含微氣泡的單相連續(xù)介質(zhì),通過引進氣相體積分數(shù)α 來定義混合相的密度ρm,僅求解一套控制方程?;诨旌厦芏鹊木|(zhì)混合流的連續(xù)方程為

    動量方程為

    式中:下標l、v 分別指液相和氣相;μm為混合流的動力粘性系數(shù),定義與混合密度類似;μt,m為湍流引起的渦粘性系數(shù),需要對混合流采用湍流模型進行求解。本文采用RNG k-ε 兩方程模型求解渦粘系數(shù),該模型在求解旋轉流動和剪切流時有優(yōu)勢,壁面第一層網(wǎng)格節(jié)點的無量綱距離y+要求在30~100 之間。

    關于蒸汽質(zhì)量分數(shù)f 的控制方程為

    式中:Re和Rc分別為由汽化和凝結引起的相變率;f 按下式定義:

    根據(jù)對(5)式右端源項的建模方式,可以區(qū)分出不同的空化模型。常用的有Zwart-Gerber-Belamri模型和Schnerr-Sauer 模型。本文采用Zwart 模型:

    當p≤pv時,

    當p>pv時,

    式中:pv為水的飽和蒸氣壓,與溫度相關;氣核半徑RB默認取值1 μm;αnuc為水中所含氣核的體積分數(shù),默認取值5×10-4;Ce和Cc分別為汽化系數(shù)和凝結系數(shù),默認取值50 和0.01。

    2 研究對象與網(wǎng)格

    PPTC 槳是指德國波茲坦水池(SVA)設計的可調(diào)螺距槳VP1304,也是第二屆和第四屆船舶推進器國際研討會(SMP’11 和SMP’15)的算例槳,其在均流和斜流情形下的敞水性能、伴流分布、空化性能及壓力脈動等試驗數(shù)據(jù)通過互聯(lián)網(wǎng)為全世界的水動力學研究人員所共享。本文即以SMP’11 的算例槳PPTC’11 為研究對象,對研討會上發(fā)布的空化案例Case 2.3 進行數(shù)值模擬,以校驗空化模型及相應的空化流場計算方法。

    圖1 研究對象與計算域Fig.1 Computed object and computational field

    PPTC 為直徑D=250 mm 的5 葉右旋可調(diào)螺距槳,這使得槳葉在導邊和隨邊處與槳轂之間存在0.3 mm 的間隙,數(shù)值模擬時予以忽略。SMP’11 的空泡測試在SVA 的空泡筒K15A 中完成,軸無傾角。分別在15、20、25 r/s 三個轉速下對敞水性能進行了測試,以探討雷諾數(shù)對敞水性能的影響;敞水試驗在高壓時進行,以避免空泡。圖1 為SVA 提供的PPTC 槳模與空泡筒試驗段幾何模型。為了更好地模擬均勻定常來流下槳的水動力性能,按SMP’11 的要求,計算時將試驗段變換成等截面的圓柱段,其它如槳模安裝位置和試驗段長度保持不變。

    依據(jù)流動最大雷諾數(shù)及湍流模型的適用范圍估算了壁面第一層網(wǎng)格節(jié)點的距離,約為2.4D‰,據(jù)此在近壁區(qū)進行邊界層網(wǎng)格加密。尤其是對槳葉梢部的梢渦脫出區(qū)域進行了區(qū)域劃分及域內(nèi)網(wǎng)格加密,以保證渦核區(qū)域內(nèi)徑向15 個網(wǎng)格節(jié)點的密度(如圖2 所示)。生成全流道計算域網(wǎng)格,共8 560 925 個混合型單元。

    圖2 梢部網(wǎng)格加密Fig.2 Local grid refinement of the tip vicinity

    3 無空化流場計算

    PPTC 槳被設計專門用于產(chǎn)生梢渦,并展開廣泛的模型試驗,用以校驗螺旋槳的各種勢流或粘流分析程序。在無空化情形下,SVA 在空泡筒對3 個轉速下的PPTC 敞水性能進行了測試。本章選擇3條敞水性能曲線中進速系數(shù)等于1 附近的3 個試驗點進行數(shù)值模擬。計算條件同試驗條件(見表1),表中pA為環(huán)境壓力。旨在獲得可靠的空化流場初始解,并對網(wǎng)格劃分及計算方法進行無空化計算校驗。

    表1 無空化計算條件Tab.1 Computational conditions w/o CAV

    計算在單相流框架下進行。湍流模型采用旋轉流適用的RNG k-ε 模型,近壁區(qū)的處理采用標準壁面函數(shù)。壓力速度耦合采用PISO(Pressure-Implicit with Splitting of Operators)算法,即隱式算子分割算法,該算法適用于高傾斜度網(wǎng)格??臻g離散均采用二階以上高精度離散方式:梯度擴散項采用基于單元體的格林—高斯離散;壓力采用高旋流適用的PRESTO 格式;動量方程和湍流方程均采用適用于四面體網(wǎng)格的二階迎風格式。表2 給出了推力系數(shù)KT和扭矩系數(shù)KQ的計算結果。由表可見:KT和KQ的計算誤差在3.0%~5.4%以內(nèi);隨著轉速提高(n=15~25 r/s),雷諾數(shù)增大(Re0.7=1.0×106~1.6×106,以0.7R 半徑處的合速度及槳葉弦長為特征量),計算精度略有降低。

    表2 無空化計算結果與誤差分析Tab.2 Results and error analysis w/o CAV

    為了準確定義并描述旋轉的流場,需要對漩渦結構進行定義。關于漩渦的判別法則有許多,比較流行的是渦量準則、Q 準則、Δ 準則和λ2準則。以上這些判別方式,都定義了一個可以逐點計算數(shù)值的函數(shù),然后判定每一點位于漩渦內(nèi)或者漩渦外。盡管判定方式不少,但是對渦進行準確而嚴格的定義卻是較為困難的。通常,漩渦被認為是渦量較高的區(qū)域,但是對該區(qū)域的準確閾值卻沒有嚴格的定義,隨著流動性質(zhì)的不同,閾值的判定將會嚴重影響結果。而且,在沒有漩渦的強剪切流動中渦量也可能較大,所以,利用渦量準則來判定漩渦是不夠準確的。利用流線來定義渦并進行描述也是不可靠的,因為即使在簡單的勻速平移運動中,流線的拓撲結構也容易被改變。而后三種方法(Q 準則、Δ 準則和λ2準則)則在漩渦判斷的方式上具有等價性。本文采用Q 準則對漩渦進行定義和直觀的表達。

    Q 準則由Hunt 等在1988 年提出[12],它定義流場中速度梯度張量▽V 的第二不變量Q 具有正值的區(qū)域為漩渦。另外,它要求漩渦區(qū)域的壓強低于周圍的壓強。Q 定義為:

    式中:Ω 為渦量,S 為變形率,皆為二階張量。符號“||||”表示張量的二范數(shù)。S 和Ω 分別代表了流場中一點的變形和旋轉。Q 準則反映了流場中一個流體微團旋轉和變形之間的一種平衡,Q>0 則反映了旋轉在流動中占據(jù)統(tǒng)治地位,同時也體現(xiàn)了該旋轉區(qū)域的壓力為附近區(qū)域的最小值(極小值)。圖3 給出了以轉速為特征量的無量綱等值面圖。由圖可見,采取本文所述網(wǎng)格生成方法得到的梢渦流場較之采用傳統(tǒng)的不對梢渦進行區(qū)域劃分的網(wǎng)格計算結果要顯著得多,所捕獲的梢渦長度要長得多:梢渦從導邊后緣拖出,并向螺旋槳后方移動并逐漸耗散;當值增大,梢渦區(qū)域向渦核中心線收縮。

    圖3 等值面Fig.3 ISO surfaces of

    4 空化流場計算

    螺旋槳空化性能預報對于分析螺旋槳的設計工況和非設計工況都意義重大。SMP’11 對工作狀態(tài)下的螺旋槳空化性能進行了觀測,均在轉速n=25 r/s 下進行,按不同的推力系數(shù)和空化數(shù)分為3 個案例,分別標注為Case 2.3.1、Case 2.3.2 和Case 2.3.3,試驗條件見表3。其中,Case 2.3.1 和Case 2.3.3 為非設計工況,Case 2.3.2 為設計工況。本章將對Case 2.3.1 和Case 2.3.2 進行數(shù)值模擬。計算條件同試驗條件。

    表3 空化計算條件Tab.3 Computational conditions with CAV

    表3 中σn為以旋轉線速度為特征速度的空化數(shù),定義如下:由此可以換算出環(huán)境壓力pA。

    先計算無空化情形,此時環(huán)境壓力可以取作一個大氣壓。校驗推力系數(shù)(扭矩系數(shù)缺失),并作為空化流場的初始解。然后將計算切換至均質(zhì)混合流框架下,改變環(huán)境壓力后,開啟空化模型并設置空化參數(shù)。此時,與無空化計算相比,增添了一個氣相體積分數(shù)的控制方程。該方程的空間離散選用QUICK(Quadratic upwind interpolation)格式。此格式用于六面體網(wǎng)格具有三階精度,用于四面體網(wǎng)格或混合網(wǎng)格時只具有二階精度。其他在無空化計算時已有的設置則保持不變。

    4.1 算例1(Case 2.3.1)

    該算例屬于非設計工況,載荷偏高情形。先進行無空化計算(此時pA=101.325 kPa),將計算得到的推力系數(shù)與試驗值進行比較。由表4 可見,推力系數(shù)的計算誤差為2.9%。重新設置環(huán)境壓力(此時pA=42.207 kPa),開啟空化計算。由表4 可見,空化將引起推力和扭矩系數(shù)下降;KT和KQ的計算誤差均在3.0%以內(nèi)。

    表4 計算結果與誤差分析(Case 2.3.1)Tab.4 Results and error analysis for Case 2.3.1

    空泡形態(tài)一般采用汽相體積分數(shù)琢的等值面來表示(0<琢≤1),因此空泡的大小長短與琢的取值相關。圖4 給出了空泡形態(tài)的計算結果(琢=0.2)與試驗結果(EFD)對比。其中,“SS”表示吸力面;“PS”表示壓力面;“SVSS”表示吸力面?zhèn)纫晥D。圖4(a)的試驗結果顯示,此工況下,除了發(fā)生顯著的梢渦空化外,吸力面葉根處也發(fā)生了片空化。圖4(b)~(d)的3 個視圖全方位地展示了梢渦空泡、吸力面葉根處的片空泡及轂渦空泡,計算與試驗符合較好。

    圖4 空泡形態(tài)計算與試驗比較Fig.4 Observed and calculated cavitation(Case 2.3.1)

    圖5 等值面Fig.5 ISO surfaces of (Case 2.3.1)

    4.2 算例2(Case 2.3.2)

    此工況為設計工況,中等載荷情形。先進行無空化計算(此時pA=101.325 kPa),將計算得到的推力系數(shù)與試驗值進行比較。由表5 可見,推力系數(shù)的計算誤差為4.5%。重新設置環(huán)境壓力(此時pA=30.528 kPa),開啟空化計算。由表5 可見:空化將引起推力和扭矩系數(shù)下降;KT的計算誤差為2.4%,KQ為4.5%。

    表5 計算結果與誤差分析(Case 2.3.2)Tab.5 Results and error analysis for Case 2.3.2

    圖6 給出了空泡形態(tài)的計算結果與試驗結果對比。圖6(a)的試驗結果顯示,此工況下的梢渦空泡在靠近隨邊處發(fā)生并拖出,吸力面葉根處發(fā)生了較大面積的片空化,并伴隨有泡狀空泡。圖6(b)~(d)分3 個視圖給出了數(shù)值模擬結果(琢=0.2)。由圖可見,除了梢渦空泡、吸力面葉根處的片空泡,壓力面葉根處也發(fā)生了較小的片空泡,計算與試驗符合較好。此工況下,沒有出現(xiàn)轂渦空泡。

    圖6 空泡形態(tài)計算與試驗比較Fig.6 Observed and calculated cavitation(Case 2.3.2)

    圖7 等值面Fig.7 ISO surfaces of (Case 2.3.2)

    4.3 水中含氣率對空化計算的影響

    數(shù)值實驗表明,當水中氣核體積分數(shù)琢nuc在一定范圍內(nèi)取不同值時,對空泡形態(tài)的影響并不顯著,卻會影響宏觀量—推力和扭矩系數(shù)的計算。表6 給出了算例Case 2.3.2 中琢nuc對KT和KQ的影響。由表可見,琢nuc越大,KT和KQ有降低的趨勢,從而更接近試驗值,但此時的計算容易發(fā)散。如何讓琢nuc的取值與模型試驗中的水中含氣量等環(huán)境變量相關聯(lián),或者說,進一步探究或改良現(xiàn)有空化模型是下一步研究需要解決的問題。

    表6 含氣率對KT 和KQ 的影響Tab.6 Effects of αnuc on KT&KQ

    5 結 語

    本文采用梢渦區(qū)域劃分的網(wǎng)格形式,對PPTC 槳進行了梢渦及梢渦空化模擬,在重載及中載工況的空化模擬中取得了初步的成效。通過數(shù)值模擬能得到顯著的梢渦及梢渦空泡,同時獲得較為準確的推力和扭矩系數(shù)。研究表明:梢渦空泡的捕捉很大程度上取決于梢渦的捕捉,梢渦空化只發(fā)生在渦量較大、壓力較低的渦核區(qū)域。

    猜你喜歡
    空泡空化螺旋槳
    功率超聲作用下鋼液中空化泡尺寸的演變特性
    鋼鐵釩鈦(2023年5期)2023-11-17 08:48:34
    水下航行體雙空泡相互作用數(shù)值模擬研究
    基于CFD的螺旋槳拉力確定方法
    三維扭曲水翼空化現(xiàn)象CFD模擬
    不同運動形式下水物相互作用空化數(shù)值模擬
    基于LPV的超空泡航行體H∞抗飽和控制
    基于CFD的對轉槳無空泡噪聲的仿真預報
    船海工程(2015年4期)2016-01-05 15:53:28
    3800DWT加油船螺旋槳諧鳴分析及消除方法
    廣東造船(2015年6期)2015-02-27 10:52:46
    螺旋槳轂帽鰭節(jié)能性能的數(shù)值模擬
    SPH在水下高速物體空泡發(fā)展模擬中的應用
    計算物理(2014年1期)2014-03-11 17:00:22
    99久久精品国产亚洲精品| 亚洲专区中文字幕在线 | 香蕉国产在线看| 建设人人有责人人尽责人人享有的| 丝瓜视频免费看黄片| 精品少妇久久久久久888优播| www日本在线高清视频| 晚上一个人看的免费电影| 欧美日本中文国产一区发布| 亚洲美女黄色视频免费看| 国产亚洲午夜精品一区二区久久| 丝袜脚勾引网站| 亚洲欧美色中文字幕在线| 色94色欧美一区二区| 超色免费av| 美女扒开内裤让男人捅视频| 国产精品99久久99久久久不卡 | 男女边吃奶边做爰视频| 在现免费观看毛片| 日韩电影二区| 亚洲精品在线美女| 欧美变态另类bdsm刘玥| 亚洲激情五月婷婷啪啪| 国产精品一区二区在线不卡| 亚洲欧美激情在线| 亚洲成av片中文字幕在线观看| 亚洲激情五月婷婷啪啪| 成人影院久久| 亚洲一区中文字幕在线| 日日啪夜夜爽| 狠狠婷婷综合久久久久久88av| 精品国产一区二区久久| 我的亚洲天堂| 国产欧美日韩一区二区三区在线| 成人亚洲精品一区在线观看| 国产精品秋霞免费鲁丝片| 一区二区三区四区激情视频| 国产精品久久久久久精品古装| 午夜免费男女啪啪视频观看| 国产黄频视频在线观看| 国产高清国产精品国产三级| 国产无遮挡羞羞视频在线观看| 亚洲欧美成人精品一区二区| 一区在线观看完整版| 两个人看的免费小视频| 国产成人欧美| av有码第一页| 最新在线观看一区二区三区 | 丝袜脚勾引网站| 日韩 欧美 亚洲 中文字幕| 国产伦人伦偷精品视频| 熟女少妇亚洲综合色aaa.| 中文字幕精品免费在线观看视频| 日韩 欧美 亚洲 中文字幕| 亚洲国产精品一区二区三区在线| 欧美中文综合在线视频| 亚洲熟女毛片儿| 国产精品久久久人人做人人爽| 永久免费av网站大全| 久久久久国产一级毛片高清牌| 美国免费a级毛片| 久久精品国产a三级三级三级| 久久久久久久久久久免费av| 亚洲精品乱久久久久久| 欧美精品亚洲一区二区| 欧美黑人精品巨大| 久久性视频一级片| 国产1区2区3区精品| 悠悠久久av| 欧美人与善性xxx| 国产黄频视频在线观看| 亚洲国产欧美网| 久久久久久久国产电影| 久久国产精品大桥未久av| 久久久精品国产亚洲av高清涩受| 一本大道久久a久久精品| 欧美av亚洲av综合av国产av | 青春草国产在线视频| 国产一卡二卡三卡精品 | 老司机影院成人| 美女中出高潮动态图| 亚洲一码二码三码区别大吗| 国产精品一区二区在线不卡| 亚洲美女搞黄在线观看| 亚洲欧美一区二区三区黑人| av一本久久久久| 久久久久国产精品人妻一区二区| 捣出白浆h1v1| av片东京热男人的天堂| 亚洲精品国产av成人精品| 免费黄频网站在线观看国产| 国产成人精品无人区| 午夜激情久久久久久久| 日韩 亚洲 欧美在线| 欧美变态另类bdsm刘玥| 国产极品粉嫩免费观看在线| 我要看黄色一级片免费的| 欧美精品高潮呻吟av久久| 国产女主播在线喷水免费视频网站| 大片电影免费在线观看免费| 国产日韩欧美视频二区| 国产精品二区激情视频| 成人国语在线视频| 一级黄片播放器| 久久亚洲国产成人精品v| 91aial.com中文字幕在线观看| 综合色丁香网| 涩涩av久久男人的天堂| 亚洲精品日韩在线中文字幕| 黄片播放在线免费| 伊人久久国产一区二区| 亚洲一区二区三区欧美精品| 美女福利国产在线| 国产 一区精品| 校园人妻丝袜中文字幕| 午夜精品国产一区二区电影| 国产免费福利视频在线观看| 9色porny在线观看| 国产一级毛片在线| 一级毛片黄色毛片免费观看视频| av一本久久久久| 妹子高潮喷水视频| 女人久久www免费人成看片| 七月丁香在线播放| 日韩中文字幕欧美一区二区 | 国产精品无大码| 免费观看性生交大片5| 在现免费观看毛片| 国产伦理片在线播放av一区| 国产成人免费观看mmmm| 亚洲一区中文字幕在线| 久久天躁狠狠躁夜夜2o2o | 日韩一区二区三区影片| 国产精品一区二区在线观看99| 日本wwww免费看| 亚洲精品国产av成人精品| 精品一品国产午夜福利视频| 成人黄色视频免费在线看| 少妇猛男粗大的猛烈进出视频| 国产乱人偷精品视频| 久久人妻熟女aⅴ| 男女边吃奶边做爰视频| 晚上一个人看的免费电影| 大陆偷拍与自拍| av又黄又爽大尺度在线免费看| 97精品久久久久久久久久精品| 不卡视频在线观看欧美| 久久久久视频综合| 国产精品一二三区在线看| 亚洲精品美女久久av网站| 精品国产一区二区久久| 国产精品99久久99久久久不卡 | 亚洲精品国产av成人精品| 欧美黑人精品巨大| 日本wwww免费看| 亚洲国产精品成人久久小说| 午夜福利一区二区在线看| 一区二区三区精品91| 美女福利国产在线| 久久天堂一区二区三区四区| 女人精品久久久久毛片| 麻豆精品久久久久久蜜桃| 性高湖久久久久久久久免费观看| 可以免费在线观看a视频的电影网站 | 蜜桃在线观看..| 男人操女人黄网站| 国产精品嫩草影院av在线观看| 热99久久久久精品小说推荐| 久久精品国产亚洲av涩爱| 日本色播在线视频| 老司机影院成人| 高清视频免费观看一区二区| 国产成人欧美在线观看 | 久久狼人影院| 久久精品久久久久久噜噜老黄| 国产精品熟女久久久久浪| 热99久久久久精品小说推荐| 国产精品亚洲av一区麻豆 | 亚洲欧洲国产日韩| 最近的中文字幕免费完整| 精品人妻熟女毛片av久久网站| 精品久久久久久电影网| 国产日韩欧美视频二区| 成年美女黄网站色视频大全免费| 两个人看的免费小视频| 精品一区二区三卡| 人人妻人人澡人人看| 日本欧美国产在线视频| 看免费成人av毛片| 国产亚洲最大av| 在现免费观看毛片| videosex国产| 成人18禁高潮啪啪吃奶动态图| av线在线观看网站| 久久性视频一级片| 久久久久精品性色| 国产精品一区二区在线观看99| 国产成人啪精品午夜网站| 啦啦啦啦在线视频资源| 亚洲国产精品一区二区三区在线| 国产av国产精品国产| 少妇精品久久久久久久| 一区二区三区精品91| 咕卡用的链子| 日本黄色日本黄色录像| 夫妻午夜视频| 亚洲国产毛片av蜜桃av| 国产精品偷伦视频观看了| 美国免费a级毛片| a级片在线免费高清观看视频| 欧美日韩福利视频一区二区| 一级黄片播放器| 久久国产精品大桥未久av| 国产片内射在线| 自线自在国产av| 久久天躁狠狠躁夜夜2o2o | 伊人久久国产一区二区| 日本色播在线视频| 麻豆乱淫一区二区| 制服丝袜香蕉在线| 精品卡一卡二卡四卡免费| 一区二区日韩欧美中文字幕| 七月丁香在线播放| 18禁观看日本| 一区二区三区精品91| 亚洲精品国产av成人精品| 午夜福利,免费看| 桃花免费在线播放| 精品国产超薄肉色丝袜足j| 国产有黄有色有爽视频| 国产av国产精品国产| 在线观看免费视频网站a站| 最黄视频免费看| 欧美日韩av久久| 成人免费观看视频高清| 天堂8中文在线网| 午夜福利视频在线观看免费| 少妇被粗大猛烈的视频| 天堂俺去俺来也www色官网| 美女脱内裤让男人舔精品视频| 制服人妻中文乱码| 涩涩av久久男人的天堂| 亚洲av成人精品一二三区| 亚洲欧美一区二区三区国产| 人体艺术视频欧美日本| 高清不卡的av网站| 亚洲精品一二三| 国产精品亚洲av一区麻豆 | 国产成人精品无人区| 免费观看性生交大片5| h视频一区二区三区| 又大又黄又爽视频免费| 亚洲精品成人av观看孕妇| 国产男女超爽视频在线观看| 黄片播放在线免费| 久久天堂一区二区三区四区| 亚洲国产中文字幕在线视频| 18禁动态无遮挡网站| 交换朋友夫妻互换小说| 亚洲精品国产av成人精品| 久久久久精品人妻al黑| 国产片内射在线| 一个人免费看片子| 99热网站在线观看| 999久久久国产精品视频| 视频在线观看一区二区三区| 午夜福利影视在线免费观看| 成人免费观看视频高清| 一区二区三区四区激情视频| 精品亚洲成国产av| 亚洲av日韩在线播放| 日日啪夜夜爽| 亚洲精品一二三| 考比视频在线观看| 午夜日本视频在线| 国产免费视频播放在线视频| 在线观看一区二区三区激情| 妹子高潮喷水视频| 亚洲在久久综合| 午夜av观看不卡| 高清不卡的av网站| 丝瓜视频免费看黄片| 最近中文字幕2019免费版| 亚洲国产欧美在线一区| 亚洲国产av影院在线观看| 一边摸一边抽搐一进一出视频| 亚洲第一青青草原| 免费在线观看黄色视频的| 免费观看性生交大片5| 伊人久久国产一区二区| 精品福利永久在线观看| 久久久久久久久免费视频了| 国产乱人偷精品视频| 亚洲精品视频女| 国产成人a∨麻豆精品| 亚洲成人免费av在线播放| 美女脱内裤让男人舔精品视频| 一级毛片电影观看| 悠悠久久av| 国产成人精品久久二区二区91 | 91精品伊人久久大香线蕉| 亚洲国产毛片av蜜桃av| 亚洲精品视频女| 国产人伦9x9x在线观看| 一边摸一边做爽爽视频免费| 19禁男女啪啪无遮挡网站| 成人国产av品久久久| 精品少妇内射三级| av线在线观看网站| 亚洲精品aⅴ在线观看| av一本久久久久| 丝瓜视频免费看黄片| 国产精品久久久久久人妻精品电影 | 777米奇影视久久| 热99国产精品久久久久久7| 国产精品香港三级国产av潘金莲 | 美女扒开内裤让男人捅视频| 欧美日韩亚洲高清精品| 亚洲av电影在线进入| 中文字幕色久视频| 一二三四在线观看免费中文在| av有码第一页| 激情五月婷婷亚洲| 亚洲精品国产av蜜桃| 免费高清在线观看日韩| 国产一区亚洲一区在线观看| 亚洲成人手机| 久久天堂一区二区三区四区| 国产黄频视频在线观看| 丝袜美足系列| 老汉色∧v一级毛片| 交换朋友夫妻互换小说| 精品少妇久久久久久888优播| 天天影视国产精品| 最新的欧美精品一区二区| 中文字幕av电影在线播放| 一区二区日韩欧美中文字幕| 交换朋友夫妻互换小说| 最近中文字幕高清免费大全6| 欧美97在线视频| 国产女主播在线喷水免费视频网站| 男男h啪啪无遮挡| 国产深夜福利视频在线观看| 久久精品人人爽人人爽视色| 亚洲综合色网址| 亚洲精品av麻豆狂野| 亚洲三区欧美一区| 国产成人一区二区在线| 亚洲国产av影院在线观看| 国产黄色视频一区二区在线观看| 亚洲成国产人片在线观看| 在现免费观看毛片| 啦啦啦 在线观看视频| 一级片'在线观看视频| 999久久久国产精品视频| 亚洲激情五月婷婷啪啪| 卡戴珊不雅视频在线播放| 国产野战对白在线观看| 国产一区二区三区av在线| 亚洲一区中文字幕在线| av有码第一页| 操出白浆在线播放| 男人爽女人下面视频在线观看| 欧美黄色片欧美黄色片| av片东京热男人的天堂| 亚洲成人一二三区av| 天天操日日干夜夜撸| 国产成人av激情在线播放| 新久久久久国产一级毛片| 欧美日韩成人在线一区二区| 精品久久久精品久久久| 高清视频免费观看一区二区| 中国三级夫妇交换| 亚洲精品成人av观看孕妇| 欧美日韩综合久久久久久| 高清不卡的av网站| 人妻 亚洲 视频| 菩萨蛮人人尽说江南好唐韦庄| 免费黄色在线免费观看| 亚洲国产精品一区二区三区在线| 91精品伊人久久大香线蕉| 蜜桃在线观看..| 国产一区二区在线观看av| 国产日韩欧美视频二区| 亚洲第一区二区三区不卡| 成人影院久久| 国产日韩欧美在线精品| 在线观看免费日韩欧美大片| 中文字幕高清在线视频| 丁香六月天网| 欧美精品一区二区免费开放| 老司机在亚洲福利影院| 精品酒店卫生间| 国产精品一区二区在线不卡| 亚洲欧美色中文字幕在线| 中文欧美无线码| 中文乱码字字幕精品一区二区三区| 久久久久久久久久久免费av| 亚洲色图综合在线观看| 精品视频人人做人人爽| 看非洲黑人一级黄片| 妹子高潮喷水视频| 午夜日本视频在线| 国产精品免费视频内射| 最近中文字幕2019免费版| 中文字幕最新亚洲高清| 精品亚洲成国产av| 成人黄色视频免费在线看| 国产精品免费视频内射| 国产伦人伦偷精品视频| 91精品伊人久久大香线蕉| a级片在线免费高清观看视频| 午夜激情久久久久久久| 高清在线视频一区二区三区| 看非洲黑人一级黄片| 久久人人爽人人片av| av.在线天堂| 啦啦啦视频在线资源免费观看| 免费在线观看完整版高清| 国产精品三级大全| 巨乳人妻的诱惑在线观看| 又大又黄又爽视频免费| 一边摸一边做爽爽视频免费| 国产精品一区二区在线观看99| 国产亚洲午夜精品一区二区久久| 97在线人人人人妻| 亚洲精品久久久久久婷婷小说| 老司机在亚洲福利影院| 国产av精品麻豆| 国产免费一区二区三区四区乱码| 日韩 欧美 亚洲 中文字幕| 久久久精品国产亚洲av高清涩受| 下体分泌物呈黄色| 美女脱内裤让男人舔精品视频| 搡老乐熟女国产| 97在线人人人人妻| 亚洲天堂av无毛| 欧美人与性动交α欧美软件| 久久女婷五月综合色啪小说| 啦啦啦中文免费视频观看日本| 国产伦理片在线播放av一区| 日韩大码丰满熟妇| 美女福利国产在线| 精品少妇黑人巨大在线播放| 男女高潮啪啪啪动态图| 亚洲国产精品国产精品| 香蕉国产在线看| 欧美日韩精品网址| 国产av码专区亚洲av| 久久99热这里只频精品6学生| 国产野战对白在线观看| 久久99一区二区三区| 日韩制服骚丝袜av| 国产爽快片一区二区三区| 国产福利在线免费观看视频| 制服人妻中文乱码| 一区福利在线观看| 成人国产av品久久久| 精品一区在线观看国产| 十八禁高潮呻吟视频| 最新的欧美精品一区二区| 亚洲精品日本国产第一区| 亚洲欧美精品自产自拍| 欧美乱码精品一区二区三区| 啦啦啦 在线观看视频| 精品一区在线观看国产| 午夜影院在线不卡| 久久这里只有精品19| 天天躁日日躁夜夜躁夜夜| 在线天堂中文资源库| 久久人人爽av亚洲精品天堂| 精品一区二区三卡| 日韩大片免费观看网站| 国产在线一区二区三区精| 精品国产乱码久久久久久小说| 亚洲精品国产av蜜桃| 99香蕉大伊视频| 超碰97精品在线观看| 国产精品无大码| 久久久久国产精品人妻一区二区| 国产精品av久久久久免费| 午夜福利视频精品| 久久av网站| 久久人人爽av亚洲精品天堂| 亚洲欧美成人精品一区二区| 中文字幕最新亚洲高清| 在线观看免费视频网站a站| 国产一区亚洲一区在线观看| 美女午夜性视频免费| 亚洲av福利一区| 无遮挡黄片免费观看| 视频区图区小说| 成年av动漫网址| 精品一区二区免费观看| 蜜桃国产av成人99| 久久性视频一级片| 一级黄片播放器| 一级片免费观看大全| 日本欧美国产在线视频| 肉色欧美久久久久久久蜜桃| 男男h啪啪无遮挡| videosex国产| 毛片一级片免费看久久久久| 午夜福利免费观看在线| 午夜福利视频精品| www.av在线官网国产| 久久热在线av| 午夜福利,免费看| 青青草视频在线视频观看| 国产一卡二卡三卡精品 | 最近2019中文字幕mv第一页| 男女边摸边吃奶| 97精品久久久久久久久久精品| 交换朋友夫妻互换小说| 欧美日韩av久久| 一级片免费观看大全| 亚洲国产精品成人久久小说| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久成人av| 国产成人精品久久久久久| 亚洲国产最新在线播放| 欧美最新免费一区二区三区| 亚洲国产欧美网| 老司机影院毛片| 女人爽到高潮嗷嗷叫在线视频| 国产免费现黄频在线看| 一二三四中文在线观看免费高清| 久久久久精品人妻al黑| 9色porny在线观看| 熟女av电影| 黄片播放在线免费| 99re6热这里在线精品视频| 天美传媒精品一区二区| 中文字幕人妻丝袜一区二区 | 国产精品.久久久| 亚洲婷婷狠狠爱综合网| 色网站视频免费| av国产久精品久网站免费入址| 欧美最新免费一区二区三区| 久久精品亚洲熟妇少妇任你| 亚洲天堂av无毛| 国产免费又黄又爽又色| 久久97久久精品| 天天操日日干夜夜撸| 三上悠亚av全集在线观看| 中国三级夫妇交换| 欧美人与善性xxx| 亚洲四区av| 免费观看a级毛片全部| 欧美精品高潮呻吟av久久| 久久久久久久精品精品| 纵有疾风起免费观看全集完整版| 天天躁夜夜躁狠狠躁躁| 一级毛片电影观看| 成人亚洲欧美一区二区av| 欧美成人午夜精品| 国产精品亚洲av一区麻豆 | 一级毛片我不卡| 大香蕉久久成人网| 丝袜喷水一区| 久久精品久久精品一区二区三区| 欧美精品一区二区大全| 狠狠精品人妻久久久久久综合| svipshipincom国产片| 夫妻性生交免费视频一级片| 一区二区三区激情视频| 一本—道久久a久久精品蜜桃钙片| 人妻一区二区av| 色婷婷久久久亚洲欧美| 99香蕉大伊视频| 久热这里只有精品99| 99久久99久久久精品蜜桃| 只有这里有精品99| 午夜福利视频在线观看免费| 日本av手机在线免费观看| 欧美日韩一区二区视频在线观看视频在线| 婷婷成人精品国产| 国产激情久久老熟女| 国产在线视频一区二区| 最近2019中文字幕mv第一页| 亚洲av男天堂| 精品人妻熟女毛片av久久网站| 欧美日韩av久久| 色婷婷av一区二区三区视频| 91精品三级在线观看| 黄色 视频免费看| 一边亲一边摸免费视频| 夫妻性生交免费视频一级片| 精品国产超薄肉色丝袜足j| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品美女久久av网站| 黄色视频在线播放观看不卡| 亚洲av男天堂| 老汉色av国产亚洲站长工具| 国产精品久久久久久精品电影小说| 久久久精品国产亚洲av高清涩受| 免费在线观看完整版高清| 五月天丁香电影| 国产一区二区在线观看av| 人人妻人人爽人人添夜夜欢视频| 男人爽女人下面视频在线观看| 日韩制服骚丝袜av| 男女床上黄色一级片免费看| 女性被躁到高潮视频| 日韩制服骚丝袜av| 国产精品二区激情视频| a 毛片基地| 国产成人免费观看mmmm| 国产一区二区三区综合在线观看| 亚洲欧美中文字幕日韩二区| 精品久久久久久电影网| 亚洲在久久综合| 一二三四在线观看免费中文在| 午夜老司机福利片| 高清视频免费观看一区二区|