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

    吊艙式CRP推進器水動力性能數(shù)值模擬

    2012-03-08 06:41:46立,熊鷹,楊
    艦船科學(xué)技術(shù) 2012年10期
    關(guān)鍵詞:吊艙推進器槳葉

    盛 立,熊 鷹,楊 勇

    (海軍工程大學(xué) 船舶與海洋工程系,湖北 武漢 430033)

    吊艙式CRP推進器水動力性能數(shù)值模擬

    盛 立,熊 鷹,楊 勇

    (海軍工程大學(xué) 船舶與海洋工程系,湖北 武漢 430033)

    建立吊艙式CRP推進器數(shù)值模型,結(jié)合RANS方程和SST k-ω湍流模型,運用滑移網(wǎng)格方法對吊艙式CRP推進器在均勻流場中水動力性能進行非定常數(shù)值預(yù)報。將數(shù)值預(yù)報所得的敞水性能結(jié)果與在真實空泡水洞內(nèi)利用吊艙動力儀及長軸動力儀對吊艙式CRP推進器進行敞水試驗得到的試驗數(shù)據(jù)進行比較;同時得到了吊艙式CRP推進器前后槳葉面及葉背壓力系數(shù)分布與前后槳及吊艙的非定常性能,將計算結(jié)果和不附帶吊艙相同對轉(zhuǎn)槳計算結(jié)果進行比較分析。結(jié)果表明,本文所用數(shù)值計算方法對吊艙式CRP推進器水動力性能的預(yù)報具有較高的可信度,能達到工程應(yīng)用的要求。

    吊艙式CRP推進器;水動力性能;數(shù)值模擬;RANS方程;滑移網(wǎng)格

    1 概 述

    近年來,由于特種推進器在提高船舶推進效率,減少艦船振動以及降低噪聲,節(jié)省燃油消耗等方面發(fā)揮了巨大優(yōu)勢,越來越受到國內(nèi)外研究機構(gòu)和學(xué)者的重視,并且創(chuàng)造了不少新的推進形式,相繼進行了一系列的理論和試驗研究[1-18]。吊艙式電力推進系統(tǒng)于1989年由ABB公司首先推出,獲得了巨大成功。在此基礎(chǔ)上,2000年推出了吊艙式CRP系統(tǒng)[15],如圖1所示,前槳就是傳統(tǒng)推進系統(tǒng)的螺旋槳,而在傳統(tǒng)舵葉所在的位置,安裝1個可360°旋轉(zhuǎn)的吊艙式推進器。2個螺旋槳位于同一軸線上,轉(zhuǎn)向相反。同時后面的吊艙式推進器,還起著舵的作用,稱其為主動舵,而傳統(tǒng)的舵則稱為被動舵。吊艙式CRP推進器綜合了吊艙推進器和CRP推進器的優(yōu)點,除了具有很高的推進效率、良好的操作性能和空泡性能、節(jié)能等特點外,還具備系統(tǒng)可靠性高、空間體積小、機械結(jié)構(gòu)相對簡單、機動性和靈活性好的優(yōu)點[16]。

    圖1 吊艙式CRP推進器Fig.1 Podded contra-rotating propulsor

    吊艙式CRP推進器模型具有特殊的結(jié)構(gòu)形式,因此敞水、自航試驗需要開發(fā)專門的試驗設(shè)備,研究相應(yīng)的試驗方法;在各種吊艙外殼阻力的尺度修正方法中,ITTC推薦的是以數(shù)值模擬為基礎(chǔ)的方法,也即吊艙式CRP推進器的試驗測量對數(shù)值模擬形成明顯的借助關(guān)系[19-20]。這是吊艙式CRP推進器與常規(guī)槳試驗測量相比一個明顯的特點。試驗方法和數(shù)值模擬是研究吊艙式CRP推進器水動力性能的重要手段,通過試驗和計算所得到的數(shù)據(jù)是螺旋槳、艙體以及支架的設(shè)計與優(yōu)化、實船性能預(yù)報的基礎(chǔ)。

    本文針對吊艙式CRP推進器建立數(shù)學(xué)模型,結(jié)合RANS方程和SST k-ω湍流模型,運用滑移網(wǎng)格方法對吊艙式CRP推進器在均勻流場中水動力性能進行了非定常數(shù)值預(yù)報,將數(shù)值預(yù)報所得的敞水性能數(shù)據(jù)與在真實空泡水洞內(nèi)利用吊艙動力儀及長軸動力儀對吊艙式CRP推進器進行敞水試驗得到的試驗數(shù)據(jù)進行比較,同時將上述方法得到的葉面及葉背壓力系數(shù)和非定常性能進行分析。

    2 數(shù)值模擬方法

    2.1 RANS方程和湍流模型

    RANS方程是粘性流體運動學(xué)和動力學(xué)的普適性控制方程,是求解吊艙推進器水動力性能計算的基本方程。RANS方程是通過將瞬時N-S方程中的速度、壓強、質(zhì)量力、密度等流體變量進行時歷平均化后得到,其形式如下:

    式中:ρ為流體密度;p為靜壓;fi為單位質(zhì)量的質(zhì)量力;ui和uj為速度分量。雷諾平均N-S方程與瞬時N-S方程在形式上基本一致,只是方程中的變量是時歷平均值。另外,還出現(xiàn)了代表湍流效應(yīng)的雷諾應(yīng)力項-,湍流應(yīng)力項的出現(xiàn)增加了方程的未知數(shù)個數(shù),為了使方程封閉,須引入湍流模式。

    式中α1為特定常數(shù)。

    SST k-ω模型,綜合了k-ω模型在近壁區(qū)計算的優(yōu)點和k-ε模型在遠場計算的優(yōu)點,將k-ω模型和標(biāo)準(zhǔn)k-ε都乘以一個混合函數(shù)后再相加就得到這個模型。在近壁區(qū),混合函數(shù)的值等于1,因此在近壁區(qū)等價于k-ω模型。在遠離壁面的區(qū)域混合函數(shù)的值則等于0,因此自動轉(zhuǎn)換為標(biāo)準(zhǔn)k-ε模型。

    2.2 旋轉(zhuǎn)區(qū)域處理

    滑移網(wǎng)格技術(shù)的基本原理是將幾何模型網(wǎng)格劃分成幾個區(qū)域,交界面兩側(cè)網(wǎng)格相互滑動,而不要求交接面兩側(cè)的網(wǎng)格結(jié)點相互重合,但要計算交界面兩側(cè)的通量,使其相等。為了計算交界面的通量,首先在每一個新的時間步確定出交界面兩邊交界區(qū)的重合面。基本上,通過網(wǎng)格重合面的通量由交界面兩邊交界區(qū)的重合面計算(見圖2)。交界面區(qū)域是由A-B,B-C和D-E,E-F組成(見圖3)。這2個區(qū)域的相交產(chǎn)生d-b,b-e和e-c,2個網(wǎng)格單元區(qū)塊在d-b,b-e和e-c上的重疊構(gòu)成了內(nèi)部區(qū)域。為計算通過單元Ⅲ的通量(D-E上),在計算過程中將不考慮D-E,而是由d-b和b-e來代替,通過d-b和b-e分別由單元Ⅰ和單元Ⅱ把流場信息帶入到單元Ⅲ中。Rhinoceros軟件生成吊艙艙體模型,見圖4,再導(dǎo)入FLUENT前處理軟件GAMBIT導(dǎo)入螺旋槳型值生成吊艙式CRP推進器,如圖5所示。

    圖4 吊艙艙體數(shù)值模型Fig.4 Numberical model of pod

    3 研究對象建模及網(wǎng)格劃分

    3.1 模型主參數(shù)

    吊艙式CRP推進器由吊艙、前槳及后槳等構(gòu)成。吊艙艙體參數(shù)見表1。螺旋槳模型分為前槳和后槳。前槳選用 DTMB3686,后槳選用DTMB3849,兩槳間距為前槳直徑的0.1415倍[17-18],模 型 主 要 參 數(shù) 見 表 2。 建 模 中 采 用

    表1 吊艙模型參數(shù)Tab.1 Main parameters of the pod

    表2 對轉(zhuǎn)槳模型主參數(shù)Tab.2 Main parameters of the CRP propellers

    圖5 吊艙式CRP推進器數(shù)值模型Fig.5 Numberical model of Poded CRP proplsor

    3.2 網(wǎng)格劃分

    為更真實地模擬空泡水洞內(nèi)吊艙式CRP推進器敞水性能以及前后槳和艙體、支架、尾鰭等之間的相互作用,建模中將吊艙式CRP推進器放置于數(shù)值空泡水洞內(nèi),空泡水洞工作段截面尺寸為0.6 m×0.6 m。前后螺旋槳采用滑移網(wǎng)格技術(shù)來實現(xiàn)螺旋槳的旋轉(zhuǎn)效應(yīng),滑移面將整體網(wǎng)格分成固定區(qū)域和動區(qū)域,動區(qū)域之間以及動區(qū)域和固定區(qū)域之間接觸面上的網(wǎng)格盡量一致,接觸面附近區(qū)域的網(wǎng)格要求盡量精細。建立計算域網(wǎng)格模型如圖6所示,總網(wǎng)格數(shù)220萬。固定區(qū)域網(wǎng)格數(shù)為115萬;動區(qū)域分為2部分,都為圓柱體,直徑1.18D,前槳動區(qū)域網(wǎng)格數(shù)約為60萬,后槳動區(qū)域網(wǎng)格數(shù)約為55萬。整個流場區(qū)域網(wǎng)格最差畸變度為0.843 3,如圖6所示。上游邊界距離前槳槳盤約為3倍的前槳直徑,下游邊界距離前槳槳盤約為15倍的前槳直徑,由于葉梢、導(dǎo)邊、隨邊、槳轂上流動變化比較劇烈,對上述區(qū)域的網(wǎng)格進行了加密,如圖7和圖8所示。

    4 水動力性能分析

    4.1 水動力性能數(shù)據(jù)

    對于吊艙式CRP推進器需要計算:前槳推力TF,吊艙推進器推力TP,前槳扭矩QF,后槳推力TA,后槳扭矩QA,螺旋槳轉(zhuǎn)速n,水流進速VA。

    為了便于比較分析,通常均以前槳直徑無因次化。

    4.2 敞水性能分析

    利用吊艙動力儀和長軸動力儀針對吊艙式CRP推進器在海軍工程大學(xué)空泡水洞中利用吊艙動力儀和長軸動力儀在相同工況下進行了3組重復(fù)性敞水試驗,每組試驗結(jié)果之間吻合很好,將3組試驗結(jié)果在對應(yīng)進速系數(shù)下進行平均,如表3所示。

    表3 吊艙式CRP推進器敞水性能試驗數(shù)據(jù)Tab.3 The experimental data of poded contra-rotating propulsor

    在數(shù)值空泡水洞內(nèi)對吊艙式CRP推進器數(shù)值模型敞水性能進行了計算,敞水性能數(shù)值預(yù)報結(jié)果見表4。試驗值和計算值之間的誤差見表5。

    表4 吊艙式CRP推進器敞水性能數(shù)值預(yù)報結(jié)果Tab.4 Numberical predicted results of poded contra-rotating propulsor

    表5 試驗值和計算值之間的誤差Tab.5 Comparison of experimental results and predicted results

    從結(jié)果可知,前槳的誤差要比后槳的誤差總體來說偏低,計算值與試驗數(shù)據(jù)吻合較好,最大誤差為7.81%,絕大部分誤差控制在6%以內(nèi),說明本文所用數(shù)值計算方法適用于吊艙式CRP推進器敞水性能預(yù)報,且都得到較好的預(yù)報精度和穩(wěn)定性,滿足工程應(yīng)用的要求。

    4.3 槳葉表面壓力分析

    本文運用滑移網(wǎng)格方法對吊艙式CRP推進器每個葉片的不同半徑上的壓力分布進行了預(yù)報,并與不附帶吊艙相同對轉(zhuǎn)槳壓力分布計算結(jié)果進行了比較,圖9和圖10分別為吊艙式CRP推進器和CRP推進器在進速系數(shù)J=1.1時,r/R=0.3,0.7,0.9這3個不同半徑上各個槳葉壓力系數(shù)的分布。壓力系數(shù)計算公式如下:

    式中:(p-p0)為相對壓力;VA為遠前方來流速度。

    從圖中可以看出,不論吊艙式CRP推進器還是CRP推進器,前后槳葉數(shù)不同,而對轉(zhuǎn)速度一致,前槳與后槳槳葉之間的相位產(chǎn)生周期性變化,從而導(dǎo)致每個槳葉上壓力系數(shù)分布不同。但是由于吊艙的阻塞效應(yīng),吊艙式CRP推進器前后槳不同槳葉上壓力分布的差別要比CRP推進器的差別明顯,尤其是距離吊艙較近的后槳;吊艙式CRP推進器前槳在對應(yīng)半徑上葉背與葉面壓力差值(即升力)和CRP推進器計算結(jié)果差別不大,而后槳在對應(yīng)半徑上葉背與葉面壓力差值相比CRP推進器計算結(jié)果偏大。

    4.4 非定常性能數(shù)值分析

    圖9 吊艙式CRP推進器在J=1.1時前后槳葉上壓力系數(shù)的分布Fig.9 Pressure coefficient distribution of poded contra-rotating propulsor

    圖10 CRP推進器在J=1.1時前后槳葉上壓力系數(shù)的分布Fig.10 Pressure coefficient distribution of contra-rotating propellers

    圖11 吊艙式CRP推進器前槳推力系數(shù)和扭矩系數(shù)變化Fig.11 Variation of thrust and moment coefficient of poded contra-rotating propulsor’s forward propeller

    本文運用滑移網(wǎng)格方法對吊艙式CRP推進器非定常性能進行了預(yù)報,圖11~圖13是吊艙式CRP推進器在進速系數(shù)J=1.1時,前槳與后槳相對旋轉(zhuǎn)時推力系數(shù)和扭矩系數(shù)以及吊艙推力系數(shù)的變化規(guī)律。本文計算了前后槳及吊艙在一個轉(zhuǎn)動周期(360°)推力系數(shù)和扭矩系數(shù)的變化。同時與不附帶吊艙相同對轉(zhuǎn)槳非定常性能計算數(shù)據(jù)進行了比較,圖14和圖15為CRP推進器前后槳推力系數(shù)和扭矩系數(shù)在單個葉片轉(zhuǎn)動周期內(nèi) (前槳90°,后槳72°)推力系數(shù)和扭矩系數(shù)變化規(guī)律。

    由圖11~圖15可以看出,對于吊艙式CRP推進器來說,前后槳在單個槳葉轉(zhuǎn)動周期內(nèi)推力系數(shù)和扭矩系數(shù)成周期性變化,且前槳和后槳推力系數(shù)與扭矩系數(shù)在單個槳葉轉(zhuǎn)動周期內(nèi)變化頻率和不附帶吊艙相同對轉(zhuǎn)槳一樣,分別為軸頻10倍和8倍,由此可得不論有無吊艙,前后槳轉(zhuǎn)動1周時推力系數(shù)和扭矩系數(shù)變化頻率都為軸頻40倍;但由于吊艙的影響,在單個槳葉轉(zhuǎn)動周期內(nèi)這種變化相較于CRP推進器來說,是不規(guī)律的,但在1個轉(zhuǎn)動周期內(nèi) (360°)總體來看,推力系數(shù)和扭矩系數(shù)的變化是有一定規(guī)律的,其峰值變化頻率和吊艙推力系數(shù)變化頻率一致,都為軸頻的9倍。

    5 結(jié) 語

    本文建立吊艙式CRP推進器數(shù)學(xué)模型,結(jié)合RANS方程和SST湍流模型,運用滑移網(wǎng)格方法對吊艙式CRP推進器在均勻流場中水動力性能進行了非定常數(shù)值預(yù)報。結(jié)果表明:

    1)本文所用數(shù)值計算方法適用于吊艙式CRP推進器敞水性能預(yù)報,且得到較好的預(yù)報精度和穩(wěn)定性,滿足工程應(yīng)用的要求。

    2)由于吊艙的阻塞效應(yīng),吊艙式CRP推進器前后槳不同槳葉上壓力分布的差別要比CRP推進器的差別明顯,尤其是距離吊艙較近的后槳;吊艙式CRP推進器前槳在對應(yīng)半徑上葉背與葉面壓力差值(即升力)和CRP推進器計算結(jié)果差別不大,而后槳在對應(yīng)半徑上葉背與葉面壓力差值相比CRP推進器計算結(jié)果偏大。

    3)不論有無吊艙,前后槳轉(zhuǎn)動一周時推力系數(shù)和扭矩系數(shù)變化頻率都為軸頻40倍,但由于吊艙的影響,在單個槳葉轉(zhuǎn)動周期內(nèi)吊艙式CRP推進器前后槳推力系數(shù)和扭矩系數(shù)變化相較于CRP推進器是不規(guī)律的,但在1個轉(zhuǎn)動周期內(nèi)總體來看,吊艙式CRP推進器前后槳推力系數(shù)和扭矩系數(shù)的變化是有一定規(guī)律的 ,其峰值變化頻率和吊艙推力系數(shù)變化頻率一致,都為軸頻的9倍。

    [1]熊鷹,葉金銘.吊艙電力推進系統(tǒng)性能評估及設(shè)計方法[J].海軍工程大學(xué)學(xué)報,2002,14(1):23-26.

    XIONG Ying,YE Jin-ming.Evalution of poded eletric propulsion system performance and design method of poded eletric propulsion [J].Journal of Naval University of Engineering,2002,14(1):23-26.

    [2]The Propulsion Committee.Final report and recommendations to the 22nd ITTC[R].Proceeding of 22nd ITTC,1999.

    [3]Poded propulsor tests and extrapodation[R].ITTC Recommended Procedures,Proceedings of the 23nd ITTC,2002.

    [4]The specialist committee on azimuthing poded propulsion[R].Final Report and Recommendations to the 24th ITTC,2005,2.

    [5]The Propulsion Committee[R].Final Report and Recommendations to the 24th ITTC,Proceedings of the 24th ITTC,2005,1.

    [6]ISLAM M,TAYLOR R, QUINTON J.Numerical investigation of propulsive characteristics of podded propel1ers[C].The 1st Intemational Conference on Technological Advances in Podded Propulsion.Newcastle:The University of Newcastle,2004.

    [7]SAKIR B,HAKAN A,MESUT G.Preliminary results of a numerical method for podded propulsors[C].TPOD 2006,October 3-5,L`Aber-Wrac`h,F(xiàn)rance.

    [8]ISLAM M F,MOLLOY S,HE M,et al.Hydrodynamic study of podded propulsors with systematically varied geometry[C].T-POD 2006,October 3-5,L'Aber-Wrac'h,F(xiàn)rance.

    [9]ANTONIO S C,JAAKKO V.P.RANS predictions for flow patterns around a compact azipod[C].T-POD 2006,October 3-5,L'Aber-Wrac'h,F(xiàn)rance.

    [10]VAN R M,HOLTROP J.Investigations on a pod open water test set-up[C].1st T-POD,University of Newcastle,UK,2004.

    [11]HEINKE H J.Investigation about the forces and moments at poded drives[C].1st T-POD,University of Newcastle,UK,2004.

    [12]NAKATAKE K,ANDO J,YOSHITAKE A,SATO Y,TAMASHIMA M.On propulsive performance of a small bulk-carrier model with twin poded propellers[C].1st TPOD,University of Newcastle,UK,2004.

    [13]楊晨俊,錢正芳,馬騁.吊艙對螺旋槳水動力性能的影響[J].上海交通大學(xué)學(xué)報,2003,37(8):1229-1233.

    YANG Chen-jun,QIAN Zheng-fang,MA Cheng.Influences of pod on the propeller performance[C].Journal of Shanghai Jiaotong University,2003,37(8):1229-1233.

    [14]熊鷹,葉金銘.用面元法預(yù)報吊艙推進器水動力性能[J].2007,19(2):41-45.

    XIONG Ying,YE Jin-ming.Prediction of poded propulsion hydrodynamic performance by surface-panel method[J].Journal of Naval University of Engineering,2007,19(2):41-45.

    [15]OSKAR L.Advanced machinery with CRP propulsion for fast RoPax vessels[C].Proceedings of the 24th Motorship Marine Propulsion Conference.Copenhagen.2002.

    [16]聶延生,韓學(xué)勝,等.對轉(zhuǎn)螺旋槳的結(jié)構(gòu)原理及特點分析[J].船電技術(shù),2005,(2):50-52.

    NIE Yan-sheng,HAN Xue-sheng,et al.Analysis about the construction principle and characteristics of contra rotating propeller[J].Marine Electric & Electronic Engineering,2005,(2):50-52.

    [17]MARLIN L M.Experimental determination of unsteady forces on contrarotating propellers in uniform flow[R].DAVID W TAYLOR Naval Ship Research and Development Center U.S.1976.

    [18]劉小龍,唐登海,侯櫻.對轉(zhuǎn)槳定常面元法水動力性能預(yù)估[J].中國造船,2009,50(3):1-8.

    LIU Xiao-long,TANG Deng-h(huán)ai,HOU Ying.Prediction of steady performance of contra-rotating propellers by potential based panel method[J].Shipbuilding of China,2009,50(3):1-8.

    [19]Podded propulsor tests and extrapodation[R].ITTC Recommended Procedures,Proceedings of the 23nd ITTC,2002.

    [20]The specialist committee on azimuthing podded propulsion[R].Final Report and Recommendations to the 24th ITTC,2005,2.

    Experimental investigation and numberical simulation on the open-water performance of poded contra-rotating propulsor

    SHENG Li,XIONG Ying,YANG Yong
    (Department of Naval Architecture and Ocean Engineering,Naval Unirevsity of Engineering,Wuhan 430033,China)

    Numberical model of poded contra-rotating propulsor is estabished,the unsteady hydrodynamics performance of the poded contra-rotating propulsor in uniform flow are predicted by the RANS formula with SST k-ωturbulence model based on sliding mesh method;the prediction results of the thrust coefficients,torque coefficients are compared with the experimental value of the same poded contra-rotating propulsor in the real cavitation tunnel by pod dynamical instrument and long-axis dynamical instrument.And then analysis the unsteady performance and the pressure coefficients of face side and back side of blade.It is shown that,the numerical method presented in this paper has good precision in the prediction of hydrodynamics performance of poded contrarotating propulsor,can achieve the requirements of engineering application.

    poded contra-rotating propulsor;hydrodynamics performance;experimental investigations;numberical simulation;RANS formula;sliding mesh

    TB535

    A

    1672-7649(2012)10-0009-08

    10.3404/j.issn.1672-7649.2012.10.002

    2011-09-07

    盛立(1984-),男,博士,工程師,研究方向為船舶水動力學(xué)。

    猜你喜歡
    吊艙推進器槳葉
    探究奇偶旋翼對雷達回波的影響
    基于有限元的吊艙內(nèi)部溫度場仿真分析與優(yōu)化
    戰(zhàn)機吊艙火力賦能“神器”
    基于CFD扇翼推進器敞水性能預(yù)報分析
    立式捏合機槳葉結(jié)構(gòu)與槳葉變形量的CFD仿真*
    發(fā)揮考核“指揮棒”“推進器”作用
    NEO—2X光電吊艙
    航空模型(2016年10期)2017-05-09 06:30:58
    讓黨建成為信仰播種機和工作推進器
    直升機槳葉/吸振器系統(tǒng)的組合共振研究
    吊艙式電力推進船舶螺旋槳匹配設(shè)計仿真研究
    色婷婷av一区二区三区视频| 亚洲精品一二三| 欧美日韩视频高清一区二区三区二| 亚洲av国产av综合av卡| 一级毛片黄色毛片免费观看视频| 亚洲国产精品一区二区三区在线| 亚洲精品乱码久久久久久按摩| 久久久亚洲精品成人影院| 一级毛片黄色毛片免费观看视频| 久久人人97超碰香蕉20202| 一边摸一边做爽爽视频免费| 咕卡用的链子| 国产成人a∨麻豆精品| 免费av中文字幕在线| 国产精品三级大全| 激情五月婷婷亚洲| 亚洲欧洲日产国产| 亚洲美女视频黄频| 国产精品久久久久久久电影| 国产不卡av网站在线观看| 国产亚洲av片在线观看秒播厂| 超碰97精品在线观看| 黄网站色视频无遮挡免费观看| 天堂俺去俺来也www色官网| 国产乱人偷精品视频| 两个人免费观看高清视频| 少妇高潮的动态图| 天天躁夜夜躁狠狠久久av| 国产极品粉嫩免费观看在线| 搡老乐熟女国产| videossex国产| 秋霞在线观看毛片| 国产精品麻豆人妻色哟哟久久| 国产伦理片在线播放av一区| 久久国产精品大桥未久av| 亚洲激情五月婷婷啪啪| av又黄又爽大尺度在线免费看| 免费av中文字幕在线| 美女福利国产在线| 欧美3d第一页| a级片在线免费高清观看视频| 久久国产精品大桥未久av| 欧美日韩成人在线一区二区| 亚洲av福利一区| 久久青草综合色| 下体分泌物呈黄色| 久久精品国产a三级三级三级| 80岁老熟妇乱子伦牲交| 寂寞人妻少妇视频99o| 中文欧美无线码| 大片电影免费在线观看免费| 美女内射精品一级片tv| 久久人人97超碰香蕉20202| 99re6热这里在线精品视频| 国产在线视频一区二区| 久久人妻熟女aⅴ| 日本午夜av视频| 午夜激情久久久久久久| 亚洲精品av麻豆狂野| 欧美 亚洲 国产 日韩一| 国产精品一二三区在线看| 免费黄频网站在线观看国产| 亚洲中文av在线| av网站免费在线观看视频| 国产熟女午夜一区二区三区| 精品第一国产精品| 狂野欧美激情性bbbbbb| 一区二区av电影网| 最近中文字幕高清免费大全6| 在现免费观看毛片| 老女人水多毛片| 亚洲在久久综合| 国产成人精品无人区| 久久这里只有精品19| 久久热在线av| 久久久久国产精品人妻一区二区| av在线app专区| 宅男免费午夜| 日韩一本色道免费dvd| 日韩中文字幕视频在线看片| 精品少妇内射三级| 国产精品国产三级国产专区5o| 亚洲国产毛片av蜜桃av| 少妇高潮的动态图| 日本午夜av视频| 伦理电影免费视频| av卡一久久| 亚洲欧美一区二区三区国产| 青春草视频在线免费观看| 你懂的网址亚洲精品在线观看| 制服丝袜香蕉在线| 精品午夜福利在线看| 美女内射精品一级片tv| 97在线人人人人妻| 国产探花极品一区二区| 1024视频免费在线观看| 欧美另类一区| 国产精品一区www在线观看| 看非洲黑人一级黄片| 中文欧美无线码| 精品国产乱码久久久久久小说| 波多野结衣一区麻豆| 亚洲熟女精品中文字幕| 精品国产一区二区三区久久久樱花| 国产成人精品在线电影| 欧美精品人与动牲交sv欧美| 这个男人来自地球电影免费观看 | 欧美性感艳星| 一级毛片电影观看| 搡老乐熟女国产| 校园人妻丝袜中文字幕| av国产久精品久网站免费入址| 欧美丝袜亚洲另类| a级片在线免费高清观看视频| 亚洲色图 男人天堂 中文字幕 | 秋霞在线观看毛片| 成年女人在线观看亚洲视频| 日本色播在线视频| 91精品国产国语对白视频| 丝袜在线中文字幕| 国产一区有黄有色的免费视频| 狠狠婷婷综合久久久久久88av| 久久久久久伊人网av| 大片免费播放器 马上看| 国产老妇伦熟女老妇高清| 久久这里有精品视频免费| 大话2 男鬼变身卡| 男女边吃奶边做爰视频| 麻豆乱淫一区二区| 欧美最新免费一区二区三区| 亚洲欧美成人综合另类久久久| 精品国产国语对白av| 日韩视频在线欧美| 大话2 男鬼变身卡| 一区二区三区四区激情视频| 久热久热在线精品观看| 插逼视频在线观看| 国产精品一区二区在线不卡| 成人二区视频| 女性生殖器流出的白浆| 99热国产这里只有精品6| 一区二区三区四区激情视频| 国产亚洲av片在线观看秒播厂| 欧美日韩一区二区视频在线观看视频在线| 最新的欧美精品一区二区| 亚洲国产av新网站| 性色avwww在线观看| 国产69精品久久久久777片| 新久久久久国产一级毛片| 女性生殖器流出的白浆| 免费不卡的大黄色大毛片视频在线观看| 国产av国产精品国产| 国产一级毛片在线| 欧美激情国产日韩精品一区| 极品人妻少妇av视频| 啦啦啦视频在线资源免费观看| 久热久热在线精品观看| 亚洲人成77777在线视频| 久久精品久久久久久噜噜老黄| av播播在线观看一区| 午夜激情久久久久久久| 热re99久久精品国产66热6| 国产精品不卡视频一区二区| 成人黄色视频免费在线看| 中国国产av一级| 一边摸一边做爽爽视频免费| www.av在线官网国产| 成年美女黄网站色视频大全免费| 亚洲综合色惰| 高清不卡的av网站| 日日啪夜夜爽| 精品熟女少妇av免费看| xxx大片免费视频| av电影中文网址| 国产亚洲av片在线观看秒播厂| 日本-黄色视频高清免费观看| 成人国语在线视频| 欧美 亚洲 国产 日韩一| 久久国产精品男人的天堂亚洲 | 国产精品成人在线| 午夜久久久在线观看| 欧美另类一区| 熟女电影av网| 丰满乱子伦码专区| 久久久久国产精品人妻一区二区| 久久久久精品人妻al黑| 美女国产视频在线观看| 国产1区2区3区精品| 久久久久久久亚洲中文字幕| 亚洲精品乱码久久久久久按摩| 久久ye,这里只有精品| 欧美3d第一页| 国产白丝娇喘喷水9色精品| 精品国产一区二区久久| 一区二区三区乱码不卡18| 啦啦啦啦在线视频资源| 宅男免费午夜| 纵有疾风起免费观看全集完整版| 99热国产这里只有精品6| 咕卡用的链子| 人妻 亚洲 视频| 国产极品天堂在线| 国产免费视频播放在线视频| 久久精品久久久久久久性| 少妇人妻精品综合一区二区| 免费看不卡的av| 制服诱惑二区| 亚洲,一卡二卡三卡| 亚洲av综合色区一区| videosex国产| 国产成人免费观看mmmm| 精品少妇久久久久久888优播| 国产精品久久久久久精品电影小说| 欧美97在线视频| 成人毛片60女人毛片免费| 午夜久久久在线观看| 中文字幕精品免费在线观看视频 | 免费观看在线日韩| 各种免费的搞黄视频| 一区二区三区四区激情视频| tube8黄色片| 久久久国产精品麻豆| 丝袜人妻中文字幕| av免费观看日本| 一区二区三区精品91| 国产一区二区在线观看av| 免费久久久久久久精品成人欧美视频 | 国产av一区二区精品久久| 香蕉国产在线看| 国产成人免费无遮挡视频| 日韩精品免费视频一区二区三区 | 99久久综合免费| 国产亚洲av片在线观看秒播厂| 亚洲av欧美aⅴ国产| 成人漫画全彩无遮挡| 免费黄网站久久成人精品| 国产免费一级a男人的天堂| 丰满迷人的少妇在线观看| 男人舔女人的私密视频| 国产成人精品无人区| 黄片播放在线免费| 伊人久久国产一区二区| 午夜日本视频在线| 91精品国产国语对白视频| 乱码一卡2卡4卡精品| 91国产中文字幕| 久久精品国产综合久久久 | 亚洲av福利一区| 综合色丁香网| 国产精品麻豆人妻色哟哟久久| 母亲3免费完整高清在线观看 | a级片在线免费高清观看视频| 满18在线观看网站| 一区二区三区四区激情视频| 男人添女人高潮全过程视频| 久久久久精品久久久久真实原创| 黄色怎么调成土黄色| 亚洲av在线观看美女高潮| 国产高清三级在线| 久久99精品国语久久久| 精品卡一卡二卡四卡免费| 天天影视国产精品| 香蕉国产在线看| 高清欧美精品videossex| 国产又色又爽无遮挡免| 我的女老师完整版在线观看| 亚洲av综合色区一区| 国产av码专区亚洲av| 欧美少妇被猛烈插入视频| 日日爽夜夜爽网站| 日韩不卡一区二区三区视频在线| 亚洲成人av在线免费| 亚洲第一区二区三区不卡| 国产一区二区三区av在线| 2018国产大陆天天弄谢| 久久久久久人人人人人| 国产成人精品久久久久久| 又大又黄又爽视频免费| 中国国产av一级| 亚洲精品一区蜜桃| 天天操日日干夜夜撸| 色视频在线一区二区三区| 免费少妇av软件| 母亲3免费完整高清在线观看 | 欧美激情 高清一区二区三区| 国产精品一国产av| 好男人视频免费观看在线| 亚洲精品一区蜜桃| 国产免费现黄频在线看| 中文字幕最新亚洲高清| av在线app专区| 91午夜精品亚洲一区二区三区| 久久97久久精品| 男女无遮挡免费网站观看| 亚洲天堂av无毛| 欧美人与性动交α欧美软件 | 精品一区二区三区四区五区乱码 | 亚洲成国产人片在线观看| 亚洲精品第二区| 80岁老熟妇乱子伦牲交| 亚洲av.av天堂| 建设人人有责人人尽责人人享有的| 亚洲人成77777在线视频| 一本色道久久久久久精品综合| 最近2019中文字幕mv第一页| 久久久久久久久久人人人人人人| 久久久久人妻精品一区果冻| 国产又爽黄色视频| 男人操女人黄网站| 在线观看一区二区三区激情| 免费在线观看黄色视频的| 成年人午夜在线观看视频| 成年av动漫网址| 亚洲精品视频女| 校园人妻丝袜中文字幕| 亚洲丝袜综合中文字幕| 18禁裸乳无遮挡动漫免费视频| 欧美丝袜亚洲另类| 精品人妻熟女毛片av久久网站| 免费不卡的大黄色大毛片视频在线观看| 黄网站色视频无遮挡免费观看| 插逼视频在线观看| 国产av国产精品国产| 久热久热在线精品观看| 少妇的逼水好多| 国产麻豆69| 国产精品国产三级国产专区5o| 老司机影院毛片| 国产精品一区二区在线不卡| 亚洲国产看品久久| 欧美精品一区二区免费开放| 国产黄色视频一区二区在线观看| 18禁国产床啪视频网站| 精品久久久久久电影网| a级毛色黄片| 91精品三级在线观看| 中文精品一卡2卡3卡4更新| 中文字幕亚洲精品专区| 汤姆久久久久久久影院中文字幕| 国产在视频线精品| 九草在线视频观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲情色 制服丝袜| 精品国产一区二区三区四区第35| 日韩av在线免费看完整版不卡| 日本av手机在线免费观看| 久热这里只有精品99| 男的添女的下面高潮视频| 一二三四中文在线观看免费高清| 免费看不卡的av| 亚洲精品aⅴ在线观看| 国产黄色免费在线视频| 日韩制服骚丝袜av| 久久女婷五月综合色啪小说| 久久精品国产亚洲av天美| 国产精品久久久久久精品电影小说| 亚洲精品一区蜜桃| 99国产精品免费福利视频| 18在线观看网站| xxxhd国产人妻xxx| 亚洲国产欧美日韩在线播放| 看免费成人av毛片| 伦精品一区二区三区| 国产在线免费精品| 插逼视频在线观看| 多毛熟女@视频| 免费看av在线观看网站| 另类亚洲欧美激情| 90打野战视频偷拍视频| 国国产精品蜜臀av免费| 亚洲欧美色中文字幕在线| 久久久久精品久久久久真实原创| 热re99久久精品国产66热6| 欧美精品人与动牲交sv欧美| 国产黄色免费在线视频| 亚洲av.av天堂| 亚洲精品456在线播放app| 久久久亚洲精品成人影院| 日韩欧美精品免费久久| 免费在线观看完整版高清| 精品久久国产蜜桃| 美女中出高潮动态图| 婷婷色综合大香蕉| 妹子高潮喷水视频| 在线观看免费日韩欧美大片| 97精品久久久久久久久久精品| 国产一区二区三区综合在线观看 | 美女国产高潮福利片在线看| 99视频精品全部免费 在线| 亚洲成国产人片在线观看| 亚洲av综合色区一区| 欧美xxⅹ黑人| 又黄又粗又硬又大视频| 久久女婷五月综合色啪小说| 91久久精品国产一区二区三区| 丝袜喷水一区| 精品一区二区免费观看| 18禁动态无遮挡网站| 99视频精品全部免费 在线| 成人国产麻豆网| av女优亚洲男人天堂| 在线天堂最新版资源| 久久韩国三级中文字幕| 亚洲,欧美,日韩| 黑丝袜美女国产一区| 中文字幕最新亚洲高清| 国产又色又爽无遮挡免| 精品久久国产蜜桃| 精品一品国产午夜福利视频| 十分钟在线观看高清视频www| 中文字幕免费在线视频6| 女性生殖器流出的白浆| 一级a做视频免费观看| 亚洲第一av免费看| 女人精品久久久久毛片| 女性被躁到高潮视频| 日本免费在线观看一区| 又黄又粗又硬又大视频| 免费观看av网站的网址| 国精品久久久久久国模美| 午夜福利视频在线观看免费| 大话2 男鬼变身卡| 高清欧美精品videossex| 亚洲av电影在线进入| 久久99热这里只频精品6学生| 尾随美女入室| 久久午夜综合久久蜜桃| 亚洲国产精品999| 街头女战士在线观看网站| 少妇人妻 视频| av一本久久久久| 美女中出高潮动态图| 色婷婷av一区二区三区视频| 精品人妻一区二区三区麻豆| 国产成人免费观看mmmm| 亚洲欧美中文字幕日韩二区| 制服丝袜香蕉在线| 国产精品久久久久成人av| 精品国产露脸久久av麻豆| 国产免费一区二区三区四区乱码| 亚洲成人av在线免费| 成人无遮挡网站| 最近中文字幕2019免费版| 91在线精品国自产拍蜜月| 97人妻天天添夜夜摸| 又黄又爽又刺激的免费视频.| 国产极品粉嫩免费观看在线| 狠狠精品人妻久久久久久综合| 大片免费播放器 马上看| 午夜福利乱码中文字幕| 精品国产一区二区三区四区第35| 国产一区二区三区综合在线观看 | 久久久久久久久久成人| 一区二区av电影网| 精品一区在线观看国产| 黑人猛操日本美女一级片| 赤兔流量卡办理| av一本久久久久| 国产午夜精品一二区理论片| 亚洲一级一片aⅴ在线观看| 精品99又大又爽又粗少妇毛片| 日本色播在线视频| 国产欧美日韩一区二区三区在线| 亚洲成色77777| 9191精品国产免费久久| 欧美老熟妇乱子伦牲交| 亚洲国产精品一区三区| 国产又色又爽无遮挡免| 欧美日韩视频高清一区二区三区二| 久久av网站| 伦理电影免费视频| 一级爰片在线观看| 精品久久国产蜜桃| av黄色大香蕉| 各种免费的搞黄视频| 欧美丝袜亚洲另类| 这个男人来自地球电影免费观看 | 日韩一区二区三区影片| av播播在线观看一区| 久久99热6这里只有精品| 香蕉国产在线看| 亚洲精品第二区| 国产熟女欧美一区二区| 美女内射精品一级片tv| www.av在线官网国产| 大香蕉97超碰在线| 人成视频在线观看免费观看| 午夜老司机福利剧场| 熟女电影av网| 高清欧美精品videossex| 亚洲av中文av极速乱| 一二三四在线观看免费中文在 | 中文字幕人妻丝袜制服| 亚洲国产色片| 午夜免费观看性视频| 91精品三级在线观看| 国产精品无大码| 美国免费a级毛片| 国产1区2区3区精品| 黄色 视频免费看| 国产免费现黄频在线看| 亚洲av.av天堂| 99久国产av精品国产电影| 亚洲人成网站在线观看播放| 美女主播在线视频| 国产黄色视频一区二区在线观看| 不卡视频在线观看欧美| 日韩中文字幕视频在线看片| 亚洲精品国产av蜜桃| 精品一品国产午夜福利视频| 18在线观看网站| 亚洲四区av| 精品酒店卫生间| 26uuu在线亚洲综合色| 久久久久久人人人人人| 少妇人妻久久综合中文| 美女xxoo啪啪120秒动态图| 精品人妻在线不人妻| 国产一区二区在线观看日韩| 亚洲精品一区蜜桃| 亚洲国产欧美在线一区| h视频一区二区三区| 亚洲一级一片aⅴ在线观看| 国产精品久久久av美女十八| 国产激情久久老熟女| 熟女人妻精品中文字幕| 亚洲av日韩在线播放| 曰老女人黄片| 一级,二级,三级黄色视频| 欧美成人午夜精品| 国产极品粉嫩免费观看在线| 久久狼人影院| 成人影院久久| 精品少妇黑人巨大在线播放| 我要看黄色一级片免费的| 国产成人精品久久久久久| 少妇 在线观看| 婷婷色麻豆天堂久久| 97在线人人人人妻| 一级a做视频免费观看| 欧美激情 高清一区二区三区| 人人妻人人添人人爽欧美一区卜| 在线天堂最新版资源| 黄色视频在线播放观看不卡| 亚洲色图 男人天堂 中文字幕 | av在线播放精品| 亚洲国产精品国产精品| 中文精品一卡2卡3卡4更新| 国产成人午夜福利电影在线观看| 黑人高潮一二区| 成人午夜精彩视频在线观看| 日日撸夜夜添| 久久精品国产a三级三级三级| 五月玫瑰六月丁香| 国产精品久久久久久久久免| 99久久中文字幕三级久久日本| 精品一区二区三卡| 成人漫画全彩无遮挡| 岛国毛片在线播放| 欧美 亚洲 国产 日韩一| 男女下面插进去视频免费观看 | videosex国产| www日本在线高清视频| 国产精品成人在线| 久久精品国产亚洲av天美| 欧美另类一区| av有码第一页| 秋霞在线观看毛片| 日韩精品有码人妻一区| 老熟女久久久| 最近最新中文字幕免费大全7| 国产在线一区二区三区精| 少妇的逼好多水| 色5月婷婷丁香| 青春草亚洲视频在线观看| 亚洲国产日韩一区二区| 丰满少妇做爰视频| 高清黄色对白视频在线免费看| 黑人欧美特级aaaaaa片| 午夜激情av网站| av免费观看日本| 成人二区视频| 亚洲av.av天堂| 久热这里只有精品99| 黄色怎么调成土黄色| 九九在线视频观看精品| 日韩一区二区三区影片| 极品人妻少妇av视频| 精品人妻在线不人妻| 汤姆久久久久久久影院中文字幕| 免费av中文字幕在线| 热99久久久久精品小说推荐| 九色成人免费人妻av| 久久精品国产亚洲av涩爱| 少妇猛男粗大的猛烈进出视频| 男女免费视频国产| 国产精品一二三区在线看| 精品久久久久久电影网| 最后的刺客免费高清国语| av免费观看日本| 免费少妇av软件| 视频中文字幕在线观看| 90打野战视频偷拍视频| 51国产日韩欧美| 日韩中字成人| 成人二区视频| 精品福利永久在线观看| 天堂俺去俺来也www色官网| av有码第一页| 成人无遮挡网站| 亚洲国产毛片av蜜桃av| 如日韩欧美国产精品一区二区三区| 亚洲国产精品国产精品| 男女高潮啪啪啪动态图| 国产成人一区二区在线|