• <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在线观看美女高潮| 欧美变态另类bdsm刘玥| 亚洲激情五月婷婷啪啪| 国产精品人妻久久久影院| 好男人视频免费观看在线| 国产精品久久久久久久久免| 亚洲在久久综合| 日韩人妻高清精品专区| av一本久久久久| 我要看黄色一级片免费的| 老熟女久久久| 色94色欧美一区二区| 亚洲熟女精品中文字幕| 久久 成人 亚洲| 免费看日本二区| 亚洲精品一二三| 亚洲欧美一区二区三区黑人 | 精品人妻熟女av久视频| 国产在视频线精品| 久久精品国产亚洲av天美| 国产成人freesex在线| 免费人妻精品一区二区三区视频| 久久人人爽av亚洲精品天堂| 久久精品国产亚洲网站| 秋霞伦理黄片| 亚洲av国产av综合av卡| 男女国产视频网站| 麻豆乱淫一区二区| 国产精品99久久久久久久久| 精品国产乱码久久久久久小说| 在线观看美女被高潮喷水网站| 国产黄色视频一区二区在线观看| 国产av一区二区精品久久| 在线看a的网站| 亚洲一级一片aⅴ在线观看| 精品国产国语对白av| 国产91av在线免费观看| 美女脱内裤让男人舔精品视频| 亚洲欧洲国产日韩| 午夜老司机福利剧场| 高清毛片免费看| 黄色一级大片看看| 极品人妻少妇av视频| 国产成人a∨麻豆精品| 极品教师在线视频| 成人免费观看视频高清| 久久久久久久久大av| 六月丁香七月| 一级爰片在线观看| 制服丝袜香蕉在线| 一边亲一边摸免费视频| 国产日韩欧美视频二区| 日本黄色日本黄色录像| 亚洲美女视频黄频| 日韩欧美精品免费久久| 91成人精品电影| 肉色欧美久久久久久久蜜桃| 2018国产大陆天天弄谢| 久久精品国产亚洲网站| 黄色一级大片看看| 日本vs欧美在线观看视频 | av网站免费在线观看视频| 国产成人精品婷婷| 久久精品国产亚洲网站| 大香蕉久久网| 女人精品久久久久毛片| 麻豆成人午夜福利视频| 国产伦在线观看视频一区| 自拍偷自拍亚洲精品老妇| 国产综合精华液| 亚洲av综合色区一区| 麻豆成人av视频| 国产一区二区三区av在线| 在线观看美女被高潮喷水网站| 成人二区视频| av福利片在线观看| 大片电影免费在线观看免费| 国产日韩欧美视频二区| 97超碰精品成人国产| 成人毛片60女人毛片免费| 久久99蜜桃精品久久| 久久99蜜桃精品久久| 亚洲美女视频黄频| 成人18禁高潮啪啪吃奶动态图 | 99热网站在线观看| 国产男女内射视频| 久久 成人 亚洲| 日韩 亚洲 欧美在线| 建设人人有责人人尽责人人享有的| 香蕉精品网在线| 亚洲国产毛片av蜜桃av| 在线播放无遮挡| 国产男女超爽视频在线观看| 午夜日本视频在线| 亚州av有码| 女性被躁到高潮视频| 亚洲丝袜综合中文字幕| 蜜臀久久99精品久久宅男| 一级二级三级毛片免费看| 国产精品久久久久久久久免| 国产白丝娇喘喷水9色精品| 国产欧美日韩综合在线一区二区 | 亚洲国产精品999| 街头女战士在线观看网站| 亚洲精品久久午夜乱码| 91成人精品电影| 美女脱内裤让男人舔精品视频| 亚洲精品乱码久久久久久按摩| 日本猛色少妇xxxxx猛交久久| 日本av手机在线免费观看| 国产精品一二三区在线看| 国产精品99久久久久久久久| 一本色道久久久久久精品综合| 亚洲av免费高清在线观看| av专区在线播放| 蜜臀久久99精品久久宅男| 久久99热这里只频精品6学生| 18禁动态无遮挡网站| av福利片在线观看| 免费久久久久久久精品成人欧美视频 | 亚洲精品自拍成人| 啦啦啦视频在线资源免费观看| 日日啪夜夜爽| 一级,二级,三级黄色视频| 欧美性感艳星| 91成人精品电影| 纵有疾风起免费观看全集完整版| 色视频在线一区二区三区| 亚洲av不卡在线观看| 最近中文字幕2019免费版| 中文精品一卡2卡3卡4更新| 久久精品夜色国产| 久久久国产精品麻豆| 我的女老师完整版在线观看| 最近的中文字幕免费完整| 99国产精品免费福利视频| 精品国产一区二区久久| 午夜激情福利司机影院| 欧美最新免费一区二区三区| 精品一区在线观看国产| 亚洲精品国产av蜜桃| 精品国产露脸久久av麻豆| 国产亚洲午夜精品一区二区久久| 亚洲国产欧美日韩在线播放 | 日韩 亚洲 欧美在线| 欧美激情国产日韩精品一区| 精品午夜福利在线看| 久久精品国产亚洲网站| 日本免费在线观看一区| 婷婷色麻豆天堂久久| 亚洲情色 制服丝袜| 超碰97精品在线观看| 国产黄片美女视频| 久久久久人妻精品一区果冻| 亚洲人成网站在线播| av.在线天堂| 久久久久人妻精品一区果冻| 久久久国产欧美日韩av| 深夜a级毛片| 好男人视频免费观看在线| 男人添女人高潮全过程视频| 视频区图区小说| 国产精品熟女久久久久浪| 亚洲精品一二三| 一区二区av电影网| 麻豆乱淫一区二区| 亚洲欧美精品自产自拍| 亚洲精品久久午夜乱码| 亚洲丝袜综合中文字幕| 亚洲av成人精品一二三区| 亚洲国产色片| 女人久久www免费人成看片| 丰满人妻一区二区三区视频av| 国产中年淑女户外野战色| 欧美日本中文国产一区发布| 一级a做视频免费观看| 丰满饥渴人妻一区二区三| 国产黄片美女视频| 搡老乐熟女国产| 九九爱精品视频在线观看| 天天操日日干夜夜撸| 久久青草综合色| 老女人水多毛片| 久久久午夜欧美精品| 91久久精品国产一区二区成人| 在线观看一区二区三区激情| 欧美日韩精品成人综合77777| av不卡在线播放| 久久av网站| 亚洲伊人久久精品综合| 久久精品久久久久久久性| 国产男女超爽视频在线观看| 久久人人爽av亚洲精品天堂| 免费黄频网站在线观看国产| 亚洲精品久久久久久婷婷小说| 久久韩国三级中文字幕| 精品一区二区三卡| 国产成人91sexporn| 51国产日韩欧美| h日本视频在线播放| 2018国产大陆天天弄谢| 新久久久久国产一级毛片| 亚洲av中文av极速乱| 久久午夜福利片| 18+在线观看网站| 极品人妻少妇av视频| 亚洲精品成人av观看孕妇| 久久久久久久久久久免费av| 免费观看a级毛片全部| 97超视频在线观看视频| 久久人妻熟女aⅴ| 国产精品国产三级专区第一集| 日韩中文字幕视频在线看片| 久久精品久久久久久久性| 日本免费在线观看一区| 久久韩国三级中文字幕| 亚洲久久久国产精品| 日韩三级伦理在线观看| 91精品国产九色| 国产精品一区二区性色av| 校园人妻丝袜中文字幕| av视频免费观看在线观看| 最新中文字幕久久久久| 99热这里只有是精品在线观看| 大陆偷拍与自拍| 亚洲国产欧美在线一区| 啦啦啦中文免费视频观看日本| 日本vs欧美在线观看视频 | 精品亚洲成a人片在线观看| 日本黄色日本黄色录像| 精品久久久久久久久亚洲| 国产成人免费无遮挡视频| 久久精品国产a三级三级三级| 国产伦精品一区二区三区四那| 久久久国产精品麻豆| 日韩伦理黄色片| 天堂8中文在线网| 一级av片app| 91精品一卡2卡3卡4卡| 久久99精品国语久久久| 视频中文字幕在线观看| 久久久久久久大尺度免费视频| 欧美激情极品国产一区二区三区 | av女优亚洲男人天堂| 亚洲人成网站在线播| 午夜免费男女啪啪视频观看| 3wmmmm亚洲av在线观看| 久久精品夜色国产| 国产女主播在线喷水免费视频网站| 91精品伊人久久大香线蕉| 国产在线男女| 啦啦啦在线观看免费高清www| 亚洲国产成人一精品久久久| 久久99精品国语久久久| 亚洲四区av| 久久鲁丝午夜福利片| 精品久久国产蜜桃| 欧美精品高潮呻吟av久久| 国产欧美日韩综合在线一区二区 | 午夜福利网站1000一区二区三区| 卡戴珊不雅视频在线播放| 亚洲人成网站在线观看播放| 一级a做视频免费观看| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩视频精品一区| a级毛色黄片| 大话2 男鬼变身卡| 久久国内精品自在自线图片| 色网站视频免费| 免费观看性生交大片5| 国产午夜精品一二区理论片| 夫妻午夜视频| 一级毛片aaaaaa免费看小| 国产精品秋霞免费鲁丝片| 国产视频内射| 男人舔奶头视频| 成人无遮挡网站| 国产精品99久久久久久久久| 亚洲情色 制服丝袜| 在线观看免费视频网站a站| 一级毛片 在线播放| 丝袜在线中文字幕| 黄片无遮挡物在线观看| 人妻系列 视频| 男女免费视频国产| 性色avwww在线观看| 黄色一级大片看看| 日日爽夜夜爽网站| 免费看光身美女| 午夜久久久在线观看| 观看免费一级毛片| 一级毛片 在线播放| 狠狠精品人妻久久久久久综合| 最新的欧美精品一区二区| 热99国产精品久久久久久7| 夫妻性生交免费视频一级片| 热re99久久国产66热| 亚洲,欧美,日韩| 亚洲精品国产成人久久av| 大又大粗又爽又黄少妇毛片口| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久久久免| 日本色播在线视频| 亚洲国产欧美日韩在线播放 | 热re99久久精品国产66热6| 久久99一区二区三区| 国产免费一区二区三区四区乱码| 久久女婷五月综合色啪小说| 嫩草影院入口| 国产精品一区二区在线观看99| 黄色怎么调成土黄色| 欧美老熟妇乱子伦牲交| 日韩强制内射视频| 国产成人午夜福利电影在线观看| av一本久久久久| 日韩精品免费视频一区二区三区 | 男女免费视频国产| 色婷婷久久久亚洲欧美| 永久免费av网站大全| 一本久久精品| 欧美最新免费一区二区三区| videossex国产| 看十八女毛片水多多多| 免费观看无遮挡的男女| 久久久久久久久久久免费av| 美女内射精品一级片tv| 色94色欧美一区二区| 国精品久久久久久国模美| 午夜福利,免费看| 国产精品一区二区在线不卡| 日韩在线高清观看一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 伦理电影大哥的女人| 国产综合精华液| 亚洲国产精品999| 成人影院久久| 91精品一卡2卡3卡4卡| 成人特级av手机在线观看| 国产欧美日韩精品一区二区| 多毛熟女@视频| 中文字幕久久专区| 国产亚洲午夜精品一区二区久久| h视频一区二区三区| 在线观看三级黄色| h视频一区二区三区| 久久人人爽人人爽人人片va| 日韩一区二区视频免费看| 女人精品久久久久毛片| 99热这里只有是精品在线观看| 天堂中文最新版在线下载| 精品久久久噜噜| 22中文网久久字幕| 午夜免费鲁丝| 乱码一卡2卡4卡精品| 天天躁夜夜躁狠狠久久av| 啦啦啦视频在线资源免费观看| 精品久久久噜噜| 国产一区二区在线观看日韩| 夜夜骑夜夜射夜夜干| 国产又色又爽无遮挡免| 亚洲怡红院男人天堂| 麻豆成人午夜福利视频| 看免费成人av毛片| av免费在线看不卡| 性色av一级| av.在线天堂| 日本爱情动作片www.在线观看| 国产精品麻豆人妻色哟哟久久| 只有这里有精品99| 亚洲欧美清纯卡通| 日韩成人伦理影院| 亚洲伊人久久精品综合| a级毛色黄片| 三级国产精品欧美在线观看| 91在线精品国自产拍蜜月| 五月天丁香电影| 大香蕉97超碰在线| 在线观看免费视频网站a站| 26uuu在线亚洲综合色| 在现免费观看毛片| 亚洲性久久影院| 国产男女内射视频| 国产爽快片一区二区三区| 22中文网久久字幕| 美女内射精品一级片tv| 成人美女网站在线观看视频| 亚洲精品日韩av片在线观看| 五月天丁香电影| 亚洲成色77777| 18禁在线无遮挡免费观看视频| 久久久久久久久久人人人人人人| 中文字幕精品免费在线观看视频 | 在线观看免费视频网站a站| 97在线视频观看| 免费人成在线观看视频色| 亚洲内射少妇av| 国产午夜精品一二区理论片| 日韩成人av中文字幕在线观看| 精品少妇黑人巨大在线播放| 久久青草综合色| 国产在线男女| 成人黄色视频免费在线看| 天天躁夜夜躁狠狠久久av| 亚洲精品乱久久久久久| 成人二区视频| 国产极品天堂在线| 亚洲四区av| 成年人免费黄色播放视频 | 青青草视频在线视频观看| 久久精品国产亚洲av涩爱| 久久 成人 亚洲| 国产男女超爽视频在线观看| 色婷婷久久久亚洲欧美| 汤姆久久久久久久影院中文字幕| 亚洲成人一二三区av| 少妇丰满av| 久久精品国产亚洲av天美| 我的女老师完整版在线观看| 秋霞伦理黄片| 日本猛色少妇xxxxx猛交久久| 亚洲精品一区蜜桃| 七月丁香在线播放| 在线观看免费视频网站a站| 丰满人妻一区二区三区视频av| 中文字幕制服av| 欧美亚洲 丝袜 人妻 在线| 国产成人a∨麻豆精品| 久久 成人 亚洲| 国产又色又爽无遮挡免| 女人精品久久久久毛片| 九九爱精品视频在线观看| 国产精品成人在线| 国产精品一区www在线观看| 久久国产精品大桥未久av | 99久国产av精品国产电影| 一个人看视频在线观看www免费| 在线观看三级黄色| 亚洲精品乱码久久久v下载方式| 久久99蜜桃精品久久| 国产黄色视频一区二区在线观看| 免费观看在线日韩| 国产成人a∨麻豆精品| 99re6热这里在线精品视频| 97在线视频观看| 欧美一级a爱片免费观看看| 亚洲欧洲国产日韩| 香蕉精品网在线| 在线观看免费日韩欧美大片 | 免费在线观看成人毛片| 在线看a的网站| 最近手机中文字幕大全| 成人免费观看视频高清| av在线播放精品| 中国美白少妇内射xxxbb| 伦理电影免费视频| 国产伦精品一区二区三区四那| √禁漫天堂资源中文www| 午夜影院在线不卡| 最近最新中文字幕免费大全7| 亚洲一区二区三区欧美精品| 另类亚洲欧美激情| 99久久精品热视频| 国产视频内射| 只有这里有精品99| a级毛色黄片| 免费人妻精品一区二区三区视频| 亚洲精品日本国产第一区| 男女免费视频国产| 久久久久精品性色| 欧美精品人与动牲交sv欧美| 日日爽夜夜爽网站| 久久久久国产网址| 在线观看免费高清a一片| 女人久久www免费人成看片| 久久韩国三级中文字幕| 欧美精品高潮呻吟av久久| 亚洲欧美成人综合另类久久久| av线在线观看网站| 丰满人妻一区二区三区视频av| 老女人水多毛片| 久久人人爽人人爽人人片va| 99热这里只有是精品50| 国产精品人妻久久久影院| 久久国产精品男人的天堂亚洲 | 99re6热这里在线精品视频| 久久午夜综合久久蜜桃| 在线亚洲精品国产二区图片欧美 | 亚洲天堂av无毛| 男人狂女人下面高潮的视频| 汤姆久久久久久久影院中文字幕| 欧美3d第一页| 国产精品福利在线免费观看| 妹子高潮喷水视频| 日本av免费视频播放| 啦啦啦啦在线视频资源| 免费观看a级毛片全部| 免费在线观看成人毛片| 亚洲精品亚洲一区二区| 少妇熟女欧美另类| 亚洲成人av在线免费| 性色avwww在线观看| 亚洲av日韩在线播放| 国产黄色视频一区二区在线观看| 亚洲激情五月婷婷啪啪| 99热6这里只有精品| 国产精品偷伦视频观看了| 久久久久国产精品人妻一区二区| 欧美区成人在线视频| 美女内射精品一级片tv| 色网站视频免费| 我要看黄色一级片免费的| 寂寞人妻少妇视频99o| 国产美女午夜福利| 一区二区三区精品91| 亚洲国产欧美日韩在线播放 | 亚洲精品乱码久久久v下载方式| 国产精品成人在线| 另类精品久久| 日韩中字成人| 国产91av在线免费观看| 七月丁香在线播放| 亚洲三级黄色毛片| 欧美另类一区| 免费看不卡的av| 国产精品女同一区二区软件| 久久毛片免费看一区二区三区| 大片电影免费在线观看免费| 在线观看免费视频网站a站| 黑人猛操日本美女一级片| 久热这里只有精品99| 天美传媒精品一区二区| 欧美激情极品国产一区二区三区 | 草草在线视频免费看| 亚洲国产精品国产精品| 3wmmmm亚洲av在线观看| 99久久人妻综合| 69精品国产乱码久久久| 十八禁高潮呻吟视频 | 日本欧美视频一区| 男女边吃奶边做爰视频| 国产探花极品一区二区| 男男h啪啪无遮挡| 欧美最新免费一区二区三区| 嫩草影院新地址| .国产精品久久| 男男h啪啪无遮挡| 中文字幕人妻丝袜制服| 亚洲精品日本国产第一区| 22中文网久久字幕| 国产中年淑女户外野战色| 哪个播放器可以免费观看大片| 国产深夜福利视频在线观看| 777米奇影视久久| 久久av网站| 嫩草影院入口| 交换朋友夫妻互换小说| 建设人人有责人人尽责人人享有的| 搡女人真爽免费视频火全软件| 国产精品麻豆人妻色哟哟久久| 多毛熟女@视频| 久久鲁丝午夜福利片| 我要看黄色一级片免费的| 久久久午夜欧美精品| 久久久久久久国产电影| 各种免费的搞黄视频| 青春草视频在线免费观看| 亚洲精品国产色婷婷电影| 少妇被粗大的猛进出69影院 | 久久这里有精品视频免费| 亚洲欧洲精品一区二区精品久久久 | 69精品国产乱码久久久| 久久精品国产a三级三级三级| 爱豆传媒免费全集在线观看| 日韩电影二区| 国产精品女同一区二区软件| 精品国产一区二区三区久久久樱花| 99久国产av精品国产电影| freevideosex欧美| 在线亚洲精品国产二区图片欧美 | 成人漫画全彩无遮挡| 国产成人精品福利久久| 99久久精品国产国产毛片| 婷婷色麻豆天堂久久| 熟妇人妻不卡中文字幕| 丰满人妻一区二区三区视频av| 国产午夜精品一二区理论片| 天美传媒精品一区二区| 99re6热这里在线精品视频| 亚洲,一卡二卡三卡| 高清午夜精品一区二区三区| 国产精品一区www在线观看| 日产精品乱码卡一卡2卡三| 青春草国产在线视频| 精品一区二区免费观看| 一级毛片aaaaaa免费看小| 嫩草影院新地址| 肉色欧美久久久久久久蜜桃| 最近2019中文字幕mv第一页| 少妇人妻一区二区三区视频| 乱码一卡2卡4卡精品| 国产视频首页在线观看| 久久久久久久久久久免费av| 男人爽女人下面视频在线观看| 日日啪夜夜撸| 丝袜在线中文字幕| 精品午夜福利在线看| 久久精品国产亚洲av涩爱| 国内精品宾馆在线| 国产精品偷伦视频观看了| 在线观看人妻少妇|