牛躍聽，穆希輝，姜志保，楊振海

（1. 總裝備部 軍械技術(shù)研究所，石家莊 050000；2. 北京工業(yè)大學(xué)，北京 100022）

自然貯存環(huán)境下某型加速度計(jì)貯存壽命評(píng)估

牛躍聽1，穆希輝1，姜志保1，楊振海2

（1. 總裝備部 軍械技術(shù)研究所，石家莊 050000；2. 北京工業(yè)大學(xué)，北京 100022）

某型加速度計(jì)貯存于我國亞濕熱、亞干熱、溫和、干燥4個(gè)典型氣候區(qū)。針對(duì)歷年統(tǒng)計(jì)的加速度計(jì)的成敗型、不完全故障數(shù)據(jù)，假設(shè)其服從指數(shù)分布族、Weibull分布族、極值分布族、對(duì)數(shù)正態(tài)分布族，結(jié)合工程實(shí)際處理異常數(shù)據(jù)、“倒掛”數(shù)據(jù)，采用極小χ2估計(jì)對(duì)分布函數(shù)的參數(shù)進(jìn)行估計(jì)，采用極小χ2檢驗(yàn)對(duì)各分布函數(shù)的合理性進(jìn)行了驗(yàn)證，計(jì)算“服從不同自由度下的χ2分布隨機(jī)變量”檢驗(yàn)的擬合優(yōu)度，得到了加速度計(jì)貯存可靠性分布函數(shù)。通過統(tǒng)計(jì)的故障數(shù)據(jù)，驗(yàn)證了所提出的數(shù)學(xué)模型評(píng)估方法的適用性和正確性，得到了置信度0.90、可靠度0.95條件下加速度計(jì)自然貯存壽命。

加速度計(jì)；自然貯存；壽命評(píng)估

加速度計(jì)是為武器慣導(dǎo)平臺(tái)提供運(yùn)動(dòng)載體加速度信息的關(guān)鍵部件，其顯著特點(diǎn)是長期貯存，一次使用。圍繞武器系統(tǒng)關(guān)鍵部組件貯存壽命的評(píng)估問題是近年來研究熱點(diǎn)之一。如果沒有長期貯存期間的加速度計(jì)檢測數(shù)據(jù)，通常在一些假設(shè)前提條件下，采用加速壽命試驗(yàn)的方法，預(yù)估其壽命[1-5]。

某型加速度計(jì)貯存于我國亞濕熱、亞干熱、溫和、干燥4個(gè)典型氣候區(qū)，至今已經(jīng)自然貯存接近10年時(shí)間，期間積累了大量的檢測數(shù)據(jù)。自然貯存使加速度計(jì)受到各種環(huán)境因素的綜合作用，可以真實(shí)、直觀地反映其在多環(huán)境因素作用下的性能變化規(guī)律[6]。那么，加速度計(jì)的貯存壽命到底有多長？這是生產(chǎn)和使用單位都十分關(guān)注的問題。本文旨在評(píng)估加速度計(jì)可靠壽命指標(biāo)，為決策部門提供加速度計(jì)訂購生產(chǎn)、儲(chǔ)備布局、使用維護(hù)決策的依據(jù)。

1 自然貯存加速度計(jì)故障統(tǒng)計(jì)

分別抽取自然環(huán)境中貯存9年、8年、7年、6年、5年、4年、3年的加速度計(jì)樣本進(jìn)行性能檢測，檢測結(jié)果如表1。

在貯存期間加速度計(jì)每2年檢測一次，檢測數(shù)據(jù)為成敗型數(shù)據(jù)，并且不知道加速度計(jì)在檢測周期內(nèi)出現(xiàn)故障的具體時(shí)間。因此，加速度計(jì)抽樣檢測數(shù)據(jù)屬于成敗型不完全數(shù)據(jù)。

表1 加速度計(jì)故障統(tǒng)計(jì)表Tab.1 Accelerometer fault statistics

根據(jù)加速度計(jì)故障統(tǒng)計(jì)結(jié)果，結(jié)合工程實(shí)際，運(yùn)用數(shù)理統(tǒng)計(jì)的方法評(píng)估其自然貯存壽命。

2 貯存可靠性分布函數(shù)類型選擇及參數(shù)確定

為了合理地評(píng)估加速度計(jì)在自然貯存環(huán)境下的貯存壽命，首先需要研究其貯存可靠性的變化規(guī)律，確定貯存可靠性分布函數(shù)。

下面給出問題的統(tǒng)計(jì)模型：

設(shè)產(chǎn)品的貯存壽命為X，X的分布函數(shù)為F(X,θ)。在時(shí)刻ti對(duì)X作了ni次觀察，其中有Yi次觀察值大于或等于ti，ni-Yi次小于ti，但未觀察到X的準(zhǔn)確值。將X≥ti視為成功，X＜ti視為失敗，對(duì)連續(xù)變量X，觀察數(shù)據(jù)(ni, Yi, ti)是成敗型不完全數(shù)據(jù)，易見：

貯存可靠性分布函數(shù)的確定歸結(jié)為以下兩個(gè)基本問題：

2.1 分布函數(shù)類型的選擇

即檢驗(yàn)假設(shè)H0: F∈p0，p0={F( · ;θ)，θ∈Θ}是分布族，Θ是參數(shù)空間。對(duì)于貯存可靠性，常見的分布族有指數(shù)分布族、威布爾分布族、極值分布族、對(duì)數(shù)正態(tài)分布族等，這一問題稱為模型檢驗(yàn)。

指數(shù)分布族：

威布爾分布族：

極值分布族：

對(duì)數(shù)正態(tài)分布族：

2.2 分布函數(shù)中參數(shù)的確定

在F∈p0的假設(shè)下，基于成敗型不完全數(shù)據(jù)(ni, Yi, ti)，1≤i≤m，估計(jì) F(X, θ)中的未知參數(shù)θ。

記時(shí)刻ti的第j個(gè)樣品的貯存壽命為X（iji=1,2,…; j=1,2,…，Xij)獨(dú)立同分布，其共同的分布函數(shù)為F( · ;θ)，θ=(θ1,θ2,…,θs)′∈Θ?RS，θ為參數(shù)，Θ為參數(shù)空間。令：

則：Yi～B(ni, Pi)，Pi= Pi(θ)=1-F(ti, θ)，1≤i≤m。

1) 構(gòu)造χ2統(tǒng)計(jì)量

假定貯存壽命X的分布函數(shù)F(t, θ)為某一分布族，基于數(shù)據(jù)(ni, Yi, ti)，1≤i≤m，可以得到在ti時(shí)刻：

將以上諸式代入式(1)，即可算得χ2統(tǒng)計(jì)量。

2) 分布函數(shù)中參數(shù)的估計(jì)

易見，若真實(shí)參數(shù)為θ0，當(dāng)ni→ ∞，1≤i≤m，χ2(θ0)漸近于自由度為m的χ2分布。若滿足：則稱是θ 的極小χ2估計(jì)。

② 當(dāng)n→∞時(shí)，

③ χ2()的極限分布（n→∞）為自由度為m-s的分布，其中m是觀測樣本的次數(shù)，s是分布函數(shù)中未知參數(shù)的個(gè)數(shù)。

利用上述定理就可以選擇較為合適的貯存可靠度分布類型，并對(duì)其中的未知參數(shù)進(jìn)行估計(jì)。

4) χ2統(tǒng)計(jì)量的最小值χ2()及檢驗(yàn)

對(duì)上述χ2統(tǒng)計(jì)量關(guān)于θ求最小值，得到極小χ2估計(jì)和極小χ2統(tǒng)計(jì)量χ2()。利用定理的第三條結(jié)論進(jìn)行極小χ2檢驗(yàn)，對(duì)于給定的檢驗(yàn)水平α，確定拒絕域?yàn)椋?/p>

為明確檢驗(yàn)的置信度，可以在作出拒絕或接受原假設(shè)的結(jié)論的基礎(chǔ)上，進(jìn)一步計(jì)算檢驗(yàn)的p值。記服從自由度為m-s的分布的χ2隨機(jī)變量為χ2(m-s)，則：P=P(χ2(m-s)＞χ2())，p為檢驗(yàn)的擬合優(yōu)度，衡量了假定的分布與數(shù)據(jù)的擬合程度。一般說來，p值越大，擬合程度越好。

對(duì)于貯存可靠性中的四種常見分布族：指數(shù)分布族、威布爾分布族、極值分布族和對(duì)數(shù)正態(tài)分布族，可以逐一按上述方法進(jìn)行檢驗(yàn)，從中選擇一種最為合適的分布類型。

3 貯存壽命評(píng)估的理論與方法

記R(t)為某一貯存對(duì)象的貯存可靠度，RL為給定的可靠度下限，1-α為置信度，欲求貯存年限T，使得P(R(t)≥RL)=1-α，易見：

其中，N(0,1)是均值為0，方差為1的標(biāo)準(zhǔn)正態(tài)分布。

由于其滿足正態(tài)分布，則可以推出：

再通過求解不等式可以得到R(t)的下限：

式中，Uα為標(biāo)準(zhǔn)正態(tài)分布的 α的上分位數(shù)。上面的式子中假定ni全相等。但在本問題中，由于ni不全相等，為了得到較為保守的R(t)下限，可用下式求得：

從式(2)中解出R(t)，即可得到所要的T。

4 故障統(tǒng)計(jì)數(shù)據(jù)預(yù)處理

4.1 異常的統(tǒng)計(jì)數(shù)據(jù)

分析表1發(fā)現(xiàn)：由于抽樣誤差、檢測誤差、加速度計(jì)批次質(zhì)量等方面的原因，現(xiàn)有的性能檢測數(shù)據(jù)存在一些問題，主要表現(xiàn)在：

圖1 加速度計(jì)可靠度分布Fig.1 Accelerometer reliability distribution

為了確定合理的貯存可靠度分布，進(jìn)而對(duì)貯存壽命進(jìn)行比較可靠的評(píng)估，在計(jì)算過程中我們針對(duì)不同情況對(duì)所需的檢測數(shù)據(jù)進(jìn)行了必要的工程處理。

4.2 異常的統(tǒng)計(jì)數(shù)據(jù)的處理

1）剔除某些異常值。

對(duì)于加速度計(jì)，個(gè)別年份的故障數(shù)過高或過低會(huì)給分布的選擇、確定及壽命預(yù)測帶來嚴(yán)重影響。對(duì)明顯的異常數(shù)據(jù)，能確定產(chǎn)生原因的應(yīng)予以剔除[7]。

由表1可知：貯存6年的加速度計(jì)故障數(shù)量為17。通過對(duì)貯存 6年的不同批次加速度計(jì)故障進(jìn)行統(tǒng)計(jì)（表2）。發(fā)現(xiàn)第g、i批次的部件A故障數(shù)量明顯高于其他批次，說明該批次質(zhì)量存在問題，處理數(shù)據(jù)時(shí)應(yīng)剔除該故障數(shù)據(jù)。同時(shí)，也要在貯存6年檢測樣本總量中去除這兩個(gè)批次數(shù)量（300個(gè)）。

表2 貯存6年的不同批次部件A故障統(tǒng)計(jì)表Tab.2 Fault statistics of part A in different batch after storage for 6 years

表3 剔除異常值后的加速度計(jì)部組件故障統(tǒng)計(jì)表Tab.3 Accelerometer fault statistics after eliminating outliers

剔除異常值后的統(tǒng)計(jì)表如表3所示。

2）盡量不采用數(shù)據(jù)倒掛年份數(shù)據(jù)。

3）基于連續(xù)年份的數(shù)據(jù)進(jìn)行計(jì)算。

在計(jì)算過程中還遵循了這樣的原則：盡量使用連續(xù)幾個(gè)年份的數(shù)據(jù)，以減少由于外推產(chǎn)生的誤差。

5 計(jì)算結(jié)果

各假設(shè)分布類型、假設(shè)檢驗(yàn)、擬合優(yōu)度、貯存壽命評(píng)估結(jié)果如表4所示。比較表4計(jì)算結(jié)果可知：極值分布的χ2統(tǒng)計(jì)量中，0.2569最小，p值最高，并且極小χ2檢驗(yàn)結(jié)果為接收，所以判定加速度計(jì)符合極值分布。取置信度 1-α=0.90，RL=0.95，利用式(1)算得加速度計(jì)的貯存壽命為12.94年。

表4 加速度計(jì)備選分布類型相關(guān)計(jì)算情況Tab.4 Related calculation situations of accelerometer’s alternative distribution types

6 驗(yàn)證與評(píng)價(jià)

我們從以下幾個(gè)方面進(jìn)行了分析和考察，以驗(yàn)證和評(píng)價(jià)計(jì)算結(jié)果的合理性。

6.1 擬合優(yōu)度評(píng)價(jià)

在應(yīng)用極小χ2檢驗(yàn)選擇貯存可靠度分布類型時(shí)，計(jì)算了檢驗(yàn)的p值，即擬合優(yōu)度[8-9]，其p值在0.3065至0.880之間，擬合得比較好[10]。

6.2 可靠度分布函數(shù)計(jì)算與年度統(tǒng)計(jì)數(shù)據(jù)可靠度比較

利用確定的分布函數(shù)對(duì)各年份的可靠度進(jìn)行了預(yù)測，與檢測數(shù)據(jù)得到的可靠度估計(jì)進(jìn)行了比較。

以貯存3年的加速度計(jì)為例，其年度統(tǒng)計(jì)數(shù)據(jù)可靠度為：

加速度計(jì)服從極值分布，由其可靠度分布函數(shù)計(jì)算得到貯存3年時(shí)的可靠度為：

極值分布族：

兩者的相對(duì)誤差為：

同理，可求得加速度計(jì)的年度統(tǒng)計(jì)數(shù)據(jù)可靠度、可靠度分布函數(shù)年度可靠度，以及兩者相對(duì)誤差，見表5?？梢妰烧呦鄬?duì)誤差較小，預(yù)測精度較高。

表5 加速度計(jì)各年份可靠度預(yù)測值與估計(jì)值的比較Tab.5 Comparison of accelerometer reliability between the predicted values and the estimated value each year

（References）：

[1] 袁宏杰，李樓德，段剛，吳浩. 加速度計(jì)貯存壽命與可靠性的步進(jìn)應(yīng)力加速退化試驗(yàn)評(píng)估方法[J]. 中國慣性技術(shù)學(xué)報(bào)，2012，20(1)：113-116.

YUAN Hong-jie, LI Lou-de, DUAN Gang, WU Hao. Storage life and reliability evaluation of accelerometer by step stress accelerated degradation testing[J]. Journal of Chinese Inertial Technology, 2012, 20(1): 113-116.

[2] Lee Jinsuk, Pan Rong. Analyzing step-stress accelerated life testing data using generalized linear models[J]. IIE Transactions, 2010, 42(8): 589-598.

[3] Bae S J, Kim S-J, Park J I, et al. Lifetime prediction through accelerated degradation testing of membrane electrode assemblies in direct methanol fuel cells[J]. International Journal of Hydrogen Energy, 2010, 35: 9166-9167.

[4] Wohlgemuth J H, Kurtz S. Using accelerated testing to predict module reliability[C]//IEEE Photovoltaic Specialists Conference, 2011: 3601-3605.

[5] LI Xiao-yang, JIANG Tong-min. Optimal design for stepstress accelerated degradation with competing failure modes [C]//Annual Reliability and Maintainability Symposium, 2009: 64-68.

[6] 楊振海，安保社. 基于成敗型不完全數(shù)據(jù)的參數(shù)估計(jì)[J].應(yīng)用概率統(tǒng)計(jì)，1994, 10(2): 142-147. YANG Zhen-hai, AN Bo-she. Parameter Estimation Based on Binary Date[J]. Chinese Journal of Applied Probabi1ity and Statistics, 1994, 10(2): 142-147.

[7] 李東陽. 彈藥儲(chǔ)存可靠性分析設(shè)計(jì)與試驗(yàn)評(píng)估[M].北京：國防工業(yè)出版社，2013. LI Dong-yang. Ammunition storage reliability analysis design and test assessment[J]. Beijing: Defense Industry Press, 2013.

[8] Su Y, Shi H F. Goodness-of-fit analysis for uniformity in a bounded solid region[C]//Proceeding of the Ninth International Conference on Machine Learning and Cybernetics, 2010(1): 231-235.

[9] ZHAG Feng, GAO Yuan-Ning, HUO Lei. A viable method for goodness-of-fit test in maximum likeihood fit[J]. Chinese Physics C, 2011, 35(6): 580-584.

[10] 楊振海，程維虎，張軍艦. 擬合優(yōu)度檢驗(yàn)[M]. 北京：科學(xué)出版社，2011.

Storage life assessment of an accelerometer under natural storage environment

NIU Yue-ting1, MU Xi-hui1, JIANG Zhi-bao1, YANG Zhen-hai2
(1. General Armament Ordnance Technology Institute, Shijiazhuang 050000, China; 2. Beijing University of Technology, Beijing 100022, China)

In our country, a certain type of accelerometer is stored in four typical climate zones which are sub humid, Asia hot, mild and dry zones. In view of the historical statistics (accelerometer's success or failure, and incomplete failure data), assuming that they are subject to exponential distribution, Weibull distribution, extreme value distribution and logarithmic normal distribution, and by combining with the engineering data to handle the exceptions and "upside down" data, a minimum χ2estimation method is used to estimate the parameters of the distribution function, and the fit goodness of “obey distributed random variables under different degrees of freedom χ2” test is calculated. In this way the storage reliability distribution function of an accelerometer is obtained. The fault data statistics verify the applicability and validity of the proposed mathematical model assessment method, and the accelerometer natural storage life is obtained under the condition of confidence 0.90 and reliability 0.95.

accelerometer; natural storage; life evaluation

1005-6734(2014)04-0552-05

10.13695/j.cnki.12-1222/o3.2014.04.024

TJ760.6

A

2014-04-18；

2014-07-28

中國博士后科學(xué)基金資助項(xiàng)目（2013M532181）

牛躍聽（1978—），男，博士后，主要研究方向?yàn)樾畔⒒瘡椝帀勖u(píng)估與延壽技術(shù)。E-mail：nyt369@sina.com