故障率のソースを表示

'''故障率'''（こしょうりつ、{{Lang-en-short|failure rate}}）とは、システムや部品が故障する頻度で、単位時間当たりの故障数で表される。通常、[[ギリシア文字]]の [[λ]]（ラムダ）で表され、[[信頼性工学]]でよく用いられる。

故障率は通常、時間に依存している。例えば、自動車の5年目の故障率は1年目のそれの何倍にもなる可能性がある。

保険、金融、商業においても安全なシステムを設計するための基本的なパラメータである。

1年間の使用でシステムや部品が故障する確率を特に「年間故障率」と表現することがある。

== 平均故障間隔 ==
{{See Also|平均故障間隔}}
故障率の代わりに[[平均故障間隔]]（MTBF、1/λ）が報告されることがよくある。MTBFは、故障率が一定の場合に有効であり、複雑なシステムや電子機器によく用いられる。一部の信頼性基準（軍事および航空宇宙）では一般的となっている。故障率が一定ということは、[[バスタブ曲線]]の平坦な領域、つまり摩耗型の故障が起こる前の時間領域での寿命を示す。このため、MTBFから耐用年数を推定することは正しくない。バスタブ曲線の「寿命末期の磨耗型」の故障率ははるかに高いため、一般的に耐用年数はMTBFよりもはるかに短寿命になるためである。

MTBFが故障率より好まれる理由は、大きな正の数値（2000時間など）の方が、非常に小さな数値（1時間あたり0.0005など）よりも直感的で覚えやすいためである。

MTBFは故障率を管理する必要があるシステム、特に安全系において重要なパラメータである。設計要件に頻繁に登場し、必要なシステムの保守・点検の頻度を定めるのに役立つ。故障からの回復時間が無視でき、かつ故障率が時間に対して一定である場合、MTBFは故障率の逆数（1/λ）になる。

[[運輸業|運輸業界]]、特に[[鉄道]]や[[貨物自動車|トラック]]輸送で使われている類似の指標に「平均故障距離間隔」（mean distance between failures、MDBF）がある。

== 故障率の推定 ==
故障率はいくつかの方法で推定できる。故障時間の母集団を正しく推定できれば良い。一般的な手段は次のとおりである。

=== 寿命試験 ===
最も正確な方法は寿命試験（life testing）であり、機器やシステムのサンプルをテストして実際の故障データを得る。故障時間が全くばらつかない場合、1サンプルのみを試験すれば良いが、現実には個体毎にばらつく場合が多い。このばらつきを把握する（故障時間の母集団を推定する）ため、サンプル規模は大きくながちで非常に費用がかかる。また、故障時間が非常に長い場合、現実的な時間で試験を終了させることも困難になる。

=== 市場実績 ===
市場で実際に発生した故障データから、統計分析で故障率を推定することができる。正確な故障率を得るためには、分析者は機器の動作、データ収集の手順、故障率に影響を与える主要な環境変数、システムレベルでの機器、部品の使用方法などを十分に理解している必要がある。

=== 生産履歴 ===
多くの組織では、製造している機器やシステムの故障情報を内部データベースとして保持しており、故障率を算出するために使用することができる。新しい機器やシステムでデータが少ない場合には、過去のデータが有用な事前情報として役立つ。

=== 政府および商用の故障率データ ===
さまざまな部品の故障率のハンドブックが政府や民間から入手できる。「電子機器の信頼度予測（''Reliability Prediction of Electronic Equipment''）」（MIL-HDBK-217F）は、多くの軍用電子部品の故障率データを提供する[[軍用規格]]である。また、非電子部品を含む商用部品に焦点を当てたいくつかの故障率データソースが市販されている。

=== 予測法 ===
故障率推定の重大な欠点の1つに時間がかかることがある。故障率データが集まる頃には、対象となる機器が旧式になっていることがよくある。これを克服するために以下のサイクル試験などが故障率予測法として開発されている。

=== サイクル試験 ===
機械的な動作は機器の摩耗を引き起こし、主な故障メカニズムとなる。この消耗故障点は、サイクル試験（cycle testing）により、機器が故障するまでに実行されたサイクル数で定量化される。サイクル試験は機器が故障するまで可能な限り迅速に繰り返される。複数のサンプルを試験し、例えばその内の10%の個体が故障するまで試験が行われる。
<!-- 英語版記事の「FMEDA」の節は翻訳から外した。（理由：商品宣伝であると考えた） -->

== 離散的な故障率 ==
故障率は次のように定義される：

{{Indent|規定の条件下における特定の測定間隔の間に、あるアイテムの母集団内で起こった故障総数を、その集団が費やした合計時間で割ったもの。(MacDiarmid, et al.)}}

故障率 <math>\lambda (t)</math> は、しばしば時間 <math>t</math> 以前に故障がない場合に特定の間隔で故障が発生する[[確率]]と考えられがちだが、1を超えることもあるので実際には確率ではない。故障率を誤ってパーセント（%）で表現すると、特に修理可能なシステム、故障率一定型でないシステム、または動作時間が異なる複数のシステムについて測定する場合に、この尺度を正しく認識できない可能性がある。故障率 <math>\lambda (t)</math> は、[[信頼性関数]] <math>R(t)=1-F(t)</math> （時刻 <math>t</math> 以前に故障が発生しない確率。'''生存関数'''とも呼ばれる）を用いて次のように定義される。

:: <math>\lambda(t) = \frac{f(t)}{R(t)}</math>,

ここに <math>f(t)</math> は（最初の）故障までの時間分布（すなわち故障密度関数）である。<math>t_1</math>（または <math>t</math>）から <math>t_2</math> までの時間区間 <math>\Delta t</math> = <math>(t_2-t_1)</math> において、

:: <math>\lambda(t) = \frac{R(t_1)-R(t_2)}{(t_2-t_1) \cdot R(t_1)}
 = \frac{R(t)-R(t+\Delta t)}{\Delta t \cdot R(t)} \!</math>

ただし、これは[[条件付き確率]]であり、時間 <math>t</math> 以前に故障が発生していないことが条件である。そのため、分母には <math>R(t)</math> が含まれている。

ハザード率（後述）とROCOF（rate of occurrence of failures、故障発生率）は、しばしば故障率と同じものと誤解されることがある{{clarify|date=April 2015}}。違いを明確にするならば、アイテムの修理が早ければ早いほど、またすぐに壊れるので、ROCOFは高くなる。しかし、ハザード率は、修復時間や物流遅延時間には依存しない。

== 連続的な故障率 ==
[[ファイル:Loglogistichaz.svg|右|サムネイル|300x300ピクセル|ハザード関数 <math>h(t)</math> を{{仮リンク|対数ロジスティック分布|en|Log-logistic distribution}}を選択してプロットした。]]
故障率をより小さな時間間隔で計算すると、'''{{visible anchor|ハザード関数}}'''（hazard function、'''ハザード率'''（hazard rate）とも呼ばれる） <math>h(t)</math> が得られる。これは、<math>\Delta t </math> がゼロに近づくにつれて、瞬間故障率（instantaneous failure rate）あるいは瞬間ハザード率（instantaneous hazard rate）と呼ばれるものになる。

: <math>h(t)=\lim_{\Delta t \to 0} \frac{R(t)-R(t+\Delta t)}{\Delta t \cdot R(t)}.</math>

連続故障率は、時刻 ''<math>t</math>'' まで（少なくとも）の故障確率 <math>\operatorname{Pr}</math> を表す[[累積分布関数|累積分布型]]の'''故障分布関数''' <math>F(t)</math> に依存しており、

: <math>\operatorname{Pr}(T\le t)=F(t)=1-R(t),\quad t\ge 0 \!</math>

と表される。ここに <math>{T}</math> は故障時間である。

この故障分布関数は、故障密度関数 <math>f(t)</math> の積分で、

: <math>F(t)=\int_{0}^{t} f(\tau)\, d\tau \!</math>

である。これによりハザード関数は、

: <math>h(t)=\frac{f(t)}{1-F(t)}=\frac{f(t)}{R(t)}</math>

と定義できる。

[[ファイル:Exponential_pdf.svg|右|サムネイル|300x300ピクセル|指数型の故障密度関数。これらはそれぞれ（異なる）一定型のハザード関数を持っている（本文参照）。]]
故障分布のモデル化においては、多くの確率分布を用いることができる（{{仮リンク|確率分布のリスト|en|List of probability distributions}}を参照）。

一般的なモデルは、[[指数分布|指数密度関数]]に基づく'''指数故障分布'''、

: <math>F(t)=\int_{0}^{t} \lambda e^{-\lambda \tau}\, d\tau = 1 - e^{-\lambda t} \!</math>

である。

これに対するハザード率関数は、

: <math>h(t) = \frac{f(t)}{R(t)} = \frac{\lambda e^{-\lambda t}}{e^{-\lambda t}} = \lambda </math>

である。このように、指数故障分布では、ハザード率は時間に対して一定である（つまり「{{仮リンク|無記憶性|en|Memorylessness}}」分布）。[[ワイブル分布]]や[[対数正規分布]]のような他の分布では、ハザード関数は時間に対して一定ではない場合がある。[[退化分布|確定的分布]]などの一部では[[単調増加]]であり（「[[摩耗]]」に類似）、[[パレート分布]]などの他の分布では単調減少であるが（「{{仮リンク|バーンイン|en|Burn-in}}」に類似）、多くの場合は単調ではない。

微分方程式

: <math>h(t)=\frac{f(t)}{1-F(t)}=\frac{F'(t)}{1-F(t)}</math>

を <math>F(t)</math> について解くと、

: <math>F(t) = 1 - \exp{\left(-\int_0^t h(t) dt \right)}</math>

であることがわかる。

== 故障率減少型 ==
{{Main|故障率曲線}}
故障率減少型（decreasing failure rate、DFR）とは、ある事象が将来の一定の時間間隔で発生する確率が、時間の経過とともに減少していく現象を表す。故障率減少型は、初期に起こる故障が解消または修正される「初期故障期間」を表すことができ<ref>{{Cite book|doi=10.1007/978-1-84800-986-8_1|chapter=Introduction|first=Maxim|last=Finkelstein|title=Failure Rate Modelling for Reliability and Risk|series=Springer Series in Reliability Engineering|pages=1–84|year=2008|isbn=978-1-84800-985-1}}</ref>、λ(t) が[[減少|減少関数]]である状況に対応する。

DFR変数の混合<!-- mixtures -->はDFRである<ref name="brown1980">{{Cite journal|last1=Brown|first1=M.|year=1980|title=Bounds, Inequalities, and Monotonicity Properties for Some Specialized Renewal Processes|journal=The Annals of Probability|volume=8|issue=2|pages=227–240|doi=10.1214/aop/1176994773|jstor=2243267|doi-access=free}}</ref>。[[指数分布]]確率変数の混合は、{{仮リンク|超指数分布|en|Hyperexponential distribution}}である。

=== 再生過程 ===
DFR再生関数を持った{{仮リンク|再生理論|en|Renewal theory|label=再生過程}}（renewal processes）では、再生間時間<!--inter-renewal times-->は凹になる<ref name="brown1980" /><ref name="shanthikumar">{{Cite journal|last1=Shanthikumar|first1=J. G.|year=1988|title=DFR Property of First-Passage Times and its Preservation Under Geometric Compounding|journal=The Annals of Probability|volume=16|issue=1|pages=397–406|doi=10.1214/aop/1176991910|jstor=2243910|doi-access=free}}</ref>{{訳語疑問点|date=2021年10月}}。Brownは逆に、再生間時間が凹になるためにはDFRが必要であると推測したが<ref>{{Cite journal|last1=Brown|first1=M.|year=1981|title=Further Monotonicity Properties for Specialized Renewal Processes|journal=The Annals of Probability|volume=9|issue=5|pages=891–895|doi=10.1214/aop/1176994317|jstor=2243747|doi-access=free}}</ref>、この推測は離散的な場合<ref name="shanthikumar" />にも連続的な場合にも成り立たないことが示されている<ref>{{Cite journal|last1=Yu|first1=Y.|year=2011|title=Concave renewal functions do not imply DFR interrenewal times|journal=Journal of Applied Probability|volume=48|issue=2|pages=583–588|arxiv=1009.2463|doi=10.1239/jap/1308662647}}</ref>。

=== 用途 ===
故障率増加型（Increasing failure rate、IFR）は、部品が消耗することによって起こる直感的な概念である。故障率減少型（DFR）は、経年変化によって改善されるシステムを表す<ref name="proschan" />。宇宙船の寿命においても故障率減少型が見られ、Baker and Bakerは『最後に残ったこれらの宇宙船は、延々と続く。』（"These spacecraft that last, last on and on."）とコメントしている<ref>{{Cite journal|last1=Baker|first1=J. C.|last2=Baker|first2=G. A. S. .|year=1980|title=Impact of the space environment on spacecraft lifetimes|journal=Journal of Spacecraft and Rockets|volume=17|issue=5|pages=479|bibcode=1980JSpRo..17..479B|doi=10.2514/3.28040}}</ref><ref>{{Cite book|doi=10.1002/9781119994077.ch1|chapter=On Time, Reliability, and Spacecraft|first1=Joseph Homer|last1=Saleh|first2=Jean-François|last2=Castet|title=Spacecraft Reliability and Multi-State Failures|pages=1|year=2011|isbn=9781119994077}}</ref>。航空機の空調システムの信頼性は、個々に[[指数分布]]を持つことが分かっており、プールされた母集団ではDFRとなる<ref name="proschan">{{Cite journal|last1=Proschan|first1=F.|year=1963|title=Theoretical Explanation of Observed Decreasing Failure Rate|journal=Technometrics|volume=5|issue=3|pages=375–383|doi=10.1080/00401706.1963.10490105|jstor=1266340}}</ref>。

=== 変動係数 ===
故障率が減少している場合、[[変動係数]] ⩾&nbsp;1 であり、故障率が増加している場合、変動係数 ⩽&nbsp;1 である<ref>{{Cite journal|last1=Wierman|first1=A.|last2=Bansal|first2=N.|last3=Harchol-Balter|first3=M.|year=2004|title=A note on comparing response times in the M/GI/1/FB and M/GI/1/PS queues|url=http://users.cms.caltech.edu/~adamw/papers/fbnote.pdf|journal=Operations Research Letters|volume=32|pages=73–76|doi=10.1016/S0167-6377(03)00061-0|author-link1=Adam Wierman|author3-link=Mor Harchol-Balter}}</ref>。この結果は、すべての t&nbsp;⩾&nbsp;0 <ref>{{cite book|title=Analysis of Queues: Methods and Applications|first=Natarajan|last=Gautam|publisher=CRC Press|year=2012|page=703|isbn=978-1439806586}}</ref>に対して故障率が定義されている場合にのみ成立し、逆の結果（変動係数が故障率の性質を決める）は成立しない。

=== 単位 ===
故障率は任意の時間の尺度で表すことができるが、実際には'''[[時間 (単位)|時間]]'''（hours）が最も一般的な単位である。時間の代わりに、マイルや回転数などの他の単位を使用することもできる。

故障率は非常に低いことが多いため、特に個々の部品の故障では、100万個あたりの故障数または{{仮リンク|工学的記数法|en|Engineering notation|label=工学表記}}（10<sup>-6</sup>）<!-- 10E-6 -->で表されることがよくある。

[[電子部品|デバイス]]の'''故障率'''（Failures In Time、'''FIT'''）は、10億（10<sup>9</sup>）<!-- 10E+9 -->デバイス時間の動作で予想される故障の数である<ref>
Xin Li; Michael C. Huang; Kai Shen; Lingkun Chu.
[http://www.cs.rochester.edu/~kshen/papers/usenix2010-li.pdf "A Realistic Evaluation of Memory Hardware Errors and Software System Susceptibility"].
2010.
p. 6.
</ref>。たとえば、1000個のデバイスで100万時間や、100万個のデバイスで1000時間など、さまざまな組み合わせが考えられる。この用語は、特に[[半導体]]業界で使用されている。

FITとMTBFの関係は次のように表される。

: MTBF = 1,000,000,000 × 1/FIT

=== 加法性 ===
ある種の工学的な仮定（たとえば、故障率一定型に関する上記の仮定に加えて、考慮されるシステムには関連する[[冗長性]]がないという仮定）の下では、複雑な[[システム]]の故障率は、単位が一貫している限り、その構成要素の個々の故障率の単純な合計となる（たとえば、100万時間当たりの故障数）。これにより、個々の構成要素またはサブシステムをテストすることが可能になり、それらの故障率を加算してシステム全体の故障率を得ることができる<ref>
[http://www.weibull.com/hotwire/issue108/relbasics108.htm "Reliability Basics"].
2010.
</ref><ref>
Vita Faraci.
[http://src.alionscience.com/pdf/1Q2006.pdf "Calculating Failure Rates of Series/Parallel Networks"].
2006.
</ref>。

[[単一故障点]]をなくすために「冗長」部品を追加すると、ミッション故障率<!-- mission failure rate -->は改善するが、直列故障率<!-- series failure rate -->（ロジスティクス故障率とも呼ばれる）は悪化する。つまり、追加の部品は平均重大故障間隔（mean time between critical failures、MTBCF）を改善させる反面、何かが故障するまでの平均時間は悪化する<ref>
[https://www.quanterion.com/mission-reliability-and-logistics-reliability-a-design-paradox/ "Mission Reliability and Logistics Reliability: A Design Paradox"].
</ref>。

=== 計算例 ===
特定の部品の故障率を推定する必要があると仮定する。その故障率を推定するために、テストを行うことができる。同一の部品10個は、それぞれが故障するか1000時間に達するまでテストされ、その時点でその部品のテストを終了する。(この例では、統計的な[[信頼区間]]は考慮されていない)。その結果は次のとおりである。

推定故障率は、

: <math>\frac{6\text{ failures}}{7502\text{ hours}} = 0.0007998\, \frac{\text{failures}}{\text{hour}} = 799.8 \times 10^{-6}\, \frac{\text{failures}}{\text{hour}}, </math>

となり、100万時間稼働するごとに799.8件の故障が発生する。

== 参照項目 ==
{{Portal|数学}}
{{Div col|colwidth=22em}}
* {{ill2|年間故障率|en|Annualized failure rate}}
* {{ill2|バーンイン|en|Burn-in}}
* [[故障]] [[:en:Failure|Failure]]
* [[故障モード]]
* [[死力]]
* {{ill2|超過頻度|en|Frequency of exceedance}}
* [[信頼性工学]]
* [[生存分析]]
* [[ワイブル分布]]
* [[故障率曲線]]
{{div col end}}

== 脚注 ==
{{Reflist|30em}}

== 推薦文献 ==

* {{Citation|last=Goble|first=William M.|year=2018|title=Safety Instrumented System Design: Techniques and Design Verification|publisher=International Society of Automation|location=Research Triangle Park, NC 27709}}
* {{cite book|last=Blanchard|first=Benjamin S.|year=1992|title=Logistics Engineering and Management|edition=Fourth|pages=26–32|publisher=Prentice-Hall|location=Englewood Cliffs, New Jersey|isbn=0135241170}}
* {{cite book|last=Ebeling|first=Charles E.|year=1997|title=An Introduction to Reliability and Maintainability Engineering|pages=23–32|publisher=McGraw-Hill|location=Boston|isbn=0070188521}}
* [[:en:Federal Standard 1037C|Federal Standard 1037C]]
* {{cite book|last=Kapur|first=K. C.|last2=Lamberson|first2=L. R.|year=1977|title=Reliability in Engineering Design|pages=8–30|publisher=John Wiley & Sons|location=New York|isbn=0471511919}}
* {{cite journal|last=Knowles|first=D. I.|year=1995|title=Should We Move Away From 'Acceptable Failure Rate'?|journal=Communications in Reliability Maintainability and Supportability|volume=2|issue=1|page=23|publisher=International RMS Committee, USA}}
* {{cite book|last=MacDiarmid|first=Preston|last2=Morris|first2=Seymour|date=n.d.|title=Reliability Toolkit|edition=Commercial Practices|pages=35–39|publisher=Reliability Analysis Center and Rome Laboratory|location=Rome, New York}}
* {{cite book|last=Modarres|first=M.|author-link=Mohammad Modarres|last2=Kaminskiy|first2=M.|last3=Krivtsov|first3=V.|year=2010|title=[[:en:Reliability Engineering and Risk Analysis: A Practical Guide|Reliability Engineering and Risk Analysis: A Practical Guide]]|edition=2nd|publisher=CRC Press|isbn=9780849392474}}
* {{cite journal|last=Mondro|first=Mitchell J.|date=June 2002|title=Approximation of Mean Time Between Failure When a System has Periodic Maintenance|url=http://www.mitre.org/work/best_papers/02/mondro_approx/mondro_approx.pdf|journal=IEEE Transactions on Reliability|volume=51|issue=2|pages=166–167|doi=10.1109/TR.2002.1011521}}
* {{cite book|last=Rausand|first=M.|last2=Hoyland|first2=A.|year=2004|title=System Reliability Theory; Models, Statistical methods, and Applications|publisher=John Wiley & Sons|location=New York|isbn=047147133X}}
* {{cite book|last=Turner|first=T.|last2=Hockley|first2=C.|last3=Burdaky|first3=R.|year=1997|title=The Customer Needs A Maintenance-Free Operating Period|work=1997 Avionics Conference and Exhibition, No. 97-0819, P. 2.2|publisher=[[:en:ERA Technology Ltd|ERA Technology Ltd]].|location=Leatherhead, Surrey, UK}}
* U.S. Department of Defense, (1991) ''Military Handbook, “Reliability Prediction of Electronic Equipment, MIL-HDBK-217F, 2''

== 外部リンク ==
* [http://www.asq.org/reliability/articles/bathtub.html Bathtub curve issues], ASQC
* [http://lamspeople.epfl.ch/kirrmann/Pubs/FaultTolerance/Fault_Tolerance_Tutorial_HK.pdf ''Fault Tolerant Computing in Industrial Automation''] by Hubert Kirrmann, ABB Research Center, Switzerland

{{Statistics|analysis}}

{{DEFAULTSORT:こしようりつ}}
[[Category:保険数理]]
[[Category:工学の失敗]]
[[Category:信頼性工学]]
[[Category:生存分析]]
[[Category:保守]]
[[Category:統計比率]]
[[Category:誤差測度]]
[[Category:率・割合]]