Динамика вибрационного сверления устройством с электромагнитным вибровозбудителем

ДИНАМИКА, ПРОЧНОСТЬ, НАДЕЖНОСТЬ \

УДК 621.9.01:517.91

А. М. Гуськов, В. В. Захаров

ДИНАМИКА ВИБРАЦИОННОГО СВЕРЛЕНИЯ УСТРОЙСТВОМ С ЭЛЕКТРОМАГНИТНЫМ ВИБРОВОЗБУДИТЕЛЕМ

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

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

В настоящей работе рассматривается расчетная схема устройства с электромагнитным вибровозбудителем для вибрационного сверления (рис. 1). Для возбуждения колебаний в системе используется двух-зазорный электромагнитный вибровозбудитель переменного тока, состоящий из двух электромагнитов с П-образными сердечниками 1 и двух жестко связанных между собой якорей 2 из ферромагнитного материала. Якоря крепятся на балках 3 постоянного сечения, которые связаны с дополнительной массой 4. Связь инструмента 5 с дополнительной массой моделируется в виде пружины с линейной жесткостью и демпфера.

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

Рис. 1. Расчетная схема устройства для вибрационного сверления

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

Математическая модель. Математическая модель электродинамического вибратора и инструмента при сверлении состоит из:

• системы нелинейных дифференциальных уравнений электромеханических колебаний;

• двухпараметрического закона резания;

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

Уравнения электромеханических колебаний. Рассматриваемая система является электромеханической. Динамика электромеханических систем описывается уравнениями Лагранжа-Максвелла [6, 7]. Для систем, в которых токи замкнуты, эти уравнения можно записать в виде

Ф r +

m

R dW = E

brs ^ ж r j

д Фя

s=1

1дТ_дГ + дП_0 | dw

dt dqk dqk dqk dqk'

r = 1,

k = 1,

n.

Здесь Фг — магнитные потоки; Ягз — активные сопротивления в контуре; Ж — энергия магнитного поля системы; Ег — сумма сторонних ЭДС в контуре г; т — число замкнутых неразветвленных контуров; qk — обобщенная координата системы; П, Т — потенциальная и кинетическая энергии системы; (к — непотенциальные силы; п — число степеней свободы системы. Первая группа уравнений описывает динамику электромагнитной составляющей системы, а вторая — механической.

На рис. 2 приведена расчетная схема вибровозбудителя. Предположим, что магнитное сопротивление общей части двух П-образных сердечников мало по сравнению с суммой магнитных сопротивлений

m;

воздушных промежутков и остальных участков сердечников. Это позволяет считать, что силовые линии замыкаются так, как показано на рис. 2, и не учитывать часть потока (штриховая прямая, см. рис. 2).

Предположим, что электромагниты одинаковые, размеры сечения сердечника велики по сравнению с расстоянием между якорем и сердечником, но малы по сравнению с длиной силовых линий. Это позволяет считать скалярный потенциал постоянным по сечению, а поле в промежутках между якорем и сердечником — однородным; поле же в ферромагнетике можно учесть, введя магнитное сопротивление Ям сердечника и якоря.

Энергия магнитного поля в системе определяется по формуле [12]

Рис. 2. Расчетная схема вибровозбудителя

w=2

дм(Ф2 + Ф2) +

2

Ф2

Цо3 Цо3

где Ф — магнитный поток в ферромагнетике; З — площадь зазора; ц0 — магнитная проницаемость воздуха; ( — зазор между якорями и сердечниками в невозбужденной системе; х — перемещение верхнего якоря, направленное в сторону увеличения зазора. Выражение (1) составлено в предположении, что магнитный поток через любое сечение магнитопровода и промежутки между сердечниками и якорями один и тот же. Величины Ям считаются постоянными, т.е. не учитывается гистерезис и насыщение ферромагнетика.

Обобщенными импульсами в системе являются не потоки Ф1 и Ф2, а потокосцепления пФ1 и пФ2, где п — число витков в обмотках электромагнитов. Поэтому первая группа уравнений Лагранжа-Максвелла для рассматриваемой системы будет иметь вид

пФк + Ягк = и

где Я — активное сопротивление обмотки электромагнита; 1к — ток в обмотке; и(£) = и0 эт(2п П ¿) — напряжение, подаваемое на обмотку; П — частота сети, Гц.

Связь между потоками и токами можно представить в виде [7]

2 (

Uli + U3 13 =

Ul2 - U3 13 =

Rm +

Rm +

^0 S 2 d S

1+x 1 - x

Ф1;

ф2

Рис. 3. Расчетная схема механической составляющей системы

где п3 — число витков обмоток в цепи подмагничивания; ¿3 — ток в цепи подмагничивания. Соотношения (3) — это закон Ома для магнитных цепей электромагнитов, причем параметры в левых частях равны значениям магнитодвижущих сил.

После подстановки соотношений (3) в уравнения (2) получаем

Ф i +

Ф 2 +

R Rm

n2

R RM

n

2

+ 2

+ 2

Rd n2^o S

Rd n2^0 S

1+x 1 - x

7->n3 • U0 . ,n ~ ,4

Ф1 - R—г3 = — sin(2n П t);

n2

n

Ф2 + Rn|гз = — sln(2n П t).

n2

n

Ток ¿3 в цепи подмагничивания определяется независимо от динамики системы, поэтому в дальнейшем считается известным.

Используя условие жесткого соединения якорей, расчетную схему механической составляющей системы можно привести к виду, показанному на рис. 3. Полученная система имеет две степени свободы: поперечные х и продольное ы перемещения якорей и инструмента. Остальные обозначения на рис. 3 имеют следующий смысл: т ь т2, т3 — суммарная масса якорей, масса дополнительного элемента и инструмента соответственно; I, Е,]х — длина и жесткость балки; Ц — пондеромоторная сила; Рс — осевая составляющая силы резания.

Продольные перемещения и1 якорей и дополнительной массы и2 выражаются через поперечное перемещение х якоря. Для получения этой связи представим форму изогнутой оси балок в виде уравнения = x(t)f (в), где 8 € [0,1] — безразмерная продольная координата сечений балки. Функция f (в) должна удовлетворять следующим граничным условиям:

s = 0

s = 1

v(t, 0) = 0, v(t, 1) = 1,

dv

ds dv

= 0;

s=0

ds

= 0.

s=1

Тогда перемещения щ^) и и2^) выражаются через поперечное перемещение якоря х следующим образом:

1 С1 ( ду\ 2

и^) = — ! ( — ) U2(t) = 2ul(t).

Подставляя выражение для в), получаем

х2 х2

« = 21 ¿1, и2 = ТЛ

- ¿1=I («в )2 «в -1.

Кинетическая энергия системы

Т = 1 [т1 х2 + т1 «2 + т2 «2 + т3 Ш2] = 2

I 1 I J1 2 \ -2 | -2

mi ( 1 + — р-21 x ) x + тз w

_ 1 = 2

где ^21 = 1 + 4 т^/ть

Потенциальная энергия системы

П = + 2 IV,

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

и =12Е3Х Г1 ()2 = Е7х¿х2

иизг=~2~1^1, «в =-¿г-^ x,

г1 («2/)2 , Л

где 32 = I «в. Окончательно потенциальную энергию си-

стемы можно записать как

2 Е ¿х ¿2 2 к [ ¿1 2 4 2

П = —13— х + х -'

Рассеяние энергии в механической системе учитывается с помощью диссипативной функции Рэлея:

в1х2 в2(и2 — Ш)2 в1х2 в2 {2 71 хх

Ф = ---+ —-— = ---+ — -;--Ш

2 2 2 2 \ I

где в1, в2 — коэффициенты демпфирования. Обобщенные силы вязкого трения определяются по уравнениям:

(х* = -= -Ах + 2 в¿гхШ - 4 в242"х2 х; дх 1 12

(ш* = -дф = 2 в2"Тх х - в Ш. дш 1

Пондеромоторная сила [7]

А = -^ = Л (Ф*2 - Ф.2

дх ^о Ь

После подстановки полученных выражений во вторую группу уравнений Лагранжа-Максвелла получаем

Ji2 2\ •• п  4 E 7Ж J2

ш1 + "ргх ^ X + +---х + N (ж, X, ад, ад) =

= ф22 - ф12; (5) £

кЛ 2 2 в2 Л . р

ш3 + р2 w + кад--— ж---— ж ж = Рс.

I I

Нелинейные слагаемые в первом из уравнений (5) имеют вид

лт, . .4 Ш1^21 ^12-2 2 к 71 2 в2 Л. ,

N (ж, ж, w, w) =-—-ж ж---— w ж---— w ж+

I2 1 1

4 в2 л2 . 2 2 к 72 3

+ жж + ж .

12 12

Уравнения (4) и (5) описывают динамику рассматриваемой электромеханической системы.

Сила резания. Кинематика резания. Для описания осевой силы резания воспользуемся следующим двухконстантным представлением [8]:

Т-, ксо Н0 -д

Рс =-Пд,

д

где кс0 = даь Лд [Н0/(пс Л)]д-1 — статическая жесткость резания; д — постоянная формы режущей кромки; аь — характерное напряжение обрабатываемого материала; Л — радиус сверла; д — параметр нелинейности закона резания, 0 < д ^ 1; Н0 — номинальная подача на оборот; пс — число режущих кромок; п — приведенная толщина снимаемого слоя, определяется мгновенными значениями толщины Н^ каждой из снимаемых режущих кромок:

П =

Пс ^ V c ho

i/q

При непрерывном резании без вибраций n = 1; сила резания совпадает с «квазистатическим» законом резания.

Рассмотрим обрабатываемую поверхность как сигнал, поступающий под режущую кромку в момент времени t и сформированный предыдущей режущей кромкой в момент времени t — T/nc, где T = 1/ш — период вращения детали, с; ш — частота вращения детали, Гц. Под j-й режущей кромкой в момент времени t величина этого сигнала равна расстоянию от поверхности торца до поверхности дна под кромкой с номером J = (j — 2modnc) + 1 в момент времени t — T/nc. Тогда

q

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

Dj(t) = -w(t) - H + vt - LJ(t - T/nc); hj (t) = max [0, Dj (t)] ; Lj (t) = Lj (t - T/nc) + hj (t); Lj (t) = Loj (t), t < 0.

Здесь ^ (t) — расстояние между з-й режущей кромкой и дном отверстия; Н — расстояние от инструмента до средней поверхности торца детали в начальный момент времени; V — постоянная скорость осевой подачи; Ь] ^ — Т/пс) — отклонение поверхности торца от плоскости под з-й кромкой; Ь^ (^ — характеристика обработанной поверхности (глубина отверстия вдоль образующей напротив з-й режущей кромки в текущий момент времени); Ь0^ (^ — начальная характеристика поверхности торца детали. В случае плоского торца, перпендикулярного осевому движению детали, следует предположить Ь0^ (^ = 0 при t ^ 0. В этом случае все режущие кромки имеют одинаковые геометрические характеристики.

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

Нормирование уравнений. Для удобства дальнейших вычислений полученные уравнения приводятся к безразмерному виду. В качестве масштаба времени берется время Т одного оборота детали, масштаба длины — номинальная подача на оборот Л,0 = vT, а для масштаба магнитных потоков вводится обозначение ф :

т = T; {фl, ф2> = -1{Фъ Ф2};

С S} = -1 {x w d}.

ho

Уравнения динамики системы (4) и (5) в безразмерном виде принимают вид:

ф1 + Р

ф2 + р

Рм + Р ( 1 + |

Рм + Р ( 1 - £

ф1 - г*| = u sln(2n vt); Ф2 + г*\ = u sln(2n vt);

(1 + P21 f) + 4 п aipi^' + 4 п2 pi2 e + N (e, e', Z, Z') =

=

= Ф2 - Ф1;

z " + 4 n «2 p2 z ' + 4 п2 p22 z + n (e, e', z, z о = 4 n2p2

2 2 K T

q

Здесь штрих обозначает дифференцирование по новому времени т и введены следующие безразмерные параметры электромагнитной составляющей системы:

Л 2¿Т

РР = ~2 -С '

П2 р0 6

R RmT

Р Рм

n2

. = R пз гзТ

Р г* 2

n2

Ф.

u=

UoT

n Ф* ,

v = П T.

Параметр р характеризует отношение активного сопротивления электромагнита к индуктивному, обычно р ^ 1 [7].

Безразмерные параметры механической составляющей системы можно представить как:

Т

Р1 =

ho

£ = Ji~l ,

4 J2 EJx

a1 =

А

2W mr

4 J2 EJx

«2 =

2 п V ш1 /3

в2

2>/ k тз'

P2 =

к =

2п kc 0

k ho

k Ш3

/3

Смысл параметров р1 и р2 — число свободных колебаний соответственно якоря и инструмента за один оборот детали.

Нелинейные слагаемые в последних двух уравнениях системы (7) имеют вид

N1 (е , е', С , С') = — £ (8 П2 Р31 Р22 се + 8 П Р31 «2 Р2 С' 0 +

+£2 (р21 е е'2 + 16 п рз1 «2 Р2 е2 е' + 8 п2 рз1 Р22 е3), N (е, е ',с,с') = —£ (4 п2Р22е2+8 п«2 р ее').

Здесь р31 = ш3/ш1. Также добавлено уравнение для определения масштаба магнитных потоков:

Т2 Ф2

* = 1.

ш1^0р06

Запишем безразмерные характеристики обрабатываемой поверхно-

сти как:

(Д,-, Л,-, п, Н} = -1 (^, Ь3, Н3, Н}. Я-0

Соответственно, уравнения образования новых поверхностей в процессе резания можно представить в безразмерном виде:

Aj(т) = -Z(т) - H + т - Aj(т - 1/nc);

n(т) = max[0, Aj(т)];

Aj(т) = Aj(т - 1/nc)+ nj(т);

Aj (т) = Aj0 (т), т < 0; J = (j - 2modnc) + 1.

Полученные уравнения электромеханических колебаний (7) и образования новой поверхности в процессе резания (8) являются нелинейными. Источниками нелинейности в системе служат:

• пондеромоторные силы;

• связь продольных перемещений якорей и дополнительной массы с поперечными перемещениями якорей;

• зависимость толщины срезаемого слоя от перемещений инструмента (прерывистость резания);

• нелинейность закона резания.

Аналитическое решение. Система дифференциальных уравнений (7) естественно содержит малый параметр е = 7 Ло// ~ ~ (1... 10)10-4. Поэтому в случае, когда сила резания равна нулю, можно получить приближенное аналитическое решение одним из асимптотических методов. Воспользуемся методом многомасштабных разложений. Этот метод позволяет получить равномерно пригодное решение, которое можно использовать для любых моментов времени т > 0.

Введем новые временные переменные [9]:

¿о = т, ¿1 = ет, ¿2 = е2т, ... , tj = еj т,

где ¿0 — основное, "быстрое" время; tj при ] > 0 — "медленное" время. Величины tj представляют собой различные масштабы времени. Производные по времени т преобразуются следующим образом:

= до + е д1 + е2д2 + ... ;

от

^ = до2 + 2 едо д1 + е2 (д!2 + 2 до д2) + ... ,

где дj = д/дtj. Решения уравнений (7) будем искать в виде рядов по е:

те

д = д(^ ¿ъ ...,е) = ^ еj ^j,

j=о

где д = {ф1, ф2, С, С}. Подставляя ряды в уравнения (7) и собирая слагаемые при соответствующих степенях е, получаем систему рек-курентных уравнений для определения составляющих искомых функ-ций1. Система для определения функций с индексом "0" называется порождающей и в рассматриваемом случае имеет вид

до ф1о = и эт(2п до ф2о = и эт(2п^о);

до2 & + 4п«1 р1 Со + 4п2р12 Со = Ф2о2 - Фю2; до2 Со + 4 п«2 Р2 Со + 4 п2 р22 Со = 0.

1Для дальнейших выкладок заменим ^ = ест, где ст > 0.

Решение этой системы (для перемещений используем только частные решения дифференциальных уравнений) можно представить следующим образом:

и

фю = А1 - -— шз(2п^о); и

Ф2о = А2 - ^^ шз(2п^о);

2nv

Со = -kci(Ai2 - A22) +-(Ai - A2) cos(2nvto + ^1);

nv

где

Co = 0,

ko =

1

ki =

4n2pi2i

i

v vPi

- ) - 1 I + 4ai2 ( -

Величины А1, А2 являются неизвестными функциями медленных масштабов времени. Для их определения используем уравнения для потоков, образованные коэффициентами при е в первой степени,

до фи = -di^io - о доф2! = -дlф2o - о

I I 1 I

Рм + И 1 + у

I 1

Рм + Р ( 1 - у

фЮ - i* ( ; ф20 + i*

Для обеспечения равномерной пригодности решений необходимо и достаточно, чтобы правые части уравнений (10) не содержали секуляр-ных слагаемых, т.е. слагаемых, которые приведут к появлению частных решений, неограниченно возрастающих с увеличением времени. Секулярных слагаемых не будет, если в разложения правых частей уравнений (10) в ряд Фурье не войдут постоянные составляющие:

г1/и г1/и

у\ Н1 = 0; у\ Н2 = 0.

оо

Здесь Н1, Н2 — обозначения для правых частей уравнений (10). Подставляя в условия (11) выражения для Hj и учитывая, что

д1фю = д1 Аь д1ф2о = д1 А2,

получаем систему уравнений для определения неизвестных функций Aj. Эти уравнения называются уравнениями установления и имеют

2

2

2

вид

diAi = Pi(Ai, A2) = а { Ai [а - b (Ai - A|)] - c (Ai - A2) - e};

diA2 = P2(Ai, A2) = а {A2 [а + b (A? - A2)] + c (Ai - A2) + e}.

Здесь введены следующие обозначения:

Г i/v

e = v г* dto, Jo

а = Рм + P,

b = kop,

2

u2P L / c = . „ 2 ki cos .

4n2v

Параметр е — это среднее за период возбуждения значение тока в цепи подмагничивания.

Стационарные движения в системе определяем, положив д?А? = = 5^2 = 0:

А? [а - Ь (А? - А2)] - с (А? - А2) - е = 0;

A2 [а + b (Ai - A2)] + c (Ai - A2) + e = 0.

Система нелинейных алгебраических уравнений (13) допускает решение вида А2 = -А?, называемое симметричным [7]. Это решение будет

А? =

а - 2с

Сложив уравнения (13) и затем сократив на А? + А2, что для несимметричных решений допустимо, получаем

а - b (Ai - A2)2 = 0.

С помощью уравнения (14) можно найти несимметричные решения в явном виде:

Ai = -

/а а - 2c ^ e * b \l b + v/Ob

A2 = -

I ±J + 2 V=

b V b -\7äb

Ai = -

A2 = -

a - 2c e -:- — 2

vVb

a - 2c e -:- — 2

v^öb

Устойчивость решения можно определить по собственным числам матрицы Якоби системы

/ др др \

а - с - Ь(3А2 - А2) с + 2 6А1 А2

с + 2 Ь А1 А2 а - с + Ь(А2 - 3А2),

J

dAi дА2 дР2 дР2

V^Ai дА2/

2

Решение будет устойчиво, если корни характеристического уравнения

Л2 + 81 Л + 82 = 0

матрицы (15) будут иметь отрицательные действительные части. Здесь коэффициенты 81, 82 имеют вид

81 = 2 а(Ь (А2 + А2) + с - а), 82 = а2 (а2 - Ь (А1 - А2)2 (-2 с + 3 Ь (А1 + А2)2) --2а (с + Ь (А12 + А22))).

Чтобы действительные части корней характеристического уравнения (16) были отрицательными, необходимо и достаточно, чтобы были положительными его коэффициенты:

81 > 0, 82 > 0.

Для несимметричных решений, используя уравнение (14), из условия для 82 получаем невыполняемое неравенство

Ь2 а2 (А12 - А22)2 < 0.

Отсюда делаем вывод, что несимметричные решения неустойчивы.

Для симметричного решения условия устойчивости решения имеют вид

(а - 2 с)2 (с - а) + 2Ье2 > 0;

ч2 4 Ье2

а (а - 2 с)2--> 0.

v ' а - 2с

В симметричном режиме первые гармоники сил, действующих на якорь со стороны обоих электромагнитов, складываются, остальные — уничтожаются, поэтому

ик1^ 2е

Co =

cos(2n vto + ф),

п V а - 2с

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

Со = АI cos(2пvtо + ф),

где

A =

e u

2 n2ppi3v

Рм + 1 P

\

-1

+4ai2 —

v

u2 cos p1 J 8 n4v 2p12i

"р = " +

Угол равен нулю, если знаменатель выражения для амплитуды больше нуля, и равен п — в противном случае.

Как видно, амплитуда колебаний якорей Ар линейно зависит от среднего за период возбуждения значения тока в цепи подмагничива-ния е и нелинейно зависит от амплитуды подаваемого напряжения и. Причем, существуют такие пары значений параметров {V, и}, которые обращают в нуль знаменатель выражения для амплитуды перемещений якоря. При этих значениях параметров в системе будут колебания, сопровождающиеся соударениями якорей и сердечников электромагнита, что подтверждается экспериментально [7].

Используя полученное решение для перемещения якоря, получаем дифференциальное уравнение для определения первого приближения перемещений инструмента:

С1' + 4 П«2 Р2 Z1 +4 П2 Р22 Zl =

= 2 п2 р22 a

v

1 + cos (2 T) - 4a2— sin (2 T P2

где Т = 2 п + "р. Частное решение уравнения (17) имеет вид Ар

Z1 = [1 + k cos(4n vto + 2 ^ + )], 2

где

k =

\

1 + 4а22 (

V2J

2А2 -1

Р2

+ 4а22 ( -

Р2

= arctg

2а2 I — )

Р2 )

2v

-

. — (1 - 4«22) - 1

На рис. 4 приведены амплитудно-частотные характеристики (АЧХ) перемещений якоря и инструмента (сплошные линии) для системы с

2

2

2

2

2

Рис. 4. Амплитудно-частотные характеристики:

а — якоря; 6 — инструмента; сплошные линии — аналитическое решение, точки — результаты численного интегрирования

S = 50, а2 = 0,05.

2

параметрами2:

и = 80, е = 5, р =2,8,

Р1 = 1, Р2 = 1,4, ах = 0,05,

Нелинейность системы дифференциальных уравнений не привела к качественному изменению АЧХ по сравнению с АЧХ линейной системы. На АЧХ для перемещения инструмента еще один локальный максимум расположен вблизи значения частоты возбуждения V, равного половине частоты собственных колебаний инструмента р2. Для сравнения на том же рисунке приведены значения амплитуд перемещений, полученные в результате численного интегрирования исходных уравнений (7) при Рс = 0 (амплитуды определялись как полуразмах установившихся вынужденных колебаний).

По результатам проведенного анализа динамики системы без резания можно сделать следующие выводы:

• уравнения для потоков допускают несколько стационарных решений: симметричное и несимметричные. Несимметричные решения либо не существуют, либо неустойчивы. Для симметричного решения получены условия устойчивости;

• для возбуждения в системе колебаний необходимо, чтобы среднее за период возбуждения значение тока е в цепи подмагничивания было отлично от нуля;

• колебания якорей происходят на частоте возбуждения, а колебания инструмента — на удвоенной частоте возбуждения;

• амплитуда колебаний якорей линейно зависит от среднего значения тока е в цепи подмагничивания и нелинейно — от амплитуды подаваемого напряжения и;

2При вычислении угла фх по формуле (9) необходимо контролировать непрерывность функции aгctg.

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

Численное моделирование системы в процессе резания. Для исследования динамики рассматриваемой технологической системы наиболее адекватным методом является использование отображения Пуанкаре в виде последовательности экстремумов перемещений якоря и инструмента или толщины слоя, снимаемого одной из режущих кромок (или приведенной толщины снимаемого слоя) — £ж£г[£], ], £х£г[п]. При этом на диаграмме для каждого значения параметра показываются все экстремумы функции, попавшие в интервал наблюдения.

На рис. 5, а показана диаграмма отображения Пуанкаре приведенной толщины снимаемого слоя для системы, в которой отсутствует ток в цепи подмагничивания, т.е. возбуждений колебаний инструмента нет. Изменяемым параметром для построения диаграммы является безразмерная величина р2 — число свободных колебаний инструмента за один оборот детали. На диаграмме видны зоны значений изменяемого параметра, в которых имеет место прерывистое резание с регулярной стружкой (авторезонансное вибрационное сверление [8]). В этих зонах происходит потеря устойчивости процессом резания с постоянной толщиной снимаемого слоя [8].

На рис. 5, б приведена диаграмма Пуанкаре (для системы с тем же набором безразмерных параметров) при вынужденных колебаниях инструмента, возбуждаемых вибропитателем. Безразмерная частота подаваемого напряжения V = ПТ = 1,1. В результате дополнительного возбуждения колебаний инструмента в ранее полученных зонах прерывистого резания изменяется поведение системы в процессе резания — стружка становится нерегулярной; а вблизи значения параметра р2 = 2, 2, равного удвоенному значению частоты возбуждения,

Рис. 5. Отображение Пуанкаре приведенной толщины снимаемого слоя для системы без возбуждения (а) и с возбуждением (б) вибраций инструмента; изменяемый параметр р

Рис. 6. Численное моделирование системы в процессе резания (р = 1,5):

а — без возбуждения вибраций инструмента; б — с возбуждением вибраций инструмента

на диаграмме появляется еще одна зона с прерывистым резанием и регулярной стружкой — зона синхронизации.

На рис. 6 приведены решения для обеих систем при значении параметра р2 = 1, 5. Показаны зависимости перемещений якорей и инструмента и приведенной толщины срезаемого слоя от времени, и построен аналог фазового портрета для толщины срезаемого слоя: по оси абсцисс откладывается текущее значение ), а по оси ординат — значение в некоторый предыдущий момент времени п(т — Т). В данном случае, для Т взято время запаздывания Т = 1/пс. По такому портрету можно судить о стационарности и регулярности значения толщины снимаемого слоя в процессе резания.

На рис. 7 приведены зоны регулярного прерывистого резания на плоскости двух параметров системы: р — число свободных колебаний

инструмента за один оборот детали; к — безразмерная статическая жесткость резания, включающая в себя технологические характеристики процесса резания. Точки на диаграмме, находящиеся в зонах 1, соответствуют значениям параметров, при которых имеет место авторезонансное вибрационное сверление с выходом инструмента. Рис. 7. Зоны регулярного прерывисто- _ ....

го резания на плоскости параметров При этом вибровозбУдитель вы-{p2, к} ключен. Точки, находящиеся в

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

Работа поддерживается грантами РФФИ № 05-01-08062, 07-08-00592-а и 07-08-00253-а.

СПИСОК ЛИТЕРАТУРЫ

1.Altintas Y.,Weck M. Chatter stability of metal cutting and grinding // CIRP Annals — Manufacturing Technology, 53 (2004), pp. 619-642.

2. Гуськов А. М. Динамика двух резцового точения // Станки и инструменты. -2004. -№ 11, 12.

3. Gouskov A. M., Voronov S. A., Paris H., Batzer S. A. Nonlinear dynamics of a machining system with two interdependent delays // Communications in Nonlinear Science and Numerical Simulation. December 2002. Vol. 7. No. 3, pp. 207-221.

4. T s e n g C. Y., T u n g P. C. Dynamics of a flexible beam with active nonlinear magnetic force // ASME Journal of vibration and acoustics. January 1998. Vol. 120. No. 1, pp. 39-46.

5. Tariq A. Nayfeh, Edward Emaci, Alexander F. Vakakis. Application of nonlinear localization to the optimization of a vibration isolation system // AIAA Journal. August 1997. Vol. 35. No. 8, pp. 1378-1386.

6. Вибрации в технике: Справочник. В 6 т. / Ред. совет: В. Н. Челомей (пред.). — М.: Машиностроение, 1979 — Т.2. Колебания нелинейных механических систем / Под. ред. И. И. Блехмана. 1979. — 351 с.

7. С к у б о в Д. Ю., Ходжаев К. Ш. Нелинейная электромеханика. — М.: Физматлит, 2003. — 360 с.

8. Г у с ь к о в А. М. Нелинейная динамика вибрационного сверления. Роль уравнений образования новых поверхностей // 4-й Международный конгресс "Конструкторско-технологическая информатика". - МГТУ-Станкин, 2000. -C. 123-130.

9. Н а й ф э А. Введение в методы возмущений. — М.: Мир, 1984. — 536 с.

Статья поступила в редакцию 15.06.2007