• <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é)處理中的“吸聲”
    午夜91福利影院| 国产精品国产三级国产专区5o| 两个人看的免费小视频| 欧美中文综合在线视频| 可以免费在线观看a视频的电影网站 | 伊人久久大香线蕉亚洲五| tube8黄色片| 日韩,欧美,国产一区二区三区| 久久99一区二区三区| 香蕉国产在线看| 女人精品久久久久毛片| 好男人视频免费观看在线| 国产黄色视频一区二区在线观看| 最近2019中文字幕mv第一页| 在线天堂最新版资源| 亚洲国产精品国产精品| 一区二区三区激情视频| 男人爽女人下面视频在线观看| av女优亚洲男人天堂| 欧美日韩国产mv在线观看视频| xxx大片免费视频| 街头女战士在线观看网站| 成年美女黄网站色视频大全免费| 免费在线观看完整版高清| 男女无遮挡免费网站观看| 日日啪夜夜爽| 三上悠亚av全集在线观看| 久久久久久免费高清国产稀缺| 日韩 亚洲 欧美在线| 欧美国产精品va在线观看不卡| 春色校园在线视频观看| 久久精品国产综合久久久| 欧美亚洲日本最大视频资源| 两个人免费观看高清视频| 国产在线一区二区三区精| 国产深夜福利视频在线观看| 少妇的逼水好多| 久久久国产精品麻豆| 中文天堂在线官网| 自拍欧美九色日韩亚洲蝌蚪91| 人人澡人人妻人| 国产精品熟女久久久久浪| 久久久久久人妻| 在线观看免费视频网站a站| 色94色欧美一区二区| 汤姆久久久久久久影院中文字幕| 国产一区亚洲一区在线观看| 美女福利国产在线| 免费看av在线观看网站| kizo精华| 天天躁夜夜躁狠狠久久av| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽av亚洲精品天堂| 下体分泌物呈黄色| 亚洲国产看品久久| 色网站视频免费| 最新的欧美精品一区二区| 爱豆传媒免费全集在线观看| 咕卡用的链子| 丝瓜视频免费看黄片| 久久久a久久爽久久v久久| 久久鲁丝午夜福利片| 永久免费av网站大全| 99香蕉大伊视频| 夫妻午夜视频| 热99国产精品久久久久久7| 伊人久久大香线蕉亚洲五| 亚洲第一青青草原| 1024视频免费在线观看| 天堂8中文在线网| 麻豆乱淫一区二区| 十八禁网站网址无遮挡| 国产精品秋霞免费鲁丝片| 久久午夜福利片| 久久精品国产鲁丝片午夜精品| 亚洲经典国产精华液单| 一本—道久久a久久精品蜜桃钙片| 丰满迷人的少妇在线观看| 少妇猛男粗大的猛烈进出视频| 午夜老司机福利剧场| 亚洲av免费高清在线观看| 日本爱情动作片www.在线观看| 国产黄频视频在线观看| 一区二区三区精品91| 男人操女人黄网站| av在线观看视频网站免费| 亚洲天堂av无毛| 五月伊人婷婷丁香| 又大又黄又爽视频免费| 精品久久久久久电影网| 性色avwww在线观看| 亚洲国产最新在线播放| a级毛片黄视频| 亚洲精品中文字幕在线视频| 纯流量卡能插随身wifi吗| 中文字幕制服av| 亚洲国产精品国产精品| 亚洲国产毛片av蜜桃av| 亚洲国产精品999| 欧美最新免费一区二区三区| 国产成人精品福利久久| av天堂久久9| 午夜激情久久久久久久| 日本av免费视频播放| 大码成人一级视频| 观看美女的网站| 麻豆精品久久久久久蜜桃| 9191精品国产免费久久| 韩国精品一区二区三区| 国产成人精品在线电影| 国产免费福利视频在线观看| 色婷婷av一区二区三区视频| 欧美国产精品一级二级三级| 啦啦啦在线免费观看视频4| 午夜福利视频精品| 香蕉精品网在线| 热99久久久久精品小说推荐| 免费观看在线日韩| 亚洲欧洲国产日韩| 国产高清国产精品国产三级| 夜夜骑夜夜射夜夜干| 国产成人欧美| 在线观看免费视频网站a站| 欧美日韩综合久久久久久| 日韩一本色道免费dvd| 亚洲精品乱久久久久久| 少妇猛男粗大的猛烈进出视频| 亚洲欧美日韩另类电影网站| 久久久a久久爽久久v久久| 岛国毛片在线播放| 最新的欧美精品一区二区| 中文精品一卡2卡3卡4更新| 久久午夜福利片| 亚洲图色成人| 欧美人与善性xxx| 天堂中文最新版在线下载| 只有这里有精品99| 亚洲综合色网址| 久久午夜福利片| 18禁动态无遮挡网站| 熟女av电影| 欧美精品一区二区大全| 亚洲国产av新网站| 亚洲精品久久午夜乱码| 国产精品99久久99久久久不卡 | 亚洲四区av| 熟妇人妻不卡中文字幕| 亚洲,欧美精品.| 天天操日日干夜夜撸| 美女主播在线视频| 久久午夜综合久久蜜桃| a级毛片黄视频| 国产精品麻豆人妻色哟哟久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩大片免费观看网站| 国产1区2区3区精品| 韩国高清视频一区二区三区| 欧美变态另类bdsm刘玥| 亚洲精品aⅴ在线观看| 久久久久久伊人网av| 亚洲 欧美一区二区三区| 亚洲精品第二区| 秋霞伦理黄片| 最近最新中文字幕大全免费视频 | 久久精品国产亚洲av天美| 欧美人与善性xxx| 国产精品 国内视频| 亚洲成国产人片在线观看| 精品国产国语对白av| 亚洲国产成人一精品久久久| 尾随美女入室| freevideosex欧美| 一级毛片电影观看| 老司机影院成人| 波多野结衣av一区二区av| 国产免费现黄频在线看| 亚洲色图 男人天堂 中文字幕| 搡老乐熟女国产| 9热在线视频观看99| 亚洲少妇的诱惑av| 亚洲成色77777| 99精国产麻豆久久婷婷| 精品国产乱码久久久久久小说| 搡老乐熟女国产| 精品少妇一区二区三区视频日本电影 | 国产成人精品婷婷| 嫩草影院入口| 婷婷色av中文字幕| 国产成人午夜福利电影在线观看| 亚洲av日韩在线播放| 最近最新中文字幕免费大全7| 一级片免费观看大全| 成人手机av| 熟女少妇亚洲综合色aaa.| 毛片一级片免费看久久久久| 婷婷色av中文字幕| 一区二区日韩欧美中文字幕| 成年女人在线观看亚洲视频| 国产色婷婷99| 黑丝袜美女国产一区| 永久免费av网站大全| 成人亚洲精品一区在线观看| 亚洲欧美成人精品一区二区| 国产成人精品在线电影| 国产极品天堂在线| 亚洲av.av天堂| 亚洲,一卡二卡三卡| 少妇人妻精品综合一区二区| www.自偷自拍.com| 久久av网站| 亚洲av中文av极速乱| 日本猛色少妇xxxxx猛交久久| 午夜老司机福利剧场| 狠狠精品人妻久久久久久综合| 亚洲精品国产av成人精品| 国产精品欧美亚洲77777| 午夜福利视频精品| 亚洲国产精品成人久久小说| 精品亚洲乱码少妇综合久久| 大香蕉久久网| av在线观看视频网站免费| 国产精品一国产av| 爱豆传媒免费全集在线观看| 一级毛片我不卡| 999久久久国产精品视频| 欧美老熟妇乱子伦牲交| 国产女主播在线喷水免费视频网站| 免费播放大片免费观看视频在线观看| 亚洲综合色惰| 久久国产亚洲av麻豆专区| 国产精品久久久久久久久免| 久久久欧美国产精品| 亚洲欧美精品自产自拍| 国产一区二区在线观看av| 日本av手机在线免费观看| 久久久久国产精品人妻一区二区| 精品一品国产午夜福利视频| 99精国产麻豆久久婷婷| 一区二区av电影网| 在线观看美女被高潮喷水网站| 国产黄色视频一区二区在线观看| 男女无遮挡免费网站观看| 高清黄色对白视频在线免费看| 久久精品久久久久久噜噜老黄| 人人妻人人添人人爽欧美一区卜| 日日啪夜夜爽| 久久久亚洲精品成人影院| 午夜免费鲁丝| 精品视频人人做人人爽| 欧美人与性动交α欧美精品济南到 | 青草久久国产| 午夜91福利影院| 国产又色又爽无遮挡免| 国产精品秋霞免费鲁丝片| 丝袜脚勾引网站| 久久久欧美国产精品| 人人妻人人澡人人看| 国产老妇伦熟女老妇高清| 国产在线一区二区三区精| 久久久国产精品麻豆| 一边亲一边摸免费视频| 999久久久国产精品视频| 男女下面插进去视频免费观看| 久久久久久久久免费视频了| 免费少妇av软件| 日本av免费视频播放| 婷婷色综合www| 国产精品国产av在线观看| 久久久久精品久久久久真实原创| 精品少妇一区二区三区视频日本电影 | 国产1区2区3区精品| 熟妇人妻不卡中文字幕| 日韩一区二区三区影片| 一二三四中文在线观看免费高清| 777久久人妻少妇嫩草av网站| 涩涩av久久男人的天堂| 赤兔流量卡办理| 亚洲精品aⅴ在线观看| 亚洲综合色网址| 97人妻天天添夜夜摸| 国产午夜精品一二区理论片| 久久久久精品人妻al黑| 日本vs欧美在线观看视频| 美女xxoo啪啪120秒动态图| 国产免费视频播放在线视频| 久久精品国产亚洲av天美| 久久精品国产亚洲av高清一级| 精品第一国产精品| 性色av一级| 国产一区二区激情短视频 | 高清av免费在线| 少妇人妻精品综合一区二区| 捣出白浆h1v1| 国产国语露脸激情在线看| 丰满少妇做爰视频| 日韩电影二区| 哪个播放器可以免费观看大片| 美女福利国产在线| 午夜福利视频在线观看免费| 日韩av在线免费看完整版不卡| 91在线精品国自产拍蜜月| 国产乱人偷精品视频| 中文乱码字字幕精品一区二区三区| 天堂俺去俺来也www色官网| 久久精品国产亚洲av涩爱| 亚洲熟女精品中文字幕| 嫩草影院入口| 精品一品国产午夜福利视频| 国产精品一国产av| av国产久精品久网站免费入址| 美女午夜性视频免费| 一级片免费观看大全| 亚洲,一卡二卡三卡| 欧美亚洲日本最大视频资源| 国语对白做爰xxxⅹ性视频网站| 自线自在国产av| 亚洲欧美精品自产自拍| 亚洲av电影在线观看一区二区三区| 一级黄片播放器| 久久久a久久爽久久v久久| 99re6热这里在线精品视频| 国产精品成人在线| 黑人欧美特级aaaaaa片| 18禁动态无遮挡网站| 大陆偷拍与自拍| xxx大片免费视频| 三级国产精品片| 99精国产麻豆久久婷婷| 大片电影免费在线观看免费| 免费黄色在线免费观看| 中文字幕人妻丝袜一区二区 | 老鸭窝网址在线观看| 久久久精品94久久精品| 在线观看www视频免费| av有码第一页| 亚洲综合色网址| 亚洲五月色婷婷综合| 人妻一区二区av| 99热网站在线观看| 日韩一区二区三区影片| videosex国产| 精品一区二区免费观看| 亚洲色图 男人天堂 中文字幕| 久久韩国三级中文字幕| 国产亚洲av片在线观看秒播厂| 美女大奶头黄色视频| 天美传媒精品一区二区| 亚洲av综合色区一区| 日韩熟女老妇一区二区性免费视频| 国产成人免费观看mmmm| 亚洲精品日本国产第一区| a级片在线免费高清观看视频| 免费观看a级毛片全部| 王馨瑶露胸无遮挡在线观看| 美女中出高潮动态图| 两性夫妻黄色片| 亚洲国产最新在线播放| 国产淫语在线视频| 如日韩欧美国产精品一区二区三区| 卡戴珊不雅视频在线播放| 国产精品久久久av美女十八| 精品福利永久在线观看| 国产成人av激情在线播放| 十八禁高潮呻吟视频| 最近最新中文字幕大全免费视频 | 久久av网站| 国产精品亚洲av一区麻豆 | 秋霞在线观看毛片| 精品酒店卫生间| 在线观看三级黄色| av女优亚洲男人天堂| 国产伦理片在线播放av一区| 99久久中文字幕三级久久日本| 一二三四在线观看免费中文在| 欧美+日韩+精品| 国产成人免费观看mmmm| 亚洲精品成人av观看孕妇| 国产一区二区激情短视频 | 老司机亚洲免费影院| 尾随美女入室| 精品久久久久久电影网| 亚洲欧美精品综合一区二区三区 | 少妇人妻久久综合中文| 中文乱码字字幕精品一区二区三区| 美女中出高潮动态图| 男人操女人黄网站| 一区二区三区四区激情视频| 久久久欧美国产精品| 电影成人av| 一级毛片我不卡| 亚洲人成77777在线视频| 成人18禁高潮啪啪吃奶动态图| 国产精品香港三级国产av潘金莲 | av视频免费观看在线观看| 成人国产av品久久久| 国产精品一区二区在线不卡| 啦啦啦啦在线视频资源| 女性生殖器流出的白浆| 国产精品免费大片| 国产无遮挡羞羞视频在线观看| 狠狠婷婷综合久久久久久88av| 久久久久久久久免费视频了| 少妇精品久久久久久久| 在线观看免费日韩欧美大片| 精品卡一卡二卡四卡免费| 精品视频人人做人人爽| 男女边摸边吃奶| 男男h啪啪无遮挡| 一级爰片在线观看| 丝袜人妻中文字幕| 中文字幕精品免费在线观看视频| 国产精品国产av在线观看| 国产成人91sexporn| 狂野欧美激情性bbbbbb| 精品人妻熟女毛片av久久网站| 国产男女内射视频| 熟妇人妻不卡中文字幕| 美女大奶头黄色视频| 高清欧美精品videossex| 亚洲人成77777在线视频| 99久久综合免费| 日本黄色日本黄色录像| 青春草亚洲视频在线观看| 久久狼人影院| 国产高清不卡午夜福利| 国产精品偷伦视频观看了| 欧美 亚洲 国产 日韩一| 国产在线免费精品| 久久久精品94久久精品| av线在线观看网站| 久久久亚洲精品成人影院| 午夜福利视频精品| a级片在线免费高清观看视频| 爱豆传媒免费全集在线观看| 久久人人爽人人片av| 最近中文字幕2019免费版| 人人澡人人妻人| 日韩视频在线欧美| 欧美日韩视频精品一区| 亚洲精品国产色婷婷电影| 人人澡人人妻人| 国产淫语在线视频| 国产精品免费大片| 在线观看免费视频网站a站| 免费久久久久久久精品成人欧美视频| 国产不卡av网站在线观看| 人人妻人人添人人爽欧美一区卜| 欧美最新免费一区二区三区| 婷婷色综合www| 久久久久久久久久久久大奶| 成人二区视频| 免费在线观看黄色视频的| 日本91视频免费播放| 久久精品久久久久久噜噜老黄| 久久人人97超碰香蕉20202| 亚洲欧洲国产日韩| 国产精品久久久久久精品古装| 老司机影院毛片| 欧美xxⅹ黑人| 国精品久久久久久国模美| 精品一品国产午夜福利视频| 18禁裸乳无遮挡动漫免费视频| 飞空精品影院首页| 亚洲美女黄色视频免费看| 国产老妇伦熟女老妇高清| 十八禁网站网址无遮挡| 亚洲国产欧美在线一区| 久久人人爽av亚洲精品天堂| a级毛片在线看网站| 日韩免费高清中文字幕av| 亚洲在久久综合| 国产亚洲午夜精品一区二区久久| 美女中出高潮动态图| 色吧在线观看| 免费播放大片免费观看视频在线观看| 国产成人免费无遮挡视频| 国产伦理片在线播放av一区| 男女国产视频网站| 日韩一区二区视频免费看| 国产精品亚洲av一区麻豆 | 最近的中文字幕免费完整| 人成视频在线观看免费观看| 日本欧美国产在线视频| 狠狠婷婷综合久久久久久88av| 丰满少妇做爰视频| 国产高清不卡午夜福利| 精品国产国语对白av| 波野结衣二区三区在线| 精品少妇内射三级| 大香蕉久久网| 久久国产亚洲av麻豆专区| 免费大片黄手机在线观看| 91在线精品国自产拍蜜月| 国产在线视频一区二区| 日本av手机在线免费观看| 尾随美女入室| 两个人看的免费小视频| 99国产精品免费福利视频| 青青草视频在线视频观看| av国产精品久久久久影院| 国产午夜精品一二区理论片| 欧美另类一区| 国产淫语在线视频| 国产精品人妻久久久影院| 久久久国产欧美日韩av| 日日撸夜夜添| 韩国精品一区二区三区| 观看av在线不卡| 999精品在线视频| 国产精品一国产av| a级毛片黄视频| 久久久久久久国产电影| a级毛片在线看网站| 亚洲欧美一区二区三区国产| 制服丝袜香蕉在线| 免费日韩欧美在线观看| 纵有疾风起免费观看全集完整版| 精品少妇内射三级| 老汉色av国产亚洲站长工具| 亚洲内射少妇av| 美女午夜性视频免费| 一级毛片 在线播放| 99久久人妻综合| 两性夫妻黄色片| 麻豆乱淫一区二区| 亚洲一区二区三区欧美精品| 人妻少妇偷人精品九色| 老鸭窝网址在线观看| 最近的中文字幕免费完整| 母亲3免费完整高清在线观看 | 一级黄片播放器| 最新的欧美精品一区二区| 亚洲三区欧美一区| 久久99蜜桃精品久久| 日韩中文字幕欧美一区二区 | 丝瓜视频免费看黄片| 亚洲精品乱久久久久久| 欧美+日韩+精品| 亚洲综合精品二区| 午夜福利在线观看免费完整高清在| 老司机亚洲免费影院| 日韩熟女老妇一区二区性免费视频| 久久免费观看电影| 亚洲欧美清纯卡通| 亚洲成人av在线免费| 9191精品国产免费久久| 大陆偷拍与自拍| 又黄又粗又硬又大视频| 亚洲av欧美aⅴ国产| 亚洲欧美日韩另类电影网站| av在线观看视频网站免费| 婷婷成人精品国产| 99香蕉大伊视频| 99久久人妻综合| 18禁裸乳无遮挡动漫免费视频| 91成人精品电影| 亚洲一级一片aⅴ在线观看| 日韩电影二区| 亚洲成人手机| 欧美日韩成人在线一区二区| 国产精品久久久av美女十八| 亚洲,一卡二卡三卡| 日韩人妻精品一区2区三区| 免费黄网站久久成人精品| 国产亚洲最大av| 香蕉丝袜av| 精品少妇一区二区三区视频日本电影 | 80岁老熟妇乱子伦牲交| 亚洲国产最新在线播放| 国产精品无大码| 美女中出高潮动态图| 毛片一级片免费看久久久久| 免费高清在线观看视频在线观看| 欧美av亚洲av综合av国产av | 好男人视频免费观看在线| 欧美精品一区二区大全| 超碰97精品在线观看| 久久久久久伊人网av| 国产成人精品一,二区| 欧美少妇被猛烈插入视频| 欧美+日韩+精品| 日韩欧美精品免费久久| 欧美精品一区二区大全| 激情视频va一区二区三区| 免费在线观看黄色视频的| 久久久a久久爽久久v久久| 国产精品欧美亚洲77777| 深夜精品福利| av国产久精品久网站免费入址| videos熟女内射| 婷婷色麻豆天堂久久| 欧美日韩一级在线毛片| 欧美成人精品欧美一级黄| 国产亚洲午夜精品一区二区久久| 99精国产麻豆久久婷婷| 欧美亚洲日本最大视频资源| 男女边吃奶边做爰视频| 国产av码专区亚洲av| 欧美日韩亚洲高清精品| 97在线视频观看| 午夜福利,免费看| av国产久精品久网站免费入址| 老鸭窝网址在线观看| 国产野战对白在线观看| xxxhd国产人妻xxx| av有码第一页| 久久久久网色| 国产一区二区三区av在线| 亚洲国产精品国产精品| 国产精品一区二区在线不卡| 最近最新中文字幕大全免费视频 | 国产精品人妻久久久影院| 女人久久www免费人成看片| 亚洲国产精品成人久久小说| 69精品国产乱码久久久|