Устойчивость линейных дифференциальных уравнений с постоянным запаздыванием, описывающих процесс фрезерования

Устойчивость линейных дифференциальных уравнений с постоянным запаздыванием, описывающих процесс фрезерования

АЯ КрасильниковКЮ, Кравченко

В ходе изучения процессов механической обработки был обнаружен эффект наложения волн при возникновении автоколебаний [1, 2]. При этом поведение системы описывается дифференциальными уравнениями с запаздывающим аргументом (ДУЗА). Исследование устойчивости ДУЗА позволило построить диаграммы устойчивости на плоскости параметров резания. Графики, отображающие положение кривой границы устойчивости на плоскости «глубина резания - скорость вращения шпинделя», уникальны и однозначно определяются динамическими характеристиками отдельно взятой системы деталь-инструмент-приспособление-станок (ДИПС). В то же время, диаграмма устойчивости позволяет определить оптимальные режимы обработки, обеспечивающие устойчивость резания при высокой скорости снятия материала. В связи с этим вопрос об эффективном быстром построении таких диаграмм является актуальным.

Автоколебания всегда существуют при резании. Традиционно, во избежание неустойчивого режима автоколебаний параметры обработки значительно занижаются. Для проектирования эффективной обработки необходимо оптимизировать режимы резания по критерию асимптотической устойчивости.

В ряде работ, среди которых [3, 4, 5], представлена и исследована на устойчивость одномерная модель фрезерования. При механической обработке, в частности, при фрезеровании на обработанной поверхности остается след от колебаний режущей кромки, вызванных силой резания. При последующем резании этой поверхности режущая кромка режет уже с переменной глубиной резания, зависящей от колебаний предыдущего зуба

фрезы. Сила резания при этом, зависящая от глубины резания, также осциллирует.

Ниже представлена схема концевого фрезерования в ортогональной плоскости с одной степенью свободы в случае попутного направления фрезерования (рис.1).

X

Рис. 1. - Схема фрезерования с одной степенью свободы Представленная схема позволит описать модель фрезерования с учетом вращения инструмента. Уравнение движения рассматриваемой системы имеет вид (1).

^ + 2?Юсоб ^ +-L )=^ F

dt dt kynp

где q - относительный коэффициент демпфирования; юсоб - частота собственных колебаний, рад/с; купр - коэффициент упругости, Н/м; Т -постоянный период запаздывания, с; F(t) - сила резания, Н.

Мы рассматриваем систему с одним возможным смещением вдоль оси OX. Величина проекции силы F(t) на выбранную ось определяется (2):

Fx

j=1

где Ks - удельная сила резания, МПа; акр - осевая глубина резания, мм; f -подача на зуб, мм/зуб; фу - угловое положение j зуба относительно положительной полуоси OY; z - эффективное число зубьев, одновременно участвующих в резании; T - период запаздывания между двумя

последовательными резами, Т = 60/(zO6N), где N - скорость шпинделя, об/мин; zo6 - общее количество зубьев; в = л/2-y, где у - передний угол инструмента.

Постоянное значение К5акр/8тф/-8т(ф/-+Р) не влияет на устойчивость системы, поэтому его можно отбросить при дальнейшем анализе, исследуя лишь динамическую составляющую силы резания. В проекции на ось OX получим уравнения для компонента Fx (3).

Fx

j=1

где Ax = x(t-T) - x(t).

Параметр А(ф) изменяется во времени, т.к. ф = ф(t), причем Ае£2([ф5Ь ф^]), где фа - угол входа фрезы в контакт с заготовкой, фел. - угол выхода фрезы из контакта с заготовкой. Перейти к независимому от ф виду возможно путем рассмотрения первого члена соответствующего ряда Фурье

(4).

а +х>

А

2 r =1

где коэффициенты Фурье функции:

І ф^ І фex І фex

— J A

— — —

ф st Vst Vst

Рассматривая только первый член ряда Фурье мы накладываем ограничение на модель в виде малости периода контакта каждого зуба фрезы с заготовкой. Иными словами, рассматриваем случай, когда радиальная глубина резания составляет всего несколько процентов от диаметра фрезы. Однако, как показано в [6], рассмотрение только первого члена ряда Фурье возможно для всего диапазона значений радиальной глубины резания. Таким образом, переходя к инвариантному во времени виду, получим (6).

a І ф“ z |^ex

Ao = -f = ao = — J A

2 2<, 8п іф „

В рассматриваемом исходном уравнении (І) произведем следующие замены, согласно (7).

dt

dT

dT

Тогда уравнение (1) примет вид (8):

Уравнение (8) описывает рассматриваемую систему резания с одной степенью свободы в унифицированном виде с безразмерными коэффициентами a и k. Заметим, что коэффициенты a и k являются положительными вещественными числами.

Поставим задачу определить границу устойчивости уравнения (8) на

плоскости параметров T и k при фиксированном значении параметра a. Для этого исследуем нулевое состояние равновесия на устойчивость. Критерием асимптотической устойчивости является отрицательность вещественной части всех характеристических показателей дифференциального уравнения.

Другими словами, существует такое значение T, при котором корни характеристического уравнения на комплексной плоскости переходят из правой полуплоскости в левую полуплоскость и наоборот.

Теорема! При фиксированном значении параметра a уравнения

каждому значению параметра T соответствует такое значение параметра ^р, что для любых k < ^р решение уравнения асимптотически устойчиво при условии

Доказательство. Определим решение соответствующего однородного

уравнения как х( t) = Хехр[/ю t - у-]. Тогда х( t - Т) = х( £ )ехр[/0] согласно [5]. Подставим указанные выражения в уравнение (8).

где 0 = ю Т -єу; єу - сдвиг между фазой уі функции без запаздывания и фазой ут функции запаздывания; і - мнимая единица. Результатом решения уравнения (10) относительно ю является (11).

ю±

Юд І Юд — 2к — 1

где ю0 = +1 - а2/2 .

Уравнение (11) при рассмотрении полуоси ю± > 0, а также условие положительности параметров а и k приводит к набору ограничений, который гарантируют дальнейшую структуру области неустойчивости.

(а, к )а є(0,л/2 ) к

( 2 а

~2

м

+ а,+ю

у

Как показано в [5, 7] исходное уравнение (1) на мнимой оси может иметь корни ±ю+ и ±ю_ для всех ю± > 0, исходя из критерия устойчивости

Найквиста, при 0 = 2пЬ, где Ь = 0, 1, 2, ... В то же время значение фазового сдвига определяется как (12) согласно [7]. Исходя из вышеизложенного, получим (13).

єу = 2п — 2 аг^ап

Л 2 Л

1 — Ю±

т±

2п — 2 аг^ап

1 — ю±

: (к) ,

аю.

+ 2пЬ

Ь є N

Критическое значение ккр численно равно минимальному значению k, удовлетворяющему условию (13) для фиксированного значения Т.

Теорема доказана.

Зависимость (13) является уравнением кривой устойчивости. Ранее [1, 3, 5, 7, 8] кривая устойчивости задавалась параметрически. Очевидно, что

граница устойчивости на плоскости параметров Т и k состоит из счетного числа подобных элементов (Ь - натуральное число). Каждый элемент границы будет складываться из двух несимметричных половин - Т- -(к) и

Т+ -(к). Кривые Т_-(к) и Т+ -(к) сближаются, только если 0_-0+^0, что возможно в единственном случае, когда к ^ к0, где к0 = а2/2 + а.

Сделаем ряд замечаний. Во-первых, параметры а и к являются безразмерными. Безразмерной величиной, следовательно, также является и относительная частота колебаний ю±. В действительности, значение ю± есть отношение частоты осцилляций к собственной частоте системы юсоб. Во-вторых, ограничения (9) носят условный характер. Эти ограничения определяют существование вещественных корней на полуоси ю >0. Как показано, на этой полуоси ю имеет всего два корня - ю+ и ю-. Значение параметра к = к0 задает предел, ниже которого расположена зона устойчивых решений и нет участков неустойчивых решений. Выше предела к = к0 появляются зоны неустойчивости.

Одним из параметров резания является величина относительного радиального врезания Я/ = Ь/Бфр, где Ь - ширина фрезерования, мм; Бфр -диаметр фрезы, мм. Этот параметр характеризует глубину резания в ортогональной плоскости и также влияет на устойчивость обработки. При попутном фрезеровании Я/(фд) = 0,5(1-соБфд), при встречном фрезеровании Я!(фд) = 0,5(1+соБфд), где фд = фех - ф^.

Теорема! Функция акр(Я!) на отрезке Ше(0,1) имеет хотя бы один экстремум, причем при к=сот1

купр к

КА

Доказательство. Исходя из сделанной замены (7), получим (15).

а = купр к

кр КА

Очевидно, что поведение функции акр(Я/) при постоянных прочих параметрах будет зависеть от поведения функции А0(Я1). Выразим А0 из (6) через фд, причем фех = п для попутного фрезерования и ф5Г = 0 для встречного фрезерования (16).

«,(Я/) = 7Г^ (14>

А0

4п

где «-» - для случая попутного фрезерования, «+» - для случая встречного фрезерования, фде(-п, 0)П(0, п) исходя из Я/е(0,1). Диапазон фде(-п, 0) соответствует попутному фрезерования, фде(0, п) - встречному, исходя из условия положительности акр.

= +тЧс08Р_ с0з

dфд 4п

Таким образом, функция А0(фд) испытывает максимум при фд = (-п+Р) для попутного фрезерования и максимум при фд = в в случае встречного направления фрезерования. Соответствующая функция акр(Я/) достигнет минимума при Я/ = 0,5(1-собР) при попутном фрезеровании и минимума при Я/ = 0,5(1+собР) при встречном фрезеровании.

Теорема доказана.

Сделаем замечание. В формулировке теоремы 2 нет указания на вид экстремума, поскольку вид функции А0(фд) зависит от выбранной оси координат, в нашем случае это ось ОХ. При выборе основной оси ОУ функция А0(фд) также будет иметь экстремумы, но не обязательно максимумы.

С помощью изложенных теорем разработан метод оптимизации режимов резания по критерию асимптотической устойчивости. Ранее метод оптимизации параметров резания был предложен в [9], однако предложенный метод является последовательным. Кратко, метод, изложенный в [9], сводится к последовательности шагов:

1. для произвольного значения Я1=1 построить график границы устойчивости акр(Аг);

2. на графике задать пару (акр, Ы);

3. варьируя Я/ в пределах от 0 до 1, для выбранного на шаге два значения N определить множество значений акр;

4. на плоскости акр(Я1) расположить пары (акр, Я1).

Новый метод с параллельной оптимизацией параметров резания сводится к следующей последовательности вычислений:

1. задать пару (к, T), удовлетворяющую условию к < ккр и T = T (к), согласно теореме 1;

2а. вычислить N = 60юсоб/(2об T);

2б. для диапазона Л/е(0,1) вычислить значения акр при выбранном значении k.

3. построить график акр

Шаги 2а и 2б не являются зависимыми и могут быть выполнены параллельно. Таким образом, оптимизация параметров резания разделяется на оптимизацию скорости шпинделя и оптимизацию пары параметров -осевая глубина резания и относительное радиальное врезание. Также, как следует из (7), изменение жесткости системы влечет изменение параметра акр, а варьирование частоты собственных колебаний юсоб влияет на изменение величины скорости вращения шпинделя.

Разделение оптимизируемых параметров в значительной степени упрощает программное решение задачи, поскольку для реализации первого метода требуется перебирать массивы в поисках акр. В то же время новый

метод позволяет работать локально, поскольку изначально задана пара (к, T).

Применение программных методов является неотъемлемым аспектом математического моделирования [10, 11]. Для решения задачи оптимизации применялся язык объектно-ориентированного программирования Python 2.6.

На рис.2 изображена граница устойчивости для группы системы с коэффициентом относительного демпфирования 1,65%.

°'°0 10 20 30 40 50 60

Г

Рис. 2. - Граница устойчивости для а = 0,033 (серым цветом показана область неустойчивых решений) Соответствующая граница устойчивости для параметров Юсоб = 1308,11 Гц и кУпр = 4,459 105 Н/м, К = 1253 Мпа, в = 61°, г = 1, гоб = 4, ф5Г = -п/2, фех = -п приведена на рис. 3.

' 3000 4000 5000 6000 7000 8000 9000 10000

Скорость вращения шпинделя, об/мин

Рис. 3. - Граница устойчивости для а=0,0165, юсоб = 1308,11 Гц, купр = 4,459-105 Н/м, К = 1253 Мпа, Р = 61°, г = 1, гоб = 4, ф5Г = -п/2, фех = -п Для скорости шпинделя 9000 об/мин, что соответствует Т = 13,7 на рис.4 представлен график зависимости акр(Я1).

1.2

1.0

Р) 0.8

ё °'6

S'

S 0-4

0.2

N=9000.0

0.7 0.4 Об 0.8

Относительное радиальное врезание

1.0

Рис. 4. - График зависимости акр(ЯТ) при N = 9000 об/мин Очевидно, что кривая на рис.4 имеет минимум. По графику можно задать пару (акр, ЯТ), соответствующую устойчивому режиму автоколебаний. В результате оптимизации определяются три параметра резания, соответствующие устойчивому режимы автоколебаний.

В статье рассмотрена модель концевого фрезерования с одной степенью подвижности. Составлено соответствующее уравнение движения. Получено унифицированное линейное дифференциально уравнение с постоянным запаздыванием. Сформулирована теорема об асимптотической устойчивости унифицированного уравнения. В работе описан метод оптимизации параметров обработки по критерию асимптотической устойчивости.

Литература:

1. Tobias, S. Machine tool vibration / S. Tobias. - New York: Wiley, 1961. - 352 p.

2. Tlusty J., Polacek M. The stability of the machine tools against selfexcited vibrations in machining // International Research in Production Engineering, 1963, 1(1), P. 465-474.

3. Budak E,Altintas Y. Analytical prediction of chatter stability in milling -Part I: general formulation // Journal of Dynamic Systems, Measurement and Control, Transactions of ASME, 1998, 6(1), P. 22-30.

4. Stepan G. Modeling nonlinear regenerative effects in metal cutting // Philosophical transactions of the Royal Society of London, Series A., 2001, 359(1781), P. 739-757.

5. Красильников А.Я., Кравченко К.Ю. Исследование устойчивости систем с запаздыванием, описывающих процесс фрезерования в случае с одной степенью свободы // Вестник машиностроения, 2013. - №9. - С. 67-75.

6. Budak E. Mechanics and dynamics of milling thin walled structures: Ph.D thesis, the University of British Columbia, 1994. - 284 p.

7. Красильников А.Я., Кравченко К.Ю. Аналитические методы исследования устойчивости систем с запаздыванием, описывающие процесс фрезерования // Справочник. Инженерный журнал с приложением, 2013. -№9 - С. 23-31.

8. Insperger T., Stepan G., Turi J. State-dependent delay in regenerative turning processes // Nonlinear Dynamics, 2007, 47(1), P. 275-283.

9. Budak E., Tekeli A. Maximizing Chatter Free Material Removal Rate in Milling through Optimal Selection of Axial and Radial Depth of Cut Pairs // CIRP Annals - Manufacturing Technology, 2005, 54(1), P. 353-356.

10. Литвинов С.В., Клименко Е.С., Кулинич И.И. и др. Расчет на устойчивость стержней из ЭДТ-10 при различных вариантах закрепления [Электронный ресурс] // «Инженерный вестник Дона», 2011, №2. - Режим доступа: http://ivdon.ru/magazine/archive/n2y2011/415 - Загл. с экрана. - Яз.

рус.

11. Зотова Е.В., Панасюк Л.Н. Численное моделирование динамических систем с большим числом степеней свободы на импульсные воздействия [Электронный ресурс] // «Инженерный вестник Дона», 2012, №3. - Режим доступа: http://ivdon.ru/magazine/archive/n3y2012/933 - Загл. с экрана. - Яз.

рус.