截尾試驗(yàn)（Censored Test）數(shù)據(jù)處理與MTBF推算

在可靠性測(cè)試中，由于時(shí)間、成本或?qū)嶋H條件的限制，我們往往無(wú)法等到所有樣品都失效才結(jié)束試驗(yàn)。這種包含未失效樣品（截尾數(shù)據(jù)）的試驗(yàn)稱(chēng)為截尾試驗(yàn)。正確處理截尾數(shù)據(jù)，并從中推算出可靠的MTBF（平均無(wú)故障時(shí)間），是可靠性工程的核心技能之一。

本文將詳細(xì)介紹截尾試驗(yàn)的基本概念、不同類(lèi)型截尾試驗(yàn)的數(shù)據(jù)處理方法、MTBF推算公式以及置信區(qū)間計(jì)算方法。

一、截尾試驗(yàn)的基本概念

1.1 什么是截尾數(shù)據(jù)？

截尾數(shù)據(jù)（Censored Data）是指在試驗(yàn)結(jié)束時(shí)，部分樣品仍未失效，我們只知道它們的壽命大于當(dāng)前試驗(yàn)時(shí)間，而不知道確切的失效時(shí)間。

完全數(shù)據(jù)與截尾數(shù)據(jù)的對(duì)比：

1.2 為什么會(huì)出現(xiàn)截尾數(shù)據(jù)？

1.3 截尾試驗(yàn)的類(lèi)型

二、定時(shí)截尾試驗(yàn)（Type I Censoring）

2.1 定義

定時(shí)截尾試驗(yàn)是指預(yù)先設(shè)定一個(gè)測(cè)試時(shí)間T，投入n個(gè)樣品進(jìn)行試驗(yàn)，當(dāng)達(dá)到時(shí)間T時(shí)立即停止，記錄發(fā)生的故障數(shù)r和每個(gè)故障的發(fā)生時(shí)間t?, t?, ..., t?。

特點(diǎn)：

• 試驗(yàn)時(shí)間T固定

• 故障數(shù)r隨機(jī)

• 未失效樣品在時(shí)間T截尾

2.2 MTBF點(diǎn)估計(jì)

對(duì)于指數(shù)分布，MTBF的點(diǎn)估計(jì)為：

$$\widehat{\theta }=\frac{T_{total}}{r}$$

其中總試驗(yàn)時(shí)間 $T_{total}={\sum }_{i=1}^r{t}_i+(n?r)\times T$

計(jì)算示例：

• 10臺(tái)設(shè)備測(cè)試1000小時(shí)

• 發(fā)生3次故障，故障時(shí)間分別為：200h、450h、800h

• 總試驗(yàn)時(shí)間 = 200 + 450 + 800 + 7 × 1000 = 8450小時(shí)

• MTBF點(diǎn)估計(jì) = 8450 / 3 = 2817小時(shí)

2.3 MTBF置信區(qū)間

單側(cè)置信下限（最常用）：

$${\theta }_L=\frac{2T_{total}}{{\chi }^2(2r+2,1?\alpha )}$$

雙側(cè)置信區(qū)間：

$$[\frac{2T_{total}}{{\chi }_{1?\alpha \mathrm{/}2}^2(2r+2)},\frac{2T_{total}}{{\chi }_{\alpha \mathrm{/}2}^2(2r)}]$$

示例（續(xù)上）：

• T_total = 8450小時(shí)

• r = 3

• 90%置信下限（α=0.1）：
  χ2(2×3+2=8, 0.9) = 13.36
  θ_L = 2×8450 / 13.36 = 1265小時(shí)

2.4 零失效特殊情況

當(dāng)r=0時(shí)，上述公式不能用，使用：

$${\theta }_L=\frac{T_{total}}{?\ln ?(1?C)}$$

其中C為置信度。

示例：

• T_total = 10000小時(shí)

• 無(wú)故障

• 90%置信度：θ_L = 10000 / 2.3026 = 4343小時(shí)

三、定數(shù)截尾試驗(yàn)（Type II Censoring）

3.1 定義

定數(shù)截尾試驗(yàn)是指預(yù)先設(shè)定一個(gè)故障數(shù)r，投入n個(gè)樣品進(jìn)行試驗(yàn)，當(dāng)達(dá)到第r個(gè)故障時(shí)立即停止，記錄每個(gè)故障的時(shí)間t?, t?, ..., t?（按順序排列）。

特點(diǎn)：

• 故障數(shù)r固定

• 試驗(yàn)時(shí)間T隨機(jī)

• 未失效樣品在第r個(gè)故障時(shí)間截尾

3.2 MTBF點(diǎn)估計(jì)

對(duì)于指數(shù)分布，MTBF的點(diǎn)估計(jì)為：

$$\widehat{\theta }=\frac{T_{total}}{r}$$

其中總試驗(yàn)時(shí)間 $T_{total}={\sum }_{i=1}^r{t}_i+(n?r)\times t_r$

計(jì)算示例：

• 10臺(tái)設(shè)備測(cè)試到第5次故障停止

• 故障時(shí)間：100h、250h、400h、600h、900h

• 總試驗(yàn)時(shí)間 = 100+250+400+600+900 + 5×900 = 6750小時(shí)

• MTBF點(diǎn)估計(jì) = 6750 / 5 = 1350小時(shí)

3.3 MTBF置信區(qū)間

單側(cè)置信下限：

$${\theta }_L=\frac{2T_{total}}{{\chi }^2(2r,1?\alpha )}$$

雙側(cè)置信區(qū)間：

$$[\frac{2T_{total}}{{\chi }_{1?\alpha \mathrm{/}2}^2(2r)},\frac{2T_{total}}{{\chi }_{\alpha \mathrm{/}2}^2(2r)}]$$

示例（續(xù)上）：

• T_total = 6750小時(shí)

• r = 5

• 90%置信下限（α=0.1）：
  χ2(2×5=10, 0.9) = 15.99
  θ_L = 2×6750 / 15.99 = 844小時(shí)

四、兩種截尾試驗(yàn)的對(duì)比

五、隨機(jī)截尾試驗(yàn)

5.1 定義

隨機(jī)截尾是指由于各種原因（如樣品退出、試驗(yàn)中斷等），樣品的截尾時(shí)間各不相同。

5.2 數(shù)據(jù)處理方法

對(duì)于隨機(jī)截尾數(shù)據(jù)，通常采用Kaplan-Meier估計(jì)法（乘積極限法）或Nelson-Aalen估計(jì)法。

Kaplan-Meier估計(jì)公式：

$$\widehat{R}(t)=\prod _{i:t_i\le t}(1?\frac{d_i}{n_i})$$

其中：

• d_i：在時(shí)間t_i的失效數(shù)

• n_i：在t_i之前仍處于風(fēng)險(xiǎn)中的樣品數(shù)

六、不同分布的MTBF推算

6.1 指數(shù)分布

指數(shù)分布是最常用的壽命分布，特點(diǎn)是失效率恒定。

MTBF = θ（尺度參數(shù)）

6.2 威布爾分布

威布爾分布可以描述失效率隨時(shí)間變化的情況。

概率密度函數(shù)：

$$f(t)=\frac{\beta }{\eta }{\left(\frac{t}{\eta }\right)}^{\beta ?1}{e}^{?(t\mathrm{/}\eta {)}^{\beta }}$$

MTBF公式：

$$MTBF=\eta ?\mathrm{\Gamma }(1+\frac{1}{\beta })$$

其中Γ為伽馬函數(shù)。

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

MTBF = exp(μ + σ2/2)

七、數(shù)據(jù)處理流程

7.1 數(shù)據(jù)整理

7.2 總試驗(yàn)時(shí)間計(jì)算

通用公式：

$$T_{total}=\sum _{i=1}^n{t}_i^{?}$$

其中t_i^*是每個(gè)樣品的最終時(shí)間（失效時(shí)間或截尾時(shí)間）。

7.3 MTBF計(jì)算

點(diǎn)估計(jì)： θ? = T_total / r
置信區(qū)間： 根據(jù)截尾類(lèi)型選擇公式

7.4 示例：混合數(shù)據(jù)

數(shù)據(jù)：

• 樣品1：失效@100h

• 樣品2：失效@200h

• 樣品3：截尾@300h

• 樣品4：失效@350h

• 樣品5：截尾@500h

• 樣品6：失效@600h

• 樣品7：截尾@600h

• 樣品8~10：截尾@600h（定時(shí)截尾）

計(jì)算：

• r = 4（樣品1、2、4、6）

• T_total = 100 + 200 + 300 + 350 + 500 + 600 + 600 + 600 + 600 + 600 = 4450小時(shí)

• θ? = 4450 / 4 = 1112.5小時(shí)

90%置信下限（定時(shí)截尾）：
χ2(2×4+2=10, 0.9) = 15.99
θ_L = 2×4450 / 15.99 = 556小時(shí)

八、軟件工具

8.1 常用軟件

8.2 Excel計(jì)算公式

九、常見(jiàn)問(wèn)題與解答

Q1: 為什么定時(shí)截尾用χ2(2r+2)而定數(shù)截尾用χ2(2r)？

A: 這是因?yàn)槎〞r(shí)截尾在試驗(yàn)結(jié)束時(shí)還可能發(fā)生一次故障，統(tǒng)計(jì)上相當(dāng)于多了一個(gè)自由度。

Q2: 故障數(shù)很少時(shí)怎么辦？

A: 故障數(shù)少時(shí)置信區(qū)間會(huì)很寬，這是正常現(xiàn)象。可以：

• 增加樣品數(shù)

• 延長(zhǎng)測(cè)試時(shí)間

• 接受較寬的置信區(qū)間

Q3: 零失效時(shí)如何報(bào)告？

A: 零失效時(shí)只能報(bào)告單側(cè)置信下限，不能報(bào)告點(diǎn)估計(jì)。常用公式：θ_L = T_total / (-ln(1-C))

Q4: 如何選擇截尾方式？

A:

• 時(shí)間嚴(yán)格受限：選擇定時(shí)截尾

• 精度要求高：選擇定數(shù)截尾

• 一般情況：定時(shí)截尾更常用

十、小結(jié)

截尾試驗(yàn)的數(shù)據(jù)處理是可靠性測(cè)試的核心技能：

掌握截尾試驗(yàn)的數(shù)據(jù)處理方法，能夠從有限的測(cè)試數(shù)據(jù)中科學(xué)地推算出產(chǎn)品的MTBF，為可靠性評(píng)估提供依據(jù)。

訊科標(biāo)準(zhǔn)檢測(cè)
ISTA認(rèn)可實(shí)驗(yàn)室 | CMA | CNAS
地址：深圳寶安

訊科標(biāo)準(zhǔn)檢測(cè)是一家專(zhuān)業(yè)的第三方檢測(cè)機(jī)構(gòu)，已獲得CNAS、CMA及ISTA等多項(xiàng)資質(zhì)認(rèn)可。實(shí)驗(yàn)室可提供可靠性測(cè)試方案設(shè)計(jì)、截尾試驗(yàn)數(shù)據(jù)處理、MTBF推算分析等技術(shù)服務(wù)。檢測(cè)報(bào)告可用于產(chǎn)品質(zhì)檢、市場(chǎng)準(zhǔn)入及客戶(hù)驗(yàn)證等場(chǎng)景。

?? 咨詢(xún)熱線(xiàn)：0755-27909791 / 15017918025（同微）
?? 郵箱：cs@xktest.cn
地址：深圳市寶安區(qū)航城街道