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

Основы теории надежности. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

проверить гипотезу о виде распределения;

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

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

Решение

Построение интервального вариационного ряда. Оптимальный интервал h вычисляется по формуле Стэрджеса [4]:

h (lmax lmin )(1 3,322lg N ) ,

где lmax = 312,5 тыс. км, lmin = 66,3 тыс. км – соответственно максимальное и минимальное значение ресурса (табл. 2.2); N = 66 – общее число наблюдений.

Тогда

h = (312,5 – 66,3) / (1 + 3,322 lg N) = 34,9 тыс. км.

Окончательно принимаем ближайшее целое число h = 35 тыс. км. За начало первого интервала принимаем величину

l1 lmin h2 = 66,3 – 35/2 = 48,8 тыс. км;

начало второго интервала совпадает с концом первого интервала и

равно

 

 

 

 

 

l2 l1 h = 48,8 + 35 = 83,8 тыс. км;

 

далее

 

 

 

 

l3

= 83,8 + 35 = 118,8 тыс. км;

l4

= 118,8 + 35 = 153,8

тыс. км;

l5

= 153,8 + 35 = 188,8 тыс. км;

l6

= 188,8 + 35 = 223,8

тыс. км;

l7

= 223,8 + 35 = 258,8 тыс. км;

l8

= 258,8 + 35 = 293,8

тыс. км;

l9

= 293,8 + 35 = 328,8 тыс. км.

 

 

 

Шкала интервалов и группировка результатов наблюдений приведены в табл. 2.2.

51

 

 

 

Таблица 2.2

Шкала интервалов и частота попадания в интервал

 

 

 

Частота попада-

Номер

Границы интерва-

Середины интерва-

интервала

лов (тыс. км)

лов lr, (тыс. км)

ния в интервал ni

1

48,8…83,8

66,3

1

2

83,8…118,8

101,3

4

3

118,8…153,8

136,3

6

4

153,8…188,8

171,3

21

5

188,8…223,8

206,3

20

6

223,8…258,8

241,3

12

7

258,8…293,8

276,3

1

8

293,8…328,8

311,3

1

Определение параметров и характеристик нормального распре-

деления. Плотность вероятности имеет вид

 

 

 

 

(l m )2

 

 

 

1

 

i

l

 

 

 

 

2 2

 

f (l)

 

 

e

 

 

.

 

2

 

 

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

математическое ожидание (для интервального вариационного ряда, табл. 2.2)

 

 

 

 

r

 

r

n 66,3

1 101,3 4 136,3 6 171,3 21

m

l

 

l

n

 

l

 

 

 

r

i

i

1 4 6 21

 

 

 

i 1

 

i 1

206,3 20 241,3 12

276,3 1 311,3 1 12495,8 189,3 тыс. км;

 

 

 

 

20 12

1 1

66

среднеквадратическое отклонение (для интервального вариационного ряда)

r

(lr ml )2 ni

 

i 1

 

44,8 тыс. км;

 

N 1

 

 

 

52

значения эмпирической плотности распределения вероятностей fэ(lr ) по интервалам наработки (табл. 2.3):

 

 

f

э

(l )

 

ni

;

 

 

 

 

 

 

 

 

 

 

 

r

 

Nh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.3

 

Расчетные параметры распределения

 

 

 

 

 

 

 

 

 

 

 

ni

yr

fэ(lr )

 

 

 

f0 ( yr )

fт(lr )

n1

–2,746

0,0004

 

 

 

0,0092

0,0002

 

n2

–1,964

0,0017

 

 

 

0,0579

0,0013

 

n3

–1,183

0,0026

 

 

 

0,1982

0,0044

 

n4

–0,402

0,0091

 

 

 

0,3680

0,0082

 

n5

0,379

0,0087

 

 

 

0,3710

0,0083

 

n6

1,161

0,0052

 

 

 

0,2036

0,0045

 

n7

1,942

0,0004

 

 

 

0,0603

0,0013

 

n8

2,723

0,0004

 

 

 

0,0099

0,0002

 

нормированные и центрированные отклонения середины интервалов yr (см. табл. 2.3):

 

 

 

 

 

yr (lr ml ) ;

 

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

fт(lr )

(см. табл. 2.3):

 

 

 

 

 

 

fт(lr ) (1/ ) f0 ( yr ) ,

где

f0

( yr )

1

e yr2 2

– плотность вероятности нормального распре-

2

деления (см. табл. П2 приложения);по результатам расчета строятся гистограмма середины интерва-

лов fn (l), кривая распределения плотностей вероятностей fэ(lr ) и теоретическая кривая распределения fт(lr ) (рис. 2.5).

Проверка согласия между эмпирическим и теоретическим (нормальным) распределениями по критерию Пирсона. Как следует из анализа интервального вариационного ряда (см. табл. 2.2), для применения критерия Пирсона необходимо объединить интервалы 1 и 2, интервалы 7 и 8 с интервалом 6, так как в интервалах 1, 2, 7, 8 количество

53

наблюдений менее 5. В результате объединения получим следующий

ряд распределения:

 

 

 

 

 

 

 

 

 

 

 

Интервалы

48,8…

118,8…

 

 

153,8…

 

188,8…

223,8…

mi = ni

118,8

 

153,8

 

 

 

188,8

 

223,8

328,8

5

 

 

 

6

 

 

 

21

 

20

14

f(l)

 

 

n

 

 

 

 

 

 

 

 

 

 

0,010

 

25

 

 

 

 

 

fэ(l)

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

15

 

 

 

 

 

 

 

fn(l)

 

0,005

 

 

 

 

fт(l)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

101,3

171,3

241,3

l,

 

 

 

 

 

 

 

 

 

 

311,3 тыс. км

 

 

 

66,3

136,3

 

 

206,3

276,3

 

 

 

 

 

 

 

Рис. 2.5

 

 

 

 

Определим вероятность

pi (i 1,5) ;

р1

выражает вероятность того,

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

принимает значение, принадлежащее интервалу (48,8…118,8), т. е.

p P(48,8 X

118,8) P

48,8 ml

X ml

118,8 ml

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

118,8 m

 

1

 

 

48,8 m

 

 

 

2

Ф

 

l

 

2

Ф

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

54

 

 

 

 

 

 

 

 

1

118,8 189,3

 

 

1

 

48,8 189,3

 

 

2

Ф

44,8

 

2

Ф

44,8

 

 

 

 

 

 

 

 

12 Ф( 1,57) 12 Ф( 3,14)

12 Ф(3,14) 12 Ф(1,57) 0,4992 0,4418 0,0574,

где Ф(z) определяется по табл. П4 приложения. Аналогично получаем:

 

 

p2 0,1565 ;

p3 0,2813 ;

 

p4 0,2834 ;

p5

0, 2197 .

Для нахождения статистики 2

составим табл. 2.4.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.4

 

 

 

 

Расчет статистики 2

 

 

 

 

 

 

 

Интервал после

mi

 

рi

 

i

 

 

mi npi

 

 

(mi npi )2

 

 

п/п

 

объединения

 

 

 

 

 

 

1,21

 

 

npi

1

 

48,8…

118,8

5

 

0,0574

 

3,79

 

 

 

0,385

2

 

118,8…153,8

6

 

0,1565

 

10,33

 

4,33

 

 

1,815

3

 

153,8…188,8

21

 

0,2813

 

18,57

 

2,43

 

 

0,318

4

 

188,8…223,8

20

 

0,2834

 

18,70

 

1,30

 

 

0,090

5

 

223,8…328,8

14

 

0,2197

 

14,50

 

0,50

 

 

0,017

 

 

 

66

 

0,9983

 

65,89

 

 

 

 

 

2,625 = χ2

В

результате имеем:

количество интервалов

r

= 5, количество

параметров распределения (ml , ) s = 2, число степеней свободы k = r – s – 1 = 5 – 2 – 1 = 2. В табл. П3 приложения для выбранного

уровня значимости 0,01 и k = 2 соответствует 22;0,01 9,2 .

В нашем случае 2,625 < 9,2, следовательно, нет оснований отвергнуть гипотезу о нормальности распределения ресурса.

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

ность безотказной работы детали по интервалам наработки (см. табл. 2.2) по формуле

55

 

1

r

 

r

 

 

P(l ) 1

n

N

 

n

N ;

 

 

i

 

i

 

i

 

 

N i 1

 

i 1

 

 

P(l1) = (66 – 1)/66 = 0,98;

 

 

P(l2) = (66

– 5)/66 =

0,92;

P(l3) = (66 – 11)/66 = 0,83;

 

 

P(l4) = (66

– 32)/66 = 0,52;

P(l5) = (66 – 52)/66 = 0,21;

 

 

P(l6) = (66

– 64)/66 = 0,03;

P(l7) = (66 – 65)/66 = 0,015;

 

 

P(l8) = (66

– 66)/66 = 0.

P(l)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,0

66,3

101,3 136,3 171,3 206,3 241,3 276,3 311,3

l

 

тыс. км

Рис. 2.6.

На рис. 2.6 построена кривая вероятности безотказной работы детали P(l) в зависимости от наработки l.

56

3.ПОТОКИ ОТКАЗОВ И ВОССТАНОВЛЕНИЙ

ВТЕОРИИ НАДЕЖНОСТИ

3.1.Характеристики потоков отказов

ивосстановлений

Влюбой момент времени t объект может находиться в одном из двух состояний: отказа и работоспособности. Обозначим:

Z (t) – неизвестное текущее состояние объекта в момент времени t; Z*– состояние отказа; Z – состояние работоспособности.

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

ξ1, η1, ξ2, η2, … ξi, ηi (рис. 3.1), где ξi, ηi – длительность i-го по счету периода работоспособности и периода отказа соответственно. Поэтому в теории надежности принято пользоваться понятием потока событий.

 

ξ1

 

η1

 

ξ2

η2

ξ3

 

η3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

t

Рис. 3.1

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

Основными потоками событий в теории надежности являются потоки отказов и восстановлений. По своему характеру они могут быть различными. Наибольшее применение в теории надежности получили простейший поток и потоки Эрланга [5]. Наиболее важными характеристиками потоков отказов и восстановлений являются математическое ожидание числа отказов (восстановлений) Ω(t) на интервале (0, t), интенсивность λ(t) или μ(t) и ω(t) – параметр потока отказов.

Под интенсивностью потока отказов (восстановлений) понимает-

ся математическое ожидание числа отказов (восстановлений) в единицу времени.

Простейший поток. Простейшим потоком событий называется поток, обладающий свойствами стационарности, ординарности и отсутствия последействия. Название «простейший» обусловлено тем, что

57

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

Стационарность потока отказов – вероятность появления опреде-

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

Ординарность – вероятность возникновения на элементарном отрезке времени t двух или более отказов пренебрежимо мала по сравнению с вероятностью возникновения одного отказа. Применительно к описанию надёжности ординарность означает, что одновременное возникновение двух разных отказов, например у автомобиля, практически маловероятно.

Отсутствие последействия – вероятность появления фиксирован-

ного числа отказов на интервале времени (t, t + τ) не зависит от того,

сколько отказов появилось до момента τ. Практически отсутствие по-

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

ляются в те или другие моменты вре-

Pk

 

 

 

 

мени независимо друг от друга.

a = 1

 

 

 

На практике суммирование не менее

 

 

 

0,3

 

 

 

 

6…8 элементарных потоков приводит к

a = 3

 

 

 

образованию простейшего или близко-

 

 

 

 

0,2

 

 

 

 

го к нему потока. Для простейшего по-

 

 

a = 6

 

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

 

 

 

 

0,1

 

 

 

 

ния определенного числа отказов в те-

 

 

 

 

чение времени

определяется законом

0

 

 

 

 

Пуассона (рис. 3.2):

4

6

8

k

P (t) ( t)k e t ,

2

 

Рис. 3.2

 

 

k

k!

 

 

 

 

 

 

 

 

 

 

 

где k = 0, 1, 2, – число отказов, возникающих за время t;

– пара-

метр потока отказов.

 

 

 

 

 

58

В реальных условиях производства обычно фиксируют значение t,

например, 1 ч, 1 смена, 1 неделя и так далее, т. е. t = 1, а t 0 a – среднее число отказов, возникающих за время t. В этом случае

Pk ak e a . k!

Так, поступление автомобиля в зону ремонта для устранения неисправности принято называть требованием. В реальных условиях требование может включать комбинацию неисправностей агрегатов и автомобилей. Используя последнюю формулу, можно установить вероятность появления определенного числа требований Pk при известном значении а. Например, при значении а = 3 вероятность отсутствия

требований Pk 0 30 e 3 0,05 , или 5 %; вероятность появления одно- 0!

го требования 0,15; двух 0,22; трех тоже 0,22; четырех 0,16 и т. д. (рис. 3.2). Таким образом, загрузка постов и оборудования носит вероятностный характер: 22 % от всех смен будет иметь фактическое число требований, совпадающее со средним, у 42 % – (5 + 15 + 22) загрузка будет меньше, а в 36 % имеем (100 – 22 – 42) случаев больше средней.

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

неполной загрузке зон и участков;

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

В зависимости от стоимости простоя автомобилей в ожидании ре-

монта (Са), а также оборудования и рабочих в ожидании автомобилей

(Со.р), требующих ремонта, определяют оптимальную пропускную

способность зон, участков, постов ТО и ремонта. Эта задача решается с использованием теории массового обслуживания [1, 8], и из условия минимизации выражения, называемого целевой функцией,

u = Са + Со.р.

Характерным признаком закона Пуассона является равенство дисперсии среднему значению, поэтому коэффициент вариации потока

требований υ а 0,5 . Это означает, что с увеличением программы работы вариация ее фактического значения сокращается.

59

Средняя программа Коэффициент вариации

1…………………………………1

2…………………………………0,71

3…………………………………0,58

4…………………………………0,50

5…………………………………0,45

9…………………………………0,30

25…………………………………0,20

Закон распределения становится более симметричным с увеличением программы (рис. 3.2, при а = 6), что благоприятно сказывается на организации технологического процесса технического обслуживания (ТО) и ремонта. Поэтому укрупнение предприятий, централизация и кооперирование ТО и ремонта, приводящие к увеличению программы работы, являются одним из направлений совершенствования технической эксплуатации автомобилей.

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

x 1/ , x , υ 1 .

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

следующим образом:

t

a (t)dt .

t

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

Потоки Эрланга. Нарушение свойств стационарности или наличие последействия приводит к тому, что поток становится непростейшим. К непростейшим потокам относятся потоки Эрланга, которые возникают «просеиванием» простейшего потока отказов.

60

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