• <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ù)值模擬
    亚洲,一卡二卡三卡| 久久久久性生活片| 亚洲欧美中文字幕日韩二区| 国产人妻一区二区三区在| 国产 精品1| 插阴视频在线观看视频| 久久久久网色| 91久久精品国产一区二区成人| 偷拍熟女少妇极品色| 国产精品一区二区在线不卡| 亚洲精品色激情综合| 欧美区成人在线视频| 91精品国产九色| 久热久热在线精品观看| 久久久久久人妻| 国产高清三级在线| 亚洲av成人精品一区久久| 国产在线一区二区三区精| 99热这里只有精品一区| 亚洲美女搞黄在线观看| 久久久久久人妻| 中文字幕制服av| 精品亚洲乱码少妇综合久久| 亚洲av在线观看美女高潮| 大又大粗又爽又黄少妇毛片口| 精品一品国产午夜福利视频| 亚洲精品一区蜜桃| videos熟女内射| 亚洲精品国产av蜜桃| 中文字幕亚洲精品专区| 精品人妻一区二区三区麻豆| 97在线人人人人妻| 老师上课跳d突然被开到最大视频| av线在线观看网站| av在线老鸭窝| 青春草视频在线免费观看| a级一级毛片免费在线观看| 精品一区二区三卡| 极品教师在线视频| 国产深夜福利视频在线观看| 欧美一级a爱片免费观看看| 国产欧美亚洲国产| av福利片在线观看| 不卡视频在线观看欧美| 久久久久久久久久久免费av| 搡女人真爽免费视频火全软件| 久久精品久久精品一区二区三区| 久久久久久人妻| 伦理电影大哥的女人| 精华霜和精华液先用哪个| 18禁在线播放成人免费| 久久99热6这里只有精品| 日韩一区二区视频免费看| 夫妻性生交免费视频一级片| 日韩三级伦理在线观看| 美女cb高潮喷水在线观看| 99久久人妻综合| av天堂中文字幕网| 亚洲国产欧美在线一区| 免费黄色在线免费观看| 内射极品少妇av片p| 日本与韩国留学比较| 我要看日韩黄色一级片| 极品教师在线视频| 三级经典国产精品| 国产成人精品婷婷| 韩国av在线不卡| 欧美高清成人免费视频www| 日韩av在线免费看完整版不卡| 麻豆国产97在线/欧美| 亚洲国产精品成人久久小说| 蜜桃在线观看..| 婷婷色麻豆天堂久久| 我要看日韩黄色一级片| 蜜桃久久精品国产亚洲av| 老司机影院毛片| 成人国产av品久久久| 2022亚洲国产成人精品| 国产视频首页在线观看| 蜜桃亚洲精品一区二区三区| 老司机影院毛片| 一级毛片 在线播放| 日韩免费高清中文字幕av| 国产 一区 欧美 日韩| 国产精品久久久久久久电影| 一级a做视频免费观看| 国产精品99久久99久久久不卡 | 老熟女久久久| 国产av码专区亚洲av| av免费在线看不卡| 乱系列少妇在线播放| 国产av一区二区精品久久 | 精品久久久精品久久久| h视频一区二区三区| 中文字幕av成人在线电影| 黑人猛操日本美女一级片| 国产乱来视频区| 麻豆成人av视频| 哪个播放器可以免费观看大片| 91久久精品国产一区二区成人| 久久久久久久久久久丰满| 大陆偷拍与自拍| 免费不卡的大黄色大毛片视频在线观看| 中文字幕精品免费在线观看视频 | 国产精品女同一区二区软件| 老熟女久久久| 亚洲国产高清在线一区二区三| 亚洲自偷自拍三级| 啦啦啦啦在线视频资源| 成人亚洲欧美一区二区av| av免费在线看不卡| 国产视频内射| 亚洲成人av在线免费| 亚洲一级一片aⅴ在线观看| 亚洲四区av| 亚洲丝袜综合中文字幕| 亚洲天堂av无毛| 少妇猛男粗大的猛烈进出视频| 老熟女久久久| 亚洲精品日韩在线中文字幕| 国产精品久久久久久av不卡| 成人一区二区视频在线观看| 欧美日韩视频精品一区| 夜夜骑夜夜射夜夜干| 亚洲三级黄色毛片| 亚洲熟女精品中文字幕| 国产精品久久久久成人av| 日韩视频在线欧美| 精品久久久久久久久亚洲| 国精品久久久久久国模美| 国产深夜福利视频在线观看| 人妻夜夜爽99麻豆av| 成人影院久久| 直男gayav资源| av.在线天堂| 五月玫瑰六月丁香| 三级国产精品片| 精品酒店卫生间| 我的女老师完整版在线观看| 秋霞在线观看毛片| 国产男女超爽视频在线观看| 舔av片在线| 国产爽快片一区二区三区| 一级毛片黄色毛片免费观看视频| 99九九线精品视频在线观看视频| 亚洲精品一区蜜桃| 免费人妻精品一区二区三区视频| 国产精品精品国产色婷婷| 美女视频免费永久观看网站| 亚洲欧美精品自产自拍| 色视频在线一区二区三区| 欧美高清性xxxxhd video| 精品熟女少妇av免费看| 中国国产av一级| 欧美成人午夜免费资源| 久久热精品热| 美女内射精品一级片tv| 欧美日本视频| 视频区图区小说| 男人和女人高潮做爰伦理| 国产成人午夜福利电影在线观看| 亚洲怡红院男人天堂| 亚洲va在线va天堂va国产| 一级毛片aaaaaa免费看小| av在线观看视频网站免费| 日日摸夜夜添夜夜爱| 国产日韩欧美在线精品| 久久久午夜欧美精品| 国产av国产精品国产| 在线观看人妻少妇| 亚洲精品亚洲一区二区| 王馨瑶露胸无遮挡在线观看| h视频一区二区三区| 久久国产乱子免费精品| 老司机影院成人| 伦精品一区二区三区| 亚洲综合色惰| 妹子高潮喷水视频| 欧美精品一区二区大全| 一二三四中文在线观看免费高清| 成年女人在线观看亚洲视频| 免费播放大片免费观看视频在线观看| 狂野欧美白嫩少妇大欣赏| 久久鲁丝午夜福利片| 一本—道久久a久久精品蜜桃钙片| 日韩 亚洲 欧美在线| 一个人看视频在线观看www免费| 精品一区二区免费观看| 国产视频内射| 久久久久久久久久人人人人人人| 久久人人爽av亚洲精品天堂 | 在线观看三级黄色| 国产欧美日韩一区二区三区在线 | 看十八女毛片水多多多| av一本久久久久| 香蕉精品网在线| 欧美三级亚洲精品| 久久6这里有精品| 久久久久久久亚洲中文字幕| 精品久久久久久久久亚洲| 国产精品嫩草影院av在线观看| 熟女av电影| 精品一区二区免费观看| 亚洲精品456在线播放app| 日本-黄色视频高清免费观看| 国产精品女同一区二区软件| 亚洲三级黄色毛片| 一区二区av电影网| 三级国产精品片| 熟女电影av网| 久久国产亚洲av麻豆专区| 联通29元200g的流量卡| 精品国产三级普通话版| 性色av一级| 免费av不卡在线播放| 国产视频首页在线观看| 免费少妇av软件| 亚洲欧美日韩无卡精品| 欧美精品亚洲一区二区| 国产av码专区亚洲av| 国产精品伦人一区二区| 午夜福利高清视频| 亚洲综合精品二区| 免费大片黄手机在线观看| 国产国拍精品亚洲av在线观看| 美女视频免费永久观看网站| 国产乱人视频| 一级毛片久久久久久久久女| 免费黄网站久久成人精品| 身体一侧抽搐| 一级毛片电影观看| 亚洲美女搞黄在线观看| 国产 一区精品| 亚洲av免费高清在线观看| 久久久久久久久久成人| 激情五月婷婷亚洲| 免费少妇av软件| 男女下面进入的视频免费午夜| 丝袜脚勾引网站| 黄色日韩在线| 天美传媒精品一区二区| 18+在线观看网站| 国产成人免费观看mmmm| 国产欧美日韩一区二区三区在线 | 亚洲欧美精品自产自拍| 日日摸夜夜添夜夜爱| av视频免费观看在线观看| 久久久久久久久大av| 人妻系列 视频| 午夜福利在线观看免费完整高清在| 人人妻人人添人人爽欧美一区卜 | 亚洲av.av天堂| 干丝袜人妻中文字幕| 1000部很黄的大片| 亚洲综合精品二区| av天堂中文字幕网| 少妇猛男粗大的猛烈进出视频| 久久久欧美国产精品| 一本色道久久久久久精品综合| 亚洲精品中文字幕在线视频 | 深夜a级毛片| 女性被躁到高潮视频| 伦理电影大哥的女人| 久久ye,这里只有精品| 亚洲av成人精品一区久久| 99热国产这里只有精品6| 国产美女午夜福利| 熟女电影av网| 成人无遮挡网站| 亚洲在久久综合| 黑人猛操日本美女一级片| 最近最新中文字幕大全电影3| 欧美日韩在线观看h| 成年人午夜在线观看视频| 网址你懂的国产日韩在线| 日韩强制内射视频| av专区在线播放| 色婷婷av一区二区三区视频| 亚洲精品乱码久久久v下载方式| 欧美亚洲 丝袜 人妻 在线| 最近手机中文字幕大全| 成人亚洲欧美一区二区av| 最近中文字幕高清免费大全6| 麻豆国产97在线/欧美| 国产在线一区二区三区精| 欧美日韩视频精品一区| 国产69精品久久久久777片| 日韩欧美 国产精品| 国产免费一区二区三区四区乱码| 性高湖久久久久久久久免费观看| 久久97久久精品| 日日摸夜夜添夜夜添av毛片| 777米奇影视久久| 亚洲人成网站在线播| 亚洲真实伦在线观看| 老司机影院成人| 欧美区成人在线视频| 免费观看无遮挡的男女| www.av在线官网国产| 一级二级三级毛片免费看| 日韩电影二区| 国产男人的电影天堂91| 国产高清三级在线| 观看av在线不卡| 午夜激情福利司机影院| 国产永久视频网站| 在线免费十八禁| 人体艺术视频欧美日本| 好男人视频免费观看在线| 亚洲最大成人中文| 一区二区三区四区激情视频| 91狼人影院| 男女边吃奶边做爰视频| 午夜免费观看性视频| 国产一区亚洲一区在线观看| 亚洲国产毛片av蜜桃av| 五月玫瑰六月丁香| 成人综合一区亚洲| 亚洲av福利一区| 最近最新中文字幕大全电影3| 亚洲美女视频黄频| 欧美高清性xxxxhd video| av国产免费在线观看| 久久人人爽人人片av| 永久免费av网站大全| 毛片女人毛片| 亚洲国产高清在线一区二区三| 精品久久久久久久末码| 亚洲成人av在线免费| 久久6这里有精品| 十八禁网站网址无遮挡 | 天堂俺去俺来也www色官网| 91午夜精品亚洲一区二区三区| 日本欧美国产在线视频| 成年美女黄网站色视频大全免费 | 久久99热这里只频精品6学生| 久久久欧美国产精品| 亚洲欧洲日产国产| 久久毛片免费看一区二区三区| 精品久久国产蜜桃| 少妇的逼水好多| 一级片'在线观看视频| 最近最新中文字幕大全电影3| 国产成人91sexporn| 欧美高清性xxxxhd video| 精品国产乱码久久久久久小说| 久久精品久久久久久久性| 2022亚洲国产成人精品| 51国产日韩欧美| 日韩欧美 国产精品| 女性被躁到高潮视频| 永久网站在线| 下体分泌物呈黄色| 在线观看av片永久免费下载| 狂野欧美激情性bbbbbb| av网站免费在线观看视频| 国产精品久久久久久精品电影小说 | 美女主播在线视频| 熟妇人妻不卡中文字幕| 久久99热这里只有精品18| av线在线观看网站| 一本—道久久a久久精品蜜桃钙片| 涩涩av久久男人的天堂| 五月玫瑰六月丁香| 国产欧美日韩精品一区二区| 日韩强制内射视频| 欧美人与善性xxx| 青青草视频在线视频观看| 亚洲国产精品999| 国产精品精品国产色婷婷| 久久亚洲国产成人精品v| 噜噜噜噜噜久久久久久91| 欧美变态另类bdsm刘玥| 国产高清国产精品国产三级 | 六月丁香七月| 国产精品久久久久久久电影| 久久久久网色| 亚洲精品日韩在线中文字幕| 国产69精品久久久久777片| 国产午夜精品久久久久久一区二区三区| 久久人人爽人人片av| 国产极品天堂在线| 亚洲精品456在线播放app| 精华霜和精华液先用哪个| 18禁动态无遮挡网站| 精品久久久久久久久亚洲| 午夜精品国产一区二区电影| 国产成人91sexporn| 乱码一卡2卡4卡精品| 日日撸夜夜添| 男人添女人高潮全过程视频| av一本久久久久| 91精品一卡2卡3卡4卡| 搡老乐熟女国产| 亚洲成人手机| 男男h啪啪无遮挡| 国产91av在线免费观看| 看免费成人av毛片| 国产精品秋霞免费鲁丝片| 欧美一区二区亚洲| 精品亚洲成a人片在线观看 | 超碰av人人做人人爽久久| 久久婷婷青草| 国产精品av视频在线免费观看| 啦啦啦啦在线视频资源| 中文天堂在线官网| av卡一久久| 91狼人影院| 亚洲精品国产色婷婷电影| tube8黄色片| 99热国产这里只有精品6| 水蜜桃什么品种好| 大陆偷拍与自拍| 日本av手机在线免费观看| 国产在线男女| 亚州av有码| 亚洲成人一二三区av| 在线观看av片永久免费下载| 国产精品熟女久久久久浪| 亚洲,欧美,日韩| 国产精品欧美亚洲77777| av线在线观看网站| 久久久久久伊人网av| 99热这里只有是精品50| 丰满乱子伦码专区| 日本黄色日本黄色录像| 午夜老司机福利剧场| 日本黄大片高清| 在线 av 中文字幕| 午夜日本视频在线| 欧美成人精品欧美一级黄| 日韩伦理黄色片| 新久久久久国产一级毛片| 国产欧美另类精品又又久久亚洲欧美| 联通29元200g的流量卡| 噜噜噜噜噜久久久久久91| 亚洲在久久综合| 日日摸夜夜添夜夜爱| 免费av不卡在线播放| 午夜福利影视在线免费观看| 日韩欧美 国产精品| 日韩 亚洲 欧美在线| 亚洲精品久久久久久婷婷小说| 天堂俺去俺来也www色官网| 深夜a级毛片| 国产欧美日韩一区二区三区在线 | 一个人看的www免费观看视频| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在| 男人和女人高潮做爰伦理| 国产女主播在线喷水免费视频网站| 国产精品99久久99久久久不卡 | 国产无遮挡羞羞视频在线观看| av播播在线观看一区| 亚洲国产成人一精品久久久| 亚洲av男天堂| 成年女人在线观看亚洲视频| 少妇人妻 视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲怡红院男人天堂| tube8黄色片| a级毛片免费高清观看在线播放| 欧美激情国产日韩精品一区| 亚洲婷婷狠狠爱综合网| 亚洲欧美精品自产自拍| 亚洲精品第二区| 欧美日韩视频精品一区| 伊人久久国产一区二区| 日本猛色少妇xxxxx猛交久久| 亚洲人成网站在线观看播放| 国产av国产精品国产| 久久久久久久久久人人人人人人| 亚洲精品日本国产第一区| 久久久午夜欧美精品| 国产精品久久久久久久久免| 亚洲在久久综合| 99九九线精品视频在线观看视频| 十八禁网站网址无遮挡 | 国产精品熟女久久久久浪| 交换朋友夫妻互换小说| 99热这里只有是精品在线观看| 日日啪夜夜撸| 日本黄色片子视频| 国产精品av视频在线免费观看| 亚洲人成网站高清观看| 狂野欧美激情性xxxx在线观看| 人人妻人人爽人人添夜夜欢视频 | 久久人人爽av亚洲精品天堂 | 亚洲四区av| 国产永久视频网站| 免费在线观看成人毛片| 视频区图区小说| 久久久久久久亚洲中文字幕| 国产精品免费大片| 中文资源天堂在线| 亚洲aⅴ乱码一区二区在线播放| 天堂中文最新版在线下载| 久久久久精品性色| 免费观看a级毛片全部| 免费黄网站久久成人精品| 色吧在线观看| 国产日韩欧美亚洲二区| 欧美国产精品一级二级三级 | 三级国产精品欧美在线观看| 最近2019中文字幕mv第一页| 亚洲天堂av无毛| 内地一区二区视频在线| 午夜激情久久久久久久| 国产av一区二区精品久久 | 青春草视频在线免费观看| 美女xxoo啪啪120秒动态图| 三级经典国产精品| 在线免费十八禁| 日本午夜av视频| 国产亚洲av片在线观看秒播厂| 大片电影免费在线观看免费| 在线精品无人区一区二区三 | 久久久久久久久久人人人人人人| 欧美精品亚洲一区二区| 欧美日韩精品成人综合77777| 丝袜喷水一区| 国产视频首页在线观看| 日韩,欧美,国产一区二区三区| 美女内射精品一级片tv| av在线蜜桃| 日韩 亚洲 欧美在线| 亚洲精品日韩av片在线观看| 国产高清三级在线| 一边亲一边摸免费视频| 久久精品国产亚洲av涩爱| 亚洲欧美清纯卡通| 国产一区二区三区av在线| 国产一区有黄有色的免费视频| 一级毛片久久久久久久久女| 日本-黄色视频高清免费观看| 国产无遮挡羞羞视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 中国美白少妇内射xxxbb| 一本久久精品| 蜜臀久久99精品久久宅男| 老熟女久久久| 舔av片在线| 一级毛片黄色毛片免费观看视频| 干丝袜人妻中文字幕| 日韩中文字幕视频在线看片 | 日本一二三区视频观看| 自拍偷自拍亚洲精品老妇| 国产色爽女视频免费观看| 天堂8中文在线网| 天堂中文最新版在线下载| 91精品国产九色| 九色成人免费人妻av| 国产精品一区二区三区四区免费观看| 亚洲国产色片| 亚洲美女视频黄频| 国产精品不卡视频一区二区| 午夜老司机福利剧场| 国产精品福利在线免费观看| 一级爰片在线观看| 久久6这里有精品| 中文字幕av成人在线电影| 精品亚洲乱码少妇综合久久| 伦理电影免费视频| 老司机影院毛片| 成年免费大片在线观看| 亚州av有码| 成人高潮视频无遮挡免费网站| 下体分泌物呈黄色| 日韩制服骚丝袜av| 色5月婷婷丁香| 精品久久久精品久久久| 精品人妻视频免费看| 免费黄频网站在线观看国产| 国产黄片视频在线免费观看| 22中文网久久字幕| 又爽又黄a免费视频| 天堂8中文在线网| 两个人的视频大全免费| 国产久久久一区二区三区| 成人特级av手机在线观看| 99视频精品全部免费 在线| 国产精品久久久久久精品古装| 亚洲精品色激情综合| 黄色一级大片看看| 亚洲精品成人av观看孕妇| 一本色道久久久久久精品综合| 亚洲最大成人中文| 国产精品一区www在线观看| 亚洲av成人精品一二三区| 中国美白少妇内射xxxbb| 日韩亚洲欧美综合| 又大又黄又爽视频免费| 久久毛片免费看一区二区三区| 国产在视频线精品| 黄片wwwwww| 欧美高清成人免费视频www| 国产成人a∨麻豆精品| 身体一侧抽搐| 性色av一级| 大码成人一级视频| 黄色配什么色好看| 亚洲国产色片| 久久国产乱子免费精品| 乱码一卡2卡4卡精品| 超碰av人人做人人爽久久| 成人毛片a级毛片在线播放| 99久久人妻综合| 欧美日韩视频高清一区二区三区二| 边亲边吃奶的免费视频| 精品人妻一区二区三区麻豆| 日本爱情动作片www.在线观看| 日韩一本色道免费dvd| av在线老鸭窝| 在线观看免费高清a一片| 国产精品成人在线|