CC BY
46
УДК 62-83:621.313.2:681.513.68
ОПТИМИЗАЦИЯ НЕСТАЦИОНАРНЫХ ЭЛЕКТРОМЕХАНИЧЕСКИХ СИСТЕМ С ДИСКРЕТНЫМ УПРАВЛЕНИЕМ
В.Г. Букреев, И.Ю. Краснов, А.К. Чащин, С.К. Соснин
Томский политехнический университет E-mail: [email protected]
Рассматривается процедура синтеза структуры системы и законов управления на основе функции Гамильтона, изменение которой на оптимальной траектории движения позволяет организовать адаптивные алгоритмы. Приводится пример адаптивного управления электроприводом постоянного тока при существенном изменении динамического момента нагрузки.
Введение
Для обеспечения требуемых показателей качества технологического процесса в условиях нестационарности параметров внешней среды и неконтролируемости переменных состояния многосвязной электромеханической системы (ЭМС) необходимо организовать адаптивное управление протекающими процессами в системе [1, 2]. На этапе синтеза таких ЭМС принимают определенные допущения, которые позволяют на основе методов аналитического конструирования оптимальных регуляторов, синтезировать структуру и параметры адаптивных регуляторов. Количественные значения допущений обычно определяются из идеализированных условий протекания технологических процессов и функционирования электромеханической системы. Таким образом, уже на стадии проектирования ЭМС предполагается отличие назначенного функционала качества от действительного значения. Дальнейшее снижение качества в ЭМС с жесткой и неизменяемой структурой регулятора будет наблюдаться при увеличении числа неконтролируемых переменных, неизмеряемых возмущений или превышении определенных граничных значений измеряемых возмущений.
Одним из вариантов снижения эффекта влияния нестационарности возмущений на качественные показатели процесса управления ЭМС является применение методов алгоритмического конструирования нестационарных систем. Принципиальное отличие методов алгоритмического конструирования нестационарных систем от методов аналитического конструирования оптимальных регуляторов заключается в том, что первые позволяют целенаправленно создать некоторое множество структур адаптивных регуляторов, включение которых в структуру ЭМС определяется в процессе функционирования системы условиями оптимизации функционала качества. Для синтеза структуры системы и законов управления нестационарными объектами можно использовать функцию Гамильтона, изменение которой на оптимальной траектории движения ЭМС позволяет организовать адаптивные алгоритмы [3].
Постановка задачи
Динамика электромеханического объекта (ЭМО) с учетом стационарности его параметров на интервалах дискретности (/О+уТ^+уТ+О и ОО+уТ+^+Гтх)
управляющих сигналов (у - целое действительное число) может быть представлена дифференциальными уравнениями [4]:
x(t) = A1 x(t) + bU (t) + mH
при t £ (t0 + jT, to + jT + tj),
x (t) = A2 x(t) + bU (t) + mH
пРи t e (to + jT +11,to + ТЩД
(1)
где х=х(/) - вектор состояния ЭМО в момент времени /; А1, А2 - матрицы параметров объекта на интервалах подключения (от /0 до /1) и отключения (от /1 до Т^) управляющих сигналов объекта соответственно; =, Т, 2Т,... - моменты времени; /0 - время, характеризующее начальное состояние электромеханического объекта х(/=/0)=х0.
В результате интегрирования ур. (1) на интервалах времени от /0 до /1 и от /1 до Т^, где ТПх - максимальный период дискретности управляющего сигнала, записывается дискретное векторно-матрич-ное уравнение, характеризующее движение ЭМС:
х,+1 = рх, + (4 X + 6*)7 +
где ^ехр^Тпах); ^хр^^М-А*]; 0*= =ехр(А2 Тпах) ЪиЩ+А Тпах)( ¿0 Ц0+тн); 0*=Т^Ь Ц0+тн); Ьи, Ь0, тн - векторы соответствующей размерности, характеризующие параметры силовых цепей импульсного преобразователя и звеньев механической нагрузки объекта, у - управляющее воздействие, в качестве которого рассматривается относительная продолжительность включения импульсного преобразователя:
[| к (г) и( }Т )| при \н (]Т )| < Т / к (г)
[ Т при |и()Т)| > Т/к(г),
и(г) при г0 + ]Т < г < г0 + ]Т + у(и(]Т))
7(u(jT)) =
U (t) =
U2(t) при
to + jT + Y(u( jT)) < t < to + (j +1) T,
где и(уТ), к(/) - входной сигнал и коэффициент передачи импульсного преобразователя электромеханического модуля соответственно; у(и(у'Т)) - длительность включенного состояния импульсного преобразователя; Т - период импульсного преобразователя; Щ /), и*( /) - выходные напряжения импульсного преобразователя, модулируемые по требуемому закону на соответствующих интервалах времени.
Вектор измеряемых переменных в системе управления ЭМО формируется следующим образом:
Л = Сх,,
где С - тхн - матрица, состоящая из нулей и единиц и характеризующая включение в состав технически измеряемых компонент вектора х, состояний объекта.
Для решения задачи конструирования системы управления ЭМО требования к качественным показателям формируются в виде функционала от переменных состояния, управления и времени. Распространенной формой записи функционала качества является квадратичная форма, представленная для дискретных ЭМС с одномерным управлением в виде:
I (х, ,и,) = £ (хТах + и,2Я), (2)
,= 0
где компоненты матрицы 2 и постоянной Я характеризуют вклад соответствующих переменных и управления и. Условия минимума этого функционала, с учетом 2=2г>0 и Я>0, позволяют построить управление, стабилизирующее переходные процессы в замкнутой системе регулирования на рассматриваемом интервале времени N. Техническая сложность или практическая невозможность измерения полного вектора х состояния электромеханических объектов и их исполнительных приводов обуславливает необходимость построения адаптивного управления ими с наблюдателем состояния в условиях неконтролируемых возмущений и неполной информации о протекающих процессах.
Синтез алгоритма адаптации
Представляя функционал качества (2) с учетом использования оценки х динамического процесса X в виде суммы двух функционалов 11(а=х-х) и 12(Х, и), построение структуры СУ ЭМО и синтез параметров контуров адаптации можно осуществить в три этапа [5]:
Первый - заключается в синтезе структуры системы управления при решении двух подзадач: минимизация функционала Це) качества оценки X динамического процесса х, по наблюдениям у, и построение управления и, доставляющего минимум функционалу 12(х,, и).
Второй - предполагает построение алгоритмов адаптации в структурном пространстве системы управления с помощью минимизации функционалов /1(е,) и 1г(х„ и) в пространстве параметров ЭМО:
Л ^ V?,/? ,ЛД Г.
где ур, у/1, \уД - подвектор строки структурных параметров регулятора, наблюдающего устройства и датчиков соответственно.
Третий этап - заключается в выборе алгоритмов настройки ур, у/1, до5:
/,P =/ (01P . [е,01P O(x,,Ut)P],
/(н =/ tMH • [e,MH O(Xt,Ut)H],
/tfl =/ 101Д • 01Д O(x, ,Ut)д ],
где Др, Дн, ДД - матрицы параметров контуров адаптации для соответствующих составных элементов системы управления; O(x,,U)р, O(x,,U)н, O(x,,U)A -условия оптимальности функционалов I1(st) и I2(x,,U). Здесь определяются такие компоненты матриц Др, Дн, ДД, которые обеспечивают асимптотический переход периферийных значений функционала I2(x, U) к его экстремальным значениям.
Предлагаемая структура (рис. 1) системы управления ЭМО позволяет реализовать множество алгоритмов с возможностью анализа условий декомпозиции эталонной модели (наблюдателя состояний). На структурной схеме дополнительно обозначены векторы: xM0 - начальных состояний эталонной модели; xMt - состояний эталонной модели; АКЙ - приращений перенастраиваемых параметров регулятора.
Анализатор формирует необходимые значения логических сигналов pM, pr, pOC, которые осуществляют переключение структур регулятора, эталонной модели, устройств обратной связи при выполнении условий устойчивости многосвязного ЭМО.
Раздельное решение задачи оптимизации движения ЭМО на первом этапе предполагает включение блока задержки, которое на интервале времени [ t0,N1] отключает управляющее воздействие, поступающее с выхода широтно-импульсного преобразователя на исполнительные приводы ЭМО. На этом промежутке времени в наблюдателе состояния происходит формирование наилучшей оценки xt состояния xt по измерениям у,, и на интервале [N1,N| подключается управляющее воздействие U на входы исполнительных приводов и наблюдателя состояний ЭМО, который исполняет роль эталонной модели.
Следовательно, для уменьшения или устранения нежелательных переходных процессов в ЭМО, вызываемых неадекватностью начальных условий необходимо ввести блок задержки. Длительность интервала N1 при формировании наилучшей оценки вектора xt состояния по измерениям у, можно определять в результате экспериментальных исследований ЭМО, либо на основе сравнений компонент векторов yt и у==Cxt. Причем второй способ предполагает наличие в у переменной с максимальным временем переходного процесса. Данный способ отыскания N1 имеет ограниченные возможности при значительных изменениях параметров ЭМО. Поэтому актуальной задачей становится задача построения алгоритма функционирования блока задержки, использующей значения измеряемых переменных, параметров наблюдателя и регулятора.
Рассмотрим упрощенную модель ЭМС вида
X+i = Fx, + GU , (3)
и функции Ляпунова квадратичной формы
Рис. 1. Структура адаптивной системы управления ЭМО с разделением процессов наблюдения и управления
¥1(е,) = еХ е<,
л л Т л
У2( X ,) = X , К р х,,
где Кн, КР - матрицы коэффициентов усиления наблюдателя состояний и регулятора ЭМС.
Предположим, что наблюдатель состояния характеризуется матрицей Ь параметров соответствующей размерности. Тогда, записывая первую разность Д^(е,) и ДГ2(х() соответствующих функций Ляпунова:
Таким образом, асимптотические свойства переходных процессов в наблюдателе состояния выполняются при следующих неравенствах:
(^ - О Кр)-1 > (^ - О Кр)Т, ХхТ, [Кр - (^ - О КР)Т Кр(^ - ОКр)] Хг --(у, - сХ,)ТЬТКр Ь (у, - сХг) > 0.
Последнее неравенство можно записать в следующем виде:
т
АГ1(е,) = е[+1 + Кн в,+1 - вт, Кн в,
(4)
А V, (х2) = хмКр х,+1 - х, Кр х,. (5)
Учитывая ур. (3), уравнения (4, 5) запишем в следующем виде:
ДV(в,) = -вТ[Кн - (^ -ЬС)ТКн(^ -ЬС)] е,
х, [Кр - (^ - О Кр )Т Кр (^ - О Кр)]) хх, (у - Сх1 )ТЬТКр Ь(у - Сх).
(7)
(6)
АК, (х,) = -х, [Кр - - ОКр)Т Кр (^ - О Кр)] х + +вТСТЬТК„ЬС в.
Из ур. (6) получаем условие асимптотической сходимости процесса восстановления:
(^ - ЬС)-1 > (^ - ЬС )Т.
Неравенство (7) является алгоритмом работы (рис. 1) устройства переключения, при выполнении которого осуществляется подключение управляющего воздействия и исполнительным приводом в момент минимальной оценки значений вектора состояний ЭМС.
Пример моделирования
В качестве объекта исследования для имитационного моделирования алгоритмов управления использовался электропривод постоянного тока с
Переходные процессы тока / и скорости п в электроприводе при скачкообразном уменьшении неконтролируемого /Д на 50 % от /Д НОМ без подключения контура адаптации
Относ. величины 40
32
24
16
8
0
-8
-16
0,000
"зад
. изменение
0,188
0,375
0,563
0,750 1.с
Рис. 3. Переходные процессы тока / и скорости п в электроприводе при скачкообразном уменьшении неконтролируемого /Д в 15 раз от /Д НОМ с адаптивным регулятором
з
электродвигателем типа ДК-1-2.3 и Г-образным силовым фильтром на входе полупроводникового преобразователя, параметры которого приведены в таблице.
Таблица. Параметры электропривода
Параметры электродвигателя ДК-1-2.3 Значение
Индуктивность якорной обмотки, мГн 2,24
Активное сопротивление цепи якоря, Ом 0,25
Конструктивная постоянная 0,035
Момент инерции якоря, кг/см2 2,3
Магнитный поток, мВб 1,75...2,20
ЭДС двигателя при пНОм, В 32...41
Напряжение питания, В 45
Номинальный ток якоря, А 7,5
Номинальная частота вращения пНОМ, об/мин 1000
Индуктивность силового фильтра, мГн 2
Активное сопротивление силового фильтра, Ом 0,02
Емкость конденсатора силового фильтра, мкФ 4000
Период коммутации преобразователя, мс 0,33
На примере электропривода постоянного тока синтезированный алгоритм перенастройки коэффициентов передачи по току К и скорости Хш дви-
гателя адаптивного регулятора является высокоэффективным способом компенсации влияния такого возмущения, как момент инерции /Д нагрузки двигателя. Так, при скачкообразном уменьшении /д в 15 раз от номинального значения /Д НОМ практически полностью исключается колебательный характер процессов стабилизации тока I и скорости п двигателя (рис. 2, 3).
При моделировании алгоритмов управляющее воздействие электропривода и значения перенастраиваемых параметров адаптивного регулятора для будущего интервала дискретности импульсного преобразователя вычисляются на текущем интервале.
Выводы
Предложен метод синтеза структуры системы управления электромеханическим объектом с использованием функции Гамильтона. За счет организации контура адаптации в несколько раз повышается эффективность регулирования: точность управления с параметрическими и внешними возмущениями возрастает в 4 раза.
СПИСОК ЛИТЕРАТУРЫ
1. Афанасьев В.Н., Колмановский В.Б., Носов В.Р. Математическая теория конструирования систем управления. - М.: Высшая школа, 1998. - 574 с.
2. Красовский А.А. Алгоритмические основы оптимальных адаптивных регуляторов нового класса // Автоматика и телемеханика. - 1995. - № 9. - С. 104-106.
3. Борцов Ю.А., Поляхов Н.Д., Путов В.В. Электромеханические системы с адаптивным и модальным управлением. - Л.: Энер-гоатомиздат, 1984. - 216 с.
4. Борцов Ю.А., Федоров С.В. Адаптивные электроприводы и следящие системы // Электротехника. - 1993. - № 7. - С. 4-8.
5. Букреев В.Г., Параев Ю.И. Адаптивные регуляторы в дискретных системах управления сложными электромеханическими объектами. - Томск: Изд-во Том. ун-та, 2000. - 278 с.
УДК 621.314:658.512
РЕШЕНИЕ ЗАДАЧ ДИНАМИКИ ЭЛЕКТРОМЕХАНИЧЕСКИХ СИСТЕМ В СРЕДЕ АВТОМАТИЗИРОВАННОГО МОДЕЛИРОВАНИЯ МАРС
Т.Н. Зайченко
Томский государственный университет систем управления и радиоэлектроники E-mail: [email protected]
Рассматриваются вопросы применения отечественной универсальной среды автоматизированного моделирования МАРС для решения задач динамики электромеханических систем. Приведены формализованное представление электрической и механической частей электромеханической системы, примеры структурного и схемотехнического моделирования.
Введение и постановка задачи
На кафедре теоретических основ электротехники Томского государственного университета систем управления и радиоэлектроники ведутся работы по созданию комплекса программ автоматизации функционального проектирования электротехнических устройств и систем [1]. Их основой являются:
- метод компонентных цепей (КЦ) Е.А. Арайса и В.М. Дмитриева [2] как теоретическая база автоматизации моделирования неоднородных технических устройств;
- универсальная среда автоматизированного моделирования МАРС (Моделирование и Автоматический Расчет Систем) для Windows [3], являющаяся программной реализацией теоретических основ.
В общем случае математическая модель КЦ состоит из трех типов уравнений относительно переменных связей КЦ: топологических (для потоковых переменных), компонентных и базового узла (для потенциальной переменной; присутствует только в моделях КЦ с энергетическими связями) [1, 2]. В свою очередь среда (система) МАРС включает две части: инвариантную - методы формирования и решения модели КЦ и объектно-ориентированную -библиотеки моделей компонентов соответствующей предметной области. Таким образом, среда МАРС позволяет реализовать моделирование электромеханических систем (ЭМС) при наличии в ее составе моделей элементов и функциональных узлов ЭМС.
Решение задач динамики в системе МАРС рассматривалось в работах Е.А. Арайса, В.М. Дмитриева, А.В. Шутенкова, Л.А. Арайс [2, 4, 5], но применительно к механическим системам. В них использова-
лись модели электроприводов на уровне идеализированных моделей: источников скорости и силы (момента) либо постоянной величины, либо изменяющейся во времени по определенному закону (трапецеидальному, синусоидальному). Такой подход неприменим для решения задач динамики ЭМС, где закон изменения момента и скорости электрической машины (ЭМ) заранее неизвестен, а его определение является одной из целей моделирования.
Целью настоящей работы является развитие объектно-ориентированной составляющей (библиотеки моделей компонентов) системы МАРС для решения задач динамики ЭМС. Разработка модели компонента включает выбор математических моделей функционального и визуального аспектов элемента ЭМС и их программную реализацию в системе МАРС (язык Microsoft Visual C++, среда разработки Visual Studio.NET). При этом математическая модель функционального аспекта отражает процесс функционирования элемента, а модель визуального аспекта - модель условного графического обозначения (УГО) на чертеже КЦ и в дереве раздела библиотеки моделей компонентов среды МАРС.
Структурное моделирование ЭМС
Структурное моделирование ЭМС базируется на математических моделях ЭМ и ЭМС, формализме структурных схем и методах теории автоматического управления. Для структурного моделирования ЭМС необходима реализация моделей элементов структурных схем и математических блоков. Модели компонентов линейных систем автоматического управления (интегрирующее и пропорциональное звенья, сумматор, регуляторы, фильтр) были разработаны и представлены в статье [6]. Для
