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

    實尺度噴水推進(jìn)船拖泵工況數(shù)值模擬與分析

    2016-05-04 18:39:20靳栓寶魏應(yīng)三胡鵬飛
    船舶力學(xué) 2016年11期
    關(guān)鍵詞:推進(jìn)器船體流場

    靳栓寶,沈 洋,王 東,劉 濤,魏應(yīng)三,胡鵬飛,祝 昊

    (海軍工程大學(xué) 艦船綜合電力技術(shù)國防科技重點實驗室,武漢430033)

    實尺度噴水推進(jìn)船拖泵工況數(shù)值模擬與分析

    靳栓寶,沈 洋,王 東,劉 濤,魏應(yīng)三,胡鵬飛,祝 昊

    (海軍工程大學(xué) 艦船綜合電力技術(shù)國防科技重點實驗室,武漢430033)

    對于四泵推進(jìn)的噴水推進(jìn)船,在巡航工況時中間加速泵通常處于鎖軸狀態(tài),其拖曳阻力的大小對噴水推進(jìn)器的選型以及船泵機(jī)的最優(yōu)匹配有著重要影響。然而,拖泵阻力很難通過船模試驗的方法獲得。為此,該研究在驗證均勻和非均勻條件噴水推進(jìn)器數(shù)值模型的準(zhǔn)確性基礎(chǔ)上,采用數(shù)值試驗的方法對18節(jié)航速下某雙泵推進(jìn)噴水推進(jìn)船的實尺度“船體+兩臺噴水推進(jìn)器”系統(tǒng)帶自由液面的流場進(jìn)行了數(shù)值模擬,計算此時噴泵拖曳阻力及其所占船體阻力的百分比。以此噴泵拖曳阻力作為參考,對尺寸與上述噴水推進(jìn)泵相近的某四泵推進(jìn)噴水推進(jìn)船的噴泵進(jìn)行了選型和設(shè)計,并對該船在18節(jié)航速下加速泵拖曳阻力的大小進(jìn)行了計算,進(jìn)一步驗證選型時拖曳阻力取值的合理性。為消除尺度效應(yīng)的影響采用實尺度模型對“船體+四臺噴水推進(jìn)器”系統(tǒng)帶自由液面的非定常流場進(jìn)行計算,并探索了大尺度條件下船泵系統(tǒng)考慮自由液面和重力影響的非定常計算方法。

    船舶;噴水推進(jìn)器;實尺度;拖泵阻力;數(shù)值模擬;非定常

    0 引 言

    噴水推進(jìn)技術(shù)應(yīng)用于高速高性能艦船已成為國內(nèi)外公認(rèn)的發(fā)展趨勢,其應(yīng)用范圍正不斷地向大功率、大噸位艦艇延伸[1]。與常規(guī)螺旋槳推進(jìn)相比,噴水推進(jìn)具有機(jī)動性和操縱性優(yōu)越、高航速時具有極佳的抗空泡能力和更高的推進(jìn)效率、水動力噪聲小等眾多優(yōu)點,在世界海軍水面艦船上有著廣泛的應(yīng)用[2-3]。船舶噴水推進(jìn)器設(shè)計前需要根據(jù)噴水推進(jìn)基本原理以及船體阻力等參數(shù)對噴水推進(jìn)器主要設(shè)計參數(shù)進(jìn)行優(yōu)化選取,使“船體—噴水推進(jìn)器—原動機(jī)”三者之間達(dá)到最佳匹配。另外,對于多泵推進(jìn)的噴水推進(jìn)船,在巡航工況下加速泵處于鎖軸工況,其將產(chǎn)生額外的阻力,該阻力連同船體阻力構(gòu)成了整個船泵系統(tǒng)在巡航工況下的阻力,巡航所用噴泵的選型需要依據(jù)總阻力值進(jìn)行主要參數(shù)的確定。然而,國外有關(guān)拖泵阻力計算鮮有報道,國內(nèi)關(guān)于多泵推進(jìn)噴水推進(jìn)船舶的設(shè)計處于起步階段,拖泵阻力的經(jīng)驗數(shù)據(jù)不足。

    第21屆至24屆國際拖曳水池會議(ITTC)成立的噴水推進(jìn)專家委員會指出,隨著高性能計算機(jī)的出現(xiàn)和數(shù)值計算方法的改進(jìn),CFD技術(shù)在噴水推進(jìn)研究中的應(yīng)用將越來越廣泛[4-7]。將噴水推進(jìn)器與整個船體作為整體進(jìn)行船舶粘性流場流體動力性能計算來分析船體與噴水推進(jìn)器之間的相互作用,目前正成為噴水推進(jìn)CFD研究的一個方向,其計算已達(dá)到較高的精度[8]。

    本文采用數(shù)值計算方法對實尺度雙泵推進(jìn)噴水推進(jìn)船的拖泵阻力進(jìn)行計算,分析噴泵拖曳阻力特性,并以此拖曳阻力作為參考,對某四泵推進(jìn)的噴水推進(jìn)船所采用的噴泵進(jìn)行最優(yōu)化選型,進(jìn)而完成噴水推進(jìn)器的設(shè)計。為驗證該四泵推進(jìn)噴水推進(jìn)船在18節(jié)工況下加速泵的拖曳阻力特性,對實尺度“船體+四臺噴水推進(jìn)器”系統(tǒng)帶自由液面的非定常流場進(jìn)行計算和分析,并探索了大尺度條件下船泵系統(tǒng)考慮自由液面和重力影響的瞬態(tài)計算方法。

    1 噴水推進(jìn)器數(shù)值計算模型校驗

    CFD方法已經(jīng)廣泛應(yīng)用于船舶推進(jìn)領(lǐng)域的設(shè)計和優(yōu)化問題,計算精度已遠(yuǎn)達(dá)到工程應(yīng)用許用誤差。本研究主要基于CFD方法完成實尺度噴泵拖泵阻力的計算和分析,在噴泵和船體性能計算分析之前,首先對本研究所采用的計算模型進(jìn)行校驗。以KaMeWa 71SII混流式噴水推進(jìn)泵為研究對象[9],利用CAD軟件(UG)構(gòu)造該泵的幾何模型(噴泵進(jìn)口直徑710 mm),圖1(a)為該混流式噴水推進(jìn)泵的幾何模型,在實際工作中該泵通過泵進(jìn)口前安裝的進(jìn)水流道從船底吸水,圖1(b)顯示了模擬該泵實際工作的幾何模型,其主要由混流泵、進(jìn)水流道、船底(用于模擬船底邊界層對噴泵進(jìn)流的影響)和水下計算域等組成?;炝魇絿娝七M(jìn)泵、進(jìn)水流道及水下控制體各區(qū)域均采用六面體結(jié)構(gòu)化網(wǎng)格進(jìn)行離散,如圖2所示。

    圖1 KaMeWa 71SII混流式噴水推進(jìn)泵數(shù)值計算幾何模型Fig.1 Geometry model of KaMeWa 71SII

    圖2 KaMeWa 71SII混流式噴水推進(jìn)泵數(shù)值計算網(wǎng)格示意圖Fig.2 Surface mesh of KaMeWa 71SII model

    采用SST湍流模型進(jìn)行計算求解,其融合了k-ε和k-ω兩種湍流模型,在自由流動區(qū)域使用k-ε模式,而在近壁面區(qū)域(y+<2.5)使用k-ω模型中的低雷諾數(shù)公式,兩者之間通過混合函數(shù)來過渡,這樣可以不需要使用壁面函數(shù),能夠較好地模擬粘性底層的流動[10]。

    湍動能方程:

    湍流耗散率方程:

    SST模型中各系數(shù)是k-ε模型和k-ω模型中相應(yīng)系數(shù)的線性組合,即模型中封閉系數(shù)選取如下。各系數(shù)取值為:

    圖3顯示了該混流式噴水推進(jìn)泵在各個轉(zhuǎn)速下功率的計算結(jié)果,其中,基于SST湍流模型定常模擬計算了該混流泵7個轉(zhuǎn)速下的流動特性,計算得到的葉輪功率與試驗數(shù)據(jù)最大誤差為2.9%,該混流泵在不均勻進(jìn)流條件下三個工作轉(zhuǎn)速的葉輪功率計算結(jié)果與廠商提供數(shù)據(jù)的最大誤差為2.1%,驗證了實尺度條件下噴水推進(jìn)器在均勻和非均勻條件下數(shù)值計算模型的準(zhǔn)確性和有效性。

    圖3 KaMeWa 71SII混流式噴水推進(jìn)泵數(shù)值計算結(jié)果與試驗值對比Fig.3 Comparison of power performance calculation of KaMeWa 71SII with test data

    2 某雙泵噴水推進(jìn)船拖泵特性數(shù)值模擬

    為研究在既定航速下噴泵鎖軸時在來流沖壓作用下,流體經(jīng)進(jìn)水流道、泵過流通道后從噴口流出所產(chǎn)生的拖泵阻力特性,本研究對某雙泵推進(jìn)的噴水推進(jìn)船在18節(jié)航速下雙泵處于鎖軸狀態(tài)下的流場進(jìn)行了實尺度數(shù)值計算,該船長132 m,寬6 m,正常排水量船體吃水8 m,噴泵進(jìn)口直徑1.95 m,噴口直徑1.2 m。為消除尺度效應(yīng)的影響,準(zhǔn)確模擬噴水推進(jìn)器的實際不均勻進(jìn)流,全面反映實船邊界層對噴水推進(jìn)泵進(jìn)流特性和推進(jìn)性能的真實影響,采用足尺1:1模型和SST湍流模型VOF模擬方法來求取“船體+噴水推進(jìn)泵+進(jìn)水流道”系統(tǒng)的帶自由液面的粘性流場。應(yīng)用不可壓縮的三維N-S方程模擬船泵系統(tǒng)的流體性能,采用有限體積法離散控制方程,對流項采用一階迎風(fēng)格式,擴(kuò)散項采用二階中心差分格式,基于SIMPLEC算法實現(xiàn)速度和壓力之間的耦合求解。

    2.1 控制體的選取和邊界條件設(shè)置

    由于該計算模型沿船體縱中剖面呈鏡像分布,只對縱中剖面一側(cè)的流場進(jìn)行建模與計算,以節(jié)省數(shù)值計算的時間。圖1所示為“船體+噴水推進(jìn)泵+進(jìn)水流道”系統(tǒng)的數(shù)值計算模型和邊界條件設(shè)置。計算域總長取5倍船長(船長L為132 m),寬度為1.5倍船長,船底至計算域最下端長度取0.5倍船長。采用了速度進(jìn)口和壓力出口邊界條件,進(jìn)口速度為18節(jié),出口壓力在水面以上為大氣壓力,水面以下為按照水深變化的水壓,船體縱中剖面設(shè)置為對稱面。

    圖4 “船體+噴泵+進(jìn)水流道”數(shù)值計算模型和邊界條件Fig.4 Numerical model and boundary conditions of the‘hull+waterjet+inlet duct’

    圖5 船艉板附近及噴泵表面網(wǎng)格Fig.5 Surface mesh of the ship stern and pumps

    圖6 拖泵工況下船泵系統(tǒng)自由液面分布圖Fig.6 Free surface of ship and pump at towing condition

    2.2 網(wǎng)格劃分

    “船體+噴水推進(jìn)泵+進(jìn)水流道”系統(tǒng)計算域采用全結(jié)構(gòu)化六面體網(wǎng)格進(jìn)行離散,總網(wǎng)格節(jié)點數(shù)約2100萬,圖5為船尾進(jìn)水口附近及噴水推進(jìn)泵網(wǎng)格分布圖。

    2.3 計算結(jié)果分析

    由于此時噴泵轉(zhuǎn)速為零,數(shù)值計算采用了基于SST湍流模型的穩(wěn)態(tài)計算方法來求取拖泵工況的“船體+噴水推進(jìn)泵+進(jìn)水流道”系統(tǒng)在18節(jié)下的粘性流場。圖6(a)顯示了數(shù)值計算得到的船體周圍的興波示意圖,圖6(b)顯示了拖泵工況下噴水推進(jìn)泵的射流。我們可以看出,雖然噴泵處于鎖軸狀態(tài),但是由于來流的沖壓作用,有流體經(jīng)進(jìn)水流道、泵過流通道后從噴口流出。數(shù)值計算穩(wěn)定后,采用壁面積分法計算了此時噴泵的拖曳阻力[9],計算結(jié)果表明在18節(jié)下該噴泵的拖曳阻力占船體阻力的1.80%。

    3 某四泵噴水推進(jìn)船拖泵特性的數(shù)值模擬

    以第一章中兩泵推進(jìn)噴水推進(jìn)船的拖泵阻力特性作為參考,對尺寸與之相近的某四泵推進(jìn)噴水推進(jìn)船的拖泵阻力進(jìn)行了預(yù)估(船長121 m,寬7.2 m,吃水6 m,大泵進(jìn)口直徑2.1 m,小泵進(jìn)口直徑1.6 m),然后結(jié)合該船阻力對其四臺噴水推進(jìn)器進(jìn)行了選型(其舷側(cè)的兩臺相同型號噴水推進(jìn)器用于18節(jié)航速下工況,中間兩臺相同型號噴水推進(jìn)器在高航速時使用)。進(jìn)而,采用泵的三元設(shè)計方法與數(shù)值試驗相互結(jié)合的方法完成噴泵的設(shè)計[11],運(yùn)用參數(shù)設(shè)計方法和數(shù)值模擬分析相互結(jié)合的方法完成進(jìn)水流道的設(shè)計。

    3.1 網(wǎng)格劃分

    “船體+噴水推進(jìn)泵+進(jìn)水流道”系統(tǒng)計算域采用全結(jié)構(gòu)化六面體網(wǎng)格進(jìn)行離散,船體外流場總網(wǎng)格節(jié)點數(shù)約2010萬,噴水推進(jìn)器流場總網(wǎng)格節(jié)點數(shù)約720萬,總網(wǎng)格節(jié)點數(shù)約2730萬。圖7(a)為進(jìn)水流道和船艉板附近網(wǎng)格分布圖,圖7(b)為兩類噴水推進(jìn)泵葉輪和導(dǎo)葉表面網(wǎng)格示意圖。

    圖7 計算域各個部件表面網(wǎng)格Fig.7 Surface mesh of parts of the model

    3.2 計算結(jié)果分析

    采用足尺1:1模型和VOF模擬方法來求取“船體+四臺噴水推進(jìn)器”系統(tǒng)的帶自由液面的粘性流場,運(yùn)用與第一章中雙泵推進(jìn)噴水推進(jìn)船相同的數(shù)值計算域大小和邊界條件。該四泵推進(jìn)噴水推進(jìn)船拖泵工況與兩泵噴水推進(jìn)泵拖泵工況有所區(qū)別,四泵推進(jìn)噴水推進(jìn)船在18節(jié)的拖泵工況下其兩臺加速泵處于鎖軸,另外兩臺弦側(cè)泵處于正常工作,因此四泵推進(jìn)噴水推進(jìn)船的拖泵工況的流場更為復(fù)雜,其船泵系統(tǒng)的數(shù)值模擬需要分兩步進(jìn)行。第一步采用穩(wěn)態(tài)計算來求取拖泵工況的“船體+四臺噴水推進(jìn)泵+四個進(jìn)水流道”的全船流場。此時,噴泵轉(zhuǎn)速設(shè)定為零,計算的時間步長適度加大,以加速模擬船體周圍的興波運(yùn)動,直至整個計算域流場穩(wěn)定。然后,以上一步拖泵工況的穩(wěn)態(tài)流場為初始值,采用瞬態(tài)計算方法求取噴水推進(jìn)泵開始轉(zhuǎn)動后的流場,以精確求取進(jìn)水流道進(jìn)水口附近、噴水推進(jìn)泵內(nèi)部以及噴口射流場的流動特性。在瞬態(tài)計算時噴水推進(jìn)器轉(zhuǎn)速需要由小到大漸進(jìn)的增加,直至舷側(cè)泵轉(zhuǎn)速增加到設(shè)計值。值得說明的是噴水推進(jìn)器每一次轉(zhuǎn)速增加后的瞬態(tài)計算需要以上一步計算的穩(wěn)定流場為初始值,轉(zhuǎn)速增加過快或流場未穩(wěn)定急加轉(zhuǎn)速都可能使得計算因為自由液面的劇烈變化而發(fā)散。圖8(a)顯示了巡航工況下加速泵處于鎖軸、弦側(cè)泵正常工作時的船體興波流場,圖8(b)為四泵均處于鎖軸的數(shù)值計算中間過渡狀態(tài)的噴泵射流情況,圖8(c)顯示了巡航工況下加速泵處于鎖軸、弦側(cè)泵正常工作時的四臺噴泵射流情況,圖8(d)顯示了拖泵工況下弦側(cè)泵內(nèi)部流線分布圖。從圖中我們可以看出在拖帶工況下弦側(cè)泵(額定轉(zhuǎn)速)和中間加速泵(鎖軸)時噴射水流的區(qū)別,數(shù)值計算結(jié)果表明在18節(jié)下該加速泵的拖曳阻力占船體阻力的1.63%,進(jìn)一步驗證了選型時所采用的拖泵阻力預(yù)估值的合理性。

    4 結(jié) 論

    本研究采用實尺度模型對某雙泵推進(jìn)噴水推進(jìn)船的船泵系統(tǒng)帶自由液面流場和拖曳阻力特性進(jìn)行了數(shù)值計算,并以此拖泵阻力為參考對尺寸與之相當(dāng)?shù)哪乘谋猛七M(jìn)的噴水推進(jìn)船的加速泵拖泵阻力進(jìn)行了預(yù)估。進(jìn)而,完成了該四泵推進(jìn)噴水推進(jìn)船所采用噴泵的選型和設(shè)計,并對該四泵推進(jìn)的噴水推進(jìn)船的“船體+四臺噴水推進(jìn)器”系統(tǒng)帶自由液面的非定常流場以及拖泵阻力進(jìn)行了計算,驗證了拖泵阻力預(yù)估的合理性,也為噴水推進(jìn)船拖泵阻力特性的研究提供有益參考。另外,本文探索了大尺度條件下船泵系統(tǒng)考慮自由液面和重力影響的瞬態(tài)計算方法,研究認(rèn)為對于該尺度下拖泵數(shù)值模型復(fù)雜的非定常流場計算需要采用分步計算的方法進(jìn)行,首先對該船在18節(jié)航速下全部泵處于鎖軸工況時的穩(wěn)態(tài)流場進(jìn)行數(shù)值計算,計算的時間步長適度加大,以加速模擬船體附近的興波流場;然后以穩(wěn)態(tài)計算結(jié)果為初始值進(jìn)行該工況下瞬態(tài)流場計算,以提高非定常計算的收斂速度和穩(wěn)定性;最后,以上一步瞬態(tài)流場為初始值,逐步增加泵轉(zhuǎn)速直到達(dá)到既定轉(zhuǎn)速。

    [1]Thor Peter Andersen.Design of rim driven waterjet pump for small rescue vessel[D].Gothenburg,Sweden:Chalmers University of Technology,2014.

    [2]Ronald O’Rourke.Navy Littoral Combat Ship(LCS)program:Background,issues,and options for congress[R].America Congressional Research Service Report for congress,2011.

    [3]靳栓寶.高性能噴水推進(jìn)器與泵噴推進(jìn)器的設(shè)計與性能分析[D].武漢:海軍工程大學(xué),2013. Jin Shuanbao.Design and analysis of high performance waterjet and pumpjet[D].Wuhan:College of Naval Architecture and Power,2013

    [4]Kruppa C.Report of the specialist committee on waterjets[R].Paris:21st ITTC,1996.

    [5]Hoyt III J G.Report of the specialist committee on waterjets[R].Netherlands:22nd ITTC,1999.

    [6]Van T.The specialist committee on validation of waterjet test procedure[R].Italy:23rd ITTC,2002.

    [7]Van T.The specialist committee on validation of waterjet test procedure[R].UK:24th ITTC,2005.

    [8]易文彬,王永生,靳栓寶等.雷諾數(shù)相似噴水推進(jìn)船模預(yù)報實船推進(jìn)性能[J].武漢理工大學(xué)學(xué)報,2014,38(5):1074-1078. Yi Wenbin,Wang Yongsheng,Jin Shuanbao,et al.Prediction of ship propulsion performance using Reynolds similarity waterjet ship model[J].Journal of Wuhan University of Techonlogy,2014,38(5):1074-1078.

    [9]董新國.混流式噴水推進(jìn)泵動葉輪、導(dǎo)葉體結(jié)構(gòu)強(qiáng)度分析[D].上海:上海交通大學(xué),2011. Dong Xinguo.Structure strength analysis of impeller and guide vane chamber of mixed flow waterjet pump[D].Shanghai:Shanghai Jiaotong University,2011.

    [10]孫存樓,王永生.噴水推進(jìn)器推力預(yù)報的兩種不同方法比較[J].船舶力學(xué),2008,14(11):1208-1212. Sun Cunlou,Wang Yongsheng.Comparison of two different methods for thrust prediction of marine waterjet[J].Journal of Ship Mechanics,2010,14(11):1208-1212.

    [11]靳栓寶,王永生,丁江明,蘇永生.混流式噴水推進(jìn)泵三元設(shè)計及數(shù)值試驗[J].哈爾濱工程大學(xué)學(xué)報,2012,33(10): 1223-1227. Jin Shuanbao,Wang Yongsheng,Ding Jiangming,et al.Three-dimensional design and numerical experiment of mixedflow waterjet[J].Journal of Harbin Engineering University,2012,33(10):1223-1227.

    Research on the pump towing drag of full scale waterjet ship with CFD

    JIN Shuan-bao,SHEN Yang,WANG Dong,LIU Tao,WEI Ying-san,HU Peng-fei,ZHU Hao
    (National Key Laboratory of Science and Technology on Vessel Integrated Power System,Naval University of Engineering,Wuhan 430033,China)

    Accelerating pump towing drag of the ship propelled by four waterjet pumps at cruising condition has great influence on the process of waterjet selection and optimal match between hull,waterjets and the main engine.However,the pump towing drag is difficult to be obtained from the ship model test,and so far there is little record about the pump towing drag at home and abroad.For this reason,CFD method of waterjet in uniform and non-uniform flow were validated,and CFD method was used for the simulation of the viscous flow of the full scale ship propelled by two waterjets,and the pump towing drag and its percentage on the ship drag could be calculated.Taking this pump towing drag as a reference,the waterjet pumps of the ship propelled by four waterjet were selctected and designed.And then the towing drag of accelerating pumps of the ship could be calculated with CFD simulation,the results indicate that the percentage of the pump towing drag on the ship drag were roughly the same as hypothetical.In order to eliminate the scale effect,the method for calculating the unsteady viscous flow of full scale‘hull+four waterjets’model with gravity and free surface was researched and introduced briefly.

    ship;waterjet;full scale;pump towing drag;numerical simulation;unsteady

    U664.34

    A

    10.3969/j.issn.1007-7294.2016.07.004

    1007-7294(2016)11-1381-07

    2016-04-03

    國家自然科學(xué)基金資助項目(51309229,51307177,51409256)

    靳栓寶(1983-),男,博士后,E-mail:hjgcjin@163.com;沈 洋(1983-),男,講師。

    猜你喜歡
    推進(jìn)器船體流場
    船體行駛過程中的壓力監(jiān)測方法
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場計算
    基于CFD扇翼推進(jìn)器敞水性能預(yù)報分析
    發(fā)揮考核“指揮棒”“推進(jìn)器”作用
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計分析
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    焊接殘余應(yīng)力對船體結(jié)構(gòu)疲勞強(qiáng)度的影響分析
    焊接(2015年9期)2015-07-18 11:03:51
    讓黨建成為信仰播種機(jī)和工作推進(jìn)器
    赴美軍“仁慈”號醫(yī)院船駐船體會
    日本与韩国留学比较| xxx大片免费视频| 亚洲av综合色区一区| av不卡在线播放| 少妇 在线观看| 在线观看三级黄色| 亚洲一区二区三区欧美精品| 巨乳人妻的诱惑在线观看| 毛片一级片免费看久久久久| 国产麻豆69| 91午夜精品亚洲一区二区三区| 成人亚洲精品一区在线观看| 国产又爽黄色视频| 婷婷色av中文字幕| 国语对白做爰xxxⅹ性视频网站| 免费观看性生交大片5| 精品久久国产蜜桃| 有码 亚洲区| 看免费成人av毛片| 国产日韩欧美在线精品| 精品一区二区三卡| 久久久久久人妻| 午夜福利视频在线观看免费| 精品国产一区二区三区久久久樱花| 99香蕉大伊视频| 美女xxoo啪啪120秒动态图| 亚洲欧美中文字幕日韩二区| 午夜福利乱码中文字幕| 高清黄色对白视频在线免费看| 精品一区二区三区四区五区乱码 | 免费少妇av软件| 欧美日本中文国产一区发布| freevideosex欧美| 久久久精品94久久精品| 97精品久久久久久久久久精品| 亚洲国产毛片av蜜桃av| 寂寞人妻少妇视频99o| 精品一区二区免费观看| 国语对白做爰xxxⅹ性视频网站| 亚洲图色成人| 最近中文字幕2019免费版| 久久精品国产自在天天线| 两个人看的免费小视频| 欧美xxxx性猛交bbbb| 美女中出高潮动态图| 久久久久久久久久久免费av| 最新中文字幕久久久久| 丝袜美足系列| 99久久人妻综合| 亚洲性久久影院| 亚洲精品乱码久久久久久按摩| 国产国拍精品亚洲av在线观看| 一本大道久久a久久精品| av线在线观看网站| 香蕉丝袜av| 美女国产视频在线观看| 99久久精品国产国产毛片| 内地一区二区视频在线| 国产熟女欧美一区二区| 中文精品一卡2卡3卡4更新| 性色avwww在线观看| 久久精品久久久久久久性| 美女中出高潮动态图| 在线看a的网站| 99热网站在线观看| 欧美人与性动交α欧美软件 | av又黄又爽大尺度在线免费看| 久久午夜福利片| 久久久久久久大尺度免费视频| 91精品三级在线观看| 日本av免费视频播放| 精品第一国产精品| 另类精品久久| 伦理电影大哥的女人| 国产一区二区在线观看日韩| 欧美日韩国产mv在线观看视频| 欧美成人午夜免费资源| 最近最新中文字幕大全免费视频 | 久久久欧美国产精品| 欧美精品一区二区大全| 中文字幕最新亚洲高清| 妹子高潮喷水视频| 久久鲁丝午夜福利片| 男人操女人黄网站| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩一区二区三区在线| 欧美日韩综合久久久久久| 少妇人妻久久综合中文| 国产精品熟女久久久久浪| 亚洲欧美色中文字幕在线| 日本免费在线观看一区| 免费不卡的大黄色大毛片视频在线观看| av在线播放精品| 欧美3d第一页| 成年人免费黄色播放视频| 免费观看a级毛片全部| 日本av手机在线免费观看| 亚洲成人一二三区av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久女婷五月综合色啪小说| 国产一区亚洲一区在线观看| 少妇 在线观看| 蜜桃在线观看..| 日产精品乱码卡一卡2卡三| 久久久久人妻精品一区果冻| 女的被弄到高潮叫床怎么办| 国产av码专区亚洲av| 最近最新中文字幕免费大全7| 我要看黄色一级片免费的| 午夜影院在线不卡| 超色免费av| 黑人欧美特级aaaaaa片| 久久久久久人妻| 成人亚洲欧美一区二区av| 久久这里只有精品19| 好男人视频免费观看在线| 国产亚洲av片在线观看秒播厂| 全区人妻精品视频| 黑人猛操日本美女一级片| 秋霞在线观看毛片| videosex国产| 亚洲精品色激情综合| 亚洲av成人精品一二三区| 91国产中文字幕| 亚洲精品,欧美精品| 久久人妻熟女aⅴ| 永久免费av网站大全| 欧美人与善性xxx| 日韩电影二区| 人妻 亚洲 视频| kizo精华| www.av在线官网国产| 97超碰精品成人国产| 热99久久久久精品小说推荐| 亚洲在久久综合| 亚洲精品美女久久久久99蜜臀 | 两性夫妻黄色片 | 午夜免费男女啪啪视频观看| 免费人妻精品一区二区三区视频| 校园人妻丝袜中文字幕| 丝袜喷水一区| 欧美日韩成人在线一区二区| 99热6这里只有精品| 在线精品无人区一区二区三| 18禁裸乳无遮挡动漫免费视频| 亚洲av免费高清在线观看| 午夜免费观看性视频| 亚洲精品,欧美精品| 国产1区2区3区精品| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 国产精品一区二区在线不卡| 韩国精品一区二区三区 | 秋霞在线观看毛片| 国产高清不卡午夜福利| 有码 亚洲区| 在线观看www视频免费| 伦理电影大哥的女人| 成人国产av品久久久| 夫妻午夜视频| 久久久久久人人人人人| 亚洲av日韩在线播放| 丰满乱子伦码专区| 久久精品国产a三级三级三级| 国产熟女欧美一区二区| 在线观看美女被高潮喷水网站| 韩国精品一区二区三区 | 欧美3d第一页| 制服丝袜香蕉在线| 久久99热这里只频精品6学生| 欧美成人午夜免费资源| 母亲3免费完整高清在线观看 | 好男人视频免费观看在线| 久久久国产欧美日韩av| 最新的欧美精品一区二区| 看免费av毛片| 日韩视频在线欧美| 免费在线观看黄色视频的| 99九九在线精品视频| 欧美97在线视频| 日产精品乱码卡一卡2卡三| 久久 成人 亚洲| 97在线人人人人妻| 欧美亚洲 丝袜 人妻 在线| 观看美女的网站| 九九爱精品视频在线观看| 欧美激情国产日韩精品一区| 日韩av在线免费看完整版不卡| 亚洲色图综合在线观看| 亚洲精品456在线播放app| kizo精华| 国产一区亚洲一区在线观看| 免费看av在线观看网站| 亚洲欧美一区二区三区黑人 | 久久人人爽人人片av| 五月开心婷婷网| 啦啦啦在线观看免费高清www| 韩国av在线不卡| 亚洲中文av在线| 大片电影免费在线观看免费| 亚洲国产成人一精品久久久| 日韩av在线免费看完整版不卡| 国产免费视频播放在线视频| 免费在线观看完整版高清| 亚洲av男天堂| 美女视频免费永久观看网站| 成人毛片a级毛片在线播放| 看免费av毛片| 两个人免费观看高清视频| 观看av在线不卡| 2021少妇久久久久久久久久久| 亚洲精华国产精华液的使用体验| 黑人高潮一二区| 免费看av在线观看网站| 日韩人妻精品一区2区三区| 啦啦啦中文免费视频观看日本| 又黄又爽又刺激的免费视频.| 国产精品无大码| 最近手机中文字幕大全| 亚洲中文av在线| 亚洲av福利一区| 亚洲成人一二三区av| 免费女性裸体啪啪无遮挡网站| 9色porny在线观看| 一级爰片在线观看| 人人妻人人澡人人看| 午夜激情av网站| 精品少妇黑人巨大在线播放| av又黄又爽大尺度在线免费看| 亚洲国产欧美日韩在线播放| 亚洲精品自拍成人| 一级毛片我不卡| 亚洲精华国产精华液的使用体验| 亚洲av福利一区| 晚上一个人看的免费电影| 国产福利在线免费观看视频| videossex国产| 亚洲精品一二三| 桃花免费在线播放| 国产爽快片一区二区三区| 精品福利永久在线观看| 最后的刺客免费高清国语| 亚洲精华国产精华液的使用体验| 九色成人免费人妻av| 男人舔女人的私密视频| av网站免费在线观看视频| 久久 成人 亚洲| 欧美变态另类bdsm刘玥| 日韩中文字幕视频在线看片| 久久久久久人妻| 午夜日本视频在线| 五月天丁香电影| 少妇高潮的动态图| 一级毛片电影观看| 国产成人精品一,二区| 免费观看a级毛片全部| 热99国产精品久久久久久7| 黑人巨大精品欧美一区二区蜜桃 | 免费人妻精品一区二区三区视频| 大码成人一级视频| 久久人人97超碰香蕉20202| 视频在线观看一区二区三区| 亚洲国产色片| 日本免费在线观看一区| 日韩欧美精品免费久久| 黄色配什么色好看| 亚洲国产精品一区三区| 日韩欧美一区视频在线观看| 亚洲成国产人片在线观看| 中文字幕制服av| 亚洲av中文av极速乱| 亚洲综合精品二区| 日韩熟女老妇一区二区性免费视频| 精品一区在线观看国产| 精品人妻在线不人妻| 美女国产视频在线观看| 九九爱精品视频在线观看| av在线观看视频网站免费| 又黄又粗又硬又大视频| 久久精品人人爽人人爽视色| 日本av手机在线免费观看| 日韩一区二区视频免费看| 草草在线视频免费看| 亚洲精品中文字幕在线视频| 超色免费av| 丝袜脚勾引网站| 99精国产麻豆久久婷婷| 嫩草影院入口| 成人亚洲精品一区在线观看| 91精品国产国语对白视频| 色5月婷婷丁香| 一边亲一边摸免费视频| 91成人精品电影| 中文字幕人妻熟女乱码| 亚洲精品中文字幕在线视频| 国产黄频视频在线观看| 建设人人有责人人尽责人人享有的| 男女免费视频国产| 中文字幕精品免费在线观看视频 | av不卡在线播放| 看十八女毛片水多多多| 日日爽夜夜爽网站| 中国美白少妇内射xxxbb| 免费久久久久久久精品成人欧美视频 | 久久久久久人人人人人| 日本午夜av视频| 人妻系列 视频| 午夜福利在线观看免费完整高清在| 免费观看无遮挡的男女| 国产成人精品福利久久| 黄色毛片三级朝国网站| 在线天堂中文资源库| 黄片播放在线免费| 欧美人与善性xxx| 精品一区在线观看国产| 亚洲丝袜综合中文字幕| 国产国拍精品亚洲av在线观看| 日本91视频免费播放| 久久精品久久久久久噜噜老黄| 建设人人有责人人尽责人人享有的| 伊人亚洲综合成人网| 亚洲四区av| 国产男女超爽视频在线观看| 亚洲婷婷狠狠爱综合网| 看十八女毛片水多多多| 美女福利国产在线| 伦理电影大哥的女人| 日本av手机在线免费观看| 日韩成人伦理影院| 国产视频首页在线观看| 久久亚洲国产成人精品v| av一本久久久久| 男人添女人高潮全过程视频| 国产激情久久老熟女| 成人综合一区亚洲| 久久久国产一区二区| 亚洲国产毛片av蜜桃av| 自拍欧美九色日韩亚洲蝌蚪91| 国产深夜福利视频在线观看| www.色视频.com| 国产精品蜜桃在线观看| 日本色播在线视频| 欧美亚洲日本最大视频资源| 美女主播在线视频| 午夜精品国产一区二区电影| 亚洲美女视频黄频| av黄色大香蕉| 国产欧美另类精品又又久久亚洲欧美| 日本91视频免费播放| 久久久久久伊人网av| 国产爽快片一区二区三区| 欧美bdsm另类| 在线亚洲精品国产二区图片欧美| 国产欧美日韩综合在线一区二区| 男人舔女人的私密视频| 日本wwww免费看| 三上悠亚av全集在线观看| 男人操女人黄网站| 午夜精品国产一区二区电影| 亚洲国产av影院在线观看| 韩国精品一区二区三区 | 日韩电影二区| 97在线视频观看| 欧美97在线视频| xxx大片免费视频| 国产又爽黄色视频| 国国产精品蜜臀av免费| 蜜桃在线观看..| 中国国产av一级| 亚洲精品一二三| 丰满迷人的少妇在线观看| 蜜桃在线观看..| 制服人妻中文乱码| 久久久久精品人妻al黑| 91精品国产国语对白视频| 看十八女毛片水多多多| 毛片一级片免费看久久久久| 国产av码专区亚洲av| 性色avwww在线观看| 国产有黄有色有爽视频| 国产精品熟女久久久久浪| 满18在线观看网站| 中文字幕精品免费在线观看视频 | 三级国产精品片| 超色免费av| 啦啦啦啦在线视频资源| 五月伊人婷婷丁香| 丝瓜视频免费看黄片| 少妇人妻久久综合中文| 日本欧美视频一区| 多毛熟女@视频| 午夜免费男女啪啪视频观看| 色视频在线一区二区三区| 久久鲁丝午夜福利片| av国产精品久久久久影院| 亚洲中文av在线| 大片免费播放器 马上看| 岛国毛片在线播放| 丰满乱子伦码专区| 久久久久网色| 欧美xxxx性猛交bbbb| 国产欧美亚洲国产| 777米奇影视久久| 国产成人欧美| 国产免费福利视频在线观看| 一区二区三区四区激情视频| av视频免费观看在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美丝袜亚洲另类| 狠狠婷婷综合久久久久久88av| 国产成人av激情在线播放| 免费黄网站久久成人精品| 国产免费现黄频在线看| 久久久久人妻精品一区果冻| 丝袜喷水一区| 免费在线观看黄色视频的| 久久久久国产精品人妻一区二区| 国产精品一区二区在线观看99| 又粗又硬又长又爽又黄的视频| 黄色 视频免费看| 亚洲精品乱码久久久久久按摩| 一本久久精品| 看免费成人av毛片| 成人18禁高潮啪啪吃奶动态图| h视频一区二区三区| 色婷婷av一区二区三区视频| 欧美成人午夜精品| 亚洲精品aⅴ在线观看| 国产精品嫩草影院av在线观看| 男的添女的下面高潮视频| 免费黄色在线免费观看| 99香蕉大伊视频| 亚洲伊人色综图| 日韩中文字幕视频在线看片| 午夜免费男女啪啪视频观看| av在线播放精品| 大陆偷拍与自拍| www日本在线高清视频| 欧美精品人与动牲交sv欧美| 国产无遮挡羞羞视频在线观看| 国产极品粉嫩免费观看在线| 丝袜人妻中文字幕| 亚洲精华国产精华液的使用体验| 91国产中文字幕| 九色亚洲精品在线播放| 日本色播在线视频| 国产在线一区二区三区精| 精品午夜福利在线看| 国产精品国产av在线观看| 国产综合精华液| 亚洲经典国产精华液单| 美女国产高潮福利片在线看| 国产黄频视频在线观看| 久久久久久人妻| 亚洲精品色激情综合| 国产午夜精品一二区理论片| 天堂8中文在线网| videosex国产| 国产片特级美女逼逼视频| 国产精品国产三级专区第一集| 超碰97精品在线观看| 国产麻豆69| 亚洲国产看品久久| 人人妻人人添人人爽欧美一区卜| 国产高清不卡午夜福利| 亚洲欧美一区二区三区国产| 91aial.com中文字幕在线观看| 久久久久国产网址| 精品卡一卡二卡四卡免费| 制服丝袜香蕉在线| 啦啦啦啦在线视频资源| 99热网站在线观看| 亚洲人成网站在线观看播放| 免费观看性生交大片5| 亚洲av电影在线进入| av在线播放精品| 久久久久久久久久久免费av| 春色校园在线视频观看| 午夜福利在线观看免费完整高清在| 国产成人aa在线观看| 日本午夜av视频| av在线老鸭窝| 制服人妻中文乱码| 夫妻性生交免费视频一级片| 丝袜脚勾引网站| 日本色播在线视频| 国产亚洲av片在线观看秒播厂| 久久久精品94久久精品| 看十八女毛片水多多多| 人妻少妇偷人精品九色| 99热国产这里只有精品6| 国产一区有黄有色的免费视频| 丰满乱子伦码专区| 成人亚洲精品一区在线观看| 熟女人妻精品中文字幕| 99热国产这里只有精品6| 日韩视频在线欧美| xxxhd国产人妻xxx| 免费观看av网站的网址| 国产亚洲一区二区精品| 十分钟在线观看高清视频www| 午夜视频国产福利| 在线观看免费高清a一片| 全区人妻精品视频| 国产亚洲一区二区精品| 中国三级夫妇交换| 亚洲国产精品999| 午夜91福利影院| 天天躁夜夜躁狠狠久久av| 国产欧美日韩一区二区三区在线| 少妇的逼好多水| 欧美日本中文国产一区发布| 老司机影院成人| 亚洲美女黄色视频免费看| 精品亚洲成a人片在线观看| 国产亚洲最大av| 伦理电影大哥的女人| 亚洲国产色片| 国产欧美日韩综合在线一区二区| 美女中出高潮动态图| 91在线精品国自产拍蜜月| 亚洲欧美一区二区三区黑人 | av在线老鸭窝| 欧美日韩视频高清一区二区三区二| 全区人妻精品视频| 热re99久久精品国产66热6| 最近最新中文字幕免费大全7| 热99久久久久精品小说推荐| 午夜91福利影院| 18在线观看网站| 高清在线视频一区二区三区| 精品一区在线观看国产| 欧美日韩综合久久久久久| 国产乱来视频区| 国产成人免费观看mmmm| 欧美激情 高清一区二区三区| 精品卡一卡二卡四卡免费| 男女边摸边吃奶| 精品亚洲乱码少妇综合久久| 精品一区在线观看国产| 99精国产麻豆久久婷婷| 久久久精品区二区三区| 亚洲精品色激情综合| 亚洲精品美女久久久久99蜜臀 | 男人操女人黄网站| 日韩视频在线欧美| 啦啦啦视频在线资源免费观看| 伊人亚洲综合成人网| 亚洲美女搞黄在线观看| 成人影院久久| 色网站视频免费| 亚洲欧美日韩另类电影网站| 亚洲国产精品专区欧美| 日韩在线高清观看一区二区三区| 亚洲综合精品二区| 制服诱惑二区| 成年女人在线观看亚洲视频| 久久 成人 亚洲| 亚洲精品乱码久久久久久按摩| 亚洲国产精品国产精品| 两个人看的免费小视频| 寂寞人妻少妇视频99o| 国产欧美日韩一区二区三区在线| 免费av中文字幕在线| 久久久国产精品麻豆| 亚洲人与动物交配视频| 菩萨蛮人人尽说江南好唐韦庄| 18+在线观看网站| 日韩欧美精品免费久久| 国产一区二区三区av在线| 青春草视频在线免费观看| 卡戴珊不雅视频在线播放| 天天躁夜夜躁狠狠久久av| 国产乱人偷精品视频| 久久人人爽av亚洲精品天堂| 成年美女黄网站色视频大全免费| 欧美成人精品欧美一级黄| 97在线视频观看| 熟女电影av网| 亚洲成人一二三区av| 免费久久久久久久精品成人欧美视频 | 波多野结衣一区麻豆| 成人综合一区亚洲| 国产精品久久久久久av不卡| 婷婷色麻豆天堂久久| videos熟女内射| 最近中文字幕高清免费大全6| 精品人妻偷拍中文字幕| 精品亚洲乱码少妇综合久久| 视频中文字幕在线观看| 久久久久久久久久久久大奶| 亚洲伊人色综图| 视频中文字幕在线观看| 美女大奶头黄色视频| 精品亚洲成a人片在线观看| xxxhd国产人妻xxx| 91午夜精品亚洲一区二区三区| 午夜福利视频在线观看免费| 一区二区三区四区激情视频| 午夜日本视频在线| 在线观看一区二区三区激情| 卡戴珊不雅视频在线播放| 久久99精品国语久久久| 国产成人精品一,二区| 久久韩国三级中文字幕| 日本与韩国留学比较| 精品一品国产午夜福利视频| 18禁在线无遮挡免费观看视频| 制服诱惑二区| 赤兔流量卡办理| av不卡在线播放| 亚洲美女视频黄频| 亚洲成av片中文字幕在线观看 | 久久久久久人妻| 久久人人爽av亚洲精品天堂|