что такое вероятность безотказной работы

21.12.202323.09.2023 admin 0 Comments

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

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

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

P*( t )=P<Т>t>

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

где : N₀ – число образцов аппаратуры в начале испытаний;

n(t) –число отказавших образцов за время t.

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

Точность определения Р* (t) зависит от числа N₀. Функция вероятности безотказной работы Р (t) является не возрастающей функцией времени и обладав следующими очевидными свойствами:

в период от до ;

Считаем, что изделия при испытании (эксплуатации) не восстанавливаются и

не заменяются новыми.

Типичная кривая изменения частоты отказов изделий в соответствии с

зависимостью f * (t) показана на рис. 2.2. На этой кривой можно отметить три характерных участка.

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

К этой группе отказов можно отнести также эксплуатационные отказы, вызванные слабым знанием правил эксплуатации (или отсутствием необходимого опыта). Поэтому первый период называется периодом приработки изделий (элементов).

Второй участок характерен сравнительно постоянным значением частоты отказов и называется периодом нормальной эксплуатации.

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

Длительность вышерассмотренных участков различна: для первого она составляет величину порядка нескольких десятков, нескольких сотен часов в зависимости от сложности изделия; второй участок самый продолжительный и при расчетах принимается равным 2/3 технического ресурса; третий период зачастую не достигается, так как раньше наступает «моральное старение», при котором данный ТУ списывается.

Рисунок 2.2. Кривая изменения частоты отказов от наработки

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

где — среднее число элементов, продолжающих исправно работать на интервале Dt.

Рисунок 2.3. Зависимость интенсивности отказов от наработки

Число определяют с помощью выражения:

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

Зависимость l* от t представлена на рисунке 2.3.

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

Источник

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

Понятиям 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/ λ.

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

Источник

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

Вероятностью безотказной работы называется вероятность того, что при определенных условиях эксплуатации, в пределах заданной наработки не произойдет ни одного отказа.
Вероятность безотказной работы обозначается как P(l), которая определяется по формуле (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 тыс. км., в локомотивном депо произведен сбор информации по отказам ТЭД. При этом количество исправных ТЭД в начале периода эксплуатации составляло N 0 = 180 шт. Суммарное количество отказавших ТЭД за анализируемый период составило ∑r(600000) = 60. Интервал пробега принять равным 100 тыс. км. При этом количество отказавших ТЭД по каждому участку составило: 2, 12, 16, 10, 14, 6.

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

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

Первоначально по уравнению (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.5.). Необходимо отметить, что первой точкой на графике, т.е. при пробеге равном 0, величина вероятности безотказной работы примет максимальное значение – 1.

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

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

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

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

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

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

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

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

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

Исходные данные.
В результате испытания десяти топливных насосов высокого давления получены наработки их до отказа: 400, 440, 500, 600, 670, 700, 800, 1200, 1600, 1800 ч. Предполагая, что наработка до отказа топливных насосов подчиняется экспоненциальному закону распределения.

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

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

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

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

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

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

Исходные данные.
Время простоя тепловозов в неплановых ремонтах по вине вспомогательного оборудования подчиняется закону распределения Вейбулла-Гнеденко с параметрами b=2 и a=46.

Требуется.
Необходимо определить вероятность выхода тепловозов из неплановых ремонтов после 24 ч. простоя и время простоя, в течение которого работоспособность будет восстановлена с вероятностью 0,95.

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

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

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

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

Исходные данные.
Известно, что наработки контакторов до отказа по параметрам старения изоляции катушек можно описать функцией распределения Рэлея с параметром S = 260 тыс.км.

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

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

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

Исходные данные.
Нерезервированная система состоит из 5 элементов. Интенсивности их отказов соответственно равны 0,00007; 0,00005; 0,00004; 0,00006; 0,00004 ч-1

Требуется.
Необходимо определить показатели надежности системы: интенсивность отказов, среднее время наработки до отказа, вероятность безотказной работы, частота отказов. Показатели надежности P(l) и a(l) получить в интервале от 0 до 1000 часов с шагом в 100 часов.

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

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

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

Графическая иллюстрация 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). Далее рассчитаем надежность каждого блока в отдельности, а затем найдем надежность всей системы.