Управление током статора синхронного электродвигателя с возбуждением постоянными магнитами при частотном пуске
Анатолий Коршунов
В статье рассмотрен контур регулирования тока статора синхронного двигателя с постоянными магнитами, получено условие устойчивости этого контура и исследована точность поддержания заданного тока с учетом ЭДС вращения, наводимой в обмотке статора.
Введение
В предыдущей статье автора [1] был рассмотрен частотный пуск синхронного двигателя (СД) с возбуждением постоянными магнитами при синусоидальных токах статора с неизменной амплитудой. Как известно, на ток статора оказывают существенное влияние многие факторы, например, увеличение ЭДС вращения и сопротивления обмоток вследствие повышения частоты при пуске СД, изменение величины прикладываемого напряжения. Следовательно, для поддержания тока статора неизменным необходимы специальные меры.
В электроприводе постоянного тока, построенном по принципу подчиненного управления, задачу ограничения тока и поддержания его на этом уровне в переходных режимах успешно решает контур регулирования тока. Поэтому логично использовать его и в электроприводе переменного тока. Поскольку при частотном управлении двигателями переменного тока используют обычно трехфазный автономный инвертор (АИН), построенный на транзисторах типа MOSFET или IGBT, весьма чувствительных к перегрузке по току, требования к динамическим свойствам регулятора тока статора достаточно высоки. С другой стороны, весьма жесткими оказываются и требования к ослаблению влияния возмущающих воздействий на величину тока. Как известно, при жестких требованиях к качеству управления возможно нежелательное проявление свойств широтно-импульсного управления инвертором, например, недопустимых субгармонических автоколебаний. Таким образом, необходим учет и этого фактора.
1. Математическое описание силовой части трехфазного контура регулирования тока
Рассмотрим АИН, нагруженный на обмотку статора СД, соединенную звездой (рис. 1) при следующих допущениях:
- в силу симметрии трехфазной обмотки статора ra = rb = rc = r, La = Lb = Lc = L, Mab = Mba = Mac = Mca = Mbc = Mcb = –M;
- транзисторы инвертора считаем идеальными ключами;
- в любой момент времени включены три ключа — по одному из каждого плеча;
- ЭДС вращения ea, eb, ec, наводимые в обмотке статора вращающимся ротором, образуют трехфазную симметричную систему.
Принятые допущения позволяют записать систему уравнений электрического равновесия:
Упрощения в двух последних уравнениях дают
Подставляя в уравнение (5) ток ic, выраженный из уравнения (1), и складывая с удвоенным уравнением (4), получаем
С учетом равенств
первое из которых вытекает из трехфазной симметрии ЭДС вращения, а второе записано по второму закону Кирхгофа для внешнего контура статорных обмоток СД, уравнение (6) приводится к виду
По симметрии можно записать следующие два уравнения
Разумеется, любое из уравнений (8, 9, 10) является следствием двух других, поскольку из трех внутренних контуров статора СД независимые только два.
2. Построение регулятора фазных токов статора
Состояние ключей в каждом плече инвертора (рис. 1) будем изменять только один раз в течение периода коммутации Т. Длительность проводимости верхнего ключа обозначим τ1i, а нижнего — τ2i, где i = a, b, c. Очевидно,
Применение разрабатываемых в последнее время алгоритмов векторной широтно-импульсной модуляции (ШИМ) в трехфазном АИН [2, 3], обладающих определенными преимуществами для инверторов напряжения, в рассматриваемом случае нерационально. Причина этого — сложность реализации и применение законов коммутации ключей АИН, использующих более трех (4 или 6) переключений за период коммутации. При регулировании тока гораздо проще и рациональнее изменять длительности включения транзисторных ключей АИН не программным способом, как при векторной ШИМ, а в зависимости от ошибок управления регуляторов фазных токов. Вместо рекомендуемого числа переключений за период коммутации, равного 4 или 6, лучше увеличить частоту коммутации в 4/3 или в 2 раза.
При τ1 = τ2 = Т/2 в каждом плече АИН линейные напряжения статорной обмотки uab, ubc, uca равны нулю, поскольку все три конца статорной обмотки (а, b и с) подключены к одной точке. В течение первой половины периода они соединены с «+» источника питания, а в оставшийся полупериод — с его «–».
Изменение τ1 и τ2 только в одном плече моста, например, связанном с фазой а, вызывает изменение сразу двух линейных напряжений статорной обмотки: uab и uca. При этом очевидно, что
где 1(t) = 0 при t < 0 и 1(t) = 1 при t > 0 — единичная функция.
В случае τ1а > T/2 uab представляет собой положительный прямоугольный импульс, а при τ1а < T/2 — отрицательный. С учетом (11) в правых частях дифференциальных уравнений (8–10) получаем соответственно
Из выражений (12) очевидно, что увеличение времени подключения одной из фаз обмотки статора к «+» источника питания 1(τ1j > T/2) приводит не только к появлению положительного импульса в правой части дифференциального уравнения для тока этой фазы, но и к отрицательным импульсам вдвое меньшей амплитуды в правых частях дифференциальных уравнений для токов двух других фаз. При (τ1j < T/2) полярности всех трех упомянутых выше импульсов изменяются на противоположную. Меняется и положение импульсов на временной оси. В первом случае (τ1j > T/2) они начинаются в середине периода, а во втором (τ1j < T/2) в середине периода заканчиваются.
Таким образом, для увеличения тока в j-й фазе обмотки статора следует увеличивать τ1j (уменьшать τ2j) в j-м плече инвертора, а для уменьшения тока (с учетом его знака) следует уменьшать τ1j (увеличивать τ2j).
Следовательно, регулятор тока должен в случае положительной ошибки управления, увеличивать τ1j, а при отрицательной ошибке — уменьшать τ1j.
где Iзj и ij – заданное и действительное значение тока.
Очевидно, что для надежного ограничения токов транзисторных ключей АИН в переходных режимах регулятор тока не должен быть слишком инерционным. В простейшем случае регулятор может быть пропорциональным.
Из изложенного выше ясно, что регуляторы не будут независимы, поскольку увеличение тока в одной фазе вызывает уменьшение токов в двух других фазах. Это взаимное влияние в общем случае нежелательно. Действительно, полагая равными ошибки управления всех трех контуров регулирования
и учитывая вытекающее из этого равенство
получаем из схемы силовой части рис. 1
Таким, образом, равные ошибки управления в контурах регулирования тока нарушают процесс регулирования токов статора СД.
Однако в рассматриваемом случае заданные токи: Iза (t), Iзb (t), Iзс (t) образуют симметричную трехфазную систему. Вследствие этого с учетом тождества (1) получаем
Следовательно, в рассматриваемом случае ошибки управления в контурах регулирования тока могут быть равными между собой только в случае равенства каждой из них нулю.
3. Разностные уравнения трехфазного регулятора тока при пропорциональных регуляторах и их анализ
Разностное уравнение необходимо для анализа устойчивости контура с учетом нелинейного дискретного характера управления и для определения ошибок управления.
А. Вначале для простоты положим r = 0, и будем считать постоянными ЭДС вращения
a также и заданные значения токов
Допущения (17) и (18) тем более верны, чем ниже частота заданных токов и чем быстрее протекают электромагнитные переходные процессы по сравнению с механическими.
Пренебрегая дискретностью значений ошибок управления, вычисляемых в микропроцессоре, определяемой ограниченной разрядной сеткой, будем полагать
где τ1jn – значение τ1j в n-ом периоде коммутации,
δm – граница линейной зоны широтно-импульсного модулятора, Кр – коэффициент усиления пропорционального регулятора.
В линейной зоне ШИМ, то есть при
выполняется условие
Для определенности положим
Это позволяет построить временные диаграммы линейных напряжений: uab (t), ubc (t) и uca (t) (рис. 2) и определить импульсы напряжения в правых частях уравнений (8–10) в n-ом периоде коммутации:
Структура полученных выражений позволяет предположить, что и при других соотношениях между τ1an, τ1bn, τ1cn, отличных от принятых (22), выражения импульсов (23) останутся неизменными. Это утверждение несложно проверить.
Интегрирование уравнений (8–10) с учетом принятых допущений дает выражения (24).
В полученных разностных уравнениях значения τ1jn, j = a, b, c определяются уравнением модуляции (19).
Очевидно, разностные уравнения (24) сохраняют силу и вне линейной зоны ШИМ, когда хотя бы при одном значении j
При выполнении же условия (20) контур регулирования тока можно описать системой линейных разностных уравнений
где β = EКрT/[6(L + M)δm], ? = T/(L + M).
С учетом тождества (16) система разностных уравнений (26) распадается на три независимых разностных уравнения
где α = 1 – 3β = 1 – EКрT/[2δm (L + M)].
Установившееся решение уравнений (27) имеет вид
где
Решения (28) разностных уравнений (27) устанавливаются только при их асимптотической устойчивости. Для асимптотической устойчивости необходимо и достаточно выполнения неравенства
Несложно проверить, что неравенство (30) ограничивает коэффициент усиления регулятора тока условием
Заметим, что 0 < α α 1 соответствует монотонному затуханию переходного процесса установления δj∞, а–1 < α < 0 – колебательному с частотой вдвое меньшей частоты коммутации. Особое место занимает значение α = 0, при котором переходный процесс в пределах линейной зоны ШИМ заканчивается за 1 период коммутации. Поэтому рационально с точки зрения быстродействия выбирать соответствующее значение Кр
Последнее значение Кр обеспечивает двукратный запас устойчивости контура регулирования тока, что следует из (31) и (32).
Величина установившейся ошибки при выборе Кр из условия максимального быстродействия в линейной зоне ШИМ (32) имеет значение
Из последнего выражения видно, что уменьшение периода коммутации при сохранении двукратного запаса устойчивости позволяет пропорционально уменьшать ошибки управления контура регулирования тока, вызываемые ЭДС вращения, наводимыми в статорной обмотке.
Разумеется, формула (33) справедлива в случае значения установившихся ошибок управления, удовлетворяющих условию линейности ШИМ (20)
Из последнего выражения с учетом (32) получаем минимально необходимую величину напряжения питания инвертора
где em — амплитуда ЭДС, наводимой полюсами вращающегося ротора в фазовых обмотках статора.
Согласно (33) при постоянном задающем воздействии (18) установившаяся ошибка вызвана только возмущающим воздействием (17). Причина этого в допущении о равенстве нулю активного сопротивления статорной обмотки и нулевом падении напряжения, создаваемом постоянным током на индуктивности.
При изменяющемся задающем воздействии неизбежна составляющая ошибки управления, создающая на чисто индуктивной нагрузке напряжение, необходимое для изменения тока в ней. Кроме того, в действительности и возмущающее воздействие не удовлетворяет условию (17).
Б. Учтем гармоническое изменение задающего воздействия
где φa = 0, φb = 2π/3,φc = 4π/3, и возмущающего воздействия
Поскольку при пуске СД скорость повышения частоты напряжения на его статорных обмотках ограничена механической инерционностью СД и его нагрузки, при анализе значительно более быстрых электромагнитных переходных процессов в контуре регулирования тока можно считать Ω = const.
Интегрирование уравнений (8–10) дает при r = 0 выражения (37):
Полагая работу контура регулирования тока в линейной зоне ШИМ (20), с учетом очевидных выражений
систему разностных уравнений (37) несложно привести к виду
Линейные разностные уравнения (39) имеют две составляющие установившихся решений, то есть
Первая составляющая δej∞(nT) вызвана изменением задающих воздействий (35), а вторая — δзj∞(nT) — возмущающими воздействиями (36).
При выборе Кр согласно (32), что обеспечивает окончание переходных процессов при работе в линейной зоне ШИМ в первом периоде коммутации, из (39) легко получить
Из формул (41–44) очевидно, что при ΩT, стремящемся к нулю, (ΩT?0) ошибки, вызванные изменением задающего воздействия, обращаются в 0, а ошибки, вызванные возмущающими воздействиями, стремятся к значениям, определяемым формулой (33), в которой ej = –Emsin(φj + ?), j = a, b, c. Таким образом, предыдущий результат представляет собой частный случай последнего.
Напряжение питания инвертора, необходимое для работы контура регулирования тока в линейной зоне ШИМ, определим для наиболее тяжелого режима, в котором ошибки δзj∞ и δej∞(nT) совпадают по фазе. Для этого случая необходимо выполнить неравенство
Учитывая рост ЭДС вращения при разгоне СД, требуемое напряжение питания АИН оказывается наибольшим в конце разгона. Его минимальное допустимое значение согласно (45) равно удвоенной сумме амплитуд ЭДС вращения и падения напряжения на индуктивном сопротивлении фазы статорной обмотки СД.
4. Учет влияния активного сопротивления статорных обмоток
Интегрирование дифференциальных уравнений (8–10) при r ≠ 0 и гармонических ЭДС вращения ej (t) (36) для известных значений τ1j, j = a, b, c дает
где j = a, b, c, Tc = (L+M)/r — постоянная времени статорных обмоток, φ = arctgΩTc,
C учетом (38) разностное уравнение (46) приводится к виду
Систему разностных уравнений удобно записать в векторно-матричной форме
Зависимость матрицы D(Xn) и скаляра σ(Xn) от вектора Xn нелинейна даже в линейной зоне ШИМ (20), поскольку τ1jn, зависящие от элементов вектора Xn, входят в показатель степени экспонент, являющихся нелинейными функциями. Таким образом, при учете активного сопротивления статорных обмоток контур регулирования тока представляет собой нелинейную импульсную систему.
Декомпозиции системы нелинейных разностных уравнений (48) на три независимых разностных уравнения препятствует наличие в каждом из системы трех уравнений (48) члена σ(nT), нелинейно зависящего от δa (nT), δb (nT) и δс (nT). Нелинейность разностных уравнений и невозможность их декомпозиции существенно усложняет анализ и синтез контура регулирования тока. Поэтому желательно аппроксимировать нелинейные разностные уравнения контура линейными уравнениями. Один из вариантов заключается в пренебрежении активным сопротивлением обмоток статора, что и было сделано выше. Другой подход основан на линеаризации экспоненциальной функции
Полагая работу контура регулирования тока в линейной зоне ШИМ (20), в результате подстановки (50) с учетом (16) получаем
Выражения (51, 52) свидетельствуют о декомпозиции и линеаризации системы разностных уравнений (48). Подстановка (52) в (48) дает
Из выражения (54) находим критическое значение коэффициента усиления регуляторов фазных токов, соответствующее границе устойчивости контура (α1 = –1):
Сравнивая (55) с (31), видим, что с учетом активного сопротивления обмоток статора критический коэффициент усиления регулятора оказывается в (1 + еT/Tc)/2 ≈ 1 + rT/[2(L + M)] раз больше.
Соответствующее α1 = 0 значение Кр, полученное из формулы (54), совпадает со значением, определяемым формулой (32) и найденным из условия α = 0 в уравнении (27). Выбрав указанное значение Кр, соответствующее максимальному быстродействию контура в линейной зоне ШИМ, получаем установившийся режим изменения ошибок управления фазных регуляторов
В наиболее неблагоприятном случае θ ≈ φ и согласно выражению (56) составляющие ошибки управления, вызванные задающим и возмущающим воздействиями, складываются. Поскольку в реальных случаях выполняются неравенства
из выражения (56) можно найти максимальную ошибку управления для фазовых контуров регулирования тока. Приравняв ее значению δm/Кр, соответствующему границе линейности ШИМ, несложно получить неравенство, определяющее выбор напряжения питания инвертора
Принимая во внимание малость r по сравнению с Ω(L+M) при максимальной частоте Ω и векторное суммирование падений напряжения на активном и индуктивном сопротивлениях статорной обмотки, величина Е, определенная с учетом активного сопротивления (57), незначительно отличается от значения (45), определяемого без учета активного сопротивления.
Из выражения (56) очевидно, что при уменьшении периода коммутации Т и сохранении неизменного запаса устойчивости ошибка управления контура регулирования тока уменьшается практически пропорционально, поскольку
5. Повышение точности контура регулирования тока
Если по энергетическим соображениям выбран период коммутации, недостаточный для получения допустимых ошибок управления, для повышения точности контура регулирования тока можно использовать принцип комбинированного управления, сочетающий управление по отклонению с управлением по задающему и возмущающему воздействию.
Для полного устранения составляющей ошибки, вызванной задающим воздействием, необходимо подавать на входы однофазных контуров регулирования тока задающие воздействия Iзj (t) (35), увеличенные на значение соответствующей ошибки. Очевидно, что, пренебрегая для простоты активным сопротивлением обмоток статора (r = 0), можно определять амплитуду добавки к задающему воздействию по формуле (42), а ее фазовый сдвиг по отношению к задающему воздействию принять равным ?/2 в сторону опережения. Физическая причина опережающего фазового сдвига на ?/2 состоит в том, что ток в индуктивности отстает по фазе от напряжения на ней на ?/2.
Для компенсации ошибки, вызываемой возмущающим воздействием, необходима полная информация о нем, то есть необходимо измерять ЭДС, наводимые в обмотках статора полюсами вращающегося ротора. Добавка к задающему воздействию при r = 0 должна иметь амплитуду, определяемую выражением (43), а по фазе должна практически совпадать с ЭДС вращения. Совпадение фаз объясняется тем, что для компенсации действия ЭДС вращения напряжение на обмотке необходимо повысить на величину этой ЭДС. Однако компенсация ошибки, вызванной возмущающим воздействием, может привести к образованию положительной обратной связи с коэффициентом больше критического и к самовозбуждению привода, поскольку ток статора влияет на скорость ротора и, следовательно, на величину ЭДС вращения. Указанная опасность в условиях нестабильных параметров и нелинейности системы допускает только частичную компенсацию ошибок управления, вызываемых ЭДС вращения.
6. Проверка полученных результатов
Проверка проводилась на математической модели контура регулирования тока, построенной в среде Matlab 6.5, Simulink 5, SimPower Systems, представленной на рис. 3.
С целью количественной проверки результатов в модели инвертора использованы идеальные ключи (Ideal Switch), а не модели MOSFET или IGBT. Индуктивную нагрузку инвертора моделирует блок 3-Phase RLC Branch1. Генератор заданных значений трехфазного тока (Iзj (t), j = a, b, c) построен на программируемом трехфазном генераторе напряжения (3-Phase Programmable Voltage Source 1), трехфазной активной нагрузке (3-Phase Series Branch) и трехфазном измерителе (3 – Phase V1 Measurement), моделирующем датчик трехфазного тока. Демультиплексоры Demux 1 и Demux 2 вместе с тремя сумматорами представляют трехфазный измеритель ошибок управления в контуре регулирования тока. Пропорциональные регуляторы тока моделируются тремя усилителями Gain 1, 2, и 3. Широтно-импульсные преобразователи, описываемые уравнением (19), представлены в модели фиксаторами нулевого порядка (Zero-Order Hold 1, 2 и 3) и блоками Subsystem 1, 2 и 3. Последние представляют собой широтноимпульсные модуляторы с синхронными генераторами линейно возрастающего опорного сигнала.
В модели установлены следующие гипотетические значения параметров:
Результаты моделирования при практически постоянных, полученных при очень низкой частоте Ω генератора 3-Phase Programmable Voltage, рассмотрены ниже. На рис. 4 представлен переходный процесс установления постоянных токов при значении
обеспечивающем наибольшее быстродействие контура, то есть окончание переходного процесса в первом периоде коммутации. Для протекания переходного процесса в линейной зоне ШИМ задано значение Im√3/2 = 0,56 А. Увеличение Кр вдвое соответствует границе устойчивости контура регулирования тока. При тех же значениях Iзj переходный процесс начинается в зоне насыщения ШИМ. Со второго периода коммутации наблюдаются незатухающие колебания половинной частоты в линейной зоне ШИМ (рис. 5), соответствующие значению α = –1. Однако при критическом значении Кр контур регулирования тока является «негрубой» системой. Это означает, что достаточно ничтожного изменения любого из параметров, влияющих на его значение (31), чтобы контур регулирования тока сместился с колебательной границы устойчивости в область устойчивости или в область неустойчивости. В первом случае колебания со временем затухнут. Во втором же случае амплитуда колебаний возрастает. При вхождении колебаний в зону насыщения ШИМ рост амплитуды колебаний прекращается и устанавливаются автоколебания определенной амплитуды.
На рис. 6 представлен переходный процесс установления субгармонических автоколебаний при Кр = 3,4 (на 0,2 больше критического). Из рис. 6 видно, что частоты наблюдаемых автоколебаний оказываются разными в разных фазах. В фазах «b» и «c» наблюдаются автоколебания половинной частоты, а в фазе «a» автоколебания имеют частоту коммутации. Интересно, что в режиме автоколебаний постоянные составляющие токов ia (t), ib (t), ic (t) отличаются от заданных значений Iза, Iзb, Iзc. При учете сопротивления статорных обмоток r = 6 Ом при том же значении Кр = 3,4 устанавливаются автоколебания половинной частоты во всех фазах (рис. 7). Происходят они в линейной зоне ШИМ. Объясняется это уже отмеченной выше нелинейностью контура регулирования тока при отличном от нуля активном сопротивлении обмоток статора.
Результаты моделирования процесса отработки гармонических токов (35) при нулевых начальных условиях, r = 0, Кр = 1,6, частоте Ω = 2? Ѕ 108 рад/с (108 Гц) и той же величине Im, при которой исследовали отработку постоянных заданных токов (Im√3/2 = 0,56 А) в линейной зоне, представлены на рис. 8. Контур регулирования тока выходит на установившийся режим уже в конце первого периода коммутации. Величины Ω и (L + M) взяты не из реальных условий, а из соображений наглядности временных диаграмм. Из рис. 8 видно, что отрабатываемые токи ij (t) отличаются от заданных токов Iзj (t), j = a, b, c практически только по фазе. Это вызывает ошибки управления, показанные на рис. 9.
Для компенсации ошибок управления, вызываемых изменяющимися задающими воздействиями (35), достаточно прибавить к задающим воздействиям их производные с соответствующими коэффициентами. Дифференцируя (35), получаем
и сравнивая с выражением (41), записанным для реальных условий (Ωt << 1),
находим, что выходной сигнал дифференциатора надо умножить на коэффициент, равный периоду коммутации Т.
Математическая модель контура регулирования тока с компенсацией ошибки, вызываемой изменяющимся задающим воздействием, представлена на рис. 10. Дифференцирование задающего воздействия осуществляют блоки Transfer Function 1, 2, 3 c передаточной функцией 0,00025S/(0,000001S + 1), выходные сигналы которых складываются с задающими воздействиями.
На рис. 11 представлен переходный процесс отработки тех же, что и выше, гармонических задающих воздействий при компенсации ошибки их отработки. Из рис. 12, на котором показаны кривые изменения ошибок управления (13), видно, что после окончания переходного процесса ошибки управления содержат только пульсационную составляющую, как и ожидалось.
Выводы
- Трехфазный контур регулирования тока при чисто индуктивном сопротивлении обмоток статора синхронного двигателя с постоянными магнитами на роторе в линейной зоне широтно-импульсного модулятора описывается тремя независимыми линейными разностными уравнениями первого порядка.
- Учет активного сопротивления обмоток делает разностные уравнения нелинейными и зависимыми.
- Превышение критического коэффициента усиления регуляторов тока вызывает неустойчивость контура регулирования тока, что приводит к субгармоническим автоколебаниям.
- При чисто индуктивном сопротивлении статорных обмоток субгармонические автоколебания происходят в зоне насыщения ШИМ, а при активно-индуктивном сопротивлении могут существовать и в линейной зоне ШИМ.
- При заданном запасе устойчивости в контуре регулирования тока его ошибки управления пропорциональны периоду коммутации автономного инвертора.
- Точность контура регулирования тока можно существенно повысить за счет использования принципа комбинированного управления.
- Напряжение питания автономного инвертора, обеспечивающее отработку заданных синусоидальных токов должно в 2 раза и более превышать сумму амплитуд ЭДС вращения и падения напряжения, создаваемого заданным током на сопротивлении обмотки статора синхронного двигателя.
Литература
- Коршунов А. Равноускоренный частотный пуск синхронного двигателя с постоянными магнитами на роторе // Силовая электроника. 2007. № 1.
- Изосимов Д. Б., Байда С. В. Алгоритмы векторной широтно-импульсной модуляции трехфазного автономного инвертора // Электротехника. 2004. № 4.
- Крутиков К. К., Рожков В. В., Петрухин Ю. В. Анализ и моделирование симплексных алгоритмов управления трехфазным автономным инвертором для частотного асинхронного электропривода // Электротехника. 2006. № 5.