
книги из ГПНТБ / Кацев, П. Г. Статистические методы исследования режущего инструмента
.pdfФункция распределения
|
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 ~ h° . |
(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-V2л |
2 |
\ s { b} t |
s {b) |
||
при tx |
< t |
< |
12. |
|
|
|
Введем |
обозначение |
|
|
|
||
|
|
|
hk — ftp. |
a == |
_b__ |
|
|
|
|
P = s{b) ’ |
|
s {b} |
|
Таким образом, распределение примет вид:
f(T) |
с0_ |
|
t*V 2я |
(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 2л |
|
|
|
|
|
||||||||
хехр |
|
|
|
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 мм можно считать нормальным.
Рассчитаем параметры нормального распределения, исходя из полученных анализом реализаций результатов.
7р |
г |
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