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

    螺旋槳敞水性能計(jì)算及堵塞效應(yīng)研究

    2016-12-12 03:07:24瞿沐淋孫江龍黃本燊馬超龍
    艦船科學(xué)技術(shù) 2016年11期
    關(guān)鍵詞:螺旋槳扭矩坐標(biāo)系

    瞿沐淋,孫江龍,2,3,黃本燊,鐘 誠,馬超龍

    (1. 華中科技大學(xué) 船舶與海洋工程學(xué)院,湖北 武漢 430074;2. 船舶與海洋工程水動力湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430074;3. 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海 200240)

    螺旋槳敞水性能計(jì)算及堵塞效應(yīng)研究

    瞿沐淋1,孫江龍1,2,3,黃本燊1,鐘 誠1,馬超龍1

    (1. 華中科技大學(xué) 船舶與海洋工程學(xué)院,湖北 武漢 430074;2. 船舶與海洋工程水動力湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430074;3. 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海 200240)

    根據(jù)螺旋槳的投影原理以及其幾何參數(shù),用三維建模軟件 CATIA 建立三維螺旋槳數(shù)值模型。根據(jù)計(jì)算流體動力學(xué)(CFD)原理,使用流體動力學(xué)軟件 Fluent 對螺旋槳數(shù)值模型進(jìn)行分析計(jì)算。采用 RANS 方法結(jié)合RSM 湍流模型求解螺旋槳三維粘性流場,計(jì)算域的離散采用非結(jié)構(gòu)網(wǎng)格方法,運(yùn)用相對旋轉(zhuǎn)坐標(biāo)方法(MRF)來模擬螺旋槳的運(yùn)動,以此求出該螺旋槳在常態(tài)以及堵塞效應(yīng)下的流場特性,并將螺旋槳的數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行對比以確定該方法的適用性。最后研究堵塞效應(yīng)的相關(guān)性質(zhì)并將螺旋槳普通敞水性能與螺旋槳在堵塞效應(yīng)下的敞水性能進(jìn)行對比,得出堵塞效應(yīng)對螺旋槳敞水性能的影響。

    螺旋槳;CFD;MRF;敞水性能;堵塞效應(yīng)

    0 引 言

    由于數(shù)值模擬方法相對于實(shí)驗(yàn)方法擁有成本低、周期短、操作方便等一系列優(yōu)點(diǎn),故而計(jì)算流體動力學(xué)(CFD)應(yīng)用十分廣泛。

    本文利用 CATIA 軟件,以螺旋槳的投影原理和螺旋槳的相關(guān)型值參數(shù)為基礎(chǔ),建立三維螺旋槳模型。用 Ansys icem 建立相應(yīng)流域、計(jì)算域,再劃分非結(jié)構(gòu)性網(wǎng)格。應(yīng)用 CFD 軟件(Ansys,F(xiàn)luent),運(yùn)用 MRF技術(shù),首先對螺旋槳的敞水性能進(jìn)行數(shù)值模擬并將得到的結(jié)果與實(shí)驗(yàn)值對比,發(fā)現(xiàn)計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果測

    量值吻合良好。進(jìn)而分析螺旋槳在敞水中的壓力、速度分布。通過計(jì)算螺旋槳在不同情況的堵塞效應(yīng)下的敞水性能,分析、研究堵塞效應(yīng)的性質(zhì)。并比較堵塞效應(yīng)下與常態(tài)下螺旋槳的敞水性能,得到堵塞效應(yīng)對于螺旋槳的敞水性能的影響,這對于研究螺旋槳的敞水性能具有一定的參考價(jià)值。

    1 計(jì)算基本公式

    1.1 RANS 方程

    RANS 方程是粘性流體運(yùn)動學(xué)和動力學(xué)的普適性控制方程,本文用該方程作為求解螺旋槳水動力性能計(jì)算的基本方程。其具體方程式為:

    式中:ρ 為流體密度;P 為靜壓;fi為單位質(zhì)量的質(zhì)量力;ui,uj分別為速度分量。

    1.2 湍流脈動動能方程

    湍流脈動動能方程即 k 方程:

    1.3 雷諾應(yīng)力模型方程

    雷諾應(yīng)力模型是求解雷諾應(yīng)力張量的各個(gè)分量的輸運(yùn)方程。具體形式如下:

    式中:Cij為對流項(xiàng);Dijr為湍流擴(kuò)散項(xiàng);Dijl為分子擴(kuò)散項(xiàng);Pij為應(yīng)力產(chǎn)生項(xiàng);Gij為浮力產(chǎn)生項(xiàng);?ij為壓力應(yīng)變項(xiàng);εij為耗散項(xiàng) ,F(xiàn)ij為系統(tǒng)旋轉(zhuǎn)產(chǎn)生項(xiàng)。

    2 模型建立

    2.1 螺旋槳投影原理

    DTMB4119 槳是一種無側(cè)斜無后傾分布的 3 葉螺旋槳,被 ITTC 選為考證數(shù)值方法預(yù)報(bào)精度的標(biāo)準(zhǔn)。其直徑為 0.304 8 m,螺距比為 1.084,轂徑比為 0.2,剖面為 NACA-66mod 型。

    根據(jù)螺旋槳投影原理來建立模型(見圖 1),首先將坐標(biāo)系 O1X1Y1Z1中的型值點(diǎn)經(jīng)過一次旋轉(zhuǎn)得到坐標(biāo)系 O1UVW 中坐標(biāo),然后由坐標(biāo)系 O′X′Y′Z′ 與O1UVW 的關(guān)系和圖1 中平面坐標(biāo)與圖中柱坐標(biāo)的關(guān)系,可推導(dǎo)出螺旋槳葉切面處局部坐標(biāo)系 O1X1Y1Z1到全局坐標(biāo)系 OXYZ 的坐標(biāo)變換公式為:

    式中:θ 為縱傾角,φ 為螺距角;L 為長度;r 為切面所在圓柱面的半徑。全局坐標(biāo)系 OXYZ 的 OXY 平面與漿轂斷面平行,坐標(biāo)系與 OXYZ 平行,坐標(biāo)系的軸經(jīng)過葉切面的最厚處,坐標(biāo)系經(jīng)過一次旋轉(zhuǎn)得到。

    圖1 螺旋槳投影原理Fig.1 Propeller projection principle

    2.2 模型建造

    將以型值表形式表達(dá)的螺旋槳葉切面的局部坐標(biāo)轉(zhuǎn)換為全局坐標(biāo)后,在 CAD 軟件(本文使用 CATIA)中進(jìn)行三維建模,坐標(biāo)軸的定義為:Z 軸與螺旋槳旋轉(zhuǎn)軸一致,來流方向指向 z 軸負(fù)向;X 軸與槳葉參考線一致;Y 軸服從右手系(見圖 2)。計(jì)算時(shí)采用全尺模型,轉(zhuǎn)速為 20 r/s,通過改變進(jìn)流速度來實(shí)現(xiàn)不同的進(jìn)速系數(shù)下螺旋槳的計(jì)算。

    對于堵塞效應(yīng)而言,由于需要研究起堵塞效應(yīng)而導(dǎo)致螺旋槳敞水性能的變化,故而需要在螺旋槳前方、計(jì)算域內(nèi)建立一個(gè)或者數(shù)個(gè)圓柱體來模擬堵塞固體塊。本文建立以下 2 種對比模型:

    1)建一個(gè)半徑為 25 mm、高 10 mm 的圓柱體軸線平行 z 軸,其 z 軸正向圓柱表面圓心為(0,70,172.4),這是加 1 個(gè)堵塞物的模型(見圖 3)。

    圖3 添加 1 個(gè)堵塞物模型Fig.3 Add a blockage model

    2)與上一模型類似,建立 5 個(gè)圓柱,半徑為 20 mm,高 15 mm。其 z 軸正向圓柱表面圓心分別為(130,0,146.4),(80,100,250),(30,30,200.4),(80,80,200.4),(0,–80,200.4)。這是加 5 個(gè)堵塞物的模型(見圖 4)。

    圖4 添加 5 個(gè)堵塞物模型Fig.4 Add five blockage model

    3 網(wǎng)格劃分及計(jì)算方法

    對于 CFD 軟件數(shù)值模擬來說,網(wǎng)格的劃分是最為費(fèi)時(shí)、費(fèi)力的環(huán)節(jié)。因?yàn)樵摥h(huán)節(jié)會影響到后來的計(jì)算,包括模擬計(jì)算的精度以及效率,網(wǎng)格過密或是過稀都會導(dǎo)致計(jì)算環(huán)節(jié)出問題。如果網(wǎng)格過于稀疏,則模擬精度不夠,得到一個(gè)不準(zhǔn)切的值,甚至是錯(cuò)誤的結(jié)果。但若是網(wǎng)格過密,則會導(dǎo)致計(jì)算量增大,耗時(shí)過長,甚至無法收斂。

    本文采用 Ansys icem 對螺旋槳劃網(wǎng)格。局部加密,其余區(qū)域網(wǎng)格相對較稀,以便于計(jì)算。為模擬無限流暢,所建立的流域必須足夠大,為方便計(jì)算,流域尺寸選取為螺旋槳直徑的整數(shù)倍,故而建立一個(gè)直徑為 3 倍的螺旋槳直徑,長 5 倍的螺旋槳直徑的圓柱體作為流域,其中螺旋槳距離流域來流部分為 1 倍直徑,距尾部 4 倍直徑。為模擬螺旋槳轉(zhuǎn)動,貼近螺旋槳建立一個(gè)小的計(jì)算域,并運(yùn)用 MRF 技術(shù)即相對旋轉(zhuǎn)的原理,通過旋轉(zhuǎn)坐標(biāo)系來模擬螺旋槳的轉(zhuǎn)動。其中小的計(jì)算域?yàn)橹睆綖?605.6 mm,長 509.6 mm 的圓柱。

    為減小網(wǎng)格數(shù)目,故而設(shè)計(jì)流域網(wǎng)格密度較稀而計(jì)算域的網(wǎng)格相對較密,且在螺旋槳表面處加密。本文建立的都是非結(jié)構(gòu)性網(wǎng)格。建立流域與計(jì)算域之間的區(qū)域?yàn)?YU1,建立計(jì)算域與螺旋槳之間的區(qū)域?yàn)閅U2。通過 y-plus 計(jì)算可得,螺旋槳第 1 層網(wǎng)格大小在 1.6 ~ 6.3 mm 之間。故而其他網(wǎng)格參數(shù)設(shè)定為 YU1為流域,其最大網(wǎng)格尺寸為 30 mm;YU2 即為計(jì)算域,其最大網(wǎng)格尺寸為 15 mm。而螺旋槳表面網(wǎng)格更小,設(shè)定最大網(wǎng)格尺寸為 2 mm。

    計(jì)算策略為:瞬態(tài)計(jì)算,計(jì)算模型選擇湍流模型,入口以及流域壁面設(shè)定為速度入口,出口設(shè)定為水流出口。

    4 計(jì)算結(jié)果分析

    4.1 數(shù)值驗(yàn)證

    首先我們對 Fluent 軟件計(jì)算螺旋槳敞水性能的準(zhǔn)確性進(jìn)行驗(yàn)證,表 1 分別為通過 Fluent 軟件計(jì)算以及船模試驗(yàn)所得的螺旋槳敞水實(shí)驗(yàn)數(shù)據(jù)。

    表1 螺旋槳敞水實(shí)驗(yàn)與計(jì)算數(shù)據(jù)對比Tab.1 Data comparison of propeller open water test and calculation

    本文使用 MRF 動坐標(biāo)系方法來模擬流動,用 Fluent 軟件來進(jìn)行流體計(jì)算。由表 1 可知,槳敞水性能的數(shù)值模擬計(jì)算結(jié)果和試驗(yàn)結(jié)果相差不大。對于扭矩系數(shù)而言無論是單個(gè)還是總體誤差均不大,其誤差均在10% 之內(nèi);推力系數(shù)誤差最大其次為推進(jìn)效率,除此之外誤差在均不大??偟膩碚f在所計(jì)算的 case 中,進(jìn)速系數(shù)在 0.5~1.1 之間,F(xiàn)luent 數(shù)值模模擬的結(jié)果與實(shí)驗(yàn)值 Kt,Kq,η 分別在 J 等于 0.833、0.9 時(shí)吻合較好;當(dāng) J = 0.5、0.7 時(shí), Kt 和 η 的誤差略大,但總的

    來說平均誤差在 10% 左右,在允許范圍內(nèi)。

    由此可以得出結(jié)論,在誤差允許范圍內(nèi),F(xiàn)luent 數(shù)值模擬得到的螺旋槳敞水性能相應(yīng)數(shù)值與實(shí)驗(yàn)值相同,即用數(shù)值模擬來求解螺旋槳的敞水性能的方法是可靠。

    4.2 堵塞效應(yīng)分析

    4.2.1 固體塊在垂直于螺旋槳軸線的面的運(yùn)動的影響

    為研究這個(gè)問題,本文特地在加 1 個(gè)和 5 個(gè)固體塊的模型中,對固體塊設(shè)置邊界條件為 moving wall,且分別設(shè)置 60 rad/s 和 125.6 rad/s 速度旋轉(zhuǎn),進(jìn)行 Fluent 計(jì)算。最后得出計(jì)算結(jié)果見表 2。在 5 個(gè)不同的case 中,不論是推力還是扭矩相差均不大,推力相差在 1% 以內(nèi),而扭矩相差也很小。即轉(zhuǎn)速的大小對其影響不大。故而可得出結(jié)論:堵塞固體塊在垂直于軸線面上速度大小在一定范圍內(nèi)對螺旋槳的水動力性能影響不大。

    表2 固體塊不同轉(zhuǎn)速下螺旋槳推力、扭矩Tab.2 Solid block at different speeds propeller thrust, torque

    4.2.2 堵塞固體塊的數(shù)量對螺旋槳的性能影響

    理論上而言,對于流場,堵塞固體塊越多對流場的影響越大,進(jìn)而對螺旋槳的推力和扭矩影響越大。本文設(shè)定 2 種不同數(shù)量的堵塞固體塊量,但堵塞固體塊大小均相同。模型 1 代指只有 1 個(gè)堵塞固體塊,模型2 代指有 5 個(gè)堵塞固體塊(具體建模見上文),其他條件保持一致,在 Fluent 中計(jì)算,具體結(jié)果如表 3 所示。

    由表 3 可知,整體呈整體呈增大趨勢,雖然計(jì)算的 case 有限,但趨勢確定。即我們可確定在一定范圍內(nèi)隨著堵塞固體塊的增多,堵塞效應(yīng)越嚴(yán)重,進(jìn)而螺旋槳所受推力以及扭矩會增大。

    表3 堵塞效應(yīng)下 2 種模型的推力以及扭矩Tab.3 Thrust and torque at blocking effect of the two models

    4.2.3 堵塞效應(yīng)對螺旋槳敞水性能的影響

    本文用模型 2(即堵塞物較多的模型)對比螺旋槳敞水下數(shù)值,推力、扭矩對比如表 4 所示。

    表4 螺旋槳在敞水及堵塞效應(yīng)下的推力、扭矩Tab.4 Thrust and torque of propeller in open water and clogging effects

    由表 4 數(shù)據(jù)已看出加了固體快的模型在其他條件相同的情況下對應(yīng)的推力和扭矩相對于普通的螺旋槳而言更大。即由 4 個(gè)計(jì)算的 case 可知堵塞作用下,螺旋槳的水動力載荷即推力和扭矩均增大,平均增量約為 15%。

    圖5 敞水螺旋槳壓力云圖Fig.5 Propeller open water stress nephogram

    圖6 堵塞作用下螺旋槳壓力云圖Fig.6 Stress nephogram of propeller under the plugging effect

    圖 5 和圖 6 分別為進(jìn)速系數(shù)為 0.7 時(shí)敞水螺旋槳的壓力云圖,由圖易得由于堵塞作用螺旋槳的槳葉葉面壓力普遍高于敞水螺旋槳的葉面壓力。與表中得出的結(jié)論一致,即堵塞效應(yīng)下的螺旋槳對應(yīng)推力、扭矩均

    有一定程度增大。本文計(jì)算的模型增大約為 15%。

    5 結(jié) 語

    1)根據(jù)螺旋槳投影原理以及原始數(shù)據(jù),利用CATIA 軟件建立 4119 槳的三維模型,實(shí)現(xiàn)了數(shù)據(jù)到模型的轉(zhuǎn)化。

    2)分析對比了數(shù)值模擬得到的數(shù)據(jù)以及實(shí)驗(yàn)得到的數(shù)據(jù),發(fā)現(xiàn)在誤差允許范圍內(nèi)數(shù)值模擬的結(jié)果與實(shí)驗(yàn)結(jié)果相同,驗(yàn)證了 CFD 數(shù)值模擬方法的可行性。并得到相應(yīng)的誤差范圍。

    3)對堵塞作用對于螺旋槳的影響做出來研究,分別研究了堵塞固體塊的速度以及數(shù)量對螺旋槳敞水性能的影響,通過對比發(fā)現(xiàn)固體塊在垂直螺旋槳旋轉(zhuǎn)軸平面的速度對螺旋槳所受推力、扭矩影響不大;固體塊的數(shù)量在一定范圍內(nèi)與螺旋槳的推力、扭矩成正比。為螺旋槳的研究提供了依據(jù)。

    [1]沈閱. 螺旋槳非正常工作狀態(tài)時(shí)的水動力性能研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2007. SHEN Yue. The research on the characteristic of propeller in abnormal condition[D]. Harbin: Harbin Engineering University, 2007.

    [2]高富東, 潘存云, 蔡汶珊, 等. 基于CFD的螺旋槳敞水性能數(shù)值分析與驗(yàn)證[J]. 機(jī)械工程學(xué)報(bào), 2010, 46(8): 133–139. GAO Fu-dong, PAN Cun-yun, CAI Wen-shan, et al. Numerical analysis and validation of propeller open-water performance based on CFD[J]. Journal of Mechanical Engineering, 2010, 46(8): 133–139.

    [3]繆宇躍, 孫江龍. CFD敞水螺旋槳性能計(jì)算分析[J]. 中國艦船研究, 2011, 6(5): 63–68. MIAO Yu-yue, SUN Jiang-long. CFD Analysis of hydrodynamic performance of propeller in open water[J]. Chinese Journal of Ships Research, 2011, 6(5): 63–68.

    [4]蔡榮泉, 陳鳳明, 馮學(xué)梅. 使用Fluent軟件的螺旋槳敞水性能計(jì)算分析[J]. 船舶力學(xué), 2006, 10(5): 41–48. CAI Rong-quan, CHEN Feng-ming, FENG Xue-mei. Calculation and analysis of the open water performance of propeller by CFD software Fluent[J]. Journal of Ship Mechanics, 2006, 10(5): 41–48.

    [5]劉帥. 潛艇操縱運(yùn)動水動力數(shù)值研究[D] 上海: 上海交通大學(xué), 2011. LIU Shuai. Numerical study of hydrodynamic forces on a maneuvering submarine[D] Shanghai: Shanghai Jiaotong University, 2011.

    [6]孫存樓, 王永生, 李堅(jiān)波. 基于CFD的調(diào)距槳調(diào)距過程槳葉水動力性能研究[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2010, 31(2): 149–153. SUN Cun-lou, WANG Yong-sheng, LI Jian-bo. Using CFD to analyze the hydrodynamic performance of controllable pitch propeller blades during pitch adjustments[J]. Journal of Harbin Engineering University, 2010, 31(2): 149–153.

    [7]王國棟. 螺旋槳水動力、空泡和噪聲性能預(yù)報(bào)方法研究[D].武漢: 華中科技大學(xué), 2013. WANG Guo-dong. Investigation on the numerical simulation of propeller hydrodynamics、cavitation and noise[D]. Wuhan: Huazhong University of Science and Technology, 2013.

    [8]王超. 螺旋槳水動力性能、空泡及噪聲性能的數(shù)值預(yù)報(bào)研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2010. WANG Chao. The research on performance of propeller’ hydrodynamics, cavitation and noise[D] Harbin: Harbin Engineering University, 2010.

    [9]SHAZLY M, PRAKASH V, LERCH B A. High strain rate compression testing of ice[M]. USA: National Aeronautics and Space Administration, Glenn Research Center, 2006.

    [10]KASTEN M. Contmllahle pitch propellers[M]. UK: Kasten Marine Design, Inc., 1998.

    Propeller open water performance calculation and blockage effect research

    QU Mu-lin1, SUN Jiang-long1,2,3, HUANG Ben-shen1, ZHONG Cheng1, MA Chao-long1
    (1. School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, China; 2. Hubei Key Laboratory of Naval Architecture & Ocean Engineering Hydrodynamics, Wuhan 430074, China; 3. Collaboration Innovation Center for Advanced Ship and Deep-Sea Exploration, Shanghai 200240, China)

    Based on the projection principle of propeller and its geometric parameters, using three-dimensional modeling software CATIA to establish a three-dimensional numerical model of propeller. Then according to the principle of computational fluid dynamics (CFD), the use of fluid dynamics software FLUENT numerical model for analysis and calculation of propeller. Using RANS method and combining with the RSM turbulence model propellers, three dimensional viscous flow field computation domain of discrete using unstructured grid method, using the relative rotating coordinates method (MRF) to simulate the movement of the propeller, in order to calculate the propeller under different into the coefficient of the flow field characteristics, and compared the numerical calculation results and experimental results of the propeller to determine the applicability of the method. Finally to study the properties of blockage effect and the ordinary propeller open water performance and propeller under the blockage effect of open water performance comparison, draw a blockage effect on propeller open water performance impact.

    propelle;CFD;MRF;open water performance;the arctic route

    U664.3

    A

    1672 – 7619(2016)11 – 0039 – 05

    10.3404/j.issn.1672 – 7619.2016.11.007

    2016 – 03 – 02;

    2016 – 04 – 05

    國家自然科學(xué)基金資助項(xiàng)目(51474109)

    瞿沐淋(1994 – ),男,碩士研究生,研究方向?yàn)榕灤畡恿W(xué)。

    猜你喜歡
    螺旋槳扭矩坐標(biāo)系
    電動汽車扭矩濾波控制策略研究
    水垂比對管柱摩阻和扭矩的影響實(shí)驗(yàn)
    中國煤層氣(2021年5期)2021-03-02 05:53:14
    基于CFD的螺旋槳拉力確定方法
    解密坐標(biāo)系中的平移變換
    坐標(biāo)系背后的故事
    液壓扭矩扳手的不確定度評定
    基于重心坐標(biāo)系的平面幾何證明的探討
    氣密封接頭最佳上扣扭矩計(jì)算與分析
    3800DWT加油船螺旋槳諧鳴分析及消除方法
    廣東造船(2015年6期)2015-02-27 10:52:46
    螺旋槳轂帽鰭節(jié)能性能的數(shù)值模擬
    亚洲,欧美,日韩| a级片在线免费高清观看视频| 一区二区三区四区激情视频| 欧美日韩福利视频一区二区| 国产女主播在线喷水免费视频网站| 精品国产一区二区三区四区第35| www.自偷自拍.com| 国产亚洲av片在线观看秒播厂| 丝瓜视频免费看黄片| 在线亚洲精品国产二区图片欧美| 悠悠久久av| av线在线观看网站| 少妇精品久久久久久久| 国产野战对白在线观看| 欧美av亚洲av综合av国产av| 久久亚洲精品不卡| 亚洲成人国产一区在线观看 | 中文字幕人妻熟女乱码| 亚洲免费av在线视频| 男人操女人黄网站| 亚洲七黄色美女视频| 国产精品av久久久久免费| 国产有黄有色有爽视频| 欧美在线一区亚洲| 一级,二级,三级黄色视频| 国产女主播在线喷水免费视频网站| 交换朋友夫妻互换小说| 欧美激情高清一区二区三区| 国产高清不卡午夜福利| 黄频高清免费视频| 后天国语完整版免费观看| 国产人伦9x9x在线观看| 黄网站色视频无遮挡免费观看| 欧美变态另类bdsm刘玥| 亚洲五月婷婷丁香| 高潮久久久久久久久久久不卡| 久久人人爽人人片av| 日韩欧美一区视频在线观看| 亚洲一区二区三区欧美精品| 男女边吃奶边做爰视频| 亚洲五月色婷婷综合| 日韩精品免费视频一区二区三区| 丝袜人妻中文字幕| 校园人妻丝袜中文字幕| 婷婷成人精品国产| 老汉色av国产亚洲站长工具| 日本午夜av视频| 色视频在线一区二区三区| 一边摸一边做爽爽视频免费| 国产又色又爽无遮挡免| 青草久久国产| 中国国产av一级| 久久人人爽人人片av| 热99久久久久精品小说推荐| 午夜av观看不卡| 国产精品一区二区免费欧美 | 两个人免费观看高清视频| 大话2 男鬼变身卡| 男女高潮啪啪啪动态图| 美女国产高潮福利片在线看| 欧美成人午夜精品| 亚洲av综合色区一区| 麻豆乱淫一区二区| 免费在线观看视频国产中文字幕亚洲 | 亚洲激情五月婷婷啪啪| 欧美日韩成人在线一区二区| 嫁个100分男人电影在线观看 | 每晚都被弄得嗷嗷叫到高潮| 另类精品久久| 欧美 日韩 精品 国产| 性色av一级| 国产野战对白在线观看| 国产91精品成人一区二区三区 | 午夜91福利影院| 亚洲精品第二区| 欧美97在线视频| 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区| 99热全是精品| 久久精品国产亚洲av涩爱| 日本午夜av视频| 啦啦啦视频在线资源免费观看| 三上悠亚av全集在线观看| 日本黄色日本黄色录像| 桃花免费在线播放| 午夜激情av网站| 国产午夜精品一二区理论片| 观看av在线不卡| 男男h啪啪无遮挡| 十八禁高潮呻吟视频| 人体艺术视频欧美日本| 精品欧美一区二区三区在线| 日韩一区二区三区影片| 免费在线观看完整版高清| 亚洲熟女精品中文字幕| 真人做人爱边吃奶动态| 国产野战对白在线观看| 婷婷色麻豆天堂久久| 亚洲av电影在线观看一区二区三区| 超碰97精品在线观看| 亚洲av美国av| 电影成人av| 99久久精品国产亚洲精品| 大片电影免费在线观看免费| 激情五月婷婷亚洲| 精品一区二区三区av网在线观看 | 亚洲专区国产一区二区| 亚洲av在线观看美女高潮| 亚洲视频免费观看视频| 久久女婷五月综合色啪小说| 亚洲精品自拍成人| 亚洲欧洲国产日韩| 久久久精品94久久精品| 亚洲自偷自拍图片 自拍| 久热爱精品视频在线9| 国产真人三级小视频在线观看| a级毛片在线看网站| 成人18禁高潮啪啪吃奶动态图| 性色av一级| 五月开心婷婷网| 色婷婷av一区二区三区视频| 国产片特级美女逼逼视频| 两个人免费观看高清视频| 国产色视频综合| 亚洲 欧美一区二区三区| 99久久人妻综合| 性色av乱码一区二区三区2| 成人三级做爰电影| 久久精品久久久久久久性| 九草在线视频观看| av又黄又爽大尺度在线免费看| 欧美日韩精品网址| 看十八女毛片水多多多| 亚洲欧洲日产国产| 18禁观看日本| 久久久国产精品麻豆| 日本午夜av视频| 亚洲国产欧美在线一区| av国产精品久久久久影院| 最近手机中文字幕大全| 在线观看www视频免费| 美女中出高潮动态图| 看十八女毛片水多多多| 婷婷色麻豆天堂久久| 国产精品人妻久久久影院| 日日摸夜夜添夜夜爱| 狂野欧美激情性bbbbbb| 男女午夜视频在线观看| 久久久久精品人妻al黑| 亚洲黑人精品在线| 在线av久久热| 国产精品久久久久久精品古装| 久久精品人人爽人人爽视色| 国产高清videossex| 999精品在线视频| 久久国产精品大桥未久av| xxx大片免费视频| 亚洲精品第二区| 国产在线免费精品| svipshipincom国产片| 色婷婷久久久亚洲欧美| 亚洲精品国产av蜜桃| 免费看av在线观看网站| 国产精品麻豆人妻色哟哟久久| 一二三四社区在线视频社区8| 交换朋友夫妻互换小说| 久久久久久亚洲精品国产蜜桃av| 看免费av毛片| 狂野欧美激情性bbbbbb| 搡老岳熟女国产| 成在线人永久免费视频| 亚洲美女黄色视频免费看| 男女边吃奶边做爰视频| 欧美 亚洲 国产 日韩一| 欧美+亚洲+日韩+国产| 国产一区二区三区综合在线观看| 狂野欧美激情性xxxx| 午夜两性在线视频| 欧美激情高清一区二区三区| 看免费av毛片| 欧美变态另类bdsm刘玥| 国产成人91sexporn| 亚洲欧美中文字幕日韩二区| 久久女婷五月综合色啪小说| 两个人看的免费小视频| 午夜影院在线不卡| 99精品久久久久人妻精品| 超碰成人久久| 性高湖久久久久久久久免费观看| 久久午夜综合久久蜜桃| 欧美av亚洲av综合av国产av| 日本一区二区免费在线视频| 国产亚洲av高清不卡| 亚洲成色77777| 国产激情久久老熟女| 999精品在线视频| 久久人妻熟女aⅴ| 国产精品一区二区免费欧美 | 午夜福利在线免费观看网站| 欧美久久黑人一区二区| 久久久国产欧美日韩av| 婷婷色综合大香蕉| 一级黄片播放器| 久久精品人人爽人人爽视色| 黑人欧美特级aaaaaa片| 精品亚洲成国产av| 国产黄色免费在线视频| 亚洲精品日本国产第一区| 亚洲国产日韩一区二区| 国产片内射在线| av线在线观看网站| 国产福利在线免费观看视频| 可以免费在线观看a视频的电影网站| 久久国产亚洲av麻豆专区| 丝袜脚勾引网站| 日韩欧美一区视频在线观看| 国产极品粉嫩免费观看在线| 97人妻天天添夜夜摸| 国产欧美亚洲国产| 国产片内射在线| 精品亚洲乱码少妇综合久久| 汤姆久久久久久久影院中文字幕| 老司机影院成人| 一二三四社区在线视频社区8| 国产一区亚洲一区在线观看| 久久精品国产a三级三级三级| 日韩av在线免费看完整版不卡| 好男人视频免费观看在线| 欧美变态另类bdsm刘玥| 亚洲av男天堂| 久久久久久久国产电影| 久久精品亚洲av国产电影网| 免费在线观看视频国产中文字幕亚洲 | 啦啦啦 在线观看视频| 亚洲国产欧美在线一区| 午夜精品国产一区二区电影| 视频在线观看一区二区三区| 国产精品熟女久久久久浪| 日韩,欧美,国产一区二区三区| 一级片'在线观看视频| 中文字幕制服av| 一级片免费观看大全| 在线观看免费高清a一片| av视频免费观看在线观看| 日韩一卡2卡3卡4卡2021年| 日本91视频免费播放| 亚洲色图 男人天堂 中文字幕| 免费在线观看视频国产中文字幕亚洲 | 午夜福利在线免费观看网站| 啦啦啦中文免费视频观看日本| 午夜福利乱码中文字幕| av又黄又爽大尺度在线免费看| 中文字幕制服av| 视频在线观看一区二区三区| 下体分泌物呈黄色| 国产高清不卡午夜福利| 亚洲精品久久午夜乱码| 亚洲,欧美精品.| 婷婷丁香在线五月| 一个人免费看片子| 一区二区av电影网| 免费在线观看影片大全网站 | 欧美 亚洲 国产 日韩一| 亚洲国产日韩一区二区| 天堂8中文在线网| 国产真人三级小视频在线观看| 国产极品粉嫩免费观看在线| 黄色 视频免费看| 国产精品久久久久久人妻精品电影 | 日韩伦理黄色片| 久久久久国产精品人妻一区二区| 男女免费视频国产| 亚洲精品国产av成人精品| 国产精品九九99| 欧美精品av麻豆av| 成人午夜精彩视频在线观看| 日韩大片免费观看网站| 十八禁高潮呻吟视频| 波多野结衣一区麻豆| 国产精品久久久久久精品电影小说| 天天添夜夜摸| 看免费av毛片| 欧美精品高潮呻吟av久久| 亚洲国产欧美一区二区综合| 国产精品三级大全| 久热爱精品视频在线9| 人妻 亚洲 视频| 亚洲av国产av综合av卡| 高潮久久久久久久久久久不卡| 嫩草影视91久久| 成人免费观看视频高清| 免费不卡黄色视频| 亚洲精品日本国产第一区| 精品视频人人做人人爽| 国产精品成人在线| 国产在线视频一区二区| 精品国产乱码久久久久久小说| av国产久精品久网站免费入址| 久久国产精品男人的天堂亚洲| 麻豆av在线久日| 2021少妇久久久久久久久久久| av有码第一页| 啦啦啦视频在线资源免费观看| 亚洲国产精品999| www.精华液| 亚洲成人免费av在线播放| 中文字幕人妻熟女乱码| 亚洲欧美中文字幕日韩二区| 宅男免费午夜| 亚洲图色成人| 国产av精品麻豆| 欧美少妇被猛烈插入视频| 国产精品久久久久久精品古装| 女性生殖器流出的白浆| 菩萨蛮人人尽说江南好唐韦庄| 国产精品熟女久久久久浪| 18禁国产床啪视频网站| 视频区图区小说| 日本欧美国产在线视频| 国产不卡av网站在线观看| 中文字幕人妻熟女乱码| 不卡av一区二区三区| 欧美老熟妇乱子伦牲交| 少妇人妻 视频| 亚洲欧美中文字幕日韩二区| 久久久久精品人妻al黑| 一级毛片 在线播放| 中文乱码字字幕精品一区二区三区| 亚洲男人天堂网一区| 一级毛片电影观看| 国产欧美日韩一区二区三 | 美女主播在线视频| 国产亚洲一区二区精品| 久久精品亚洲av国产电影网| 超碰97精品在线观看| 欧美日韩成人在线一区二区| 性色av一级| 超色免费av| 又紧又爽又黄一区二区| 嫁个100分男人电影在线观看 | 国产成人av教育| 国产高清视频在线播放一区 | 久久99热这里只频精品6学生| 国产成人一区二区在线| 亚洲精品久久午夜乱码| 黄色视频不卡| 50天的宝宝边吃奶边哭怎么回事| 一二三四社区在线视频社区8| 久久人人爽av亚洲精品天堂| 国产欧美日韩精品亚洲av| 两个人看的免费小视频| 一二三四社区在线视频社区8| 久久综合国产亚洲精品| 中文字幕人妻丝袜一区二区| av国产精品久久久久影院| 久久久久精品人妻al黑| 真人做人爱边吃奶动态| 国产高清不卡午夜福利| 欧美人与性动交α欧美软件| 国产成人精品久久二区二区91| 色婷婷久久久亚洲欧美| 国产精品 国内视频| 午夜激情av网站| 交换朋友夫妻互换小说| 侵犯人妻中文字幕一二三四区| 欧美精品人与动牲交sv欧美| 一边摸一边抽搐一进一出视频| 男女下面插进去视频免费观看| 亚洲自偷自拍图片 自拍| 天堂中文最新版在线下载| 中文字幕精品免费在线观看视频| 欧美xxⅹ黑人| 欧美精品亚洲一区二区| 国产成人影院久久av| 亚洲欧美激情在线| 国产黄频视频在线观看| 999精品在线视频| 这个男人来自地球电影免费观看| 亚洲欧美清纯卡通| 日本91视频免费播放| 久久99一区二区三区| 成人影院久久| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品久久久久久| 99热网站在线观看| 中文乱码字字幕精品一区二区三区| 午夜免费观看性视频| 99国产精品一区二区蜜桃av | 中文字幕亚洲精品专区| 纵有疾风起免费观看全集完整版| 国产1区2区3区精品| 午夜福利乱码中文字幕| 老司机午夜十八禁免费视频| 国产视频一区二区在线看| 超碰97精品在线观看| 国产爽快片一区二区三区| 亚洲精品自拍成人| 91精品伊人久久大香线蕉| kizo精华| 爱豆传媒免费全集在线观看| 久久久久久人人人人人| 久久久亚洲精品成人影院| 一本久久精品| 在线精品无人区一区二区三| 欧美日韩亚洲国产一区二区在线观看 | 久久热在线av| 熟女av电影| 91老司机精品| 日韩熟女老妇一区二区性免费视频| 亚洲av片天天在线观看| 中文字幕另类日韩欧美亚洲嫩草| 下体分泌物呈黄色| 国产成人精品无人区| 国产免费一区二区三区四区乱码| 亚洲国产日韩一区二区| 国产精品人妻久久久影院| 久久久久久久大尺度免费视频| 视频在线观看一区二区三区| 免费一级毛片在线播放高清视频 | av福利片在线| 午夜两性在线视频| 日本a在线网址| 亚洲中文日韩欧美视频| 中文字幕人妻丝袜一区二区| 国产精品欧美亚洲77777| 国产亚洲av片在线观看秒播厂| 欧美日韩一级在线毛片| 国产视频一区二区在线看| 亚洲七黄色美女视频| 大片免费播放器 马上看| 亚洲,欧美,日韩| 国产男女内射视频| 欧美av亚洲av综合av国产av| 国产一区亚洲一区在线观看| 久久久欧美国产精品| 91老司机精品| 免费高清在线观看日韩| 男人爽女人下面视频在线观看| 免费看十八禁软件| 久久久久久久大尺度免费视频| 国产一区二区在线观看av| 国产亚洲av片在线观看秒播厂| 蜜桃在线观看..| 精品一区在线观看国产| 欧美日本中文国产一区发布| 久久精品国产综合久久久| 亚洲国产精品一区二区三区在线| 精品少妇一区二区三区视频日本电影| 日本午夜av视频| 人人妻,人人澡人人爽秒播 | 女人高潮潮喷娇喘18禁视频| 美女扒开内裤让男人捅视频| 黄片小视频在线播放| 一级毛片女人18水好多 | 精品福利观看| 午夜av观看不卡| 乱人伦中国视频| 国产男人的电影天堂91| 久久久久久久精品精品| 亚洲av电影在线进入| 日本欧美国产在线视频| 久久99精品国语久久久| 亚洲av成人不卡在线观看播放网 | 亚洲国产欧美网| 精品免费久久久久久久清纯 | 日本一区二区免费在线视频| 搡老岳熟女国产| 亚洲精品国产色婷婷电影| 精品少妇黑人巨大在线播放| 好男人视频免费观看在线| 男女免费视频国产| 久久精品久久久久久久性| 国产成人av教育| 免费少妇av软件| 制服人妻中文乱码| 波野结衣二区三区在线| 欧美日韩精品网址| 69精品国产乱码久久久| 丁香六月天网| 国产高清国产精品国产三级| 美女高潮到喷水免费观看| 美国免费a级毛片| 日韩人妻精品一区2区三区| 婷婷成人精品国产| 七月丁香在线播放| 精品一品国产午夜福利视频| 久久精品久久精品一区二区三区| 久久人妻福利社区极品人妻图片 | 99热全是精品| 日本一区二区免费在线视频| 精品国产一区二区三区四区第35| 黄片播放在线免费| 国产成人精品在线电影| a 毛片基地| 久9热在线精品视频| 国产亚洲精品第一综合不卡| 国产人伦9x9x在线观看| 国产免费现黄频在线看| 午夜激情av网站| 国产熟女午夜一区二区三区| tube8黄色片| 搡老乐熟女国产| 国产欧美日韩精品亚洲av| 亚洲欧美清纯卡通| 久久99精品国语久久久| 黑人巨大精品欧美一区二区蜜桃| 亚洲伊人久久精品综合| 午夜福利视频在线观看免费| 多毛熟女@视频| 国产三级黄色录像| 国产精品久久久久久精品电影小说| 亚洲熟女精品中文字幕| 国产精品香港三级国产av潘金莲 | 国产一区二区三区av在线| 国产精品久久久久成人av| 乱人伦中国视频| 国产欧美日韩一区二区三 | 蜜桃在线观看..| 亚洲国产中文字幕在线视频| 久久久久久人人人人人| 母亲3免费完整高清在线观看| 中文乱码字字幕精品一区二区三区| cao死你这个sao货| 国产在视频线精品| 美女脱内裤让男人舔精品视频| av电影中文网址| 你懂的网址亚洲精品在线观看| 一本综合久久免费| 丝袜美腿诱惑在线| 王馨瑶露胸无遮挡在线观看| 久久国产精品大桥未久av| 啦啦啦在线免费观看视频4| 1024香蕉在线观看| 亚洲欧洲日产国产| 丝瓜视频免费看黄片| 日韩av不卡免费在线播放| 97精品久久久久久久久久精品| 色精品久久人妻99蜜桃| 一级毛片我不卡| 亚洲av日韩精品久久久久久密 | 在线亚洲精品国产二区图片欧美| 成年人黄色毛片网站| 亚洲精品第二区| 亚洲国产欧美一区二区综合| 一二三四在线观看免费中文在| 这个男人来自地球电影免费观看| 欧美国产精品va在线观看不卡| 久久人妻熟女aⅴ| 日日爽夜夜爽网站| 黑丝袜美女国产一区| 午夜91福利影院| 免费观看av网站的网址| 五月天丁香电影| 日本a在线网址| 精品卡一卡二卡四卡免费| 啦啦啦 在线观看视频| 免费在线观看日本一区| 午夜免费男女啪啪视频观看| 欧美中文综合在线视频| 日韩欧美一区视频在线观看| 国产精品三级大全| 成年动漫av网址| 在线观看国产h片| 男人添女人高潮全过程视频| 亚洲人成电影免费在线| 亚洲精品中文字幕在线视频| 天天躁夜夜躁狠狠久久av| 中文字幕精品免费在线观看视频| 黄片播放在线免费| 亚洲国产欧美在线一区| 日本91视频免费播放| 在线观看免费午夜福利视频| 飞空精品影院首页| 观看av在线不卡| 999久久久国产精品视频| 亚洲成国产人片在线观看| 在线观看一区二区三区激情| 在线av久久热| 男女无遮挡免费网站观看| 大码成人一级视频| 丝袜喷水一区| 黄色片一级片一级黄色片| 黑人巨大精品欧美一区二区蜜桃| 少妇的丰满在线观看| 久久久久久久大尺度免费视频| 免费不卡黄色视频| 在线观看一区二区三区激情| 久久久精品94久久精品| 国产一区二区三区综合在线观看| 丝袜美足系列| 亚洲一码二码三码区别大吗| 亚洲欧美日韩另类电影网站| 亚洲色图 男人天堂 中文字幕| 国产深夜福利视频在线观看| 中文字幕色久视频| av欧美777| 50天的宝宝边吃奶边哭怎么回事| 国产三级黄色录像| 老司机亚洲免费影院| 一区二区三区精品91| 在线亚洲精品国产二区图片欧美| 亚洲成国产人片在线观看| 久久久欧美国产精品| 看免费av毛片| 一级毛片黄色毛片免费观看视频| 人妻人人澡人人爽人人| 亚洲九九香蕉| 欧美 亚洲 国产 日韩一| 电影成人av| 国产片特级美女逼逼视频| 国产精品 国内视频| 黄色 视频免费看| 国产成人欧美在线观看 | 精品人妻在线不人妻| 男女国产视频网站|