Weibull 可靠度分析

💡 支援單組或雙組數據分析與比對

分析模式 (Analysis Mode):
組別 A
序號 壽命 狀態 刪除
組別 B
序號 壽命 狀態 刪除

分析結果

$$ R(t) = e^{-(t/\eta)^\beta} $$

組別 A

β -
η -
-

組別 B

β -
η -
-

機率圖 (Probability Plot)

可靠度圖 (Reliability Plot)

📋 批次輸入數據
使用說明:
• 從 Excel 或 CSV 複製數據,直接貼上到下方文字框
• 格式:每行一筆數據,用 Tab、逗號或空格分隔
• 第一欄:壽命 (t),第二欄:狀態 (F=失效, S=右設限)
• 範例:100 F200,S300 F
📘 Weibull 分析理論指南
基礎概念
參數判讀
B-Life 指南
樣本數建議
雙組比對

1. 數據類型:失效與右設限

在可靠度測試中,我們會遇到兩種數據類型:

失效數據 (Failure, F)

產品在測試期間確實發生失效的數據點。例如:產品在運作 300 小時後損壞。這是最直接、最有價值的數據。

右設限數據 (Suspended/Censored, S)

右設限是指產品在測試結束時仍然正常運作,我們只知道它「至少」存活了這麼久,但不知道實際壽命。

為什麼會有右設限數據?

  • 測試時間限制:產品測試 500 小時後結束,但有些樣品還沒壞 → 記錄為「500 小時,設限 (S)」
  • 樣品移除:測試中途因其他原因(非失效)移除樣品,例如送去做其他檢測
  • 成本考量:長時間測試成本高,無法等到所有樣品都失效

重要:右設限數據雖然沒有失效,但仍然提供了「產品至少能活這麼久」的資訊,對於提升分析準確度非常有幫助!

實際範例

假設測試 5 個產品,測試時間上限 600 小時:

  • 產品 A:300 小時失效 → 300, F
  • 產品 B:450 小時失效 → 450, F
  • 產品 C:600 小時仍正常 → 600, S (右設限)
  • 產品 D:550 小時失效 → 550, F
  • 產品 E:600 小時仍正常 → 600, S (右設限)

這樣的數據組合可以進行 Weibull 分析,右設限數據會被正確納入計算。

1.5 可靠度與不可靠度的關係

在工程實務中,產品的狀態是一體兩面的:不是「存活」就是「已失效」。因此我們定義了互補的兩個指標:

  • \(R(t)\) 可靠度 (Reliability)
    通俗名稱:存活率 (Survival Rate)
    定義:代表產品運作到時間 \(t\) 仍未損壞的機率。
  • \(F(t)\) 不可靠度 (Unreliability)
    通俗名稱:累積失效機率 (Cumulative Probability of Failure)
    定義:代表產品從開始運作直到時間 \(t\) 已經損壞的總比例(即 CDF)。
$$ R(t) + F(t) = 1 $$

(存活率) + (累積失效機率) = 100%

2. 核心公式

Weibull 分佈的 可靠度函數 (Reliability Function) \(R(t)\) 公式如下:

$$ R(t) = e^{-(t/\eta)^\beta} $$
  • \(t\):壽命(或行駛里程、循環次數)。
  • \(R(t)\):在時間 \(t\) 時,產品尚未失效的機率(可靠度)。
  • \(\beta\) (Beta):形狀參數 (Shape Parameter),決定失效模式。
  • \(\eta\) (Eta):尺度參數 (Scale Parameter),又稱特徵壽命。

3. 理解關鍵參數

在代入公式前,必須先透過數據擬合求出 \(\beta\) 和 \(\eta\)。

\(\beta\) (形狀參數):決定失效類型

  • \(\beta < 1\):早期失效期 (Infant Mortality)。產品剛出廠就壞,通常是生產瑕疵。
  • \(\beta = 1\):隨機失效期 (Random Failures)。失效律是常數,類似指數分佈。
  • \(\beta > 1\):磨耗失效期 (Wear-out)。產品老化、疲勞,失效律隨時間上升。

\(\eta\) (尺度參數):特徵壽命

代表當 63.2% 的產品已經失效時的時間點(此時 \(R(t) = 36.8\%\))。這是衡量產品壽命長短的主要指標。

\(R^2\) (擬合度):數據可信度指標

擬合度 (R-squared, R²) 是用來評估 Weibull 模型與實際數據的吻合程度,數值範圍為 0 到 1。

  • \(R^2 \geq 0.95\):優秀擬合,數據非常符合 Weibull 分佈,分析結果高度可信。
  • \(0.90 \leq R^2 < 0.95\):良好擬合,結果可接受,但可能存在少量異常數據。
  • \(0.80 \leq R^2 < 0.90\):尚可擬合,建議檢查數據品質或考慮其他分佈模型。
  • \(R^2 < 0.80\):擬合不佳,數據可能不符合 Weibull 分佈,或存在嚴重異常值。

注意:如果 \(R^2\) 過低,建議重新檢視數據來源、測試條件是否一致,或考慮使用對數常態分佈等其他模型。

4. 計算步驟實例

假設你分析一組數據後得到參數:\(\beta = 2.5\) (磨耗型) 與 \(\eta = 1000\) 小時

情境 A:計算特定時間的可靠度

問題:產品運作 500 小時 後,還存活的機率是多少?

代入公式:

$$ R(500) = e^{-(500/1000)^{2.5}} $$ $$ R(500) = e^{-(0.5)^{2.5}} = e^{-0.1768} \approx 0.838 $$

答案: 可靠度為 83.8%


情境 B:計算 B10 壽命 (B10 Life)

問題:多少時間後,會有 10% 的產品失效(即 90% 可靠度)?

公式變形: \( t = \eta \times (-\ln(R))^{1/\beta} \)

$$ t = 1000 \times (-\ln(0.9))^{1/2.5} $$ $$ t = 1000 \times (0.10536)^{0.4} $$ $$ t \approx 1000 \times 0.406 = 406 \text{ 小時} $$

答案: B10 壽命為 406 小時

4. B-Life 概念與工程意義

在可靠度工程(特別是 Weibull 分析)中,「B」 代表 Basic Life(基本壽命),後面的數字代表 累積失效百分比。簡單來說,Bx 壽命就是:「當 x% 的產品失效時,所經過的時間」。以下是 B1、B10、B63.2 的詳細定義與它們在工程上的意義:

1. B10 壽命 (B10 Life) —— 行業通用標準

定義:當 10% 的產品失效(也就是還有 90% 存活)的時間點。

數學意義:可靠度 \(R(t) = 90\%\) 或 \(0.9\)。

應用場景:這是汽車業、軸承業、馬達業最常用的指標。通常用來設定 保固期 (Warranty)

實際例子:如果車廠算出某零件的 B10 壽命是 3年,代表 3年內預計會有 10% 的車回來維修,這有助於計算保固成本。

2. B1 壽命 (B1 Life) —— 高安全標準

定義:當 1% 的產品失效(也就是還有 99% 存活)的時間點。

數學意義:可靠度 \(R(t) = 99\%\) 或 \(0.99\)。

應用場景:用於高安全性、高可靠度的產品,如 航太、醫療設備、汽車安全氣囊。對於這些產品,10% 的失效是不可接受的,必須確保極少量的產品會在該時間點前壞掉。

3. B63.2 (特徵壽命 / \(\eta\)) —— 數學基準

定義:當 63.2% 的產品失效(只剩 36.8% 存活)的時間點。

為什麼是 63.2%?這來自 Weibull 的公式結構。當時間 \(t\) 正好等於尺度參數 \(\eta\) 時:

$$ R(\eta) = e^{-(\eta/\eta)^\beta} = e^{-1} \approx 0.3678 $$ $$ F(t) = 1 - 0.3678 = 0.6322 \text{ (即 63.2%)} $$

意義:B63.2 就是 \(\eta\) (Eta, 尺度參數)。它代表這批產品的「特徵壽命」(Characteristic Life)。雖然 63% 的失效率聽起來很高(大半都壞了),但它在數學上是一個非常穩定的基準点,用來衡量整體壽命的長短。如果 \(\eta\) 越大,代表產品整體撐得越久。

總結比較表

指標 失效比例 存活率 (可靠度) 意義與用途
B1 1% 99% 極低風險:用於飛機、醫療、安全件。
B10 10% 90% 工業標準:用於制定保固期、馬達壽命、一般電子產品。
B50 50% 50% 中位數壽命:一半壞了一半好的時間點 (平均壽命)。
B63.2 63.2% 36.8% 特徵壽命 (\(\eta\)):Weibull 的尺度參數,數學上的基準點。

如何透過我們剛做的網頁計算?

在我們的網頁工具算出 \(\beta\) 和 \(\eta\) 後,您可以手動用以下公式反推這些 B-Life:

$$ t = \eta \times (-\ln(R))^{1/\beta} $$
  • 算 B1:將 R 設為 0.99
  • 算 B10:將 R 設為 0.90
  • 算 B63.2:答案就是網頁上顯示的 Eta (\(\eta\)) 值。

5. 樣本數量建議:多少顆才夠?

進行可靠度分析(特別是 Weibull 分析)時,樣本數量的選擇永遠是「統計精確度」與「測試成本/時間」之間的權衡。雖然數學上只需要 2 個失效點就能畫出一條線(求出 \(\beta\) 和 \(\eta\)),但在工程實務上,這是遠遠不夠的。

建議的樣本數量 (Rule of Thumb)

樣本數量 評價 適用情境
3 - 5 顆 極低 (高風險) 僅適用於非常昂貴的產品(如衛星、大型模具)或初步摸底測試。誤差範圍極大,僅能看大致趨勢。
6 - 10 顆 勉強接受 適用於開發階段的快速驗證。可以算出參數,但對於 B10 壽命(早期失效)的預測仍不夠準確。
10 - 20 顆 工程推薦 (甜蜜點) 大多數工業標準的選擇。在成本與精確度之間取得最佳平衡,足以判斷失效模式 (\(\beta\)) 並估算 B10 壽命。
20 - 30 顆+ 統計顯著 適用於高可靠度要求(如車規、醫療、航太)或低成本元件。能提供較窄的信賴區間。

為什麼?(背後的統計原因)

A. 信賴區間 (Confidence Bounds) 的寬度

Weibull 分析算出的只是「估計值」。樣本越少,不確定性(信賴區間)就越寬。

例子:

  • 5 顆樣品:算出 B10 壽命是 500 小時,但 90% 信賴區間可能是 [100 小時, 2500 小時]。這範圍大到幾乎無法做決策。
  • 20 顆樣品:算出 B10 壽命是 500 小時,信賴區間可能縮小至 [400 小時, 650 小時]。這數據就很有參考價值。

B. 為了抓到「早期失效」(B1/B10)

工程師通常不關心平均壽命 (\(\eta\), B63.2),而是關心「最早壞的那幾顆」 (B10 或 B1)。

  • Weibull 分佈的尾端(早期失效區)對數據非常敏感。
  • 如果樣本太少,你很難確定最早失效的那顆是「隨機的特例」還是「整批產品的常態」。
  • 數據量大才能穩定描繪出曲線的尾端。

C. 確認失效模式 (\(\beta\) 值)

\(\beta\) 值決定了產品是「早夭」、「隨機失效」還是「磨耗」。

樣本少時,隨便一顆異常值(Outlier)都會讓直線斜率劇烈波動,導致你誤判失效模式(例如把磨耗誤判為隨機失效),進而選錯改善方向。

如果預算/時間有限,樣本很少怎麼辦?(Mouldex 實務建議)

在模具或高單價產品分析中,往往不可能測 20 顆。這時有幾種策略:

策略一:使用「Weibayes 分析」(1參數 Weibull)

如果你只有很少的樣本(例如 2-3 顆),甚至沒有失效(全都是 Suspended),但你有過去類似產品的經驗。

方法:固定 \(\beta\) 值(根據歷史經驗,例如金屬疲勞通常 \(\beta \approx 3\)),只計算 \(\eta\)。

優點:需要的樣本數極少,且結果比標準 Weibull 更可信(因為利用了歷史知識)。

策略二:增加「設限數據」(Suspended Data)

方法:你可以只測到 2-3 顆失效,但同時放入 10 顆樣品一起跑。當前 3 顆失效後,停止測試。

原因:那 7 顆未失效的數據(右設限)雖然沒壞,但它們提供了「這段時間內沒壞」的資訊,這能顯著提升分析的準確度,比只測 3 顆全部測到壞要好得多。

總結建議

對於 Mouldex 的應用場景

  • 理想情況:盡量爭取 10 顆以上的樣本。
  • 最低底線:至少要有 5-6 顆
  • 測試策略:如果只能測到少數幾顆壞掉,請務必保留那些「還沒壞」的數據作為右設限 (Suspended) 輸入到我們剛做好的網頁工具中,這樣能大幅提升準確度。

6. 如何取得參數 (\(\beta, \eta\))

在實際工程中,我們不會手算參數,通常使用圖解法或軟體(如上方計算機):

圖解法 (Excel 原理)

利用 Weibull 線性化公式製作散佈圖:

$$ \ln(-\ln(1 - F(t))) = \beta \ln(t) - \beta \ln(\eta) $$

在圖表中,趨勢線的斜率即為 \(\beta\)

Python 程式範例

如果你熟悉程式,使用 Python 的 lifelines 套件是最快的方法:

from lifelines import WeibullFitter

# 1. 輸入失效時間數據
T = [100, 150, 200, 280, 400]

# 2. 擬合數據
wf = WeibullFitter()
wf.fit(T)

# 3. 輸出參數
print(f"Beta (rho_): {wf.rho_}")
print(f"Eta (lambda_): {wf.lambda_}")

7. 為什麼需要較多樣本?

A. 信賴區間的寬度

樣本越少,不確定性越高。5 顆樣品的 B10 信賴區間可能非常寬;20 顆則明顯收斂。

B. 早期失效尾端的敏感度

B1/B10 屬於分佈尾端,對少量異常值極為敏感,需較多樣本穩定估計。

C. 確認失效模式 (\(\beta\) 值)

樣本太少時,斜率容易受單點影響而誤判失效型態,導致錯誤改善方向。

8. 實務策略

策略一:Weibayes (1 參數 Weibull)

樣本極少或多為設限時,依歷史經驗固定 \(\beta\),僅估 \(\eta\)。

策略二:增加設限數據

少量失效搭配更多設限樣本,可顯著提升擬合穩定度。

10. 差異分析的理論補充(雙組比較)

\(\beta\) 差異的意義

  • \(\beta_B > \beta_A\):B 組失效率隨時間上升更快,屬磨耗加劇,尾端失效更集中
  • \(\beta_B \approx \beta_A\):兩組失效型態相近,適合用 \(\eta\) 或 B‑Life 直接比較壽命
  • \(\beta_B < \beta_A\):B 組較傾向早期失效,需檢討製程/品質一致性

\(\eta\) 改善率的判讀

$$ \text{\% 改善} = \left(\frac{\eta_B}{\eta_A} - 1\right) \times 100\% $$
  • 適用情境:\(|\beta_B - \beta_A| \le 0.2\)(失效型態近似)
  • 解讀:\(\eta_B/\eta_A = 1.30\) 表示壽命延長約 30%

當 \(\beta\) 不同時的比較方式

  • 固定可靠度比較:選定 \(R\)(如 95%、90%),比較兩組的 \(t(R)\)
$$ t(R) = \eta \times \big(-\ln(R)\big)^{1/\beta} $$
  • B‑Life 比較:以 \(B10\)、\(B50\) 或 \(B1\) 為指標,比較 \(t_{B,\,x} / t_{A,\,x}\)
  • 優點:避免僅以 \(\eta\) 判斷而忽略 \(\beta\) 對尾端的影響

判讀範例(文字化)

例一(型態相近):\(\beta_A = 1.2,\ \eta_A = 500;\ \beta_B = 1.2,\ \eta_B = 650\)。則 \(\eta\) 改善 ≈ 30%,B10 比較亦約 30% 提升。

例二(型態不同):\(\beta_A = 1.0,\ \eta_A = 520;\ \beta_B = 1.6,\ \eta_B = 650\)。以 \(R=90\%\) 比較:計算 \(t_A(0.9)\) 與 \(t_B(0.9)\);若 \(t_B/t_A > 1\),表示在該可靠度下 B 組更耐久。

判讀優先順序

步驟 重點 說明
1 檢查 \(R^2\) 擬合度不足時避免深入比較
2 比較 \(\beta\) 型態相近者可直接用 \(\eta\)/B‑Life;型態不同用 \(t(R)\)
3 選擇比較指標 相近:\(\eta\) 改善率;不同:固定 \(R\) 或 B‑Life
4 形成結論 用明確百分比或時間比值呈現差異