УДК 621.914
АНАЛИТИЧЕСКОЕ МОДЕЛИРОВАНИЕ СИЛОВОГО РЕЖИМА ВЫСОКОСКОРОСТНОГО ФРЕЗЕРОВАНИЯ МАТЕРИАЛА С ВЯЗКОПЛАСТИЧЕСКИМ УПРОЧНЕНИЕМ
© 2012 А. И. Хаймович, А. В. Кузнецов
Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)
В настоящей статье рассматривается аналитическое моделирование силового режима высокоскоростной обработки (ВСО) фрезерованием жаропрочных сплавов с учетом их вязкопластического упрочнения. В модели применено экспериментально полученное определяющее соотношение Джонсона-Кука, характеризующее динамическое упрочнение и температурное разупрочнение обрабатываемого материала при ВСО.
Высокоскоростная обработка, жаропрочные сплавы, реологические свойства, напряженно-деформированное состояние, модель Джонсона-Кука.
В современных авиационных ГТД используются детали, изготавливаемые из вязких жаропрочных сплавов, обработка которых вызывает ряд трудностей.
Целью настоящего исследования является аналитическое моделирование силового режима высокоскоростной обработки фрезерованием жаропрочных сплавов с учетом их вязкопластического упрочнения.
Мощность резания в очаге пластической деформации (ОПД) объемом W, локализованной у режущей кромки, определяется зависимостью Р = (е2,е2,в)е2^ , (1)
где
(£2А,в)
эквивалентное
напряжение, е2 - интенсивность скоростей
деформаций. В свою очередь
е2 = V),е2 = р2Ж .
(2)
моделью пластического течения Джонсона-Кука [1,2]:
а, = [а0 + аее2Па 111 + Ь 1пе
1 -в'
(3)
где
е2
'21а,=а,.
относительная
интенсивность деформации, е2 = —
"2а,=а,.
относительная
деформации,
интенсивность в в
скорости
в =■
в
в
тв -
где У0 - скорость резания, у1 - виртуальное кинематически допустимое поле скоростей (КДПС) в объеме W, е2 - интенсивность
деформаций, в -температура у режущей кромки. Для того, чтобы определить функциональную связь между
напряжениями с одной стороны и е2, е2, в с
другой стороны в форме а, = а,(е2,е2,в) , воспользуемся феноменологической
относительная температура, па , Ь эмпирически определяемые коэффициенты, характеризующие свойства материала при высоких скоростях течения и температурах деформирования, а, 0 = а0 + ае - предел текучести при в = 1,е2 = 1,е2 = 1 и скорости деформирования V 0| = .
0|а, а, о
При высоких скоростях резания, в силу локализации ОПД в малом объеме, целесообразно полагать
аее2 а = а
(4)
Модель Джонсона-Кука
характеризует температурное разупрочнение деформируемого материала (коэффициент тв ), возникающее в ОПД при высоких скоростях резания, одновременно связкопластичным (коэффициент Ь ) и
в
п
V
01ст. =а
V
о
деформационным (коэффициент ае ) упрочнением.
Оптимальное сочетание скорости резания, глубины и ширины резания обеспечивает минимальное значение мощности резания (1) за счет разупрочнения обрабатываемого материала при резании с высокими скоростями. Таким образом, экспериментально определяемые
постоянные коэффициенты, входящие в модель (3,4) и аналитические зависимости, характеризующие интенсивность скорости деформации в ОПД при резании, позволяют найти оптимальные режимы
высокоскоростной обработки лезвийным инструментом согласно зависимости:
ГР Л
(5)
V st у где ОПД,
{у,а,ги }-іпуа
-»штп
V - действительное поле скоростей в
шщ
0
глобальный
минимум,
(V,,,, г, ^ } - технологические изменяемые оптимизируемые параметры, (у, а, ги} -
фиксированные параметры, определяющие наряду с (,, г} геометрию ОПД, у, а, ги -соответственно передний, задний углы и радиус при вершине режущего клина, ,, г -относительные ширина и глубина резания. Для фрезерования ,, г - две ортогональные величины, измеряемые в плоскости, перпендикулярной вектору скорости резания, ^ - подача на зуб.
Произведение относительных
величин - ширины s и глубины резания t -определяют производительность лезвийной обработки. Поэтому согласно (5) минимум
:э
минимальной диссипации мощности резания, затрачиваемую на пластическую деформацию при максимальной
производительности.
Поскольку действительная картина течения материала в ОПД априори не известна, то для практических целей поиска глобального минимума min0 достаточно
определить значения а, = а, (е2,е2,в) в
характерных областях ОПД. Докажем последнее утверждение.
Условно разделим ОПД на N
характерных областей, геометрия которых определяется положением условной
плоскости сдвига, задаваемой
коэффициентом усадки стружки Л6, и
геометрией режущего клина. В этом случае (1) запишется как
2 я Є2 і ,в] )Є2 № . (6)
р=Ё|(є
і=1 Wj
В каждой из этих областей зададим vi - виртуальное КДПС течения металла при
резании. Согласно экстремальным принципам механики сплошной среды, если КДПС vi приближается к действительному
полю скоростей, то пластический потенциал Р стремится к своему локальному минимуму шіпх при фиксированных
значениях v0,в . В этом случае из выражения для мощности резания (6) в соответствии с экстремальными принципами механики сплошных сред и по теореме о наибольшем значении функции имеем:
:£ [СТ, (є2тах і )1Р,
Р '£
----------> ШІП, ,
vr ^v 1 '
і=1
Рі = Є
2шах і
(8)
где є 2
(vi) = шах{є2і(vi)} - максимальное
' 2шах і V і )
значение интенсивности скорости деформации на j -участке для КДПС vi ,
'2е^
(Vi) = ■
Wi
средняя
по
функционала
соответствует локальному объему Wj интенсивность
скорости деформаций на j -участке для КДПС V,,
р. имеет физический смысл относительной мощности резания при а, = 1 , которая
W■
рассеивается в единичном объеме — для j -
ой области ОПД.
Таким образом, для того, чтобы иметь верхнюю оценку относительной мощности
t-соті
резания
Р
достаточно в соответствии с
зависимостью (8) определить значения эквивалентных напряжений (є2шахі ^г.), ві)
в характерных областях ОПД, а также вычислить значения средних
интенсивностей скоростей деформаций є2 (уі) в этих областях.
Предположим, что значения р. и
є2шах і для условий ортогонального резания
можно определить, если ОПД разделить на 4 характерные области.
Анализ этих областей ОПД позволяет сделать вывод, что их геометрия хорошо описывается в цилиндрических координатах.
На рисунке 1 представлены характерные области ОПД с сегментными границами, описанными в 2-х
противонаправленных цилиндрических
системах координат с разнесенными началами.
Введем безразмерные величины, характеризующие степень деформации в областях 1, 2, 3 (рис. 2):
Л = = —,Л = ^, (9)
Г Г Г
Г1 2 ГЪ
где 11 = Л -1 - степень деформации в 1
области и коэффициент л = — - усадка
Я
(уширение) стружки.
ортогональном резании с неотрицательным углом у
Область 1 и 2 рисунка 1
соответствуют области первичной интенсивной деформации в зоне резания, область 3 (область стружки) - зона вторичной деформации, образованной трением по передней грани режущего клина. Области 1,2,3 являются пластичными,
область 4 - жесткой. Интенсивные
деформации сдвига, вызывающие
относительные смещения областей из-за разрыва скоростей движения материала, наблюдаются по границам областей 1-4, 2-4, и 2-3.
В результате получаем систему из 15 уравнений, которые связывают между собой 14 размеров ОПД и 4 относительных параметра ЛЛЛЛ , характеризующие силовой режим резания.
Чтобы уровнять количество
уравнений и неизвестных, выберем из числа рассмотренных 3 независимых размера, которые будем определять методами теории пластического деформирования. В качестве такой тройки параметров целесообразно использовать множества |Д1, Л,, Щ или
ЛЛЛ} . Действительно, Л1 определяет деформацию в заготовке, Л3 - деформацию в стружке, Л, - соотношения деформаций
между областямив заготовке и стружке.
В итоге диссипация мощности пластической деформации и максимальная интенсивность скоростей деформаций для рассмотренных областей определяется
зависимостями:
-< 2^ . Л
3 л-1
Л
>/2
3
'2шах1
Р2 < 2
Л -1
=< 2^
3 Л-1
л/2
>2шах1
Г1
р3 < 0,62v(
3 л-1
а
КАЯ - й)1п К є
2шах2
= є
Л ссеу
Л
1—
л
2,633
В
2шах1 ’
(10)
є ? = к •
2шах3
Здесь
cсs уЯ3
-В_
1
v
о
B = nbB + (1 -nb )B ■ ,
b max V b S mm ’
в =
• 2 Q3
ЛП ^3
= [1 + A(A ^3)max (A(P, ^3)max + 2)sm ^ ] Bmin = [1 + A(P,^3)mm(A(P,^3)m,n + 2)s1^ ^ Й
A(P, Ф)max =
С08 уЯ3
n(4 -1)k
1
Л
- (24 -1) • —
R3 j
24 - (24-1) • R-^4 -
V r1 R3 r1
1 - k2
• 0,, sin^3/
A(P,^)min =
С08 у^3
• Q3/
sin ад/
- + -
n(4 - 1)k
1zkL k2 ’
R
+1 - 4
R
1 v R R
P1-4 =;/, —-1[ rL(R> - R2) - R'!n rl]cos n01,
P2-4 =
'2max1-4
R
yfin 4 -1
cos nQ1 - cos nQ
2R
V R2
' 2 max 1: 1 - 1
n
t
+ sin nQ1 - sin nQ I,
= t
2max2-4 2max2^
A,:
1 + -
'2max1,2-HH
V3n2 V 4 -1 :-HH = max(t
(1 - cos nQ),
2 max1 2 max2
Pi-ин
2v0 R3
k +1
cosy
■ С0Б
Qi
+
P2-3 =
1 24-1 .
+ —:-------TsinQ ^
n 4-1 J
t2max3-HH t2max3 '
л/34 (k + 1)cosy
С0Б
Q
3/ —
J
Vn
• cos nQ(R1 ln 42 - R2 + r1),
л/3л, (Л-1)
е2тах2-3 = тах(е2тах 2 , 82тах3}.
В зависимостях (10) значения
n = 0,612, k = 2.
2
= 1,634
Температура резания оказывается достаточно сложной функцией большого количества технологических параметров, характеризующих свойства обрабатываемого и инструментального материалов, геометрии инструмента, режимов резания
Для оценки температуры резания был применен аналитический метод, изложенный в работах [3,4].
В качестве объекта
экспериментальных исследований был выбран сплав ХН45МВТЮБР-ИД, широко используемый для производства
компрессорных лопаток ГТД.
Полуэмпирическая зависимость
эквивалентного напряжения (интенсивности напряжений) от скоростного и
температурного режима фрезерования была получена авторами на основе обработки экспериментальных данных и
формализована в виде закона, подобного модели Джонсона-Кука:
о. =600-10е
0,85(1 + 1,3t20,86)+ 0,061ln
V - vn
Vmax - V0
1 -1,347
6-6
rn
6 - 450
(11)
V,™,, = 600 м / мин, V, = 50 м / мин,6„п = 20000 С.
^таХ ^ 0 ^ ПЛ
На рисунке 2 представлены
результаты аналитического моделирования силового режима фрезерования по
изложенной выше методике в форме
поверхности относительной удельной мощности фрезерования р = р(у, /). На рисунке 3 изображена поверхность относительного усилия при вершине радиуса вдоль режущей кромки зуба фрезы в форме
зависимости / = —г— = / (V / ) , где Р -
Аа, 2
усилие на режущей кромке вдоль контактной поверхности Аг.
Анализ геометрической формы поверхности р = р^, /) показывает, что
увеличение скорости резания при
фиксированной подаче на зуб приводит к увеличению усилия резания по логарифмическому закону, что связано с повышенными вязкопластическими
свойствами обрабатываемого материала.
1
2
k
1
2
1
V
0
V
0
V
0
При постоянной скорости резания в
интервале подач в области / = 0,04...0,08
мм/зуб, процесс резания характеризуется локальным минимумом по мощности резания (рисунок 2), что объясняется оптимальной формой ОПД в зоне резания.
Относительное удельное усилие резания
Расчетная скорость резания,
Подача на зуб, мм/зуб
обрабатывающем центре МікгопЦСР800 Durс. Технологические параметры обработки: боковой съем £ = 0,1мм, съем по высоте і = 4 мм, подача на оборот фрезы /об =
0,32 мм/об, фрезерование попутное, обработка с СОЖ.
Относительное удельное усилие резания
"Ж. - экспериментальные значения.
-■п- - аналитич. значения
резания в зависимости скорости резания и величины подачи на зуб р = р(у, / ) при ВСО сплава ХН45МВТЮБР-ИД
При скоростях резания выше
300м/мин контактные напряжения на
радиусе режущей кромки (рисунок 3) составляют более 0,7а £ , что существенно снижает стойкость режущего инструмента.
, -,
0,95
0,90
0,85
0,80
0,75
0,70
0,65
0,60
0,55
Рис.З.Расчетное относительное удельное усилие на радиусе режущей кромки фрезы в зависимости скорости резания и величины подачи на зуб при ВСО сплава ХН45МВТЮБР-ИД
С целью проверки состоятельности аналитической модели были проведены экспериментальные исследования по измерению усилия резания при высокоскоростном фрезеровании образцов из сплава ХН45МВТЮБР-ИД.
Сила резания измерялась с помощью 3-х осевого динамометрического стола Kistler9257B с регистрацией составляющих усилия резания Fx,Fy,Fz в реальном времени. Обработку образцов производили концевой цилиндрической фрезой диаметром 12мм на
200 220 240 260 280 300 320 340 360 380 400 420 440 460
скорость резания м/с
Рис. 4. - Зависимость относительной силы резания от скорости резания при высокоскоростном фрезеровании образца из ХН45МВТЮБР-ИД (ЭП718 ИД). Подача 0,32 мм на оборот, фреза диаметром 12 мм, боковой съем 0,1 мм, съем по высоте 4 мм
На рисунке 4 приведены теоретические и экспериментальные кривые зависимости усилия резания от скорости резания. Экспериментальная кривая получена аппроксимацией опытных значений методом минимизации
среднеквадратичных отклонений.
Теоретическая кривая является сечением поверхности р = р(у, /) плоскостью
/ = 0,08мм / зуб .
Анализ приведенных данных свидетельствует, что представленная
аналитическая модель процесса
фрезерования вязкопластического материала является верхнеграничной оценкой силового режима фрезерования. Точность
аналитической модели в большой степени определяется точностью определяющих соотношений а, = а, (е2,е2,в) ,
устанавливающих связь между
инвариантами напряженного и
деформированного состояния для
конкретного обрабатываемого материала в определенных условиях резания.
Библиографический список
1. Johnson, G.R. Cook, W.H. A constitutive model and data for metals subjected to larges trains, high strainrates and high temperatures. ^.neKTpoHHbiHpecypc]: /Johnson G.R. ,Cook W.H.//Proceedings of the 7th International Symposium on Ballistics.-1983.- p 541-547/-
http://www.lajss.org/HistoricalArticles/A%20co
nstitutive%20model%20and%20data%20for%2
0metals.pdf.
2. Steinberg, D.J. Cochran, S.G.Guinan,
M.W. A constitutive model for metals applicable at high-strain
rate. ^.neKTpoHHbiHpecypc] :/Steinberg, D.J.
Cochran, S.G. Guinan, M.W.
//JournalofAppliedPhysics. -51: 1498,1980.http://link.aip.org/link/?JAPIAU/ 51/1498/1
3. Силин, С.С. Теоретическое
определение параметров процесса резания.[Текст]: /Силин С.С.
//Производительная обработка и
технологическая надежность деталей машин.Межвузовский сб. научных трудов.-Ярославль: ЯПИ, № 6, 1977. - С. 3—16.
4. Силин, С. С. Метод подобия при резании материалов [Текст] / С. С. Силин. М.: Машиностроение, 1979. - 152 с.
ANALYSIS ОF THE CUTTING FORCE OF HEAT RESISTING ALLOYS DURING
HIGH-SPEED MILLING
© 2012 A. I. Khaimovich, A. V. Kuznetsov
Samara State Aerospace University named after Academician S.P. Korolyov (National Research University)
In this article the analytical model of high speed milling (HSM) of the heat resisting alloys taking into account their hardening properties is considered. Johnson-Cook's experimentally received constitutive model, characterizing dynamic hardening of worked material is applied.
High speed milling, heat resisting alloys, rheological properties, Johnson-Cook constitutive model.
Информация об авторах Хаймович Александр Исаакович, кандидат технических наук, доцент кафедры производства двигателей летательных аппаратов, Самарский государственный
аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Е-таП: ЬегШ [email protected]. Область научных интересов: исследование поликристаллических материалов, подвергнутых интенсивной деформации.
Кузнецов Антон Владимирович, аспирант, Самарский государственный
аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Е-та11: а.у [email protected]. Область научных интересов: программирование многокоординатного оборудования с ЧПУ, оптимизация параметров лезвийной обработки резанием.
Khaimovich A.I., Candidate of Technical Sciences, associate professor at production of aircraft engines department, Samara State Aerospace University named after Academician S.P. Korolyov (National Research University). E-mail: berill [email protected]. Area of research: research of the polycrystalline materials subjected to intensive deformation.
Kuznetsov A.V., graduate student, Samara State Aerospace University named after Academician S.P. Korolyov (National Research University). E-mail: a.v [email protected]. Area of research: programming NC multi-axis equipment, optimization of cutting processing parameters.