蘇興榮,尹聚文
(空軍裝備研究院,北京 100085)
某型GaAs MMIC功率放大器是用于某型裝備的新研大功率芯片,要求在正常結溫150℃的條件下,失效率>10-6/h,可靠壽命≥15年。為考核其可靠性指標是否達到要求,采用恒定溫度應力加速壽命試驗的方法進行驗證。根據(jù)理論分析,失效機理主要為柵下沉[1-2],滿足最弱環(huán)節(jié)壽命特征模型。因此作以下假定。
1)器件在正常應力水平θ0和加速應力水平θi下的壽命都服從威布爾分布,其分布函數(shù)為[3-4]:
式中:m為形狀參數(shù),mi>0;ηi為特征壽命,ηi>0。此假定表明,改變應力水平不會改變壽命分布類型。
2)在θ0和θi下,產品的失效機理不變,由于威布爾分布的形狀參數(shù)反映失效機理,故此假定意味著:m0=mi。
3)產品的特征壽命ηi與所加的溫度應力水平θi間有如下的加速模型,它們之間的關系符合阿倫
尼斯方程,即:ηi=ea+b/θi或 l nηi=a+b/θi。
根據(jù)國內外相關資料和該芯片的生產工藝,確定在結溫225,240,255,270℃等4個加速溫度應力下進行1200 h的定時截尾試驗,每個溫度應力下投入10個樣品,在試驗過程中實時監(jiān)測漏源電流變化量(ΔIDS≤20%),定時測試輸出功率(下降≤1 dB)。在225℃的溫度下無失效,各溫度應力下的失效情況見表1。
表1 各溫度應力下失效情況Table 1 Failure time at each temperature level
首先利用表1的數(shù)據(jù),用Matlab繪制威布爾分布圖[5-6](如圖 1所示)。從圖 1中可以看出,θ2(240℃)和θ3(255℃)溫度應力下2條威布爾分布直線基本平行,θ4(270℃)時直線斜率出現(xiàn)顯著變化??梢猿醪脚袛喔鳒囟葢ο缕骷膲勖耐紶柗植?,但在高溫度應力270℃時失效機理發(fā)生了變化。
進一步利用范-蒙特福特檢驗法檢驗各溫度應力下器件的壽命是否服從威布爾分布[7],這種方法比圖估法更為精確。假設H0:器件的壽命分布為威布爾分布W(m,η),令:
圖1 各溫度應力下的威布爾分布曲線Fig.1 Weibull distribution curve at each temperature level
表2 溫度應力為240℃下范-蒙特福特檢驗計算表Table 2 Calculation table of Van-Montfort test at 240℃
1)240℃時的分布檢驗。溫度應力為240℃時,Gi的值見表2。
表2中的E(Zi)可通過查《可靠性試驗用表》E(Zr,n)得到[8]。把Gi均分為2組,則統(tǒng)計量為:
對于給定的顯著水平α=0.1,如果上述F值滿足以下條件,則認為假設成立:
由于F=0.8386,介于F0.95(2,4)和F0.05(2,4)之間,所以不能拒絕H0,即可以認為該溫度應力下的壽命分布是威布爾分布。
2)255,270℃時的分布檢驗。同理,溫度應力為255℃時,F(xiàn)=3.8525,F(xiàn)0.9(56,6)=0.2336,F(xiàn)0.0(56,6)=4.28;溫度應力為 270 ℃時,F(xiàn)=2.023,F(xiàn)0.95(8,8)=0.290 7,F(xiàn)0.0(58,8)=3.44。因此,在255,270℃溫度應力下的壽命分布都是威布爾分布。
利用阿倫尼烏斯方程估算壽命曲線,需要3個溫度應力下的參數(shù)。在試驗中,雖然270℃下的失效數(shù)據(jù)服從威布爾分布,但其曲線斜率與前2個溫度有明顯變化,故認為失效機理發(fā)生了變化,不能使用該溫度下的失效數(shù)據(jù),而225℃下又無失效。為此,先用經典方法對θ2,θ3參數(shù)進行估計,再使用貝葉斯公式對225℃時的參數(shù)進行估計。
在240,255℃溫度下,用以下公式求得μi和σi的近似無偏估計[1]:
式中:k(ri,ni)和Ec(Zri,ni)的數(shù)值由樣本總量n、失效數(shù)r查近似無偏估計系數(shù)表(極值分布)得到。
溫度為 240 ℃時,k(4,10)=4.496 7,Ec(Z4,10)=-0.966 385;溫度為255 ℃時,k(7,10)=9.135 356,Ec(Z7,10)=-0.013 997。
下面用貝葉斯原理估算θ1(溫度應力225℃)水平下參數(shù)[9-15]。
設有n1個產品進行壽命試驗,截尾時間為τ1,在時刻t1前有r1個產品失效,其中每個產品的失效概率均設為p1,其似然函數(shù)可用二項分布表示:
其中p1為產品在τ1前失效的概率,即:
試驗中,10個產品在時刻1 200 h前無一失效,即r1=0,則似然函數(shù)為:
L(0|p1) =(1-p1)10
現(xiàn)確定參數(shù)p1的先驗分布,p1總界于0與1之間,由于加速壽命試驗產品的失效概率不會很大,因此,p1靠近于0。另外,低溫度應力下的失效概率不可能超過高溫度應力的失效概率,所以p1的一個上限λ可以由240℃下的試驗信息確定:。
若?。?,λ)上的均勻分布作為p1的先驗分布,再聯(lián)系到上述似然函數(shù),由Bayes定理可得p1的后驗分布為:
其中,0<p2<0.3744,利用數(shù)值積分法得 ln的后驗期望為:
利用3種估計的參數(shù),在溫度225,240,255℃等3個應力水平下,加速壽命方程中a和b的計算按表3進行。
按下列公式計算a和b:
1)在正常工作情況下,結溫Tj=150℃時的特征壽命:
2)如果按10年計算,在溫度150℃的正常工作溫度下,平均失效率為:
3)在正常工作溫度為150℃的條件下,平均壽命:
表3 計算加速壽命試驗方程中的系數(shù)a和bTable 3 Calculation of coefficients a and b in the accelerate life test equation
4)產品在某溫度150℃條件下,可靠度為0.9時的可靠壽命為:
從試驗結果及以上分析和計算結果,可以得出以下結論。
1)該GaAs MMIC功率放大器在結溫270℃時的失效機理與正常溫度應力下相比已經發(fā)生改變,不能使用該溫度應力下的失效數(shù)據(jù)對芯片的可靠性指標進行估計。
2)利用器件在低溫度應力下失效概率應小于高溫度應力下的物理經驗,采用均勻分布作為先驗布,利用貝葉斯公式,可估算出低溫度應力下無失效數(shù)據(jù)的威布爾分布參數(shù)。
3)在150℃正常溫度應力下,該功率放大器使用10年的平均失效率小于10-6/h,可靠度0.9時的可靠壽命達到30年以上,滿足使用指標要求。
參考文獻:
[1] 崔曉英,許燕,黃云.GaAs PHEMT器件高溫加速壽命試驗及物理分析[J].電子器件,2010,33(1):22-26.CUI Xiao-ying,XU Yan,HUANG Yun.High Temperature Accelerant Life Test and Physical Analysis of GaAs PHEMT Device[J].Chinese Journal of Electron Device,2010,33(1):22-26.
[2] 王蘊輝,莫郁薇,吳海東,等.砷化鎵微波單片集成電路可靠性預計模型研究[J].電子產品環(huán)境與可靠性試驗,2007,25(2):43-48.WANG Yun-hui,MO Yu-wei,WU Hai-dong,et al.The Research on Reliability Prediction Model for MMIC[J].Electronic Product Reliability and Environmental Testing,2007,25(2):43-48.
[3] 茆詩松,王玲玲.可靠性統(tǒng)計[M].上海:華東師范大學出版社,1984.MAO Shi-song,WANG Ling-ling.The Statistics of Reliability[M].Shanghai:East China Normal University Publishing House,1984.
[4] 劉海濤.氦氖激光器加速壽命試驗理論研究[D].西安:西安電子科技大學,2008.LIU Hai-tao.The Theory Study of He-Ne Lasers Accelerated Life Test[D].Xi'an:Xi'an Electronic Technology University,2008.
[5] GB 2689.1-1981,壽命試驗和加速壽命試驗方法總則[S].GB 2689.1-1981,The General Life Test and AcceleratedLife Test Method[S].
[6] GB 2689.2-1981,壽命試驗和加速壽命試驗的圖估計法(用于威布爾分布)[S].GB 2689.2-1981,The Probability Paper Estimation for Life Test and Accelerated Life Test(for Weibull Distribution)[S].
[7] 茆詩松,王玲玲.加速壽命試驗[M].上海:華東師范大學出版社,1984.MAO Shi-song,WANG Ling-ling.The Acceleration Life Test[M].Shanghai:East China Normal University Publishing House,1984.
[8] 中國電子標準化研究所.可靠性試驗用表(增訂本)[M].北京:國防工業(yè)出版社,1987.The Chinese Electronics Standardization Institute.Tables for Reliability Test(Revised and Enlarged Book)[M].Beijing:National Defense Industry Press,1987.
[9] 茆詩松,李億民,陸淑蘭,等.恒定應力加速壽命試驗中無失效數(shù)據(jù)的處理[J].應用概率統(tǒng)計,1993,9(2):216-218.MAO Shi-song,LI Yi-min,LU Shu-lan,et al.Data Processing for Accelerated Life Test in Constant Stress[J].Chinese Journal of Applied Probability and Statistics,1993,9(2):216-218.
[10] GB 2689.4-1981,壽命試驗和加速壽命試驗的最好線性無偏估計[S].GB 2689.4-1981,The Best Linear Unbiased Estimation for Life Test and Accelerated Life Test[S].
[11] 徐健.定時截尾恒加壽命試驗下無失效數(shù)據(jù)的Bayes處理方法[J].上饒師范學院學報,2007,27(3):1-7.XU Jian.Bayesian Statistical Method for Non-failure Data under Accelerated Life Test Based on Type I Censoring[J].Journal of Shangrao Normal College,2007,27(3):1-7.
[12] 王乃生,王玲玲.恒定應力加速壽命試驗數(shù)據(jù)缺失時的統(tǒng)計分析[J].華東師范大學學報(自然科學版),2002(1):35-44.WANG Nai-sheng,WANG Ling-ling.Statistical Analysis of Constant Stress Accelerated Life Testing under Multiply Type-II Censoring[J].Journal of East China Normal University(Natural Science),2002(1):35-44.
[13] 李億民.恒加壽命試驗中無失效數(shù)據(jù)的處理[J].山東工程學院學報,1996,10(4):55-59.LI Yi-min.A Method for Processing Non-failure Data in Constant Stress Accelerating Life Test[J].Journal of Shandong Institute of ENG,1996,10(4):55-59.
[14] 嚴栓航,師義民.恒應力加速壽命試驗中無失效數(shù)據(jù)的處理[J].純粹數(shù)學與應用數(shù)學,2006,22(2):145-148.YAN Shuan-hang,SHI Yi-min.Zero-failure Data Process of Constant Stress Accelerated Life Testing[J].Pure and Applied Mathematics,2006,22(2):145-148.
[15] 唐彥峰,王森,劉祥凱,等.基于Bayes方法的無失效數(shù)據(jù)壽命分布研究[J].軍事交通學院學報,2011,13(9):58-60.TANG Yan-feng,WANG Sen,LIU Xiang-kai,et al.On Bayesian Method for Life Distribution of Zero-failure Data[J].Journal of Military Transportation University ,2011,13(9):58-60.