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

    欠頻響壓力探針測量壓氣機動態(tài)流場的結(jié)果分析

    2016-09-23 03:37:58馬宏偉
    航空發(fā)動機 2016年2期
    關(guān)鍵詞:頻響總壓壓氣機

    馬宏偉,魏 巍,張 良,官 輝

    (1.北京航空航天大學(xué)能源與動力工程學(xué)院;2.先進(jìn)航空發(fā)動機協(xié)同創(chuàng)新中心:北京100191)

    欠頻響壓力探針測量壓氣機動態(tài)流場的結(jié)果分析

    馬宏偉1,2,魏巍1,2,張良1,2,官輝1,2

    (1.北京航空航天大學(xué)能源與動力工程學(xué)院;2.先進(jìn)航空發(fā)動機協(xié)同創(chuàng)新中心:北京100191)

    壓氣機內(nèi)部流動具有很強的非定常特征,尤其是小尺寸發(fā)動機葉片通過頻率高達(dá)30~80 kH z,而動態(tài)壓力探針的響應(yīng)頻率只能與葉片通過頻率相當(dāng),采用其測量不易獲得具體的高頻流場結(jié)構(gòu)。為了研究這種欠頻響情況下測得流場與真實流場的差異,以某低速壓氣機為對象進(jìn)行了模擬研究。采用頻響為420 H z的動態(tài)壓力探針及頻響大于33 kH z的單斜絲熱線探針,測量了葉片通過頻率為303 H z的轉(zhuǎn)子出口流場。修正了動態(tài)壓力探針的容腔效應(yīng)對測得流場的影響,校驗了在測試系統(tǒng)響應(yīng)頻率較低時,修正后流場與熱線測得的參考流場之間的差異。對比結(jié)果表明:原始流場壓力信號的部分頻率成分產(chǎn)生了較大偏差,出口截面流場的周向位置有較大誤差。經(jīng)修正后,周向位置較為準(zhǔn)確,泄漏堵塞區(qū)域形態(tài)以及強度都與參考流場的更為接近,但由于尾跡區(qū)存在部分高頻,修正后結(jié)果與參考流場的差異仍較明顯。

    欠頻響壓力探針;動態(tài)測量;管傳遞函數(shù);動態(tài)修正;壓氣機;航空發(fā)動機

    0 引言

    準(zhǔn)確測量動態(tài)壓力信號對描述航空發(fā)動機壓氣機部件中的非定常流動特征起著非常重要的作用,但相對于壓氣機內(nèi)部流場極高頻的動態(tài)特征,高頻壓力探針的頻響往往不能滿足測量要求。目前,常規(guī)尺寸的航空發(fā)動機壓氣機部件的典型葉片通過頻率為5~20 kHz,小尺寸發(fā)動機的葉片通過頻率則高達(dá)30~80 kHz,而多數(shù)高頻壓力探針的頻響只有30~50 kHz[1-8]??梢哉J(rèn)為在多數(shù)情況下,這些動態(tài)壓力探針都是在欠頻響即接近流場特征頻率的狀態(tài)下工作。

    因此,為了擴大測量系統(tǒng)的頻響范圍、提高動態(tài)測量的準(zhǔn)確度,必須進(jìn)行仔細(xì)的動態(tài)標(biāo)定和修正[9]。此前,很多研究者曾給出過展寬測量系統(tǒng)頻響或動態(tài)修正的方法[10-13]。李繼超等[10]標(biāo)定了容腔效應(yīng)對高頻動態(tài)探針的影響,通過改進(jìn)容腔尺寸來提高高頻動態(tài)探針的頻響,從而滿足試驗需求;楊林等[11]采用半無限引壓管測量動態(tài)壓力,并進(jìn)行了動態(tài)標(biāo)定,有效地減小動態(tài)壓力信號的能量衰減;郭勇等[12]利用動態(tài)補償數(shù)字濾波器技術(shù),展寬了高頻壓力探針的工作頻帶。

    綜上所述,現(xiàn)有的研究集中在探針的結(jié)構(gòu)設(shè)計、動態(tài)修正等方面,以提高動態(tài)壓力探針的頻響,但在動態(tài)探針頻響不能遠(yuǎn)大于流場特征頻率情況下,少有文獻(xiàn)剖析測得的流場形態(tài)是否準(zhǔn)確。

    在真實發(fā)動機試驗環(huán)境中,現(xiàn)有的試驗技術(shù)很難提供標(biāo)準(zhǔn)的動態(tài)流場來與動態(tài)探針測得的流場進(jìn)行對比。本文以某低速壓氣機試驗設(shè)備為對象進(jìn)行了模擬研究。采用頻響接近流場特征頻率的動態(tài)壓力探針,研究了測量流場與參考流場的差異,并探討了該情況下旋轉(zhuǎn)探針法的適用性。

    1 試驗系統(tǒng)介紹

    1.1測試設(shè)備及探針動態(tài)標(biāo)定

    本文采用圓柱狀動態(tài)壓力探針,側(cè)面的取壓孔直徑為2.2 mm,探針直徑為4 mm,引壓管長400 mm,管徑為2.2 mm。傳感器選用SMI5882硅微壓傳感器,其穩(wěn)態(tài)標(biāo)定精度為0.1%。由于做了封裝處理,傳感器的標(biāo)稱頻響為500 Hz,由溫漂和重復(fù)性造成的誤差限小于8 Pa。

    另外,參考動態(tài)流場采用單斜絲熱線探針進(jìn)行測量,其熱絲直徑為5 μm,方波測試結(jié)果顯示整套熱線系統(tǒng)的動態(tài)頻響大于33 kHz。試驗中使用1套NI-PXIe高速數(shù)據(jù)采集系統(tǒng)進(jìn)行數(shù)據(jù)采集,所有的動態(tài)采集的采樣率均為120 kHz,采樣點數(shù)為600 kHz。

    動態(tài)標(biāo)定設(shè)備如圖1所示。由旋轉(zhuǎn)盤提供周期性的氣流激勵,旋轉(zhuǎn)速度可調(diào);轉(zhuǎn)盤在邊緣沿圓周方向均布了120個直徑為2 mm的圓孔,改變單脈沖信號的頻率可通過改變圓盤轉(zhuǎn)速來實現(xiàn)。

    圖1 探針動態(tài)標(biāo)定裝置

    動態(tài)標(biāo)定的主要目的是標(biāo)定該段引壓管的容腔效應(yīng)對動態(tài)響應(yīng)的影響。在標(biāo)定過程中,裸傳感器和探針頭部分別被脈沖氣流激勵,然后根據(jù)測得數(shù)據(jù)建立管傳遞函數(shù)。經(jīng)FFT分析,激勵信號為5 ms的脈沖氣流激勵包含0~500 Hz的頻率成分。

    管傳遞函數(shù)的建立采取2種方法。第1種是經(jīng)典的FFT算法,即直接對裸傳感器得到的“標(biāo)準(zhǔn)信號”和探針測到的信號進(jìn)行FFT分析,將頻域中輸出量和輸入量幅值的比值和相位差作為此總壓探針的管傳遞函數(shù),文獻(xiàn)[9]對該方法進(jìn)行了比較詳細(xì)的描述;在修正數(shù)據(jù)的過程中,需要先將數(shù)據(jù)變換成頻域,然后利用管傳遞函數(shù)反變換成時域;第2種方法則利用Matlab軟件中的ARMA模型辨識方法(自回歸滑動平均模型)對2組數(shù)據(jù)進(jìn)行識別,建立系統(tǒng)的動態(tài)模型,此分析方法為時域建模法,考慮了測量探針系統(tǒng)輸入和輸出信號的測量噪聲,其頻域特性曲線更光滑。本文采用第2種方法。2種方法的函數(shù)圖像對比如圖2所示。從圖中可見,在0~420 Hz內(nèi),采用2種方法得到的結(jié)果很接近,在更高的頻率區(qū)間內(nèi),由于存在測量噪聲、有限的激勵頻率成分以及有限的傳感器頻響等因素,2種方法產(chǎn)生了明顯差距。

    為了校驗得到的管傳遞函數(shù),在旋轉(zhuǎn)盤上打開了4個連續(xù)的孔,產(chǎn)生了4組頻率約為100 Hz的周期性脈沖。其中“標(biāo)準(zhǔn)信號”、探針原始信號、修正容腔效應(yīng)后信號三者的圖像對比如圖3所示。從圖中可見,在此激勵作用下,管傳遞函數(shù)模型的修正還是十分可靠的。

    圖3 管傳遞函數(shù)驗證

    圖4 單轉(zhuǎn)子壓氣機試驗臺

    1.2試驗臺及測試方案試驗在航空發(fā)動機氣動熱力國家級重點實驗室的大尺寸壓氣機試驗臺上進(jìn)行。研究對象如圖4所示,為1個單轉(zhuǎn)子,含17個葉片,葉片造型為C4葉型,葉頂間隙設(shè)計為3 mm。壓氣機外徑為1 m,輪轂比為0.6,葉片彎角為26.5°,安裝角為33.4°。壓氣機運轉(zhuǎn)的折合轉(zhuǎn)速為1100 r/min,轉(zhuǎn)速的波動約為±1 r/min?;谶M(jìn)口流速和葉片弦長的雷諾數(shù)約為7.5E5,試驗選取了2個典型的工作點:流量系數(shù)Cva為0.55的近失速點和0.65的設(shè)計點。

    測試方案的布置情況如圖5所示。其中,T3表示位移機構(gòu)安裝座,位于轉(zhuǎn)子尾緣下游約10%軸向弦長處,用于安裝位移機構(gòu)和探針,徑向測量站數(shù)為45站。為保證在每個徑向站探針頭部都正對來流方向,探針的安裝角度由5孔針測得的氣流偏轉(zhuǎn)角給定。

    圖5 測試方案布置

    轉(zhuǎn)子試驗件的幾何偏差如圖6所示。其中,葉型最大厚度與平均值的差別如圖6(a)所示,轉(zhuǎn)子進(jìn)、出口柵距的周向不均勻分布如圖6(b)所示。周向的非均勻特征為校驗非定常壓力場測量方法提供了1個良好的對象。

    在試驗中,13#葉片葉頂貼有錫箔紙,用于觸發(fā)光學(xué)探頭。在葉輪機測試過程中,硬件鎖相采樣平均技術(shù)是1種常見方法,然而在持續(xù)多周的測試中,葉輪機的旋轉(zhuǎn)角速度并不是絕對均勻的,既定的單周采樣點數(shù)可能導(dǎo)致轉(zhuǎn)子在單周測量中轉(zhuǎn)過的角位移有所不同;如果只對單次測量或系綜平均結(jié)果進(jìn)行頻域分析,其采樣點數(shù)少,頻率分辨率很低,而且間斷的單周測量系綜平均的數(shù)據(jù)會產(chǎn)生相位誤差。為保證數(shù)據(jù)的連續(xù)性,動態(tài)測量技術(shù)均使用軟件鎖相采樣平均技術(shù),即在試驗過程中,同時記錄連續(xù)的流場信號和鎖相觸發(fā)信號,在后處理中,在觸發(fā)信號的上升沿分離出每個單周的流場信號,進(jìn)行插值和系綜平均。

    2 試驗結(jié)果分析

    2.1總壓測量及修正結(jié)果

    在不同狀態(tài)、葉頂間隙下轉(zhuǎn)子出口的總壓系數(shù)分布情況如圖7所示??倝合禂?shù)定義為當(dāng)?shù)乜倝号c用葉尖切線速度計算的動壓頭之比。對比圖7(a)、(b)可知,在設(shè)計工況下,修正后的尾跡的周向分布變化較大,這主要是由于引壓管引起的遲滯作用造成的。主流區(qū)域的總壓分布的偏差較大,這主要是由于從尾跡到主流的周向突躍的頻率較高,傳感器并不能準(zhǔn)確識別;但在尖區(qū),由于泄漏流動的影響,總壓的周向變化呈近似正弦波動的特征,近似于轉(zhuǎn)子通過頻率,反而被捕捉的較為準(zhǔn)確,泄漏的形態(tài)和強度都比較一致。

    在設(shè)計工況下,轉(zhuǎn)子葉尖間隙分別為3 mm和1 mm時的總壓系數(shù)分布情況如圖7(b)、(c)所示。從圖中可見,在葉尖間隙為1 mm的情況下,泄漏流動的尺度和總壓虧損都有所降低,而在主流區(qū)域,流動的重復(fù)性較好,總壓分布較為一致。即總壓探針能夠分辨出總壓分布的不同,且能夠保持一直重復(fù)下去。

    圖7 轉(zhuǎn)子出口總壓系數(shù)分布

    在近失速工況(NS)下的總壓系數(shù)分布情況如圖7(d)、(e)、(f)所示。從圖中可見,修正后的主流區(qū)域流動更為均勻,這符合使用標(biāo)準(zhǔn)動態(tài)傳感器測得的轉(zhuǎn)子出口總壓分布[14-15]。類似地,泄漏流動的尺度和強度差異也能夠清晰地識別。

    2.2與參考動態(tài)流場的對比分析

    在不同狀態(tài)、不同葉頂間隙下單斜絲熱線測得的轉(zhuǎn)子出口的軸向速度系數(shù)分布如圖8所示。軸向速度系數(shù)定義為軸向速度與葉尖切向速度之比。對比圖7、8可見,由于修正后的總壓分布減弱了引壓管的遲滯效應(yīng),尾跡的周向位置更符合熱線測到的結(jié)果。同時,在設(shè)計點和近失速點,泄漏流的所在位置的對應(yīng)關(guān)系較好;在近失速點,泄漏流與葉尖吸力面尾跡有明顯的摻混現(xiàn)象。在葉片根部,轉(zhuǎn)子和輪轂之間的縫隙所產(chǎn)生的泄漏流動也可以被捕捉到。

    圖8 轉(zhuǎn)子出口軸向速度系數(shù)分布

    由于主流部分的測量結(jié)果受限于探針的頻響特性,不能很好地表現(xiàn)流動結(jié)構(gòu)。而在近失速工況下,中間葉高處各葉片尾跡寬度卻有所增加,表明測量站得到的動態(tài)數(shù)據(jù)所包含的頻率成分在向低頻成分移動。因此,經(jīng)過修正后,主流區(qū)域的流動更接近真實值。對比圖7(b)、8(a),以及圖7(e)、8(c),可以證明以上推論。

    2.3周向非均勻流場的識別

    在設(shè)計點、葉尖間隙為3 mm時轉(zhuǎn)子出口整周的總壓系數(shù)和軸向速度系數(shù)如圖9所示,其中總壓系數(shù)為修正后結(jié)果。對比(a)、(b)2圖可知,總壓探針對尾跡位置和泄漏流動的捕捉在整周上都與熱線測得的結(jié)果比較接近。由轉(zhuǎn)子周向非均勻?qū)е碌牧鲃硬町愐材芎芎玫乇辉撎结槻蹲降?,且二者和圖3所給出的幾何參數(shù)對應(yīng)較好,2#、6#、12#等葉片的尾跡相對相鄰葉片都較寬,故總壓虧損和軸向速度虧損都較大。

    從整周來看,欠頻響狀態(tài)下工作的動態(tài)壓力探針基本能夠反映出典型的泄漏流、尾跡等二次流動。但由于缺少高頻部分的流場信息,其周向的波動強度被弱化,這種規(guī)律與在真實壓氣機的試驗中也相同。

    圖9 3 mm葉尖間隙時轉(zhuǎn)子周向非對稱導(dǎo)致的周向流場不均勻分布

    圖9(a)測量結(jié)果FFT分析沿徑向分布的瀑布圖如圖10(a)所示。其中轉(zhuǎn)子葉片通過頻率為303 Hz。由于探針的固有頻率為110 Hz左右,總壓測量結(jié)果的頻譜在此頻率成分下被放大。由于受轉(zhuǎn)子葉片周向非均勻幾何分布的影響,在轉(zhuǎn)子頻率的整數(shù)倍成分上(1~16倍轉(zhuǎn)子頻率),均存在較清晰的峰值。在經(jīng)過修正后,110 Hz左右的頻率成分均被抑制,如圖10 (b)所示,頻域分析結(jié)果與熱線所得結(jié)果(圖11)類似。這種頻域中的偏差在高速壓氣機試驗中也存在,一方面是容腔效應(yīng)導(dǎo)致的,另一方面取決于傳感器本身的幅頻特性。與前文的修正方法類似,二者均可以通過動態(tài)補償?shù)玫揭欢ǖ男拚?,但修正后的結(jié)果仍將受限于探針的頻率特征,高頻成分的信息會有不同程度的偏差。

    圖10 動態(tài)探針測量結(jié)果的FFT分析

    2.4對旋轉(zhuǎn)探針法適用性的探討

    在葉輪機械的級間動態(tài)測量中,多孔動態(tài)探針由于尺寸原因,并不一定適用。研究者更多地使用單孔探針或雙孔探針配合旋轉(zhuǎn)探針法來獲得級間動態(tài)流場。旋轉(zhuǎn)探針法的應(yīng)用方法如下:在某一徑向高度,利用多次(3~11次)旋轉(zhuǎn)探針獲得的擬合曲線,與穩(wěn)態(tài)的氣動標(biāo)定數(shù)據(jù)進(jìn)行對比,從而得到氣流的偏轉(zhuǎn)角、俯仰角、總壓、靜壓、速度等參數(shù)。在這個過程中,以觸發(fā)信號為基準(zhǔn),提取每一時刻多次旋轉(zhuǎn)測得地原始信號作為當(dāng)前時刻的擬合曲線。前提是需在每個旋轉(zhuǎn)位置下進(jìn)行提取,測孔對氣流中的所有頻率成分的響應(yīng)的幅值和相位均一致,但從前文的分析來看,這顯然是不現(xiàn)實的。

    圖11 熱線探針測得軸向速度系數(shù)的FFT分析

    目前研究中探針的管傳遞函數(shù)的標(biāo)定主要集中在當(dāng)測壓孔正對來流時的頻率響應(yīng)。但當(dāng)測壓孔側(cè)對來流時頻率響應(yīng)的研究并不常見。利用總壓探針和旋轉(zhuǎn)探針法測得的偏轉(zhuǎn)角分布情況如圖12(a)所示,對三個旋轉(zhuǎn)位置(-30°,0°,30°)測得的壓力信號進(jìn)行修正,并計算所得的偏轉(zhuǎn)角分布情況如圖12(b)所示,熱線得到的偏轉(zhuǎn)角分布情況如圖12(c)所示。

    圖12 轉(zhuǎn)子出口氣流角分布

    對比圖12(a)、(b)以及標(biāo)準(zhǔn)分布圖12(c),雖然主流區(qū)域的氣流偏轉(zhuǎn)角略有改善,但尾跡區(qū)并沒有明顯改善,甚至產(chǎn)生了更大的誤差,在尖區(qū)和角區(qū)的泄漏流動只是在周向位置上有所改善,但區(qū)域和大小并沒有改善。因此,本文認(rèn)為,在使用旋轉(zhuǎn)探針法時,熱線和齊平安裝的動態(tài)壓力傳感器可以不考慮多個旋轉(zhuǎn)位置下的不同動態(tài)標(biāo)定問題,但有引壓管存在時,則必須謹(jǐn)慎處理,應(yīng)在動態(tài)標(biāo)定設(shè)備中對多個旋轉(zhuǎn)位置進(jìn)行動態(tài)標(biāo)定,并分別修正各個位置的數(shù)據(jù),然后再使用旋轉(zhuǎn)探針法處理氣流角數(shù)據(jù)。

    3 結(jié)論

    (1)通過動態(tài)修正,減小了引壓管的容腔效應(yīng)的影響。修正后,流場的周向位置與參考流場基本一致。在尾跡區(qū)域,由于尾跡到主流的周向突躍變化較快,測得結(jié)果偏差較大;相對的,由于受泄漏流動的影響,尖區(qū)總壓的周向變化呈現(xiàn)出近似正弦波動的特征,頻率成分集中于轉(zhuǎn)子通過頻率,因此被捕捉得較為準(zhǔn)確,泄漏的形態(tài)和強度都和標(biāo)準(zhǔn)場比較一致。

    (2)使用旋轉(zhuǎn)探針法測量葉輪機械中的動態(tài)流場時,如果傳感器前有引壓管,應(yīng)盡可能考慮對每個旋轉(zhuǎn)位置下都進(jìn)行動態(tài)標(biāo)定并分別利用標(biāo)定數(shù)據(jù)對原始數(shù)據(jù)進(jìn)行修正,以還原真實的流場信息。

    (3)本文的研究對象為低速壓氣機,但在高速設(shè)備中,以上結(jié)論也適用。衡量二次流損失的尾跡寬度、泄漏尺度等參數(shù)的選擇均應(yīng)在動態(tài)修正的基礎(chǔ)上審慎考慮。

    [1]KupferschmiedP,KoppelP,GizziW,etal.Time-resolvedflow measurements with Fast-Response Aerodynamic Probes in turbo machines[J].MeasurementScienceTechnology,2000,11(7):1036-1054.

    [2]Ainsworth R W,Miller R J,Moss R W,et al.Unsteady pressure measurement[J].Measurement Science Technology,2000,11(7):1055-1076.

    [3]Sieverding C,Arts T,Denos R,et al.Measurement techniques for unsteady flows in turbo machines[J].Experiments in Fluids,2000,28 (4):285-321.

    [4]Persico G,Gaetani P,Guardone A.Design and analysis of new concept Fast-ResponsePressureProbes[J].MeasurementScienceand Technology,2005,16(9):1741-1750.

    [5]Georgakis C,Bennett I,Ivey P C.Fast response probes measuring unsteady flows in high-speed research compressors[R].ASME 2003-GT-38160.

    [6]Brouckaert J F.Fast Response Aerodynamic Probes for measurements in turbo machines[J].Journal of Power and Energy,2007,221(6):801-813.

    [7]蔣浩康,李雨春,張洪,等.研究轉(zhuǎn)子內(nèi)流動的大尺寸軸流壓氣機試驗裝置和動態(tài)測試技術(shù)[J].航空動力學(xué)報,1992,7(1):1-8.

    JIANG Haokang,LI Yuchun,ZHANG Hong,et al.A large-scale axial flow compressor facility and dynamic measurement techniques for rotor flow study[J].Journal of Aerospace Power,1992,7(1):1-8.(in Chinese)

    [8]李井洋,馬宏偉,賀象.楔頂雙孔探針測量跨聲多級壓氣機轉(zhuǎn)子出口三維動態(tài)流場的方法 [J].航空動力學(xué)報,2012,27(10):2262-2268.

    LI Jingyang,MA Hongwei,HE Xiang.Method of measuring 3-D unsteady flow at exits of transonic compressor rotor passages using a two-hole tip-wedge pressure probe[J].Journal of Aerospace Power,2012,27(10):2262-2268.(in Chinese)

    [9]楊慧,鄭赟.測量葉輪機振動葉片表面非定常氣動響應(yīng)的實驗技術(shù)-影響系數(shù)法[J].實驗流體力學(xué),2011,25(3):71-77.

    YANG Hui,ZHENG Yun.Experimental technique for measuring unsteady aerodynamics of vibrating blade in turbo machinery influence coefficient method[J].Journal of Experiments in Fluid Mechanics,2011,25(3):71-77.(in Chinese)

    [10]李繼超,王偲臣,林峰,等.一種容腔效應(yīng)標(biāo)定技術(shù)及其在高頻響動態(tài)探針中的應(yīng)用[J].航空動力學(xué)報,2011,26(12):2749-2756.

    LI Jichao,WANG Sichen,LIN Feng,et al.A technique to calibrate cavity effect in unsteady pressure probes with high frequency response [J].Journal of Aerospace Power,2011,26(12):2749-2756.(in Chinese)

    [11]楊林,王偲臣,林峰,等.基于半無限長管效應(yīng)的動態(tài)壓力測量方法[J].航空動力學(xué)報,2014,29(10):2457-2463.

    YANG Lin,WANG Sichen,LIN Feng,et al.Method ofdynamic pressure measurement based on semi-infinite pressure tube effect[J]. Journal of Aerospace Power,2014,9(10):2457-2463.(in Chinese)

    [12]郭勇,王效葵,王洪偉.高頻壓力探針動態(tài)性能的補償方法[J].航空動力學(xué)報,2005,20(2):273-277.

    GUOYong,WANGXiaokui,WANGHongwei.Improvementon dynamic characteristics of high response pressure probe[J].Journal of Aerospace Power,2005,20(2):273-277.(in Chinese)

    [13]勞力云,陳珙,張宏建,等.引壓管路對動態(tài)差壓檢測影響的分析[J].計量學(xué)報,1999,20(4):286-291.

    LAO Liyun,CHEN Gong,ZHANG Hongjian,et al.The effect of pressure ducts on dynamic characteristics of pressure drop measuring system[J].Acta Metrologica Sinica,1999,20(4):286-291.(in Chinese)

    [14]馬宏偉,蔣浩康.壓氣機轉(zhuǎn)子出口流場的發(fā)展及三維紊流特性[J].航空動力學(xué)報,2000,15(2):113-118.

    MA Hongwei,JIANG Haokang.Development of three-dimensional turbulent flow field at exit of an axial compressor rotor passage[J]. Journal of Aerospace Power,2000,15(2):113-118.(in Chinese)

    [15]馬宏偉,蔣浩康.壓氣機不同狀態(tài)下轉(zhuǎn)子出口三維紊流流場[J].航空動力學(xué)報,1997,12(3):268-272.

    MA Hongwei,JIANG Haokang.3D turbulent flow at compressor rotor exit under different flow conditions[J].Journal of Aerospace Power,1997,12(3):268-272.(in Chinese)

    (編輯:栗樞)

    Analysis of Measured Unsteady Flow Field Using a Quasi Fast Response Pressure Probe

    MA Hong-wei1,2,WEI Wei1,2,ZHANG Liang1,2,GUAN Hui1,2
    (1.School of Energy and Power Engineering,Beihang University;2.Collaborative Innovation Center of Advanced Aero-Engine,Beijing 100191,China)

    It is well known that the flow filed inside a compressor has strong dynamic characteristic,especially in the small-scale ones whose blade passing frequency could reach 30~80 kHz.In most circumstances,because the Fast Response Aerodynamic Probe(FRAP)can only meet the blade passing frequency,the dynamic characteristic of the probe can not satisfy the requirement of reproducing a detailed high-frequency flow structure.In order to investigate the difference between measured flow field and real flow field using this quasi FRAP,a dynamic pressure probe with a response frequency of 420 Hz and a slanted single hot wire probe with a response frequency more than 33 kHz was utilized to measure the unsteady flow filed at the outlet of an isolated compressor rotor which had a blade passing frequency of 303 Hz.The frequency response of the probe was measured to correct the influence of square-cavity effect on the measured flow field.The corresponding corrected ones were compared with the reference flow field measured with hot wire probe as the response frequency of measurement system was low.The comparison shows that some frequency components of the original signals were amplified.After correction,the circumferential position of the unsteady field,the blockage and pattern of the tip leakage area are obviously more close to the reference flow field.Nevertheless,due to the high frequency component in the wake,the correction still cannot get a good result.

    quasi FRAP;dynamic measurement;tube transfer function;dynamic correction;compressor;areoengine

    V 231.3

    A

    10.13477/j.cnki.aeroengine.2016.02.013

    2015-06-12基金項目:國家自然科學(xué)基金(51136003)資助

    馬宏偉(1968),男,博士,博士生導(dǎo)師,研究方向為葉輪機復(fù)雜流動,發(fā)動機氣動測試技術(shù);E-mail:mahw@buaa.edu.cn。

    引用格式:馬宏偉,魏巍,張良,等.欠頻響壓力探針測量壓氣機動態(tài)流場的結(jié)果分析[J].航空發(fā)動機,2016,42(2):67-72.MAHongwei,WEIWei,ZHANG Liang,etal.AnalysisofmeasuredunsteadyflowfieldusingaquasiFastResponsePressureProbe[J].Aeroengine,2016,42(2):67-72.

    猜你喜歡
    頻響總壓壓氣機
    總壓探針性能結(jié)構(gòu)敏感性分析
    軸流壓氣機效率評定方法
    重型燃?xì)廨啓C壓氣機第一級轉(zhuǎn)子葉片斷裂分析
    基于分塊化頻響函數(shù)曲率比的砌體房屋模型損傷識別研究
    可調(diào)式總壓耙設(shè)計及應(yīng)用
    壓氣機緊湊S形過渡段內(nèi)周向彎靜子性能數(shù)值計算
    亞聲速條件下總壓探針臨壁效應(yīng)的數(shù)值研究
    2 m超聲速風(fēng)洞流場變速壓控制方法研究
    美團外賣哥
    頻響函數(shù)殘差法在有限元模型修正中的應(yīng)用
    亚洲熟妇熟女久久| 亚洲欧美激情综合另类| 婷婷丁香在线五月| 亚洲精品久久午夜乱码| 久久久精品区二区三区| 他把我摸到了高潮在线观看| 少妇裸体淫交视频免费看高清 | 18在线观看网站| 国产精品亚洲av一区麻豆| 日日夜夜操网爽| 一二三四在线观看免费中文在| 飞空精品影院首页| 老汉色av国产亚洲站长工具| 国产精华一区二区三区| 欧美精品av麻豆av| 国产乱人伦免费视频| 国产亚洲精品第一综合不卡| 老司机影院毛片| 欧美日韩一级在线毛片| 一区二区三区国产精品乱码| 欧美+亚洲+日韩+国产| 黄色女人牲交| 国产精品 国内视频| 最新在线观看一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 在线播放国产精品三级| 欧美日韩精品网址| 午夜免费观看网址| 亚洲五月婷婷丁香| 免费观看精品视频网站| 中文字幕制服av| 国产精品二区激情视频| 免费看十八禁软件| 国产精品免费大片| 久久国产精品影院| 欧美老熟妇乱子伦牲交| 久久精品91无色码中文字幕| 亚洲欧洲精品一区二区精品久久久| www日本在线高清视频| 国产亚洲精品第一综合不卡| 69av精品久久久久久| 中文字幕色久视频| 下体分泌物呈黄色| 国产不卡av网站在线观看| 电影成人av| 欧美另类亚洲清纯唯美| 日韩中文字幕欧美一区二区| 99热只有精品国产| 别揉我奶头~嗯~啊~动态视频| 亚洲在线自拍视频| 老汉色∧v一级毛片| 黄片小视频在线播放| 又大又爽又粗| 精品少妇久久久久久888优播| 久久久久久免费高清国产稀缺| 一本一本久久a久久精品综合妖精| 亚洲五月婷婷丁香| 真人做人爱边吃奶动态| 老熟女久久久| 亚洲七黄色美女视频| 久久久久视频综合| 又黄又爽又免费观看的视频| 国产欧美日韩一区二区三| 国产精品久久视频播放| 国产在线精品亚洲第一网站| 国产极品粉嫩免费观看在线| 欧美日韩一级在线毛片| 十八禁人妻一区二区| 欧美精品人与动牲交sv欧美| 亚洲国产欧美一区二区综合| 久9热在线精品视频| 熟女少妇亚洲综合色aaa.| 国产aⅴ精品一区二区三区波| 日日夜夜操网爽| 大片电影免费在线观看免费| 无限看片的www在线观看| 在线视频色国产色| 性色av乱码一区二区三区2| 少妇的丰满在线观看| 欧美人与性动交α欧美软件| 午夜福利一区二区在线看| 久久亚洲真实| 亚洲精品国产一区二区精华液| 亚洲一区中文字幕在线| 国产精品免费一区二区三区在线 | 国产激情欧美一区二区| 亚洲成人免费av在线播放| 国内毛片毛片毛片毛片毛片| av一本久久久久| 国产视频一区二区在线看| 欧美老熟妇乱子伦牲交| 美女 人体艺术 gogo| 亚洲av片天天在线观看| 成年动漫av网址| av超薄肉色丝袜交足视频| 国产精品永久免费网站| 国产一区二区三区视频了| 亚洲精品粉嫩美女一区| 国产精品秋霞免费鲁丝片| 99久久99久久久精品蜜桃| 村上凉子中文字幕在线| 久久天躁狠狠躁夜夜2o2o| 成人18禁在线播放| 超碰97精品在线观看| 99在线人妻在线中文字幕 | 精品一区二区三卡| 中国美女看黄片| 国产aⅴ精品一区二区三区波| 欧美久久黑人一区二区| 国产精品久久久av美女十八| 韩国精品一区二区三区| 成年动漫av网址| 午夜福利免费观看在线| 久久亚洲真实| 精品一区二区三区四区五区乱码| 老司机午夜福利在线观看视频| 人人澡人人妻人| 两性午夜刺激爽爽歪歪视频在线观看 | 19禁男女啪啪无遮挡网站| 久久久久久久久久久久大奶| 亚洲人成电影免费在线| 婷婷丁香在线五月| 成人精品一区二区免费| 一边摸一边做爽爽视频免费| 日韩制服丝袜自拍偷拍| 久久精品国产清高在天天线| 国产亚洲欧美在线一区二区| 男人的好看免费观看在线视频 | 成年人免费黄色播放视频| 女人久久www免费人成看片| а√天堂www在线а√下载 | av有码第一页| aaaaa片日本免费| 99国产精品一区二区三区| 久久精品国产a三级三级三级| 天天躁夜夜躁狠狠躁躁| 两人在一起打扑克的视频| 性少妇av在线| 久久人妻av系列| 精品高清国产在线一区| 久久精品亚洲熟妇少妇任你| 亚洲色图综合在线观看| 国产精品免费大片| 伦理电影免费视频| 欧美不卡视频在线免费观看 | 国产午夜精品久久久久久| 69精品国产乱码久久久| 色婷婷久久久亚洲欧美| 亚洲av成人av| 亚洲avbb在线观看| 国产深夜福利视频在线观看| 最近最新免费中文字幕在线| av国产精品久久久久影院| 亚洲美女黄片视频| 亚洲一区二区三区欧美精品| 精品卡一卡二卡四卡免费| av欧美777| 国产精品偷伦视频观看了| 日日爽夜夜爽网站| 亚洲欧美色中文字幕在线| 中亚洲国语对白在线视频| av国产精品久久久久影院| tube8黄色片| 色精品久久人妻99蜜桃| 国产成人免费无遮挡视频| 午夜影院日韩av| 中国美女看黄片| 真人做人爱边吃奶动态| 91麻豆精品激情在线观看国产 | 成人精品一区二区免费| 国产欧美日韩一区二区三| 国产精品美女特级片免费视频播放器 | 最新在线观看一区二区三区| 9热在线视频观看99| 高清视频免费观看一区二区| 可以免费在线观看a视频的电影网站| 麻豆乱淫一区二区| 亚洲精品中文字幕一二三四区| 久久久国产成人免费| 一级毛片精品| 国产真人三级小视频在线观看| 女人久久www免费人成看片| 欧美日韩视频精品一区| 欧美日韩国产mv在线观看视频| 日韩 欧美 亚洲 中文字幕| 黄片大片在线免费观看| 中文字幕高清在线视频| 他把我摸到了高潮在线观看| 无人区码免费观看不卡| 欧美成人午夜精品| 高清黄色对白视频在线免费看| xxx96com| 精品亚洲成国产av| 中国美女看黄片| 丝袜美足系列| 亚洲视频免费观看视频| 在线av久久热| 成人av一区二区三区在线看| 这个男人来自地球电影免费观看| 亚洲专区字幕在线| 国产午夜精品久久久久久| 人妻久久中文字幕网| 妹子高潮喷水视频| 久久天躁狠狠躁夜夜2o2o| 看黄色毛片网站| 亚洲精品久久成人aⅴ小说| 在线观看日韩欧美| 精品福利观看| 大香蕉久久网| 色尼玛亚洲综合影院| 亚洲专区中文字幕在线| 国产精品成人在线| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲 | 亚洲av美国av| 18禁裸乳无遮挡动漫免费视频| 精品国产超薄肉色丝袜足j| 午夜亚洲福利在线播放| 国产99久久九九免费精品| 一级a爱片免费观看的视频| 成人手机av| 中文字幕精品免费在线观看视频| 中文字幕人妻丝袜一区二区| 成人三级做爰电影| 纯流量卡能插随身wifi吗| 国产精品久久久久久人妻精品电影| 免费在线观看日本一区| 他把我摸到了高潮在线观看| 国产亚洲欧美98| 亚洲成人免费电影在线观看| 久久中文字幕人妻熟女| 一级毛片高清免费大全| 亚洲男人天堂网一区| 女人被狂操c到高潮| 在线观看www视频免费| 一夜夜www| 丁香六月欧美| 久久精品国产亚洲av香蕉五月 | 精品人妻在线不人妻| 黄色片一级片一级黄色片| 黄色成人免费大全| 国产av精品麻豆| 精品第一国产精品| 日韩欧美一区视频在线观看| 色婷婷av一区二区三区视频| 欧美日韩亚洲高清精品| 嫁个100分男人电影在线观看| av电影中文网址| 天堂√8在线中文| 夫妻午夜视频| 搡老岳熟女国产| 69精品国产乱码久久久| 在线观看www视频免费| 亚洲av日韩在线播放| av国产精品久久久久影院| 国产蜜桃级精品一区二区三区 | 国产av一区二区精品久久| 久久这里只有精品19| 国产精品亚洲一级av第二区| 大型av网站在线播放| 91精品国产国语对白视频| 国产免费av片在线观看野外av| 十分钟在线观看高清视频www| 身体一侧抽搐| 91国产中文字幕| 下体分泌物呈黄色| 亚洲av熟女| 99国产精品99久久久久| 人人妻人人爽人人添夜夜欢视频| 成人影院久久| 国产又爽黄色视频| 国产成人欧美在线观看 | 久久久久国产一级毛片高清牌| 18禁观看日本| 一级作爱视频免费观看| 久久久精品区二区三区| 国产精品久久视频播放| 水蜜桃什么品种好| 免费人成视频x8x8入口观看| 超碰成人久久| 免费女性裸体啪啪无遮挡网站| 久久青草综合色| 国产不卡av网站在线观看| 精品午夜福利视频在线观看一区| 国产男靠女视频免费网站| 一级作爱视频免费观看| 亚洲国产欧美日韩在线播放| 久久性视频一级片| 丰满饥渴人妻一区二区三| 大型av网站在线播放| 看免费av毛片| www.自偷自拍.com| av视频免费观看在线观看| 极品教师在线免费播放| tocl精华| 又黄又粗又硬又大视频| 国产精品久久久人人做人人爽| 99香蕉大伊视频| 色婷婷久久久亚洲欧美| 18禁国产床啪视频网站| 老鸭窝网址在线观看| 久久久国产一区二区| 午夜视频精品福利| 午夜精品久久久久久毛片777| 麻豆成人av在线观看| 男女下面插进去视频免费观看| 黑人操中国人逼视频| 岛国毛片在线播放| 脱女人内裤的视频| 日韩三级视频一区二区三区| 国产成人av激情在线播放| 他把我摸到了高潮在线观看| 丰满人妻熟妇乱又伦精品不卡| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成国产人片在线观看| 国产精品影院久久| 美女 人体艺术 gogo| 飞空精品影院首页| 亚洲人成电影免费在线| 999久久久国产精品视频| 成年人黄色毛片网站| 在线观看66精品国产| 亚洲av熟女| 国产精品久久久久成人av| 久久久国产精品麻豆| 夜夜躁狠狠躁天天躁| 日韩欧美一区二区三区在线观看 | 亚洲专区中文字幕在线| 三级毛片av免费| 丝袜在线中文字幕| 一区在线观看完整版| 麻豆乱淫一区二区| 成人手机av| 亚洲三区欧美一区| 成年人黄色毛片网站| 一区二区三区激情视频| 久久九九热精品免费| 侵犯人妻中文字幕一二三四区| 最近最新中文字幕大全免费视频| 国产精品乱码一区二三区的特点 | 国产精品九九99| 成年动漫av网址| 欧美成狂野欧美在线观看| 少妇 在线观看| 天天操日日干夜夜撸| 在线十欧美十亚洲十日本专区| 国产亚洲精品久久久久5区| 亚洲av片天天在线观看| 久久中文字幕一级| 无限看片的www在线观看| 久久中文字幕一级| 免费日韩欧美在线观看| 亚洲黑人精品在线| 亚洲五月色婷婷综合| 夜夜爽天天搞| 欧美成人午夜精品| 如日韩欧美国产精品一区二区三区| 亚洲三区欧美一区| 亚洲熟女精品中文字幕| 悠悠久久av| 午夜精品久久久久久毛片777| 婷婷成人精品国产| 国产精品美女特级片免费视频播放器 | 精品少妇久久久久久888优播| 男女床上黄色一级片免费看| 亚洲情色 制服丝袜| 1024香蕉在线观看| 精品国产乱子伦一区二区三区| 一二三四社区在线视频社区8| 高潮久久久久久久久久久不卡| 日韩欧美三级三区| 午夜免费鲁丝| 欧美日韩中文字幕国产精品一区二区三区 | 黄频高清免费视频| 曰老女人黄片| 亚洲精华国产精华精| 久久国产精品影院| 午夜精品在线福利| 一进一出抽搐gif免费好疼 | 亚洲午夜理论影院| 亚洲九九香蕉| 一边摸一边抽搐一进一出视频| 亚洲av第一区精品v没综合| 免费日韩欧美在线观看| 国产精品久久视频播放| 妹子高潮喷水视频| 国产欧美日韩综合在线一区二区| 极品少妇高潮喷水抽搐| 亚洲综合色网址| 精品一品国产午夜福利视频| 黄色丝袜av网址大全| a级片在线免费高清观看视频| 亚洲欧洲精品一区二区精品久久久| 久久精品国产清高在天天线| 免费观看精品视频网站| 一个人免费在线观看的高清视频| 亚洲自偷自拍图片 自拍| 最近最新中文字幕大全免费视频| 手机成人av网站| 久久天躁狠狠躁夜夜2o2o| 一边摸一边做爽爽视频免费| 国产色视频综合| 欧美日韩亚洲高清精品| 亚洲成a人片在线一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 午夜成年电影在线免费观看| 人妻一区二区av| 亚洲色图综合在线观看| 叶爱在线成人免费视频播放| 下体分泌物呈黄色| 国产精品成人在线| 嫩草影视91久久| x7x7x7水蜜桃| 国产在线观看jvid| √禁漫天堂资源中文www| 1024香蕉在线观看| 亚洲精品av麻豆狂野| 啦啦啦 在线观看视频| 久久国产亚洲av麻豆专区| 成人亚洲精品一区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲黑人精品在线| 久久久久久免费高清国产稀缺| 别揉我奶头~嗯~啊~动态视频| 老熟女久久久| 欧美成人免费av一区二区三区 | 18禁裸乳无遮挡动漫免费视频| 国产一区在线观看成人免费| 久久午夜综合久久蜜桃| 国产三级黄色录像| 欧美 亚洲 国产 日韩一| 国产成人精品无人区| 国产精品98久久久久久宅男小说| 久久性视频一级片| 日韩欧美国产一区二区入口| 久久精品国产99精品国产亚洲性色 | 久久热在线av| 亚洲精品美女久久av网站| 亚洲国产欧美网| 51午夜福利影视在线观看| 下体分泌物呈黄色| 亚洲色图 男人天堂 中文字幕| 欧美日韩亚洲高清精品| 亚洲成国产人片在线观看| 久久亚洲精品不卡| 国产在线精品亚洲第一网站| 午夜福利在线免费观看网站| 老司机亚洲免费影院| 亚洲国产欧美一区二区综合| 黄色 视频免费看| 人人妻人人澡人人看| 亚洲国产毛片av蜜桃av| 国产精品一区二区在线不卡| 熟女少妇亚洲综合色aaa.| 搡老岳熟女国产| 久久久久精品国产欧美久久久| 国产精品免费一区二区三区在线 | 久久中文看片网| 啦啦啦视频在线资源免费观看| 91av网站免费观看| 欧美成人午夜精品| 三级毛片av免费| 国产又爽黄色视频| 国产在线精品亚洲第一网站| 国产精品自产拍在线观看55亚洲 | 久9热在线精品视频| 日本五十路高清| 成熟少妇高潮喷水视频| 欧美老熟妇乱子伦牲交| 在线国产一区二区在线| 夜夜爽天天搞| 午夜视频精品福利| 99久久综合精品五月天人人| 亚洲男人天堂网一区| 午夜日韩欧美国产| 电影成人av| 老司机靠b影院| 老司机午夜十八禁免费视频| 久久久精品免费免费高清| 久久久久久久久免费视频了| videosex国产| 真人做人爱边吃奶动态| 午夜福利免费观看在线| 久久国产亚洲av麻豆专区| 久久中文看片网| av网站在线播放免费| 亚洲欧美色中文字幕在线| 99热国产这里只有精品6| 变态另类成人亚洲欧美熟女 | 午夜精品国产一区二区电影| 亚洲色图 男人天堂 中文字幕| 天堂俺去俺来也www色官网| 最近最新免费中文字幕在线| 欧美不卡视频在线免费观看 | 国产精品免费一区二区三区在线 | 女人被躁到高潮嗷嗷叫费观| 国产精品偷伦视频观看了| 亚洲欧美激情在线| 中国美女看黄片| 又紧又爽又黄一区二区| 最新的欧美精品一区二区| 一二三四在线观看免费中文在| 日韩人妻精品一区2区三区| 国产精品99久久99久久久不卡| 夜夜爽天天搞| 精品国产超薄肉色丝袜足j| 操美女的视频在线观看| 国产精品久久久人人做人人爽| 国产成人免费无遮挡视频| 国产成人精品无人区| 婷婷成人精品国产| 大码成人一级视频| 一本一本久久a久久精品综合妖精| 亚洲中文日韩欧美视频| 国产精品久久久久久人妻精品电影| 悠悠久久av| 欧美成人免费av一区二区三区 | 正在播放国产对白刺激| 久久久久久久午夜电影 | 一二三四在线观看免费中文在| 亚洲人成电影免费在线| 一区福利在线观看| av国产精品久久久久影院| 男女免费视频国产| 成人av一区二区三区在线看| 人妻久久中文字幕网| 日本黄色日本黄色录像| 两性午夜刺激爽爽歪歪视频在线观看 | 很黄的视频免费| 黄色女人牲交| 久久精品国产综合久久久| 久久久久国内视频| videos熟女内射| 久久人妻av系列| 91av网站免费观看| 身体一侧抽搐| 人妻久久中文字幕网| 99在线人妻在线中文字幕 | 亚洲熟女精品中文字幕| 99re在线观看精品视频| 啦啦啦视频在线资源免费观看| 夜夜爽天天搞| av线在线观看网站| 国产精品 欧美亚洲| 露出奶头的视频| 久久精品亚洲熟妇少妇任你| 久久久国产成人免费| 精品视频人人做人人爽| 一级毛片高清免费大全| 国产黄色免费在线视频| 婷婷精品国产亚洲av在线 | 日韩视频一区二区在线观看| 麻豆乱淫一区二区| 18禁裸乳无遮挡动漫免费视频| 男女床上黄色一级片免费看| 色婷婷久久久亚洲欧美| 免费在线观看黄色视频的| 精品视频人人做人人爽| 热re99久久国产66热| 又大又爽又粗| 12—13女人毛片做爰片一| 在线观看免费午夜福利视频| 国产在线精品亚洲第一网站| 日韩欧美国产一区二区入口| 午夜老司机福利片| 女人被躁到高潮嗷嗷叫费观| 波多野结衣一区麻豆| 精品电影一区二区在线| 三级毛片av免费| 国产野战对白在线观看| av网站在线播放免费| 午夜免费观看网址| xxxhd国产人妻xxx| 亚洲人成77777在线视频| 国产亚洲精品一区二区www | 国产aⅴ精品一区二区三区波| 啦啦啦 在线观看视频| 搡老熟女国产l中国老女人| 国产成人精品无人区| 麻豆乱淫一区二区| 亚洲人成电影免费在线| 久久国产乱子伦精品免费另类| 亚洲精品国产一区二区精华液| 亚洲第一青青草原| 真人做人爱边吃奶动态| 欧美日韩乱码在线| 曰老女人黄片| 制服诱惑二区| 国产日韩欧美亚洲二区| 日韩欧美在线二视频 | 久久久久久久久免费视频了| 久久久国产精品麻豆| 久久久久国产精品人妻aⅴ院 | 亚洲精品久久成人aⅴ小说| 王馨瑶露胸无遮挡在线观看| 欧美成狂野欧美在线观看| 亚洲精品久久成人aⅴ小说| 欧美乱码精品一区二区三区| 欧美色视频一区免费| 精品一品国产午夜福利视频| 咕卡用的链子| 午夜精品久久久久久毛片777| 日韩欧美一区二区三区在线观看 | 视频区欧美日本亚洲| 天堂俺去俺来也www色官网| 一级片'在线观看视频| 黑丝袜美女国产一区| 亚洲国产欧美日韩在线播放| av一本久久久久| а√天堂www在线а√下载 | 多毛熟女@视频| 中文字幕另类日韩欧美亚洲嫩草| 黄色 视频免费看| 久久ye,这里只有精品| 性色av乱码一区二区三区2|