【摘要】文章從產(chǎn)品壽命是衡量質(zhì)量的重要指標入手，通過三種方法來闡述產(chǎn)品壽命試驗的損傷失效率的數(shù)學模型.在分析過程中，將產(chǎn)品壽命假定為服從幾何分布，借助相應(yīng)的概率論與數(shù)理統(tǒng)計知識，將在步加試驗中，對服從幾何分布的產(chǎn)品壽命的損傷失效率數(shù)學模型中的相關(guān)參數(shù)做出極大似然估計，最后通過一個簡單的例題加以印證.

【關(guān)鍵詞】損傷失效率；數(shù)學模型；步加試驗；損傷因子；極大似然估計

一、產(chǎn)品壽命試驗的損傷失效率數(shù)學模型

隨著我國經(jīng)濟融入全球經(jīng)濟一體化進程，產(chǎn)品質(zhì)量成為衡量商品進出口的重要指標.為提高我國產(chǎn)品與服務(wù)質(zhì)量的總體水平，推動中國制造走向世界，為我國經(jīng)濟全面、協(xié)調(diào)、可持續(xù)發(fā)展，我國產(chǎn)品質(zhì)量協(xié)會始終不渝地堅持服務(wù)政府、服務(wù)社會、服務(wù)企業(yè)和廣大消費者（用戶）的宗旨.而壽命是衡量產(chǎn)品質(zhì)量好壞的一個重要指標.

產(chǎn)品進行壽命試驗，就目前情況來看，一般有三種不同的試驗方法：第一種稱恒定應(yīng)加速壽命試驗（簡稱恒加試驗），指的是產(chǎn)品在進行試驗的過程中應(yīng)力保持不變；第二種稱為步進應(yīng)力加速壽命試驗（簡稱步加試驗），指的是產(chǎn)品在試驗過程中，應(yīng)力呈階梯狀上升；第三種稱為序進應(yīng)力加速壽命試驗（簡稱序加試驗），指的是產(chǎn)品在進行試驗過程中，應(yīng)力則是連續(xù)上升的（通常是指線性上升），它是一種加于受試產(chǎn)品上的應(yīng)力隨時間連續(xù)增加的一種壽命試驗，試驗一直持續(xù)到某一固定時間或受試產(chǎn)品有部分或全部失效為止.對于恒加試驗，當應(yīng)力水平較低時，試驗需要的時間較長.而步加試驗和序加試驗卻能縮短試驗時間，節(jié)省大量人力、物力和財力.本文以步加試驗產(chǎn)品壽命來加以研究.

在步加試驗中，設(shè)S1

損傷失效率數(shù)學模型，考慮的是一批產(chǎn)品在步進應(yīng)力加速試驗下，從時刻t=0開始直到一個固定的時刻t1都遭受到一個應(yīng)力S1，在時刻t1未失效的產(chǎn)品受到應(yīng)力S2(>S1)，且試驗到產(chǎn)品都失效為止.假設(shè)這種應(yīng)力變化的結(jié)果是導(dǎo)致開始時失效率函數(shù)λ1(y)乘上與變化點t1有關(guān)的一個未知因子α(>1).記步進應(yīng)力壽命時間Y*的失效函數(shù)為λ*(y)，所提議的損傷失效率數(shù)學模型為：

λ*(y)=λ1(y)，當y≤t1時，αλ1(y)，當y>t1時.

因子α將由S1和S2確定，而且有可能和時間t1也有關(guān).于是α一般記為α(t1)，在此稱為損傷因子.

二、幾何分布產(chǎn)品壽命試驗損傷失效率數(shù)學模型下的統(tǒng)計分析

設(shè)應(yīng)力S1下產(chǎn)品壽命服從參數(shù)為p的幾何分布，即P(X=k)=pqk-1，k=1，2，3，….

在應(yīng)力S1下將n個產(chǎn)品投入試驗，試驗進行到第k0次后將應(yīng)力提高到S2(>S1)繼續(xù)做試驗，試驗持續(xù)到所有產(chǎn)品均失效為止.在應(yīng)力S1下有r個產(chǎn)品失效，其次序失效時間為X(1)≤X(2)≤…≤X(r).

在應(yīng)力S2下有n-r個產(chǎn)品失效，其次序失效“時間”（從0算起）X(r+1)≤X(r+2)≤…≤X(n)，而X(r)≤k0

在應(yīng)力S1下產(chǎn)品的失效率

λ1(k)=P(X=k)P(X≥k)=pqk-1∑∞l=kpql-1=qk-1qk-11-q=p.

當在k0時應(yīng)力從S1提高到S2，此時，假定產(chǎn)品失效率服從損傷籌集資金率（TFR）模型，也即失效率λ(k)=αλ1(k)=αp，k>k0.

而其中損傷因子α>1，其值將由S1和S2確定，而且有可能和“時間k0”也有關(guān)，在此α也可記作α(k0).

當k≥k0+1時，P(X=k)=αpP(X≥k)=αpP(X=k)+αpP(X≥k+1)，即(1-αp)P(X=k)=αpP(X≥k+1).

也即P(X=k)=αp1-αpP(X≥k+1)

=αp1-αp［P(X=k+1)+P(X≥k+2)］

=αp1-αpαp1-αpP(X≥k+2)+P(X≥k+2)

=αp(1-αp)2P(X≥k+2).



由此，一般地，對K≥k≥k0+1，有

P(X=k)=αp(1-αp)K-kP(X≥K).

從而有P(X=k)=αp(1-αp)kP(X≥K)(1-αp)K，K≥k≥k0+1.

令θ=limK→∞P(X≥K)(1-αp)K（θ應(yīng)為α，p和k0的函數(shù)）.

又由于P(X≥k0)=1-P(X≤k0)=1-∑k0i=1pqi-1=qk0，

由此得θ=(1-p)k0(1-αp)k0+1.對k≥k0+1，有

P(X=k)=αp(1-αp)k(1-p)k0(1-αp)k0+1

=αp(1-p)k0(1-αp)k-k0-1.



特別地，P(X=k0+1)=αp(1-p)k0.

值得一提的是，若k0→0，此時可看作產(chǎn)品一開始便在應(yīng)力S2下做試驗，于是有：

P(X=k)=limk0→0［α(k0)p(1-p)k0(1-α(k0)p)k-k0-1］

=α(0)p(1-α(0)p)k-1.



也就是說在恒應(yīng)力S2下產(chǎn)品壽命仍服從幾何分布，其參數(shù)為α(0)p.下面研究參數(shù)的極大似然估計：

似然函數(shù)為：（其中A為下常數(shù)）

L［α，p］=A∏ri=1［pqx(i)-1］∏ni=r+1［αp(1-p)x0#8226; (1-αp)x(j)-k0-1］

=Apr(1-p)∑ri=1x(i)-rαn-rpn-r(1-p)k0(n-r)#8226; (1-αp)∑NJ=R+1x(j)-(k0+1)(n-r)

=Apn(1-p)∑ri=1x(i)-rαn-r(1-p)k0(n-r)#8226; (1-αp)∑NJ=R+1x(j)-(k0+1)(n-r).



lnL［α，p］=lnA+nlnp+∑ri=1x(i)+k0n-(k0+1)r#8226; ln(1-p)+(n-r)lnα+ ∑NJ=R+1x(j)-(k0+1)(n-r)ln(1-αp).



lnL［α，p］α=n-rα-p∑nj=r+1x(j)-(k0+1)(n-r)1-αp.

lnL［α，p］p=np-∑ri=1x(i)+nk0-(k0+1)r1-p-

α∑nj=r+1x(j)-(k0+1)(n-r)1-αp.

令lnL［α，p］α=0，lnL［α，p］p=0，

化簡得α1-αp∑nj=r+1x(j)-(k0+1)(n-r)=n-rp，

np-∑ri=1x(i)+nk0-(k0+1)r1-p-n-rp=0，

rp=∑ri=1x(i)+nk0-(k0+1)r1-p，

r-rp=p∑ri=1x(i)+nk0-(k0+1)r，

r=p∑ri=1x(i)+(n-r)k0.

由此參數(shù)p的極大似然估計為：

=r∑ri=1x(i)+(n-r)k0.

又由于1-αpαp

=∑nj=r+1x(j)-(k0+1)(n-r)n-r

1αp=∑nj=r+1x(j)-(n-r)k0n-r，



進而得損傷因子α的極大似然估計為：

=n-rr#8226;∑ri=1x(i)+(n-r)k0∑nj=r+1x(j)-(n-r)k0.

例如，設(shè)在應(yīng)力S1下將10個產(chǎn)品投入試驗，當有5個產(chǎn)品失效時（次序失效“時間”為：120，223，296，321，386），在k0=400時將應(yīng)力提高至S2繼續(xù)做試驗，試驗持續(xù)到所有產(chǎn)品失效為止（次序失效“時間”為：411，458，466，496，518），運用本文方法得參數(shù)的極大似然估計為：=14943×10-3，=104304.