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

    4級(jí)重復(fù)級(jí)低速模擬壓氣機(jī)的試驗(yàn)與計(jì)算

    2018-09-17 06:09:00劉寶杰于賢君
    航空發(fā)動(dòng)機(jī) 2018年5期
    關(guān)鍵詞:靜子總壓壓氣機(jī)

    劉寶杰 ,張 帥 ,于賢君

    (1.北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院北京100191;2.北京航空航天大學(xué)航空發(fā)動(dòng)機(jī)氣動(dòng)熱力國(guó)防重點(diǎn)試驗(yàn)室,北京100191;3.先進(jìn)航空發(fā)動(dòng)機(jī)協(xié)同創(chuàng)新中心,北京100191)

    0 引言

    高壓壓氣機(jī)的研制一直是現(xiàn)代航空發(fā)動(dòng)機(jī)設(shè)計(jì)中的關(guān)鍵問(wèn)題之一[1]。在當(dāng)前的設(shè)計(jì)體系下,高性能高壓壓氣機(jī)設(shè)計(jì)很大程度上依賴于對(duì)其內(nèi)部復(fù)雜流場(chǎng)的認(rèn)識(shí),由于受其轉(zhuǎn)速高、葉片相對(duì)尺寸小、結(jié)構(gòu)緊湊以及高溫高壓的環(huán)境、昂貴的試驗(yàn)費(fèi)用等的限制,采用常規(guī)方法很難詳細(xì)測(cè)量高壓壓氣機(jī)內(nèi)部流場(chǎng)。這也直接限制了壓氣機(jī)設(shè)計(jì)體系的發(fā)展和完善。

    20世紀(jì)70年代中后期發(fā)展起來(lái)的低速模擬試驗(yàn)研究技術(shù),實(shí)現(xiàn)了在低速環(huán)境下模擬高壓壓氣機(jī)后面級(jí)高速環(huán)境流動(dòng),試驗(yàn)成本大大降低,而且便于進(jìn)行內(nèi)部流場(chǎng)的詳細(xì)測(cè)量。D.C.Wisler等[2-3]最早在GE公司E3計(jì)劃支持下,首次提出保持雷諾數(shù)、稠度、葉表無(wú)量綱速度分布等參數(shù)相同的低速模型,建立了4級(jí)重復(fù)級(jí)低速大尺寸壓氣機(jī)試驗(yàn)臺(tái)并以此為依托發(fā)展了定制裁剪葉型和3維造型技術(shù);PW公司[4-5]、NASA[6]以及RR公司聯(lián)合Cambridge大學(xué)和Cranfield大學(xué)[7-9]等,也開(kāi)展了低速模擬技術(shù)的研究。低速模擬技術(shù)的飛速發(fā)展對(duì)20世紀(jì)七八十年代壓氣機(jī)設(shè)計(jì)水平的進(jìn)步提供了關(guān)鍵支持[10]。

    國(guó)內(nèi)對(duì)于低速模擬試驗(yàn)技術(shù)也做了大量工作,北京航空航天大學(xué)在20世紀(jì)80年代建立了國(guó)內(nèi)首個(gè)低速大尺寸壓氣機(jī)試驗(yàn)臺(tái)[11],并研發(fā)了1整套轉(zhuǎn)子流場(chǎng)動(dòng)態(tài)測(cè)試技術(shù);南京航空航天大學(xué)在2010年前后,按照GE公司的思路建立了國(guó)內(nèi)首個(gè)4級(jí)重復(fù)級(jí)低速大尺寸壓氣機(jī)試驗(yàn)臺(tái)[12],并發(fā)展了相關(guān)的試驗(yàn)測(cè)量技術(shù)[13];中科院工程熱物理所在近年也進(jìn)行了低速模擬技術(shù)的研究并設(shè)計(jì)和建造了低速大尺寸壓氣機(jī)試驗(yàn)臺(tái)[14];上海交通大學(xué)和中國(guó)航發(fā)沈陽(yáng)發(fā)動(dòng)機(jī)研究所合作開(kāi)展了高壓壓氣機(jī)進(jìn)口級(jí)低速模擬設(shè)計(jì)的數(shù)值分析[15]。國(guó)內(nèi)對(duì)于壓氣機(jī)低速模擬技術(shù)的發(fā)展比國(guó)外落后一些,除了需要建立完備的低速模擬技術(shù)設(shè)計(jì)體系外,完善準(zhǔn)確的試驗(yàn)測(cè)量也迫在眉睫。此外,當(dāng)前對(duì)于壓氣機(jī)設(shè)計(jì)人員而言,在對(duì)設(shè)計(jì)方案進(jìn)行最終試驗(yàn)驗(yàn)證前往往通過(guò)商用軟件校核,但數(shù)值模擬軟件的可靠性有待商榷,而完整的低速模擬試驗(yàn)數(shù)據(jù)提供了軟件校核的細(xì)致數(shù)據(jù)。

    本文在北京航空航天大學(xué)新建立的4級(jí)重復(fù)級(jí)低速模擬壓氣機(jī)試驗(yàn)臺(tái)[16]上搭建了全面的測(cè)試系統(tǒng),對(duì)試驗(yàn)臺(tái)做了細(xì)致測(cè)量,同時(shí)應(yīng)用試驗(yàn)數(shù)據(jù)對(duì)常用的數(shù)值模擬軟件進(jìn)行了校驗(yàn)分析。

    1 試驗(yàn)設(shè)備及試驗(yàn)方案

    1.1 試驗(yàn)設(shè)備簡(jiǎn)介

    為了系統(tǒng)地開(kāi)展4級(jí)重復(fù)級(jí)低速模擬研究,于賢君、劉寶杰等完成了模擬某先進(jìn)發(fā)動(dòng)機(jī)后面級(jí)典型流動(dòng)的氣動(dòng)方案設(shè)計(jì)[16]的同時(shí),與某研究所共同完成了4級(jí)重復(fù)級(jí)試驗(yàn)臺(tái)的結(jié)構(gòu)設(shè)計(jì)、加工和裝配以及試驗(yàn)測(cè)量。試驗(yàn)臺(tái)的結(jié)構(gòu)模型如圖1所示。試驗(yàn)段結(jié)構(gòu)由導(dǎo)葉和4級(jí)幾何參數(shù)相同的葉片組成,試驗(yàn)葉片加工采用樹(shù)脂加玻璃纖維強(qiáng)化的方式,試驗(yàn)臺(tái)(下文簡(jiǎn)稱4Stages)的主要設(shè)計(jì)參數(shù)見(jiàn)表1。

    1.2 試驗(yàn)測(cè)量方案

    本文系統(tǒng)測(cè)量了4級(jí)重復(fù)級(jí)壓氣機(jī)的特性、級(jí)間氣動(dòng)參數(shù)的展向分布以及靜子后的流場(chǎng)。試驗(yàn)臺(tái)的所有測(cè)量截面如圖2所示。

    圖1 4級(jí)重復(fù)級(jí)低速模擬壓氣機(jī)結(jié)構(gòu)模型

    表1 4級(jí)重復(fù)級(jí)低速模擬壓氣機(jī)主要設(shè)計(jì)參數(shù)

    圖2 4級(jí)重復(fù)級(jí)測(cè)量截面

    本文測(cè)量的壓氣機(jī)的特性是總壓升系數(shù)和效率隨流量系數(shù)的變化曲線,通過(guò)調(diào)節(jié)壓氣機(jī)出口機(jī)械閥門的開(kāi)度來(lái)控制壓氣機(jī)的流量狀態(tài)。為了測(cè)量試驗(yàn)臺(tái)的流量系數(shù),通過(guò)在進(jìn)口Inlet截面均勻布置的4個(gè)壁面靜壓孔測(cè)量到的平均靜壓Ps,in,結(jié)合通過(guò)振筒氣壓計(jì)和高精度溫度傳感器測(cè)量的壓氣機(jī)進(jìn)口的大氣壓強(qiáng) Pt,in和溫度 Tt,in,通過(guò)式(1)求得

    式中:Umid為轉(zhuǎn)子葉中切線速度。

    總壓升系數(shù)定義為

    式中:Pt,ave為試驗(yàn)臺(tái)出口截面的平均總壓,通過(guò)在出口截面布置的8個(gè)總壓耙通過(guò)面積加權(quán)的方式求得。

    本文測(cè)量的壓氣機(jī)效率為扭矩效率

    式中:π為壓氣機(jī)壓比;Win為壓氣機(jī)進(jìn)口物理流量;n為物理轉(zhuǎn)速;k為定比熱;R為氣體參數(shù);M為扭矩。

    本文通過(guò)控制位移機(jī)構(gòu)帶動(dòng)5孔氣動(dòng)探針移動(dòng)到待測(cè)點(diǎn),并設(shè)置等待時(shí)間,待流場(chǎng)穩(wěn)定后采集5個(gè)孔壓力數(shù)據(jù),經(jīng)過(guò)數(shù)據(jù)處理獲取級(jí)間各截面的氣動(dòng)參數(shù),5孔針數(shù)據(jù)處理方法見(jiàn)文獻(xiàn)[17]。測(cè)量用的5孔探針?lè)糯笳掌娃D(zhuǎn)子后測(cè)量用的位移機(jī)構(gòu)如圖3所示,轉(zhuǎn)子后測(cè)點(diǎn)沿徑向的布置如圖4所示。近機(jī)匣第1個(gè)測(cè)點(diǎn)距離機(jī)匣2 mm,近輪轂第1個(gè)測(cè)點(diǎn)距離輪轂3 mm,主流區(qū)域的測(cè)點(diǎn)間隔約5%葉高,在葉根和葉尖區(qū)域加密,共布置24個(gè)測(cè)點(diǎn)。

    圖3 試驗(yàn)用5孔探針?lè)糯蠹稗D(zhuǎn)子后測(cè)量用位移機(jī)構(gòu)

    圖4 轉(zhuǎn)子后測(cè)點(diǎn)沿徑向的布置

    靜子后的流場(chǎng)測(cè)量采用2維位移機(jī)構(gòu)控制5孔探針掃描靜子出口扇形截面的方式進(jìn)行。靜子后測(cè)量用的位移機(jī)構(gòu)和測(cè)量網(wǎng)格如圖5所示。測(cè)量節(jié)點(diǎn)共18×17(306)個(gè),覆蓋了1.1倍靜子周期,并在尾跡區(qū)域和端壁進(jìn)行加密。與輪轂最近的測(cè)點(diǎn)距離為4 mm,與機(jī)匣最近的測(cè)點(diǎn)距離為2 mm。靜子出口氣動(dòng)參數(shù)的展向分布通過(guò)相同半徑的測(cè)點(diǎn)數(shù)據(jù)進(jìn)行加權(quán)得到。

    1.3 5孔針重復(fù)性測(cè)量和精度

    為了考察5孔探針測(cè)量的重復(fù)性,對(duì)第2級(jí)轉(zhuǎn)子(R2)后的設(shè)計(jì)狀態(tài)(流量系數(shù)為0.62)和第2級(jí)靜子(S2)后的堵塞狀態(tài)(流量系數(shù)為0.68)進(jìn)行重復(fù)測(cè)量。第2級(jí)轉(zhuǎn)子后氣流角和流量系數(shù)(Vz/Umid)沿展向分布的對(duì)比如圖6所示。從圖中可見(jiàn),前后2次試驗(yàn)數(shù)據(jù)具有高度的一致性。

    圖5 靜子后測(cè)量用位移機(jī)構(gòu)及測(cè)量網(wǎng)格

    圖6 設(shè)計(jì)狀態(tài)下第2級(jí)轉(zhuǎn)子展向分布的重復(fù)性對(duì)比

    圖7 堵塞狀態(tài)下第2級(jí)靜子展向分布的重復(fù)性對(duì)比

    堵塞狀態(tài)(流量系數(shù)為0.68)下第2級(jí)靜子后氣流角和流量系數(shù)展向參數(shù)以及流場(chǎng)的重復(fù)性對(duì)比如圖7、8所示。從圖中可見(jiàn),靜子出口流場(chǎng)周期性良好,流場(chǎng)結(jié)構(gòu)和氣動(dòng)參數(shù)展向分布都具有很強(qiáng)的重復(fù)性。5孔針測(cè)量的精度見(jiàn)表2。

    圖8 堵塞狀態(tài)下第2級(jí)靜子后無(wú)量綱軸向速度重復(fù)性對(duì)比

    表2 5孔針測(cè)量精度

    圖9 數(shù)值模擬采用的計(jì)算網(wǎng)格

    2 數(shù)值模擬方法

    當(dāng)前的設(shè)計(jì)體系的初步方案校驗(yàn)主要采用3維數(shù)值模擬軟件。本文利用試驗(yàn)數(shù)據(jù)校驗(yàn)了3種常用的計(jì)算程序:Turbo程序和商用軟件CFX、Numeca。

    Turbo程序最初由Denton[18]開(kāi)發(fā),湍流效應(yīng)通過(guò)采用混合長(zhǎng)度模型考慮,氣體的黏性通過(guò)薄剪切層近似模擬,計(jì)算程序可以實(shí)現(xiàn)在可接受的精度下盡量占用少的計(jì)算資源,迅速得到計(jì)算結(jié)果,從而幫助設(shè)計(jì)者很快地調(diào)整設(shè)計(jì)方案。本文采用的計(jì)算程序是劉寶杰等[19]通過(guò)對(duì)大量試驗(yàn)結(jié)果校驗(yàn)改進(jìn)的版本TURBO-N7。Turbo采用簡(jiǎn)單的H型網(wǎng)格,如圖9上圖所示。這里只給出導(dǎo)葉和第1級(jí)網(wǎng)格,因?yàn)楹竺?級(jí)網(wǎng)格與第1級(jí)網(wǎng)格參數(shù)完全一致,IGV網(wǎng)格塊在流向、周向和徑向分別有175、37和45個(gè)節(jié)點(diǎn),Stator和Rotor的網(wǎng)格節(jié)點(diǎn)數(shù)是138、37和45個(gè),第4級(jí)靜子的網(wǎng)格節(jié)點(diǎn)數(shù)為178、37和45個(gè)。計(jì)算網(wǎng)格總數(shù)約217.1萬(wàn)。網(wǎng)格在壁面進(jìn)行適當(dāng)加密y+=11左右,以滿足湍流模型的計(jì)算要求。

    CFX與Numeca計(jì)算采用同一套網(wǎng)格,如圖9的下圖所示。使用Autogrid5模塊生成,葉片主流通道采用典型的O6H型網(wǎng)格拓?fù)浣Y(jié)構(gòu),間隙內(nèi)采用OH蝶型網(wǎng)格,從而保證網(wǎng)格的正交性,同時(shí)在葉片前尾緣、葉片表面、間隙、壁面適當(dāng)加密,網(wǎng)格的Y+絕大部分控制在1左右,總網(wǎng)格數(shù)為767萬(wàn),保證了對(duì)流場(chǎng)特征的捕捉。Numeca計(jì)算采用S-A湍流模型,CFX采用kω湍流模型。

    3個(gè)計(jì)算程序的計(jì)算求解都是給定了進(jìn)口的總壓、總溫和氣流角,計(jì)算域出口給定的是反壓。

    3 試驗(yàn)結(jié)果分析

    3.1 壓氣機(jī)特性

    試驗(yàn)測(cè)量得到的壓氣機(jī)特性如圖10所示。設(shè)計(jì)點(diǎn)流量系數(shù)為0.62,效率為88.7%,總壓升系數(shù)為3.06,峰值效率約89%,綜合裕度為39.04%,這些性能指標(biāo)對(duì)于當(dāng)前的壓氣機(jī)設(shè)計(jì)水平來(lái)說(shuō)具有典型性,并且具有較高水平,證明該低速模擬壓氣機(jī)的高壓原型同樣具有良好的性能指標(biāo)。本文沿特性線選取4個(gè)狀態(tài)點(diǎn)分別為近堵點(diǎn)(NC)、設(shè)計(jì)點(diǎn)(DE)、中間狀態(tài)點(diǎn)(MID)和近失速點(diǎn)(NS),流量系數(shù)分別對(duì)應(yīng)0.68、0.62、0.54、0.46,進(jìn)行級(jí)間參數(shù)和流場(chǎng)的測(cè)量。

    圖10 4級(jí)重復(fù)級(jí)試驗(yàn)特性

    3.2 展向氣動(dòng)參數(shù)分布

    各級(jí)轉(zhuǎn)子出口流量系數(shù)展向分布在4個(gè)流動(dòng)狀態(tài)下的對(duì)比如圖11所示。首先在NC和DE流量較大的狀態(tài)下第1級(jí)轉(zhuǎn)子葉尖區(qū)域的端壁堵塞較弱,出口流量的展向分布比較均勻,經(jīng)過(guò)第1級(jí)的充分發(fā)展后端壁的流通能力在第2~4級(jí)轉(zhuǎn)子出截面尖部區(qū)域基本相似;在MID和NS流量較小的狀態(tài)下,4級(jí)轉(zhuǎn)子出口端壁的流通能力基本相當(dāng);在所有工況下,4級(jí)轉(zhuǎn)子的根部流通能力逐漸降低。盡管各級(jí)轉(zhuǎn)子的流量展向分布有些許差異,但是總體來(lái)看,在所有工況下,第2~4級(jí)轉(zhuǎn)子出口的流量沿展向的分布非常近似。第2~4級(jí)轉(zhuǎn)子的氣流折轉(zhuǎn)角(如圖12所示)和D因子(如圖13所示)的展向分布也是一致的,說(shuō)明2~4級(jí)轉(zhuǎn)子內(nèi)的流動(dòng)也基本一致。

    圖11 不同流量狀態(tài)下各級(jí)轉(zhuǎn)子出口流量系數(shù)展向分布對(duì)比

    圖12 不同流量狀態(tài)下各級(jí)轉(zhuǎn)子氣流折轉(zhuǎn)角展向分布對(duì)比

    圖13 不同流量狀態(tài)下各級(jí)轉(zhuǎn)子D因子展向分布對(duì)比

    前3級(jí)靜子在機(jī)匣輪轂處都沒(méi)有間隙,第4級(jí)靜子是懸臂安裝,在輪轂一側(cè)有1 mm的間隙。所以由流量展向分布來(lái)看,4個(gè)狀態(tài)下第4級(jí)靜子根部的流通能力較弱,這是由于泄漏造成的。不同流量狀態(tài)下各級(jí)靜子出口流量系數(shù)、D因子展向分布對(duì)比分別如圖14、15所示。從圖中可見(jiàn),第2、3級(jí)靜子出口截面流量的展向分布有很高的一致性,靜子D因子沿展向的分布也是相同的規(guī)律,第2、3級(jí)靜子基本一致。結(jié)合轉(zhuǎn)子的流動(dòng)特點(diǎn),本文的4級(jí)重復(fù)級(jí)在第2、3級(jí)流場(chǎng)有很好的重復(fù)性。

    圖14 不同流量狀態(tài)下各級(jí)靜子出口流量系數(shù)展向分布對(duì)比

    圖15 不同流量狀態(tài)下各級(jí)靜子D因子展向分布對(duì)比

    圖16 近堵點(diǎn)各級(jí)靜子出口無(wú)量綱軸向速度

    3.3 靜子出口流場(chǎng)

    4種狀態(tài)下各級(jí)靜子出口的流場(chǎng)如圖16~19所示。在每個(gè)流動(dòng)狀態(tài)都清晰可見(jiàn),第4級(jí)靜子葉根間隙泄漏所形成流動(dòng)堵塞,在每個(gè)狀態(tài)下靜子出口流場(chǎng)結(jié)構(gòu)清晰可辨且結(jié)構(gòu)比較一致;隨著流量系數(shù)的減小,靜子的尾跡均勻增厚,但是并未發(fā)生明顯的角區(qū)分離,可以看出本文4級(jí)重復(fù)級(jí)壓氣機(jī)靜子的流場(chǎng)組織合理,是比較理想的壓氣機(jī)設(shè)計(jì),也是高效率高壓壓氣機(jī)的必要特征。

    圖17 設(shè)計(jì)點(diǎn)各級(jí)靜子出口無(wú)量綱軸向速度

    圖18 中間狀態(tài)點(diǎn)各級(jí)靜子出口無(wú)量綱軸向速度

    圖19 近失速點(diǎn)各級(jí)靜子出口無(wú)量綱軸向速度

    3.4 低速模擬與高速原型對(duì)比

    4級(jí)重復(fù)級(jí)第3級(jí)轉(zhuǎn)靜子進(jìn)、出口的氣流角與高速原型的對(duì)比如圖20所示。從圖中可見(jiàn),轉(zhuǎn)靜子的進(jìn)口氣流角比較吻合,出口氣流角有所偏差,但是分布規(guī)律較相似,這個(gè)偏差是由于低速情況下氣流是不可壓的,增加氣流的折轉(zhuǎn)來(lái)補(bǔ)償高速氣流的壓縮性造成的。

    圖20 低速模擬試驗(yàn)進(jìn)、出口氣流角與高速原型對(duì)比

    圖21 低速模擬試驗(yàn)轉(zhuǎn)、靜子D因子與高速原型對(duì)比

    4級(jí)重復(fù)級(jí)第3級(jí)轉(zhuǎn)靜子D因子與高速原型的對(duì)比如圖21所示。從圖中可見(jiàn),轉(zhuǎn)靜子的D因子與原型的分布規(guī)律基本相似,低速模擬試驗(yàn)結(jié)果D因子稍偏大。總體來(lái)看重復(fù)級(jí)的設(shè)計(jì)較好地模擬出高壓原型的氣動(dòng)參數(shù)的分布情況,其負(fù)荷水平和內(nèi)部流動(dòng)特點(diǎn)相對(duì)于當(dāng)前的壓氣機(jī)設(shè)計(jì)水平來(lái)說(shuō)具有典型性。

    綜上所述,4級(jí)重復(fù)級(jí)試驗(yàn)臺(tái)是為了研究高壓壓氣機(jī)后面級(jí)中的復(fù)雜流動(dòng),導(dǎo)葉的作用是為了給定下游轉(zhuǎn)子的工作條件,前2級(jí)是最為重要的試驗(yàn)研究級(jí),第4級(jí)是為了提供下游流場(chǎng)環(huán)境,減少設(shè)計(jì)和加工成本,通常選用重復(fù)級(jí)的形式。重復(fù)級(jí)最為核心的設(shè)計(jì)準(zhǔn)則是在合適的雷諾數(shù)和馬赫數(shù)條件下,保證進(jìn)口來(lái)流條件相似的情況下,保證壓氣機(jī)級(jí)的進(jìn)出口氣流角和擴(kuò)散因子等氣動(dòng)參數(shù)相似。在保證上述相關(guān)參數(shù)相似的情況下,其它氣動(dòng)參數(shù)基本相似。需要指出的是,這里無(wú)法保證速度三角形完全相同,只能保證級(jí)進(jìn)出口氣流角相同。為了保證負(fù)荷相似,葉片的彎角增大,因此轉(zhuǎn)子/靜子的出口/進(jìn)口氣流角必然不同,但這種不同是可以預(yù)知的,并且只有這樣才能保證高速和低速情況的流動(dòng)基本相似。

    圖22 4級(jí)重復(fù)級(jí)數(shù)值計(jì)算與試驗(yàn)的特性對(duì)比

    4 數(shù)值軟件的校驗(yàn)

    4.1 特性

    3維數(shù)值模擬結(jié)果與試驗(yàn)特性的對(duì)比如圖22所示。從圖中可見(jiàn),對(duì)于效率特性,在大流量狀態(tài)下,CFX、Numeca計(jì)算值比較吻合,Turbo計(jì)算值偏低;在小流量狀態(tài)下只有Numeca計(jì)算值相差較小,CFX計(jì)算的效率整體趨勢(shì)比較吻合,但是效率預(yù)估的太過(guò)理想,Turbo計(jì)算值也偏理想。對(duì)于總壓升特性,在大流量狀態(tài)下,Turbo和Numeca計(jì)算值比較吻合,但CFX計(jì)算值偏小;在小流量狀態(tài)下,與試驗(yàn)值相比,3種軟件計(jì)算值都偏高,Numeca計(jì)算值偏離得最遠(yuǎn),CFX的最接近,Turbo的居中。3種計(jì)算軟件在試驗(yàn)的近失速狀態(tài)的算例均不能收斂,Numeca計(jì)算得到的穩(wěn)定工作的流量范圍最窄,Turbo的稍有增加,CFX計(jì)算的穩(wěn)定工作的流量范圍最寬但距離試驗(yàn)值還有一定差距。計(jì)算得到的重復(fù)級(jí)的綜合裕度,Numeca的最低為21.47%,Turbo和CFX的相當(dāng),分別為 31.17%和32.94%,但是均比試驗(yàn)值39.04%的低。

    由于是重復(fù)級(jí)設(shè)計(jì),由前述試驗(yàn)數(shù)據(jù)分析,實(shí)際上各級(jí)的工作狀態(tài)基本相似,在中間狀態(tài)下Numeca和CFX計(jì)算的單級(jí)特性對(duì)比如圖23所示。相對(duì)于Numeca而言,CFX計(jì)算的總特性效率偏高、總壓升偏低的特點(diǎn)在每一級(jí)呈現(xiàn)相同的規(guī)律。下面進(jìn)行第3級(jí)在設(shè)計(jì)狀態(tài)和中間狀態(tài)下氣動(dòng)參數(shù)展向分布和流場(chǎng)的3種軟件計(jì)算值與試驗(yàn)結(jié)果的對(duì)比分析。

    圖23 中間狀態(tài)下單級(jí)特性Numeca與CFX計(jì)算值對(duì)比

    圖24 設(shè)計(jì)狀態(tài)下展向參數(shù)對(duì)比

    4.2 設(shè)計(jì)狀態(tài)

    從圖22中可知,在設(shè)計(jì)狀態(tài)下的壓氣機(jī)總壓升特性,Turbo計(jì)算值與試驗(yàn)值相近,CFX的偏小,Numeca的偏大。轉(zhuǎn)子的總壓升系數(shù)展向分布及靜子的損失沿展向分布的對(duì)比如圖24所示。從圖中可見(jiàn),對(duì)于轉(zhuǎn)子總壓升系數(shù),Numeca計(jì)算值比試驗(yàn)值整體偏高,靜子的損失水平相當(dāng);而Turbo轉(zhuǎn)子總壓升系數(shù)計(jì)算值在尖部比試驗(yàn)值低,但靜子在尖部的總壓損失也比試驗(yàn)值低,所以總壓升特性與試驗(yàn)值相近;而CFX計(jì)算的轉(zhuǎn)子總壓升偏小,靜子損失水平基本相當(dāng),故而計(jì)算得到的總壓升特性比試驗(yàn)值偏低。

    在亞聲速轉(zhuǎn)子中,轉(zhuǎn)子的總壓升與氣流的折轉(zhuǎn)角成正比,計(jì)算和試驗(yàn)得到的轉(zhuǎn)子攻角、落后角沿展向分布的對(duì)比如圖25所示,在設(shè)計(jì)狀態(tài)下轉(zhuǎn)子效率和轉(zhuǎn)子氣流折轉(zhuǎn)角的展向分布對(duì)比如圖26所示。從圖中可見(jiàn),轉(zhuǎn)子攻角落后角計(jì)算值比試驗(yàn)值都有不同程度的偏差,Numeca計(jì)算值最接近,但是根部氣流折轉(zhuǎn)角比試驗(yàn)值大,CFX與Turbo計(jì)算值相當(dāng)?shù)仍囼?yàn)值偏小。也就是說(shuō)轉(zhuǎn)子壓升系數(shù)的差異本質(zhì)上是由于數(shù)值模擬軟件對(duì)于落后角計(jì)算的偏差使得葉片的工作狀態(tài)與實(shí)際工作狀態(tài)不一致造成的。

    圖25 設(shè)計(jì)狀態(tài)下試驗(yàn)值和計(jì)算值對(duì)比

    圖26 設(shè)計(jì)狀態(tài)下展向分布對(duì)比

    對(duì)于在設(shè)計(jì)狀態(tài)下壓氣機(jī)的效率特性,CFX計(jì)算值與試驗(yàn)值較接近,Numeca和Turbo計(jì)算值偏低,由于試驗(yàn)沒(méi)有得到轉(zhuǎn)子效率的展向分布,這里只給出3種軟件得到的轉(zhuǎn)子效率展向分布對(duì)比(圖26左圖),效率的差異主要集中在葉尖區(qū)域,Turbo計(jì)算的損失最大,Numeca的次之,CFX計(jì)算的尖區(qū)損失最小效率最高,正好與特性線上Turbo的最低,CFX的最高的計(jì)算結(jié)果對(duì)應(yīng)。尖區(qū)效率不一致反映了各軟件對(duì)于葉尖泄漏渦的模擬存在差異。

    4.3 中間狀態(tài)

    在中間狀態(tài)下轉(zhuǎn)子總壓升系數(shù)和靜子損失展向參數(shù)對(duì)比如圖27所示。從圖中可見(jiàn),對(duì)于在中間狀態(tài)下的壓氣機(jī)總壓升特性,3種軟件計(jì)算值都偏高(圖27左圖),Numeca的最高,Turbo的次之。對(duì)于總壓升系數(shù),Numeca計(jì)算值明顯偏高,CFX的稍高,Turbo的主流和根部稍高,70%~90%葉高偏低;但是對(duì)于靜子的總壓損失(圖27右圖),Turbo計(jì)算值明顯偏低,所以總壓升比試驗(yàn)值高,Numeca計(jì)算的根部損失偏大,尖部偏小,綜合起來(lái)壓氣機(jī)總壓升特性偏高,CFX計(jì)算值也呈相同的規(guī)律,只是偏高的程度沒(méi)有Numeca的大。對(duì)比DE和MID2種狀態(tài)下轉(zhuǎn)子總壓升系數(shù)展向分布的差別可見(jiàn),隨著流量的減小,CFD計(jì)算的總壓升的增大程度要比試驗(yàn)值的增大程度大,也意味著轉(zhuǎn)子的負(fù)荷隨流量的減小增大得較快,這也是流量邊界估計(jì)不準(zhǔn)的原因。由于Numeca計(jì)算的負(fù)荷在相同的流量系數(shù)下最大,所以Numeca失速邊界的流量最大,CFX與Turbo計(jì)算的失速邊界相近,但也不能與試驗(yàn)值相吻合。

    圖27 中間狀態(tài)下展向參數(shù)對(duì)比

    對(duì)于在中間狀態(tài)下壓氣機(jī)的效率,Numeca計(jì)算值與試驗(yàn)值比較接近,Turbo的偏高,CFX的最高。效率展向分布對(duì)比如圖28所示。Turbo計(jì)算值仍然是尖區(qū)的效率最低,但是不同的是Numeca計(jì)算的根部的效率偏低,再加上Turbo對(duì)于靜子的損失估計(jì)的理想,所以使得Numeca計(jì)算的效率比Turbo的低,CFX對(duì)于葉尖區(qū)域的計(jì)算過(guò)于理想,效率偏高,根部沒(méi)有明顯虧損,所以CFX計(jì)算的效率最高。

    綜上所述,數(shù)值計(jì)算主要在葉片的攻角落后角估計(jì)、葉尖泄漏的模擬及靜子損失估計(jì)等方面存在不足,導(dǎo)致壓氣機(jī)特性模擬產(chǎn)生偏差,同時(shí)對(duì)于失速邊界的估計(jì)也存在偏差。

    圖28 中間狀態(tài)下轉(zhuǎn)子效率展向分布計(jì)算結(jié)果對(duì)比

    5 結(jié)論

    本文對(duì)某高壓壓氣機(jī)后面級(jí)的4級(jí)重復(fù)級(jí)低速大尺寸軸流壓氣機(jī)進(jìn)行了詳細(xì)的試驗(yàn)測(cè)量和數(shù)值模擬對(duì)比分析,得到如下結(jié)論:

    (1)獲得4級(jí)重復(fù)級(jí)低速大尺寸壓氣機(jī)完整、詳細(xì)的試驗(yàn)數(shù)據(jù),利用小尺寸高精度5孔探針獲得靜子高分辨率的出口流場(chǎng),可以清晰地分辨尾跡和角區(qū)等流動(dòng)細(xì)節(jié);

    (2)流場(chǎng)測(cè)量結(jié)果表明,該4級(jí)低速模擬壓氣機(jī)實(shí)現(xiàn)了典型的重復(fù)級(jí)流動(dòng),較好地模擬了高壓壓氣機(jī)后面級(jí)原型的流動(dòng)特征,驗(yàn)證了原型壓氣機(jī)具有氣動(dòng)特征效率高和綜合裕度大的性能優(yōu)勢(shì);

    (3)利用詳細(xì)流場(chǎng)測(cè)量結(jié)果校驗(yàn)了常用3維計(jì)算軟件Numeca、CFX和Turbo,校驗(yàn)結(jié)果表明:3種軟件在估計(jì)葉片攻角、落后角、葉尖泄漏流等方面存在不足,無(wú)法準(zhǔn)確估計(jì)壓氣機(jī)的特性,特別是近失速工況;

    (4)4級(jí)重復(fù)級(jí)低速模擬壓氣機(jī)試驗(yàn)臺(tái)為高壓壓氣機(jī)后面級(jí)的關(guān)鍵技術(shù)驗(yàn)證和設(shè)計(jì)參數(shù)選擇以及設(shè)計(jì)計(jì)算軟件的校驗(yàn)提供了良好的試驗(yàn)手段及平臺(tái)。

    猜你喜歡
    靜子總壓壓氣機(jī)
    總壓探針性能結(jié)構(gòu)敏感性分析
    軸流壓氣機(jī)效率評(píng)定方法
    重型燃?xì)廨啓C(jī)壓氣機(jī)第一級(jí)轉(zhuǎn)子葉片斷裂分析
    可調(diào)式總壓耙設(shè)計(jì)及應(yīng)用
    壓氣機(jī)緊湊S形過(guò)渡段內(nèi)周向彎靜子性能數(shù)值計(jì)算
    何必喧囂慰寂寥
    亞聲速條件下總壓探針臨壁效應(yīng)的數(shù)值研究
    堅(jiān)持了十年的書(shū)信
    2 m超聲速風(fēng)洞流場(chǎng)變速壓控制方法研究
    山鄉(xiāng)一瞥
    国产免费福利视频在线观看| 亚洲中文av在线| 国产精品.久久久| 亚洲精品视频女| 亚洲自偷自拍三级| 丰满少妇做爰视频| 老女人水多毛片| 十八禁网站网址无遮挡 | 欧美日韩亚洲高清精品| 国产在线免费精品| 亚州av有码| 九九久久精品国产亚洲av麻豆| 国产高清国产精品国产三级| 97精品久久久久久久久久精品| 日本欧美视频一区| 插逼视频在线观看| 99久久精品热视频| 日韩欧美精品免费久久| 在线观看人妻少妇| 一级毛片黄色毛片免费观看视频| 婷婷色综合大香蕉| 91精品国产国语对白视频| 国产精品伦人一区二区| 久久久久国产网址| 免费播放大片免费观看视频在线观看| 国产成人免费无遮挡视频| 国产69精品久久久久777片| 少妇的逼水好多| 纵有疾风起免费观看全集完整版| 免费av中文字幕在线| 老司机影院成人| 观看免费一级毛片| 99视频精品全部免费 在线| 免费播放大片免费观看视频在线观看| 美女福利国产在线| av线在线观看网站| 久久女婷五月综合色啪小说| 十分钟在线观看高清视频www | 久久亚洲国产成人精品v| 精品人妻偷拍中文字幕| 啦啦啦在线观看免费高清www| 夜夜骑夜夜射夜夜干| 中文在线观看免费www的网站| 亚洲成色77777| 精品少妇黑人巨大在线播放| 91精品一卡2卡3卡4卡| 在线观看av片永久免费下载| 久久午夜综合久久蜜桃| 成人二区视频| 女人久久www免费人成看片| 亚洲天堂av无毛| 久久久久人妻精品一区果冻| 国产无遮挡羞羞视频在线观看| 亚洲精品一二三| 日韩av在线免费看完整版不卡| 日韩av免费高清视频| 亚洲四区av| 老女人水多毛片| 肉色欧美久久久久久久蜜桃| 精品久久国产蜜桃| 国产精品国产三级国产专区5o| 国产午夜精品一二区理论片| 22中文网久久字幕| 成人漫画全彩无遮挡| 少妇熟女欧美另类| 久久精品国产亚洲av涩爱| 亚洲天堂av无毛| 一区二区三区四区激情视频| 特大巨黑吊av在线直播| 插逼视频在线观看| 国产亚洲91精品色在线| 午夜激情福利司机影院| 99国产精品免费福利视频| 欧美日韩在线观看h| 一级二级三级毛片免费看| 日韩一区二区三区影片| 在线亚洲精品国产二区图片欧美 | 香蕉精品网在线| 国产在视频线精品| 精品一品国产午夜福利视频| 亚洲精品一二三| 少妇高潮的动态图| av福利片在线观看| 永久网站在线| 黄片无遮挡物在线观看| 国产精品国产三级国产av玫瑰| 亚洲av电影在线观看一区二区三区| 好男人视频免费观看在线| 亚洲人与动物交配视频| 丝瓜视频免费看黄片| 免费看不卡的av| 国产av精品麻豆| 免费高清在线观看视频在线观看| 黑人高潮一二区| 在线观看免费高清a一片| 亚洲av综合色区一区| 亚洲人成网站在线观看播放| 看十八女毛片水多多多| 亚洲精品乱久久久久久| 国产精品熟女久久久久浪| 精品久久久噜噜| 永久免费av网站大全| 欧美国产精品一级二级三级 | 91午夜精品亚洲一区二区三区| a级毛片免费高清观看在线播放| 国产老妇伦熟女老妇高清| 日韩一区二区三区影片| 男女啪啪激烈高潮av片| 亚洲婷婷狠狠爱综合网| a 毛片基地| 中文字幕制服av| 人妻 亚洲 视频| 在线天堂最新版资源| 亚洲av综合色区一区| 一级毛片aaaaaa免费看小| 亚洲中文av在线| 免费高清在线观看视频在线观看| 国产成人精品婷婷| 大片免费播放器 马上看| 99热6这里只有精品| 精品久久久噜噜| 看十八女毛片水多多多| 欧美人与善性xxx| 青春草视频在线免费观看| 80岁老熟妇乱子伦牲交| 精品人妻偷拍中文字幕| 91久久精品国产一区二区三区| 美女福利国产在线| 18禁动态无遮挡网站| 久久久久国产网址| 欧美日韩视频精品一区| 亚洲国产av新网站| 99久久综合免费| 如何舔出高潮| 国产亚洲av片在线观看秒播厂| 国产在线男女| av天堂中文字幕网| 亚洲美女黄色视频免费看| 亚洲内射少妇av| 一级毛片我不卡| 亚洲人成网站在线播| 亚洲精华国产精华液的使用体验| 日本猛色少妇xxxxx猛交久久| 看十八女毛片水多多多| 国产精品福利在线免费观看| 午夜久久久在线观看| 国精品久久久久久国模美| 精品亚洲乱码少妇综合久久| 精品人妻偷拍中文字幕| 嫩草影院入口| 亚洲av成人精品一区久久| 亚洲av欧美aⅴ国产| 日韩中文字幕视频在线看片| 亚洲国产最新在线播放| 午夜福利在线观看免费完整高清在| 欧美最新免费一区二区三区| 国产精品偷伦视频观看了| 中国三级夫妇交换| 国语对白做爰xxxⅹ性视频网站| 男男h啪啪无遮挡| 亚洲国产精品一区三区| 五月天丁香电影| 亚洲国产精品国产精品| 国产午夜精品一二区理论片| 性高湖久久久久久久久免费观看| 蜜臀久久99精品久久宅男| 一边亲一边摸免费视频| 中文字幕精品免费在线观看视频 | 少妇人妻一区二区三区视频| 亚洲欧洲精品一区二区精品久久久 | 国产精品久久久久久精品电影小说| 精品亚洲成国产av| 黑丝袜美女国产一区| 精品国产一区二区久久| 国产av国产精品国产| 亚洲av免费高清在线观看| 久久精品久久精品一区二区三区| 女性被躁到高潮视频| 如何舔出高潮| 色5月婷婷丁香| 成人毛片60女人毛片免费| 草草在线视频免费看| 日本爱情动作片www.在线观看| 亚洲欧美日韩另类电影网站| 亚洲精品一区蜜桃| 中国美白少妇内射xxxbb| 一区二区三区精品91| 亚洲,一卡二卡三卡| 男人狂女人下面高潮的视频| 少妇的逼水好多| 亚洲精品日韩在线中文字幕| 一级毛片电影观看| 插阴视频在线观看视频| 熟妇人妻不卡中文字幕| av福利片在线观看| 日本wwww免费看| 人人妻人人添人人爽欧美一区卜| 黄色视频在线播放观看不卡| 欧美精品高潮呻吟av久久| 肉色欧美久久久久久久蜜桃| 大码成人一级视频| av国产久精品久网站免费入址| 久热久热在线精品观看| 欧美精品一区二区大全| 午夜福利视频精品| 伦精品一区二区三区| 七月丁香在线播放| 亚洲精品国产av蜜桃| 亚洲欧美一区二区三区黑人 | 成年女人在线观看亚洲视频| 在线观看三级黄色| 少妇人妻精品综合一区二区| 你懂的网址亚洲精品在线观看| a级毛片免费高清观看在线播放| 色视频www国产| 久久亚洲国产成人精品v| 99视频精品全部免费 在线| 国产精品国产三级专区第一集| 亚洲av.av天堂| 国产精品三级大全| 春色校园在线视频观看| av线在线观看网站| 久久99蜜桃精品久久| 黄色毛片三级朝国网站 | 尾随美女入室| 久久久久久久国产电影| 国产成人一区二区在线| 自拍偷自拍亚洲精品老妇| 国产女主播在线喷水免费视频网站| 日产精品乱码卡一卡2卡三| 久久综合国产亚洲精品| 内地一区二区视频在线| 大话2 男鬼变身卡| 99九九线精品视频在线观看视频| 国产无遮挡羞羞视频在线观看| 亚洲人成网站在线观看播放| 亚洲真实伦在线观看| 青春草亚洲视频在线观看| 超碰97精品在线观看| 国产在线免费精品| 亚洲av男天堂| 99久久综合免费| 老女人水多毛片| 边亲边吃奶的免费视频| 成人特级av手机在线观看| 欧美丝袜亚洲另类| 精品人妻熟女毛片av久久网站| 如日韩欧美国产精品一区二区三区 | 少妇被粗大的猛进出69影院 | 韩国高清视频一区二区三区| 日韩,欧美,国产一区二区三区| 久久久久精品久久久久真实原创| 国模一区二区三区四区视频| 久久久国产欧美日韩av| 边亲边吃奶的免费视频| 麻豆成人av视频| 欧美+日韩+精品| 亚洲美女视频黄频| 性色avwww在线观看| 极品教师在线视频| 亚洲无线观看免费| 精华霜和精华液先用哪个| 热re99久久国产66热| 丝瓜视频免费看黄片| 熟女电影av网| 国产亚洲欧美精品永久| 简卡轻食公司| 亚洲精品一二三| 精品国产乱码久久久久久小说| 在线观看免费高清a一片| 91成人精品电影| 国产免费一级a男人的天堂| 一级片'在线观看视频| 国产 精品1| 久久亚洲国产成人精品v| 国内揄拍国产精品人妻在线| 日本欧美国产在线视频| 日韩成人伦理影院| 日本wwww免费看| 中文乱码字字幕精品一区二区三区| 黄色一级大片看看| 在线观看美女被高潮喷水网站| 黄色欧美视频在线观看| 最近最新中文字幕免费大全7| 在线免费观看不下载黄p国产| 国产精品99久久久久久久久| 观看av在线不卡| 免费观看无遮挡的男女| 日韩精品免费视频一区二区三区 | 免费观看无遮挡的男女| 大又大粗又爽又黄少妇毛片口| 国产极品天堂在线| 亚洲成人手机| 日日啪夜夜爽| 人体艺术视频欧美日本| 欧美性感艳星| 日日摸夜夜添夜夜添av毛片| 男女国产视频网站| 国产精品国产av在线观看| 97超视频在线观看视频| 丝袜喷水一区| 97超视频在线观看视频| videos熟女内射| 日韩中字成人| 丰满饥渴人妻一区二区三| 国产免费一级a男人的天堂| 精品久久久精品久久久| 人体艺术视频欧美日本| 性色avwww在线观看| 成年人免费黄色播放视频 | 狂野欧美激情性xxxx在线观看| 又粗又硬又长又爽又黄的视频| 蜜桃久久精品国产亚洲av| 亚洲精品,欧美精品| 最近2019中文字幕mv第一页| 大话2 男鬼变身卡| 国产老妇伦熟女老妇高清| 美女脱内裤让男人舔精品视频| 91精品伊人久久大香线蕉| 高清视频免费观看一区二区| 狂野欧美激情性xxxx在线观看| 久久婷婷青草| 97在线视频观看| videos熟女内射| 男女边摸边吃奶| 日本爱情动作片www.在线观看| 国产爽快片一区二区三区| 22中文网久久字幕| 一区二区三区四区激情视频| 熟妇人妻不卡中文字幕| 全区人妻精品视频| 午夜久久久在线观看| 亚洲av不卡在线观看| 你懂的网址亚洲精品在线观看| 欧美性感艳星| 久久久久久久大尺度免费视频| 欧美亚洲 丝袜 人妻 在线| 欧美xxxx性猛交bbbb| 男女免费视频国产| 亚洲av免费高清在线观看| 一级毛片黄色毛片免费观看视频| 99久久精品国产国产毛片| 精品国产露脸久久av麻豆| av卡一久久| 九草在线视频观看| 热re99久久国产66热| 成人免费观看视频高清| videossex国产| 精品酒店卫生间| 天美传媒精品一区二区| 高清视频免费观看一区二区| 免费观看性生交大片5| 在线精品无人区一区二区三| 天美传媒精品一区二区| av有码第一页| 午夜91福利影院| 大陆偷拍与自拍| 日本黄色日本黄色录像| 欧美日韩视频高清一区二区三区二| 日日摸夜夜添夜夜添av毛片| 午夜免费鲁丝| 日韩,欧美,国产一区二区三区| 嫩草影院新地址| 国产成人一区二区在线| av国产精品久久久久影院| 亚洲在久久综合| 18禁在线播放成人免费| 日韩在线高清观看一区二区三区| 七月丁香在线播放| 国产日韩一区二区三区精品不卡 | 亚洲不卡免费看| 国产精品一区二区三区四区免费观看| 视频区图区小说| 午夜老司机福利剧场| 另类亚洲欧美激情| 我的女老师完整版在线观看| 韩国高清视频一区二区三区| 欧美区成人在线视频| 丰满迷人的少妇在线观看| 一个人免费看片子| 亚洲天堂av无毛| 国产精品久久久久久久久免| 精品国产乱码久久久久久小说| 午夜精品国产一区二区电影| 亚洲欧美成人综合另类久久久| 黑人巨大精品欧美一区二区蜜桃 | 亚洲真实伦在线观看| 韩国高清视频一区二区三区| 少妇被粗大猛烈的视频| 又爽又黄a免费视频| 亚洲综合色惰| 国产高清三级在线| 国产国拍精品亚洲av在线观看| 成人综合一区亚洲| 久久ye,这里只有精品| 久久亚洲国产成人精品v| 久久人人爽av亚洲精品天堂| 内射极品少妇av片p| 观看美女的网站| 国产毛片在线视频| 91久久精品国产一区二区成人| 九草在线视频观看| 人妻人人澡人人爽人人| 两个人的视频大全免费| 狂野欧美激情性bbbbbb| 如日韩欧美国产精品一区二区三区 | 性色avwww在线观看| 中文字幕人妻丝袜制服| 亚洲国产色片| 国产黄频视频在线观看| 亚洲人成网站在线播| 少妇人妻久久综合中文| 亚洲三级黄色毛片| 成人18禁高潮啪啪吃奶动态图 | 全区人妻精品视频| 人妻少妇偷人精品九色| 蜜臀久久99精品久久宅男| 免费看av在线观看网站| 免费观看性生交大片5| av有码第一页| 插阴视频在线观看视频| 成人毛片a级毛片在线播放| 国产中年淑女户外野战色| 中文欧美无线码| 纯流量卡能插随身wifi吗| 最新中文字幕久久久久| 又黄又爽又刺激的免费视频.| 又粗又硬又长又爽又黄的视频| 久久久精品免费免费高清| 成人漫画全彩无遮挡| 日韩,欧美,国产一区二区三区| 亚洲成人手机| 国产免费一级a男人的天堂| 久久精品夜色国产| 精华霜和精华液先用哪个| 久久精品国产亚洲av天美| 国产成人免费观看mmmm| 青春草视频在线免费观看| 亚洲精品国产成人久久av| 99re6热这里在线精品视频| 卡戴珊不雅视频在线播放| 黑丝袜美女国产一区| 老司机亚洲免费影院| 久久久久久久久久久久大奶| 亚洲精品中文字幕在线视频 | a级毛色黄片| 欧美日韩精品成人综合77777| 国产乱人偷精品视频| 精品午夜福利在线看| 少妇的逼水好多| 精品少妇久久久久久888优播| 看免费成人av毛片| 免费观看的影片在线观看| 日韩av不卡免费在线播放| 男人狂女人下面高潮的视频| 久久久久久伊人网av| 美女中出高潮动态图| 亚洲国产精品一区三区| 波野结衣二区三区在线| 日韩亚洲欧美综合| 蜜臀久久99精品久久宅男| 亚洲国产成人一精品久久久| 免费观看av网站的网址| 日本欧美国产在线视频| 伦理电影大哥的女人| 国产精品欧美亚洲77777| 97在线人人人人妻| 久久热精品热| 午夜福利影视在线免费观看| 欧美激情极品国产一区二区三区 | 青春草视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 中文在线观看免费www的网站| 国产一区亚洲一区在线观看| 另类亚洲欧美激情| 久久韩国三级中文字幕| 日韩欧美 国产精品| 国产美女午夜福利| 纵有疾风起免费观看全集完整版| 在线 av 中文字幕| 久久精品国产亚洲av天美| 日韩熟女老妇一区二区性免费视频| 国产日韩一区二区三区精品不卡 | 国模一区二区三区四区视频| 人妻制服诱惑在线中文字幕| 婷婷色综合www| 国产黄色免费在线视频| 亚洲图色成人| 免费av中文字幕在线| 免费大片黄手机在线观看| 一本久久精品| 亚洲精品视频女| 蜜桃久久精品国产亚洲av| 午夜激情久久久久久久| 欧美+日韩+精品| 人妻人人澡人人爽人人| 蜜臀久久99精品久久宅男| 乱人伦中国视频| 97在线视频观看| 岛国毛片在线播放| 国产淫片久久久久久久久| 日韩人妻高清精品专区| 国产美女午夜福利| 这个男人来自地球电影免费观看 | 男女无遮挡免费网站观看| 亚洲三级黄色毛片| 97超碰精品成人国产| 国产亚洲欧美精品永久| 国产一级毛片在线| 日本av免费视频播放| 97精品久久久久久久久久精品| 寂寞人妻少妇视频99o| 国模一区二区三区四区视频| 人人妻人人澡人人爽人人夜夜| 久久精品久久精品一区二区三区| 肉色欧美久久久久久久蜜桃| 欧美丝袜亚洲另类| 黑人高潮一二区| 欧美+日韩+精品| 中文字幕精品免费在线观看视频 | 色网站视频免费| 久久久久国产精品人妻一区二区| 少妇丰满av| 男人添女人高潮全过程视频| 日本av手机在线免费观看| 欧美日韩一区二区视频在线观看视频在线| 久久人人爽人人片av| 亚洲伊人久久精品综合| kizo精华| 婷婷色综合www| 街头女战士在线观看网站| 日本91视频免费播放| 18禁在线播放成人免费| 亚洲第一区二区三区不卡| a级毛色黄片| 十分钟在线观看高清视频www | 多毛熟女@视频| 三级经典国产精品| 亚洲精品成人av观看孕妇| 色婷婷av一区二区三区视频| 久久毛片免费看一区二区三区| 久久国产亚洲av麻豆专区| 美女大奶头黄色视频| 亚洲美女搞黄在线观看| 日韩强制内射视频| 国产男女超爽视频在线观看| 久久人人爽人人片av| 我的女老师完整版在线观看| 国产伦在线观看视频一区| 大香蕉97超碰在线| 婷婷色综合大香蕉| 99久久精品一区二区三区| 久久 成人 亚洲| 久久热精品热| 97在线人人人人妻| 天堂8中文在线网| 肉色欧美久久久久久久蜜桃| 91精品国产九色| 免费久久久久久久精品成人欧美视频 | 久久 成人 亚洲| 男女啪啪激烈高潮av片| 老熟女久久久| 丰满迷人的少妇在线观看| 久久人妻熟女aⅴ| 特大巨黑吊av在线直播| 多毛熟女@视频| 欧美日韩综合久久久久久| 欧美丝袜亚洲另类| 一边亲一边摸免费视频| 亚洲经典国产精华液单| 99热这里只有是精品在线观看| 日韩三级伦理在线观看| 欧美激情极品国产一区二区三区 | 亚洲精品日本国产第一区| 久久久久精品久久久久真实原创| 国产又色又爽无遮挡免| 99久久精品国产国产毛片| 观看av在线不卡| 少妇的逼水好多| 最近手机中文字幕大全| 深夜a级毛片| av专区在线播放| 嫩草影院入口| 街头女战士在线观看网站| 青春草视频在线免费观看| 精品少妇内射三级| 中文欧美无线码| 成年女人在线观看亚洲视频| 精品国产国语对白av| 国产真实伦视频高清在线观看| 免费少妇av软件| 久久6这里有精品| 97超碰精品成人国产| 国产亚洲5aaaaa淫片| 精品久久久噜噜| 国产成人a∨麻豆精品| 亚洲精品,欧美精品| 亚洲三级黄色毛片| 亚洲婷婷狠狠爱综合网| 亚洲美女黄色视频免费看| 成人美女网站在线观看视频| 亚洲色图综合在线观看| 亚洲精品亚洲一区二区| 在线观看一区二区三区激情| 少妇猛男粗大的猛烈进出视频| 久久99精品国语久久久| 美女中出高潮动态图| 国模一区二区三区四区视频| 国产精品无大码| 成人午夜精彩视频在线观看| 一本久久精品| 久久精品久久久久久噜噜老黄| 尾随美女入室|