邱云峰，班元朗，陳 林

（貴州航天計(jì)量測(cè)試技術(shù)研究所,貴州 貴陽 550025）

引言

航空、航天等高可靠裝備領(lǐng)域要求所選用的元器件具備高性能、高可靠、長(zhǎng)壽命，甚至“零缺陷”。電子元器件壽命與可靠性評(píng)估的傳統(tǒng)方法是基于器件的故障或失效進(jìn)行，以失效時(shí)間、失效模式和故障物理作為統(tǒng)計(jì)分析對(duì)象，通過對(duì)參數(shù)退化、失效數(shù)據(jù)的統(tǒng)計(jì)，推斷得到器件的壽命分布模型（例如威布爾分布），從而進(jìn)行可靠性評(píng)定。但傳統(tǒng)方法是通過加速壽命試驗(yàn)方式獲取大量失效模式和數(shù)據(jù)，對(duì)失效壽命的推算的，而在某些特殊應(yīng)用領(lǐng)域，例如導(dǎo)彈武器系統(tǒng)，使用某一類型元器件較少且器件本身可靠性較高，能夠得到失效數(shù)據(jù)太少或無失效數(shù)據(jù)，尤其是對(duì)高可靠、長(zhǎng)壽命和昂貴的產(chǎn)品，應(yīng)用量少，難以獲得失效數(shù)據(jù)，也不能進(jìn)行完全壽命實(shí)驗(yàn)，甚至不能進(jìn)行壽命實(shí)驗(yàn)，因此難以有效評(píng)估出裝機(jī)元器件的剩余使用壽命（Remaining Useful Life，RUL）[1-2]?；跀?shù)據(jù)驅(qū)動(dòng)的壽命預(yù)測(cè)方法以元器件生成過程數(shù)據(jù)作為基礎(chǔ)，不需要故障演化過程或壽命退化過程的精確解析模型，直接通過大數(shù)據(jù)分析方法對(duì)各類可用數(shù)據(jù)進(jìn)行處理，通過人工神經(jīng)網(wǎng)絡(luò)算法，挖掘?qū)ο笙到y(tǒng)數(shù)據(jù)中隱含的健康狀態(tài)或退化特征信息。此方法能夠在失效數(shù)據(jù)太少或無失效數(shù)據(jù)的情況，以電子元器件部分初態(tài)數(shù)據(jù)為基礎(chǔ)進(jìn)行壽命預(yù)測(cè)，能夠有效解決裝機(jī)元器件壽命評(píng)估難的問題[3]。以雙極型晶體管為對(duì)象，研究常見失效模式與機(jī)理，識(shí)別性能退化敏感參數(shù)，采集雙極型晶體管檢測(cè)試驗(yàn)過程中敏感參數(shù)測(cè)試數(shù)據(jù)，構(gòu)建數(shù)據(jù)集，構(gòu)建特征提取模型對(duì)數(shù)據(jù)集進(jìn)行特征提取，并利用BP 神經(jīng)網(wǎng)絡(luò)算法對(duì)提取到的晶體管壽命特征信息進(jìn)行建模分析[4]，構(gòu)建特征變量與剩余使用壽命之間的映射關(guān)系，從而實(shí)現(xiàn)雙極型晶體管剩余使用壽命的預(yù)測(cè)。

1 晶體管壽命特征參數(shù)提取

雙極晶體管在實(shí)際應(yīng)用過程中受到的環(huán)境應(yīng)力是復(fù)雜的，同時(shí)受到溫度、濕度、電、力、磁等多應(yīng)力的復(fù)合影響，實(shí)際上也是由于這些應(yīng)力的綜合效果影響了產(chǎn)品的壽命，因此在加速試驗(yàn)中引入綜合應(yīng)力的加速模型，可以更精確的模擬實(shí)際環(huán)境條件，從而利用加速試驗(yàn)中得到的數(shù)據(jù)，預(yù)計(jì)正常使用狀態(tài)下的壽命，然而每一種應(yīng)力引起器件失效的模式不一樣，綜合應(yīng)力下故障特性與壽命關(guān)系模型很難確定。雙極型晶體管包含集電極- 發(fā)射極擊穿電壓V（BR）CEO、電流傳輸比hFE等20 多個(gè)性能參數(shù)，每個(gè)參數(shù)又有多種應(yīng)力環(huán)境的數(shù)據(jù)，若將這些數(shù)據(jù)全部作為數(shù)據(jù)集，使用神經(jīng)網(wǎng)絡(luò)來預(yù)測(cè)是難以實(shí)現(xiàn)的。業(yè)內(nèi)常用方法是通過對(duì)雙極型晶體管失效模式和失效機(jī)理進(jìn)行研究，找出幾個(gè)失效敏感參數(shù)或者退化明顯參數(shù)作為壽命預(yù)測(cè)的關(guān)鍵特征參數(shù)。晶體管參數(shù)的變化幾乎都與溫度有關(guān)。PN 結(jié)反向電流是由于勢(shì)壘區(qū)的產(chǎn)生電流及表面漏電流產(chǎn)生的，由于半導(dǎo)體材料中存在雜質(zhì)缺陷和制造過程摻雜工藝的不均勻性，在高溫的作用下，由于一些潛在的工藝缺陷被激活，導(dǎo)致PN 結(jié)反向電流大幅度增加。相比于發(fā)射極，晶體管的集電極摻雜濃度低且表面面積大，更容易受外界應(yīng)力的影響，使得表面漏電增加。隨著晶體管工作時(shí)間積累，漏電流會(huì)逐漸增加，當(dāng)漏電流超出器件正常使用范圍時(shí)，器件失效，達(dá)到壽命終點(diǎn)，所以加速壽命試驗(yàn)中晶體管的ICBO參數(shù)退化明顯，因此常做為表征晶體管壽命特征的主要參數(shù)之一[5]。同時(shí)，由于晶體管Si-SiO2界面氧化層中存在的可動(dòng)離子電荷，既可以在界面的垂直方向，也可在其水平方向移動(dòng)帶電粒子。在溫度的作用下，硅表面感應(yīng)負(fù)電荷的數(shù)量及其位置也隨之發(fā)生變化，從而引起晶體管電流增益hFE參數(shù)的漂移，隨著使用時(shí)間增加Si-SiO2界面氧化層退化，可動(dòng)離子電荷增多，hFE的漂移也就越來越大，hFE參數(shù)漂移和晶體管壽命之間就建立了關(guān)聯(lián)關(guān)系[6]。另外在對(duì)加速壽命試驗(yàn)過程采集的數(shù)據(jù)分析發(fā)現(xiàn)VCE（sat）也有一定的退化趨勢(shì)。因此將ICBO、hFE、VCE（sat）作為壽命預(yù)測(cè)的三個(gè)關(guān)鍵特征參數(shù)。晶體管特性參數(shù)的退化直接反映出了壽命，而晶體管特征參數(shù)眾多，各參數(shù)直接具有一定的相關(guān)性，直接使用全部特征參數(shù)進(jìn)行壽命預(yù)測(cè)實(shí)現(xiàn)難度較大，使用主成分分析方法對(duì)特征參數(shù)分析，提取影響壽命的主成分特征序列，即對(duì)特征參數(shù)降維，降低預(yù)測(cè)模型實(shí)現(xiàn)難度。

2 基于BP 神經(jīng)網(wǎng)絡(luò)的壽命預(yù)測(cè)模型構(gòu)建

以何種形式構(gòu)建預(yù)測(cè)模型模型是影響模型預(yù)測(cè)精度的關(guān)鍵因素。常見的預(yù)測(cè)模型包括多項(xiàng)式回歸模型、移動(dòng)最小二乘模型、Kriging 模型、徑向基函數(shù)模型、支持向量回歸模型、ANFIS 自適應(yīng)模糊神經(jīng)網(wǎng)絡(luò)系統(tǒng)和BP 神經(jīng)網(wǎng)絡(luò)模型等。元器件的敏感性能參數(shù)、壽命和發(fā)生故障概率之間的關(guān)系是一個(gè)非線性時(shí)間序列，對(duì)元器件的壽命預(yù)測(cè)是一個(gè)典型的多元非線性向量回歸問題，BP 神經(jīng)網(wǎng)絡(luò)采用信號(hào)正向傳播和誤差反向傳播算法，通過監(jiān)督學(xué)習(xí)能夠生成一個(gè)具有泛化和容錯(cuò)能力的模型，它結(jié)構(gòu)簡(jiǎn)單，且可操作性強(qiáng)，故本次選用神經(jīng)網(wǎng)絡(luò)的方法進(jìn)行元器件壽命預(yù)測(cè)。本研究以雙極型晶體管為對(duì)象，研究常見失效模式與機(jī)理，識(shí)別性能退化敏感參數(shù)，采集雙極型晶體管加速老化過程中敏感參數(shù)測(cè)試數(shù)據(jù)，構(gòu)建數(shù)據(jù)集，構(gòu)建特征提取模型對(duì)數(shù)據(jù)集進(jìn)行特征提取，并利用BP 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)算法對(duì)提取到的特征信息進(jìn)行時(shí)序建模，構(gòu)建特征變量與剩余使用壽命之間的映射關(guān)系，完成雙極型晶體管剩余使用壽命預(yù)測(cè)，并嘗試將壽命預(yù)測(cè)方法應(yīng)用到航天元器件可靠性評(píng)估中，降低元器件質(zhì)量與可靠性保障成本。

BP 神經(jīng)網(wǎng)絡(luò)是一種具有反向傳播功能的人工神經(jīng)網(wǎng)絡(luò)，是壽命預(yù)測(cè)領(lǐng)域常用的數(shù)學(xué)模型。通常BP 神經(jīng)網(wǎng)絡(luò)由三層以上網(wǎng)絡(luò)構(gòu)成。網(wǎng)絡(luò)由輸入層、輸出層和隱藏層組成，隱藏層往往有多層。網(wǎng)絡(luò)中的每個(gè)神經(jīng)元是一個(gè)節(jié)點(diǎn)，連接權(quán)重系數(shù)連接前層和后層，是網(wǎng)絡(luò)形成一個(gè)整體。BP 神經(jīng)網(wǎng)絡(luò)由前向傳播鏈和反向傳播鏈兩個(gè)部分構(gòu)成。前向傳播過程中，如果輸出達(dá)到了預(yù)期值，則學(xué)習(xí)過程終止；否則，將進(jìn)行反向傳播。反向傳播根據(jù)原始連接路徑，反方向計(jì)量誤差(預(yù)期輸出和實(shí)際輸出之間的差值)。為了減少誤差，可以采用梯度下降法來調(diào)整每個(gè)神經(jīng)元的權(quán)重，從而提供網(wǎng)絡(luò)計(jì)算速度，以使網(wǎng)絡(luò)的推斷值盡可能，更快的近似于目標(biāo)輸出值。圖1 為三層網(wǎng)絡(luò)BP 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)。

圖1 三層BP 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

BP 神經(jīng)網(wǎng)絡(luò)中隱藏層的神經(jīng)元輸出可以通過式（1）、（2）表示

其中yi是隱藏層的輸出，neti是第i 個(gè)節(jié)點(diǎn)的激活值，fH為節(jié)點(diǎn)的激活函數(shù)，BP 神經(jīng)網(wǎng)絡(luò)中常用的激活函數(shù)包括Sigmoid 函數(shù)、tanh 函數(shù)、ReLU 函數(shù)以及Leaky ReLU 函數(shù)等。

基于BP 神經(jīng)網(wǎng)絡(luò)的壽命預(yù)測(cè)模型分為數(shù)據(jù)預(yù)處理（數(shù)據(jù)分析、歸一化等）、網(wǎng)絡(luò)模型構(gòu)建、模型訓(xùn)練及優(yōu)化、模型驗(yàn)證四部分，具體流程見圖2。

圖2 基于BP 神經(jīng)網(wǎng)絡(luò)的壽命預(yù)測(cè)模型原理框圖

2.1 數(shù)據(jù)處理

（1） 首先將收到的不同文檔內(nèi)的電子元器件數(shù)據(jù)整合到同一個(gè)文檔中，將壽命數(shù)據(jù)添加進(jìn)去，刪除無用數(shù)據(jù)。

（2） 對(duì)晶體管的失效機(jī)理進(jìn)行研究，找出與元器件失效關(guān)系較大的敏感參數(shù)，對(duì)這些敏感參數(shù)予以保留，刪除其他參數(shù)，并為相關(guān)數(shù)據(jù)加入數(shù)據(jù)標(biāo)簽。

（3） 將壽命特征參數(shù)數(shù)據(jù)集劃分訓(xùn)練集和測(cè)試集，隨機(jī)劃分?jǐn)?shù)據(jù)集中的80%作為訓(xùn)練數(shù)據(jù)，剩下數(shù)據(jù)作為測(cè)試集。

（4） 為了避免不同參數(shù)數(shù)據(jù)值大小對(duì)模型訓(xùn)練及預(yù)測(cè)結(jié)果的影響，對(duì)數(shù)據(jù)進(jìn)行歸一化。

2.2 模型構(gòu)建

考慮到此壽命預(yù)測(cè)問題為一個(gè)典型的非線性回歸問題，建立了基于BP 神經(jīng)網(wǎng)絡(luò)的非線性回歸模型。根據(jù)數(shù)據(jù)量和數(shù)據(jù)輸出情況，結(jié)合數(shù)據(jù)特點(diǎn)等對(duì)BP神經(jīng)網(wǎng)絡(luò)模型的層數(shù)、隱藏層神經(jīng)元個(gè)數(shù)、激活函數(shù)以及優(yōu)化器等進(jìn)行優(yōu)化設(shè)置。

由于晶體管壽命特征數(shù)據(jù)為非線性數(shù)據(jù)，是一個(gè)隨時(shí)間變化的，時(shí)間序列，而單層神經(jīng)網(wǎng)絡(luò)只能用于表示線性可分離函數(shù)，為了更準(zhǔn)確的實(shí)現(xiàn)壽命預(yù)測(cè)，選擇了三層隱含層作為模型的層數(shù)。在神經(jīng)元個(gè)數(shù)選擇上，根據(jù)BT5551 晶體管數(shù)據(jù)集數(shù)據(jù)量選擇25 個(gè)輸入層神經(jīng)元，15 個(gè)隱藏層神經(jīng)元；激活函數(shù)均為relu 函數(shù)，1 作為輸出層神經(jīng)元個(gè)數(shù)，表示輸出一個(gè)預(yù)測(cè)數(shù)據(jù)。

在進(jìn)行模型優(yōu)化時(shí)，選擇RMSprop 模型優(yōu)化器、學(xué)習(xí)率設(shè)為0.001，模型損失函數(shù)為均方誤差mse，評(píng)測(cè)標(biāo)準(zhǔn)為均方誤差mse 和平均絕對(duì)值誤差mae。

2.3 模型訓(xùn)練及優(yōu)化

在數(shù)據(jù)訓(xùn)練過程中，通過每完成一次訓(xùn)練打印一個(gè)點(diǎn)來顯示訓(xùn)練進(jìn)度。選擇20 000次作為模型的最大數(shù)據(jù)擬合次數(shù)，按20%比例從訓(xùn)練集中取出一部分?jǐn)?shù)據(jù)作為驗(yàn)證、查看模型效果。為了防止過擬合，定義回調(diào)函數(shù)callbacks 的操作，每經(jīng)過10 次訓(xùn)練，會(huì)調(diào)用該函數(shù)進(jìn)行一次判斷，判斷依據(jù)是當(dāng)驗(yàn)證數(shù)據(jù)的val_loss 不再下降時(shí)停止訓(xùn)練。訓(xùn)練過程中會(huì)返回一個(gè)history 對(duì)象，儲(chǔ)存訓(xùn)練過程中的loss、mae、mse 信 息，后 續(xù) 通 過 查 看history 中 的mae、val_mae 和mse、val_mse 數(shù)據(jù)和測(cè)試集的預(yù)測(cè)數(shù)據(jù)來觀察訓(xùn)練效果。如果訓(xùn)練效果較差，則返回調(diào)整輸入層或者隱藏層神經(jīng)元個(gè)數(shù)或者增加隱藏層層數(shù)，直至訓(xùn)練效果達(dá)到相應(yīng)標(biāo)準(zhǔn)。

2.4 模型驗(yàn)證

利用一開始劃分的測(cè)試集數(shù)據(jù)對(duì)訓(xùn)練好的模型進(jìn)行驗(yàn)證、獲得預(yù)測(cè)數(shù)據(jù)。將預(yù)測(cè)數(shù)據(jù)與測(cè)試數(shù)據(jù)進(jìn)行比較，查看誤差分布及大小。

3 BT5551 晶體管壽命預(yù)測(cè)結(jié)果

選取BT5551 晶體管壽命試驗(yàn)后測(cè)試數(shù)據(jù)共計(jì)2 098 條，根據(jù)前文分析漏電流ICEO 是影響晶體管壽命主要因素之一，因此選取ICEO 值異常數(shù)據(jù)20 條、正常數(shù)據(jù)30 條，共計(jì)50 條數(shù)據(jù)進(jìn)行壽命的仿真和計(jì)算，隨機(jī)選取其中40 條數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)集訓(xùn)練模型，剩下10 條作為測(cè)試數(shù)據(jù)集。測(cè)試集理論值數(shù)據(jù)和預(yù)測(cè)結(jié)果數(shù)據(jù)見圖3。

圖3 BT5551 晶體管壽命預(yù)測(cè)結(jié)果

根據(jù)式（3）對(duì)預(yù)測(cè)結(jié)果的置信度進(jìn)行分析。對(duì)于BT5551 系列晶體管，選出的數(shù)據(jù)量一共為50 條，壽命理論值的平均值μ=202 428.9，標(biāo)準(zhǔn)誤差σ=17 619.321 34。0.8 置信度下的置信區(qū)間為（201 794.18，203 063.61）。壽命預(yù)測(cè)值的平均值μ=207 185.30，標(biāo)準(zhǔn)誤差σ=11 007.08。0.8 置信度下的置信區(qū)間（206 788.79, 207 581.82）。預(yù)測(cè)值的置信區(qū)間上限誤差為（207 581.82-203 063.61）/203 063.61 ≈2.225%，預(yù)測(cè)值的置信區(qū)間下限誤差為（206 788.79-201 794.18）/201 794.18≈2.475%。

4 結(jié)論

信息系統(tǒng)裝備元器件大多具有“壽命較長(zhǎng)、貯存時(shí)間較久”的特點(diǎn)，在正常環(huán)境影響應(yīng)力作用下，元器件一般很難出現(xiàn)故障，但是對(duì)于處于長(zhǎng)時(shí)間惡劣環(huán)境中工作的元器件，卻容易受到溫度、濕度、電、磁等環(huán)境應(yīng)力的影響，造成其性能下降，甚至失效。所以這些元器件在投入使用前，為了驗(yàn)證其可靠性，會(huì)進(jìn)行大量的可靠性測(cè)試或者加速壽命試驗(yàn)，然而這些測(cè)試或試驗(yàn)往往會(huì)消耗大量的人力和物力。本研究通過對(duì)元器件質(zhì)量數(shù)據(jù)的分析，構(gòu)建了基于深度學(xué)習(xí)的神經(jīng)網(wǎng)絡(luò)模型，運(yùn)用所建立的BP 神經(jīng)網(wǎng)絡(luò)模型對(duì)晶體管進(jìn)行了剩余壽命預(yù)測(cè)。選取漏電流ICEO 作為壽命敏感參數(shù)，選取50 只BT5551 晶體管可靠性壽命試驗(yàn)后的測(cè)試數(shù)據(jù)作為BP 神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練數(shù)據(jù)集和測(cè)試數(shù)據(jù)集。對(duì)壽命預(yù)測(cè)模型的輸出結(jié)果進(jìn)行分析，最后的預(yù)測(cè)結(jié)果誤差＜3%，說明用BP 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)晶體管的剩余壽命是可行的，也驗(yàn)證了所建立的神經(jīng)網(wǎng)絡(luò)模型的準(zhǔn)確性和良好的泛化能力。