Аннотация: приведена методика синтеза системы разрывного полеориентированного управления асинхронным электроприводом с активным демпфированием колебаний упругого передаточного устройства.
Ключевые слова: синтез, математическая модель, передаточное устройство, регулятор, упругие колебания, демпфирование, контур регулирования, полеориентированное управление, скользящий режим.
Технические науки
УДК 62-83: 621.313.333
Клименко Юрий Михайлович
кандидат технических наук, доцент,
доцент кафедры электротехники и
электромеханики
Днепродзержинского государственного
технического университета
Садовой Александр Валентинович
доктор технических наук, профессор,
проректор по научной работе
Днепродзержинского государственного
технического университета
Klimenko Yuri Mikhailovich
candidate of technical Sciences,
associate Professor the Department of
electrical engineering and electromechanics
Dneprodzerzhinsk state technical
University
Sadovoy Alexander Valentinovich
doctor of technical Sciences,
Professor, Vicerector on scientific work
Dneprodzerzhinsk state technical
University
СИНТЕЗ СИСТЕМЫ ПОЛЕОРИЕНТИРОВАННОГО УПРАВЛЕНИЯ АСИНХРОННЫМ ЭЛЕКТРОПРИВОДОМ С ДЕМПФИРОВАНИЕМ КОЛЕБАНИЙ УПРУГОГО ПЕРЕДАТОЧНОГО УСТРОЙСТВА
SYNTHESIS OF FIELD-ORIENTED CONTROL OF ASYNCHRONOUS ELECTRIC DRIVE WITH DAMPING VIBRATION OF AN ELASTIC TRANSMISSION DEVICES
Аннотация: приведена методика синтеза системы разрывного полеориентированного управления асинхронным электроприводом с активным демпфированием колебаний упругого передаточного устройства.
Ключевые слова: синтез, математическая модель, передаточное устройство, регулятор, упругие колебания, демпфирование, контур регулирования, полеориентированное управление, скользящий режим.
Summary: the method of synthesis of a discontinuous field-oriented control of asynchronous electric drive with active damping vibration of an elastic transmission devices
Key words: synthesis, mathematical model, transfer device controller, elastic vibrations, damping, contour regulation, field oriented control, sliding mode.
Интенсивное развитие современных технологий и оборудования выдвигает повышенные требования к электромеханическим системам (ЭМС) точного воспроизведения сложных движений (ТВСД). К ним относятся системы управления (СУ): антенными установками; радиолокационными станциями; зеркалами радиотелескопов, радиотехнических комплексов навигации, зондирования и наблюдения; специализированными устройствами военной техники, которые осуществляют поиск, наведение и автосопровождение подвижных объектов. Достижение высокой надежности, требуемого ресурса длительного функционирования без профилактического обслуживания и оперативных ремонтных работ при традиционно используемых в ЭМС ТВСД электроприводах (ЭП) на основе машин постоянного тока усложнено или невозможно. По этой причине актуальным и бурно развивающимся направлением в настоящее время для большинства механизмов высококачественного управления движением является применение бесконтактных глубокорегулируемых асинхронных электроприводов (АЭП) с векторным полеориентированным управлением (ВПУ), реализуемым быстродействующими, работающими в ключевом режиме преобразователями на силовых IGBT модулях.
Сложность задач создания высококачественных АЭП для систем ТВСД усугубляется и тем, что короткозамкнутый асинхронный двигатель (КАД) как ОУ характеризуется совокупностью звеньев, накапливающих электромагнитную и механическую энергию; имеет сложную многомерную структуру с внутренними перекрестными связями (ВПС); отличается значительной нелинейностью, вызванной флуктуациями характеристик цепи намагничивания и неидеальными свойствами передаточных устройств (ПУ); обладает нестационарностью параметров в виде изменений количественных связей между ними в статических и динамических режимах работы АЭП при нагреве электродвигателя, изменениях насыщения магнитопровода и вытеснении тока; подвержен значительным возмущениям со стороны питающих электрических цепей и ПУ; выделяется сложностью прямого измерения потокосцепления и электромагнитного момента (ЭМ).
Традиционно в системах ВПУ используются П, ПИ и ПИД регуляторы, которые обладают недостаточной робастностью, плохо противодействуют изменяющимся внешним воздействиям, упругим свойствам ПУ, люфтам и нелинейному трению.
Изложенные особенности не позволяют при синтезе СУ АЭП использовать принципы линейной теории управления по причине резкого снижения качества управления при отклонениях параметров ОУ от расчетных, нарушения автономности каналов управления потоком и электромагнитным моментом, низкого запаса устойчивости СУ вплоть до потери ее работоспособности.
Цель работы – синтез системы ВПУ КАД с высокими требованиями к качеству отработки управляющих воздействий в режимах слежения и позиционирования, обладающей низкой чувствительностью к параметрическим и координатным возмущениям, обеспечивающей активное демпфирование колебаний упругого ПУ и устойчивую работу на основе информации от наблюдателей трудноизмеряемых координат АЭП и ПУ.
Достижение поставленной цели осуществлено путем синтеза АЭП с ВПУ в классе нелинейных систем с преднамеренно организованными в контурах регулирования (КР) многомерными скользящими режимами (МСР) при замыкании последних как по непосредственно измеренной, так и по полученной методами идентификации информации о векторах состояния АЭП, необходимой для реализации синтезированных алгоритмов управления (АУ).
Для достижения поставленной цели составим математическую модель (ММ) комплекса «КАД - ПУ с упругими свойствами», которая будет использована в качестве инструмента при синтезе и исследованиях разрабатываемых АЭП.
Анализ известных и подробно описанных в научно-технической литературе математических моделей ПУ ЭМС, построенных с учетом упруго-вязких свойств механических передач, приводит к выводу о том, что в большинстве случаев при анализе и синтезе ЭМС с повышенными требованиями к качеству управления целесообразным является использование описания ее моделью эквивалентной двухмассовой упругой системы, которая несмотря на предельное упрощение, отражает физические особенности упругих ПУ с достаточной для инженерных расчетов точностью.
Математическую модель ПУ механической части ЭМС с упругими свойствами составим с учетом следующих известных допущений:
В соответствии с принятыми допущениями ПУ ЭМС с наиболее распространенной на практике кинематической схемой, изображенной на рис. 1, представим в виде обобщенной двухмассовой упругодиссипативной системы, описываемой уравнениями:
где φ1,2, ω1,2; Jд,1.,2; - угловые перемещения, угловые скорости, моменты инерции двигателя (д), первой (1) и второй (2) масс двухмассовой ЭМС:
J1 = Jд + Jд / i2 ;
i - передаточное число;
M12 и МC2- моменты упругого взаимодействия масс и нагрузки на
выходном валу;
Mf1,2 - моменты внешнего вязкого трения:
Mf1= β1 ω1 (при Mfд = 0), Mf 2 = β2 ω2 (2)
β1,2 - коэффициенты внешнего вязкого трения валов 1 и 2;
С12 - коэффициент жесткости;
b12 - коэффициент внутреннего трения в деформируемой передаче.
Рис.1 Кинематическая схема механической передачи ЭМС с упругостью первого рода (составлена автором)
В (1) и (2) параметры без индексов (•)' получены приведением к валу электродвигателя величин, обозначенных на кинематической схеме J'2; ω'1,2; φ'1,2; C'12; b'1; M'12; M'c2; M'f, 1,2. Приведение осуществлено методом эквивалентных преобразований [1], обеспечивающим выполнение закона сохранения энергии.
С учетом выражений (2) и базовых величин ωб и Мб введем в (1) следующие обозначения:
TМ1 = J1 ωб Mб-1, TМ2 = J2 ωб Mб-1, TС = Mб (C12 ωб)-1,
kf 1 = β1 ωб Mб-1, kf 2 = β2 ωб Mб-1, kс = b12 ωб Mб-1,
где ТМ1 ,ТМ2 - механические постоянные времени сосредоточенных
масс;
ТС - механическая постоянная жесткости;
k f1, k f2, kc - безразмерные коэффициенты внешнего и внутреннего
вязкого трения.
Представив третье - пятое уравнения системы (1) с учетом обозначений (3), получим:
dω°1 /dt = (M° - M°12 - kf 1 ω°1 )/ ТM1
dω°2 /dt = (M°12 - M°c2 - kf 2 ω° 2)/ ТM2 (4)
M°12 = kс (ω°1 - ω°2 ) + (φ°1 - φ°2 )/ ТС ,
где: (•)° = (•)/(•)б - параметры в относительных единицах.
В дальнейшем индексы ° при (•) упущены.
На рис. 2 представлена структурная схема нормированной системы уравнений (4), применяемая многими авторами при синтезе СУ и приводящая к алгоритмам, построенным на использовании операций дифференцирования. Для получения алгоритмов управления (АУ) с простой практической реализацией упростим модель (4), представив ее в пространстве состояний. С этой целью продифференцируем уравнение для вычисления M12 из (4) при постоянстве параметров Тс, kс и запишем результат с учетом уравнений (1)
M12 /dt = С12( ω1 - ω2) + b12 ( dω1 / dt - dω2 / dt) (5)
Рис. 2. Структурная схема упрощенной модели двухмассовой ЭМС (составлена автором)
Подставим первое и второе уравнения системы (4) в выражение (5) и запишем полученный результат совместно с первым и вторым уравнениями системы (4):
где М – электромагнитный момент КАД (7)
где С1=m Zp /2; m, Zp - числа фаз и пар полюсов КАД.
Полученная нормированная модель ПУ механической части ЭМС с упругостью первого рода описывает динамические процессы в упругой передаче с учетом диссипативных сил. Уровень информативности полученной модели позволяет использовать ее при анализе и синтезе АУ с активным демпфированием колебаний упругих ПУ, наблюдателя координат (НК) КАД и упругого ПУ при их совместной работе.
Для синтеза алгоритмов ВПУ, обеспечивающих активное демпфирование колебаний в АЭП с упругими ПУ, в качестве расчетной используем математическую модель (6), дополнив ее дифференциальным уравнением для вычисления электромагнитного момента. Указанное уравнение получим путем дифференцирования выражения (7) при ΨοΑ=const, C1=const и выполнении условий ориентации координатного базиса по вектору ΨοΑ: Ψou = ΨοΑ, Ψv=0, dΨv/dt =0.
d M/dt ≡ Сm d isv/dt , (8)
где Сm = C1∙ΨοΑ при ΨοΑ=const - постоянный коэффициент.
Используя третье уравнение системы (20) [2, с.521]
(9)
соотношение (8) и равенство М = Сm isv получим:
dМ /dt = - C4М + СmUsy / L's +Сm (C2 L's-1 ω1 - ωΨο ) isu - Сmω1 L's-1 Ψ , (10)
где C2 = (1-a)Lm ; C4 = 1-C5/ L's ; C5 = аLm-LL
.
Записав совместно (6) и (10), получим расчетную ММ в развернутой векторно-матричной форме записи:
где а11= - kf2 / Tм2 ; a12 = 1/ Tм2 ; a21 = -1/ Tc + kc kf2 / Tм2 ; a22 = -kc (Tм1+ Tм2) / Tм1 Tм2; a24 = -1/ Tc + kc kf1 / Tм1 ; a33 = -СmС4; а42 = - 1/ Tм1 ; а34 =Сm L's-1 (С2isu - Ψ); а43 = 1 / Tм1 ; а44= - kf 1 / Tм1; b2 = Сm / L's ; μ1 = -1 / Tм2 ; μ2 = kc / Tм2 ; μ3 = -СmωΨο isu.
Математическая модель (11) и соответствующая ей структурная схема (рис.3) отражают многомерность, нелинейный характер и динамическую взаимосвязь электромагнитных, электромеханических параметров, управляющих и возмущающих воздействий в ЭМС рассматриваемой структуры. Однако структурная и параметрическая избыточность такой математической модели затрудняет использование ее при синтезе СУ с минимальным количеством контурных регуляторов и сложностью информационно - датчиковой системы.
Рис.3. Структурная схема исходной математической модели (11) АЭП с ВПУ и упругим ПУ (составлена автором).
Упрощение модели и понижение ее порядка осуществим с использованием модального метода редуцирования [3, с.65÷70], предусматривающего выделение доминирующих координат ХI , непосредственное управление которыми определяет требуемые динамические характеристики, и не доминирующих ХII , которые можно свести к возмущающим воздействиям, компенсируемым контурными регуляторами координат Z I.
Вектор состояния системы Х «неособым» преобразованием Z = ΤX, det T≠0 преобразуем в новый координатный базис векторов Ζ модальных координат [4]:
d Z(t) /dt = TАT-1 Z(t) + TB U(t)+Tμ , Y = CT-1Z(t), (12)
Обозначив TAT-1 = Λ; T = V; Tμ = m; CT-1=C , представим (12) в "новом" координатном базисе:
d Z(t)/dt= Λ X(t) + VBU(t)+ m, Y = C Z(t). (13)
C учетом разбивки на группы координат ZI и Z II представим (13) в развернутой форме:
Редуцирование порядка при этом сведем к исключению дифференциального уравнения вычисления ZII . Выбор преобразующих матриц в (14) осуществлен таким образом, чтобы кроме устранения взаимосвязи координат ZI и ZII, получить матрицу Λ Ι коэффициентов Z Ι с верхнетреугольной формой:
где λ11=a21 ; λ12=a12a21 (a11)-1; λ22= a33a42(a43)-1; λ23= a242a33(a43)-1(a24(a44)-1- a23 ( a43 )-1) (a22 - a12 a21 ( a11 )-1 - a23 a42 ( a43 )-1)-1 ; λ33= a42 (a34 ( a44 )-1-a33 ( a43 )-1 - λ23
Выбор формы представления редуцированной модели в виде (15) позволяет синтезировать алгоритмы ВПУ и технически просто реализуемые структуры систем управления.
Используя концепцию Ляпунова о возмущенном - невозмущенном движении [5], представим (15) в фазовом пространстве координат возмущенного движения, переход к которому позволяет исключить из рассмотрения возмущающие воздействия μ:
где ηi = Z i – Z i* , i=1,2,3 - координаты возмущенного движения;
Z i и Z i* - значения истинных и программно-заданных координат;
U∆sv=Usv - U*sv - дополнительное управление, однозначно определяемое компонентами векторов ηi .
Задача синтеза сводится к выбору ограниченных по модулю управляющих воздействий, обеспечивающих минимум интегральных критериев качества, задаваемых для каждого из контурных регуляторов в виде функционалов [ 5] :
где V01=V03V13/V33, V02=V03V23 /V33, V03 = - λ11 λ22 λ33 - весовые коэффициенты.
Минимизация функционалов (17) гарантирует апериодический характер переходных процессов с минимально возможной для данных параметров и ограничений управляющих воздействий постоянной времени. Коэффициенты функционалов (17) определим в соответствии с функцией А.М. Ляпунова, выбранной в виде:
где J=1,2,3 – индексы переменных ω2, М12 и М, управляемых регуляторами Р1, Р2 и Р3.
Оптимальные управляющие воздействия на выходах регуляторов сформируем в виде:
где Upj – амплитуда напряжения релейного элемента;
U(pj) - напряжение на выходе регулятора координаты j.
Определим коэффициенты функций А.М.Ляпунова, входящие в АУ (19):
для Р1 V13(1) = λ22 λ33, V23(1) = - λ12 λ33 , V33(1) = λ12 λ23;
для Р2 V13(2) = 0 , V23(2) = - λ11 λ33 , V33(2) = λ11λ23 (20)
для Р3 V13(3) = 0 , V23(3) = 0 , V33(3) = -λ11λ22.
В соответствии с (20) алгоритмы (19) представим в виде:
или после перехода в исходный координатный базис фазового пространства отклонений систему (21) запишем в виде
На рис. 4 представлена функциональная схема системы ВПУ [7], в которой реализованы АУ (22). В состав схемы входят: контур управления потоком; контуры управления с релейными регуляторами Р1, Р2, Р3, осуществляющими управление координатами механического движения и активное демпфирование упругих колебаний УПУ; наблюдатель координат
Рис.4. Функциональная схема системы ВПУ КАД с активным
демпфированием упругих колебаний УПУ по АУ (22)
асинхронного двигателя и упругого передаточного устройства НК АД УПУ; типовые для систем ВПУ преобразователи координат ПК1-ПК3; контур формирования фазных токов КФФТ с преобразователем на силовых транзисторных модулях ПСТМ, датчиками ДФТ и релейными регуляторами фазных токов РФТ; вычислитель «эквивалентного» управления ВЭУ, обеспечивающий МСР в СУ; короткозамкнутый асинхронный двигатель М с тахогенератором ДС на валу, упругим передаточным устройством УПУ и датчиком положения ДП2 на выходном валу.
Наблюдатели координат НК АД и УПУ синтезированы и исследованы авторами в [ 6, с 358÷360].
Организация МСР в системе управления осуществлена на основе применения метода «эквивалентного» управления [7, с.82÷89], позволяющего разделить разнотемповые движения во внешних и внутренних контурах управления путем включения между ними выделителя ВЭУ «эквивалентного» управления, выполненного в виде двухканальной, замкнутой по сигналам на ее выходах, модели контура тока, содержащей релейные регуляторы активного и реактивного тока, модели вычисления тока с параметрами управляемого КАД.
Результаты моделирования синтезированных АУ (22) при ОКБ по |Ψr| для режима слежения представлены на рис. 5. Графики изменений задающего воздействия , координат ω1(t), ω2(t), Μс(t), Му(t), М(t), ε(t)=
при |Ψr| =const, |Ψr|(t), isu(t), Ψr α,β(t), is α,β(t), приведены для режимов: возбуждения машины до уровня Ψr ном. (0÷t1); работы с (0÷t4) при Μс = 0 (0÷t2), Μс = - Мном (t2÷t3), Μс = Мном (t3÷t4); отработки задающего воздействия
при Μс=Мном (t4÷t5), при Μс= -Мном (t5÷t6), Μс=0 (t6÷t).
Для наглядности изображений на рис.5 сигналы ω2,(ω1) и ω*2 показаны в разных масштабах.
С целью оценки чувствительности скоростной подсистемы с АУ (22) к неточности задания при расчете или изменениям в процессе эксплуатации коэффициента жесткости Сж и момента инерции J2 проведены исследования ошибок регулирования ε(t)= для режимов: работы с
(0÷t3) при Μс=0 (0÷t1), Μс= -Мном (t1÷t2), Μс=Мном (t2÷t3); отработки задающего воздействия
при Μс=Мном (t3÷t4), при Μс= -Мном (t4÷t5), Μс=0 (t5÷t).
На рис.6 представлены результаты исследования ошибок ε(t), полуполученные для случаев ε0(t) - при расчетных значениях ; ε1(t) - при Сж= 0.5
и
; ε2(t) - при Сж= 1.5
и
; ε3(t) - при J2=0.5
и
; ε4(t) - при J2=1.5
и
. Величины ошибок не превышают значений ±5%.
Рис.5. Результаты моделирования системы ВПУ АЭП с АУ (22) (получены авторами)
Выводы. Синтезированная система полеориентированного управления АД с активным демпфированием колебаний упругого передаточного устройства по алгоритмам (22) с регуляторами скорости, электромагнитногомомента АД и упругого момента ПУ при работе с МСР, организованны-
Рис.6 Результаты моделирования чувствительности ЭМС управления скоростью с АУ ( 22) к изменениям Сж и J2 от расчетных (получены авторами).
ми на основе информации от наблюдателей координат АД и УПУ, реализованных в виде работающих в реальном масштабе времени прямых замкнутых динамических моделей КАД и УПУ, снабженных работающими в скользящем режиме контурами слежения, обеспечивает высококачественное управление АЭП, активное демпфирование упругих колебаний УПУ при низкой чувствительности к параметрическим и координатным возмущениям. Разработанная СУ АЭП удовлетворяет высоким требованиям, предъявляемым к ЭМС ТВСД.
Список литературы