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

    浮式振蕩立軸葉輪的三維效應(yīng)影響

    2019-06-19 02:09:10張錦嵐張萬(wàn)超周亞輝
    關(guān)鍵詞:展弦比立軸水輪機(jī)

    張錦嵐,王 凱,張萬(wàn)超,周亞輝

    (1. 武漢第二船舶設(shè)計(jì)研究所, 湖北 武漢 430064;2. 江蘇科技大學(xué) 船舶與海洋工程學(xué)院, 江蘇 鎮(zhèn)江 212003;3. 江蘇科技大學(xué) 能源與動(dòng)力學(xué)院, 江蘇 鎮(zhèn)江 212003)

    海洋可再生能源儲(chǔ)量豐富,總儲(chǔ)量可達(dá)7.66億千瓦,其開發(fā)利用逐漸受到世界各國(guó)的重視[1-3]。其中潮流能因其分布相對(duì)集中、能流密度高、有規(guī)律可循和相對(duì)穩(wěn)定的屬性特點(diǎn)[4-6],受到更多的重視。目前潮流能開發(fā)利用,已經(jīng)進(jìn)入復(fù)雜海況下的定常流研究階段。通過(guò)限定水輪機(jī)運(yùn)動(dòng)特性,探究復(fù)雜特殊工況下的水輪機(jī)水動(dòng)力特性及能量轉(zhuǎn)換是目前急需解決的任務(wù)。立軸葉輪強(qiáng)迫晃蕩耦合運(yùn)動(dòng)是解決上述問(wèn)題的關(guān)鍵,但相關(guān)研究成果還比較有限[7-9]。

    目前為止,利用計(jì)算流體力學(xué)(Computational Fluid Dynamics,CFD)方法模擬葉輪的運(yùn)行工況已經(jīng)比較成熟。2008年,Nabavi[10]應(yīng)用FLUENT軟件模擬了固定偏角水輪機(jī)的性能,并采用Spalart-Ailmaras湍流模型和滑移網(wǎng)格技術(shù)模擬水輪機(jī)二維及三維運(yùn)動(dòng),此外還分析了時(shí)間步、邊界距離、葉片處網(wǎng)格厚度對(duì)模擬結(jié)果的影響。2014年,Danao等[11]應(yīng)用FLUENT軟件分析了二維和三維垂直軸風(fēng)機(jī)在穩(wěn)態(tài)風(fēng)作用下的性能,湍流度、葉片表面節(jié)點(diǎn)數(shù)、邊界距離及最小時(shí)間步等因素被討論,算例采用滑移網(wǎng)格技術(shù),提出時(shí)間步應(yīng)取風(fēng)機(jī)轉(zhuǎn)過(guò)0.5°所需時(shí)間。2015年,Alaimo等[12]應(yīng)用FLUENT軟件模擬了二維及三維垂直軸水輪機(jī)在不同速比下的運(yùn)動(dòng)過(guò)程,著重分析了模擬時(shí)間步、模型網(wǎng)格數(shù)、葉片處網(wǎng)格尺寸等因素對(duì)模擬結(jié)果的影響。模型采用滑移網(wǎng)格技術(shù)且各算例均采用k-ε湍流模型。這些學(xué)者研究了葉輪三維和二維模擬下的葉輪水動(dòng)力性能影響。對(duì)于立軸葉輪在浪流作用下發(fā)生振蕩運(yùn)動(dòng)的三維和二維水動(dòng)力性能差別還少有學(xué)者進(jìn)行研究。

    近幾年,一些學(xué)者對(duì)浪流作用下的水平軸水輪機(jī)的水動(dòng)力性能做了一些研究[13-16],這些學(xué)者認(rèn)為浪流作用下水平軸水輪機(jī)的水動(dòng)力會(huì)發(fā)生波動(dòng),力和力矩的峰值會(huì)提高,但是平均值基本不變。哈爾濱工程大學(xué)張亮團(tuán)隊(duì)從潮流能水平軸水輪機(jī)實(shí)際運(yùn)行中也發(fā)現(xiàn)浪流作用對(duì)水輪機(jī)的水動(dòng)力性能有明顯影響。

    1 數(shù)值模擬

    1.1 理論基礎(chǔ)

    表1為立軸水輪機(jī)振蕩運(yùn)動(dòng)下的各參數(shù)定義。表1除去常規(guī)的葉輪參數(shù)定義外,特別定義了葉輪振蕩運(yùn)動(dòng)下的振蕩幅值、振蕩頻率和振蕩速度。

    表1 立軸水輪機(jī)葉輪參數(shù)

    通過(guò)表1,定義無(wú)量綱參數(shù)如下:

    推力系數(shù)

    CFX=FX/(0.5ρ·V2·D·H)

    (1)

    側(cè)向力系數(shù)

    CFY=FY/(0.5ρ·V2·D·H)

    (2)

    能量效率系數(shù)

    Cp=Q·ω/(0.5ρ·V3·D·H)

    (3)

    葉尖速比

    λ=R·ω/V

    (4)

    無(wú)因次振蕩速度

    (5)

    無(wú)因次振蕩加速度

    (6)

    為了簡(jiǎn)化描述,使用變量Φ來(lái)表示坐標(biāo)軸方向,其中Φ為X時(shí)表示X軸方向,Φ為Y時(shí)表示Y軸方向。葉輪的受力可以分解為三類,即均勻水動(dòng)力、阻尼力和附加質(zhì)量力,則葉輪受力表示為

    (7)

    對(duì)阻尼力和附加質(zhì)量力進(jìn)行分解,表示為

    (8)

    (9)

    其中,ψ表示葉輪位置角和各項(xiàng)展開項(xiàng)的錯(cuò)位角。將式(8)和式(9)代入式(7),則葉輪振蕩運(yùn)動(dòng)下的受力可表示為

    (10)

    利用ANSYS-CFX軟件對(duì)立軸葉輪振蕩運(yùn)動(dòng)進(jìn)行模擬仿真,并將得到的葉輪振蕩運(yùn)動(dòng)下的推力和側(cè)向力時(shí)歷曲線利用最小二乘法擬合,就可以得到式(10)中的各項(xiàng)系數(shù),總結(jié)對(duì)比二維模擬和三維模擬各參數(shù)對(duì)立軸葉輪振蕩運(yùn)動(dòng)的影響規(guī)律。

    1.2 立軸水輪機(jī)模型

    選用雙葉片葉輪模型,其葉片采用NACA0018對(duì)稱翼型,葉輪直徑為0.8 m,弦長(zhǎng)為0.12 m,來(lái)流速度為1.0 m/s。

    圖1和圖2分別為立軸振蕩葉輪數(shù)值計(jì)算的二維和三維網(wǎng)格模型。本文主要計(jì)算葉輪旋轉(zhuǎn)運(yùn)動(dòng)和振蕩運(yùn)動(dòng)相結(jié)合工況下的二維和三維結(jié)果,其中振蕩運(yùn)動(dòng)包括橫蕩運(yùn)動(dòng)和縱蕩運(yùn)動(dòng)。兩幅圖中分別給出了二維和三維工況下的各運(yùn)動(dòng)形式計(jì)算域,主要包括旋轉(zhuǎn)域、振蕩域和隨振蕩運(yùn)動(dòng)網(wǎng)格發(fā)生微幅變形的變形域。

    圖3為葉輪振蕩運(yùn)動(dòng)模擬設(shè)置的邊界條件。邊界條件的設(shè)置和常規(guī)旋轉(zhuǎn)葉輪的邊界條件設(shè)置類似,所不同的是本文中將旋轉(zhuǎn)域以外的流場(chǎng)域分為振蕩域和振蕩變形域,其中振蕩域在計(jì)算中發(fā)生橫蕩或者縱蕩運(yùn)動(dòng),運(yùn)動(dòng)過(guò)程中網(wǎng)格只存在滑移運(yùn)動(dòng),網(wǎng)格不發(fā)生變形;振蕩變形域在計(jì)算中網(wǎng)格發(fā)生微幅拉伸和壓縮,但是對(duì)網(wǎng)格質(zhì)量幾乎沒(méi)有影響。

    (a) 葉片(a) Blade (b) 旋轉(zhuǎn)域(b) Rotation domain

    (c) 橫蕩域(c) Sway region (d) 橫蕩變形域(d) Sway deformation region

    (e) 縱蕩域(e) Surge region (f) 縱蕩變形域(f) Surge deformation region 圖1 二維網(wǎng)格模型Fig.1 2D grid model

    (a) 葉片(a) Blade (b) 旋轉(zhuǎn)域(b) Rotation domain

    (c) 橫蕩域(c) Sway region (d) 橫蕩變形域(d) Sway deformation region

    (e) 縱蕩域(e) Surge region (f) 縱蕩變形域(f) Surge deformation region圖2 三維網(wǎng)格模型Fig.2 3D grid model

    (a) 二維橫蕩(a) Two-dimensional sway motion

    (b) 二維縱蕩(b) Two-dimensional surge motion

    (c) 三維橫蕩(c) Three-dimensional sway motion

    (d) 三維縱蕩(d) Three-dimensional surge motion圖3 邊界條件Fig.3 Boundary conditions

    2 結(jié)果分析

    2.1 CFD方法的可行性證明

    在計(jì)算模擬前,首先驗(yàn)證立軸葉輪仿真計(jì)算的有效性,分別從二維模擬和三維模擬方面對(duì)計(jì)算值和試驗(yàn)值進(jìn)行對(duì)比分析。

    二維計(jì)算驗(yàn)證:采用本文的計(jì)算方法對(duì)美國(guó)Strikland教授做的關(guān)于立軸葉輪的二維試驗(yàn)[17]進(jìn)行模擬,并將計(jì)算得到的單個(gè)葉片切向力、法向力與試驗(yàn)結(jié)果對(duì)比。圖4為試驗(yàn)值和本文計(jì)算結(jié)果的對(duì)比,從兩幅圖中可以看出,計(jì)算結(jié)果和試驗(yàn)值相近。教授也對(duì)該葉輪進(jìn)行了數(shù)值模擬計(jì)算,本文結(jié)果和其計(jì)算結(jié)果類似,在90°~270°的上游盤面,計(jì)算值和試驗(yàn)值誤差較小,而270°~450°的下游盤面誤差較大。Strikland教授分析可能是:一方面試驗(yàn)時(shí)葉片安裝有問(wèn)題,造成下游盤面試驗(yàn)值存在凸起現(xiàn)象;另一方面下游盤面水流受到上游盤面影響,水流復(fù)雜,導(dǎo)致試驗(yàn)值和計(jì)算值存在一定差異。

    為了驗(yàn)證本文計(jì)算方法三維模擬的有效性,對(duì)比了哈爾濱工程大學(xué)“HEU”型兩葉片立軸水輪機(jī)的水動(dòng)力性能,結(jié)果如圖5所示。圖5(a)和圖5(b)分別為葉輪效率和推力的試驗(yàn)值和三維計(jì)算值對(duì)比。不論是效率還是推力,三維仿真計(jì)算結(jié)果都略大于試驗(yàn)值,這是因?yàn)榉抡嬗?jì)算中未考慮葉輪運(yùn)行過(guò)程中機(jī)械摩擦、壁效應(yīng)的影響,這些影響都會(huì)造成葉輪效率和推力的降低。綜合來(lái)看,立軸葉輪三維模擬計(jì)算精度較高,能夠模擬立軸葉輪真實(shí)運(yùn)行情況。

    (a) 切向力(a) Tangential force

    (b) 法向力(b) Normal force圖4 二維計(jì)算對(duì)比圖Fig.4 Comparisons between calculated and experimental values

    (a) 效率對(duì)比(a) Comparison of efficiency

    (b) 推力系數(shù)對(duì)比(b) Comparison of thrust coefficient圖5 三維計(jì)算對(duì)比Fig.5 Results comparison between the experimental and 3D numerical simulation

    2.2 立軸葉輪振蕩運(yùn)動(dòng)效率

    當(dāng)展弦比較小時(shí),葉輪振蕩運(yùn)動(dòng)下二維計(jì)算和三維計(jì)算結(jié)果相差較大。為了尋找不同展弦比下二維結(jié)果和三維結(jié)果的變化規(guī)律,首先計(jì)算了不同展弦比的立軸葉輪振蕩運(yùn)動(dòng),并將效率作為對(duì)比對(duì)象和二維結(jié)果進(jìn)行比較。如圖6所示,兩圖為橫蕩運(yùn)動(dòng)和縱蕩運(yùn)動(dòng)下不同展弦比三維葉輪效率和二維效率對(duì)比的結(jié)果。隨著展弦比增大,三維計(jì)算結(jié)果逐漸增大,當(dāng)展弦比達(dá)到10左右時(shí),三維計(jì)算結(jié)果和二維計(jì)算結(jié)果相近。葉輪三維計(jì)算結(jié)果和二維計(jì)算結(jié)果存在差異的最主要原因是三維計(jì)算中存在三維效應(yīng),葉片展長(zhǎng)兩側(cè)存在機(jī)翼繞流,葉梢處高壓區(qū)流體會(huì)流向低壓區(qū),造成葉梢處尾渦,帶走一部分能量,因此葉輪三維效率計(jì)算結(jié)果小于二維計(jì)算結(jié)果。

    (a) 橫蕩運(yùn)動(dòng)(a) Sway motion

    (b) 縱蕩運(yùn)動(dòng)(b) Surge motion圖6 不同展弦比的水輪機(jī)效率曲線Fig.6 Comparison of power output curves of vertical axis turbine with different span-chord ratios

    2.3 立軸葉輪振蕩運(yùn)動(dòng)推力和側(cè)向力

    圖7所示為葉輪橫蕩和縱蕩運(yùn)動(dòng)下的不同展弦比的推力系數(shù)和側(cè)向力系數(shù)的對(duì)比。其中圖7(a)和圖7(c)為推力系數(shù)對(duì)比,圖7(b)和圖7(d)為側(cè)向力系數(shù)對(duì)比。四幅圖中分別對(duì)比了展長(zhǎng)為0.2 m、0.6 m和1.0 m的三維計(jì)算結(jié)果和二維計(jì)算結(jié)果。從四幅圖中可以看出,隨著展長(zhǎng)增大,相對(duì)應(yīng)的推力和側(cè)向力系數(shù)逐漸接近二維計(jì)算結(jié)果。推力和側(cè)向力的波峰和波谷都隨著展長(zhǎng)增大而增大。這是因?yàn)檎归L(zhǎng)越長(zhǎng),葉輪的三維效應(yīng)影響越小,葉梢處尾渦帶走的能量越少,計(jì)算結(jié)果越接近二維計(jì)算值。

    (a) 橫蕩推力(a) Thrust of sway

    (b) 橫蕩側(cè)向力(b) Lateral force of sway

    (c) 縱蕩推力(c) Thrust of surge

    (d) 縱蕩側(cè)向力(d) Lateral force of surge圖7 不同展長(zhǎng)的葉輪振蕩狀態(tài)下的推力和側(cè)向力系數(shù)對(duì)比Fig.7 Comparison of thrust and lateral force coefficients under swaying motion and surging motion of different elongation

    圖8為三個(gè)不同展長(zhǎng)的葉輪流場(chǎng)圖,從三幅圖中可以看出上游葉片葉梢兩端產(chǎn)生流體分離,帶走了一部分能量。圖8(a)所示展長(zhǎng)0.2 m的葉輪流體分離幾乎布滿整個(gè)葉片,明顯大于圖8(b)和圖8(c)。結(jié)合圖9所示的不同展長(zhǎng)的葉輪渦量場(chǎng),在同一渦量指標(biāo)下,展長(zhǎng)0.2 m的葉輪兩端渦量強(qiáng)度明顯大于0.4 m和0.6 m展長(zhǎng)的渦量強(qiáng)度,其帶走的能量也越多,因此導(dǎo)致葉輪水動(dòng)力性能降低。

    從圖8和圖9可以看出,展弦比越大,葉輪由于三維效應(yīng)在葉片兩端產(chǎn)生的尾渦越小,帶走的能量也越小,因此葉輪效率、推力等水動(dòng)力系數(shù)越高。結(jié)合2.2節(jié),當(dāng)葉輪的展弦比達(dá)到10左右時(shí),葉輪尾渦帶走的能量可以忽略不計(jì),水動(dòng)力性能也與二維葉輪的水動(dòng)力性能相近。

    (a) H=0.2 m (b) H=0.4 m (c) H=0.6 m圖8 不同展長(zhǎng)的葉輪流場(chǎng)流線圖Fig.8 Streamlines of the flow field under different elongations

    (a) H=0.2 m (b) H=0.4 m (c) H=0.6 m圖9 不同展長(zhǎng)的葉輪渦量場(chǎng)Fig.9 Vorticity fields under different elongations

    3 水動(dòng)力系數(shù)

    根據(jù)式(10)將三維模擬計(jì)算得到的不同展長(zhǎng)的立軸葉輪推力和側(cè)向力系數(shù)進(jìn)行擬合分析,得到各受力系數(shù)的均勻水動(dòng)力、阻尼力和附加質(zhì)量力系數(shù)。如表2~5所示,對(duì)比了0.2 m、0.4 m、0.6 m、0.8 m四個(gè)展長(zhǎng)和2D模擬橫蕩與縱蕩運(yùn)動(dòng)下的各項(xiàng)系數(shù)。

    表2 橫蕩運(yùn)動(dòng)CFX展開系數(shù)

    表3 縱蕩運(yùn)動(dòng)CFX展開系數(shù)

    表4 橫蕩運(yùn)動(dòng)CFY展開系數(shù)

    表5 縱蕩運(yùn)動(dòng)CFY展開系數(shù)

    表4和表5為振蕩運(yùn)動(dòng)下側(cè)向力的展開系數(shù),規(guī)律和表2、表3類似。隨著展長(zhǎng)增大,振蕩葉輪的側(cè)向力展開項(xiàng)中附加質(zhì)量力所占比例較小,且另外兩項(xiàng)系數(shù)絕對(duì)值隨著展長(zhǎng)增大而增大,逐漸接近二維結(jié)果。

    4 結(jié)論

    利用流體計(jì)算軟件ANSYS-CFX計(jì)算了不同展弦比的葉輪三維振蕩運(yùn)動(dòng),包括橫蕩和縱蕩運(yùn)動(dòng),并將計(jì)算結(jié)果和二維模擬結(jié)果做了對(duì)比,結(jié)論如下:

    1)立軸葉輪三維模擬結(jié)果中水動(dòng)力系數(shù)和二維結(jié)果存在一定差異,當(dāng)展弦比接近10時(shí),三維結(jié)果和二維結(jié)果相差不大。

    2)隨著展弦比增大,三維振蕩葉輪的輸出效率、推力和側(cè)向力系數(shù)逐漸增大。

    3)三維振蕩葉輪的水動(dòng)力性能計(jì)算結(jié)果低于二維計(jì)算結(jié)果,最主要原因是三維葉片兩端流體從高壓區(qū)流向低壓區(qū),造成能量損失。

    4)三維振蕩葉輪的推力和側(cè)向力各項(xiàng)中,附加質(zhì)量力和均勻水動(dòng)力、阻尼力相比為小量。

    猜你喜歡
    展弦比立軸水輪機(jī)
    水輪機(jī)過(guò)流面非金屬材料的修復(fù)及防護(hù)
    江山如畫
    ——遼寧省博物館藏中國(guó)古代立軸山水畫展
    中華書畫家(2022年4期)2022-04-21 09:31:06
    八大山人的大尺幅行草立軸書法探析
    大中型水斗式水輪機(jī)的關(guān)鍵技術(shù)
    不同展弦比下扭轉(zhuǎn)葉片振動(dòng)特性分析
    水輪機(jī)虛擬仿真動(dòng)畫制作的研究
    大展弦比復(fù)合材料機(jī)翼結(jié)構(gòu)細(xì)節(jié)抗疲勞優(yōu)化
    名家藝苑
    雅昌在線拍 最新拍品
    優(yōu)雅(2018年2期)2018-01-30 12:26:26
    氫動(dòng)力無(wú)人機(jī)大展弦比機(jī)翼靜氣彈特性分析
    777久久人妻少妇嫩草av网站| 欧美性猛交黑人性爽| 午夜福利免费观看在线| 岛国在线免费视频观看| 欧美另类亚洲清纯唯美| 日韩高清综合在线| 欧美久久黑人一区二区| 精品国产超薄肉色丝袜足j| 国产一区二区三区视频了| 99久久精品热视频| 九九热线精品视视频播放| 国产欧美日韩一区二区精品| 久久草成人影院| 国产精品永久免费网站| 三级男女做爰猛烈吃奶摸视频| 少妇粗大呻吟视频| avwww免费| 亚洲av第一区精品v没综合| 欧美色欧美亚洲另类二区| 国产av麻豆久久久久久久| 国产精品久久久久久久电影 | 色综合婷婷激情| 身体一侧抽搐| 亚洲精品国产精品久久久不卡| 少妇裸体淫交视频免费看高清 | 亚洲成人免费电影在线观看| 丰满的人妻完整版| 国产亚洲av嫩草精品影院| 高清在线国产一区| 国产99白浆流出| 美女大奶头视频| 听说在线观看完整版免费高清| 日韩有码中文字幕| 黑人操中国人逼视频| 精品电影一区二区在线| 色综合站精品国产| 女人被狂操c到高潮| 韩国av一区二区三区四区| 嫩草影视91久久| 亚洲国产欧美人成| 亚洲成a人片在线一区二区| 精品高清国产在线一区| 国产三级在线视频| 听说在线观看完整版免费高清| 亚洲人成伊人成综合网2020| 亚洲专区字幕在线| 日韩精品免费视频一区二区三区| 亚洲一码二码三码区别大吗| 国产精品永久免费网站| 女同久久另类99精品国产91| 丰满人妻一区二区三区视频av | 亚洲男人天堂网一区| 国产成人影院久久av| 级片在线观看| 精品不卡国产一区二区三区| 十八禁网站免费在线| 午夜福利视频1000在线观看| 啦啦啦观看免费观看视频高清| 欧美日韩黄片免| 国产精品一区二区三区四区久久| 可以在线观看的亚洲视频| 一级毛片精品| 亚洲乱码一区二区免费版| 久久人妻福利社区极品人妻图片| 成年女人毛片免费观看观看9| 久久精品亚洲精品国产色婷小说| 日韩精品免费视频一区二区三区| 国模一区二区三区四区视频 | 亚洲国产高清在线一区二区三| 18禁国产床啪视频网站| 99精品久久久久人妻精品| 脱女人内裤的视频| 一区二区三区高清视频在线| 久久久久久国产a免费观看| 女警被强在线播放| 欧美性长视频在线观看| 久久天堂一区二区三区四区| 在线观看日韩欧美| 亚洲av中文字字幕乱码综合| 久久人人精品亚洲av| 99国产综合亚洲精品| 欧美黄色淫秽网站| av免费在线观看网站| 国产精品久久久人人做人人爽| 18禁裸乳无遮挡免费网站照片| 久久午夜综合久久蜜桃| 999久久久国产精品视频| 床上黄色一级片| 国产麻豆成人av免费视频| 中文字幕高清在线视频| 国产精品久久久久久亚洲av鲁大| 叶爱在线成人免费视频播放| 久久久久国产一级毛片高清牌| 少妇熟女aⅴ在线视频| 啦啦啦免费观看视频1| 欧美在线一区亚洲| 日韩欧美国产在线观看| 天天躁夜夜躁狠狠躁躁| 国产精品av久久久久免费| 首页视频小说图片口味搜索| 国产乱人伦免费视频| 一进一出好大好爽视频| 午夜精品久久久久久毛片777| 真人做人爱边吃奶动态| 国产av在哪里看| 免费观看人在逋| 午夜福利成人在线免费观看| 国产一级毛片七仙女欲春2| 久久久水蜜桃国产精品网| 欧美日韩福利视频一区二区| 日本熟妇午夜| 99精品欧美一区二区三区四区| 好男人在线观看高清免费视频| 麻豆一二三区av精品| 亚洲人成网站在线播放欧美日韩| 欧美 亚洲 国产 日韩一| 18禁美女被吸乳视频| 欧美日韩乱码在线| 高清毛片免费观看视频网站| 亚洲男人的天堂狠狠| 欧美日韩精品网址| 小说图片视频综合网站| 草草在线视频免费看| 岛国在线免费视频观看| 色综合欧美亚洲国产小说| 国产免费男女视频| 午夜影院日韩av| 丰满人妻一区二区三区视频av | 久久亚洲真实| 精品久久蜜臀av无| 法律面前人人平等表现在哪些方面| 欧美绝顶高潮抽搐喷水| 国产精品亚洲一级av第二区| 国产精品香港三级国产av潘金莲| 18禁美女被吸乳视频| 欧美一级毛片孕妇| 国产精品亚洲美女久久久| 亚洲av五月六月丁香网| ponron亚洲| 久久天堂一区二区三区四区| 在线观看日韩欧美| 在线观看舔阴道视频| 巨乳人妻的诱惑在线观看| 色精品久久人妻99蜜桃| 精品久久久久久久人妻蜜臀av| 久久人人精品亚洲av| av超薄肉色丝袜交足视频| 国产精品一区二区免费欧美| 亚洲自拍偷在线| 日本 欧美在线| 国内精品久久久久精免费| 国产精品影院久久| 国产精品久久视频播放| 老司机福利观看| 蜜桃久久精品国产亚洲av| 18禁美女被吸乳视频| 国产欧美日韩一区二区精品| 舔av片在线| 嫩草影院精品99| 青草久久国产| 美女扒开内裤让男人捅视频| 91在线观看av| 男人舔奶头视频| 床上黄色一级片| 免费看日本二区| 久久国产精品影院| 在线观看美女被高潮喷水网站 | 国产视频一区二区在线看| 久久精品影院6| 又爽又黄无遮挡网站| 一级黄色大片毛片| 琪琪午夜伦伦电影理论片6080| 中出人妻视频一区二区| 人妻夜夜爽99麻豆av| 欧美性猛交╳xxx乱大交人| 欧美日韩亚洲国产一区二区在线观看| 妹子高潮喷水视频| 亚洲人成77777在线视频| 精品午夜福利视频在线观看一区| 午夜福利高清视频| 日本免费a在线| 欧美国产日韩亚洲一区| av中文乱码字幕在线| 国产精品久久久久久人妻精品电影| 国产精品一区二区免费欧美| 亚洲专区中文字幕在线| av有码第一页| 亚洲第一电影网av| 亚洲专区国产一区二区| 草草在线视频免费看| 99精品久久久久人妻精品| 亚洲最大成人中文| 久久精品综合一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 神马国产精品三级电影在线观看 | 老汉色∧v一级毛片| 少妇裸体淫交视频免费看高清 | 久久亚洲精品不卡| x7x7x7水蜜桃| 亚洲中文字幕日韩| 亚洲一码二码三码区别大吗| 日日爽夜夜爽网站| 美女午夜性视频免费| 制服诱惑二区| 亚洲欧美日韩东京热| 美女高潮喷水抽搐中文字幕| 最新在线观看一区二区三区| 一区二区三区激情视频| 亚洲国产精品成人综合色| 免费观看精品视频网站| 国产不卡一卡二| 午夜精品久久久久久毛片777| 国内精品一区二区在线观看| 亚洲av电影在线进入| 久久伊人香网站| 伊人久久大香线蕉亚洲五| 男女视频在线观看网站免费 | 狠狠狠狠99中文字幕| 欧美一级a爱片免费观看看 | 国产1区2区3区精品| 美女午夜性视频免费| 99精品久久久久人妻精品| 亚洲色图av天堂| 欧美性猛交黑人性爽| 国产高清videossex| 欧美又色又爽又黄视频| 亚洲色图 男人天堂 中文字幕| 亚洲专区字幕在线| 淫秽高清视频在线观看| 人妻丰满熟妇av一区二区三区| 日韩 欧美 亚洲 中文字幕| 麻豆国产av国片精品| 亚洲欧美日韩东京热| 麻豆国产97在线/欧美 | 久久精品国产亚洲av香蕉五月| 精品熟女少妇八av免费久了| 变态另类成人亚洲欧美熟女| 国产黄色小视频在线观看| 成人一区二区视频在线观看| 在线永久观看黄色视频| 精品欧美一区二区三区在线| 久久欧美精品欧美久久欧美| 欧美人与性动交α欧美精品济南到| 欧美激情久久久久久爽电影| 免费电影在线观看免费观看| 欧美乱妇无乱码| 麻豆一二三区av精品| 色在线成人网| 色播亚洲综合网| 欧美大码av| 一进一出抽搐gif免费好疼| 欧美黑人巨大hd| 亚洲中文日韩欧美视频| 免费高清视频大片| 国内久久婷婷六月综合欲色啪| 国产av在哪里看| 久久午夜亚洲精品久久| 在线观看免费视频日本深夜| 在线十欧美十亚洲十日本专区| 一二三四社区在线视频社区8| 精品午夜福利视频在线观看一区| 村上凉子中文字幕在线| 热99re8久久精品国产| 国产亚洲av高清不卡| 一夜夜www| 一级毛片精品| 国产av不卡久久| 日本三级黄在线观看| 两个人的视频大全免费| 男女视频在线观看网站免费 | av超薄肉色丝袜交足视频| 搡老熟女国产l中国老女人| 五月玫瑰六月丁香| 午夜久久久久精精品| 无遮挡黄片免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久免费高清国产稀缺| 免费搜索国产男女视频| 国产精品亚洲美女久久久| 女生性感内裤真人,穿戴方法视频| 国产人伦9x9x在线观看| 国产精品久久久人人做人人爽| 成人国语在线视频| 亚洲乱码一区二区免费版| 精品熟女少妇八av免费久了| 久久婷婷成人综合色麻豆| 亚洲国产欧美一区二区综合| 欧美日韩精品网址| 亚洲av第一区精品v没综合| 国产成人精品久久二区二区免费| 国产午夜精品久久久久久| 国产高清videossex| 美女午夜性视频免费| 亚洲成a人片在线一区二区| 日本免费一区二区三区高清不卡| 毛片女人毛片| 妹子高潮喷水视频| 欧美成狂野欧美在线观看| 久热爱精品视频在线9| 日日干狠狠操夜夜爽| 9191精品国产免费久久| 一个人免费在线观看的高清视频| 成人精品一区二区免费| 免费高清视频大片| 脱女人内裤的视频| 长腿黑丝高跟| 国内精品久久久久精免费| 国产91精品成人一区二区三区| 精品日产1卡2卡| av中文乱码字幕在线| 亚洲成人久久爱视频| 欧美三级亚洲精品| 毛片女人毛片| 身体一侧抽搐| 久久香蕉激情| 日韩免费av在线播放| 成年人黄色毛片网站| 精品不卡国产一区二区三区| 俺也久久电影网| 国产精品1区2区在线观看.| 制服丝袜大香蕉在线| 亚洲av片天天在线观看| 国产精品自产拍在线观看55亚洲| 色综合婷婷激情| 女生性感内裤真人,穿戴方法视频| 50天的宝宝边吃奶边哭怎么回事| 国产高清视频在线观看网站| 啦啦啦观看免费观看视频高清| 国产精品久久久久久精品电影| 免费看日本二区| 2021天堂中文幕一二区在线观| 亚洲成人久久爱视频| 成人av一区二区三区在线看| 日韩有码中文字幕| 嫩草影视91久久| 欧美日韩一级在线毛片| 精品电影一区二区在线| 美女 人体艺术 gogo| 少妇粗大呻吟视频| 成人av在线播放网站| 身体一侧抽搐| 久99久视频精品免费| 悠悠久久av| 亚洲色图 男人天堂 中文字幕| 成人特级黄色片久久久久久久| 最近最新免费中文字幕在线| 亚洲欧美日韩东京热| 亚洲九九香蕉| 国产午夜精品论理片| 亚洲性夜色夜夜综合| 欧美最黄视频在线播放免费| 亚洲国产欧美人成| 黑人操中国人逼视频| 中文字幕熟女人妻在线| 国产亚洲精品av在线| 亚洲狠狠婷婷综合久久图片| 亚洲国产日韩欧美精品在线观看 | 18禁国产床啪视频网站| 精品第一国产精品| 国产精品一区二区三区四区久久| 日日夜夜操网爽| 变态另类丝袜制服| 免费一级毛片在线播放高清视频| 久久香蕉国产精品| 午夜日韩欧美国产| 俄罗斯特黄特色一大片| 日韩国内少妇激情av| 我要搜黄色片| 欧美性长视频在线观看| 99久久无色码亚洲精品果冻| 国产精品国产高清国产av| 黄频高清免费视频| 婷婷亚洲欧美| 欧美成人性av电影在线观看| 欧美绝顶高潮抽搐喷水| 一个人观看的视频www高清免费观看 | 亚洲性夜色夜夜综合| 亚洲国产欧美网| 亚洲人成77777在线视频| 91在线观看av| 国内毛片毛片毛片毛片毛片| 不卡av一区二区三区| 一个人免费在线观看电影 | 美女免费视频网站| 99久久综合精品五月天人人| 在线观看一区二区三区| 一本精品99久久精品77| 亚洲av美国av| 亚洲 欧美一区二区三区| av有码第一页| 欧美乱码精品一区二区三区| 亚洲欧美一区二区三区黑人| 手机成人av网站| 老鸭窝网址在线观看| 又紧又爽又黄一区二区| 午夜精品在线福利| 国产一区二区三区视频了| 一本一本综合久久| www.自偷自拍.com| 国产精品亚洲av一区麻豆| 久久久国产欧美日韩av| 精品久久久久久成人av| 国产精品一及| 可以在线观看毛片的网站| 国产一区二区激情短视频| 一区二区三区国产精品乱码| 午夜亚洲福利在线播放| av免费在线观看网站| www.www免费av| 动漫黄色视频在线观看| 国产成人av教育| 国产男靠女视频免费网站| 性欧美人与动物交配| 啦啦啦免费观看视频1| 日韩欧美在线二视频| 欧美午夜高清在线| 成人欧美大片| 少妇裸体淫交视频免费看高清 | 亚洲中文av在线| 日韩欧美在线乱码| 国产高清有码在线观看视频 | 少妇裸体淫交视频免费看高清 | 日韩欧美免费精品| 白带黄色成豆腐渣| 两性夫妻黄色片| 男人舔女人的私密视频| 欧美3d第一页| 国产免费av片在线观看野外av| 欧美久久黑人一区二区| 久久香蕉国产精品| 熟妇人妻久久中文字幕3abv| 亚洲黑人精品在线| 国产99白浆流出| 日本一区二区免费在线视频| 亚洲国产精品sss在线观看| 欧美极品一区二区三区四区| 在线免费观看的www视频| 国产激情欧美一区二区| 91av网站免费观看| 亚洲精品一卡2卡三卡4卡5卡| 日韩大码丰满熟妇| 久久久久久九九精品二区国产 | 淫妇啪啪啪对白视频| 人人妻,人人澡人人爽秒播| 最新在线观看一区二区三区| 两个人免费观看高清视频| 久久午夜综合久久蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 无限看片的www在线观看| 一个人观看的视频www高清免费观看 | 中文资源天堂在线| 欧美黑人巨大hd| 动漫黄色视频在线观看| 国产麻豆成人av免费视频| 国产野战对白在线观看| 欧美中文综合在线视频| 亚洲国产欧洲综合997久久,| 日韩欧美在线乱码| 亚洲国产精品成人综合色| 久久香蕉精品热| 久久久国产成人免费| 免费看a级黄色片| 亚洲男人的天堂狠狠| 亚洲av美国av| 国产亚洲精品av在线| 欧美在线一区亚洲| 久久婷婷成人综合色麻豆| 国产日本99.免费观看| 日韩精品免费视频一区二区三区| 国产成人av激情在线播放| 一个人免费在线观看的高清视频| 亚洲午夜精品一区,二区,三区| 男人的好看免费观看在线视频 | 精品久久久久久成人av| 国产欧美日韩一区二区三| 中文字幕久久专区| 黄频高清免费视频| 一区二区三区激情视频| xxxwww97欧美| 欧美黑人精品巨大| 国产av麻豆久久久久久久| 听说在线观看完整版免费高清| 亚洲色图av天堂| 好男人在线观看高清免费视频| 国产男靠女视频免费网站| √禁漫天堂资源中文www| www.999成人在线观看| 亚洲国产欧美网| 亚洲国产精品成人综合色| 999精品在线视频| e午夜精品久久久久久久| 男女做爰动态图高潮gif福利片| 看免费av毛片| 日本a在线网址| 久久精品91蜜桃| 精品一区二区三区四区五区乱码| 亚洲精品久久成人aⅴ小说| 国产精品一区二区免费欧美| 叶爱在线成人免费视频播放| 欧美成人午夜精品| 三级毛片av免费| 一区福利在线观看| 日本一本二区三区精品| 亚洲性夜色夜夜综合| 国产精品1区2区在线观看.| 欧美日韩福利视频一区二区| 午夜福利欧美成人| 啦啦啦观看免费观看视频高清| 在线观看美女被高潮喷水网站 | 亚洲国产欧美一区二区综合| 日本免费一区二区三区高清不卡| 欧美日韩福利视频一区二区| 国内少妇人妻偷人精品xxx网站 | 亚洲国产精品久久男人天堂| 黄色 视频免费看| 午夜精品在线福利| 日本熟妇午夜| 亚洲七黄色美女视频| 婷婷精品国产亚洲av| 午夜福利在线观看吧| 亚洲中文字幕日韩| 好男人在线观看高清免费视频| 亚洲专区中文字幕在线| av片东京热男人的天堂| 国产在线观看jvid| 欧美一级毛片孕妇| 免费av毛片视频| 久久久久久久久免费视频了| 在线永久观看黄色视频| 亚洲人成伊人成综合网2020| 色在线成人网| 成人国语在线视频| 美女午夜性视频免费| 亚洲avbb在线观看| 国产又色又爽无遮挡免费看| 国产成人啪精品午夜网站| 久久中文看片网| 久久久久精品国产欧美久久久| 1024香蕉在线观看| 黄色a级毛片大全视频| 国产午夜精品久久久久久| 亚洲av片天天在线观看| 免费看美女性在线毛片视频| 欧美性长视频在线观看| 色播亚洲综合网| 国产精品一区二区精品视频观看| 一级片免费观看大全| 全区人妻精品视频| 国产真实乱freesex| 成人国产一区最新在线观看| 性欧美人与动物交配| 国产高清激情床上av| 精品午夜福利视频在线观看一区| 国产蜜桃级精品一区二区三区| 久久久久精品国产欧美久久久| 黄片小视频在线播放| 成人欧美大片| 黄色视频,在线免费观看| 欧美乱色亚洲激情| 欧美在线一区亚洲| 又紧又爽又黄一区二区| 一进一出抽搐动态| 午夜福利18| 欧美中文日本在线观看视频| 一级a爱片免费观看的视频| 久久草成人影院| 五月伊人婷婷丁香| 丰满人妻熟妇乱又伦精品不卡| 天堂√8在线中文| 看片在线看免费视频| 欧美性猛交╳xxx乱大交人| 日韩国内少妇激情av| 1024视频免费在线观看| 久久婷婷成人综合色麻豆| 久久这里只有精品19| 午夜福利18| www.精华液| 国产成年人精品一区二区| 看免费av毛片| 听说在线观看完整版免费高清| 老熟妇乱子伦视频在线观看| 久久久久久久精品吃奶| 日韩中文字幕欧美一区二区| 99在线视频只有这里精品首页| 久久国产乱子伦精品免费另类| 久久久久精品国产欧美久久久| 久久精品综合一区二区三区| 日韩精品中文字幕看吧| 男插女下体视频免费在线播放| ponron亚洲| 成人国语在线视频| 他把我摸到了高潮在线观看| 真人做人爱边吃奶动态| 欧美日韩乱码在线| 久热爱精品视频在线9| 欧美成狂野欧美在线观看| 99久久久亚洲精品蜜臀av| 少妇的丰满在线观看| 日韩 欧美 亚洲 中文字幕| 正在播放国产对白刺激| 亚洲精品粉嫩美女一区| 韩国av一区二区三区四区| 中亚洲国语对白在线视频| 久久久久久免费高清国产稀缺| 热99re8久久精品国产| 两个人免费观看高清视频| 50天的宝宝边吃奶边哭怎么回事| 久久中文字幕人妻熟女| 亚洲精品一卡2卡三卡4卡5卡| 成年人黄色毛片网站| 欧美一区二区精品小视频在线| 精品久久蜜臀av无| 久久99热这里只有精品18| 国产成人精品久久二区二区免费| 欧美成人性av电影在线观看|