Задача синтеза систем управления является одной из важнейших в теории управления. Различным способам синтеза алгоритмов управления посвящено большое количество работ. Зачастую они посвящены синтезу законов управления для линейных систем [1,2], что не всегда приемлемо для нелинейных объектов.
Одним из рассматриваемых в данной работе способов разработки алгоритма управления, является метод синтеза системы управления в пространстве состояний.
Расчет параметров настройки регулятора при помощи методов модального управления предполагает изменение собственных чисел матрицы объекта (мод) для достижения желаемой цели управления. Суть такого метода заключается в приравнивании желаемого и реального характеристического уравнения замкнутой системы для вычисления из полученного соотношения параметров
регулятора [3]. |
̇= + , |
(1) |
|
Объект управления рассматривается в форме: |
|
||
где x — n-мерный вектор |
= , |
|
|
|
|
состояния; y — l-мерный вектор выходных (ре- |
|
гулируемых) переменных; u — m-мерный вектор управляющих воздействий; A
— матрица, определяющая динамические свойства объекта управления, размерности n × n; B — матрица входа управляющих воздействий размерности n × m; C — матрица выхода размерности l × n; n — порядок объекта управления.
Для решения поставленной задачи необходимо перейти от исходного |
||||||||
полинома: |
|
|
|
|
̅= ̅+ |
|
||
описания объекта управления к преобразованной системе |
(2) |
|||||||
Матрица А̅ |
является |
сопровождающей матрицей характеристического |
||||||
|
|
|
̇ |
̅ |
|
|
||
|
|
|
φА(λ)=λn+an-1 |
λn-1+ an-2 λn-2+…+a1λ+a0 |
(3) |
|||
Связь между описаниями (1) и (2) определяется матрицей перехода: |
||||||||
|
|
|
|
|
Pu= ̅MyMy-1, |
(4) |
||
где ̅My – матрица управляемости преобразованной системы, |
My – матри- |
|||||||
ца управляемости исходной системы. |
|
|
|
|||||
Управляющее воздействие формируется как линейная функция всехпе- |
||||||||
где k- |
|
|
= −( 1 1 |
+ + ) = − ̅ |
|
|||
ременных состояния: |
|
|
|
|
, |
(5) |
||
ра. |
матрица размерностью 1 × n параметров настройки регулято- |
|||||||
|
|
|
|
|
|
|
|
|
Дифференциальные уравнения системы получаются в результате подстановки выражения (5) в выражение̇= ((1):+ )
(6)
Желаемая сопровождающая матрица А̅полинома φA(λ) имеет вид:
150
̅ |
0 |
1 |
|
0 |
|
… |
0 |
, |
|
0 |
0 |
|
1 |
|
… |
0 |
(7) |
||
= |
0 |
0 |
|
… |
|
… |
1 |
|
|
|
|
0 |
|
|
|
||||
|
− |
− − … − |
|
||||||
в которой элементами0 |
последней строки являются коэффициенты харак- |
||||||||
|
1 |
|
2 |
|
−1 |
|
|
||
теристического полинома φA(λ) с противоположными знаками. |
|
|||||||||||
Матрица-столбец В̅состоит из нулей, кроме единицы в последней строке, |
||||||||||||
и запишется: |
|
|
|
= 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(8) |
||
|
|
|
0 |
1 |
0 |
0 |
|
… |
0 |
|
||
|
|
|
1 |
|
|
|
||||||
Матрица системы А-ВК будет |
иметь вид: |
|
|
|
|
|||||||
… |
1 |
|
… |
0 |
|
(9) |
||||||
А + ВК = |
0 |
0 |
|
|
|
|||||||
0 |
0 |
|
|
… |
|
… |
1 |
|
||||
|
|
|
|
|
0 |
+ 2 |
|
|
||||
Искомые |
|
|
0 + 1 |
1 + 2 |
2 |
−1 |
+ −1 |
|
|
|||
А+ВК: |
коэффициенты регулятора находятся из равенства матриц А̅ и |
|||||||||||
|
|
|
ki=ai-1-a̅i-1 . |
|
|
|
(10) |
|||||
|
|
|
|
|
|
|
||||||
Параметр, определяющий радиус распределения корней характеристиче- |
||||||||||||
ского полинома ω0 |
находится из соотношения: |
|
|
|
|
|||||||
|
|
|
|
|
|
0 |
, |
|
|
|
(11) |
|
|
|
|
|
время0 |
регулирования в системе с эталонным |
|||||||
где τ0 – безразмерное |
= |
р |
|
|
|
|
|
|||||
характеристическим полиномом; τр – время регулирования в реальной системе. Таким образом, система управления объектом (1) будет включать в себя
управляющее воздействие в виде пропорционального регулятора. Пропорциональные регуляторы, как правило, не обеспечивают отслежи-
вания установок или возмущающих воздействий по нагрузке в системе. Таким образом, желательно переформулировать постановку задачи синтеза таким образом, чтобы в управлении появилась интегральная составляющая, снижающая ошибки регулирования [4].
Существует много способов введения интегральной составляющей. Один из них основан на расширении состояния, путем добавления к n первоначаль-
ным p новых переменных z(t), где |
|
|
|
(12) |
||
|
|
|
|
|||
Расширенное состояние |
будет являться (n+p) - мерным вектором: |
|
||||
|
̇= |
|
|
|||
|
= −1 |
|
. |
(13) |
||
|
̅ − ∫ |
|
|
|||
Соответствующий |
оптимальный регулятор имеет вид: |
|
||||
= … |
|
|
||||
|
|
|
|
2 |
|
|
151
Для существования такого регулятора необходимо, чтобы число переменных состояния, для которых вводится интегральное управляющее воздействие, не превышало числа переменных управления. [4]
Так как химические реакторы, как правило, являются многомерными, многосвязными и нелинейными объектами управления, то системы управления, разработанные методами в линейной постановке, не всегда обеспечивают поддержание оптимального режима работы. Поэтому в данной работе рассматривается синтез системы управления на основе аналитического конструирования агрегированных регуляторов (АКАР) [5,6].
Суть синергетической концепции управления заключается в формировании в фазовом пространстве объекта искусственных аттракторов – притягивающих инвариантных многообразий. Все траектории движения замкнутой системы «объект-регулятор» стягиваются на инвариантное многообразие за счет соответствующих законов управления [7].
Для решения поставленной задачи требуется найти закон управления:
u(ψ)=u(x1,…xn), |
(14) |
который обеспечивает перевод изображающей точки системы из произвольного |
|
начального состояния сначала в окрестность многообразия |
|
ψ(x1,…xn)=0, |
(15) |
в пространстве координат x1,…xn , а затем дальнейшее движение вдоль этого многообразия в некоторую стационарную точку [8].
При создании систем управления также часто используются астатические законы управления, включающие в себя интегральные составляющие, что улучшает точностные свойства систем в установившемся режиме движения. [8] В данной работе на примере химического реактора решается задача синтеза и анализа алгоритмов управления концентрацией целевого компонента для
системы управления в пространстве состояний, а так же при помощи синергетической теории управления, с использованием методов компьютерного моделирования.
Описание объекта и постановка задачи управления |
|||||
В химическом реакторе емкостного типа с мешалкой протекает двухста- |
|||||
+ 1 |
→ 1 |
, |
1 |
→ 2 |
, |
дийная экзотермическая реакция: |
1 |
|
2 |
||
|
|
|
|
|
|
где A, B – исходные реагенты; P , P2 – продукты реакции; k1, k2 – константы скоростей. Целевой продукт – P1.
152
Рис. 1. Принципиальная схема химического реактора
Аппарат функционирует в изотермическом режиме. Назначением объекта управления является проведение химической реакции. Цель функционирования
– получение продукта Р1 с заданной концентрацией.
На рис. 1 введены обозначения: ν1, ν2 – расходы потоков исходных компонентов; САвх, СВвх – концентрации исходных компонентов; ν – расход смеси на выходе из реактора; t – температура смеси в реакторе; СА, СВ, СР1, СР2 – концентрации продуктов реакции; V – объем смеси в реакторе.
В соответствии с допущениями математическая модель химического реактора примет вид:
|
|
dC |
A |
|
|
|
|
вх −(ν |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
|
|
|
=ν |
C |
1 |
+ν |
2 |
)C |
A |
+V (−K C C |
|
) |
|
||||||||||
|
dτ |
|
|
|||||||||||||||||||||
|
|
|
|
|
1 |
A |
|
|
|
|
|
|
|
1 |
A B |
|
|
|||||||
|
|
dC |
B |
|
|
|
|
вх |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(16) |
V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
dτ |
|
=ν |
2CВ |
+V (−K1CACB )−(ν1+ν2)CB |
|
|||||||||||||||||||
|
|
dC |
Р1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
V |
|
=V (K C C |
|
−K C |
Р1 |
)−(ν |
1 |
+ν |
2 |
)C |
Р1 |
|
||||||||||||
|
|
|
|
|||||||||||||||||||||
|
|
dτ |
|
|
|
|
1 |
A B |
|
|
2 |
|
|
|
|
|||||||||
Система (16) решается при соответствующих начальных условиях: CA =CA0,CB =CB0,Cр1=CР10 , где CA0,CB0,Cр10 - решение уравнений статики.
Анализ математической модели показывает, что объект управления является нелинейным, многомерным и многосвязным.
Задача системы управления следует из цели функционирования аппарата
– стабилизация концентрации целевого компонента на заданных значениях в условиях действия возмущений.
В качестве управляющего воздействия для регулирования концентрации целевого компонента целесообразно использовать расход входного потока ν1.
Исследование системных свойств объекта
Важным этапом исследования химического реактора как объекта управления является исследование таких свойств, как устойчивость, управляемость, наблюдаемость и построение статических и динамических характеристик.
153
Анализ статических и динамических характеристик показал, что объект нелинеен по большинству каналов, в том числе по каналу управления.
Анализ устойчивости показал, что объект обладает свойством устойчиво-
сти свободного движения, так как корни характеристического уравнения (λi) |
|||||
|
|
|
|
|
− . . |
имеют отрицательные вещественные части: |
|||||
Были исследованы свойства |
|
= −. |
|||
варианта, где измеряется только |
|
|
−. |
||
имеют полный ранг. |
|
|
управляемости и наблюдаемости для такого |
||
|
С = ( ) |
|
|||
|
|
|
концентрация целевого вещества, и матрица |
||
измерений имеет вид: |
|
|
. Матрицы управляемости и наблюдаемости |
||
|
Таким образом, объект полностью наблюдаем и полностью |
||||
управляем в пространстве состояний.
Синтез системы управления в пространстве состояния
Поскольку математическая модель объекта нелинейна, была проведена
процедура линеаризации в окрестности рабочей точки для получения модели в |
|||||||||||||||||
|
|
|
|
|
−. |
− . |
|
|
|
|
|
|
|
|
|||
пространстве состояний (1) и получены матрицы А и В: |
|
|
|
|
|
|
|||||||||||
|
|
А = |
|
− . |
− . |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
. |
. |
|
|
− . |
|
|
|
|
|
|||
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
Управляющее |
|
|
В = − . |
|
|
|
|
|
|
|
|
|||||
|
Зададимся |
|
|
|
|
−. |
|
|
|
= − |
|
|
|
|
(17) |
||
|
|
= −( + + ) |
|
|
|
|
|
||||||||||
|
|
|
воздействие (6) будет иметь вид: |
|
|
|
|
|
|
||||||||
|
|
желаемым характеристическим полиномом (3) в силу того, что |
|||||||||||||||
порядок объекта n=3: |
D(s)=s3+3s2ω0+3sω02+ω03=0 |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
(18) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Матрица A̅(7) и B̅(8) будет иметь вид: |
|
|
|
|
|
|
|
|
||||||||
|
А = |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
−. − |
− . − |
|
− . − |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
. |
|
|
|||||
|
|
|
. |
|
|
|
|||||||||||
|
Получили матрицу |
преобразований Pu из соотношения (4): |
|
|
|
|
|||||||||||
|
|
|
В = |
|
|
|
|
|
|
|
|
|
|
||||
|
= |
|
|
|
|
|
|
|
|
− . |
|
|
|||||
|
Из соотношения−.(11) получили−параметр. ω |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
||||
τ |
=610 мин. |
|
. |
|
|
0= |
|
при τ0=6 мин и |
|||||||||
р |
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|||
154