- •С.П. Никитин
- •Никитин С.П.
- •ВВЕДЕНИЕ
- •1. РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ
- •1.1. Основы метода прямой аналогии
- •1.1.1. Сущность метода
- •1.1.2. Основные принципы метода прямой аналогии
- •1.1.3. Выделение в исходном объекте однородных физических подсистем
- •1.1.5. Установление связей между подсистемами
- •1.2.2. Проверка корректности
- •1.2.4. Линеаризация нелинейных уравнений
- •I.2.S. Построение линейной системы уравнений
- •2. МЕТОДЫ АНАЛИЗА ФИЗИЧЕСКОЙ СИСТЕМЫ ПО МАТЕМАТИЧЕСКОЙ МОДЕЛИ
- •2.1. Схема анализа физической системы по математической модели
- •2.3.6. Расчет частотных характеристик по передаточной функции системы
- •3. МАТЕМАТИЧЕСКИЕ МОДЕЛИ ТИПОВЫХ УЗЛОВ ФИЗИЧЕСКИХ СИСТЕМ
- •3.1. Моделирование рычажной системы
- •3.2. Моделирование взаимодействия твердых тел
- •3.2.1. Особенности моделирования динамики твердых тел
- •3.2.4. Моделирование взаимодействия двух твердых тел
- •4.1. Пример моделирования шпиндельного узла
- •4.5. Разработка математической модели плоскодоводочного станка «Растр»
- •4.6. Разработка математической модели тепловых процессов при резании
- •5. АНАЛИЗ ВЛИЯНИЯ ПРОЦЕССА РЕЗАНИЯ НА ДИНАМИКУ ТЕХНОЛОГИЧЕСКОЙ СИСТЕМЫ
- •5.1. Влияние параметров процесса резания на устойчивость системы с одной степенью свободы
- •5.2. Влияние параметров процесса резания на вынужденные колебания динамической системы с одной степенью свободы
- •ЗАКЛЮЧЕНИЕ
- •Список литературы
- •Оглавление
5. АНАЛИЗ ВЛИЯНИЯ ПРОЦЕССА РЕЗАНИЯ НА ДИНАМИКУ ТЕХНОЛОГИЧЕСКОЙ СИСТЕМЫ
5.1. Влияние параметров процесса резания на устойчивость системы с одной степенью свободы
На основе математических моделей, полученных методом прямой анало гии, для динамической системы токарной обработки с врезной подачей рас смотрим колебания системы по оси у. Метод прямой аналогии позволяет по следовательным отражением в модели физических явлений, имеющих место в динамической системе станка, выяснить их влияние на динамику и устойчи вость, определить оптимальную сложность математической модели процесса резания.
Для анализа использован пример прорезания канавки на токарном станке (см. рис. 4.17) [13]. Жесткость системы в точке резания ку= 108 Н/м, масса т = = 98 кг, постоянная времени демпфирования 7/, = 10'5 с, удельная сила резания к>л= 2 1 09 Н/м2, коэффициент силы резания кР= 10Об7 Н/м, постоянная времени процесса резания 7> = 3 1 0-4 с.
Для расчета динамических характеристик используется программа PAN, которая позволяет определять собственные значения системы, график переход ного процесса, амплитудно-частотные характеристики для анализа вынужден ных колебаний в системе, амплитудно-фазовые характеристики для анализа ус тойчивости замкнутой системы при резании.
Собственные значения в общем виде представляют собой комплексно
сопряженные корни |
|
Su+\ =-а±7со. |
(5.1) |
Анализ собственных значений позволяет судить о динамике разомкнутой динамической системы. Отрицательность вещественных частей (а) собствен ных значений говорит об устойчивости системы, только при этом условии, пе реходные процессы в разомкнутой динамической системе будут затухать, то есть с течением времени все отклонения от возмущений будут стремиться к ну лю. Время затухания переходного процесса определяется величиной вещест венной части корней. Мнимая часть корней дает информацию о собственных частотах системы
(0/_ |
(5.2) |
v; = 2я |
Об устойчивости разомкнутой системы можно судить также на основе ряда критериев (Гурвица, Рауса), которые можно получить без решения харак теристического уравнения. В настоящее время использование их все менее оп равданно, так как постоянно расширяющиеся возможности компьютеров по-
зволяют решать сравнительно слоэкные уравнения И получать собственные зна чения, которые несут многогранную информацию.
Устойчивость замкнутой системы чаще всего определяется на основе критерия Найквиста (см. п. 2.3.3). О степени устойчивости при этом судят по расстоянию А от точки пересечения графиком АФЧХ разомкнутой системы действительной оси до точки (-1,0) (см. рис. 2.15).
При анализе динамической системы станка необходимо всякий раз учи тывать возможную неустойчивость замкнутой системы при резании, так как из вестно, что одно апериодическое звено, в качестве которого обычно рассматри вается процесс резания, может при некотором значении коэффициента усиле ния (коэффициента силы резания кР) приводить к неустойчивости замкнутую систему.
На основе этих динамических характеристик рассмотрим влияние на ус тойчивость динамической системы различной сложности математической мо дели и полноты учета параметров процесса резания.
Динамическая система токарной обработки без учета резания представля ет собой колебательное звено с одной степенью свободы [12]:
ту+ hy у+ Суу = F. |
(5.3) |
Известно, что такое звено устойчиво при любых условиях.
Решение уравнения (5.3) для приведенных выше данных дает один ком плексно-сопряженный корень
5 12 = - 5 ± у-999,987.
Отсюда собственная частота системы
v = 159,2 Гц.
Если учесть в динамической системе процесс резайИя как апериодическое звено, то есть демпфированием процесса резания пренебречь, то математиче ская модель будет иметь вид (4.15). Решение этой модели дает три собственных значения:
5, =-3014,912,
523 =+2,456 ± у-1022,155.
Следовательно, учет процесса резания He3Ha4HfenbH0 увеличивает собст венную частоту системы до v = 163 Гц, что подтверждается практикой. Поло жительное значение действительной части корня а = 2,456 свидетельствует о неустойчивом состоянии разомкнутой динамической системы, включающей в себя процесс резания с указанными выше параметрами.
При учете процесса резания в виде полной линейной модели, то есть с учетом демпфирования процесса резания, математическая модель будет иметь вид:
Му у+ hy у+ Суу + куР = О,
-^-Р+Р = hPy у+ кРуу.
1р
Решение такой математической модели дает следующие собственные значения:
Sx = -3001,497,
S2,з =-4,25± 7-1024,43.
Учет демпфирования процесса резания существенно влияет на динамиче ские характеристики, прежде всего на вещественную часть комплексно сопряженного корня, то есть на устойчивость разомкнутой системы. Значение а = - 4,25 свидетельствует об устойчивом состоянии системы. На основе этого можно констатировать, что приближенное значение собственной частоты мо жет быть получено уже по динамической модели без учета процесса резания. О колебательных свойствах системы, устойчивости ее состояния можно судить только по динамической модели с учетом отражения в ней процесса резания в виде полной линейной модели.
На рис. 5.1 приведены АФЧХ разомкнутой динамической системы при учете процесса резания как апериодического звена и в виде полной линейной модели. Для того и другого случая графики существенно различаются. Для процесса резания как апериодического звена в замкнутом состоянии мы имеем устойчивое состояние, а для полной линейной модели процесса резания
Рис. 5.1. АФЧХ динамической системы токарной обработки при учете про цесса резания: 1 - как апериодического звена; 2 - в виде полной линейной модели
система потенциально неустойчива, то есть график АФЧХ пересекает отрицательную часть действительной оси правее точки (-1,0).
Далее рассмотрено влияние отдельных параметров динамической системы на ее устойчивость в замкнутом и разомкнутом состоянии, а также другие динамические характеристики на основе модели, представленной на рис. 4.17. Об устойчивости системы в разомкнутом состоянии будем судить по величине вещественной части а собственных значений, а в замкнутом состоянии - по величине отрезка Д от точки пересечения графиком АФЧХ действительной оси до координаты (-1,0) (см. рис.2.15).
Наиболее сильное влияние на устойчивость процесса обработки оказыва ет коэффициент резания кР, который определяется выражением
кР =кух-Ь, |
(5-5) |
где куЛ- удельная сила резания, приходящаяся на единицу площади срезаемо го слоя;
b - ширина срезаемого слоя.
На рис 5.2 представлены зависимости значений А и а от значения кРдля рассматриваемого примера. Они показывают, что с увеличением коэффициента резания устойчивость динамической системы в разомкнутом состоянии и в замкнутом состоянии снижается. Причем лимитирующим критерием является устойчивость в замкнутом состоянии, так как точка пересечения графика А с
осью абсцисс находится левее (уже при кР = |
11,3*106 система в замкнутом со |
||||||
стоянии оказывается на границе устойчивости). |
|
||||||
|
|
1 |
1 |
1 |
1 |
1 |
|
|
|
устойчиво |
|
||||
|
|
|
|
|
|
|
[J а |
|
|
|
|
|
|
|
А |
—' — |
|
|
|
|
|
А |
г |
|
Г Кь |
|
|
|
|
|
|
11 •10 |
* 1 К « |
|
|
|
|
1 HI |
|
|
|
1 |
|
|
|
|
|
|
|
■ч |
|
|
|
|
|
|
|
ч ' ч |
> |
|
|
1 |
Л |
|
|
|
|
|
|||
неустойчиво |
|
|
1 |
|
и |
||
|
|
|
|
|
|||
|
|
|
|
\ |
|
|
|
1 1 |
I I 1 1 |
|
|
|
\ |
|
|
|
|
|
|
|
|
Рис.5.2. Влияние кРна устойчивость системы: а - в разомкнутом состоянии; А - в замкнутом состоянии
Значение собственной частоты с изменением коэффициента резания не связано, оно практически остается постоянным.
Другой характеристикой процесса резания является постоянная времени стружкообразования 7>, которая характеризует запаздывание изменения силы резания при изменении сечения стружки. При увеличении постоянной времени стружкообразования устойчивость системы в разомкнутом состоянии (а) и в замкнутом состоянии (Д) снижается (рис. 5.3).
Рис. 5.3. Влияние постоянной времени стружкообразования (7>) на устойчивость системы: а - в разомкнутом состоянии; А- в замкнутом состоянии
Лимитирующим критерием является устойчивость в замкнутом состоя нии. Значение собственной частоты при изменении постоянной времени струж кообразования Тр практически остается неизменным.
Большое значение имеет коэффициент демпфирования резания /*/>, кото рый зависит от жесткости процесса резания, постоянной времени демпфирова ния и определяется выражением:
|
hp = Г/, - к],, |
(5.6) |
где Г/, - |
постоянная времени демпфирования резания; |
|
к,} |
- установившийся коэффициент резания. |
|
Увеличение значения коэффициента демпфирования за счет изменения постоянной времени демпфирования (кР= const) приводит к повышению устоичиврсти динамической системы как в разомкнутом, так и замкнутом состоянии (рис. 5.4). Значение собственной частоты практически не зависит от коэффици ента демпфирования (постоянной времени демпфирования) резания.
Еще одной характеристикой, которая связана с условиями процесса реза ния, является величина коэффициента передачи К. Значение этого коэффициен-
га определяется проекцией силы резания на соответствующую координатную ось:
1&
>
1
К = c o s (a ) . |
|
|
|
|
|
(5.7) |
|
|
У f |
1 |
1 |
1 |
1 |
1 |
1 |
a |
устойчиво |
|
|||||
|
7 □ |
|
|
|
|||
v\ |
|
|
|
|
|||
< |
|
) |
|
дь |
|
|
• |
|
|
|
|
|
|
||
Л |
4 |
V 0^ |
|
|
|
||
/ |
|
|
|
|
|
|
|
JГ |
|
|
|
|
|
|
|
/й |
неустойчиво |
|
|
|
|
|
-t
Рис.5.4. Влияние постоянной времени демпфирования резания (7а) на устой чивость системы: a - в разомкнутом состоянии; Д - в замкнутом состоянии
Таким образом, величина К может лежать в пределах от 0 до 1. Увеличе ние коэффициента передачи приводит к снижению устойчивости системы в ра зомкнутом и замкнутом состояниях. Определяющим является критерий устой чивости в замкнутом состоянии. При заданных параметрах система становится неустойчивой при превышении К = 0,7 (рис. 5.5).
Рис. 5.5. Влияние коэффициента передачи (К) на устойчивость системы: a - в разомкнутом состоянии, Д - в замкнутом состоянии
Параметры процесса резания определяются условиями обработки, поэто му влиять через них на устойчивость динамической системы не представляется возможным. Предельные значения этих параметров, обеспечивающие устойчи вую работу станка, являются критериями динамического качества конструкции. Так как значения постоянных времени не определены, то в качестве такого кри терия чаще всего выступает критическая величина коэффициента резания кР или критическая ширина срезаемого слоя для данного обрабатываемого мате риала.
Для конструктора важно знать влияние конструктивных параметров на устойчивость динамической системы. К ним относятся жесткость, демпфирова ние и инерционность элементов динамической системы. Как показывают ре зультаты исследования, такое влияние наблюдается и имеется возможность управлять устойчивостью.
На рис. 5.6 представлены зависимости критериев устойчивости Л,а от податливости упругой механической системы. Устойчивость в разомкнутом со стоянии (а) несущественно зависит от податливости. Только в области малых значений это влияние значительно и снижение податливости увеличивает ус тойчивость. Зависимость устойчивости в замкнутом состоянии от податливости представляет собой кривую, которая имеет минимум при 1/С = 5ТО"8 м/Н.
А,а
1 1 L.1 1 1
устойчиво
й
т |
т |
о |
ш |
А |
L |
|
|||||
|
|
|
|
неустойчиво |
м/Н |
|
|
|
|
|
Рис. 5.6. Влияние податливости (/,) на устойчивость системы: а - в разомкну том состоянии; А- в замкнутом состоянии
Величина податливости несущей системы определяет величину собст венной частоты, наряду с инерционностью. Влияние массы несущей системы на устойчивость динамической системы представлено на рис. 5.7.
т п Г г г т
устойчиво
1 |
а |
|
|
|
|
|
|
|
Д |
* |
—■ |
6 |
8 |
-1 - : JL )— |
|
г |
И " Т |
||
|
|
|
. |
неустойчиво |
Рис. 5.7. Влияние массы (т) на устойчивость системы: а - в разомкнутом состоянии; Л - в замкнутом состоянии
Устойчивость разомкнутой динамической системы с ростом массы не сущей системы уменьшается, асимптотически приближаясь к нулю. Устойчи вость в замкнутом состоянии также снижается, достигая критического значения уже при т = 500 кг. При дальнейшем увеличении массы устойчивость в замк нутом состоянии вновь повышается, однако колебательность системы, устой чивость разомкнутой системы продолжает снижаться. Поэтому устойчивая ра бота для заданных параметров возможна при массе менее 500 кг.
Влияние коэффициента демпфирования элементов несущей системы И\ на устойчивость системы представлено на рис. 5.8.
I |
I |
7 |
"рг |
1 1 1 |
ГТ |
|
а |
/ |
|
устойчиво |
|
|
r*N |
/ |
|
|
|
|
|
L |
|
|
|
А/ / |
|
_ i Л — / / . - — - |
|
|
|
/ |
Ъ |
■ 4 |
|
; |
' л |
|
|
/с
>/
неустойчиво
/1
Рис. 5.8. Влияние коэффициента демпфирования (h\) на устойчивость системы: а - в разомкнутом состоянии; Д - в замкнутом состоянии