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

    壓氣機氣動失穩(wěn)故障的聲學(xué)診斷技術(shù)研究進展

    2024-01-12 13:54:38陳禹西楊明綏楊國安
    測控技術(shù) 2023年12期
    關(guān)鍵詞:壓氣機聲學(xué)氣動

    陳禹西, 楊明綏, 楊國安

    (1.中國航發(fā)沈陽發(fā)動機研究所,遼寧 沈陽 110015; 2.北京化工大學(xué) 機電工程學(xué)院,北京 100029)

    作為航空發(fā)動機的重要部件,壓氣機處于不穩(wěn)定狀態(tài)時,不僅其性能會下降,還可能對自身結(jié)構(gòu)產(chǎn)生額外載荷,導(dǎo)致其壽命減少甚至結(jié)構(gòu)失效,對發(fā)動機正常工作甚至整機飛行安全帶來嚴重影響。圖1展示了氣動失穩(wěn)對壓氣機葉片造成的結(jié)構(gòu)破壞。為改善壓氣機的穩(wěn)定性,研究人員提出了各種故障檢測與診斷的方案,這些方案一般通過壓氣機上安裝的傳感器監(jiān)測壓氣機失穩(wěn)時的異常信號,然后向控制系統(tǒng)發(fā)出指令,調(diào)節(jié)壓氣機的穩(wěn)定裕度或使壓氣機的工作點遠離不穩(wěn)定邊界。這種失穩(wěn)控制手段已經(jīng)在整機上得到應(yīng)用[1-2]。但目前這種方法僅能在失穩(wěn)發(fā)生時發(fā)揮作用,并不能提前預(yù)防失穩(wěn)的發(fā)生。為了更加及時地捕捉到壓氣機的不穩(wěn)定工作狀態(tài),減少壓氣機不穩(wěn)定工作的時間,有必要探究基于各種傳感器方案的氣動失穩(wěn)監(jiān)測和預(yù)警技術(shù)。

    圖1 喘振導(dǎo)致的壓氣機葉片結(jié)構(gòu)破壞[3]

    已有經(jīng)驗表明,發(fā)動機整機及部件在運行時往往會出現(xiàn)一些反映其運行狀態(tài)和運行品質(zhì)的異常聲響或振動,而這些異常聲振現(xiàn)象發(fā)生時均在其聲音信號上有顯著體現(xiàn)[4],可以被試車員、飛行員等依據(jù)“響聲”“異響”主觀感知到,這為使用聲學(xué)原理和聲學(xué)信號研究壓氣機內(nèi)的流動和振動現(xiàn)象奠定了基礎(chǔ)。隨著研究的深入,人們逐漸認識到聲學(xué)信號不僅能夠有效表征葉片氣動故障與振動故障信息,而且聲信號獨有的波動傳播理論在葉片流固耦合機理研究方面具有天然優(yōu)勢,是分析葉片非同步共振、聲腔共振、氣動失速、喘振、顫振等故障機理的優(yōu)異分析手段。Sundstr?m等[5]提出,在低質(zhì)量流量下,壓氣機發(fā)出的噪聲主要來源于作用在葉片表面的非定常壓力載荷。Cumpsty[6]也提出,氣動失穩(wěn)時的壓氣機具有特定的聲學(xué)特征。聲學(xué)測量已經(jīng)在壓氣機故障研究中獲得了一些應(yīng)用[4,7],基于聲學(xué)信號的發(fā)動機葉片故障監(jiān)測與診斷技術(shù)具有非接觸、強便攜性、便于安裝、無需對壓氣機進行過多測試改裝等優(yōu)點[8],且在航空發(fā)動機氣動運行狀態(tài)、旋轉(zhuǎn)機械運行品質(zhì)檢測方面,聲學(xué)信號獨有的“可有效與人耳聽覺的主觀評價相關(guān)聯(lián),且可量化”的優(yōu)點,最具潛力發(fā)展形成“人機合一”綜合檢測與診斷能力的物理信號。通過對聲信號的深入研究分析,研究人員有望研制出更加準確、靈敏的壓氣機故障診斷、監(jiān)測和預(yù)警方法。

    因此,本文主要關(guān)注國內(nèi)外基于聲學(xué)手段的壓氣機氣動失穩(wěn)故障診斷研究的進展情況,對相關(guān)的研究工作進行梳理、提煉和闡述,理清氣動失穩(wěn)現(xiàn)象與聲信號表征的機理分析、試驗方案、測試布局與方法、數(shù)據(jù)處理手段等技術(shù)發(fā)展脈絡(luò)。在此基礎(chǔ)上,提煉分析并給出了壓氣機氣動失穩(wěn)故障的聲學(xué)監(jiān)測與預(yù)警方法的發(fā)展趨勢與關(guān)鍵技術(shù),為壓氣機失穩(wěn)檢測診斷手段的發(fā)展與完善提供一定的指導(dǎo)。

    1 氣動失穩(wěn)的發(fā)聲機理與信號特征

    旋轉(zhuǎn)失速與喘振是典型的壓氣機氣動失穩(wěn)現(xiàn)象,隨著壓氣機內(nèi)部流動條件惡化,一些初始擾動出現(xiàn),流道內(nèi)會出現(xiàn)流動分離,壓氣機發(fā)生旋轉(zhuǎn)失速,若流動狀況進一步惡化,壓氣機則可能發(fā)生喘振。當轉(zhuǎn)速足夠高時,流動條件惡化也可能直接造成喘振。不同轉(zhuǎn)速下的失穩(wěn)形式如圖2所示。旋轉(zhuǎn)失速發(fā)生時,某一個或幾個葉片流道首先發(fā)生氣流分離,形成分離團,如圖3所示,從整體上看,分離團在轉(zhuǎn)子上逆轉(zhuǎn)子旋轉(zhuǎn)方向轉(zhuǎn)動。在分離區(qū)域內(nèi),葉片的流動分離與渦系結(jié)構(gòu)十分復(fù)雜,具有強烈的非定常性,形成了獨特的聲源,圖4展示了失穩(wěn)過程中聲學(xué)信號的變化。一些情況下,在失穩(wěn)發(fā)生前,會首先出現(xiàn)模態(tài)波的特征頻率,并伴有“異?!甭暡ǖ漠a(chǎn)生;另一些情況下失穩(wěn)過程發(fā)生得更突然,未觀察到明顯的先兆現(xiàn)象即直接進入失速或喘振狀態(tài)。進入失穩(wěn)狀態(tài)后,聲信號會發(fā)生明顯變化,包括幅值的顯著波動和頻率成分的改變。

    圖2 不同轉(zhuǎn)速下的失穩(wěn)形式

    圖3 基元中的流動分離

    圖4 失穩(wěn)過程中聲學(xué)信號變化

    旋轉(zhuǎn)不穩(wěn)定性是另一種出現(xiàn)在壓氣機高負荷工況時的動失穩(wěn)現(xiàn)象,目前普遍認為其與葉尖泄露而形成的漩渦有關(guān),不同于失速或喘振,在信號特征上,旋轉(zhuǎn)不穩(wěn)定性可以視為繞轉(zhuǎn)子旋轉(zhuǎn)的脈動聲源。如圖5所示,脈動聲源的周向分布可以看成多個連續(xù)階數(shù)模態(tài)空間傅里葉分量的疊加,導(dǎo)致旋轉(zhuǎn)不穩(wěn)定性的頻率特性一般表現(xiàn)為離散頻率峰值及其諧波。如圖6所示,RI頻率通常伴隨著BPF+RI和BPF-RI的出現(xiàn),而且RI 頻帶通常表現(xiàn)為多重離散峰值的疊加,相鄰2個離散峰的間隔頻率幾乎相等。

    圖5 旋轉(zhuǎn)不穩(wěn)定聲源特征

    圖6 旋轉(zhuǎn)不穩(wěn)定的頻率特征

    2 氣動失穩(wěn)的聲學(xué)表征方法

    借助各種信號分析手段,獲取氣動失穩(wěn)的聲學(xué)表征是發(fā)展相應(yīng)診斷技術(shù)的必要途徑。為此,國內(nèi)外研究人員圍繞著聲信號的采集與分析處理開展了大量工作。隨著技術(shù)的發(fā)展,在采集方法方面,手段逐漸趨于復(fù)雜,通道數(shù)量、采樣率等參數(shù)不斷提升;在分析處理方面,從簡單時頻表征手段逐漸發(fā)展到高階參數(shù)表征、模態(tài)分量表征、空間特征表征。其中信號分析處理技術(shù)的發(fā)展起到了重要的推動作用。因此本節(jié)主要從信號處理技術(shù)的角度,對氣動失穩(wěn)的聲學(xué)表征方法進行梳理。

    2.1 基于傳統(tǒng)時頻理論的表征方法

    早期的失穩(wěn)聲學(xué)測量研究使用的傳聲器數(shù)目較少,信號處理分析的手段也僅集中于基本的時頻分析,盡管如此,研究人員已經(jīng)能夠通過比較正常與失穩(wěn)狀態(tài)下聲信號的差異,從一些時域統(tǒng)計量和頻率成分變化的角度對失穩(wěn)信號的特征進行表征。

    國內(nèi)可以檢索到的最早使用聲學(xué)手段研究壓氣機失穩(wěn)的工作為20世紀末胡宗安等[9-10]在單、雙級壓氣機試驗臺上進行的旋轉(zhuǎn)失速與顫振現(xiàn)象研究,他們提出了根據(jù)聲信號中滿足倍頻程關(guān)系的頻率確定失速團傳播速度的方法。同一時期,Lawless等[11-12]也對某離心壓氣機失穩(wěn)現(xiàn)象的聲信號特性進行了研究,描述了隨著流量減少,聲信號中失速特征頻率和喘振特征頻率的變化,還開創(chuàng)性地使用了環(huán)形傳聲器陣列獲得了壓氣機的聲模態(tài)信息,發(fā)現(xiàn)部分模態(tài)在時域信號的幅值出現(xiàn)異常前就已經(jīng)發(fā)生了變化。Aretakis等[13-14]使用如圖7所示的測試布局研究了某離心壓氣機在幾種結(jié)構(gòu)故障下,振動、脈動壓力與聲信號的變化,發(fā)現(xiàn)故障發(fā)生時,信號功率譜上一階與二階BPF處能量將發(fā)生變化;于2004年對失速以及喘振現(xiàn)象展開了試驗研究[15],發(fā)現(xiàn)壓氣機進入失速狀態(tài)后,聲信號的幅值增加主要來源于0~0.9BPF范圍的能量增強,如圖8、圖9所示,進入喘振狀態(tài)后,時域上信號的平穩(wěn)性消失,在頻譜上可以觀察到對應(yīng)低頻的位置出現(xiàn)了大量能量聚集,同時還發(fā)現(xiàn),相比于振動傳感器,傳聲器對失速引起的氣流擾動更加敏感。受到Aretakis以及Lawless等工作的啟發(fā),2007年,Morini等[16]在費拉拉大學(xué)的軸流-離心混合式的多級壓氣機試驗臺開展研究,使用如圖10所示的測試布局,捕捉到了喘振發(fā)生前后,聲信號中類似的低頻成分變化。

    圖7 Aretakis等的測量方案[15]

    圖8 聲信號隨壓氣機工作點的變化[15]

    圖9 聲信號頻譜隨工作點的變化[15]

    圖10 Morini等的測量方案[16]

    2.2 基于現(xiàn)代信號分析理論的表征方法

    僅從有限的時域參數(shù)和特征頻率角度對失穩(wěn)聲信號進行描述并不全面:一方面,不同型號的壓氣機甚至同一壓氣機在不同轉(zhuǎn)速下的信號幅值和失穩(wěn)特征頻率都是不一致的,難以發(fā)展出可靠的故障診斷方法;另一方面,對于失穩(wěn)診斷的需求正逐漸從診斷失穩(wěn)本身過渡到失穩(wěn)先兆診斷,失穩(wěn)先兆的特征一般十分微弱,難以直接從時頻分析中獲得。一些先進的信號分析手段,特別是各種高階統(tǒng)計量、信號分解與降噪方法,能夠多角度、更準確地捕捉微弱的信號特征。有望對失穩(wěn)先兆進行更加深入的表征。

    Munari等[17]使用循環(huán)平穩(wěn)分析的方法分析了喘振的聲信號。結(jié)果表明,正常運轉(zhuǎn)的壓氣機信號具有較強的隨機性,當壓氣機進入喘振狀態(tài)時,信號隨機性消失,轉(zhuǎn)變?yōu)橛蓡我活l率主導(dǎo)。Ferrari等[18]對T100微型燃機進行了失速喘振的信號測量,通過定義用于分析并量化振動與聲信號的不平穩(wěn)性的包絡(luò)與方差函數(shù),有效識別了喘振信號。在之后的研究中[19],進一步運用自相關(guān)函數(shù)與高階譜方法對失穩(wěn)狀態(tài)的聲信號進行了分析。對自相關(guān)函數(shù)進行分析發(fā)現(xiàn),接近失穩(wěn)狀態(tài)時,傳聲器可以檢測到聲共振激發(fā)的壓氣機強迫響應(yīng)現(xiàn)象,其表現(xiàn)為信號自相關(guān)函數(shù)外側(cè)趨勢的改變。對信號雙譜分析的結(jié)果表明,不穩(wěn)定狀態(tài)下低頻部分能量相較于穩(wěn)態(tài)有明顯增強,雙譜分析還揭示了系統(tǒng)中存在調(diào)制現(xiàn)象。

    曹昳劼等[20]基于快速傅里葉變換和小波分析等方法,實現(xiàn)了喘振特征信號的提取,得到了與Lawless、Aretakis等類似的結(jié)論。金剛[21]比較了壓力脈動信號、聲信號和振動信號在失穩(wěn)檢測效果上的差異,發(fā)現(xiàn)動態(tài)壓力和聲壓信號對氣動失穩(wěn)的反應(yīng)更加靈敏。Sun等[22]也通過判斷特征頻率是否出現(xiàn)在聲壓信號頻譜上,分析了聲學(xué)測量監(jiān)測壓氣機不穩(wěn)定狀態(tài)的能力,指出傳聲器的位置因素對能否檢測到氣動失穩(wěn)影響較大。Guan等[23]研究了聲壓信號的IMF分量和POD能量占比,發(fā)現(xiàn)壓氣機失穩(wěn)時,部分內(nèi)涵模態(tài)分量會發(fā)生大幅階躍。閆思齊等[24]使用如圖11所示的測試布局采集了臺架喘振試驗時壓氣機左右兩側(cè)的聲壓信號,應(yīng)用小波低頻重構(gòu)聲壓信號實現(xiàn)了對某型渦軸發(fā)動機喘振的識別與特征提取。

    圖11 閆思齊等的傳感器布置方案[24]

    2.3 基于空間模態(tài)的表征方法

    除了上述特征,空間模態(tài)信息對壓氣機的一些失穩(wěn)現(xiàn)象研究也尤為重要[25],王良鋒等[26]對某單級風(fēng)扇的管道周向聲模態(tài)進行了試驗測量,比較了信號互相關(guān)和傳統(tǒng)的均方根兩種模態(tài)分解方法的結(jié)果。梁東等[27]發(fā)展了當測點數(shù)目不滿足條件或測點損壞時,求解出聲模態(tài)的通用分析方法。Sack等[28]提出了基于神經(jīng)網(wǎng)絡(luò)的管道模態(tài)分解方案。程禮等[25]提出了聲模態(tài)時序圖的繪制方法,發(fā)現(xiàn)后行聲模態(tài)波與壓氣機效率下降故障有關(guān),驗證了聲模態(tài)診斷方法的可行性。西安交通大學(xué)與中國航發(fā)四川燃氣渦輪研究院基于圖12所示的測試系統(tǒng)開展了基于聲學(xué)測量的故障診斷、先兆識別和信號處理的研究。在聲模態(tài)方面,提出了基于壓縮感知的聲模態(tài)重構(gòu)方法[29],實現(xiàn)了少數(shù)傳聲器對主導(dǎo)聲模態(tài)的階數(shù)識別和幅值重構(gòu)。在故障先兆識別方面,開展了多級風(fēng)扇試驗器由穩(wěn)定旋轉(zhuǎn)至喘振的瞬態(tài)試驗[30],通過在壓氣機進口機匣安裝的周向傳聲器陣列獲得故障信號,觀測到了風(fēng)扇在進入喘振前,強烈單音噪聲和失速氣團的特征頻率,獲得了與Enrico等觀察到的現(xiàn)象相吻合的結(jié)果。近期,研究人員還使用類似手段對氣動彈性問題開展了研究[31],通過對聲信號的分析,提取到了強迫響應(yīng)、聲共振導(dǎo)致的顫振等故障的特征參數(shù)。

    圖12 測試系統(tǒng)示意圖[29]

    3 氣動失穩(wěn)及其先兆的診斷與預(yù)警方法

    壓氣機的氣動失穩(wěn)及其先兆常常表現(xiàn)為特定規(guī)律的壓力波動。這種波動可以被脈動壓力傳感器或傳聲器捕捉到,基于這兩類信號的氣動失穩(wěn)診斷方法具有一定的相似性。不同于脈動壓力傳感器,由于傳聲器測得的是壓力的動態(tài)量,無法獲得壓力本身的變化,因此基于聲信號的診斷方法更加關(guān)注時域的波動,如信號的變化率、方差、相關(guān)系數(shù)等。另一些方法則基于信號的頻域特點[32],如中心頻率法、能量法等。不同方法的差異主要體現(xiàn)在診斷參數(shù)和報警邏輯上。一般基于時域特征參數(shù)的方法也會通過濾波等預(yù)處理手段引入頻率信息。在已發(fā)表的文獻中,還有很多工作使用不同信號處理分析手段提取出了失穩(wěn)信號的特征,但尚未形成完整的診斷方法,例如自相關(guān)函數(shù)[19,33]、小波分析[34-37]、信號熵[38-39]、聲學(xué)模態(tài)分解[25,31],以及高階譜方法[19],這些方法參見對應(yīng)的文獻。

    3.1 基于信號變化率的診斷方法

    壓氣機進入失穩(wěn)狀態(tài)時氣流參數(shù)的急劇變化可以在信號時域被直觀地捕捉到。雷杰等[40]提出了一種基于脈動壓力變化率的喘振檢測方法,這種方法首先使用式(1)提取出脈動壓力中的動態(tài)量,獲得類似聲壓的結(jié)果ΔP1。喘振發(fā)生時,ΔP1會發(fā)生急劇變化,可以通過式(2)來監(jiān)測這種變化并給出檢測結(jié)果。

    (1)

    |d(ΔP1)/dt|>A

    (2)

    式中:P1(i)為傳感器采集到的壓力信號;N為計算均值所使用的點數(shù);k為當前時刻;A為檢測閾值。這種方法計算簡單,適用于在線診斷,但需要開展大量試驗工作,確定各種狀態(tài)下脈動壓力變化率的分布規(guī)律以及診斷閾值。

    3.2 基于信號方差的診斷方法

    方差可以反映一段時間內(nèi)信號圍繞其平均值的波動量,而波動量是壓氣機失穩(wěn)時信號時域上最顯著的特征。失穩(wěn)狀態(tài)下的聲壓波動也可以使用統(tǒng)計方差[41-42]、延遲方差[39]、自協(xié)方差[21]等來度量,劉勃等[43]提出了一種改進的方差方法。

    (3)

    (4)

    (5)

    自協(xié)方差是描述平穩(wěn)隨機信號統(tǒng)計特征最常用的二階統(tǒng)計量。對于平穩(wěn)隨機信號,自協(xié)方差函數(shù)等于去除了均值部分的自相關(guān)函數(shù),因此自協(xié)方差受壓氣機工況變化影響較小。對于k時刻前采集到的長度為N的信號P(i),可以使用式(6)計算時差為m時的自協(xié)方差,并根據(jù)方差范圍給出失穩(wěn)診斷結(jié)果。

    (6)

    (7)

    3.3 基于信號相關(guān)性的診斷方法

    變化率與方差主要關(guān)注失穩(wěn)信號的突變與波動范圍是否在閾值之內(nèi),但單一的閾值往往缺乏通用性,并且這2個指標對于特征微弱的失穩(wěn)先兆響應(yīng)并不靈敏。信號相似度方法將2個子段信號的相關(guān)性指標作為檢測參數(shù),在一定程度上解決了該問題。相似度越大,說明傳感器采集到的壓力波信號具有越好的時間周期性,壓氣機工作穩(wěn)定。李琳等[45-46]較早地驗證了這種方法在失速先兆檢測方面的可行性。李繼超等[47]提出了比較不同位置信號相似度的方法,進一步獲取了失穩(wěn)起始點的位置信息。毛發(fā)金[48]使用概率分布函數(shù)改進了這一方法。Li等[49]對這種方法背后的原理做出了詳細的闡述。Dong等[50]基于這種方法評估了處理機匣對壓氣機穩(wěn)定性的提升。

    對于單個傳感器采集到的聲壓信號P,t時刻的信號相似度R(t)可以使用式(8)計算:

    (8)

    式中:dwnd為子段長度;tshaft為所要比較的2個子段的時間間隔。余弦相似度方法對信號的周期性較為敏感,李琳等[46]建議使用一個轉(zhuǎn)子周期作為時間間隔,子段長度可以根據(jù)計算量要求從一個葉片周期到一個轉(zhuǎn)子周期之間選擇[51]。

    互相關(guān)方法使用在壓氣機周向布置的多個傳感器同一時刻的余弦相似度作為檢測參數(shù)。當壓氣機進入失穩(wěn)狀態(tài)時,失速先兆會導(dǎo)致不同位置傳感器信號的差異。對于來自通道X與Y的2個信號,余弦相似度為

    (9)

    為了降低正常狀態(tài)下相似度波動導(dǎo)致的誤報,可以使用累計分布函數(shù),將一段時間內(nèi)相似度超過閾值的概率作為檢測指標,如式(10)所示,相較于R值,接近失穩(wěn)狀態(tài)時,F(Rth)上升得更加平穩(wěn)。

    F(Rth)=P(R≤Rth)

    (10)

    基于相關(guān)性的方法計算方便,通用性更強。據(jù)報道,其已經(jīng)在多次高速和低速壓氣試驗中得到成功應(yīng)用[48,52-53],并且由于失穩(wěn)先兆出現(xiàn)時也會引發(fā)相關(guān)性的波動,使得這種方法在預(yù)警速度上具有明顯優(yōu)勢。但是當子段間隔或長度選擇不當,導(dǎo)致相位發(fā)生變化時,相似度可能發(fā)生波動。

    3.4 基于統(tǒng)計分布特征的診斷方法

    對于一些方法閾值難以確定的問題,李長征等[54]發(fā)現(xiàn)壓氣機出口總壓脈動的有效值服從正態(tài)分布,并據(jù)此提出了基于統(tǒng)計特征的喘振檢測方法,其流程如圖13所示,首先根據(jù)正常運行時的數(shù)據(jù),估計信號的正態(tài)分布參數(shù)μ、σ,通過判斷當前壓力值是否超過μ±Cthσ給出預(yù)警。申世才等[55]發(fā)現(xiàn)轉(zhuǎn)速變化過程中的壓力參數(shù)波動會導(dǎo)致分布參數(shù)估計失敗,減小用于參數(shù)估計的樣本容量可以解決這一問題。隨后李長征等[56]又提出了一種聯(lián)合時域頻域特征喘振檢測方法,在時域上計算信號的方差,在頻域上計算信號的中心頻率,借助二維正態(tài)分布函數(shù)對這兩個參數(shù)開展檢測。

    圖13 基于統(tǒng)計分布特征的失穩(wěn)檢測方法[56]

    3.5 基于頻域能量分布的診斷方法

    借助頻譜分析信號的周期性或頻率特性進行檢測是故障檢測常用的手段。Aretakis等[15]發(fā)現(xiàn)0~0.9BPF頻率范圍的聲壓信號均方根幅值σ0.9X隨著壓氣機流量的降低呈現(xiàn)單調(diào)增加的趨勢。因此可以確定出區(qū)分穩(wěn)定與失速、失速與喘振的閾值σL1與σL2,并建立起轉(zhuǎn)速與σL1、σL2的關(guān)系。通過判斷σ0.9X與當前轉(zhuǎn)速對應(yīng)的σL1、σL2的大小,給出失穩(wěn)預(yù)警,流程如圖14所示。李長征等[57]提出了基于低頻成分瞬時能量的失穩(wěn)檢測方法。杜軍等[58]提出根據(jù)低頻能量與其余能量的比例檢測失穩(wěn)的方法,如圖15所示,實現(xiàn)了對不同換算轉(zhuǎn)速下的軸流壓氣機喘振故障的診斷。這種方法無須設(shè)置閾值,更加適用于工程應(yīng)用。

    圖14 基于RMS的壓氣機故障診斷方法[15]

    圖15 基于聲壓能量比的壓氣機故障診斷方法[58]

    3.6 基于數(shù)據(jù)驅(qū)動的診斷方法

    基于數(shù)據(jù)驅(qū)動的各種故障識別方法已經(jīng)在很多場景取得了良好的效果,并有望應(yīng)用于壓氣機氣動失穩(wěn)故障的監(jiān)測甚至預(yù)警,特別是失穩(wěn)先兆監(jiān)測這種微小特征識別任務(wù)。但由于試驗成本高昂、數(shù)據(jù)難以積累等原因,現(xiàn)有的研究大多基于低速風(fēng)扇/壓氣機試驗器開展。Schulze等[59]基于離心壓氣機噪聲信號的轉(zhuǎn)頻能量、葉片通過頻率能量等特征,使用人工免疫算法對氣動失穩(wěn)故障進行了診斷,實現(xiàn)了對于氣動失穩(wěn)的提前預(yù)警。

    Wang等[60-63]進行了數(shù)據(jù)驅(qū)動方法的一系列代表性研究工作,對包含外部擾動的微小且快速變化的失速動態(tài)信號進行了RBF神經(jīng)網(wǎng)絡(luò)建模和動態(tài)模式快速識別,通過比較當前信號與不同失穩(wěn)模式庫的殘差,實現(xiàn)了對短尺度Spike型失速的提前預(yù)警。另一些研究工作更聚焦于端到端學(xué)習(xí),即直接將原始信號作為診斷模型的輸入,將診斷結(jié)果作為輸出。Hipple等[64]使用LSTM網(wǎng)絡(luò)基于多種傳感器信號開展了失速預(yù)測研究,指出使用回歸模型預(yù)測壓氣機的失穩(wěn),將失穩(wěn)概率作為網(wǎng)絡(luò)的輸出獲得的結(jié)果更加準確。Quan等[65-66]提出了一種用于失穩(wěn)預(yù)測的LR-WaveNet模型,模型的輸入為時序數(shù)據(jù)和時域特征,輸出為當前時刻對應(yīng)的失穩(wěn)概率,并通過與LR模型、WaveNet模型、LSTM模型進行對比驗證了所提出方法的預(yù)測精度與穩(wěn)定性較好。

    4 結(jié)束語

    聲信號中包含了大量關(guān)于壓氣機工作狀態(tài)的信息。通過對聲信號的深入研究分析,有望研制出更加準確、靈敏的壓氣機故障診斷、監(jiān)測和預(yù)警方法。通過對以上研究工作的介紹,主要得到以下結(jié)論。

    ① 壓氣機的各種氣動失穩(wěn)形式都可以在聲信號上得到表征。相比于機匣外側(cè),安裝在進氣道靠近氣流位置的傳聲器對氣動失穩(wěn)故障的響應(yīng)更加靈敏,信噪比更高。壓氣機進入氣動失穩(wěn)時通常伴隨著聲信號幅值的增大,而這種幅值增大通常有氣動失穩(wěn)相關(guān)的先兆特征,如模態(tài)波等。

    ② 基于上述信號變化特征,可以使用聲信號對壓氣機的失穩(wěn)故障進行監(jiān)測與預(yù)測,目前常用的方法主要關(guān)注聲信號時域幅值特征、頻段能量特征和信號周期性特征,基本實現(xiàn)了對于失穩(wěn)的監(jiān)測,但是對于特征更加微弱的失穩(wěn)先兆監(jiān)測能力有所不足。先進的信號處理方法,例如各種模態(tài)分解算法、智能方法有望為失穩(wěn)預(yù)警提供有效手段。

    ③ 目前對于壓氣機聲學(xué)測量的研究仍存在一些局限,限制了聲學(xué)方法作為能夠在線診斷、健康監(jiān)測的數(shù)據(jù)來源的可能。這種局限主要來源于:很多聲學(xué)試驗都是在實驗室或露天試車臺環(huán)境中開展,伴隨著嚴格的消聲措施與混響控制,而實際飛行條件下測得的聲信號可能包含了大量的噪聲與干擾。為解決這一問題,可能需要在傳感器的選擇與布局、聲源定位、降噪方法等方面開展進一步研究,提高信號的信噪比。

    ④ 對于數(shù)據(jù)驅(qū)動的失穩(wěn)診斷方法數(shù)據(jù)量不足的問題,可能需要在充分認識氣動失穩(wěn)機理的基礎(chǔ)上,探究借助數(shù)字孿生等手段,人工構(gòu)造失穩(wěn)信號數(shù)據(jù)。盡管廣泛收集來自各種型號的數(shù)據(jù),構(gòu)造通用的失穩(wěn)監(jiān)測預(yù)測模型較為困難,但對于某些重點型號,針對性地大量搜集試驗數(shù)據(jù),構(gòu)造數(shù)據(jù)集并訓(xùn)練準確率更高的失穩(wěn)監(jiān)測與預(yù)警模型是有可能實現(xiàn)的。

    猜你喜歡
    壓氣機聲學(xué)氣動
    中寰氣動執(zhí)行機構(gòu)
    軸流壓氣機效率評定方法
    基于NACA0030的波紋狀翼型氣動特性探索
    重型燃氣輪機壓氣機第一級轉(zhuǎn)子葉片斷裂分析
    愛的就是這股Hi-Fi味 Davis Acoustics(戴維斯聲學(xué))Balthus 70
    壓氣機緊湊S形過渡段內(nèi)周向彎靜子性能數(shù)值計算
    Acoustical Treatment Primer:Diffusion談?wù)劼晫W(xué)處理中的“擴散”
    Acoustical Treatment Primer:Absorption談?wù)劼晫W(xué)處理中的“吸聲”(二)
    基于反饋線性化的RLV氣動控制一體化設(shè)計
    Acoustical Treatment Primer:Absorption 談?wù)劼晫W(xué)處理中的“吸聲”
    黄片无遮挡物在线观看| 婷婷色av中文字幕| 亚洲av在线观看美女高潮| 丝袜脚勾引网站| 亚洲情色 制服丝袜| 男男h啪啪无遮挡| a级毛片黄视频| 少妇熟女欧美另类| 国精品久久久久久国模美| 欧美日韩国产mv在线观看视频| 国产亚洲精品久久久com| 五月天丁香电影| 国产成人精品婷婷| 亚洲综合色惰| 欧美精品高潮呻吟av久久| 精品视频人人做人人爽| 丝瓜视频免费看黄片| 久久精品久久久久久久性| 黄色 视频免费看| 成人亚洲欧美一区二区av| 乱码一卡2卡4卡精品| 在线亚洲精品国产二区图片欧美| 国产精品一二三区在线看| 日韩av免费高清视频| 99热全是精品| 在现免费观看毛片| a级毛片黄视频| 亚洲欧美中文字幕日韩二区| 亚洲高清免费不卡视频| 国产免费一区二区三区四区乱码| 欧美日韩视频高清一区二区三区二| 国产精品免费大片| 啦啦啦中文免费视频观看日本| 国产成人av激情在线播放| 成人国产av品久久久| 婷婷成人精品国产| 国产爽快片一区二区三区| 人妻 亚洲 视频| 看免费成人av毛片| 免费高清在线观看日韩| 男女啪啪激烈高潮av片| 18在线观看网站| 哪个播放器可以免费观看大片| 男女免费视频国产| 大话2 男鬼变身卡| 亚洲激情五月婷婷啪啪| 亚洲第一区二区三区不卡| 亚洲av日韩在线播放| 精品国产露脸久久av麻豆| 中文字幕人妻丝袜制服| 国产永久视频网站| 大香蕉久久成人网| 亚洲精品美女久久av网站| 伦精品一区二区三区| 涩涩av久久男人的天堂| 美女中出高潮动态图| 97超碰精品成人国产| 亚洲成av片中文字幕在线观看 | 婷婷色综合大香蕉| 国产精品一区二区在线不卡| 久久久亚洲精品成人影院| 亚洲,欧美,日韩| 18禁动态无遮挡网站| 97人妻天天添夜夜摸| 国产在视频线精品| 国产亚洲精品久久久com| 午夜激情av网站| 99九九在线精品视频| 亚洲综合色惰| 国产精品熟女久久久久浪| 我的女老师完整版在线观看| 中文字幕人妻熟女乱码| 亚洲国产av影院在线观看| 亚洲精品视频女| 少妇被粗大猛烈的视频| 制服丝袜香蕉在线| 这个男人来自地球电影免费观看 | 久久久国产一区二区| 国产熟女午夜一区二区三区| 免费高清在线观看视频在线观看| 国产黄色免费在线视频| 观看av在线不卡| 亚洲美女视频黄频| 国产精品一国产av| 五月玫瑰六月丁香| 99热网站在线观看| 亚洲色图 男人天堂 中文字幕 | 亚洲av中文av极速乱| 免费少妇av软件| 少妇被粗大猛烈的视频| 女人久久www免费人成看片| 波野结衣二区三区在线| 精品一区二区三卡| 丰满饥渴人妻一区二区三| 国产片特级美女逼逼视频| 亚洲成色77777| 男人舔女人的私密视频| 免费少妇av软件| 欧美人与善性xxx| 久久久国产欧美日韩av| 日韩电影二区| 中文字幕制服av| 少妇人妻精品综合一区二区| 黄片播放在线免费| 啦啦啦中文免费视频观看日本| 日韩av在线免费看完整版不卡| 亚洲激情五月婷婷啪啪| 国产欧美另类精品又又久久亚洲欧美| 国产成人av激情在线播放| 大香蕉97超碰在线| 丝袜在线中文字幕| 国产亚洲精品第一综合不卡 | 婷婷色综合www| 亚洲欧美清纯卡通| 国产av国产精品国产| 宅男免费午夜| 色94色欧美一区二区| 国产成人精品无人区| 亚洲熟女精品中文字幕| 欧美日韩精品成人综合77777| 日韩av在线免费看完整版不卡| 丝袜人妻中文字幕| 亚洲av欧美aⅴ国产| 精品福利永久在线观看| 最新中文字幕久久久久| 亚洲av.av天堂| 国产亚洲欧美精品永久| 免费久久久久久久精品成人欧美视频 | 免费观看性生交大片5| 亚洲 欧美一区二区三区| 亚洲综合精品二区| 最后的刺客免费高清国语| 欧美3d第一页| 美女视频免费永久观看网站| 国产不卡av网站在线观看| 国产av码专区亚洲av| 69精品国产乱码久久久| 亚洲欧洲日产国产| 午夜免费鲁丝| 午夜老司机福利剧场| 18禁在线无遮挡免费观看视频| 一本大道久久a久久精品| 丰满饥渴人妻一区二区三| 久久久久久久大尺度免费视频| 久久99热6这里只有精品| 国产精品国产三级国产专区5o| 黄色 视频免费看| 亚洲国产成人一精品久久久| 国产老妇伦熟女老妇高清| 成年人午夜在线观看视频| 美国免费a级毛片| 不卡视频在线观看欧美| 91国产中文字幕| 搡女人真爽免费视频火全软件| 国产熟女午夜一区二区三区| 中国三级夫妇交换| 色哟哟·www| 亚洲av电影在线观看一区二区三区| 午夜免费鲁丝| 黄色视频在线播放观看不卡| 国产欧美日韩一区二区三区在线| 亚洲欧美日韩另类电影网站| freevideosex欧美| 免费不卡的大黄色大毛片视频在线观看| 国产精品欧美亚洲77777| 中文字幕av电影在线播放| 女人被躁到高潮嗷嗷叫费观| 在现免费观看毛片| 亚洲成人手机| 久久久久久人人人人人| 涩涩av久久男人的天堂| 久久影院123| 精品国产国语对白av| 国产高清三级在线| 寂寞人妻少妇视频99o| 国产一区二区在线观看日韩| 校园人妻丝袜中文字幕| 色婷婷久久久亚洲欧美| 天天影视国产精品| 午夜福利网站1000一区二区三区| 国产精品秋霞免费鲁丝片| 黄色毛片三级朝国网站| 内地一区二区视频在线| 亚洲精品美女久久av网站| 日韩制服骚丝袜av| 超色免费av| 男人操女人黄网站| 少妇的逼水好多| 亚洲人成77777在线视频| 日本av免费视频播放| 国产精品免费大片| 亚洲欧洲精品一区二区精品久久久 | 一级爰片在线观看| 国产午夜精品一二区理论片| 18禁裸乳无遮挡动漫免费视频| 永久免费av网站大全| 久久久久久伊人网av| 久久 成人 亚洲| 免费黄网站久久成人精品| 欧美激情极品国产一区二区三区 | 久久热在线av| 成年人午夜在线观看视频| 日本91视频免费播放| 久久久久国产网址| 国产亚洲午夜精品一区二区久久| 亚洲美女黄色视频免费看| 免费不卡的大黄色大毛片视频在线观看| 日本爱情动作片www.在线观看| 国产 精品1| 午夜影院在线不卡| 国产精品一区二区在线观看99| 国产极品天堂在线| 国产黄色免费在线视频| 51国产日韩欧美| 国产成人av激情在线播放| 母亲3免费完整高清在线观看 | 欧美 亚洲 国产 日韩一| 飞空精品影院首页| 免费黄网站久久成人精品| 在线天堂中文资源库| 亚洲美女黄色视频免费看| 久久久精品94久久精品| 美女福利国产在线| 青春草国产在线视频| 日本欧美视频一区| 免费观看在线日韩| 国产精品偷伦视频观看了| 色吧在线观看| 国产精品三级大全| 亚洲精华国产精华液的使用体验| 久久av网站| 超色免费av| 国产精品免费大片| 亚洲国产色片| 亚洲国产看品久久| 成年美女黄网站色视频大全免费| 超色免费av| 日日摸夜夜添夜夜爱| 日韩制服丝袜自拍偷拍| 亚洲经典国产精华液单| 久久国产精品大桥未久av| 汤姆久久久久久久影院中文字幕| 日本猛色少妇xxxxx猛交久久| 国产综合精华液| 亚洲精品久久午夜乱码| 成年人午夜在线观看视频| 久久ye,这里只有精品| 中文字幕免费在线视频6| 国产精品国产三级国产专区5o| 久久狼人影院| 母亲3免费完整高清在线观看 | 免费黄网站久久成人精品| 久久久精品免费免费高清| 亚洲第一区二区三区不卡| 亚洲综合色惰| av卡一久久| 欧美 日韩 精品 国产| 男女午夜视频在线观看 | 亚洲 欧美一区二区三区| 欧美 亚洲 国产 日韩一| 韩国精品一区二区三区 | 国产成人91sexporn| 国产av国产精品国产| 黄色 视频免费看| 在线观看人妻少妇| 国产亚洲欧美精品永久| 国产在视频线精品| 国产精品国产av在线观看| 一区在线观看完整版| 一级毛片我不卡| 国产一级毛片在线| 成年女人在线观看亚洲视频| 欧美精品一区二区大全| 美女国产视频在线观看| 亚洲精品国产av成人精品| 国产高清三级在线| 亚洲欧美一区二区三区黑人 | 在线天堂中文资源库| 日韩中字成人| 国产日韩一区二区三区精品不卡| 另类精品久久| 久久久久久久国产电影| 女人久久www免费人成看片| 中文字幕人妻熟女乱码| 人妻人人澡人人爽人人| 在线观看一区二区三区激情| 免费观看在线日韩| 中国三级夫妇交换| 女性被躁到高潮视频| 亚洲三级黄色毛片| 亚洲av日韩在线播放| 交换朋友夫妻互换小说| 90打野战视频偷拍视频| 亚洲精品,欧美精品| 日韩中字成人| 777米奇影视久久| 黄色视频在线播放观看不卡| 国产又色又爽无遮挡免| av网站免费在线观看视频| 免费不卡的大黄色大毛片视频在线观看| 高清欧美精品videossex| 一个人免费看片子| 国产免费福利视频在线观看| 丰满迷人的少妇在线观看| 99久久人妻综合| 你懂的网址亚洲精品在线观看| 亚洲av福利一区| 久久精品夜色国产| 免费看不卡的av| 久久精品久久久久久噜噜老黄| 日本av手机在线免费观看| 老熟女久久久| av.在线天堂| 欧美人与性动交α欧美精品济南到 | 夫妻性生交免费视频一级片| 极品人妻少妇av视频| 最近中文字幕高清免费大全6| 啦啦啦啦在线视频资源| 国产成人午夜福利电影在线观看| 一区二区三区精品91| 国产视频首页在线观看| 免费看av在线观看网站| 午夜福利乱码中文字幕| 国产成人精品无人区| 一级毛片 在线播放| 女性被躁到高潮视频| av在线老鸭窝| 你懂的网址亚洲精品在线观看| 国产一区亚洲一区在线观看| 黄片播放在线免费| 丝袜在线中文字幕| 亚洲国产日韩一区二区| 日本午夜av视频| 午夜老司机福利剧场| 亚洲精品一区蜜桃| 一级毛片黄色毛片免费观看视频| av免费在线看不卡| 下体分泌物呈黄色| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美3d第一页| 人人妻人人澡人人爽人人夜夜| 国产精品一国产av| 观看美女的网站| 黄色 视频免费看| 亚洲少妇的诱惑av| 国产亚洲精品久久久com| 天天影视国产精品| 成年美女黄网站色视频大全免费| 亚洲伊人久久精品综合| 亚洲国产精品成人久久小说| 久久精品aⅴ一区二区三区四区 | 亚洲精品日本国产第一区| 蜜桃在线观看..| 美女国产视频在线观看| 五月天丁香电影| 另类亚洲欧美激情| 一区在线观看完整版| 免费av中文字幕在线| 看非洲黑人一级黄片| 制服诱惑二区| 美女大奶头黄色视频| 男女高潮啪啪啪动态图| 欧美精品亚洲一区二区| 久久亚洲国产成人精品v| av网站免费在线观看视频| 青春草亚洲视频在线观看| 欧美亚洲日本最大视频资源| 90打野战视频偷拍视频| 少妇的逼好多水| 麻豆精品久久久久久蜜桃| 国产欧美亚洲国产| 国产精品免费大片| 日韩一区二区三区影片| 赤兔流量卡办理| 亚洲成人手机| 免费少妇av软件| av播播在线观看一区| 99九九在线精品视频| 久久精品aⅴ一区二区三区四区 | 九色亚洲精品在线播放| 久久精品国产亚洲av涩爱| 成人国产av品久久久| 麻豆精品久久久久久蜜桃| 九九爱精品视频在线观看| 亚洲经典国产精华液单| 99久久中文字幕三级久久日本| 久久久久网色| 精品一品国产午夜福利视频| 日韩人妻精品一区2区三区| 中文字幕最新亚洲高清| av在线app专区| 狠狠婷婷综合久久久久久88av| 夫妻午夜视频| 侵犯人妻中文字幕一二三四区| 欧美成人精品欧美一级黄| 日本wwww免费看| 精品卡一卡二卡四卡免费| 精品视频人人做人人爽| 男女边摸边吃奶| 午夜免费鲁丝| 午夜福利视频在线观看免费| 人成视频在线观看免费观看| 国产精品一区二区在线不卡| 九九爱精品视频在线观看| 成人亚洲精品一区在线观看| 久久这里只有精品19| 最近手机中文字幕大全| 91成人精品电影| 国产一区二区激情短视频 | 久久国产精品男人的天堂亚洲 | 最近最新中文字幕免费大全7| 咕卡用的链子| 高清在线视频一区二区三区| 制服丝袜香蕉在线| 精品久久国产蜜桃| 日韩 亚洲 欧美在线| 99久久人妻综合| 制服丝袜香蕉在线| 国产1区2区3区精品| 十八禁高潮呻吟视频| 在线看a的网站| 国产国拍精品亚洲av在线观看| 亚洲精品国产av成人精品| 国产又爽黄色视频| 国产精品 国内视频| 啦啦啦视频在线资源免费观看| 超碰97精品在线观看| 日韩,欧美,国产一区二区三区| 精品一区二区三区四区五区乱码 | 亚洲精品视频女| 咕卡用的链子| 国产在线视频一区二区| 亚洲精品色激情综合| 熟妇人妻不卡中文字幕| 日韩中字成人| 欧美 亚洲 国产 日韩一| 女人久久www免费人成看片| 9色porny在线观看| 99久久综合免费| 如何舔出高潮| 国产在线视频一区二区| 少妇人妻久久综合中文| 免费观看性生交大片5| 青春草亚洲视频在线观看| 国产精品蜜桃在线观看| 免费av不卡在线播放| 成人国语在线视频| 校园人妻丝袜中文字幕| 在现免费观看毛片| 捣出白浆h1v1| 少妇的逼好多水| 女的被弄到高潮叫床怎么办| 最近最新中文字幕大全免费视频 | 超碰97精品在线观看| 日韩一区二区视频免费看| 香蕉国产在线看| 在线 av 中文字幕| 国产一区二区激情短视频 | 夫妻午夜视频| av电影中文网址| 最近的中文字幕免费完整| 久久午夜综合久久蜜桃| 99热全是精品| 国产精品久久久久久久电影| 国产亚洲午夜精品一区二区久久| 天天躁夜夜躁狠狠躁躁| 亚洲av中文av极速乱| 中文字幕av电影在线播放| av黄色大香蕉| 免费在线观看黄色视频的| 熟女人妻精品中文字幕| 亚洲成色77777| 搡老乐熟女国产| 久久精品aⅴ一区二区三区四区 | 日本91视频免费播放| 大话2 男鬼变身卡| 高清视频免费观看一区二区| 日韩欧美精品免费久久| 男人舔女人的私密视频| 在线看a的网站| 精品卡一卡二卡四卡免费| 免费在线观看完整版高清| 高清黄色对白视频在线免费看| 91国产中文字幕| 国产黄色视频一区二区在线观看| 黄色配什么色好看| 中文字幕亚洲精品专区| 免费看光身美女| 国产男人的电影天堂91| 丝袜在线中文字幕| 一区二区三区精品91| 亚洲欧美中文字幕日韩二区| 亚洲,欧美精品.| 久久久精品94久久精品| 99久久人妻综合| 日韩精品有码人妻一区| 韩国高清视频一区二区三区| 午夜日本视频在线| 精品一区二区三区视频在线| 黑丝袜美女国产一区| 91精品国产国语对白视频| 丝袜脚勾引网站| 在线免费观看不下载黄p国产| 高清视频免费观看一区二区| 一级片'在线观看视频| av线在线观看网站| av福利片在线| 黄网站色视频无遮挡免费观看| 寂寞人妻少妇视频99o| 欧美精品av麻豆av| 97在线人人人人妻| 国产精品秋霞免费鲁丝片| 天堂俺去俺来也www色官网| 插逼视频在线观看| 狂野欧美激情性xxxx在线观看| 天堂俺去俺来也www色官网| 日本免费在线观看一区| 99视频精品全部免费 在线| 中文精品一卡2卡3卡4更新| 精品人妻一区二区三区麻豆| 黄片无遮挡物在线观看| 春色校园在线视频观看| 久久久久精品人妻al黑| 国产免费福利视频在线观看| 最后的刺客免费高清国语| 在线精品无人区一区二区三| 欧美另类一区| 在线观看免费视频网站a站| 国产男人的电影天堂91| 中文字幕免费在线视频6| 免费观看av网站的网址| 男女无遮挡免费网站观看| 精品久久久久久电影网| 内地一区二区视频在线| 亚洲激情五月婷婷啪啪| 一级毛片我不卡| 色网站视频免费| 成年人午夜在线观看视频| 一级,二级,三级黄色视频| 久久久精品区二区三区| 美女主播在线视频| 三上悠亚av全集在线观看| 精品久久久久久电影网| 国产亚洲av片在线观看秒播厂| 免费观看a级毛片全部| 人妻系列 视频| 香蕉国产在线看| 欧美亚洲日本最大视频资源| 婷婷色综合大香蕉| 国产乱来视频区| 久久久亚洲精品成人影院| 最近的中文字幕免费完整| 国产精品久久久久久久电影| xxxhd国产人妻xxx| 亚洲久久久国产精品| 视频中文字幕在线观看| 不卡视频在线观看欧美| 男女国产视频网站| 蜜桃在线观看..| 中文字幕最新亚洲高清| 免费看不卡的av| 中文字幕最新亚洲高清| 伦精品一区二区三区| 国产日韩一区二区三区精品不卡| 秋霞伦理黄片| 最近最新中文字幕免费大全7| 亚洲国产成人一精品久久久| 国产黄色免费在线视频| 国产成人aa在线观看| 亚洲精品乱码久久久久久按摩| 国产在线免费精品| 亚洲av男天堂| 精品亚洲成国产av| 桃花免费在线播放| tube8黄色片| 国产国拍精品亚洲av在线观看| 久久精品国产鲁丝片午夜精品| 搡女人真爽免费视频火全软件| 999精品在线视频| 亚洲美女视频黄频| 韩国精品一区二区三区 | 一区二区三区四区激情视频| 不卡视频在线观看欧美| 亚洲精品国产色婷婷电影| 国产日韩欧美在线精品| 免费久久久久久久精品成人欧美视频 | 美女国产高潮福利片在线看| 精品福利永久在线观看| 一二三四在线观看免费中文在 | √禁漫天堂资源中文www| 成年动漫av网址| 亚洲伊人久久精品综合| 一区二区av电影网| 三上悠亚av全集在线观看| www.色视频.com| 久久人人97超碰香蕉20202| 色吧在线观看| 免费少妇av软件| 国产精品一区二区在线不卡| 国产成人一区二区在线| 久久热在线av| 捣出白浆h1v1| 精品国产国语对白av| 午夜福利在线观看免费完整高清在| 婷婷色综合www| 边亲边吃奶的免费视频| 精品一区二区三区视频在线| 纵有疾风起免费观看全集完整版| a 毛片基地| 亚洲美女黄色视频免费看| 超碰97精品在线观看| 最新的欧美精品一区二区| 亚洲国产精品成人久久小说| 18在线观看网站| 自线自在国产av|