マイクロ波無線装置の補修用シート保有に関する信頼性の検討

Transcription

1 マイクロ波無線装置の故障と信頼性維持について 補修用ユニットの保有に関する検討. 概説通信回線は設備が直列に接続されるため 一箇所でも切断すると系統全体が機能しなくなるという 信頼度上特異な構成となっている このため 冗長構成 ダイバシティー化 ルート化などの対策を施して所要の信頼度を確保しているが 設備の構成要素 ( ユニット ) が故障すると回線信頼度は直接又は間接的に脅かされることになる 特にマイクロ波無線装置のような多重通信機器が故障するとその影響は大きい このため送受信装置は冗長構成となっているが その片系が故障すると迅速に補修をしない限り所定の信頼度を確保できないことは当然である 最近の送受信装置は回路が複雑で現場での修理はほとんど不可能であり 故障ユニットをメーカに送って修理を依頼するが 近年この修理期間が問題となっている 中通における調査結果では 短いもので 56 日 長いものでは 3 日余りもの修理期間を要している この間 片系故障のまま放置することはできないので メーカから代替ユニットを借用して対策しているのが現状である このような折 近年における無線施設数の増大と機種の多様化により メーカ側も代替品の確保と供給に責任を持つことが困難になってきたとして あるメーカから代替ユニット供給に関する保守契約を結ぶよう申し入れがあった この問題に対処するためには 冗長系に故障が生じたとき 何時間以内に修理又は代替品を供給しなければならないかを技術的に検討する必要がある 回線切断に至る故障としては 共通部故障と冗長系二重故障があるが これらの故障に対する修理時間 (MTTR) についての検討を行った いくつかの条件を設定した上で信頼性理論に基づく計算を行ったが 結論的として 共通部故障のときは 直ちに修理を行う必要があり また冗長部故障の場合でも 速やかなる代替品の供給が不可欠であることを明らかにした 該当文書は 予備品等管理検討作業会 の資料として作成したものであるが これを技術研修資料として使用するために 次のような変更を加えた 文書のタイトルを マイクロ波無線装置の補修用シート保有に関する信頼性の検討 から本タイトルに変更した 概説 結論 及び 参考文献 を追加した 3 補修用シート を一般的な信頼性用語として 補修用ユニット に変更した 4 計算結果を若干修正したが 結論に影響するような修正ではない 5 本文にあった MTBF の区間推定 は付録とした また付録に 信頼性に関する基礎理論 を追加した 6 その他 研修資料として必要と思われる注書き等を追加した

2 . 適用理論. 信頼度関数と MTBF(Mean Time Between Failure) F t () 信頼度関数 R ( t) と不信頼度関数 ( ) 時間 ~t の間に故障する確率 R( t) と故障しない確率 F( t) ( t) + F ( t) = () 故障密度関数 f ( t) 時間 t における故障確率密度を f ( t) で表わす d F ( t) d R ( t) の間に次の関係が成立する R () f ( t) = = () dt dt ( 注 ) 確率密度関数と累積分布関数の累積分布関数は 次式で与えられる 一般に確率密度関数 f ( x) F x ( x) = f ( x)dx 確率密度関数 (Probability Denity Funtion) は P.D.F とも呼ばれ 関数記号は小文字で表わされる 累積分布関数 (Cumulative Ditribution Funtion) は C.D.F とも呼ばれ 関数記号は大文字で表わされる (3) 平均故障間隔 MTBF E t が MTBF となる f ( t) の期待値 ( ) E ( t) = MTBF = t f ( t) dt ( 注 ) 期待値の定義 変数 x の確率密度関数が f ( x) E ( x) x f ( x) + = dx 式 (3) に次の部分積分公式を適用する ( x) v ( x) dx = v ( x) u( x) u ( x) v( x)dx の場合 その期待値 (Expetation) は次式で与えられる u (4) ( x) = t, v ( x) = f ( t), u ( x) =, v( x) F ( t) u = F t f ( t) dt = tf( t) F( t) dt ( t) R( t) = であるから ( t) dt = tf( t) { R ( t) } dt = t F( t) t R ( t) dt t f + E ( t) = t f ( t) dt = [ t F( t) t] + R( t) dt ( t) lim F = t であるから [ F( t) t] = t [ F( t) ] t となる = (3)

3 E ( t) = MTBF = R( t) dt (5). 故障数の分布 () ポアソン分布 ある時間内に平均 µ 台故障する系において その時間内に x 台が故障する確率は一般に次 のポアソン分布に従うとされている x µ µ p( x) = e (6) x! 単一装置の故障率を とすると t 時間内の平均故障数は µ = t ( 台 ) となる 更に同一機種 が n 台ある場合の全体の故障数は p ( n t ) ( ) ( nt) x e x µ = n t となり ポアソン分布は次式となる = (7) x! 時間 t の間に故障しない確率は 式 (7) において x = とすれば得られ 指数分布となる ( ここで x! =! = である ) n t ( ) = e p (8) () 指数型信頼度関数式 (8) において 時間 t を変数にとると e n t は ~t 時間に故障しない確率 すなわち指数分布に従う信頼度関数を与える 指数型故障モデルは故障率が一定 ( 対象期間内に装置が劣化傾向を持たない ) モデルであり 本検討ではこのモデルを使用する 信頼度関数は次式となる R n t ( t) e = (9) 故障密度関数 f ( t) は式 () より dr ( ) ( t) n t f t = n e dt この分布の期待値 E( t) E = () が MTBF となる ( t) = MTBF = t f ( t) dt = t n e n t dt = n ().3 不通率 ( 保全度係数 ) 平均修復時間を MTTR (Mean Time to Repair) とし t 時間内に m 台の装置が故障したと仮定すると その間の不通時間は次式で与えられる 3

4 ( m) ( n t) m t Uhm = m MTTR p = m MTTR e () m! Uhm の期待値は その系の平均不通時間 Uhn となる p( ) + MTTR p( ) + + n MTTR p( m) + + n MTTR p( n) { p( ) + p( ) + np( n) } (3) Uhn = MTTR = MTTR + 式 (3) の { } 内はポアソン分布の期待値 n t となるので Uhn = n t MTTR (4) となる 指数型故障モデルでは = であるから t 時間内における不通時間は次式となる MTBF MTTR Uhn = n t MTTR = t (5) MTBF n t = とおくと 単位時間当たりの不通時間となり 保全度係数と呼ばれる ρ = MTTR n MTBF (6) n 不通率 ( 不稼働率 )U は次式で与えられるが 一般に MTTR MTBF であるから実際上 ρ が不通率となる MTTR MTTR U = = ρ n (7) MTBF MTBF + MTTR n n m n.4 システムの故障率 () 直列システム複数のブロックからなるシステムにおいて いずれかのブロックが故障すると全体の系がシステムダウンとなる系は信頼度的に直列なシステムと呼ばれる R ( t ) R ( t ) R 3 3 ( t ) R n( t ) n 図 直列システム 直列システムにおいては 各部分の信頼度関数の積が全体の信頼度関数となる 4

5 R ( t) R ( t) R ( t) R ( t) R ( t) = 3 n (8) 指数型故障モデルでは次に示すように 各部分の故障率の和が全体の故障率となる R t t 3t n t ( n) t ( t) = e e e e = e = e (9) = n () 総合の MTBF は次式となる = MTBF () () 並列システム並列システムにおいては 各部分の不信頼度の積が全体の不信頼度関数となる 図 に示すような二つのブロックからなる並列システムの総合 MTBF を求める F( t) = F ( t) F ( t) = { R ( t) }{ R ( t) } R ( t) = R ( t) = R ( t) = = の場合 F ( t) = { R ( t) } = R ( t) + R ( t) 総合信頼度関数 R ( t) は次式となる ( t) F( t) = R ( t) R ( t) R = () t ( t) = e R となるから R 総合 MTBF は t t ( t) = e e t t 3 MTBF = R( t) dt = = = e + e (3) ( t) R, ( t) R, 図 並列システム 式 (3) は修理を伴わない並列システムの MTBF であり 片系が故障しても修理をしないで放置しておくと 並列システムの MTBF は片系のみの場合の MTBF( = / ) の.5 5

6 倍にしかならないことになる (3) 修理を伴う並列システム単位時間内に修復できる確率を保全度 (M) という 信頼度関数 保全度関数とも指数分布 t µτ に従うものとし それぞれ R ( t) = e M ( τ ) e = とすると 状態推移確率行列を 解くことにより 結論として次式を得る [ 文献 [] p34] 3 + µ MTBF = (4) ここで µ は単位時間内に修復される確率であり 修復率と呼ばれる 式 (4) において µ = とおくと MTBF = 3 となり 式 (3) に示す修理をしない並列システムの MTBF と等しくなる 一般に µ であり 式 (4) の分子を + µ µ µ 3 とすると MTBF となる この 式は 比較的短期間に修理される並列システムの信頼度に関しては 故障率の乗法則が成立することを示唆している 片系故障時の保全度を µ 両系故障時の保全度を µ とし 両系故障による回線停止に対す る平均修復時間は MTTR = となるので 不通率は次式となる MTTR µ u = µ = = MTBF e 3 + µ µ ( 3 + µ ) (5) 式 (5) の MTBF は両系故障に対する等価平均故障間隔である e 以上は完全並列システムの場合であり 実際のシステムには共通部があるので これを考慮する必要がある ( 第 4 章参照 ) 3. 検討モデル 3. 回線系統モデル北向マイクロが 区間 局 西向マイクロが 区間 3 局であるが モデル回線として次の 区間 局の系統を想定する 3 区間 局 図 3 回線系統モデル 6

7 3. 計算モデル () 単純故障 MTBF(MTBF) 回線停止に至るかどうかに関係なく 修理を要する故障が発生した場合の平均故障間隔を単純 MTBF と呼ぶことにし これを MTBFとする この故障は発生頻度が高く データの取得が容易である () 等価 MTBF(MTBFe) 回線停止に直接結びつく故障 すなわち共通部故障 冗長系の同時故障などを対象とする平均呼称間隔を等価 MTBF をと呼ぶことにし これを MTBFe とする これは発生頻度が少なく 有意なデータの取得が事実上不可能である (3) 換算係数 ( k ) 単純 MTBF と等価 MTBF の比を換算係数 k とする MTBFe k = (6) MTBF マイクロ波無線装置の場合 k の値は約 5 とされているので本検討ではこの値を使用する 3.3 装置モデル () 検討条件 (a) 直列システム : 故障率の和とする 式 () 参照 (b) 並列システム : 等価 MTBF は式 (4) より次式となる 3 + µ MTBFe = (7) 送受信部の両系同時故障 ( 注 3) による等価故障率 は次式となる = = MTBFe 3 + µ TR (8) TR ( 注 3) 両系同時故障率装置の片系が故障し その修理が終わらないうちに他の片系が故障する確率を指す () 重故障までを対象とする () 送受信装置故障モデル装置の信頼度的構成を図 4のように設定する 電源部は便宜上共通部に含める 7

8 頼度的に並列共通部送信部受信部信 t t r r 信頼度的に直列 tr 図 4 送受信装置の故障モデル 4. 信頼度の計算 4. 冗長系を含む等価故障率計算の単純化のために 号 号の送信部 号 号の受信部の全ての故障率を等しい ものとし これを とする ' () 単純故障 + tr = (9) MTBF tr ' t + t + r + r 4 = (3) = () 回線停止故障 冗長系送受信部の両系同時故障率を とすると 回線停止故障率 は次式となる TR = + TR = (3) MTBFe 送受信部の故障率を等しいとして ' t t = r = r = トの両系同時故障率は 式 (4) より次式となる ' ' = とすると 送受信各ユニッ ' TR = (3) 3 + µ ここで µ は片系故障時の保全度である 8

9 (3) 計算上の故障率マイクロ波送受信装置の単純 MTBF は文献 [5] により約 万時間とされている 回線停止を対象とする等価 MTBF は単純 MTBF の約 /5 とされているのでこの値を使用する 送受信装置 台当りの故障率は次のように表される MTBF ' = = = (33) 5 回線停止対象故障率は次式となる 7 + TR = = = (34) 5 5 MTBF 5 一般に であるから TR 7 (35) 送受信各ユニット単体当りの故障率は であるから 4 ' 5 = = =.5 (36) 4 送受各ユニットの冗長部同時故障率は式 (3) より ' 5 ' (.5 ) = = ' µ µ TR (37) 4. 電源開発 無線装置の MTBF 及び MTTR 実績 () 運用実績と故障数電源開発 にて昭和 63 年度から平成 3 年度までの 4 年間について 稼動中のマイクロ波無線装置のコンポーネントアワーと故障数について調べた結果は表 の通りであった 表 無線装置の故障実績 年度コンポーネントアワー ( 注 4) 故障数 昭和 63 年度 74 台 876 時間 =,4,4 時間 平成 年度 7 台 876 時間 =,373,96 時間 平成 年度 75 台 876 時間 =,49, 時間 7 平成 3 年度 9 台 876 時間 =,549,6 時間 8 合計 9,73,36 時間 49 ( 注 4) コンポーネントアワーの算出に当たっては 年度途中で運用開始した装置についても年度初めから運用しているものとして計算した 9

10 () MTBF 点推定値 コンポーネントアワーがT = 時間で その間の故障数の合計が 49 件であったので MTBF の点推定値は以下の通りとなる MTBF = = 986 ( 時間 ) (38) 49 (3) MTBF 区間推定値 χ 分布による 9% 両側信頼区間を求める MTBF 区間推定表 ( 注 5) にて r = 5 から求めると次のようになる の欄 上限値 98,6.35=63, ( 時間 ) 下限値 98,6.768=53, ( 時間 ) (39) ( 注 5)MTBF 区間推定表式 (3.9) の係数.35 及び.768 は MTBF 区間推定表から求めた値である この表は例えば文献 [4] p9 に与えられている ( 付録 Ⅱ MTBF 区間推定 参照 ) 下限値と比較しても モデル計算における MTBF =. 時間よりは若干良い値となっている 4.3 片系故障時の信頼度低下特性 () モデルケースの場合 5 7 =.45 = MTBFe = 片系故障時の保全度を µ 両系故障時 の保全度を µ とすると 式 (3) より 無線装置 台当りの不通率は次式で求められる MTTR 4 u = = = + (4) MTBFe µ µ µ ( 3 + µ ) 前述の 4. (3) 計算上の故障率 において 無線装置 台を送受 現用予備の 4 台に分解したので 区間 (4 装置 ) 合計台数は 6 台となる なお 共通部故障率は無線装置 台当りの数値として扱う 従って合計不通率は次式となる u = u + u ( 合計故障率 ) r u = 4 ( 共通部不通率 ) µ (4)

11 ' u r = 6 ( 冗長部不通率 ) µ ' ( 3 + µ ) モデルケースの不通率について計算した結果を表 に示す 保全度については µ = µ = µ として計算した 片系故障時と両系故障時では緊急度の相違から修理時間に差があると考えられるが そのようなデータは得られていないので 同じとして計算した = 表 モデルケースの不通率 µ µ µ = 保全度 µ MTTR 時間 共通部不通率 冗長部不通率 () 実測 MTBF の場合電源開発 における MTBF 実績値の区間推定を求める 式 (39) に MTBF の区間推定値 が与えられているので MTBFe=MTBF 5 = = ' = と MTBF MTBFe 4 して故障率を求めると表 3のようになる 表 3 実測 MTBF による故障率推定値 MTBF MTBFe 単純故障率 共通部故障 ユニット単純故 ( 時間 ) ( 時間 ) 率 障率 ' 上限値 63, 下限値 53,

12 共通部不通率 = 4 冗長部不通率 = 6 µ µ ' ( 3 ' + µ ) として上限及び下限の不通率 を求めると表 4 のようになる 表 4 実測 MTBF による不通率 u MTTR( 時間 ) 上限共通部 上限冗長部下限共通部 下限冗長部 3.4E-6.888E- 5.4E E E-6.55E-9.5E E E-9.6E-5.56E E-8 5.4E E E-7.5E E E-7.6E-4.5E E-6 5.4E E E-5.5E E E-5.6E E E E E-.836 (3) MTTR の実績故障した送受信ユニットをメーカに修理に出して返ってくるまでの日数を 平成 年 4 月から平成 4 年 月までの 4 件について調べた結果によると 最短日数が 7 日 最長日数が 3 日 平均は 56 日であった 別紙の 装置故障に伴う不通率推定曲線 にはこのデータを記入している (4) 目標信頼度モデル回線に対する目標信頼度としては 文献 [5] で示されている無線装置 台当りの信頼度 - 6 を採用した 回線系統の目標信頼度は 4 台 - 6 = 8-5 となる 別紙の目標信頼度にはこの値を与えている 4.4 片系故障時の伝搬路信頼度低下特性 6QAM 方式の無線装置は 号 号受信機ともダイバシティー用受信機を備えているが 4PSK 方式の無線装置は 号受信機を SD 用受信機として使用しているため 号のいずれかが故障して片肺運用になるとダイバシティー効果がなくなる この影響による信頼度の低下はさほど大きいとは考えられないが その概要を把握するために以下の検討を

13 行う ここでは計算方法のみを示し 具体的計算は行っていない () 検討モデル 中継区関数 : 区間 受信装置の単純故障率 : 台当り.5-6 式(3.7) より 回線系統の伝搬路信頼度基準値 : 5-5 ( 電協研 6kmモデルより ) 故障した区間には改善度 A ( 倍 ) のダイバシティーが採用されているものとする () 計算方法次の二つの計算方法を示す 年間にある 区間が故障した場合の条件付確率で瞬断率を求める 全受信装置 n 台のうち r 台が故障し その故障区間でレーレーフェージングが発生したとして 瞬断率を求める (3) 条件付確率による瞬断率の求め方 ダイバシティーのある 区間の信頼度 = 5-5 /=.5-6 片系が故障した 区間の信頼度 =.5-6 A A はダイバシティー改善度 ( 倍 ) である ( 注 ) ダイバシティー効果が無くなった 区間の伝搬路信頼度ダイバシティー区間は 所要伝搬路信頼度.5-6 をダイバシティーにより満足させているはずである 従って その 区間の伝搬路信頼度に関して次の関係が成立する P R [ ] A M =.5 ここで P : レーレーフェージング発生確率 R A : ダイバシテイ改善度 ( 倍 ) [ M ]: フェージングマージン ( 真値 ) ダイバシティー効果がなくなると その区間の伝搬路信頼度は次式となる P R =.5 A [ M ] 年間瞬断時間 Dy は以下のとおりになる Dy =.5 A MTTR +.5 ( 876 MTTR) ( 時間 ) 故障した 区間の修 故障した 区間のその 故障していない 9 区間 理中の断時間 他の断時間 の瞬断時間 3

14 瞬断率 = Dy/ 876 (4) 上の式にダイバシティー改善度 A ( 例えば 倍 ) 及び MTTR の値を入れれば瞬断率を計 算することができる (4) 故障した r 台の受信区間でフェージングが発生した場合の瞬断率 受信装置台数 : n = 区間 4 台 = 8 台 台の故障率が p =.5 の場合に k 台が同時に故障する確率 p k は次の二項分布で与えれる n k n k p =! k p ( p) k! n k! (43) ( ) k 台の受信装置が故障し その修理期間中にフェージングが発生して瞬断が生じる確率は次式で求められる 瞬断率 = p.5 + p {(.5 A MTTR )/ ( 876 MTTR) 876} / p {(.5 A MTTR )/ ( 876 MTTR) / 876} p n {(.5 A MTTR) / n.5 ( 876 MTTR) / 876} (44) 上の式にダイバシティー改善度 A ( 例えば 倍 ) 及び MTTR の値を入れれば瞬断率を計算することができる 以上は計算方法を示しただけで 具体的には算出していないが 片系の受信装置が故障し その区間のダイバシティーが不能となっても系全体に与える影響は少ない この問題は 海上区間など ダイバシティー機能が不可欠な区間について個別に検討すべきであると考える 5. 結論無線装置の代替ユニットに関する保守契約検討時の資料を基に研修用資料を作成した 研修目的であるから全面的に作り変えることも可能であったが 通信設備保全に関する工学的アプローチの具体例として なるべく元の文書の内容を残すよう考慮した 4

15 本検討による計算結果は別紙のグラフに示してある モデル回線の共通部不通率カーブによると 信頼度目標を満たすためにはMTTR= すなわち 時間以内に修理を完了しなければならない事を示している これは年間の許容不通時間を示しているのではなく 共通部故障率に基づく故障が発生した場合にこの時間内に修理を施す必要があることを示している 共通部故障率 台は モデル回線において 4 年に 回の割合で発生する確率である モデル回線の冗長部不通率では 目標信頼度を満たす MTTR は 時間ほどになっている 中通の予備ユニットに対する MTTR 実績はこの期間を大きく超過しており 代替ユニットの供給が不可欠であることを示している グラフには実測 MTBF に基づく不通率の上下限カーブも点線で示した 実測 MTBF は文献 [5] の調査結果よりは少し向上しているが 昭和 5 年頃のアナログ無線装置に比べて大きく改善されたとは言いがたい 元の文書は 代替ユニットの必要性を概略的に示すのが目的であり 条件設定等はあまり厳密ではない 従ってこの検討結果を他の目的に利用する場合は注意を要する 例えば 故障が発生した場合の駆けつけ時間に関する検討が目的であれば それに合わせたモデルと条件設定を行う必要があり 当然別の方法論を考える必要がある このような理由から この資料は研修目的のみに使用することを希望する 参考文献 [] 佐藤喜代蔵 信頼性工学 オーム社 [] 塩見弘 信頼性の基礎 電気 電子工学体系 コロナ社 [3] 塩見弘著 高木昇監修 信頼性概論 信頼性工学体系 東京電機大学 [4] 塩見弘著 高木昇監修 信頼性の基礎数学 信頼性工学体系 東京電機大学 [5] 電気協同研究第 33 巻第 号 電力用マイクロ波無線装置信頼度向上対策 電力用マイクロ波無線装置信頼度向上専門委員会 5

16 付録 Ⅰ MTBF の信頼区間. 少数サンプルの MTBF 少ない故障数から求めた実測 MTBF を母集団の MTBF の推定値と考えるとき その推定 誤差が問題となる 以下には χ ( カイ二乗 ) 分布を使用する MTBF の信頼区間の求め方 を述べる ( 付録 Ⅱ MTBF の区間推定 参照 ) ( 注 ) 参考文献 文献 [3]p3 p89 (MTBF の区間推定 ) 文献 []p45( ポアソン分布の部分和と χ 分布の関係 ) ( 注 ) 信頼区間 統計標本から母数を推定するとき ( 例えば 標本平均から母平均を推定するとき ) 母数が例えば 9% の確率である区間内に存在すると推定されたとき これを母数の 9% 信頼区間という. ポアソン分布による信頼区間指数型故障モデルにおいて 一定時間における平均故障数を m とすると 故障数 r の分布は次のポアソン分布となる f ( r) r m m = e (a) r! 累積分布は次式となる F ( r) = f ( r) r= (a) 故障数の実測値が r であった場合 母平均 m の推定モデルは付録図 のようになる p p m L : 母平均下限推定値 m U : 母平均上限推定値 p : 危険率 m L r m U 付録図 故障数の区間推定図 ( 注 5) ポアソン分布は離散型分布であるが 付録図 では便宜上連続型分布で図示している

17 付録図 から m, m を求めることは 付録図 に示すように 平均 r のポアソン分布の p L U p 及び の値を求めることと同じである p p m L r m U 付録図 平均 r の分布からの推定 MTBF の推定値は次式から求めることができる MTBF ( 上限値 ) = T m L T MTBF ( 下限値 ) = T : コンポーネントアワー m U (a3) ( 注 ) コンポーネントアワー (Component hour) 系 機器 部品などについて 測定された個々の動作時間のまたは試験時間の総計値をいう (JIS Z 85 信頼性用語 より) ただし ポアソン分布は離散型分布であるため 故障数は整数であり 丁度 p に相当する 位置に整数値がくるとは限らない 故障数 r が大きい場合 ( r > ) は 正規近似によって求 めることができるが 少ない場合は誤差が大きくなる 少数故障を含む MTBF の区間推定法としては次の χ ( カイ二乗 ) 分布による推定法が適用される 3. χ 分布による MTBF の区間推定 () ポアソン分布の部分和と χ 分布の関係 + x x 不定積分の漸化式 e x dx = e ( x + x +! ) を使うと ポアソン分布の部分和 は次式で表される F r m m m ( ) = e = ( m + m + +! ) = r= r!! Γ m ( + ) e x x dx (a4)

18 Γ + = ただし ( )! 一方 自由度 φ の χ 分布は次式で与えられる 密度関数 ( χ ) 分布関数 F( ) φ χ χ φ = e ( ) (a5) φ Γ χ / φ ( ) x φ / χ = e x dx (a6) Γ 次に F( χ ) の補数 F( χ ) ( ) を確率 p に等しくするような χ の値を考える φ x F χ = e x dx = p (a8) φ Γ( ) χ 上式と式 (a4) を比較すると χ = m φ = ( +) とすると両式が一致することがわか る すなわち となる m χ φ, p χ 分布の (%) = m = χ φ, p p 値を χ φ, p とすると () MTBF の区間推定 ( 定時打切方式 ) 前述のように観測値が r の場合の故障数の上 下限値は r を平均値とするポアソン分布の p (%) 及び ( p) (%) に等しい 従って 平均 m のポアソン分布の (%) p, とすると m =, p = χ { ( r + ) p} p を m p χ φ となる p p の値を MTBF の下限に の値を上限対応させると 定時打切方式の場合の MTBF の区間推定値は次式となる T T MTBF ( 下限 ) = = p p χ ( r + ), χ ( r + ), (a9) T T MTBF ( 上限 ) = = p p χ r, χ r, 上限と下限で自由度 ( 注 ) が異なるのは 定時打切方式の場合における打切り時点の想定 の相違によるもので 打切られた時点が次の故障の直前であったとすると自由度は ( r +) 故障の直後であったとすると r となる 定数打切方式の場合は 故障直後に打切るので 上下限とも自由度は r となる 3

19 ( 注 ) χ 分布の自由度と上側確率 χ φ, p χ 分布の変数が ( ) で表されるとき φ は自由度であり p は上側確率である χ 分布の確率密度関数を f ( χ ) 累積分布関数を F( χ ) 次の図のようになる とすると これらの関係は f ( χ ) ( ) p = F χ ( ) p = F χ χ 分布関係図 ( p) χ φ, χ 以上 4

20 付録 Ⅱ 信頼性に関する基礎理論 この資料は電源開発 殿より開発電子技術 が受託した 通信機器劣化診断に関する技術調査業務 の報告書( 平成 年 3 月 ) の第 3 章を抜粋したものであり 電発殿の許可を得て研修資料として使用するものである 3. 信頼性の基礎 3-3. ワイブル確率紙 修理を施されるシステムの信頼性 システムの信頼性 MTBF の区間推定 3-6 参考文献 3-