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

    航空發(fā)動機振動監(jiān)測與故障診斷技術(shù)研究進展

    2024-04-07 02:21:58胡明輝高金吉江志農(nóng)王維民鄒利民周濤凡云峰王越馮家欣李晨陽
    航空學報 2024年4期
    關(guān)鍵詞:故障診斷發(fā)動機振動

    胡明輝,高金吉,江志農(nóng),王維民,鄒利民,周濤,凡云峰,王越,馮家欣,李晨陽

    1.北京化工大學 高端壓縮機及系統(tǒng)技術(shù)全國重點實驗室,北京 100029

    2.北京化工大學 發(fā)動機健康監(jiān)控及網(wǎng)絡(luò)化教育部重點實驗室,北京 100029

    3.北京化工大學 高端機械裝備健康監(jiān)控與自愈化北京市重點實驗室,北京 100029

    航空發(fā)動機是一個國家科技、工業(yè)和國防實力的重要體現(xiàn),作為“工業(yè)皇冠上的明珠”,航空發(fā)動機在結(jié)構(gòu)設(shè)計、材料和制造工藝及試驗方面具有極高的要求。為保障國家安全和鞏固大國地位,中國自2015年相繼啟動“中國制造2025”計劃、“兩機專項”、國家“十四五規(guī)劃綱要”和“2035年遠景”目標綱要,均明確了加快航空發(fā)動機核心技術(shù)發(fā)展的急迫需求[1]。

    先進航空發(fā)動機不斷朝著高效率、高推重比的方向發(fā)展,使得結(jié)構(gòu)更為復雜,且逐漸趨于輕量化,同時發(fā)動機長期服役于高速、高溫、高壓的極端環(huán)境,導致其關(guān)鍵部件的疲勞損傷、突發(fā)故障等問題頻發(fā)[2]。發(fā)動機故障一旦發(fā)生,會直接危及飛機的飛行安全,極易導致機毀人亡等惡性事故,如圖1[3]所示。因此,因發(fā)動機故障導致飛機動力嚴重不足之前,如果能及時、準確地監(jiān)測診斷出故障原因,對保障飛行安全、降低事故發(fā)生率具有重大意義。

    圖1 疑似發(fā)動機故障導致飛機損毀[3]Fig.1 Aircraft damage caused by suspected engine failure[3]

    從20 世紀50年代開始,歐美各國持續(xù)開展了狀態(tài)監(jiān)測和故障診斷技術(shù)的研究[4-5],發(fā)動機健康管理(EHM)系統(tǒng)也逐漸得到了廣泛關(guān)注。EHM 系統(tǒng)可在發(fā)動機運行過程中對關(guān)鍵部件狀態(tài)進行實時監(jiān)測,通過振動分析、滑油監(jiān)測、性能分析等技術(shù)實現(xiàn)故障預測、診斷和壽命管理,進而向飛行員報告危及飛行安全的早期故障信息,避免二次損傷和飛機事故。

    美國普惠公司首席科學家Volponi[6]于2014年綜述了航空發(fā)動機健康管理的發(fā)展,強調(diào)了預測與健康管理(PHM)技術(shù)對航空發(fā)動機運行安全、事故預防以及運維成本控制的重要作用。從1988年實施的IHPTET 計劃中,美國已將狀態(tài)監(jiān)測和故障診斷系統(tǒng)列為亟待突破的關(guān)鍵技術(shù)。美國1997年設(shè)立的航空安全中長期研究計劃[7],2006年實施的VAATE 計劃、目前正在實施的ATTAM 計劃,PHM 技術(shù)均被列為重要研究領(lǐng)域,將顯著提高發(fā)動機運行安全性和可靠性、改善維修模式作為主要研究目標[8]。

    航空發(fā)達國家發(fā)動機PHM 技術(shù)的研究與應用均已取得了顯著突破。PHM 技術(shù)的最早成功應用可追溯到1977年,TF-41 發(fā)動機監(jiān)測系統(tǒng),將A-7E 飛機的故障率降低了91%,且維修效率平均提高了64%[9];1979年,通用電氣公司為F414-GE-400 型發(fā)動機配裝了IECMS 發(fā)動機狀態(tài)監(jiān)測系統(tǒng),經(jīng)美國海軍應用評估,發(fā)動機異常檢測率為100%、虛警率為0%[10];1987年,普惠公司為F100-PW-200 發(fā)動機研制并安裝了健康管理系統(tǒng),實現(xiàn)了所有同型發(fā)動機的健康監(jiān)測和管理,并增加了發(fā)動機監(jiān)測數(shù)據(jù)趨勢分析及故障診斷功能[11];1990年,渦扇發(fā)動機EJ200 配備的EHM 系統(tǒng)采用機載系統(tǒng)和地面支持系統(tǒng)的綜合健康管理架構(gòu)方案,已具備完善的狀態(tài)監(jiān)測和故障診斷功能[12]。F135 發(fā)動機配裝的健康管理系統(tǒng)延續(xù)了F119 的雙FADEC 和CEDU 模式,并增加了故障預測功能,提出了分層區(qū)域預測健康管理系統(tǒng)的概念,集成先進傳感器、先進算法和智能模型來實時監(jiān)視、診斷、預測和管理發(fā)動機的健康狀態(tài),顯著提高了發(fā)動機使用安全性[13]。近年來,美國又提出了地空一體化的航空發(fā)動機健康管理系統(tǒng)構(gòu)架,如圖2[6]所示。

    圖2 美國提出的EHM 系統(tǒng)框架[6]Fig.2 EHM system framework proposed by the United States[6]

    航空發(fā)達國家的發(fā)動機故障預測與健康管理技術(shù)已日趨成熟并趨向智能化,已成為確保飛機安全的關(guān)鍵支撐技術(shù)。同時,航空工業(yè)發(fā)達國家亦對發(fā)動機PHM 系統(tǒng)的規(guī)范化、標準化提出了一系列要求。1981年,美國汽車工程師協(xié)會(SAE)制定了ARP 1587《航空燃氣渦輪發(fā)動機監(jiān)視系統(tǒng)指南》,基于當時的技術(shù)水平提出了對航空發(fā)動機健康管理的要求[14]。隨著電子和計算機技術(shù)水平的不斷發(fā)展,2007年SAE 修訂出版了ARP 1587B《航空燃氣渦輪發(fā)動機健康管理系統(tǒng)指南》[15],對發(fā)動機健康管理系統(tǒng)進行了全面的定義,同時已將最新的健康管理系統(tǒng)具備的故障預測能力囊括在內(nèi)。2008年頒布的AIR 4061B《典型發(fā)動機健康管理系統(tǒng)功能與飛機系統(tǒng)集成指南》,是發(fā)動機健康管理系統(tǒng)功能與飛機系統(tǒng)集成的實行準則。至今,SAE已陸續(xù)發(fā)布了31 項關(guān)于航空燃氣渦輪發(fā)動機健康管理的系列標準,以指導航空發(fā)動機的狀態(tài)監(jiān)測、故障診斷及健康管理系統(tǒng)的設(shè)計、研制、使用和維護[11]。

    表1[12]總結(jié)了美國發(fā)動機監(jiān)測系統(tǒng)、維修策略與相關(guān)技術(shù)文件的發(fā)展沿革。由此可見,PHM是美國保障飛行安全和新型維修模式變革的主要實現(xiàn)技術(shù),且日趨規(guī)范化[16]。

    與發(fā)達國家相比,中國對航空發(fā)動機PHM技術(shù)的研究與應用起步較晚。在役航空發(fā)動機仍主要采用以機上超限告警監(jiān)測、機下離線飛參數(shù)據(jù)判讀為主的監(jiān)測診斷技術(shù),往往在發(fā)動機部件損傷或結(jié)構(gòu)故障發(fā)展到較嚴重程度才能發(fā)揮作用[12]。

    雖然近年來一些科研院所取得了部分技術(shù)成果,但鮮見實現(xiàn)機載的廣泛應用。由北京飛機維修工程公司、中國東方航空集團有限公司等聯(lián)合研制的發(fā)動機監(jiān)測與診斷系統(tǒng)在40 多臺發(fā)動機上試用成功,可以實現(xiàn)氣路故障的粗略分析,但無法實現(xiàn)故障部件的精確定位[17]。清華大學研制了高速數(shù)據(jù)采集系統(tǒng)、信號處理系統(tǒng)及基于模糊聚類的故障診斷系統(tǒng),通過地面實驗進行了監(jiān)測和故障診斷方法的探索[18]。西北工業(yè)大學的張永峰[19]開發(fā)了適用于飛行試驗的發(fā)動機振動監(jiān)測系統(tǒng),能夠有效實現(xiàn)穩(wěn)態(tài)和過渡態(tài)振動參數(shù)的趨勢分析和譜分析。北京航空工程技術(shù)研究所研制了一套發(fā)動機振動監(jiān)測與分析系統(tǒng),綜合振動總量評估與頻譜分析手段,對某型發(fā)動機開展多點振動監(jiān)測分析,在試車臺、外場通過多次實驗驗證了其主要功能[20]。

    振動監(jiān)測診斷是航空發(fā)動機健康管理最重要的方面之一,本文將從振動監(jiān)測振動系統(tǒng)研制及應用、典型故障特征及診斷方法研究現(xiàn)狀進行具體綜述,并對整機振動類故障診斷技術(shù)的發(fā)展方向進行分析。

    1 整機振動監(jiān)測診斷系統(tǒng)研制及應用

    氣路、振動、滑油PHM 技術(shù)是航空發(fā)動機預測與健康管理系統(tǒng)的關(guān)鍵,這3 種技術(shù)可以基本涵蓋發(fā)動機的主要氣路、轉(zhuǎn)動和傳動部件監(jiān)測分析需求[21]。

    其中,振動類故障是航空發(fā)動機的主要故障模式。國外先進發(fā)動機EJ200、F119、F135 配裝的健康管理系統(tǒng)均包含振動診斷技術(shù),并且隨著發(fā)動機健康管理系統(tǒng)的發(fā)展,振動診斷在整個系統(tǒng)中的比重隨之加大[12]。同時,根據(jù)中國某保障單位近3年的數(shù)據(jù)統(tǒng)計,振動超標、性能故障、滑油消耗量大、金屬含量超標、漏油、滑油壓力低等是某型現(xiàn)役發(fā)動機地面試車常出現(xiàn)的故障現(xiàn)象[12],其中振動超標占70%以上。因此,提升振動PHM 技術(shù)是中國航空發(fā)動機健康監(jiān)測發(fā)展的重中之重。

    國外航空發(fā)動機振動健康監(jiān)測能力和水平已達到成熟的應用狀態(tài),在航空發(fā)動機振動監(jiān)測與診斷技術(shù)及應用方面,具有代表性的是瑞士Vibro-Meter 公司的UEVM 系統(tǒng)[22],羅羅(R-R)公司的QUICK 系統(tǒng)[23],美國GE 公司的HUMS系統(tǒng)[24],霍尼韋爾(HUNEYWELL)的Zing 系統(tǒng)[25]等,這些用于渦扇發(fā)動機、燃氣輪機等的狀態(tài)監(jiān)測預警與故障診斷系統(tǒng)已相對成熟。

    瑞士Vibro-Meter 公司研制的通用發(fā)動機振動監(jiān)測儀(UEVM)已取得了大量應用,可實時監(jiān)測整機振動狀況,包括振動強度、頻率等內(nèi)容,同時實現(xiàn)遠程監(jiān)測及平衡調(diào)整,其系統(tǒng)架構(gòu)如圖3[22]所示。牛津大學與羅羅(R-R)公司聯(lián)合開發(fā)了QUICK 先進機載監(jiān)測系統(tǒng),能夠?qū)崟r監(jiān)測航空發(fā)動機在飛行過程中的各項關(guān)鍵參數(shù),包括振動、溫度、壓力等,并將振動數(shù)據(jù)與瞬時性能測量值進行信息融合,針對航空發(fā)動機提供預測性維護方案。通用電氣(GE)航空部門使用自研的HUMS 振動監(jiān)測系統(tǒng),該系統(tǒng)在振動實時監(jiān)測和振動故障診斷的基礎(chǔ)上,可進一步實現(xiàn)自適應學習,能夠根據(jù)設(shè)備的運行情況和環(huán)境變化進行動態(tài)調(diào)整,提高系統(tǒng)的故障識別準確性和效率。霍尼韋爾近年來推出了代號為Zing 的健康狀態(tài)與使用監(jiān)控系統(tǒng)(HUMS),如圖4[26]所示,能夠?qū)u軸發(fā)動機實現(xiàn)振動監(jiān)測與故障診斷,并幫助維護人員進行航空發(fā)動機風險管理和安全評估。

    可以看出,國外已經(jīng)應用的系統(tǒng)實現(xiàn)了對航空發(fā)動機的振動健康監(jiān)測。而中國的整機振動監(jiān)測診斷系統(tǒng)仍處于起步階段,機載主要功能為狀態(tài)監(jiān)視、嚴重故障告警,虛警率、漏警率均較高,導致發(fā)動機故障難以實現(xiàn)及時、有效的分析診斷。同時,地面振動檢測設(shè)備也較為欠缺,主要依賴發(fā)動機地面檢驗操縱臺、綜合調(diào)節(jié)器地面檢查操縱臺、孔探儀、渦流檢測儀等,導致發(fā)動機排故能力仍有待提高。因此,近年來中國航空發(fā)動機振動故障監(jiān)測診斷技術(shù)研究受到廣泛關(guān)注,但目前主要研究均集中在理論基礎(chǔ)層面,急需進一步開展技術(shù)落地應用。

    中國現(xiàn)役三代機發(fā)動機沿用的是俄制磁電速度式振動傳感器和整機振動監(jiān)控標準,實踐發(fā)現(xiàn),原俄制振動傳感器存在頻響低、數(shù)量少兩大問題。頻響低,主要滿足轉(zhuǎn)子轉(zhuǎn)速跟蹤用,導致軸承、葉片、齒輪等高頻振動特征難以感知;數(shù)量少,現(xiàn)役發(fā)動機大多監(jiān)測一個機載振動測點,導致測試信號難以直接、有效地監(jiān)測主要旋轉(zhuǎn)部件[11]。

    2 典型故障特征及診斷方法研究現(xiàn)狀

    2.1 轉(zhuǎn)子不平衡故障

    轉(zhuǎn)子不平衡故障是由于轉(zhuǎn)子部件質(zhì)量偏心或轉(zhuǎn)子部件出現(xiàn)缺損造成的。在實際生產(chǎn)中,不平衡質(zhì)量是不可避免的,因此轉(zhuǎn)子不平衡是航空發(fā)動機乃至旋轉(zhuǎn)機械中常見的故障之一。據(jù)某發(fā)動機維修廠統(tǒng)計,每年因不平衡振動故障而造成的提前返廠量占返廠發(fā)動機總數(shù)的60%以上,嚴重影響發(fā)動機出勤率。

    航空發(fā)動機多采用雙轉(zhuǎn)子結(jié)構(gòu),如圖5所示。國內(nèi)外學者通常將剛度、阻尼等線性化處理,以便分析雙轉(zhuǎn)子系統(tǒng)動力學特性。國外學者Gupta等[27]將傳遞矩陣法和試驗研究相結(jié)合,分析了雙轉(zhuǎn)子模型的不平衡響應,盡管該雙轉(zhuǎn)子模型與實際航空發(fā)動機轉(zhuǎn)子系統(tǒng)存在差異,但仍能說明內(nèi)外轉(zhuǎn)子之間存在相互影響。Ferraris 等[28]采用動力學建模的方法,研究了雙轉(zhuǎn)子系統(tǒng)同轉(zhuǎn)以及對轉(zhuǎn)時的響應特征。國內(nèi)學者同樣對轉(zhuǎn)子不平衡故障開展了深入研究。胡絢[29]、蔣云帆[30]等研究雙轉(zhuǎn)子系統(tǒng)不平衡振動響應特性時,采用傳遞矩陣法與試驗相結(jié)合的方式探究了雙轉(zhuǎn)子系統(tǒng)同轉(zhuǎn)以及對轉(zhuǎn)時,臨界轉(zhuǎn)速與不平衡響應的差異。韓軍等[31-32]利用傳遞矩陣法建立雙轉(zhuǎn)子模型,分析了雙轉(zhuǎn)子不同測點的不平衡響應與轉(zhuǎn)速、不平衡量之間的靈敏度關(guān)系,并且通過對雙轉(zhuǎn)子同轉(zhuǎn)或?qū)D(zhuǎn)的數(shù)值仿真分析,發(fā)現(xiàn)當轉(zhuǎn)速比較接近時發(fā)動機會出現(xiàn)拍振現(xiàn)象,引起振動強度過大。廖明夫等[33-34]通過理論公式推導,分析了中介軸承剛度以及阻尼對雙轉(zhuǎn)子模型中高、低壓轉(zhuǎn)子不平衡響應的影響,同時也指出當高、低壓轉(zhuǎn)速接近時,會出現(xiàn)明顯的拍振現(xiàn)象。馬平平等[35]針對某型航空發(fā)動機結(jié)構(gòu),建立了四盤五支點雙轉(zhuǎn)子模型,研究了在不平衡激勵下,轉(zhuǎn)速、不平衡位置、不平衡量大小、不平衡相位、中介支承剛度等因素對雙轉(zhuǎn)子系統(tǒng)振動響應的影響。

    圖5 航空發(fā)動機復雜轉(zhuǎn)子系統(tǒng)Fig.5 Complex rotor system of the aero-engine

    航空發(fā)動機轉(zhuǎn)子系統(tǒng)中,不同轉(zhuǎn)子間常采用中介軸承連接,受聯(lián)軸器、軸承、擠壓油膜阻尼器等非線性部件的影響,航空發(fā)動機轉(zhuǎn)子系統(tǒng)的動力特性更加復雜。國外學者Hibner[36]基于傳遞矩陣法分析對比了有/無非線性黏性阻尼時,多軸航空發(fā)動機轉(zhuǎn)子系統(tǒng)的不平衡振動響應特點。Gunter 等[37]、Bonello 和Hai[38]考慮在擠壓油膜阻尼器的非線性影響下航空發(fā)動機轉(zhuǎn)子的不平衡響應。國內(nèi)學者同樣對此展開了研究,陳曦等[39]針對大涵道比渦扇發(fā)動機低壓轉(zhuǎn)子中套齒式剛性聯(lián)軸器結(jié)構(gòu),利用理論分析與試驗驗證,并通過改變軸向擰緊力矩來改變聯(lián)軸器的非線性連接剛度,發(fā)現(xiàn)對應的不平衡響應幅值會有不同程度地改變,且改變不平衡量的位置會導致不平衡響應幅值大小及其變化率有所不同。路振勇[40]、鄧四二[41]等考慮中介支承非線性影響,建立雙轉(zhuǎn)子系統(tǒng)動力學模型,研究了支承結(jié)構(gòu)參數(shù)對雙轉(zhuǎn)子系統(tǒng)不平衡振動響應的影響。羅貴火等[42]考慮中介軸承非線性影響,建立了航空發(fā)動機雙轉(zhuǎn)子-滾動軸承-機匣耦合動力學模型,采用數(shù)值分析與試驗驗證相結(jié)合的方式,對同轉(zhuǎn)和對轉(zhuǎn)雙轉(zhuǎn)子系統(tǒng)的不平衡響應進行了分析,指出對于同轉(zhuǎn)和對轉(zhuǎn)的雙轉(zhuǎn)子系統(tǒng),當高、低壓轉(zhuǎn)子的轉(zhuǎn)速差較小時,系統(tǒng)的拍振響應明顯,振動響應較大。陳果[43]針對典型雙轉(zhuǎn)子航空發(fā)動機結(jié)構(gòu),建立了雙轉(zhuǎn)子-滾動軸承-機匣耦合動力學模型,分析了擠壓油膜阻尼器的減振性能以及系統(tǒng)對不平衡故障的靈敏度。楊喜關(guān)等[44]綜合考慮了軸承與擠壓油膜阻尼器支承系統(tǒng)的非線性對雙轉(zhuǎn)子系統(tǒng)不平衡振動的影響。

    國內(nèi)外有關(guān)航空發(fā)動機轉(zhuǎn)子不平衡故障的研究成果十分豐富。針對航空發(fā)動機轉(zhuǎn)子不平衡的研究多集中在航空發(fā)動機復雜模型的建立、部件結(jié)構(gòu)參數(shù)的影響、不同位置的響應特點、轉(zhuǎn)子同轉(zhuǎn)及對轉(zhuǎn)的影響等方面。多利用工頻幅值變化規(guī)律并結(jié)合航空發(fā)動機不平衡故障振動特性,綜合診斷轉(zhuǎn)子不平衡故障。目前,針對航空發(fā)動機轉(zhuǎn)子不平衡故障問題,仍需要對以下幾個方面進行進一步研究。

    1) 考慮聯(lián)軸器、軸承、擠壓油膜阻尼器等多部件的耦合,建立起更加復雜全面的航空發(fā)動機模型。當前建立的模型對航空發(fā)動機各部件結(jié)構(gòu)進行了不同程度的簡化,且多數(shù)文獻只考慮某一部件結(jié)構(gòu)的影響,忽略了其他部件結(jié)構(gòu)的影響,仍需進一步研究。

    2) 結(jié)合航空發(fā)動機實測數(shù)據(jù),研究航空發(fā)動機不同位置發(fā)生不平衡,以及各部件結(jié)構(gòu)參數(shù)變化時的不平衡響應特征。目前國內(nèi)外已有大量文獻研究了轉(zhuǎn)子不同位置發(fā)生不平衡故障以及部件結(jié)構(gòu)參數(shù)變化時的響應特征,且多通過搭建轉(zhuǎn)子實驗臺進行實驗驗證,但轉(zhuǎn)子實驗臺實驗數(shù)據(jù)與航空發(fā)動機實測數(shù)據(jù)仍有很大差距。因此結(jié)合航空發(fā)動機實測數(shù)據(jù)進行分析及驗證的問題仍需進一步深入開展。

    3) 研究航空發(fā)動機轉(zhuǎn)子不平衡的故障定位以及定量診斷技術(shù)。航空發(fā)動機的復雜結(jié)構(gòu)、傳遞路徑、多部件耦合等多種因素,以及工作條件惡劣、傳感器安裝位置受限等問題,使得航空發(fā)動機轉(zhuǎn)子不平衡故障問題變得更加復雜。準確識別航空發(fā)動機轉(zhuǎn)子不平衡故障位置以及故障程度,對航空發(fā)動機的安全穩(wěn)定運行具有重大意義。

    2.2 支承不同心故障

    由于現(xiàn)代航空發(fā)動機對高效率、高性能等需求的提升,其轉(zhuǎn)子系統(tǒng)多采用雙轉(zhuǎn)子結(jié)構(gòu)設(shè)計。其中,低壓轉(zhuǎn)子系統(tǒng)作為大跨度柔性轉(zhuǎn)子系統(tǒng),為防止撓曲變形過大,常采用多支承的結(jié)構(gòu)設(shè)計。這種結(jié)構(gòu)設(shè)計雖然能夠滿足航空發(fā)動機的設(shè)計需求,但是受制造、裝配和長時間工作的影響,轉(zhuǎn)子支承結(jié)構(gòu)的中心往往存在一定程度的相對偏移,各支承的中心很難保持在同一條軸線上,從而出現(xiàn)“支承不同心”問題,如圖6所示。當轉(zhuǎn)子系統(tǒng)發(fā)生支承不同心故障后,將直接影響發(fā)動機的正常運行,輕則造成轉(zhuǎn)子振動加劇,軸承磨損,嚴重時甚至會發(fā)生動靜碰摩故障。

    圖6 支承不同心故障示意圖Fig.6 Schematic diagram of misalignment fault

    針對此類故障,目前主要研究集中在以下2大方面:一是通過動力學建模,研究支承不同心故障下轉(zhuǎn)子系統(tǒng)的振動響應特征,為航空發(fā)動機支承不同心故障的診斷提供理論依據(jù)。這一方面研究的關(guān)鍵在于動力學建模方法。常見的航空發(fā)動機動力學建模方法有質(zhì)心集中質(zhì)量法、傳遞矩陣法、有限元法等。目前針對航空發(fā)動機支承不同心故障的動力學建模主要集中于質(zhì)心集中質(zhì)量法和有限元法,傳遞矩陣法鮮有提及。在利用質(zhì)心集中質(zhì)量法建模方面,國內(nèi)于2012年,由張振波和Li 等[45-46]利用Lagrange 能量法建立了適用于航空發(fā)動機轉(zhuǎn)子支承不同心的動力學模型,該模型以增加非線性剛度項的方式引入支承不同心故障。通過研究模型中不同心量和不平衡量對振動響應的影響,獲得了轉(zhuǎn)子支承不同心時的頻率響應特征。在此基礎(chǔ)上,張振波和Ren 等[47-48]建立了不同連接結(jié)構(gòu)的運動微分方程,得出了不同連接結(jié)構(gòu)出現(xiàn)支承不同心后,不同心激振力各不相同的結(jié)論。進一步地,劉永泉等[49]將當量剛度引入動力學建模,實現(xiàn)對航空發(fā)動機轉(zhuǎn)子支承不同心時附加非線性剛度的定量描述。

    在利用有限元法建模方面,國內(nèi)的吳英祥等[50]則利用支承軸承間隙模型描述轉(zhuǎn)子支承不同心問題,研究發(fā)現(xiàn)支承不同心需要軸承的間隙補償,同時過大的間隙會使轉(zhuǎn)子運動呈現(xiàn)混沌現(xiàn)象。馮國全等[51]將支承不同心看作軸承不對中問題,并將軸承不對中等效為不對中彎矩,基于此建立了轉(zhuǎn)子支承軸承不對中有限元模型,利用模型分析發(fā)現(xiàn),振動響應中轉(zhuǎn)子轉(zhuǎn)頻的2 倍頻成分占據(jù)主導。王美令、張宏獻和肖森等[52-54]將支承不同心引起的附加載荷引入轉(zhuǎn)子系統(tǒng)有限元模型中,得出了出現(xiàn)支承不同心故障后振動頻率更加復雜以及振動響應特征的變化規(guī)律。柏長青等[55]基于多體動力學建立了轉(zhuǎn)子-軸承-支承系統(tǒng)的有限元模型,利用附加驅(qū)動轉(zhuǎn)速的方式模擬支承不同心,分析發(fā)現(xiàn)軸向振動工頻和倍頻成分是轉(zhuǎn)子支承不同心的典型特征。除此之外,Han 等[56-57]基于彈性鼠籠支承的剛性轉(zhuǎn)子系統(tǒng)建立了有限元模型,對支承不同心時模型的應力、應變能分布進行了分析,發(fā)現(xiàn)不同心故障后鼠籠的應力水平顯著增加,甚至具有損壞的風險。

    以上是國內(nèi)外研究學者針對航空發(fā)動機支承不同心故障的動力學建模方面的研究。通過比較國內(nèi)外研究學者的研究現(xiàn)狀,發(fā)現(xiàn)相較于國外學者,國內(nèi)針對航空發(fā)動機支承不同心故障的動力學建模研究更為深入、全面;而國外研究學者普遍將支承不同心故障視為轉(zhuǎn)子不對中故障的一部分,因此鮮少有國外學者將其單獨考慮。

    目前航空發(fā)動機支承不同心故障的振動特性可以總結(jié)為:

    1) 振動響應中轉(zhuǎn)子轉(zhuǎn)頻的2 倍頻響應占據(jù)主導地位。

    2) 軸向和徑向2 倍頻幅值相對于工頻幅值的大小可以表征支承不同心故障的嚴重程度。

    3) 軸心軌跡包含了軸承錯位的位置和方向信息,可用于支承不同心故障的診斷。

    二是通過開展航空發(fā)動機轉(zhuǎn)子系統(tǒng)支承不同心故障模擬實驗,進行航空發(fā)動機支承不同心故障的診斷方法研究。在該方向,國內(nèi)的王美令等[52,58]設(shè)計了具有不同程度支承不同心故障的實驗裝置,通過調(diào)整墊片的厚度實現(xiàn)故障程度的模擬,實驗分析發(fā)現(xiàn)支承不同心引起了明顯的軸向振動,頻譜圖中以工頻及其倍頻成分為主導。馬平平[59]基于低壓渦輪后支點不對中現(xiàn)象設(shè)計了相關(guān)模擬實驗,通過實驗分析發(fā)現(xiàn)故障后頻譜圖中出現(xiàn)顯著的轉(zhuǎn)頻2 倍頻成分,同時軸心軌跡為水滴形。

    目前航空發(fā)動機支承不同心故障的診斷方法可以總結(jié)為:采用頻譜分析轉(zhuǎn)子轉(zhuǎn)頻及其倍頻成分(尤其是2 倍頻成分),同時結(jié)合軸心軌跡綜合診斷支承不同心故障。

    現(xiàn)有研究工作中對具有支承不同心故障的航空發(fā)動機整體建模研究較少,對航空發(fā)動機支承不同心故障的診斷方法考慮較少,使得目前故障診斷方法無法滿足航空發(fā)動機支承不同心故障的診斷需求,未來的航空發(fā)動機支承不同心故障診斷的主要研究方向是:

    2) 在故障診斷方法方面,現(xiàn)有支承不同心故障特征與轉(zhuǎn)子不對中故障特征相似,未來需要基于支承不同心故障理論研究建立針對航空發(fā)動機支承不同心故障的診斷方法。

    3) 航空發(fā)動機支承不同心故障的定量診斷。對于航空發(fā)動機來說,支承不同心故障的定量分析是一個難題,現(xiàn)有相關(guān)研究較少。故障的定量診斷對提高故障診斷準確率和保障設(shè)備安全運行具有重要意義,未來仍需要對航空發(fā)動機支承不同心故障的定量診斷進行深入研究。

    2.3 轉(zhuǎn)子不對中故障

    轉(zhuǎn)子不對中同樣是航空發(fā)動機乃至旋轉(zhuǎn)機械中典型故障之一,通常是指由于聯(lián)軸器結(jié)構(gòu)的制造誤差、裝配不當、熱變形等因素,導致應在同一軸線的相鄰轉(zhuǎn)子存在一定的偏移或傾斜量的現(xiàn)象。因此聯(lián)軸器不對中通常包括平行不對中、角度不對中以及綜合不對中,如圖7[60]所示。轉(zhuǎn)子不對中故障會引起轉(zhuǎn)子系統(tǒng)的軸向振動和徑向振動,嚴重時會影響整機穩(wěn)定運行,導致發(fā)生嚴重事故。國外的統(tǒng)計表明,航空發(fā)動機整機振動與不對中相關(guān)的故障占到系統(tǒng)故障的50%~60%;據(jù)美國海軍飛機維護庫數(shù)據(jù)調(diào)查顯示,40%的固定翼飛機和70%旋轉(zhuǎn)翼飛機存在的附件花鍵連接損傷問題主要是與不對中故障有關(guān)[60-62]。

    圖7 不對中故障形式[60]Fig.7 Form of Misalignment fault[60]

    國內(nèi)外學者一般通過3 種方法獲得聯(lián)軸器不對中的動力學模型[60,63]。第1 種是根據(jù)聯(lián)軸器結(jié)構(gòu)特點,基于聯(lián)軸器的變形幾何關(guān)系和受力分析建模;第2 種是對于自由度較少的轉(zhuǎn)子系統(tǒng),采用系統(tǒng)整體拉格朗日能量法建模;第3 種是采用有限元或集中參數(shù)法,將聯(lián)軸器看作轉(zhuǎn)子系統(tǒng)的一部分,用等效軸段建模。

    國內(nèi)外學者對航空發(fā)動機轉(zhuǎn)子不對中故障的研究取得了一定成果。國外學者Leen 等[64]通過模擬試飛條件并結(jié)合三維有限元模型分析,研究了航空發(fā)動機聯(lián)軸器的疲勞與微動磨損。航空發(fā)動機轉(zhuǎn)子不對中問題在中國也受到了越來越多的關(guān)注。國內(nèi)學者李全坤等[65]根據(jù)幾何關(guān)系與受力分析,建立了考慮不對中故障的雙轉(zhuǎn)子動力學模型,通過解析求解、數(shù)值分析和試驗驗證,研究了雙轉(zhuǎn)子系統(tǒng)存在不對中故障時,高、低壓轉(zhuǎn)子的振動響應特性。結(jié)果表明低壓轉(zhuǎn)子發(fā)生不對中故障時,由于高、低壓轉(zhuǎn)子間的耦合作用,低壓轉(zhuǎn)子不對中故障會影響高壓轉(zhuǎn)子的振動響應。陳果和李興陽[66]主要以航空發(fā)動機高壓轉(zhuǎn)子及后機匣為研究對象,考慮套齒聯(lián)軸器綜合不對中,并根據(jù)幾何關(guān)系與受力分析,建立了轉(zhuǎn)子-滾動軸承-機匣耦合動力學模型,通過數(shù)值仿真分析,發(fā)現(xiàn)不對中故障將引起系統(tǒng)的超諧共振,其中基頻在臨界轉(zhuǎn)速下會引起共振,2 倍頻會在1/2 臨界轉(zhuǎn)速下引起共振。李俊慧等[67]利用有限元法對航空發(fā)動機常用的套齒聯(lián)軸器結(jié)構(gòu)進行了動力學分析,并提出了套齒結(jié)構(gòu)動力學設(shè)計方法。蔣玲莉等[68-69]以某型航空發(fā)動機雙轉(zhuǎn)子系統(tǒng)為研究對象,建立了航空發(fā)動機雙轉(zhuǎn)子系統(tǒng)有限元模型,對比了正常狀態(tài)與平行不對中故障對固有頻率的影響,并且考慮了聯(lián)軸器的影響,將聯(lián)軸器?;癁橐粋€當量軸段,之后通過數(shù)值仿真分析,研究了存在不平衡激勵與聯(lián)軸器發(fā)生角度不對中情況下不同支點處的振動響應與軸心軌跡規(guī)律,結(jié)果表明低壓轉(zhuǎn)子聯(lián)軸器不對中振動特性可通過中介軸承傳遞到高壓轉(zhuǎn)子,進而影響高壓轉(zhuǎn)子的振動。

    國內(nèi)外有關(guān)轉(zhuǎn)子不對中的成果十分豐富,但關(guān)于航空發(fā)動機轉(zhuǎn)子不對中的相關(guān)文獻仍舊較少。針對航空發(fā)動機轉(zhuǎn)子不對中的研究多集中在聯(lián)軸器的動力學建模、聯(lián)軸器不對中與轉(zhuǎn)子系統(tǒng)的耦合作用、轉(zhuǎn)子不對中傳遞規(guī)律等方面。多采用頻譜分析2 倍頻幅值變化,并結(jié)合航空發(fā)動機轉(zhuǎn)子不對中故障振動特性,綜合診斷轉(zhuǎn)子不對中故障。針對航空發(fā)動機轉(zhuǎn)子不對中問題,以下幾方面仍需進一步研究。

    1) 建立航空發(fā)動機聯(lián)軸器動力學模型。國內(nèi)外關(guān)于聯(lián)軸器動力學建模已取得大量成果,但航空發(fā)動機的復雜結(jié)構(gòu)以及工作環(huán)境等因素,使得航空發(fā)動機聯(lián)軸器動力學建模問題更加復雜,仍需進一步研究。

    2) 結(jié)合航空發(fā)動機實測數(shù)據(jù)以及動力學模型,研究航空發(fā)動機轉(zhuǎn)子不對中傳遞規(guī)律。中介軸承耦合作用以及其他復雜結(jié)構(gòu)使得航空發(fā)動機的振動信號傳遞十分復雜,使得診斷分析變得困難,仍需要進一步研究。

    黃河作為世界第五長河,世界含沙量最大的黃色河流,孕育了五千年華夏文明,而現(xiàn)今黃河的旅游開發(fā)狀況與其國際地位嚴重不符。山西黃河全域旅游建設(shè)應放眼國際,改變單一、保守的國內(nèi)客源市場開拓為主的觀念,利用國內(nèi)外強勢媒體進行宣傳推介,打造“黃河旅游”這一國家名片。

    2.4 動靜碰摩故障

    碰摩故障是旋轉(zhuǎn)機械中常見的故障類型,航空發(fā)動機作為一種高精度的氣動熱力旋轉(zhuǎn)機械,出現(xiàn)動靜碰摩故障的概率極大[70]。目前,在航空發(fā)動機設(shè)計時,為了增加喘振裕度和提高壓氣機的效率,普遍將航空發(fā)動機設(shè)計成雙轉(zhuǎn)子-軸承-機匣機構(gòu)的復雜耦合系統(tǒng),這使得發(fā)動機的動力學特性變得復雜多變。同時,為了提高發(fā)動機的推重比和結(jié)構(gòu)效率,轉(zhuǎn)子與靜子之間間隙設(shè)計的越來越小。這也導致了發(fā)動機碰摩故障發(fā)生的風險急劇增加[71]。葉片與機匣碰摩示意圖,如圖8所示。當發(fā)動機發(fā)生碰摩故障時,會引起整機劇烈振動,嚴重時甚至會造成動靜葉片斷裂,轉(zhuǎn)軸永久性彎曲等事故,從而嚴重影響生命以及財產(chǎn)安全。因此,針對航空發(fā)動機碰摩故障振動特性以及故障預測的研究有著極其重要的意義。

    圖8 葉片與機匣的碰摩示意圖Fig.8 Schematic diagram of rub-impact between blades and casing

    針對航空發(fā)動機碰摩故障,國外學者也在進行不斷的研究。Prabith 和Praveen krishna[72]為代表的國外學者針對碰摩故障建立了多種碰摩模型,通過模型仿真研究了碰摩故障的復雜振動特性。針對不同碰摩類型、系統(tǒng)參數(shù)變化以及碰摩引發(fā)的不同耦合故障機理進行著不斷仿真和實驗分析。

    在航空發(fā)動機碰摩故障類型模擬以及各參數(shù)變化對碰摩故障的影響方面,國內(nèi)學者做出了巨大貢獻。林學森等[73]建立了帶柔性靜子部件的軸承共腔結(jié)構(gòu)的渦軸發(fā)動機的碰摩動力學模型。結(jié)果表明該結(jié)構(gòu)與傳統(tǒng)雙轉(zhuǎn)子結(jié)構(gòu)中中介軸承具有相似的動力學特征。袁惠群等[74]建立了雙轉(zhuǎn)子-機匣簡化耦合系統(tǒng)的動力學模型,得出當碰摩剛度增加時,共振區(qū)范圍以及振動幅值明顯增加。張?zhí)斐痰龋?5]建立了航空發(fā)動機雙轉(zhuǎn)子系統(tǒng)含碰摩力的彎扭耦合動力學模型,得出扭轉(zhuǎn)振動信號頻譜特性更適用于碰摩故障診斷的結(jié)論。羅貴火等[76]提出了一種解決含局部非線性雙轉(zhuǎn)子動力特性的高效方法,并得出碰摩故障會嚴重影響內(nèi)轉(zhuǎn)子系統(tǒng)同步響應頻率幅值。聶日敏等[77]通過模型得出高壓渦輪碰摩會導致系統(tǒng)發(fā)生反向渦動、低壓渦輪碰摩會抑制高壓碰摩以及高低壓渦輪同時碰摩相比于高壓渦輪碰摩呈現(xiàn)出轉(zhuǎn)速區(qū)間更窄、振動幅值更低的特點。丁小飛等[78]進行了高壓渦輪葉片和機匣的碰摩仿真分析,仿真結(jié)果表明碰摩故障導致發(fā)動機振動幅值明顯增大。秦海勤等[79]得出可由轉(zhuǎn)子和機匣頻譜圖判斷是否發(fā)生局部碰摩,并提出高壓轉(zhuǎn)子的碰摩幅值隨偏心量和碰摩剛度的增加愈加明顯的結(jié)論。陳松霆和吳志強[80]認為碰摩剛度和支承彈簧剛度的增大有利于保持系統(tǒng)的穩(wěn)定性。張婭等[81]通過模型分析得出僅由不平衡量導致的軸向碰摩系統(tǒng)具有較好的穩(wěn)定性的結(jié)論。

    航空發(fā)動機碰摩故障信號的處理方法也在不斷的更新改進,以提高故障診斷的準確性。在碰摩信號處理方法方面。何田等[82]提出了基于奇異值分解(SVD)的故障特征提取方法以適用于碰摩特征提取,能夠在強干擾背景下診斷,并能夠精確的檢測不同位置碰摩發(fā)生先后。劉洋等[83]認為基于傅里葉分解方法(FDM)的轉(zhuǎn)子碰摩故障診斷方法,表現(xiàn)了出更高的可靠性和有效性。邊杰等[84]提出一種參數(shù)自適應的變分模態(tài)分解方法(APVMD)對雙轉(zhuǎn)子航空發(fā)動機振動狀態(tài)實測彈性支承應變信號進行分析。陳果等[85]提出了一種基于倒譜分析的碰摩部位識別方法,可準確識別碰摩部位樣本特征。左紅艷等[86]提出的小波變換和Hilbert-Huang 變換方法在密度分布不均勻的多類型混合數(shù)據(jù)表現(xiàn)為更高的精度。張志禹等[87]采用廣義S 變換方法檢測和提取故障特征,該方法相對于小波變換的結(jié)果具有更精確的碰摩位置檢測和良好的抗干擾能力。周海侖等[88]利用支持向量機從大量樣本中獲取碰摩知識并進行判斷,結(jié)果表明該方法針對碰摩故障具有很高的識別率。

    現(xiàn)階段,國內(nèi)外學者都通過碰摩故障模型的搭建與改進,來不斷完善理論碰摩故障模型。通過模型來不斷分析不同碰摩類型以及不同參數(shù)下的碰摩故障特征。并且在碰摩故障振動信號的處理方面也在不斷的更新,傅里葉變換、小波變換等處理方法的應用使得碰摩故障的診斷更加快速準確。但如今,國內(nèi)外大多數(shù)學者關(guān)于碰摩故障的研究主要集中于對所建模型的數(shù)值仿真模擬進行計算,而缺乏真實的實際工程驗證。在航空發(fā)動機運轉(zhuǎn)過程中,復雜的運行機理可能會導致碰摩故障的振動特征難以凸顯,難以驗證所建模型及仿真分析結(jié)果的有效性。同時還應該進一步建立更靠近真實工況的碰摩模型進行碰摩分析并發(fā)展更為完善的碰摩損傷程度評估方法。

    2.5 葉片斷裂故障

    航空發(fā)動機是一種高速旋轉(zhuǎn)的動力機械,葉片作為航空發(fā)動機的關(guān)鍵部件,因為其工作環(huán)境惡劣,如高溫、高壓和高速,導致其成為故障頻發(fā)的部件[89]。據(jù)統(tǒng)計,葉片故障是航空發(fā)動機中常見的重大問題,約占整機故障的42%。葉片斷裂會導致轉(zhuǎn)子的平衡性被破壞,甚至可能導致災難性事故和發(fā)動機的徹底毀壞[90],如圖9所示。目前國內(nèi)外學者主要從以下幾個方面對葉片斷裂故障進行了研究。

    圖9 葉片斷裂故障示意圖Fig.9 Schematic diagram of blade fracture fault

    2.5.1 基于動力學模型的葉片斷裂故障

    通過建立動力學模型來分析葉片斷裂引起的系統(tǒng)響應特性的變化,分析模型中不同部件的耦合特性,發(fā)展了基于轉(zhuǎn)軸振動、葉片振動和機匣振動跟蹤分析的葉片故障診斷方法。

    2013年,Sinha[91]考慮風扇葉片斷裂后轉(zhuǎn)子的非對稱特征。推導了具有可變形葉片的非對稱轉(zhuǎn)子的轉(zhuǎn)子動力學方程。結(jié)果表明轉(zhuǎn)盤中心和機匣響應中呈現(xiàn)出了非線性現(xiàn)象。Hong 等[92]提出了葉片斷裂后轉(zhuǎn)子載荷傳遞系統(tǒng)的瞬時微分方程。得出葉片斷裂后突加的沖擊載荷在靜子傳遞過程中會嚴重衰減,阻尼會降低沖擊響應的極值。Hong 等[93-94]、Xie 等[95]研究了葉片丟失下航空發(fā)動機整機動力學響應,得到了葉片丟失使轉(zhuǎn)子轉(zhuǎn)速降低,并在復雜載荷激勵下呈現(xiàn)出非線性特性的結(jié)論。當葉片丟失導致碰摩發(fā)生時,在轉(zhuǎn)子系統(tǒng)響應中伴隨有轉(zhuǎn)子橫向固有振動,持續(xù)碰摩還會改變系統(tǒng)的臨界轉(zhuǎn)速,同時提出了平均頻率和帶寬指標作為葉片損傷監(jiān)測指標的一種方法。Wang 等[96-97]建立了葉片斷裂下與機殼發(fā)生碰摩故障的雙轉(zhuǎn)子-葉片-機匣系統(tǒng)有限元模型。葉片斷裂會在機匣上產(chǎn)生沖擊頻率為葉片通過頻率的振動,并且在沖擊頻率附近激勵起雙轉(zhuǎn)子轉(zhuǎn)速差的倍頻成分;由于陀螺力矩的影響,也會在轉(zhuǎn)軸位移信號中激勵出前2 階轉(zhuǎn)子進動頻率。Wang 等[98]基于有限元方法和自由度減縮,建立了航空發(fā)動機含葉片脫落的轉(zhuǎn)子-葉片系統(tǒng)模型,通過采用分段擬合的方法對轉(zhuǎn)子和葉片結(jié)構(gòu)進行建模,通過試驗臺驗證了模型的有效性,結(jié)果表明,葉片脫落的瞬間會產(chǎn)生較大的瞬態(tài)沖擊力以及較大的系統(tǒng)阻尼,同時葉片會激起轉(zhuǎn)子的低階固有振型。

    基于葉片動力學模型的葉片故障診斷方法是通過建立動力學方程,分析葉片斷裂引起的響應變化,從而實現(xiàn)對葉片故障的診斷。在基于動力學模型分析葉片固有振動特性變化的研究中,尚需提高模型預測的準確性、增強葉片故障特征的解釋性以及解決方法在工程實際應用中的可行性與適用性問題。

    2.5.2 基于信號處理的葉片斷裂故障

    基于葉片斷裂故障的產(chǎn)生機理,分析葉片斷裂前后葉片振動信號、位移信號、轉(zhuǎn)軸振動信號、轉(zhuǎn)速信號、聲壓信號、溫度信號和機匣振動信號中的故障敏感特征,基于先進的信號處理算法形成不同的葉片斷裂故障診斷方法。

    2013年,F(xiàn)orbes 和Randall[99]首次根據(jù)航空發(fā)動機機匣振動信號估計葉片固有頻率,并推導了旋轉(zhuǎn)葉片的壓力信號以及機匣振動響應模型(SDOF)。得出在機匣振動信號中包含以某個轉(zhuǎn)頻倍頻±葉片固有頻率間隔分布的窄帶峰值。Liska 等[100]采用FFT、濾波、加權(quán)平均等算法提取軸振信號中倒頻譜的葉片特征分量,可實時自動跟蹤葉片的固有頻率及其振幅,并在250 MV汽輪機中進行了商業(yè)驗證。胡明輝等[101-103]研究了基于機匣振動的航空發(fā)動機壓氣機動葉斷裂故障機理,建立了葉片斷裂下機匣振動響應的數(shù)學模型。分析了振動信號中葉片通過頻率幅值、工頻幅值、工頻相位等斷裂故障敏感特征,采用先進信號處理算法建立了基于機匣振動信號的葉片斷裂故障診斷方法,并在實際燃氣輪機葉片斷裂故障案例中得到有效應用。Feng 等[104]提出了一種基于稀疏諧波積譜(SHPS)的燃氣輪機葉片狀態(tài)監(jiān)測模型。該模型可實現(xiàn)葉片相關(guān)振動分量與機匣振動分量的準確分離,比傳統(tǒng)的閾值法更早、更準確地挖掘潛在地葉片斷裂等故障類型。Shan 等[105]以時域跟蹤法和時頻轉(zhuǎn)換法計算了葉片的非線性響應,結(jié)果表明凸肩結(jié)構(gòu)對于航空發(fā)動機風扇葉片的一階彎曲振動有很好的減振效果,并且可以降低葉片發(fā)生斷裂故障時的應力。Gubran 和Sinha[106]通過軸的瞬時角速度(IAS)信號來監(jiān)測葉片的狀態(tài)。葉片故障會導致在葉片固有頻率附近出現(xiàn)2 個或多個峰值及高次諧波,可為葉片的健康監(jiān)測提供基礎(chǔ)。Zhang等[107]提出了一種基于自適應噪聲完全集成經(jīng)驗模態(tài)分解(CEEMADAN)和變步長歸一化最小均方(VSS-NLMS)算法的聲壓信號降噪算法,對含強噪聲干擾下不同轉(zhuǎn)速、不同狀態(tài)(正常葉片、斷裂葉片)產(chǎn)生的聲壓信號具有良好的降噪效果,并且在斷裂葉片頻譜中,葉片通過頻率及諧波頻率被完整保留,可用于葉片斷裂故障診斷。Abdelrhma 等[108]提出了2 種具有高頻分辨率和高時間分辨率的小波算法,用于識別多級轉(zhuǎn)子系統(tǒng)中微弱的葉片故障特征。

    基于信號處理算法的葉片故障診斷主要通過對與葉片相關(guān)的多種振動信號進行分析和處理,提取出故障敏感特征,并基于這些特征形成葉片斷裂故障的診斷方法。主要的故障特征包括葉片固有頻率、振動幅值、相位、工頻幅值、工頻相位以及諧波頻率等。然而,基于信號處理的葉片斷裂故障診斷方法仍然面臨一些待解決的問題,主要體現(xiàn)在算法優(yōu)化、特征解釋性和工程實際應用等方面。

    2.5.3 基于葉間定時技術(shù)的葉片故障

    轉(zhuǎn)子葉片振動測量技術(shù)分為接觸式測量和非接觸式測量2 種。其中,非接觸式測量相較而言具有很大的優(yōu)勢。葉尖定時(BTT)是一種新興的葉片振動監(jiān)測技術(shù),用于檢測、測量和分析旋轉(zhuǎn)葉片組件中的葉片振動,通過分析BTT 信號中與葉片故障相關(guān)的特征信息,可實現(xiàn)葉片故障診斷。

    國外對葉尖定時方法的研究始于20 世紀60~70年代,德國MTU(Motoren-und Turbine-Union)研究人員(2005)基于自主研發(fā)的BSSM 葉尖定時系統(tǒng),MTU 的研究人員通過監(jiān)測葉片二彎模態(tài)振動的最大應力,在共振狀態(tài)下識別了高階模態(tài)振動導致的高周裂紋[109]。NASA Glenn 研究中心和克利夫蘭州立大學合作,開展了葉尖定時系統(tǒng)在轉(zhuǎn)子輪盤裂紋故障檢測中的探索[110-111],GE 公司則將葉尖定時系統(tǒng)和其他測量方法一起應用于TF41-A1B 發(fā)動機的振動故障檢測[112-113]。韓國學者針對BTT 信號欠采樣的問題,定義了一個稱為馬氏距離的統(tǒng)計指標,并將其用于具有裂紋的旋轉(zhuǎn)多葉片系統(tǒng)的故障診斷。并利用模擬模型獲得的BTT 信號,研究了裂紋存在和信噪比對所提出方法可靠性的影響[114]。

    中國學者對葉尖定時技術(shù)的研究起步較晚,天津大學段發(fā)階團隊通過建立BTT 仿真模型,設(shè)計并搭建BTT 測量系統(tǒng),實現(xiàn)葉片的狀態(tài)監(jiān)測[115]。中國航發(fā)沈陽發(fā)動機研究所和湖南動力機械研究所引進了HOOD 公司的BTT 測量設(shè)備,并在相關(guān)的試車臺上進行了測試實驗,以獲取振動數(shù)據(jù)用于葉片的強度分析[116]。此外,西安交通大學、北京化工大學等對BTT 技術(shù)也有較多文獻報道,主要研究重點集中在高速數(shù)據(jù)采集、振動參數(shù)提取、優(yōu)化探頭布局及測振系統(tǒng)集成開發(fā)等方面。利用BTT 信號實現(xiàn)葉片健康檢測的主要難點之一在于準確獲取葉片的振動參數(shù),由于BTT 信號具有嚴重的欠采樣特性,近幾年學者從不同的角度開發(fā)了許多分析方法來研究如何從BTT 采樣信息中有效地表征葉片的真實振動。Chen 及其團隊提出了一種有效且高效的基于時延的頻譜重構(gòu)方法,適用于大多數(shù)BTT方法的參數(shù)識別或頻譜重構(gòu),為BTT 測量提供了指導[117]。

    綜上所述,基于BTT 技術(shù)的葉片故障診斷依賴于對BTT 信號特征信息的分析,以識別葉片振動參數(shù)和共振頻率等特征,從而實現(xiàn)葉片健康狀態(tài)的檢測和故障診斷。但是,應用BTT 技術(shù)實現(xiàn)葉片故障診斷仍需實現(xiàn)數(shù)據(jù)處理算法的優(yōu)化、BTT 系統(tǒng)的穩(wěn)定性與可靠性的提高、工程實際應用以及故障診斷準確性的提高等。

    2.5.4 基于人工智能的葉片斷裂故障

    基于人工智能的分類系統(tǒng)可有效識別航空發(fā)動機葉片故障,最大限度地減少人工干預。人工智能方法能從大量數(shù)據(jù)中提取故障敏感特征,并將其作為故障診斷分類器的輸入,對葉片故障進行分類診斷。

    Kuo 等[118]利用非線性系統(tǒng)的特性將完整渦輪振動信號加載到統(tǒng)一的混沌系統(tǒng)中,開發(fā)了一種渦輪葉片的故障診斷儀,根據(jù)可拓學理論對信號進行人工智能故障診斷。Li 等[119]提出一種結(jié)合歷史監(jiān)測數(shù)據(jù)的相似性傳播聚類的故障診斷方法,并在航空發(fā)動機轉(zhuǎn)子故障試驗臺上模擬葉片斷裂故障并驗證了方法的有效性。Shapiai 等[120-121]提出了基于小波變換的人工神經(jīng)網(wǎng)絡(luò)模型來識別航空發(fā)動機振動信號中的葉片故障特征,根據(jù)小波系數(shù)統(tǒng)計參數(shù)特征作為(ANN)的輸入,結(jié)果表明,該模型對葉片故障分類的準確率為88.43%。Wang 等[122]將多通道卷積神經(jīng)網(wǎng)絡(luò)(MCNN)算法應用于風力渦輪機葉片角度異常、葉片表面損傷、葉片斷裂等故障診斷中,開發(fā)了一種葉片狀態(tài)監(jiān)測的故障診斷系統(tǒng)。Tan 等[123]提出了一種基于本征模態(tài)函數(shù)的經(jīng)驗模態(tài)分解 和平均統(tǒng)計參數(shù)特征提取的葉片故障診斷方法。

    綜合來看,針對航空發(fā)動機中的葉片斷裂故障,未來研究的重點將集中在深入挖掘?qū)崪y數(shù)據(jù),進一步完善葉片動力學模型、發(fā)展葉片的振動監(jiān)測技術(shù)、結(jié)合信號處理算法和人工智能技術(shù),實現(xiàn)對葉片斷裂故障的精準診斷和預測,以提高葉片的安全性和可靠性。

    2.6 結(jié)構(gòu)共振類故障

    振動問題一直是航空發(fā)動機設(shè)計、制造和使用中無法避開的議題,但引起航空發(fā)動機振動超標的原因很復雜[124]。其中結(jié)構(gòu)共振導致的振動超標雖不經(jīng)常發(fā)生,但一有發(fā)生便易引起嚴重故障,因此針對航空發(fā)動機結(jié)構(gòu)共振類故障的研究十分必要。航空發(fā)動機作為振動系統(tǒng),擁有多子系統(tǒng)、多階固有頻率。當其出現(xiàn)某些振動故障后將在系統(tǒng)內(nèi)引入激勵,此時如果高、低壓轉(zhuǎn)頻等激振頻率接近某一階固有頻率,系統(tǒng)會出現(xiàn)結(jié)構(gòu)共振,振幅急劇增大[125]。這種由某些振動故障引起的故障共振稱為結(jié)構(gòu)共振類故障。

    針對此類故障,目前主要研究集中在以下2大方面:一是通過動力學建模,研究結(jié)構(gòu)共振故障下轉(zhuǎn)子系統(tǒng)的振動響應特征,為航空發(fā)動機結(jié)構(gòu)共振類故障的診斷提供理論依據(jù)。這一方面研究的關(guān)鍵在于故障轉(zhuǎn)子系統(tǒng)的建立。常見的故障轉(zhuǎn)子系統(tǒng)有裂紋故障轉(zhuǎn)子系統(tǒng)、碰摩故障轉(zhuǎn)子系統(tǒng)等。首先,對于裂紋故障轉(zhuǎn)子系統(tǒng),國外的Spagnol 等[126]將呼吸裂紋模型中重力主導的假設(shè)去除,并引入不平衡研究兩者的耦合效應,研究結(jié)果發(fā)現(xiàn)較大的不平衡量和180°的偏角會導致超諧共振頻率分量的消失。國內(nèi)的Hou[127-130]在呼吸裂紋模型的基礎(chǔ)上,引入了機動載荷[127-128]、慣性激勵[129]和不平衡激勵[130],通過分析系統(tǒng)的非線性響應,發(fā)現(xiàn)轉(zhuǎn)子系統(tǒng)超諧共振的產(chǎn)生受到呼吸裂紋、慣性激勵、機動載荷和不平衡激勵的顯著影響。牛和強[131]建立了含斜裂紋、直裂紋和半拋物線斜裂紋的轉(zhuǎn)子系統(tǒng)有限元模型,研究發(fā)現(xiàn)半拋物線斜裂紋和斜裂紋在扭轉(zhuǎn)方向會有明顯的共振響應而直裂紋沒有。

    對于碰摩故障轉(zhuǎn)子系統(tǒng),國內(nèi)的聶日敏[132]建立了帶有高、低壓渦輪碰摩的航空發(fā)動機雙轉(zhuǎn)子模型,分析雙轉(zhuǎn)子的振動響應,研究結(jié)果發(fā)現(xiàn)轉(zhuǎn)子碰摩會產(chǎn)生激勵頻率,引發(fā)組合共振。侯磊[133]建立了水平盤旋飛行環(huán)境下的轉(zhuǎn)子系統(tǒng)模型,對系統(tǒng)的全局分叉進行分析,發(fā)現(xiàn)當系統(tǒng)接近2 倍臨界轉(zhuǎn)速時,機動載荷會誘發(fā)系統(tǒng)的亞諧共振,進而導致動靜碰摩。Hou 等[134]在碰摩轉(zhuǎn)子系統(tǒng)的基礎(chǔ)上引入了機動載荷,通過模型分析發(fā)現(xiàn),由于機動載荷作用,系統(tǒng)會發(fā)生次諧波共振,從而導致碰摩。

    第二方面是通過實驗分析,進行航空發(fā)動機結(jié)構(gòu)共振類故障的診斷方法研究。在該方向,國內(nèi)的馬會防等[135]基于某風扇增壓級振動故障,研究發(fā)現(xiàn)動靜碰摩造成了行波共振,使得振動突增。王海霞等[136]基于航空發(fā)動機附件機匣振動超標問題,分析了具體的超標原因為固有頻率共振,通過改變固有頻率的措施,降低了振動。鄭旭東和張連祥[124]介紹了因甩油孔位置不當引發(fā)的自激振動故障,說明了共振故障對發(fā)動機的危害之大,在此基礎(chǔ)上提出了改進方案,證明了有效性。

    目前針對航空發(fā)動機結(jié)構(gòu)共振類故障的研究多集中在故障轉(zhuǎn)子系統(tǒng)的動力學建模、結(jié)構(gòu)共振類故障的振動特性、結(jié)構(gòu)共振類故障診斷方法研究等方面。有關(guān)航空發(fā)動機結(jié)構(gòu)共振類故障的診斷方法可以總結(jié)為:利用臨界轉(zhuǎn)速狀態(tài)的振動趨勢圖,結(jié)合頻譜分析綜合診斷結(jié)構(gòu)共振類故障。

    現(xiàn)有研究工作中對航空發(fā)動機結(jié)構(gòu)共振類故障的診斷方法考慮較少,使得目前故障診斷方法無法滿足航空發(fā)動機結(jié)構(gòu)共振類故障的診斷需求,未來的航空發(fā)動機結(jié)構(gòu)共振類故障診斷的主要研究方向是:

    1) 目前,針對航空發(fā)動機結(jié)構(gòu)共振類故障建立的動力學模型多以轉(zhuǎn)子系統(tǒng)為研究對象,但真實情況下,信號多采集自航空發(fā)動機機匣測點。而轉(zhuǎn)子系統(tǒng)響應信號與機匣采集信號間存在明顯差異。未來可在現(xiàn)有轉(zhuǎn)子系統(tǒng)響應信號的基礎(chǔ)上,研究響應信號與機匣采集信號間的傳遞關(guān)系,為后續(xù)故障診斷研究奠定基礎(chǔ)。

    2) 鑒于轉(zhuǎn)子系統(tǒng)響應信號與機匣采集信號間的差異,未來需要建立考慮薄壁機匣結(jié)構(gòu)特性的航空發(fā)動機動力學模型,為后續(xù)故障特征研究奠定基礎(chǔ)。

    3) 現(xiàn)有結(jié)構(gòu)共振類故障案例積累較少,因此為實現(xiàn)對該故障的診斷,未來可將結(jié)構(gòu)共振類故障案例與小樣本學習相結(jié)合,通過小樣本訓練模型,實現(xiàn)對航空發(fā)動機結(jié)構(gòu)共振類故障的分析和故障診斷。

    2.7 主軸承損傷

    航空發(fā)動機主軸承是發(fā)動機的關(guān)鍵部件,主軸承運行在極端苛刻工況及復雜環(huán)境條件下,承受著高溫、高速、重載、貧油、斷油等極端工況,失效形式十分復雜且可能出現(xiàn)多種失效并存,主要失效模式包括疲勞、磨損、熱損傷、打滑蹭傷、保持架失效等。如圖10[137]所示。在振動信號檢測方面,由于航空發(fā)動機結(jié)構(gòu)限制,振動傳感器布置嚴重受限,通常只能安裝在外機匣表面。軸承振動信號的傳遞路徑長且傳遞部件多為薄壁復雜結(jié)構(gòu),導致其頻率成分復雜且信噪比極低。在相關(guān)研究中,陳果等[138]通過對比分析機匣測點信號與軸承座測點信號,揭示了在機匣與軸承連接剛度較小時,機匣信號在傳遞路徑中會發(fā)生較大的衰減。但通過選擇合適的方法依舊可以進行較為準確的診斷。

    圖10 主軸承主要失效模式[137]Fig.10 Main failure modes of bearings[137]

    目前,主軸承故障診斷的主要方法有4 類:面向故障特征頻率提取的振動信號分析、多傳感器協(xié)同分析診斷、基于人工智能的故障辨識以及軸承剩余壽命預測。

    通過提取振動信號的故障特征頻率,計算時域特征、頻域特征以及時頻域特征可以對軸承狀態(tài)進行診斷。2015年在法國羅安理工學院舉辦的賽峰競賽提供了民用飛機發(fā)動機的2 個損壞軸承振動和轉(zhuǎn)速信號[139]。為準確診斷軸承故障,普遍采用的方法是首先在角域?qū)π盘栠M行重采樣,然后采用不同濾波技術(shù)削弱背景噪聲對軸承故障特征信號的影響。圍繞航空發(fā)動機主軸承中的中介軸承,廖明夫等[140]提出了轉(zhuǎn)差域頻譜和轉(zhuǎn)差域包絡(luò)譜的概念,并以此分析了中介軸承外圈故障的信號特征。北京化工大學診斷與自愈工程研究中心(DSE 中心)也建立了相應的雙轉(zhuǎn)子故障模擬實驗臺,并基于主軸承(含中介軸承)信號傳遞特性,建立了一系列主軸承故障診斷方法,包括預白化、共振解調(diào)、轉(zhuǎn)差域階次跟蹤等預處理[2],變分模態(tài)分解[141]、沖擊脈沖法[142]等特征分析,卷積神經(jīng)網(wǎng)絡(luò)[143]、支持向量機[144]等智能診斷方法。

    此外,Zhang 等[145]提出了一種軸承故障診斷方法,在保持相似性以及特征結(jié)構(gòu)完整性的同時,提高了高轉(zhuǎn)速下主軸承故障的診斷精度。Rzadkowski等[146]通過葉尖計時和葉尖間隙技術(shù)分析了渦噴發(fā)動機中介軸承故障,為渦噴發(fā)動機軸承故障分析提供了新的思路。西安交通大學的陳雪峰等[147]提出了匹配同步壓縮變換、匹配解調(diào)變換[148]、瞬時頻帶與同步壓縮方法[149]方法,顯著提升了時頻分析方法的聚集性與航空發(fā)動機振動監(jiān)測的有效性,并經(jīng)過某型航空發(fā)動機整機試車驗證,驗證了方法在振動監(jiān)測與振動突跳溯源排故中具有不可替代的作用。沙云東等[150]提出了一種基于閾值參數(shù)篩選的航空發(fā)動機主軸承故障特征提取方法。西安交通大學Wang 等[151]提出了非凸正則化稀疏診斷、脊加權(quán)稀疏診斷方法[152],顯著提升了航空發(fā)動機軸承剝落故障特征提取與診斷能力。

    由于振動信號采集方便,國內(nèi)外學者也提出了許多基于振動信號的航空發(fā)動機故障診斷方法,并且通過實驗證實了這些方法的可靠性。目前該方法對于機械系統(tǒng)進行故障診斷的理論已較為成熟,然而面對航空發(fā)動機高溫、高應力、油霧環(huán)境、工況復雜、高背景噪聲等條件,EMD、FFT、沖擊脈沖法等單一方法不能有效提取故障特征,需要根據(jù)應用環(huán)境進行改進。

    多傳感器信息融合在主軸承故障診斷中亦受到了廣泛關(guān)注。主軸承故障特征弱、背景噪聲強,常規(guī)監(jiān)測手段僅通過單個傳感器收集狀態(tài),而多傳感器系統(tǒng)可將不同特征進行合理支配,將冗余互補信息進行融合,最終輸出更有用的信息。鑒于上述思路,林桐等[153]提出了基于標準化歐氏距離的多特征融合評估方法,該方法對于航空發(fā)動機軸承故障的不同類型、不同程度診斷有較高靈敏度。Ma 等[154]提出了一種基于弱磁檢測的發(fā)動機軸承故障檢測方法,可通過無損非接觸方式實現(xiàn)振動信號和內(nèi)圈、保持架、滾動元件旋轉(zhuǎn)頻率的同步采集,進而實現(xiàn)多傳感器協(xié)同的軸承故障檢測。目前多傳感器信息融合的方法已經(jīng)成為設(shè)備狀態(tài)檢測的一種趨勢,主要的融合方式有數(shù)據(jù)融合、特征融合、決策融合3 種方式。但是在其增強了振動信號有效性的同時,也帶來了運算量較大,原始數(shù)據(jù)信息損失等問題。

    基于人工智能的航空發(fā)動機主軸承故障辨識也正在一個主要的研究方向,該類方法主要包括神經(jīng)網(wǎng)絡(luò)、遺傳算法、支持向量機等。Wang等[155]提出了一種基于支持矢量機(SVM)的智能診斷方法,用于定量診斷高速航空發(fā)動機軸承的早期微弱故障。張向陽等[156]提出了基于卷積神經(jīng)網(wǎng)絡(luò)(CNN)的故障診斷方法,利用峭度圖等方法將原始信號轉(zhuǎn)變?yōu)閳D像信號,通過卷積神經(jīng)網(wǎng)絡(luò)對故障類型進行識別。為了滿足航空發(fā)動機滾動軸承在線監(jiān)測的要求,Lin 等[157]提出了一種新的超球面距離判別方法(HDD),將振動加速度信號中提取的原始多維特征轉(zhuǎn)換為分布在超球面的同維重構(gòu)特征,利用距離判別分析模型實現(xiàn)軸承故障檢測和退化評估。人工智能的監(jiān)測方法主要依賴于觀測數(shù)據(jù),與系統(tǒng)的復雜程度無關(guān),但其診斷精度卻依賴于樣本數(shù)據(jù)的完整性,需要對各種狀態(tài)下的監(jiān)測數(shù)據(jù)進行學習。

    對軸承進行剩余壽命預測估計軸承狀態(tài)亦是航空發(fā)動機主軸承故障狀態(tài)識別的未來發(fā)展方向。西安交通大學李乃鵬等[158]為了對軸承非線性隨機退化行為進行表征,對傳統(tǒng)的指數(shù)隨機退化模型進行改進,引入貝葉斯更新和粒子濾波方法對退化模型參數(shù)和軸承損傷狀態(tài)進行聯(lián)合更新,降低了退化數(shù)據(jù)的隨機波動對壽命預測的干擾;之后,進一步考慮了運行工況時變對軸承退化速率和監(jiān)測數(shù)據(jù)幅值的影響,構(gòu)建了雙尺度工況轉(zhuǎn)換的壽命預測方法,解決了時變工況下的滾動軸承壽命預測問題[159];同時針對工程實際中很難獲取完整的全壽命數(shù)據(jù)的問題,提出了一種自數(shù)據(jù)驅(qū)動的壽命預測方法框架,基于軸承自身部分退化數(shù)據(jù)從多樣化模型庫中自適應優(yōu)選最佳模型用于退化過程的描述和預測[160]。

    綜上,總結(jié)第2 節(jié)內(nèi)容,目前研究進展匯總?cè)绫?所示。

    表2 故障特征及診斷方法研究現(xiàn)狀匯總Table 2 Summary of current research status on fault characteristics and diagnostic methods

    3 整機振動故障診斷技術(shù)的發(fā)展方向

    航空發(fā)動機這一復雜對象的實際特點,決定了常規(guī)的轉(zhuǎn)子系統(tǒng)振動故障診斷方法并不適用。雖然針對上述常見故障的振動機理、診斷方法研究較多,但大多通過嚴重簡化的實驗器進行驗證,鮮見于成果在實際數(shù)據(jù)中驗證及應用的公開報道,因此,航空發(fā)動機的振動故障診斷還面臨寬頻高精度振動采集、振動機理與傳遞路徑分析、機理與案例數(shù)據(jù)深度融合、機載智能診斷與地面大數(shù)據(jù)挖掘協(xié)同等諸多挑戰(zhàn)。在總結(jié)上述工程應用、理論研究現(xiàn)狀的基礎(chǔ)上,借鑒國外發(fā)展經(jīng)驗,梳理航空發(fā)動機整機振動故障診斷技術(shù)的發(fā)展方向及已取得典型進展如下。

    3.1 寬頻高精度振動采集是技術(shù)發(fā)展的根基

    以主軸承損傷、葉片斷裂為代表的故障特征主要集中在數(shù)千甚至上萬赫茲的高頻帶,而不平衡、不同心等轉(zhuǎn)子類故障特征信號主要集中在幾十至數(shù)百赫茲的中低頻,因此需采集寬頻振動信號,分析特定故障特征分量的變化,以實現(xiàn)航空發(fā)動機典型振動故障診斷。

    同時,由于航空發(fā)動機振動測點多布置在外機匣上,轉(zhuǎn)子系統(tǒng)及支承發(fā)生故障時,通過復雜路徑(如圖11[161]所示)傳遞到測點位置時已極其微弱;并且,由于航空發(fā)動機的頻繁、大幅變工況特點,其振動信號通常為具有強時變非平穩(wěn)特性的快變信號,需要對快變過程每個瞬態(tài)的振動信號進行精準采集分析。

    圖11 某型雙轉(zhuǎn)子航空發(fā)動機支承結(jié)構(gòu)示意圖[161]Fig.11 Schematic diagram of support structure for a certain type of the aero-engine[161]

    綜上,發(fā)展航空發(fā)動機振動故障診斷技術(shù)以實現(xiàn)寬頻高精度振動感知、快變條件下微弱特征分量的精準提取為根基。國內(nèi)相關(guān)單位亦開展了相關(guān)探索研究,但機載的推廣應用還需進一步加強:如針對C919客機自主研發(fā)的PHM 系統(tǒng)[162]、廣州航新科技公司針對直升機開發(fā)的HUMS 系統(tǒng)[163]等。北京化工大學DSE 中心亦致力于整機振動監(jiān)測診斷技術(shù)的研究和應用工作,目前已取得豐碩的成果。PHM 系統(tǒng)要求有足夠的傳感器才能獲得實現(xiàn)預測與健康管理所需的底層信息[164-166],DSE 中心與某發(fā)動機保障單位合作,綜合考慮振動敏感性與工程可行性[167-168],實現(xiàn)了某型航空發(fā)動機前機匣、中介機匣、后機匣共7 個測點的寬頻高精度振動狀態(tài)監(jiān)測,已積累航空發(fā)動機寬頻高質(zhì)量數(shù)據(jù)2 000 余臺次。

    3.2 整機試驗及案例是技術(shù)發(fā)展的前提

    航空發(fā)動機是試驗科學的產(chǎn)物,故障診斷更需要實踐和試驗來指導。整機振動故障診斷技術(shù)研究主要涵蓋4 個方面:模型、專家知識、信號處理和數(shù)據(jù),如圖12[168]所示。

    圖12 故障診斷的主要方法[168]Fig.12 Main methods of fault diagnosis[168]

    基于模型的故障診斷,即通過建立發(fā)動機模型用以描述系統(tǒng)或部件的力學、數(shù)學特性,輸入動態(tài)激勵后求解系統(tǒng)故障響應,以此指導故障診斷[169-170],而發(fā)動機的故障案例即是驗證模型仿真結(jié)果準確性的重要依據(jù)。

    基于專家知識的故障診斷,即以發(fā)動機故障診斷和維修保障等領(lǐng)域?qū)<覍W者多年積累的系統(tǒng)性知識和經(jīng)驗為基礎(chǔ),構(gòu)建典型故障診斷規(guī)則庫,進而指導故障診斷[171-172],發(fā)動機的故障案例是專家知識形成的主要依據(jù)。

    基于信號處理的故障診斷,即對振動信號進行時域、頻域、時頻特征分析,對各個振動分量進行有效分離和特征提取,結(jié)合診斷知識對提取出的特征進行物理解釋,用于故障診斷[173-174],發(fā)動機故障案例是信號處理的主要對象。

    基于數(shù)據(jù)的診斷方法,即通過學習算法挖掘隱藏在案例數(shù)據(jù)中的規(guī)律和特征,利用人工智能實現(xiàn)發(fā)動機狀態(tài)評估與故障識別[175-176],發(fā)動機的故障案例是機器學習的基礎(chǔ)。

    綜上所述,發(fā)動機的故障案例在這4 方面中均發(fā)揮著重要作用,是航空發(fā)動機整機振動故障診斷技術(shù)發(fā)展的前提。然而,發(fā)動機的故障案例存在樣本少且分布極其不均的問題。因此,需要開展航空發(fā)動機整機試驗,在最大程度還原發(fā)動機故障狀態(tài)的情況下,補充故障案例,支撐振動診斷技術(shù)發(fā)展。

    為彌補整機試驗及案例高質(zhì)量數(shù)據(jù)的缺失,相關(guān)單位已初步開展了2 方面工作:① 整機故障模擬。北京化工大學DSE 中心、某研究所分別以某型渦扇發(fā)動機為原型,建設(shè)了電驅(qū)真機故障模擬試驗裝置,主要用于在整機狀態(tài)下開展軸承損傷、轉(zhuǎn)子不平衡、動靜碰摩、葉片損傷等典型振動故障模擬[177];② 真機案例的高質(zhì)量數(shù)據(jù)積累。北京化工大學DSE 中心與某保障單位、中國航發(fā)沈陽發(fā)動機研究所、中國航發(fā)株洲動力機械研究所、中國航發(fā)沈陽黎明航空發(fā)動機公司等單位合作,成功捕獲了500 余起典型振動故障案例;西安交通大學陳雪峰等[147]通過與中國燃氣渦輪研究院、中國航發(fā)西安航空發(fā)動機有限公司等單位合作,亦積累了航空發(fā)動機試車故障案例。上述數(shù)據(jù)的積累,可為診斷技術(shù)研究奠定良好的數(shù)據(jù)基礎(chǔ)。

    3.3 振動機理與傳遞路徑分析是技術(shù)發(fā)展的基礎(chǔ)

    發(fā)動機整機振動機理研究是航空領(lǐng)域的重大課題[178]。發(fā)動機轉(zhuǎn)子-支承-機匣系統(tǒng)高度耦合,且受發(fā)動機結(jié)構(gòu)和工作條件的制約,振動傳感器大多分布于外機匣,因此分析典型故障振動機理、信號傳遞機理格外重要[179]。以某型雙轉(zhuǎn)子航空發(fā)動機低壓壓氣機工作葉片和盤產(chǎn)生的不平衡、不同心、葉片斷裂等故障激勵為例,該振動載荷通過鼓筒傳遞至前、后軸頸,再分別通過1 號軸承和2 號軸承傳遞到低壓壓氣機靜子機匣、中介機匣,信號經(jīng)復雜薄壁結(jié)構(gòu)、彈性支承等部件傳遞,振動響應必然會發(fā)生顯著改變。因此,若想實現(xiàn)航空發(fā)動機整機振動監(jiān)測與故障診斷,首先必須以整機振動的激振源為出發(fā)點,結(jié)合發(fā)動機整機結(jié)構(gòu)動力學,分析故障的發(fā)生發(fā)展機理、振動響應的傳遞規(guī)律。

    關(guān)于振動機理相關(guān)研究,本文第2 節(jié)已進行總結(jié),此處不再贅述。在傳遞路徑分析方面,天津大學[180]、哈爾濱工業(yè)大學[181]等都已形成了部分成果,但鮮見在實際航空發(fā)動機上應用的公開報道。故障機理及傳遞路徑的研究,還需結(jié)合實際發(fā)動機結(jié)構(gòu)及運行特點,通過整機試驗、測試數(shù)據(jù)來繼續(xù)深入。

    3.4 機理與案例數(shù)據(jù)深度融合是技術(shù)發(fā)展的主要途徑

    航空發(fā)動機的復雜結(jié)構(gòu)與多變工況,使得其與傳統(tǒng)旋轉(zhuǎn)機械差異顯著,傳統(tǒng)診斷理論與方法難以直接應用。數(shù)據(jù)驅(qū)動的診斷模型需要大量高質(zhì)量的帶標簽故障樣本作為訓練條件,而積累了高質(zhì)量數(shù)據(jù)的發(fā)動機故障案例樣本匱乏,且人為標記的案例標簽未必精確,導致其難以在航空發(fā)動機振動故障診斷中直接應用;同時,僅基于數(shù)據(jù)訓練得到的模型由于其處于“黑箱”狀態(tài),其結(jié)果難以令人信服,也不利于進行故障源分析及方法優(yōu)化改進。因此,結(jié)合機理知識的可解釋性與實際案例數(shù)據(jù)的真實性,構(gòu)建機理與數(shù)據(jù)深度融合模型,可實現(xiàn)準確、可靠的發(fā)動機整機振動故障診斷。

    該深度融合方式可被大致概括為2 個方向:基于機理信息的嵌入算法與數(shù)字孿生?;跈C理信息的嵌入算法是將物理知識融入機器學習過程中,利用機理知識對數(shù)據(jù)驅(qū)動的隨機性起到限制作用[182],西安交通大學Wang 等[183-185]提出了模型驅(qū)動的算法展開深度網(wǎng)絡(luò)構(gòu)造方法,不僅網(wǎng)絡(luò)設(shè)計具有事前可解釋性,網(wǎng)絡(luò)學習結(jié)果亦具有事后可解釋性,學習出了與故障機理相對應的網(wǎng)絡(luò)卷積參數(shù)以及故障特征。而數(shù)字孿生則是一個更為廣闊的方向,不止包含有效物理信息的輸入,也包括實際案例數(shù)據(jù)對模型的參數(shù)識別,通過數(shù)字孿生模型將能為故障辨識提供與物理對象相同、甚至更多的信息,輔助數(shù)據(jù)驅(qū)動模型得到更好的學習結(jié)果。

    為建立適用于航空發(fā)動機的故障診斷技術(shù),依托國家重大項目支持,北京化工大學、西安交通大學、西北工業(yè)大學、中國航發(fā)等單位均開展了大量研究,但鮮見應用于真機試車、飛行中的實時分析[19]。

    3.5 機載智能診斷與地面大數(shù)據(jù)挖掘協(xié)同是技術(shù)發(fā)展的目標

    隨著現(xiàn)代振動信號處理、特征提取工程、人工智能技術(shù)的不斷發(fā)展,航空發(fā)動機振動故障診斷技術(shù)也不斷拓展、繁榮。傳統(tǒng)事后維修、定期檢修方式,費時費力,耗資居高不下,難以保障發(fā)動機長期有效工作,已無法滿足現(xiàn)代飛機發(fā)動機故障診斷和維護需求,自動化、信息化、智能化已成為一種剛需,不只是地面檢修維護,更亟需在空中飛行過程便對發(fā)動機振動故障進行智能告警與診斷。然而,受限于機載計算資源有限、難以實現(xiàn)多臺次發(fā)動機的綜合對比分析等條件,航空發(fā)動機的PHM 體系通常還包括地面站的大數(shù)據(jù)分析平臺。

    針對于此,綜合利用各項信號分析處理、機理模型構(gòu)建及診斷技術(shù),機載端實現(xiàn)典型故障的智能告警與快速智能診斷,地面站實現(xiàn)機群間對比分析、大數(shù)據(jù)挖掘及溯源診斷,兩者協(xié)同,實現(xiàn)航空發(fā)動機飛行過程中實時監(jiān)測與快準診斷,及時準確地檢測故障發(fā)生并提供相關(guān)告警,以輔助進行飛行決策,降低重大事故發(fā)生概率,保障飛行人員生命安全。同時,地面站可進行故障溯源與定位,指導后續(xù)靶向保障,轉(zhuǎn)被動、事后維修為狀態(tài)維修,縮短航空發(fā)動機維護時間,減少定期維修,提高發(fā)動機出勤率。

    因此,機載智能診斷與地面大數(shù)據(jù)挖掘協(xié)同是航空發(fā)動機整機振動故障監(jiān)測診斷的發(fā)展目標,也是未來智能發(fā)動機的迫切需求。針對這一難題,中國航發(fā)沈陽發(fā)動機研究所、中國航發(fā)株洲動力機械研究所、中國航發(fā)貴陽發(fā)動機設(shè)計研究所等工業(yè)部門正在聯(lián)合高校開展相關(guān)理論研究與應用探索。

    3.6 故障隔離與整機抑制是振動故障的快準解決策略

    故障隔離是指當航空發(fā)動機某部分(或某臺)發(fā)生故障時,采取停機/降速運行等技術(shù)手段,將故障隔離起來,使其對其他部分(或雙發(fā)的另一臺發(fā)動機)盡可能不產(chǎn)生/少產(chǎn)生影響,進而降低故障的危害性。

    整機抑制包括微創(chuàng)抑振與自愈調(diào)控。微創(chuàng)抑振指發(fā)動機在地面保障過程中,采取整機動平衡、局部分解更換部件等方式,在不對發(fā)動機進行大范圍拆解的條件下治愈振動故障;自愈調(diào)控指發(fā)動機在實際運行中,采取自動平衡、損傷自修復等技術(shù),在運行中“自行”抑制或消除故障。

    顯然,上述技術(shù)是振動故障得到精準診斷后的快準解決策略,可極大保障發(fā)動機及飛機的安全,并顯著提高維修效率。針對這一重大難題,北京化工大學聯(lián)合中國航發(fā)沈陽發(fā)動機研究所、中國航發(fā)株洲動力機械研究所形成了2 項成果:整機動平衡技術(shù),真機試車驗證了其有效性;積極探索了自動平衡方法,在部件試驗器上驗證了其抑振效果。

    猜你喜歡
    故障診斷發(fā)動機振動
    振動的思考
    科學大眾(2023年17期)2023-10-26 07:39:14
    振動與頻率
    天天愛科學(2020年6期)2020-09-10 07:22:44
    發(fā)動機空中起動包線擴展試飛組織與實施
    中立型Emden-Fowler微分方程的振動性
    因果圖定性分析法及其在故障診斷中的應用
    新一代MTU2000發(fā)動機系列
    UF6振動激發(fā)態(tài)分子的振動-振動馳豫
    計算物理(2014年2期)2014-03-11 17:01:44
    基于LCD和排列熵的滾動軸承故障診斷
    基于WPD-HHT的滾動軸承故障診斷
    機械與電子(2014年1期)2014-02-28 02:07:31
    高速泵的故障診斷
    河南科技(2014年3期)2014-02-27 14:05:48
    如何舔出高潮| 在线观看三级黄色| av免费观看日本| 精品一区二区三卡| 在线免费观看不下载黄p国产| 成人一区二区视频在线观看| 国产精品嫩草影院av在线观看| 欧美xxxx黑人xx丫x性爽| 久久久久久久午夜电影| 色哟哟·www| 中文精品一卡2卡3卡4更新| 成人黄色视频免费在线看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲美女搞黄在线观看| 欧美成人午夜免费资源| 国产精品久久久久久久久免| 久久久久精品性色| 69av精品久久久久久| videossex国产| 午夜老司机福利剧场| 亚洲最大成人手机在线| 国产精品av视频在线免费观看| 身体一侧抽搐| 国产又色又爽无遮挡免| 男男h啪啪无遮挡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 永久网站在线| 国产欧美亚洲国产| 男人爽女人下面视频在线观看| 精华霜和精华液先用哪个| 国产日韩欧美在线精品| 51国产日韩欧美| 日本爱情动作片www.在线观看| 精品久久久久久电影网| 亚洲国产精品专区欧美| 欧美精品人与动牲交sv欧美| 亚洲自偷自拍三级| 亚洲高清免费不卡视频| 内地一区二区视频在线| 国产精品一区二区在线观看99| 午夜日本视频在线| 少妇的逼好多水| 国语对白做爰xxxⅹ性视频网站| 欧美成人精品欧美一级黄| 97超碰精品成人国产| 日韩欧美 国产精品| 免费电影在线观看免费观看| 日韩欧美精品免费久久| 老师上课跳d突然被开到最大视频| www.av在线官网国产| 2022亚洲国产成人精品| 欧美精品国产亚洲| 老司机影院成人| 伊人久久国产一区二区| av天堂中文字幕网| 国产精品久久久久久精品电影| 久久精品综合一区二区三区| 午夜爱爱视频在线播放| 小蜜桃在线观看免费完整版高清| 交换朋友夫妻互换小说| 国产成年人精品一区二区| 亚洲精品久久午夜乱码| 晚上一个人看的免费电影| 免费看不卡的av| 身体一侧抽搐| 亚洲欧洲日产国产| 久久精品夜色国产| 人妻 亚洲 视频| 高清午夜精品一区二区三区| 天天躁夜夜躁狠狠久久av| 乱码一卡2卡4卡精品| 国产伦精品一区二区三区四那| 亚洲国产高清在线一区二区三| 国产成人91sexporn| 国产成人a∨麻豆精品| 久久久久久久大尺度免费视频| 人人妻人人爽人人添夜夜欢视频 | 联通29元200g的流量卡| 亚洲av欧美aⅴ国产| 免费不卡的大黄色大毛片视频在线观看| 女人十人毛片免费观看3o分钟| 一区二区三区四区激情视频| 久久精品久久久久久久性| 特大巨黑吊av在线直播| 国产又色又爽无遮挡免| 国产极品天堂在线| 青青草视频在线视频观看| 国产乱来视频区| 久久久精品免费免费高清| 国产爱豆传媒在线观看| 国产精品爽爽va在线观看网站| 日韩一区二区三区影片| 大香蕉久久网| 久久久久久伊人网av| 深夜a级毛片| 99九九线精品视频在线观看视频| 欧美 日韩 精品 国产| 国产探花极品一区二区| 成人国产麻豆网| 国产午夜福利久久久久久| 麻豆精品久久久久久蜜桃| 在现免费观看毛片| 国产精品国产三级专区第一集| 久久久午夜欧美精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久6这里有精品| 在现免费观看毛片| 亚洲精品乱码久久久v下载方式| 成人午夜精彩视频在线观看| 十八禁网站网址无遮挡 | 综合色丁香网| 久久久精品欧美日韩精品| 久久人人爽av亚洲精品天堂 | 亚洲av福利一区| 又粗又硬又长又爽又黄的视频| 中文字幕av成人在线电影| 成人毛片60女人毛片免费| 国产午夜福利久久久久久| 亚洲成人av在线免费| 色综合色国产| 天天一区二区日本电影三级| 国产免费一级a男人的天堂| 99热网站在线观看| 日本欧美国产在线视频| 亚洲精华国产精华液的使用体验| 欧美日韩精品成人综合77777| 丝袜美腿在线中文| 又黄又爽又刺激的免费视频.| 99视频精品全部免费 在线| 亚洲精品色激情综合| 成年人午夜在线观看视频| 看黄色毛片网站| 2022亚洲国产成人精品| 国产精品久久久久久av不卡| 国产免费又黄又爽又色| 国产免费一级a男人的天堂| 综合色av麻豆| 亚洲精品国产av蜜桃| 一个人看的www免费观看视频| 亚洲综合色惰| 国产精品99久久99久久久不卡 | 久久人人爽人人爽人人片va| 午夜激情福利司机影院| 五月开心婷婷网| 久久久久久久久大av| 亚洲婷婷狠狠爱综合网| 成人一区二区视频在线观看| 日本午夜av视频| 国产成人a∨麻豆精品| 久久精品国产a三级三级三级| 精品少妇黑人巨大在线播放| 麻豆久久精品国产亚洲av| 欧美日韩精品成人综合77777| 亚洲欧美成人综合另类久久久| 精品午夜福利在线看| 春色校园在线视频观看| 久久精品国产亚洲av天美| 黄色一级大片看看| 亚洲欧洲国产日韩| 久久综合国产亚洲精品| 亚洲欧美清纯卡通| 极品少妇高潮喷水抽搐| 久久久久国产网址| 亚洲精品日韩在线中文字幕| 亚洲自拍偷在线| 欧美日韩精品成人综合77777| 视频区图区小说| 亚洲精品色激情综合| 国产高清国产精品国产三级 | 狂野欧美白嫩少妇大欣赏| 在线观看一区二区三区激情| 一个人看的www免费观看视频| 久久精品国产自在天天线| 午夜免费鲁丝| 人妻系列 视频| 成人美女网站在线观看视频| 欧美区成人在线视频| 国产爱豆传媒在线观看| 日韩中字成人| 国产国拍精品亚洲av在线观看| 中国三级夫妇交换| 最近手机中文字幕大全| 国产欧美日韩一区二区三区在线 | 夜夜看夜夜爽夜夜摸| 少妇人妻久久综合中文| 老司机影院成人| 国产色婷婷99| 国产爱豆传媒在线观看| 亚洲国产欧美人成| 久久99热6这里只有精品| 国产成人一区二区在线| 亚洲一区二区三区欧美精品 | 婷婷色av中文字幕| tube8黄色片| 熟女电影av网| 18禁裸乳无遮挡动漫免费视频 | 蜜桃久久精品国产亚洲av| 日本一本二区三区精品| 伊人久久精品亚洲午夜| av在线观看视频网站免费| 一边亲一边摸免费视频| 欧美区成人在线视频| 欧美一级a爱片免费观看看| 国产白丝娇喘喷水9色精品| 91aial.com中文字幕在线观看| 国产女主播在线喷水免费视频网站| 免费观看a级毛片全部| 色吧在线观看| 嫩草影院新地址| 18禁动态无遮挡网站| 久久久久国产网址| 国产精品无大码| 午夜激情福利司机影院| 男人狂女人下面高潮的视频| 中文字幕亚洲精品专区| 联通29元200g的流量卡| 久久久亚洲精品成人影院| 九色成人免费人妻av| 99久国产av精品国产电影| 一级爰片在线观看| 亚洲人成网站在线观看播放| 全区人妻精品视频| 日韩中字成人| 国产成人免费无遮挡视频| 国产69精品久久久久777片| av在线app专区| 久久鲁丝午夜福利片| 深夜a级毛片| 国国产精品蜜臀av免费| 国产精品嫩草影院av在线观看| 国产乱人偷精品视频| 91狼人影院| 蜜臀久久99精品久久宅男| 大码成人一级视频| 日韩成人av中文字幕在线观看| 国产片特级美女逼逼视频| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久久成人| 神马国产精品三级电影在线观看| 激情 狠狠 欧美| 日本黄色片子视频| 国产精品不卡视频一区二区| 亚洲真实伦在线观看| xxx大片免费视频| 国产美女午夜福利| 亚洲精华国产精华液的使用体验| 一级毛片aaaaaa免费看小| 亚洲欧美一区二区三区黑人 | 精品久久久久久久久av| 国产黄片视频在线免费观看| 久久久久国产精品人妻一区二区| 国产黄色视频一区二区在线观看| 精品久久久噜噜| 色网站视频免费| 色吧在线观看| 一级黄片播放器| 特级一级黄色大片| 亚洲精品一二三| 男女边摸边吃奶| 国产精品女同一区二区软件| 一区二区三区乱码不卡18| 亚洲av在线观看美女高潮| .国产精品久久| 在线观看三级黄色| 亚洲欧美日韩卡通动漫| 99热6这里只有精品| 欧美 日韩 精品 国产| 国产综合懂色| 寂寞人妻少妇视频99o| 亚洲怡红院男人天堂| 少妇裸体淫交视频免费看高清| 婷婷色av中文字幕| 国产白丝娇喘喷水9色精品| 日韩欧美精品v在线| 一级二级三级毛片免费看| 最近中文字幕高清免费大全6| 国产淫片久久久久久久久| 一级爰片在线观看| 黄色一级大片看看| 免费av毛片视频| 亚洲高清免费不卡视频| 中文字幕人妻熟人妻熟丝袜美| 波野结衣二区三区在线| 菩萨蛮人人尽说江南好唐韦庄| 男人和女人高潮做爰伦理| 免费看不卡的av| 韩国av在线不卡| 欧美变态另类bdsm刘玥| av在线蜜桃| 青青草视频在线视频观看| 尤物成人国产欧美一区二区三区| 99热这里只有是精品50| 寂寞人妻少妇视频99o| 一级毛片 在线播放| 春色校园在线视频观看| 18+在线观看网站| 亚洲伊人久久精品综合| 精品久久久久久久末码| 亚洲在线观看片| 欧美亚洲 丝袜 人妻 在线| 久久ye,这里只有精品| 欧美 日韩 精品 国产| 国内精品美女久久久久久| 久久综合国产亚洲精品| 男女那种视频在线观看| 日韩欧美精品v在线| 国产亚洲av嫩草精品影院| 人人妻人人爽人人添夜夜欢视频 | 九草在线视频观看| 人妻夜夜爽99麻豆av| 日本猛色少妇xxxxx猛交久久| 亚洲人与动物交配视频| 女人久久www免费人成看片| 99久久九九国产精品国产免费| 国产精品av视频在线免费观看| 午夜爱爱视频在线播放| 亚洲色图av天堂| 国产一区二区亚洲精品在线观看| 亚洲人成网站在线播| 男女啪啪激烈高潮av片| 一个人观看的视频www高清免费观看| 成年女人在线观看亚洲视频 | 一级av片app| 插逼视频在线观看| 国产亚洲最大av| 中文在线观看免费www的网站| 一区二区三区精品91| 欧美日韩国产mv在线观看视频 | 欧美区成人在线视频| 精品一区在线观看国产| av黄色大香蕉| 91aial.com中文字幕在线观看| 插阴视频在线观看视频| 亚洲在线观看片| 黄片wwwwww| 亚洲国产精品成人久久小说| 麻豆成人av视频| 最新中文字幕久久久久| 成年人午夜在线观看视频| 亚洲欧美成人精品一区二区| 日本午夜av视频| 六月丁香七月| av在线app专区| 欧美三级亚洲精品| 日韩伦理黄色片| 超碰av人人做人人爽久久| 午夜激情福利司机影院| 少妇的逼水好多| a级一级毛片免费在线观看| 欧美97在线视频| 久久99精品国语久久久| 久久精品综合一区二区三区| 精品人妻熟女av久视频| 国产精品人妻久久久久久| 亚洲精品日韩在线中文字幕| 一级毛片电影观看| 国产在线一区二区三区精| 大陆偷拍与自拍| 亚洲精品乱久久久久久| 大香蕉97超碰在线| 免费av观看视频| 亚洲精品久久久久久婷婷小说| 夫妻午夜视频| 午夜激情福利司机影院| 免费看不卡的av| 男女下面进入的视频免费午夜| 久久久久久久精品精品| 久久精品国产鲁丝片午夜精品| 汤姆久久久久久久影院中文字幕| 91久久精品国产一区二区成人| 看黄色毛片网站| 国产欧美另类精品又又久久亚洲欧美| 中国三级夫妇交换| 国产精品久久久久久精品电影小说 | 亚洲欧美日韩无卡精品| 美女cb高潮喷水在线观看| 国产毛片a区久久久久| 日韩强制内射视频| 男人添女人高潮全过程视频| www.色视频.com| 久久精品久久久久久久性| 亚洲精品aⅴ在线观看| 一区二区三区四区激情视频| 亚洲丝袜综合中文字幕| www.av在线官网国产| 亚洲国产欧美在线一区| 亚洲国产色片| h日本视频在线播放| 大片电影免费在线观看免费| 国产亚洲精品久久久com| 欧美激情国产日韩精品一区| 日韩av免费高清视频| 免费看日本二区| 汤姆久久久久久久影院中文字幕| 能在线免费看毛片的网站| 国产毛片a区久久久久| 国产亚洲午夜精品一区二区久久 | www.色视频.com| 色视频www国产| 看免费成人av毛片| 欧美一区二区亚洲| 免费看不卡的av| 97在线视频观看| 亚洲人成网站在线播| 亚洲精品视频女| 下体分泌物呈黄色| 全区人妻精品视频| 91久久精品国产一区二区成人| 在线播放无遮挡| 日韩三级伦理在线观看| 久久久久久久国产电影| 亚洲在久久综合| 久久精品综合一区二区三区| 美女视频免费永久观看网站| 久久久成人免费电影| 国产成人免费无遮挡视频| 国产在线一区二区三区精| 久久99蜜桃精品久久| av卡一久久| 国产综合懂色| 国产精品国产三级专区第一集| 欧美国产精品一级二级三级 | 黄色视频在线播放观看不卡| 午夜免费男女啪啪视频观看| 亚洲av免费在线观看| 日日啪夜夜爽| 亚洲精品成人久久久久久| 久久久久性生活片| 日韩欧美 国产精品| 一二三四中文在线观看免费高清| 人妻 亚洲 视频| 久久久久精品久久久久真实原创| 国产精品99久久99久久久不卡 | 亚洲真实伦在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲av成人精品一区久久| 国产真实伦视频高清在线观看| 一边亲一边摸免费视频| 国产在线一区二区三区精| 亚洲精品亚洲一区二区| 亚洲欧美日韩卡通动漫| 成人午夜精彩视频在线观看| 日本午夜av视频| 青春草国产在线视频| 精品国产一区二区三区久久久樱花 | 一级av片app| 亚洲怡红院男人天堂| 久久影院123| 久热久热在线精品观看| 九九在线视频观看精品| 久久人人爽人人爽人人片va| 香蕉精品网在线| 国产熟女欧美一区二区| 美女xxoo啪啪120秒动态图| 国产69精品久久久久777片| 亚洲欧美一区二区三区黑人 | 人人妻人人爽人人添夜夜欢视频 | 午夜福利网站1000一区二区三区| 成人特级av手机在线观看| 亚洲精品456在线播放app| 天堂中文最新版在线下载 | av在线天堂中文字幕| 夜夜爽夜夜爽视频| 精品一区二区免费观看| 欧美xxxx性猛交bbbb| 岛国毛片在线播放| 欧美少妇被猛烈插入视频| 又黄又爽又刺激的免费视频.| 精品国产三级普通话版| 丰满少妇做爰视频| 青春草视频在线免费观看| 国产午夜精品久久久久久一区二区三区| 亚洲最大成人手机在线| 夜夜看夜夜爽夜夜摸| 欧美xxxx黑人xx丫x性爽| 搞女人的毛片| 亚洲欧洲日产国产| 老师上课跳d突然被开到最大视频| 97人妻精品一区二区三区麻豆| 日韩 亚洲 欧美在线| 国产免费福利视频在线观看| 卡戴珊不雅视频在线播放| 看免费成人av毛片| 69人妻影院| 色5月婷婷丁香| 久久久久久久午夜电影| 国内精品美女久久久久久| 亚洲综合色惰| 如何舔出高潮| av在线观看视频网站免费| 国产精品99久久99久久久不卡 | 成人毛片a级毛片在线播放| 大片免费播放器 马上看| 亚洲经典国产精华液单| 久久久久九九精品影院| 九九在线视频观看精品| 91在线精品国自产拍蜜月| 国产又色又爽无遮挡免| 久久久久久久午夜电影| 国产精品国产av在线观看| av又黄又爽大尺度在线免费看| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 18禁在线无遮挡免费观看视频| 国产成人精品一,二区| 久久国产乱子免费精品| 性色av一级| 免费观看a级毛片全部| 亚洲精品中文字幕在线视频 | 国产成人a∨麻豆精品| 少妇高潮的动态图| 免费观看性生交大片5| 黄色一级大片看看| 男男h啪啪无遮挡| 久久女婷五月综合色啪小说 | 成人美女网站在线观看视频| 国产欧美另类精品又又久久亚洲欧美| 国产高清三级在线| 大码成人一级视频| 日韩视频在线欧美| 亚洲国产欧美人成| 亚洲人成网站在线观看播放| 搡老乐熟女国产| 亚洲性久久影院| 69人妻影院| 亚洲成色77777| 最后的刺客免费高清国语| 国产av码专区亚洲av| 久久久久网色| 22中文网久久字幕| 亚洲精品成人av观看孕妇| 欧美日韩精品成人综合77777| 免费少妇av软件| 精品少妇久久久久久888优播| 欧美xxⅹ黑人| 女的被弄到高潮叫床怎么办| 亚洲国产欧美在线一区| 免费大片黄手机在线观看| 男人爽女人下面视频在线观看| 亚洲欧美日韩无卡精品| 精品一区二区三卡| 丝袜美腿在线中文| 久久精品人妻少妇| 又爽又黄a免费视频| 中国国产av一级| 高清午夜精品一区二区三区| 国产69精品久久久久777片| 亚州av有码| 国产成人freesex在线| 国产精品嫩草影院av在线观看| 日本午夜av视频| 日韩电影二区| av卡一久久| 2022亚洲国产成人精品| videos熟女内射| 成人漫画全彩无遮挡| 夫妻性生交免费视频一级片| 国产乱来视频区| 亚洲精品乱久久久久久| 亚洲最大成人中文| 水蜜桃什么品种好| 国产成人91sexporn| av在线亚洲专区| 热99国产精品久久久久久7| 精华霜和精华液先用哪个| 欧美潮喷喷水| 欧美成人一区二区免费高清观看| 2021天堂中文幕一二区在线观| 日本一本二区三区精品| 久久99热这里只频精品6学生| 国产午夜福利久久久久久| 晚上一个人看的免费电影| 九色成人免费人妻av| 国产精品一区二区在线观看99| 97精品久久久久久久久久精品| 成人国产麻豆网| 赤兔流量卡办理| videossex国产| 91aial.com中文字幕在线观看| 欧美97在线视频| 色播亚洲综合网| 蜜臀久久99精品久久宅男| 日本一本二区三区精品| 久久6这里有精品| 欧美日韩一区二区视频在线观看视频在线 | 久久精品国产亚洲av天美| 欧美另类一区| 欧美区成人在线视频| a级毛色黄片| 啦啦啦在线观看免费高清www| 亚洲国产最新在线播放| 久久99精品国语久久久| 精品一区二区三卡| 精品久久久久久久末码| 国产精品无大码| 一级毛片aaaaaa免费看小| 一区二区三区免费毛片| 国产精品偷伦视频观看了| 黄色欧美视频在线观看| 69人妻影院| 成年女人看的毛片在线观看| 一级黄片播放器| 男女下面进入的视频免费午夜| 又粗又硬又长又爽又黄的视频| 夜夜爽夜夜爽视频| 国产精品爽爽va在线观看网站| 亚洲无线观看免费| 国产高清有码在线观看视频| 熟女电影av网| 少妇人妻精品综合一区二区| 哪个播放器可以免费观看大片| 国产乱人偷精品视频| 久久久久久国产a免费观看| 久久鲁丝午夜福利片| 免费观看的影片在线观看|