УДК 621.91
ИССЛЕДОВАНИЕ СИЛ РЕЗАНИЯ, ВОЗНИКАЮЩИХ ПРИ ОРТОГОНАЛЬНОМ СВОБОДНОМ ТОЧЕНИИ СПЛАВА Д16, С ПОМОЩЬЮ МЕТОДИК АНАЛИТИЧЕСКОГО
И КОНЕЧНО-ЭЛЕМЕНТНОГО АНАЛИЗА
© 2017 Д.В. Евдокимов, М.А. Олейник
Самарский национальный исследовательский университет имени академика С.П. Королёва
Статья поступила в редакцию 31.08.2017
В статье приводится сравнение результатов вычисления составляющих силы резания при ортогональном точении алюминиевого сплава Д16. Исследования построены на основе методики аналитического определения силы резания при точении, а также на основе созданной компьютерной модели базирующейся на методе конечных элементов. Задействованная аналитическая методика определения силы резания при точении была дополнена рекомендациями А.Н. Резникова. Ключевые слова: точение, ортогональное точение, сила резания, составляющие силы резания, феноменологическая модель Джонсона-Кука, плоскость сдвига, касательные напряжения, срезаемый слой, стружка.
ВВЕДЕНИЕ
Процесс механической обработки материалов лезвийным инструментом сопровождается противодействующими силами, знание которых необходимо для выполнения расчетов на прочность, жесткость и виброустойчивость элементов конструкции станков и режущих инструментов.
Кроме того существует подход к оценке действующей температуры в зоне резания при точении, основывающийся на знании величин плотностей тепловых потоков, определить которые возможно при условии известных величин составляющих силы резания. Как известно, на формирование состояния поверхностного слоя деталей оказывают влияние силовой и температурный факторы. Причем, если температура в зоне резания достигает значений, соответствующих значениям критических точек аллотропических превращений в поверхностном слое материале заготовки, то влияние температуры на формирование поверхностного слоя существенно возрастает.
Остаточные напряжения, формируемые в поверхностном слое заготовок при резании и являющиеся одним из основных параметров, определяющих качество обработки, возникают в результате воздействия силового и температурных полей в зоне резания [1, 2, 3]. Также известно, что остаточные напряжения в значительной степени влияют на сопротивление усталости деталей. Существуют методики оценки остаточных напряжений в поверхностном слое заготовок наведенных после их
Евдокимов Дмитрий Викторович, аспирант. E-mail: dmitry.evd.ssau@gmail.com Олейник Максим Андреевич, студент. E-mail: oleynik1997@mail.ru
обработки. Такие методики обычно требуют знание величины силы резания и температуры, действовавшей в зоне резания [4].
Поэтому необходимо проводить исследования силы резания, применять различные методики для ее оценки, находить их недостатки и искать пути совершенствования существующих методик.
МЕТОДИКА АНАЛИТИЧЕСКОГО ОПРЕДЕЛЕНИЯ СИЛЫ РЕЗАНИЯ
При создании методики авторы работы [5] использовали схему распределения сил в зоне резания, которая была составлена Юджином Мёрчантом (M. Eugene Merchant). Для исследования, представленного в данной работе, в схему распределения сил были внесены изменения в виде добавленных составляющих силы резания, таких как N3 и F3 [6, 7]. Необходимость их присутствия объясняется наличием у неидеального режущего клина площадки его контакта по задней поверхности с заготовкой. Распределение сил показано на рис. 1.
В соответствии с методикой, изложенной в работе [5], касательная сила (сила сдвига), возникающая в плоскости сдвига АВ (рис. 1), определяется по формуле:
Рх = AB -b- тсд,
где b - ширина срезаемого слоя, Тсд. - касательные напряжения, возникающие в зоне сдвига, АВ - длина плоскости сдвига, которую просто вычислить, зная глубину резания t и угол сдвига ß;:
t
Угол сдвига можно определить, используя зависимость А.Н. Резникова [7]:
Машиностроение и машиноведение
Рис. 1. Схема распределения сил в зоне резания
sin ^
cosy
Jk2 - 2/csiny+l'
(1)
Мёрчанта [5], что стало возможным благодаря известному углу сдвига (1):
где k - коэффициент усадки стружки.
Определив силу сдвига Р , по зависимости, представленной в работе [8] возможно вычислить величину нормальной силы Pn , перпендикулярной к Р%. Равнодействующей этих двух сил является сила R, а её реакцией, направленной в режущий клин - сила R1 (R = R1). Зная силу R1 можно определить главную, тангенциальную Pz и радиальную Ру составляющую силы резания, а также нормальную силу Nn и силу трения Fn , действующие на передней поверхности режущего клина.
В методике, изложенной в работе [5] определение тангенциальной Pz и радиальной силы резания Ру сведено к следующим формулам:
р _ fi-со5(рц - у) р _ fi-sin(p^ - У)
Z cos (рх + р^ - у) ' У cos (рх + р^ - у) '
где р^ - угол действия силы давления [9].
Но как уже было сказано выше, к данным формулам были добавлены слагаемые в виде силы трения задней поверхности режущего клина о заготовку F3 и силы упруго-пластического взаимодействия задней поверхности режущего клина с заготовкой N3. Оценить величину данных сил помогли рекомендации А.Н. Резникова [7], согласно которым сила F3 принимается примерно равной силе N3 и может быть определена по зависимости:
F3 « 0,252 ■ св • b • 13,
где aB - предел прочности обрабатываемого материала, b - ширина срезаемого слоя, l3 - длина участка контакта задней поверхности режущего клина с заготовкой.
При проведении исследования угол р^ был получен с использованием формулы Юджина
П
Р
В =---* + -.
И1 А 11
У
Далее, для определения всех представленных на рис. 1 сил, необходимо рассчитать значение силы сдвига Р%. Ее определение сводится к умножению величины касательных напряжений Тсд. в плоскости сдвига АВ на площадь плоскости сдвига. Для того чтобы определить величину касательных напряжений в плоскости сдвига АВ авторами работы [5, 8] было использовано уравнение Джонсона-Кука:
1
VI
хсд
А + В
1 + С1п —
e0J
[1 - Тт],
0 при Т <ТС
Т =
Т-Тг
при ТС<Т< Гп,
1 при Т > Гп,
где 8СД - деформация сдвига; 8СД - скорость деформации сдвига; £0 - справочная скорость пластической деформации (обычно, перед тем как определять коэффициенты зависимости, данное значение принимается равным 1 с-1), что упрощает вычисление напряжений пластического течения по уравнению; Г - справочная температура (принимается обычно равной комнатной температуре); Тл - температура плавления обрабатываемого материала; А - предел текучести материала, соответствующий справочным скоростям деформации и температуры; В - модуль деформационного упрочнения; п -показатель степени деформационного упрочнения; т - экспонента термопластикации; С - показатель чувствительности к скорости деформирования.
РАЗРАБОТКА КОМПЬЮТЕРНОЙ МОДЕЛИ ДЛЯ ПРОВЕДЕНИЯ ЧИСЛЕННЫХ ЭКСПЕРИМЕНТОВ
В универсальной программной системе конечно-элементного анализа ЛЫБУБ была разработана компьютерная модель ортогонального свободного точения алюминиевого сплава Д16, применительно к установившемуся процессу резания. Модель представлена в виде двух виртуальных тел заготовки и режущего клина, для каждого из которых задано свое кинематическое описание движения и свое описание материала внутри компьютерной программы.
ОПИСАНИЕ ГЕОМЕТРИИ ВИРТУАЛЬНЫХ ТЕЛ КОМПЬЮТЕРНОЙ ПРОГРАММЫ
Заготовка имеет следующие геометрические параметры: высота 5 мм, ширина 2 мм, длина 10 мм. Режущий клин имеет передний угол у = 12°, ширина режущего клина 4 мм, а радиус скругле-ния режущей кромки составляет 0,025 мм. Две другие кромки виртуального тела режущего клина имеют скругление 0,25 мм. Это связано с задачей корректного наложения сетки конечных элементов в данных областях режущего клина (рис. 2).
ОПИСАНИЕ МАТЕРИАЛА ВИРТУАЛЬНЫХ ТЕЛ
Помимо геометрических данных заготовки и режущего клина, в роли исходных данных для вычислений выступили также описания материалов заготовки и режущего клина при использовании модели пластичности Джонсона-Кука и модели разрушения Джонсона-Кука, уравнение которой приведено ниже:
£/ =
1 + ДЛп( —
(1 + о5т),
где стт является средней величиной по трем главным напряжениям, О - эквивалентные напряжения, полученные по критерию пластично-
сти фон Мизеса, Е^ - эквивалентная деформация при разрушении, а Di - это коэффициенты, получаемые по результату экспериментов.
КИНЕМАТИЧЕСКОЕ ОПИСАНИЕ ДВИЖЕНИЯ ВИРТУАЛЬНЫХ ТЕЛ
Выполненные расчеты, что при задействовании аналитической методики определения силы резания, что при использовании компьютерной модели, были проведены для режимов обработки сплава Д16, соответствующих подаче 0,47 мм/об и 1 мм/об при глубине резания составляющей 2 мм и скорости резания менявшейся в диапазоне от 15 до 75 м/мин.
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ
Прежде всего, стоит отметить, что на величину силы резания при использовании аналитической методики, оказывает большее влияние изменение коэффициента усадки стружки, нежели изменение скорости резания. Изменение скорости резания влечет к изменению касательных напряжений, в плоскости сдвига. В свою очередь, это изменение, при всех исходных данных представленного исследования оказалось небольшим и для увеличения точности прогнозируемой величины силы резания необходимо знание изменения усадки стружки в зависимости от скорости и глубины резания. В представленной работе коэффициент усадки стружки был принят постоянным для всех скоростей резания, т.е. не учитывалось его изменение в зависимости от скорости резания.
Ниже представлены графики изменения составляющих силы резания и их результирующей, полученные посредством компьютерной модели (рис. 3 и 4).
Далее приведена таблица со средними значениями сил, вычисленными с помощью аналитической методики и сравненными с аналогичными значениями, полученными с помощью компьютерной модели (табл. 1).
Рис. 2. Кинематическая схема компьютерной модели и изометрия виртуальных тел заготовки и режущего клина
Машиностроение и машиноведение
Рис. 3. Зависимость составляющих силы резания от скорости резания при подаче £ = 0,47 мм/об
Рис. 4. Зависимость составляющих силы резания от скорости резания при подаче £ = 1 мм/об
Таблица 1. Величины средних значений сил при ортогональном свободном точении сплава Д16, при глубине резания £ = 2 мм на диапазоне скоростей резания 15-75 м/мин
5 = 0,47 мм/об
сила, Н комп. модель мат. модель А, %
Ру 138 175 21
Рг 141 167 15
Я 198 242 18
Б = 1 мм/об
Ру 89 98 10
Рг 419 461 9
Я 428 471 9
ВЫВОДЫ
Были получены зависимости составляющих силы резания при ортогональном свободном точении сплава Д16 с помощью разработанной компьютерной модели.
Величины средних значений составляющих силы резания, полученные с помощью компьютерной модели, были сравнены с аналогичными значениями, полученными с помощью аналитической методики, использование которой при исследовании сопровождалось внедрением некоторых рекомендаций и зависимостей А.Н. Резникова. Разница полученных результатов лежит в диапазоне 9-21%.
СПИСОК ЛИТЕРАТУРЫ
1. Кравченко Б.А. Митряев К.Ф. Обработка и выносливость высокопрочных материалов. Куйбышев: Куйбышевское книжное издательство, 1968. 132 с.
2. Evdokimov D.V., Skuratov D.L., Fedorov D.G. Thermal Fields in the End Milling of VT6 Titanium Alloy with Cooling // Russian Engineering Research, 2014. V. 35. No. 10. P. 773-776.
3. EvdokimovD.V., Skuratov D.L., Fedorov D.G. Influence
of Tool Wear on the Heat Flux Distribution and Temperature at the Contact Surfaces in the End Milling of OT4 Titanium Alloy // Russian Engineering Research, 2016. V. 36, No. 4. P. 324-327.
4. Skuratov D.L., Evdokimov D.V., Fedorov D.G. Thermal Stress Resarch of Processing and Formation of Residual Stress When End Milling of a Workpiece // World Applied Sciences Journal, 2014. V. 31 (1). P. 51-55.
5. Binglin L., Xuelin W., Yujin H., Chenggang L. Analytical prediction of cutting forces in orthogonal cutting using unequal division shear-zone model // International Journal of Advanced Manufacturing Technology, 2011. V. 54 (5). P. 431-443.
6. Евдокимов Д.В., Скуратов Д.Л. Методика расчета составляющих силы резания при концевом фрезеровании на базе феноменологической модели Джонсона-Кука // СТИН. 2017. № 1. С. 29-33.
7. Резников А.Н. Теплофизика резания // М.: Машиностроение, 1969. 288 с.
8. Zhongtao, F., Wenyu Y., Xuelin W. and Jürgen L. Analytical modeling of milling forces for helical end milling based on a predictive machining theory // Procedia CIRP, 2015. V. 31. P. 258-263.
9. Ozlu E., Budak E., Molinari A. Analytical and experimental investigation of rake contact and friction behavior in metal cutting // International Journal of Machine Tools and Manufacture, 2009. V. 49 (11). P. 865-875.
RESEARCH OF CUTTING FORCES AT ORTHOGONAL FREE TURNING OF D16 ALLOY BY MEANS OF TECHNIQUES ANALYTICAL AND A TERMINATING ULTIMATE ANALYSIS
© 2017 D.V. Evdokimov, M.A. Oleynik
Samara National Research University named after Academician S.P. Korolyov
The article compares the results of calculating the cutting force components for the orthogonal turning of an aluminum alloy D16. The investigations are based on the method of analytical determination of the cutting force at turning, and also on the basis of the computer model based on the finite element method. The analytical method used to determine the cutting force during turning was supplemented by the recommendations of A.N. Reznikov.
Keywords: turning, orthogonal turning, cutting force, cutting force components, Johnson-Cook phenomenological model, shear plane, tangential stresses, shear layer, chips.
Dmitry Evdokimov, Postgraduate Student.
E-mail: dmitry.evd.ssau@gmail.com
Maxim Oleynik, Student. E-mail: oleynik1997@mail.ru