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

    徑向總壓畸變對(duì)風(fēng)扇前后流場(chǎng)的影響研究

    2017-06-05 14:16:08高緒萬吳宇豪
    關(guān)鍵詞:總壓畸變風(fēng)扇

    鄧 甜,雷 攀,高緒萬,吳宇豪

    (中國(guó)民航大學(xué)中歐航空工程師學(xué)院,天津 300300)

    徑向總壓畸變對(duì)風(fēng)扇前后流場(chǎng)的影響研究

    鄧 甜,雷 攀,高緒萬,吳宇豪

    (中國(guó)民航大學(xué)中歐航空工程師學(xué)院,天津 300300)

    采用數(shù)值模擬的方法,對(duì)某大涵道比小型渦扇發(fā)動(dòng)機(jī)的風(fēng)扇部分進(jìn)行建模,并分別設(shè)置等效低壓區(qū)為扇形角度60°、90°、120°及150°這4種不同程度的總壓畸變進(jìn)口條件,研究進(jìn)氣道徑向總壓畸變進(jìn)口條件下,風(fēng)扇前后流場(chǎng)的變換規(guī)律。研究發(fā)現(xiàn),風(fēng)扇后流場(chǎng)的總壓畸變指數(shù)比風(fēng)扇前大,說明風(fēng)扇后流場(chǎng)受到風(fēng)扇轉(zhuǎn)動(dòng)作用,其不穩(wěn)定度程度比風(fēng)扇前大;相對(duì)于切向速度和軸向速度,徑向速度畸變指數(shù)受總壓畸變的影響更大;總壓畸變指數(shù)在風(fēng)扇前呈增長(zhǎng)趨勢(shì),且不同徑向總壓畸變進(jìn)口條件的變化規(guī)律十分相似,總壓畸變指數(shù)在風(fēng)扇后存在最大值區(qū)域。

    渦扇發(fā)動(dòng)機(jī);徑向總壓畸變;畸變指數(shù);數(shù)值模擬

    徑向總壓畸變屬于進(jìn)口總壓畸變中的一種,其存在能在一定程度上影響發(fā)動(dòng)機(jī)的穩(wěn)定性,近現(xiàn)代由于發(fā)動(dòng)機(jī)轉(zhuǎn)速的提升,使得總壓畸變的影響程度也相應(yīng)增加到不得不考慮的范圍,國(guó)內(nèi)外關(guān)于總壓畸變的研究如下。

    WAWalter等[1]基于壓力和溫度方面的考慮對(duì)F100型發(fā)動(dòng)機(jī)進(jìn)行了實(shí)驗(yàn)研究。Awatef Hamed[2]和Milt Davis等[3]為發(fā)動(dòng)機(jī)地面測(cè)試提供了方案,主要對(duì)畸變發(fā)生位置進(jìn)行了深入研究。文獻(xiàn)[4]中采用三維非穩(wěn)態(tài)流場(chǎng)進(jìn)行計(jì)算,對(duì)壓氣機(jī)進(jìn)氣道畸變的響應(yīng)進(jìn)行了分析。William T Cousins等[5-6]利用實(shí)際實(shí)驗(yàn),通過安裝畸變擾流板人為地創(chuàng)造畸變條件,并對(duì)某型高涵道比渦扇發(fā)動(dòng)機(jī)進(jìn)行相關(guān)實(shí)驗(yàn)研究。Jose Rodrigue等[7]為改善進(jìn)口流場(chǎng)的穩(wěn)定度,針對(duì)畸變條件專門設(shè)計(jì)了新型壓氣機(jī)進(jìn)口。

    張靖煊[8]于2007年研究了進(jìn)氣道總壓畸變對(duì)簡(jiǎn)單葉型的壓氣機(jī)葉頂間隙流的穩(wěn)定性,并提出了相應(yīng)的改善方法。孫鵬[9]通過數(shù)值模擬,對(duì)跨聲速風(fēng)扇流場(chǎng)結(jié)構(gòu)對(duì)畸變的響應(yīng)進(jìn)行了研究。鄭寧等[10]從三維模型出發(fā),進(jìn)行了風(fēng)扇進(jìn)氣畸變的技術(shù)研究。張環(huán)等[11]進(jìn)行了帶旋轉(zhuǎn)的總壓畸變進(jìn)口條件對(duì)壓氣機(jī)穩(wěn)定性影響的研究。2010年任偉峰[12]針對(duì)小流量的渦噴發(fā)動(dòng)機(jī),對(duì)總壓畸變的評(píng)定技術(shù)進(jìn)行了分析。王春利等[13]進(jìn)行了進(jìn)氣畸變對(duì)兩級(jí)風(fēng)扇穩(wěn)定性影響的數(shù)值模擬研究。朱愛迪等[14]通過總壓變化圖譜,基于畸變擾流板進(jìn)行實(shí)際實(shí)驗(yàn),以觀察畸變流場(chǎng)的變化規(guī)律。此外,劉大響等[15-19]一直在進(jìn)行包括總壓畸變測(cè)評(píng)技術(shù)、發(fā)動(dòng)機(jī)關(guān)于總壓畸變的穩(wěn)定性測(cè)評(píng)方法、總壓畸變實(shí)驗(yàn)方法等的相關(guān)研究和總結(jié)。

    以上國(guó)內(nèi)外的研究大多只針對(duì)風(fēng)扇前流場(chǎng),而對(duì)風(fēng)扇后流場(chǎng)對(duì)畸變響應(yīng)的考慮十分少。特別是,針對(duì)具體某型號(hào)發(fā)動(dòng)機(jī),尤其當(dāng)涉及到不同類型和大小的壓氣機(jī)時(shí),此類型畸變研究的參考價(jià)值更是十分有限。

    本文基于進(jìn)氣道徑向總壓畸變進(jìn)口條件,通過完整的風(fēng)扇流場(chǎng)模型,包括進(jìn)氣道和延展的風(fēng)扇后流場(chǎng),對(duì)某型渦扇發(fā)動(dòng)機(jī)DGEN380進(jìn)行進(jìn)氣道徑向總壓畸變對(duì)其風(fēng)扇前后流場(chǎng)影響的研究。本文選取60°、90°、120°和150°這4種等效低壓區(qū)面積的局部進(jìn)口低壓區(qū),低壓區(qū)總壓設(shè)定依據(jù)總壓差值與均值之比的畸變指數(shù)計(jì)算方法,初始條件參照南京航空航天大學(xué)的進(jìn)口總壓畸變穩(wěn)定性評(píng)定方法研究[20],以25%畸變指數(shù)換算得到?;冎笖?shù)定義如下

    其中:DC(pt)為總壓畸變指數(shù);pt.max為截面最大總壓;pt.min為截面最小總壓;pt.av為截面總壓均值。

    1 數(shù)值模擬模型及計(jì)算條件

    本文運(yùn)用NUMECA軟件,利用某小型渦扇發(fā)動(dòng)機(jī)采集風(fēng)扇設(shè)計(jì)參數(shù)進(jìn)行葉型設(shè)計(jì)和風(fēng)扇模型建立。模型的建立基于跟DGEN380配套的數(shù)值模擬平臺(tái)內(nèi)風(fēng)扇相關(guān)的數(shù)據(jù),借助于NUMECA軟件的AUTOBLADE模塊進(jìn)行葉型建立,利用IGG模塊進(jìn)行網(wǎng)格劃分,通過Fine模塊進(jìn)行數(shù)值計(jì)算,最后借助于CFVIEW模塊進(jìn)行后處理和數(shù)據(jù)分析,再利用Matlab進(jìn)行數(shù)據(jù)處理。整個(gè)模型將具有14個(gè)葉片的風(fēng)扇置于流道的前中部,進(jìn)口和出口截面都距風(fēng)扇有一定距離,這樣便于觀察總壓畸變進(jìn)口的風(fēng)扇前后流場(chǎng)的變化趨勢(shì)。

    如圖1所示,模型流道長(zhǎng)30 cm,并設(shè)定流場(chǎng)進(jìn)口為軸向原點(diǎn),出口為軸向距離0.3 m,風(fēng)扇葉片弦長(zhǎng)0.08 m,進(jìn)口處位于軸向0.045 m處,出口處位于軸向0.125 m處。為較好觀察風(fēng)扇出口流場(chǎng)各物理量隨總壓畸變的變化趨勢(shì),設(shè)計(jì)選取了較長(zhǎng)的流場(chǎng)出口區(qū)域。

    全風(fēng)扇共有葉片14片,其結(jié)構(gòu)與網(wǎng)格劃分如圖2所示。

    圖1 風(fēng)扇單通道模型及網(wǎng)格Fig.1 Single channel model and grid

    圖2 全風(fēng)扇流道3D和網(wǎng)格化Fig.2 3D grid of fan flow

    數(shù)值模擬首先進(jìn)行了無畸變地面工況下的風(fēng)扇部件穩(wěn)態(tài)流場(chǎng)計(jì)算,總壓為1 bar,總溫為288 K,轉(zhuǎn)速為13 150 r/m。隨后分別進(jìn)行不同進(jìn)口區(qū)遮擋面積周向總壓畸變的非穩(wěn)態(tài)流場(chǎng)計(jì)算。

    2 徑向總壓畸變的影響

    徑向總壓畸變是進(jìn)氣道總壓畸變研究中常見的一種總壓畸變形式。其低壓區(qū)的定義為周向全覆蓋,徑向由葉尖向葉根部位進(jìn)行覆蓋,其覆蓋面積的變化代表不同程度的徑向總壓畸變。

    本文選取與60°~150°這4種低壓區(qū)覆蓋面積不同的徑向總壓畸變作為進(jìn)口條件。整個(gè)模型的進(jìn)口葉根部位半徑為0.046 m,葉尖部位為0.175 m,為簡(jiǎn)化模型,建模時(shí)未考慮葉頂間隙的存在。經(jīng)換算得到其低壓區(qū)等效面積由葉尖向葉根部位覆蓋的半徑如表1所示。

    表1 徑向總壓畸變進(jìn)口低壓區(qū)半徑覆蓋位置Tab.1 Radius coverage of total radial pressure distortion in inlet low pressure zone

    除了進(jìn)口總壓低壓區(qū)覆蓋角度變化外,其他計(jì)算條件均不作任何改變,且與無畸變時(shí)一樣。以等效低壓區(qū)覆蓋面積為60°和150°扇形面積的徑向總壓畸變?yōu)槔?,進(jìn)口條件為60°和150°徑向總壓畸變進(jìn)口總壓分布云圖和數(shù)值圖如圖3和圖4所示。

    圖3 徑向60°進(jìn)氣道總壓畸變進(jìn)口Fig.3 Radial 60°total inlet pressure distortion

    圖4 徑向150°進(jìn)氣道總壓畸變進(jìn)口Fig.4 Radial 150°total inlet pressure distortion

    圖3(a)和圖4(a)分別為60°和150°徑向總壓畸變進(jìn)口條件的進(jìn)口總壓分布情況,淺色區(qū)域?yàn)榈蛪簠^(qū),總壓數(shù)值為77 660 Pa,深色區(qū)域?yàn)檎4髿鈮簭?qiáng)101 300 Pa;圖3(b)和圖4(b)為對(duì)應(yīng)的沿徑向的總壓數(shù)值變化圖。

    2.1 風(fēng)扇前總壓畸變指數(shù)的變化趨勢(shì)

    圖5(a)和圖5(b)分別為60°和150°等效低壓區(qū)面積的徑向總壓畸變進(jìn)口條件下風(fēng)扇前0.03 m處的總壓云圖。從圖中可以看出,對(duì)于徑向總壓畸變進(jìn)口條件,其總壓數(shù)值會(huì)梯度化沿徑向分布,周向呈規(guī)律性變化。

    圖5 徑向總壓畸變進(jìn)口條件下風(fēng)扇前0.03 m處總壓云圖Fig.5 Total pressure contour of total radial pressure distortion inlet at 0.03 m

    對(duì)于徑向總壓畸變來說,其風(fēng)扇前流場(chǎng)總壓畸變指數(shù)的計(jì)算基于上面的公式,可得到如圖6所示的畸變指數(shù)隨軸向位置變化的數(shù)值圖,RD代表radial distortion,括號(hào)內(nèi)不同數(shù)字代表不同的等效低壓區(qū)面積。

    徑向總壓畸變的低壓區(qū)采用周向全覆蓋,徑向取不同數(shù)值大小的方式。從圖6中也可以看出,風(fēng)扇前的徑向總壓畸變指數(shù)從60°~150°都呈現(xiàn)穩(wěn)定增長(zhǎng)的趨勢(shì),且隨著畸變低壓區(qū)面積的增大而呈現(xiàn)上升趨勢(shì)。

    2.2 風(fēng)扇后總壓畸變指數(shù)的變化趨勢(shì)

    圖7(a)~圖7(d)分別為60°徑向總壓畸變進(jìn)口條件下,風(fēng)扇后流場(chǎng)軸向位置為0.15 m、0.2 m、0.25 m和0.29 m處的總壓云圖,圖8(a)~圖8(d)分別為150°徑向總壓畸變進(jìn)口條件下,風(fēng)扇后流場(chǎng)軸向位置為0.15m、0.2m、0.25 m和0.29 m處的總壓云圖。

    圖6 風(fēng)扇進(jìn)口總壓畸變指數(shù)隨軸向位置的數(shù)值圖Fig.6 Numerical figure of fan inlet total pressure distortion index with axial position

    圖7 60°徑向總壓畸變進(jìn)口條件下的風(fēng)扇后流場(chǎng)總壓云圖Fig.7 Total pressure contour of fan with 60°total radial pressure distortion

    圖8 150°徑向總壓畸變進(jìn)口條件下的風(fēng)扇后流場(chǎng)總壓云圖Fig.8 Total pressure contour of fan with 150°total radial pressure distortion

    從圖8中可以看出,風(fēng)扇后流場(chǎng)受風(fēng)扇轉(zhuǎn)動(dòng)作用的影響還很大,由于葉片的帶動(dòng)作用,總壓分布也會(huì)呈周期性分布,但這種分布會(huì)隨著軸向位置遠(yuǎn)離葉片而慢慢消失。此外,由于靠近葉尖區(qū)域存在低壓區(qū),故在葉尖區(qū)域其總壓數(shù)值沒有明顯的周期性變化規(guī)律。

    相對(duì)于風(fēng)扇前總壓畸變指數(shù)呈現(xiàn)出穩(wěn)定增長(zhǎng)的趨勢(shì),風(fēng)扇后流場(chǎng)的總壓畸變指數(shù)則呈現(xiàn)出無規(guī)則的變化趨勢(shì),其畸變指數(shù)隨軸向位置的數(shù)值變化如圖9所示。

    從圖9可以看出,靠近風(fēng)扇葉片部分的徑向總壓畸變指數(shù)受葉片附面層影響而相對(duì)較大,之后的正常區(qū)域呈先增大后減小的趨勢(shì),且在0.2 m處附近會(huì)取得最大值,之后慢慢趨于穩(wěn)定。值得注意的是,在趨于穩(wěn)定的時(shí)候,其不同程度的徑向總壓畸變指數(shù)的穩(wěn)定值會(huì)隨著畸變區(qū)面積的減小而增大。

    圖9 風(fēng)扇出口總壓畸變指數(shù)隨軸向位置的數(shù)值圖Fig.9 Numerical figure of fan outlet total pressure distortion index with axial position

    為進(jìn)一步對(duì)比徑向總壓畸變對(duì)風(fēng)扇前后流場(chǎng)的影響,將風(fēng)扇前后徑向總壓畸變指數(shù)隨軸向位置變化的數(shù)值圖放在一幅圖中,可得到如圖10所示的數(shù)值變化圖。

    圖10 徑向總壓畸變進(jìn)口條件下風(fēng)扇前后畸變指數(shù)的數(shù)值圖Fig.10 Numerical figure of distortion index of forward and backward fans with radial total pressure distortion

    從圖10中可以看出,風(fēng)扇出口的徑向總壓畸變指數(shù)相對(duì)于風(fēng)扇進(jìn)口整體上會(huì)更大。由此說明,在經(jīng)過風(fēng)扇轉(zhuǎn)動(dòng)作用后的總壓畸變指數(shù)會(huì)變得更大,即徑向總壓畸變進(jìn)口條件對(duì)風(fēng)扇前流場(chǎng)固然有一定的影響,然而風(fēng)扇前的不穩(wěn)定流場(chǎng)在經(jīng)過風(fēng)扇轉(zhuǎn)動(dòng)作用后會(huì)變得更加不穩(wěn)定。因此,在進(jìn)行進(jìn)氣道總壓畸變相關(guān)研究中,鑒于風(fēng)扇后流場(chǎng)相對(duì)于風(fēng)扇前更加不穩(wěn)定,且風(fēng)扇后流場(chǎng)直接會(huì)影響發(fā)動(dòng)機(jī)后續(xù)部件的工作效率,更應(yīng)該進(jìn)行進(jìn)氣道總壓畸變對(duì)風(fēng)扇后流場(chǎng)的相關(guān)研究。

    2.3 進(jìn)口總壓畸變對(duì)風(fēng)扇后速度場(chǎng)的影響

    速度畸變計(jì)算公式類似于總壓畸變指數(shù),采用某截面速度最大值與最小值的差值與均值之比,能夠得到如圖11所示的徑向總壓畸變進(jìn)口條件下的風(fēng)扇后流場(chǎng)徑向速度畸變指數(shù)隨軸向位置的數(shù)值變化圖。其中,不同符號(hào)標(biāo)記的曲線代表不同程度的徑向總壓畸變進(jìn)口條件,如RSD(60)(RSD,radial speed distortion),代表等效遮擋面積為60°的徑向總壓畸變進(jìn)口條件下的風(fēng)扇后徑向速度畸變指數(shù)隨軸向位置的變化。

    圖11 徑向速度畸變指數(shù)隨軸向位置的數(shù)值圖Fig.11 Numerical figure of radial velocity distortion index with axial position

    從圖11中可以看出,徑向速度畸變指數(shù)隨軸向位置的變化趨勢(shì)。先維持穩(wěn)定,后呈現(xiàn)先增大后減小的趨勢(shì),且基于不同低壓區(qū)等效面積的徑向總壓畸變進(jìn)口條件,其徑向總壓畸變指數(shù)在變化過程中均在點(diǎn)0.22 m處附近取得最大值。除此之外,不同低壓區(qū)等效面積的徑向總壓畸變指數(shù)會(huì)隨著低壓區(qū)面積的增大而呈現(xiàn)減小的趨勢(shì)。

    基于同樣的速度畸變指數(shù)計(jì)算方法,可以得到如圖12和圖13所示的切向速度畸變指數(shù)和軸向速度畸變指數(shù)隨軸向位置變化的數(shù)值圖,CSD代表circumferential speed distortion,ASD代表 axial speed distortion。相對(duì)于徑向速度畸變指數(shù),切向速度畸變指數(shù)和軸向速度畸變指數(shù)相對(duì)來說要小一些。

    從圖12中可以看出,切向速度畸變指數(shù)在60°和90°兩種情況在靠近風(fēng)扇部分會(huì)有小幅度上升趨勢(shì),隨后一直處于下降趨勢(shì),120°和150°兩種情況則在風(fēng)扇后一直處于下降趨勢(shì),并慢慢趨于穩(wěn)定,且不同程度的徑向總壓畸變進(jìn)口條件,切向速度畸變指數(shù)會(huì)隨著等效遮擋面積的增大而減小。從圖13中的軸向速度畸變指數(shù)可以看出,4種不同程度的徑向總壓畸變進(jìn)口條件下的風(fēng)扇后流場(chǎng)的軸向速度畸變指數(shù)會(huì)在靠近風(fēng)扇部分呈明顯下降趨勢(shì),隨后在0.17 m處有小幅上升,并逐漸趨于穩(wěn)定。比較不同的徑向總壓畸變進(jìn)口條件,軸向速度會(huì)隨著等效遮擋面積的增加而上升,這一點(diǎn)與切向速度畸變指數(shù)正好相反。

    為進(jìn)一步比較相同進(jìn)口條件下各方向上速度畸變指數(shù)的相對(duì)大小,將同一種進(jìn)口條件下的3種不同的速度畸變指數(shù)放到一幅圖中,圖14和圖15分別為低壓區(qū)等效面積為60°和150°扇形面積進(jìn)口徑向總壓畸變條件下3種速度的速度畸變指數(shù)數(shù)值圖,其中RD(60)為radial distortion,等效面積為60°扇形面積時(shí)的徑向總壓畸變進(jìn)口條件。

    對(duì)于這兩種不同的徑向總壓畸變進(jìn)口條件,比較其不同方向上的速度畸變指數(shù)發(fā)現(xiàn),徑向速度畸變指數(shù)要明顯大于切向速度畸變指數(shù)和軸向速度畸變指數(shù),且切向速度畸變指數(shù)和軸向速度畸變指數(shù)基本維持相同的數(shù)值范圍,具有相似的變化趨勢(shì)。與周向總壓畸變進(jìn)口不同的是,徑向速度畸變指數(shù)不會(huì)一直處于增長(zhǎng)狀態(tài),在經(jīng)過一段距離的增長(zhǎng)后,在0.22 m(也就是總壓畸變指數(shù)開始趨于穩(wěn)定的點(diǎn))到0.26 m區(qū)間趨于平穩(wěn),在0.26 m后,徑向速度畸變指數(shù)總體上又呈下降趨勢(shì)。

    圖12 切向速度畸變指數(shù)隨軸向位置的數(shù)值圖Fig.12 Numerical figure of tangential velocity distortion index with axial position

    圖13 軸向速度畸變指數(shù)隨軸向位置的數(shù)值圖Fig.13 Numerical figure of axial velocity distortion index with axial position

    圖14 60°徑向總壓畸變進(jìn)口條件下的速度畸變指數(shù)數(shù)值圖Fig.14 Velocity distortion exponent of 60°total radial pressure distortion

    圖15 150°徑向總壓畸變進(jìn)口條件下的速度畸變指數(shù)數(shù)值圖Fig.15 Velocity distortion exponent of 150°total radial pressure distortion

    3 結(jié)語(yǔ)

    1)徑向總壓畸變進(jìn)口條件下的風(fēng)扇后流場(chǎng)的總壓畸變指數(shù)比風(fēng)扇前流場(chǎng)的總壓畸變指數(shù)大,總的來說發(fā)動(dòng)機(jī)風(fēng)扇前的畸變流場(chǎng)受風(fēng)扇轉(zhuǎn)動(dòng)作用影響到達(dá)風(fēng)扇后會(huì)有更高的不穩(wěn)定度。

    2)徑向總壓畸變進(jìn)口條件下風(fēng)扇前總壓畸變指數(shù)的變化在不同程度的進(jìn)口條件下,其變化規(guī)律十分類似;風(fēng)扇出口正常區(qū)域流場(chǎng)的總壓畸變指數(shù)整體變化先增大后減小,但不同程度的進(jìn)口條件,增長(zhǎng)規(guī)律不同。

    3)徑向總壓畸變進(jìn)口條件對(duì)速度場(chǎng)有一定的影響,且徑向速度畸變指數(shù)明顯大于軸向速度畸變指數(shù)和切向速度畸變指數(shù),但3個(gè)方向的速度畸變指數(shù)都會(huì)穩(wěn)定在一定范圍。

    總的來說,本文對(duì)于此型號(hào)以及類似的大涵道比小型渦扇發(fā)動(dòng)機(jī)的總壓畸變數(shù)值模擬研究工作,一方面提供了風(fēng)扇前后流場(chǎng)的畸變程度的有效對(duì)比,另一方面對(duì)此型號(hào)發(fā)動(dòng)機(jī)即將進(jìn)行的試車臺(tái)總壓畸變實(shí)驗(yàn)提供了參考數(shù)據(jù)。此外,對(duì)此型號(hào)發(fā)動(dòng)機(jī)的適航工作,特別是針對(duì)進(jìn)氣道總壓畸變程度的測(cè)試,無論是測(cè)試方法還是需要考慮的畸變可能發(fā)生的位置等都具有一定的參考意義。

    [1]WALTER W,SHAW M.Predicted F100 Engine Response to Circumferential Pressure and Temperature Distortion[C]//Awatef Hamed.Inlet Distortion Considerations for High Cycle Fatigue in Gas Turbine Engine,AIAA,1997.

    [2]HAMED A,NUMBERS K,HAMEDA,et al.Inlet Distortion ConsiderationsforHighCycleFatigueinGas Turbine Engine[C]//Joint Propulsion Conference and Exhibit,2013.

    [3]DAVIS M,HALE A,BEALE D.An argument for enhancement of the current inlet distortion ground test practice for aircraft gas turbine engines[J].Journal of Turbomachinery,2002,124(2):V001T01A011.

    [4]AIAA.Unsteady Three-Dimensional Analysis of Inlet Distortion in Turbomachinery[C]//Joint Propulsion Conference and Exhibit,1997:35-38.

    [5]COUSINS W T.History,Philosophy,Physics,and Future Directions of Aircraft Propulsion System/Inlet Integration[C]//ASME Turbo Expo 2004:Power for Land,Sea and Air,2004:305-320.

    [6]COUSINS W T,GEORGES M J,REZAEI H.Inlet Distortion Testing andAnalysisofaHigh-BypassRatioTurbofanEngine[C]//ISABE,2003.

    [7]RODRIGUEZ J,KLUMPP S,BIESINGER T,et al.A New Inlet Distortion and Pressure Loss Based Design of an Intake System for Stationary Gas Turbines[C]//ASME Turbo Expo 2013:Turbine Technical Conference and Exposition,2013:V06BT38A002.

    [8]張靖煊.畸變條件下軸流壓氣機(jī)葉頂間隙流對(duì)流動(dòng)失穩(wěn)影響的非定常機(jī)制[D].北京:中國(guó)科學(xué)院研究生院(工程熱物理研究所),2007.

    [9]孫 鵬.總壓畸變對(duì)跨聲速風(fēng)扇流場(chǎng)結(jié)構(gòu)影響的全流道數(shù)值模擬[D].哈爾濱:哈爾濱工業(yè)大學(xué),2007.

    [10]鄭 寧,鄒正平,徐力平.風(fēng)扇進(jìn)氣畸變?nèi)S非定常數(shù)值模擬技術(shù)研究[J].航空動(dòng)力學(xué)報(bào),2007,22(1):60-65.

    [11]張 環(huán).旋轉(zhuǎn)總壓畸變對(duì)壓氣機(jī)穩(wěn)定性影響的研究[D].南京:南京航空航天大學(xué),2009.

    [12]任偉峰.小流量渦噴發(fā)動(dòng)機(jī)進(jìn)口總壓畸變?cè)u(píng)定技術(shù)研究[D].長(zhǎng)沙:湖南大學(xué),2010.

    [13]王春利,喬渭陽(yáng).畸變進(jìn)氣對(duì)兩級(jí)風(fēng)扇穩(wěn)定性影響的數(shù)值模擬[J].燃?xì)鉁u輪試驗(yàn)與研究,2014(4):1-6.

    [14]朱愛迪,鐘易成.基于進(jìn)氣道總壓畸變?cè)囼?yàn)的插板數(shù)值模擬研究[J].航空動(dòng)力學(xué)報(bào),2015(3):143-145,157.

    [15]劉大響,葉培梁,胡 駿,等.航空燃?xì)鉁u輪發(fā)動(dòng)機(jī)穩(wěn)定性設(shè)計(jì)與評(píng)定技術(shù)[M].北京:北京航空工業(yè)出版社,2004:186-278.

    [16]李 陽(yáng),胡 駿.航空發(fā)動(dòng)機(jī)進(jìn)氣總壓周向畸變數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào),2005,31(2):11-13.

    [17]黃順洲,胡 駿.進(jìn)氣畸變對(duì)發(fā)動(dòng)機(jī)穩(wěn)定性影響的分析模型[J].航空動(dòng)力學(xué)報(bào),2006,27(5):426-430.

    [18]王志強(qiáng),胡 駿,李 亮,等.進(jìn)口總壓隨機(jī)脈動(dòng)的數(shù)值模擬方法[J].航空動(dòng)力學(xué)報(bào),2010,25(1):87-91.

    [19]胡 駿,趙運(yùn)生,丁 寧,等.進(jìn)氣畸變對(duì)大涵道比渦扇發(fā)動(dòng)機(jī)穩(wěn)定性的影響[J].航空動(dòng)力學(xué)報(bào),2013,39(6):7-12.

    [20]李 亮.進(jìn)口總壓畸變穩(wěn)定性評(píng)定方法研究[D].南京:南京航空航天大學(xué),2008.

    (責(zé)任編輯:楊媛媛)

    Research on influence of radial total pressure distortion on forward and backward fan flow

    DENG Tian,LEI Pan,GAO Xuwan,WU Yuhao
    (Sino-European Institute of Aviation Engineering,CAUC,Tianjin 300300,China)

    Modelling and simulation are used to study a small turbofan engine with high bypass ratio.With four situations of inlet total pressure distortion of equivalent surface of sector angles 60°,90°,120°and 150°,the influence of inlet radial total pressure distortion on forward and backward fan flow of turbofan engine is studied.It is found that the total pressure distortion index is bigger for backward fan part,and it is proved that the unstable level of backward fan part is bigger than forward fan part because of fan rotation;compared with tangential speed and axial speed, the radial speed distortion support is influenced more by inlet radial total pressure distortion;for the flow field of forward fan part,the total pressure distortion has the same type of increase for the four different inlet situations, but it is not the same for backward fan part which has the biggest value area.

    turbofan-engine;radial total pressure distortion;distortion index;modelling and simulation

    V231.3

    A

    1674-5590(2017)02-0024-07

    2016-09-09;

    2016-10-16

    國(guó)家自然科學(xué)基金項(xiàng)目(51506216);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)(ZXH2012H001)

    鄧甜(1982—),女,天津人,講師,博士,研究方向?yàn)橥牧鳌上嗔?

    猜你喜歡
    總壓畸變風(fēng)扇
    總壓探針性能結(jié)構(gòu)敏感性分析
    可調(diào)式總壓耙設(shè)計(jì)及應(yīng)用
    亞聲速條件下總壓探針臨壁效應(yīng)的數(shù)值研究
    2 m超聲速風(fēng)洞流場(chǎng)變速壓控制方法研究
    在Lightroom中校正鏡頭與透視畸變
    電風(fēng)扇
    故事大王(2017年4期)2017-05-08 07:53:40
    基于智能手機(jī)控制風(fēng)扇運(yùn)行的實(shí)現(xiàn)
    新蒙迪歐車?yán)鋮s風(fēng)扇常高速運(yùn)轉(zhuǎn)
    輻射誘導(dǎo)染色體畸變的快速FISH方法的建立
    奇瑞A5車散熱風(fēng)扇無低速擋
    大又大粗又爽又黄少妇毛片口| 嫩草影院精品99| 尾随美女入室| 国产aⅴ精品一区二区三区波| 久久99热6这里只有精品| 在线免费观看不下载黄p国产| 成年免费大片在线观看| 欧美极品一区二区三区四区| 久久久国产成人精品二区| 老女人水多毛片| 免费无遮挡裸体视频| 日本与韩国留学比较| 一级毛片电影观看 | 美女大奶头视频| 国产成人福利小说| 1024手机看黄色片| 99精品在免费线老司机午夜| 亚洲一区二区三区色噜噜| 国产精品一及| 在现免费观看毛片| 成人一区二区视频在线观看| 欧美+亚洲+日韩+国产| 亚洲熟妇熟女久久| 国产精品美女特级片免费视频播放器| 赤兔流量卡办理| av在线蜜桃| 国产成人91sexporn| 岛国在线免费视频观看| 亚洲不卡免费看| 熟妇人妻久久中文字幕3abv| 国产一区二区三区av在线 | 狂野欧美白嫩少妇大欣赏| 国产探花在线观看一区二区| 国产一级毛片七仙女欲春2| 我的女老师完整版在线观看| 一级毛片我不卡| 欧美国产日韩亚洲一区| 欧美三级亚洲精品| 天堂√8在线中文| 亚洲四区av| 十八禁国产超污无遮挡网站| 人妻夜夜爽99麻豆av| 国产成年人精品一区二区| 久久久久九九精品影院| 国产高清三级在线| 成人av一区二区三区在线看| 国产精品爽爽va在线观看网站| 国产黄色视频一区二区在线观看 | 亚洲人成网站在线播| 国产探花极品一区二区| 久久亚洲精品不卡| 99久久精品国产国产毛片| 一本精品99久久精品77| 国产精品一区二区性色av| 五月玫瑰六月丁香| 精品人妻熟女av久视频| 久久久久免费精品人妻一区二区| 18禁黄网站禁片免费观看直播| 综合色丁香网| 日本三级黄在线观看| 亚洲av五月六月丁香网| 免费电影在线观看免费观看| 欧美色视频一区免费| 欧美一级a爱片免费观看看| 亚洲第一区二区三区不卡| 美女黄网站色视频| 人妻少妇偷人精品九色| 国产伦一二天堂av在线观看| 成人国产麻豆网| 少妇猛男粗大的猛烈进出视频 | 久久精品国产亚洲av涩爱 | 日本爱情动作片www.在线观看 | 97超视频在线观看视频| 午夜免费激情av| 一本精品99久久精品77| 一进一出抽搐gif免费好疼| 国产国拍精品亚洲av在线观看| 精品欧美国产一区二区三| 午夜精品一区二区三区免费看| 永久网站在线| 18禁在线播放成人免费| 久久精品国产鲁丝片午夜精品| 人人妻人人澡人人爽人人夜夜 | 干丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 久久99热这里只有精品18| 99国产精品一区二区蜜桃av| 男女那种视频在线观看| 日韩一本色道免费dvd| 九九热线精品视视频播放| 久久久久性生活片| 婷婷亚洲欧美| 亚洲18禁久久av| 麻豆成人午夜福利视频| 在线国产一区二区在线| 精华霜和精华液先用哪个| 欧美激情在线99| 亚洲不卡免费看| 九九在线视频观看精品| 久久久精品94久久精品| 精品日产1卡2卡| 免费无遮挡裸体视频| 99热网站在线观看| 熟女人妻精品中文字幕| 亚洲最大成人中文| 少妇丰满av| 嫩草影院入口| 天堂网av新在线| 中文字幕av成人在线电影| 久久久久久大精品| 99久国产av精品| 看十八女毛片水多多多| 床上黄色一级片| 亚洲人与动物交配视频| 国产精品国产高清国产av| 婷婷六月久久综合丁香| 国产人妻一区二区三区在| 亚洲成av人片在线播放无| 村上凉子中文字幕在线| 岛国在线免费视频观看| 成人一区二区视频在线观看| 五月伊人婷婷丁香| 老司机影院成人| 天堂av国产一区二区熟女人妻| 国产男靠女视频免费网站| 搡老妇女老女人老熟妇| 亚洲人成网站高清观看| 一区二区三区高清视频在线| 最近中文字幕高清免费大全6| 欧美高清性xxxxhd video| 亚洲专区国产一区二区| a级毛片免费高清观看在线播放| 97超级碰碰碰精品色视频在线观看| 熟妇人妻久久中文字幕3abv| 菩萨蛮人人尽说江南好唐韦庄 | 国产av不卡久久| 免费在线观看成人毛片| 一个人看的www免费观看视频| av免费在线看不卡| 性色avwww在线观看| 中文字幕免费在线视频6| 亚洲av成人精品一区久久| 国产成人一区二区在线| 一本久久中文字幕| 夜夜夜夜夜久久久久| 男女那种视频在线观看| 欧美日韩国产亚洲二区| 免费无遮挡裸体视频| 日日啪夜夜撸| 国产av不卡久久| 亚洲va在线va天堂va国产| 国产91av在线免费观看| 色5月婷婷丁香| 欧美+亚洲+日韩+国产| 国产成年人精品一区二区| 中国美女看黄片| 91精品国产九色| 国产亚洲精品av在线| 午夜福利高清视频| 亚洲精品色激情综合| 国产精品久久视频播放| 午夜精品国产一区二区电影 | 中出人妻视频一区二区| 两个人的视频大全免费| 亚洲欧美日韩卡通动漫| 亚洲av中文字字幕乱码综合| 婷婷六月久久综合丁香| 日本撒尿小便嘘嘘汇集6| 看十八女毛片水多多多| 中文字幕熟女人妻在线| 丝袜喷水一区| 日日啪夜夜撸| 禁无遮挡网站| 熟妇人妻久久中文字幕3abv| 免费电影在线观看免费观看| 国产精品久久久久久久久免| 国内少妇人妻偷人精品xxx网站| 一夜夜www| 亚洲天堂国产精品一区在线| 亚洲国产精品sss在线观看| 亚洲va在线va天堂va国产| 中文字幕av在线有码专区| 精品久久久久久久久亚洲| 男人的好看免费观看在线视频| 波多野结衣巨乳人妻| 国产不卡一卡二| 精品欧美国产一区二区三| 亚洲,欧美,日韩| 久久久久久久久大av| 国产精品美女特级片免费视频播放器| 国产高清不卡午夜福利| 偷拍熟女少妇极品色| 久久久久久久久久久丰满| ponron亚洲| 亚洲av第一区精品v没综合| 真实男女啪啪啪动态图| www日本黄色视频网| 久久国内精品自在自线图片| 精品少妇黑人巨大在线播放 | 美女免费视频网站| 国产精品美女特级片免费视频播放器| 九九在线视频观看精品| 日本爱情动作片www.在线观看 | av在线天堂中文字幕| 直男gayav资源| 国产乱人视频| 俺也久久电影网| 久久亚洲国产成人精品v| 精品一区二区免费观看| 欧美色欧美亚洲另类二区| 波野结衣二区三区在线| 欧美一区二区亚洲| 十八禁国产超污无遮挡网站| 国产一区二区在线观看日韩| 可以在线观看的亚洲视频| 婷婷亚洲欧美| 美女xxoo啪啪120秒动态图| 国产又黄又爽又无遮挡在线| 精华霜和精华液先用哪个| 精品欧美国产一区二区三| 国产精品乱码一区二三区的特点| 精品乱码久久久久久99久播| 久久久久九九精品影院| 欧美最新免费一区二区三区| 亚洲精品乱码久久久v下载方式| 最近2019中文字幕mv第一页| 久久久a久久爽久久v久久| 99久久九九国产精品国产免费| 美女免费视频网站| 69人妻影院| 九九爱精品视频在线观看| 亚洲av二区三区四区| 久久亚洲精品不卡| 欧美成人a在线观看| 亚洲不卡免费看| 一级黄色大片毛片| 国产私拍福利视频在线观看| 久久九九热精品免费| 免费人成在线观看视频色| 亚洲精品亚洲一区二区| 美女xxoo啪啪120秒动态图| 久久精品国产99精品国产亚洲性色| 亚洲性久久影院| 国产极品精品免费视频能看的| 婷婷色综合大香蕉| 丰满人妻一区二区三区视频av| 亚洲欧美日韩卡通动漫| 99久久精品国产国产毛片| 午夜福利成人在线免费观看| 最近最新中文字幕大全电影3| 伊人久久精品亚洲午夜| 午夜福利高清视频| www.色视频.com| 免费看美女性在线毛片视频| 午夜激情福利司机影院| 亚洲人成网站高清观看| 三级经典国产精品| 99国产极品粉嫩在线观看| 99热只有精品国产| 亚洲国产欧美人成| 久久久久国内视频| 欧美成人精品欧美一级黄| 成人一区二区视频在线观看| 日韩,欧美,国产一区二区三区 | 国产成人91sexporn| 99九九线精品视频在线观看视频| 18禁在线无遮挡免费观看视频 | 欧美日韩一区二区视频在线观看视频在线 | 色5月婷婷丁香| 人妻夜夜爽99麻豆av| 久久人人精品亚洲av| АⅤ资源中文在线天堂| 久久韩国三级中文字幕| 夜夜夜夜夜久久久久| 免费观看精品视频网站| 女同久久另类99精品国产91| 成人三级黄色视频| 露出奶头的视频| 国产黄色视频一区二区在线观看 | 久久亚洲精品不卡| 日韩欧美免费精品| 国产精品美女特级片免费视频播放器| 久久精品国产鲁丝片午夜精品| 色5月婷婷丁香| 三级经典国产精品| 国产精品无大码| 欧美性感艳星| 欧美xxxx性猛交bbbb| 久久久精品大字幕| 少妇高潮的动态图| 日本与韩国留学比较| 久久久久久久午夜电影| 在现免费观看毛片| 国产 一区 欧美 日韩| 日韩国内少妇激情av| 亚洲精品在线观看二区| 国产精品一区二区三区四区久久| 国产精品久久视频播放| 国产又黄又爽又无遮挡在线| 内射极品少妇av片p| 欧美一区二区精品小视频在线| 日韩欧美三级三区| 性插视频无遮挡在线免费观看| 69人妻影院| 热99re8久久精品国产| 午夜激情福利司机影院| 亚洲不卡免费看| 内射极品少妇av片p| 神马国产精品三级电影在线观看| 激情 狠狠 欧美| 秋霞在线观看毛片| 大又大粗又爽又黄少妇毛片口| 亚洲av五月六月丁香网| 婷婷六月久久综合丁香| а√天堂www在线а√下载| 久久午夜亚洲精品久久| 成年av动漫网址| 一个人看视频在线观看www免费| 最好的美女福利视频网| 成人漫画全彩无遮挡| 免费看av在线观看网站| 赤兔流量卡办理| 日韩精品青青久久久久久| 波多野结衣高清作品| 久久国内精品自在自线图片| 欧美一级a爱片免费观看看| 国产精品一区二区性色av| 久久久久久久午夜电影| 亚洲人与动物交配视频| 麻豆一二三区av精品| 一区二区三区高清视频在线| 天堂动漫精品| 亚洲四区av| 国产探花在线观看一区二区| 日本黄大片高清| 女生性感内裤真人,穿戴方法视频| 春色校园在线视频观看| 最近最新中文字幕大全电影3| 高清毛片免费看| 18+在线观看网站| 91午夜精品亚洲一区二区三区| 天堂√8在线中文| 成人av在线播放网站| 免费观看精品视频网站| 麻豆久久精品国产亚洲av| 亚洲精品色激情综合| 毛片女人毛片| 男女那种视频在线观看| 成人特级黄色片久久久久久久| 综合色av麻豆| 秋霞在线观看毛片| 麻豆一二三区av精品| 中文字幕av成人在线电影| 亚洲国产欧美人成| 给我免费播放毛片高清在线观看| 一本精品99久久精品77| 伦理电影大哥的女人| 一a级毛片在线观看| 岛国在线免费视频观看| 日本黄色片子视频| 欧美激情国产日韩精品一区| av在线观看视频网站免费| 亚洲av一区综合| 久久草成人影院| 人人妻人人澡人人爽人人夜夜 | 级片在线观看| 99久久久亚洲精品蜜臀av| 舔av片在线| 男人狂女人下面高潮的视频| 国产伦在线观看视频一区| 欧美日韩在线观看h| 美女xxoo啪啪120秒动态图| 亚洲国产精品成人久久小说 | 欧美另类亚洲清纯唯美| 日韩国内少妇激情av| 精华霜和精华液先用哪个| 成人午夜高清在线视频| 精品午夜福利视频在线观看一区| 大香蕉久久网| 久久亚洲精品不卡| 日韩高清综合在线| 观看免费一级毛片| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久久久黄片| 国产精华一区二区三区| 国产探花在线观看一区二区| 免费在线观看影片大全网站| 久久久久久国产a免费观看| 三级毛片av免费| 国产成人一区二区在线| 成人三级黄色视频| 国产高清有码在线观看视频| 九九爱精品视频在线观看| 日本成人三级电影网站| 国产高清视频在线观看网站| 亚洲国产精品成人久久小说 | 国产在视频线在精品| 高清午夜精品一区二区三区 | 国产一区亚洲一区在线观看| 国产麻豆成人av免费视频| 精品少妇黑人巨大在线播放 | 99热这里只有精品一区| 你懂的网址亚洲精品在线观看 | 可以在线观看毛片的网站| 此物有八面人人有两片| 亚洲av一区综合| 亚洲久久久久久中文字幕| 国产熟女欧美一区二区| 成年女人永久免费观看视频| 夜夜夜夜夜久久久久| 欧美不卡视频在线免费观看| 国产一区亚洲一区在线观看| 三级毛片av免费| 国产精品乱码一区二三区的特点| 国产精品久久久久久精品电影| 男人和女人高潮做爰伦理| 国产激情偷乱视频一区二区| 少妇人妻一区二区三区视频| 午夜福利在线观看吧| 99精品在免费线老司机午夜| 亚洲精品日韩av片在线观看| 成人无遮挡网站| 啦啦啦观看免费观看视频高清| 亚洲乱码一区二区免费版| 亚洲婷婷狠狠爱综合网| 成人一区二区视频在线观看| 亚洲精品一区av在线观看| 两个人视频免费观看高清| 又爽又黄a免费视频| 中文字幕av成人在线电影| 一级毛片久久久久久久久女| 午夜久久久久精精品| 亚洲精品日韩在线中文字幕 | 国产精品女同一区二区软件| 欧美中文日本在线观看视频| 九九在线视频观看精品| 亚洲精品粉嫩美女一区| 久久精品国产99精品国产亚洲性色| h日本视频在线播放| 国产黄a三级三级三级人| 性欧美人与动物交配| av福利片在线观看| 色哟哟·www| 女人被狂操c到高潮| 女同久久另类99精品国产91| 精品欧美国产一区二区三| 欧美另类亚洲清纯唯美| 永久网站在线| 欧美性感艳星| 乱码一卡2卡4卡精品| 嫩草影视91久久| 女人十人毛片免费观看3o分钟| 看黄色毛片网站| 一进一出好大好爽视频| 国产伦精品一区二区三区视频9| 久久午夜福利片| 天堂av国产一区二区熟女人妻| 成人国产麻豆网| 久久久国产成人免费| 特级一级黄色大片| 深夜精品福利| 欧美极品一区二区三区四区| 国产爱豆传媒在线观看| 日韩av在线大香蕉| 国产精品人妻久久久久久| 久久久成人免费电影| 少妇熟女欧美另类| 成人特级黄色片久久久久久久| 免费观看人在逋| 欧美潮喷喷水| 久久精品久久久久久噜噜老黄 | 99视频精品全部免费 在线| 久久久久久大精品| 国内久久婷婷六月综合欲色啪| 蜜臀久久99精品久久宅男| 春色校园在线视频观看| 久久亚洲国产成人精品v| 搡女人真爽免费视频火全软件 | 亚洲精品乱码久久久v下载方式| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品成人综合色| 中文字幕久久专区| 成人国产麻豆网| 三级毛片av免费| 国产熟女欧美一区二区| 日韩欧美免费精品| 淫妇啪啪啪对白视频| 免费看a级黄色片| 给我免费播放毛片高清在线观看| 中国美白少妇内射xxxbb| 国产高清视频在线观看网站| 色av中文字幕| 亚洲真实伦在线观看| 一本久久中文字幕| 精品少妇黑人巨大在线播放 | 国产 一区 欧美 日韩| 色在线成人网| 久久久国产成人免费| 免费无遮挡裸体视频| 免费大片18禁| 国产欧美日韩精品一区二区| 久久精品国产亚洲av香蕉五月| 热99re8久久精品国产| 国产高清不卡午夜福利| 深爱激情五月婷婷| 小蜜桃在线观看免费完整版高清| 久久久久久伊人网av| 国产精品一二三区在线看| 夜夜看夜夜爽夜夜摸| 欧美zozozo另类| 免费av观看视频| 国产精品国产三级国产av玫瑰| 老司机福利观看| 此物有八面人人有两片| 亚洲精品国产av成人精品 | 十八禁网站免费在线| 人妻久久中文字幕网| 日韩一本色道免费dvd| 日本一二三区视频观看| 成人高潮视频无遮挡免费网站| 综合色av麻豆| 国产精品一区二区免费欧美| av国产免费在线观看| 少妇熟女欧美另类| 我的老师免费观看完整版| av卡一久久| 国产精品亚洲美女久久久| 久久久久久国产a免费观看| 久久久久久久久久久丰满| 国产精品无大码| 秋霞在线观看毛片| 我的老师免费观看完整版| 亚洲婷婷狠狠爱综合网| 精品少妇黑人巨大在线播放 | 久久久久久久亚洲中文字幕| 村上凉子中文字幕在线| 国产淫片久久久久久久久| 国产亚洲精品久久久久久毛片| 少妇裸体淫交视频免费看高清| 免费一级毛片在线播放高清视频| 99热网站在线观看| 亚洲欧美精品综合久久99| 男人舔奶头视频| 给我免费播放毛片高清在线观看| 欧美最新免费一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 成人永久免费在线观看视频| 综合色av麻豆| 99久久精品一区二区三区| 五月玫瑰六月丁香| 又爽又黄a免费视频| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利成人在线免费观看| 亚洲精品日韩av片在线观看| 国产视频内射| 亚洲精品在线观看二区| 亚洲人与动物交配视频| 国产高潮美女av| 欧美色视频一区免费| 亚洲中文字幕日韩| 亚洲va在线va天堂va国产| 国产精品1区2区在线观看.| 99久久无色码亚洲精品果冻| 欧美激情久久久久久爽电影| 精品免费久久久久久久清纯| 久久久久国产精品人妻aⅴ院| 国产黄色小视频在线观看| 91狼人影院| 国产精品精品国产色婷婷| 国产精品人妻久久久久久| 97在线视频观看| 国产欧美日韩精品一区二区| 亚洲精品粉嫩美女一区| 黄色欧美视频在线观看| 亚洲av免费在线观看| 一本精品99久久精品77| 97超视频在线观看视频| 男女那种视频在线观看| 亚洲一级一片aⅴ在线观看| av黄色大香蕉| 国产精品久久久久久精品电影| 尤物成人国产欧美一区二区三区| 99热网站在线观看| 尤物成人国产欧美一区二区三区| 国产蜜桃级精品一区二区三区| 久久精品国产鲁丝片午夜精品| 国产伦精品一区二区三区视频9| 国产乱人偷精品视频| 亚洲av美国av| 亚洲av成人精品一区久久| 深夜精品福利| 精品国内亚洲2022精品成人| 少妇的逼好多水| 午夜视频国产福利| 久久久精品94久久精品| 人人妻人人看人人澡| 色综合站精品国产| 天堂av国产一区二区熟女人妻| 欧美一区二区国产精品久久精品| 成人精品一区二区免费| 最近视频中文字幕2019在线8| av在线亚洲专区| 男女之事视频高清在线观看| 久久精品国产亚洲av涩爱 | 99久久中文字幕三级久久日本| 国产又黄又爽又无遮挡在线| 变态另类丝袜制服| 精品久久久久久久人妻蜜臀av| 欧美性猛交黑人性爽| 亚洲成av人片在线播放无| 俄罗斯特黄特色一大片| 国产精品亚洲美女久久久| 一夜夜www| 伦理电影大哥的女人| 一区二区三区免费毛片|