• <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
    or卡值多少钱| 在线免费观看的www视频| 日本黄大片高清| 亚洲第一区二区三区不卡| 亚洲中文字幕日韩| 亚洲av二区三区四区| 亚洲专区国产一区二区| 在线a可以看的网站| 国产黄片美女视频| 成人特级黄色片久久久久久久| 亚洲电影在线观看av| 成人毛片a级毛片在线播放| 国产精品综合久久久久久久免费| 久久精品人妻少妇| 大型黄色视频在线免费观看| 自拍偷自拍亚洲精品老妇| 男人舔女人下体高潮全视频| 美女黄网站色视频| 91麻豆精品激情在线观看国产| 亚洲精品成人久久久久久| 少妇被粗大猛烈的视频| 日本五十路高清| 国产亚洲精品久久久com| 一级黄色大片毛片| 日日啪夜夜撸| 欧美成人a在线观看| 波多野结衣巨乳人妻| 有码 亚洲区| 18禁黄网站禁片免费观看直播| 天天一区二区日本电影三级| 搞女人的毛片| 国产成人freesex在线 | 小蜜桃在线观看免费完整版高清| 中文资源天堂在线| h日本视频在线播放| 欧美日韩一区二区视频在线观看视频在线 | 国产一区二区亚洲精品在线观看| 舔av片在线| 日本爱情动作片www.在线观看 | 成人高潮视频无遮挡免费网站| 亚洲av五月六月丁香网| 国产老妇女一区| 你懂的网址亚洲精品在线观看 | 六月丁香七月| 性色avwww在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲人与动物交配视频| 中文字幕免费在线视频6| 99久久成人亚洲精品观看| 国产精品野战在线观看| 伦精品一区二区三区| 欧美激情久久久久久爽电影| 色5月婷婷丁香| 三级经典国产精品| 成年版毛片免费区| 在线a可以看的网站| 国内少妇人妻偷人精品xxx网站| 亚洲精品一区av在线观看| 欧美日韩在线观看h| 九九久久精品国产亚洲av麻豆| 少妇的逼好多水| 综合色丁香网| 日韩av在线大香蕉| 精品欧美国产一区二区三| 麻豆国产av国片精品| 欧美极品一区二区三区四区| 亚洲av成人精品一区久久| 91狼人影院| 亚洲电影在线观看av| av国产免费在线观看| av免费在线看不卡| av在线天堂中文字幕| 久久九九热精品免费| 久久人人爽人人片av| 天堂网av新在线| 黄色一级大片看看| ponron亚洲| 18禁黄网站禁片免费观看直播| 我的女老师完整版在线观看| 激情 狠狠 欧美| 性插视频无遮挡在线免费观看| 国产人妻一区二区三区在| 亚洲欧美精品综合久久99| 天堂影院成人在线观看| 精品久久久久久久人妻蜜臀av| 亚洲av五月六月丁香网| 日韩av不卡免费在线播放| 久久精品国产清高在天天线| 女生性感内裤真人,穿戴方法视频| 亚洲av熟女| 能在线免费观看的黄片| 别揉我奶头 嗯啊视频| 国产亚洲精品av在线| 国产亚洲精品久久久久久毛片| 久久精品国产自在天天线| 久久久久久国产a免费观看| av视频在线观看入口| 国产成人aa在线观看| 午夜久久久久精精品| 在线国产一区二区在线| 九九在线视频观看精品| 国产真实伦视频高清在线观看| 国产一区二区在线av高清观看| 午夜亚洲福利在线播放| 狂野欧美白嫩少妇大欣赏| 久久草成人影院| 人妻丰满熟妇av一区二区三区| 97超碰精品成人国产| 夜夜看夜夜爽夜夜摸| 欧美+亚洲+日韩+国产| 久久久久久久久大av| 99九九线精品视频在线观看视频| videossex国产| 蜜桃亚洲精品一区二区三区| 啦啦啦观看免费观看视频高清| 久久久欧美国产精品| 国产精品亚洲美女久久久| 亚洲国产高清在线一区二区三| 午夜福利18| 丰满人妻一区二区三区视频av| 九九热线精品视视频播放| 日韩欧美精品免费久久| 亚洲国产精品成人久久小说 | 精品日产1卡2卡| 国产乱人偷精品视频| 人人妻人人看人人澡| 精品少妇黑人巨大在线播放 | 国产人妻一区二区三区在| 久久精品国产鲁丝片午夜精品| 国产精品人妻久久久久久| 亚洲欧美精品自产自拍| 搡老熟女国产l中国老女人| a级一级毛片免费在线观看| 国产白丝娇喘喷水9色精品| 日产精品乱码卡一卡2卡三| 五月玫瑰六月丁香| 国产精品久久久久久av不卡| 又爽又黄无遮挡网站| 免费一级毛片在线播放高清视频| 免费看光身美女| 在线观看av片永久免费下载| 一区福利在线观看| 舔av片在线| 一进一出抽搐动态| 成人美女网站在线观看视频| 老熟妇乱子伦视频在线观看| 精品人妻偷拍中文字幕| 99热网站在线观看| 精品久久国产蜜桃| 男人和女人高潮做爰伦理| 中出人妻视频一区二区| 国产视频内射| 啦啦啦啦在线视频资源| 欧美最黄视频在线播放免费| 国语自产精品视频在线第100页| 2021天堂中文幕一二区在线观| 久久久久久久久大av| 日本a在线网址| a级毛色黄片| 又黄又爽又刺激的免费视频.| 亚洲最大成人中文| 人人妻人人澡欧美一区二区| 久久久久久国产a免费观看| 春色校园在线视频观看| 久久精品国产99精品国产亚洲性色| 日韩欧美国产在线观看| 一个人看视频在线观看www免费| 久久精品91蜜桃| 久久精品国产自在天天线| 久久精品综合一区二区三区| 性色avwww在线观看| 成年av动漫网址| 国产探花极品一区二区| 亚洲一区二区三区色噜噜| 国产蜜桃级精品一区二区三区| 超碰av人人做人人爽久久| 熟女人妻精品中文字幕| 色尼玛亚洲综合影院| 亚洲国产日韩欧美精品在线观看| 欧美三级亚洲精品| 久久这里只有精品中国| 国产黄色小视频在线观看| 99在线视频只有这里精品首页| 国产 一区 欧美 日韩| 熟妇人妻久久中文字幕3abv| 在线天堂最新版资源| eeuss影院久久| 女生性感内裤真人,穿戴方法视频| 亚洲最大成人手机在线| 免费看日本二区| 亚洲av中文字字幕乱码综合| 欧美不卡视频在线免费观看| 国产三级在线视频| АⅤ资源中文在线天堂| 夜夜爽天天搞| 久久久国产成人精品二区| 午夜免费男女啪啪视频观看 | videossex国产| 成熟少妇高潮喷水视频| 赤兔流量卡办理| 国产成人精品久久久久久| 黄片wwwwww| 床上黄色一级片| 日本一二三区视频观看| 校园春色视频在线观看| 久久鲁丝午夜福利片| 亚洲成人中文字幕在线播放| 在线观看美女被高潮喷水网站| 久久99热6这里只有精品| 国产亚洲精品综合一区在线观看| 午夜福利高清视频| 国产成人福利小说| 人妻少妇偷人精品九色| 精品久久久噜噜| 日本成人三级电影网站| www日本黄色视频网| 男女之事视频高清在线观看| 日韩 亚洲 欧美在线| 男女那种视频在线观看| 成人亚洲欧美一区二区av| 搡老岳熟女国产| 一个人观看的视频www高清免费观看| 久久精品国产99精品国产亚洲性色| 中文字幕av在线有码专区| 男人舔奶头视频| 精品久久久久久久久av| 一级毛片我不卡| 久久精品国产鲁丝片午夜精品| 亚洲欧美日韩东京热| 免费看a级黄色片| 露出奶头的视频| 51国产日韩欧美| 久久久精品大字幕| 欧美日韩一区二区视频在线观看视频在线 | 国产精品日韩av在线免费观看| 亚洲国产精品成人久久小说 | 亚洲无线观看免费| 国产精品久久久久久亚洲av鲁大| 欧美性猛交黑人性爽| 三级毛片av免费| 亚洲人成网站在线观看播放| 欧美一区二区精品小视频在线| 五月玫瑰六月丁香| 麻豆国产av国片精品| 日韩一区二区视频免费看| 两个人的视频大全免费| 亚洲av美国av| 亚洲va在线va天堂va国产| 亚洲人成网站高清观看| 黄色一级大片看看| 欧美日韩精品成人综合77777| 久久久精品欧美日韩精品| 国产精品一二三区在线看| 久久久欧美国产精品| 又粗又爽又猛毛片免费看| 床上黄色一级片| 国产精品美女特级片免费视频播放器| 99热精品在线国产| 日韩欧美 国产精品| 国产高清视频在线观看网站| 毛片女人毛片| 国产 一区 欧美 日韩| 亚洲第一区二区三区不卡| 亚洲av免费在线观看| 麻豆国产av国片精品| 99热精品在线国产| 国内少妇人妻偷人精品xxx网站| 亚洲国产欧美人成| 国产人妻一区二区三区在| 联通29元200g的流量卡| 日韩一本色道免费dvd| 国产一区亚洲一区在线观看| 色吧在线观看| а√天堂www在线а√下载| 在线a可以看的网站| 禁无遮挡网站| 最近最新中文字幕大全电影3| 一进一出好大好爽视频| 99久久精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 免费高清视频大片| 成人三级黄色视频| 国产精品一区二区性色av| 成年av动漫网址| 亚洲va在线va天堂va国产| 久久人人爽人人爽人人片va| 天天躁夜夜躁狠狠久久av| 色5月婷婷丁香| 欧美丝袜亚洲另类| 久久婷婷人人爽人人干人人爱| 午夜福利高清视频| 三级国产精品欧美在线观看| 亚洲性久久影院| 国产成人a区在线观看| 色综合色国产| 亚洲成人精品中文字幕电影| 久久亚洲精品不卡| 色综合色国产| 一区二区三区高清视频在线| 色综合亚洲欧美另类图片| 亚洲国产色片| 久久鲁丝午夜福利片| 国产成年人精品一区二区| 国产片特级美女逼逼视频| 欧美潮喷喷水| 欧美+亚洲+日韩+国产| 国产一区二区亚洲精品在线观看| 午夜爱爱视频在线播放| 日日摸夜夜添夜夜添av毛片| 亚洲欧美精品自产自拍| 亚洲内射少妇av| 老司机福利观看| 国产成人freesex在线 | 美女cb高潮喷水在线观看| 99久国产av精品| 国产黄片美女视频| 男女做爰动态图高潮gif福利片| 国产精品女同一区二区软件| 亚洲成人久久爱视频| 国产探花在线观看一区二区| 日韩一本色道免费dvd| 亚洲熟妇熟女久久| 亚洲欧美精品综合久久99| 免费观看在线日韩| 观看美女的网站| 高清毛片免费看| 人妻夜夜爽99麻豆av| 偷拍熟女少妇极品色| 99热这里只有是精品在线观看| 日日摸夜夜添夜夜添小说| 村上凉子中文字幕在线| 久久午夜福利片| 国产精品人妻久久久久久| 日本免费一区二区三区高清不卡| avwww免费| 午夜福利在线观看吧| 校园春色视频在线观看| 在线天堂最新版资源| 久久久久久大精品| 内地一区二区视频在线| 国产极品精品免费视频能看的| 一级毛片久久久久久久久女| 国产精品永久免费网站| 久久精品人妻少妇| 日韩av在线大香蕉| 中文字幕人妻熟人妻熟丝袜美| 日本黄大片高清| 久久精品综合一区二区三区| 干丝袜人妻中文字幕| 欧美日韩综合久久久久久| 九九热线精品视视频播放| 久久99热6这里只有精品| av福利片在线观看| 国产一区二区激情短视频| 如何舔出高潮| 亚洲一区二区三区色噜噜| 亚洲av一区综合| 男女那种视频在线观看| av在线观看视频网站免费| 久久久久九九精品影院| 国产免费男女视频| 国产一区亚洲一区在线观看| 免费看日本二区| 小说图片视频综合网站| 最近视频中文字幕2019在线8| 别揉我奶头 嗯啊视频| 国产视频一区二区在线看| 国产精品精品国产色婷婷| 九九爱精品视频在线观看| 日韩人妻高清精品专区| 国产成人a∨麻豆精品| 国产亚洲精品av在线| 免费av毛片视频| 99热这里只有是精品在线观看| 国产女主播在线喷水免费视频网站 | 精品熟女少妇av免费看| 淫秽高清视频在线观看| 日本三级黄在线观看| 寂寞人妻少妇视频99o| 99热只有精品国产| 一夜夜www| 在线免费观看不下载黄p国产| 日韩精品中文字幕看吧| 女生性感内裤真人,穿戴方法视频| 精品99又大又爽又粗少妇毛片| 久久精品人妻少妇| 亚洲av免费高清在线观看| 亚洲性久久影院| 国产国拍精品亚洲av在线观看| 成熟少妇高潮喷水视频| 国产一级毛片七仙女欲春2| 亚洲aⅴ乱码一区二区在线播放| 淫妇啪啪啪对白视频| 91久久精品电影网| 九色成人免费人妻av| 亚洲最大成人手机在线| 午夜久久久久精精品| 大香蕉久久网| av国产免费在线观看| 色播亚洲综合网| 麻豆一二三区av精品| 久久人人爽人人爽人人片va| 别揉我奶头 嗯啊视频| 亚洲自拍偷在线| 国产亚洲91精品色在线| av视频在线观看入口| 国产精品美女特级片免费视频播放器| 精品久久久久久成人av| 日韩一区二区视频免费看| 少妇高潮的动态图| 身体一侧抽搐| 国产不卡一卡二| 在线免费观看的www视频| 我要看日韩黄色一级片| 卡戴珊不雅视频在线播放| 精品人妻熟女av久视频| 亚洲国产精品国产精品| 国产成人一区二区在线| 久久精品91蜜桃| 少妇人妻一区二区三区视频| 精品午夜福利在线看| 卡戴珊不雅视频在线播放| 久久精品国产亚洲网站| 久久人妻av系列| 少妇熟女aⅴ在线视频| 午夜福利高清视频| 香蕉av资源在线| 免费观看人在逋| 国产黄色视频一区二区在线观看 | 国产精品1区2区在线观看.| 最好的美女福利视频网| 亚洲18禁久久av| 久久韩国三级中文字幕| 99热6这里只有精品| 国产又黄又爽又无遮挡在线| 不卡视频在线观看欧美| 男女啪啪激烈高潮av片| 久久精品国产清高在天天线| 亚洲av一区综合| 国产黄色视频一区二区在线观看 | 欧美一区二区精品小视频在线| 久久亚洲精品不卡| 超碰av人人做人人爽久久| 国产成人aa在线观看| 丝袜喷水一区| 亚洲欧美成人综合另类久久久 | 国产91av在线免费观看| 日韩欧美在线乱码| 国产女主播在线喷水免费视频网站 | 欧美精品国产亚洲| 日韩制服骚丝袜av| 99热这里只有是精品在线观看| 黄色配什么色好看| 青春草视频在线免费观看| 变态另类成人亚洲欧美熟女| 美女大奶头视频| 我要看日韩黄色一级片| 国产麻豆成人av免费视频| 亚洲在线观看片| 99视频精品全部免费 在线| 久久亚洲精品不卡| 日日摸夜夜添夜夜爱| 精品日产1卡2卡| 亚洲欧美日韩东京热| 午夜福利18| 国产一区二区三区在线臀色熟女| 日韩 亚洲 欧美在线| 成人综合一区亚洲| 你懂的网址亚洲精品在线观看 | 欧美一区二区亚洲| 久久人人爽人人片av| 日韩欧美在线乱码| 此物有八面人人有两片| 九色成人免费人妻av| 成人性生交大片免费视频hd| 国产aⅴ精品一区二区三区波| 网址你懂的国产日韩在线| 热99在线观看视频| 免费人成视频x8x8入口观看| 久久99热6这里只有精品| 搡老熟女国产l中国老女人| 欧美日本亚洲视频在线播放| 国产精品女同一区二区软件| 又爽又黄无遮挡网站| 亚洲国产精品合色在线| 91av网一区二区| 丝袜美腿在线中文| 免费黄网站久久成人精品| 白带黄色成豆腐渣| 亚洲欧美成人综合另类久久久 | 99热这里只有是精品50| 人妻夜夜爽99麻豆av| 久久久久久久午夜电影| 久久精品夜色国产| 欧美日韩国产亚洲二区| 色在线成人网| 亚洲aⅴ乱码一区二区在线播放| 亚洲美女搞黄在线观看 | 综合色av麻豆| 欧美丝袜亚洲另类| 不卡视频在线观看欧美| 精品国产三级普通话版| 亚洲精品国产av成人精品 | 亚洲自拍偷在线| 色播亚洲综合网| 观看免费一级毛片| 老熟妇仑乱视频hdxx| 中国国产av一级| 中文字幕免费在线视频6| 国产精品1区2区在线观看.| 97超视频在线观看视频| 可以在线观看的亚洲视频| 久久久久久久久中文| 色在线成人网| 欧美成人a在线观看| 国产 一区 欧美 日韩| 欧美一区二区精品小视频在线| 嫩草影院精品99| 成人二区视频| 国产视频一区二区在线看| 少妇的逼好多水| 啦啦啦啦在线视频资源| 亚洲av二区三区四区| 亚洲av美国av| 最近的中文字幕免费完整| 国产精品亚洲美女久久久| 久久精品国产亚洲av香蕉五月| 久久午夜亚洲精品久久| 久久久欧美国产精品| 91久久精品国产一区二区三区| 小蜜桃在线观看免费完整版高清| 有码 亚洲区| 在线国产一区二区在线| 国产成年人精品一区二区| 色尼玛亚洲综合影院| 五月伊人婷婷丁香| 美女 人体艺术 gogo| 国产亚洲精品久久久久久毛片| 嫩草影院新地址| 午夜日韩欧美国产| 十八禁网站免费在线| 青春草视频在线免费观看| 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 国内少妇人妻偷人精品xxx网站| а√天堂www在线а√下载| 又黄又爽又免费观看的视频| eeuss影院久久| 别揉我奶头~嗯~啊~动态视频| 国产成人一区二区在线| 美女黄网站色视频| 久久人人爽人人片av| 精品不卡国产一区二区三区| 激情 狠狠 欧美| 久久久a久久爽久久v久久| 日韩欧美精品v在线| 男插女下体视频免费在线播放| 麻豆一二三区av精品| 亚洲美女视频黄频| 狂野欧美白嫩少妇大欣赏| av天堂中文字幕网| 欧美性猛交╳xxx乱大交人| 国产 一区 欧美 日韩| 丰满的人妻完整版| 午夜久久久久精精品| 久久精品影院6| 亚洲人成网站在线观看播放| 男女下面进入的视频免费午夜| 精品日产1卡2卡| 日本色播在线视频| 国产欧美日韩一区二区精品| 欧美+日韩+精品| www.色视频.com| 极品教师在线视频| 97人妻精品一区二区三区麻豆| 久久精品国产亚洲网站| 免费观看在线日韩| 免费搜索国产男女视频| 日日摸夜夜添夜夜添小说| 色噜噜av男人的天堂激情| 六月丁香七月| 久久久欧美国产精品| 日本精品一区二区三区蜜桃| 亚洲久久久久久中文字幕| 亚洲国产高清在线一区二区三| 日韩成人av中文字幕在线观看 | 黄片wwwwww| 少妇熟女欧美另类| 全区人妻精品视频| 卡戴珊不雅视频在线播放| 亚洲国产日韩欧美精品在线观看| 国产精品一区二区三区四区久久| 国产久久久一区二区三区| 国产真实伦视频高清在线观看| 久久久国产成人免费| 在线观看免费视频日本深夜| 我要搜黄色片| 日韩欧美在线乱码| 亚洲精品粉嫩美女一区| 国产久久久一区二区三区| 美女大奶头视频| 国产日本99.免费观看| 日本黄大片高清| 高清毛片免费观看视频网站| 黄色配什么色好看| 美女xxoo啪啪120秒动态图| 欧美国产日韩亚洲一区| av中文乱码字幕在线| 1024手机看黄色片| 亚洲国产高清在线一区二区三| 国产精品久久久久久av不卡| 亚洲欧美日韩东京热| 超碰av人人做人人爽久久| 国产精品人妻久久久久久| 国产精品三级大全|