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

    軸流壓氣機轉(zhuǎn)子端壁造型多目標(biāo)優(yōu)化的數(shù)值研究

    2017-04-17 01:38:40晁曉亮賈一哲
    流體機械 2017年3期
    關(guān)鍵詞:總壓壓氣機輪轂

    高 翔,晁曉亮,賈一哲

    (中國飛行試驗研究院,陜西西安 710089)

    軸流壓氣機轉(zhuǎn)子端壁造型多目標(biāo)優(yōu)化的數(shù)值研究

    高 翔,晁曉亮,賈一哲

    (中國飛行試驗研究院,陜西西安 710089)

    針對NASA跨音速轉(zhuǎn)子67,建立了可靠的數(shù)值計算模型,對3種不同的目標(biāo)函數(shù)開展了非軸對稱端壁造型的優(yōu)化工作。3種優(yōu)化方案均使得轉(zhuǎn)子性能得到了不同程度改善,效率提升最大為0.41%。優(yōu)化結(jié)果表明,由于相同流量的約束,轉(zhuǎn)子通道中流量分布沿徑向發(fā)生了改變,對于近輪轂的端壁區(qū)域來說,造型后靜壓梯度的方向?qū)α鲌鼍哂兄匾绊?。不同造型結(jié)果中葉片尾緣三維角區(qū)分離的位置與范圍也有所不同,當(dāng)三維角區(qū)分離后移且沿徑向發(fā)生遷移時,會使得效率得到的增益明顯降低。

    跨音速軸流壓氣機轉(zhuǎn)子;非軸對稱端壁造型;二次流;多目標(biāo)優(yōu)化

    author_info_translate: Reliable numerical computation method was built on rotor67.Based on this numerical computation method,optimization of non-axisymmetric endwall contouring focusing on three different objective functions were conducted.The three optimal results all improved the performances of the rotor and the maximum improvement was 0.41%.Because of the restriction of the mass flow rate during optimization,the flow fields in different span height changed to different extent.As for the regions near hub,validity of the gradient direction influences the flow field significantly.In different optimal results,the location and scope of the three-dimensional corner separation near the trailing edge of the blade altered obviously,and the gain in efficiency drops when the separation region moves backward and migrates towards spanwise direction.

    1 前言

    非軸對稱端壁造型是對葉輪機械的端壁引入凹凸變化結(jié)構(gòu)來達到調(diào)整局部流動的目的。近年來,非軸對稱造型在渦輪中的應(yīng)用已經(jīng)較為成熟[1~4],但是在壓氣機中,由于級負荷比渦輪小,以及其獨特的逆壓梯度環(huán)境,使得非軸對稱端壁造型的應(yīng)用比渦輪復(fù)雜的多。

    目前,國內(nèi)外針對壓氣機展開的非軸對稱端壁造型研究,主要是在壓氣機葉柵上展開,而直接針對壓氣機轉(zhuǎn)子進行的非軸對稱端壁造型優(yōu)化較少[5~7]。因此,本文以跨音速壓氣機轉(zhuǎn)子Rotor67為研究對象,采用多目標(biāo)優(yōu)化的方法,對其轉(zhuǎn)子輪轂進行非軸對稱端壁造型優(yōu)化,通過優(yōu)化結(jié)果的對比分析,闡述非軸對稱端壁造型在壓氣機轉(zhuǎn)子中的直接應(yīng)用與葉柵中應(yīng)用的不同點,同時具體分析其流動控制機理,進而對該造型技術(shù)在壓氣機中的實際應(yīng)用提供參考。

    2 研究對象

    本文采用NASA的Rotor67作為研究對象,其詳細的試驗測量數(shù)據(jù)已在NASA報告中進行了公布[8~12]。Rotor67的基本設(shè)計參數(shù)見表1。

    表1 Rotor67設(shè)計參數(shù)

    3 數(shù)值模型建立與驗證

    采用NUMECA軟件進行網(wǎng)格塊的繪制,數(shù)值計算采用Jameson的有限體積差分格式,結(jié)合Spalart-Allmaras湍流模型求解相對坐標(biāo)系下的三維雷諾平均Navier-Stokes方程,采用四階Runge-Kutta法時間推進獲得定常解。

    圖1給出了數(shù)值計算結(jié)果與試驗結(jié)果的對比,可以看出效率和壓比特性線整體趨勢與試驗值均吻合較好,壓比計算結(jié)果相對試驗值有一定程度偏低,但趨勢基本一致,可以作為后續(xù)計算的基礎(chǔ)。

    (a) 效率特性

    (b) 壓比特性

    4 優(yōu)化參數(shù)設(shè)置

    在本文的參數(shù)化過程中,Rotor67轉(zhuǎn)子輪轂曲面采用NURBS曲面生成的非軸對稱端壁造型,如圖2所示,該曲面由平行于葉片中弧線的5條貝塞爾曲線構(gòu)成,每條控制曲線均勻分布著4個控制點,造型區(qū)域從距前緣點10%軸向弦長開始,至距尾緣點10%軸向弦長結(jié)束,即整個造型型面完全處于轉(zhuǎn)子通道內(nèi),因此所有20個控制點均可以在給定范圍內(nèi)自由擾動,同時可以滿足造型曲面與上下游原始端壁曲面的一階導(dǎo)數(shù)連續(xù)性,20個控制點可沿徑向在±5 mm(3.5%平均葉高)的范圍內(nèi)自由變化。優(yōu)化工況為Rotor67轉(zhuǎn)子的設(shè)計點(m=33.25 kg/s),所有邊界條件均與原始計算設(shè)置相同。針對20個自由變量,最終通過數(shù)據(jù)庫生成了70個樣文用于后續(xù)的優(yōu)化工作。

    圖2 非軸對稱端壁造型示意

    優(yōu)化的過程中,除了將壓氣機的總體性能指標(biāo)——效率與總壓比作為優(yōu)化目標(biāo)參數(shù)外,還選取了圖2中距轉(zhuǎn)子尾緣10%軸向弦長處平面的二次流動能大小作為優(yōu)化目標(biāo)參數(shù)。本文采用以下方法來定義二次流速度:

    (1)選定二次流優(yōu)化平面。二次流在通道中充分發(fā)展累積,在通道出口平面附近往往會達到最大值,因此本文選取葉片尾緣點上游10%軸向弦長處的平面作為二次流優(yōu)化的目標(biāo)平面。

    (2)在圓柱坐標(biāo)系下,通過周向平均的方法,求出選定的平面上每一個徑向網(wǎng)格上的速度Vz(r)、Vr(r)、Vθ(r),則此時每層徑向網(wǎng)格上的主流速度方向可定義為:

    (1)

    (2)

    (3)

    (5)二次流動能SKE

    (4)

    為了保證造型過程中流道面積變化不至于過大,將喉道區(qū)域的面積作為一個優(yōu)化目標(biāo)參數(shù),使所得優(yōu)化轉(zhuǎn)子造型的喉道面積盡可能的與原始轉(zhuǎn)子喉道面積接近。綜上所述,目標(biāo)函數(shù)包含了4個部分:效率、壓比、出口平面SKE以及喉道面積,即:

    OF=QW1·Pη+QW2·Ppr+QW3·PSKE+QW4·PThroat Area

    (5)

    式中OF——目標(biāo)函數(shù)QWi——不同目標(biāo)參數(shù)給定的權(quán)重系數(shù),通過改變其值的大小,來實現(xiàn)不同的優(yōu)化方案,i=1,2,3,4

    采用了3種方案進行優(yōu)化,分別以效率、二次流動能、效率及二次流動能為主要優(yōu)化目標(biāo)。設(shè)置目標(biāo)函數(shù)時,總壓比采取了限定最小值的方式進行優(yōu)化,即當(dāng)優(yōu)化產(chǎn)生的樣本總壓比低于原始轉(zhuǎn)子設(shè)計工況的計算壓比1.6587時,目標(biāo)函數(shù)罰值會迅速增加至很大,使得該樣本在進一步的遺傳算法學(xué)習(xí)中被淘汰。因此,在有效樣本中,對目標(biāo)函數(shù)具有貢獻的只有效率、SKE以及喉道面積3個參數(shù)。表2~4分別給出這3個參數(shù)的目標(biāo)值、權(quán)重系數(shù)以及其罰值對目標(biāo)函數(shù)貢獻的比例。

    表2 效率為主要優(yōu)化目標(biāo)

    表3 二次流動能為主要優(yōu)化目標(biāo)

    表4 效率與二次流動能為主要優(yōu)化目標(biāo)

    5 優(yōu)化結(jié)果分析

    5.1 總體性能分析

    將原始Rotor67轉(zhuǎn)子記為origin,3種方案最終優(yōu)化得出的結(jié)果分別記為opt-eff、opt-ske、opt-eff&ske,則設(shè)計工況下(m=33.25 kg/s)的相關(guān)參數(shù)優(yōu)化結(jié)果如表5所示??梢钥闯觯?種優(yōu)化結(jié)果均獲得了一定程度的效率提升,但是出口平面二次流動能并非全部降低,對于獲得最大效率提升的opt-eff來說,其出口二次流動能反而有所增加,其原因?qū)⒃诤笪姆治鲋羞M一步給出。

    表5 目標(biāo)參數(shù)優(yōu)化結(jié)果

    圖3給出了3種優(yōu)化方案所得到的輪轂半徑變化云圖。opt-eff在通道內(nèi)靠近葉片壓力面的區(qū)域形成了一定程度的上凸,且在通道前半部分區(qū)域的上凸幅值較大,最大幅值約為4 mm,在通道尾部,上凸結(jié)構(gòu)從壓力面向通道中部有所延伸但并未達到葉片吸力面;在30%~50%軸向弦長區(qū)域內(nèi),靠近葉片吸力面的位置處形成了較為明顯的下凹結(jié)構(gòu),下凹最大幅值約為3 mm。opt-ske其半徑變化云圖十分明顯的一個特點就是在轉(zhuǎn)子通道中間出現(xiàn)了大面積的下凹變化,僅僅在50%~60%軸向弦長區(qū)域內(nèi)的靠近葉片吸力面位置出現(xiàn)了幅值約為2 mm的上凸變化。對于opt-eff&ske,其變化幅值比前2種方案小,輪轂面的幾何光順程度最高,在30%~80%軸向弦長區(qū)域靠近葉片壓力面的位置形成了最大幅值約為3 mm的下凹,在葉片通道中部形成了幅值僅約為1 mm的微小上凸。

    (a) opt-eff

    (b) opt-ske

    (c) opteff&ske

    5.2 流場結(jié)果分析

    非軸對稱端壁造型的最初目的是改善近端壁區(qū)域的流動狀況,減弱端區(qū)二次流流動,從而帶來損失的減小與效率的提升。圖4給出了轉(zhuǎn)子出口平面(NASA報告中Station2位置)的相對總壓損失系數(shù)Cp的徑向分布。Cp定義為:

    (6)

    式中Pt-rel-inlet——進口的相對總壓Pt-rel——出口的當(dāng)?shù)叵鄬倝害裪nlet——進口的密度winlet——進口的相對速度

    (a) 0~30%葉高

    (b) 30%~100%葉高

    為了更加清晰地進行比較,圖4中將整個葉展區(qū)域的相對總壓損失系數(shù)分為了兩部分。從圖4(a)可以看出,在5%葉高范圍內(nèi),opt-eff的相對總壓損失系數(shù)比原始轉(zhuǎn)子有所增加,盡管這與非軸對稱端壁造型改善近端壁區(qū)域流動的目的相違背。從圖4(b)中可以發(fā)現(xiàn),在50%以上的葉高部分,opt-eff的相對總壓損失系數(shù)比原始轉(zhuǎn)子降低,兩部分綜合作用的結(jié)果使整個轉(zhuǎn)子通道內(nèi)的相對總壓損失減小,達到了提升效率的目的;而對于opt-ske來說,其情況與opt-eff剛好相反,在近端壁區(qū)域的5%葉高內(nèi),相對總壓損失得到了有效的降低,而50%以上葉高相對總壓損失有一定程度增加,但是增加的幅值并沒有opt-eff中減小的幅值那么大;對于opt-eff&ske來說,其相對總壓損失系數(shù)減小的區(qū)域集中在3%~15%葉高范圍內(nèi),在特別靠近端壁的區(qū)域以及15%以上的葉高部分,其總壓損失系數(shù)幾乎沒有變化。從以上對損失的分析中可以發(fā)現(xiàn),對壓氣機轉(zhuǎn)子的非軸對稱端壁造型,與葉柵端壁造型明顯不同,由于優(yōu)化過程中相同設(shè)計流量(m=33.25 kg/s)的限制,使得流量在葉高范圍內(nèi)重新分布,最終影響范圍延伸至整個葉片通道。

    雖然非軸對稱端壁造型所帶來的是全葉高范圍內(nèi)的作用,但是近內(nèi)輪轂面的低葉展區(qū)域仍然是受到造型影響的最直接區(qū)域,因此后續(xù)著重分析了近轉(zhuǎn)子輪轂面的低葉展區(qū)域的流場變化。

    圖5給出了轉(zhuǎn)子輪轂面的極限流線以及靜壓分布。從圖中可以看出,非軸對稱端壁造型引起輪轂面的凹凸變化,其所帶來的最直接的影響就是輪轂面靜壓分布的改變。對于opt-eff,葉片壓力面前緣以及尾緣的高靜壓區(qū)域的范圍均有明顯增大;對于opt-ske,原始轉(zhuǎn)子在葉片壓力面前緣與尾緣的小面積高靜壓區(qū)域均消失,取代的是在葉片壓力面的中部位置出現(xiàn)了很大范圍的高靜壓區(qū)域,這與圖3中該位置很大范圍的下凹型面直接相關(guān);對于opt-eff&ske,由于葉片吸力面附近的凹凸變化不大,其吸力面部分靜壓幾乎沒有改變,而在壓力面25%、75%兩個位置分別形成了兩小塊對稱的高靜壓區(qū)域。

    (a) 原型轉(zhuǎn)子 (b) opt-eff

    (c) opt-ske

    (d) opt-eff&ske

    圖5 輪轂表面極限流線及靜壓分布

    圖5中箭頭標(biāo)出了轉(zhuǎn)子通道中不同區(qū)域的靜壓梯度方向,這些靜壓梯度的方向所代表的就是通道壓差對氣流驅(qū)動力的有效方向。在非軸對稱端壁造型中,并非通道橫向壓差的改變越大,對流場的作用效果就越強,只有當(dāng)壓差的有效作用方向與原始流線接近正交時,通過壓差作用才能更有效地影響到流動的偏轉(zhuǎn)方向。與原始轉(zhuǎn)子相比,opt-eff與opt-ske在尾緣的有效壓差方向都有很大改變,產(chǎn)生了較強的逆壓梯度。除此之外,opt-ske在前半部分的有效壓差方向也有所改變,由原來的“順流向”壓差變?yōu)榱恕澳媪飨颉眽翰睢?/p>

    圖5中各圖的流線出發(fā)位置相同,比較同位置的流線可以發(fā)現(xiàn),對于opt-ske與opt-eff&ske來說,通道橫向流動達到相鄰葉片吸力面的位置均有所推后,即延緩了通道中的橫向流動的發(fā)展,而這一作用也使得葉片尾緣的分離渦更靠近葉片吸力面,如圖中C、D位置所示。而對于opt-eff,其通道中橫向流動到達相鄰葉片吸力面的位置略有所提前,尾緣分離渦位置稍遠離葉片吸力面,但是效果不是十分明顯。

    圖6給出了葉片吸力面的極限流線以及相對馬赫數(shù)分布。可以看出opt-ske與opt-eff&ske中端壁二次流與葉片吸力面的交匯位置明顯后移,且尾緣的分離區(qū)變得更加“細長”;而opt-eff中端壁二次流與吸力面交匯位置稍有提前,但不是十分明顯,這都與前述圖5中的分析相一致。值得注意的是在端壁二次流與吸力面交匯位置前,對于opt-eff以及opt-eff&ske來說,近端壁區(qū)域的相對馬赫數(shù)均有了明顯增加,即近端壁區(qū)域的通流能力得到了一定程度改善,如圖中A、C位置所示;而在opt-ske中,雖然交匯位置前的相對馬赫數(shù)也有一定程度提高,但在前緣20%的位置處則出現(xiàn)了一小塊面積不大的堵塞區(qū)域,如圖中B位置所示,這意味著在通道內(nèi)部出現(xiàn)了一定范圍的堵塞加重,損失增加的情況,進而使得全局效率的提升不是十分明顯。

    (a) 原型轉(zhuǎn)子

    (b) opt-eff

    (c) opt-ske

    (d) opt-eff&ske

    5.3 非設(shè)計工況分析

    雖然非軸對稱端壁造型優(yōu)化是對Rotor67轉(zhuǎn)子的設(shè)計工況(m=33.25 kg/s)進行的,但是輪轂型面的改變必然也會引起到其它工況點流場的變化,而其它工況點性能的改變,轉(zhuǎn)子工作裕度的改變等對整臺壓氣機也具有較為重要的影響。因此,在得到3種優(yōu)化結(jié)果opt-eff、opt-ske及opt-eff&ske后,對其分別進行了全工況的數(shù)值模擬,得到全工況范圍內(nèi)的性能特性線如圖7所示,需要指出的是,圖中所得曲線在近失速點逼喘時,采取的最小背壓提升間隔為100 Pa。

    (a) 效率特性

    (b) 壓比特性

    從圖7中可以看出,雖然在優(yōu)化過程中限制了轉(zhuǎn)子喉道面積的變化,但是對于opt-ske來說,堵塞點流量還是發(fā)生了較為明顯的改變(增加約1.0%),而opt-eff與opt-eff&ske的堵塞點流量基本沒有改變。仔細考慮這種現(xiàn)象產(chǎn)生的原因,可能是由于壓氣機轉(zhuǎn)子在設(shè)計時設(shè)計點流量的給定無法考慮到三維角區(qū)分離引起的堵塞。近輪轂區(qū)域的流場均為亞音速流場,其出口平面的流場變化可以傳播到流場上游,因此,雖然轉(zhuǎn)子流量在設(shè)計時是以喉道面積進行限制的,但是下游堵塞區(qū)域嚴重時對流量的限制比喉道的限制更加嚴格,此現(xiàn)象發(fā)生時相當(dāng)于喉道位置“后移”至堵塞最嚴重的區(qū)域,流量限制的關(guān)鍵因素變成了堵塞位置能通過的流量大小。opt-ske通過非軸對稱端壁造型,消除了堵塞工況近輪轂位置的堵塞,使得流通能力比與原始轉(zhuǎn)子增強,因此其堵塞工況的流量增加。而opt-eff與opt-eff&ske沒有消除堵塞工況時近輪轂區(qū)域的堵塞,因此它們的流量限制因素與原始轉(zhuǎn)子相同,均被該堵塞位置的通流能力限制。

    從以上分析可以看出,在非軸對稱端壁造型技術(shù)應(yīng)用于實際的壓氣機中時,需要將其對堵塞點流量的改變考慮進去,因為堵塞點流量的大小會關(guān)系到壓氣機在系統(tǒng)中工作時的匹配問題。

    6 結(jié)論

    (1)在設(shè)計工況,3種優(yōu)化方案的效率均得到了提升,其中opt-eff與opt-eff&ske的效率提升較為明顯,分別達到了0.41%與0.32%,而opt-ske的效率提升僅為0.03%。

    (2)在壓氣機轉(zhuǎn)子中采用非軸對稱端壁造型的初衷是改善近輪轂區(qū)域的流動,但由于優(yōu)化過程中同流量點的限制,非軸對稱端壁造型引起了轉(zhuǎn)子通道中徑向流量的重新分布,從而使得全葉高范圍內(nèi)的流場均受到了一定影響。

    (3)非軸對稱端壁造型通過輪轂面的凹凸變化直接改變其附近的靜壓分布,造型后靜壓梯度的有效方向十分關(guān)鍵,當(dāng)其與原始二次流流向接近垂直時,作用效果最為明顯。

    (4)非軸對稱端壁造型會引起葉片尾緣三維角區(qū)分離的徑向擴展,同時會使得轉(zhuǎn)子通道內(nèi)堵塞情況有所改變。

    (5) 將非軸對稱端壁造型技術(shù)應(yīng)用于實際的壓氣機中時,需要對堵塞點流量的改變進行考慮,以滿足發(fā)動機部件之間的匹配問題。

    [1] Harvey N W,Rose M G.Non-axisymmetric turbine endwall design :three2dimensional linear design system[C].Proceedings of ASME Turbo Expo 1999 Power for Land ,Sea & Air ,New York ,USA ,1999.

    [2] Hartland J C,Gregory-Smith D.Non-axisymmetric turbine end wall design:experimental validation[C].Proceedings of ASME Turbo Expo 1999 Power for Land,Sea& Air,New York ,USA ,1999.

    [3] Brennan G,Harvey N W,Rose M G,et al.Improving the efficiency of the trent 500 HP turbine using non-axisymmetric endwalls part I:turbine design[R].ASME Paper 2001GT-0444,2001.

    [4] 黃洪雁,王仲奇,馮國泰.上端壁翹曲對渦輪葉柵流場的影響[J].推進技術(shù),2002,23(1):36-39.

    [5] Hoeger M,Cardamone P,Fottner L.Influence of endwall contouring on the transonic flow in a compressor blade[R].ASME Paper GT2002-30440,2002.

    [6] 盧家玲,楚武利,吳艷輝.利用靜子非軸對稱輪轂造型提高軸流壓縮機效率[J].風(fēng)機技術(shù),2009(2):3-6.

    [7] Harvey N W.Some effects of non axisymmetric endwall profiling on axial flow compressor aero dynamics Part II:Multi-Stage HPC CFD Study[R].ASME Paper GT2008-50991,2008.

    [8] 鄭覃,盧少鵬,羌曉青,等.動葉局部彎曲對壓氣機流場和性能的影響[J].流體機械,2015,43(7):27-31.

    [9] 劉力歌,馮文卓,高歡,等.一種新型網(wǎng)殼儲罐的支座結(jié)構(gòu)及其優(yōu)化[J].壓力容器,2016,33(2):30-34.

    [10] 祝華云,徐志明,胡國軍,等.尾緣厚度對壓氣機葉柵氣動性能的影響研究[J].機電工程,2016,33(2):187-190.

    [11] 王晉華,同利平,趙珊珊.離心壓縮機組管道優(yōu)化布置及應(yīng)力分析計算[J].化工設(shè)備與管道,2014,51(4):50-53.

    [12] Strazisar A J.Laser Anemometer Measurement in a Transon Axial-Flow Fan Rotor[M].NASA-TP-2879,1989.

    Numerical Investigation on a Transonic Axia-flow Compressor Rotor with Multi-objective Endwall Profiling Optimization

    GAO Xiang,CHAO Xiao-liang,JIA Yi-zhe

    (Chinese Flight Test Establishment,Xi′an 710089,China)

    transonic axial-flow compressor rotor;non-axisymmetric endwall;secondary flow;multi-objective optimization

    1005-0329(2017)03-0026-07

    2016-05-19

    2016-06-15

    TH41;V231.3

    A

    10.3969/j.issn.1005-0329.2017.03.006

    高翔(1990-),男,碩士,主要從事航空發(fā)動機性能特性試飛、葉輪機械氣動熱力學(xué)等方面的研究,通訊地址:710000 陜西西安市閻良區(qū)凌云路5號 中國飛行試驗研究院發(fā)動機所,E-mail:397180893@qq.com。

    猜你喜歡
    總壓壓氣機輪轂
    總壓探針性能結(jié)構(gòu)敏感性分析
    軸流壓氣機效率評定方法
    抗壓痕透明粉在精車鋁輪轂上的應(yīng)用研究
    上海涂料(2021年5期)2022-01-15 06:09:26
    重型燃氣輪機壓氣機第一級轉(zhuǎn)子葉片斷裂分析
    可調(diào)式總壓耙設(shè)計及應(yīng)用
    壓氣機緊湊S形過渡段內(nèi)周向彎靜子性能數(shù)值計算
    亞聲速條件下總壓探針臨壁效應(yīng)的數(shù)值研究
    2 m超聲速風(fēng)洞流場變速壓控制方法研究
    基于CPS 的汽車輪轂制造系統(tǒng)設(shè)計
    基于ANSYS的輪轂支架結(jié)構(gòu)設(shè)計
    国产乱来视频区| 亚洲欧洲日产国产| 欧美人与性动交α欧美精品济南到| 看免费成人av毛片| 免费黄网站久久成人精品| netflix在线观看网站| 一区二区av电影网| 色94色欧美一区二区| 男女床上黄色一级片免费看| 在线亚洲精品国产二区图片欧美| 亚洲国产最新在线播放| 大码成人一级视频| 婷婷色麻豆天堂久久| 国产成人精品在线电影| 亚洲精品久久成人aⅴ小说| 老司机靠b影院| 在现免费观看毛片| 黑人欧美特级aaaaaa片| 亚洲国产av影院在线观看| 伊人亚洲综合成人网| 99re6热这里在线精品视频| 日本91视频免费播放| 亚洲欧美日韩另类电影网站| 日日爽夜夜爽网站| 国产精品三级大全| 欧美日韩福利视频一区二区| 欧美xxⅹ黑人| 精品少妇黑人巨大在线播放| 无遮挡黄片免费观看| 日韩一区二区三区影片| 在线观看免费高清a一片| 在线观看三级黄色| 久久久久精品久久久久真实原创| 老鸭窝网址在线观看| 久久久久久久久久久免费av| 国产午夜精品一二区理论片| 涩涩av久久男人的天堂| 丰满饥渴人妻一区二区三| av电影中文网址| 久久综合国产亚洲精品| 亚洲国产欧美在线一区| 看非洲黑人一级黄片| 国产精品 国内视频| 最新的欧美精品一区二区| 汤姆久久久久久久影院中文字幕| 丝袜人妻中文字幕| 99精国产麻豆久久婷婷| 久久久久久久久久久免费av| 国产精品久久久人人做人人爽| 日韩制服丝袜自拍偷拍| 国产精品一区二区精品视频观看| 国产成人a∨麻豆精品| av在线播放精品| 少妇被粗大猛烈的视频| 母亲3免费完整高清在线观看| 国产一区二区三区综合在线观看| 国产成人系列免费观看| 中文字幕制服av| 久久人人爽人人片av| 亚洲欧美清纯卡通| 国产99久久九九免费精品| 国产激情久久老熟女| 性少妇av在线| 国产精品一国产av| 男女边摸边吃奶| 欧美日韩综合久久久久久| 亚洲成人国产一区在线观看 | 91老司机精品| 自拍欧美九色日韩亚洲蝌蚪91| 观看av在线不卡| 午夜日本视频在线| 男女之事视频高清在线观看 | 国产精品嫩草影院av在线观看| 最近中文字幕高清免费大全6| www.熟女人妻精品国产| 99精品久久久久人妻精品| 亚洲一区二区三区欧美精品| 人体艺术视频欧美日本| a级片在线免费高清观看视频| 欧美av亚洲av综合av国产av | a级毛片在线看网站| 日韩熟女老妇一区二区性免费视频| 十八禁网站网址无遮挡| 女性生殖器流出的白浆| 叶爱在线成人免费视频播放| 中文天堂在线官网| 欧美xxⅹ黑人| a级毛片在线看网站| av有码第一页| 91精品三级在线观看| 亚洲精品久久成人aⅴ小说| 91成人精品电影| 悠悠久久av| 麻豆精品久久久久久蜜桃| 成年人午夜在线观看视频| 午夜福利视频精品| 亚洲国产精品一区三区| 亚洲国产中文字幕在线视频| 桃花免费在线播放| 国产精品秋霞免费鲁丝片| 婷婷色综合大香蕉| 如日韩欧美国产精品一区二区三区| 国产成人精品久久二区二区91 | 亚洲精品自拍成人| 大话2 男鬼变身卡| 成人亚洲精品一区在线观看| 欧美日本中文国产一区发布| 亚洲国产精品999| 两个人免费观看高清视频| 亚洲国产看品久久| 丁香六月天网| 999久久久国产精品视频| 妹子高潮喷水视频| 一区二区三区激情视频| 麻豆乱淫一区二区| 亚洲图色成人| 日日啪夜夜爽| 欧美97在线视频| 国产精品麻豆人妻色哟哟久久| 新久久久久国产一级毛片| 日韩制服丝袜自拍偷拍| 悠悠久久av| 黄色 视频免费看| 激情视频va一区二区三区| 久久久久国产精品人妻一区二区| 黄色视频在线播放观看不卡| 9191精品国产免费久久| 美女福利国产在线| 少妇猛男粗大的猛烈进出视频| 综合色丁香网| 久久影院123| 久久午夜综合久久蜜桃| 日韩av免费高清视频| 大码成人一级视频| 美女扒开内裤让男人捅视频| 亚洲精品乱久久久久久| svipshipincom国产片| 高清不卡的av网站| 久久久久网色| 欧美 日韩 精品 国产| 日韩一区二区视频免费看| 国产免费福利视频在线观看| 一区在线观看完整版| 爱豆传媒免费全集在线观看| 成人影院久久| 日本欧美国产在线视频| 一边亲一边摸免费视频| 亚洲欧美一区二区三区久久| 国产成人欧美| 女性生殖器流出的白浆| 国产av精品麻豆| 国产精品亚洲av一区麻豆 | 捣出白浆h1v1| 日韩人妻精品一区2区三区| 久久久国产精品麻豆| 成人黄色视频免费在线看| 亚洲一区二区三区欧美精品| 久久久亚洲精品成人影院| 制服诱惑二区| 18禁裸乳无遮挡动漫免费视频| 女人精品久久久久毛片| 一本色道久久久久久精品综合| xxx大片免费视频| 一边亲一边摸免费视频| 欧美国产精品一级二级三级| 中国国产av一级| 91老司机精品| 国产99久久九九免费精品| 国产成人a∨麻豆精品| 国语对白做爰xxxⅹ性视频网站| 亚洲四区av| 下体分泌物呈黄色| 久久人人爽人人片av| 999精品在线视频| 男人舔女人的私密视频| 亚洲精品一区蜜桃| 极品少妇高潮喷水抽搐| tube8黄色片| 久久久精品免费免费高清| 午夜福利,免费看| 捣出白浆h1v1| 黑丝袜美女国产一区| 最近最新中文字幕大全免费视频 | 亚洲成人av在线免费| 亚洲欧美一区二区三区黑人| 日韩一本色道免费dvd| 久久精品亚洲av国产电影网| 街头女战士在线观看网站| 久久久久精品久久久久真实原创| 七月丁香在线播放| 午夜激情久久久久久久| 丰满乱子伦码专区| 最近中文字幕高清免费大全6| 9191精品国产免费久久| 女性被躁到高潮视频| 一级黄片播放器| 成年av动漫网址| 日韩 欧美 亚洲 中文字幕| 亚洲欧美成人精品一区二区| 久久鲁丝午夜福利片| 夫妻性生交免费视频一级片| 久久精品国产亚洲av高清一级| 亚洲av国产av综合av卡| 国产精品人妻久久久影院| 美女福利国产在线| 亚洲精品在线美女| 日韩欧美一区视频在线观看| 色网站视频免费| 亚洲熟女精品中文字幕| 一级毛片 在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 免费av中文字幕在线| 亚洲情色 制服丝袜| 国精品久久久久久国模美| 精品卡一卡二卡四卡免费| 夜夜骑夜夜射夜夜干| 老司机深夜福利视频在线观看 | 丝袜美腿诱惑在线| 一级,二级,三级黄色视频| 看非洲黑人一级黄片| 久久久久国产一级毛片高清牌| 国产高清国产精品国产三级| 一二三四中文在线观看免费高清| 在现免费观看毛片| 一区二区三区激情视频| 高清在线视频一区二区三区| 午夜免费观看性视频| 在线天堂中文资源库| 天堂中文最新版在线下载| 久久性视频一级片| 国产97色在线日韩免费| 黑人欧美特级aaaaaa片| 最近手机中文字幕大全| 亚洲精品成人av观看孕妇| www.自偷自拍.com| h视频一区二区三区| 日韩熟女老妇一区二区性免费视频| 成人毛片60女人毛片免费| 亚洲熟女精品中文字幕| 大码成人一级视频| 久久久久久久精品精品| 成年人午夜在线观看视频| 女的被弄到高潮叫床怎么办| 久久精品熟女亚洲av麻豆精品| 国产一区二区三区综合在线观看| 高清黄色对白视频在线免费看| 超色免费av| 狂野欧美激情性xxxx| 日本爱情动作片www.在线观看| 老司机影院成人| 999久久久国产精品视频| 午夜福利一区二区在线看| 侵犯人妻中文字幕一二三四区| 亚洲精品第二区| 十八禁高潮呻吟视频| 亚洲欧美成人精品一区二区| 精品国产露脸久久av麻豆| 香蕉国产在线看| 一边亲一边摸免费视频| 精品少妇久久久久久888优播| 老司机影院成人| 久久久久久人妻| 久久毛片免费看一区二区三区| 色网站视频免费| 国产一区二区三区av在线| av又黄又爽大尺度在线免费看| 亚洲色图综合在线观看| 亚洲欧美精品综合一区二区三区| 亚洲av成人不卡在线观看播放网 | 久久热在线av| 欧美成人午夜精品| 亚洲精品一二三| 亚洲精品日韩在线中文字幕| 2018国产大陆天天弄谢| 一本大道久久a久久精品| 婷婷色av中文字幕| 久久久久久久久久久免费av| 日本猛色少妇xxxxx猛交久久| 少妇人妻精品综合一区二区| 精品酒店卫生间| 久久精品久久久久久噜噜老黄| 我要看黄色一级片免费的| 亚洲人成网站在线观看播放| 69精品国产乱码久久久| 啦啦啦视频在线资源免费观看| 亚洲少妇的诱惑av| 久久人妻熟女aⅴ| 少妇人妻 视频| 久久国产亚洲av麻豆专区| 国产爽快片一区二区三区| 黄色一级大片看看| 99精品久久久久人妻精品| 日日爽夜夜爽网站| 18禁动态无遮挡网站| 天堂中文最新版在线下载| 青青草视频在线视频观看| 亚洲av在线观看美女高潮| 伊人亚洲综合成人网| 丁香六月欧美| 精品一区二区三卡| 亚洲精品自拍成人| 人人妻人人添人人爽欧美一区卜| 婷婷成人精品国产| 国产精品.久久久| 日韩大片免费观看网站| 欧美日韩精品网址| 免费在线观看完整版高清| 99热国产这里只有精品6| 国产精品熟女久久久久浪| 国产不卡av网站在线观看| 观看av在线不卡| 无限看片的www在线观看| 另类亚洲欧美激情| 美女主播在线视频| 精品国产国语对白av| 国产精品香港三级国产av潘金莲 | 18在线观看网站| 亚洲婷婷狠狠爱综合网| 在线观看免费视频网站a站| 一区在线观看完整版| 2021少妇久久久久久久久久久| 国产不卡av网站在线观看| 日韩制服骚丝袜av| 777久久人妻少妇嫩草av网站| 亚洲欧美成人精品一区二区| 男人添女人高潮全过程视频| 精品福利永久在线观看| 99久久99久久久精品蜜桃| www.自偷自拍.com| 欧美在线一区亚洲| 亚洲成人手机| 国产日韩欧美在线精品| 看非洲黑人一级黄片| 久久久久网色| 国产高清不卡午夜福利| 国产片内射在线| 中文字幕制服av| 午夜福利在线免费观看网站| 亚洲欧洲日产国产| 晚上一个人看的免费电影| 肉色欧美久久久久久久蜜桃| 亚洲精品国产色婷婷电影| 18禁动态无遮挡网站| 免费观看人在逋| 国产精品一区二区在线不卡| 九九爱精品视频在线观看| 国产成人精品无人区| 国产99久久九九免费精品| 在线观看一区二区三区激情| 午夜日本视频在线| 亚洲成人手机| 亚洲国产成人一精品久久久| 国产一区二区激情短视频 | 免费日韩欧美在线观看| 亚洲一码二码三码区别大吗| av女优亚洲男人天堂| 欧美xxⅹ黑人| 下体分泌物呈黄色| 午夜日本视频在线| 热re99久久精品国产66热6| 最新的欧美精品一区二区| 亚洲五月色婷婷综合| 亚洲精品日韩在线中文字幕| 久久精品久久久久久久性| 18在线观看网站| 韩国av在线不卡| 日韩av在线免费看完整版不卡| 黄色视频在线播放观看不卡| 亚洲精品久久久久久婷婷小说| 久久人人爽人人片av| 这个男人来自地球电影免费观看 | 中文字幕人妻熟女乱码| 老熟女久久久| 各种免费的搞黄视频| 国产精品久久久av美女十八| 国产亚洲欧美精品永久| 日韩欧美一区视频在线观看| 免费av中文字幕在线| 赤兔流量卡办理| 国产亚洲午夜精品一区二区久久| 亚洲成人国产一区在线观看 | 亚洲精品中文字幕在线视频| 操出白浆在线播放| av线在线观看网站| 日韩大片免费观看网站| avwww免费| 欧美少妇被猛烈插入视频| 少妇被粗大的猛进出69影院| 黄片播放在线免费| 女性生殖器流出的白浆| 国语对白做爰xxxⅹ性视频网站| 精品国产露脸久久av麻豆| 欧美人与性动交α欧美软件| 欧美精品人与动牲交sv欧美| 亚洲三区欧美一区| 香蕉国产在线看| 欧美 日韩 精品 国产| 久久久久精品久久久久真实原创| 丁香六月天网| 美女扒开内裤让男人捅视频| 亚洲成人国产一区在线观看 | xxxhd国产人妻xxx| 久久精品人人爽人人爽视色| netflix在线观看网站| 一级毛片黄色毛片免费观看视频| 亚洲国产精品成人久久小说| 亚洲人成电影观看| 狠狠精品人妻久久久久久综合| 精品国产一区二区久久| 在线观看一区二区三区激情| 欧美少妇被猛烈插入视频| 宅男免费午夜| av国产精品久久久久影院| 女人久久www免费人成看片| 9191精品国产免费久久| 国产免费福利视频在线观看| 免费黄色在线免费观看| 如日韩欧美国产精品一区二区三区| 晚上一个人看的免费电影| 精品国产一区二区久久| 国产免费又黄又爽又色| 在线观看国产h片| 久久久久精品人妻al黑| 欧美乱码精品一区二区三区| 午夜福利网站1000一区二区三区| 欧美成人午夜精品| 精品国产乱码久久久久久小说| 国产精品久久久久久久久免| 99久久精品国产亚洲精品| 午夜影院在线不卡| 中文字幕最新亚洲高清| 美女高潮到喷水免费观看| 好男人视频免费观看在线| 中文字幕人妻丝袜制服| 一二三四在线观看免费中文在| 免费人妻精品一区二区三区视频| 纯流量卡能插随身wifi吗| av片东京热男人的天堂| 久久婷婷青草| 亚洲欧美清纯卡通| 99精国产麻豆久久婷婷| a 毛片基地| 久久性视频一级片| 亚洲精品一区蜜桃| xxx大片免费视频| 赤兔流量卡办理| 久久久久久久精品精品| 国产精品av久久久久免费| 美女国产高潮福利片在线看| 两性夫妻黄色片| 狠狠精品人妻久久久久久综合| 一级毛片 在线播放| 亚洲精品久久成人aⅴ小说| 欧美久久黑人一区二区| 夫妻性生交免费视频一级片| 下体分泌物呈黄色| 亚洲熟女精品中文字幕| 丁香六月天网| 久久久久精品国产欧美久久久 | 亚洲综合精品二区| 人人妻人人澡人人爽人人夜夜| 精品一区二区三区av网在线观看 | 在线观看人妻少妇| 日韩伦理黄色片| 欧美亚洲日本最大视频资源| 只有这里有精品99| 99热全是精品| 久久久亚洲精品成人影院| www日本在线高清视频| 美女脱内裤让男人舔精品视频| 午夜日韩欧美国产| 国产成人精品福利久久| 天美传媒精品一区二区| 黄片小视频在线播放| 男女国产视频网站| 久久精品熟女亚洲av麻豆精品| 各种免费的搞黄视频| 亚洲少妇的诱惑av| 菩萨蛮人人尽说江南好唐韦庄| 三上悠亚av全集在线观看| 国产精品亚洲av一区麻豆 | 国产成人a∨麻豆精品| 亚洲欧美一区二区三区黑人| 久久精品久久久久久久性| 日韩制服骚丝袜av| 欧美变态另类bdsm刘玥| 美女福利国产在线| 欧美精品人与动牲交sv欧美| 国产精品蜜桃在线观看| 在线观看三级黄色| 欧美人与善性xxx| 亚洲欧美一区二区三区久久| 亚洲激情五月婷婷啪啪| 自线自在国产av| 又粗又硬又长又爽又黄的视频| 人妻一区二区av| 午夜日韩欧美国产| 一级,二级,三级黄色视频| 三上悠亚av全集在线观看| xxx大片免费视频| 人妻一区二区av| 丰满少妇做爰视频| 国产黄色免费在线视频| 人人妻,人人澡人人爽秒播 | 美女福利国产在线| 欧美精品一区二区大全| 中文字幕亚洲精品专区| 大片电影免费在线观看免费| 欧美日本中文国产一区发布| 老司机亚洲免费影院| 91aial.com中文字幕在线观看| 国产日韩欧美亚洲二区| 久久影院123| 色网站视频免费| 亚洲精品第二区| 精品国产一区二区久久| 80岁老熟妇乱子伦牲交| 丝袜美足系列| 亚洲av中文av极速乱| 国产深夜福利视频在线观看| 男女高潮啪啪啪动态图| 久久女婷五月综合色啪小说| 国产精品免费视频内射| 亚洲国产av新网站| 国产又色又爽无遮挡免| 久久精品国产综合久久久| 亚洲av中文av极速乱| 亚洲伊人色综图| 久久久久精品国产欧美久久久 | 亚洲少妇的诱惑av| 在线观看国产h片| 人体艺术视频欧美日本| 国产亚洲一区二区精品| 久久久久视频综合| 91精品三级在线观看| 久久精品亚洲熟妇少妇任你| 免费观看a级毛片全部| 精品一区二区三卡| 久久 成人 亚洲| 一级片免费观看大全| 亚洲美女搞黄在线观看| 国产亚洲最大av| 色网站视频免费| 天堂俺去俺来也www色官网| 日韩一卡2卡3卡4卡2021年| 免费看不卡的av| 一区二区三区乱码不卡18| 亚洲精品,欧美精品| 欧美av亚洲av综合av国产av | 狠狠婷婷综合久久久久久88av| 性色av一级| 日韩中文字幕视频在线看片| 久久精品国产亚洲av涩爱| 国产老妇伦熟女老妇高清| av福利片在线| av视频免费观看在线观看| netflix在线观看网站| 免费人妻精品一区二区三区视频| 国产精品麻豆人妻色哟哟久久| 校园人妻丝袜中文字幕| 亚洲精品aⅴ在线观看| 自线自在国产av| 中文字幕最新亚洲高清| 你懂的网址亚洲精品在线观看| 亚洲精品在线美女| 久久久国产一区二区| videosex国产| 日韩中文字幕欧美一区二区 | a级毛片黄视频| 亚洲欧洲国产日韩| 国产不卡av网站在线观看| 亚洲欧美清纯卡通| 欧美人与性动交α欧美精品济南到| 精品一品国产午夜福利视频| 99久久99久久久精品蜜桃| 青草久久国产| 精品视频人人做人人爽| 王馨瑶露胸无遮挡在线观看| 亚洲精品,欧美精品| 日韩视频在线欧美| 嫩草影院入口| 91国产中文字幕| 久久狼人影院| 色播在线永久视频| 两性夫妻黄色片| 男女无遮挡免费网站观看| av卡一久久| 亚洲精品久久久久久婷婷小说| 国产精品麻豆人妻色哟哟久久| 亚洲色图综合在线观看| 亚洲激情五月婷婷啪啪| 波多野结衣av一区二区av| 久久97久久精品| 这个男人来自地球电影免费观看 | 国产精品女同一区二区软件| 国产成人精品无人区| 99久久人妻综合| 国产成人av激情在线播放| 中文天堂在线官网| 涩涩av久久男人的天堂| 久久精品久久精品一区二区三区| 日韩一本色道免费dvd| 丝袜喷水一区| 国产成人免费无遮挡视频| 欧美亚洲 丝袜 人妻 在线| 在线观看一区二区三区激情| 91精品伊人久久大香线蕉| 黄片无遮挡物在线观看| h视频一区二区三区| 美国免费a级毛片| 美女主播在线视频| kizo精华| 男人爽女人下面视频在线观看| 色婷婷av一区二区三区视频| 国产在线视频一区二区|