CC BY
669
УДК 621.313.3
ОПРЕДЕЛЕНИЕ УГЛОВОГО ПОЛОЖЕНИЯ РОТОРА СИНХРОННОГО ДВИГАТЕЛЯ С ВОЗБУЖДЕНИЕМ ОТ ПОСТОЯННЫХ МАГНИТОВ*
Н.В. Клиначев, Н.Ю. Кулёва, С.Г. Воронин
Описан наблюдатель электрического угла поворота вала синхронного двигателя с постоянными магнитами, состоящий из вычислителя проекций вектора потокосцепления на оси неподвижной системы координат а-в и следящего контура, сопровождающего этот вектор по фазе (контур фазовой автоподстройки частоты, далее ФАПЧ). Показано, что в режимах движения электропривода с постоянным ускорением ошибка следящего контура имеет постоянную составляющую. Для её устранения предложено использовать метод комбинированного управления. В качестве опорного сигнала на интегратор угла следящего контура (контур ФАПЧ) подается выходная координата наблюдателя скорости вала эквивалентного двигателя постоянного тока. ПИ-регулятор компенсирует малые отклонения в контуре. Модификация наблюдателя не препятствует функционированию следящего контура в скользящем режиме. Сторонний наблюдатель скорости вала позволяет замкнуть контур скорости привода с сохранением предельной для двигателя динамики.
Ключевые слова: электропривод, бездатчиковый, векторный, релейно-векторный, синхронный двигатель с постоянными магнитами, наблюдатель, потокосцепление, ФАПЧ, моделирование в Jigrein4WEB.
Введение
Для обеспечения работы синхронного двигателя с возбуждением от постоянных магнитов (СДПМ) в режиме вентильного двигателя (ВД) требуется информация об угловом положении ротора относительно статора. При этом во многих случаях использование для этой цели отдельного датчика углового положения в виде вращающегося трансформатора, редуктосина или энкодера нежелательно или недопустимо. Отсюда возникает задача косвенного определения положения ротора относительно статора по наблюдаемым координатам, мгновенным значениям фазных напряжений и токов, и формирования на этой основе одной из перечисленных координат. Имеются технические решения, позволяющие решить такую задачу [1]. К ним относится, в частности, метод определения положения ротора по наблюдаемым значениям потокосцепления статорной обмотки в неподвижной системе координат по ортогональным осям а и в [2]. При всей очевидности и кажущейся простоте этого метода при его практической реализации возникает ряд сложностей. Во-первых, возникает проблема численного интегрирования наблюдаемых производных от потокосцепления [3], во-вторых, возможны затруднения при использовании полученной информации для формирования фазных напряжений или токов ВД и, что особенно важно параметры двигателя - индуктивность, активное сопротивление обмотки и намагниченность ротора могут изменяться. В настоящей статье сделана попытка решения перечисленных проблем с
подтверждением эффективности полученных технических решений путём численного моделирования электропривода.
Наблюдатель электрического угла поворота вала СДПМ для цифровой системы управления
Обратимся к известному уравнению напряжений статорной обмотки электродвигателя переменного тока, записанному в векторной форме:
и = /Д + <«у^ , (1)
где и, I, Фх - мгновенные значения векторов напряжения, тока и потокосцепления статора, R - активное сопротивление статорной обмотки.
Потокосцепление статорной обмотки можно записать следующим уравнением:
ф = • <П/<И + Ф, (2)
где - индуктивное сопротивление статорной обмотки, а Ф - потокосцепление статорной обмотки с магнитным полем ротора.
Тогда уравнение напряжения обмотки статора можно представить выражением
П = /Д + Ь5 • <Л/<И + Ф. (3)
Перейдём от векторной формы записи к уравнениям напряжения, записанным в неподвижных статорных координатах аир, совместив ось а с магнитной осью первой фазы. В результате получим систему уравнений:
иа = 1аЯ8 + Ь8 ■ ¿Иа/СИ + йфа/йЬ,
Щ = ¿рДх + ■ й1р/(И + ^фр/^£, (4)
где иа, щ, ¿и, ¿р, фи и фр - скалярные значения напряжений, токов и потокосцепления обмотки ротора по осям а и в, которые могут быть получены путём
* Работа выполнена при поддержке ФГБОУ «Российский фонд фундаментальных исследований» в рамках договора № НК 14-07-00514\14 от 10.02.2014.
известных координатных преобразований по измеренным мгновенным значениям фазных напряжений и токов, Rs - активное сопротивление статора.
Разрешим уравнения (4) относительно пото-косцепления обмотки ротора, получим мгновенные значения его проекций на оси а и в:
Фа _ / (^а — ia.Rs)d.t — laLs,
Фр = / Op - ipRs)dt - ipLs. (5)
Примем форму ЭДС двигателя синусоидальной, тогда мгновенные значения проекций пото-косцепления содержат информацию о действительном угловом положении оси намагниченности ротора фг относительно статорной оси а и могут быть записаны выражениями:
Фа = Фт COS фг, фр = Фт Sin фг, (6)
где Фт - модуль потокосцепления обмотки ротора с магнитным полем статора, фг - действительное угловое положение оси намагниченности ротора.
Учитывая, что Фт есть величина постоянная и определяется обмоточными данными статорной обмотки и намагниченностью ротора, оценку углового положения последнего определим из выражения
ф = ЭТС^р^Д (7)
где ф - оценка углового положения оси намагниченности ротора.
Для исключения длительного переходного процесса при включении бездатчикового электропривода, интеграторы вычислителя мгновенных значений проекций вектора потокосцепления (5) охватывают неглубокой отрицательной обратной связью [3] (замещают их апериодическими звеньями первого порядка с большой постоянной времени):
фа = (^а — ^a^s)/^ — ^a^s *
* (иа - iaRs) ■ Тф/(1 + V) - iaLs, (8)
Фр = (мр - i$Rs)/s - ipLs *
* (ир - ¿p^s) ■ 7’ф/(1 + 7’ф5) - ipl^
где 7ф - постоянная времени апериодического звена первого порядка, s - оператор Лапласа.
Логарифмическая фазо-частотная характеристика (ЛФЧХ) апериодического звена меняется от 0° до требуемых -90° на протяжении двух декад. Поэтому при низких частотах вращения наблюдатель электрического угла поворота вала выдает оценку с соответствующим упреждением. Для первой итерации рекомендуется выбрать постоянную времени в соответствии с формулой
где шн - максимальная частота тока статора, гр - число пар полюсов, D - нижний предел диапазона регулирования скорости вала, который определен стремящимися к бесконечности в начале шкалы погрешностями измерения тока и регулирования питающего статор напряжения и равен (30...20), 10 - декада от частоты сопряжения участков.
Наблюдатель электрического угла поворота вала СДПМ для аналоговой системы управления
В случае отсутствия возможности вычислить функцию арктангенса [4] угол идентифицируют с помощью следящего контура фазовой автоподстройки частоты (ФАПЧ) - специализированной системы автоматического регулирования, которая способна вести сопровождение входного сигнала (вектора потокосцепления) по фазе, подстраивая до совпадения фазу собственного управляемого генератора. Следящий контур ФАПЧ состоит из четырехквадрантного фазового детектора (преобразователя Парка), ПИ-регулятора и генератора управляемого напряжением. Мгновенные значения проекций вектора потокосцепления фи и фр являются аргументами для преобразования Парка:
Фт = Re(vP) = фи cos ф + фр sin9, (10)
где Фт - модуль потокосцепления.
Для обеспечения максимального момента необходимо, чтобы векторы статора и ротора были ортогональны (у = 0). Для оценки отклонения от ортогональности векторов введём переменную х, которая оценивает величину проекции вектора статора и определяется по выражению
0 * х = Im(vP) = фр cos ф — фи sin ф. (11)
Модуль потокосцепления будет вычисляться безошибочно, если проекция вектора на мнимую ось вращающейся системы координат (d-q) будет равна нулю (что выполниться при совпадении действительного значения угла фг и его оценки ф). Выполнение условия (х =0) отслеживает ПИ-регулятор, формирующий управляющий сигнал частоты w для интегратора угла по уравнению
v = x-Kp(1 + Ts-s)/s, (12)
Ts = ^, Kv> -.
5 Rs> V Ts
Тогда оценка углового положения оси намагниченности ротора определяется следующим уравнением
ф = w ■ (1/s). (13)
Наблюдатель скорости и электрического угла поворота вала СДПМ с прецизионным следящим контуром ФАПЧ
Очевидно, что при линейном увеличении или снижении скорости, ошибка определения углового положения будет ненулевой. Предложим решение, позволяющее исправить это. Дополним систему управления наблюдателем скорости вала эквивалентного двигателя постоянного тока (ДПТ), который получим из уравнения напряжения по оси q вращающейся системы координат двигателя, при условии у = 0 в виде
uq*i^Rs + L% + e, (14)
где uq = UpCosф — и^тф,
iq = ip cos ф — ia sin ф, e = Фт • О).
Откуда
ш = ——----------,
где Т = Ls/Rs - постоянная времени двигателя, 5 - оператор Лапласа.
Выходной сигнал наблюдателя скорости необходимо подать на апериодическое звено с постоянной времени равной электромагнитной постоянной двигателя Т8, для получения устойчивой к шумам разностной вычислительной схемы. Оценка скорости исказится пренебрежимо мало, а операция дифференцирования измеряемого тока статора будет исключена:
ш[п] - (ич/(1 + 7^-5)- 1чИ5)/(Ке/гр), (16)
где ш[п] - вычисленная угловая скорость с задержкой на период дискретизации.
Вычисленную угловую скорость с задержкой на период дискретизации контура тока (16) в соответствии с методом комбинированного управления, введём в контур идентификации электрического угла поворота вала, т. е. просуммируем с выходным сигналом ПИ-регулятора, которая определяется по уравнению:
Лш = х^Хр(1 + Г^5)/5, (17)
где Дш - сигнал на выходе регулятора, Кр - коэффициент пропорционального канала ПИ-регулятора.
Тогда, оценка углового положения оси намагниченности ротора определяется следующим уравнением:
ф[п] « (ш[п — 1] + Дш) • (1/5). (18)
Очевидно, что чем точнее будет идентифицирована скорость (15, 16), тем меньше будет сигнал на выходе регулятора (17). Данное решение позволяет устранить постоянную составляющую ошибки при линейном увеличении и уменьшении скорости. Для уменьшения величины ошибки слежения контур ФАПЧ может быть переведен в скользящий режим [5, 6]. Поэтому для замыкания контура скорости бездатчикового электропривода рекомендуется использовать выходную координату того же наблюдателя (16).
Восстановление и кэширование аргументов наблюдателя электрического угла поворота вала для релейновекторного электропривода
Релейно-векторный электропривод (с прямым управлением моментом) отличается способом формирования питающего статор напряжения [2]. В соответствующей системе управления нет координат ua, Up, Uq. Их можно вычислить, зная напряжения на стойках моста иА, ив и ис: ип = (иА + ив+ ис)/3, ии = иА- ип,
Uv Uq ^n, ^С ,
Ua = Uu, Up = (uv — uw)/sqrt(3) = (19)
= (uu + uv + uv)/sqrt(3),
Uq = Up cos ф — ua sin ф.
Поскольку в релейно-векторной системе управления вектор питающего статор напряжения принимает лишь шесть ненулевых значений фиксированной амплитуды, программа для микроконтроллера может кэшировать ортогональные координаты иа, Up в относительных единицах, в виде таблицы констант и приводить их к абсолютным значениям, используя в качестве базовой величины измеряемое напряжение шины постоянного тока Udc.
Результаты численного эксперимента
Для проверки предложенных решений были построены математические модели (в программе Jigrein [7]) бездатчиковых электроприводов. Параметры двигателей приведены в таблице. На рис. 1, а представлены наблюдатель потокосцепления, на выходе которого оценка углового положения оси намагниченности ротора, ПИД-регулятор (s-PID), наблюдатель скорости эквивалентного ДТП, на выходе которого оценка скорости (Speed) и углового положения оси намагниченности ротора (Omega). Вышеупомянутые системы управления настроены так, чтобы привод отрабатывал типовую циклограмму (разгон с постоянным ускорением, стабилизация скорости (0,1... 1)шн, реверс с постоянным ускорением) с предельной для двига-
Параметры двигателей, использованных в численном эксперименте
Тип СДПМ ДВМ100.021 TY4H-105
Число пар полюсов, 1„ 13 4
Моментный коэффициент, К, Н м/А 1.50.62 1.54.00
Момент инерции ротора двигателя, J, кгм 2.0e-3 150e-3
Активное сопротивление двух секций статора, включенных последовательно, 2Я, Ом 2.5 6.4
Реактивное сопротивление двух секций статора включенных последовательно, 2L, мГн 5.0 40
Номинальный ток, 1н, А 2.5 7.3
Номинальный момент, Мн, Нм 2.0 38
Номинальная мощность, Рн, Вт 125 3k0
Отношение механической постоянной времени к электромагнитной постоянной времени, Тт/Та 2.16 3.20
Рис. 1. Блок-схема модифицированного наблюдателя угла поворота вала СДПМ, дополненная наблюдателем скорости вала эквивалентного ДПТ
теля динамикой. Критерием качества служила осциллограмма момента, искажения которой не превышали 10 % от номинальной величины.
Результаты моделирования представлены на рис. 1, б, где 1 - оценка скорости наблюдателя скорости вала эквивалентного ДТП, 2 - оценка скорости в следящем контуре идентификации электрического угла, 3 - осциллограмма ошибки (х), 4 - постоянная составляющая ошибки.
Осциллограмма, представленная на рис. 1, в, демонстрирует устранение постоянной составляющей ошибки (4) в осциллограмме ошибки (3) контура ФАПЧ после включения канала комбинированного управления (линия связи 5).
Выводы
1. Включение канала комбинированного управления (подмешивание сигнала наблюдателя
скорости эквивалентного ДПТ) устраняет постоянную составляющую ошибки определения вектора потокосцепления.
2. Представленный наблюдатель электрического угла поворота вала электродвигателя (ф) по методу комбинированного управления без использования функции arctg ф, позволяет реализовать векторное управление вентильного двигателя с меньшими по сравнению с другими известными методами объемом вычислительных ресурсов.
Литература
1. Driss Yousfi. Efficient Sensorless PMSM Drive for Electric Vehicle Traction Systems / Driss Yousfi, Elbacha Abdelhadi, Abdellah Ait Ouahman // Electric Vehicles - Modelling and Simulations / Dr. Seref Soylu (Ed.). - 2011. - http://www.intechopen.com/books/ electric-vehicles-modelling-andsimulations/ efficient-
sensorless-pmsm-drive-for-electric-vehicle-traction-systems.
2. Воронин, С.Г. Векторное управление синхронными двигателями с возбуждением от постоянных магнитов / С.Г. Воронин, Д.А. Курносов, А.С. Кульмухаметова // Электротехника. - 2013. -№ 10. - С. 50-54.
3. Виноградов, А.Б. Векторное управление электроприводами переменного тока / А.Б Виноградов. - Иваново: Иванов. гос. энергет. ун-т им. В. И. Ленина. - 2008. - 298 с.
4. Atmel AVR1636: Configurable PMSMSensor-less Field Oriented Control using the XMEGA / [APPLICATIONNOTE]Rev.: 42061A-AVR-01/2013 // Atmel web site. - www.atmel.com/images/doc42061.pdf (дата обращения: 01.02.2014).
5. AN1078, Sensorless Field Oriented Control of
a PMSM (DS01078B) / Jorge Zambada, Debraj Deb // Microchip web site. - http://ww1.microchip.com/ downloads/en/AppNotes/01078B.pdf (дата обращения: 01.02.2014).
6. Sensorless PMSM Vector Control with a Sliding Mode Observer for Compressors Using MC56F8013 / Document Number: DRM099, Rev.:2-09/2008 //Freescale Semiconductor web site. - http://cache.freescale.com/ files/microcontrollers/ doc/ref_manual/DRM099.pdf (дата обращения: 01.02.2014).
7. Клиначёв, Н.В. Цифровые управляющие системы для электропривода: фрагменты учебного курса. 2011-2013 гг. / Н.В. Клиначёв //Моделирующая программа Jigrein: Теория, программа, руководство, модели. 2006-2014 гг. - http://model. exponenta.ru/k2/Jigrein/dcs_knv.htm (дата обращения: 01.02.2014).
Клиначев Николай Васильевич, канд. техн. наук, доцент кафедры «Электромеханики и электромеханических систем», Южно-Уральский государственный университет, г. Челябинск; klinacherv_nv@mail.ru.
Кулева Надежда Юрьевна, аспирант кафедры «Электромеханики и электромеханических систем», Южно-Уральский государственный университет, г. Челябинск; nadezhda.kuleva@mail.ru.
Воронин Сергей Григорьевич, д-р техн. наук, профессор кафедры «Электромеханики и электромеханических систем», Южно-Уральский государственный университет, г. Челябинск; voroninsg@susu.ac.ru.
Поступила в редакцию 26 марта 2014 г.
Bulletin of the South Ural State University Series “Power Engineering” _____________2014, vol. 14, no. 2, pp. 49-54
ROTOR POSITION ESTIMATION FOR PERMANENT MAGNET SYNCHRONOUS MOTOR
N. V. Klinachev, South Ural State University, Chelyabinsk, Russian Federation, klinacherv_nv@mail. ru,
N. Yu. Kuleva, South Ural State University, Chelyabinsk, Russian Federation, nadezhda.kuleva@mail. ru,
S.G. Voronin, South Ural State University, Chelyabinsk, Russian Federation, voroninsg@susu. ac. ru
The paper describes rotor position observer for permanent magnet synchronous motor, consisting of calculator of flux vector projections on the axis of the fixed coordinate system a-P and loop tracking this vector by phase (phase-locked loop, referred to as PLL). It is shown that in the electric drive motion modes with constant acceleration the loop tracking error has a DC component. To eliminate it the usage of combined control method is proposed. The reference signal for integrator circuit (in the PLL circuit for VCO) is provided by output coordinate of the equivalent DC motor speed estimator. PI controller compensates small deviations in the tracking loop. Modified observer does not disturb the sliding mode of tracking loop. Outside shaft speed estimator allows to close drive speed loop while preserving the motor dynamics on its maximum.
Keywords: electric drive, sensorless, vector, relay-vector, estimator, PLL, FOC, DTC, PMSM, BLAC Motor, modeling in Jigrein4WEB.
References
1. Driss Yousfi, Abdelhadi Elbacha and Abdellah Ait Ouahman (2011). Efficient Sensorless PMSM Drive for Electric Vehicle Traction Systems, Electric Vehicles - Modelling and Simulations, Dr. Seref Soylu (Ed.), InTech, Available at: http://www.intechopen.com/books/electric-vehicles-modelling-andsimulations/ efficient-sensorless-pmsm-drive-for-electric-vehicle-traction-systems.
2. Voronin S.G., Kurnosov D.A., Kul’mukhametova A.S. Vector Control of Permanent-Magnet Synchronous Motors, Russian Electrical Engineering, 2013, vol. 84, iss. 10, pp. 581-585.
3. Vinogradov A.B. Vektornoe upravlenie jelektroprivodami peremennogo toka [Field oriented control for alternating current electrical drives]. Ivanovo, Ivanovo State Power University named after V.I. Lenin Publ., 2008. 298 p.
4. Atmel AVR1636: Configurable PMSM Sensorless Field Oriented Control using the XMEGA [APPLICATION NOTE] Rev.: 42061A-AVR-01/2013 Available at: www.atmel.com/images/doc42061.pdf. Accessed: 01.02.2014.
5. AN1078, Sensorless Field Oriented Control of a PMSM (DS01078B) Jorge Zambada, Debraj Deb. Microchip web site. Available at: http://ww1.microchip.com/ downloads/en/AppNotes/01078B.pdf. Accessed: 01.02.2014.
6. Sensorless PMSM Vector Control with a Sliding Mode Observer for Compressors Using MC56F8013 Document Number: DRM099, Rev.:2-09/2008 // Freescale Semiconductor web site. Available at: http://cache. freescale.com/files/microcontrollers/ doc/ref_manual/DRM099.pdf. Accessed: 01.02.2014.
7. Klinachev N.V. Digital control systems for electric drive: fragments of the training course, 20112013. Simulation program Jigrein: Theory, program, user’s guide, models, 2006-2014. Available at: http://model.exponenta.ru/k2/Jigrein/ dcs_knv.htm. Accessed: 01.02.2014.
Received 26 March 2014
