999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

部分葉高串列擴(kuò)壓器對離心壓氣機(jī)性能的影響

2017-02-08 03:00:15黃月晴楊愛玲宋文杰陳二云
動力工程學(xué)報 2017年1期
關(guān)鍵詞:效率

黃月晴,楊愛玲,宋文杰,陳二云,戴 韌

(1.上海理工大學(xué) 能源與動力工程學(xué)院,上海市動力工程多相流動與傳熱重點(diǎn)實(shí)驗(yàn)室,上海 200093)

部分葉高串列擴(kuò)壓器對離心壓氣機(jī)性能的影響

黃月晴1,2,楊愛玲1,2,宋文杰1,2,陳二云1,2,戴 韌1,2

(1.上海理工大學(xué) 能源與動力工程學(xué)院,上海市動力工程多相流動與傳熱重點(diǎn)實(shí)驗(yàn)室,上海 200093)

為了保證離心壓氣機(jī)壓比同時拓寬其穩(wěn)定工作范圍,將部分葉高概念引入串列擴(kuò)壓器設(shè)計(jì)中,獲得部分葉高串列擴(kuò)壓器,并基于數(shù)值方法研究串列葉柵前后排葉片相對周向位置以及前排葉片相對高度對離心壓氣機(jī)與串列擴(kuò)壓器性能的影響.在相對周向位置為30%的串列擴(kuò)壓器基礎(chǔ)上,利用不同部分葉高葉片設(shè)計(jì)串列擴(kuò)壓器前排葉片,獲得不同的部分葉高串列擴(kuò)壓器.結(jié)果表明:串列葉柵前后排葉片相對周向位置對擴(kuò)壓器的擴(kuò)壓能力以及穩(wěn)定工作范圍有很大影響,當(dāng)相對周向位置處于20%~30%內(nèi)時離心壓氣機(jī)的整體性能達(dá)到最佳;與原離心壓氣機(jī)相比,當(dāng)前排葉片相對高度h/B分別為40%和50%時,離心壓氣機(jī)的喘振裕度可分別增加21%和25%,總壓比和等熵效率僅下降1%左右.

離心壓氣機(jī);串列擴(kuò)壓器;相對周向位置;部分葉高葉片;數(shù)值模擬

串列擴(kuò)壓器因其擴(kuò)壓能力強(qiáng)和高效率得到了國內(nèi)外學(xué)者的關(guān)注.Pampreen[1]設(shè)計(jì)了一種由3排葉柵組成的串列擴(kuò)壓器,并與通道型擴(kuò)壓器的性能進(jìn)行比較,發(fā)現(xiàn)串列擴(kuò)壓器不論是穩(wěn)定運(yùn)行范圍還是效率,都優(yōu)于通道型擴(kuò)壓器.Japikse[2]引用Pampreen未公開發(fā)表的實(shí)驗(yàn)數(shù)據(jù),經(jīng)分析發(fā)現(xiàn)當(dāng)?shù)诙湃~柵吸力面前緣相對于第一排葉柵壓力面尾緣相對周向位置(Relative Circumferential Position,RCP)為節(jié)距的10%時,串列葉柵的運(yùn)行效率最高.Seleznev等[3]通過對比單排葉片的葉柵擴(kuò)壓器和串列擴(kuò)壓器的性能,發(fā)現(xiàn)串列擴(kuò)壓器的效率高于單排葉片的葉柵擴(kuò)壓器,且當(dāng)相對周向位置為10%時,其效率最高.李紹斌等[4]通過數(shù)值方法研究了串列葉柵后排靜葉周向位置對壓氣機(jī)性能的影響,結(jié)果表明串列葉柵后排靜葉周向位置合理有利于抑制流動分離,降低損失.吳東坡等[5]對分別采用無葉擴(kuò)壓器和串列擴(kuò)壓器時的某高能頭半開式離心壓縮機(jī)的級進(jìn)行了流動模擬和性能分析,結(jié)果表明在高能頭系數(shù)的級中采用串列擴(kuò)壓器可以有效地改善葉輪出口流場,降低流動損失,提高壓力恢復(fù)系數(shù)和級效率.

另一方面,由于受到離心力與科氏力的作用,附面層內(nèi)低速氣流逐漸向葉片吸力面與輪蓋側(cè)堆積,導(dǎo)致此區(qū)域速度較低,而輪盤與壓力面?zhèn)人俣容^高,因此離心葉輪出口氣流沿周向和軸向呈非均勻分布.非均勻的氣流與下游擴(kuò)壓器葉片之間的動靜干涉使得擴(kuò)壓器進(jìn)口流場產(chǎn)生了非定常現(xiàn)象,不僅影響葉片擴(kuò)壓器的性能,甚至還可能影響離心壓氣機(jī)的整機(jī)性能.朱營康等[6]對安裝在輪蓋側(cè)的半葉高擴(kuò)壓器進(jìn)行了研究,發(fā)現(xiàn)其可以在較大的工況范圍內(nèi)提高擴(kuò)壓器的壓力恢復(fù)系數(shù),并且提高整個壓氣機(jī)級的總壓和多變效率.為了降低離心葉輪尾跡與葉片擴(kuò)壓器的相互干涉,筆者在串列擴(kuò)壓器的前排葉片引入部分葉高設(shè)計(jì).部分葉高概念最早由Yoshinaga等[7]于1987年提出,利用部分葉高葉片可以提高無葉擴(kuò)壓器的壓力恢復(fù)系數(shù),同時使沿葉高的速度分布更均勻.

筆者將安裝在輪蓋側(cè)的部分葉高葉片引入串列擴(kuò)壓器前排,相關(guān)公開發(fā)表的研究表明,裝在輪蓋側(cè)的半葉高擴(kuò)壓器性能稍優(yōu)于其他形式的半葉高擴(kuò)壓器[8],因此筆者將研究安裝在輪蓋側(cè)的半葉高擴(kuò)壓器,通過數(shù)值模擬分析部分葉高對串列擴(kuò)壓器性能的影響以及影響機(jī)理.

1 部分葉高串列擴(kuò)壓器模型

本文的離心葉輪型線來源于發(fā)動機(jī)及渦輪機(jī)聯(lián)盟弗里的希哈芬股份有限公司(MTU)設(shè)計(jì)的開式后彎離心葉輪,其設(shè)計(jì)轉(zhuǎn)速為35 200 r/min,葉片數(shù)為15;擴(kuò)壓器則是Rothstein[8]設(shè)計(jì)完成的通道型擴(kuò)壓器,葉片數(shù)為23,該離心壓氣機(jī)設(shè)計(jì)質(zhì)量流量qm,d=2.405 kg/s.宋文杰[9]首先利用數(shù)值模擬獲得了該離心葉輪出口氣流參數(shù)(如出口氣流角、馬赫數(shù)等)的分布,然后根據(jù)串列葉柵設(shè)計(jì)方法,參考NACA平面葉柵實(shí)驗(yàn)數(shù)據(jù)選擇NACA 65-(4A10)06作為葉型,為該離心葉輪設(shè)計(jì)了串列擴(kuò)壓器.表1給出了r3/r2=1.14條件下該串列擴(kuò)壓器前后排葉柵參數(shù),r2、r3分別指離心葉輪出口和串列擴(kuò)壓器葉片進(jìn)口半徑(見圖1).

表1 串列葉柵幾何參數(shù)(r3/r2=1.14)

圖1 部分葉高串列擴(kuò)壓器示意圖

前后排葉片周向位置示意圖如圖2所示.為了獲得串列擴(kuò)壓器前后排葉片的最佳相對位置,首先令前后排葉片的相對周向位置LRCP分別為10%、20%、…、80%和90%,從而獲得9種串列擴(kuò)壓器模型,通過數(shù)值方法分析相對周向位置對離心壓氣機(jī)性能的影響,獲得性能最佳的串列擴(kuò)壓器.然后將部分葉高設(shè)計(jì)引入串列擴(kuò)壓器前排葉片,令串列擴(kuò)壓器前排葉片高度h與串列擴(kuò)壓器高度B的比值分別為20%、30%、…、90%和100%,形成不同的部分葉高串列擴(kuò)壓器,如圖1所示,h/B為100%時對應(yīng)采用全葉高串列擴(kuò)壓器的原離心壓氣機(jī).最后,通過數(shù)值模擬方法分析h/B對離心壓氣機(jī)性能及擴(kuò)壓器內(nèi)流動的影響.

圖2 串列擴(kuò)壓器幾何示意圖

2 數(shù)值模擬方法與驗(yàn)證

2.1 流動控制方程及邊界條件

數(shù)值模擬基于計(jì)算流體力學(xué)商用軟件CFX完成.流動控制方程為N-S方程,湍流場通過SST湍流模型進(jìn)行模擬,近壁面使用自動壁面處理函數(shù).控制方程組采用基于有限元的有限體積法進(jìn)行離散.

為了減小計(jì)算量,計(jì)算時僅考慮了旋轉(zhuǎn)葉輪域和擴(kuò)壓器域,離心葉輪的進(jìn)口和擴(kuò)壓器出口有適當(dāng)?shù)难娱L段,計(jì)算域包含2個離心葉輪流道和3個擴(kuò)壓器流道,如圖3(a)所示.圖3(b)和圖3(c)分別為單葉片流道和單擴(kuò)壓器流道圖.計(jì)算域采用結(jié)構(gòu)化網(wǎng)格進(jìn)行離散,為了保證近壁區(qū)的模擬精度,預(yù)設(shè)第一層壁面網(wǎng)格的y+等于1,計(jì)算域總網(wǎng)格數(shù)為102萬.數(shù)值模擬涉及進(jìn)口、出口、固體壁面、動靜交界面以及周期性5類邊界條件.離心葉輪(擴(kuò)壓器)流道兩側(cè)采用節(jié)點(diǎn)一一對應(yīng)的周期性邊界條件,葉輪域與擴(kuò)壓器域的動靜交界面則采用混合平面法處理,固體壁面給定無滑移條件,進(jìn)口指定進(jìn)口總溫、總壓(湍流強(qiáng)度設(shè)置為5%),出口給出質(zhì)量流量條件.計(jì)算時間步長為1/Ω,Ω為離心葉輪旋轉(zhuǎn)角速度.

當(dāng)離心壓氣機(jī)運(yùn)行接近喘振時,其內(nèi)部流動將出現(xiàn)明顯的非定常流動特征,基于定常N-S方程很難獲得收斂的數(shù)值解.因此,在近喘振點(diǎn),求解的控制方程為非定常N-S方程.另外,從離心壓氣機(jī)性能曲線可知,當(dāng)其運(yùn)行接近喘振工況時隨著質(zhì)量流量的減小壓比增長梯度非常大,此時再采用質(zhì)量流量出口邊界條件不易收斂,改為壓力出口邊界條件.其他邊界條件和數(shù)值離散格式與定常流場求解相似,只是時間步長從1/Ω改為5.839 55×10-6s,離心葉輪與擴(kuò)壓器的動靜交界面處理改為瞬態(tài)轉(zhuǎn)子靜子法.

(a) 計(jì)算域整體

(b) 離心葉輪流道

(c) 擴(kuò)壓器流道

2.2 數(shù)值模擬方法驗(yàn)證

圖4給出了80%離心葉輪設(shè)計(jì)轉(zhuǎn)速(28 541 r/min)下、r3/r2=1.14時,原離心壓氣機(jī)的總壓比π*及等熵效率η隨質(zhì)量流量qm的分布,圖中實(shí)驗(yàn)值來自文獻(xiàn)[10],且P1、P2、M和S1分別對應(yīng)近喘振點(diǎn)、稍遠(yuǎn)離喘振點(diǎn)、氣動設(shè)計(jì)點(diǎn)和近堵塞點(diǎn).從圖4可以看出,無論是總壓比還是等熵效率,模擬值與實(shí)驗(yàn)值基本吻合,其中最大誤差出現(xiàn)在P2工況,總壓比的相對誤差值為3%.圖5為P1工況下,離心葉輪出口絕對速度和子午速度云圖.從圖5(a)可以看出,數(shù)值方法獲得的絕對速度稍大于實(shí)驗(yàn)值,與實(shí)驗(yàn)值相比,高速區(qū)的范圍更接近于流道中間而不是偏向吸力面.從圖5(b)可以看出,在流動核心區(qū),模擬值與實(shí)驗(yàn)值分布基本接近.但近輪轂側(cè)的模擬值的低速范圍要大于實(shí)驗(yàn)值.圖4和圖5表明,所采用的方法在離心壓氣機(jī)級的性能及內(nèi)部流場模擬方面均具有較好的精度,說明本文數(shù)值模擬方法是可靠的.

(a) 總壓比

(b) 等熵效率

(a) 葉輪出口絕對速度

(b) 葉輪出口子午速度

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

對離心壓氣機(jī)而言,最重要的性能參數(shù)包括總壓比、等熵效率、喘振裕度以及穩(wěn)定工作范圍,因此首先分析前后排葉片相對周向位置的變化對離心壓氣機(jī)性能參數(shù)的影響.圖6給出了設(shè)計(jì)轉(zhuǎn)速下,不同相對周向位置串列擴(kuò)壓器的離心壓氣機(jī)總壓比和等熵效率曲線.由圖6可知,當(dāng)LRCP增大時,離心壓氣機(jī)級的總壓比逐漸增大,等熵效率逐漸升高,在LRCP=30%時達(dá)到最大值,隨著LRCP繼續(xù)增大,總壓比逐漸減小,等熵效率降低.不過LRCP從10%增加到30%,級的總壓比和等熵效率的增長幅度相對較小,當(dāng)質(zhì)量流量為2.36 kg/s時,總壓比增大了0.6%,等熵效率提高了0.44%.

(a) 總壓比

(b) 等熵效率

圖6 串列擴(kuò)壓器相對周向位置對離心壓氣機(jī)總壓比和等熵效率的影響

Fig.6 Influence of RCP on the total pressure ratio and isentropic efficiency

圖7給出了喘振裕度隨LRCP的變化趨勢,圖中縱坐標(biāo)為離心壓氣機(jī)的喘振裕度,其定義為

(1)

圖7 串列擴(kuò)壓器相對周向位置對離心壓氣機(jī)喘振裕度的影響

可見,當(dāng)LRCP=10%、20%和30%時,離心壓氣機(jī)喘振裕度相對較大,分別為7.04%、7.28%和7.15%;當(dāng)LRCP=40%和50%時,離心壓氣機(jī)喘振裕度迅速減小.這表明串列擴(kuò)壓器前后排葉片的相對周向位置對離心壓氣機(jī)的穩(wěn)定工作范圍有較大影響,LRCP在20%~30%內(nèi)較為有利.

筆者在LRCP=30%的串列擴(kuò)壓器基礎(chǔ)上,將擴(kuò)壓器前排葉片設(shè)計(jì)為部分葉高形式,并對離心壓氣機(jī)流場與性能進(jìn)行分析.圖8為設(shè)計(jì)轉(zhuǎn)速下、前排葉片相對高度h/B從20%增加到90%時,離心壓氣機(jī)的總壓比和等熵效率.由圖8可知,雖然h/B由20%增加到90%時,離心壓氣機(jī)的總壓比和等熵效率持續(xù)小幅下降,這表明串列擴(kuò)壓器前排葉片采用部分葉高設(shè)計(jì)將不利于離心壓氣機(jī)總壓比的增大和等熵效率的提高.從圖8還可以看出,采用部分葉高串列擴(kuò)壓器的離心壓氣機(jī)的性能曲線變得更為平坦,這說明其穩(wěn)定工作范圍增加.圖9為離心壓氣機(jī)喘振裕度S隨前排葉片相對高度h/B的變化趨勢.從圖9可以看出,喘振裕度隨h/B的增加呈先增加后減小的趨勢,當(dāng)h/B增加到80%和90%時,離心壓氣機(jī)在工作質(zhì)量流量大于設(shè)計(jì)質(zhì)量流量時已發(fā)生了喘振.當(dāng)擴(kuò)壓器前排葉片高度等于40%和50%的擴(kuò)壓器高度時,喘振裕度為8.6%和8.9%,相比全葉高串列擴(kuò)壓器的喘振裕度(7.1%),分別增加了21%和25%;而由圖8可知,此時h/B為40%和50%的離心壓氣機(jī)的總壓比和等熵效率與全葉高時相比下降量都在0.9%~1.1%.這表明,在串列擴(kuò)壓器前排葉片采用部分葉高技術(shù)后,離心壓氣機(jī)的總壓比和等熵效率雖略有下降,但是合適的部分葉高串列擴(kuò)壓器能明顯增加離心壓氣機(jī)的喘振裕度和穩(wěn)定工作范圍.

(a) 總壓比

(b) 等熵效率

圖9 前排葉片相對高度對離心壓氣機(jī)喘振裕度的影響

圖10給出了額定工況下具有不同前排葉片相對高度的串列擴(kuò)壓器內(nèi)不同流向位置處的渦量云圖.圖10中,首先分析前排葉片所在流道部分,發(fā)現(xiàn)前排葉片尾緣處的尾跡范圍隨著前排葉片相對高度的增加而增加,且與h/B=100%時相比,由于葉根存在間隙,所以壓力面?zhèn)纫约叭~根的低能流體會在壓力面與吸力面?zhèn)鹊膲翰钭饔孟孪蛭γ鎮(zhèn)葯M向流動,在葉根間隙形成泄漏渦.從圖10還可以看出,沿著流向,擴(kuò)壓器前排葉片輪轂側(cè)上的渦量較大區(qū)域不斷增加,且?guī)缀跽紦?jù)了整個流道的寬度,而當(dāng)h/B=100%時,此渦量較大區(qū)域只存在于前排葉片壓力面?zhèn)雀浇γ鎮(zhèn)葞缀醪淮嬖?由于受到前排葉片的影響,后排葉片吸力面?zhèn)鹊牧鞯烂黠@存在前排葉片的尾跡,且尾跡區(qū)域隨著前排葉片相對葉高的增加而增大.

圖10 氣動設(shè)計(jì)點(diǎn)下擴(kuò)壓器內(nèi)的渦量云圖

圖11給出了前后排葉片流道的總壓損失系數(shù)ω隨h/B的變化曲線,圖中FVP表示前排葉片,AVP表示后排葉片.由圖11可知,隨著前排葉片相對高度h/B從10%逐漸增加到90%,前排葉片所在流道的總壓損失一直呈上升趨勢.這種分布特點(diǎn)與圖10給出的渦量分布云圖吻合,也是采用部分葉高串列擴(kuò)壓器后離心壓氣機(jī)等熵效率與總壓比略降的原因.圖11表明后排葉片流道的總壓損失系數(shù)隨h/B的變化趨勢與前排葉片所在流道的總壓損失系數(shù)變化趨勢大致相同.

圖11 部分葉高串列擴(kuò)壓器內(nèi)總壓損失系數(shù)

4 結(jié) 論

(1) 串列葉柵相對周向位置對擴(kuò)壓器的擴(kuò)壓能力以及穩(wěn)定工作范圍有很大影響,前后排葉片相對周向位置LRCP在20%~30%內(nèi)較為有利.

(2) 在串列葉柵相對周向位置為30%的基礎(chǔ)上,對前排葉片使用部分葉高串列擴(kuò)壓器,并通過數(shù)值模擬方法分析了不同前排葉片相對葉高條件下的離心壓氣機(jī)性能.數(shù)值結(jié)果表明,在串列擴(kuò)壓器前排葉片采用部分葉高技術(shù)后,當(dāng)h/B為40%和50%時,離心壓氣機(jī)的總壓比和等熵效率雖略有下降,但是喘振裕度和穩(wěn)定工作范圍明顯增加.

[1] PAMPREEN R C.The use of cascade technology in centrifugal compressor vaned diffuser design[J].Journal of Engineering for Power,1972,94(3):187-192.

[2] JAPIKSE D.Centrifugal compressor design and performance[M].Vermont,USA:Concepts ETI Inc.,1996.

[3] SELEZNEV K P,GALERKIN I B.Centrifugal compressors[M].Leningred,Russia:Mashinostroenie,1982.

[4] 李紹斌,王松濤,馮國泰,等.串列葉柵后排靜葉周向位置對壓氣機(jī)性能影響的數(shù)值研究[J].工程熱物理學(xué)報,2004,25(6):943-945.

LI Shaobin,WANG Songtao,FENG Guotai,etal.Numerical investigation on the compressor performances under different circumferential positions of the rear blades of tandem stator[J].Journal of Engineering Thermophysics,2004,25(6):943-945.

[5] 吳東坡,劉長勝.離心壓縮機(jī)級半開式葉輪采用無葉和串列擴(kuò)壓器的性能分析[J].風(fēng)機(jī)技術(shù),2009(3):8-10,14.

WU Dongpo,LIU Changsheng.Performance analy-sis of vaneless and tandem diffuser on unshrouded impeller of centrifugal compressor stage[J].Compressor,Blower & Fan Technology,2009(3):8-10,14.

[6] 朱營康,楊彥.半高導(dǎo)葉擴(kuò)壓器的研究[J].流體機(jī)械,1992,20(8):4-10.

ZHU Yingkang,YANG Yan.Study on the half vane diffuser[J].Fluid Machinery,1992,20(8):4-10.

[7] YOSHINAGA Y,KANEKI T,KOBAYASHI H,etal.A study of performance improvement for high specific speed centrifugal compressors by using diffusers with half guide vanes[J].Journal of Fluids Engineering,1987,109(4):359-366.

[8] ROTHSTEIN M.Entwicklung und experimentelle untersuchung eines radialverdichters mit beschaufeltem diffusor untersuchung eines radialverdichters mit beschaufeltem diffusor variabler geometrie[D].Aachen,Germany:RWTH Aachen,1993.

[9] 宋文杰.基于HALF GUIDE VANES串列葉柵擴(kuò)壓器設(shè)計(jì)[J].能源研究與信息,2015,31(4):234-241.

SONG Wenjie.Influence of relative circumferential position of tandem cascade on the performance of centrifugal compressor[J].Energy Research and Information,2015,31(4):234-241.

[10] ZIEGLER K U,GALLUS H E,NIEHUIS R.A study on impeller-diffuser interaction—part I:influence on the performance[J].Journal of Turbomachinery,2003,125(1):173-182.

Effect of Partial-height Tandem Diffuser on Performance of the Centrifugal Compressor

HUANGYueqing1,2,YANGAiling1,2,SONGWenjie1,2,CHENEryun1,2,DAIRen1,2

(1.School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China; 2.Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering,Shanghai 200093,China)

To ensure the pressure ratio of a centrifugal compressor and simultaneously broaden its stable working range,the concept of partial-height blade was introduced to the design of tandem diffuser,and subsequently the effects of relative circumference position (RCP) and relative height of forward blade (RHFB) on the performance of the centrifugal compressor and tandem diffuser were researched by numerical simulation.Based on the tandem diffuser with 30% RCP,different partial-height tandem diffusers were designed with partial-height blades used on the front row.Results show that RCP significantly affects the diffusion ability and stable working range of the diffuser,and the best performance appears as RCP ranges from 20% to 30%.Compared to the whole blade height,the surge margin of centrifugal compressor with RHFB lying in 40%-50% may be increased by 21%-25%,when the total pressure ratio and isentropic efficiency would be reduced only by 1%.

centrifugal compressor; tandem diffuser; relative circumference position; partial-height blade; numerical simulation

2015-11-13

2016-03-26

上海市科委基地建設(shè)資助項(xiàng)目(13DZ2260900)

黃月晴(1991-),女,安徽宣城人,碩士研究生,研究方向?yàn)閯恿C(jī)械.電話(Tel.):18110351223;E-mail:huangyueqing1991@outlook.com.

1674-7607(2017)01-0026-07

TH452

A 學(xué)科分類號:470.30

猜你喜歡
效率
你在咖啡館學(xué)習(xí)會更有創(chuàng)意和效率嗎?
提升朗讀教學(xué)效率的幾點(diǎn)思考
甘肅教育(2020年14期)2020-09-11 07:57:42
注意實(shí)驗(yàn)拓展,提高復(fù)習(xí)效率
效率的價值
商周刊(2017年9期)2017-08-22 02:57:49
引入“倒逼機(jī)制”提高治霾效率
質(zhì)量與效率的爭論
跟蹤導(dǎo)練(一)2
提高食品行業(yè)清潔操作的效率
OptiMOSTM 300V提高硬開關(guān)應(yīng)用的效率,支持新型設(shè)計(jì)
“錢”、“事”脫節(jié)效率低
主站蜘蛛池模板: 欧美在线一级片| 国产第一页屁屁影院| 国产福利免费在线观看| 国产精品视频导航| 中文字幕天无码久久精品视频免费| 99精品视频在线观看免费播放| 无码一区中文字幕| 日韩毛片基地| 久久青草视频| 色偷偷综合网| 九色国产在线| 欧美日韩动态图| 日本午夜影院| 亚洲毛片一级带毛片基地| 青草视频网站在线观看| 五月丁香伊人啪啪手机免费观看| 精品国产免费人成在线观看| 国产欧美日韩va| 亚洲天堂免费在线视频| 亚洲av日韩av制服丝袜| 试看120秒男女啪啪免费| …亚洲 欧洲 另类 春色| 亚洲色图欧美激情| 精品撒尿视频一区二区三区| 国产成人1024精品下载| 黄色福利在线| 97se亚洲综合在线天天| 日本免费精品| a级毛片一区二区免费视频| 亚洲三级成人| 精品1区2区3区| 最近最新中文字幕在线第一页| 91视频区| 99久久性生片| 亚洲最大在线观看| AV网站中文| 久久免费视频6| 欧美在线导航| 亚洲另类色| www.99在线观看| 亚洲天堂日韩av电影| 99人妻碰碰碰久久久久禁片| 亚洲黄色片免费看| 欧美啪啪网| 丰满人妻久久中文字幕| 综合社区亚洲熟妇p| 在线观看av永久| 久久一级电影| 国产91色| 99久视频| 欧美伊人色综合久久天天| 中文一区二区视频| 久久久亚洲色| 国产二级毛片| 怡春院欧美一区二区三区免费| 少妇露出福利视频| 久久天天躁狠狠躁夜夜躁| 欧美一区二区三区欧美日韩亚洲 | 69av在线| 黄色网站在线观看无码| 女人18一级毛片免费观看| 精品国产一区二区三区在线观看| 成人免费视频一区二区三区 | 香蕉综合在线视频91| 国产白浆视频| 国产麻豆精品久久一二三| 久久婷婷色综合老司机| 精品成人一区二区| 四虎精品黑人视频| 免费看久久精品99| 成人毛片免费在线观看| 免费人成在线观看成人片| 色噜噜在线观看| 激情综合五月网| 成人国产精品2021| 中文字幕丝袜一区二区| 五月婷婷综合网| 国产一级片网址| 福利国产微拍广场一区视频在线| 国产精品不卡永久免费| 波多野结衣亚洲一区| 欧美不卡视频一区发布|