梁曉峰，鄧熠，司國雷，陳君輝，陳文秀

(1.中國人民解放軍海軍裝備部，四川成都 610100；2.四川航天烽火伺服控制技術有限公司，四川成都 611130)

0 前言

燃油調節(jié)器具備航空發(fā)動機啟停機、供油、按需調節(jié)油量等功能，是發(fā)動機的重要組成部件。而弱狀態(tài)是指燃油調節(jié)器處于內部部件磨損但不影響正常使用的工作狀態(tài)。燃油調節(jié)器的故障具有隱蔽性、長期潛伏性及突發(fā)性，若想提早發(fā)現(xiàn)其隱藏故障，建立能夠描述監(jiān)測信號本質屬性的本征特征空間，以實現(xiàn)對燃油調節(jié)器物理狀態(tài)的監(jiān)測，則尤為重要。然而，因流體固有的物理特性及機-液耦合所帶來的影響，使得人們很難通過常規(guī)手段獲取用于描述其狀態(tài)的信號本征特征，導致其自身物理狀態(tài)呈現(xiàn)較強的模糊性且將長期處于弱狀態(tài)，存在極大的安全隱患。傳統(tǒng)的信號處理方法所提取的信號本征特征在特征空間呈現(xiàn)交錯性分布，且受測試硬件和時間的限制，在一定測試時間內獲取的工作狀態(tài)的信息量有限，導致提取的特征信息不完備。而當利用未知信息判斷燃油調節(jié)器自身物理狀態(tài)時，會因特征信息的缺陷，導致難以進行準確的信息匹配，從而影響狀態(tài)監(jiān)測系統(tǒng)推理機的決策精度。因此，獲取較為完備的狀態(tài)信息并構建互不交錯的特征空間，對燃油調節(jié)器狀態(tài)監(jiān)測系統(tǒng)有重要意義。

在液壓系統(tǒng)狀態(tài)特征提取中，許多學者介紹了LMD、領域粗糙集、經(jīng)驗模態(tài)、徑向基、貝葉斯神經(jīng)網(wǎng)絡、深度神經(jīng)網(wǎng)絡及小波包等算法，將它們作用于壓力、流量、振動等信號中，并結合決策級/特征級數(shù)據(jù)融合技術，取得了良好的特征數(shù)據(jù),但存在效率較低、硬件和時間成本較高的缺點?；谏鲜龇治?，本文作者針對弱狀態(tài)下燃油調節(jié)器自身物理狀態(tài)難以描述的問題，以燃油調節(jié)器出口組件的壓力脈動為目標，利用無監(jiān)督聚類算法和多源信息融合技術對壓力脈動信號進行特征初次提取并設計尺度變化準則，進行特征重構。利用Java編寫特征提取和-NN(-Nearest Neighbor)故障分類器的腳本，進行算法測試。

1 弱狀態(tài)下狀態(tài)表征現(xiàn)象及特征提取策略

1.1 燃油調節(jié)器的弱狀態(tài)表征現(xiàn)象

圖1所示為燃油調節(jié)器工作資源的分配方式。燃油泵向發(fā)動機供給的燃油進入燃油調節(jié)器后，在計量活門2油口開度的控制下流經(jīng)出口組件5，發(fā)動機油路為發(fā)動機正常工作提供所需要的燃油。壓差活門1則在任意時刻下保證流入發(fā)動機的燃油壓力恒定，占空比活門3、數(shù)字閥4及電子控制器6則根據(jù)發(fā)動機的工作需要，對計量活門的油口開度進行控制，以實現(xiàn)對燃油供給量的調整。由于燃油調節(jié)器長期工作在弱狀態(tài)下，當其內部零部件出現(xiàn)一定程度磨損(未達到完全失效的狀態(tài))時，其壓力-流量特性與正常狀態(tài)下相差不大。但是，若故障突然發(fā)生，將會造成不可逆轉的后果。而壓差活門1是影響發(fā)動機供油壓力穩(wěn)定的關鍵部件且閥芯磨損是最主要的失效方式。為獲得不同弱狀態(tài)下，燃油調節(jié)器的壓力脈動信號，將壓差活門的閥芯人為磨損至不同程度，并對出口組件5出油口處的壓力脈動進行監(jiān)測。此次測試將圖1中所示的壓差活門1的閥芯進行人工失效，使它與內部腔體的單側配合間隙分別為0.02、0.03、0.04和0.05 mm(完全失效的臨界值)，而單側間隙0.01 mm為設計尺寸且其余液壓元件計量活門2、占空比活門3、數(shù)字閥4、出口組件5、燃油泵及發(fā)動機油路保持初始設計狀態(tài)。測試過程從啟動持續(xù)至平穩(wěn)運行，總時長為20 s，采集數(shù)據(jù)時的采樣步長為0.01 s，每種弱狀態(tài)下的樣本量共計2 000個。出口組件出油口的壓力脈動時域狀態(tài)如圖2所示。

圖1 燃油調節(jié)器資源分配

圖2 不同壓差活門狀態(tài)下的出口組件壓力脈動

由圖2可知：當壓差活門的閥芯與腔體單側配合間隙(側隙)從0.01 mm(設計值)增加至0.05 mm(臨界失效值)時，出口組件出油口的壓力脈動最終將穩(wěn)定在12 MPa左右；當其內部元件出現(xiàn)不同程度的失效(側隙為0.02～0.05 mm)時，壓力脈動變化不明顯，且其曲線分布間隔較為均勻，表明此時的燃油調節(jié)器工作在弱狀態(tài)下，難以利用直接觀察或常規(guī)的時頻域分析方法提取特征以描述其當前的物理狀態(tài)；當壓差活門的閥芯磨損至臨界失效值(0.05 mm)時，出口組件的壓力脈動曲線出現(xiàn)一定程度劇烈波動，線條分布相較于其他情況更為緊密，表明此時的壓差活門已瀕臨失效，流向發(fā)動機供油路的燃油壓力出現(xiàn)劇烈變化。綜上分析，當燃油調節(jié)器內部元件出現(xiàn)磨損但未到達臨界失效值時，其供給發(fā)動機油路的燃油壓力基本處于穩(wěn)定狀態(tài)，但此時的燃油調節(jié)器內部參數(shù)已發(fā)生變化，長此以往會使燃油調節(jié)器由弱狀態(tài)逐步退化至全面失效，若無法在退化形成前識別其物理狀態(tài)并采取相應措施，則會出現(xiàn)航空事故。

1.2 弱狀態(tài)下故障特征提取策略

弱狀態(tài)下燃油調節(jié)器中出口組件的壓力脈動信號存在較強的模糊性，難以直接分辨它所處的物理狀態(tài)。加之數(shù)據(jù)采集過程中硬件條件和時間成本的限制，無法長時間利用物理實驗方法進行高效的數(shù)據(jù)處理及分析工作。故本文作者提出如圖3所示的弱狀態(tài)故障特征提取策略及算法測試流程。

圖3 弱狀態(tài)燃油調節(jié)器狀態(tài)特征提取策略

2 特征提取、重構及K-NN算法設計

2.1 特征提取算法及參數(shù)設計

(1)-means無監(jiān)督聚類算法設計

-means無監(jiān)督聚類算法為利用給定的包含個數(shù)據(jù)的數(shù)據(jù)集尋找數(shù)據(jù)間相似度并構建個簇的方法。該方法將總體數(shù)據(jù)劃分為個內部相似度極高的簇，原則上每個簇至少包含一個數(shù)據(jù)對象，每個數(shù)據(jù)對象必須屬于且只能屬于一個簇。其算法流程如圖4所示。

圖4 k-means算法流程

①初始值的選擇

根據(jù)現(xiàn)有的文獻及-means在模式識別領域的應用研究，誤差平方和所有樣本的聚類誤差可作為值的有效選擇指標，如下所示：

(1)

式中:為第個簇;為中的樣本點;為的質心(中所有樣本的均值);為所有樣本的聚類誤差;為每個簇中的元素。

②平均損失函數(shù)及其閾值

(2)

圖5 不同k值下的fLoss和聚類時間

由圖5(a)可知：隨著的增加，在0.1～0.5的范圍內浮動且隨著的增大，不同弱狀態(tài)聚類結果的均呈現(xiàn)局部上升。進一步由一維搜索原理，可確定各弱狀態(tài)下的最優(yōu)聚類個數(shù)=10。由圖5(b)可知：當=8時，算法完成聚類所需時間最小，平均耗時1.548 s;當=10時，平均耗時2.977 s，但其聚類結果的相較于=8時的普遍表現(xiàn)較優(yōu);而當>10時，平均耗時超過10 s。綜上可知，當前數(shù)據(jù)源下的聚類個數(shù)=10時，平均損失函數(shù)的收斂閾值為0.1。

(2)信息融合算法設計

信息融合又稱多傳感器信息融合，此處的傳感器不僅包括物理意義上的各種傳感器系統(tǒng)，也包括與觀測環(huán)境匹配的各種信息獲取系統(tǒng)。本文作者將提取的時頻域特征視作多個不同傳感器獲取的多源信息，把聚類后各子類的每種時頻域信息都視為來自一個傳感器，通過信息融合進行多時頻域特征變換。利用多種時頻域信息進行原始信號的描述，以有效彌補單特征下信息不完備的問題。

通常情況下，信息融合后新特征具有極高的維度，使得新特征空間的交錯現(xiàn)象更加明顯，增加狀態(tài)監(jiān)測系統(tǒng)使用時的時間復雜度。而主成分分析(Principal Components Analysis，PCA)算法通過線性變換提取數(shù)據(jù)的主要特征分量，將高維數(shù)據(jù)變換為一組各維度線性無關的表示，是提取高維數(shù)據(jù)關鍵信息及降維的最有效方法之一。PCA算法流程如圖6所示，通過PCA獲得的降維數(shù)據(jù)在累計貢獻率的限制下進行篩選。而累積貢獻率的大小表示當前選擇的所有主成分攜帶原數(shù)據(jù)的信息的比例，故文中PCA中的累計貢獻率=0.99。

圖6 PCA算法流程

2.2 空間質心尺度變化方法

利用所提策略對原始信號數(shù)據(jù)進行一次和二次特征提取，獲得的特征數(shù)據(jù)集中，部分弱狀態(tài)的特征數(shù)據(jù)集的特征空間分布呈現(xiàn)較高的位置交錯性。空間質心尺度變化是指在歐氏空間中，將數(shù)據(jù)集中遠離其質心的數(shù)據(jù)點，利用空間坐標變化產(chǎn)生新的數(shù)據(jù)點并使它們聚集在以質心為圓心、為平均半徑的包絡曲面內。通過質心尺度變化算法所產(chǎn)生的新數(shù)據(jù)點，在空間中的分布具有較低的位置交錯性，能有效地提升狀態(tài)監(jiān)測系統(tǒng)的決策精度。原始特征數(shù)據(jù)的空間分布、空間質心尺度變化原理如圖7所示。

圖7 原始特征數(shù)據(jù)分布及質心尺度變化原理

如圖7(a)所示，原始特征數(shù)據(jù)分布在以第一主成分、第二主成分和第三主成分構造的歐氏空間中。記弱狀態(tài)原始特征數(shù)據(jù)集為，此時={∪|≠且,=0.01,…,0.05}，其中和為不同的弱狀態(tài)特征集。圖中可清晰地看出，各弱狀態(tài)下的原始特征數(shù)據(jù)集并未呈現(xiàn)相互獨立的分布，即當≠且,=0.01～0.05時，存在∩≠?，特別是虛線框內所示的交錯現(xiàn)象最為嚴重。當有未知信號的特征集輸入時，狀態(tài)監(jiān)測系統(tǒng)將采取最大相似度原則進行決策，此時會有多個弱狀態(tài)的特征數(shù)據(jù)與未知數(shù)據(jù)存在較高相似性，導致難以獲得準確的決策結果。因此，設計空間質心尺度變化算法，以降低不同弱狀態(tài)下特征數(shù)據(jù)的相似性，使每個相鄰數(shù)據(jù)集間都滿足∩=0，在提升狀態(tài)監(jiān)測系統(tǒng)的決策精度方面具有重要意義。

(3)

2.3 K-NN算法參數(shù)設計

-NN算法基本思路：在訓練集中數(shù)據(jù)和標簽已知的情況下，輸入測試數(shù)據(jù)，將測試數(shù)據(jù)的特征與訓練集中對應的特征進行比較，找到訓練集中與之最為相似的前個數(shù)據(jù)，則將測試數(shù)據(jù)所屬類別對應為個數(shù)據(jù)中出現(xiàn)次數(shù)最多的類別。-NN算法具有易實現(xiàn)、對異常點不敏感、無訓練時間、擅長處理多分類、預測及準確度高的優(yōu)點，廣泛應用于各類型決策系統(tǒng)、推薦系統(tǒng)。其算法描述如下：

輸入：歷史數(shù)據(jù)集、歷史數(shù)據(jù)集對應標簽；

輸出：測試數(shù)據(jù)集所屬標簽。

具體步驟：

(1)計算相似度。計算測試數(shù)據(jù)集每個樣本與所有隸屬數(shù)據(jù)的距離；

(2)形成決策數(shù)據(jù)集。設計超參數(shù)，并以此產(chǎn)生決策數(shù)據(jù)集；

(3)選擇決策法則。將測試數(shù)據(jù)集歸類為個最近鄰樣本中具有最多樣本個數(shù)的類別。

根據(jù)圖2所示數(shù)據(jù)及其數(shù)據(jù)量的特點，將選擇歐氏距離函數(shù)為-NN故障分類器的相似度函數(shù)。同時，為增強-NN算法的健壯性，保證每次算法都能得出結果，將其超參數(shù)設置為奇數(shù)且在決策法則中選用少數(shù)服從多數(shù)的投票準則。

3 算法測試及結果分析

(1)特征提取與重構算法測試分析

利用圖2所示的出口組件各弱狀態(tài)下的壓力脈動數(shù)據(jù)，在Linux環(huán)境下利用Java和Python平臺分別設計燃油調節(jié)器狀態(tài)特征提取算法及數(shù)據(jù)可視化功能。其中，-means算法的聚類個數(shù)=10，損失函數(shù)的閾值為0.1，該算法一次提取的特征包含各子類的均值、標準差、均方根、方根幅值、峭度、偏度、峰值、裕度及波形，共9類時頻域特征。另外，設定PCA信息融合算法的累計貢獻率=0.99，空間質心尺度變化準則的相對尺度參數(shù)分別取為二次提取后的數(shù)據(jù)的標準差、為二次提取后的數(shù)據(jù)的標準差峰值及利用人為調參經(jīng)驗獲取線性尺度=100。運行燃油調節(jié)器狀態(tài)特征提取算法并讀取最終的數(shù)據(jù)可視化結果，如圖8所示。

圖8 fLoss=0.02時，不同尺度下的重構特征分布

由圖8可知:各弱狀態(tài)在標準差尺度下的重構特征數(shù)據(jù)分布的交錯性比峰值尺度變化、線性尺度變化的分布明顯，表明此時的重構特征數(shù)據(jù)質量較差，且標準差不是一個穩(wěn)定的尺度變化參數(shù)；而經(jīng)峰值尺度和線性尺度變化后的重構特征數(shù)據(jù)的分布處于穩(wěn)定的相互獨立狀態(tài)，表明這兩種尺度參數(shù)下的重構特征數(shù)據(jù)分布交錯性較弱，有效地抑制了特征空間的模糊性。進一步分析產(chǎn)生上述現(xiàn)象的機制可知，如式(3)所示的數(shù)學模型將→′的點空間位置變化轉化為以質心為核心，在尺度變化因子限制下的小步長移動。當數(shù)據(jù)點與質心的相對尺度越大時，其尺度變化因子越小，此時質心的移動步長越小，其重構數(shù)據(jù)點′的位置越靠近于質心，更容易實現(xiàn)重構特征在空間中的獨立分布。故圖(d)所示的線性尺度下重構特征數(shù)據(jù)的分布相較于峰值尺度下的重構特征數(shù)據(jù)分布更集中于質心。

綜上可知，參數(shù)的選取將直接影響重構特征。為實現(xiàn)對狀態(tài)特征模糊性的抑制，在使用空間質心尺度變化算法進行特征數(shù)據(jù)重構時，其尺度變化參數(shù)可選擇峰值或線性參數(shù)。進一步從算法感知角度看，用峰值能直接描述出一組數(shù)據(jù)中的最大數(shù)據(jù)且對數(shù)據(jù)自身的完備性和數(shù)據(jù)量大小不具有依賴性，能準確地反映數(shù)據(jù)自身情況，因此峰值尺度下的重構特征能更好地體現(xiàn)出燃油調節(jié)器狀態(tài)特征提取算法在數(shù)據(jù)感知階段所具備的智能性。

(2)-NN分類算法測試分析

利用圖2所示的5組數(shù)據(jù)，在Linux環(huán)境中利用Java按照如表1所示的占比，采用隨機法從各弱狀態(tài)數(shù)據(jù)中抽取相應數(shù)量的數(shù)據(jù)，并在隨機位置拼接成5組不同類別的測試數(shù)據(jù)。

表1 測試數(shù)據(jù)集中各弱狀態(tài)數(shù)據(jù)占比

進一步,利用Java搭建-NN分類算法模型，并分別將未經(jīng)特征提取與重構的測試數(shù)據(jù)集和已進行特征提取與重構(峰值尺度)的測試數(shù)據(jù)集作為-NN分類算法的輸入數(shù)據(jù)，并將其超參數(shù)設定為0.7，即選擇前70%的數(shù)據(jù)作為決策數(shù)據(jù)集。根據(jù)少數(shù)服從多數(shù)的投票法則，-NN分類算法對2種數(shù)據(jù)集決策結果分別如表2、表3所示。

表2 K-NN分類算法對未經(jīng)特征提取與重構的測試數(shù)據(jù)集T的識別率(結果保留小數(shù)點后兩位)

表3 K-NN分類算法對峰值尺度重構的測試數(shù)據(jù)集T的識別率(結果保留小數(shù)點后兩位)

由表2可知：當使用未經(jīng)特征提取與重構的測試數(shù)據(jù)時，5組測試數(shù)據(jù)平均耗時1.251 s，-NN算法基本做出了正確的決策，但其最高識別率均未超過0.5，并且在對第3組數(shù)據(jù)決策時出現(xiàn)模糊決策；而當-NN算法使用經(jīng)過特征提取與重構后測試數(shù)據(jù)時，5組測試數(shù)據(jù)平均耗時0.472 s，-NN算法識別率均高于0.5且其決策結果完全正確。

深入分析-NN分類算法模型可知，由于算法運行時需要搜尋整個數(shù)據(jù)空間，直接將未經(jīng)特征提取與重構的測試數(shù)據(jù)(原始數(shù)據(jù))作為-NN算法的數(shù)據(jù)輸入，會因輸入數(shù)據(jù)量過大，導致運行時間較長，降低識別效率。同時，原始數(shù)據(jù)內含較多的冗余數(shù)據(jù)且數(shù)據(jù)間關聯(lián)性較強，算法產(chǎn)生的決策數(shù)據(jù)集品質不高，難以對各弱狀態(tài)進行準確描述，導致-NN分類算法的決策效率較低且精度較差。而經(jīng)由特征提取與重構(峰值尺度)的數(shù)據(jù)，其數(shù)據(jù)量和數(shù)據(jù)間的關聯(lián)性均得到大幅度削弱，這使得-NN分類算法能夠在保持較少搜尋數(shù)據(jù)量的基礎上，產(chǎn)生具有良好品質的決策數(shù)據(jù)，從而提升-NN分類算法的整體性能。

4 結論

針對弱狀態(tài)下燃油調節(jié)器的狀態(tài)信號存在較強模糊性，難以進行特征提取和故障診斷的問題，在結合無監(jiān)督聚類、多源信息融合與空間質心尺度變化的基礎上，設計特征提取和重構算法。經(jīng)算法測試，驗證了重構后的特征數(shù)據(jù)在特征空間中分布特性良好。進一步利用Java搭建-NN分類算法，將經(jīng)特征提取和重構后的數(shù)據(jù)作為-NN算法的輸入，此時-NN分類算法能夠有效地對測試數(shù)據(jù)進行識別，且在效率與準確率方面均表現(xiàn)良好。綜上可知，本文作者所設計的特征提取與重構方法，不僅能夠有效地提取弱狀態(tài)下燃油調節(jié)器的狀態(tài)特征，還能夠幫助-NN分類算法更加高效且準確地進行故障診斷診斷，為燃油調節(jié)器的后續(xù)健康管理系統(tǒng)搭建提供參考。