PosWeibull對數(shù)正態(tài)分布計算結(jié)果準確性驗證及使用注意事項
一、對數(shù)正態(tài)分布模型及應用范圍
對數(shù)正態(tài)分布是可靠性工程中較為常用的一個壽命分布模型。對數(shù)正態(tài)分布模型的兩個參數(shù)分別為均值mu和標準差σ。對數(shù)正態(tài)分布與正態(tài)分布類似，但是，對數(shù)正態(tài)分布是假設隨機變量值的對數(shù)值呈正態(tài)分布而不是隨機變量本身呈正態(tài)分布。因此，對數(shù)正態(tài)分布的所有隨機變量值都是正的，并且分布形狀向左傾斜。
在開展壽命數(shù)據(jù)分析時，通常會遇到到底選擇威布爾分布還是對數(shù)正態(tài)分布更合適的問題。對數(shù)正態(tài)分布模型和威布爾分布模型均適用于較多類型的壽命數(shù)據(jù)。但是，對數(shù)正態(tài)分布模型，它在工程中適用于描述失效模式具有疲勞應力特性的產(chǎn)品壽命分布規(guī)律（大部分的機械產(chǎn)品都具有該特性），經(jīng)常用于金屬疲勞試驗數(shù)據(jù)、維修數(shù)據(jù)（維修時間）、化學工藝設備的故障和維修數(shù)據(jù)、材料特性以及非線性加速退化數(shù)據(jù)分析。當失效時間是由效應相乘（疊加）得到時，建議使用對數(shù)正態(tài)分布。例如，設備隨著環(huán)境應力不斷惡化而出現(xiàn)了裂紋，且裂紋隨著應力增加不斷增大，類似這樣的數(shù)據(jù)可以使用對數(shù)正態(tài)分布模型進行分析。
`f(t')=\frac{1}{\sigma'\sqrt{2\pi}}e^{\-1/2(\frac{t'-u'}{\sigma'})^2}`
`t'`是故障時間t的自然對數(shù)值，即`t'=ln?(t)`。需要記住這個特征，后續(xù)的換算需要用到。
同樣，`u'`、`σ'`分別為故障時間t的自然對數(shù)的均值、標準差。
二、PosWeibull的對數(shù)正態(tài)分布計算準確性測試
為驗證PosWeibull軟件的計算精度，選擇國外同類的Re軟件、國外的M軟件與PosWeibull的22版本進行對比驗證。
（1）測試案例1
某產(chǎn)品可靠性試驗的故障數(shù)據(jù)如下：2，5，11，23，29，37，43，59，單位是小時。
現(xiàn)選擇對數(shù)正態(tài)分布作為擬合的分布類型，計算方法分別選擇極大似然法，置信度取0.95（雙側(cè)），分別使用PosWeibull的壽命分析模塊、Re軟件的Weibull模塊、M軟件對上述數(shù)據(jù)進行計算。計算結(jié)果如下表所示。
結(jié)果分析：選擇極大似然法時，三個工具計算得到的均值參數(shù)`u`、標準差`σ`參數(shù)`η`基本一致，誤差小于0.01%。其中，PosWeibull的計算結(jié)果與M軟件的完全一致。
表1 案例1計算結(jié)果對比（極大似然法）
| 輸入數(shù)據(jù) | ||||
|---|---|---|---|---|
| 2，5，11，23，29，37，43，59 | 是否刪失 | 無刪失 | ||
| 計算設置 | 計算軟件 | PosWeibull | Re軟件 | M軟件 |
| 計算方法 | 極大似然 | 極大似然 | 極大似然 | |
| 置信度 | 0.95（雙側(cè)） | 0.95（雙側(cè)） | 0.95（雙側(cè)） | |
| 分布類型 | 對數(shù)正態(tài)分布 | 對數(shù)正態(tài)分布 | 對數(shù)正態(tài)分布 | |
| 計算結(jié)果 | 均值`u` | 2.83162 | 2.83（四舍五入） | 2.83162 |
| 標準差`σ` | 1.09875 | 1.10（四舍五入） | 1.09875 | |
果.png)
圖1-PosWeibull軟件計算結(jié)果
果.png)
圖2-M軟件計算結(jié)果
（2）測試案例2
使用Kececioglu 所編的《Reliability & Life Testing Handbook》資料里面的的案例進行對比驗證。
故障數(shù)據(jù)點分別為：62.5、91.9、100.3、117.4、141.1、146.8、172.7、192.5、201.6、235.8、249.2、297.5、318.3、410.6、550.5。將該數(shù)據(jù)集的數(shù)據(jù)分別錄入到PosWeibull的壽命分析模塊、Re軟件的Weibull模塊、M軟件進行計算。計算方法選擇RRX方法。
《Reliability & Life Testing Handbook》資料通過概率圖方法給出的結(jié)果為`u`=5.22575，`σ`=0.62048。
Re軟件計算結(jié)果為`u`=5.2303，σ=0.6283。
M軟件計算結(jié)果為`u`=5.23028，σ=0.630225。
PosWeibull軟件的計算結(jié)果為`u`=5.23028，`σ`=0.630225。
可以看出，PosWeibull軟件的計算結(jié)果與M軟、Re軟件的基本一致，與《Reliability & Life Testing Handbook》資料通過概率圖方法給出的結(jié)果也基本一致，存在細微差別的原因是《Reliability & Life Testing Handbook》資料的計算方法是通過手動繪制擬合線的，得到的結(jié)果精度相對較差一些。
表2 案例2計算結(jié)果對比（最小二乘法）
| 輸入數(shù)據(jù) | ||||
|---|---|---|---|---|
| 62.5、91.9、100.3、117.4、141.1、146.8、172.7、192.5、201.6、235.8、249.2、297.5、318.3、410.6、550.5 | 是否刪失 | 無刪失 | ||
| 計算設置 | 計算軟件 | PosWeibull | Re軟件 | M軟件 |
| 計算方法 | RRX | RRX | RRX | |
| 置信度 | 0.95（雙側(cè)） | 0.95（雙側(cè)） | 0.95（雙側(cè)） | |
| 分布類型 | 對數(shù)正態(tài)分布 | 對數(shù)正態(tài)分布 | 對數(shù)正態(tài)分布 | |
| 計算結(jié)果 | 均值`u` | 5.23028 | 5.2303 | 5.23028 |
| 標準差`σ` | 0.630225 | 0.6283 | 0.630225 | |
果.png)
圖3-PosWeibull軟件計算結(jié)果
果.png)
圖4-M軟件計算結(jié)果
（3）測試案例3：較多刪失數(shù)據(jù)
以Nelson的《Applied Life Data Analysis》的機車控制設備的故障數(shù)據(jù)為例，分別使用Re軟件、PosWeibull軟件、M軟件進行計算。一共96個樣本，其中37個樣本出現(xiàn)故障，59個樣本未故障（刪失）。
正態(tài)故障刪失.png)
Re軟件的計算結(jié)果為`u`=2.2223，`σ`=0.3064，均為10為底的對數(shù)的結(jié)果。
M軟件的計算結(jié)果為`u`=5.11692，`σ`=0.70549，均為自然對數(shù)的結(jié)果。
PosWeibull的計算結(jié)果為`u`=5.11692，`σ`=0.705494，均為自然對數(shù)的結(jié)果。
由于Re軟件給出的是10為底的對數(shù)的計算結(jié)果，將PosWeibull軟件的計算結(jié)果轉(zhuǎn)為10為底的對數(shù)的計算結(jié)果，轉(zhuǎn)換后的PosWeibull的計算結(jié)果（10為底）為`u'`=2.22225，`σ'`=0.306351?；九cRe軟件的計算結(jié)果一致。表明在使用對數(shù)正態(tài)分布處理包含較多的刪失數(shù)據(jù)的時候，PosWeibull所計算得到的結(jié)果與國外Re等軟件的計算結(jié)果基本一致，滿足工程需求。但是，使用的時候需要特別注意，PosWeibull軟件對數(shù)正態(tài)分布計算得到的是取對數(shù)后的變量的計算結(jié)果??！
表3 案例3計算結(jié)果對比（極大似然法）
| 輸入數(shù)據(jù) | ||||
|---|---|---|---|---|
| 22.5,37.5,46,48.5,51.5,53,54.5,57.5,66.5, 68,69.5,76.5,77,78.5,80,81.5, 82,83,84,91.5,93.5,102.5,107,108.5,112.5,113.5,116, 117,118.5,119,120,122.5,123,127.5,131,132.5,134,135(+,59) | 是否刪失 | 存在刪失，刪失標記(+) | ||
| 計算設置 | 計算軟件 | PosWeibull | Re軟件 | M軟件 |
| 計算方法 | 極大似然 | 極大似然 | 極大似然 | |
| 置信度 | 0.95（雙側(cè)） | 0.95（雙側(cè)） | 0.95（雙側(cè)） | |
| 分布類型 | 對數(shù)正態(tài)分布 | 對數(shù)正態(tài)分布 | 對數(shù)正態(tài)分布 | |
| 計算結(jié)果 | 均值`u` | 5.11692 | 2.2223 | 5.11692 |
| 標準差`σ` | 0.705494 | 0.3064 | 0.70549 | |
圖5-機車設備故障概率圖
三、使用注意事項
（1）在使用PosWeibull的對數(shù)正態(tài)分布模型進行壽命數(shù)據(jù)分析的時候，需要注意的是，所計算得到的均值`u'`、標準差`σ'`均為自然對數(shù)（e為底）時的結(jié)果。如果要利用PosWeibull軟件計算得到的均值、標準差進行故障時間預測時，可采用如下方法：a、利用PosWeibull軟件計算所得到的均值`u'`、標準差`σ'`，計算出對應的故障時間`t'`。b、根據(jù)`t'=ln(t)`關系，使用`t=exp(t')`換算出`t`的值。
（2）使用對數(shù)正態(tài)分布時，需要注意與威布爾分布的相似性和差異性。