bilim.techО проекте

Надежность и отказы

Почему новая машина ломается чаще обкатанной, а старая снова чаще? Запусти парк машин, смотри, как гаснут отказавшие, и строй кривую ванны для подшипников, шин и гидронасосов.

01 Что это такое

Надежность это вероятность, что узел проработает время t и не откажет.

Когда откажет конкретный подшипник, сказать нельзя. Но если поставить в работу сотню одинаковых подшипников, доля исправных со временем падает по вполне определенной кривой R(t). Скорость, с которой отказывают еще работающие узлы, называют интенсивностью отказов. У большинства машин она меняется по «кривой ванны»: в начале высокая из-за брака сборки, потом долго низкая и ровная, а к концу ресурса снова растет из-за износа. От того, на каком участке кривой работает узел, зависит, поможет ли плановая замена.

приработкаβ < 1нормальная работаβ = 1, λ постояннаизносβ > 1брак сборки,ошибки монтажаудары, перегрузки, случайностиусталость, истираниенаработка, чинтенсивность отказов h(t)
Кривая ванны. Плановая замена узла помогает только на участке износа. На участке случайных отказов новый узел ломается так же часто, как старый.
R(t)
Надежность. Доля узлов, которые проработали время t без отказа. В начале 1, со временем падает к 0.
h(t), λ
Интенсивность отказов. Какая доля еще работающих узлов откажет за следующий час. Если она постоянна, ее обозначают λ.
MTBF
Средняя наработка на отказ. Сколько часов в среднем работает узел до отказа. При постоянной λ равна 1/λ.
02 Опыт
Пусти время. Смотри, как гаснут машины и как это повторяет кривая R(t). РАСЧЕТ В РЕАЛЬНОМ ВРЕМЕНИ
ВЕРОЯТНОСТЬ И СТАТИСТИКА

Сколько машин останется в строю?

Рис. 7.6
СЛУЧАЙ
Машин в строю
Надежность R(t)
Средняя наработка, MTBF

ⓘ Наработка каждого узла до отказа это случайное число по закону Вейбулла: t = η · (−ln u)1/β, где u равномерно от 0 до 1. Отказавшие машины в опыте не ремонтируют. Средняя наработка считается через гамма-функцию: MTBF = η · Γ(1 + 1/β).

03 Формулы по буквам

Экспонента, Вейбулл и схемы соединения

Первая формула описывает случайные отказы с постоянной интенсивностью. Вторая добавляет приработку и износ. Третья собирает надежность машины из надежностей узлов. Под каждой формулой числа из опыта.

Сколько узлов доживет до t, если отказы случайные?

R(t) = e−λt,  MTBF = 1λ

Каждый час отказывает одна и та же доля λ еще работающих узлов. Из такого правила получается экспонента. Узел не стареет: проработавший 10 000 часов откажет в следующий час так же вероятно, как новый. Средняя наработка обратна λ.

λ
интенсивность отказов, 1/ч. Например, 0,05 на 1000 ч
t
наработка, ч
MTBF
средняя наработка на отказ, ч. При t = MTBF в строю e⁻¹ ≈ 37 %

А если узел прирабатывается или изнашивается?

R(t) = e−(t/η)β,  h(t) = βη · (t/η)β − 1

Закон Вейбулла. Параметр β задает форму, η задает масштаб времени. При β = 1 это снова экспонента с λ = 1/η. При β больше 1 интенсивность h(t) растет со временем: так ведут себя подшипники и шины. При β меньше 1 она падает: так проявляется брак, который отсеивается в первые часы.

β
параметр формы: 0,5…1 приработка, 1 случайные отказы, 2…4 износ
η
характерный ресурс, ч: к нему откажет 63 % узлов
h(t)
интенсивность отказов в момент t

Какая надежность у машины из нескольких узлов?

Rцеп = R₁ · R₂ · … · Rk
Rрез = 1 − (1 − R)²

Цепочка работает, только если работают все узлы, поэтому надежности перемножаются, и цепочка всегда хуже самого слабого узла. Резерв отказывает, только если откажут оба узла. Вероятность этого (1 − R)², все остальное машина работает.

R₁…Rk
надежности узлов в цепочке в один и тот же момент t
1 − R
вероятность, что узел уже отказал
Rрез
два одинаковых узла, хватает одного

Откуда берется экспонента

Пусть каждый час отказывает четверть еще работающих узлов, λ = 0,25 1/ч. Из 100 узлов через час останется 75, через два часа 75 · 0,75 ≈ 56, через три 42. Каждый час число умножается на одно и то же 0,75. Такое убывание и есть экспонента: доля в строю R = (1 − λ)t, а при малых шагах времени R = e−λt.

В законе Вейбулла доля отказов за час не постоянна, а меняется как h(t). Тогда в показатель экспоненты попадает не λt, а накопленная интенсивность (t/η)β. Средняя наработка это площадь под кривой R(t): для экспоненты она ровно 1/λ.

dR/dt = −h(t) · R
убыль работающих пропорциональна их числу
R = e−∫h dt
решение этого уравнения
MTBF = ∫₀∞ R(t) dt
средняя наработка: площадь под R(t)
1000751562423324245186137108часов работыR = 0,75t = e−0,29tкаждый час ×0,75
Числа над столбиками: сколько узлов из 100 еще работают. Каждый час их становится в одно и то же число раз меньше.
04 Где и зачем это применяют

Везде, где планируют ремонты и запчасти

Допущения модели+

Узлы в опыте отказывают независимо, все одинаковые, и отказавшие не ремонтируют. В жизни машины чинят, и тогда считают уже коэффициент готовности, где важно и время ремонта. Нагрузка в модели постоянная, а на деле смена режима, температура и пыль сдвигают η. Резервный узел работает все время вместе с основным, поэтому изнашивается так же. Холодный резерв, который включается только при отказе, надежнее. Кривая ванны на рисунке нарисована схематично: закон Вейбулла с одним β описывает только один ее участок.

Еще в направлении «Математика»

Все темы направления «Математика» →