CC BY
69
Электронное научно-техническое издание
НАУКА и ОБРАЗОВАНИЕ
_Эл № ФС 77 - 30569. Государственная регистрация №0420900025. ISSN 1994-0408_
Исследование динамики процесса резания при обработке жаропрочных материалов
# 08, август 2009 автор: Яковлев М. Г.
Первые публикации
Программы и программные системы
Учебные программы
Студенческая весна
Общие проблемы инженерного образования
Инженер в современной России
Экобионика
Зарубежное образование
История технического прогресса
Будущий инженер
Вне рубрик
Расширеный поиск Подписаться на новости
ПОИСК
Ред.совет Специальности Рецензентам Авторам Архив
УДК 621.91.01
cntl@salut.ru
НИИД ФГУП ММПП "Салют"
Введение
Исследования показывают, что на период стойкости режущего инструмента в значительной степени влияют вибрации при резании [1]. Основные характеристики вибраций, частота и амплитуда, зависящие от жесткости, массы и демпфирования технологической системы резания, также являются функциями режимов резания и свойств обрабатываемого материала. На динамическое поведение инструмента при обработке резанием труднообрабатываемых материалов большое влияние оказывают возбуждаемые автоколебания инструмента. Существуют следующие основные гипотезы причин возбуждения автоколебаний инструмента:
1) Нелинейность характеристик сил резания [2]. При этом сила резания считается связанной с режимами обработки и геометрией инструмента аналитическим соотношением следующего вида:
F= Кр Sf* V™ ■ b
rKv
■м
где Б0 - подача, V - скорость резания, Ь - глубина резания, Кр,Ку,КБ,Ю - эмпирические коэффициенты.
F(t)=Kp*h(t-Tp)
где Tp - время запаздывания.
F(t)=Kp[h(t)+Kih2(t)+ K2h2(t)+ K3 h 2 (t)],
(1)
Сразу отметим, что такая статическая характеристика силы резания не применима для анализа процессов обработки материалов с вибрациями, так как коэффициенты в формуле (1) при вибрациях изменяются с изменением режимов обработки.
2) Падающая зависимость коэффициента трения от относительной скорости между стружкой и инструментом, а также между деталью и инструментом [3].
3) Регенерация колебаний при движении по следу, образующемуся на поверхности резания в процессе обработки [4].
В реальных технологических системах при обработке труднообрабатываемых материалов возможно возбуждение автоколебаний по нескольким причинам одновременно. Значит необходимо учесть все выше описанные причины, вызывающие вибрации инструмента.
силы резания (1) могут быть рассмотрены следующие динамические
Кроме статической характеристики характеристики:
1) Линейная динамическая характеристика, представляющая собой зависимость силы резания от толщины срезаемого слоя [4]:
Р(Ц=Кр*М(Ц (2)
где Кр - эмпирический коэффициент, - толщина срезаемого слоя в произвольный момент времени t.
2) Линейная динамическая характеристика с запаздыванием, когда сила резания в настоящий момент времени t пропорциональна толщине срезаемого слоя, измеренного в момент времени ^Тр [4]
(3)
3) Нелинейная динамическая характеристика, представляющая собой зависимость силы резания от толщины срезаемого слоя в следующем виде [6]:
(4)
где Ki - эмпирические коэффициенты.
Определение эмпирических коэффициентов и получение решения нелинейной динамической характеристики вида (4) является очень трудоемкой задачей. Для построения вибрационной модели поведения режущей кромки инструмента при обработке труднообрабатываемых сплавов будем использовать линейную динамическую характеристику (2). Характеристика по гипотезе регенеративного возбуждения колебаний при движении по следу от предыдущего прохода хорошо и быстро решается численными методами с помощью ЭВМ.
Экспериментальные данные для построения вибрационной модели поведения режущей кромки инструмента получены на диагностическом стенде. Стенд разработан на базе токарного станка фирмы Jesco Machinery 1650ENC с ЧПУ системы Fagor 800TGI и укомплектован трехкомпонентным динамометром «9257 ВА» со встроенным усилителем заряда, акселерометром «Kistler 8614A1000M1» с усилителем и фильтром сигнала «Kistler 5127В» и трех компонентным акселерометром "ГлобалТест АР2043-50" с высокочастотным фильтром сигнала. Стенд обеспечивает автоматизированный прием информации о динамических и вибрационных процессах при резании в реальном масштабе времени и обработку этой информации с целью определения физических критериев, характеризующих состояние технологической операции.
СОБЫТИЯ
Международная
научно-практическая
конференция
"Использование технологий дистанционного обучения в высшей школе "
III Международная
научно-практическая
конференция
"Информационная среда вуза XXI века "
НОВОСТНАЯ ЛЕНТА
12.10.2009
Всероссийская научная школа для участников программы «Участник Молодежного Научно-Инновационного Конкурса» (У.М.Н.И.К)
14.09.2009
Абитуриенты смогут поступать только в пять вузов - проект приказа
24.08.2009
Комсомольская правда > 24.08.2009 > Евгений СТРИГУНОВ > Толковым студентам -платникам государство будет помогать
13.07.2009
РИА "Новости " > 13.07.2009 > Около 28 тысяч студентов вузов переведены с платных на бюджетные места в 2009 году
8.07.2009
Письмо Рособразования от 29.05.2009 N 20-353 "О Конкурсе 2009 года на соискание медалей Российской академии наук с премиями для молодых ученых РАН..."
Логин
ВХОД
регистрация забыли пароль?
Пресс-релизы
Библиотека
Конференции
Выставки
Доска объявлений
Архив
Ассоциация технических Университетов Информация о проекте
Ниже мы рассмотрим вибрационную модель с линейной системой уравнений, основанной на линейной динамической характеристике силы резания, и вибрационную модель с нелинейной системой уравнений, учитывающих фрикционные явления в зоне резания. Кроме того, мы сделаем вывод о возможности применения данных моделей, сравнив расчетные данные с экспериментальными.
Расчетные схемы и основные уравнения, необходимые для построения моделей.
При точении инструментом с большим вылетом наименее жестким элементом в технологической системе является державка. Поэтому при расчете колебаний системы можно ограничиться только рассмотрением вибраций державки. Для упрощения расчетов, возьмем за основу схему свободного прямоугольного резания рис. 1( а).
а)
б)
Рис. 1. а) Принципиальная схема процесса ортогонального резания, б) Расчетная схема при свободном
прямоугольном резани/
Согласно схеме рис. 1(б), режущий инструмент представляется в виде точечной массы с двумя степенями свободы, совершающей колебания в направлениях Х и Z. В каждом направлении точечная масса подкреплена упругими и демпфирующими элементами. На массу действует также возмущающая сила, являющаяся результатом взаимодействия режущих кромок инструмента с материалом заготовки (сила резания).
По такой же схеме можно вести расчет и при использовании абсолютно жесткого инструмента и маложесткой детали, то есть при обычном точении деталей типа тонкостенных труб и дисков, деталей с большим вылетом из патрона без поджатия задней бабкой. При этом определяются динамические характеристики детали, а не инструмента, но алгоритм расчетов остается таким же.
В условиях динамического равновесия державки с учетом действующих на нее сил резания получим следующие линейные уравнения колебаний режущей кромки инструмента:
Мх ■ х + сх ■ х + кх ■ х =
. Мх X М_ -2 в -
Где л, , ¿1 - силы инерции колебательной системы,
Mx, Mz- массы колебательной системы,
X г % - скорость перемещений колебательной системы, сх, ^ - коэффициенты трения в системе (коэффициент демпфирования), кх, ^ - коэффициенты жесткости колебательной системы, Рх, Fz - составляющие силы резания, действующие на систему.
(5)
Рис. 2. Схема расчета ширины резания при свободном прямоугольном резании.
При свободном прямоугольном резании, например (рис.2), когда резец перемещается параллельно оси трубы, во
Авторы
Координационный совет
время первого оборота детали инструмента оставляет на гладкой поверхности трубы волнистые следы в результате колебаний державки в направлении подачи Х. Когда начинается второй оборот детали, поверхность уже имеет волнистость и по внутреннему контуру, по которому режет инструмент х(Ц, и по внешней поверхности резания в результате вибраций во время предыдущего прохода х(^ТТогда динамическую (мгновенную) толщину стружки можно представить используя линейную динамическую характеристику с запаздыванием:
М(Ц=80+х(Ц-х(Ъ- Т
Так как в рассматриваемом случае державка колеблется и в плоскости Е, то вибрационная составляющая перемещений в этом сечении Дс1 на следующем обороте будет равняться разности перемещений на этом обороте и на предыдущем Е(^Т0)
ДС= ЕГО^-То)
(7)
Тогда время запаздывания Т в плоскости Х с учетом вибрационной составляющей в плоскости Е будет определяться следующим образом:
,
(8)
где
К - радиус заготовки, шд-частота вращения заготовки. Преобразуя (8) с учетом (7), получим:
(9)
Тогда динамическая составляющая толщины стружки с учетом вибраций в системе с двумя степенями свободы в двух перпендикулярных плоскостях может быть выражена как:
(10)
Выражая (5) через (10) получим зависимость составляющих силы резания от динамической толщины срезаемого слоя Ь:
(11)
Тогда система уравнений перемещений инструмента, вызванных вибрациями, для случая с двумя степенями свободы будет выглядеть следующим образом:
Мх л +с, л +*,■ л = -К ■ Л(0-Л(г-(То+т~2(1~То)))1
М-а+с ¿+к 1= к АДр+л(0--(7Ъ+
2 6} К
(12)
Решение линейной системы уравнений движения резца численным методом в среде Ма^аЬ.
Для нахождения решения дифференциального уравнения п-го порядка с использованием функций пакета Ма^аЬ следует привести его к эквивалентной системе дифференциальных уравнений первого порядка.
Пусть
, тогда запишем систему уравнений (12) так:
(13)
Эта линейная систему уравнений интегрируется численно с использованием функций пакета Ма^аЬ. Для решения этой системы уравнений необходимо определить:
1) эмпирические коэффициенты модели сил резания,
2) коэффициенты жесткости, массы колебательной системы, коэффициенты демпфирования.
При проведении экспериментов по определению эмпирических коэффициентов модели силы резания использовались специальные образцы заготовок дисков с канавками по торцу. Выступы между канавками стачивались инструментом с заданными углами заточки, с постоянной глубиной и толщиной резания, с углом в плане ]=90°, что обеспечивает условия ортогонального резания. При этом измерялись значения составляющих силы резания в окружном направлении Лг, и в направлении подачи Лх, при различных значениях подачи вд. Для рассматриваемого случая обработки жаропрочного никелевого сплава ХН73МБТЮ были получены следующие коэффициенты:
Крх=1,6*109 Н/мм2, Кр2=5,2*109 Н/мм2
Динамические характеристики системы (коэффициенты жесткости, массы колебательной системы, коэффициенты демпфирования) были идентифицированы при помощи частотного анализа технологической системы и ее отдельных элементов. Этот метод заключается в ударном возбуждении и измерении вибраций тестируемой технологической системы. Полученные в результате измерений данные обрабатывались с помощью специальных программ обработки сигналов, которые используют Фурье-преобразование и описаны в работе [9]. Данная методика полностью представлена в работе [8]. Динамические параметры технологической системы, на которой проводились эксперименты по токарной обработке жаропрочного никелевого сплава ХН73 МБТЮ, следующие:
Жесткость кх=1.3*107Н/м, к2=1.7*107Н/м;
Масса колебательной системы: тх=0.063 кг, т2=0.096 кг;
Коэффициенты демпфирования в системе: сх=9%, с2=9%.
Рис.3. Экспериментально полученные амплитудо-частотные характеристики процесса резания сплава ХН73МБТЮ на режиме резания Ь =1 мм, V = 20 м/мин, вр = 0,1 мм/об: амплитуда 5мкм, частота 2478Гц.
Анализ экспериментальных (рис. 3) и расчетных (рис. 4) данных показал, что исходная система дифференциальных уравнений (13) не полностью отражает характер вибраций, так как линейные системы данного вида базируются на регенеративном механизме возбуждения вибраций и в системах с большой жесткостью и демпфированием дают результат с быстрым затуханием вибраций. Поэтому для более точных расчетов необходимо рассмотреть другую модель, более общую, которая представляет собой нелинейную систему дифференциальных уравнений, учитывающих еще и фрикционные явления между инструментом и заготовкой. Эта модель более точно отражает вид вибрационных перемещений и позволяет прогнозировать изменение амплитуды вибраций в зависимости от
скорости резания.
Рис.4. График изменения рассчитанных перемещений режущей кромки инструмента во времени при точении сплава ХН73МБТЮ на режиме резания Ь =1 мм, V = 20 м/мин, Эд = 0,1 мм/об.
Учет влияния фрикционных явлений на систему резания.
В процессе резания всегда присутствует сила трения, которая возникает между контактирующими поверхностями инструмента и заготовки. Тогда составляющие силы трения в направлениях X и 1 зависят от относительной скорости скольжения между инструментом и заготовкой и описываются нелинейной функцией, например полиномом третьей степени:
-У?)
где у-эмпирический коэффициент, (мин/м )2.
(14)
Тогда
динамическая
модель
автоколебаний
технологической
примет
вид:
Мх х - схх ■ (1 - )+кхх = -кгь ■ [&+^) - ^ - (ть + Мг2 -Сг1-(\-р2)+кг1- -К Ь ■ [Яэ + х(1) - х{1 - (То +
2(£)-2(1-То)
щК
2(£)-2(ь-То)
))]
(15)
Эта нелинейная система уравнений интегрируется численно с использованием функций пакета Ма^аЬ, как и
7
система (12). Но для этого необходимо наити зависимость между эмпирическим фрикционным коэффициентом ' и скоростью резания. Для установления этой зависимости были проведены испытания точением никелевого сплава ХН73 МБТЮ со скоростью резания 10-30 м/мин с замером амплитуды вибраций кромки инструмента. Геометрия режущего инструмента: а = 6°, Л = 7°, ф = 95°, Ф1 = 5°, у = 7°, К = 0.8 мм.
Так как автоматический поиск нужного значения коэффициента у в системе с нелинейными дифференциальными уравнениями занимает очень много времени и требует больших вычислительных мощностей, выполним его в ручном
У
режиме. То есть будем подставлять значения коэффициента * вручную и решать систему (15) до совпадения величины амплитуды экспериментальной с расчетной, требуемый результат достигается быстро и просто. Остальные коэффициента мы берем из ранее решенной линейной системы дифференциальных уравнений (13):
тх=0.063 кг, кх=1.3*107Н/м, сх=9%, Крх=1,6*109 Н/мм2;
т2=0.096 кг, к2=1.7*107Н/м, с2=9%, Кр2=5,2*109 Н/мм2.
Получим зависимость коэффициента * от скорости резания (рис.5). Аппроксимируем полученную зависимость к Г)
экспоненциальной функцией У 10325-5 (рис.5)
как наиболее подходящей по виду и имеющей наименьшее среднеквадратичное отклонение р=0,9974 по сравнению с другими видами аппроксимирующих функций: степенной (р =0,986), линейной (р=0,7) или полиномом второй степени (р=0,95) .
системы
Рис. 5. Аппроксимация зависимости значения коэффициента * от скорости резания V.
у
В итоге, заменяя в выражении (16) ' экспоненциальной функцией / " , мы получаем нелинейную
динамическую модель процесса токарной обработки описывающую автоколебания инструмента как системы с двумя степенями свободы, под действием нелинейных сил взаимодействия при резании, учитывающую регенеративный механизм возбуждения колебаний при движении инструмента по поверхности от предыдущего оборота и нелинейную зависимость силы трения от относительной скорости между инструментом и деталью:
¿ = Х, 1
V =-Г—■ + *(*) - х(1-(То +
к
1 г
£" =-■ + *(0 - х(1 - (То +
Л/.
(16)
Рассчитаем по системе дифференциальных уравнений (16) вибрационные перемещения кромки инструмента (рис.6) и амплитудно -частотную характеристику при токарной обработке сплава ХН73 МБТЮ на режиме резания Ь=1 мм, V = 20 м/мин, Б0 = 0,1 мм/об.
Рис.6. Пример численного решения системы нелинейных уравнений описывающих вибрационные перемещения кромки инструмента: Амплитуда 4,2мкм, частота 2400Гц.
Выводы:
Разработанная нелинейная система, учитывающая фрикционные явления, моделирует вибрации режущей кромки инструмента при обработке никелевых сплавов лучше, чем линейная, точнее отражая вид вибраций и достаточно точно описывая их параметры.
Используя систему (16), можно моделировать вибрации для технологической систем с различными жесткостными
характеристиками на разных режимах резания при обработке различных никелевых сплавов.
По результатам расчетов вибрационных характеристик процессов резания жаропрочных сплавов при получистовом точении, был сделан вывод о том, что жесткость технологических систем не должна быть меньше 107Н/м, т. к. возникающие вибрации снижают стойкость режущего инструмента.
Литература.
1. Подураев В.П. Обработка резанием с вибрациями . - М.: 1970. - 351 с.
2. Заковоротный В.А. Исследование динамической характеристики резания при автоколебаниях инструмента // Известия Северо-Кавказского научного центра высшей школы. Технические науки. - 1978. - №2. - С. 37-41.
3. Каширин А. И. Исследование вибраций при резании металлов. - М. - Л.: Изд-во АН СССР, 1944. - 237 с.
4. Меррит. Теория автоколебаний металлорежущих станков //
Конструирование и технология машиностроения. - 1965. - Т.87, №4. - С. 62-72.
5. Бобров В.Ф. Основы теории резания металлов. - М.: Машиностроение, 1975. - 344 с.
6. Сараванья-Фабис, Д"Суза. Нелинейный анализ устойчивости автоколебаний при резании // Конструирование и технология машиностроения . - 1974. - Т. 96, №2. - С. 292-299.
7. Мурашкин Л.С., Мурашкин С.Л. Прикладная нелинейная механика станков. - Л.: Машиностроение, 1977. - 192 с.
8. Altintas Y. Manufacturing automation. - Cambridge: Cambridge University Press, 2000. - 286 c.
9. Вибрации в технике. Т.2 // под редакцией М.Д. Генкина. - М.: Машиностроение, 1981. - 497с.
Публикации с ключевыми словами: точение, вибрации, нелинейная динамическая модель Публикации со словами: точение, вибрации, нелинейная динамическая модель Смотри так же:
• Учебное пособие «Обработка заготовок на станках токарной группы »
Тематические рубрики:
• Наука в образовании: Электронное научное издание
■ УЧАСТНИКАМИ
т Ассоциация технических Университетов Координационный совет Вузы Новости Информационное агентство УМО Вузов Jmaqazine@xware.ru т^п^фон (ядчч) 263-68-67 0 RSS П STOCK GROUP © 2003-2009 «Наука и образование: электронное научно-техническое издание»
