Упрощенная математическая модель синхронного электродвигателя с возбуждением постоянными магнитами
Анатолий Коршунов
Для синхронного двигателя с постоянными магнитами на роторе выведены дифференциальные уравнения, в которых учитываются главные факторы, определяющие свойства машины, — зависимость электромагнитных процессов и вращающего момента от скорости и углового положения ротора.
Введение
Благодаря своим высоким эксплуатационным характеристикам синхронный двигатель (СД) с постоянными магнитами на роторе является перспективной электрической машиной для приводов малой и средней мощности: от стиральной машины до металлорежущего станка и электровоза [1]. Отсутствие обмотки возбуждения и успокоительной обмотки несколько упрощает процессы в СД и их описание.
Однако, как известно, сложность явлений, происходящих в электрических машинах переменного тока при переходных процессах, делает их математическое описание и исследование без ряда упрощающих допущений практически невозможным [2]. Стремление учесть главные факторы, определяющие свойства машины, и пренебречь второстепенными факторами приводит к рассмотрению идеализированной электрической машины. Такая машина характеризуется:
- отсутствием насыщения, гистерезиса и вихревых токов в магнитной цепи;
- отсутствием вытеснения тока в проводниках обмоток;
- синусоидальным распределением в воздушном зазоре магнитодвижущей силы и магнитной индукции;
- полной симметрией обмоток статора;
- независимостью индуктивностей рассеяния обмоток статора от положения ротора.
Обычно идеализируют и источник питания, считая его выходное сопротивление нулевым, а ЭДС — образующими трехфазную симметричную систему.
Имея конечной целью использование математической модели синхронного двигателя (СД) в анализе и синтезе замкнутых частотно управляемых электроприводов (ЭП), желательно получить наиболее простую и удобную модель. Для этого примем дополнительные упрощающие допущения. Учитывая необходимое для защиты силовых ключей автономно инвертора ограничение тока СД, будем считать, что статорная обмотка питается от трехфазного генератора тока. При наличии в ЭП контура регулирования тока это допущение весьма близко к реальности и позволяет рассматривать электромеханические процессы независимо от электромагнитных процессов. Кроме этого, пренебрежем индуктивностью рассеяния обмоток статора СД и учтем отсутствие успокоительной обмотки.
Учитывая очень большую коэрцетивную силу современных редкоземельных магнитов (Nd–Fe–B или Sm–Co), будем считать индукцию магнитного поля, создаваемого постоянными магнитами ротора в зазоре СД, неизменной.
Электромагнитный вращающий момент СД
Рассмотрим поперечное сечение СД с одной парой полюсов (2р = 2), представленное на рис. 1. Буквами AX, BY, CZ обозначены секторы фазных обмоток статора. Фазовые токи образуют трехфазную симметричную систему:
Направления токов указаны для момента t = 0. Ось магнитного поля статора совпадает с осью обмотки фазы А, причем направление силовых линий поля соответствует указанным на рис. 1 магнитным полюсам статора Nc и Sc. Нетрудно понять из рис. 1 и формул (1), что поле статора вращается по часовой стрелке.
Для определения амплитуды электромагнитного вращающего момента СД положим, что в момент t = 0 ось поля ротора совпадает с плоскостью обмотки фазы А (то есть перпендикулярна к оси обмотки), а направление его силовых линий соответствует показанным на рис. 1 полюсам поля ротора: Np и Sp.
По закону Ампера, на проводники обмотки статора действуют силы, направленные согласно правилу левой руки в сторону, противоположную направлению вращения поля статора. Таким образом, на статор СД со стороны его ротора действует вращающий момент, направленный против направления вращения поля статора.
Согласно третьему закону Ньютона, точно такой же вращающий момент действует со стороны статора СД на его ротор, но направлен он в противоположную сторону, то есть в сторону вращения поля статора.
Для вычисления вращающего момента примем, что ток статорной обмотки не сосредоточен в ее проводниках, а распределен вдоль окружности статора, как показано на рис. 2.
Плотность тока в секторах фаз А, В и С определяется выражением:
где ij —мгновенное значение фазного тока, W — число витков фазовой обмотки, α — относительная длина дуги окружности в радианах, iA = –ix = Im, iB = –iy = Im cos [–(2π/3)], ic = = –iz = Imcos [–(4π/3)]. Элементарный момент dМ, создаваемый элементарным током статора di=(di/da) da, определяется выражением
где В (α)—индукция в данной точке воздушного зазора, R — расстояние от оси машины до проводников статора, l — длина активной части ротора.
Отсчитывая α от плоскости обмотки фазы А в направлении вращения поля, можно записать выражение для величины вращающего момента
Выполняя интегрирование по фазным секторам обмотки, на которых плотность тока статора постоянна, с учетом
где Вm — амплитуда магнитной индукции, создаваемой в зазоре полюсами ротора, получаем
Рассмотрим далее общий случай, когда в момент t = 0 токи в фазах имеют значения
что означает поворот поля статора на угол φ в направлении его вращения (φ > 0). Ротор сместим на угол
в том же направлении (
> 0). В этом случае полюсы статора и ротора сдвинуты уже не на π/2, как на рис. 1, а на угол, равный π/2 + φ –
.
На рис. 3 представлены законы изменения плотности тока статора в относительных единицах (то есть без множителя (3ImW)/π и магнитной индукции, создаваемой ротором в воздушном зазоре СД). Аналогично предыдущему определяется вращающий электромагнитный момент двигателя, что с учетом (6) дает (8).
где Мm = (9/π) BmWRlIm =CмIm, Cм = (9/π) WBmRl.
При смещении точки отсчета угла поворота ротора на π/2 в сторону вращения поля статора, то есть отсчитывая
от оси обмотки фазы А, в выражение (8) следует подставить
+ (π/2), что дает
где φ –
— угловое смещение оси полюса ротора Sp (Np) относительно оси полюса статора Nc (Sc), отсчитываемое от Sp до Nc в направлении вращения поля. Очевидно, что при отставании поля ротора от поля статора на угол до 180° на ротор действует вращающий момент, а при опережении на такой же угол – тормозящий момент. При совпадении осей полюсов Sp и Nc (Sc и Np) электромагнитный момент равен нулю. Как известно, вращающий и тормозящий моменты обращаются в нуль и при φ –
= ±π, что в случае одной пары полюсов соответствует одному и тому же положению ротора относительно статора. Однако это положение неустойчиво, поскольку любое отклонение от него вызывает увеличение отклонения под действием тормозящего или разгоняющего момента. Разумеется, приведенные рассуждения относятся к случаю отсутствия момента нагрузки.
ЭДС, наводимая постоянными магнитами вращающегося ротора в обмотке статора СД
При вращении ротора в каждом проводнике обмотки наводится ЭДС вращения
где Bi — магнитная индукция, создаваемая ротором в месте i-го проводника обмотки, ν — линейная скорость проводника относительно поля ротора.
Мгновенное значение полной ЭДС, наводимой в проводнике фазовой обмотки, равно сумме ЭДС, наводимых в данный момент в каждом проводнике обмотки.
Определим положение ротора углом
между плоскостью обмотки фазы А и осью S полюса ротора, отсчитываемым в направлении вращения поля. Распределение магнитной индукции в зазоре у сектора фазы А обмотки статора можно описать выражением
где α — угол, отсчитываемый от плоскости обмотки фазы А в направлении вращения поля статора.
Полагая угловое расстояние между соседними проводниками обмотки (рис. 1) равным
получаем величину суммарной ЭДС в фазе А:
Множитель 2 перед суммой объясняется суммированием ЭДС проводников только в одном секторе фазной обмотки.
Ряд косинусов легко сворачивается, если воспользоваться формулой Эйлера
Подставляя выражение (14) в ряд (13), получаем с учетом Δα << 1:
Подстановка полученного результата (15) в формулу (13) дает
Выразив линейную скорость υ через угловую (d
/dt), что дает
С учетом пространственного сдвига фазных обмоток статора СД имеем
Для определения направления еА рассмотрим рис. 1 при
= 0. При вращении ротора в сторону вращения поля статора, то есть по часовой стрелке (ч. с.), проводники ротора вращаются относительно поля ротора в противоположную сторону, то есть против ч. с. Согласно правилу правой руки, еА направлена в сторону, противоположную указанному на рис. 1 направлению тока фазы А, совпадающему с выбранным положительным направлением тока в этой фазе статорной обмотки СД. По симметрии направим ЭДС остальных фаз: еВ и еС противоположно соответствующим фазным токам iB и iC.
Учитывая симметрию фазовых обмоток и угол сдвига между их осями 2π/3, можно считать взаимоиндуктивности одинаковыми и равными –М.
Уравнения электрического равновесия для фазовых обмоток статора СД (рис. 4) имеют вид
где L и r — индуктивность и активное сопротивление фазы статора, uА, uВ, uС — фазные напряжения.
Установившийся режим работы СД
Положим электромагнитный момент СД равным постоянному моменту нагрузки:
а фазные токи — определенными формулами (1).
Из выражения (21) получаем
Учитывая, что в выражениях (18), (19) угол поворота ротора отсчитывается от плоскости обмотки фазы А, а в формуле (22) — от оси этой обмотки, подстановкой в эти выражения
+ π/2 вместо
, а также
где ωр — скорость вращения ротора, равная в установившемся режиме скорости вращения поля статора ω (ωр = ω), получаем
Подставив выражения (1), (22) и (24) в систему уравнений (20), находим фазовые напряжения
где Em= Ceω — амплитуда ЭДС вращения.
Переписав первое из уравнений системы (25) c помощью элементарных тригонометрических преобразований в виде:
легко построить векторную диаграмму для фазового напряжения, представленную на рис. 5. Из выражения (26) и рис. 5 следует, что фаза ЭДС вращения относительно фазового тока равна π/2 –
0, где
0 — угол отставания ротора от вращающегося поля статора. Следовательно, на векторной диаграмме (см. рис. 5) положение вектора Em определяет угол
0, отсчитываемый от направления, перпендикулярного вектору тока Im, по часовой стрелке.
Из рис. 5 очевидно, что уменьшение отставания поля ротора от поля статора
0, необходимое при неизменных фазовых токах в случае уменьшения момента нагрузки Mн, требует увеличения фазных напряжений. Поэтому рациональнее уменьшать фазные токи при сохранении
0, близким к π/2. Это обеспечит меньшие фазные напряжения и меньший сдвиг по фазе между током и напряжением φ (рис. 5) и, следовательно, высокий коэффициент мощности cosφ.
Для проверки баланса мощностей умножим фазовое напряжение uА (ωt) на ток iA (ωt) и вычислим среднюю мощность фазы А, что дает:
Суммируя мощности трех фаз, получаем потребляемую СД электрическую мощность:
Отдаваемая механическая мощность с учетом формул (6), (18), (21), (22) определяется выражением:
Из сравнения выражений (28) и (29) следует, что потребляемая СД электрическая мощность превышает отдаваемую механическую мощность на величину потерь в активном сопротивлении обмотки статора. Это при принятых допущениях соответствует закону сохранения энергии.
Дифференциальные уравнения СД, управляемого от источника напряжения
При управлении от источника напряжения даже при наличии контуров регулирования фазных токов в переходных режимах условие (1) может не выполняться, что следует учитывать.
Положим известными мгновенные значения токов в фазах обмотки статора: iA, iB и iC и угол сдвига оси поля ротора относительно плоскости обмотки фазы A статора
, отсчитываемый в направлении вращения поля статора.
Подставляя в формулу (4) значения плотностей тока в секторах фаз (2):
с учетом закона изменения индукции в зазоре (11) (рис. 6):
получаем (см. ниже для Мэм).
Можно получить выражение для Мэм в другом виде, используя формулу (8). Для этого представим мгновенные значения токов фаз СД в виде (7):
Очевидно, система уравнений не противоречит первому уравнению Кирхгофа для статорных обмоток СД, соединенных звездой (рис. 4).
Представив второе уравнение системы (7) в виде:
и решив его совместно с первым уравнением, получаем
Из полученных формул следует, что любые токи статора СД, удовлетворяющие условию (31), можно рассматривать как мгновенные значения трехфазной системы токов (7), амплитуды которых определяются формулой (32), а текущая фаза — формулой (33). Подставляя Im и φ в формулу (8), получаем мгновенное значение Мэм.
Воспользовавшись дифференциальными уравнениями непрерывной модели автономного инвертора, работающего на статорную обмотку СД [2], и используя уравнения электромагнитного момента (30) и ЭДС вращения (18, 19), получаем систему дифференциальных уравнений электрического и механического равновесия (34), где LΣ = L + M, Iзj, j = A, B, C — заданные значения фазных токов, Кр — коэффициент усиления пропорциональных регуляторов фазных токов, L, M и r — индуктивность, модуль взаимоиндуктивности и активное сопротивление фазных обмоток статора, Е — напряжение питания инвертора, δm — граница линейной зоны широтно-импульсных модуляторов (ШИМ), управляющих ключами инвертора, sat (x) — функция насыщения, J — приведенный к ротору СД суммарный момент инерции двигателя и его нагрузки, Мн — момент нагрузки, в общем случае зависящий от угла поворота ротора СД
и его скорости d
/dt,
отсчитывается от оси обмотки фазы А до оси полюса S ротора.
В линейной зоне широтно-импульсного модулятора система уравнений электрического равновесия допускает декомпозицию и полная система (32) имеет вид (35).
Уравнения (34) и (35) получены в предположении одной пары полюсов у СД (2р = 2). В случае большего числа пар полюсов угол
измеряется в электрических радианах:
эл = р
. В этом случае можно сохранить уравнения неизменными, полагая в них
измеряемым в электрических радианах, но несколько изменив значения коэффициентов в дифференциальных уравнениях. Коэффициент Се в первых трех уравнениях следует разделить на число пар полюсов р, а суммарный момент инерции J в последнем уравнении следует разделить на р, оставив без изменения Се в его правой части.
Экспериментальное определение параметров СД с постоянными магнитами технически проще выполнить в генераторном режиме, особенно при значительной его мощности. Это позволяет обойтись без мощного источника переменного напряжения регулируемой частоты и напряжения [4], а также снимает проблему нагрузки двигателя. Из опыта холостого хода легко определить значение Ce, а из опыта короткого замыкания — значение LΣ. Величину активного сопротивления фазной обмотки статора r можно измерить непосредственно, например, методом амперметра и вольтметра.
Выводы
Учет только основных особенностей синхронного двигателя с постоянными магнитами, определяющих его свойства, позволил получить относительно простые дифференциальные уравнения электромагнитных и механических процессов в нем.
Литература
- Панкратов В. Вентильный электропривод: от стиральной машины до металлорежущего станка и электровоза // Электронные компоненты. 2007. № 2.
- Сипайлов Г. А., Кононенко Е. В., Хорьков К. А. Электрические машины (специальный курс). Высшая школа, М.: 1987.
- Коршунов А. И. Непрерывная модель трехфазного контура регулирования тока синхронного двигателя с постоянными магнитами на роторе // Силовая электроника. 2007. № 4.
- Дробкин Б. З., Карзунов Р. А., Крутяков Е. А., Павлов П. А., Пронин М. В. Высоковольтные тиристорные преобразователи частоты ОАО «Электросила» // Электротехника. 2003. № 5.