CC BY
36
СПИСОК ЛИТЕРАТУРЫ
1. Беллуян З.А. Анализ выпускаемых синхронных генераторов мощностью до 100 кВт // ИТУ (сборник), Ереван, 2002. -№4-1. - С. 101-105.
2. Налимов В.В., Чернова Н.А. Статистические методы планирования экспериментальных исследований. - М.: Наука, 1965.
3. ОСТ 16 0.801.218-84. Машины электрические вращающиеся от 63 до 355 мм габарита включительно. Генераторы синхронные
явнополюсные высокоскоростные. Методика ускоренных испытаний на надежность.
4. ОСТ 16 0.800.882-81. Генераторы синхронные с высотой оси вращения от 56 до 355 мм (мощностью от 0,5 до 100 кВт). Надежность. Методика ускоренных испытаний.
5. Беллуян З.А. Ускоренные испытания на надежность обмоток статора и ротора синхронных генераторов // Известия Томского политехнического университета. - 2004. - Т. 307. - № 2. -С. 141-143.
УДК 621.31-5
ОЦЕНКА РОБАСТНОСТИ АЛГОРИТМОВ УПРАВЛЕНИЯ НЕСТАЦИОНАРНЫМИ ЭЛЕКТРОМЕХАНИЧЕСКИМИ ОБЪЕКТАМИ
В.Г. Букреев, И.Ю. Краснов
Томский политехнический университет E-mail: [email protected]
Осуществлен синтез робастных алгоритмов управления нестационарными электромеханическими объектами, произведена оценка их чувствительности.
Введение
Принципиально важным свойством адаптивных методов обработки информации является робаст-ность, под которой понимается статистическая надежность метода и нечувствительность результатов к изменениям условий наблюдения. Основной особенностью рассматриваемых объектов является их нестационарность: параметры объектов с течением времени изменяются в определённом диапазоне. Нестационарность вносит принципиальные трудности как в изучение структурных свойств объекта (устойчивости, управляемости и наблюдаемости), так и в разработку алгоритмов оценивания и управления.
1. Описание модели объекта управления
Пусть математическая модель функционирования электромеханического объекта описана системой линейных нестационарных дифференциальных уравнений вида [1]:
x(t) = A(t ) x(t ) + b(t) u (t), x(t 0) = x0, (1)
где x(t) - «-мерный вектор, компоненты которого определяют состояние электромеханического объекта в момент времени t; u(t)__- m-мерный вектор управляющих воздействий; A(t) - матрица параметров объектов управления размерности (nxn), b (t) - матрица влияния управляющих воздействий размерности (nxm); x0 - начальное состояние электромеханического объекта в момент времени t0.
Так как сигнал с ПЭВМ поступает в дискретные моменты времени, дискретная модель, соответствующая непрерывной модели (1), имеет вид:
x(k +1) = A (k ) x(k ) + b(k ) u (k ), x(0) = x0, (2)
где матрицы A(k) и b(k) рассчитываются по формулам [2]:
A (k) _£
Atm A{th )m m!
~ _^A.mA(tk )mb(tk ) b(k (m-1)! •
(3)
Если обозначить через А(к) и Ь(к) суммы N первых членов рядов (3) соответственно, то матрицы Ад(к) и Ь1(к) аппроксимируют А(к) и ¿(к) с погреш-ностью~порядка 0(Д/Л). При этом число слагаемых в А(к) и ¿(к) можно задавать заранее (вычисление с фиксированной точностью) или определять автоматически с помощью соотношения:
AN (к )|| -I An-i(k )|| b (k )|| -|\bN-i(k )|||
A N (k
bN (k
<s, (4)
где е выбирается из условия обеспечения максимальной точности вычислений на ПЭВМ конкретного типа.
Такт к в (2-4) соответствует моменту времени 4=0+кД/, Д/- период дискретности по времени.
2. Структура широтно-импульсного модулятора (ШИМ)
В качестве управляющего воздействия рассматривается выходное напряжение импульсного преобразователя, модулированного по широтно-им-пульсному закону и поступающее на исполнительный элемент электромеханического объекта. Ши-ротно-импульсный модулятор (ШИМ) описан следующими уравнениями:
m_0
u(k +1) =
ЭТU(k) при U(k) < —,
T
при u(k) > —;
ЭТ
U (k) =
\Un при t0 + kT < k < t0 + kT + u(k), lo при t0 + kT + u(k) < k < t0 + (k +1)T,
где и(к) - входной сигнал в такт времени /0+кД/, ЭТ и Т- коэффициент передачи и период импульсного преобразователя электромеханического модуля соответственно, Пп - напряжение питания исполнительной части объекта управления.
3. Синтез регуляторов функционирования
системы управления объектом
Для увеличения качества процесса управления объектом (увеличения быстродействия и точности достижения заданного значения) синтезируем несколько типов регуляторов.
3.1. Пропорциональный регулятор
на основе функций Ляпунова 2-го порядка
Параметры регулятора будем определять, исходя из минимизации квадратичного интегрального критерия:
т
I = | ет (<) бе(<) <И,
<0
где в(1)=уш1(1)-у(1) - ошибка регулирования, Т - знак транспонирования, 2 - диагональная матрица (их и), элементы которой определяются аналогично с п. 3.2.
Для получения необходимой плавности движения в окрестности заданной точки позиционирования и асимптотической устойчивости возмущённой системы управления синтез управления приводом может быть осуществлён на основе второго метода Ляпунова.
Используя дискретное представление объекта (2) и функцию Ляпунова квадратичной формы:
V(x(k)) = x(k)T Px(k),
(5)
Матрица Р определяется из условия минимума критерия (6):
Е(к)т РЕ(к) - Р = -б.
Построенный закон управления обеспечивает наиболее быстрое убывание функции Ляпунова (5), определённой на решениях системы (2).
3.2. Пропорциональный регулятор
на основе функций Беллмана
Пусть модель объекта задана в виде (2). Будем считать, что вектор состояния х измеряем полностью.
Функционал качества имеет вид [3]:
N -1
3 = Х (х(к )т б (к) х(к) +
к =0
+и (к)тЯи (к)) + х(N)т Л х(N), (6)
где Л - неотрицательно определенные симметрические матрицы размера (ихи), Я - положительно определенная симметрическая матрица (дхд). Требуется найти управление и(к) с полной обратной связью, минимизирующее функционал (13). Уравнение Беллмана принимают вид:
~х(к )тбх (к) + й(к )тК (к) и (к) + +В (к +1, А(к) х + В (к) и (к))
J (k, x) = min
(7)
где Р - положительно-определённая (ихи) - матрица, управление и(к) будем находить из условия минимума первой разности функции У(х(к)):
д V (х(к)) = V (х(к +1)) - V (х(к)),
вычисленной на решении системы (2). Минимизируя (5) по и(к) получим: И (к) = -(( Е (к )тРВ(к ))т + +В (к )т РЕ (к)) (В(к )т Р В (к))-1,
где Дк) и В(к) - линеаризованные по методу Эйлера в первом приближении матрицы:
Е (к) = I + А (к) т, В(к) = Е (к) Ь(к).
Функция Беллмана J(k, x) определяется по формуле:
J(k, x) = xT (k)P(k) x(k), (8)
где P(k) - неизвестная неотрицательно определенная симметрическая матрица размера (их и). Получаем, что в задаче (6-8) оптимальное управление определяется соотношением:
U (k) = K (k) x(k),
k = 0,1,...,(N-1),
где K(k) - матрица коэффициентов усиления регулятора размера (qxn):
K (k) = [Л + B(k )TP(k +1) B(k )]-1 х
xB(k)TP(k +1) A(k), k = 0,1,...,(N -1),
а матрица P(k) размера (их и) удовлетворяет уравнению:
P(k) = Q + K (k )TRK (k) + +(A(k) - B(k) K(k))T P(k +1) (A(k) - B(k) K(k)), k = (N -1),...,1,0, P(N) = Л.
Матрица Q оценивает отклонения реальных переменных состояния от желаемых и представляет собой своеобразный "штраф" за большие ошибки при любом t0<t<tk. Матрица R, будучи всегда положительной, оценивает стоимость управления. Физически она характеризует затрагиваемую энергию на управление. Матрица Л характеризует ошибку управления в конечный момент времени.
Предварительный выбор значений элементов весовых матриц 2 и Я может быть осуществлен с помощью следующих рекомендаций:
- обычно матрицы 2 и Я назначаются постоянными и диагональными, т.е. матрица 2 содержит п ненулевых элементов qn, /=1,2,...,п, а матрица Я—г содержит г элементов у /'=1,2,...,п;
- принимаем, что максимально допустимые отклонения переменных состояния х(к) в любой момент времени вносят в функционал качества одинаковый вклад.
Распространяя аналогичные рассуждения и на отклонения сигналов управления и(к), можно записать:
Л
Чгг =
у
(
у
И
чи Л тах у
(9)
X ЧиХ21 тах = X
Г ] тах, Л ^ '
(10)
Л=1
й(!) = К
е($) +11 е($ )Л +
(11)
где Хр, Ти, Тй - пропорциональный, интегральный, дифференциальный коэффициенты регулятора соответственно, е(/) - ошибка рассогласования.
Запишем (11) в конечных разностях путем замены: =кТк, Д/ - период дискретности во времени:
й(к) = К
е(к) +1£ е(1)Тк +1
е(к) - е(к -1)
Тк
Здесь тах и х тах - максимально допустимые отклонение, определяемое техническим заданием, /-ой переменной состояния (/=2,3,...,п) и /'-ого сигнала управления (/=2,3,...,г) соответственно. Общий вклад максимально допустимых отклонений переменных состояния должен приблизительно соответствовать общему вкладу максимально допустимых отклонений сигналов управления.
- произвольно выбираем значение г11 (например, ги=1) и по формулам (9, 10) вычисляем значения остальных коэффициентов. Полученные значения весовых коэффициентов следует рассматривать как начальные оценки. Если отклонения по всем переменным достигают своего максимума не одновременно, формулы (9, 10) неверно отражают требования, предъявляемые к системе. Поэтому окончательный выбор весовых коэффициентов целесообразно производить после нескольких пробных процедур синтеза и моделирования системы управления.
3.3. Пропорционально-интегрально-дифференциальный регулятор (ПИД)
Наиболее распространенными алгоритмами являются пропорционально-интегральный (ПИ) и ПИД алгоритмы цифрового управления. При правильной настройке эти алгоритмы обеспечивают достаточно хорошее качество управления для большинства объектов промышленной технологии.
Рассмотрим процедуру вывода алгоритма цифрового ПИД-регулятора из соответствующего непрерывного закона, имеющего вид [4]:
( Т , , Л
где к=1,2,3,...,(Л+1), Тк - номер и величина периода квантования.
Отметим, что при достаточно малых периодах квантования цифровой ПИД закон управления обеспечивает почти такое же качество процессов управления, что и исходный непрерывный закон (11). Для того, чтобы эффект квантования по времени мало сказывался на динамику системы цифрового регулирования, рекомендуется выбирать период квантования из соотношения [5]:
Т Т
191 < т <121
15 <1к < 5 '
где Т95 - время достижения выходным сигналом уровня 95 % от установившегося значения при подаче на вход объекта сигнала ступенчатой формы.
Настройку ПИД-регулятора (нахождения параметров Кр, Ти, Т), будем осуществлять, руководствуясь заданной относительной степенью затухания ф [5]. Определяем значение логарифмического декремента затухания т:
- 1п(1 -ф)
т =--——.
2п
Система автоматического управления будет обладать требуемой относительной степенью затухания, если амплитудно-фазовая характеристика разомкнутой системы автоматического управления будет проходить через точку на комплексной плоскости (-1,/0) [6]: Кр=0,409, Ти=8,696, Т,=0,035.
4. Оценка робастности
Проведём анализ синтезированных регуляторов на робастность (реакцию на изменение параметров объекта). Будем считать, что объект является нестационарным, т.е. его параметры изменяются во времени. Положим, что объектом управления является двигатель постоянного тока ДК-1. Переходные процессы по току двигателя без вынужденных возмущений и при увеличении момента сопротивления в момент времени 4=/0+к-Д/, к=200 (заданное значение по скорости равно 1) изображены на рисунке.
Таблица. Характеристики синтезированных пропорциональных регуляторов
Тип регулятора Перерегулирование, % Установившееся значение Время, такты
Без возмущений С на-грузкой Без воз-муще-ний С на-грузкой Без воз-муще-ний С на-грузкой
К=1 1,532 - 0,999 - 144 -
По функциям Ляпунова 0 0 0,995 0,983 126 300
ПИД 0 0 0,994 0,987 140 261
По функциям Беллмана 0 0 0,995 0,984 122 295
и
1 тах
1 тах
г,, =
г
х
V I тах J
1=1
- По фушсциям Белмана
Рисунок. Переходные процессы по току на обмотке двигателя, обеспечиваемые синтезированными регуляторами (без вынужденных возмущений - слева и при увеличении нагрузки на валу в 3 раза - справа)
По результатам, приведённым в таблице, можно заключить, что робастность регуляторов, синтезированных по предложенным алгоритмам, определяется в соответствии с конкретно поставленными задачами. В зависимости от требований, предъявляемых к объекту управления, поставленными перед ним целями и возможными затратами, может быть выбран тот или иной алгоритм управления.
В нашем случае наиболее подходящим является ПИД-регулятор.
Заключение
В современной теории управления нестационарными электромеханическими объектами большое внимание уделяется робастности алгоритмов к действующим возмущениям. В данной статье были исследованы на нечувствительность несколько типов регуляторов, приведены результаты исследования, сделаны рекомендации по поводу использования того или иного регулятора.
СПИСОК ЛИТЕРАТУРЫ
1. Сю Д., Мейер А. Современная теория автоматического управления и ее применение. - М.: Машиностроение, 1972. - 183 с.
2. Григорьев В.В., Дроздов В.Н., Лаврентьев В.В., Ушаков А.В. Синтез дискретных регуляторов при помощи ЭВМ. - Л.: Машиностроение, 1983. - 357 с.
3. Колосов С.П., Калмыков И.В., Нефедова В.И. Элементы автоматики. - М.: Машиностроение, 1970. - 285 с.
4. Семенов В.В., Пантелеев А.В., Бортаковский А.С. Математическая теория управления в примерах и задачах. - М.: МАИ, 1997. - 514 с.
5. Атанс М., Фабл П. Оптимальное управление. - М.: Машиностроение, 1968. - 764 с.
6. Мирошник И.В., Никифоров В.О., Фрадков А.Л. Нелинейное и адаптивное управление сложными динамическими системами. - СПб.: Наука, 2000. - 549 с.
