陳娟 劉浩 吳建國 李志強 李佳,3 隗靖昆 姜博文 向剛 吳昶霖

（1 北京航空航天大學(xué)，北京 100191；2 北京環(huán)境強度研究所，北京 100076；3 米蘭理工大學(xué)；493160部隊，北京 100076；5 北京航天自動控制研究所，北京 100854；6 宇航智能控制技術(shù)國家級重點實驗室，北京 100854）

0 引言

航天器作為一個國家機械、電子、控制等領(lǐng)域的重要體現(xiàn)，代表著國家的科技實力和綜合國力。航天器的研究、設(shè)計和貯存等環(huán)節(jié)影響著其使用可靠性。一般情況下，航天器沒有過多的實際使用機會，因此，航天器的貯存環(huán)節(jié)占了其生命周期的絕大部分。在航天器長期的貯存環(huán)節(jié)中，存在著諸多影響其貯存和使用可靠性的因素，如振動、溫度、濕度、電應(yīng)力、機械應(yīng)力，等等。航天器各模塊各單機的失效模式和失效機理不盡相同，但一般能夠?qū)⑵錃w結(jié)為器件老化、機械性能或電性能退化等物理、化學(xué)退化過程。同時，針對航天器的部分特定單機而言，其長時間的加電工作狀態(tài)也是健康評估和壽命預(yù)測的重要研究對象?，F(xiàn)階段判斷航天器貯存和使用壽命的主流方法是對實體航天器整機或單機進行加載加速實驗，獲取相關(guān)數(shù)據(jù)進行分析評估，此類方法成本較高，難以獲得大量的數(shù)據(jù)樣本；目前有學(xué)者提出了基于數(shù)字孿生模型的航天器貯存可靠性評估技術(shù)[1]，運用虛擬模型建模仿真技術(shù)，通過原理仿真和故障仿真，能夠使用較低成本獲得大量數(shù)據(jù)，數(shù)據(jù)集能夠有效支撐基于機器學(xué)習(xí)算法的航天器的可靠性評估和健康狀態(tài)評估。

1 研究現(xiàn)狀

1.1 航天器供配電系統(tǒng)

供配電系統(tǒng)肩負著在航天器發(fā)射后為其它系統(tǒng)和有效載荷進行供電的重要任務(wù)，是航天器系統(tǒng)用電設(shè)備的核心，其可靠性直接影響了航天器的工作狀態(tài)及質(zhì)量。航天器的供配電系統(tǒng)主要組成為電源系統(tǒng)（一次電源、二次電源）和電路系統(tǒng)（配電器和電纜網(wǎng)）。航天器的總體電路系統(tǒng)由配電器、電纜網(wǎng)等組成，其中最重要的部件是配電器。配電器的功能是接收單體電池供電，并實現(xiàn)對彈上其他設(shè)備的電壓電流分配，完成對彈體系統(tǒng)各設(shè)備的一次母線供電、集中二次供電以及電壓的變換分配，根據(jù)航天器飛行任務(wù)的要求和地面指令，對航天器負載等進行加載、卸載。配電器具體模塊組成如表1所示。

表1 配電器模塊組成Table.1 The component of power distribution system

1.2 航天器健康狀態(tài)評估

針對航天器系統(tǒng)而言，故障預(yù)測與健康管理（Prognostics and Health Management, PHM）是系統(tǒng)健康評估的關(guān)鍵技術(shù)之一，PHM的原理是利用盡量少的傳感器采集系統(tǒng)運行狀態(tài)參數(shù)，從傳感器返回的數(shù)據(jù)中提取能夠反映航天器系統(tǒng)健康特征的信息，并利用相關(guān)算法將傳感器返回的多維數(shù)據(jù)轉(zhuǎn)化為一維的健康狀態(tài)指數(shù)（或健康置信度）[2]，最終目的是實現(xiàn)航天器系統(tǒng)健康狀況的評價和分級?？紤]到目前以及未來戰(zhàn)爭的復(fù)雜性和嚴酷性，為提高航天器全生命周期的裝備故障及健康狀態(tài)檢測與管理能力，需要開展更為全面的航天器健康狀況認知管理系統(tǒng)建設(shè)。目前較為常用的健康狀態(tài)評估方法可以大致分為基于知識經(jīng)驗的方法、基于模型驅(qū)動的方法以及基于數(shù)據(jù)驅(qū)動的方法?；谥R經(jīng)驗的方法一般是在缺少傳感系統(tǒng)和模型仿真的情況下，通過專家經(jīng)驗進行判斷評估；基于模型驅(qū)動的方法是通過對設(shè)備運行狀態(tài)進行充分把握，通過構(gòu)建物理模型和仿真分析，實現(xiàn)對系統(tǒng)的健康評估；基于數(shù)據(jù)驅(qū)動的方法是采取大量傳感器進行數(shù)據(jù)采集、監(jiān)測，并對系統(tǒng)各個運行參數(shù)進行綜合分析，以達到健康狀態(tài)評估的目的[3][4]。

2 航天器配電系統(tǒng)機理建模

2.1 配電器工作原理

本課題方案中的配電器單機的主要元器件為繼電器（包括線圈電阻、接觸開關(guān)等）、穩(wěn)壓二極管及降壓電阻，其中，繼電器的作用為接收控制信號，控制配電器啟動；穩(wěn)壓管的作用是利用二極管的反向特性使配電器輸出電壓穩(wěn)定在要求范圍內(nèi)；降壓電阻實現(xiàn)270V電源電壓降壓至28V輸出。配電器的工作邏輯如圖1所示。

圖1 配電器工作邏輯Fig.1 The working logic of power distribution system

2.2 配電器功能建模

在對配電器單機功能及相關(guān)參數(shù)分析完成后，使用Simulink建立了配電器的電路模型。模塊化封裝的配電器單機模型如圖2所示。

圖2 配電系統(tǒng)功能Fig.2 Function of power distribution system

本課題將配電器模塊細分為以下子模塊—電源切換模塊、穩(wěn)壓模塊、繼電器模塊以及降壓回路。其中，電源切換模塊由地面電源和箭載電源組成，實現(xiàn)發(fā)射時的電源切換；穩(wěn)壓模塊由穩(wěn)壓二極管和保護電阻構(gòu)成，實現(xiàn)穩(wěn)定電壓270V輸出；繼電器模塊由繼電器線圈和觸點開關(guān)等組成，實現(xiàn)控制回路邏輯判斷、觸點動作等功能；降壓回路由降壓電阻和保護二極管構(gòu)成，實現(xiàn)將270V電壓降壓成三路28V電壓輸出。

2.3 配電器模型仿真結(jié)果

配電器單機仿真模型中的各個模塊有機聯(lián)系、相互配合，實現(xiàn)了仿真要求的功能及輸出，仿真結(jié)果如圖3所示，實現(xiàn)了高壓270V一路輸出（UH）以及低壓28V三路輸出（UL1, UL2 & UL3）。同時保證了輸出回路電流穩(wěn)定在2A。

圖3 配電系統(tǒng)整體模型Fig.3 Virtual model of power distribution system

3 配電系統(tǒng)故障注入及退化仿真

3.1 配電器故障注入

針對配電器單機而言，常見的失效模式主要有以下幾類——單機無輸出、單機接收控制信號但無響應(yīng)、單機響應(yīng)時間超差，等等。除偶發(fā)因素導(dǎo)致的失效外，上述失效模式的主要影響因素是溫濕度等環(huán)境應(yīng)力、電應(yīng)力、沖擊和振動應(yīng)力。下面將以幾種常見失效模式闡述各類應(yīng)力的影響機理[5]。

針對配電器單機無輸出這一失效模式而言，其內(nèi)因主要有以下幾點：繼電器線圈或觸點退化導(dǎo)致開路、高溫導(dǎo)致繼電器動作簧片蠕變斷裂、觸點粘接、銜鐵卡死等，主要影響因素為環(huán)境溫度和電應(yīng)力。在配電器單機貯存和工作過程中，繼電器某些部位受環(huán)境溫度和時間綜合作用發(fā)生退化甚至失效，特別是繼電器線圈、觸點和簧片這類較為“脆弱”的點位，退化或故障模式最終的外在表現(xiàn)為繼電器開路。

針對配電器單機無響應(yīng)這一失效模式而言，其內(nèi)因主要是觸點動作電壓的改變，如觸點間存在雜質(zhì)和汽結(jié)霜、繼電器線圈漆包層老化脫落等，主要影響因素為環(huán)境溫度和環(huán)境濕度。受工作或貯存環(huán)境因素的影響，配電器內(nèi)部會隨時間積累產(chǎn)生雜質(zhì)累積、結(jié)霜等現(xiàn)象，此外，繼電器線圈等部位也會隨工作時間逐漸退化，此類退化失效因素最終會導(dǎo)致配電器單機動作電壓的改變，對外表現(xiàn)為配電器單機接收控制信號但不進行正常響應(yīng)。

針對配電器單機響應(yīng)時間超差這一失效模式而言，其內(nèi)因主要是動作簧片因老化或侵蝕導(dǎo)致的動作不同步、觸點熔焊或冷焊，主要影響因素是環(huán)境應(yīng)力和電應(yīng)力。在繼電器的長期貯存和工作過程中，內(nèi)部線圈的匝數(shù)、阻抗等參數(shù)認為是幾乎不發(fā)生退化的，而在環(huán)境應(yīng)力和時間的綜合影響下，其內(nèi)部的金屬簧片剛度逐漸降低，其自身吸合力減弱，從而導(dǎo)致其動作時間參數(shù)發(fā)生變化，反映在單機層面的失效形式上，形成了單機響應(yīng)時間超差這一故障表現(xiàn)。

除上述幾類典型失效模式外，配電器單機還有許多較為常見的失效表現(xiàn)形式，邵旭東等通過經(jīng)驗分析和理論推導(dǎo)等多種方法研究并總結(jié)了沖擊應(yīng)力對航天產(chǎn)品失效的影響[6]，王鴻鑫研究并分析了振動應(yīng)力對航空航天產(chǎn)品的影響[7]。本文對配電器單機工作狀態(tài)及工作貯存環(huán)境因素進行細致分析和梳理，整理了配電器單機的FMEA表。

潛在的失效模式及后果分析（Failure Mode and Effects Analysis, FMEA）是指對構(gòu)成產(chǎn)品的子系統(tǒng)、零件、構(gòu)成過程的各個工序逐一分析，找出所有潛在的失效模式，并分析其可能的后果，從而預(yù)先采取必要的措施，以提高產(chǎn)品的質(zhì)量和可靠性的一種系統(tǒng)化活動[8]。

將FMEA應(yīng)用在航天器配電器單機上，結(jié)合具體型號及廠商提供的相關(guān)數(shù)據(jù)，本課題對兩類單機的各元器件，如繼電器、二極管、電阻絲、電池等進行了細致分析，梳理完成了配電器和熱電池單機的FMEA表。FMEA表中的故障模式、故障參數(shù)及故障表現(xiàn)等參數(shù)，對應(yīng)了第二章中建立的單機模型，指導(dǎo)了后續(xù)單機故障注入模型的建立。

表2 配電器FMEA表Table2 The FMEA table of power distribution system

結(jié)合梳理完成的配電器FMEA表，將表內(nèi)的故障參數(shù)注入到配電器機理模型中，得到了不同故障參數(shù)對應(yīng)的配電器故障模型，對故障模型進行仿真，能夠得到不同故障參數(shù)對應(yīng)的單機故障輸出特性曲線。下面展示了幾類較為典型的故障輸出曲線。

從圖4可以看到，改變動作開關(guān)的觸點電阻可以影響單機輸出特性，當觸點阻值為無窮大時，配電器單機輸出電壓和輸出電流為0。

圖4 配電器機理模型仿真結(jié)果Fig.4 Output of power distribution system model simulation

根據(jù)圖5可以看出，當繼電器觸點動作電壓升高后，接受同樣的觸點閉合指令，配電器不再正常工作，輸出電壓和輸出電流接近0；根據(jù)圖6可以看到，當設(shè)置觸點動作時間延長這一故障后，配電器單機的響應(yīng)時間也隨之延長。

圖5 注入故障信號（觸點閉合阻值為無窮大）的輸出特性曲線Fig.5 Output of the injected fault signal (contact resistance is infinite)

圖6 注入故障信號（觸點動作電壓升高）的輸出特性曲線 Fig.6 Output of injected fault signal (contact action voltage rise)

圖7 注入故障信號（觸點動作時間延長）的輸出特性曲線Fig.7 Output of injected fault signal (contact action time extension)

3.2 配電器退化仿真

在故障注入模型完成的基礎(chǔ)上，總結(jié)部分元器件級的退化規(guī)律，將退化模式注入模型中，得到了航天器典型單機的在特定條件下的退化規(guī)律。需要指出的是，退化規(guī)律仿真與故障注入仿真的不同之處在于：故障注入仿真的結(jié)果與時間無關(guān)，只是反映了單機模型在某一時刻發(fā)生故障是的工作參數(shù)變化情況；而退化仿真則是加入了強的時間關(guān)聯(lián)，將元器件的退化規(guī)律按時間關(guān)系注入仿真模型中，得到的仿真結(jié)果表征了單機模型在整個工作過程中隨時間的退化關(guān)系。

對配電器中的降壓電阻而言，由于未指定具體型號，本文選用了由典型金屬薄膜電阻的退化曲線擬合而成的對數(shù)型退化公式[9][10]

這一公式給出了確定材料的金屬薄膜電阻在某工作溫度下的電阻值隨工作時間變化的退化規(guī)律。選取一定量的數(shù)據(jù)點注入配電器的機理模型中，可以得到配電器單機的降壓電阻退化輸出特性曲線。通過仿真結(jié)果可以看出，在降壓電阻退化的影響下，配電器輸出特性會呈現(xiàn)變化率逐漸減小的類似對數(shù)型的退化趨勢。

借助本節(jié)獲得的退化仿真數(shù)據(jù)集，對輸出特性退化曲線進行擬合分析，得到對應(yīng)的退化函數(shù)關(guān)系，再設(shè)置合適的評價閾值，可以實現(xiàn)針對不同退化模式影響的單機健康狀態(tài)評估。

3.3 基于BP神經(jīng)網(wǎng)絡(luò)的健康狀態(tài)評估

BP（Back Propagation）神經(jīng)網(wǎng)絡(luò)算法是由Rumelhart在1985年基于誤差反向傳播理論提出的深度學(xué)習(xí)算法，是在多層神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)上構(gòu)建的包含輸入層、隱含層、輸出層的神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，能夠通過誤差反饋學(xué)習(xí)實現(xiàn)最終輸出誤差最小化[11]。在多層神經(jīng)網(wǎng)絡(luò)中，前一層的輸出作為下一層的輸入，同一層的神經(jīng)元之間沒有連接，不同層神經(jīng)元的連接對應(yīng)著各自的權(quán)值，具體結(jié)構(gòu)如圖8(a)所示。從數(shù)學(xué)角度看，神經(jīng)網(wǎng)絡(luò)就是一個函數(shù)f(x)，作用就是將輸入x映射為輸出y，而多層神經(jīng)網(wǎng)絡(luò)則相當于符合函數(shù)y=fn(…f1（f0(x)）…)；從具體應(yīng)用角度看，神經(jīng)網(wǎng)絡(luò)就是一個模型，各層、各節(jié)點及連接為模型內(nèi)部結(jié)構(gòu)，各連接對應(yīng)的權(quán)值就是模型的參數(shù)。在多層神經(jīng)網(wǎng)絡(luò)中加入損失反饋結(jié)構(gòu)就構(gòu)成了BP神經(jīng)網(wǎng)絡(luò)模型，如圖8(b)所示，該模型最大的特點時人為確定的損失函數(shù)會沿各層各節(jié)點的連接進行反向傳播，通過反向傳播算法計算各層的優(yōu)化梯度，結(jié)合梯度下降算法最終實現(xiàn)損失函數(shù)的優(yōu)化，即誤差的最小化。

圖8 配電器輸出電壓在降壓電阻退化影響下的退化曲線Fig.8 Degradation curve of distributor output voltage under the influence of step-down resistance degradation

首先，本文將未進行數(shù)據(jù)集擴充的退化數(shù)據(jù)輸入到BP神經(jīng)網(wǎng)絡(luò)模型中，將大部分數(shù)據(jù)作為訓(xùn)練集，將最后三組數(shù)據(jù)作為檢驗集以檢驗輸出精度，得到的模型預(yù)測結(jié)果如圖9所示。可以看到，在樣本較少的情況下，模型預(yù)測精度較差，無法較好的實現(xiàn)單機健康狀態(tài)評估和壽命預(yù)測功能。樣條插值是一種較為常用的數(shù)據(jù)集擴充方式，其實質(zhì)是用一個近似函數(shù)逼近未知函數(shù)，然后在這個近似函數(shù)曲線上取不同的值實現(xiàn)數(shù)據(jù)集擴充。使用MATLAB庫中自帶的樣條插值函數(shù)能夠?qū)崿F(xiàn)此功能，分別將兩類退化模式的退化數(shù)據(jù)擴充為50組。使用擴充后的訓(xùn)練集進行BP神經(jīng)網(wǎng)絡(luò)算法預(yù)測，將擴充前的原始數(shù)據(jù)作為測試集，預(yù)測的輸出結(jié)果如圖10所示，通過與上圖對比能夠十分明顯地看出預(yù)測精度地改善。根據(jù)BP神經(jīng)網(wǎng)絡(luò)算法預(yù)測的降壓電阻退化模式下的單機壽命為950.40h（以95%的規(guī)定輸出特性為閾值）。

圖9 多層神經(jīng)網(wǎng)絡(luò)的基本結(jié)構(gòu)Fig.9 The structure of Multilayer neural network

圖10 訓(xùn)練集不足情況下的降壓電阻退化BP神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果Fig.10 Prediction results of BP neural network for step-down resistor degradation under insufficient training set

圖11 神經(jīng)網(wǎng)絡(luò)壽命預(yù)測結(jié)果Fig.11 Prediction results of BP neural network for step-down resistor degradation

4 結(jié)論

本文基于航天器實際貯存環(huán)境下故障數(shù)據(jù)不足、壽命預(yù)測和健康評估支持技術(shù)缺少等實際問題，主要進行了以下幾部分內(nèi)容的研究和探索：

1）充分調(diào)研了航天器配電系統(tǒng)的工作原理，在此基礎(chǔ)上運用Simulink仿真平臺建立了單機的機理模型，實現(xiàn)了預(yù)定的工作模式和單機模型仿真輸出。

2）細致總結(jié)和梳理了航天器配電器單機的故障模式，建立了基本覆蓋所有元器件的FMEA表，基于FMEA表和建立的單機機理模型，進行了單機故障注入和故障建模仿真，獲得了單機故障注入模型的仿真輸出特性，總結(jié)出了兩類單機的故障表現(xiàn)形式。

3）基于單機故障模型和典型元器件的退化規(guī)律，進行單機性能退化仿真，獲得了幾類較為典型的單機性能退化數(shù)據(jù)集。在此基礎(chǔ)上，基于BP神經(jīng)網(wǎng)絡(luò)算法完成了配電系統(tǒng)的壽命預(yù)測。

通過對上述內(nèi)容的研究，本文得到了以配電器為主的航天器控制系統(tǒng)核心單機的故障機理和故障模型的具體表現(xiàn)形式，在此基礎(chǔ)上獲得了單機退化數(shù)據(jù)，進一步地，可以借助基于理論退化公式的推導(dǎo)、基于多項式曲線回歸和基于BP神經(jīng)網(wǎng)絡(luò)等多種方式和算法以完成核心單機的健康狀態(tài)評估和壽命預(yù)測。本文的研究內(nèi)容能夠為解決實際貯存環(huán)境下航天器故障數(shù)據(jù)缺乏等問題提供較為有效的解決方案，為航天器健康評估、航天器退化規(guī)律探索、航天器生命周期管理分析及壽命預(yù)測提供數(shù)據(jù)和技術(shù)支持。