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

    相對流速對非球形鋁滴旋轉(zhuǎn)受力的影響①

    2015-04-24 08:35:52劉叢林康文武
    固體火箭技術(shù) 2015年2期
    關(guān)鍵詞:來流雷諾數(shù)升力

    劉叢林,郜 冶,賀 征,康文武

    (1.哈爾濱工程大學(xué) 航天與建筑工程學(xué)院,哈爾濱 150001;2.北京航空航天大學(xué) 能源與動力工程學(xué)院,北京 100191)

    ?

    相對流速對非球形鋁滴旋轉(zhuǎn)受力的影響①

    劉叢林1,郜 冶1,賀 征1,康文武2

    (1.哈爾濱工程大學(xué) 航天與建筑工程學(xué)院,哈爾濱 150001;2.北京航空航天大學(xué) 能源與動力工程學(xué)院,北京 100191)

    針對固體火箭發(fā)動機(jī)中非球形鋁滴的旋轉(zhuǎn)受力問題進(jìn)行數(shù)值模擬,詳細(xì)分析流速變化對顆粒的受力影響。計算結(jié)果表明,非球顆粒在流場中的升力系數(shù)與無量綱轉(zhuǎn)速之間已經(jīng)不再滿足經(jīng)典的理論關(guān)系式。雷諾數(shù)Re<1時,周向升力系數(shù)存在著交替出現(xiàn)的2組極大值與極小值;升力系數(shù)均為正值,且周向變化較為平緩。隨來流速度增加,阻力系數(shù)變得均衡,系數(shù)值向0值逼近;升力系數(shù)也急速降低。當(dāng)相對流速由0.5 m/s提高到5 m/s時,顆粒所受阻力銳降89%。顆粒的阻力系數(shù)與雷諾數(shù)之間滿足新的關(guān)系式。

    非球形鋁滴;相對流速;阻力系數(shù);升力系數(shù)

    0 引言

    在含顆粒的多相流計算中,一般視顆粒為離散相,顆粒在流場中的運動軌跡及空間分布取決于作用于其上的力和力矩,對流場中的離散相進(jìn)行受力分析是獲得準(zhǔn)確計算結(jié)果的基礎(chǔ)。尤其固體火箭發(fā)動機(jī)中,金屬鋁滴在流場中相變?nèi)紵龝r,會因部分氧化鋁凝相回落到顆粒表面而轉(zhuǎn)變?yōu)椴灰?guī)則的非球體[1-2],導(dǎo)致動力學(xué)特性發(fā)生一系列改變,從而影響整個多相流場的流動特性。但經(jīng)典的顆粒受力理論是以規(guī)則球形為基礎(chǔ)進(jìn)行推導(dǎo)的,實際應(yīng)用中將不可避免地產(chǎn)生偏差[3]。另外,實際流場中,顆粒的運動狀態(tài)往往是轉(zhuǎn)動伴隨平動同時存在,但為簡化計算,數(shù)值模擬中,前者往往被忽略,理論上講,顆粒在發(fā)動機(jī)中所受阻力與相對流速度呈平方反比關(guān)系。所以,流速變化對顆粒的受力影響很大。固體火箭發(fā)動機(jī)中,在推進(jìn)劑剛開始燃燒時,燃?xì)馑俣容^小,但幾十毫秒后,燃?xì)鈯A帶著顆粒向噴管加速流動,由于慣性作用,顆粒相速度滯后,與主燃?xì)忾g存在著速度差,再加上自身旋轉(zhuǎn)作用,受力狀態(tài)較為復(fù)雜。在固體火箭發(fā)動機(jī)內(nèi)流場中,金屬顆粒的運動特性對整個工作過程起至關(guān)重要的作用,綜合考慮實際環(huán)境中可能存在的顆粒形態(tài)及運動狀態(tài),針對旋轉(zhuǎn)非球形單顆粒的受力問題進(jìn)行詳細(xì)數(shù)值分析,對進(jìn)一步分析顆粒間凝聚、碰撞等運動特性以及顆粒與流體間相互作用有一定的參考作用,對進(jìn)一步認(rèn)清多相流的本質(zhì)特征、完善多相流模型有重要意義。

    1 物理模型

    金屬鋁極易氧化,鋁顆粒存放期間,其表面就包裹著氧化物,初始進(jìn)入發(fā)動機(jī)時,顆粒形狀已經(jīng)發(fā)生變化,不再是規(guī)則的球形[4],受力狀態(tài)與運動特性亦發(fā)生相應(yīng)改變。Merrill[5]研究了鋁顆粒在固體火箭發(fā)動機(jī)環(huán)境中的相變?nèi)紵^程,對其結(jié)構(gòu)變化過程做了仔細(xì)地推導(dǎo)和計算。結(jié)果表明,顆粒將在90 ms內(nèi)完成燃燒過程,其間主要以Al和Al2O3兩球形相切的方式出現(xiàn),隨燃燒的發(fā)展,外形不斷變化,如圖1所示。數(shù)值模擬中,通常以等體積球體代替非球體進(jìn)行受力分析。對Merrill教授的研究結(jié)果進(jìn)行統(tǒng)計可知,隨鋁滴在燃燒室中狀態(tài)的改變,等體積顆粒在來流方向投影的橫截面面積偏離真實顆粒的水平較大,最大相差49.9%,即便在顆粒初始進(jìn)入流場時,也有1.1%的區(qū)別,這將對顆粒的受力分析產(chǎn)生極大誤差。

    (a)t=0 ms (b)t=30 ms

    (c)t=60 ms (d)t=90 ms

    以Merrill的研究為基礎(chǔ),利用Fluent商用軟件中的滑移網(wǎng)格技術(shù),對典型階段下旋轉(zhuǎn)鋁滴的受力狀態(tài)進(jìn)行分析。鋁滴顆粒初始半徑R=50 μm,計算區(qū)域取得足夠大,以保證顆粒的運動特性不受邊界流動影響:平行于來流方向的長度l=3 mm,垂直于來流方向的寬度d為2 mm。顆粒置于流場中靠近來流方向處,中心與入口邊界相距1.1 mm。采用非正交網(wǎng)格進(jìn)行計算,近壁面處做邊界層處理,顆粒附近區(qū)域進(jìn)行了加密處理,最小網(wǎng)格尺寸為5 μm,最大為170 μm,總網(wǎng)格數(shù)約為35萬,如圖2所示。為便于分析,記來流v的方向為x,顆粒繞z軸順時針旋轉(zhuǎn),軸線與來流方向瞬時針夾角為θ,如圖3所示。

    (a)三維模型 (b)中心截面圖

    圖3 θ的表示方法Fig.3 Representation of θ

    2 控制方程

    假設(shè)來流是穩(wěn)態(tài)的,其物性在運動過程中不發(fā)生變化,計算的控制方程為

    ·U=0

    (1)

    (2)

    為簡化模擬條件,不考慮氣固兩相間的傳熱,以空氣代替發(fā)動機(jī)內(nèi)的燃?xì)?,顆粒所受x方向阻力系數(shù)CD由顆粒所受到的阻力FD定義[6]:

    (3)

    同樣,顆粒所受到的與來流垂直方向的阻力系數(shù),即升力系數(shù)CL,由升力FL定義:

    (4)

    計算中,顆粒僅有轉(zhuǎn)動,無平動,顆粒雷諾數(shù)Rep小于100,屬層流不可壓縮流動。采用三階精度Quick差分離散格式對方程組離散,利用Simplec算法求解控制方程組。收斂判據(jù)為殘差和小于10-6,迭代時間步長為10-3。

    3 結(jié)果與分析

    3.1 模型驗證

    對于規(guī)則球形顆粒的旋轉(zhuǎn)受力情況,已有學(xué)者進(jìn)行了大量研究,并總結(jié)出相應(yīng)的升力經(jīng)驗公式CL=2Γ[7],Γ為無量綱轉(zhuǎn)速Γ=rω/vr,ω為顆粒旋轉(zhuǎn)角速度。無量綱轉(zhuǎn)速是顆粒角速度與來流速度的比值,可更好地反映顆粒的旋轉(zhuǎn)狀態(tài)。為驗證模型的適用性,針對典型工況進(jìn)行計算,并與經(jīng)驗公式對比。取初始半徑R=50 μm的球形顆粒進(jìn)行驗證計算,無量綱轉(zhuǎn)速Γ分別為0.1、1、2、3。模擬得顆粒的阻力系數(shù)對比結(jié)果如圖4所示。當(dāng)雷諾數(shù)較小時,計算結(jié)果與理論解的最大誤差不超過5%,二者符合良好,證明文中所用模型是正確的。以此模型為計算基礎(chǔ),對其他狀態(tài)下的旋轉(zhuǎn)顆粒進(jìn)行了受力分析。

    圖4 升力系數(shù)模擬結(jié)果與理論值的比較Fig.4 Simulation result vs theoretical value

    3.2 相對流速變化對顆粒受力的影響

    顆粒旋轉(zhuǎn)速度發(fā)生改變時,會影響因旋轉(zhuǎn)而產(chǎn)生的升力,也導(dǎo)致顆粒表面壓力梯度的變化。選取t=90 ms的典型狀態(tài),對比計算轉(zhuǎn)速分別為5、50、100、200 rad/s時顆粒的受力情況。計算結(jié)果如圖5所示。

    計算表明,旋轉(zhuǎn)速度對顆粒受力的影響較小,各大轉(zhuǎn)速下,顆粒的升力系數(shù)與阻力系數(shù)相差較小。因此,以鋁滴燃盡狀態(tài)(t=90 ms)為基礎(chǔ),取50 rad/s轉(zhuǎn)速狀態(tài),模擬計算相對流速vr分別為0.5、1.0、1.5、5.0 m/s時顆粒的受力情況,對比分析流場對顆粒受力的影響。4種工況中,顆粒雷諾數(shù)Re分別為0.035、0.07、0.14、0.35,無量綱轉(zhuǎn)速Γ分別為0.007、0.003 5、0.002 3、0.000 7。計算結(jié)果表明,非球顆粒在流場中的升力系數(shù)與無量綱轉(zhuǎn)速之間已經(jīng)不再滿足CL=2Γ的關(guān)系。隨相對流速增加,顆粒的無量綱轉(zhuǎn)速減小,阻力系數(shù)也隨之減小。圖6~圖8詳細(xì)反映了顆粒在一個旋轉(zhuǎn)周期內(nèi)不同方位角處的受力情況。

    圖6反映了當(dāng)相對流速由0.5 m/s逐漸變化到5 m/s時,顆粒在y方向所受升力系數(shù)的變化過程。當(dāng)來流速度相對較小時,在一個旋轉(zhuǎn)周期內(nèi),顆粒的升力系數(shù)變化較為明顯。周向不同位置處存在著2組極大值與極小值,它們交替出現(xiàn),每對極值所在位置相差約180°,2組極值的連線交叉近似“十”字形,且阻力系數(shù)在負(fù)向的最大振幅大于正向,這一點在圖6(a)表現(xiàn)得尤為明顯。對比圖6(a) ~(d)可看出,各工況下,阻力系數(shù)極值出現(xiàn)的位置固定,與相對流速無關(guān)。隨來流速度的增加,正負(fù)向極值分別迅速衰減和增加,周向上阻力系數(shù)變得更加均衡,均逐漸逼近于零。由圖6(d)可看出,當(dāng)流速增加到5 m/s時,周向上各位置處的阻力系數(shù)均在0值環(huán)線附近,說明此時流場對顆粒在x向所產(chǎn)生的阻力作用已經(jīng)很弱。

    (a) 升力系數(shù)對比

    (b) 阻力系數(shù)對比

    圖7反映了當(dāng)相對流速由0.5 m/s逐漸變化到5 m/s時,顆粒在x方向所受阻力系數(shù)的變化過程。與x向變化規(guī)律類似,在一個旋轉(zhuǎn)周期內(nèi),周向不同位置處存在著2組交替出現(xiàn)的極值,每對極值所處位置也相差約180°。與x向不同的是顆粒所受到的y向阻力系數(shù)均為正值,且周向變化較為平緩,曲線近似橢圓形。相對而言,水平位置處極值稍大于垂直位置,但二者差別不大。流速對顆粒阻力影響很大,由圖7(a)~(d) 明顯可見,阻力系數(shù)迅速減小,當(dāng)相對流速由0.5 m/s增加至5 m/s時,阻力系數(shù)均值由614.3銳降到64.5。流速越大,對顆粒產(chǎn)生的阻力作用越弱。

    圖8(a)、(b)分別直觀反映了顆粒在周向不同方位時所受x方向阻力和y方向升力的波動情況。顆粒在x向的阻力系數(shù)變化接近“正弦”曲線,數(shù)值正負(fù)交替出現(xiàn),說明顆粒在不同位置時所受到x向的阻力方向是變化的,且相對流速越大時,相對應(yīng)的“振幅”越小,曲線越加光滑。當(dāng)相對流速增加到5 m/s時,曲線已接近于數(shù)值為零的直線,比流速為0.5 m/s時的平均值銳減了89%。在顆粒主軸與來流約成90°和270°兩位置處,顆粒所受x向的阻力系數(shù)均為零。顆粒在y向的阻力系數(shù)曲線波動較小,相對流速愈大,曲線愈加平滑,當(dāng)流速增加到5 m/s時,曲線已經(jīng)接近于直線。與x向阻力變化類似,y向阻力均值也銳減了89%。

    (a) vr=0.5 m/s

    (b) vr=1 m/s

    (c) vr=2 m/s

    (d) vr=5 m/s

    (a) vr=0.5 m/s

    (b) vr=1 m/s

    (c) vr=2 m/s

    (d) vr=5 m/s

    3.3 顆粒升力系數(shù)和阻力系數(shù)與顆粒雷諾數(shù)的關(guān)系

    計算表明,顆粒的阻力系數(shù)與來流相對速度成正比,當(dāng)物性條件確定時,顆粒的雷諾數(shù)Re只與相對流速有關(guān)。當(dāng)旋轉(zhuǎn)速度確定時,顆粒的無量綱轉(zhuǎn)速Γ也只與相對來流速度有關(guān),即與雷諾數(shù)Re相關(guān),由此可得到阻力系數(shù)與雷諾數(shù)Re的函數(shù)關(guān)系,即CD=f(Re)。為探索二者的關(guān)系,首先統(tǒng)計顆x向無量綱相對阻力系數(shù)與相對雷諾數(shù)的比值,如表1所示。

    (a) 升力系數(shù)對比

    (b) 阻力系數(shù)對比

    表1 顆粒相對阻力系數(shù)與相對雷諾數(shù)的比值統(tǒng)計Table 1 Statistics of the numberless drag coefficient and Reynolds number

    由表1可知,當(dāng)旋轉(zhuǎn)角速度恒定時,顆粒的相對阻力系數(shù)與顆粒雷諾數(shù)之比值滿足:

    (5)

    式中CD1和CD2與Re1和Re2代表2種不同工況下的阻力系數(shù)和雷諾數(shù);α為一常數(shù)。

    由此可推斷,顆粒阻力系數(shù)與雷諾數(shù)滿足以下關(guān)系:

    CD=βRe-1

    (6)

    式中β為常數(shù)。

    引入球形度(Sphericity)的概念[8-9],ψ=(dv/ds),可得到顆粒的球形度為0.912。Haider和Chhabra[10]曾提出球形度與阻力系數(shù)間存在著一定的關(guān)系:

    (7)

    其中

    A=exp(2.328 8-6.458 1ψ+2.448 6ψ2)

    B=0.096 4+0.556 5ψ

    C=exp(4.905-13.894 4ψ+18.422 2ψ2-10.259 9ψ3)

    D=exp(1.468 1+12.258 4ψ-20.732 2ψ2+15.885 5ψ3)

    當(dāng)顆粒的球形度ψ=0.912時,可求得顆粒的阻力系數(shù)與雷諾數(shù)間關(guān)系為

    CD=24.48Re-1

    (8)

    此式與式(6)完全符合,此時β取為24.48,這也再次印證了本文計算所得結(jié)果可信。

    統(tǒng)計顆粒y向升力系數(shù)計算結(jié)果,得到表2。與表1的統(tǒng)計結(jié)果相比,表2中旋轉(zhuǎn)周向不同位置處,顆粒的無量綱升力系數(shù)與雷諾數(shù)比值存在著差別,并不是常數(shù)關(guān)系。

    表2 顆粒相對升力系數(shù)與相對雷諾數(shù)的比值統(tǒng)計Table2 Statistics of the numberless lift coefficient and Reynolds number

    式(6)應(yīng)修正為

    (9)

    即比例系數(shù)α是方位角θ的函數(shù),α=f(θ)。統(tǒng)計表明,當(dāng)顆粒雷諾數(shù)Re=0.14~0.35時,α值為0.8~1.17。此時,顆粒升力系數(shù)與雷諾數(shù)Re之間的關(guān)系也應(yīng)修正為

    CL=C(θ)Re-1

    (10)

    4 結(jié)論

    (1)固體火箭發(fā)動機(jī)燃燒室環(huán)境中,非球形鋁滴的旋轉(zhuǎn)受力與經(jīng)典球形顆粒的理論解有很大區(qū)別,顆粒所受x向阻力系數(shù)與y向升力系數(shù)變化規(guī)律差別較大,阻力系數(shù)明顯大于升力系數(shù)。小雷諾數(shù)條件下,一個旋轉(zhuǎn)周期內(nèi),周向不同位置處的升力系數(shù)方向發(fā)生變化,存在著交替出現(xiàn)的兩組極大值與極小值。其中,負(fù)向振幅大于正向,阻力系數(shù)極值出現(xiàn)的位置固定,與相對流速無關(guān);阻力系數(shù)均為正,但周向變化較為平緩,各方位角處數(shù)值相差不大。隨來流速度增加,升力系數(shù)變得均衡,系數(shù)值逐漸逼近于零;阻力系數(shù)也急速降低。當(dāng)相對流速由0.5 m/s提高到5 m/s時,顆粒所受阻力銳降89%。流速越大,對顆粒產(chǎn)生的阻力作用越弱。

    (2)顆粒的相對升力系數(shù)與顆粒雷諾數(shù)之比值滿足:CL1/CL2=α(θ)Re2/Re1,比例系數(shù)α是方位角θ的函數(shù)。當(dāng)顆粒雷諾數(shù)Re=0.14~0.35時,升力系數(shù)中α值為0.8~1.17,阻力系數(shù)中α=1。顆粒阻力系數(shù)與雷諾數(shù)滿足關(guān)系:CD=β(θ)Re-1。其中,x向阻力系數(shù)中,β為常數(shù);y向升力系數(shù)中,β是方位角θ的函數(shù)。

    [1] 方丁酉.兩相流體力學(xué)[M]. 長沙:國防科技大學(xué)出版社,1988.

    [2] Olsen S E, Beckstead M W. Burn time measurements of single aluminum particles in steam and CO2mixtures[J]. Journal of Propulsion and Power, 1996, 12(4):662-671.

    [3] 林建中. 超常顆粒多相流體動力學(xué)——圓柱狀顆粒兩相流[M]. 北京:科學(xué)出版社, 2008.

    [4] Robert Geisler. A global view of the use of aluminum fuel in solid rocket motors[R].AIAA 2002-3748.

    [5] Merrill K King. Aluminum combustion in a solid rocket motor environment[J]. Proceedings of the Combustion Institute, 2009, 32(2):2107-2114.

    [6] Ounis H, Ahmadi G, McLaughlin J B. Brownian diffusion of submicrometer particles in the viscous sublayer[J]. Journal of Colloid and Interface Science, 1991, 143(1):266-277.

    [7] Rubinow S I, Keller J B. The transerverse force on a spinning sphere moving in a viscous fluid[J]. Fluid Mechanics, 1961, 11: 447-459.

    [8] Lu Hui-lin , Liu Wen-tie, Zhao Guang-bo. Computational modeling of dense gas particle flow in a pipe: kinetic theory approach of granular flow[J]. Journal of Chemical Industry and Engineering ( China), 2000, 5(1):31-38.

    [9] Gan Lin, Xu Mao-sheng, Zhu Bing-chen. Heat transfer parameters of packed bed with ring pellet[J]. Journal of Chemical Industry and Engineering ( China), 2000, 5(6):778-783.

    [10] Chhabra R P, Agarwal L, Sinha N Kl. Drag on non spherical particles: an evaluation of available methods[J].Powder Technology,1999,101:288-295.

    (編輯:崔賢彬)

    Influence of relative velocity on the force of non-spherical rotating aluminum droplet

    LIU Cong-lin1, GAO Ye1, HE Zheng1, KANG Wen-wu2

    (1. Aerospace and Civil Engineering, Harbin Engineering University, Harbin 150001, China; 2. School of Energy and Power Engineering, Beihang University, Beijing 100191,China)

    The rotation characteristics of non-sphere aluminum in different kinds of fluid velocity were simulated in order to analyze the effect of flow rate on the forces acting on particles. Results show that, the relationship between lift coefficient and dimensionless rotational speed doesn't obey the typical law. When Reynolds number is less than 1, two groups of maximum and minimum values appear alternately on certain circumferential location for drag coefficient. Lift coefficient is always positive and smooth. With the increasing of the inflow velocity, drag coefficient becomes even, and the value is approaching to zero, while the lift coefficient falls quickly. If the relative velocity increase from 0.5 m/s to 5 m/s, drag coefficient of the droplet decreases by 89% sharply. Drag coefficient and Reynolds number obey a new law.

    non-spherical aluminum droplet; relative velocity; drag coefficient; lift coefficient

    2014-05-06;

    :2014-06-18。

    國家自然科學(xué)基金(11372079);中央高?;究蒲袠I(yè)務(wù)費專項資金(HEUCF 130203)。

    劉叢林(1981—),女,助理研究員,研究方向為多相流內(nèi)流場燃燒。E-mail:liuconglin2006@126.com

    V435

    A

    1006-2793(2015)02-0214-06

    10.7673/j.issn.1006-2793.2015.02.012

    猜你喜歡
    來流雷諾數(shù)升力
    高速列車車頂–升力翼組合體氣動特性
    兩種典型來流條件下風(fēng)力機(jī)尾跡特性的數(shù)值研究
    能源工程(2022年2期)2022-05-23 13:51:48
    無人機(jī)升力測試裝置設(shè)計及誤差因素分析
    不同來流條件對溢洪道過流能力的影響
    基于自適應(yīng)偽譜法的升力式飛行器火星進(jìn)入段快速軌跡優(yōu)化
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    升力式再入飛行器體襟翼姿態(tài)控制方法
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計研究
    民機(jī)高速風(fēng)洞試驗的阻力雷諾數(shù)效應(yīng)修正
    丰满饥渴人妻一区二区三| 桃花免费在线播放| 国产熟女午夜一区二区三区 | 久久人妻熟女aⅴ| 国产免费一级a男人的天堂| 交换朋友夫妻互换小说| 人体艺术视频欧美日本| 亚洲丝袜综合中文字幕| 成年女人在线观看亚洲视频| 99久久人妻综合| 丰满饥渴人妻一区二区三| 国产乱来视频区| 婷婷色综合www| 人人妻人人添人人爽欧美一区卜| 欧美另类一区| 国产亚洲欧美精品永久| 亚州av有码| 日韩中文字幕视频在线看片| 亚洲av.av天堂| 亚洲怡红院男人天堂| 看非洲黑人一级黄片| 男女啪啪激烈高潮av片| 亚洲精品aⅴ在线观看| 亚洲成人av在线免费| 国产精品.久久久| 伦理电影大哥的女人| 黄色配什么色好看| 国产精品久久久久久久久免| 免费日韩欧美在线观看| 国产免费现黄频在线看| 国产伦理片在线播放av一区| 久久久久久久亚洲中文字幕| 99久久人妻综合| 新久久久久国产一级毛片| 午夜福利视频在线观看免费| 涩涩av久久男人的天堂| 亚洲第一av免费看| 高清欧美精品videossex| 国产黄色免费在线视频| 这个男人来自地球电影免费观看 | 涩涩av久久男人的天堂| 国产亚洲一区二区精品| 麻豆精品久久久久久蜜桃| 国产亚洲一区二区精品| 在线观看www视频免费| 特大巨黑吊av在线直播| 国产无遮挡羞羞视频在线观看| 欧美丝袜亚洲另类| 亚洲av不卡在线观看| 久久精品国产鲁丝片午夜精品| 一边摸一边做爽爽视频免费| 简卡轻食公司| 亚洲国产精品成人久久小说| 我要看黄色一级片免费的| 国产在视频线精品| 高清欧美精品videossex| 久久国内精品自在自线图片| 女性被躁到高潮视频| 一本大道久久a久久精品| 久久青草综合色| 免费观看无遮挡的男女| 91精品伊人久久大香线蕉| 久久国内精品自在自线图片| 制服人妻中文乱码| 欧美少妇被猛烈插入视频| 黑人猛操日本美女一级片| 精品久久国产蜜桃| 国产亚洲最大av| av天堂久久9| 在线观看国产h片| 欧美少妇被猛烈插入视频| 校园人妻丝袜中文字幕| 成人漫画全彩无遮挡| 中国国产av一级| av卡一久久| 最近手机中文字幕大全| 色婷婷久久久亚洲欧美| 国产黄色视频一区二区在线观看| 亚洲国产成人一精品久久久| av线在线观看网站| 97在线人人人人妻| 亚洲国产精品999| 高清午夜精品一区二区三区| 亚洲精品久久久久久婷婷小说| 国产极品粉嫩免费观看在线 | 黑人猛操日本美女一级片| 最近2019中文字幕mv第一页| 91精品三级在线观看| 日本av免费视频播放| 精品少妇久久久久久888优播| 欧美 日韩 精品 国产| 日韩熟女老妇一区二区性免费视频| 少妇丰满av| 国产精品久久久久久av不卡| 亚洲国产毛片av蜜桃av| 亚洲av国产av综合av卡| 亚洲欧美日韩卡通动漫| 国产精品麻豆人妻色哟哟久久| 国产欧美亚洲国产| 天天操日日干夜夜撸| 最近中文字幕2019免费版| 亚洲丝袜综合中文字幕| 肉色欧美久久久久久久蜜桃| 国产成人91sexporn| 亚洲成人手机| 春色校园在线视频观看| 亚洲国产毛片av蜜桃av| h视频一区二区三区| 寂寞人妻少妇视频99o| 日本av免费视频播放| 日本猛色少妇xxxxx猛交久久| 亚洲高清免费不卡视频| 日韩一区二区三区影片| 老熟女久久久| 高清视频免费观看一区二区| 亚洲激情五月婷婷啪啪| 2018国产大陆天天弄谢| 老司机亚洲免费影院| 乱码一卡2卡4卡精品| 国产成人精品福利久久| 亚洲美女黄色视频免费看| 国产精品一区二区在线不卡| 亚洲欧美一区二区三区国产| 国产色爽女视频免费观看| 亚洲欧美成人精品一区二区| 99国产精品免费福利视频| 亚洲av国产av综合av卡| 日韩中文字幕视频在线看片| 啦啦啦视频在线资源免费观看| 一区二区三区免费毛片| 亚洲人与动物交配视频| 精品99又大又爽又粗少妇毛片| tube8黄色片| 欧美日韩一区二区视频在线观看视频在线| 美女中出高潮动态图| 精品人妻熟女毛片av久久网站| 欧美精品人与动牲交sv欧美| 99热网站在线观看| 99久久综合免费| 亚洲av不卡在线观看| 午夜视频国产福利| 纵有疾风起免费观看全集完整版| 国产精品秋霞免费鲁丝片| 91在线精品国自产拍蜜月| 五月天丁香电影| 午夜激情福利司机影院| 丝袜在线中文字幕| 热99久久久久精品小说推荐| 永久网站在线| 亚洲怡红院男人天堂| 啦啦啦中文免费视频观看日本| 久久人人爽人人片av| 日韩视频在线欧美| 国产熟女欧美一区二区| 99热国产这里只有精品6| 丰满迷人的少妇在线观看| 简卡轻食公司| 免费人成在线观看视频色| 成人毛片60女人毛片免费| 亚洲国产欧美在线一区| 久久人人爽人人片av| 成人免费观看视频高清| 少妇的逼好多水| 人妻少妇偷人精品九色| 国产淫语在线视频| 欧美激情 高清一区二区三区| 夜夜看夜夜爽夜夜摸| av网站免费在线观看视频| av免费观看日本| 精品久久蜜臀av无| 色视频在线一区二区三区| 久久 成人 亚洲| 99热6这里只有精品| 中文字幕免费在线视频6| 国产成人精品久久久久久| 在线观看国产h片| 国产深夜福利视频在线观看| 久久久久久久久久久久大奶| 99九九在线精品视频| h视频一区二区三区| 亚洲天堂av无毛| 春色校园在线视频观看| 国产不卡av网站在线观看| 日韩伦理黄色片| 一本一本综合久久| 伦理电影大哥的女人| 大片电影免费在线观看免费| 国产爽快片一区二区三区| 亚洲国产精品一区二区三区在线| 国产成人免费无遮挡视频| 久久精品熟女亚洲av麻豆精品| 国产精品偷伦视频观看了| 大片免费播放器 马上看| 又黄又爽又刺激的免费视频.| 国产精品国产三级专区第一集| 亚洲欧美日韩卡通动漫| 国产精品国产av在线观看| 视频区图区小说| 亚洲国产精品成人久久小说| 国产男女超爽视频在线观看| 97在线视频观看| 成年av动漫网址| 精品人妻偷拍中文字幕| av卡一久久| 另类亚洲欧美激情| 人人妻人人澡人人看| 久久久久精品性色| videosex国产| 9色porny在线观看| 日韩精品免费视频一区二区三区 | 午夜福利在线观看免费完整高清在| 十分钟在线观看高清视频www| 国产免费一区二区三区四区乱码| 各种免费的搞黄视频| 国产亚洲欧美精品永久| 男男h啪啪无遮挡| 一边摸一边做爽爽视频免费| 久久99热6这里只有精品| 亚洲国产欧美日韩在线播放| 亚洲av成人精品一区久久| 考比视频在线观看| 免费观看av网站的网址| 精品国产一区二区久久| 丰满少妇做爰视频| 涩涩av久久男人的天堂| 考比视频在线观看| 老司机影院毛片| 最新的欧美精品一区二区| 国产成人91sexporn| 日韩成人av中文字幕在线观看| 美女cb高潮喷水在线观看| 久久久久国产精品人妻一区二区| 在线亚洲精品国产二区图片欧美 | 久久久午夜欧美精品| 美女国产视频在线观看| xxx大片免费视频| 国产亚洲欧美精品永久| 欧美性感艳星| 亚洲av福利一区| 王馨瑶露胸无遮挡在线观看| 免费黄网站久久成人精品| 国产男女超爽视频在线观看| 日韩大片免费观看网站| av电影中文网址| 欧美97在线视频| 男的添女的下面高潮视频| 国产亚洲午夜精品一区二区久久| 中文乱码字字幕精品一区二区三区| 亚洲丝袜综合中文字幕| 亚洲av日韩在线播放| 一级片'在线观看视频| 国产成人aa在线观看| 大陆偷拍与自拍| 国产成人av激情在线播放 | 国产熟女午夜一区二区三区 | 精品国产乱码久久久久久小说| 黄色视频在线播放观看不卡| 日本黄大片高清| 91精品伊人久久大香线蕉| 亚洲精品一二三| 女人精品久久久久毛片| 狠狠精品人妻久久久久久综合| 大话2 男鬼变身卡| 中文字幕av电影在线播放| 亚洲人成网站在线观看播放| 亚洲内射少妇av| 国产乱人偷精品视频| 亚洲色图 男人天堂 中文字幕 | 国产一级毛片在线| 欧美精品人与动牲交sv欧美| 日韩亚洲欧美综合| 亚洲av二区三区四区| 色哟哟·www| 新久久久久国产一级毛片| 日本午夜av视频| 51国产日韩欧美| 日本免费在线观看一区| 最后的刺客免费高清国语| 国产毛片在线视频| 在线观看国产h片| 日韩电影二区| 国产成人精品婷婷| 欧美一级a爱片免费观看看| 国产伦精品一区二区三区视频9| 另类精品久久| 只有这里有精品99| 一边摸一边做爽爽视频免费| 亚洲精品视频女| 色婷婷久久久亚洲欧美| 夜夜看夜夜爽夜夜摸| 91精品国产国语对白视频| 人妻一区二区av| 夫妻午夜视频| 久久这里有精品视频免费| 免费不卡的大黄色大毛片视频在线观看| 精品99又大又爽又粗少妇毛片| 亚洲精品,欧美精品| 人妻少妇偷人精品九色| 中文天堂在线官网| 高清欧美精品videossex| 成人亚洲精品一区在线观看| 久久久久久久精品精品| 国产探花极品一区二区| 国产综合精华液| tube8黄色片| 2022亚洲国产成人精品| 久久人人爽av亚洲精品天堂| 秋霞在线观看毛片| 麻豆成人av视频| 少妇被粗大猛烈的视频| 嘟嘟电影网在线观看| 2022亚洲国产成人精品| 国产日韩欧美视频二区| 精品人妻偷拍中文字幕| 最近2019中文字幕mv第一页| 91午夜精品亚洲一区二区三区| 色婷婷久久久亚洲欧美| 国产成人精品久久久久久| 日韩欧美精品免费久久| 亚洲成人av在线免费| 久久 成人 亚洲| 免费看光身美女| 亚洲情色 制服丝袜| 考比视频在线观看| 日韩强制内射视频| 国产精品不卡视频一区二区| av天堂久久9| 少妇的逼水好多| 久久久国产一区二区| 欧美精品高潮呻吟av久久| 国产成人91sexporn| 熟妇人妻不卡中文字幕| 又大又黄又爽视频免费| 亚洲图色成人| 精品一区二区免费观看| 十八禁高潮呻吟视频| 亚洲av成人精品一二三区| 亚洲国产欧美日韩在线播放| 日韩伦理黄色片| 亚洲国产欧美日韩在线播放| 免费观看a级毛片全部| 国产乱来视频区| 婷婷色麻豆天堂久久| 亚洲人与动物交配视频| 2022亚洲国产成人精品| 日本黄色片子视频| 免费高清在线观看视频在线观看| 青春草国产在线视频| 最新的欧美精品一区二区| 久久精品夜色国产| 日产精品乱码卡一卡2卡三| 少妇 在线观看| 91午夜精品亚洲一区二区三区| 水蜜桃什么品种好| 最新的欧美精品一区二区| 97精品久久久久久久久久精品| a级毛片在线看网站| 黄色视频在线播放观看不卡| 亚洲精品日韩在线中文字幕| 丰满饥渴人妻一区二区三| av在线观看视频网站免费| 日韩一区二区三区影片| 欧美精品一区二区免费开放| 香蕉精品网在线| 欧美日韩视频精品一区| 久久久久久久国产电影| 精品一区二区免费观看| 国产高清三级在线| 久久人人爽人人爽人人片va| av视频免费观看在线观看| 国产男女内射视频| 久久久欧美国产精品| 国产高清国产精品国产三级| 久久久久久久久久成人| 黄色毛片三级朝国网站| 夫妻性生交免费视频一级片| 一本大道久久a久久精品| 中文字幕人妻丝袜制服| 成年人午夜在线观看视频| 国产女主播在线喷水免费视频网站| 中文乱码字字幕精品一区二区三区| 午夜福利网站1000一区二区三区| 亚洲欧美一区二区三区国产| 在线观看免费日韩欧美大片 | 我的女老师完整版在线观看| 免费少妇av软件| 久久久国产一区二区| 免费人成在线观看视频色| 一级a做视频免费观看| 考比视频在线观看| 青青草视频在线视频观看| 91精品三级在线观看| 亚洲色图综合在线观看| 免费av不卡在线播放| 能在线免费看毛片的网站| 欧美日韩一区二区视频在线观看视频在线| 女性被躁到高潮视频| 人妻夜夜爽99麻豆av| 黄色视频在线播放观看不卡| 超色免费av| 婷婷成人精品国产| 亚洲怡红院男人天堂| 国产亚洲精品第一综合不卡 | 中文乱码字字幕精品一区二区三区| 天天操日日干夜夜撸| 国产黄色视频一区二区在线观看| av播播在线观看一区| 亚洲精品456在线播放app| av免费观看日本| 自线自在国产av| av免费在线看不卡| 成年女人在线观看亚洲视频| 成人漫画全彩无遮挡| 精品人妻熟女毛片av久久网站| 国精品久久久久久国模美| .国产精品久久| 亚洲精品国产av成人精品| 亚洲av日韩在线播放| 水蜜桃什么品种好| 午夜av观看不卡| 国产爽快片一区二区三区| 久久久久久伊人网av| 国产高清三级在线| a级毛片免费高清观看在线播放| 另类亚洲欧美激情| 亚洲,一卡二卡三卡| 亚洲人成77777在线视频| 中文字幕精品免费在线观看视频 | 免费av中文字幕在线| 尾随美女入室| 国产精品久久久久久av不卡| 亚洲欧美精品自产自拍| 黄色视频在线播放观看不卡| 国产一区二区三区av在线| 国产精品蜜桃在线观看| 最近最新中文字幕免费大全7| 国产色婷婷99| 在线观看国产h片| 国产欧美亚洲国产| 高清欧美精品videossex| 美女大奶头黄色视频| 天堂俺去俺来也www色官网| av福利片在线| av在线老鸭窝| 日韩强制内射视频| 五月开心婷婷网| 一本久久精品| 亚洲精品自拍成人| 国产欧美日韩一区二区三区在线 | 国产精品国产av在线观看| 国产成人av激情在线播放 | 黄色怎么调成土黄色| 国产精品人妻久久久影院| 黑人巨大精品欧美一区二区蜜桃 | 国产精品熟女久久久久浪| 一级,二级,三级黄色视频| 亚洲精品中文字幕在线视频| 亚洲av男天堂| 久热这里只有精品99| 国产综合精华液| 成人亚洲精品一区在线观看| 国产片特级美女逼逼视频| av在线观看视频网站免费| 国产爽快片一区二区三区| 一本—道久久a久久精品蜜桃钙片| av专区在线播放| 午夜免费鲁丝| 精品国产露脸久久av麻豆| a 毛片基地| 欧美日韩视频精品一区| 精品亚洲成国产av| 搡老乐熟女国产| 国产欧美亚洲国产| 在现免费观看毛片| 91aial.com中文字幕在线观看| 久久精品国产亚洲av天美| xxx大片免费视频| a级片在线免费高清观看视频| 中文字幕av电影在线播放| 欧美日韩一区二区视频在线观看视频在线| 亚洲五月色婷婷综合| 少妇被粗大猛烈的视频| 黄色欧美视频在线观看| 中文字幕最新亚洲高清| 色94色欧美一区二区| 久久久久网色| 午夜av观看不卡| 中文字幕制服av| 精品国产乱码久久久久久小说| 91精品国产九色| 91午夜精品亚洲一区二区三区| 国产精品熟女久久久久浪| 久久久久久久久久人人人人人人| 黄色视频在线播放观看不卡| 3wmmmm亚洲av在线观看| 国产免费又黄又爽又色| 日韩成人av中文字幕在线观看| 人体艺术视频欧美日本| 一级黄片播放器| 26uuu在线亚洲综合色| 丝袜喷水一区| 9色porny在线观看| 全区人妻精品视频| 久久精品国产亚洲网站| 亚洲综合色惰| 亚洲精品一二三| 亚洲婷婷狠狠爱综合网| 欧美精品人与动牲交sv欧美| 在线亚洲精品国产二区图片欧美 | 亚洲性久久影院| 最后的刺客免费高清国语| 亚洲一区二区三区欧美精品| av线在线观看网站| 激情五月婷婷亚洲| 性高湖久久久久久久久免费观看| 精品一区二区三区视频在线| 婷婷色综合大香蕉| 日韩电影二区| 黄色毛片三级朝国网站| av网站免费在线观看视频| 久久久久久久大尺度免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 在线精品无人区一区二区三| 欧美 亚洲 国产 日韩一| 免费人妻精品一区二区三区视频| 91成人精品电影| 丰满迷人的少妇在线观看| 精品国产一区二区三区久久久樱花| 亚洲精品第二区| 精品久久久久久电影网| 蜜桃在线观看..| 丰满少妇做爰视频| 大香蕉久久成人网| 国产高清三级在线| 色视频在线一区二区三区| 免费观看在线日韩| 国产精品久久久久久精品古装| 黑人高潮一二区| 久久精品国产亚洲av涩爱| 高清欧美精品videossex| 欧美精品一区二区免费开放| 亚洲第一av免费看| 一区二区三区四区激情视频| 男的添女的下面高潮视频| 99热这里只有精品一区| 色哟哟·www| 成人亚洲精品一区在线观看| 一级黄片播放器| 女人精品久久久久毛片| 亚洲性久久影院| av不卡在线播放| 国产熟女午夜一区二区三区 | 日本爱情动作片www.在线观看| 大话2 男鬼变身卡| 中文字幕免费在线视频6| 多毛熟女@视频| 欧美最新免费一区二区三区| 秋霞在线观看毛片| 一本色道久久久久久精品综合| 亚洲久久久国产精品| 老司机影院成人| 少妇丰满av| 久久毛片免费看一区二区三区| 欧美日韩国产mv在线观看视频| 日本与韩国留学比较| 亚洲第一av免费看| 久久 成人 亚洲| 亚洲精品自拍成人| 免费黄频网站在线观看国产| 最近中文字幕高清免费大全6| 天天操日日干夜夜撸| 欧美成人午夜免费资源| av天堂久久9| 欧美日韩一区二区视频在线观看视频在线| 蜜桃国产av成人99| 精品久久国产蜜桃| 欧美日韩综合久久久久久| 热re99久久精品国产66热6| 欧美精品国产亚洲| 成人免费观看视频高清| 国产日韩欧美视频二区| 日韩中文字幕视频在线看片| 黄色一级大片看看| 秋霞在线观看毛片| 精品久久蜜臀av无| 国产精品.久久久| 麻豆乱淫一区二区| 日韩三级伦理在线观看| 久久鲁丝午夜福利片| 狂野欧美白嫩少妇大欣赏| 一级毛片 在线播放| 国产在线视频一区二区| 日韩 亚洲 欧美在线| av免费在线看不卡| 久久鲁丝午夜福利片| 少妇 在线观看| 各种免费的搞黄视频| 国产成人精品福利久久| 麻豆乱淫一区二区| 亚洲综合色网址| 亚洲av二区三区四区| 久久久精品区二区三区| 国产国拍精品亚洲av在线观看| 亚洲经典国产精华液单| 啦啦啦啦在线视频资源| 国产精品免费大片| 插逼视频在线观看| 亚洲av不卡在线观看| 中文字幕最新亚洲高清| 寂寞人妻少妇视频99o| 丝瓜视频免费看黄片| 综合色丁香网| 精品人妻偷拍中文字幕| 国产色爽女视频免费观看| 久久 成人 亚洲|