CC BY
37
УДК 621.95.01
СТАТИСТИЧЕСКОЕ ИССЛЕДОВАНИЕ СИЛОВЫХ ХАРАКТЕРИСТИК ПРИ ОБРАБОТКЕ СИНТЕГРАНА СВЕРЛЕНИЕМ
В.А. Рогов, А.К. Велис, В.В. Копылов
Кафедра технологии машиностроения, металлорежущих станков и инструментов
Инженерный факультет Российский университет дружбы народов ул. Орджоникидзе, 3, Москва, Россия, 115419
В статье рассматриваются исследования проведенных на синтегране для оценки влияния режимов резания сверлением на выходные силовые характеристики станка: осевая сила Pz и крутящий момент Mz для минимизирования напряжений узлов станка. Исследования проведены с помощью спиральных твердосплавных сверл из сплава BK8, поверхность которых была обработана магнетронным напылением сплава NbHfTi для получения градиентного по содержанию легирующих элементов промежуточного слоя перед нанесением износостойкого покрытия (TiAl)N. Исследования проводились на обрабатывающем центре, оснащенном ЧПУ, с возможностью диагностирования силовых характеристик станка в режиме реального времени. Исследования показали оптимальные значения частоты вращения шпинделя станка и осевой подачи для сверления син-теграна.
Ключевые слова: синтегран, сверление синтеграна, сверла ВК8, модифицирование поверхности, компьютерная диагностика процесса резания, износостойкие покрытия, осевая подача, крутящий момент.
Целью настоящей работы является анализ работы спиральных твердосплавных сверл из сплава ВК8 (К20 по ISO), прошедших комплексную поверхностную обработку магнетронным напылением сплавом NbHfTi, и получение математических уравнений регрессии, описывающих зависимости выходных параметров Pz, Mz от входных параметров — частоты вращения шпинделя n и осевой подачи S.
График для осевой силы Pz, полученный с помощью сигналов системы ЧПУ станка, показан на рис. 1, где четко видны 8 разных отрезков от разных опытов при сверлении до глубины 20 мм на различных режимах.
Pz, % 22
17 12
7
2
0 5 10 15 20 25 t, сек
Рис. 1. График Pz — t, полученный с помощью ЧПУ станка
»4I [l \
-1..........
График крутящего момента Мг (рис. 2), также полученный с помощью системы ЧПУ станка, показывает неустойчивые значения крутящего момента, меняющегося при увеличении глубины отверстия. Это изменение может объясняться чувствительностью датчиков к этому показателю и неравномерностью свойств синтеграна. Значения момента М2 вдоль отверстия будем считать нормально распределенными, следовательно, их можно изучать статистически.
мг, % 27
22
17
12
7
2
0 5 10 15 20 25 t, сек
Рис. 2. График Рг — t, полученный с помощью ЧПУ станка
Запишем матрицу планирования для нашего полнофакторного эксперимента для определения математической модели шероховатости и силовых характеристик при обработке синтеграна.
Результаты опытов, проведенных на синтегране, приведены в табл. 1.
Таблица 1
Результаты опытов
Опыт п в Р* Pz2 Уu М21 Уu
Х1 Х2 Уи1 Уи2 Уи1 Уи2
1—5 1 600 40 13,70 14,75 14,22 10,99 10,08 10,54
2—6 400 40 16,07 15,93 16,00 11,25 11,72 11,49
3—7 1 600 5 15,19 16,16 15,67 6,81 5,69 6,25
4—8 400 5 17,55 17,48 17,51 8,72 7,68 8,20
Полнофакторный эксперимент типа N = 2к проведен на основе однофакторных экспериментов, проведенных в работе [1]. Суть его состоит в проведении эксперимента по заранее заданной схеме (матрица планирования) [2—5].
Матрица планирования полнофакторного эксперимента типа ПФЭ 22 для обработки синтеграна сверлением, учитывающая эффект взаимодействия, приведена в табл. 2.
Интервал варьирования и основной уровень выбрали из области определения [2; 1]. Для исключения систематических ошибок проведены два эксперимента с рандомизацией каждого опыта.
|| 1
/VI
Таблица 2
2
Матрица планирования эксперимента 2 с параллельными опытами и с эффектом взаимодействия факторов
Опыт Фактор Х1Х2 Переменная состояния
Х 0 Х1 Х2 Уи1 Уи2
1—5 + 1 + 1 + 1 У11 У21
2—6 + 1 -1 + 1 -1 У12 У22
3—7 + 1 + 1 -1 -1 У13 У23
4—8 + 1 -1 -1 + 1 У14 У24
Предполагаемая математическая модель принята линейной с эффектом взаимодействия между факторами, уравнение (1).
У = Ь0 + Ь1 *1 + Ь2*2 + Ь12 *1 ^ (1)
где у — функция отклика факторов x0, x1, x2; Ь0, Ь1, Ь2 — коэффициенты уравнения регрессии; Ь12 — коэффициент взаимодействия между факторами х1 и х2.
Таблица 3
Соответствие уровней и факторов для сверления синтеграна
Уровень Фактор
п, мм/об в, мм/мин
Х1 Х2
Верхний + 1 1 600 40
Основной 0 1 000 22,5
Нижний -1 400 5
С помощью матрицы планирования (см. табл. 2) вычисляем коэффициенты Ь0, Ь1, Ь2 по формулам (2), (3), (4):
1 Ы
Ь = -ЕХшуи; (и = 1,2,..., N), (2)
^ и=1
1 N _
Ьу = NЕхиХ]иУи; (и =12,..., N);(г Ф у), (3)
^ и=1
1 М
где уи =—Е Уик; (4)
М к=1
( Х,ы — кодированное значение фактора, уи — среднее значение по параллельным опытам к-х строк матрицы планирования, и — порядковый номер строки матрицы или номер опыта, г, у — номер колонки, г Ф у, М — число повторных опытов).
После расчета коэффициентов уравнения регрессии по формулам (1) и (2) запишем все значения в табл. 4.
Уравнения регрессии для Pz, Mz будут иметь следующий вид:
Pz ^ у = 15,85 - 0,9Ц - 0,74х2 + 0,02х1х2, (5)
Mz ^ у = 9,12 - 0,73^ +1,90х2 + 0,25X1 х2. (6)
Таблица 4
Коэффициенты уравнений регрессии для Рг, Мг
Р? Мг
Ьо 15,85 9,12
ь, -0,91 -0,73
Ь2 -0,74 1,90
Ь1, 2 0,02 0,25
Проведем статистический анализ математической модели, который включает
1) проверку однородности дисперсии или оценки дисперсии воспроизводимости;
2) проверку коэффициентов регрессии; 3) проверку адекватности математической модели.
Для статистического анализа нам нужно рассчитать все построчные дисперсии параллельных опытов (уравнение (7)) и записать в табл. 5.
1 и
S2u =ХГ-ГI(Ук -Уи)2, (7)
И 1 к=1
где: (И - 1) — число степеней свободы; М — число повторных опытов; Б^ — дисперсия выборки; и — порядковый номер строки параллельного опыта.
Таблица 5
Значения дисперсии для Рг, Мг
и ви, - для Р1 М
1 в2 = 1/(2 - 1)(13,70 - 14,22)2 + (14,75 - 14,22)2 = 0,55 0,41
2 вЦ = 1/(2 - 1)(16,07 - 16,00)2 + (15,93 - 16,00)2 = 0,01 0,11
3 в2 = 1/(2 - 1)(15,19 - 16,67)2 + (16,16 - 15,67)2 = 0,47 0,63
4 в2 = 1/(2 - 1)(17,55 - 17,51 )2 + (17,48 - 17,51)2 = 0,00 0,54
N 2А2 и=1 1,03 1,70
Проверка однородности дисперсий
Используем критерий Кохрена для проверки однородности, который является соотношением между максимальной дисперсией и суммой всех дисперсий. Математически выражается таким образом:
(Би ) 0 55
Срасч = = Щ = 0,53. (8)
1 и=1Би
Таблица 6
Расчетный критерий Кохрена для Рг, Мг
с расч
Рz М1
0,53 0,24
Дисперсия является однородной, если расчетные критерии Кохрена не превышают табличных.
Если Срасч < СТабл, (д,/1,/2), дисперсии однородны.
Если Срасч > Стабл, имеются отклонения в результатах опытов вследствие систематических ошибок и неучтенных факторов.
/ — число степеней свободы (/ = М - 1 = 2 - 1 = 1); /2 — число измерений переменной у /2 = N = 4); М — число параллельных опытов (М = 2). Так как / = 1 и /2 = N = 4 ^ из табл. 6 Стабл = 0,91.
Табл. 7 показывает Срасч и Стабл, рассчитана по формуле (8) для Р2, М2 и Я2. Р2 ^ 0,53 < 0,91; М2 ^ 0,24 < 0,91 ^ дисперсия однородна.
Таблица 7
Сравнение между Срасч и Ста6л для Р1, М2
Стабл с расч
Pz М2 Р2 М2
0,91 0,91 0,53 0,37
Вычислены значения ошибки опыта или дисперсий воспроизводимости по формуле (9) и занесены в табл. 8.
1 N 1
502 =1Е Б2и = -(0,55 + 0,01 + 0,47 + 0,00) = 0,26. (9)
N п=1 4
Таблица 8
Значения ошибок опыта для Р , М
Э2 — ошибка опыта или дисперсия воспроизводимости
Р2 М2
0,26 0,42
Расчет дисперсии коэффициентов уравнения регрессии Б^ дан в табл. 9.
4 = (10)
где Б ы = л[52 — среднеквадратичное отклонение; Б^ = 026 = 0,065.
Таблица 9
Значения ошибок опыта для Рг, Мг
я2 Я»
Р. Мг Рг Мг
0,065 0,106 0,25 0,33
Проверка значимости коэффициентов уравнений регрессии
Проверим значимость коэффициентов уравнений для Р2, М2 из табл. 4. Для проверки равенства средних значений в двух выборках (двух повторных опытов) используем критерий Стьюдента или ¿-критерии.
При нулевой гипотезе предполагается, что средние значения равны.
I = Ы(М - 1) = 4 — число степеней свободы;
q = 10% ^ ¿табл = 2,13 — уровень значимости;
(¿г)расч > ¿табл^,I) — условие значимости коэффициентов.
Ъ )расч = ^, (11)
"Ы
где (¿¿)расч — расчетное значение критерия Стьюдента, табл. 10. Бы — среднеквадратич-
Таблица 10
"и расч
ное отклонение.
Расчетное значение критерий Стьюдента для Р2, М2
Опыт ^расч
Р. Мг
£ расч 62,39 27,98
£ расч 3,56 2,22
£ расч 2,91 5,81
^расч 0,06 0,77
Если (¿¿)расч > ¿табл ^ коэффициент значим; если (¿¿)расч < ¿табл ^ коэффициент незначим.
Форма уравнений регрессии с учетом значимости коэффициентов регрессии показана ниже.
Уравнения регрессии с учетом значимости для Р2, М2. Р2 у = 15,85 - 0,9Ц - 0,74х2, (12)
М2 у = 9,12 - 0,73х1 +1,90 х2. (13)
Расчет дисперсной адекватности
Проводим проверку адекватности математической модели про критерий Фишера.
РоговВ.А., Велис А.К., КопыловВ.В. Статистическое исследование силовых характеристик... Остаточная дисперсия, или дисперсия адекватности:
-> и ч и N о
=-М-=-М-у 2, (14)
" N - И - 1 п N - И -1 и
ад
где И — число параллельных опытов в и-й строке матрицы; N — число различных опытов (число строк матрицы); Б^ — средняя квадратическая ошибка прогноза
( N N \
Б2 = у Б2 = у (у - у )2
п ¿ш^ ип ¿ш^ и у и' и=1 и=1
Из (12), (13) у = 15,85 - 0,9Ц - 0,74Х2,
~ = Хг - Хг.0
Х' = Дх ,
где Х{ — кодированное значение г-го фактора; х1 — натуральное значение г-го фактора;
хг 0 — натуральное значение основного (нулевого) уровня г-го фактора; Дхг- — интервал варьирования г-го фактора; г — номер фактора.
Х1 - Х1,0 , л п А Х2 - Х2,0
у = 15,85 - 0,91—-10 + 0,74
Дх1 Дх2
так как х10 = 1000; а х20 = 22,5; Дх1 = 600; Дх2 = 17,5.
Средние значения экспериментальных данных и результаты для расчетных значений, полученных по уравнениям (12) и (14) для Р2, И2, приведены в табл. 11.
Таблица 11
Средние и расчетные значения опытов, по формулам (12) и (13)
Опыт р2 м
Уи Уи Я2 Уи Уи Я2
1 14,22 14.21 0,0001 10,54 10,29 0,06
2 16,00 16.02 0,0004 11,49 11,74 0,06
3 15,67 15.69 0,0004 6,25 6,50 0,06
4 17,51 17.50 0,0001 8,20 7,95 0,06
У (Уи - Уи )2 и=1 0,0010 N У (Уи- и=1 -Уи )2 0,24
С помощью значений табл. 11 и формулы (14) рассчитаем дисперсию адекватности табл. 12.
Таблица 12
Дисперсия адекватности уравнений регрессии для Рг, Мг
Pz М2
С2 йад С2 йад
0,002 0,48
Проверка адекватности модели
Ерасч < Етабл; при (д; /1; /2) — условие адекватности модели;
Ерасч — расчетное значение критерия Фишера;
Етабл — табличное значение критерия Фишера.
Из табл. 13 и табл. 8 получим Ерасч и значения занесем в табл. 13.
Расчетное и табличное значения критерия Фишера уравнений регрессии для Рг и Мг
Таблица 13
Pz М2
^асч ^табл ^асч ^табл
0,0077 7,71 1,14 7,71
Е_
= = 0002 = 0 0077
расч 02 « ~ ' ''
S2
0,26
Значения Е берем из таблицы (критерий Фишера) и запишем в табл. 13. /1 = N - М - 1 = 1; /2 = N (М - 1) = 4; Етабл = 7,71, Ерасч < Етабл ^ модель адекватна.
Так как Ерасч < Етабл, из табл. 13 следует, что математические модели для Р2 и М2 адекватны.
Х - Х " Х2 - Х2,0
у = 15,85 - 0,91—-10 + 0,74
Ах1
где у = Р ; х, = п; х2 = S; у = М ; х = п; х2 = S.
Ах„
После математических преобразований уравнения регрессии приобретут следующий вид:
Ря = 15,85 - 0,9^-П—1000 + 0,74£-22,5, 2 600 17,5
Р = 16,42 - 0,00152 • п - 0,0423 • £.
Р2
м.2
Уравнения регрессии для Рг, М2
Р2 = 16,42 - 0,00152 • п - 0,0423 • £, (15)
М2 = 7,89 - 0,00122 • п + 0,109 • £. (16)
Значения осевых усилий, возникающих на сверле, показаны на рис. 3, 4. Как ожидалось, максимальные усилия Р2, возникающие на сверле, получаются при минимальных частотах вращения п и минимальной осевой подачи £, в то время как минимальные значения Р2 — при максимальных значениях п (рис. 3). Линии одного цвета представляют уровни равного усилия Р2.
Рис. 3. Поверхность отклика осевого усилия Р2 по уравнению (15)
Из анализа поверхности оклика (см. рис. 3) видно, что максимальный крутящий момент, возникший при обработке синтеграна сверлением, получен при минимальных значениях вращения шпинделя и максимальных осевых подачах, в то время как минимальные значения крутящего момента возникают при максимальных значениях вращения шпинделя и минимальных значениях осевой подачи. Линии одного тона представляют уровни равного крутящего момента Ы2.
Рис. 4. Поверхность отклика крутящего момента Мг по уравнению (16)
Из вышесказанного следует, что оптимальным значениям для п и $ при обработке синтеграна (при осевой подаче 25—30 мм/мин и значениях числа оборотов 1300—1600 об/мин) соответствует светло-серый район на рис. 3 и 4.
После обработки экспериментальных данных можно сделать следующие выводы:
— минимальные осевое усилие получается при максимальных значениях n и S, в то время как крутящий момент минимальный, когда частота вращения максимальна и осевая подача минимальна, но этот режим не производителен;
— оптимальный режим, обеспечивающий станку нормальную производительность работы, выбран для n — 25—30 мм/мин и S — 1300—1600 об/мин;
— при сверлении отверстий в синтегране с оптимальными режимами обработки можно получить 5—6 класс шероховатости;
— из проведенных опытов рекомендуем обрабатывать синтегран сверлами с напаиваемыми твердосплавными пластинами с износостойким покрытием (TiAl)N, с дополнительным специальным покрытием сплавом NbHfTi, так как время, затраченное на заточку, уменьшает производительность.
ЛИТЕРАТУРА
[1] Рогов В.А., Кокарев В.И., Велис А.К. Исследование процесса обработки синтеграна твердосплавными сверлами. Т. 1. — М.: Фундаментальные исследования, 2013. [Rogov V.A., Kokarev V.I., Veliz A.C. Issledovanie protsessa obrabotki sintegrana tverdosplavnyimi sverlami. T. 1. — M.: Fundamentalnyie Issledovaniya, 2013.]
[2] Кокарев В.И. Применение статистических методов планирования эксперимента при идентификации процесса эксплуатации режущего инструмента. — М.: Крук, 2013. [Kokarev V.I. Primenenie statisticheskih metodov planirovaniya eksperimenta pri identifikatsii protsessa ekspluatatsii rezhuschego instrumenta. — Moskva: Kruk, 2013.]
[3] Адлер Ю.П.М.Е.В.Г.Ю.В. Планирование эксперимента при поиске оптимальных условий. — М.: Наука, 1976. [Adler Yu.P.M.E.V.G.Yu.V. Planirovanie eksperimenta pri poiske optimalnyih usloviy. — M.: Nauka, 1976.]
[4] Douglas C.M. Design and analisys of experiments. — Arizona: John Wiley & Sons Inc., 2001.
[5] Ашихмин В.Н., Гитман М.Б., Келлер И.Э. Введение в математическое моделирование. — М.: ЛОГОС, 2005. [Ashihmin V.N., Gitman M.B., Keller I.E. Vvedenie v matematicheskoe modelirovanie. — M.: LOGOS, 2005.]
[6] Тюкпиеков В.Н. Повышение эффективности обработки синтеграна на основе физ. и мат. моделирования: Диссертация. — М., 2002. [Tyukpiekov V.N. Povyishenie effektivnosti obrabotki sintegrana na osnove fiz. i mat. modelirovaniya: Dissertatsiya. — M., 2002.]
[7] Валид Махмуд Мохамед Э. С.Ш. Мат. и физ. моделирование динамики процесса резания композиционных структурно-неоднородных материалов (на примере синтеграна): Диссертация. — М., 2005. [Valid Mahmud Mohamed E.S.Sh. Mat. i fiz. modelirovanie dinamiki protsessa rezaniya kompozitsionnyih strukturno-neodnorodnyih materialov (na primere sin-tegrana): Dissertatsiya. — M., 2005.
[8] Кокарев В.И., Федоров В.И., Велис А.C. Исследование процесса сверления синтеграна твердосплавными сверлами c комплексным модифицированием режущей части // Вести СТАНКИНА. — 2013. [Kokarev V.I., Fedorov V.I., Veliz A.C. Issledovanie protsessa sverleniya sintegrana tverdosplavnyimi sverlami c kompleksnyim modifitsirovaniem rezhuschey chasti // Vesti STANKINA. — 2013.
[9] Елин А.В. Повышение эффективности и качества обработки полимербетонов шлифованием на примере синтеграна: Диссертация. —М., 2008. [Elin A.V. Povyishenie effektivnosti i kachestva obrabotki polimerbetonov shlifovaniem na primere sintegrana: Dissertatsiya. — M., 2008.
STATISTICAL STUDY ON THE MECHANICAL CHARACTERISTICS WHEN DRILLING SINTETHYC GRANITE
V.A. Rogov, A.C. Veliz, V.V. Kopylov
Department of Mechanical Engineering, Metal-cutting and Tool Engineering Faculty People's Friendship University of Russia Ordzhonikidze str., 3, Moscow Russia, 115419
This paper is a research on sintegran to determine the effects of drilling on the power characteristics of machine: axial force — Pz, and torque — Mz to minimize the stresses on the drilling process. Investigations were carried out using spiral carbide drills alloy BK8, the surface of which was treated with an electron beam by magnetron sputtering alloy NbHfTi gradient for the content of alloying elements intermediate layer before hardfacing (TiAl) N. Studies were conducted on a machining center equipped with CNC (computer numerical control), with the possibility of diagnosing power machine performance in real time. Studies have shown the optimal values of spindle speed and axial feed for machining sintegran.
Key words: Sintegran, drilling of sintegran, carbide drills, surface modification, computer diagnostics of mechanical tests, wear resistant coatings, axial feed, spindle speed.
