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

книги из ГПНТБ / Эксплуатационная надежность сельскохозяйственных машин

..pdf
Скачиваний:
14
Добавлен:
22.10.2023
Размер:
10.69 Mб
Скачать

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

Близость законов определяется величиной Хм (кри­ терий Колмогорова), равной максимальной, по абсолют­ ной величине, разности

Dm — max | Psuaif) ^теор(01 >

(3.39)

умноженной на ] / N, где N — число изделий (или общий объем выборки).

Хм — DmV N .

(3.40)

Если Ля^1,0, то говорят о хорошем согласии между эмпирическим и теоретическим законом распределения.

Более точно оценка выполняется по функции вероят­ ности Р(ХN).

При P(XN) >0,3—0,4 —согласие хорошее.

До сих пор рассматривалась полная выборка. Часто при испытаниях приходится встречаться с так называе­ мыми усеченными выборками. Усеченную выборку полу­ чают например в такой ситуации: поставлены на испы­ тания N неремонтируемых изделий; из них за время ia отказало т изделий при наработках tu t% ...tm; для остальных N—т объектов наработки неизвестны, извест­ но лишь, что все они больше чем t а. Усеченные выборки могут быть и при оценке ремонтируемых изделий.

■ Если испытывается No изделий в течение времени ta и отказало из них то, при этом общее число отказов было По, то рассматривая время между отказами как время безотказной работы, можно условно считать, что испытывалось n0+Noт0 элементов, из которых N0т0

не отказало, и их наработка не менее чем / ta• То есть, как видим, задача аналогична предыдущей. Для усечен­ ной выборки нельзя определить параметры распределе­ ния Гср и ot по формулам (3.3) и (3.5). В этом случае поступают следующим образом. Вычисляют вероятности

отказов на отрезке времени ta по формулам

(рассматри­

ваем для конкретности неремонтируемые изделия)

т =

(3.41)

где ka — число интервалов At на отрезке времени 0—t а-

Затем находят преобразование, которое переводит вы­ ражение предполагаемого теоретического распределения в уравнение прямой. С помощью этого преобразования пересчитывают эмпирические значения F(ti) в новые координаты. Если пересчитанные точки F(tг) хорошо укладываются на прямую, то следовательно вид теоре­ тического распределения принят правильно. После этого определяют параметры найденной прямой. Их два: угол наклона и отрезок на оси ординат, отсекаемый прямой. По значениям этих параметров определяют два парамет­ ра теоретического распределения.

Покажем метод на примере трех законов, приведен­ ных в табл. 3.6.

Экспоненциальное распределение

F(t) = 1

 

Преобразовывая, получим:

 

1

= Xt

1п-

1 - F (0

 

или

 

у= k-t.

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

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

Для удобства расчетов используют так называемую вероятностную бумагу (рис. 8). В этом случае пара­ метр X определяется по формуле

 

X =

0,023

tg ф,

(3.42)

где

 

ширина графика, мм;

Rt = 1

'^MIlIll

■размах выборки;

 

ф

■угол

наклона,

определяемый не­

 

 

посредственно по графику.

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

52

tl

^2

*3

*4

h

 

h

*8

^9

ho

t n

5

15

25

35

45

55

65

75

85

95

105

F i

f 2

F 3

f 4

f 5

F 6

f 7

f 8

F e

Fio

F n

0 ,0 9

0 ,2 8

0 ,3 7

0 ,5

| 0 ,6

0 ,6 8 0 ,7 4 0 ,8

0 ,8 4 0 ,8 7

0 ,9

 

 

 

 

1

 

 

 

 

 

 

Рис. 8. Вероятностная бумага для экшшенциального распре­ деления.

Экспоненциальный закон:

Хе~'4\ / > 0

№ =

0; f < 0

 

к — параметр распределения. Из графика:

к = 0,023 — tgi|>; L — 100 мм;

Rt

Rt = Rmn ~ Rmn = 100 ч; ф = 42°; tg* = 0,9004;

к = 0,023-1-0,9004» 0,02 1/ч;

RT — размах вариации, ч L — ширина графика, мм.

Нормальное распределение

о о

53

Вводят обозначение аргумента

t Т

U*

о

где — квантиль нормального распределения, то есть это такое значение аргумента, которое соответствует ве­ роятности отказа — а.

Тогда

U. =

t Г

о а

и уравнение прямой имеет вид

 

1

у = kt а,

где k =

Т

----;

а = -— ; y = U« .

 

а

о

Вероятностная бумага для нормального распределе­ ния приведена на рис. 9. При у = 0 имеем

(3.43)

k

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

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

о = 48,5 ■— ctgip.

(3-44)

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

tl

*2

 

U

h

*6

u

*8

*9

* 1 0

‘u

42

fl3

0,5

1,5

2,5

3 ,5

4 , 5

5 ,5

j 6 ,5

7,5

8,5

9,5

10,5

11,5

12,5

F i

F 2

F 3

f 4 F 6 F 6

f 7 f 8

f 9

F 10 F u

F M

F13

 

 

 

 

 

1

 

 

 

 

 

 

 

0,001 |о,002jo,005|0,Ог|о,09jo,29jo, 57 0,8ljo,955 0,9sjo,99^0,9 9 4 J0 ,997

54

Рис. 9. Вероятностная бумага для нормального распределения.

Нормальный закон:

(t- ту

т =

2J2

а ]/ 2я

 

Я* = Ямам — #шш — размах вариации;

L — ширина графика, мм;

а, т — параметры распределения. Из графика:

55

at = 48,5 ctg o|r, Rt = 12-103 ч; L = 120 мм;

Ф = 78°; ctgi|> = 0,2126; a = 48,5• 100-0,2126 = 1030 ч;

T = 6300 4.

Распределение Вейбулла-Гнеденко

 

F(*) = 1 — e

[io 1 .

 

Преобразуем к виду

 

 

 

 

b lg* = lg*£ +

lg In ------------- .

 

0

 

1

-

F(t)

Уравнение прямой будет

 

 

 

 

где

у = kx а,

 

 

 

 

 

 

 

У=

lg In

1

 

х = lg*;

Щ

 

1 -

 

 

 

а = lg tb0\

k =

b.

 

Тангенс угла

наклона

графика

определяет пара­

метр Ь. Абсцисса точки у = 0 на графике определяет вто­ рой параметр *о. Тогда

а

= lg*o-

х= lg * = - у

Для этого распределения также имеется вероятност­ ная бумага (рис. 10). Параметры *0 ;и Ъ с помощью ве­ роятностной бумаги определяются по выражениям

lg*0 = St -у -;

Ь = 0 , 0 1 3 ( 3 . 4 5 )

L

Rt

Здесь St — абсцисса точки пересечения прямой с осью lg*, определенная непосредственно по графику.

Пример 3. Дан ряд распределения случайной -вели­ чины, характеризующий наработку до предельного со­ стояния подшипника КПП трактора. Проверить соот­ ветствует ли распределению Вейбулла — Гнеденко и определить параметры распределения.

56

h

*2

*3

u

*5

*6

u

*8

*9

*10

tu

*12

*13

0 , 5

1 ,5

2 , 5 | 3 , 5

4 ,5

5 , 5

6 ,5

7 ,5

8 , 5

9 ,5

10,5

11,5

12,5

Fi

f2

F3

f4

F6

F6

f7

f8

f9

F10

Fu

^12

Fj.3

 

 

 

 

 

 

 

 

 

 

 

 

0 ,0 0 8

0 ,0 2

0 ,0 5

0 ,1 4

0 ,2 6 0 ,4 0 0 ,5 4

0 ,6 7 0 ,7 9

0 ,8 6

0 ,9 4

0 ,9 7

0 ,9 9 2

Закон Вейбулла — Гнеденко:

to, b — параметры распределения.

S,

L_ igt.

 

Rt

Шкала t построена по уравнению

S t = 100 lgf мм; 1 < * < 10.

Если 10</^100, значение t умножить на 10.

Если

100 5^<1000, значение t

умножить на 100 и

т. д.

 

 

 

 

Из графика:

 

 

*0

=

7,1 • 103 = 7100 ч; Ъ= 0,013 —

- tg гф;

 

 

 

Ro

 

6

=

0,013-100tg-ф; ф = 62°;

tg лр =

1,881;

 

 

Ь = 0,013-100-1,881 =2,46.

 

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

тическому.

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

57

F

Рис. 10. Вероятностная бумага для закона Вейбулла — Гнеденко.

58

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

2*< + (Лl - m ) t a

;=1

(3.46)

 

Ремонтируемые изделия:

 

 

2 */ + М> -

m0)ta

Т,ср. о

i—\

(3.47)

пгл

 

 

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

Пример 1. 17 тракторов проходили испытания в тече­ ние 2000 часов. На трех из них имели место поломки ры­ чагов муфт сцепления в моменты времени ^ = 100 ч, ^2=700 ч, = 1.200 ч, после чего эти машины были сняты с испытаний.

Определить среднюю наработку на отказ и 80-ный ресурс.

Закон распределения экспоненциальный (внезапные отказы).

Средняя наработка на отказ

_

100 + 700+ 1200 -H 4 . 2000 =

ср.о

3

 

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

 

 

 

Я ( Г ) = е

г ср = 0 ,8 .

Отсюда находим tf =2250 ч.

Таким образом, чтобы

обеспечить 80-ный ресурс

2250 часов,

необходимо среднее время безотказной ра­

боты иметь значительно больше — 10 000 часов.

Пример 2. Для исходных данных примера 1 опреде­ лить ТСР' 0 и 7, при условии, что муфты были отремонтиро­ ваны и отработали полный срок испытаний — 2000 ч.

Теперь суммарное время испытаний равно

59

fcyM= 100 + 1900 + 700 + 1300 + 1200 -f 800 +

+ 14 ■2000 = 34000 ч;

T

=' 3400 11300 ч;

1 cp. о

 

tf — определяется из уравнения е 11300 = 0,8.

Имеем tf. =2540 ч.

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

290 часов в связи с тем,

что после ремонта отказов

не было.

Полученные выше точечные

Интервальная оценка.

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

Рис. И . Закон распреде­ ления точечной оценки.

Оценка точности выборочных характеристик произво­ дится е помощью доверительных интервалов и довери­ тельных вероятностей.

В общем случае задача о доверительных вероятно­ стях и доверительных интервалах решается следующим образом. Из теоретических соображений или из опыта определяется закон распределения точечной оценки Тср (рис. 11). Отсекая от нижней и верхней границы значе-

а

 

вероятность нахождения

оценки

ния — , получим, что

2

 

 

 

 

Тср в интервале (7\,Ср, Т2.ср) равна

 

F(Ti, ср <

7ср <

Т2,ср )=

F(T2jСр) - F(TUср) =

 

=

! -----1------ 1

~ = 1— « = р.

(3.48)

60

Соседние файлы в папке книги из ГПНТБ