Средней наработкой до первого отказа называется математическое ожидание времени работы объекта до отказа.
Математическое ожидание средней наработки до отказа T 1 вычисляется через частоту отказов (плотность распределения времени безотказной
работы):
. (1.12)
Зная, что t > 0 и P(0) = 1, а P(∞) = 0, определяют T1:
. (1.13)
Средняя наработка до первого отказа, согласно статистическим данным об отказах, вычисляется по формуле
, (1.14)
где t i – время безотказной работы i-го образца; N0 – число испытываемых объектов.
Для определения средней наработки до первого отказа необходимо знать моменты выхода из строя всех испытываемых объектов. Поэтому для вычисления пользоваться данной формулой неудобно. Имея данные о количестве вышедших из строя элементов ni в каждом i-м интервале времени, среднюю наработку до первого отказа лучше определять по уравнению
, (1.15)
где t ср и m находятся по следующим формулам:
,
, (1.16)
где ti–1 – время начала i-го интервала; ti – время конца i-го интервала; tk – время, в течение которого вышли из строя все элементы; – интервал времени.
При расчетах надежности технических устройств часто применяются законы распределения: экспоненциальный, усеченный нормальный, Рэлея, гамма, Вейбулла – Гнеденко, логарифмически-нормальный. В табл. 1.1 приведены выражения для расчета количественных характеристик объектов, соответствующих перечисленным законам распределения времени их безотказной работы.
Таблица 1.1
Интенсивность отказов элементов
Закон распределения | Частота отказов (плотность распределения) | Вероятность безотказной работы | Интенсивность отказов | Средняя наработка до первого отказа |
Экспоненциальный | ![]() | ![]() | ![]() | ![]() |
Рэлея | ![]() | ![]() | ![]() | ![]() |
Гамма (при k целом) | ![]() |
![]() ![]() | ![]() | ![]() |
Вейбулла – Гнеденко | ![]() | ![]() | ![]() |
![]() |
Усеченный нормальный | ![]() | ![]() | ![]() | ![]() |
Логарифми- чески- нормальный | ![]() | ![]() | ![]() | ![]() |
Из приведенных соотношений видно, что все характеристики, кроме средней наработки до первого отказа, зависят от времени (являются функциями времени). На рис. 1.2 показаны зависимости количественных характеристик надежности объектов от времени.
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 1.2. Зависимости количественных характеристик надежности от времени:
а – экспоненциальный закон; б – усеченный нормальный закон; в – закон Рэлея;
г – гамма-распределение: д – закон Вейбулла – Гнеденко;
е – логарифмически-нормальный закон
Рассмотренные критерии надежности позволяют достаточно полно оценивать надежность невосстанавливаемых объектов, а также надежность восстанавливаемых объектов до первого отказа. Наличие нескольких критериев не означает, что всегда нужно оценивать надежность объектов по всем характеристикам.
Наиболее полной характеристикой надежности является частота отказов f(t) (плотность распределения), она содержит в себе все данные о случайном явлении – времени безотказной работы.
Средняя наработка до первого отказа является достаточно наглядной характеристикой надежности. Однако применение этого критерия для оценки надежности сложной системы ограничено в тех случаях, когда:
– время работы системы гораздо меньше среднего времени безотказной работы;
– закон распределения времени безотказной работы не однопараметрический и для достижения полной оценки требуются моменты высших порядков;
– система резервированная;
– интенсивность отказов не постоянная;
– время работы отдельных частей сложной системы разное.
Интенсивность отказа – наиболее удобная характеристика надежности простейших элементов, так как позволяет просто вычислять количественные характеристики надежности сложных систем.
Наиболее целесообразно оценивать надежность сложных систем по критерию вероятности безотказной работы, так как:
– она входит в качестве сомножителя в другие, более общие характеристики систем, например в эффективность и стоимость (цена и стоимость различаются);
– характеризует надежность с учетом изменения во времени;
– может быть получена сравнительно простыми расчетами в процессе проектирования систем и оценена в процессе испытаний.
1.4. Критерии надежности восстанавливаемых изделий
Пусть на испытании находится N изделий, и пусть отказавшие изделия немедленно заменяются исправными (новыми или отремонтированными).
Испытания считаются законченными, если число отказов достигает величины, достаточной для оценки надежности с определенной доверительной вероятностью. Если не учитывать времени, потребного на восстановление системы, то количественными характеристиками надежности могут быть параметр потока отказов μ(t) и наработка на отказ Т.
Параметром потока отказов называется отношение числа отказавших изделий в единицу времени к числу испытываемых изделий при условии, что все вышедшие из строя изделия заменяются исправными (новыми или отремонтированными).
Согласно определению,
, (1.17)
где ∆ t – малый отрезок наработки; r(t) – число отказов, наступивших от начального момента времени до достижения наработки t. Разность r(t +∆ t) – r(t) представляет собой число отказов на отрезке ∆ t.
Статистическую оценку параметра потока отказов дают по формуле
. (1.18)
Для стационарных потоков можно применять формулу
, (1.19)
где – оценка средней наработки на отказ;
, здесь t – суммарная наработка, r(t) – число отказов, наступивших в течение этой наработки,
– математическое ожидание этого числа.
Параметр потока отказов определяется также по формуле
, (1.20)
где – число отказавших образцов в интервале времени от
до
; N – число испытываемых образцов;
– интервал времени.
Формула (1.18) является статистическим определением параметра потока отказов.
Параметр потока отказов и частота отказов для ординарных потоков с ограниченным последействием связаны интегральным уравнением Вольтерра второго рода
. (1.21)
По известной f(t) можно найти все количественные характеристики надежности невосстанавливаемых изделий. Поэтому (1.21) является основным уравнением, связывающим количественные характеристики надежности невосстанавливаемых и восстанавливаемых изделий при мгновенном восстановлении.
Уравнение (1.21) можно записать в операторной форме:
,
. (1.22)
Соотношения (1.22) позволяют найти одну характеристику через другую, если существуют преобразования Лапласа функций a(s) и m(s) и обратные преобразования выражений (1.22).
Параметр потока отказов обладает следующими важными свойствами:
1) для любого момента времени независимо от закона распределения времени безотказной работы параметр потока отказов больше, чем частота отказов, т. е. μ(t) >> f(t);
2) независимо от вида функции f(t) параметр потока отказов μ( t) при
t → ∞ стремится к 1/Тср. Это важное свойство параметра потока отказов означает, что при длительной эксплуатации ремонтируемого изделия поток его отказов независимо от закона распределения времени безотказной работы становится стационарным. Однако это не означает, что интенсивность отказов есть величина постоянная;
3) если l(t) – возрастающая функция времени, то
l ( t) > μ(t) > f(t),
4) если l(t) – убывающая функция, то
f(t) > l(t) > μ(t);
5) при l(t) = const параметр потока отказов системы не равен сумме параметров потоков отказов элементов, т. е.
. (1.23)
Согласно этому свойству параметра потока отказов можно утверждать, что при вычислении количественных характеристик надежности сложной системы нельзя суммировать имеющиеся в настоящее время значения интенсивностей отказов элементов, полученные по статистическим данным об отказах изделий в условиях эксплуатации, так как указанные величины являются фактически параметрами потока отказов;
6) при l(t) = l = const параметр потока отказов равен интенсивности отказов μ(t) = l(t) = l .
Сравнение свойств интенсивности и параметра потока отказов свидетельствует, что эти характеристики различны.
В настоящее время широко используются статистические данные об отказах, полученные в условиях эксплуатации аппаратуры. При этом они часто обрабатываются таким образом, что приводимые характеристики надежности являются не интенсивностью отказов, а параметром потока отказов μ(t). Это приводит к ошибкам при расчетах надежности. В ряде случаев они могут быть значительными.
Для получения интенсивности отказов элементов из статистических данных об отказах ремонтируемых систем необходимо воспользоваться формулой (1.6), для чего следует знать предысторию каждого элемента принципиальной схемы. Это может существенно усложнить методику сбора статистических данных об отказах. Поэтому целесообразно определять l(t) по параметру потока отказов μ(t). Методика расчета сводится к следующим вычислительным операциям:
– по статистическим данным об отказах элементов ремонтируемых изделий и по формуле (1.13) вычисляется параметр потока отказов и строится гистограмма μi(t);
– гистограмма заменяется кривой, которая аппроксимируется уравнением;
– находится преобразование Лапласа μi(s) функции μi(t);
– по известной μi(s) на основании (1.15) записывается преобразование Лапласа fi(s) частоты отказов;
– по известной fi(s) находится обратное преобразование частоты отказов fi(t);
– находится аналитическое выражение для интенсивности отказов по формуле
; (1.24)
– строится график l i(t).
Если имеется участок, где l i(t) = l = const, то постоянное значение интенсивности отказов принимается для оценки вероятности безотказной работы. При этом считается справедливым экспоненциальный закон надежности.
Приведенная методика не может быть применена, если не удается найти по f(s) обратное преобразование частоты отказов f(t). В этом случае приходится применять приближенные методы решения интегрального уравнения (1.21) или машинные методы расчета.