Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математическое моделирование процессов в машиностроении..pdf
Скачиваний:
46
Добавлен:
15.11.2022
Размер:
14.87 Mб
Скачать

U - система подмножеств, выделенная в Q; р - числовая функ­

ция, определенная на U и

обладающая свойствами р(А) > О,

p(fl) = 1 для произвольной

последовательности попарно несо­

вместимых событий А„ 1,2,... (А[ГлА) = 0 при / * j, А,е U). Функ­ ция р называется вероятностью.

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

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

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

Рассмотрим единичные показатели (рис. 4.8), которые иг­ рают важную роль в проектных расчетах и отражают основные характеристики надежности технических систем.

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

1.Вероятность безотказной работы (ВБР) - это вероят­ ность того, что в пределах заданной наработки отказ объекта не возникнет. ВБР является основной количественной характе­ ристикой безотказности объекта на заданном временном интер­ вале и наиболее полно определяет его надежность. Если обозна­ чить через Т время непрерывной исправной работы объекта от начала работы до первого отказа, а через t - время, за которое необходимо определить ВБР, то ВБР записывается в виде

P(t) = P{T> t},t> 0.

(4.21)

Рис. 4.8. Единичные показатели надежности

Функция ВБР обладает следующими очевидными свойст­

вами:

1)0 </>(/) < 1;

2) функция P(t) есть невозрастающая функция своего ар­

гумента, т.е. если h> t\, то P(h) < P(t\);

3 ) Р(0)=1,Д<ю) = 0.

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

Статистически ВБР равна

 

 

IIAI

 

 

P(t) =

lim -

/=1__^

m

(4.22)

v 7

Д/—>0

дг

Nn ‘

 

 

Nt)->со

0

 

 

где N0 - число объектов в начале испытаний; и, -

число отка­

завших объектов в интервале времени A/,; t - время, для которо­

го определяется ВБР; N(t) - число объектов, исправно работаю­ щих на интервале [0, О-

Вероятность того, что отказ объекта произойдет за время,

не превышающее заданной величины t, т.е. что Т < t, равна

 

6(0 = Q{T <t} = \ - Р(0, t > о.

(423)

Функция 6 (0 представляет собой интегральную функцию

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

Если

функция 6 (0 дифференцируема, то ее производная есть диффе­ ренциальный закон (плотность) распределения случайной вели­ чины t - времени исправной работы. Тогда можно записать:

« W - Л О и М - - Ш .

(4.24)

dt

at

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

 

 

 

2 > .

N0 ~N(t)

 

 

(4.25)

г

Nn

м,->=о Jvo

 

 

Графически функции P(t), Q(t) и /0 представлены на рис. 4.9.

Рис. 4.9. Графики функций:

1 - вероятности безотказной работы;

2 - вероятности отказа;

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

Отметим некоторые общие свойства функции распределеиия Q(t), которые, как следует из выражений (4.24), противопо­ ложны свойствам функции распределения ВБР P{t):

1)0 < 2 ( 0 ^ 1;

2)функция распределения Q{t) есть неубьтающая функЦия своего аргумента, т.е. при /2 > t\ П Ь ) > Ш ,

3)в начальный момент времени функция распределения

равна нулю: (2(0) = 0,

и в

бесконечности

она равна единице

1.

 

 

 

 

 

 

Проиллюстрируем

пред­

о

 

 

 

ставленные свойства.

Для

этого

т

t0

t

будем рассматривать величину tQ

 

 

 

 

как случайную точку

Т на оси /

 

m

 

 

(рис. 4.10), которая в

результате

Рис. 4.10. Иллюстрация

опыта может занять то или иное

свойств функции Q(t)

 

положение.

 

 

 

 

 

 

Тогда функция

распределе­

 

 

 

 

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

Будем увеличивать t, т.е. перемещать точку to вправо по оси абсцисс. Очевидно, при этом вероятность того, что случайная точка Т попадает левее to, не может уменьшиться, следовательно, функция распределения Q(t) с возрастанием t убывать не может.

Из введенного в определение вероятности безотказной ра­ боты предположения, что в начальный момент времени исчисле­ ния наработки объект был работоспособен, следует, что Q(0) = 0.

Плотность вероятности fit) статистически определяется по

формуле

 

 

Д Л ') =

An(At)

(4.26)

N0At

 

где An(At) число отказов за интервал времени At.

Очевидно, что

 

(

со

6(0 = J/COrfT И P(t)=

(4-27)

где т - момент времени, предшествующий t.

Поскольку основной параметр функционирования объекта в теории надежности - время t - изменяется в пределах [0; со], то из условия нормировки (4.10) и первого свойства функции Q(t)

следует, что для функцииД/) имеет место равенство

 

оо

 

J/(0dt = 1

(4.28)

2. Средняя наработка до отказа - это момент первого порядка (математическое ожидание) наработки объекта до пер­ вого отказа. Эту величину обозначают Гер и называют также средним временем безотказной работы:

(4.29)

Статистически средняя наработка до отказа однотипных объектов равна

 

I

ср

(4-30)

*'0 7=1

где tj - время исправной работыJ-го объекта.

3. Гамма-процентная наработка до отказа (у - это нара ботка, в течение которой отказ объекта не возникнет с вероятно­ стью у, выраженной в процентах. Гамма-процентная наработка определяется из уравнения

(4-31)

При у =100% гамма-процентная наработка называется установленной безотказной наработкой, при у = 50 % гамма­ процентная наработка называется медианной наработкой.

4. Средняя наработка на отказ —отношение наработки восстанавливаемого объекта к математическому ожиданию числа его отказов в течение этой наработки. Для интервала наработки (<i; /2) средняя наработка на отказ определяется выражением

Т

*2

(4.31)

*л —

 

Q (/,;0

где То - средняя наработка на отказ; Q(/i; /2) - математическое ожидание числа отказов на интервале наработки /2 - Л.

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

Статистически средняя наработка на отказ определяется так:

< « 2)

П /=1

где tcp- время исправной работы между ( / - 1)-м и i-м отказами объекта; п число отказов объекта.

При достаточно большом числе отказов /ср будет стре­ миться к среднему времени между двумя соседними отказами. Если испытания проводятся не с одним, а с несколькими одно­ типными объектами, то среднее время между отказами можно определить из выражения

1 м

(4.33)

Т° = м Р с ’

 

где М - число объектов.

5. Интенсивность отказов - это отношение числа отка­ завших объектов в единицу времени к среднему числу объектов, продолжающих исправно работать в данный интервал времени.

Другими словами, интенсивность отказов - это условная плотность вероятности возникновения отказа невосстанавливаемого объекта, определяемая для рассматриваемого момента времени при условии, что до этого момента отказ еще не возник.

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

4 0 =

g(t;t + At)

(4.34)

 

At

где q(t; t + At) - вероятность^ отказа невосстанавливаемого объ­ екта на интервале времени от t до / + At при условии, что до мо­ мента времени t отказ не возникает.

Определение интенсивности отказов на основе статисти­ ческих данных осуществляют по формуле

& ,

(4.35)

N(t)At

где An(At) - число отказов объекта за период времени ( —— ;

t + At

2

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

Функция Х,(/) для многих объектов имеет характерный вид (рис. 4.11) для разных периодов эксплуатации.

Период I (О < / < Л ) называют периодом приработки, в процессе которого устраняются систематические отказы эле­ ментов (конструктивные, производственные, эксплуатацион­ ные), заменяются элементы, имеющие дефекты. Все это приво­ дит к повышению вероятности безотказной работы и, следова­ тельно, снижению интенсивности отказов.

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

Встречаются случаи, когда на начальной стадии прира­ ботки интенсивность отказов возрастает (участок Г). Такие зави­ симости могут быть описаны в следующем виде:

=

1.

(4-36)

кк

Эта модель описывает смесь элементов с вероятностями безотказной работы Pk(t) и относительными долями р Группа дефектных элементов имеет малую среднюю наработку до отка­ за и быстрое старение.

В периоде II (t\ < t < /2), часто называемом периодом нор­ мальной эксплуатации, интенсивность отказов остается пример­ но постоянной. Для периода нормальной эксплуатации харак­ терно наименьшее количество отказов.

Наконец, в периоде III, называемом периодом старения и износа, интенсивность отказов увеличивается из-за необрати­ мых физико-химических процессов в элементах объекта, свя­ занных с его длительным использованием, т.е. когда вследствие естественных процессов объект или его составные части при­ ближаются к предельному состоянию по условиям физического

ю (0 =

изнашивания. Отказы III периода иногда называют деградации отыми отказами.

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

Параметр потока отказов со(/) характеризует среднее число отказов, ожидаемых на малом интервале времени, и равен

dCljt)

(4.37)

dt

где Q(t) - математическое ожидание числа отказов за время /. Отметим, что средняя наработка до отказа, гамма­

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

Для того чтобы определить показатели, характеризующие надежность восстанавливаемого элемента, необходимо определить модель его функционирования. Простейшей является такая мо­ дель, при которой объект какое-то случайное время Т| работает до первого отказа, фиксируемого достоверно, затем следует мгновен­ ное и полное восстановление свойств объекта в момент Х\ = 1\, по­ сле чего элемент снова работает случайное время х2до второго от­ каза, определяемого достоверно, затем мгновенно восстанавлива­ ется до начального состояния в момент t\ = TJ + х2и т.д.

Для таких объектов моменты отказов на оси суммарной наработки или на оси непрерывного времени образуют поток отказов. В качестве характеристики потока отказов используют «ведущую функцию» Q(T) данного потока - математическое ожидание числа отказов за время t.

Статистически параметр потока отказов можно опреде­ лить по формуле

©(/) = An(At)

(4.38)

N0At

t —At

где An(At) - число отказов объекта за период времени

2 *

t + At

~7 ~ }'

4.4.2.Показатели долговечности, ремонтопригодности

исохраняемости

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

Средний ресурс Тр- это математическое ожидание ресурса.

Гамма-процентный ресурс Т^% - это наработка, в течение которой объект не достигнет предельного состояния с заданной вероятностью у, выраженной в процентах. Гамма-процентный ресурс определяется по формуле, аналогичной (4.31).

Назначенный ресурс Гр„ определяется как суммарная на­ работка объекта, при достижении которой применение по на­ значению должно быть прекращено.

Средний срок службы Тс - это математическое ожидание срока службы.

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

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

Назначенный ресурс (срок службы) - это характеристики, устанавливаемые на основании субъективных или организаци­

онных принципов и поэтому являющиеся косвенными показате­ лями надежности.

Основным показателем ремонтопригодности является ве­ роятность восстановления. Момент восстановления работо­ способности объекта после отказа является случайным событи­ ем. Поэтому интервал времени от момента отказа до момента восстановления является случайной величиной и для характери­ стики ремонтопригодности может быть использована функция распределения этой случайной величины ©. Вероятностью вос­ становления Ря называется вероятность того, что время восста­ новления работоспособного состояния объекта не превысит за­ данного:

PB(t) = P {® < l},0< t.

(4.39)

Функция Рв(0 представляет собой интегральную функцию распределения случайной величины 0. Вероятность невосста-

новления на заданном интервале (, т.е. вероятность

того, что

0 > /, равна

 

6 .(0 = Р{® > t) = 1 - Р Ш о <1.

(4.40)

Плотность вероятности момента восстановления равна:

dP

(4.50)

/«(0 = - ^ -

По аналогии со средней наработкой до отказа момент пер­ вого порядка m\{Q} (математическое ожидание) времени вос­ становления работоспособного состояния объекта называется

средним временем восстановления:

(4.51)

 

0

 

 

II

ft? 1

/""'Ч

(4.52)

 

 

 

о

Важным показателем ремонтопригодности объекта явля­ ется интенсивность восстановления р(/), которая, следуя общей

методологии, аналогична показателю безотказности —интенсив­ ности отказов.

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

4.5. АНАЛИТИЧЕСКИЕ ЗАВИСИМОСТИ МЕЖДУ

ПОКАЗАТЕЛЯМИ НАДЕЖНОСТИ

Между основными показателями надежности существуют аналитические зависимости. Приведем некоторые из них.

1.Зависимость между вероятностью безотказной работы

исредней наработкой до отказа. По выражению (4.29)

TCf = \tf(t)dt-

(4.53)

О

 

Так как/ 0 = Q\t) и Q(f) = 1 -P{t), то, выполнив интегри­

рование по частям, получим

Тср = -\tdP (t) = -tP (tf0 +)p(t)dt

(4.55)

О

О

 

Учитывая свойства функции P(t), получаем равенство ну­ лю первого слагаемого в правой части полученного равенства.

Тогда

Tcp = \P{t)dt.

(4.55)

2. Связь между ВБР и интенсивностью отказов. Согласно свойствам функций P(t) и Q(t) получим:

q(t, t + ДО = 1 - P(t-1 + At).

(4.56)

С другой стороны вероятность безотказной

работы

P(t + At) на интервале времени от 0 до t+At в соответствии с тео­ ремой умножения вероятностей рассчитывается по формуле

 

Pit + At) = P(t)P(t; t + At).

(4.57)

Графическая запись приведенных рассуждений представ­

лена на рис. 4.12.

 

 

 

 

P(t)

Pit,t+At)

 

С учетом (4.34), (4.56) и (4.57)

имеем:

 

 

 

 

 

 

0

t t+At t

 

 

j Pit + At)

 

P it+ At)

 

40 =

Pit)

(4.58)

Puc. 4.12. Иллюстрация

 

 

At

 

P it)-P it + At)

 

теоремы умножения

 

 

 

 

PiOAt

 

вероятностей

 

 

 

 

 

 

 

 

 

 

Переходя к пределу при At -> О,

получаем

 

 

 

 

 

 

 

1

dPjt)

Л О

(4.59)

 

4 0 = Pit)

dt

Pit)

После интегрирования выражения (4.59) при начальном

условии Pi0) = 1 получаем:

 

 

 

 

 

Pit) =e "

 

(4.60)

где т - момент времени, предшествующий t.

Различие между интенсивностью отказов \(t) и плотно­ стью распределения отказов fit) заключается в том, что эле­ мент условной вероятности Xit)dt характеризует вероятность отказа объекта на элементарном участке (/; t + dt), если из­ вестно, что к моменту времени / он не отказал. Элемент веро-

ятностиХО** характеризует вероятность отказа в течение эле­ ментарной наработки на участке (t; t + dt), но при этом неиз­ вестно, исправен ли объект к моменту времени t. Кроме того, из (4.59) следует, что fit) = X(t)P(t), то есть при малых t значе­ ние P{t) близко к единице и, следовательно, fit) « X(t). При

t =0 P(t) - 1

и fit) = Ц/). Во всех остальных случаях, когда

p { t)< \,fit)< m -

3.

Связь между вероятностью безотказной работы, интен­

сивностью отказов и средней наработкой до отказа.

Среднюю наработку до отказа можно вычислить через ин­

тенсивность отказов. Подставив в выражение (4.55) значение (4.60), получим

(4.61)

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

Пусть в момент t = 0 на испытании находится N0объектов. По мере выхода из строя отказавшие объекты заменяются но­ выми (выборка с возмещением), если объекты невосстанавливаемые. Тогда на основании расчета среднего числа отказавших объектов в зависимости от плотности распределения отказов^/)* параметра потока отказов со(/) и длительности рассматриваемого интервала времени можно получить уравнение Вольтерра, свя­ зывающее со(/) иfit):

(4.62)

О

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

Решение уравнения (4.62) показывает, независимо от вида функции f(t), т.е. независимо от закона распределения времени безотказной работы, что параметр потока отказов объектов стремится к постоянной величине - обратной средней наработке на отказ, т.е. справедливо выражение

lim©(0 = -Jr.

(4.63)

т0

 

5. Связь между вероятностью восстановления и интенсив ностью восстановления.

Рассмотрим условную вероятность PB(t/т) того, что вос­ становление работоспособности объекта произойдет на интервале времени следующим за интервалом х после отказа, на ко­ тором еще не удалось восстановить работоспособность:

р, « / *)=

РЛ‘+ *)-РщЮ

(4.64)

1 - а д

 

 

Предел отношения PB(tlx) при t —» 0 есть не что иное, как дифференциальная плотность вероятности восстановления в момент т при условии, что объект не был восстановлен до мо­ мента х:

PB(tk )_ d P B(x)

1___

р(т) = lirn

 

(4.65)

f->О l

dl 1 - P B( T )

Представляя (4.65) в виде

 

 

р(т) = - — 1п[1-Р(т)]

(4.66)

и учитывая, что Рв(0) = 0, получаем, аналогично (4.60),

/

- | и ( т ) А

PB(t) = \ - e "

(4.67)

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]