1. Определение коэффициентов годности и восстановления деталей
1.1 Определение технических требований к анализируемой поверхности
Проведём выкопировку эскиза указанной детали и сформируем технические требования на дефектацию заданной поверхности 6 см. [3].
Таблица 1 - Технические требования на дефектацию
Наименование
детали
|
Контролируемая
поверхность
|
Размер детали
|
Корпус коробки передач трактора
МТЗ-82
|
Поверхность
отверстия под стакан ведущей шестерни 2-й ступени редуктора
|
по
чертежу
|
допустимый в сопряжении
|
138 +0,040
|
с деталями бывшими в эксплуатации
|
с новыми деталями
|
138,07
|
138,09
|
Эскиз указанной детали приведен в приложении А.
1.2 Определение износов деталей и составление вариационного ряда
Значения размеров изношенных деталей (для отверстия – по возрастанию значений размеров) приведены в таблице 2.
Таблица 2 – Размеры изношенных деталей, мм
138,062
|
138,073
|
138,076
|
138,080
|
138,084
|
138,089
|
138,094
|
138,101
|
138,109
|
138,114
|
138,062
|
138,073
|
138,078
|
138,081
|
138,085
|
138,089
|
138,094
|
138,101
|
138,109
|
138,116
|
138,064
|
138,073
|
138,078
|
138,081
|
138,085
|
138,090
|
138,094
|
138,102
|
138,110
|
138,116
|
138,066
|
138,073
|
138,079
|
138,082
|
138,086
|
138,090
|
138,097
|
138,103
|
138,110
|
138,118
|
138,068
|
138,074
|
138,079
|
138,082
|
138,086
|
138,091
|
138,097
|
138,104
|
138,110
|
138,118
|
138,069
|
138,074
|
138,079
|
138,082
|
138,087
|
138,091
|
138,098
|
138,104
|
138,110
|
138,121
|
138,070
|
138,075
|
138,079
|
138,082
|
138,087
|
138,091
|
138,099
|
138,105
|
138,110
|
138,122
|
138,071
|
138,075
|
138,079
|
138,083
|
138,088
|
138,092
|
138,099
|
138,106
|
138,111
|
138,126
|
138,073
|
138,075
|
138,079
|
138,083
|
138,088
|
138,092
|
138,100
|
138,107
|
138,113
|
138,126
|
138,073
|
138,076
|
138,080
|
138,083
|
138,089
|
138,093
|
138,100
|
138,107
|
138,113
|
138,126
|
Вычислим износы деталей и составим их вариационный ряд в виде таблицы 3.
Износ i
-го отверстия определяют по зависимости
; (1)
где –диаметр i-го изношенного отверстия;
– наибольший конструктивный размер отверстия;
N
– число анализируемых деталей.
Пример расчета: износ 1-го отверстия:
мм.
Таблица 3 – Значения износов деталей (вариационный ряд)
Номер детали
|
Значение износа детали, мм
|
Номер детали
|
Значение износа детали, мм
|
Номер
детали
|
Значение износа детали, мм
|
Номер детали
|
Значение износа детали, мм
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
1
|
0,022
|
26
|
0,039
|
51
|
0,049
|
76
|
0,064
|
2
|
0,022
|
27
|
0,039
|
52
|
0,049
|
77
|
0,065
|
3
|
0,024
|
28
|
0,039
|
53
|
0,050
|
78
|
0,066
|
4
|
0,026
|
29
|
0,039
|
54
|
0,050
|
79
|
0,067
|
5
|
0,028
|
30
|
0,040
|
55
|
0,051
|
80
|
0,067
|
6
|
0,029
|
31
|
0,040
|
56
|
0,051
|
81
|
0,069
|
7
|
0,030
|
32
|
0,041
|
57
|
0,051
|
82
|
0,069
|
8
|
0,031
|
33
|
0,041
|
58
|
0,052
|
83
|
0,070
|
9
|
0,033
|
34
|
0,042
|
59
|
0,052
|
84
|
0,070
|
10
|
0,033
|
35
|
0,042
|
60
|
0,053
|
85
|
0,070
|
11
|
0,033
|
36
|
0,042
|
61
|
0,054
|
86
|
0,070
|
12
|
0,033
|
37
|
0,042
|
62
|
0,054
|
87
|
0,070
|
13
|
0,033
|
38
|
0,043
|
63
|
0,054
|
88
|
0,071
|
14
|
0,033
|
39
|
0,043
|
64
|
0,057
|
89
|
0,073
|
15
|
0,034
|
40
|
0,043
|
65
|
0,057
|
90
|
0,073
|
16
|
0,034
|
41
|
0,044
|
66
|
0,058
|
91
|
0,074
|
17
|
0,035
|
42
|
0,045
|
67
|
0,059
|
92
|
0,076
|
18
|
0,035
|
43
|
0,045
|
68
|
0,059
|
93
|
0,076
|
19
|
0,035
|
44
|
0,046
|
69
|
0,060
|
94
|
0,078
|
20
|
0,036
|
45
|
0,046
|
70
|
0,060
|
95
|
0,078
|
21
|
0,036
|
46
|
0,047
|
71
|
0,061
|
96
|
0,081
|
22
|
0,038
|
47
|
0,047
|
72
|
0,061
|
97
|
0,082
|
23
|
0,038
|
48
|
0,048
|
73
|
0,062
|
98
|
0,086
|
24
|
0,039
|
49
|
0,048
|
74
|
0,063
|
99
|
0,086
|
25
|
0,039
|
50
|
0,049
|
75
|
0,064
|
100
|
0,086
|
1.3 Составление статистического ряда износов
Число интервалов n
определяют по зависимости:
(2)
с последующим округлением полученного результата до целого числа
=.
Длину интервалов вычисляют по зависимости:
, (3)
где и – наибольшее и наименьшее значения СВ из вариационного ряда соответственно.
мм.
Начало t
нi
и конец t
кi
i
-го интервала вычисляют по следующим зависимостям:
t
н
1
= t
min
; t
н
i
= t
к
(i
–1)
; t
к
i
= t
н
i
+ h
(4)
Пример решения:
t
н1
= t
min
=0,022 мм;
t
к1
= t
н1
+ h
=0,022+0,0064=0,0284 мм.
Количество наблюдений (значений СВ) в i
-м интервале (i
= 1, …, n
) называется опытной частотой
. Опытная частота , отнесенная к общему числу наблюдений (объему выборки) , называется опытной вероятностью.
.
Ее значение определяется по зависимости:
, (5)
где – значение СВ в середине i
-го интервала.
Пример решения:
.
Накопленная опытная вероятность
, являющаяся статистическим аналогом функции распределения, вычисляется по зависимости:
(6)
Пример решения:
.
Таким образом, статистическим рядом распределения является таблица 4, в которой указаны границы и середины интервалов, опытные частоты, опытные и накопленные опытные вероятности.
Таблица 4 – Статистический ряд распределения износов
Границы
интервала,
мм
|
0,0220
...
0,0284
|
0,0284
...
0,0348
|
0,0348
...
0,0412
|
0,0412
...
0,0476
|
0,0476
...
0,0540
|
0,0540
...
0,0604
|
0,0604
...
0,0668
|
0,0668
...
0,0732
|
0,0732
...
0,0796
|
0,0796
…
0,0860
|
Середина интервала,
мм
|
0,025
|
0,031
|
0,038
|
0,044
|
0,050
|
0,057
|
0,063
|
0,070
|
0,076
|
0,082
|
Опытная частота
|
5
|
11
|
17
|
14
|
15,5
|
7,5
|
8
|
12
|
5
|
5
|
Границы
интервала,
мм
|
0,0220
...
0,0284
|
0,0284
...
0,0348
|
0,0348
...
0,0412
|
0,0412
...
0,0476
|
0,0476
...
0,0540
|
0,0540
...
0,0604
|
0,0604
...
0,0668
|
0,0668
...
0,0732
|
0,0732
...
0,0796
|
0,0796
…
0,0860
|
Опытная вероятность
|
0,05
|
0,11
|
0,17
|
0,14
|
0,155
|
0,075
|
0,08
|
0,12
|
0,05
|
0,05
|
Накопленная опытная вероятность
|
0,05
|
0,16
|
0,33
|
0,47
|
0,625
|
0,7
|
0,78
|
0,9
|
0,95
|
1
|
1.4 Определение числовых характеристик статистической совокупности износов
Наиболее применяемыми числовыми характеристиками совокупности значений случайной величины являются:
– среднее значение, характеризующее центр группирования случайной величины;
– среднеквадратическое отклонение и коэффициент вариации, являющиеся характеристиками рассеивания случайной величины.
Так как > 25, то характеристики вычисляются по зависимостям:
, (7)
, (8)
Анализ зависимостей для определения показывает, что его значение зависит не только от величины рассеивания, но и от абсолютных значений СВ. От этого недостатка свободен коэффициент вариации , определяемый по зависимости:
(9)
где при N
> 25 t
см
= t
н1
–0,5h
;
t
см
= t
н
1
–0,5h
=0,022 - 0,5∙0,0064= 0,0188 мм.
1.5 Проверка однородности информации об износах
Проверку на выпадающие точки проводят по критерию Ирвина , который вычисляют по зависимости:
, (10)
где и – смежные значения случайной величины вариационного ряда.
Проверку начинают с крайних значений случайной величины. Вычисленное сравнивают с табличным значением , взятом из табл. В.1 [1], при доверительной вероятности и числе наблюдений .
При переходят к проверке однородности следующего значения СВ. При проверяемое значение СВ признают выпадающим (экстремальным), и оно исключается из выборочной совокупности наблюдений.
Пример решения:
.
при N=100, значение критерия Ирвина
Вычисленные значения критерия Ирвина запишем в таблицу 5.
Таблица 5 – Значения критерия Ирвина
-
|
0
|
0
|
0
|
0,063
|
0
|
0,063
|
0,063
|
0,126
|
0,063
|
0
|
0
|
0,126
|
0,063
|
0,063
|
0
|
0
|
0
|
0
|
0,126
|
0,126
|
0
|
0
|
0
|
0
|
0,063
|
0
|
0,063
|
0,063
|
0
|
0,126
|
0
|
0,063
|
0,063
|
0,063
|
0
|
0,189
|
0,063
|
0
|
0,126
|
0,126
|
0,063
|
0
|
0
|
0
|
0,063
|
0
|
0,063
|
0
|
0
|
0,063
|
0
|
0
|
0
|
0,063
|
0
|
0,063
|
0
|
0
|
0,189
|
0,063
|
0,063
|
0
|
0
|
0
|
0
|
0,063
|
0,063
|
0
|
0,063
|
0,063
|
0
|
0
|
0,063
|
0,063
|
0,063
|
0
|
0,063
|
0,063
|
0,253
|
0,126
|
0
|
0
|
0
|
0
|
0
|
0,063
|
0,063
|
0,126
|
0
|
0
|
0,063
|
0,063
|
0
|
0,063
|
0,063
|
0
|
0
|
0
|
0
|
Вычисленные значения сравним с табличным значением
Взятом из таблицы В.1 [1] при доверительной вероятности и числе наблюдений N=100
Отсюда следует, что все точки однородны.
1.6 Графическое построение опытного распределения износов
Для наглядного представления опытного распределения, оценки качества произведенного группирования (разделения на интервалы) и более обоснованного выдвижения гипотезы о предполагаемом теоретическом распределении по данным статистического ряда строим гистограмму, полигон и график накопленной опытной вероятности (приложения Б, В, Г).
1.7
Выравнивание опытной информации теоретическим законом распределения
1.7.1 Выдвижение гипотезы о предполагаемом теоретическом законе распределения
Вычисленное значение коэффициента вариации V=0,492
При значении коэффициента вариации V=0,30…0,50 возникает неопределённость. В этой ситуации гипотезы о НЗР и ЗРВ являются равноправными, поэтому производится расчёт дифференциального и интегрального законов распределения обоих видов с последующей проверкой правдоподобия каждого из них по одному из критериев согласия и принятием соответствующего решения.
1.7.2 Расчет и построение дифференциального и интегрального ТЗР
Для нормального закона распределения
Так как при составлении статистического ряда (см. таблицу 4) были вычислены не статистические плотности функции распределения , а опытные вероятности попадания наблюдений в -й интервал , то для обеспечения сравнимости распределений вычислим теоретические вероятности этих же событий по зависимости:
, (11)
где – длина интервала, принятая при построении статистического ряда;
– квантиль нормального распределения, значение которого вычислено для середины -го интервала ;
– значение центрированной и нормированной плотности распределения из приложения Г [1] (при этом следует учесть, что );
n
- число интервалов, принятое при составлении статистического ряда.
Пример решения для середины 1-го интервала:
Значения теоретических вероятностей запишем в таблицу 6.
Таблица 6 - Значения теоретических вероятностей
Середина интервала,
мм
|
0,025
|
0,031
|
0,038
|
0,044
|
0,050
|
0,057
|
0,063
|
0,070
|
0,076
|
0,082
|
Плотность функции распределения f(z)
|
0,11
|
0,19
|
0,29
|
0,37
|
0,4
|
0,37
|
0,29
|
0,19
|
0,11
|
0,05
|
Теоретическая
вероятность
|
0,044
|
0,076
|
0,117
|
0,149
|
0,162
|
0,149
|
0,117
|
0,076
|
0,044
|
0,02
|
Вычисление функции распределения осуществляется по зависимости:
; , (12)
где – квантиль нормального распределения, значение которого вычислено для конца -го интервала ;
– значение интегральной функции нормального распределения (при этом следует учесть, что ).
Вычислим функцию распределения на 1-м интервале:
.
Значения функции распределения запишем в таблицу 7.
Таблица 7 – Значения функции распределения
Границы
интервала,
мм
|
0,0220
...
0,0284
|
0,0284
...
0,0348
|
0,0348
...
0,0412
|
0,0412
...
0,0476
|
0,0476
...
0,0540
|
0,0540
...
0,0604
|
0,0604
...
0,0668
|
0,0668
...
0,0732
|
0,0732
...
0,0796
|
0,0796
…
0,0860
|
Функция распределения
|
0,08
|
0,16
|
0,27
|
0,42
|
0,58
|
0,73
|
0,84
|
0,92
|
0,97
|
0,99
|
Используя значение функции распределения, можно определить теоретическое число интересующих нас событий (число отказов в i
-м интервале) по формуле:
(13)
Определяем теоретическое число отказов в 1-м интервале: отказов.
Определим значения теоретических чисел для каждого интервала и заполним таблицу 8.
Таблица 8 – Значения теоретических чисел для каждого интервала
Функция распределения
|
0,08
|
0,16
|
0,27
|
0,42
|
0,58
|
0,73
|
0,84
|
0,92
|
0,97
|
0,99
|
Теоретическая
частота
|
8
|
8
|
11
|
15
|
16
|
15
|
11
|
8
|
5
|
2
|
Для закона распределения Вейбулла.
Рассуждая аналогично п. 1.7.2, вычислим не , а теоретические вероятности попадания СВ в -й интервал, например, вероятность отказа объекта в -м интервале по зависимости:
; , (14)
где a
,
b
-
параметры закона распределения, причем а
параметр масштаба, имеющий размерность случайной величины t
;
b
-
параметр формы (безразмерная величина);
-
смещение зоны рассеивания случайной величины t;
значения функции приведены в таблице Е.2[1].
Параметр определяют, используя коэффициент вариации. Из этого же приложения выбирают значения коэффициентов и :
Параметр рассчитывают по одному из уравнений:
или .
Пример решения для середины 1-го интервала:
Значения теоретических вероятностей запишем в таблицу 9.
Таблица 9 – Значения теоретических вероятностей
Середина интервала,
мм
|
0,025
|
0,031
|
0,038
|
0,044
|
0,050
|
0,057
|
0,063
|
0,070
|
0,076
|
0,082
|
Плотность функции распределения f(t)
|
0,2
|
0,55
|
0,78
|
0,84
|
0,84
|
0,74
|
0,57
|
0,48
|
0,32
|
0,19
|
Теоретическая
вероятность
|
0,034
|
0,095
|
0,135
|
0,146
|
0,146
|
0,128
|
0,099
|
0,083
|
0,055
|
0,033
|
Функция распределения Вейбулла имеет вид:
(15)
Данная функция зависит от двух аргументов – от параметра и обобщенного параметра . Ее значения могут быть вычислены непосредственно по зависимости (15) или определены по таблице (приложение Ж [1]). Входами в эту таблицу являются:
– значение параметра ;
– значение обобщенного параметра ,
где – значение случайной величины на конце i
-го интервала.
Вычислим функцию распределения на 1-м интервале:
Значения функции распределения запишем в таблицу 10.
Таблица 10 – Значения функции распределения
Границы
интервала,
мм
|
0,0220
...
0,0284
|
0,0284
...
0,0348
|
0,0348
...
0,0412
|
0,0412
...
0,0476
|
0,0476
...
0,0540
|
0,0540
...
0,0604
|
0,0604
...
0,0668
|
0,0668
...
0,0732
|
0,0732
...
0,0796
|
0,0796
…
0,0860
|
Функция распределения
|
0,050
|
0,148
|
0,286
|
0,443
|
0,598
|
0,732
|
0,835
|
0,907
|
0,951
|
0,977
|
Используя значение функции распределения, можно вычислить теоретическое число интересующих нас событий, например, число отказов машин в -м интервале по формуле:
(16)
где N
– общее число испытуемых (подконтрольных) объектов.
Определяем теоретическое число отказов в 1-м интервале:
Определим значения теоретических чисел для каждого интервала и заполним таблицу 11.
Таблица 11 – Значения теоретических чисел для каждого интпрвала
Функция распределения
|
0,050
|
0,148
|
0,286
|
0,443
|
0,598
|
0,732
|
0,835
|
0,907
|
0,951
|
0,977
|
Теоретическая
частота
|
5
|
9,86
|
13,78
|
15,74
|
15,45
|
13,38
|
10,34
|
7,16
|
4,48
|
2,53
|
По вычисленным значениям и для всех интервалов строят графики и , которые приведены в приложениях В и Г.
Результаты выравнивания опытных данных теоретическими законами распределения представим в виде таблицы 12.
Таблица 12 – Результаты выравнивания опытных данных теоретическими законами распределения
Границы
интервала,
мм
|
0,0220
...
0,0284
|
0,0284
...
0,0348
|
0,0348
...
0,0412
|
0,0412
...
0,0476
|
0,0476
...
0,0540
|
0,0540
...
0,0604
|
0,0604
...
0,0668
|
0,0668
...
0,0732
|
Середина интервала,
мм
|
0,025
|
0,031
|
0,038
|
0,044
|
0,050
|
0,057
|
0,063
|
0,070
|
Опытная частота
|
5
|
11
|
17
|
14
|
15,5
|
7,5
|
8
|
12
|
Дифференциальный закон
распределения
|
Опытная вероятность
|
0,05
|
0,11
|
0,17
|
0,14
|
0,155
|
0,075
|
0,08
|
0,12
|
Теоретическая
вероятность
|
НЗР
|
0,044
|
0,076
|
0,117
|
0,149
|
0,162
|
0,149
|
0,117
|
0,076
|
ЗРВ
|
0,034
|
0,095
|
0,135
|
0,146
|
0,146
|
0,128
|
0,099
|
0,083
|
Интегральный закон
распределения
|
Накопленная опытная вероятность
|
0,05
|
0,16
|
0,33
|
0,47
|
0,625
|
0,7
|
0,78
|
0,9
|
Функция распределения
|
НЗР
|
0,08
|
0,16
|
0,27
|
0,42
|
0,58
|
0,73
|
0,84
|
0,92
|
ЗРВ
|
0,050
|
0,148
|
0,286
|
0,443
|
0,598
|
0,732
|
0,835
|
0,907
|
Теоретическая
частота
|
НЗР
|
8
|
8
|
11
|
15
|
16
|
15
|
11
|
8
|
ЗРВ
|
5
|
9,86
|
13,78
|
15,74
|
15,45
|
13,38
|
10,34
|
7,16
|
1.7.3
Проверка правдоподобия (сходимости) опытного и теоретического законов распределения
Критерий Пирсона вычисляют по зависимости:
, (17)
где – опытная частота попадания СВ в i
-й интервал статистического ряда (берется из таблицы 4);
n
– число интервалов статистического ряда;
– значение функции распределения (интегральной функции) соответственно в конце i
-го и -го интервалов;
– теоретическая частота в i
-м интервале статистического ряда.
Делаем проверку для НЗР:
Делаем проверку для ЗРВ:
Значение критерия, вычисленное по зависимости (17) для НЗР , а для ЗРВ ; число степеней свободы , где n
– число интервалов статистического ряда, а m
– число параметров ТЗР (для НЗР и ЗРВ m
= 2); приняты уровень значимости (вероятность необоснованного отклонения гипотезы) . Необходимо выбрать ТЗР, наиболее адекватный распределению статистической информации.
По таблице В.2 приложения В [1] и k=5 определяем критическое значение -критерия: .
Сравниваем с . Так как только для ЗРВ, то делаем заключение о том, что выдвинутая гипотеза о сходимости опытного с теоретическим распределением ЗРВ с вероятностью не отвергается.
Для принятия окончательного решения определим вероятность подтверждения проверяемых ТЗР. Для этого опять используем таблицу В.2 [1]. Войдя в таблицу по этим значениям с учетом интерполяции определяем, что вероятность подтверждения выдвинутой гипотезы о ЗРВ в данном примере P
=19%.
Следовательно, в этой ситуации принимается гипотеза о том, что анализируемая статистическая информация с достаточной степенью достоверности подчиняется закону распределения Вейбулла.
1.8 Интервальная оценка числовых характеристик износов
Закон распределения Вейбулла.
В этом случае доверительные границы определяют по формуле:
, (18)
где - коэффициенты распределения Вейбулла, и выбираются из таблицы В.3 приложения В[1];
Следовательно:
- нижняя граница доверительного интервала;
- верхняя граница доверительного интервала.
С вероятностью можем утверждать, что истинное значение математического ожидания попадет в интервал от 0,0482мм до 0,0540мм.
1.9
Определение относительной ошибки переноса
Более правильно характеризовать точность оценки показателя надежности относительной ошибкой, которая позволяет корректно сравнивать объекты, в том числе и по разнородным показателям.
(19)
где – верхняя граница изменения среднего значения показателя надежности, установленная с доверительной вероятностью ;
– оценка среднего значения показателя надежности.
Вычислим относительную ошибку переноса:
Максимально допустимая ошибка переноса ограничивается величиной 20%, т.е. .
1.10 Определение числа годных и требующих восстановления
деталей
1) определим допустимые износы анализируемых деталей при их сопряжении с новыми и бывшими в эксплуатации деталями.
Для отверстия:
где – допустимый размер отверстия при сопряжении его с новыми деталями;
– допустимый размер отверстия при сопряжении его с деталями, бывшими в эксплуатации;
– наибольший предельный размер отверстия.
2) вычисленное значение допустимого износа отверстия отложили по оси абсцисс (Приложение Г). Из него восстановим перпендикуляр до пересечения с теоретической кривой износов . Полученную точку спроектируем на ось ординат и снять значение вероятности того, что детали окажутся годными (их восстановление не потребуется), при условии их сборки с новыми сопрягаемыми деталями. При этом число годных деталей может быть вычислено по зависимости:
(20)
3) выполняя аналогичные графические построения для значения , определяют число годных деталей при сопряжении их с деталями, бывшими в эксплуатации:
(21)
4) число деталей, требующих восстановления , определяется как
(22)
5) следует заметить, что большее практическое значение имеют не сами числа , , , а соответствующие коэффициенты, значения которых определяются ниже.
Коэффициент годности анализируемых деталей:
Коэффициент восстановления деталей:
=1-0,53=0,47.
Вывод
По значениям вычисленных коэффициентов можно сделать вывод,что необходимо более тщательно планировать производственную программу ремонтного предприятия по анализируемой детали.
|