Файл: Кацев, П. Г. Статистические методы исследования режущего инструмента.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 166

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

принята как минимальное значение стойкости при испытаниях, равное 60 мин (см. рис. 8). Выравнивание на вероятностной бу­

маге (см. рис.

9) дает те же значения для Т и Т 0.

 

Величина

параметра

X будет

 

 

 

X

 

1

_

1

0,0059.

 

 

f

_, f o

230 — 60

 

 

 

 

 

Расчет критерия согласия дает значения %а = 0,26

и Р (х2) =

= 0,99. Значение стойкости

с вероятностью р — 0,9

подсчиты­

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

величину

порога

чувствительности Т 0.

Таким

образом, Р (Т)

=

1 — F (Т) е~ °.°059, отсюда То.д —

= 20 мин,

То,э =

То,э + То =

80 мин.

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

булла-Гнеденко

[141, [19]

применяется при анализе прочности

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

' _т]_

/(Т) = е Р ^ - Т 7-1

при Т ^ 0

(58)

0

при Т < 0.

 

Это распределение является асимметричным: в частном случае (при у = 1 ) оно превращается в экспоненциальное распределение.

Оценка параметров распределения Вейбулла-Гнеденко. Для проверки возможности использования этого распределения при описании стойкости инструмента обработаны многие испытания различных партий инструмента. Покажем на примере оценку па­ раметров этого распределения. Испытано 49 шт. сверл диаметром 3,5 мм. Испытания проводились до поломки сверл. Значения стойкости и частостей откладываем на вероятностную сетку рас­ пределения Вейбулла-Гнеденко и по полученным точкам проводим прямую I (рис. 10). Затем через точку А проводим линию, парал­

лельную прямой /, до ее пересечения с вертикальной линией Я в точке В. Ордината точки В по шкале Т дает значение у = 0,8. Для оценки (3 находим ординату /п0 = 1,2 на шкале пг точки С

пересечения прямой 1 с линией Я. Параметр (3 находится так:

(59)

где / — порядок множителя

шкалы

времени (10/)

 

В нашем случае / =» 0

 

 

 

 

в = е1.2.10°-°,8 = е1.2=

3,3201;

4 - =

0,8

= 0,24.

3,3201

 

 

р

 

 

44


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

2*0,8

}(Т) = е~~*&-0,24Т~0-2.

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

Z

т-0,8

F (Г) = J }(Т) d t = 1 — е~ зЖ.

Рис. 10. Вероятностная сетка распределения Вейбулла-Гне­ денко (стойкость сверл диаметром 3,5 мм)

Рассчитаем критерий = 7,23 (табл. 7), находим Р (%2) =

= 0,2. Рассчитаем надежность указанных сверл, используя формулу

ту

Р (Т) = 1 — F (Т) = е Э .

Определим значение

стойкости

с

вероятностью Р = 0,9

0,9 = е

_ z l

=

0,3 мин.

3|32; 7 , о ,9

Оценка параметров распределения Вейбулла-Гнеденко методом разделяющих разбиений. По результатам испытаний известны ве­ личины Т и Т 2, . . ., TN стойкости инструментов в количестве

45

Таблица 7

 

 

 

 

 

 

 

 

 

 

 

 

Расчет распределения а крит ерия %2 д л я

стойкости

сверл 3,5

м м

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кри-гернй у,2

 

Интер.

т,

Г

 

?

ту

 

т*

F( T)

f (Т)

 

 

валы

2 nk

3,32

с

3,32

(т , — (га, —т г)2

 

 

 

 

 

 

 

 

 

 

 

 

ГПп (т , —т г)

 

 

 

 

 

 

 

 

 

 

—т 2)2

т 2

0—2

13

1

13

0,26

0,30

0,74

0,26

0,26

13

 

 

 

2—4

14

3

27

0,55

0,72

0,49

0,51

0,25

12

2

4

0,33

4—6

9

5

36

0,73

1,09

0,34

0,66

0,15

7

2

4

0,57

6—8

5

7

41

0,83

1,43

0,24

0,76

0,10

5

8—10

1

9

42

0,85

1,75

0,17

0,83

0,07

3

2

4

1,33

10—12

3

11

45

0,91

2,04

0,13

0,87

0,04

2

1

1

0,50

12—14

1

13

46

0,93

2,33

0,09

0,91

0,04

2

1

1

0,50

14—16

3

15

49

1,00

2,63

0,07

0,93

0,02

1

2

4

4,00

.

7.2 = 2 = 7,23


N шт. Оценка параметров у и Р по методу разделяющих разбиений

производится следующим образом.

1. Выбирают границы разбиения ©! и ©2, при этом©! <С©2

и подсчитывается количество значений 7\, лежащих в пределах (0( ©.), / = 1,2. Обозначим эти количества т (©/)■

т (0;)

Находим отношения v (©,) = — ^— .

2.

Параметры

у

и

Р

оцениваются

по

формулам

 

 

 

 

~

 

I n i n

 

 

 

 

 

 

1

 

 

 

 

 

1 — V (0 i ) J

 

1

L 1 — v ( 0 2) J

 

 

 

 

у

 

 

 

 

 

 

 

 

 

 

In 0 ! — In 0

2

 

 

 

 

(60)

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

I

 

 

 

 

 

 

 

 

 

 

н

Г

 

1

 

 

 

 

 

 

 

 

 

 

 

 

In

 

V (0 !)

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Рассмотрим этот метод на

примере

испытаний

50

плашек

М 10X1,5 (табл.

8).

 

 

 

 

 

 

 

 

 

 

Таблица 8

 

 

 

 

 

 

 

 

 

 

 

 

 

Стойкость п ла ш ек M l О у. 1,3

 

 

 

 

 

 

 

 

Стой­

 

Стой­

 

Стой­

Л'Ь

Стой­

 

Стой­

кость

 

кость

 

кость

кость

 

кость

по пор.

в мин

по пор.

в мин

по пор.

 

в мин

по пор

в мин

по пор

в мин

1

1,30

11

 

5,00

21

 

9,70

31

 

14,60

 

41

25,20

2

1,60

12

 

6,40

22

10,00

32

 

16,52

 

42

25,50

3

2,20

13

 

6,44

23

10,50

33

 

16,81

 

43

25,60

4

2,40

14

 

6.60

24

10,60

34

 

17,30

 

44

25,70

5

3,30

15

 

7,00

25

 

11,60

35

 

17,80

 

45

28,00

6

3,60

16

 

7,70

26

 

12,40

36

 

18,50

 

46

31,00

7

4,00

17

 

8,30

27

 

12,40

37

 

18,58

 

47

31,30

8

4,80

18

 

8,60

28

 

13,40

38

 

19,20

 

48

43,80

9

4,80

19

 

8,90

29

 

14,20

39

 

22,40

 

49

47,00

10

4,90

20

 

9,20

30

 

14,60

40

 

24,20

 

50

47,00

Принимаем границы разбиений ©х = 8 мин и ©2 = 23 мин.

Соответственно

количество

плашек,

имеющих стойкость

от 0

до 8 мин и от 0 до 23 мин, будет т (0 Х)

= 15

и т (0 2)

=

39.

Находим

отношения

 

 

 

 

 

 

 

 

v(©!) =

т (®i)

15

г> о.

v (0 2)

т (02) _ 39

0,78.

N

50

U,d’

~~N

— "50

Подсчитываем

значения

параметров

у и

Р

 

 

у,

 

—‘ 1,031 — 0,415

 

 

£

81,4

51,5.

 

 

У —

2,08 — 3,13

 

— 1,4>

Р —

0,357

 

 

 

 

 

 

Подсчитаем критерий %2 — 7,57, Р (%а) = 0,11.

Гамма-распределение времени безотказной работы [19]. Рас­ смотрим идеализированную схему износа, обладающую следу­

'47


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

Плотность этого распределения имеет следующий вид:

 

 

 

f (Т) -

КТ'-'егМ- 5*0.

(61)

где

г — число

повреждений,

необходимых

для возникновения

отказа,

Г (г) — гамма-функция,

определяемая формулой

 

 

 

 

 

00

 

 

 

 

 

 

Г (/') =

|

xT~le~x clx.

 

 

 

 

 

г справедливо

о

 

 

 

Для

целых

соотношение

 

 

 

 

 

Г (г) =

( г - 1)!

 

 

Пусть hk есть предельно допустимый уровень износа, т. е.

при

11 (t) 5 »

hk наступает отказ. Отсюда следует,

что число поврежде­

ний

до

отказа

определяется

соотношением

 

 

Величина X задает среднюю скорость износа. Функция распре­ деления F (Т) — гамма-распределение задается в соответствии

с формулой (60) соотношением

т

т

 

F (Т) = ff (t) di =

- ^ r - JV - 4 - и di.

(62)

о

; 0

 

Для целых г интегрированием по частям можно получить

 

г—1

 

F (Г) = 1 —

е-лг

(63)

 

k\

 

k=Q

Математическое ожидание М (Т) и дисперсия гамма-распреде­ ления D {Т} выражены равенствами:

М[Т} = ^ ;

(64)

Параметры г и X могут принимать только положительные зна­ чения. Легко показать дифференцированием, что функция / (Т)

имеет один максимум в точке Т —■ г~^-- - , если /\> 1 . Если 0<J

48


< /- < 1, функция f (Т) является непрерывно убывающей. В част­ ном случае, когда г — 1, плотность вероятности сводится к выра­

жению

/ (Т) = Хе~ w,

т. е. к плотности экспоненциального распределения. Это согла­ суется с моделью возникновения экспоненциального распределе­ ния — мгновенный отказ при первом повреждении. С возраста­ нием г гамма-распределение стремится к нормальному. Физически

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

/■>12, или если отношение— — -----7= - > о,о, т. е. если

Vw

вариация v < -уу #=« 0,28.

Оценка параметров Я и г на основе данных о времени безот­ казной работы. Пусть имеются результаты Т ъ Т 2, . . ., TN испы­ таний на стойкость N инструментов. Рассчитываем в известном

порядке среднее значение Т и дисперсию s?. Затем, приравнивая

Т = М {Т} и D [Т] = s?, получаем соотношения, из которых можно найти X и г

Х = Т _

Т2

(65)

 

Рассмотрим выравнивание по гамма-распределению стойкости резцов Т15КЮ по результатам производственных испытаний на операции окончательной обточки колец роликоподшипника (сталь ШХ15) на 1ГПЗ. Скорость резания 135 м/мин, подача 0,4 мм/об, охлаждение — эмульсия. Критерий затупления — износ по зад­ ней поверхности 0,6 мм. Рассчитаем параметры гамма-распределе­ ния (стойкость в штуках обработанных колец).

'? __

Т

__

245,1

п

 

s2 (Т)

~

71,682 -

и ’и о -

 

7 2

 

2 4 5 ,I 2

11,7.

 

s 2 {T]

~

71.683

 

 

Получаем функцию плотности вероятности

1 1(Т) = г (11,7) -0,0511-7Г10’7е-°-05Г.

4 П. Г. Кацев

49