УДК 621.313.282.2
ОБОБЩЕННАЯ МОДЕЛЬ ДВУХКАТУШЕЧНОЙ СИНХРОННОЙ ЭЛЕКТРОМАГНИТНОЙ МАШИНЫ ДЛЯ ТЕХНОЛОГИЧЕСКИХ СИСТЕМ ВИБРОУДАРНОГО ДЕЙСТВИЯ
Нейман Людмила Андреевна1,
Нейман Владимир Юрьевич1,
1 Новосибирский государственный технический университет, Россия, 630073, г. Новосибирск, пр. Карла Маркса, 20.
Актуальность исследования обусловлена тем, что силовые электрические импульсные системы, созданные на основе линейных электромагнитных машин ударного и виброударного действия, являются перспективными для применения в технологиях по интенсификации процессов, связанных с поиском, добычей и переработкой георесурсов. В вопросах энергосбережения особый интерес вызывают линейные синхронные электромагнитные машины возвратно-поступательного движения ударного и виброударного действия. Частота ударов бойка в данных машинах синхронно связана с частотой источника промышленной сети переменного тока. К главным достоинствам оборудования, использующего линейные синхронные электромагнитные машины ударного и виброударного действия, относится снижение потребления электрической энергии при работе в резонансных режимах или близких к резонансным режимам. Исследование данных режимов связано с построением точных динамических моделей различных модификаций синхронных электромагнитных машин ударного и виброударного действия, работающих исключительно в переходных режимах.
Цель: разработать обобщенную математическую модель двухкатушечной синхронной машины с электромагнитным приводом для технологических систем виброударного действия, упрощающей процесс математического описания известных вариантов схем машин при проведении динамических расчетов.
Методы: дифференциальные уравнения электрического равновесия и механического взаимодействия, поступательно движущихся масс, уравнение Лагранжа второго рода, методы конечно-элементного моделирования двухмерного магнитного поля и структурного моделирования в программе Ма^аЬ ^тиНпк, сопоставление расчетных и экспериментальных характеристик.
Результаты. Предложена обобщенная расчетная схема построения многомассовой двухкатушечной синхронной электромагнитной машины, включающая известные модификации вариантов схем машин ударного действия. На основе расчетной схемы образована система дифференциальных уравнений обобщенной модели, позволяющая выполнять построение математических моделей, отражающих динамическое состояние известных вариантов схем синхронных ударных машин с электромагнитным приводом. Приведен пример реализации динамической модели четырехмассовой двухкатушечной синхронной электромагнитной машины с упругим реверсом бойка. Приведены результаты имитационного моделирования рабочего режима в виде диаграмм тока катушек, перемещения и скорости бойка, отражающие процессы включения и выхода на ква-зиустановившийся режим работы. На основе сопоставления значений интегральных характеристик показано, что расхождение результатов расчета и эксперимента составляет 4...6 %, что не превышает общепринятой погрешности для инженерных расчетов. Эффективность применения предложенной обобщенной модели заключается в упрощении математического описания вариантов схем синхронных машин с электромагнитным приводом и сокращении сроков реализации моделей.
Ключевые слова:
Электрическая импульсная система, линейная синхронная машина, синхронная частота, число степеней свободы, обобщенная модель, динамическая модель, моделирование.
Введение
В повышении эффективности интенсификации ряда технологических процессов и производств важную роль играют силовые электрические импульсные системы, способные создавать большие по величине ударные или виброударные воздействия на обрабатываемую среду [1-5].
Высокую долю среди них составляют силовые электромагнитные импульсные системы, содержащие в качестве основного исполнительного механизма ударный или виброударный электропривод на основе различных по конструкции линейных электромагнитных двигателей (ЛЭМД) [6-9]. Генерирование ударных нагрузок в ЛЭМД обеспечивается в результате возвратно-поступательного движения бойка. В самом общем случае подобные структуры, содержащие в своем составе ЛЭМД, устройства для преобразования
и передачи механической энергии и рабочий орган импульсной системы, образуют различные модификации электромагнитных машин ударного действия, которые находят все большее применение в прогрессивных технологиях по повышению интенсификации процессов поиска, добычи и переработки георесурсов [10-12].
Перспективным направлением считается применение электромагнитных машин ударного действия в приводах скважинных виброисточников, оказывающих силовое воздействие в обсадной колонне на уровне продуктивного пласта и повышающих интенсификацию процесса добычи нефти из малодебитных скважин [13, 14]. Их использование существенно повышает эффективность вибросейсмического метода без применения громоздкого наземного оборудования
[15].
DOI 10.18799/24131830/2021/1/3004
107
Установлена эффективность применения электромагнитных машин ударного действия различной мощности для зондирующего воздействия на грунт при сейсморазведке полезных ископаемых, что позволяет повысить точность повторения параметров ударных воздействий, упростить производство работ, обеспечить безопасность их проведения для обслуживающего персонала и окружающей среды [16-18].
Обоснована целесообразность использования линейного электропривода на основе электромагнитных машин для гидротехнического строительства морских стационарных платформ, с помощью которых ведется разработка полезных ископаемых на континентальных шельфах [19, 20].
Практика использования мобильных установок на основе электромагнитных машин ударного действия показала свою эффективность в горнодобывающей отрасли при вторичном разрушении горных пород, а также бурении ударным способом шпуров и скважин в твердых породах [21-23].
Значительным результатом следует считать серийное промышленное применение электромагнитных ударных машин в приводах ручного ударного и ударно-вращательного электроинструмента, широко используемого при монтаже в строительстве и в горном деле [24-26].
Развитие и расширение областей применения электромагнитных машин на основе ЛЭМД, прежде всего, связано с вопросами снижения энергопотребления и совершенствования методов их расчета [27, 28].
В вопросах энергосбережения среди известных вариантов схем электромагнитных машин особый интерес представляют синхронные электромагнитные машины ударного или виброударного действия (СЭМУД). В подобных машинах частота возвратно-поступательного движения бойка синхронно связана с частотой питающей промышленной сети переменного тока [29-31]. Основными преимуществами машин и оборудования на основе СЭМУД является возможность снижения потребления электрической энергии от промышленной сети переменного тока при работе в резонансных режимах или режимах близких к резонансным.
В процессе эволюционирования были созданы различные модификации вариантов схем СЭМУД, выполненные в основном по однокатушечным и двухкатушечным вариантам схем.
Разнообразие в вариантах схем, прежде всего, было вызвано необходимостью повышения их надежности, безопасности, улучшением эксплуатационных характеристик при передаче энергии ударом.
Наиболее широкую известность получили модификации вариантов схем двухкатушечных СЭМУД [29].
Существующие варианты схем имеют различную по сложности динамическую систему, способствующую оптимальному процессу формирования импульса силы и передачи его в деформируемую среду.
Сложность расчета подобных динамических систем заключается в прерывистом и многоступенчатом характере преобразования энергии. За счет импульс-
ного потребления электрической энергии работа выполняется исключительно в переходных режимах.
В цикличном процессе перемещения бойка преобладают динамические процессы, зависящие от многих факторов, таких как свойства деформируемой среды, связи в механической колебательной системе, свойства применяемых магнитных материалов, сил сопротивления движению бойка и т. д. Точность моделей для решения задач анализа и синтеза в существенной степени будет определяться точностью учета взаимосвязанных факторов.
Важным этапом расчетов является построение точных динамических моделей вариантов схем двух-катушечных СЭМУД. Учитывая их разнообразие, возникает целесообразность в построении обобщенной модели, отражающей динамическое состояние известных модификаций вариантов схем.
Цель исследования состоит в разработке обобщенной математической модели двухкатушечной синхронной машины с электромагнитным приводом для технологических систем виброударного действия, упрощающей процесс математического описания известных вариантов схем машин при проведении динамических расчетов.
Построение обобщенной математической модели
двухкатушечной СЭМУД
В общем случае в механической системе существующих вариантов схем СЭМУД возникают повторяющиеся динамические процессы, которые сопровождаются взаимным преобразованием потенциальной энергии упругих связей и кинетической энергии движущихся инерционных масс с последующей передачей энергии в деформируемую среду. Данные процессы возникают от действия вынуждающих сил электромагнитного поля, периодически возбуждаемого системой катушек.
Расчетная схема построения обобщенной модели многомассовой двухкатушечной СЭМУД, реализующая известные способы возвратно-поступательного движения бойка, представлена на рис. 1. Характер движения ударной массы бойка и управление катушками рабочего - К1 и обратного - К2 хода в обобщенной схеме устанавливается исходя из принятой модификации варианта схемы СЭМУД, согласно с заданным алгоритмом управления.
Расчетная схема (рис. 1) состоит из механической колебательной системы, имеющей пять степеней свободы. Система содержит боёк массой т1, который в процессе своего движения периодически взаимодействует с рабочим инструментом массой т2 и массивным инерционным преобразователем массой т3.
Рабочий инструмент и инерционный преобразователь посредством упругих связей связаны с магнито-проводом и закрепленными внутри катушками с общей массой т4. Также механическая колебательная система может содержать дополнительное устройство в форме кольца массой т5 для крепления СЭМУД и защиты от вибраций.
Применение дополнительного устройства в сочетании с промежуточным упругим звеном обеспечива-
ет защиту оператора от вибраций, передаваемых на корпус технологического оборудования. Величина усилия нажатия - обеспечивает устойчивый контакт динамической системы с рабочим инструментом при передаче энергии посредством удара при контакте с рабочей поверхностью деформируемой среды.
Работа механической колебательной системы сопровождается постоянными процессами поглощения энергии в упругих связях и потерями энергии за счет возникающих сил трения скольжения между соприкасающимися поверхностями.
Рассматриваемая система имеет пять степеней свободы (N=5). Выбираем в качестве обобщенных переменных координаты линейных перемещений х1... х5 центров масс ш1...ш5 (рис. 1).
тР34
/эм (1ц , Х1) _ /эм1 (гд , Х1) + /эм2 (г'д2 , Х1)
определяется координатой положения бойка х1, а также величиной намагничивающего тока катушки рабочего хода /Ц и катушки обратного хода /Ц.
Связь электрической и магнитной системы в расчетной схеме может быть учтена с помощью уравнений электромагнитного состояния, предусматривающих насыщение и потери энергии в стальном магни-топроводе на вихревые токи и гистерезис:
\ (0 = i r-
и2 (t) = i2r2 +
dt '
('ъ' X) dt ''
(1)
(2)
где м1(/), и2(/), /ь ¡1 - напряжение и токи в цепях питания катушек; ^1=/(/Ц1 ,х1), у2=/(/Ц2 ,х1) - потокосцепле-ние катушек; г1, г2 - активные сопротивления катушек рабочего и обратного хода.
Общий вид дифференциального уравнений, описывающих движение центров масс механической системы относительно обобщенных координат, может быть наиболее просто установлен через форму уравнений Лагранжа
±(дГ) дг
dtydxj dxj
(3)
///////////////////
Рис. 1. Расчетная схема обобщенной синхронной машины Fig. 1. Design diagram of a generalized synchronous machine
При выполнении расчетов полагаем, что массои-нерционные элементы механической системы рассматриваются как абсолютно твердые тела, массы которых соответственно равны m1_m5. Возникающие в системе ударные взаимодействия рассматриваются как мгновенные изменения скоростей соударяющихся тел. Силы вязкого трения пропорциональны скоростям инерционных масс. Вес бойка значительно меньше сил электромагнитного взаимодействия в механической системе.
В схеме на рис. 1 величина знакопеременной электромагнитной силы
где Т- кинетическая энергия системы; 0гп, QгФ, -обобщенные силы, соответствующие потенциальной энергии деформации упругих связей, энергиям рассеяния и внешних воздействий; х, х, - обобщенная координата и обобщенная скорость центров поступательно движущихся масс; /=1,2,3,...,N - число степеней свободы системы.
Для построения модели динамического состояния СЭМУД необходимо, обращаясь к расчетной схеме на рис. 1, выполнить операции дифференцирования левой части уравнения (3) и определить обобщенные силы системы, входящие в правую часть этого уравнения.
Рассматривая поступательные движения масс, кинетическая энергия системы
1
т = ~ X mi *<2 =
2 ,=1
1Г
= 2 ИЛ
2 + m2x22 + тъх\ + m4x24 +»75X5. (4)
Обобщенные силы, вызванные деформацией упругих связей
вП ^ = ~Ш 2 ^ (Х " Х3 )2 + к2 Х2 +
+к3 (х4 - х3)2 + к4 (х4 - х2)2 + к5 (х5 - х4)2], (5)
где П - потенциальная энергия; к1.к5 - коэффициенты жесткости упругих связей.
Обобщенные силы, соответствующие энергии рассеивания
X
еф=~=^ ■+^+
+Ъъ (х4 - х3)2 + ЪА (х4 - х2)2 + Ъ5 (х5 - Х4)2 ], (6)
где Ф - диссипативная функция (функция Рэлея); Ъ1...Ъ5 - коэффициенты линейного сопротивления, характеризующие диссипативные свойства упругих связей (вязкое трение).
Проекции обобщенной силы, вызванные трением скольжения и внешними воздействиями,
Qb — /эм , /, xi) /тЙ4 si
sign
ев - . dx2
в fTP24 Sign '
ев Г • dx3 в - /TP34sig^^;
QB — _ /эм (/ , / , Х1) + + ("/TP34 " /тР45 + /TP24 + ^4 )
QxB — /тр4! sign dХ- - FN,
dxx dt
dx.
sign
dt '
dt
m —-1+ь I —1 —31 + k (x - x) — 1 dt2 11 dt dt. lV 1 3'
— /Эм (/ V xi)- /тй 4 sign dXl •
dt
d x0 .. dx~.
m
2 dt2 2 dt
dx * dx*-.
+ b-2 - Ь -4--2 + kx^ -
dt dt
dx
(X4 - x2) —- /Ip24sign
d x •» i d^i dx.
m.
3 ^2
dt
dt dt
- - -
dx dx
dt dt
m
d x. .. I dx л dx-+ b
4 J.2 1
dt
7 I dx л dx0
-b4 — —
dx dx
b5 [ ^ - "^J + k3 (x4 - x3 ) + k4 (x4 - x2 ) -
~k5 (X5 - X4 ) = -/эм (/, / xi ) +
m.
d2 x,
—21 + b5
dt2
T - T. 1+ k5 (x - x4 ) —
,dx5 — /TP45Slgn - ^ •
(12)
(7)
где /тр14, /тр24, /тр34, /тр45 - силы сухого трения скольжения в соответствии с рис. 1.
Выражения (4)-(7) совместно с (3) позволяют образовать уравнения динамического состояния механической колебательной системы в виде следующих дифференциальных уравнений:
d х^ .. \ ^кХл
(8)
(9)
-к1 (xi - x3) - k3 (x4 - x3) — /тр34 signf; (10)
+ (^^Ф34 + /тР45 + /тР24 + /=14 ) sign ^ ^ ^
Уравнения (8)-(12) рассматривают только колебательные свойства механической системы с потерями энергии.
Дополняя уравнения электрического равновесия (1), (2) и динамического состояния механической системы (8)-(12) условиями, ограничивающими пределы возвратно-поступательного движения бойка, на границах которых возникают ударные взаимодействия, переходим к окончательному процессу формирования динамической модели с учетом рассматриваемой модификации двухкатушечной СЭМУД.
Полученную модель можно использовать для анализа рабочих процессов СЭМУД в переходных и ква-зиустановившихся режимах.
Процесс решения уравнений относительно принятого конструктивного исполнения СЭМУД должен сопровождаться предварительным расчетом функциональных зависимостей потокосцепления и электромагнитной силы у/'^хО, ^/^ХО, / эм^/^Х), /эм^Л'Цг Х0 от координаты смещения бойка и намагничивающего тока катушек. Значения указанных функциональных зависимостей в процессе моделирования определяются с помощью интерполяции табличных данных, полученных в результате решения полевой задачи в статике для различных комбинаций положения бойка и протекающего по катушкам намагничивающего тока с использованием современных методов численного расчета двухмерного магнитного поля [32-36].
Реализация динамической модели на примере
двухкатушечной СЭМУД с упругим реверсом бойка
Основные этапы использования обобщенной модели рассмотрим на примере реализации динамической модели четырехмассовой двухкатушечной СЭМУД с упругим реверсом бойка, представленной на рис. 2.
Впервые подобный вариант схемы был применен в конструкции ударного узла магнитобура СЦ-2 для бурения ударно-вращательным способом отверстий в бетоне и твердых породах, а также в ударных узлах буровых установок для разрушения горных пород [37].
Вариант схемы СЭМУД (рис. 2) содержит боек - 1, катушку рабочего хода - 2 и катушку обратного хода - 3. Обе катушки размещены внутри общего маг-нитопровода - 4, выполненного из листов электротехнической стали. При протекании тока по катушкам создается магнитное поле, с которым взаимодействует боек и совершает возвратно-поступательные движения. В процессе колебательных движений боек - 1 поочередно взаимодействует с упругим звеном - 5, выполненным в виде цилиндрической пружины, работающей на сжатие, и с торцевой поверхностью рабочего инструмента - 6. В результате этих
взаимодействий генерируется ударный импульс силы, передаваемый далее в деформируемую среду для совершения полезной работы.
Для гашения вертикальных вибраций применяется виброизолирующее устройство, состоящее из промежуточной массы - 7 и упругого звена, выполненного в виде цилиндрической пружины - 8, работающей на сжатие. Для устойчивой связи рабочего инструмента с деформируемой средой прилагается внешняя статическая сила нажатия передаваемая через промежуточную массу инструменту. Питание катушек рабочего и обратного хода выполняется импульсами тока по однополупериодной схеме выпрямления.
s+2
3
4 1
2 5-
Рис. 2. Вариант схемы двухкатушечной синхронной машины
Fig. 2. Variant of the scheme of a two-coil synchronous machine
За промежуток времени, равный одному периоду напряжения промышленной сети, осуществляется полный рабочий цикл машины. Для частоты сети f=50 Гц количество ударов бойка составляет яуд = 3000 уд/мин.
Рассматривая механическую колебательную систему рис. 2 и сопоставляя движение инерционных масс с расчетной схемой на рис. 1, следует полагать, что для реализации динамической модели следует принять: m3=0, k3=0, b3=0, x3=x4, х3=х4, =0, N=4.
Предварительно объединив на основе равенства обобщенных координат x3=x4 левые и правые части уравнений (10), (11), а также произведя некоторую
замену, получим уравнения, описывающие динамическое состояние четырехмассовой механической колебательной системы по варианту схемы на рис. 2:
d2 х.
d Xi .. I dx, dx, |т/ \
m ~0F+b1 d - tJ+ k ( x" x )-
/эм (^ ^ X1)- /ъ sign
dtj dt '
(13)
d x^ j ^—x^ .J I dx a dx.
m.
2 J,2 ' 2
+b-2 - b
+
dx.
+k2x2 - k4 (x4 - x2 ) = -/тр24 Sign ' (14)
d x„ , i dx, dx
m
4 dt2 1 k dt dt
- + b =1 -
dx a dx~
dt dt
-Ь5 ^ it - ( xi - x4 ) + k4 (x4 - x2 )-
-k5 ( x5 - x4 ) = - /эм (^ ^ x1) +
+ (/>45 + Л* + /I! ) S1
sign
dx4 dt
(15)
m.
d x5
~dF
£ - £J+k (x, - x4)-
r • dx5
(16)
Окончательный результат движения бойка необходимо рассматривать как процесс совместного действия вынуждающих электромагнитных сил и ударных импульсов сил, возникающих вследствие периодического взаимодействия бойка с рабочим инструментом.
В качестве ограничительных условий здесь выступают границы пределов поступательного перемещения бойка, гарантирующие его торможение, остановку и реверс в заданных границах. Данные пределы целесообразно рассматривать относительно промежутков рабочих воздушных зазоров 8х, 82 либо введенных фиксированных промежутков, определяемых расстояниями до ударного сечения Их, Н2, И3 и И4 (рис. 2).
Ограничения, вызванные торможением и остановкой бойка согласно рис. 2, могут быть установлены в виде следующей системы равенств:
x (t)-
0, при S+ - h4 - h2; h - h - , при о < x < h - h; h - h +|^21, при x > h - h,
(17)
Также важным условием для реализации модели является связь между кинетической энергией бойка, направленной на совершение полезной работы Тх и возвращаемой обратно в механическую систему Т2.
Связь этой энергии можно установить через соотношение скоростей ударной массы бойка до удара и после удара V в виде:
5
5
h
т = mvL (! _ к2); ^=mvL к2.
1 гу у от у 2 гу от
(18)
где кот=у1/у1 - коэффициент отскока.
Работа упругого звена при взаимодействии с бойком может быть описана следующей системой равенств, содержащих кусочно-линейные функции (рис. 2):
ЗП f0, при 0 < x < h _ h; dxx I k (X _ X), при x - h _ h
дФ dx,
о, при о < x < h _ h;
, f dx, dx,"] , ,
при х^з-Л,
) = '1 r1 +
dVi (v X1) ,
, (/) = i2 r2 +
dt
d^2 (^ X1) , dt ''
'l 'nl + 'ц1; '2 'п2 + 'ц2 ;
d2 X m —^ = A;
1 dt2
'2 dt2 d2 x
4 dt2 d2 x,
=
4 = C;
5 = D.
A =
iv1 /. . \ . dx1
1 + /эм ('«, X1 )_ /ТР14 S1gn
mv<
T
dt
dx
если —1 > 0, при 0 < x < h _ h; dt
dx
если —1 > 0, при x - h _ h; dt
/эм ('«,xi)_ ^ s1gn
dt
dx
если —1 < 0, при 0 < x < h _ h; dt
TTI
dt
C =
(19)
(20)
+ - /эм , , x j +
+ (^ _ /три ) S1gn ^
при 0 < x < К, ~ h2;
_(Лрм _ + /Тр45 ) S1gn
при x - h _ h;
dx,
dr4 dt
где
Окончательно принимая во внимание уравнения (1), (2), (13)-(16), с учетом установленных условий и ограничений (17)-(20), математическая модель варианта схемы СЭМУД на рис. 2 может быть представлена в виде системы:
D = -15+ /тр45 sign -¿~Fn>
X, = к,(х, -x.) + b, I —- - —i-1; = к. x. +b.
1 1 v 1 4' 1 а л 2 22 2 м
dt dt
dt
(21)
'5 dt2
Установленные в системе (21) переменные коэффициенты:
W1V1, Ш2У2 - импульсы сил при ударе; ти - время действия импульса; V2 - скорость приобретенная рабочим инструментом; /п, /п - токи потерь, обусловленные процессами в стальном магнитопроводе.
Полученная модель в виде системы дифференциальных уравнений (21) позволяет производить всесторонний анализ динамических процессов, происходящих в двухкатушечном СЭМУД с упругим реверсом бойка.
Решение системы (21) выполнялось в соответствии с существующими рекомендациями [38, 39] в программе Matlab Simulink.
При выполнении динамических расчетов использовались параметры модернизированного промышленного варианта СЭМУД «ИЭ-4207»: mj=0,394 кг; m2=0,32 кг; m4=3,9 кг; m5=0,2 кг А1=120^103 H/м; ¿1=8 Hx/м; А2=12^106 H/м; ¿2=0; ¿4=16^03 H/м; ¿4=35 нх/м;/тр14=3 Н;/Гр24=8 Н;/Гр45=5 Н; Fn=35 Н.
Для уменьшения потерь магнитопровод выполнен в виде пакетов из электротехнической стали 1212 толщиной 0,35 мм. Боек цельный из стали 40ХС.
Катушки выполнены медным проводником круглого сечения с использованием каркасной намотки. Параметры катушки рабочего хода: w1=1290, г1=12,9 Ом. Параметры катушки обратного хода: w2=1540, r2=23,2 Ом. Обе катушки независимо включены на переменное напряжение U=220 B через неуправляемые вентили и поочередно питаются от полуволн напряжения разной полярности.
Предварительный расчет функциональных зависимостей потокосцеплений и электромагнитных сил
V/uSO, У2=/Л,§2), /эм1=ЛМ0, /эм2=/лА) от величины воздушного рабочего зазора и тока катушек выполнялся в программе FEMM (Finite Element Method Magnetics) [40].
u
Полученные по результатам полевого расчета статических параметров трехмерные графики зависимости у1=/(/1,81) и /эм1=/(/1,81) для катушки рабочего хода, содержащие по 342 опорные точки статических параметров, представлены на рис. 3, 4.
Расчет составляющих токов /п1 и /п2 выполнялся по методике, изложенной в работе [41].
Результаты имитационного моделирования в виде диаграмм тока катушек, перемещения и скорости бойка, отражающие процесс включения СЭМУД и выхода на квазиустановившийся режим работы, рассмотрены на рис. 5-8.
Vi, Вб 4
0 -10
U, А
0 -10
i\, А
30 ' 0 &1, мм
Рис. 4. Трехмерный график расчета электромагнитного
усилия fz,Ml=f(il,8i) Fig. 4. Three-dimensional graph for calculating the electromagnetic force fem1=f(i1,S1)
Диаграммы характеризуются устойчивым периодическим режимом работы двухкатушечной СЭМУД с начальной скоростью при отскоке бойка Vi=0,5 м/с и коэффициентом отскока кот=0,11.
А
-1,0 -2,0 -3,0 -4,0 -5,0 -6,0 -7,0
0
Рис. 5. Fig. 5.
0,1 0,2 0,3 0,4 0,5 0,6 0,7 c Диаграмма тока i1(t) катушки рабочего хода Current diagram of the travel coil i1(t)
h
А 4,0 3,0 2,0 1,0 0 -1,0
0
0,1 0,2 0,3 0,4 0,5 0,6 0,7
30 s ^ 0
&i,MM
Рис. 3. Трехмерный график расчета потокосцепления W1=f(i1,S1)
Fig. 3. Three-dimensional graph for calculating flux linkage у1=f(i1,Sl)
/эм1,Н
750 500 250
Рис. 6. Диаграмма тока i2(t) катушки обратного хода Fig. 6. Current diagram of the reverse coil i2(t)
x1
мм 20 15 10 5 0
-5 0
Рис. 7. Fig. 7.
м/с 2,0 0
-2,0 -4,0 [-6,0
0,1 0,2 0,3 0,4 0,5 0,6 Диаграмма перемещения бойка x1(t) Movement diagram of the striker x1(t)
0,7
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 Рис. 8. Диаграмма скорости бойка v1(t) Fig. 8. Striker speed chart v1(t)
0
t
3
2
1
t
c
t
c
c
Также на рис. 5-8 отражена реакция модели СЭМУД по восстановлению колебательных свойств вследствие внезапного изменения упругих свойств деформируемой среды на интервале времени /=0,4...0,6 с.
С увеличением упругих свойств деформируемой среды на указанном интервале времени увеличивается начальная скорость отскока бойка Vi=2,0 м/с и коэффициент отскока кот=0,43. Заметно повышается амплитуда колебания бойка (рис. 7) при сохранении устойчивого режима работы модели. Максимальная скорость перемещения бойка (рис. 8) практически не изменяется.
Также наблюдается уменьшение амплитуды тока катушки рабочего хода (рис. 5) и катушки обратного хода (рис. 6), что вызывает снижение мощности, потребляемой от источника питания. Данное уменьшение объясняется тем, что с увеличением коэффициента отскока возрастает кинетическая энергия, возвращаемая в механическую систему при обратном ходе бойка. В то время как кинетическая энергия, передаваемая в деформируемую среду, снижается.
Таким образом, частичное снижение полезной мощности и мощности, потребляемой от источника, не нарушают колебательных свойств механической системы. Рабочий процесс характеризуются устойчивым режимом работы без применения специальных средств регулирования.
Таблица. Сравнение результатов Table. Comparison of results
Оценивание точности построения модели выполнялось сравнением интегральных значений показателей энергии удара - Ain, действующего значения тока на входе - I, активной мощности - P1, полезной мощности на выходе машины - P2 и КПД, полученных в ходе численного расчета с показателями, найденными с помощью физической модели СЭМУД с упругим реверсом бойка (рис. 2). Сравнение результатов рас-
СПИСОК ЛИТЕРАТУРЫ
1. Ашавский А.М., Вольперт А.Я., Шнейнбаум В.С. Силовые импульсные системы. - М.: Машиностроение, 1978. - 200 с.
2. Абрамов А.Д., Тюнюкова Т.К., Ижбулдин Е.А. Виброударные технологии для строительно-дорожных машин // Мир транспорта. - 2016. - Т. 14. - № 4 (65). - С. 62-69.
3. Абрамов А.Д., Галай М.С. Перспективы применения машин ударного действия при упрочнении сварных швов // Фундаментальные и прикладные вопросы горных наук. - 2016. - Т. 2. -№ 3. - С. 7-10.
4. Moshkin V.I. Pulsing linear electromagnetic motor mechanical work within the desaturated magnetic system // International Conference on Industrial Engineering, Applications and Manufacturing, ICIEAM, 2017. - P. 8076319.
5. Нейман В.Ю. Интегрированные линейные электромагнитные двигатели для импульсных технологий // Электротехника. -2003. - № 9. - С. 25-30.
чета и физического эксперимента представлено в таблице.
Расхождение результатов расчета (таблица) с результатами, полученными с помощью физической модели модернизированного варианта схемы СЭМУД с упругим реверсом бойка, составляет 4.6 %.
По результатам моделирования установлено, что точность в определении выходных характеристик моделей при выполнении динамических расчетов в значительной степени определяется точностью установки входных параметров моделей. В особенности это относится к точности выбора числа степеней свободы механической системы, обобщенных сил, связанных с энергией рассеивания, а также сил, вызванных трением скольжения и внешними воздействиями. Как правило, это параметры на входе модели, которыми разработчики чаще всего пренебрегают или пытаются существенно их упростить.
Заключение
Предложена обобщенная расчетная схема построения многомассовой двухкатушечной СЭМУД, включающая известные варианты схем машин.
На основе расчетной схемы образована система дифференциальных уравнений, позволяющая при установленных ограничениях относительно существующих модификаций вариантов схем выполнять построение математических моделей двухкатушеч-ных СЭМУД, характеризующихся многосторонними связями в механической системе, существенной нелинейностью электромагнитных тяговых характеристик и потокосцеплений, а также наличием ударной нагрузки и потерями энергии в электрической, магнитной и механической системе.
На примере четырехмассовой конструкции двух-катушечной СЭМУД с упругим реверсом бойка рассмотрен пример реализации математической модели для анализа динамических характеристик. Расхождение расчетных значений интегральных показателей с экспериментальными данными, полученными на физической модели, не превышает погрешности 4.6 %, допустимой для инженерных расчетов.
Эффективность предложенной обобщенной модели заключается в упрощении процесса математического описания вариантов схем СЭМУД при сокращении сроков реализации моделей.
6. Design and analysis of a new down-the-hole electromagnetic hammer driven by tube linear motor / T. Wu, Y. Tang, S. Tang, Y. Li, W. He, E. Chen // IET Electric Power Applications. -2017. - V. 11 (9). - P. 1558-1565.
7. Electric converters of electromagnetic strike machine with battery power / K.M. Usanov, A.V. Volgin, V.A. Kargin, A.P. Moiseev, E.A. Chetverikov // IOP Conference Series: Materials Science and Engineering. - 2018. - V. 327. - P. 052032.
8. Effect of vibroimpact treatment on the physical and mechanical properties of the surface layer of welded joints in rails / V.A. Kargin, L.B. Tikhomirova, A.D. Abramov, M.S. Galai // Welding International. - 2014. - V. 28 (3). - P. 245-247.
9. Neyman L.A., Neyman V.Yu., Shabanov A.S. Vibration dynamics of an electromagnetic drive with a half-period rectifier // The 18 international conference of young specialists on micro/nanotechnologies and electron devices: Proc. - Altai, Erlagol, 29 June - 3 July 2017. - P. 503-506.
Интегральный показатель Integral indicator Дж (J) I, A (a) P1, Вт (W) P2, Вт (W) КПД Efficiency, %
Расчет Calculation 4,21 3,54 492 211 41
Эксперимент Experiment 4,18 3,6 536 209 39
10. Yedygenov E.K., Lyashkov V. Electromagnetic rock breaker for non-explosive rocks breaking // International Symposium on planning of mining and selection of equipment of MPES 2011. -Almata, 2011. - Р. 1002-1012.
11. Павлов В.Е. Исследование режимов работы электромагнитного молота // Вестник Иркутского государственного технического университета. - 2017. - Т. 21. - № 12 (131). - С. 164-173.
12. Угаров Г.Г., Нейман В.Ю. Анализ показателей электромагнитных ударных машин // Физико-технические проблемы разработки полезных ископаемых. - 1996. - № 2. - С. 72-80.
13. Экспериментальные исследования эффективности излучения от скважинного импульсного виброисточника / Б.Ф. Симонов, В.Н. Опарин, А.О. Кордубайло, В.И. Востриков // Горный информационно-аналитический бюллетень (научно-технический журнал). - 2019. - № 8. - С. 180-189.
14. Симонов Б.Ф., Нейман В.Ю., Шабанов А.С. Импульсный линейный электромагнитный привод для скважинного виброисточника // Физико-технические проблемы разработки полезных ископаемых. - 2017. - № 1. - С. 118-126.
15. Simonov B.F., Neyman V.Y., Shabanov A.S. New conception of an electromagnetic drive for a vibration source in hole // The 18 international conference of young specialists on micro/ nanotechnologies and electron devices: Proc. - Altai, Erlagol, 29 June - 3 July 2017. - P. 507-510.
16. Ивашин В.В., Певчев В.П. Определение перемещений излучателя импульсного электромагнитного источника сейсмических сигналов // Физико-технические проблемы разработки полезных ископаемых. - 2018. - № 4. - С. 51-59.
17. Haines S.S. Design and application of an electromagnetic vibrator seismic source // Journal of Environmental and Engineering Geophysics. - 2006. - V. 11 (1). - P. 9-15.
18. Sadovskii V.M., Sadovskaya O.V., Detkov V.A. Analysis of elastic waves generated in frozen grounds by means of the electromagnetic pulse source «Yenisei» // IOP Conference Series: Earth and Environmental Science. - 2018. - V. 193 (1). -P. 012058.
19. Электромагнитные молоты для морского нефтегазопромысло-вого строительства / Б.Ф. Симонов, В.С. Вовк, Ю.В. Погар-ский, А.И. Кадышев // PROНЕФТЬ. Профессионально о нефти. - 2019. - № 2 (12). - С. 59-65.
20. Ряшенцев Н.П., Симонов Б.Ф., Ильинский Н.К. Электромагнитный молот для погружения свай // Газовая промышленность. - 1984. - № 11. - С. 11-15.
21. Манжесов В.К., Лукутина Н.А., Невенчанная Т.О. Динамика и синтез электромагнитных генераторов силовых импульсов. -Фрунзе: Изд-во Илим, 1985. - 185 с.
22. Едыгенов Е.К. Горные машины с линейным электромагнитным двигателем // Современные проблемы теории машин. -2019. - № 8. - С. 81-84.
23. Едыгенов Е.К., Васин К.А. Результаты испытаний электромагнитного молота для безвзрывного разрушения горных пород // Горный информационно-аналитический бюллетень (научно-технический журнал). - 2020. - № 5. - С. 80-90.
24. Ряшенцев Н.П., Тимошенко Е.М., Фролов А.В. Теория, расчет и конструирование электромагнитных машин ударного действия. - Новосибирск: Наука, 1970. - 259 с.
25. Ижбулдин Е.А., Абрамов А.Д. Ручной электрический ударный инструмент для реализации виброударных технологий в транспортном машиностроении и строительстве // Вестник Иркутского государственного технического университета. -2017. - Т. 21. - № 1 (120). - С. 32-39.
26. Strike action electromagnetic machine for immersion of rod elements into ground / K.M. Usanov, A.V. Volgin, E.A. Chetverikov, V.A. Kargin, A.P. Moiseev, Z.I. Ivanova // IOP Conference Series: Earth and Environmental Science. - 2017. -V. 87. - P. 032050.
27. Угаров Г.Г., Мошкин В.И. Перспективы развития силовых электромагнитных импульсных систем // Вестник Курганского государственного университета. - 2013. - № 2. - C. 88-90.
28. Гринченков В.П., Большенко И.А., Большенко А.В. Электромагнитный привод с низким энергопотреблением // Известия высших учебных заведений. Электромеханика. - 2015. -№ 5 (541). - С. 50-53.
29. Нейман Л.А., Нейман В.Ю. Линейные синхронные электромагнитные машины для низкочастотных ударных технологий // Электротехника. - 2014. - № 12. - С. 45-49.
30. Zhao, J.H., Zhang, X.F., Zhang, J.H. Field and thrust analysis of tubular permanent magnet linear synchronous motor // Electric Machines and Control - 2010. - V. 22 (3). - P. 60-63.
31. Neyman L.A., Neyman V.Yu., Obukhov K.A. New method of the synchronous vibratory electromagnetic machine mechatronic module control // The 18 international conference of young specialists on micro/nanotechnologies and electron devices: Proc. - Altai, Erlagol, 29 June - 3 July 2017. - P. 516-519.
32. Tatevosyan A.S., Tatevosyan A.A., Zaharova N.V. Calculation of non-stationary magnetic field of the polarized electromagnet with the external attracted anchor // Journal of Physics: Conference Series. - 2018. - V. 1050 (1). - P. 012086.
33. 2D Finite Element Model for computing the electric field strength-rms generated by overhead power lines / E. Lunca, S. Vornicu, A. Salceanu, O. Bejenaru // Journal of Physics: Conference Series. - 2018. - V. 1065. - P. 1-4.
34. Shevchenko V.P., Babiychuk O.B., Boltenkov V.O. Study of current transformers magnetic field by method final elements using the FEMM software complex // Applied aspects of information technology. - 2019. - V. 2 (4). - P. 317-327.
35. Cute V.I. System of automated formation of calculating models of electric machines for FEMM software environment // Technical electrodynamics. - 2018. - № 4. - P. 74-78.
36. A high torque density outer rotor claw pole stator permanent magnet synchronous motor / J. Liang, A. Parsapour, E. Cosoroaba, M. Wu, I. Boldea, B. Fahimi // IEEE Transportation Electrification Conference and Expo (ITEC). - 2018. - P. 389-393.
37. Ряшенцев Н.П., Алабужев П.М., Никишин Н.И. Ручные электрические машины ударного действия. - М.: Недра, 1970. -198 с.
38. Черных И.В. Моделирование электротехнических устройств в MATLAB, SimPowerSystems и Simulink. - М.: ДМК Пресс; СПб.: Питер, 2008. - 288 с.
39. Kiyakli A., Hamit S. Modeling of an electric vehicle with MATLAB/Simulink // International journal of automotive science and technology. - 2018. - V. 2 (4). - P. 9-15.
40. Буль О.Б. Методы расчета магнитных систем электрических аппаратов: Магнитные цепи, поля и программа FEMM. - М.: ИЦ «Академия», 2005. - 336 с.
41. Нейман Л.А., Нейман В.Ю. Моделирование динамических процессов в электромагнитных преобразователях энергии для систем генерирования силовых воздействий и низкочастотных вибраций // Известия Томского политехнического университета. Инжиниринг георесурсов. - 2015. - Т. 326. - № 4. -С. 154-162.
Поступила 31.08.2020 г.
Информация об авторах
Нейман Л.А., доктор технических наук, профессор кафедры электротехнических комплексов Новосибирского государственного технического университета.
Нейман В.Ю., доктор технических наук, профессор, заведующий кафедрой теоретических основ электротехники Новосибирского государственного технического университета.
UDC 621.313.282.2
GENERALIZED MODEL OF A TWO-COIL SYNCHRONOUS ELECTROMAGNETIC MACHINE FOR VIBRATORY IMPACT TECHNOLOGICAL SYSTEMS
Lyudmila A. Neyman1,
Vladimir Yu. Neyman1,
1 Novosibirsk State Technical University, 20, Karl Marks avenue, Novosibirsk, 630073, Russia.
Relevance of the research is explained by the perspective of electric pulse systems based on impact and vibratory impact linear electromagnetic machines application in geological exploration, mining operations and minerals processing technologies. Impact synchronous electromagnetic machines are useful for energy saving. Impact force pulse frequency of these machines is equal to or multiple of the power source frequency. The main advantage of devices and equipment based on impact and vibratory impact synchronous electromagnetic machines is the capability of power consumption decrease in near-resonant modes. These modes study is important for the design of exact dynamic models of different variants of impact synchronous electromagnetic machines operating only in transient modes.
The aim of the research is to develop the generalized mathematical model of a two-coil impact synchronous machine with an electromagnetic drive for vibratory impact technological systems. The model simplifies the mathematical description of known variants machine schemes when carrying out dynamic calculations.
The methods used in the research include differential equations of the electrical equilibrium and mechanical interaction of directly moving masses, Lagrange's equation of the second kind, finite-element methods of magnetic field modeling, methods of structured modeling in Matlab Simulink, design and experimental characteristics comparison.
Results. The authors have proposed the generalized design model of the multi-mass two-coil impact synchronous electromagnetic machine for its known and new variants. The differential equations system were derived for this model. The generalized model is the basis of dynamic models of different variants of impact synchronous electromagnetic machines. The paper considers as an example the dynamic model of the four-mass two-coil impact synchronous electromagnetic machine with the striker elastic reverse. The simulation results stated for the operating mode include the oscillograms of the current in the coils, the striker motion and velocity illustrating switching on and reaching the quasi-stationary operating mode. It follows from integral characteristics comparison that the difference between the design and experimental data is 4...6 % and it is no more than common engineering design error. The proposed generalized model helps effectively to simplify the mathematical description of the impact synchronous electromagnetic motor and decrease the time interval for model implementation.
Key words: Electrical impulse system, linear synchronous machine, synchronous frequency, number of degrees of freedom, generalized model, dynamic model, simulation.
REFERENCES
1. Ashavskiy A.M., Volpert A.Ya., Shneynbaum V.S. Silovye impulsnye sistemy [Power impulse systems]. Moscow, Mashinostroenie Publ., 1978. 200 p.
2. Abramov A.D., Tyunyukova T.K., Izhbuldin E.A. Vibroudarnye tekhnologii dlya stroitelno-dorozhnykh mashin [Vibration shock technologies for road-building machines]. Mir transporta, 2016, vol. 14, no. 4 (65), pp. 62-69.
3. Abramov A.D., Galay M.S. Perspektivy primeneniya mashin udarnogo deystviya pri uprochnenii svarnykh shvov [Prospects for the use of impact machines for strengthening welded seams]. Fundamentalnye i prikladnye voprosy gornykh nauk, 2016, vol. 2, no. 3, pp. 7-10.
4. Moshkin V.I. Pulsing linear electromagnetic motor mechanical work within the desaturated magnetic system. International Conference on Industrial Engineering, Applications and Manufacturing, ICIEAM. 2017. pp. 8076319.
5. Neiman V.Y. Integrated linear electromagnetic motors for pulsed technologies. Russian Electrical Engineering, 2003, vol. 74, no. 9, pp. 30-35. In Rus.
6. Wu T., Tang Y., Tang S., Li Y., He W., Chen E. Design and analysis of a new down-the-hole electromagnetic hammer driven by tube linear motor. IET Electric Power Applications, 2017, vol. 11, no. 9, pp. 1558-1565.
7. Usanov K.M., Volgin A.V., Kargin V.A., Moiseev A.P., Chetverikov E.A. Electric converters of electromagnetic strike machine with battery power. IOP Conference Series: Materials Science and Engineering, 2018, vol. 327, pp. 052032.
8. Kargin V.A., Tikhomirova L.B., Abramov A.D., Galai M.S. Effect of vibroimpact treatment on the physical and mechanical properties of the surface layer of welded joints in rails. Welding International, 2014, vol. 28, no. 3, pp. 245-247.
9. Neyman L.A., Neyman V.Yu., Shabanov A.S. Vibration dynamics of an electromagnetic drive with a half-period rectifier. The 18 international conference of young specialists on micro/nanotechnologies and electron devices: proc. Altai, Erlagol, 29 June - 3 July 2017. pp. 503-506.
10. Edygenov E.K., Lyashkov V. Electromagnetic rock breaker for non-explosive rocks breaking. International Symposium on planning of mining and selection of equipment of MPES 2011. Almata, 2011. pp. 1002-1012.
11. Pavlov V.E. Issledovanie rezhimov raboty elektromagnitnogo molota [Research of operating modes of the electromagnetic hammer]. Vestnik Irkutskogo gosudarstvennogo tekhnicheskogo universiteta, 2017, vol. 21, no. 12 (131), pp. 164-173.
12. Ugarov G.G., Neiman V.Y. Evaluation of operating conditions for electromagnetic impactors. Journal of Mining Science, 1996, vol. 32, no. 4, pp. 305-312. In Rus.
13. Simonov B.F., Oparin V.N., Kordubaylo A.O., Vostrikov V.I. Eksperimentalnye issledovaniya effektivnosti izlucheniya ot skvazhinnogo impulsnogo vibroistochnika [Experimental studies of the efficiency of radiation from a borehole pulsed vibration source]. Gorny informatsionno-analiticheskiy byulleten (nauchno-tekhnicheskiy zhurnal), 2019, no. 8, pp. 180-189.
14. Simonov B.F., Neiman V.Y., Shabanov A.S. Pulsed linear solenoid actuator for deep-well vibration source. Journal of Mining Science, 2017, vol. 51, no. 1, pp. 117-125. In Rus.
15. Simonov B.F., Neyman V.Y., Shabanov A.S. New conception of an electromagnetic drive for a vibration source in hole. The 18 international conference of young specialists on micro/nanotechnologies and electron devices: proc. Altai, Erlagol, 29 June - 3 July 2017. pp. 507-510.
16. Ivashin V.V., Pevchev V.P. Determining displacements of pulsed electromagnetic source of seismic signals. Journal of Mining Science, 2018, vol. 54, no. 4, pp. 575-581. In Rus.
17. Haines S.S. Design and application of an electromagnetic vibrator seismic source. Journal of Environmental and Engineering Geophysics, 2006, vol. 11 (1), pp. 9-15.
18. Sadovskii V.M., Sadovskaya O.V., Detkov V.A. Analysis of elastic waves generated in frozen grounds by means of the electromagnetic pulse source «Yenisei». IOP Conference Series: Earth and Environmental Science, 2018, vol. 193 (1), pp. 012058.
19. Simonov B.F., Vovk V.S., Pogarskiy Yu.V., Kadyshev A.I. Elektromagnitnye moloty dlya morskogo neftegazopromyslovogo stroitelstva [Electromagnetic hammers for offshore oil and gas field construction]. PRONEFT. Professionalno o nefti, 2019, no. 2 (12), pp. 59-65.
20. Ryashentsev N.P., Simonov BF., Ilinskiy N.K. Elektromagnitny molot dlya pogruzheniya svay [Electromagnetic hammer for driving piles]. Gazovaya promyshlennost, 1984, no. 11, pp. 11-15.
21. Manzhesov V.K., Lukutina N.A., Nevenchannaya T.O. Dinamika i sintez elektromagnitnykh generatorov silovykh impulsov [Dynamics and synthesis of electromagnetic power pulse generators]. Frunze, Ilim Publ., 1985. 185 p.
22. Edygenov E.K. Gornye mashiny s lineynym elektromagnitnym dvigatelem [Mining machines with linear electromagnetic motor]. Sovremennye problemy teorii mashin, 2019, no. 8, pp. 81-84.
23. Edygenov E.K., Vasin K.A. Rezultaty ispytaniy elektromagnitnogo molota dlya bezvzryvnogo razrusheniya gornykh porod [Test results of an electromagnetic hammer for non-explosive destruction of rocks]. Gorny informatsionno-analiticheskiy byulleten (nauchno-tekhnicheskiy zhurnal), 2020, no. 5, pp. 80-90.
24. Ryashentsev N.P., Timoshenko E.M., Frolov A.V. Teoriya, raschet i konstruirovanie elektromagnitnykh mashin udarnogo deystviya [Theory, calculation and design of electromagnetic impact machines]. Novosibirsk, Nauka Publ., 1970. 259 p.
25. Izhbuldin E.A., Abramov A.D. Ruchnoy elektricheskiy udarny instrument dlya realizatsii vibroudarnykh tekhnologiy v transportnom mashinostroenii i stroitelstve [Hand-held electric percussion tool for the implementation of vibration shock technologies in transport engineering and construction]. Vestnik Irkutskogo gosudarstvennogo tekhnicheskogo universiteta, 2017, vol. 21, no. 1 (120), pp. 32-39.
26. Usanov K.M., Volgin A.V., Chetverikov E.A., Kargin V.A., Moiseev A.P., Ivanova Z.I. Strike action electromagnetic machine for immersion of rod elements into ground. IOP Conference Series: Earth and Environmental Science, 2017, vol. 87, pp. 032050.
27. Ugarov G.G., Moshkin V.I. Perspektivy razvitiya silovykh elektromagnitnykh impulsnykh sistem [Prospects for the development of power electromagnetic pulse systems]. Vestnik Kurganskogo gosudarstvennogo universiteta, 2013, no. 2, pp. 88-90.
28. Grinchenkov V.P., Bolshenko I.A., Bolshenko A.V. Elektromagnitny privod s nizkim energopotrebleniem [Low power
consumption electromagnetic actuator]. Izvestiya vysshikh uchebnykh zavedeniy. Elektromekhanika, 2015, no. 5 (541), pp. 50-53.
29. Neiman L.A. Neiman V.Yu. Linear synchronous electromagnetic machines for lowfrequency impact technologies. Russian Electrical Engineering, 2014, vol. 85, no. 12, pp. 752-756. In Rus.
30. Zhao J.H., Zhang X.F., Zhang J.H. Field and thrust analysis of tubular permanent magnet linear synchronous motor. Electric Machines and Control, 2010, vol. 22 (3), pp. 60-63.
31. Neyman L.A., Neyman V.Yu., Obukhov K.A. New method of the synchronous vibratory electromagnetic machine mechatronic module control. The 18 international conference of young specialists on micro/nanotechnologies and electron device: proc. Altai, Erlagol, 29 June - 3 July 2017. pp. 516-519.
32. Tatevosyan A.S., Tatevosyan A.A., Zaharova N.V. Calculation of non-stationary magnetic field of the polarized electromagnet with the external attracted anchor. Journal of Physics: Conference Series, 2018, vol. 1050 (1), pp. 012086.
33. Lunca E., Vornicu S., Salceanu A., Bejenaru O. 2D Finite Element Model for computing the electric field strength-rms generated by overhead power lines. Journal of Physics: Conference Series, 2018, vol. 1065, pp. 1-4.
34. Shevchenko V.P., Babiychuk O.B., Boltenkov V.O. Study of current transformers magnetic field by method final elements using the FEMM software complex. Applied aspects of information technology, 2019, vol. 2 (4), pp. 317-327.
35. Cute V.I. System of automated formation of calculating models of electric machines for FEMM software environment. Technical electrodynamics, 2018, no. 4, pp. 74-78.
36. Liang J., Parsapour A., Cosoroaba E., Wu M., Boldea I., Fahimi B. A high torque density outer rotor claw pole stator permanent magnet synchronous motor. IEEE Transportation Electrification Conference and Expo (ITEC), 2018. pp. 389-393.
37. Ryashentsev N.P., Alabuzhev P.M., Nikishin N.I. Ruchnye elektricheskie mashiny udarnogo deystviya [Manual electric impact machines]. Moscow, Nedra Publ., 1970. 198 p.
38. Chernykh I.V. Modelirovanie elektrotekhnicheskikh ustroystv v MATLAB, SimPowerSystems i Simulink [Simulation of electrical devices in MATLAB, SimPowerSystems and Simulink]. Moscow, DMK Press; St. Petersburg, Piter Publ., 2008. 288 p.
39. Kiyakli A., Hamit S. Modeling of an electric vehicle with MATLAB/Simulink. International journal of automotive science and technology, 2018, vol. 2, no. 2, pp. 9-15.
40. Bul O.B. Metody rascheta magnitnykh sistem elektricheskikh apparatov: Magnitnye tsepi, polya i programma FEMM [Methods for calculating the magnetic systems of electrical devices: Magnetic circuits, fields and the FEMM program]. Moscow, Akademiya Publ., 2005. 336 p.
41. Neyman L.A., Neyman V.Yu. Simulation of dynamic processes in electromagnetic energy converters for force effects and lowvfrequency vibrations generation systems. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2015, vol. 326, no. 4, pp. 154-162. In Rus.
Received: 31 August 2020.
Information about the authors Lyudmila A. Neyman, Dr. Sc., professor, Novosibirsk State Technical University. Vladimir Yu. Neyman, Dr. Sc., professor, Novosibirsk State Technical University.