Наработка на отказ: формула расчёта MTBF

Наработка на отказ (MTBF, Mean Time Between Failures) - ключевая характеристика надёжности любого технического устройства: она показывает, сколько часов в среднем проработает изделие между последовательными отказами. Именно MTBF стоит в основе расчётов при проектировании самолётных систем, производственного оборудования, электронных блоков и медицинских приборов. Ниже разберём, откуда берётся формула, как вычислить вероятность безотказной работы за заданное время и как решить обратную задачу - найти требуемый MTBF по заданной надёжности. Чтобы сразу попробовать числа, воспользуйтесь калькулятором:

Что такое наработка на отказ и откуда берётся формула

Наработка на отказ - это математическое ожидание наработки (суммарного времени работы) изделия между двумя последовательными отказами. Если система отказывает случайно и вероятность отказа не зависит от возраста изделия, то закон распределения времени безотказной работы является показательным (экспоненциальным).

Вероятность того, что изделие проработает без отказа хотя бы $t$ часов, равна:

$P(t) = e^{-\lambda t} = e^{-t / \text{MTBF}},$

где $\lambda = 1 / \text{MTBF}$ - интенсивность отказов (число отказов в единицу времени, 1/ч). Этот закон обладает свойством отсутствия последействия: вероятность отказа в следующий час не зависит от того, сколько часов изделие уже проработало. В реальности это выполняется в «нормальной зоне» жизненного цикла - между приработкой и износом.

Рассчитать для своей задачи

Характерное значение: при $t = \text{MTBF}$ вероятность безотказной работы равна ровно $1/e \approx 0{,}368$. Иными словами, к моменту, когда наработка достигла среднего значения, больше 63% изделий уже откажет. Это противоречит интуиции («средняя наработка - значит, половина откажет»), но так работает экспоненциальный закон.

Формула и её параметры

Итоговая система формул для расчёта надёжности по экспоненциальному закону:

$\lambda = \frac{1}{\text{MTBF}},$

$P(t) = e^{-\lambda t} = e^{-t / \text{MTBF}},$

$\text{MTBF} = -\frac{t}{\ln P(t)}.$

Последняя формула - обратная задача: если известно, что за время $t$ вероятность безотказной работы составила $P$, то отсюда можно оценить MTBF. Обратите внимание: $\ln P$ при $P < 1$ отрицателен, поэтому MTBF всегда положителен.

Если за наблюдаемый период $T$ из $N$ изделий отказало $k$ штук, то оценка вероятности безотказной работы равна $\hat{P} = 1 - k/N$, и далее применяется та же обратная формула. Такой подход называют точечной оценкой MTBF по испытаниям; для получения доверительного интервала к нему прибавляют поправку на объём выборки, но в учебных задачах точечная оценка обычно достаточна.

Рассчитать для своей задачи

Пример решения: прямая задача

Разберём типовую задачу. MTBF электронного блока составляет 5000 ч. Требуется найти: (а) интенсивность отказов, (б) вероятность безотказной работы за 1000 ч, (в) ожидаемое число отказов в парке из 1000 изделий за тот же срок.

Шаг 1. Интенсивность отказов:

$\lambda = \frac{1}{5000} = 2{,}0 \cdot 10^{-4} \text{ 1/ч.}$

Шаг 2. Вероятность безотказной работы за 1000 ч:

$P(1000) = e^{-1000 / 5000} = e^{-0{,}2} \approx 0{,}819.$

Это 81,9%: около 18% изделий откажет за 1000 ч.

Шаг 3. Ожидаемое число отказов из 1000 изделий:

$N_\text{отк} = N \cdot (1 - P(t)) = 1000 \cdot (1 - 0{,}819) = 181 \text{ шт.}$

Проверка: $\lambda \cdot t = 0{,}0002 \cdot 1000 = 0{,}2$ - это «нагрузка» (число средних MTBF, вложенных в миссию). При нагрузке, значительно меньшей единицы, $P \approx 1 - \lambda t$, что даёт быструю линейную оценку: $\approx 20\%$ - грубо согласуется с точным значением 18,1%.

Пример решения: обратная задача

Задача: требуется, чтобы вероятность безотказной работы за 500 ч была не менее 0,9. Найти минимальный MTBF.

Применяем обратную формулу:

$\text{MTBF} = -\frac{t}{\ln P} = -\frac{500}{\ln 0{,}9}.$

Вычисляем: $\ln 0{,}9 \approx -0{,}10536$, значит:

$\text{MTBF} = \frac{500}{0{,}10536} \approx 4748 \text{ ч.}$

Полученный MTBF - минимальный: при более низком значении требование надёжности нарушается.

Интенсивность отказов:

$\lambda = \frac{1}{4748} \approx 2{,}11 \cdot 10^{-4} \text{ 1/ч.}$

Проверка: $P = e^{-500/4748} = e^{-0{,}1053} \approx 0{,}9$ - верно.

Связь MTBF с наработкой до отказа (MTTF)

В учебниках встречаются два близких термина: MTBF (Mean Time Between Failures) и MTTF (Mean Time To Failure). Разница принципиальна:

• MTBF используется для восстанавливаемых систем: это среднее время между двумя последовательными отказами, включает время восстановления (если оно пренебрежимо мало) или не включает, в зависимости от контекста.
• MTTF - для невосстанавливаемых изделий (лампочки, предохранители): среднее время до первого и единственного отказа.

При пренебрежимо малом времени восстановления MTBF ≈ MTTF, и формулы совпадают. В задачах ВУЗов чаще используют МТВF для любых изделий - просто как «среднее время между отказами».

Интенсивность отказов и кривая «ванны»

В реальности $\lambda(t)$ не постоянна: жизненный цикл изделия делится на три зоны. В начале - приработка (высокий начальный брак, $\lambda$ убывает), затем нормальная эксплуатация ($\lambda \approx \text{const}$, применима экспоненциальная модель), наконец износ ($\lambda$ растёт). Показательный закон с постоянной $\lambda$ - это именно модель нормальной эксплуатации.

Если задача относится к стадии износа или приработки, применяют распределение Вейбулла:

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

где $\beta$ - параметр формы ($\beta < 1$ - приработка, $\beta = 1$ - экспонента, $\beta > 1$ - износ), $\eta$ - характерная наработка. В большинстве учебных задач $\beta = 1$, то есть задачи сводятся к стандартной формуле MTBF. Чтобы проверить, применима ли экспоненциальная модель, строят вероятностную бумагу Вейбулла: если точки ложатся на прямую с наклоном примерно 45°, показательный закон адекватен и MTBF-расчёт даст корректный результат.

Частые ошибки

• Подстановка $\lambda$ в %/ч вместо 1/ч. Интенсивность отказов измеряется в 1/ч (или fit = $10^{-9}$ 1/ч). Если MTBF дан в часах, $\lambda = 1/\text{MTBF}$ - без пересчёта в проценты.
• Смешение MTBF и $T_{0,5}$ (медианной наработки). При MTBF = 5000 ч не означает, что половина изделий проработает 5000 ч. Медиана $T_{0,5} = \text{MTBF} \cdot \ln 2 \approx 0{,}693 \cdot \text{MTBF}$.
• Ошибка в знаке логарифма. В обратной формуле $\text{MTBF} = -t / \ln P$ логарифм отрицателен (так как $P < 1$), минус перед дробью делает результат положительным. Пропуск минуса даёт отрицательный MTBF.
• Использование $\log_{10}$ вместо $\ln$. Формула $P = e^{-t/\text{MTBF}}$ использует натуральный логарифм; замена на десятичный дает результат, заниженный в 2,303 раза.
• Забыть перевести годы в часы. MTBF обычно указывается в часах. Если время миссии задано в годах, нужно умножить на 8760 ч/год перед подстановкой.

FAQ

Что такое MTBF простыми словами и чему он равен? MTBF - среднее время между двумя последовательными отказами. Если MTBF = 10 000 ч, то в среднем изделие работает 10 000 ч до очередного отказа. При экспоненциальном законе именно при $t = \text{MTBF}$ вероятность безотказной работы равна $1/e \approx 36{,}8\%$.

Как найти MTBF по результатам испытаний? Если $k$ изделий отказало за суммарную наработку $T_\text{сум}$, то $\text{MTBF} = T_\text{сум} / k$. Для нескольких изделий суммарная наработка - это сумма наработок каждого до отказа (или до конца испытания). Затем из MTBF находят $\lambda = 1/\text{MTBF}$ и строят кривую надёжности.

Почему при $t = \text{MTBF}$ вероятность равна 1/e, а не 0,5? Потому что при показательном распределении математическое ожидание и параметр $1/\lambda$ совпадают, но медиана равна $\text{MTBF} \cdot \ln 2 \approx 0{,}693 \cdot \text{MTBF}$. Среднее и медиана расходятся из-за скошенности экспоненциального закона вправо: небольшая доля изделий работает очень долго, «вытягивая» среднее вверх.

Коротко

Наработка на отказ MTBF связана с вероятностью безотказной работы формулой $P(t) = e^{-t/\text{MTBF}}$; интенсивность отказов $\lambda = 1/\text{MTBF}$ постоянна в зоне нормальной эксплуатации. Обратная задача решается через $\text{MTBF} = -t / \ln P$. При $t = \text{MTBF}$ вероятность безотказной работы составляет $1/e \approx 36{,}8\%$, а медиана жизни изделия в 1,44 раза меньше MTBF.