CC BY
14
ЛВЕСТНИК
............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
DOI: 10.21821/2309-5180-2020-12-3-606-618
CONTROL OF A REACTIVE ELECTRIC MOTOR WHEN FEEDING THE WINDINGS WITH A SINUSOIDAL CURRENT
V. F. Samoseiko, E. V. Shiryaev
Admiral Makarov State University of Maritime and Inland Shipping, St. Petersburg, Russian Federation
Vector control of the electromagnetic moment of a reactive machine is considered. The control task is to eliminate pulsations of the electromagnetic moment. It can be shown that to achieve this goal, sinusoidal currents must flow through the windings, whose inductance must pulse according to the sinusoidal law. It is shown that when sinusoidal currentsflow through the phase windings, the voltage of the main and third harmonics occurs on them. At the same time, in a three-phase machine whose phase windings are connected to a star, the third harmonic disappears as part of the line voltage. In this case, a standard frequency converter can be used to synthesize the voltages applied to the stator windings using pulse-width modulation, which makes it easier and cheaper to practically control the reactive electric machine with concentrated phase windings.
The ratio of size of the stator and rotor teeth, at which the coefficient of non-sinusoidal pulsations of the phase windings inductances during the rotor rotation will be minimal, is found. Electromagnetic processes in the reactive machine are described by the linear differential equations with periodic coefficients. Lyapunov transformations for differential stress equations, which allow converting the initial equations with periodic coefficients into stress equations with constant coefficients, are found. The resulting Lyapunov equations are the basis for the synthesis of a system for controlling the dynamics of electromagnetic processes, which allows sinusoidal currents to flow through the phase windings of the machine. The control system of the reactive electric machine is represented by a structural diagram. An algorithm for controlling the electromagnetic moment of the reactive machine with a constant magnetization current is considered. The use of this algorithm allows to obtain the maximum performance of the management. A virtual scattering loop is also used to increase the speed of load current regulation. Simulation of electromechanical processes at the reactive electric machine control is carried out.
Keywords: vector control, switched reluctance motor, alternating sinusoidal current, Lyapunov transformation.
For citation:
Samoseiko, Veniamin F., and Eduard V. Shiryaev. "Control of a reactive electric motor when feeding the windings with a sinusoidal current." Vestnik Gosudarstvennogo universiteta morskogo i rechnogoflota imeni admirala S.O. Makarova 12.3 (2020): 606-618. DOI: 10.21821/2309-5180-2020-12-3-606-618.
УДК 621.3.072.6
УПРАВЛЕНИЕ РЕАКТИВНЫМ ЭЛЕКТРИЧЕСКИМ ДВИГАТЕЛЕМ ПРИ ПИТАНИИ ОБМОТОК СИНУСОИДАЛЬНЫМ ТОКОМ
со В. Ф. Самосейко, Э. В. Ширяев г
^ ФГБОУ ВО «ГУМРФ имени адмирала С. О. Макарова»,
Л Санкт-Петербург, Российская Федерация
Рассматривается векторное управление электромагнитным моментом реактивной машины, целью которого является устранение пульсаций электромагнитного момента. Для этого по фазным сосредоточенным обмоткам машины должны протекать синусоидальные токи, а их индуктивности пульсировать по синусоидальному закону. В рабате показано, что при протекании синусоидальных токов по фазным обмоткам на них возникает напряжение основной и третьей гармоники. При этом в трехфазной машине, фазные обмотки которой соединены в звезду, в составе линейных напряжений третья гармоника исчезает. В этом случае для синтеза напряжений, подводимых к обмоткам статора методом широтно-импульсной модуляции, может быть использован стандартный преобразователь частоты, что позволяет упростить и удешевить практическую реализацию управления реактивной электрической машиной с сосредоточенными фазными обмотками. Определены соотношения размеров зубцов статора и ротора, при выполнении
см о
ВЕСТНИК«
ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ЕШщ^
МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА,
которых коэффициент несинусоидальных пульсаций индуктивностей фазных обмоток при вращении ротора будет минимальным. Электромагнитные процессы в реактивной машине описаны линейными дифференциальными уравнениями с периодическими коэффициентами. Найдены преобразования Ляпунова для дифференциальных уравнений напряжений, позволяющие преобразовать исходные уравнения с периодическими коэффициентами в уравнения напряжений с постоянными коэффициентами. Полученные уравнения Ляпунова положены в основу синтеза системы управления динамикой электромагнитных процессов, позволяющей обеспечить протекание синусоидальных токов по фазным обмоткам машины. Система управления реактивной электрической машиной представлена структурной схемой. Рассмотрен алгоритм управления электромагнитным моментом реактивной машины с постоянным током намагничивания. применение которого позволяет получить максимальное быстродействие управления. Для увеличения быстродействия и робастности управления также используется контур виртуальной диссипации. Выполнено моделирование электромеханических процессов при управлении реактивной электрической машиной.
Ключевые слова: векторное управление, реактивная электрическая машина, синусоидальный ток, преобразование Ляпунова.
Для цитирования:
Самосейко В. Ф. Управление реактивным электрическим двигателем при питании обмоток синусоидальным током / В. Ф. Самосейко, Э. В. Ширяев // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2020. — Т. 12. — № 3. — С. 606-618. DOI: 10.21821/2309-5180-2020-12-3-606-618.
Введение (Introduction)
Реактивные электрические машины (РЭМ), магнитопровод которых имеет зубчатую конструкцию статора и ротора, обладают свойствами, принципиально отличающими их от других электрических машин. Это отличие обусловлено отсутствием взаимных индуктивностей между фазными
Z n
обмотками. Число зубцов статора Z и ротора Z2 всегда четное. Отношение Z - _ образует несократимую дробь, где m — целое число, называемое числом фаз статора; n — число пар зубцов ротора на полюсное деление статора, называемое также коэффициентом электрической редукции. Число зубцов статора и ротора определяются также соотношениями Z1 = 2mp; Z2 = 2np. Наибольший общий делитель числа зубцов ротора и статора является числом полюсов 2р, где р — число пар полюсов.
Пример поперечного разреза магнитной системы такой машины приведен на рис. 1. Необходимым условием функционирования РЭМ являются пульсации индуктивностей обмоток статора. На рисунке показано положение ротора, в котором индуктивность обмотки A, образованной двумя катушками, максимальна, обозначается Ld и называется продольной. При повороте ротора на 45 геометрических градусов индуктивность обмотки A будет минимальной. В этом случае она обозначается буквой Lq и называется поперечной. Очевидно, что индуктивности фазных обмоток A, B, C при вращении ротора пульсируют.
Рис. 1. Пример поперечного разреза магнитной системы реактивной машины (2 = 6, Z2= 4, т= 3, п = 2, р = 1)
ЛВЕСТНИК
............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
Сосредоточенные фазные обмотки реактивных электрических машин имеют короткие лобовые части. Это обеспечивает их преимущество перед другими электрическими машинами, а именно: низкие электрические потери мощности, высокую надежность и технологичность изготовления. Отсутствие магнитных связей между фазными обмотками позволяет алгоритмически исключить отказавшую фазу из работы и обеспечить работоспособность машины. Отсутствие обмоток на роторе и простота его конструкции обеспечивают низкий момент инерции. Кроме того, упрощается система охлаждения машины. Важнейшей особенностью реактивных машин, имеющих зубчатую конструкцию ротора и статора, является принципиальная возможность обеспечения минимальных пульсаций электромагнитного момента, что способствует существенному улучшению виброшумовых характеристик таких машин. Указанные особенности реактивных электрических машин создают перспективу их применение на объектах водного транспорта.
Проектированию электрических машин с зубчатым статором и ротором посвящено большое количество как отечественных [1]-[4], так и зарубежных [5], [6] работ. В настоящее время общепринятым является питание фазных обмоток импульсным напряжением. Это приводит к протеканию по обмоткам полигармонических токов и сопровождается существенными пульсациями электромагнитного момента [7]-[9]. В публикации [10] показано, что при синусоидальных пульсациях индуктивностей фазных обмоток электромагнитный момент двигателя создает первая гармоника тока, а также сумма нулевой и второй гармоник тока, остальные гармоники приводят лишь к пульсациям электромагнитного момента. В работах [11] и [12] рассматривается способ управления индукторной машиной, имеющей на статоре специальную обмотку возбуждения, синусоидальными напряжениями, однако последние порождают несинусоидальные токи, которые ведут к пульсациям электромагнитного момента.
В данной работе рассматривается управление синусоидальными токами реактивной машиной без специальной обмотки возбуждения на статоре. Для формирования напряжений, подводимых к обмоткам статора, может быть использован стандартный преобразователь частоты [13]. Для снижения пульсаций электромагнитного момента при синусоидальных токах в обмотках статора необходимо, чтобы индуктивности обмоток машины пульсировали по синусоидальному закону. В данной работе приводятся конструктивные соотношения зубцов статора и ротора машины, минимизирующие коэффициент несинусоидальных пульсаций индуктивностей.
Относительные единицы. Параметры РЭМ и переменные удобно представлять в относительных единицах — в виде отношения параметра к его базовому значению. Параметры, представленные в относительных единицах, помечаются верхним индексом *. Для представления параметров в относительных единицах и перехода обратно к параметрам в именованных единицах вводятся базовые величины.
Базовые величины РЭМ делятся на основные и производные от них. Основными базовыми величинами являются амплитуда первой гармоники номинального фазного напряжения — г иб = и(1), амплитуда первой гармоники номинального тока: 1б=1н 42, а также номинальная угловая
^ частота первой гармоники тока статора — юб = юн.
о Производные базовые величины находятся из основных базовых величин: сопротивле-
« ние — Я = и II' индуктивность — Lб = Яб/юб; мощность — Р = т-и-172; электромагнитный
момент — М =Р-шб; угловая скорость вращения ротора — Об = юб/р , где т — число фаз; р —
Данные значения параметров использованы при расчете параметров регуляторов и построении иллюстрационных графиков.
Электрические углы. В классических машинах переменного тока число пар зубцов ротора и пар полюсов статора совпадают. Согласно ГОСТ 27471-87, «электрическим углом называется
Методы и материалы (Methods and Materials)
расчетное число пар полюсов.
Примем за параметры среднестатистического РЭМ следующие значения:
Я* = 0,03; Ь* = 2,00; L * = 0,30.
(1)
ВЕСТНИК,
ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ВШВ
МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА,
«произведение значения геометрического угла, образованного двумя полуплоскостями, проходящими через ось вращения машины переменного тока, на число пар полюсов»1. Таким образом, в соответствии с определением, электрический угол у связан с геометрическим углом Г соотношением
у = (Х^Т^ рп-Г, (2)
где рп = 72/2 — расчетное число пар полюсов; — число зубцов ротора.
Таким образом, электрическому периоду соответствует поворот ротора на два зубцовых деления.
Угловая скорость вращения ротора О = рГ, где р — здесь и далее оператор дифференцирования по времени. Из этого выражения и формулы (2) следует связь электрической угловой скорости вращения ротора с геометрической скоростью вращения ротора:
ю = рп-а. (3)
Обеспечение синусоидальных пульсаций индуктивностей фазных обмоток. Одним из условий устранения пульсаций электромагнитного момента является синусоидальность пульсаций индуктивностей при вращении ротора. Рассмотрим геометрические соотношения зубцов машины, которые позволяют свести искажения синусоидальности пульсаций индуктивностей теоретически к нулю.
Индуктивность фазной обмотки зависит угла поворота ротора и является периодической пульсирующей функцией, имеющей Z2 максимумов на геометрическом угле поворота 2р. Размах и форма пульсаций индуктивностей зависят от конструкции зубцов статора и ротора. На графике зависимости индуктивности Lk(у¡) от электрического угла поворота ротора у фазной обмотки к = 1, 2, ..., т, приведенном на рис. 2, размерам Р* и Р* соответствуют геометрические размеры зубцов Р* = (7/2)Р и Р* = (7/2)Р .
Рис. 2. Зависимость индуктивности фазной обмотки статора от электрического угла поворота ротора
2 О 2
При проектировании РЭМ с зубчатым статором и ротором следует стремиться к синусоидальной зависимости индуктивности фазной обмотки от угла поворота ротора. Коэффициент несинусоидальных искажений зависимости индуктивности от угла поворота ротора будет минимальным, если угловые размеры зубца статора и ротора будут связаны соотношением [14]
Р5 " 3 'Р* " ъ19
(5)
где Р^ и Р^ — угловой размер зубца статора и ротора (см. рис. 2). При этом скважности зубцов статора и ротора должны удовлетворять выражениям:
1 ГОСТ 27471-87. Машины электрические вращающиеся. Введ. 1988-07-01. Взамен ГОСТ 17154-71 и ГОСТ 23375-78. С. 36.
ЛВЕСТНИК
............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
PS z Pr 1
ts = z = ^ = 2, (6)
TS 3Z 2 xR 2
где т^ и tR — геометрический угловой размер зубцового деления магнитопровода статора и ротора соответственно.
Зависимость индуктивности фазной обмотки статора от электрического угла поворота ротора, приведенная на рис. 2, построена без учета выпучивания магнитного поля из зубца статора в воздушный зазор. Выпучивание магнитного потока в паз в положении «зубец-зубец» составляет около 8 ... 12 %. С учетом выпучивания магнитного потока зависимость индуктивности фазной обмотки статора от угла поворота ротора сглаживается и становится близкой к синусоидальной. Первая гармоника индуктивности фазной обмотки может быть записана в следующем виде:
Lk (У) = Lo + Lm • cos(2Tk)' (7)
где L0 — среднее значение индуктивности; Lm — амплитуда пульсаций индуктивности; gk = g + kp — электрический угол поворота ротора относительно фазы k = 1, ..., m; p = 2-p/m — угол фазового сдвига.
Среднее значение и амплитуда пульсаций индуктивности связаны с продольной и поперечной индуктивностями фазной обмотки соотношениями:
Lo = (Ld + Lq)/2; Lm = (Ld - Lq)/2. (8)
Кроме синусоидальной составляющей индуктивности обмоток имеют также высшие гармонические составляющие, которые обусловлены влиянием магнитопровода и геометрией воздушного зазора между статором и ротором. Однако при синтезе алгоритмов управления используется описание пульсаций индуктивности в виде уравнения первой гармоники (7), так как остальные гармоники обуславливают только пульсации электромагнитного момента.
Уравнения напряжений. При составлении уравнений напряжения полагается, что реактивная машина имеет m фазных обмоток, соединенных в звезду, которые получают питание от симметричной m-фазной системы ЭДС, обеспечивающих протекание по фазным обмоткам системы симметричных синусоидальных токов:
4 (У) = ia • sin(yk + 01) = id • sin(yk) - iq • c°s(ykX (9)
где id = ia • cos (07) — продольный ток (ток намагничивания); iq = ia • sin(07) — поперечный ток (ток нагрузки); ia = (id2+iq2)1/2 — модуль вектора тока в обмотке; Э7 = arctan(iq/id) — угол токовой нагрузки; k = 1, 2, ...; m — номера обмоток.
Взаимные индуктивности между фазами отсутствуют. Каждая из фазных обмоток (k = 1, ..., m) имеет активное сопротивление R1 и собственную индуктивность Lk(g), определенную выражением (7). Уравнения напряжений на фазных обмотках k = 1, ..., m имеют следующий вид:
ик = R • ik(У) + p[Lk(у)• ik(у)}, (10)
где ik(g) — ток фазной обмотки; g = wt — электрический угол поворота ротора; w — электрическая угловая частота; p — оператор дифференцирования по времени t.
Будем рассматривать продольный и поперечный токи в выражении (9) как функции электри-щ ческого угла поворота ротора g: id = id(g); iq = iq(g). Если подставить ток в выражении (9) в уравнения "Г напряжений (10), то их можно записать в относительных единицах как сумму первой и третьей
Jj гармоники:
S щ = + ик 3). (11)
Первая и третья гармоники фазного напряжения в осях d и q имеют вид:
и*«=u* • cos(y *)+и* • sin(y *(12)
u*(3) = u*(3) • cos(3y*) + u*(3) • sin(3y*). Уравнения продольного и поперечного контуров РЭМ. Используя первую формулу уравнения (12), можно записать уравнения амплитуд первой гармоники, связывающие напряжения и токи по осям d-q:
со г
и = с • ^ • £ + рч;
и* = R^ + «*' • С + • , (13)
X* Т *
^ и ь - — относительные индуктивности продольного и поперечного контуров.
Индуктивности в уравнениях (13) связаны с продольной и поперечной индуктивностями фазной обмотки соотношениями:
г, = ^; 4 = . (14)
Используя вторую формулу уравнения (12), можно записать уравнения амплитуд третьей гармоники, связывающие напряжения и токи по осям ^-д:
=—4-(рС-• С);
и*(з) = —4 • IРЦ + 3ю* • I*), (15)
где Ьт — амплитуда пульсаций полных индуктивностей обмоток статора (5).
Вычисление амплитуд координат. Для вычисления амплитуд токов и напряжений по осям необходим датчик, который позволяет определять электрический угол положения ротора у. Преобразование наблюдаемых токов в координатах магнитных осей обмоток в токи намагничива-
• * • * х.
ния г и нагрузки г выполняется по формулам:
2 т-1 <-) т-1
т ** • к); г*=т
£ = т '2£ • сов(ук); г* = — • Xг* • sin(Ук). (16)
т к=0 т к=0
Преобразование вектора напряжений в осях координат в основную гармонику фазных напряжений выполняется по формуле
и* = и* • со<у^) + и* • мп^, (17)
где к = 1, 2, ..., т — номера фазных обмоток.
Номинальный режим работы. В номинальном режиме работы относительная скорость вращения поля статора ю* = 1, а также модули векторов относительного напряжения и тока обмотки статора:
* I *2 *2 * I *2 *2
¡а = >/ Ч + \ = 1; иа = V ud + uq = 1. (18)
Если принять, что pid = р1 = 0, то из соотношений (18) и уравнений напряжений (13) можно найти относительные номинальные значения токов намагничивания и нагрузки:
Л - I1 - ^ • ; * - К -1 (19)
1аном ~ л *2 *2 5 ^ном ~ > *2 *2 ' \1У/ 0
о
где Ь* — индуктивности продольного и поперечного контуров (14); знак приближения в дан- 2
ных равенствах обусловлен допущением 0. М
Управление электромагнитным моментом. Электромагнитная энергия, запасенная обмот- Й
ками реактивной машины, фазные обмотки которой не имеют взаимных индуктивностей, опреде- к ляется выражением
1 т
Г(у) = 2-XLt(у)• И . (20)
2 к=1
Заметим, что электромагнитная энергия зависит от электрического угла поворота ротора, который связан с геометрическим углом соотношением (2). Тогда электромагнитный момент потенциальной силы РЭМ определится выражением
_о:
Гб11
ЛВЕСТНИК
............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
1 m
м_2
8Lt (у) ;2_ Z 2 m 8Lk (у) ,-2
СО
S
см
см о
612J
_ Z2 .у QL (у) _ 2 (21)
2 Ы 8Г - _ 4 Ы 8у (21)
Если в выражение (21) подставить токи (9), то электромагнитный момент РЭМ
M = Ри- m • (LD - LQ) • id • iq = Pn-m • LD • (1 • id • iq , (22)
где рп — расчетное число пар полюсов; LD и LQ — индуктивности продольного и поперечного контуров (14); X = LqLd — коэффициент поперечного рассеяния контуров; i* и i* — соответственно продольный и поперечный относительные токи обмотки статора.
Электромагнитный момент в относительных единицах:
M* = (Ld - LQ) • i* ■ i** = Ld • (1 - • id • i* . (23)
Если принять, что ток намагничивания поддерживается системой управления на номинальном уровне id = i , то из выражения (23) следует, что управление электромагнитным моментом сводится к управлению током нагрузки i.
Уравнение движения. Скорость протекания механических процессов характеризуется механической постоянной времени:
T = J б (24)
"Mex Мб ' (24)
где J — момент инерции электропривода.
Динамика механических процессов описывается уравнением движения ротора (вторым законом Ньютона):
J • pQ = M - Mc, (25)
где W — угловая скорость вращения ротора; М — электромагнитный момент; Мс — момент сопротивления движению.
Запишем уравнение движения машины, являющееся составной частью математической модели машины, в относительных единицах:
TMex ■ рю* = M* - M* , (26)
где ю* = W/W6 — относительная угловая скорость вращения ротора.
Результаты (Results)
Синтез динамики электромагнитных и механических процессов в РЭМ поясняется структурной схемой, приведенной на рис. 3.
При синтезе желаемого динамического поведения электромагнитных процессов принимаем, что наблюдается вектор токов в осях координат магнитных осей обмоток РЭМ. Для синтеза динамических процессов наблюдаемые токи должны быть преобразованы в токи намагничивания id* и нагрузки i* в системе координат ротора d-q по формулам (16).
Управляющими воздействиями является вектор фазных напряжений на обмотках РЭМ US. Алгоритм управления синтезирует напряжения ud* и u * в системе координат ротора d-q в относительных единицах, которые должны быть преобразованы в вектор фазных напряжений на фазных обмотках РЭМ US* по формулам (17). По напряжениям US* система управления электронно-ключевого преобразователя частоты (ПЧ) синтезирует линейные модулированные напряжения US.
В трехфазной РЭМ при синусоидальных линейных напряжениях по фазным обмоткам, соединенным в звезду, потекут токи, близкие к синусоидальным. При этом фазные напряжения будут иметь третью гармонику. При числе фаз m > 3 для получения синусоидальных токов в обмотках линейные и фазные напряжения РЭМ должны иметь третью гармонику.
ВЕСТНИК«!
ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^
МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
Рис. 3. Структурная схема управления скоростью вращения ротора в первой зоне при постоянном намагничивании: 1 — преобразователь частоты; 2 и 3 — блоки преобразования координат напряжения и тока; А — регулятор тока намагничивания; Б — регулятор тока нагрузки; Г — узел компенсации влияния тока нагрузки на ток намагничивания; В — узел ограничения тока нагрузки; Ж — регулятор скорости вращения ротора; И — узел ограничения модуля напряжения на обмотке статора
Регулятор тока намагничивания. Для управления током намагничивания образован контур с пропорциональным параллельным регулятором, который далее называется контуром виртуальной диссипации (рис. 3). Коэффициент передачи пропорционального параллельного регулятора Я* далее называется параметром виртуальной диссипации [15]. Несложно показать, что при значениях параметра виртуальной диссипации:
2 • Г • Г
* * 2 * *
Ях >ю ~ 2• Гд, (27)
характер динамических процессов будет апериодическим. Значение Я* существенно превышает относительное значение сопротивления обмотки статора Я*. В этом случае можно полагать, что объектом управления контура намагничивания является апериодическое звено с передаточной функцией
1 / /?*
ъ = Т^-Рл. (28)
Постоянная времени контура намагничивания с виртуальной диссипацией:
Г*
Тп =-Ц-*. (29)
т. • Я*
■ *-г
Влияние тока нагрузки ¡^ на контур намагничивания рассматривается как возмущающее воздействие, которое ведет к возникновению ошибок управления. Для астатического управления током намагничивания целесообразно образовать второй контур (см. рис. 3) с интегральным регулятором:
_ *
го о
го
_о:
Гб1э
Я*
^ = ■ (30)
Данный регулятор настроен на технический оптимум и обеспечивает желаемый характер динамических процессов. Передаточная функция контура тока намагничивания, настроенного на технический оптимум:
ЛВЕСТНИК
............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
СО
г
сч
=. (32)
2-Т02 ■р2 + 2-Т0 .р + 1 " 2-Г0-р + \ • (31)
Для компенсации влияния тока нагрузки iна контур тока намагничивания целесообразно предусмотреть отрицательную обратную связь то току нагрузки ¡^ (см. рис. 3).
Регулятор тока нагрузки. Синтез системы управление током нагрузки выполняется путем образования контура с пропорциональным параллельным регулятором, коэффициент передачи которого Я* (см. рис. 3) и второго контура тока нагрузки с интегральным регулятором:
Щ * Р '
Постоянная времени интегрального регулятора контура нагрузки с внутренним контуром виртуальной диссипации:
Тд , (33)
где LQ — относительная поперечная индуктивность; юб — базовая угловая частота; Я* — параметр виртуальной диссипации, удовлетворяющий соотношению (27).
Несложно показать, что при значении параметра виртуальной диссипации, удовлетворяющего неравенству (27), динамические процессы в контуре нагрузки становятся апериодическими. Передаточная функция контура тока нагрузки с интегральным регулятором, настроенным на технический оптимум:
Шк" = 2Гд2 • р2 ~+2Гд • р +1" 2 тв ~р +1 • (34)
Управление электромагнитным моментом. Динамика токов намагничивания и нагрузки определяет динамику электромагнитного момента. Динамику механических процессов в электрической машине определяет электромагнитный момент. Для обеспечения максимального электромагнитного момента машины необходимо иметь максимальное намагничивание магнитопровода. Однако чрезмерное намагничивание магнитопровода требует больших токов намагничивания, что ведет к большим потерям мощности и нагреву машины. Обычно при проектировании машины номинальный режим выбирается исходя из ее рациональной намагниченности. Поэтому целесообразно выбрать ток намагничивания статора на основе номинальных данных машины. Следовательно, для обеспечения максимального быстродействия машины необходимо поддерживать ток намагничивания постоянным, а управлять моментом контуром тока нагрузки. Поэтому в данной работе рассматривается управление электромагнитным моментом при постоянном номинальном токе намагничивания по формуле (19):
а=\\ !*- %. (35)
- Ьд
В этом случае связь между электромагнитным моментом и током нагрузки определяется выражением
М* = 4-(1-9-а, - I*, (36)
сэ где X = L* / L* — коэффициент поперечного рассеяния контуров.
^ Ограничение напряжения. Узел И на структурной схеме (рис. 3) предназначен для ограни-
чения напряжения на фазных обмотках, которое происходит при и* > и0 и при ¡* > ¡0. Мультипликативный сигнал ограничения модуля вектора относительного напряжения иа вычисляется по формуле
_ 1
% 1 + ехр((щ * - щ) -X).
Если относительное значение модуля напряжение и* находится в окрестности значения u0, а относительный ток i* — в окрестности значения i0, то происходит снижение тока нагрузки. Размер окрестности задается параметром l « 100... 300.
Ограничение тока нагрузки. Узел В на структурной схеме (рис. 3) предназначен для ограничения напряжения и тока нагрузки, которое происходит в случае, если vv • cq > aM при iq > aM Сигнал задания на ток нагрузки определяется выражением
a* = max(-aM, min (aM, Vu • c*)). (38)
Регулятор скорости вращения ротора. Для построения системы управления скоростью вращения ротора используется принцип подчиненного управления. При этом контур управления скоростью является внешним, а контур управления током нагрузки — внутренним. В данной работе рассматриваются алгоритмы подчиненного управления скоростью вращения ротора в первой зоне (-1 < ю* < 1).
Передаточная функция контура тока нагрузки (31), настроенного на технический оптимум, может быть аппроксимирована апериодическим звеном первого порядка с постоянной времени T . Несложно показать, что в этом случае регулятор контура скорости, настроенный на технический оптимум, является пропорциональным с коэффициентом усиления
T
Atq
(39)
Большим значениям кср соответствует более высокое быстродействие. Однако при этом снижается запас устойчивости. Сигнал задания на контур управления скоростью вращения ротора Оз для управления в первой зоне должен принадлежать интервалу [—1, 1].
Обсуждение (Discussion)
Моделирование выполнялось в соответствии со структурной схемой, приведенной на рис.3. При моделировании динамических процессов полагалось, что производится пуск среднестатистической РЭМ с параметрами, определенными выражениями (1), с постоянной величиной номинального момента сопротивления: M* = Мном*. Моделирование проводилось с параметром виртуальной диссипации Rx* = 1 и настройкой контуров управления токами на технический оптимум. При моделировании полагалось, что предварительное намагничивание магнитопровода отсутствует: /ДО) = 0.
Результаты моделирования динамических процессов представлены рис. 4. На рис. 4, а показана динамика токов статора iq , id, ia, скорости вращения ротора ш , электромагнитного момента M* и напряжения статора ua* в относительных единицах при W3* = 1 (подъем груза), на рис. 4, б — динамика тех же переменных при W * = -1 (спуск груза):
а)
/0 = 1,5; м0 = 1,1; сэо = 1; Гмех = 0,5 с.
Рис. 4. Управление скоростью вращения ротора при постоянном намагничивании РЭМ: I*, ¡*, I* — графики токов статора; скорости вращения: ш* — ротора; М — электромагнитного момента; и* — напряжения статора (в относительных единицах): а — при П * = 1 (подъем груза); б — при П * = -1 (спуск груза)
2 О 2
_«
Гб15
ЛВЕСТНИК
............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
Выводы (Summary)
На основе проведенных в работе исследований можно сделать следующие выводы.
1. При синусоидальных пульсациях индуктивностей фазных обмоток и синусоидальных токах в фазных обмотках РЭМ с числом фаз m > 3 пульсации электромагнитного момента отсутствуют.
2. Если выполнены конструктивные соотношения между зубцами статора и ротора, определенные выражениями (5) и (6), то удается получить пульсации индуктивностей электромагнитного момента близкие к синусоидальным.
3. Преобразования Ляпунова, полученные для дифференциальных уравнений напряжений РЭМ, позволяют преобразовать исходные уравнения с периодическими коэффициентами в уравнения напряжений с постоянными коэффициентами и найти составляющие вектора тока фазной обмотки: ток намагничивания и ток нагрузки. Полученные преобразования Ляпунова позволили синтезировать динамику электромагнитных процессов путем создания контуров управления током намагничивания и током нагрузки. Выполнен синтез динамики электромагнитных процессов, позволяющий формировать синусоидальные токи в фазных обмотках.
4. Для получения максимального быстродействия управления электромагнитным моментом управление РЭМ целесообразно производить при постоянстве тока намагничивания, так как постоянная времени контура тока намагничивания выше постоянной времени контура тока нагрузки.
5. Для управления скоростью вращения ротором РЭМ целесообразно применять подчиненный принцип управления, который хорошо себя зарекомендовал в приводах с классическими типами электродвигателей. Для формирования напряжений, подводимых к обмоткам статора, может быть использован стандартный преобразователь частоты, что позволит упростить и удешевить электроприводы на базе РЭМ.
СПИСОК ЛИТЕРАТУРЫ
1. Бычков М. Г. Расчетные соотношения для определения главных размеров вентильно-индукторной машины / М. Г. Бычков, Самир Риах Сусси // Электротехника. — 2000. — № 3. — С. 15-19.
2. Ваткин В. А. Разработка вентильных индукторных электромеханических систем автотранспортного назначения: автореф. дисс. на соискание ученой степени канд. техн. наук: спец. 05.09.01 «Электромеханика и электрические аппараты» / В. А. Ваткин. — М.: МЭИ (ТУ), 2007. — 20 с.
3. Ильинский Н. Ф. Проектирование вентильно-индукторных машин общепромышленного назначения / Н. Ф. Ильинский, Й. Штайнбрунн, Ю. И. Прудникова [и др.] // Вестник Московского энергетического института. Вестник МЭИ. — 2004. — № 1. — С. 37-43.
4. Фисенко В. Г. Проектирование вентильных индукторных двигателей: методическое пособие / В. Г. Фисенко, А. Н. Попов. — М.: Изд-во МЭИ, 2005. — 56 с.
5. Krishnan R. Switched reluctance motor drives: modeling, simulation, analysis, design and applications / со R. Krishnan. — CRC Press, 2001. — 432 p.
6. Miller T. J. E. Design of a synchronous reluctance motor drive / T. J. E. Miller, A. Hutton, C. Cossar, D. A. Staton // IEEE Transactions on industry applications. — 1991. — Vol. 27. — Is. 4. — Pp. 741-749. DOI: 10.1109/28.85491.
7. Красовский А. Б. Исследование пульсаций момента вентильно-индукторного двигателя при регулировании среднего значения момента в зоне малых скоростей / А. Б. Красовский // Электротехника. — 2017. — № 5. — С. 2-8.
8. Степанов А.В. Анализ пульсаций электромагнитного момента при проектировании бесконтактного двигателя постоянного тока дискового типа / А.В. Степанов, С.И. Масленникова // Наука и образование: научное издание МГТУ им. Н.Э. Баумана. — 2015 — № 10. — С. 288-301. DOI: 10.7463/1015.0820340.
9. Шабаев В. А. Анализ источников шума вентильно-индукторного двигателя / В. А. Шабаев // Электротехника. — 2005. — № 5. — С. 62-64.
10. Самосейко В. Ф. Реактивные электрические машины. Теория и проектирование/ В. Ф. Самосей-ко. — СПб.: Изд-во ГУМРФ им. адм. С. О. Макарова, 2017. — 392 с.
г
ВЕСТНИК«!
ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^
МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
11. Жарков А. Цифровое векторное управление вентильно-индукторными двигателями с независимым возбуждением / А. Жарков, А. Анучин, А. Дроздов, В. Козаченко // Компоненты и технологии. — 2004. — № 8 (43). — С. 166-170.
12. Егоров И. Н. Векторное управление вентильно-индукторными двигателями мехатронных систем / И. Н. Егоров, В. А. Шабаев // Фундаментальные исследования. — 2015. — № 12-5. — С. 891-895.
13. Белоусов И. В. Широтно-импульсные преобразователи электрической энергии: Монография / И. В. Белоусов, Ф. А. Гельвер, В. Ф. Самосейко, В. А. Хомяк. — СПб.: ФГУП «Крыловский государственный научный центр», 2019. — 228 с.
14. Самосейко В. Ф. Реактивные электрические машины с зубчатым статором и ротором. Методика проектирования. Алгоритмы управления / В. Ф. Самосейко, Ф. А. Гельвер, В. А. Хомяк, Н. А. Лазаревский. — СПб.: ФГУП «Крыловский государственный научный центр», 2016. — 197 с.
15. Самосейко В. Ф. Адаптивный алгоритм векторного управления электроприводами с асинхронными электродвигателями / В. Ф. Самосейко // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2019. — Т. 11. — № 1. — С. 156-168. DOI: 10.21821/2309- 5180-201911-1-156-168.
REFERENCES
1. Bychkov, M. G., and Samir Riakh Sussi. "Raschetnye sootnosheniya dlya opredeleniya glavnykh razmerov ventil'no-induktornoi mashiny." Elektrotekhnika 3 (2000): 15-19.
2. Vatkin, V. A. Razrabotka ventil'nykh induktornykh elektromekhanicheskikh sistem avto-transportnogo naznacheniya. Abstract of PhD diss. M.: MEI (TU), 2007.
3. Il'inskii, N. F., et al. "Proektirovanie ventil'no-induktornykh mashin obshchepromyshlennogo naznacheniya." Vestnik Moskovskogo energeticheskogo instituta. Vestnik MEI 1 (2004): 37-43.
4. Fisenko, V. G., and A. N. Popov. Proektirovanie ventil'nykh induktornykh dvigatelei: metodicheskoepo-sobie. M.: Izd-vo MEI, 2005.
5. Krishnan, R. Switched reluctance motor drives: modeling, simulation, analysis, design and applications. CRC Press, 2001.
6. Miller, T. J. E., A. Hutton, C. Cossar, and D. A. Staton. "Design of a synchronous reluctance motor drive." IEEE Transactions on industry applications 27.4 (1991): 741-749. DOI: 10.1109/28.85491.
7. Krasovskii, A. B. "Studies of torque ripple in a switched reluctance motor under a controlled average torque value in a low-velocity mode." Russian Electrical Engineering 88.5 (2017): 247-252. DOI: 10.3103/ S1068371217050066.
8. Stepanov, A. V., and S. I. Maslennikova. "Analysis of the torque ripples in designing a disk type brushless direct current motor." Science and Education of Bauman MSTU 10 (2015): 288-301. DOI: 10.7463/1015.0820340.
9. Shabaev, V. A. "Analysis of the noise sources of the switched-reluctance motor." Elektrotekhnika 5 (2005):
62-64.
10. Samoseiko, V. F. Reaktivnye elektricheskie mashiny. Teoriya iproektirovanie. SPb.: Izd-vo GUMRF im. adm. S. O. Makarova, 2017.
11. Zharkov, A., A. Anuchin, A. Drozdov, and V. Kozachenko. "Tsifrovoe vektornoe upravlenie ventil'no-induktornymi dvigatelyami s nezavisimym vozbuzhdeniem." Komponenty i tekhnologii 8(43) (2004): 166-170.
12. Egorov, I. N., and V. A. Shabaev. "Vector control valve-inductor motor mechatronic systems." Fundamental research 12-5 (2015): 891-895.
13. Belousov, I. V., F. A. Gel'ver, V. F. Samoseiko, and V. A. Khomyak. Shirotno-impul'snye preobrazovateli elektricheskoi energii: Monografiya. SPb.: FGUP «Krylovskii gosudarstvennyi nauchnyi tsentr», 2019.
14. Samoseiko, V. F., F. A. Gel'ver, V. A. Khomyak, and N. A. Lazarevskii. Reaktivnye elektricheskie mashiny s zubchatym statorom i rotorom. Metodika proektirovaniya. Algoritmy upravleniya. SPb.: FGUP «Krylovskii gosudarstvennyi nauchnyi tsentr», 2016.
15. Samosejko, Veniamin F. "The adaptive algorithm of the vector control of electrical drives with the induction motors." Vestnik Gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S. O. Makarova 11.1 (2019): 156-168. DOI: 10.21821/2309-5180-2019-11-1-156-168.
2 О 2
_Oi
Гбп
ЛВЕСТНИК
............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
ИНФОРМАЦИЯ ОБ АВТОРАХ
INFORMATION ABOUT THE AUTHORS
Самосейко Вениамин Францевич —
доктор технических наук, профессор ФГБОУ ВО «ГУМРФ имени адмирала С. О. Макарова»
198035, Российская Федерация, Санкт-Петербург,
ул. Двинская, 5/7
e-mail: [email protected],
Ширяев Эдуард Вячеславович — аспирант
Научный руководитель:
Самосейко Вениамин Францевич
ФГБОУ ВО «ГУМРФ имени адмирала
С. О. Макарова»
198035, Российская Федерация, Санкт-Петербург, ул. Двинская, 5/7
e-mail: shiryaev. edward@yandex. ru, [email protected]
Samoseiko, Veniamin F. —
Dr. of Technical Sciences, professor Admiral Makarov State University of Maritime and Inland Shipping
5/7 Dvinskaya Str., St. Petersburg, 198035, Russian Federation e-mail: [email protected], [email protected] Shiryaev, Eduard V. — Postgraduate Supervisor:
Samoseiko, Veniamin F.
Admiral Makarov State University of Maritime
and Inland Shipping
5/7 Dvinskaya Str., St. Petersburg, 198035, Russian Federation e-mail: shiryaev. edward@yandex. ru, [email protected]
Статья поступила в редакцию 27 мая 2020 г.
Received: May 27, 2020.
