
книги из ГПНТБ / Эксплуатационная надежность сельскохозяйственных машин
..pdfнов, тем самым проверяют принятый сперва без доказа тельства теоретический закон распределения.
Близость законов определяется величиной Хм (кри терий Колмогорова), равной максимальной, по абсолют ной величине, разности
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а — квантиль нормального распределения, то есть это такое значение аргумента, которое соответствует ве роятности отказа — а.
Тогда
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