НАУЧНО-ТЕХНИЧЕСКИИ ВЕСТНИК ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ И ОПТИКИ ноябрь-декабрь 2023 Том 23 № 6 http://ntv.ifmo.ru/
SCIENTIFIC AND TECHNICAL JOURNAL OF INFORMATION TECHNOLOGIES, MECHANICS AND OPTICS
l/ITMO
N
ISSN 2226-1494 (print) ISSN 2500-0373 (online)
m hn h v I ?я N fi „„//*, , ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ Н ОПТИКИ
November-December 2023 Vol. 23 No 6 http://ntv.ifmo.ru/en/ 4 '
КРАТКИЕ СООБЩЕНИЯ BRIEF PAPERS
doi: 10.17586/2226-1494-2023-23-6-1242-1246 УДК 681.51.015
Оценивание параметров синхронного двигателя с постоянными магнитами
Антон Александрович Пыркин1®, Алексей Алексеевич Ведяков2, Антон Кириллович Голубев3
1>2>3 Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация
1 ругкт@^о.гаи, https://orcid.org/0000-0001-8806-4057
2 [email protected], https://orcid.org/0000-0003-4336-1220
3 [email protected], https://orcid.org/0009-0008-6951-2799
Аннотация
Рассмотрена задача оценивания параметров неявнополюсного синхронного двигателя с поверхностно установленными постоянными магнитами. Предложена параметризация нелинейной модели двигателя, позволяющая получить регрессионное соотношение с использованием измеренных (вычисленных) значений силы тока и напряжения в обмотках статора и углового положения ротора. С помощью метода динамического расширения и смешивания регрессора синтезирован алгоритм оценивания искомых параметров за конечное время.
Ключевые слова
синхронный двигатель с постоянными магнитами, динамическая параметризация модели, оценивание параметров Благодарности
Работа выполнена при поддержке Министерства науки и высшего образования Российской Федерации, соглашение № 075-11-2023-015 от 10.02.2023, «Создание высокотехнологичного серийного производства энергоэффективных синхронных электродвигателей со встроенным интеллектуальным датчиком положения и функциями самодиагностики для робототехники и цифровых систем автоматизации».
Ссылка для цитирования: Пыркин А.А., Ведяков А.А., Голубев А.К. Оценивание параметров синхронного двигателя с постоянными магнитами // Научно-технический вестник информационных технологий, механики и оптики. 2023. Т. 23, № 6. С. 1242-1246. doi: 10.17586/2226-1494-2023-23-6-1242-1246
Parameter estimation of permanent magnet synchronous motor Anton A. Pyrkin1H, Alexey A. Vedyakov2, Anton K. Golubev3
l,2,3 ITMO University, Saint Petersburg, 197101, Russian Federation
1 [email protected], https://orcid.org/0000-0001-8806-4057
2 [email protected], https://orcid.org/0000-0003-4336-1220
3 [email protected], https://orcid.org/0009-0008-6951-2799
Abstract
The problem of estimating the parameters of non-salient synchronous motor with surface-mounted permanent magnets is considered. A parameterization of a nonlinear motor model is proposed, which allows obtaining a linear regressor equation using measured (estimated) values of current and voltage in the stator windings and the angular rotor position. Using the method of dynamic regressor extension and mixing, an algorithm for estimating the desired parameters in finite time is designed. Keywords
synchronous motor with permanent magnets, dynamic parameterization of the model, parameter estimation Acknowledgements
The work was supported by the Ministry of Science and Higher Education of the Russian Federation, Agreement No. 07511-2023-015, 10.02.2023, "Creation of high-tech serial production of energy-efficient synchronous electric motors with integrated intelligent position sensor and self-diagnosis functions for robotics and digital automation systems".
© Пыркин А.А., Ведяков А.А., Голубев А.К., 2023
For citation: Pyrkin A.A., Vedyakov A.A., Golubev A.K. Parameter estimation of permanent magnet synchronous motor. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2023, vol. 23, no. 6, pp. 1242-1246 (in Russian). doi: 10.17586/2226-1494-2023-23-6-1242-1246
Широкое распространение синхронных двигателей с постоянными магнитами (СДПМ) в авиационной технике, автомобильном машиностроении, биомедицинской аппаратуре, бытовой технике, робототехниче-ских системах, а также в промышленности, особенно в системах регулирования скорости с большим диапазоном и высоким темпом пусков, остановок и реверса [1], привело к росту исследований в области управления данными двигателями [2, 3]. Контроль за значениями параметров математической модели (сопротивление, индуктивность, момент инерции, коэффициенты сил трения) очень важен в задаче проектирования современных двигателей с встроенной интеллектуальной системой управления и режимом самодиагностики. СДПМ с повышенной надежностью и сроком работы востребованы в современных безлюдных производствах в составе роботизированных комплексов, где участие человека минимально.
Представленная работа посвящена синтезу алгоритма оценивания параметров неявнополюсного СДПМ с поверхностной установкой постоянных магнитов, обеспечивающего глобальную сходимость ошибок оценивания параметров к нулю за конечное время. Схема подключения может быть любая конвенциональная, но это необходимо учесть при переходе к модели двухфазной машины. Алгоритм основан на динамической параметризации нелинейной модели с получением линейного регрессионного соотношения, на основе которого, с использованием метода динамического расширения и смешивания регрессора, формируются оценки всех параметров двигателя.
Типовая двухфазная оф-модель, описывающая динамику исследуемого СДПМ в неподвижной системе координат, связанной со статором, имеет следующий вид [4]:
(1)
0 = ю,
_ jЮ = -ю + ПрРЛ - хь
ъ = и - т, Ь = Ъ - ХтС(0),
где Ъ 6 Я2 — потокосцепление (полный магнитный поток); 1 6 Я2 — сила тока в обмотках статора; и 6 Я2 — напряжение на обмотках статора; Ь и Я — индуктивность и сопротивление обмоток статора; ] > 0 — момент инерции ротора; 0 6 [0, 2п] — угловое положение ротора; ю 6 Я — механическая угловая скорость; / — коэффициент вязкого трения; ть 6 Я — момент нагрузки; пр — количество пар полюсов; Ът > 0 — магнитный
и вектор C(0) =
0 -1 1 0
здесь и далее индекс «Т»
поток от постоянных магнитов; матрица J = 'cos(np0) ^ш(пр0)
означает операцию транспонирования.
На основе выражений для потокосцепления и его производной (1) получим соотношение
di
L— = -Ri - XmnpœJC(0) + u. dt
(2)
Доступными (измеряемыми или вычисленными) сигналами являются сила тока, напряжение обмоток статора и угловое положение ротора. Момент внешней нагрузки предполагается постоянным.
Выполним синтез алгоритма оценивания параметров щ = со1(Я, Ь, Ът) 6 Я3 и п2 = со1(£ j, тЬ) 6 Я3 за конечное время, при условии, что время > 0. Для / > /^получим:
||fii(0 - nill < ei, ||%(0 - %ll < £2,
(3)
где П1 = со1(Я, Ь, Ът) и п2 = со1(/, У, ТЬ) — оценки искомых параметров; е1 и е2 — величина, характеризующая точность оценивания, которая может быть принята равной нулю, если доступные сигналы не искажены шумами или ошибками измерений.
Вопрос оценивания достижимой точности при наличии помех измерения и ее зависимости от параметров алгоритма представляет собой отдельную задачу и в настоящей работе не рассмотрен.
На первом шаге получим регрессионную модель относительно неизвестных параметров Применим
1
к соотношению (2) линейный фильтр вида-, при
d Р + а
условии, что р: = ~ и а > 0:
Ь—[1] = -Я—[1] - Ът— [п„О^(0)] +
р + а
р + а
р + а
(4)
р + а
u.
Заметим, что npœJC(0) = —(C(0)), тогда получим
соотношение вида
dt
Vi = п1ФЬ
(5)
где ni = col(R, L, Xm) £ R3 и вычислимые функции Vi £ R2 и ф1 £ R3x2
1
Vi :
p + a
"[u], =
1
"[i]
[i]
p+a p+a p+a
[C(0)]
Выберем к > 2 различных значений параметра линейного фильтра а = ак в регрессионном уравнении (5), а также добавим масштабирующий коэффициент Д1 > 0, получим расширенную регрессионную модель
*i(0 = nWO
(6)
где ^i(t) = ^i row(Vi(t, ai), ..., Vi(t, ak)) £ R1x2k, Ф^) = ^i row(9i(t, ai), ..., 9i(t, ak)) £ R3x2k.
Умножим функции модели (6) на Ф1 (t), получим
+
*1(/)ФТ(0 = п!Ф1«Ф{ (t).
(7)
Применим процедуру смешивания [5, 6], получим набор скалярных моделей:
Zi(t) = Ai(t)ni,
(8)
где 2т(/) = (Т^ФЦ ««Ф^Ф^ (/)} 6 Я3, Д^) = = йеНФ^ОФ^ХО} 6 Я.
Для получения оценок параметров модели (4) за конечное время воспользуемся модифицированной версией алгоритма градиентного спуска [7, 8]:
ni(t) = , -,А(0 - îi(0)^i(t)),
1 - Wiii)
li(t) = YiAi(t)(Zi(t) - Ai(t)li(t)),
(9)
где у > 0 — коэффициент адаптации, а функция w i равна:
'Pi, wi(t) > pi,
.Wi(t), Wi(t) < Pi, ^i(t) = -YiA2(t)Wi(t), (10)
wi(0) = 1,
wi(t) =
с настроечным параметром р1 6 (0, 1). Алгоритм из выражений (9) и (10) обеспечивает глобальную сходимость ошибки оценивания к нулю в том случае, если
3ti > 0: I Ajdx > -1/Yiln(pi).
о
(11)
На втором шаге применим линейный фильтр вида
1 й
где р: = — и в > 0, ко второму уравнению си-
(P + ß)2 стемы (1):
dt
(p + ß)2
J P „ [0] = -f—[0] + n ■ (p + ß)2 J 1
- TL-[1(t)]
(p+ß):
;[i!Jî] -
где 1(0 — функция Хевисайда.
Перепишем второе уравнение системы (1) в виде линейной регрессии
V2 = "^Ф^
(12)
где n2 = colj, f tl) £ R3; ф2 £ R и ф2 £ R1x3 — функции
вида:
V2 = np
1
Ф2
P1
(P + ß)2 P
[i!Jî],
-[0] ---[0] -
(p + ß)2 (p + ß)2 (p + ß)
;[1(t)]
Функция ф2 может быть вычислена на основе доступного сигнала угла поворота 0, а ¥ зависит от недоступного сигнала потокосцепления к. Благодаря алгоритму состоящему из выражений (9), (10) и четвертому
уравнению системы (1) восстановим мгновенное значение потокосцепления и, пользуясь оценкой полного магнитного потока к при формировании регрессора, для регрессионного уравнения (12) с тремя различными параметрами в повторим алгоритм из выражений (7)-(10) с масштабирующим коэффициентом ц2 > 0 для получения оценок
П2« = . * (¿2(0 - 12(0)^2«), 1 - ю^г)
¿2« = 72^(0(22(0 - Д2«?2(^)),
V (Л =( P2, " P2,
1^(0, < Р2, ^2« = ^2(0^(0, ^2(0) = 1, У2 > 0,
21(Г) = (¥2(0Ф2(М{Ф2«Ф2(0}, Д2(t) = det{Ф2(t)Ф2(t)},
*2(0 = М^^С, в1), Р3)),
Ф2(0 = Ц2^(ф2(/, Р1), ..., ф2(/, Р3)),
¥2 = пр, 1 [iTJk], (р + а)2
Р „, 1
Ф2
;[0]
[0]
(P + ß)2 (p + ß)2 (p + ß)2
î = Li + îmC(0).
[i(t)]
(13)
Полученный алгоритм обеспечивает глобальную сходимость ошибки оценивания к нулю при отсутствии помех измерений и выполнение цели (3) при их наличии, если
1
3t2 > 0: I A2dT > — ln(p2), p2 £ (0, 1).
(14)
Ъ
Для выполнения условия (11) требуется, чтобы в регрессоре ф1 не было линейно зависимых компонент, и достаточно, если его компоненты не сходятся к нулю. Последнее справедливо в обычных режимах функционирования двигателя, кроме его полной остановки. Линейная зависимость может наблюдаться на практике при векторном управлении, если желаемое значение тока 1С1 (сила тока по продольной оси во вращающейся системе координат, связанной с ротором) будет равно нулю. В этой связи для оценивания параметров Я и Хт можно использовать подход инжекции высокочастотного сигнала [9], а именно задать желаемое значение 1с} в виде гармонической функции большой частоты и малой амплитуды. Для выполнения условий (14) требуются аналогичные условия для регрессора ф2, что может быть обеспечено при непостоянном угловом ускорении. В противном случае, в зависимости от длительности переходных процессов и параметров алгоритма, условия (14) могут быть нарушены.
Отметим, что в предложенном алгоритме использован только параметр количества пар полюсов Пр, который можно считать известным, а остальные параметры
t
t
0,2
Г 0,1
кЛ
о
t, с
Рисунок. Угловая скорость ротора (а) и нормы ошибки оценивания параметров для flj (b) и n2 (с) Figure. Rotor rotation speed (a) and norms of estimation errors for parameters ^ (b) and n2 (c)
требуют оценки. В этой связи влияние параметрических возмущений исключено. Если рассмотреть помехи измерения, то на практике для уменьшения их влияния можно использовать онлайн версию метода наименьших квадратов с аналогичной модификаций для получения сходимости за конечное время. Отметим, что используемый линейный фильтр может быть заменен на любой другой линейный оператор для лучшей фильтрации сигналов. Полученный алгоритм достаточно чувствителен к смещениям в помехах измеряемых сигналов, поэтому при их наличии можно воспользоваться подходом, предложенным в работе [10], с параметризацией и оцениванием этих смещений.
На рисунке представлены результаты численного моделирования работы двигателя BMP0701F Schneider Electric [2], согласно модели (1), со следующими значениями параметров: R = 8,875 Ом, L = 0,04003 Гн, np = 5, f = 0,05 Нмс, lm = 0,2086 Вб, j = 610-5 кгм2, tl = 0,1 Нм. При этом алгоритм (5)—(10), (12), (13) с параметрами у1 = у2 = 5, р = 0,9, ni(0) = 0, n2(0) = 0, а1 = 50, а2 = 100, ß1 = 50, ß2 = 75, ß3 = 100.
Для регулирования скорости двигателя использовано векторное управление со следующими параметрами для регуляторов скорости kp = 0,0528, ki = 23,6871 и тока kp = 0,1601, ki = 35,5. Требуемая скорость, подаваемая на вход регулятора, составила 15 рад/c. Необходимое значение тока iq (сила тока по поперечной оси, во вращающейся системе координат, связанной с
ротором) было сформировано регулятором скорости, а выбрано равным 0^ш(50/). Внешний момент Ть приложен с момента времени 0,1 с.
Дополнительно сигналы ^ и Ф1 были пропущены
50 '
через фильтр с передаточной функцией-. Кроме
/7 + 50
того, на алгоритм оценивания сигналы подавались после завершения переходных процессов на фильтрах, принятого равным 0,2 с, что реализовано с помощью изменения масштабирующих коэффициентов: Д1 = 0 для / 6 [0; 0,2) с и ц1 = 2000 для / > 0,2 с, д2 = 0 для / 6 [0; 0,4) с и д2 = 5000 для / > 0,4 с.
В измерения силы тока добавлена помеха с равномерным случайным распределением [-0,05; 0,05] А, а в значения напряжения, поступающие на вход алгоритма, помеха с равномерным случайным распределением [-0,1; 0,1] В. Важно, что напряжение, как правило, было не измерено, а вычислено, а соответствующая систематическая ошибка расчетов смоделирована подачей на алгоритм зашумленных значений напряжения, умноженных на 0,99. Сила тока после завершения переходного процесса была изменена в диапазоне [-0,85; 0,85] А, а напряжение — [-23,5; 23,5] В.
На рисунке видно, что предложенный алгоритм оценивания обеспечил сходимость норм векторов ошибок оценивания параметров в некоторую малую область при наличии помех измерения и систематической ошибки расчета напряжения.
Литература
1. Sakunthala S., Kiranmayi R., Mandadi P.N. A review on speed control of permanent magnet synchronous motor drive using different control techniques // Proc. of the 2018 International Conference on Power, Energy, Control and Transmission Systems (ICPECTS). 2018. P. 97102. https://doi.org/10.1109/icpects.2018.8521574
2. Bobtsov A.A., Pyrkin A.A., Ortega R., Vukosavic S.N., Stankovic A.M., Panteley E.V. A robust globally convergent position observer for the permanent magnet synchronous motor // Automatica. 2015. V. 61. P. 47-54. https ://doi. org/1 0. 10 1 6/j. automatica.2015.07.032
3. Ortega R., Monshizadeh N., Monshizadeh P., Bazylev D., Pyrkin A. Permanent magnet synchronous motors are globally asymptotically stabilizable with PI current control // Automatica. 2018. V. 98. P. 296301. https://doi.org/10.1016/j.automatica.2018.09.031
References
1. Sakunthala S., Kiranmayi R., Mandadi P.N. A review on speed control of permanent magnet synchronous motor drive using different control techniques. Proc. of the 2018 International Conference on Power, Energy, Control and Transmission Systems (ICPECTS), 2018, pp. 97102. https://doi.org/10.1109/icpects.2018.8521574
2. Bobtsov A.A., Pyrkin A.A., Ortega R., Vukosavic S.N., Stankovic A.M., Panteley E.V. A robust globally convergent position observer for the permanent magnet synchronous motor. Automatica, 2015, vol. 61, pp. 47-54. https://doi.org/10.1016/j. automatica.2015.07.032
3. Ortega R., Monshizadeh N., Monshizadeh P., Bazylev D., Pyrkin A. Permanent magnet synchronous motors are globally asymptotically stabilizable with PI current control. Automatica, 2018, vol. 98, pp. 296-301. https://doi.org/10.1016Zj.automatica.2018.09.031
4. Nam K.H. AC Motor Control and Electric Vehicle Applications. CRC Press, 2010. 435 p.
5. Aranovskiy S., Bobtsov A., Ortega R., Pyrkin A. Performance enhancement of parameter estimators via dynamic regressor extension and mixing // IEEE Transactions on Automatic Control. 2017. V. 62. N 7. P. 3546-3550. https://doi.org/10.1109/tac.2016.2614889
6. Ortega R., Aranovskiy S., Pyrkin A.A., Astolfi A., Bobtsov A.A. New results on parameter estimation via dynamic regressor extension and mixing: Continuous and discrete-time cases // IEEE Transactions on Automatic Control. 2021. V. 66. N 5. P. 2265-2272. https://doi. org/10.1109/tac.2020.3003651
7. Pyrkin A., Bobtsov A., Ortega R., Vedyakov A., Aranovskiy S. Adaptive state observers using dynamic regressor extension and mixing // Systems & Control Letters. 2019. V. 133. P. 104519. https:// doi.org/10.1016/j.sysconle.2019.104519
8. Pyrkin A., Bobtsov A., Ortega R., Vedyakov A., Cherginets D., Ovcharov A., Bazylev D., Petranevsky I. Robust nonlinear observer design for permanent magnet synchronous motors // IET Control Theory & Applications. 2021. V. 15. N 4. P. 604-616. https://doi. org/10.1049/cth2.12065
9. Jang J.H., Ha J.I., Ohto M., Ide K., Sul S.K. Analysis of permanentmagnet machine for sensorless control based on high-frequency signal injection // IEEE Transactions on Industry Applications. 2004. V. 40. N 6. P. 1595-1604. https://doi.org/10.1109/tia.2004.836222
10. Pyrkin A.A., Vedyakov A.A., Ortega R., Bobtsov A.A. A robust adaptive flux observer for a class of electromechanical systems // International Journal of Control. 2020. V. 93. N 7. P. 1619-1629. https://doi.org/10.1080/00207179.2018.1521995
4. Nam K.H. AC Motor Control and Electric Vehicle Applications. CRC Press, 2010, 435 p.
5. Aranovskiy S., Bobtsov A., Ortega R., Pyrkin A. Performance enhancement of parameter estimators via dynamic regressor extension and mixing. IEEE Transactions on Automatic Control, 2017, vol. 62, no. 7, pp. 3546-3550. https://doi.org/10.1109/tac.2016.2614889
6. Ortega R., Aranovskiy S., Pyrkin A.A., Astolfi A., Bobtsov A.A. New results on parameter estimation via dynamic regressor extension and mixing: Continuous and discrete-time cases. IEEE Transactions on Automatic Control, 2021, vol. 66, no. 5, pp. 2265-2272. https://doi. org/10.1109/tac.2020.3003651
7. Pyrkin A., Bobtsov A., Ortega R., Vedyakov A., Aranovskiy S. Adaptive state observers using dynamic regressor extension and mixing. Systems & Control Letters, 2019, vol. 133, pp. 104519. https://doi.org/10.1016Zj.sysconle.2019.104519
8. Pyrkin A., Bobtsov A., Ortega R., Vedyakov A., Cherginets D., Ovcharov A., Bazylev D., Petranevsky I. Robust nonlinear observer design for permanent magnet synchronous motors. IET Control Theory & Applications, 2021, vol. 15, no. 4, pp. 604-616. https://doi. org/10.1049/cth2.12065
9. Jang J.H., Ha J.I., Ohto M., Ide K., Sul S.K. Analysis of permanentmagnet machine for sensorless control based on high-frequency signal injection. IEEE Transactions on Industry Applications, 2004, vol. 40, no. 6, pp. 1595-1604. https://doi.org/10.1109/tia.2004.836222
10. Pyrkin A.A., Vedyakov A.A., Ortega R., Bobtsov A.A. A robust adaptive flux observer for a class of electromechanical systems. International Journal of Control, 2020, vol. 93, no. 7, pp. 1619-1629. https://doi.org/10.1080/00207179.2018.1521995
Авторы
Пыркин Антон Александрович — доктор технических наук, профессор, профессор, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, sc 26656070700, https://orcid.org/0000-0001-8806-4057, [email protected]
Ведяков Алексей Алексеевич — кандидат технических наук, доцент, доцент, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, sc 49664023200, https://orcid.org/0000-0003-4336-1220, [email protected]
Голубев Антон Кириллович — аспирант, ассистент, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, https://orcid. org/0009-0008-6951-2799, [email protected]
Authors
Anton A. Pyrkin — D.Sc., Full Professor, ITMO University, Saint Petersburg, 197101, Russian Federation, sc 26656070700, https://orcid. org/0000-0001-8806-4057, [email protected]
Alexey A. Vedyakov — PhD, Associate Professor, Associate Professor, ITMO University, Saint Petersburg, 197101, Russian Federation, sc 49664023200, https://orcid.org/0000-0003-4336-1220, vedyakov@ itmo.ru
Anton K. Golubev — PhD Student, Assistant, ITMO University, Saint Petersburg, 197101, Russian Federation, https://orcid.org/0009-0008-6951-2799, [email protected]
Статья поступила в редакцию 29.09.2023 Одобрена после рецензирования 27.10.2023 Принята к печати 22.11.2023
Received 29.09.2023
Approved after reviewing 27.10.2023
Accepted 22.11.2023
Работа доступна по лицензии Creative Commons «Attribution-NonCommercial»