CC BY
61
УПРАВЛЕНИЕ ЭЛЕКТРОПРИВОДОМ И СИНТЕЗ ЭЛЕКТРОТЕХНИЧЕСКИХ И ЭЛЕКТРОЭНЕРГЕТИЧЕСКИХ
СИСТЕМ
УДК 62-83:004
Е.В. Александров, д-р. техн. наук, проф., 8-910-940-03-14, [email protected] (Россия, Тула, ТулГУ), С А. Савичев, асп., 8-920-275-57-09, andrev0931 @vandex.ru (Россия, Москва, ИПУ им. В.А. Трапезникова РАН), A.A. Иванов, асп., 8-920-740-13-46, andrey0931 @yandex.ru (Россия, Москва, ИПУ им. В.А. Трапезникова РАН)
ОСОБЕННОСТИ СИНТЕЗА ЦИФРОВЫХ ЭЛЕКТРОПРИВОДОВ ПЕРЕМЕННОГО ТОКА НА БАЗЕ АСИНХРОННОГО ДВИГАТЕЛЯ
Сформулированы основные положения по организации цифрового управления и общие принципы разработки алгоритмов управления, рассмотрены математические основы анализа дискретных систем, приведен пример синтеза цифрового управления асинхронным двигателем
Кчючевые слова: цифровое управление, алгоритм, дискретная систем, анализ, синтез, разностные уравнения, асинхронный двигатель, моментный привод.
1. Основные положения по организации цифрового управления асинхронным двигателем.
С точки зрения синтеза систем управления переход к цифровой технике характеризует собой переход к дискретным системам управления с квантованием по уровню и по времени. В отличие от непрерывных сигналов сигналы, поступающие в цифровое управляющее устройство и выходящие из него, принимают дискретные значения в дискретные моменты времени [1]. Квантование по времени представляет собой периодический процесс, характеризующийся тактом дискретности или периодом дискретизации Г. Без учета этих особенностей цифровой системы управления
147
непосредственная реализация в ней аналоговых алгоритмов приводит к понижению статической точности системы, появлению колебательной компоненты с амплитудой, пропорциональной периоду дискретизации, возникновению биения.
Методы проектирования цифровых систем управления существенно отличаются от классических методов, применяемых при анализе и синтезе систем непрерывного типа. Это обусловлено следующими особенностями:
1) основой математического аппарата проектирования цифровых систем управления являются разностные уравнения, которые заменяют дифференциальные уравнения, описывающие непрерывные системы;
2) анализ и синтез системы осуществляется на периоде дискретизации Т, что позволяет осуществлять как декомпозицию управления по темпам протекания процессов, так и упрощение уравнений за счет квазипостоянства переменных на этом периоде;
3) в системе имеется память, в которой хранятся значения векторов состояния и управления в предыдущие моменты времени, которые могут быть использованы при решении задачи управления;
4) синтезированные цифровые алгоритмы реализуются чаще всего на микроконтроллерах, которые имеют свои ограничения по длительности цикла вычисления и вычислительным возможностям.
Базовой величиной при организации вычисления в микропроцессорном контроллере является длительность цикла вычисления (такта), тесно связанная с моментами времени измерений и вычислительными возможностями микропроцессора. В дальнейших рассуждениях будем предполагать, что используется синхронный принцип управления, и микропроцессор в течение одного такта, с базовой длительностью цикла, равной длительности периода ШИМ, А совпадающего с периодом дискретизации Т, обеспечивает решение поставленной задачи управления, т.е. хватает вычислительных возможностей, чтобы в течение одного такта сформировать команду управления на исполнительное устройство. Следует отметить, что при этом решение задачи управления принципиально не может быть выполнено быстрее, чем за два периода вычислений управления. На первом шаге вычислений (к, к +1) цифровой регулятор по имеющейся информации о переменных и заданиях вычисляет такие команды управления на исполнительное устройство, при которых к концу следующего (к +1, к + 2) шага гарантируется решение задачи управления. На втором шаге на исполнительное устройство подаются эти команды управления. Указанное запаздывание, если оно является существенным для данной системы, должно быть принято во внимание и соответствующим образом компенсировано. Временная диаграмма организации вычислений приведена на рис. 1.
1 .Реализация управления, рассчитанного для шага [к -1, к] 2.Расчет команд управления для шага [к, к +1] 3 .Получение задания для шага [к +1, к + 2]
4.Измерения текущих значений переменных
5.Вспомогательные вычисления
1 .Реализация управления, рассчитанного для шага [к, к +1] 2.Расчет команд управления для шага [к +1, к + 2]
3 .Получение задания для шага [к + 2, к + 3]
4.Измерения текущих значений переменных.
5.Вспомогательные вычисления
1 .Реализация управления, рассчитанного для шага [к +1, к + 2] 2.Расчет команд управления для шага [к + 2, к + 3]
3 .Получение задания для шага [к + 3, к + 4].
4.Измерения текущих значений переменных
5.Вспомогательные вычисления
Рис.1. Временная диаграмма организации вычислений
С учетом вышеизложенного для синтеза цифровых систем управления автоматизированными асинхронными ЭП с высокими показателями качества необходима разработка специальных методов проектирования таких систем, учитывающих и активно использующих вышеуказанные особенности цифровых систем управления
2. Математические основы анализа дискретных систем управления
При анализе линейных дискретных систем используются разностные ура-внения, дискретное преобразование Лапласа и 7-преобразования
Понятие о разностных уравнениях
В дискретных системах в качестве переменной функции выступает реше-тчатая функция х(кТ) = х(к).
Скорость изменения решетчатой функции характеризуется ее первой разностью (аналог первой производной для непрерывных систем)
Ах[ к ] = х[ к +1] - х[ к ]. (1)
Вторая разность (аналог второй производной для непрерывных систем)
А 2 х[ к ] = х[ к + 2] - 2 х[ к +1] + х[ к ], (2)
т.е. вторая разность есть разность первых разностей. Подставляя в (2) выражение для первых разностей получим
А2х[к] = х[к + 2] - 2х[к +1] + х(к), (3)
аналогично получим разность m -порядка
Ашх[к] = Аш-1х[к + 1] - АШ-
1х(к). (4)
Выражая конечные разности через дискреты функции, получим
Ш
. иГМ-Ш
АШх[к] = X (-1)гСгШх[к + ш - г], (5)
г=1
ш ш!
где С! =-:--коэффициент бинома Ньютона.
й(т -1)!
Если линейная стационарная непрерывная система описывается дифференциальным уравнением
ЛПО + ^ + ...о„х(0=ь0 /1+* + .Л (О (6)
йгп йгп-1 йг8 йг8-1
то соответствующая ей дискретная система описывается разностным уравнением (уравнением в конечных разностях) с постоянными коэффициентами
а0 Ах[к ] + аф п-1 х[г ] +... + апх[к ] = Ь0 А [к ] + ^А 5-1/[к ] +... + Ь8/[к ] (7)
с начальными условиями
Ах[0] = х0, Ах[0] = х1 ..., А п-1х[0] = хп-1.
Выражая конечные разности через дискретные значения функций в соответствии с формулой (4) уравнение (7) можно записать в форме рекуррентного разностного уравнения
С0х[к + г] + С1х[к + г -1] +... + Сг-1х[к +1] + Сгх[к] =
= й0 / [к + V ] + й1 / [к + V -1] +... + йу-1 / [к +1] + й1 / [к ] (8) при начальных условиях Ах[0] = х0, Ах[1] = х1 ...,х[к -1] = хк-1, где г £ п, V £ 8, г > V, а коэффициенты С1, V = 0,1,...,г, й0 выражаются через коэффициенты уравнения (7) аг, г = 0,1,...,п, Ьг, г = 0,1,...,8.
Дискретное преобразование Лапласа
На рис. 2 приведена структурная схема линейной дискретной системы.
Из структуры рис. 2 можно записать
¥
е (г) = Хе(кТ )Б (г - кТ). (9)
к=0
Так как идеальное импульсное звено (ИИЗ) генерирует последовательность мгновенных единичных импульсов (£-функций), то сигнал на выходе ИИЗ можно представить в виде
¥
е* (г) = Т$(е - кТ)е(кТ). (10)
к=0
150
Рис.2. Структурная схема линейной дискретной системы
В силу свойств ИИЗ в эквивалентной схеме имеется формирующее звено с весовой функцией w _ S(t), которой соответствует передаточная функция
Шф (р) = L{S (t)} = S (P). (11)
Для последовательности прямоугольных импульсов скважностью Y = Ти /T (рис. 3) весовая функция определяется выражением
S (t) = Ku [1(t) - 1(t - YT). (12)
Преобразовав обе части соотношения (12) по Лапласу, получим передаточную функцию ФЗ
s ( Р )
_ Kи (1 - e-YTp)
(13)
Рис.3. Формирование прямоугольных импульсов скважностью у
Учитывая (13), для Ки = 1 и у = 1 получим передаточную функцию ФЗ, которая называется экстраполятором нулевого порядка
S (p)
1 - e
- JTP
Р
(14)
Для дискретных значений выходной величины х (t) можно запи-
сать
X (t) _ I х(kR)8(t - kT). k _ 0
oo
Совокупность гармонических функций, на которые может быть
разложен 8(? - кТ), равна е-кТ. Поэтому спектр
¥
X* (уш) = I х(кТ)е-кТ . (16)
к=0
Заменив в (16) уш = р, получим формулу, определяющую прямое дискретное преобразование Лапласа для решетчатой функции
¥ г
X* (р) = I х(кТ)е-рТ . (17)
к=0
Поскольку частота квантования шт = 2 п / Т, получим е - у (ш+кшт) кТ = е - кТе - к 2 п / Т = е - ую кТ,
т.е. вся плоскость р для дискретной переменной разделена на полосы вдоль мнимой оси, одинаковые по размеру шт (рис. 4)
Рис.4. Главная полоса плоскости р, соответствующей дискретному преобразованию Лапласа
2 - преобразования
Аналогично, как с помощью преобразования Лапласа линейное дифференциальное управление преобразуется в алгебраическое, линейное разностное уравнение можно преобразовать в алгебраическое с помощью 2 - преобразований.
Обозначим в формуле (17) 2 = ерТ Подставляя р = уш и рассматривая диапазон частот (- п / Т) £ ш £ п / Т, получаем формулу 2 -преобразования
¥ 1
X(2) = I х(кТ>-к, (18)
к=о
е-уп £ 2 £ еуп. (20)
На плоскости 2 отрезок мнимой оси главной полосы плоскости р
(рис. 5) отображается в окружность единичного радиуса на плоскости z. z - преобразование применимо только для дискретной системы по времени функции / [ к ] = / (кТ), где Т - период квантования. Учитывая, что z -
комплексная переменная, обратное преобразование осуществляется по формуле
/ [ к ] = / Г (z )2к-1dz. (20)
2 щ
Рис.5.Область устойчивости линейной дискретной системы
на плоскости
3. Общие принципы разработки алгоритмов управления
При синтезе алгоритмов цифрового управления использованы следующие принципы:
- принцип декомпозиции, позволяющий проводить синтез автономно для каждой задачи управления. Возможность декомпозиции определяется существенным разделением темпов управления в отдельных контурах (процессы регулирования тока и момента, процессы оценивания величины и направления магнитного поля, процессы регулирования положения и скорости могут рассматриваться независимо);
- векторное управление асинхронным двигателем с ориентацией по вектору магнитного поля (потокосцепления) двигателя;
- независимое управление электромагнитным моментом и величиной поля при постоянной величине магнитного поля;
- построение регулятора электромагнитного момента (активного тока) во вращающейся системе координат;
- использование при синтезе алгоритма управления разностных (точных) уравнений процессов в двигателе, что соответствует (адекватно) принципу цифрового процессорного управления;
- использование автономных наблюдателей (апроксиматоров) значений регулируемых переменных и регулирование по прогнозным значениям переменных, исключающее запаздывание вычислений управления, осуществляемых в процессоре последовательно во времени;
- использование автономных оптимальных фильтров (фильтров Калмана), обеспечивающих фильтрацию измеряемых значений переменных состояния без динамической ошибки;
- использование конечношаговых алгоритмов управления (дискретных скользящих режимов), обеспечивающих максимальное быстродействие регуляторов;
- использование оптимальных векторных алгоритмов широтно-импульсного регулирования выходного напряжения автономного инвертора, питающего асинхронный двигатель.
4. Цифровые алгоритмы управления моментным приводом
4.1. Состав и структура регулятора моментного привода (рис. 6):
- регулятор тока;
- формирователь заданий тока (координатный преобразователь поворота векторов;
- формирователь заданий напряжения;
- модель потокосцепления ротора;
- автономного фильтра;
- системы приема сигналов задания момента;
4.2. Функциональные алгоритмы регулятора привода;
4.2.1. Алгоритм регулирования тока.
Рис.6. Контур управления током моментного привода
Регулятор тока осуществляет формирование заданных значений вектора напряжения на систему ШИМ так, чтобы вектор тока в двигателе был равен заданному. На вход регулятора тока поступают значения заданных токов в фазах и их измеренные значения, а выходом регулятора тока являются заданные фазные напряжения.
Уравнения в статорных обмотках имеют вид
Ж / ( = Яедг + к1пФ + к2и + /, (21)
где ¡,и,Ф - векторы тока, напряжения статора и потокосцепление ротора, п - скорость вращения ротора, ,к/ - коэффициенты, /- неконтролируемые возмущения.
Темпы изменения магнитного поля, скорости вращения ротора и возмущений, обусловленных временными задержками в передаче команд управления силовыми ключами, токовыми паузами, ошибками вычислений оценок напряжения, достаточно медленны по сравнению с темпами регулирования тока. С учетом сказанного уравнение ошибки контура тока принимает вид
йМ / (М = / + к2 Аи, (22)
й/ / Л = 0. (23)
Величина А/ измеряется, Аи - известное управление.
При переходе к разностным уравнениям получим А/к+1 = /Т + к2Т (Хк +Аик )/2, Хк+1 = Аик, (24)
/к+1 = /к .
Величина / может быть оценена по результатам измерения ошибки регулирования на двух последовательных шагах регулирования. Если имеются измерения А/к и А/к+1 и в интервале между этими измерениями управление имело значение Аи к, то
Аик = -2//к + ик+1 /2 - 1к /к2Т,
/ = - к 2( Хк + Аик) / 2 + (А/к+1 -А/к)/ Т. (25)
Поскольку оценивание и регулирование происходят независимо, значение оценки возмущения (25) может быть отфильтровано при помощи наблюдателя фильтра Калмана в соответствии с выражением
/к = /к-1 +1 (/к - /к-1), 0 < I < 1.
Все вычисления возмущения и напряжения должны быть выполнены в течение текущего интервала ШИМ.
4.2.2. Алгоритм вычисления поля статора
Уравнение магнитного поля ротора в системе координат, ориентированной по направлению вектора поля имеет вид:
(у /й = -(Яг /Ьг )(у - ) (26)
где у - величина поля; - ток намагничивания.
Направление вектора поля - угол поворота системы координат задается уравнением
= n + (27)
Lr \|/
где iq - активный ток; п - скорость ротора.
Разностные уравнения магнитного поля ротора
(28)
Lr \|/
4.2.3. Формирователь заданного тока
Формирование заданного тока заключается в координатном преобразовании токов i(j4q в неподвижную систему координат, и нахождение
проекций вектора тока на фазные орты. Преобразование координат осуществляется в «медленном» времени, в течение основного цикла расчета и задается формулами
lock = 4k cos вк ~ iqk sin вк>
ipzk = Чк sin вк - iqk cos вк. (29)
Проекции вектора тока на фазные орты задаются соотношениями
lAzk=iazfc> *Bzk = (-1/2>аzk + (л/зТ2^)/рzk ,
iCzk =(-1 ИУсвк -(V3T2)/M.
4.2.4. Формирователь заданного напряжения
Формирование заданного напряжения осуществляется в «медленном» времени по статическим соотношениям между напряжениями в цепи статора асинхронного двигателя. Как известно, область физически реализуемых фазных напряжений АД, питаемого от АИН - правильный шестиугольник, ориентированный диаметром по направляющим ортам фаз (рис. 7). Диаметр шестиугольника определяется величиной постоянного напряжения на входе инвертора. Для исключения пульсаций момента воспользуемся «гладким» ограничением, т.е. ограничением модуля напряжения по величине. Присущее этому ограничению не полное использование имеющихся ресурсов напряжения с лихвой компенсируется использованием оптимальных симплексных алгоритмов широтно-импульсной модуляции, которые позволяют увеличить использование ресурсов напряжения примерно на 13 % по сравнению с широко используемой автономной по фазной ШИМ.
С учетом разнотемповости магнитных и электрических процессов в асинхронном двигателе и постоянства в пределах двух циклов управления значения входного напряжение АИН U0, предлагается иерархическая про-
цедура проверки реализуемости рассчитанных активного напряжения и
пк\
к+1
и напряжения намагничивания иа
ограничение модуля напряжения при векторной ШИМ
ограничение фазных напряжений
граничение модуля напряжения при пофазной ШИМ
Рис. 7. Ограничения выходного напряжения инвертора
1. Вначале осуществляется проверка реализуемости напряжения
намагничивания и.
м 1
т.е.
м 1
и.
ио. Если это условие не выполняется,
то выход ограничителя - и=
к+ь
к+1
= 0: если оно вы-
полняется, то переходим к проверке реализуемости активного напряжения 2. Для проверки реализуемости активного напряжения вводится
вспомогательная величина иЛ
,к+1
МАХ
~ У > характеризую-
ли
_£+1
противном случае иа -иа
щая максимально реализуемое в данном случае активное напряжения. Ес-
£+1 £+1 _к+1 ]г+] _£+1 £+1
<ид МАХ, то выход ограничителя равен иа = иа , ¿/^ = , в
^ — ) ид_МАХ-
4.2.5. Блок векторной широтно-импульсной модуляции
Блок обеспечивает формирование закона модуляции и закона переключения силовых ключей автономного инвертора напряжения, обеспечивающих формирование в фазных обмотках асинхронного двигателя напряжений в соответствии с заданиями на напряжения, сформированными в блоке ограничения напряжения. Закон переключения силовых ключей обеспечивает минимизацию коммутационных потерь в автономном инверторе напряжения. Аналитические выражения и алгоритмы ШИМ подробно рассмотрены в [4].
Таким образом, выполненный анализ особенностей построения цифровых систем управления асинхронного электропривода с векторным принципом управления момента позволяет обосновывать математический аппарат для разработки цифровых алгоритмов привода. С учетом особенностей функционирования цифровых систем сформулированы общие принципы разработки цифровых алгоритмов, на основании которых разработаны аналитические зависимости и процедура синтеза цифровых алгоритмов управления асинхронным моментным приводом.
Список литературы
1. Микропроцессорные автоматические системы регулирования. Основы теории и элементы: учеб. пособие / В.В. Солодовников [и др.]. М.: Высшая школа, 1991. 255 с.
2. Digital vector control algorithm for induction motor torque drive A. Shipunov [at al.] // Proceedings of the 28th Annual Conference of IEEE Industrial Electronics Sosiety, IECON'02.- Sevilla, Spain, 2002. Р. 996-1001.
3. Александров Е.В. Асинхронный моментный электропривод с векторным принципом управления // Вестник ТАИИ. Тула, Изд-во ТАИИ, 2010. Вып. 1. С.334-342.
4 Рывкин С.Е. Широтно-импульсная модуляция напряжения автономных инверторов. Электричество, 1997. №6. С.
Aleksandrov E. V.
FEATURES OF SYNTHESIS OF DIGITAL AC ELECTRIC DRIVES BASED ON AN INDUCTION MOTOR
Basic provisions for organization of digital control and general principles of control algorithms development are formulated, mathematical foundations of the discrete systems analysis are considered, an example of synthesis of induction motor digital control is given.
Key words: digital control, algorithm, discrete system, analysis, synthesis, difference equations, induction motor, electric torque drive.
Получено 19.06.12
