Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы обеспечения надежности изделий машиностроения

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
10.92 Mб
Скачать

о 1 г J t

Рис. 1.4. Кривая плотности экспо­ ненциального закона распреде­ ления при Х = \

P(t) = l - U + & L - & L + « 1 - x t .

(1.39)

Плотность распределения описывают соотношением вида

/ ( 0 =

(1.40)

Графическое изображение плотности экспоненциального распределения приведено на рис. 1.4.

функция распределения выражается зависимостью

F(t) = l - P ( / ) = l _ e -*'

(1.41)

Математическое ожидание и среднее квадратическое от­ клонение для экспоненциального закона соответственно равны

M(/) = -f, а(0 = ^ .

(1.42)

Таким образом, [Л4(^)= сг(/)] для экспоненциального закона распределения, что является существенным при проверке соот­ ветствия экспериментального распределения выбранному тео­ ретическому экспоненциальному распределению.

Пример 1.3. Наработка на отказ сложной технической си­

стемы подчиняется

экспоненциальному

закону распределения

с параметром Я = 15* 10 5 ч-1.

 

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

/==100 ч, а также найти

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

Р еш ен и е . Вычислим

вероятность

безотказной работы по

формуле (1.39), так

как М = 0,015:

 

/>(/)=

1— М=1 — 15-10-®. 100= 0,985.

21

Среднее значение наработки на отказ найдем по формуле

(1.42):

 

 

M(t) = Т0 = -J- =

«6677 ч.

Л,

 

1.10. ЛОГАРИФМИЧЕСКИ НОРМАЛЬНОЕ РАСПРЕДЕЛЕНИЕ

Логарифмически нормальное распределение явля­ ется двухпараметрическим распределением случайной величины, логарифм которой распределен по нормальному закону. Как распределение положительных величин, оно несколько точнее, чем нормальное. В теориинадежности такое распределение исполь­ зуют для описания наработки до отказа деталей в период на­ ступления усталости материала, а также процессов восстанов­ ления, износовых отказов, наработки на отказ подшипников кача­ ния и электронных ламп или наработки между отказами сложных технических систем типа подвижных установок на базе шасси МАЗ-543 и ЗИЛ-135.

Плотность распределения описывается зависимостью

f(t) = (^kexp[-

v 0 при tcO .

(Int —a)‘

при /j>0;

(1.43)

2а2

 

Параметры а и а оценивают по результатам испытаний с помощью формул:

I

t;

 

a = t = i =

1

(1.44)

 

N

 

где ti — наработка до отказа /-го

изделия; N — число

изделий,

поставленных на испытания,

 

 

i

O n/,- $

(1.45)

»=л/^ i = I

 

Функция распределения имеет вид

 

In/ —а

 

 

dz

(1.46)

или

 

 

 

 

(1.47)

22

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

U P =

Int —a

(1.48)

 

Математическое ожидание и среднее квадратическое отклоне­

ние наработки до отказа соответственно равны

 

 

а2

'

 

M ( t ) = t

2

a(t) = Vexр (2а + а2) exp (а2— l).

(1-49)

Кривые логарифмически нормального распределения изо­ бражены на рис. 1.5.

Пример 1.4. Оценить вероятность безотказной работы колен­ чатого вала двигателя в течение /= 1 0 3 ч, если его наработка распределена по логарифмически нормальному закону с пара­ метрами 1п/о= 11,51, а 2,0.

Р е ш е н и е:

UP =

1п< —lnf„

InlO3 — In/,,

6,91-11,51 = -2 ,3 ;

 

"2--------

 

 

 

Ж/) = ф(цр) =

0,99.

а)

Ю

Рис. 1.5. Кривые логарифмически нормального распределения при различных параметрах а и а:

а _ плотность вероятности /(/); б — вероятность безотказной работы P(t)

23

1.11.РАСПРЕДЕЛЕНИЕ ВЕЙБУЛЛА

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

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

Плотность распределения выражается зависимостью вида

f(t) = a lta~' exp [—М“],

(1.50)

где а — параметр формы кривой распределения; А, — параметр масштаба.

Экспоненциальное распределение является частным случаем распределения Вейбулла при а = 1 .

Функция распределения описывается соотношением

F(t)= 1 - е х р [ - М “].

(1.51)

Функция надежности — величина, обратная функции распре­ деления:

P(t) = ехр[—А,/а].

(1-52)

Математическое ожидание и среднее квадратическое отклоне­ ние соответственно равны

_ г(1+т).

M(t) =

л.'/2

о(0_У:1'.^)7 г=('Щ

V Х,2/а

где Г(а) — гамма-функция.

Гамма-функция для непрерывной величины интегралом

Г ( 1 + а ) = ^ fexp( — /)d/.

(1.53)

описывается

(1.54)

24

Рис. 1.6. Кривые распределения Вейбулла при различных значениях параметров

аи X:

а— плотность распределения; б — вероятность безотказной работы

Для

вычисления значений функции Г(л + а), где

л — целое

число, а

а — дробное

при 2 ^ л ^ 6 ; рекомендуется

применять

формулу

[13]

 

 

Г(га+а) = (га — 1 + а)(л—2 + а ) . ..(1 + а)Г(1 + а). (1.55)

При га> 6 значения

Г(га-[-а) можно находить по формуле

 

 

Г (га+ 1) = л!

(1-56)

Кривые распределения Вейбулла для различных значений параметров а и X изображены на рис. 1.6.

Пример 1.5. Определить вероятность безотказной работы генератора в течение /= 1 0 3 ч, если наработка на отказ описы­ вается распределением Вейбулла с параметрами: а = 2,0 и

То= -I =1,5.10® ч. Ре ше н и е:

P(t) = е ~1'"1= e _l°3'2/(l'6’ie6) = 0,52.

1.12.ГАММА-РАС ПРЕДЕЛЕН l i t

Гамма-распределение занимает важное место в ма­ тематической статистике и теории надежности. Это распределе­ ние имеет ограничение с одной стороны Если пара­ метр формы кривой распределения а — целое число, то гамма-

25

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

Плотность вероятности гамма-распределения определяется равенствами, если А,>0 и а > 0 :

f(t) =: (

щ

(1 ' е

ПРИ

(1.57)

^ 0 при /< 0 ,

 

где

 

о

 

 

 

 

 

 

 

Г ( а ) = \

ta~'e~adt.

(1.58)

 

 

0

 

 

Функция распределения

 

 

 

 

/

r - ' e - udt при

 

F(t) =

 

0;

0

 

 

(1.59)

0

при

/< 0 .

 

 

Вид кривых свидетельствует об универсальности закона гаммараспределения. Так, при значениях параметра формы а < 1 плот­ ности распределения наработки до отказа имеют убывающий ха­ рактер, при а = 1 и X= const распределение превращается в экспо­ ненциальное, при а > 3 кривая распределения гамма-распределе­ ния приближается к нормальному.

Математическое ожидание и дисперсия соответственно равны

(1.60)

При а < 1 интенсивность отказов монотонно убывает, что соответствует ситуации, когда происходит быстрое выгорание ненадежных элементов в период приработки изделия. При а > 1 интенсивность отказов возрастает, что характеризует период из­ нашивания и старения элементов.

При а = 1 гамма-распределение совпадает с экспоненциаль­ ным распределением; при а > 1 0 гамма-распределение прибли­ жается к нормальному закону. Если а принимает значения произвольных целых положительных чисел, то такое гамма-рас­ пределение называют распределением Эрланга. Если значение а

кратно 1/2, а А, = 1/2, то гамма-распределение совпадает с рас­ пределением хи-квадрат (х2) .

26

S)

Рис. 1.7. Кривые гамма-распределения при различных значениях параметров а и А.:

а — плотность распределения; б — вероятность безотказной работы

Кривые гамма-распределения для различных значений а и X приведены на рис. 1.7.

1.13.РАСПРЕДЕЛЕНИЕ ПУАССОНА

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

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

27

Плотность вероятности дискретного распределения для цело­

численных значений

т О, 1, 2,

имеет вид

 

 

 

 

 

 

 

 

^ ( 0

= 4 г е

 

 

 

 

 

 

(1.61)

где t — фиксированный интервал времени; А,>0.

 

 

Функция распределения описывается уравнением

 

 

 

 

 

оо

 

 

I

"тГе_^

"РИ ' > 0 ;

 

 

F{t)=

£

 

Pm(t) =

(1.62)

 

 

 

 

 

 

 

т < t

 

 

 

 

 

 

 

 

 

 

m= 0

 

 

О при

/г^О.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сумма всех вероятностей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

Л, <0 = 1.

 

 

 

 

 

(1.63)

 

 

 

 

 

 

т = 0

 

 

 

 

 

 

 

 

 

Чем меньше выбрано значение параметра X, тем асиммет­

ричнее распределение (рис. 1.8).

 

 

 

 

 

 

 

 

 

 

 

Рщ

0,606

0,303

0,075

0,0125

0,0016

0,0002

 

 

 

 

 

т

 

0

1

г

3

 

 

0

 

5

 

 

0,3

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О

 

_1

 

1

1

,

 

 

 

 

 

 

 

 

 

1

2

 

3

О

5

т

 

 

 

 

 

 

 

 

 

0,0067 0,033

0,08

0,10

0,175

0,175

0,10

0,1

0,06

0,036

0,02

 

7 7 1

О

 

1

2

3

0

5

 

6

7

 

8

9

10

Рис. 1.8. Кривые распределения Пуассона для различных значений параметра к:

а — при Л.= 0,5; б — при А.= 5,0

28

Математическое ожидание и дисперсия распределения Пуас­ сона соответственно равны

M(t) = k9 D (t) = X.

(1.64)

1.14.БИНОМИАЛЬНОЕ РАСПРЕДЕЛЕНИЕ

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

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

Плотность вероятности описывается выражением

р- =

^ г Ь ) Г ^

 

т = ° ’]....... “•

с -65»

Дискретная функция распределения имеет вид

 

F(t)= X

О

при

/^ 0 ;

 

 

I

 

П\

Рт( 1 - Р )п- т при 0 < / < л ;

пг\ (п —т)!

т < t

т < t

t > n .

 

 

 

1 при

 

 

 

 

 

 

 

( 1.66)

Суммирование в формуле (1.66) распространяется на все целочисленные значения m = 0, 1, 2, я, которые меньше /, причем

X Рт= 1-

(1.67)

т = 0

 

Математическое ожидание и дисперсия соответственно равны

M(f) = nP9 D(f) = nP(l — P).

(1.68)

29

Биномиальное распределение является однопараметричес­ ким и при приближении параметра Р к значению 0,5 распре­ деление становится все более симметричным.

Пример 1.6. Производится четыре независимых выстрела в одинаковых условиях, причем вероятность попадания при вы­ стреле равна /> = 0,5.

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

Рис.. 1.9. Кривая плотности распреде­

ления

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

к примеру 1.6

Ре ш е н и е . Для

определения

вероятностей воспользуемся

формулой (1.65):

 

 

 

Р0 = (1 - Р)4 = (1 _

0,5)4 = 0,0625;

Pl = W

0’5 (l~ ° ’5)3==0’25;

Рг =

W

° ’5 2 (I_ 0 ’5)2=0’375;

Р3 = |j-0,53 (1-0,5) =0,25; Р4 = 0,54 (1 — 0,5)° =0,0625.

Плотность распределения вероятностей показана на рис. 1.9. Наиболее полное представление о различных законах рас­

пределения случайных величин дано в работах [4, 10].