Проблемы машиностроения и надежности машин, 2021, № 1, стр. 15-24
Нелинейная динамика вибрационной машины с электродинамическим возбудителем колебаний
В. К. Асташев 1, *, К. А. Пичугин 1, Е. Б. Семенова 1
1 Институт машиноведения им. А.А. Благонравова РАН
Москва, Россия
* E-mail: v_astashev@mail.ru
Поступила в редакцию 22.01.2020
Принята к публикации 22.10.2020
Аннотация
Приводятся результаты анализа динамических свойств вибрационной технологической машины с электродинамическим возбудителем колебаний при ее работе на нелинейную технологическую нагрузку. Построены амплитудно-частотные характеристики машины, как на холостом ходу, так и в рабочем режиме при питании возбудителя колебаний от источников тока и напряжения. Показана возможность и найдены условия появления областей неоднозначности амплитудно-частотных характеристик и их ветвей, отвечающих неустойчивым режимам. Выявлен ряд особенностей динамического поведения машины и ее предельные возможности.
Электродинамические возбудители колебаний [1, 2], используемые в испытательных стендах, находят применение и в разнообразных вибрационных и виброударных машинах и устройствах в качестве привода рабочих органов [3–5], особенно при взаимодействии рабочего органа с нелинейной технологической нагрузкой, возникающей при выполнении рабочего процесса. В статье выполнен анализ динамических характеристик машины при ее работе, как на холостом ходу, так и в рабочих режимах при питании возбудителя колебаний от источника тока и источника напряжения.
На рис. 1 показана схема вибрационной машины с электродинамическим возбудителем колебаний. В корпусе 1 жестко закреплена магнитная система, кольцевой магнит 2, который зажат между магнитопроводами 3 и 4. В кольцевом зазоре магнитопроводов размещена катушка 5 с возможностью перемещения вдоль оси устройства. Катушка жестко установлена на центральном стержне 6, который связан с корпусом упругими мембранами 7. На конце стержня установлен инструмент, взаимодействующий с обрабатываемым изделием или средой 8. Инструмент поджимается к изделию статической силой P, приложенной к корпусу машины.
Катушки получают питание от источника переменного тока. Проходящий через обмотку катушки переменный ток $i(t)$ создает переменную силу $f(t)$ возбуждения колебаний катушки [1, 2] и жестко связанного с ней инструмента
где B – магнитная индукция магнитного поля в зазоре; l – длина проводника обмотки катушки; i – сила тока.В то же время движение катушки в магнитном поле создает в ее обмотке противоэлектродвижущую силу
где $x(t)$ – закон движения катушки.Взаимодействие рабочего органа с обрабатываемой средой создает технологическую нагрузку на колебательную систему. Для описания технологической нагрузки используется нелинейная динамическая характеристика ${{f}_{l}} = {{f}_{l}}(x,\dot {x})$ рабочего процесса [6, 7], которая определяет зависимость силы взаимодействия инструмента с деталью от перемещения $x$ и скорости $\dot {x}$ инструмента.
Режимы холостого хода. Рассмотрим сначала динамические свойства системы при работе на холостом ходу, т.е. при отсутствии технологической нагрузки ${{f}_{l}} = 0$. Состояние электромеханической системы при установившихся гармонических колебаниях с учетом (1) и (2) описывается уравнениями
(3)
$\begin{gathered} (с - m{{\omega }^{2}} + {\text{j}}\omega b)\tilde {a} = Bl\tilde {I}, \\ (R + {\text{j}}\omega L)\tilde {I} = \tilde {U} - {\text{j}}\omega Bl\tilde {a}, \\ \end{gathered} $Если сила тока задана, т.е. питание катушек производится от источника тока, и изменяется
по гармоническому закону , где I = const – амплитуда cилы тока, из уравнений (3) находим комплексные амплитуды колебаний
инструмента и напряжения питания катушек
(4)
$\tilde {a} = \frac{{Bl}}{{W\left( {{\text{j}}\omega } \right)}}I,\quad \tilde {U} = \left[ {Z\left( {{\text{j}}\omega } \right) + \frac{{{\text{j}}\omega {{{\left( {Bl} \right)}}^{2}}}}{{W\left( {{\text{j}}\omega } \right)}}} \right]I,$![](/issues/pminm/2021/vol_2021/iss_1/PMiNM2101006Astashev/PMiNM2101006Astashev-m18.gif)
![](/issues/pminm/2021/vol_2021/iss_1/PMiNM2101006Astashev/PMiNM2101006Astashev-m19.gif)
![](/issues/pminm/2021/vol_2021/iss_1/PMiNM2101006Astashev/PMiNM2101006Astashev-m20.gif)
Из соотношений (4) после подстановки принятых обозначений получим выражения для комплексных амплитуд колебаний инструмента и напряжения питания катушек
Используя формулу Эйлера
где $\varphi $ – начальная фаза колебаний, из первого равенства (5) находим соотношения для амплитуды и начальной фазы колебаний рабочего органа(7)
$\cos \varphi = \frac{{a(c - m{{{\omega }}^{2}})}}{{BlI}},\quad \sin \varphi = \frac{{ba{\omega }}}{{BLI}}.$Выражения (6) и (7) эквивалентны соотношениям, получаемым при анализе вынужденных колебаний осциллятора при действии гармонической силы с амплитудой ${{F}_{0}} = BlI$.
Второе соотношение в (5) определяет напряжение питания. Второе слагаемое отражает влияние параметров механической части системы на импеданс электрической цепи. При питании от источника тока резонанс достигается при $\operatorname{Re} W\left( {{\text{j}}\omega } \right) = 0$, т.е. на собственной частоте ${{\omega }_{0}} = \sqrt {c{\text{/}}m} $ механической части системы. При этом достигают максимальных значений амплитуды, как колебаний инструмента, так и напряжения на выходе источника питания. Резонансные значения амплитуды колебаний инструмента и напряжения питания имеют вид
Отсюда видно, что в резонансном режиме колебания магнита создают добавочное активное электрическое сопротивление, которое тем сильнее, чем меньше коэффициент b вязкого сопротивления движению механической системы. Это объясняется возрастанием амплитуды колебаний при уменьшении b, а, следовательно, и противоэлектродвижущей силы. Заметим, что напряжение питания $U \to \infty $ при $\omega \to \infty $, что объясняется возрастанием индуктивного сопротивления катушки.
Питание электродинамического возбудителя колебаний, как правило, производится от источника
напряжения. В этом случае сила тока в проводнике катушки зависит не только от параметров
катушки. При движении в магнитном поле в проводнике катушки возникает противоэлектродвижущая
сила $e = - Bl\dot {x}$. В этом случае состояние системы описывается уравнениями (2), в которых заданным является напряжение . Решая систему уравнений (2) относительно неизвестных комплексных амплитуд $\tilde {a}$ и $\,\tilde {I}$ получим
(8)
$\tilde {a} = \frac{{BlU}}{{Z\left( {{\text{j}}\omega } \right)W\left( {{\text{j}}\omega } \right){\text{ + j}}\omega \left( {Bl} \right)}},\quad \tilde {I} = \frac{{W\left( {{\text{j}}\omega } \right)U}}{{Z\left( {{\text{j}}\omega } \right)W\left( {{\text{j}}\omega } \right) + {\text{j}}\omega {{{\left( {Bl} \right)}}^{2}}}}.$Для более наглядной физической интерпретации полученных далее решений запишем эти уравнения в виде
(9)
$\tilde {a} = \frac{{BlU}}{{Z\left( {{\text{j}}\omega } \right){{W}_{\Sigma }}({\text{j}}\omega )}},\quad \tilde {I} = \frac{{W({\text{j}}\omega )U}}{{Z({\text{j}}\omega ){{W}_{\Sigma }}({\text{j}}\omega )}},$Резонанс системы при ее питании от источника напряжения согласно (10) возникает при $\operatorname{Re} {{W}_{\Sigma }}\left( {{\text{j}}\omega } \right) = 0$, т.е. при частоте $\omega = {{\omega }_{r}}$, определяемой равенством
Поскольку коэффициент k зависит от частоты $\omega $, соотношение (11) является уравнением относительно искомой резонансной частоты ${{\omega }_{r}}$. В результате решения этого уравнения находим
Как следует из (12), при питании от источника напряжения резонансная частота системы оказывается более низкой, чем при питании от источника тока, за счет сложного сочетания параметров магнитоэлектрической и механической частей системы. В данном случае резонанс реализуется при ${{\left( {Bl} \right)}^{2}} < Lc$. При ${{\left( {Bl} \right)}^{2}} \to Lc$ резонансная частота ${{\omega }_{r}} \to 0$. Электродинамические возбудители с большими значениями ${{\left( {Bl} \right)}^{2}} \gg Lc$ используются в вибрационных стендах, позволяя получить практически постоянные амплитуды колебаний в достаточно широком диапазоне частот.
Из соотношений (9) с учетом равенств (10) и (11) найдем величины резонансных амплитуд колебаний инструмента и тока в цепи питания катушек
(13)
${{a}_{r}} = \frac{{BlU}}{{b{{{({{\omega }_{r}}L)}}^{2}} + R{{{\left( {Bl} \right)}}^{2}}}},\quad {{\tilde {I}}_{r}} = \frac{{(c - m\omega _{r}^{2} + jb{{\omega }_{r}})U}}{{b{{{({{\omega }_{r}}L)}}^{2}} + R{{{\left( {Bl} \right)}}^{2}}}}.$Рассмотрим поведение системы при двух характерных частотах. При частоте $\omega = 0$ из уравнений (8) с учетом принятых обозначений находим $I = U{\text{/}}R$, $a = BlU{\text{/}}Rc$, что вполне очевидно, так как система находится в состоянии статического равновесия. В практике эксплуатации вибрационных машин с электродинамическим возбудителем колебаний зачастую в качестве рабочей выбирают настройку на собственную частоту $\omega = {{\omega }_{0}}$ = $\sqrt {c{\text{/}}M} $ механической колебательной системы. В таком режиме амплитуда колебаний рабочего органа
Этот режим интересен тем, что амплитуда силы тока минимальна и при $b \to 0$ амплитуда силы тока $I \to 0$, т.е. механические колебания происходят, а ток в катушках отсутствует. Это происходит потому, что в этом режиме противоэлектродвижущая сила, возникающая в проводнике катушки, оказывается равной напряжению питания. По сути, в данном случае мы имеем дело с явлением антирезонанса в электрической цепи.
На рис. 2 показаны амплитудно-частотные характеристики колебаний рабочего инструмента при питании системы от источника тока (кривая 1) и источника напряжения при ${{\left( {Bl} \right)}^{2}} < Lc$ (кривая 2) и при ${{\left( {Bl} \right)}^{2}} > Lc$ (линия 3).
Нелинейная технологическая нагрузка. Перейдем к описанию работы машины с электродинамическим возбудителем колебаний при выполнении технологического процесса. Взаимодействие рабочего органа с обрабатываемым изделием или средой создает дополнительную технологическую нагрузку на рабочий орган машины. Как правило, нагрузку можно представить в виде действующей на рабочий орган силы ${{f}_{l}} = {{f}_{l}}(x,\dot {x})$, нелинейно зависящей от координаты $x$ и скорости $\dot {x}$ рабочего органа. Уравнение движения при питании от источника тока по аналогии с (3) принимает вид
Отыскивая приближенное гармоническое решение нелинейного уравнения (14), проведем гармоническую линеаризацию [6] функции ${{f}_{l}} = {{f}_{l}}(x,\dot {x})$
(15)
${{f}_{l}}\left( {x,\dot {x}} \right) \approx {{P}_{l}}(a) + {{c}_{l}}(a)x + {{b}_{l}}(a)\dot {x},$(16)
$\begin{gathered} {{P}_{l}}(a) = \frac{1}{T}\int\limits_0^T {{{f}_{l}}(x,\dot {x})} dt, \\ {{c}_{l}}(a) = \frac{2}{{Ta}}\int\limits_0^T {{{f}_{l}}(x,\dot {x})} \cos (\omega t)dt,\quad {{b}_{l}} = (a) = - \frac{2}{{Ta\omega }}\int\limits_0^T {{{f}_{l}}(x,\dot {x})} sin(\omega t)dt, \\ \end{gathered} $Поскольку вся система поджимается к обрабатываемой среде постоянной силой P (рис. 1), очевидно, что постоянная составляющая нелинейной нагрузки Pl = P, а уравнение колебаний (14) с учетом (15) принимает вид
(17)
$m\ddot {x} + (b + {{b}_{l}})\dot {x} + \left( {c + {{с}_{l}}} \right)x = BlI{{{\text{e}}}^{{j\omega t}}}.$Таким образом, получено линейное уравнение, коэффициенты которого зависят от неизвестной амплитуды. Решая уравнение (17), получаем выражение для комплексной амплитуды
(18)
$\tilde {a} = - \frac{{BlI}}{{\left( {c + {{с}_{l}}} \right) - m{{\omega }^{2}} + j\omega {{{(b + {{b}_{l}})}}_{l}}}}.$Заметим, что по форме выражение (18) совпадает с выражением (5), полученным при описании работы на холостом ходу. Поэтому проведенные выше выкладки формально оказываются справедливыми при замене величин c и b на величины $\left( {c + {{с}_{l}}} \right)$ и $\left( {b + {{b}_{l}}} \right)$ соответственно. Существенное различие равенств (5) и (18) заключается в том, что последнее является алгебраическим уравнением относительно амплитуды a, от которой зависят коэффициенты ${{с}_{l}}$ и ${{b}_{l}}$.
В качестве примера рассмотрим работу системы, предназначенной, например, для внедрения инструмента в хрупкую среду. Примем, что рабочий процесс описывается характеристикой жесткопластического материала [7]
(19)
$f\left( x \right) = \frac{1}{2}D\eta \left( {x - \Delta } \right)\left( {1 + \operatorname{sign} \dot {x}} \right),$(20)
${{P}_{l}}\left( a \right) = \frac{D}{{2\pi }}\arccos \frac{\Delta }{a},\quad {{c}_{l}}\left( a \right) = \frac{D}{{\pi a}}\sqrt {1 - {{{\left( {\frac{\Delta }{a}} \right)}}^{2}}} ,\quad {{b}_{l}}\left( a \right) = \frac{D}{{\pi a\omega }}\left( {1 - \frac{\Delta }{a}} \right).$Учитывая равенство ${{P}_{l}}\left( a \right) = P$, из первого соотношения в (20) находим ${\Delta \mathord{\left/ {\vphantom {\Delta a}} \right. \kern-0em} a}$ = = $\cos \left( {{{2\pi P} \mathord{\left/ {\vphantom {{2\pi P} D}} \right. \kern-0em} D}} \right)$ и после подстановки в (20) получаем и
(21)
${{с}_{l}}\left( a \right) = \frac{D}{{\pi a}}\sin \frac{{2\pi P}}{D},\quad {{b}_{l}}(a) = \frac{{2D}}{{\pi a\omega }}{{\sin }^{2}}\frac{{\pi P}}{D}.$Разность $a - \Delta $ определяет величину внедрения инструмента в обрабатываемую среду за один цикл колебаний. Очевидно, что средняя скорость процесса ${v} = {{\left( {a - \Delta } \right)\omega } \mathord{\left/ {\vphantom {{\left( {a - \Delta } \right)\omega } {2\pi }}} \right. \kern-0em} {2\pi }}$ и с помощью первого соотношения в (20) находим
В реальных процессах и машинах $P \ll D$. Поэтому, выражения (21) и (22) с высокой точностью можно представить в виде
(23)
${{с}_{l}}\left( a \right) = \frac{{2Р}}{a},\quad {{b}_{l}}(a) = \frac{{2\pi {{P}^{2}}}}{{a\omega D}},\quad {v} = \pi a\omega {{\left( {\frac{P}{D}} \right)}^{2}}.$Нелинейная динамика машины. Перейдем к построению динамических характеристик машины в рабочем режиме. Выше было показано, что все решения, полученные для режимов холостого хода, остаются справедливыми после замены величин c и b на величины $\left( {c + {{с}_{l}}} \right)$ и $\left( {b + {{b}_{l}}} \right)$. При питании возбудителя от источника тока по аналогии с (6) с учетом (23) находим
(24)
$а = \frac{{BlI}}{{{{{\sqrt {{{{\left( {c + \frac{{2Р}}{a} - m{{\omega }^{2}}} \right)}}^{2}} + {{\omega }^{2}}{{{\left( {b + \frac{{2\pi {{P}^{2}}}}{{a\omega D}}} \right)}}^{2}}} }}^{}}}}.$Соотношение (24) позволяет сразу найти связь резонансных частот и амплитуд с параметрами системы. Равенство нулю первого слагаемого в знаменателе (24) определяет значения резонансных частоты и амплитуды
(25)
${{a}_{{rl}}} = - \frac{{2Р}}{{c[1 - {{{\left( {{{\omega }_{{rl}}}{\text{/}}{{\omega }_{0}}} \right)}}^{2}}]}},\quad {{а}_{{rl}}} = \left( {BlI - \frac{{2\pi {{P}^{2}}}}{D}} \right){{\left( {b{{\omega }_{{rl}}}} \right)}^{{ - 1}}}.$Первое из равенств (25) является уравнением скелетной кривой амплитудно-частотной характеристики, второе – линии предельных амплитуд [6]. Решение уравнений (25) позволяет найти величины резонансных частот и амплитуд. Из (25) видно, что резонансы возникают в области частот ${{\omega }_{{rl}}} > {{\omega }_{0}}$ при условии
Из соотношений (25) следует, что резонансная частота в рабочем режиме ${{\omega }_{{rl}}} > {{\omega }_{0}}$, а резонансная амплитуда ${{a}_{{rl}}} \to \infty $ при $b \to 0$.
Уравнение (24) относительно неизвестной амплитуды колебаний рабочего органа приводится к следующему квадратному уравнению
Учитывая, что диссипация в системе существенно влияет на параметры колебаний только в окрестности резонанса, для получения ясной физической картины запишем решение уравнения, полагая b = 0
(27)
$a = \frac{{ - 2Р\left[ {1 \pm \sqrt {{{{\left( {\frac{{BlI}}{{2P}}} \right)}}^{2}} - {{{\left( {\frac{{\pi P}}{D}} \right)}}^{2}}} } \right]}}{{(c - m{{\omega }^{2}})}}.$Из выражения (27) следует, что неравенство (26) является необходимым условием существования виброударного рабочего процесса. Покажем, что вид амплитудно-частотной характеристики существенно зависит от соотношения сил возбуждения BLI и подачи P. Учитывая условие a > 0, из (27) находим, что рабочий процесс реализуется во всем частотном диапазоне при условии
На рис. 3 показаны амплитудно-частотные характеристики при различных значениях силы P. Рассмотренному случаю отвечает кривая P1. Резонансные частота и амплитуда определяются координатами точек пересечения скелетных кривых 1 и линии 2 предельных амплитуд. Ситуация радикально меняется, если увеличение силы P приводит к нарушению неравенства (28). В этом случае виброударные режимы возможны только в области $\omega > {{\omega }_{0}}$, и амплитудно-частотная характеристика имеет две ветви, разделенные скелетной кривой. На рис. 3 этому случаю отвечает кривая P2. По мере увеличения силы P подачи резонансная частота все более смещается в сторону более высоких частот. Используя соотношения (25), можно показать, что по мере приближения параметров системы к границе неравенства (26) резонансная частота ${{\omega }_{{rl}}} \to \infty $.
Показанная штриховой линией нижняя ветвь резонансной кривой P2 отвечает неустойчивым режимам. Вывод системы на устойчивые режимы с большой амплитудой можно осуществить либо внешним импульсом, либо затягиванием колебаний из области высоких частот, либо плавным увеличением силы P при заданной частоте $\omega \geqslant {{\omega }_{{rl}}}$. При этом следует иметь в виду опасность срыва колебаний при прохождении границы устойчивости.
Перейдем к построению и анализу амплитудно-частотных характеристик системы при питании от источника напряжения. По аналогии с (24), используя выражения (9), (10), (23), запишем уравнение относительно амплитуды колебаний инструмента
(29)
$a = \frac{{BlU}}{{\sqrt {({{R}^{2}} + {{\omega }^{2}}{{L}^{2}})\left\{ {{{{\left[ {\left( {c + \frac{{2Р}}{a} - {{\omega }^{2}}(m + kL)} \right)} \right]}}^{2}} + {{{\left[ {\omega b + kR + \frac{{2\pi {{P}^{2}}}}{{aD}})} \right]}}^{2}}} \right\}} }}.$и линии предельных амплитуд
(30)
${{a}_{{rl}}} = - \frac{{2Р}}{{c(1 - \omega _{{rl}}^{2}{\text{/}}\omega _{r}^{2})}},\quad {{a}_{{rl}}} = \left[ {\frac{{BlU}}{{{{R}_{{rl}}}}} - \frac{{2\pi {{P}^{2}}}}{D}} \right]{{\left( {{{\omega }_{{rl}}}b} \right)}^{{ - 1}}},$По физической сути соотношения (30) хорошо согласуются с равенствами (25). Из (30) следует, что скелетная кривая определена в частотной области $\omega > {{\omega }_{r}}$, а резонансные колебания можно реализовать при выполнении условия
Пренебрегая потерями энергии в механической (b = 0) и электрической (R = 0) частях системы, приведем уравнение (29) к квадратному
Решение уравнения имеет вид
(32)
$a = \frac{{ - 2P\left[ {1 \pm \sqrt {{{{\left( {\frac{{BlU}}{{2P{{R}_{\omega }}}}} \right)}}^{2}} - {{{\left( {\frac{{\pi P}}{D}} \right)}}^{2}}} } \right]}}{{c(1 - {{\omega }^{2}}{\text{/}}\omega _{r}^{2})}}.$Необходимое условие реализации виброударного режима имеет вид
При $\omega = {{\omega }_{{rl}}}$ неравенство (33), совпадающее с условием (31) реализации резонансных колебаний, ограничивает частотный диапазон существования виброударных процессов при заданной силе P
Заметим, что виброударный процесс в диапазоне $\omega < {{\omega }_{r}}$ возможен при
При выполнении (35) каждой частоте из диапазона (34) соответствует единственное значение амплитуды (32) колебаний. Отвечающая условию (35) амплитудно-частотная характеристика показана на рис. 4 (линия P1).
При нарушении условия (35) виброударные режимы реализуются только в частотном диапазоне ${{\omega }_{r}} < \omega < {{\omega }_{m}}$, где амплитудно-частотная характеристика (32) имеет вид петли, ветви которой разделенные скелетной кривой (30). На рис. 4 показаны амплитудно-частотные характеристики при различных значениях силы P. Тонкими линиями показаны скелетные кривые 1 и линии 2 предельных амплитуд. Штриховыми линиями показаны ветви неустойчивых режимов. По мере увеличения силы P3 > P2 > P1 частотный диапазон существования виброударного режима сужается и в пределе при P = Pm стягивается в точку, выделенную черным кружком. Графически в этой точке происходит касание скелетной кривой и линии предельных амплитуд, отвечающей предельному значению силы P.
В заключение отметим, что приведенные результаты позволяют провести полный расчет машины с электродинамическим возбудителем колебаний. Поскольку резонансные режимы являются наиболее эффективными [9] при работе вибрационной машины, понимание происходящих процессов позволяет найти лучшие пути реализации этих режимов.
Список литературы
Римский-Корсаков А.В. Электроакустика. М.: Связь, 1973. 272 с.
Харкевич А.А. Теория электроакустических преобразователей. Избранные труды. М.: Наука, 1973. Т. 1. С. 33.
Покровский В.В. Электродинамические вибровозбудители. Вибрации в технике. Справочник. М.: Машиностроение, 1981. Т. 4. С. 269.
Краснопольская Т.С. Автономное возбуждение механических колебаний электродинамическим вибратором // Прикладная механика. 1977. Т. 13. № 2. С. 108.
Astashev V.K., Babitsky V.I., Kolovsky M.Z. Dynamics and Control of Machines. Berlin: Springer, 2000. 235 p.
Бабицкий В.И. Теория виброударных систем. М.: Наука, 1967. 352 с.
Асташев В.К., Крупенин В.Л. Нелинейная динамика ультразвуковых технологических процессов. М.: МГУП им. Ивана Федорова, 2016. 372 с.
Пальмов В.А. Колебания упругопластических тел. М.: Наука, 1976. 328 с.
Антипов В.И., Асташев В.К. О принципах создания энергосберегающих вибрационных машин // Проблемы машиностроения и надежности машин. 2004. № 4. С. 3.
Дополнительные материалы отсутствуют.
Инструменты
Проблемы машиностроения и надежности машин