CC BY
71
УДК 629.423: 621.313
Е.С. Рябов, Б.Г. Любарский, Д.Ю. Зюзин, В.Л. Емельянов
К ВОПРОСУ ОПРЕДЕЛЕНИЯ ЭЛЕКТРОМАГНИТНОГО МОМЕНТА РЕАКТИВНОГО ИНДУКТОРНОГО ДВИГАТЕЛЯ С АКСИАЛЬНЫМ МАГНИТНЫМ ПОТОКОМ
Приведено математичну модель для визначення електромагнітного моменту індукторного реактивного двигуна з аксіальним магнітним потоком. Особливістю моделі є знаходження моменту шляхом інтеграції закону зміни магнітної провідності.
Приведена математическая модель для определения электромагнитного момента индукторного реактивного двигателя с аксиальным магнитным потоком. Особенностью модели является нахождение момента путем интеграции закона изменения магнитной проводимости.
ВВЕДЕНИЕ
Как известно, для «традиционных» электрических машин - постоянного тока, синхронных и асинхронных - вопрос проектирования в настоящее время сводится к следующему: по заданной мощности и частоте вращения выбирают основные геометрические параметры и электромагнитные нагрузки. Затем, используя апробированные математические модели, путём перебора определяют остальные геометрические и электрофизические параметры, которые обеспечивают требуемые рабочие и энергетические характеристики устройств. Такой подход к проектированию основан на большом опыте построения и их длительной эксплуатации, углубленных теоретических и экспериментальных исследованиях процессов, происходящих в электрических машинах.
Теория проектирования индукторных машин, к которым относиться реактивный индукторный двигатель с аксиальным магнитным потоком, находиться в стадии становления. Многочисленные публикации носят, как правило, теоретический характер, не дают полной уверенности в справедливости полученного результата [1]. В такой ситуации представляется целесообразным искать новые подходы к решению задачи синтеза индукторных машин.
Задача синтеза устройства может быть решена несколькими путями. Первый путь - это анализ (аналитический либо численный) всех возможных вариантов с последующим выбором требуемого или наиболее подходящего варианта. Однако такой подход требует больших временных затрат и ресурсов.
Вторым путём решение задач синтеза является оптимальное проектирование. В этом случае также решаются задачи анализа, но параметры задачи анализа на последующем шаге итерационного процесса однозначно определяются результатами решения задачи анализа на предыдущем шаге. Закон изменения параметров при переходе от одной итерации к другой заложен в оптимизационном методе. Критерием нахождения оптимального варианта является способность синтезируемого устройства удовлетворять функциональным требованиям с заранее установленными показателями или закономерностями.
Исходя их вышесказанного, задача синтеза реактивного индукторного двигателя с аксиальным магнитным потоком может быть сформулирована таким образом: необходимо создать устройство с заданными рабочими характеристиками (показателями) за счёт нахождения наилучшего сочетания конкретных геометрических и электромагнитных параметров уст-
ройства путём использованием процедуры оптимального проектирования.
Цель работы: разработать модель определения электромагнитного момента и магнитного напряжения индукторного двигатель с аксиальным магнитным потоком, которую можно использовать в качестве задачи анализа при оптимальном проектировании
МАТЕРИАЛ ИССЛЕДОВАНИЙ Мгновенное значение электромагнитного момента электрической машины может быть найдено из соотношения [2]
м = -
(1)
і=еотї
где дЖ' - изменении коэнергии системы, обусловленное изменением координаты ротора ду, і - токи в обмотках.
Применительно к индукторным машинам в [3] электромагнитный момент М2, реализуемый на участке зубцового деления, предложено выражать через производную от коэнергии Ж' зубцовой зоны на участке зубцового деления при условии постоянства магнитного напряжения и& зубцового слоя (рис.1):
Мг =■
(2)
С учётом того, что зубцовая зона фазы АИД образована 22 зубцовыми делениями, выражение для мгновенного значения электромагнитного момента фазы АИД имеет вид
ду
(3)
и& =еот1
Здесь 1 - число зубцов статора (ротора) по одну сторону фазного блока.
Суммарный момент на валу (в мгновенных значениях) равен
Мдв (Л У) = Е Мфк (^ у) к=1
(4)
2 - воздушный зазор, 3 - зубец ротора, 1^ - зубцовое деление статора, у - смещение ротора относительно статора, Ц& -магнитное напряжение, приложенное к зубцовому слою, Ф& - магнитный поток на участке зубцового деления
Средний момент двигателя определяется по выражению
2п
M
дв
cp 2п
1
— IM дВ (i, Y)dY =
dY
(4)
1 2л/ т \
=Е мфk(^, у)
2п 0 Iк=1 )
С учётом того, что
2п ( т ^ т (2п
| ЕМФ£ (л у)^у = Е IМФ£ ^ у)й?у
0 V £=1 ) £=1^ 0
и что среднее значение момента фазы вычисляется по выражению
2п
(5)
1 2п
M ф cp = M ф(i, Y)dY,
^cp 2п
G
выражение (2.4) примет вид
дв cp 1Mфср k .
k=1
(6)
1 T 1 Z
Mфcp = ^“ J^ф(i, Y)dY =72“Л IMмdY =
cp 2п о I ^ I G
Z-J °
_z_
2п
Ay
2n
Z
(7)
J MM dy+ J Gdy
2n
-Ay
Z - І^м -ay 2п
С учётом выражения (3) можно записать
дW'
М м = 2Z-
ду
откуда
W' z= ^м. z 2Z
Y.
(8)
(9)
Тогда приращение коэнергии на участке зубцового деления на интервале Ду составит
AWz =
M м -AY
2Z '
(1G)
Подставляя произведение Мм -Ду из (10) в (7) находим, что средний момент фазы равен
M
z 2 -aw!
<tc
(21)
^p п
и, соответственно, средний момент двигателя вычисляется по выражению
m - Z2 - AW!
(12)
"“cP n
Очевидно, что наибольшее приращение коэнергии AWz будет иметь место в случае, когда подключение фазы к источнику питания будет происходить при рассогласованном положении зубцов, а отключение - при согласованном. Это соответствует углу поворота п/Z. В этом случае приращение коэнергии равно площади криволинейной фигуры, ограниченной кривыми 1 и 2 (заштрихованная область на рис. 2):
AW'z = S(a) - S(u), (13)
где
a
S(a) = J Ф^^)dU .
G
Ua
S(u) = J Ф(u)(U)dU ,
(14)
(15)
Приняв допущение о том, что нарастание и спад тока происходят мгновенно, интервал работы фазы составляет Ду радиан по геометрическому углу поворота ротора и момент не изменяется по величине на протяжении интервала Ду и равен Мм, вычислим средний момент фазы по выражению (5):
2п
2п , г
В этих выражениях Ф(а)(и) - вебер-амперная характеристика (ВАХ) зубцового слоя в согласованном положении, Ф(м)(и) - ВАХ зубцового слоя в рассогласованном положении.
Представим зависимость Ф(и) в виде
Ф = 8(и)-и , (16)
где 8 (и) - магнитная проводимость зубцового слоя, и- магнитное напряжение.
Тогда выражения (2.14) и (2.15) примут вид
Ua
S(a) = Jga (U )U ]U,
G
Ua
S(u) = Jgu (U )U ]U .
(17)
(18)
G
Здесь ga -проводимость магнитной цепи в согласованном положении, gu - проводимость магнитной цепи в согласованном положении.
Таким образом, задача определения среднего момента двигателя сводиться к определению закона изменения магнитной проводимости зубцового слоя от магнитного напряжения для согласованного и рассогласованного положений.
Магнитная проводимость может быть определена следующими способами:
- по результатам расчета магнитного поля зубцового слоя в целом численными методами;
- на основании схемы замещения магнитной цепи, магнитные проводимости элементов которой определены на основании упрощенных аналитических ме-
G
тодов расчета магнитного поля по геометрическим размерам и свойствам материалов магнитопровода;
- путём совмещения первых двух подходов, когда магнитные проводимости части элементов схемы замещения определены по результатам расчета магнитного поля численными методами, а части - упрощенными аналитическими.
Ф
Рис. 2 Энергетическая диаграмма зубцового слоя 1 - вебер-амперная характеристика системы при согласованном положении; 2 - вебер-амперная характеристика при рассогласованном положении; Фа -магнитный поток зубца при согласованном положении зубцов и насыщенной магнитной системе; Фи - магнитный поток зубца при рассогласованном положении зубцов; иа -магнитное напряжение зубцового слоя, соответствующее номинальной магнитодвижущей силе фазной катушки
Поскольку геометрические размеры зубцового слоя определяются многими параметрами, определение магнитной проводимости зубцового слоя или его элементов по результатам решения серии полевых задач представляется нецелесообразным исходя затрат требуемого для этого времени. В связи с этим, нами принят второй подход, где магнитная проводимость зубцов статора и ротора строятся по их геометрическим размерам и кривой намагничивания материала, а магнитная проводимость воздушного зазора определяется по результатам расчёта коэффициента проводимости воздушного зазора по методу Поля [5].
Для согласованного положения, при котором индукция в частях магнитопровода превышает индукцию насыщения стали магнитопровода, схема замещения магнитной цепи будет содержать три элемента (рис. 3 а). Для рассогласованного положения в виду с низкими индукциями можно положить проводимости зубцов равными бесконечности и в схеме замещения оставить только магнитную проводимость воздушного зазора (рис. 3б).
§11 §ёа %г2 §
У у ьди
&—-Л^/ I--------1 I—I—0 &-1 I—0
а) б)
Рис. 3 Схема замещения зубцового слоя а - для согласованного положения б - для рассогласованного положения §21 - нелинейную проводимость зубца статора
- линейную проводимость воздушного зазора в согласованном положении §г2 - нелинейную проводимость зубца ротора §1и - линейная проводимость воздушного зазора в рассогласованном положении.
Тогда геометрические размеры зубцового слоя, через которые определяются вышеуказанные магнитные проводимости, и магнитная индукция в воздушном зазоре в согласованном положении будут служить параметрами, от которых зависит электромагнитный момент. Поиск оптимального сочетания этих размеров может быть сформулирован в виде задачи оптимизации АИД.
ВЫВОДЫ
Разработана математическая модель для определения электромагнитного момента индукторного реактивного двигателя с аксиальным магнитным потоком. Особенностью модели является нахождение момента путём интегрирования закона изменения магнитной проводимости.
СПИСОК ЛИТЕРАТУРЫ
1. Ильинский Н.Ф. Проектирование вентильноиндукторных машин общепромышленного назначения [Текст] / Ильинский Н.Ф., Штайнбрунн Й., Прудникова Ю.И. // Вестник МЭИ, - №1, 2004, с. 37 - 43.
2. Иванов-Смоленский А.В., Абрамкин Ю.В., Власов А.И. и др Универсальный метод расчета электромагнитных процессов в электрических машинах / Под ред. А.В. Иванова-Смоленского. - М.: Энергоатомиздат, 198б. - 213 с.
3. Птах, Г. К. Развитие методов расчёта электромагнитных процессов в электромеханических системах с индукторными машинами [Текст] : Монография / Г.К. Птах. - Юж.-Рос. гос. техн. ун-т. Новочеркасск : Ред. журн. «Изв. вузов. Электромеханика», 2003. - с. 232
4. Ткачук, В.І Електромеханотроніка [Текст] : Підручник / В Ткачук - Львів: Видавництво Національного університету «Львівська політехніка», 200б. - с. 440.
Поступила 30.09.2010
Рябов Евгений Сергеевич,
Любарский Борис Григорьевич, к.т.н., доцент,
Зюзин Дмитрий Юрьевич, к.т.н.
Национальный технический университет “Харьковский политехнический институт”
Украина, б1002, Харьков, ул. Фрунзе, 21, НТУ “ХПИ”, кафедра электрического транспорта и тепловозостроения тел. (057) 707-б5-З0, E-mail: lboris@kpi.kharkov.ua Емельянов Вадим Леонидович Национальный технический университет “Харьковский политехнический институт”
Украина, б1002, Харьков, ул. Фрунзе, 21, НТУ “ХПИ”,
кафедра электрических аппаратов
тел. (057) 707-б8-б4, E-mail: evl@kpi.kharkov.ua
E.S. Rjabov, B.G. Lubarskiy, D.Ju.Zjuzin, V.L. Emeljanov To question of determination of electromagnetic moment of inductor ramjet with axial magnetic flux.
A mathematical model is resulted for the decision of electromagnetic moment of inductor motor with the axial magnetic flow. Of model finding is a feature of moment by integration of law of change of magnetic conductivity.
Key words - mathematical model, electromagnetic moment, inductor motor, axial magnetic flow.
