Как найти вероятность отказа агрегата

1. Показатели безотказности

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

,

где N_р
– число работоспособных объектов на
момент t;

N – общее число наблюдаемых
объектов;

n(t) – число объектов,
отказавших на момент t от начала испытаний
или эксплуатации.

Вероятность
безотказной работы умень­­шается
с увеличением времени работы или
наработки объекта. Зависимость вероятности
безотказной работы от времени
характеризуется кривой убыли ресурса
объ­екта, пример которой приведен на
рисунке 9.

В начальный момент времени
для работоспособного объекта вероятность
его безотказной работы равна единице
(100 %). По мере работы объекта эта вероятность
снижается и стремится к нулю.

Вероятность отказа
характеризуется вероятностью возникновения
отказа на момент времени t
:

,

где
n(t) – число объектов, отказавших на
момент t от начала испытаний или
эксплуатации;

N – общее число наблюдаемых объектов.

Вероятность возникновения
отказа объекта возрастает с увеличе­нием
срока эксплуатации или наработки.

Пример зависимости вероятности
возникновения отказа от времени показан
на рисунке 10. Для работоспособного
объекта в начальный момент времени
вероятность отказа близка к нулю. Для
того, чтобы отказ проявился, объекту
необходимо начать работать, при этом
вероятность отказа увеличивается с
увеличением времени и стремится к
единице.

Вероятность отказа может
быть также охарактеризована плотностью
вероятности отказа

или,

где
– число отказов за промежуток времени
Δt;

N – общее число наблюдаемых
объектов.

Пример
1. После 500 часов
наработки из 56 агрегатов, поставленных
на эксплуатацию, в работоспособном
состоянии оказалось 43 агрегата. Определить
вероятность безотказной работы агрегата
в течение 500 час.

Решение:

Используем формулу для определения
вероятности безотказной работы объекта

.

Вероятность безотказной
работы агрегата в течение 500 часов
составляет 76,8 %.

Пример
2. Для предыдущего
примера определить вероятность отказа
агрегат за 500 часов работы.

Решение:

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

или

.

Таким образом, вероятность
отказа агрегата за 500 часов составляет
23,2 %.

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

1. Вероятность появления
одного из двух несовместных событий
равна сумме вероятности этих событий:

,

где A,
B
– несовместные события.

2. Вероятность совместного
появления нескольких независимых
событий равна произведению вероятностей
этих событий:

.

Первая теорема используется
для нахождения вероятности отказа при
возможности у объекта нескольких видов
несовместных отказов. С использованием
второй теоремы определяют вероятность
безотказной работы объекта, состоящего
из многих элементов, вероятность
безотказной работы которых известна.

Пример
3. Система состоит
их 4-х агрегатов. Надежность каждого
агрегата в течение времени t
характеризуется вероятностью безотказной
работы 90 %. Найти вероятность безотказной
работы системы в течение времени t
при условии независимости отказов
агрегатов.

Решение:

Используем теорему вероятности
совместного появления работоспособного
состояния всех агрегатов:

.

Следовательно, вероятность
безотказной работы системы в течение
времени t
равна 65,6 %.

Пример
4. В составе
агрегата имеются 5 узлов. Вероятность
отказа каждого узла в течение времени
t
составляет 5 %. Отказы узлов несовместны.
Определить вероятность отказа агрегата.

Решение:

Используем теорему для
вероятности хотя бы одного из нескольких
несовместных событий:

.

Таким образом, вероятность
отказа агрегата в течение времени t
составляет 25 %.

2. Интенсивность
отказов –
характеризует скорость возникновения
отказов объекта в различные моменты
времени его работы:

,

где n(t)
– число отказов за промежуток времени
t;

N_р
– число работоспособных объектов на
момент t.

Интенсивность отказов может быть найдена
теоретически

,

где f(t) – функция плотности
вероятности наработки до отказа;

P(t) – вероятность безотказной
работы,

.

Плотность распределения
f(t) наработки до отказа может быть также
определена по вероятности отказа

или.

Вероятность безотказной работы связана
с интенсивностью отказов одним из
основных уравнений теории надежности:

.

В описанных способах оценки
безотказности до первого отказа отказы
не различаются по тяжести их последствий.
В большинстве случаев при разработке
объекта необходимо установить критерий
отказа изделия по экономическим
соображением, исчерпанию ресурса или
другим характеристикам.

Критерием
отказа называют
признак или совокупность признаков
неработоспособного состояния объекта,
установленных в нормативно-технической
или конструкторской документации.

3. Средняя
наработка на отказ
– это отношение наработки восстанавливаемого
объекта к математическому ожиданию
числа его отказов в течение этой
наработки:

,

где N – общее число объектов,
поставленных на испытания или в
эксплуатацию;

t _i
– наработка i-того объекта;

m _i
– число отказов i-того объекта за весь
наблюдаемый период.

Средняя наработка на отказ
используется для характеристики
восстанавливаемых объектов.

4.
Средняя наработка до отказа –
математическое ожидание наработки
объекта до первого отказа

или,

где N_pi
– число работоспособных объектов на
интервале наработки t_i–t_i+1
;

N – общее число наблюдаемых объектов;

t
= t_i+1–t_i
– интервал времени;

k – число рассматриваемых
интервалов наработки.

Среднюю наработку до отказа можно также
определить иначе

,

где t_i
– наработка до отказа i-того объекта;

n – число объектов.

Показатель Т_ср
используется для характеристики
надежности невосстанавливаемых объектов.

5.
Средняя
наработка между отказами
– математическое ожидание наработки
объекта от окончания восстановления
его работоспособного состояния после
отказа до возникновения следующего
отказа.

Вычисляется как отношение суммарной
наработки объекта между отказами за
рассматриваемый период к числу отказов
за тот же период:

.

Показатели безотказности определяют
на разных стадиях работы объекта с целью
его совершенствования и с целью контроля
нормируемых значений при эксплуатации.

Соседние файлы в папке НАДЕЖНОСТЬ_3_ГРАДИРНИ

• #

  21.03.20168.84 Mб51НАДЕЖНОСТЬ_Ушаков_ГРАДИРНИ.pdf

• #
• #
• #

  21.03.201616.52 Mб54НАДЕЖНОСТЬ_Эдельман_ГРАДИРНИ.pdf

• #

Практика

1. Расчет показателей безотказности
1.1 Вероятность безотказной работы
1.2 Вероятность отказа
1.3 Частота отказа
1.4 Интенсивность отказа
1.5 Средняя наработка до отказа
1.6 Среднее значение параметра потока отказов
1.7 Пример расчета показателей безотказности
2. Примеры расчета показателей надежности для различных законов распределения случайных величин
2.1 Экспоненциальный закон распределения
2.2  Закон распределения Вейбулла-Гнеденко
2.3  Закон распределения Рэлея
3. Примеры расчета показателей надежности сложных систем
3.1 Основное соединение элементов
3.2 Резервное соединение

1.1 Вероятность безотказной работы

Вероятностью безотказной работы называется вероятность того, что при определенных условиях эксплуатации, в пределах заданной наработки не произойдет ни одного отказа.
Вероятность безотказной работы обозначается как P(l), которая определяется по формуле (1.1):

где N₀ – число элементов в начале испытания; r(l) – число отказов элементов к моменту наработки.Следует отметить, что чем больше величина N₀, тем с большей точностью можно рассчитать вероятность P(l).
В начале эксплуатации исправного локомотива P(0) = 1, так как при пробеге l = 0 вероятность того, что ни один элемент не откажет, принимает максимальное значение – 1. С ростом пробега l вероятность P(l) будет уменьшаться. В процессе приближения срока эксплуатации к бесконечно большой величине вероятность безотказной работы будет стремиться к нулю P(l→∞) = 0. Таким образом в процессе наработки величина вероятности безотказной работы изменяется в пределах от 1 до 0. Характер изменения вероятности безотказной работы в функции пробега показан на рис. 1.1.

Рис.2.1. График изменения вероятности безотказной работы P(l)в зависимости от наработки

Основными достоинствами использования данного показателя при расчетах является два фактора: во-первых, вероятность безотказной работы охватывает все факторы, влияющие на надежность элементов, позволяя достаточно просто судить о его надежности, т.к. чем больше величина P(l), тем выше надежность; во-вторых, вероятность безотказной работы может быть использована в расчетах надежности сложных систем, состоящих из более чем одного элемента.

1.2 Вероятность отказа

Вероятностью отказа называют вероятность того, что при определенных условиях эксплуатации, в пределах заданной наработки произойдет хотя бы один отказ.
Вероятность отказа обозначается как Q(l), которая определяется по формуле (1.2):

В начале эксплуатации исправного локомотива Q(0) = 0, так как при пробеге l = 0 вероятность того, что хотя бы один элемент откажет, принимает минимальное значение – 0. С ростом пробега l вероятность отказа Q(l) будет увеличиваться. В процессе приближения срока эксплуатации к бесконечно большой величине вероятность отказа будет стремиться к единице Q(l→∞) = 1. Таким образом в процессе наработки величина вероятности отказа изменяется в пределах от 0 до 1. Характер изменения вероятности отказа в функции пробега показан на рис. 1.2.Вероятность безотказной работы и вероятность отказа являются событиями противоположными и несовместимыми.

Рис.2.2. График изменения вероятности отказа Q(l) в зависимости от наработки

1.3 Частота отказов

Частота отказов – это отношение числа элементов в единицу времени или пробега отнесенного к первоначальному числу испытуемых элементов. Другими словами частота отказов является показателем, характеризующим скорость изменения вероятности отказов и вероятности безотказной работы по мере роста длительности работы.
Частота отказов обозначается как  и определяется по формуле (1.3):

где  –  количество отказавших элементов за промежуток пробега . 
Данный показатель позволяет судить по его величине о числе элементов, которые откажут на каком-то промежутке времени или пробега, также по его величине можно рассчитать количество требуемых запасных частей.
Характер изменения частоты отказов в функции пробега показан на рис. 1.3.

Рис. 1.3. График изменения частоты отказов в зависимости от наработки

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

Интенсивность отказов представляет собой условную плотность возникновения отказа объекта, определяемую для рассматриваемого момента времени или наработки при условии, что до этого момента отказ не возник. Иначе интенсивность отказов – это отношение числа отказавших элементов в единицу времени или пробега к числу исправно работающих элементов в данный отрезок времени.
Интенсивность отказов обозначается как  и определяется по формуле (1.4):

где 

Как правило, интенсивность отказов является неубывающей функцией времени. Интенсивность отказов обычно применяется для оценки склонности к отказам в различные моменты работы объектов.
На рис. 1.4. представлен теоретический характер изменения интенсивности отказов в функции пробега.

Рис. 1.4. График изменения интенсивности отказов в зависимости от наработки

На графике изменения интенсивности отказов, изображенном на рис. 1.4. можно выделить три основных этапа отражающих процесс экс-плуатации элемента или объекта в целом.
Первый этап, который также называется этапом приработки, характеризуется увеличением интенсивности отказов в начальный период эксплуатации. Причиной роста интенсивности отказов на данном этапе являются скрытые дефекты производственного характера.
Второй этап, или период нормальной работы, характеризуется стремлением интенсивности отказов к постоянному значению. В течение этого периода могут возникать случайные отказы, в связи с появлением внезапной концентрации нагрузки, превышающей предел прочности элемента.
Третий этап, так называемый период форсированного старения. Характеризуется возникновением износовых отказов. Дальнейшая эксплуатация элемента без его замены становится экономически не рациональной.

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

Средняя наработка до отказа – это средний пробег безотказной работы элемента до отказа. 
Средняя наработка до отказа обозначается как L₁ и определяется по формуле (1.5):

 

где l_i – наработка до отказа элемента; r_i – число отказов.
Средняя наработка до отказа может быть использована для предварительного определения сроков ремонта или замены элемента.

1.6 Среднее значение параметра потока отказов

Среднее значение параметра потока отказов характеризует среднюю плотность вероятности возникновения отказа объекта, определяемая для рассматриваемого момента времени.
Среднее значение параметра потока отказов обозначается как W_ср и определяется по формуле (1.6):

1.7 Пример расчета показателей безотказности

Исходные данные.
В течение пробега от 0 до 600 тыс. км., в локомотивном депо произведен сбор информации по отказам ТЭД. При этом количество исправных ТЭД в начале периода эксплуатации составляло N0 = 180 шт. Суммарное количество отказавших ТЭД за анализируемый период составило ∑r(600000) = 60. Интервал пробега   принять равным 100 тыс. км. При этом количество отказавших ТЭД по каждому участку составило: 2, 12, 16, 10, 14, 6.

Требуется.
Необходимо рассчитать показатели безотказности и построить их зависимости изменения во времени.

Сначала необходимо заполнить таблицу исходных данных так, как это показано в табл. 1.1.

Таблица 1.1.

Исходные данные к расчету

, тыс. км	0 – 100	100 – 200	200 – 300	300 – 400	400 – 500	500 – 600
	2	12	16	10	14	6
	2	14	30	40	54	60

 Первоначально по уравнению (1.1) определим для каждого участка пробега величину вероятности безотказной работы. Так, для участка от 0 до 100 и от 100 до 200 тыс. км. пробега вероятность безотказной работы составит:

Далее, используя зависимость (1.2) произведем расчет вероятности отказа ТЭД.

Произведем расчет частоты отказов по уравнению (1.3).

Далее по уравнению (1.4) произведем расчет интенсивности отказов ТЭД в зависимости от наработки.
Первоначально рассчитаем среднее количество работоспособных ТЭД на участке от 0 до 100 тыс. км. пробега:

Тогда интенсивность отказов на участке 0-100 тыс.км. будет равна:

Аналогичным образом определим величину интенсивности отказов для интервала 100-200 тыс. км.

По уравнениям (1.5 и 1.6) определим среднюю наработку до отказа и среднее значение параметра потока отказов.

Систематизируем полученные результаты расчета и представим их в виде таблицы (табл. 1.2.).

Таблица 1.2.

Результаты расчета показателей безотказности

, тыс.км.	0 – 100	100 – 200	200 – 300	300 – 400	400 – 500	500 – 600
	2	12	16	10	14	6
	2	14	30	40	54	60
P(l)	0,989	0,922	0,833	0,778	0,7	0,667
Q(l)	0,011	0,078	0,167	0,222	0,3	0,333
10-7, 1/км	1,111	6,667	8,889	5,556	7,778	3,333
10-7, 1/км	1,117	6,977	10,127	6,897	10,526	4,878

 Приведем характер изменения вероятности безотказной работы ТЭД в зависимости от пробега (рис. 1.5.). Необходимо отметить, что первой точкой на графике, т.е. при пробеге равном 0, величина вероятности безотказной работы примет максимальное значение – 1. 

Рис. 1.5. График изменения вероятности безотказной работы в зависимости от наработки

Приведем характер изменения вероятности отказа ТЭД в зависимости от пробега (рис. 1.6.). Необходимо отметить, что первой точкой на графике, т.е. при пробеге равном 0, величина вероятности отказа примет минимальное значение – 0. 

Рис. 1.6. График изменения вероятности отказа в зависимости от наработки

Приведем характер изменения частоты отказов ТЭД в зависимости от пробега (рис. 1.7.).

Рис. 1.7. График изменения частоты отказов в зависимости от наработки
 

На рис. 1.8. представлена зависимость изменения интенсивности отказов от наработки.

Рис. 1.8. График изменения интенсивности отказов в зависимости от наработки

2.1 Экспоненциальный закон распределения случайных величин

Во-первых, определим величину средней наработки топливных насосов до отказа по уравнению:

Затем рассчитываем величину интенсивности отказов:

Величина вероятности безотказной работы топливных насосов при наработке 500 ч составит:

Вероятность отказа в промежутке между 800 и 900 ч. работы насосов составит:

2.2 Закон распределения Вэйбулла-Гнеденко

Закон распределения Вейбулла-Гнеденко получил широкое распространение и используется применительно к системам, состоящим из рядов элементов, соединенных последовательно с точки зрения обеспечения безотказности системы. Например, системы, обслуживающие дизель-генераторную установку: смазки, охлаждения, питания топливом, воздухом и т.д.

Найдем вероятность восстановления работоспособности локомотива после простоя его в депо в течение суток по уравнению:

Для определения времени восстановления работоспособности локомотива с заданной величиной доверительной вероятности также используем выражение:

2.3 Закон распределения Рэлея

Требуется.
Для величины наработки 120 тыс.км. необходимо определить вероятность безотказной работы, интенсивность отказов и среднюю наработку до первого отказа катушки электромагнитного контактора. 

3.1 Основное соединение элементов

Система, состоящая из нескольких независимых элементов, связанных функционально таким образом, что отказ любого из них вызывает отказ системы, отображается расчетной структурной схемой безотказной работы с последовательно соединенными событиями безотказной работы элементов.

Вычислим интенсивность отказа и среднюю наработку до отказа по следующим уравнениям:

 

Значения вероятности безотказной работы и частоты отказов получим, используя уравнения приведенные к виду:

Результаты расчета P(l) и a(l) на интервале от 0 до 1000 часов работы представим в виде табл. 3.1.

Таблица 3.1.

Результаты расчета вероятности безотказной работы и частоты отказов системы на интервале времени от 0 до 1000 ч.

l, час	P(l)	a(l), час-1
0	1	0,00026
100	0,974355	0,000253
200	0,949329	0,000247
300	0,924964	0,00024
400	0,901225	0,000234
500	0,878095	0,000228
600	0,855559	0,000222
700	0,833601	0,000217
800	0,812207	0,000211
900	0,791362	0,000206
1000	0,771052	0,0002

Графическая иллюстрация P(l) и a(l) на участке до средней наработки до отказа представлена на рис. 3.1, 3.2. 

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

Рис. 3.2. Частота отказов системы.

3.2 Резервное соединение элементов

Исходные данные.
На рис. 3.3 и 3.4 показаны две структурные схемы соединения элементов: общего (рис. 3.3) и поэлементного резервирования (рис. 3.4). Вероятности безотказной работы элементов соответственно равны P1(l) = P ’1(l) = 0,95; P2(l) = P’2(l) = 0,9; P3(l) = P ’3(l) = 0,85. 

Требуется.
Необходимо рассчитать надежность двух систем. 

Рис. 3.3. Схема системы с общим резервированием.

Рис. 3.4. Схема системы с поэлементным резервированием.

Вероятность безотказной работы блока из трех элементов без резервирования рассчитаем по выражению:

Вероятность безотказной работы той же системы при общем резервировании (рис. 3.3) составит:

Вероятности безотказной работы каждого из трех блоков при поэлементном резервировании (рис. 3.4) будут равны:

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

Таким образом, поэлементное резервирование дает более существенное увеличение надежности (вероятность безотказной работы возросла с 0,925 до 0,965, т.е. на 4%).

Исходные данные.
На рис. 3.5 представлена система с комбинированным соединением элементов. При этом вероятности безотказной работы элементов имеют следующие значения: P1=0,8; Р2=0,9; Р3=0,95; Р4=0,97.

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

Рис.3.5. Схема системы при комбинированном функционировании элементов.

Для расчета в исходной системе необходимо выделить основные блоки. В представленной системе их три (рис. 3.6). Далее рассчитаем надежность каждого блока в отдельности, а затем найдем надежность всей системы.

Рис. 3.6. Сблокированная схема.

Надежность системы без резервирования составит:

Таким образом, система без резервирования является на 28% менее надежной, чем система с резервированием.

Занятие № 5

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

Цель:

ü рассмотреть основные
показатели безотказной работы (вероятность безотказной работы и вероятность
отказа).

Основные понятия:

·       
Показатели надежности

·       
Вероятность безотказной
работы

·       
Вероятность отказа

План занятия:

1.     Организационный момент: приветствие, проверяется готовность к
занятию, отмечаются в журнале отсутствующие.

2.     Проверка домашнего задания: фронтальный опрос.

3.     Актуализация знаний: сообщение темы и цели занятия.

4.     Изучение нового материала:

1.    
Показатели надежности.

2.    
Вероятность безотказной
работы.

3.    
Вероятность отказа.

4.    
Оценка вероятности
безотказной работы.

5.     Закрепление
изученного материала:

ü Что называют
вероятностью безотказной работы?

ü Что называют вероятностью
отказа?

ü Дайте определение
понятия показатели надежности.

6.     Домашнее
задание:
Яхъяев Н.Я. Основы теории надежности, стр. 39-41;

Задача:

На
испытание поставлено 1000 однотипных резисторов. За первые 10000 часов отказало
5, за последующие 5000 отказало еще 5. Определить вероятность безотказной
работы и вероятность отказа за 10000 часов, за 15000 часов и в промежутке между
10000 и 15000 часов.

7.     Подведение
итогов занятия.

– Какое
состояние называется работоспособным? (Работоспособность – это состояние
изделия, при котором оно способно выполнять заданную функцию с параметрами,
установленными требованиями технической документации, в течение расчётного
срока службы).

Отказ –
это нарушение работоспособности. Свойство элемента или системы непрерывно
сохранять работоспособность при определённых условиях эксплуатации (до первого
отказа) называется безотказностью.

Безотказность
– свойство объекта сохранять работоспособное состояние в течение некоторого
времени или наработки.

И сегодня
на занятии мы рассмотрим показатели безотказности. Тема нашего занятия: «ВЕРОЯТНОСТЬ
БЕЗОТКАЗНОЙ РАБОТЫ И ВЕРОЯТНОСТЬ ОТКАЗА».

1.     ПОКАЗАТЕЛИ
НАДЕЖНОСТИ.

Показатели
надежности – количественная характеристика одного или нескольких
свойств, составляющих надежность объекта.

Для
оценки, расчетов и исследования надежности технических устройств в процессе их
проектирования и эксплуатации используются количественные характеристики (критерии
надежности). Для показателей надежности используются две формы  представления:

v Статистическая – при
эксперементальном исследовании надежности технических систем

v Вероятностная – при
априорных аналитических расчетах надежности.

В
соответствии с ГОСТ 27.002 – 89 показатели надежности подразделяются на:

Классификация
показателей надежности

Признак	Показатель
Число характерезуемых свойств надежности	Единичный показатель надежности – показатель, характеризующий одно из свойств, составляющих надежность объекта. (Например, вероятность отказа, средний срок службы и т.п.).
Комплексный показатель надежности – показатель, характеризующий одновременно несколько свойств, составляющих надежность объекта. (Например, коэффициент готовности, удельная суммарная трудоемкость ремонтов и т.п.).
Свойство надежности	Безотказность
Долговечность
Сохраняемость
Ремонтопригодность
Метод получения	Расчетный показатель надежности – показатель, значения которого определяют расчетным методом.
Экспериментальный показатель надежности – показатель, точечную или интервальную оценку которого определяют по данным испытания.
Эксплуатационный показатель надежности – показатель, точечную или интервальную оценку которого определяют по данным эксплуатации.
Экстраполированный показатель надежности – показатель, точечную или интервальную оценку которого определяют на основании результатов расчетов, испытаний и (или) эксплуатационных данных путем экстраполирования на другую продолжительность эксплуатации и другие условия эксплуатации.
Область использования	Нормативный показатель, регламентированный в НТД
Оценочный показатель, используемый для различных сравнительных оценок при научно-исследовательских и проектно-технологических разработках
Область распространения	Групповой показатель надежности – служит для оценки надежности совокупности изделий данного типа.
Индивидуальный показатель надежности – предназначен для оценки надежности каждого изделия данного типа.

2.     ВЕРОЯТНОСТЬ
БЕЗОТКАЗНОЙ РАБОТЫ.

Пусть
испытывается некоторое число изделий N₀. По разным причинам они
будут выходить из строя, причем моменты отказов, т.е. время наработки каждого
изделия до отказов, является случайной величиной.

Вероятность
безотказной работы изделия есть вероятность того, что за
определенный рассматриваемый период времени работы (t) в заданных условиях
эксплуатации оно не откажет, т.е. вероятность того, что время наработки до
отказа (t_отк) будет больше времени работы.

Р(t) = Вер (t_отк>
t)

Если к
моменту t из поставленных на испытания N0 изделий останутся исправными N(t), то
статистическая вероятность безотказной работы изделия за время t, равно: , где  N(t) – число работоспособных
изделий на момент t; N₀ – общее число наблюдаемых изделий; n(t)  – 
число  изделий,  отказавших  на момент  t  от  начала  испытаний.

При t = 0
все изделия исправны N(0) = N₀ и P(0) = 1. Отказы изделия с течением
времени t приводят к монотонному убыванию функции Р(t). Практически для каждого
типа изделия существует наработка t*, больше которой ни одно изделие данного
типа проработать не может.

N(t) = 0,
при t ≥ t*соответственно Р(t) = 0, при t ≥ t*

Таким
образом, 0 ≤ P(t) ≤ 1.

Вероятность
безотказной работы уменьшается с увеличением времени работы или наработки
объекта. Зависимость вероятности безотказной работы от времени характеризуется
кривой убыли ресурса изделия, пример которой приведен на рисунке 1. 

Рис. 1

В
начальный момент времени для работоспособного изделия вероятность его
безотказной работы равна единице (100%). По мере работы объекта эта вероятность
снижается и стремится к нулю.

Например: После
500 часов наработки из 56 агрегатов, поставленных на эксплуатацию, в работоспособном
состоянии оказалось 43 агрегата. Определить вероятность безотказной работы
агрегата в течение 500 час.

Решение:

Используем
формулу для определения вероятности безотказной работы объекта

Вероятность
безотказной работы агрегата в течение 500 часов составляет 76,8 %.

3.     ВЕРОЯТНОСТЬ
ОТКАЗА.

Противоположным
событию безотказной работы является событие отказа

Вероятность
отказа
есть вероятность того, что время появления отказа будет меньше заданного
времени работы изделия, т.е. вероятность того, что время наработки до отказа (t_отк)
будет меньше времени работы (t).

Q(t) = Вер (t_отк<
t)

Статистическая
вероятность времени появления отказа равна: .

С течением
времени наработки число отказавших изделий непрерывно увеличивается.
Следовательно, вероятность отказов является монотонно возрастающей функцией.

Рис.2

Пример
зависимости вероятности возникновения отказа от времени показан на рисунке 2. 
Для работоспособного объекта в начальный момент времени вероятность отказа
близка к нулю. Для того, чтобы отказ проявился, объекту необходимо начать
работать, при этом вероятность отказа увеличивается с увеличением времени и стремится
к единице. 0≤Q(t)≤1

Безотказная
работа изделия и его отказ являются двумя противоположными и несовместимыми
случайными величинами, поэтому их сумма всегда равна 1.

P(t)+Q(t)=1⟹P(t)=1-Q(t)
или Q(t)=1-P(t).

Например: Для
предыдущего примера определить вероятность отказа агрегатов за 500 часов
работы.

Решение:

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

или

Таким
образом, вероятность отказа агрегата за 500 часов составляет 23,2 %.

4.     ОЦЕНКА
ВЕРОЯТНОСТИ БЕЗОТКАЗНОЙ РАБОТЫ.

Рассмотрим график наблюдения
за десятью однотипными изделиями в течение времени от 0 до t₄. Здесь
сплошной прямой линией показана продолжительность безотказной работы изделия, а
крестиком – момент возникновения отказа.

Для наглядности
разместим наработки до отказов изделий последовательно по времени их появления.

Определение
технического состояния изделий в процессе испытаний производится в моменты
времени t₁… t₄. Оценки вероятностей безотказной работы за
соответствующие интервалы времени будут иметь вид:

По
полученным данным строится ступенчатый график – гистограмма, в конце каждого
интервала времени наблюдаемое значение вероятности в данном случае снижается на
долю изделий, отказавших на данном интервале.

Полученные
значения показывают приблизительно долю изделий, которые проработают безотказно
при испытаниях другой партии таких же изделий в аналогичных условиях.

Например:
Изготовив 20 новых изделий, можно утверждать, что в течение времени t₃
приблизительно 12 изделий проработают безотказно (не проводя дополнительных
испытаний) 20*0,6=12. Это приближенная оценка будет тем точнее, чем больше
число испытанных изделий.

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

Например: Условная
вероятность безотказной работы изделия на интервале времени от t₂ до
t₃ оценивается согласно выражению:

Как вычисляется среднее время до отказа и вероятность безотказной работы?

Понятиям MTTF (Mean Time To Failure — среднее время до отказа) и другим терминам теории надежности посвящено большое количество статей, в том числе на Хабре (см., например, тут). Вместе с тем, редкие публикации «для широкого круга читателей» затрагивают вопросы математической статистики, и уж тем более они не дают ответа на вопрос о принципах расчета надежности электронной аппаратуры по известным характеристикам ее составных элементов.

В последнее время мне довольно много приходится работать с расчетами надежности и рисков, и в этой статье я постараюсь восполнить этот пробел, отталкиваясь от своего предыдущего материала (из цикла о машинном обучении) о пуассоновском случайном процессе и подкрепляя текст вычислениями в Mathcad Express, повторить которые вы сможете скачав этот редактор (подробно о нем тут, обратите внимание, что нужна последняя версия 3.1, как и для цикла по machine learning). Сами маткадовские расчеты лежат здесь (вместе с XPS- копией).

1. Теория: основные характеристики отказоустойчивости
Вроде бы, из самого определения (Mean Time To Failure) понятен его смысл: сколько (конечно, в среднем, поскольку подход вероятностный) прослужит изделие. Но на практике такой параметр не очень полезен. Действительно, информация о том, что среднее время до отказа жесткого диска составляет полмиллиона часов, может поставить в тупик. Гораздо информативнее другой параметр: вероятность поломки или вероятность безотказной работы (ВБР) за определенный период (например, за год).

Для того чтобы разобраться в том, как связаны эти параметры, и как, зная MTTF, вычислить ВБР и вероятности отказа, вспомним некоторые сведения из математической статистики.

Ключевое понятие теории надежности — это понятие отказа, измеряемое, соответственно, интервальным показателем
Q(t) = вероятность того, что изделие откажет к моменту времени t.
Соотвественно, вероятность безотказной работы (ВБР, в английской терминологии «reliability»):
P(t) = вероятность того, что изделие проработает без отказа от момента t₀=0 до момента времени t.
По определению, в момент t₀=0 изделие находится в работоспособном состоянии, т.е. Q(0)=0, а P(0)=1.

Оба параметра — это интервальные характеристики отказоустойчивости, т.к. речь идет о вероятности отказа (или наоборот, безотказной работы) на интервале (0,t). Если отказ рассматривать, как случайное событие, то, очевидно, что Q(t) — это, по определению, его функция распределения. А точечную характеристику можно определить, как
p(t)=dQ(t)/dt = плотность вероятности, т.е. значение p(t)dt равно вероятности, что отказ произойдет в малой окрестности dt момента времени t.

И, наконец, самая важная (с практической точки зрения) характеристика: λ(t)=p(t)/P(t)=интенсивность отказов.
Это (внимание!) условная плотность вероятности, т.е. плотность вероятности возникновения отказа в момент времени t при условии, что до этого рассматриваемого момента времени t изделие работало безотказно.

Измерить параметр λ(t) экспериментально можно путём испытания партии изделий. Если к моменту времени t работоспособность сохранило N изделий, то за оценку λ(t) можно принять процент отказов в единицу времени, происходящих в окрестности t. Точнее, если в период от t до t+dt откажет n изделий, то интенсивность отказов будет примерно равна
λ(t)=n/(N*dt).

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

А вот при чем.

2. Экспоненциальное распределение
В терминологии, которую мы только что использовали, пока не было никаких предположений о свойствах случайной величины — момента времени, в который происходит отказ изделия. Давайте теперь конкретизируем функцию распределения значения отказа, выбрав в качестве нее экспоненциальную функцию с единственным параметром λ=const (смысл которого будет ясен через несколько предложений).

Дифференцируя Q(t), получим выражение для плотности вероятности экспоненциального распределения:
,
а из него – функцию интенсивности отказов: λ(t)=p(t)/P(t)=const=λ.

Что мы получили? Что для экспоненциального распределения интенсивность отказов – есть величина постоянная, причем совпадающая с параметром распределения. Этот параметр и является главным показателем отказоустойчивости и его часто так и называют λ-характеристикой.

Мало того, если теперь посчитать среднее время до первого отказа – тот самый параметр MTTF (Mean Time To Failure), то мы получим, что он равен MTTF=1/ λ.

Все это замечательные свойства экспоненциального распределения. Почему мы выбрали в качестве для описания отказов именно его? Да потому что это наиболее простая модель – модель пуассоновского потока событий, которая уже была нами рассмотрена в статье про анализ конверсии сайта. Поэтому-то в теории надежности наиболее часто используется показательное (экспоненциальное) распределение, для которого, как мы выяснили:

• надежность элементов можно оценить одним числом, т.к. λ=const;
• по известной λ довольно просто оценить остальные показатели надежности (например, ВБР для любого времени t);
• λ обладает хорошей наглядностью
• λ нетрудно измерить экспериментально

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

Содержание

вероятность отказ наработка интенсивность

1. Определение вероятности отказов и вероятности безотказной работы устройств

2. Расчет средней наработки до отказа

3. Расчет интенсивности отказов

4. Расчет вероятности безотказной работы системы, состоящей из двух подсистем

Библиографический список

1. Определение вероятности отказов и вероятности безотказной работы устройств

Данные:

Массив значений наработки до отказа Т, 103ч.

10,15,7,9,6,11,13,4,15,12,12,8,5,14,

8,11,12,8,10,11,15,6,7,9,10,14,7,11,

13,5,9,8,9,15,10,9,12,14,10,12,11,8,

10,12,11,12,10,11,7,9.

Заданное значение t = 12,5Ч103ч.

Значение Т0 =4,5Ч103ч.

Объем партии = 300

Значение к = 5.

Решение:

Статистически вероятность безотказной работы устройства для наработки t определяется следующим образом по формуле (1.1)

(1.1)

где – число объектов, работоспособных на момент времени t;

;

– число объектов, неработоспособных к наработке t;

.

Вероятность отказа устройства за наработку t статистически определяется по формуле (1.2)

, (1.2)

;

Проверка:

0,18+0,82=1.

Оценка вероятности безотказной работы устройства по первым 20-ти значениям наработки до отказа обозначается как P*(t). Ее значение определяется также по формуле (1.1), но при этом N=20, и число работоспособных объектов (t) выбирается из этой совокупности.

Будем считать, что условия опыта, включающего 50 наблюдений, позволили однозначно определить вероятность безотказной работы устройства, т.е.

P(t)=1 – F(t),

где F(t) – функция распределения случайной величины “наработка до отказа”, определяющая вероятность событий .T? t при N>?.

Тогда с учетом формулы (1.1) математическое ожидание числа объектов

– работоспособных к наработке t, определяется по формуле (1.3)

, (1.3)

где N- объем партии устройств = 300;

.

2. Расчет средней наработки до отказа

Для вычисления среднего значения случайной величины Т непосредственно по ее выборочным значениям t1, t2, …..,ti,….,tN использую формулу (2.1)

, (2.1)

.

Для выполнения данного задания принимаю ?t = 3Ч103ч., а интервал m = 4.

Результаты заношу в таблицу 2.1

Таблица 2.1 – Преобразование значений наработки до отказа в статистический ряд

интервал	Число попаданий на интервал	n	Статистическая вероятность,q
№	Нижняя и верхняя границы , 103ч
1	3,5 6,5	/////	n1=5	q1= 0.1
2	6,59,5	///////////////	n2=15	q2= 0.3
3	9,5 ч 12,5	/////////////////////	n3=21	q3= 0.42
4	12,5 ч 15,5	/////////	n4=9	q4= 0.18

Статистический ряд отображаю на гистограмме (2.1).

Гистограмма 2.1

Определяю правильность расчетов используя соотношение по формуле (2.2)

, (2.2)

.

Статическую вероятность попадания случайной величины qi на i-ый интервал рассчитываю по формуле (2.3)

, (2.3)

Проверяю правильность расчетов по формуле (2.4)

, (2.4)

.

Среднее значение наработки на отказ определяю по формуле (2.5)

, (2.5)

где ;

.

Оцениваю ошибку в расчетах по формуле (2.6)

, (2.6)

3. Расчет интенсивности отказов

Рассчитаю интенсивность отказов для заданных значений t и t

Подсистема управления включает в себя k последовательно соединенных блоков (Рис.3.1).

Рисунок 3.1 – схема соединения электронных блоков

Интенсивность отказов рассчитываю по формуле (3.1).

=, (3.1)

где – статистическая вероятность отказа устройства на интервале (t, t +Дt)

P(t)-вероятность безотказной работы устройства;

Дt = 3·103 ч. принятый ранее в работу интервал наблюдения;

Определяю статистическую вероятность отказа устройства на заданном интервале (12,5·103ч ) из таблицы (2.1) и нахожу интенсивность отказов;

.

При условии, что интенсивность отказов не меняется в течении всего срока службы объекта, т.е. л=const,то наработка до отказа распределена по экспоненциальному закону и вероятность безотказной работы блока в этом случае определяется по формуле (3.2)

, (3.2)

А средняя наработка блока до отказа определяется по формуле (3.3)

, (3.3)

=1287 ч.

Интенсивность отказов подсистемы лП(t), образованной из k-последовательно включенных блоков, нахожу по формуле (3.4)

, (3.4)

Так как все блоки имеют одинаковую систему отказов, то определяю по формуле (3.5)

(3.5)

.

Вероятность безотказной работы подсистемы определяю согласно формуле (3.6)

, (3.6)

Среднюю наработку на отказ подсистемы определяю аналогично по формуле (3.3)

Результаты расчета зависимостей вероятностей безотказной работы одного блока и подсистемы от наработки заношу в таблицу 3.2

Таблица 3.2

t,ч	0	400	800	1200	1600	2000	2400
	1,000	0,733	0,537	0,393	0,288	0,211	0,154
	1,000	0,211	0,044	0,009	0,002	0,0004	0,00009
t,ч	2800	3200	3600	4000	4400	4800	5200
	0,113	0,083	0,060	0,044	0,032	0,024	0,017
	0,00001	0,0000	0,0000	0,00000	0,00000	0,0000	0,0000

Строю график зависимостей и

Рисунок 3.1 – График зависимостей и .

Для любого распределения наработки на отказ вероятность безотказной работы подсистемы, состоящей из k-последовательно соединенных блоков, связана с вероятностями безотказной работы этих блоков соотношением по формуле (3.7)

, (3.7)

Если блоки равно надежны, то вероятность безотказной работы подсистемы определяю по формуле (3.8)

, (3.8)

Рассчитываю вероятность безотказной работы подсистемы при наработке, равной по формулам (3.6) и (3.8) и сравниваю результаты:

=0,367;

=0,367.

Результаты расчета по обеим формулам одинаковы.

4. Расчет вероятности безотказной работы системы, состоящей из двух подсистем

Для наработки t= требуется рассчитать вероятность безотказной работы системы ,состоящей из двух подсистем (рис. 4.1), одна из которых является резервной.

Рисунок 4.1 – Схема резервирования элементов

Решение:

Расчет произвожу в предположении, что отказы каждой из двух систем независимы, т.е. отказ первой не нарушает работоспособность второй, и наоборот.

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

, тогда вероятность отказа одной подсистемы определяю по формуле (4.1)

, (4.1)

Вероятность отказа всей системы определяется из условия, что отказали первая и вторая подсистемы по формуле (4.2)

, (4.2)

Тогда вероятность безотказной работы системы определится формулой (4.3)

, (4.3)

Или иначе по формуле (4.4)

, (4.4)

.

Библиографический список

1. Технология производства и ремонта вагонов; под ред. К.В. Мотовилова. – 2-е изд., перераб. и доп. – М.: Маршрут, 2003. – 382 с.

2. Технология вагоностроения и ремонта вагонов; под ред. В.С. Герасимова. – М.: Транспорт, 1988. – 331 с.

3. Быков Б.В. Технология ремонта вагонов / Б.В. Быков, В.Е. Пигарев. – М.: Желдориздат, 2001. – 560 с.

4. Батюшин Т.К. Технология вагоностроения, ремонта и надёжность вагонов / Т.К. Батюшин, Д.В. Быховский, В.С. Лукашук. – М.: Машиностроение, 1990. – 360 с.

5. Воробьёв Л.Н. Технология машиностроения и ремонта машин / Л.Н. Воробьёв. – М.: Высш. шк., 1981. – 344 с.

6. Технология машиностроения. В 2 т. Т. 1. Основы технологии машиностроения; под ред. А.М. Дальского. – М.: Издательство МГТУ, 1997. – 560 с.

7. Колесов И.М. Основы технологии машиностроения / И.М. Колесов. – 2-е изд., перераб и доп. – М.: Высш. шк., 1999. – 591 с.