Надежность элемента, работающего до первого отказа

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

Если t случайная величина и известен ее закон распределения, то вероятность отказа аналогичных элементов определится из формулы

Q(t)=P{t < t}, (3.2)

где P{t < t} — вероятность отказа элемента до момента t, определенная на множестве элементов, срок жизни которых менее рассматриваемого момента.

Предполагается, что Q(t) непрерывна, и существует непрерывная плотность вероятности отказа, которая выразится первой производной

q(t)=Q'(t) (3.3)

где q(t) — плотность вероятности отказа.

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

Функция Q(t) полностью определяет надежность элемента. Однако это определение является формальным и не отражает физической сущности старения. Наряду с функцией (3.2) применяется и другая функция

Р(t)=1-Q(t)=Р{t > t}. (3.4)

Формула (3.4) определяет среднюю вероятность безотказной работы за время t . Эту функцию чаще всего называют функцией надежности. Примерный вид функций надежности и ненадежности показан на рис. 3.2.

Надежность элемента, работающего до первого отказа - student2.ru

Рис. 3.2. Вид функций: надежности – P(t), и ненадежности – Q(t)

Функция монотонно убывает Р(0) = 1, Р(t) ® 0, при t ® ¥; Q = 0 и Q(t) ® 1 при t® ¥.

Функцию Р(t) можно приближенно найти из опыта. На испытание ставится N элементов и испытывается в одинаковых условиях в течение времени t = t0. Если к моменту t0 не отказало п элементов, то по теореме Бореля

Надежность элемента, работающего до первого отказа - student2.ru при N®¥

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

Надежность элемента, работающего до первого отказа - student2.ru (3.5)

Можно функцию Р(t) найти для всех значений t£ t0. Для этого при проведении испытаний необходимо отмечать истинное время отказов элементов. По моментам отказов можно определить функцию п(t), которая численно равна числу элементов, не отказавших к моменту t. В начальный момент (при t=0) п(t)=N, а в момент каждого отказа уменьшается на единицу.

Тогда формулу (3.5) можно записать следующим образом:

Надежность элемента, работающего до первого отказа - student2.ru (3.6)

Отношение (3.6) называется эмпирической функцией надежности. С ростом N эта функция равномерно приближается к Р(t) и для больших N

Надежность элемента, работающего до первого отказа - student2.ru (3.7)

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

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

Надежность элемента, работающего до первого отказа - student2.ru (3.8)

После преобразования:

Надежность элемента, работающего до первого отказа - student2.ru (3.9)

Формула (3.9) показывает, что Т0 (среднее время безотказной работы) геометрически выражается площадью фигуры, ограниченной осями координат и кривой Р(t) (рис. 3.2).

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

Надежность элемента, работающего до первого отказа - student2.ru (3.10)

где:

Надежность элемента, работающего до первого отказа - student2.ru - среднее время жизни элемента;

τ1,2 - время жизни i-го элемента;

N – число элементов, подвергнутых испытанию.

Согласно закону больших чисел с вероятностью единица

`t®T0 при N® ¥

при большом N

Надежность элемента, работающего до первого отказа - student2.ru »T0 (3.11)

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

Надежность элемента, работающего до первого отказа - student2.ru (3.12)

n – число отказавших элементов;

T0 - среднее время безотказной работы;

ti – время жизни элемента.

Этой оценкой можно пользоваться при n, близком к N.

Вкратце упомянем о некоторых других характеристиках надежности.

1. Дисперсия времени жизни.

Теоретическая формула

Надежность элемента, работающего до первого отказа - student2.ru (3.13)

Опытная формула

Надежность элемента, работающего до первого отказа - student2.ru (3.14)

где

Надежность элемента, работающего до первого отказа - student2.ru

2. Среднее квадратичное отклонение.

Надежность элемента, работающего до первого отказа - student2.ru (3.15)

Формула (3.15) дает среднее квадратичное отклонение времени от среднего Т0. Использование дисперсии как характеристики надежности целесообразно тогда, когда s< Т0, то есть случайное время t имеет небольшой относительный разброс.

В этом случае яркое представление о надежности дает график плотности q(t), который, как правило, имеет один максимум (рис. 3.3а).

Надежность элемента, работающего до первого отказа - student2.ru

Рис. 3.3 Гистограммы плотности q(t): а) теоретическая: б) эмпирическая

Функция q(t) может быть определена из опыта. Участок оси разбивается на интервалы (рис. 3.3б).

Пусть поставлено на испытания N элементов и произошло nk отказов на интервале [(k—1)h; kh] k=1, 2, 3, ...и т.д. (целое число).

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

Надежность элемента, работающего до первого отказа - student2.ru (3.16)

если

Надежность элемента, работающего до первого отказа - student2.ru

При большом объеме испытаний и достаточно малом h

Надежность элемента, работающего до первого отказа - student2.ru (3.17)

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

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

Надежность элемента, работающего до первого отказа - student2.ru (3.18)

где

l(t) – функция интенсивности отказа;

P¢(t) – первая производная от функции надежности;

P(t) – функция надежности.

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

Применительно к теории надежности l(t) , это вероятность того, что элемент, проработавший безотказно до момента t, откажет в последующую единицу времени (при условии, что эта единица времени будет небольшой). Формула (3.18) разрешается относительно функции надежности

Надежность элемента, работающего до первого отказа - student2.ru (3.19)

Уравнение (3.19) является одним из основополагающих в стохастической теории надежности. Из него следует, что вероятность безотказной работы элемента на участке (t1;t2) будет равна

Надежность элемента, работающего до первого отказа - student2.ru

Функция интенсивности отказа может быть определена опытным путем по результатам испытаний.

Испытываются N элементов и регистрируются их отказы. Пусть n(t) -число элементов, не отказывающих к моменту t. Тогда при достаточно малом Dt и достаточно большом N

Надежность элемента, работающего до первого отказа - student2.ru

где Dn- число отказов на участке (t; t+Dt);

n(t) – функция отказов, численно равная числу элементов, не отказавших к моменту t.

Итак, если Dt достаточно мало, а Dn при этом велико, то

Надежность элемента, работающего до первого отказа - student2.ru (3.20)

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

Рассмотрим один пример.

На испытания ставятся 1000 элементов. В первый час отказывает 50 элементов. В последующие 100 часов отказывают 890 элементов, причем отказы следуют равномерно во времени. В следующий (после 100 часов) отказывают еще 20 элементов. Спрашивается: когда элемент имеет наивысшую надежность?

Абсолютные цифры не могут характеризовать этого явления. Очевидно, что наиболее надежный элемент имеет минимальную интенсивность отказов. По формуле (3.20)

а) За 1-й час

Надежность элемента, работающего до первого отказа - student2.ru ,

б) за последующие 100 часов

Надежность элемента, работающего до первого отказа - student2.ru ,

в) за последний час испытаний

Надежность элемента, работающего до первого отказа - student2.ru .

Практически кривую интенсивности отказов можно построить по методике построения кривой q(t).

Для этого:

а) ось абсцисс (времени) разбивается на интервалы длительности h,

б) из испытаний определяется число отказов, попавших в каждый интервал hk .

Если пk число отказов на k-ом интервале [(k—1)h; kh], где k=1, 2, 3 ..., то эмпирическая интенсивность отказов на этом интервале будет равна

Надежность элемента, работающего до первого отказа - student2.ru (3.21)

при

(k-1) h £ t £kh.

В формуле 3.21

N—число элементов, подвергнутых испытанию;

ni —число элементов, отказавших на предыдущих интервалах.

Подсчитав l(t) для всех интервалов, можно по полученным значениям построить график lN(t). Примерный вид построения функции lN(t) приведен на рис. 3.4.

Надежность элемента, работающего до первого отказа - student2.ru

Рис. 3.4. Порядок построения функций lN(t) по опытным данным

В общем случае кривая l(t) может иметь различный вид. Рассмотрим наиболее простой и наиболее важный частный случай кривой интенсивности отказа (рис. 3.5).

Надежность элемента, работающего до первого отказа - student2.ru

Рис. 3.5. Идеализированный вид кривой интенсивности отказов l(t).

Весь интервал по этой кривой можно разбить на 3 участка.

I—период приработки, «выжигания» дефектных элементов. В этот период выходят из строя некачественно изготовленные элементы.

II—период нормальной работы элементов.

l(t)= const или l(t)» const,

III—период форсированного «старения» в основном убыстряющегося; l(t) возрастает.

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

Поэтому во многих случаях можно принять, что

l(t)= l=const (3.22)

Тогда функция надежности (формула 3.19) примет следующий вид:

Р(t)=еxp[-lt]. (3.23)

Формула (3.23) является выражением экспоненциального закона надежности. Для него вероятность отказа за время t равна

Q(t)=1-е-lt, (3.24)

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

q(t)=lе-lt (3.25)

Среднее время жизни для этого случая

Надежность элемента, работающего до первого отказа - student2.ru (3.26)

Для экспоненциального закона интенсивность отказа обратна среднему времени. Поэтому функцию надежности можно записать и так

Р(t)=ехр Надежность элемента, работающего до первого отказа - student2.ru (3.27)

Часто интересующее нас время во много раз меньше среднего времени

t<<T0 или Надежность элемента, работающего до первого отказа - student2.ru <<1.

Тогда можно пользоваться приближенными формулами

Надежность элемента, работающего до первого отказа - student2.ru (3.28)

Надежность элемента, работающего до первого отказа - student2.ru . (3.29)

Ошибка в формулах (3.28) и (3.29) не превосходит Надежность элемента, работающего до первого отказа - student2.ru Экспоненциальный закон очень популярен в теории надежности. Это объясняется тем, что он имеет физически естественное толкование, прост и удобен для математических расчетов.

Экспоненциальный закон обладает следующим важным свойством (по Б.В. Гнеденко).

«Вероятность безотказной работы на данном интервале (t; t+t) не зависит от времени предшествующей работы t, а зависит только от длины интервала t». Популярно выражаясь, можно сказать, если элемент в данный момент исправен, то его будущая работа (безотказность) не зависит от прошлой.

Это утверждение имеет строгое математическое доказательство, с которым можно познакомиться в работах Б.В. Гнеденко.

Однако именно за это утверждение экспоненциальный закон надежности подвергается критике. Действительно, утверждение о независимости безотказной работы элемента на участке времени t; t+t от предыдущей работы элемента, не согласуется с существующими представлениями о физике старения элементов машин и деталей. Это свойство, очевидно, следует из принятого допущения о случайной природе отказов элементов и систем.

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

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

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

Нормальный закон распределения, по мнению ряда исследователей, хорошо описывает процессы, протекающие с постепенными (износовыми) отказами.

При этом предполагается, что участок tmin (рис. 3.6 а), который характеризует минимальный срок жизни элемента, установленный из опыта, представляет часть времени жизни элемента с вероятностью безотказной работы Р(t)=1 (рис. 3.6 б). Поэтому характеристики надежности могут рассчитываться в интервале времени t1; t2. Вероятность безотказной работы будет равна

Надежность элемента, работающего до первого отказа - student2.ru , (3.30)

при

t1£ t £ t2

в формуле (3.30)

Т0 — среднее время жизни элемента, определяемое по формуле

Надежность элемента, работающего до первого отказа - student2.ru (3.31)

или по формуле (3.10)

s — среднее квадратическое отклонение, равное

Надежность элемента, работающего до первого отказа - student2.ru (3.32)

Надежность элемента, работающего до первого отказа - student2.ru

Рис. 3.6. Графики нормального закона распределения

а) кривая плотности распределения q(t)

б) надежность - Р(t) и ненадежность - Q(t)

Для вычисления значения вероятности безотказной работы пользуются табулированной функцией Лапласа

Ф(z)= Надежность элемента, работающего до первого отказа - student2.ru (3.33)

где z = Надежность элемента, работающего до первого отказа - student2.ru .

Решение формулы (3.33) представляется в виде

P(t)=0,5+Ф Надежность элемента, работающего до первого отказа - student2.ru (3.34)

В формуле (3.34) применена нормированная функция Лапласа, которая изменяется в пределах

0 £ Ф(z) £0,5.

У реального элемента нередко совмещаются оба типа отказов: постепенный и внезапный. Надежность такого элемента можно определить, исходя из следующих соображений. Если P1(t)вероятность того, что за время t не произойдет внезапного отказа, а P2(t) —то же, но для постепенного отказа, то полагая:

а) что элемент работает до первого из этих отказов:

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

P(t)= Надежность элемента, работающего до первого отказа - student2.ru (3.35)

Среднее время жизни для этого случая (при s << Т0, а Т0 и Надежность элемента, работающего до первого отказа - student2.ru одного порядка)

T Надежность элемента, работающего до первого отказа - student2.ru (3.36)

Закон Вейбулла. Функция надежности имеет для него следующий вид

Р(t)= Надежность элемента, работающего до первого отказа - student2.ru . (3.37)

Среднее время жизни элемента

T 0 = Надежность элемента, работающего до первого отказа - student2.ru (3.38)

a — дополнительный параметр.

Дисперсия времени жизни

Dt= Надежность элемента, работающего до первого отказа - student2.ru Надежность элемента, работающего до первого отказа - student2.ru (3.39)

Г Надежность элемента, работающего до первого отказа - student2.ru - табулированная гамма-функция. Значения ее приводятся в таблицах. Интенсивность отказа

l(t)=alt Надежность элемента, работающего до первого отказа - student2.ru (3.40)

Если Надежность элемента, работающего до первого отказа - student2.ru , то функция Вейбулла:

Р(t) Надежность элемента, работающего до первого отказа - student2.ru (3.41)

Вид кривых l(t) при различных значениях a приведен на рис. 3.7.

Надежность элемента, работающего до первого отказа - student2.ru

Рис. 3.7. Вид кривых l(t) при различных значениях a.

При a > 1 интенсивность отказа монотонно возрастает от 0; при a < 1 интенсивность отказа убывает монотонно и не ограничена при t = 0.

Экспоненциальный закон является частным случаем закона Вейбулла при a = 1. Закон Вейбулла получил широкое распространение, так как содержит дополнительный параметр a, подбирая который можно получить лучшее соответствие функции надежности, экспериментальным данным.

ЛИТЕРАТУРА

В.А. Наумов Основы надежности и долговечности в машиностроении. Омск: 1972. 331с.

Наши рекомендации