Интенсивность отказов: расчёт и кривая ванны

Интенсивность отказов $\lambda(t)$ - центральная характеристика теории надёжности: она показывает, насколько быстро изделие «выходит из строя» в каждый конкретный момент наработки. Именно по ней строят кривую ванны, считают среднюю наработку на отказ и проектируют системы резервирования. Ниже разберём формулы, три участка жизненного цикла изделия и типовые задачи - а калькулятор ниже позволит сразу почувствовать, как меняется кривая ванны при разных параметрах.

Определение и формула интенсивности отказов

Интенсивность отказов - это условная плотность вероятности отказа в момент $t$ при условии, что до этого момента изделие работало исправно. Формально:

$$\lambda(t) = \frac{f(t)}{P(t)},$$

где $f(t)$ - плотность распределения времени до отказа, $P(t)$ - вероятность безотказной работы. На практике в задачах с постоянной интенсивностью (нормальный участок кривой ванны) используют упрощение:

$$\lambda = \frac{n_{\text{от}}}{N \cdot t},$$

где $n_{\text{от}}$ - число отказавших элементов за время $t$, $N$ - начальное число элементов. Единица измерения - $1/\text{ч}$ (или $\text{ч}^{-1}$). Типичные значения для электронных компонентов - $10^{-6} \dots 10^{-4}\ \text{ч}^{-1}$, для механических узлов - $10^{-5} \dots 10^{-3}\ \text{ч}^{-1}$.

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

Связь между $\lambda(t)$ и вероятностью безотказной работы $P(t)$ - фундаментальная:

$$P(t) = \exp\!\left(-\int_0^t \lambda(\tau)\,d\tau\right).$$

Если $\lambda = \text{const}$ (постоянная интенсивность, нормальный участок), интеграл берётся явно:

$$P(t) = e^{-\lambda t}.$$

Эта экспоненциальная модель - самая распространённая в инженерных расчётах надёжности.

Кривая ванны: три участка жизненного цикла

Реальная зависимость $\lambda(t)$ обычно имеет форму «ванны» - убывает в начале, выходит на плато, затем нарастает. Три участка:

1. Участок приработки (ранних отказов). $\lambda(t)$ убывает - выявляются производственные дефекты, скрытые нарушения технологии. Изделие проходит «технологическую обкатку». Продолжительность - сотни-тысячи часов в зависимости от сложности.

2. Участок нормальной эксплуатации. $\lambda \approx \text{const}$ - отказы случайны, не связаны с накопленным повреждением. Именно здесь действует экспоненциальная модель. Длительность - основная часть срока службы.

3. Износовый участок. $\lambda(t)$ нарастает - накапливается усталость материала, окисление, коррозия, выработка ресурса. Начало износового участка определяет назначенный ресурс изделия.

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

Математически «ванну» описывают суперпозицией трёх слагаемых:

$$\lambda(t) = \underbrace{\frac{A}{(t + 1)^\alpha}}_{\text{приработка}} + \underbrace{\lambda_0}_{\text{норм. уч.}} + \underbrace{B\,e^{\beta t}}_{\text{износ}},$$

где $A$, $\alpha$, $\lambda_0$, $B$, $\beta$ - параметры, подбираемые по статистике отказов конкретного изделия. Именно эту модель использует калькулятор выше.

Средняя наработка на отказ MTBF

Для постоянной интенсивности отказов ($\lambda = \text{const}$) средняя наработка на отказ MTBF (Mean Time Between Failures) выражается просто:

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

Это математическое ожидание времени до отказа при экспоненциальном распределении. Важное следствие: вероятность безотказной работы за время $t = \text{MTBF}$ равна

$$P(\text{MTBF}) = e^{-\lambda \cdot \frac{1}{\lambda}} = e^{-1} \approx 0{,}368,$$

то есть лишь около 37%. Это часто удивляет студентов: MTBF - не «срок гарантированной работы», а лишь среднее время, и почти две трети изделий откажут раньше него.

Для нестационарной $\lambda(t)$ MTBF вычисляется интегрированием:

$$\text{MTBF} = \int_0^\infty P(t)\,dt = \int_0^\infty \exp\!\left(-\int_0^t \lambda(\tau)\,d\tau\right)dt.$$

В этом случае аналитического решения обычно нет, и прибегают к численному интегрированию.

Надёжность системы: последовательное и параллельное соединение

В системе из $n$ последовательно соединённых элементов отказ любого из них - отказ всей системы. Если элементы независимы с постоянными интенсивностями $\lambda_i$, интенсивность системы:

$$\lambda_{\Sigma} = \sum_{i=1}^n \lambda_i, \qquad \text{MTBF}_{\Sigma} = \frac{1}{\lambda_{\Sigma}}.$$

Вероятность безотказной работы системы:

$$P_{\Sigma}(t) = \prod_{i=1}^n P_i(t) = \exp\!\left(-\lambda_{\Sigma}\,t\right).$$

Пример: три элемента с $\lambda_1 = 10^{-4}$, $\lambda_2 = 2 \times 10^{-4}$, $\lambda_3 = 5 \times 10^{-5}\ \text{ч}^{-1}$. Суммарная интенсивность $\lambda_{\Sigma} = 3{,}5 \times 10^{-4}\ \text{ч}^{-1}$, $\text{MTBF} \approx 2857$ ч.

Для параллельного резервирования (хотя бы один из $n$ элементов работает) вероятность безотказной работы системы:

$$P_{\text{пар}}(t) = 1 - \prod_{i=1}^n (1 - P_i(t)).$$

Дублирование (два одинаковых элемента параллельно, $\lambda_i = \lambda$): $P_{\text{пар}}(t) = 2e^{-\lambda t} - e^{-2\lambda t}$, а эффективная интенсивность отказов нарастает от $0$ при $t = 0$ (оба работают) до $2\lambda$ при $t \to \infty$.

Пример задачи: расчёт надёжности системы

Рассмотрим типовую задачу. Система состоит из трёх последовательных блоков с интенсивностями отказов $\lambda_1 = 1 \times 10^{-4}\ \text{ч}^{-1}$, $\lambda_2 = 2 \times 10^{-4}\ \text{ч}^{-1}$, $\lambda_3 = 5 \times 10^{-5}\ \text{ч}^{-1}$. Требуется найти MTBF системы и вероятность безотказной работы за $t = 1000$ ч.

Суммарная интенсивность:

$$\lambda_{\Sigma} = \lambda_1 + \lambda_2 + \lambda_3 = (10 + 20 + 5) \times 10^{-5} = 3{,}5 \times 10^{-4}\ \text{ч}^{-1}.$$

Средняя наработка на отказ:

$$\text{MTBF} = \frac{1}{\lambda_{\Sigma}} = \frac{1}{3{,}5 \times 10^{-4}} \approx 2857\ \text{ч}.$$

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

$$P(1000) = e^{-\lambda_{\Sigma} \cdot 1000} = e^{-0{,}35} \approx 0{,}705 \quad (70{,}5\%).$$

Если ввести горячее резервирование для блока с наибольшей $\lambda_2$, то $\lambda_2^* = 2\lambda_2 e^{-\lambda_2 t_0}/(2 - e^{-\lambda_2 t_0})$. При малых $\lambda_2 t$ это даёт примерно $\lambda_2^* \approx \lambda_2^2 \cdot t$, что при $t = 1000$ ч и $\lambda_2 = 2 \times 10^{-4}$ составляет лишь $4 \times 10^{-5}$ - на порядок меньше исходного. Резервирование критичного блока существенно улучшает характеристики всей системы.

Расчёт по статистике отказов

В реальных задачах $\lambda$ определяют из испытаний: $N_0$ изделий поставлены на стенд, за время $\Delta t$ отказало $\Delta n$ из них. Оценка интенсивности:

$$\hat{\lambda} = \frac{\Delta n}{N_{\text{ср}} \cdot \Delta t},$$

где $N_{\text{ср}}$ - среднее число работоспособных изделий за интервал. Если испытания ступенчатые (разные партии снимаются в разные моменты), применяют метод максимального правдоподобия для экспоненциального распределения:

$$\hat{\lambda} = \frac{r}{\sum_{i=1}^N t_i},$$

где $r$ - число отказов за весь эксперимент, $\sum t_i$ - суммарная наработка всех изделий (включая не отказавшие до момента прекращения испытаний - «цензурированные» наблюдения).

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

• Подмена MTBF и «гарантированного срока». MTBF - среднее время до отказа, а не гарантия. За время $t = \text{MTBF}$ откажет около 63% изделий. Не путайте это с назначенным ресурсом.
• Несоблюдение единиц. Интенсивность отказов измеряется в $1/\text{ч}$ (или $10^{-6}/\text{ч} = \text{FIT}$). Если время задано в сутках или годах, $\lambda$ нужно пересчитать: $\lambda_{\text{ч}} = \lambda_{\text{год}} / 8760$.
• Применение постоянной $\lambda$ на износовом участке. Формула $\text{MTBF} = 1/\lambda$ верна только для участка нормальной эксплуатации с постоянной интенсивностью. На участке износа $\lambda(t)$ растёт, и экспоненциальная модель даёт завышенный прогноз надёжности.
• Ошибочное суммирование вероятностей. В последовательной системе вероятность отказа суммировать нельзя - суммируется $\lambda$ (для малых $\lambda t \ll 1$) или перемножаются $P_i(t)$.
• Игнорирование цензурированных данных. При ускоренных испытаниях часть изделий снимают досрочно - их наработку нельзя выбрасывать, иначе оценка $\lambda$ будет завышена.

FAQ

Чем отличается интенсивность отказов от вероятности отказа? Вероятность отказа $F(t) = 1 - P(t)$ - кумулятивная характеристика за интервал $[0, t]$. Интенсивность отказов $\lambda(t)$ - мгновенная характеристика: с какой скоростью происходят отказы в данный момент. Формально $\lambda(t) = f(t) / P(t)$, где $f(t) = dF/dt$. На нормальном участке $\lambda = \text{const}$, тогда как $F(t)$ продолжает нарастать.

Как связаны MTBF и надёжность за заданный срок? При постоянной $\lambda$ вероятность безотказной работы за время $t$ равна $P(t) = e^{-t/\text{MTBF}}$. Чтобы обеспечить надёжность $P = 0{,}9$ за 1000 ч, нужно $\text{MTBF} \geq -1000/\ln(0{,}9) \approx 9491$ ч, то есть $\lambda \leq 1{,}05 \times 10^{-4}\ \text{ч}^{-1}$.

Как рассчитать интенсивность отказов для системы с резервированием? Для горячего дублирования (оба элемента работают одновременно, отказ системы - отказ обоих) $P_{\text{сис}}(t) = 1 - (1 - e^{-\lambda t})^2$. Эффективная интенсивность системы нестационарна: $\lambda_{\text{сис}}(t) = 2\lambda e^{-\lambda t} / (2 - e^{-\lambda t})$, что начинается от 0 и нарастает до $2\lambda$ со временем. MTBF системы: $\text{MTBF}_{\text{сис}} = \int_0^\infty P_{\text{сис}}(t)\,dt = 3/(2\lambda)$ - в 1,5 раза выше одиночного элемента.

Коротко

Интенсивность отказов $\lambda(t)$ описывает мгновенную скорость выхода изделий из строя. На нормальном участке жизненного цикла $\lambda = \text{const}$, вероятность безотказной работы $P(t) = e^{-\lambda t}$, средняя наработка на отказ $\text{MTBF} = 1/\lambda$. Реальная зависимость $\lambda(t)$ имеет форму кривой ванны с тремя участками - приработки, нормальной эксплуатации и износа. В последовательных системах интенсивности отказов элементов суммируются; резервирование снижает суммарную интенсивность и повышает надёжность.