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

книги из ГПНТБ / Кацев, П. Г. Статистические методы исследования режущего инструмента

.pdf
Скачиваний:
33
Добавлен:
20.10.2023
Размер:
8.47 Mб
Скачать

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

 

F (Г) = Ф

Т — 77,7

 

 

 

КО, 12Га + 38,4

 

 

 

 

 

 

 

Рассчитаем

критерий

х 2 и его

вероятность. Получаем х 2 =

= 3,09; р (х2)

= 0,22.

Определим

период безотказной

работы

 

 

 

 

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

 

 

 

 

0,9= 1— Ф

Т — 77,7

 

 

 

 

К 0.12Г3 +

38,4

 

 

 

 

 

Рис. 15. Эмпирические и теоретические распределения функции и плотности ве­ роятности стойкости сверл диаметром

10,5 мм:

------ —----- выравнивание

по методу разде­

ляющих р а з б и е н и й :

--------------выравнивание

по реализации износа: —0 —0 — эксперимен­ тальные данные

отсюда Т0Л = 53 мин.

Таким образом, выравни­ вание распределения по за­ кону Бернштейна на основе реализации износа дало луч­ ший критерий 'согласия, чем по методу разделяющих раз­ биений (рис. 15). Различие в оценке надежности (по зна­ чению Т0Л) оказалось несу­

щественным. Альфа-распределение стой­

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

h (0 = h0 + Bt,

(83)

где h (t) — величина

износа (например, по задней поверхности)

в момент времени t\

h0 — начальная величина износа; В — ско­

рость износа.

Начальная величина износа h0 является величиной весьма

малой по отношению к конечной величине износа. Еще меньшим является различие между величинами начального износа для различных инструментов из выборки. Исходя из этого, зависимости для различных инструментов данной выборки могут быть пред­ ставлены в виде веера (см. рис. 14) с общим началом в точке (h0, 0),

которая называется полюсом. Такие линейные случайные функ­ ции, имеющие полюс, называются веерными.

60

Математическое ожидание функции (83)

М {Л (0} = hQ+ М \ В } t,

где М {В} — математическое ожидание скорости изменения па­

раметра..

Если отложить по оси ординат величину hk, соответствующую

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

t = hk ~ .

(84)

В

 

При нормальном распределении случайной величины последняя может принимать любые значения от —оо до -foo. Поскольку возможные значения случайной величины скорости износа В

могут быть только положительными и ограничены интервалом (Ьи Ь2), то распределение В будет усеченным нормальным с плот­

ностью

1 {В) — cf (В),

(85)

где f (В) — плотность неусеченного распределения; с — норми­

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

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

f(B) = cf(B)

С

ехр

( б - б ) а

(86)

 

s {Ь} ]^2я

 

2sz {b]

 

при Ьг < b < Ь2-

Время Т достижения определенной (принятой за критерий

затупления) величины износа режущей кромки инструмента за­ висит согласно формуле (84) от скорости износа В. Распределение величины В известно, оно является усеченно-нормальным (86).

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

дует найти распределение функции от случайной переменной. Из математической статистики [14] известно следующее пра­

вило о распределении функции от случайной переменной. Если

случайная

величина

В

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

f (В),

а случайная величина

Т связана с В зависимостью Т =

ф (В)

[в нашем

случае Т — ^ h/t д ll°^ J , то плотность распределения

случайной величины

Т

имеет вид:

 

 

 

f(T) = f[B{T)]B'(T),

(87)

61

где В (Т) — функция, обратная ср (В) в нашем случае В (Т) =

-

Поэтому на основе зависимостей (84), (86) и (87) после преоб­ разований получаем: *

 

/

СП

(hk ho) с exp

1

hk — hp

b ^

 

i (6) C-V

2

\ s { b} t

s {b)

при tx

< t

<

12.

 

 

 

Введем

обозначение

 

 

 

 

 

 

hk — ftp.

a ==

_b__

 

 

 

 

P = s{b)

 

s {b}

 

Таким образом, распределение примет вид:

f(T)

с0_

 

t*V

(88)

(89)

где Р — относительный запас стойкости инструмента, имеет раз­ мерность времени; a — относительная средняя скорость износа или коэффициент однородности изменения скорости износа, без­ размерная величина.

Распределение (89) называется альфа-распределением.

Заметим, что а =

7 ^

= ■— » гДе иь — коэффициент вари­

ации скорости износа

[22].

Параметры а и р позволяют учиты­

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

В качестве примера приведем альфа-распределение стойкости сверл 4,2 мм, испытанных в количестве 108 шт. Средняя стой­

кость Т = 107 мин. Параметры распределения а — 0,48; Р =

=100; с = 1,5.

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

 

 

150

 

 

( i ^ - 0 , 4 8 ) 2]

 

 

 

](Т) = Т2 У 2л ехр

 

 

 

Расчет

критерия

согласия

 

дает

Р (%2) =

0,39.

Надежность

сверл

оценивается на основе зависимости

 

 

 

 

 

Р(Т)

Ф (а) +

Ф (z)

 

 

 

 

 

 

1 + Ф (а) ’

 

 

 

 

 

 

 

 

 

где Ф ( ) — функция

интеграла

вероятностей;

z =

О

-^---- a.

В

нашем случае

стойкость

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

Гд.д =

51

мин.

 

 

 

 

 

 

Подробнее см. работу [34].

62

Стойкость резцов, оснащенных твердым сплавом, также опи­ сывается альфа-распределением. В качестве примера приведем альфа-распределение стойкости резцов, оснащенных твердым спла­

вом Т14К8 на операции подрезки торца, материал сталь

ШХ15,

v =

134 м/мин,

s =

0 , 1 m m / o 6 ,

 

 

 

 

t

=

0,5

 

мм.

 

 

 

 

 

_

 

 

 

 

 

 

Средняя

 

стойкость

Т =

 

 

 

 

= 125 мин. Параметры рас­

 

 

 

 

пределения а — 3,2;

р =370;

 

 

 

 

с = 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(T) =

 

370

 

X

 

 

 

 

 

 

 

г- V

 

 

 

 

 

хехр

 

 

 

370

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Критерий

согласия имеет

 

 

 

 

вероятность

Р (%2) =

0,29.

 

 

 

 

Стойкость

с

 

вероятностью

 

 

 

 

Р = 0,9

будет

составлять

 

 

 

 

Тол = 82,6

мин.

приведено

 

 

 

 

 

’ На

 

рис.

16

 

 

 

 

альфа-распределение стойко­

 

 

 

 

сти

двух партий

сверл диа­

 

 

 

 

метром 5 мм: со

 

стандартной

 

 

 

 

и

увеличенной

 

толщинами

 

 

 

 

сердцевины.

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

 

 

 

 

резко

отличаются

.друг от

0

 

 

8

друга, наглядно показывают,

 

 

что увеличение толщины серд­

Рис. 16. Лльфа-распределение стойкости

цевины

 

не

только

увеличи­

 

двух партий сверл

диаметром

5 мм и

вает стойкость,

но

и увели­

толщиной сердцевины:

 

чивает надежность.

Заметим,

1 стандартной',

2

увеличенной

что

при а >

4,5

с

 

погреш­

 

 

 

 

ностью

не

более

5%

можно

 

 

 

 

среднюю

стойкость

 

принять

равной медиане

и

рассматривать

a -распределение как частный случай распределения Бернштейна. При а >> 10 альфа-распределение может быть заменено нор­

мальным.

Анализ реализаций износа

Важнейшими факторами физического содержания процесса, влияющими на закон распределения, являются модели изнашива­

ния. Рассмотрим этот вопрос применительно к режущему инстру­ менту.

Выше рассматривались случайные величины, характерные тем, что в результате опыта они принимали некоторое одно, заранее

63

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

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

врезультате опыта может принять тот или иной конкретный вид

[14].Конкретный вид, принимаемый случайной функцией в ре­ зультате опыта, называется реализацией случайной функции. Если над случайной функцией произвести группу опытов, то получим группу или семейство реализаций этой функции.

Встречающиеся на практике реализации процессов износа режущего инструмента близки к линейным: скорость износа каждого инструмента в основном постоянна, а для любого взятого инструмента скорость износа случайна — для каждого инстру­ мента своя. Для описания процессов износа инструмента поэтому целесообразно применить линейные случайные процессы, все реализации которых являются прямыми линиями. Линейные случайные процессы являются весьма удобной моделью процес­ сов износа инструмента, близки к действительным процессам из­ носа; очень просто описывают основные особенности процесса износа; требуют минимального количества экспериментальных данных для вычисления характеристик случайного процесса; дают возможность наиболее просто исследовать надежность эле­ ментов при их износе.

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

рического материала.

Рассмотрим типичные примеры реализаций износа инструмента. На рис. 17 показаны реализации износа метчиков М12 со шли­ фованным профилем резьбы, изготовленных из стали Х6ВФ с раз­

64

личной твердостью: HRC 54—55; HRC 57—58 и HRC 60—62.

Скорость резания 2,07 м/мин, длина нарезаемого отверстия 14,5 мм, обрабатываемый материал—сталь 45. По мере увеличе­ ния твердости не только повышалась стойкость метчиков, но и менялся характер реализаций — от типичной веерной до полного

переплетения.

На рис. 18 представлены реализации износа проходных рез­ цов, оснащенных сплавом Т14К8 при обработке стали ШХ15

О

100

200

300

000

500

600

700

 

Количество

нарезанных

отверстий

 

Рис. 17.

Реализации износа метчиков М 12 с различной твердостью

(случайная выборка 6 реализаций из 50 испытаний). Характерным является отсутствие заметного участка приработки и участка ускоренного износа за пределами критерия hk — 0,6 мм. Среднее

значение скорости износа является постоянной, реализации могут быть аппроксимированы прямой линией. Имеет место пере­ мешивание реализаций.

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

Рассмотрим анализ реализаций износа сверл 8 мм из быстроре­ жущей стали Р9К.5. В табл. 13 приведены данные приращений из­ носа, которые измерялись через каждые 100 шт. просверленных отверстий (14 мин работы) до 84 мин работы, т. е. число интервалов

измерений износа

было т =

6. Расчеты сведены в

табл. 14.

1.

Для каждого сверла вычисляем среднюю

величину прира

щения

износа за

интервал времени 8 = 14 мин

 

 

 

6Л(0 =

— £ 8 /i}° .

 

 

 

 

1=1

 

5

П. Г. Кацев

65

Например,

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

i =

1

W l) = 4

- ( 0 ,15 + 0,1 + 0 ,1 +

0 +

0,2 + 0,15) = 0,117мм.

2. Вычисляем общую среднюю величину приращения износа Ш

^ l==- r f i ^ <IJ =

• 2,7°9 = 0,0903.

Рис, 18. Реализации износа проходных резцов из сплава Т14К8:

1—6— данные для шести инструментов (и = 136 м/мин: t = 0,5 мм;

s— 0,4 мм/об)

3.Вычисляется дисперсия s- величин бМ°, их сумма и средн

дисперсия приращений

sht

 

 

 

 

 

2

1

T S

(бй{-° — б/г(г))2;

2 s? =

0,1231,

sr-

 

 

 

/=i

 

 

 

 

 

 

Л .

(171

^ S Si

-|

т

\\

ГвТЛП

=ёт\2.

Sfift -

-

и»_1

+

- 7 ^ - Г

А ( &г -

6h) ’

 

 

 

 

 

 

 

1=1

 

 

21 (б/г<1) — 8h)2 — 0,010035,

 

 

L—l

 

 

 

 

 

 

 

2 _ _

( 6 —

1)

0,1231

 

6 - 30 — 1

0,01 — 0,00377.

Sfift

6 - 30 — 1

 

66

Таблица 13

 

 

 

 

 

Приращение и зн о са 6Л^1)

по интервалам ( t= 1 4

мин)

 

для сверл 8 мм

 

 

 

 

 

 

 

 

Интервал /

 

51

 

t 14 мин

It

31

4<

Ы

;

 

 

 

 

 

сверл

Приращения износа 6 f t

в мм

 

 

 

 

1

0,15

0,10

0,10

0,00

0,20

0,15

2

0,20

0,10

0,05

0,05

0,10

0,15

3

0,15

0,10

0,15

0,10

0,10

0,10

14

0,15

0,10

0,15

0,20

0,00

0,10

15

0,10

0,00

0,05

0,15

0,10

0,10

16

0,10

0,10

0,05

0,05

0,10

0,00

28

0,20

0,15

0,00

0,15

0,00

0,05

29

0,10

0,10

0,00

0,10

0,10

0,05

30

0,10

0,00

0,10

0,00

0,00

0,15

Таблица 14

Анализ реализаций износа сверл 0 8 мм

сверл i

6ft<‘ )

 

■?

1S s 9c

 

 

 

 

1

0,117

 

0,00468

3,6702

2

0,108

 

0,00252

3,4014

3

0,117

 

0,00068

4,8325

14

0,117

 

0,00468

3,6702

15

0,083

 

0,00268

3,4281

16

0,067

 

0,00168

3,2253

28

0,092

 

0,00742

3,8704

29

0,075

 

0,00178

3,2504

30

0,058

 

0,00442

3,6454

6 h — 0,0903

2 s ? =

0,1231 £

lgs? =

73,4581

^

( 6ft(n

6 ft)2 =

0,010035

 

б*

67

4. Теперь можно установить, противоречат ли данные испыт ний предположению об однородности начального качества объ­ ектов.

Для этого надо убедиться, что между величинами s2, s£..........s2

и между величинами 6/i(1>, б/г(2>, . . ., <5/г("> нет существенной раз­ ницы. Для проверки гипотезы о равенстве дисперсий применим

критерий Бартлетта

[19].

 

 

Вычисляем величину х а

 

X2

 

2,3026

п (т — 1) X

1

 

 

3п (т — 1)

 

 

 

 

 

 

 

 

S -?

 

 

п

 

 

 

 

X

Ig

i=I

 

II

 

 

 

 

 

 

 

 

 

 

при

числе степеней свободы f =

п — 1

= 30 — 1

= 29.

 

 

о

 

2,3026

 

3 0 (6 — 1) X

 

 

 

Х" =

1

3 0 +

1

 

 

 

 

 

 

 

 

 

 

 

3-30 (6 — 1)

 

 

 

 

X [lg0,0041 -

- 733’04581 ]

= 19,8 .

 

По таблице

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

х а

[11]

находим

для f = 29 и

р =

0,05, что

Хкр =

42.

Так

как

%2 < Х к Р. то,

следовательно,

опытные данные не противоречат гипотезе о равенстве дисперсий.

5.

Проверка существенности различия средних величин пр

ращений

износа.

 

 

 

 

 

 

 

 

Подсчитывается сумма квадратов отклонений между сериями

испытаний сверл

 

П

 

__

 

 

 

 

 

 

 

 

 

 

 

Qi =

т 2

(Ж (0 — т у = 6 ■0,01 = 0,06

 

 

 

;=х

 

 

 

 

 

при числе степеней свободы f = п — 1 =

30 — 1 = 29.

Подсчитывается сумма квадратов отклонений внутри серий Q2

 

 

 

П

sf = (6•— 1)• 0,1231 =

 

 

Q2 = (m— 1) 2

 

0,616

при кг =

п (т — 1) =

/=1

 

 

 

 

 

30 (6 — 1) =

150.

 

 

Находим величину

критерия F:

 

 

 

 

 

F _

Qi/fti

_

0,06/29

_

0,0021

_

П с

 

Г

Q2/k2

~

0,616/150

0,0042

~

U,D-

По таблице работы [11] находим Fq — 2,2. Так как F <i Fq,

то гипотеза о равенстве средних принимается.

68

6.Так как опытные данные не противоречат предположени

об однородности

начального

качества

объектов, вычислим:

 

=

„2

п n n 'JB

 

 

jbft_=

0,0038

0,042,

 

 

ш ~

0,09

 

 

 

?

б

1

0,0081

0,213,

 

s 2

' 8 ~

0,0038-10

 

s6h

 

 

 

hk — 0,5 мм (предельно допустимый уровень износа)

h

0,5

12, 0.

У

0,042

 

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

Т-з

3025

15,4;

■> т

55

г

196

К

J96 = 0,28.

S2

 

 

 

Расчет критерия %2 согласия эмпирических данных с гамма распределением дает %2 = 1,79; Р (х2) = 0,41.

7. Можно считать, что распределение нормально, если

г/к

: 3,5.

V~rfk

В нашем случае

Следовательно, распределение стойкости сверл 0 8 мм можно считать нормальным.

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

г

12,0

сс

т = Т = 0^ 13- =

56 м и н -

=s = 16.3 шш-

Значенияи s, вычисленные по результатам стойкостных испытаний, Т — 55 мин; s = 14 мин.

Расчет критерия согласия х2 дает х2 = 0,49; Р (х2) = 0,78.

Проверка стационарности экспериментальных данных потока отказов инструмента [19]. Поток отказов является стационар­ ным, если закон распределения группы случайных величин К (ti), К (tz), . . ., К (tn) совпадает с законом распределения случайных

величин К {t1 + а) К (а); К (*,

+ а) К (а), . . .,

К (tn +

+ с) — К (а), т. е. распределение

случайных величин

К (t) не

69

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