ワイブル分布｜信頼性工学で使う代表分布

ワイブル分布

ワイブル分布は、寿命データや強度データの解析に広く用いられる連続確率分布である。形状・尺度・位置の各パラメータにより、初期故障から摩耗故障まで多様な故障様式を表現でき、信頼性工学や品質管理、材料強度、風工学などで標準的に用いられる。確率密度関数は、位置γ、尺度η（eta）、形状β（beta）を用いて f(t)= (β/η){((t-γ)/η)}^{β-1} exp{-((t-γ)/η)^β}（t≥γ）で与えられ、柔軟性の高さが特徴である。

Table Of Contents

ワイブル分布

定義とパラメータ

ワイブル分布は3母数形が最も一般的で、(β, η, γ) の三つで分布形が決まる。βは分布の形状（ハザードの増減）を制御し、ηは特性時間（尺度）として時刻のスケールを決め、γは故障が起こり始める基準（位置）を平行移動する。累積分布関数は F(t)=1-exp{-((t-γ)/η)^β}、生存関数（信頼度関数）は R(t)=exp{-((t-γ)/η)^β} である。βが同じでもηが大きいほど寿命が長い系を表すため、設計比較や耐久向上評価に有用である。

2母数形と3母数形

位置γ=0を仮定すると2母数形となり、実装が簡潔で推定も安定しやすい。一方、磨耗が始まるまでの待ち時間や規格下限が明確な系では3母数形が適合する。データサイズが小さい場合は過剰適合を避け、まず2母数形から適用し、残差や適合度で3母数形を検討する手順が実務的である。

信頼性指標とハザード

ハザード関数は h(t)=f(t)/R(t)=(β/η){((t-γ)/η)}^{β-1} で、βにより単調減少・一定・単調増加のいずれも表現できる。平均寿命は E[T]=γ+η·Γ(1+1/β)、中央値は Med[T]=γ+η(ln2)^{1/β} で与えられる。分位点は設計余寿命の見積りに直接使え、保全周期や保証期間の設定根拠となる。信頼水準qに対する寿命は t_q=γ+η{−ln(1−q)}^{1/β} で計算できる。

指数分布・レイリーとの関係

β=1で指数分布に一致し、ハザード一定（ランダム故障期）を表す。β=2ではレイリー分布となり、累積損傷や摩耗の影響が強い系を表現する。したがってワイブル分布は「指数（β=1）から摩耗（β>1）、初期欠陥（β<1）」までを連続的に包含する上位概念として理解するとよい。

形状パラメータβと故障様式

β<1は初期故障期（ハザード減少）で、製造欠陥や取付不良、表面欠陥支配の材料に対応する。β=1は偶発故障期、β>1は摩耗故障期（ハザード増加）で、摩耗・腐食・疲労亀裂進展など累積損傷を受ける機械要素に適合しやすい。設備のバスタブ曲線を分割し、区間ごとに別パラメータで当てるモデリングも行われる。

材料強度と弱リンク仮説

セラミックスや複合材では、最弱リンク（weakest-link）仮説の下でワイブル分布が破壊強度分布を良く表す。体積や試験片寸法の影響をスケーリングで扱えるため、試験から実機への外挿が比較的明瞭になる。

推定法とWeibullプロット

母数推定には最尤法（MLE）が基本で、尤度方程式を数値的に解く。打切り（右打切り・左打切り・区間打切り）を含む寿命試験でもMLEは自然に拡張できる。初期検討では線形化プロット（Weibullプロット）が有効で、横軸 ln(t−γ)、縦軸 ln{−ln(1−F̂)} をとると直線近似され、傾きがβ、切片がηに対応する。経験的累積 F̂ は中央値推定（Benard近似など）を用いる。

適合度とモデル診断

コルモゴロフ–スミルノフ検定やアンダーソン–ダーリング検定で分布仮説を検証し、残差プロットやQ–Qプロットで系統偏差を点検する。βが区間で変化しそうな場合は混合ワイブル分布や段階モデルも候補となる。

実務での使い方

保全では、信頼水準に基づく予防交換時期 t_q の設定、保証解析、予備品点検間隔の設計に直結する。設計ではワイブル分布から得たβ・ηをもとに安全率やデレーティング方針を整え、MTBFや可用度計算に反映する。試作段階では初期故障（β<1）の検出に重点を置き、スクリーニングやバーンインの要否を判断する。

打切りデータと情報活用

長寿命製品では全数故障を待てないため右打切りが常態である。MLEは打切りを自然に扱え、Weibullプロットでも打切り点を明示して視覚的に情報を活かす。逐次ベイズ更新を用いればフィールドデータで推定精度を高められる。

モデリング上の注意

3母数形で位置γを自由化すると尤度が平坦になり推定が不安定化しやすい。物理的根拠（最小作動時間、構造クリアランス、検出遅延など）を持つときに限定して採用するのがよい。環境や負荷の階層差が大きければ、共変量を導入した加速モデル（Arrhenius型、Eyring型、Cox比例ハザード）や層別解析を検討する。さらに、試験条件と実使用条件のギャップは加速係数やストレス強度モデルで橋渡しする。

ワイブル分布｜信頼性工学で使う代表分布