УДК 62-503.56 Воронин С.Г., Хафизов Г.Т.
НЕКОТОРЫЕ СПОСОБЫ РЕГУЛИРОВКИ ВЫХОДНЫХ ПАРАМЕТРОВ ЭЛЕКТРОПРИВОДА С ВЕНТИЛЬНЫМ ДВИГАТЕЛЕМ ПРИ ВЕКТОРНОМ УПРАВЛЕНИИ
В данной работе предлагается синтез нового метода регулировки параметров привода, который позволит упростить задачу векторного управления вентильным двигателем. Структурой решения данной проблемы является рассмотрение исследуемого электропривода как сложного составного объекта, получение математической модели объекта и определение задач из структуры, необходимых для решения проблемы. В ходе работы были получены аналитические выражения динамической модели электропривода, получены выражения, описывающие работу привода в псевдостатическом режиме, также точные аналитические выражения для определенного оптимального режима: максимальная скорость при известном моменте сопротивления, предложены методы использования. Проведено компьютерное моделирование синхронного двигателя с постоянными магнитами в режиме вентильного двигателя, получены соответствующие графики процессов.
Ключевые слова: вентильный двигатель, синхронный двигатель с постоянными магнитами, динамическая модель, оптимальный угол коммутации, максимальная скорость под нагрузкой.
Введение
Оптимальные режимы работы вентильных двигателей на основе синхронных двигателей с постоянными магнитами являются приоритетной задачей теории электропривода. Цель выполненной работы - получение точного аналитического выражения для оптимального угла коммутации при режиме максимальной скорости.
Решение поставленной задачи заключается в определении необходимой математической модели синхронного двигателя с возбуждением от постоянных магнитов. Далее аналитическое определение исследуемых характеристик привода в зависимости от угла коммутации, проведение компьютерного моделирования исследуемого электропривода с целью подтверждения результатов исследования. Причём, в отличие от широко распространённых методов векторного управления путём автономного регулирования токов по осям d и q, требующих привлечения значительных вычислительных ресурсов [2-4], управление предлагается осуществлять путём регулирования угла коммутации 9 [5, 10], которое при наличии датчика углового положения ротора реализуется простейшими средствами.
Физико-математическая модель
ВЕНТИЛЬНОГО ДВИГАТЕЛЯ
Математическая модель двигателя была создана для синхронного двигателя с неявнополюсной двухфазной однополюсной статорной обмоткой и постоянными магнитами на роторе и идеализированным синусоидальным управлением. За основу была взята модель с описанием физических обмоток как распределенных токовых слоев, позволяющих получить вращающееся синусоидально распределенное в пространстве магнитное поле статора. На математическую модель накладываются следующие ограничения: скорость распространения фронта электромагнитных волн считается бесконечной, магнитная проницаемость магнито-провода тоже считается бесконечной, отношение внутреннего радиуса статора к радиусу ротора равно единице, неучтены магнитные и электрические поля высших порядков, а также явление гистерезиса и вихревых
© Воронин С.Г., Хафизов Г.Т.
токов, симметричности электрической машины [1].
После всех предварительных преобразований из физической модели, основанной на величинах электродинамики, была получена модель в системе отсчета вращающегося ротора, используя величины теории электрических цепей [8]:
3и 008(0} = Я'Г + Ь'" + ш(г) + 3кш(г); " ёг "
о
- 3и 8ш(0) = Я'Г, + Ь'" ^ - Ь"Гш(г);
2 а 4
Мс = Мэм - Г
М.
ёш(г) ёг '
ёг
(1)
к
эм _ у'
= ,
где и, I - амплитуды напряжений и токов системы питания статора в системе отсчета статора, для трехфазной машины; //, // - токи статорных обмоток эквивалентного двухфазного двигателя в системе отсчета ротора dq^, Я" - активные сопротивления фазных обмоток двигателя; Ь - полная индуктивность обмотки фазы; М - механический момент, создаваемый двигателем без нагрузки на валу; Мэм - электромагнитный момент двигателя; X - момент инерции ротора двигателя; к - коэффициент пропорциональности потокос-цепления и ЭДС.
Давайте преобразуем и упростим систему уравнений (1) к значениям в относительных единицах, для этого запишем систему обозначений безразмерных величин:
и
ЦЯ'
У =
1,5и„
1,5и
Ъ=тЬ ;х=гюн;
Ь
ТЭ =Т7шн; шн = Я
1,5и
е = -
ш
(2)
к М„
ш
М М 1,5и к
м =-эм • П =-с. • М = -
Гэм , , ; Г с , , ; М П
Мп Мп
П
Я'
,ГЯ' 1,5и к
где у- относительное напряжения статора эквивалентного двухфазного двигателя в системе отсчета ротора dq^, [ ' - приведенный ток обмотки q эквивалентного
Ч
н
^м =
двухфазного двигателя в системе отсчета ротора ёд; в - относительная скорость вращения ротора; г/ - приведенный ток обмотки ё эквивалентного двухфазного двигателя в системе отсчета ротора ёд; т - безразмерное время; тэ - электромагнитная постоянная времени электродвигателя; юн - эквивалентная частота; тм -механическая постоянная времени электродвигателя; цэм - относительный электромагнитный момент эквивалентного двухфазного двигателя в системе отсчета ротора ёд; |с - относительный момент сопротивления на валу эквивалентного двухфазного двигателя в системе отсчета ротора ёд; Мп - приведенный момент.
Преобразованная система (1) с помощью обозначений (2) приводит к выражению в домене Лапласа:
у соб
(е)=г;+тэ \\р+тэ ¡а в+в;
у ^(е)=^+тэ %р -тэ Кв;
: в _
Тм 1 = 1 эм — ис;
ат
и = I .
(3)
Определение оптимального угла нагрузки е
Используя систем (3) в условиях установившихся процессов, выразим значения приведенных токов
Определим угол е, при котором наблюдается максимальная скорость вращения ротора вентильного электропривода при неком известном моменте |и> Воспользовавшись четвертой и третьей строкой выражения (3) и выражением (4), найдем выражение скорости для установившихся процессов при условии цэм=цс.
_ _у со8(е)+тэву Бт(е) — в
иэм = ис =
(6)
12
+ ТЭ 8
Выразим в из выражения (6) и получим
-_(тэ у в1п (е)—1)
8 =
(тэ у вт
21сТэ
+
+
л/(тэ у зт(е)—1)2 — 4|стэ2 (|| — у сов(е))
(7)
21сТэ
Возьмем производную в (7) по углу е при условии экстремума и получим.
— = 0 = ту сов(е)±
:е э 4 7
^ 2тэу сов(еХтэу 5Ш(е) —1)— 41стэ2у 5Ш(е)
2л1 (тэу 51п(е)— 1)2 — 4|стэ2 (1с — у сов(е))
(8)
Раскроем и упростим выражение (8), приведя к выражению вида
г: =
тэ ву сов(е)—у вт(е)— тэ в
^ =
1 + тэ 2 ш2 (т) у сов(е)+тэ ву в1п(е)—в
1 + тэ в
(4)
вт
(е)+
1сТ
с ''З
соб
(е)
= тэу.
(9)
Данное выражение (9) решается аналитически с заменой Х=соБ(е) и приводится к выражению
(|сТэ2 )— (2|сГТэ2 X + (тэ2у2 — 1^2 + ^4 = 0. (10)
Данное уравнение (10) решается методом Ферра-ри-Кардано [7]. В целях экономии места, мы запишем сразу решение уравнения (10), для этого введем необходимые сокращения.
2
2
а
я
А=тэ2у2 — 1;
В ' = А ' 3 /27 + 2П л/3/9 + 4|2тэ 2 А3 + 2|с2тэ 4у; С' = у!9В'23 —6 А В'13+12|2^2 + А' 2;
р' = V 27|4тэ8у 4 — 16|6тэб — |2тэ2 А' 4 +8|4ТЭ 4 А' 2 —36|4тэ6 у2 А+|Уу2 А' 3; Е ' = д/ А'3 +6л/3£ '—36|2ъ2 А '+54|2ъ 4у2; р, = — А • 2 С '—9СВ '23—12 А 'С'В13—1212x3 С '; / = — 8т(ли/2 — л/ 4); g = сов(ли).
Получим решение уравнения (10) в виде
0„_ = агссоБ
■ +
+
6В1/6
4Р + ^ (/^^мцг эУЕ
1/6
(12)
+ лс, с е 2,
где и - порядковый номер решения, и принимает значения 1, 2, 3, 4.
Также необходимо заметить, что для предварительной оценки можно найти оптимальный угол коммутации по формуле [9]
0
ТЭ (У-Пс).
(13)
ИССЛЕДОВАНИЕ АНАЛИТИЧЕСКОГО ВЫРАЖЕНИЯ ОПТИМАЛЬНОГО УГЛА КОММУТАЦИИ
Исследуем выражение (10) на количество положительных корней, для этого воспользуемся теоремой Декарта [6] для знаков полинома.
Согласно теореме: количество положительных корней многочлена с вещественными коэффициентами равно количеству перемен знаков в ряду его коэффициентов или на чётное число меньше этого количества, для рассмотрения количества отрицательных корней надо брать функцию с отрицательным аргументом.
Применим теорему Декарта к выражению (10):
X4 + (тэ 2у2 -1 )Л2 -
К
(2мс утэ2 МмДэ2 )=©(х)
шт(©(х) = 0 Х> 0) =
= V! 0, (тэ2у2 - ^ -(2псутэ 2 ) (мтэ 2 ) )-2к; (14)
пгг(©(-Х) = 0 X < 0=
= v(l, 0, (тэ2у2 - ^ (2Псутэ2 ) (Пстэ2 )) - 2к;
к е Я,
где V - функция количества знакоперемен.
Исходя из выражения (14) на суммарное количество рациональных корней крайне влияют знаки управляющего напряжения и момента сопротивления. Может так сложится, что оптимального угла при заданных параметрах не существует.
Также исследуем выражение (9) на запрещенные углы коммутации. При нулевом знаменателе выражение (9) теряет смысл, поэтому запишем множество запрещенных углов:
л
00=- + л^ Се2.
(15)
Исходя из выражений (11), (12), (14), (15), можно записать необходимое условие локального максимума, для условия максимальной скорости вентильного двигателя на основе синхронного с постоянными магнитами, как сложную функцию угла коммутации:
0_„ = агссоБ
'г.
Ы/ С
+
+
6 В16
Р + 81вп (/ )бл/Г2мс Тэ2 уЕ
1/6
+ лс, се 2;
(0таХ -А0) + -^ С0Б
(0)
А0е Я;
- ^э У > 0,
(16)
0 < А0 = |0П
(0тах +А0) + -Ъ^ -Тэ у< 0;
С0Б
(0)
л
0тах ^00 =~ + ЛС, Се2-
Промоделируем систему уравнений (3) для условий у=1; тэ=0,2, также для оптимального углов и без регулирования по углу 0=0 и построим механические характеристики.
1 2 3
Механические характеристики вентильного двигателя с коррекцией и без, по оптимальному углу коммутации, для режима максимальной скорости под нагрузкой (штрихованная линия - это скорректированная механическая характеристика, гладкая - скорректированная)
Как видно из рисунка, скорректированная механическая характеристика не имеет момента короткого замыкания, в этом случае происходит срыв работы привода. Но за счет коррекции значительно расширяется рабочий диапазон скоростей привода.
В силу выражений (18) и (16) и с учетом закона Мура, при нынешнем развитии науки и техники, предлагается два варианта применения оптимального угла
0тах в практике.
Введем тензор типа
®=т; (мс, Тэ, у)).
(19)
В случае малой вычислительной мощности контроллера электропривода строится трехмерный тензор-массив (19) в долгосрочной энергонезависимой памяти микроконтроллера, содержащий дискретные значения оптимального угла 0тах с шагом дискретизации т, ко-
3
личеством элементов п , для широкого ряда параметров Пс, тэ, у, предварительно рассчитанных по выражениям (16), (18) и упорядоченных согласно возрастанию
<
по индексам i, j, k.
В случае удовлетворительной вычислительной мощности микроконтроллера при псевдостатических режимах электропривода необходимо вести расчет согласно алгоритму по выражениям (16) и (18).
Заключение
Поставленные в исследовании задачи решены и получены точные аналитические выражения для регулирования угла коммутации в процессе векторного управления электроприводом, в режиме максимальной скорости при известной нагрузке на валу. Получены необходимые динамические и статические математические модели, предложены методы эксплуатации полученных выражений в электроприводе. Все предлагаемые в исследовании соотношения подкреплены конкретными расчётами и компьютерным моделированием.
СПИСОК ЛИТЕРАТУРЫ
1. Уайт Д.С., Вудсон Г.Х. Электромеханическое преобразование энергии. М.;Л.: Энергия, 1964. 528 с.
2. Model reference adaptive backstopping based direct torque control of PM synchronous motor drives / Hai, L., Yan, W.-S., Li, H., Wei, Y.// 4th IEEE Conference on Industrial Electronics and Applications, ICIEA, 2009 , art. no. 5138380, pp. 1173-1178.
Information in English
3. Speed-sensorless direct torque control system of permanent magnet synchronous motor / X. Wang, Y. Xing, Y. Liu, D Yang, // Dongbei Daxue Xuebao/Journal of Northeastern Uni-versity,2012. №33(5). Pp. 618-621.
4. Sensor less control of synchronous machines based on direct speed and position estimation in polar stator-current coordinates / D. Paulus, J.-F.Stumper, R. Kennel, IEEE Transactions on Power Electronics № 28 (5) , art. no. 6257501, 2013, pp. 2503-2513.
5. Воронин С.Г., Курносов Д.А., Кульмухаметова А.С. Некоторые схемы и алгоритмы векторного управления синхронными двигателями с возбуждением от постоянных магнитов // Электротехника. 2013. №10. С. 50-54.
6. Анго А. Математика для электро и радио инженеров с предисловием Луи Де Бройля / Андре Анго; пер с фр. под общ. ред. К.С. Шифрина. М.: Наука, 1965. 780 с.
7. Справочник по математике для инженеров и учащихся втузов / И.Н. Бронштейн, К.А. Семедяев; под ред Г. Гроше, В. Циглера. М.: Наука, 1986. 544 с.
8. Бессонов Л.А. Теоретические основы электротехники: учеб. пособие. М.: Высш. шк., 1967. 780 с.
9. Юферов Ф.М. Электрические машины автоматических устройств: учебник для студентов вузов, обучающихся по спец. «Электромеханика». 2-е изд., перераб. и доп. М.: Высш. шк., 1988. 479 с.
10. Сизякин А., Румянцев М. Без датчика положения ротора: решения компании IR для управления вентильными двигателями // Новости электроники (Электронный журнал). - 2011. - № 10. - Ст. 6.
Some Ways to Control Output Parameters of Brushless Electric Motor with Vector Handling
Voronin S.G., Hafizov G.T.
In this paper, we propose a new synthesis method of the drive parameter adjustment, which will simplify the task of vector control of a brushless motor. The structure of the solution to this problem is to consider the test actuator as a complex object, to obtain a mathematical model of the object and to determine the tasks of the structures necessary to solve the problem. in the course of the research work, analytical expressions of the dynamic electric model were obtained, which describe the operation of the drive in a pseudo-static mode, as well as accurate analytical expressions to determine the optimal mode: max speed at a certain moment of resistance. Computer simulation of a synchronous motor with permanent magnets in the ac electronic motor was carried out, besides, the authors obtained the appropriate process diagrams.
Keywords: brushless electric motor, a synchronous motor with permanent magnets, dynamic model, optimum angle switching, maximum speed under load.
References
1. White D.S, Woodson G.H. Elektromekhanicheskoe preobrazovanie energii [Electromechanical energy conversion] Moscow, Leningrad, Energy, 1964, 528 p
2. Hai, L., Yan, W.-S., Li, H., Wei, Y.Model reference adaptive backstopping based direct torque control of PM synchronous motor drives, 4th IEEE Conference on Industrial Electronics and Applications, ICIEA, 2009, art. no.5138380, pp.11731178.
3. Wang X., Xing Y., Liu Y., Yang D. Speed-sensorless direct torque control system of permanent magnet synchronous motor, Dongbei Daxue Xuebao, Journal of Northeastern University, 2012, no.33(5), pp.618-621.
4. Paulus D., Stumper J.-F., Kennel R. Sensor less control
of synchronous machines based on direct speed and position estimation in polar stator-current coordinates, IEEE Transactions on Power Electronics, no.28(5), art. no.6257501, 2013, pp.25032513.
5. Voronin S.G., Kurnosov D.A., Kulmuhametova A.S. Nekotorye skhemy i algoritmy vektornogo upravleniya sinkhronnymi dvigatelyami s vozbuzhdeniem ot postoyannykh magnitov [Some Schemes and Vector Control Algorithms for Synchronous Motors with Permanent-magnet], Electronics, 2013, no.10, pp.50-54.
6. Ango A. Matematika dlya elektro i radio inzhenerov s predisloviem Lui De Broilya [Mathematics for Electrical and Radio Engineers with a foreword Louis de Broglie], André Angot; Translated from French, under the editorship of K.S. Shifrin, Moscow, Science, 1965, 780 p.
7. Bronstein I.N., Semedyaev K.A., Spravochnik po matematike dlya inzhenerov i uchaschikhsya vtuzov [Handbook of Mathematics for Engineers and Students of Technical Colleges], Ed G. Grosche, W. Ziegler, Moscow, Science, 1986, 544 p.
8. Bessonov L.A. Teoreticheskie osnovy elektrotekhniki [Theoretical Foundations of Electrical Engineering], textbook, Moscow, Higher School, 1967, 780 p.
9. Yuferov F.M. Elektricheskie mashiny avtomaticheskikh ustroistv [Electrical Machinery of Automatic Devices]. textbook for students enrolled in the specialty "Electromechanics",2nd ed., revised and extended, Moscow, Higher School, 1988, 479 p.
10. Sizyakina A., Rumyantsev M. Bez datchika polozheniya rotora: resheniya kompanii IR dlya upravleniya ventilnymi dvigatelyami [Without Rotor Position Sensor: Solutions of IR in the Field of Brushless Motors Control], Electronics News (Electronic Journal), 2011, no.10, Article 6.