
MS_bak_220400_220700
.pdf
3.10. Синтез связанной системы управления объектом с двумя связанными параметрами по критерию минимум суммы интегрально-квадратичных ошибок и из условия автономности
Для многомерных объектов с сильными внутренними связями применяются схемы связанного управления (рис. 3.24, 3.25).
Рис. 3.24. Пример связанной системы управления концентрациями изоамиленов и изопентана на контрольных тарелках (к.т. № 41 и № 20) колонны ректификации: QE – датчики состава; QC – регуляторы состава; К12, К21 – компенсаторы перекресных связей
Расчет связанной ЦСУ может проводится двумя способами (с помощью численных методов и из условия автономности).
В первом способе определяются оптимальные настройки регуляторов по основным каналам управления (синтез одноконтурных ЦСУ (см. пункт 3.5)), а затем одновременно оптимизируются настройки цифровых компенсаторов задаваемых структур по комплексному критерию (3.84) (см. пункт 3.9).
121

Рис. 3.25. Структурная схема связанной системы управления:
WP1 ,WP2 ,WK12 ,WK21 ,WO1 ,WO2 ,WO12 ,WO21 - передаточные функции регукомпенсаторов и основных, перекрестных каналов управления; - задающие воздействия; e1 ,e2 - величины рассогласования;
u1 ,u2 ,uk 12 ,uk 21 - управляющие воздействия от регуляторов и ком-
пенсаторов; u1 ,u2 - суммарные управляющие воздействия на объект;
2 - суммарные выходы объекта
Модель связанной ЦСУ объектом с двумя связанными параметрами с последовательным включением основных и компенсирующих регуляторов может быть представлена следующей
системой разностных уравнений: |
|
|
|
|
|
|
|
|
|||||||||
u1 |
u1 |
q1 y31 |
y1 q1 |
y31 |
y1 q1 y31 |
y1 , |
|||||||||||
i |
|
i 1 |
|
0 |
i |
i |
1 |
|
i 1 |
|
|
i 1 |
2 |
i 2 |
|
i 2 |
|
u2 |
u2 |
q2 |
y32 y2 q |
2 |
y32 |
|
y2 |
q2 y32 |
|
y2 , |
|||||||
i |
|
i 1 |
|
0 |
i |
i |
1 |
i 1 |
|
i 1 |
2 |
|
i 2 |
i 2 |
|||
uk12 |
uk12 |
q12 |
u1 |
q12 |
u1 |
|
|
q12 u1 |
|
, |
|
|
|||||
|
i |
|
i 1 |
|
0 |
i |
1 |
|
i 1 |
2 |
i 2 |
|
|
|
|||
uk 21 |
uk 21 |
q21 |
u2 |
q21 u2 |
|
q21 |
u2 |
|
, |
|
|
||||||
|
i |
|
i 1 |
|
0 |
i |
1 |
|
i 1 |
2 |
i 2 |
|
|
|
|||
u1 u1 |
uk |
21 , |
|
|
|
|
|
|
|
|
|
|
|
(3.86) |
|||
i |
|
i |
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
ui2 ui2 uki12 ,
122
y1 |
|
|
a1 y1 |
|
a1 y1 |
|
|
|
b1u1 , |
|
|
|
|
|
|
|||||||
i d1 1 |
1 i d1 |
|
2 i d1 1 |
|
i |
|
|
|
|
|
|
|
|
|||||||||
y12 |
|
a12 y12 |
|
b12 u1 |
, |
|
|
|
|
|
|
|
|
|
||||||||
i d12 1 |
1 |
i d12 |
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
||||||
y2 |
|
|
a2 y2 |
|
b2 u2 |
b2 u2 |
, |
|
|
|
|
|
|
|
||||||||
i d 2 1 |
1 i d 2 |
|
1 |
i |
|
|
|
2 i 1 |
|
|
|
|
|
|
|
|
||||||
y21 |
|
a |
21 y21 |
|
a21 y |
21 |
|
a21 y21 |
|
b21u2 , |
|
|||||||||||
i d 21 1 |
1 |
i d 21 |
|
2 |
|
|
i d 21 1 |
|
|
3 |
i d 21 2 |
|
|
i |
|
|||||||
y1 |
|
y1 |
|
|
y21 |
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|||
i d1 1 |
i d1 1 |
|
|
i d1 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
y2 |
|
y2 |
|
|
y12 |
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|||
i d 2 1 |
i d 2 1 |
|
|
i d 2 1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
где |
u1 ,u |
2 |
|
- |
|
суммарные |
|
управляющие |
воздействия; |
|||||||||||||
|
|
|
i |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1 |
, y2 |
, y12 , y21 - выходы объекта по основным и перекрестным |
||||||||||||||||||||
i |
i |
|
i |
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
каналам; |
|
y1 |
, y |
2 |
|
|
- |
суммарные |
выходы |
объекта; |
||||||||||||
|
|
|
|
|
i |
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a1 |
,a1 |
,b1 ,d1,a2 |
,b2 |
,b2 ,d2 |
и a12 |
,b12 |
,d12,a21 ,a21 |
,a21 |
,b21 ,d21 - |
|||||||||||||
1 |
2 |
|
|
|
1 |
|
1 |
|
2 |
|
|
|
|
1 |
|
|
|
|
1 |
2 |
3 |
|
параметры |
моделей |
объекта; q1 |
,q1 |
,q1 |
,q2 |
,q2 |
,q |
2 - |
настройки |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
1 |
2 |
0 |
1 |
2 |
|
основных цифровых регуляторов, определяемые по основным каналам управления; q012 ,q112 ,q212 ,q021 ,q121 ,q221 - искомые настройки цифровых компенсаторов 2-го порядка.
После расчета оптимальных настроек по первому способу (табл. 3.11) выполняется исследование динамических характеристик связанной ЦСУ (рис. 3.26).
123
|
|
|
|
|
|
Таблица 3.11 |
|
|
|
Параметры оптимальной связанной ЦСУ |
|
|
|||
|
|
|
|
|
|||
|
Параметры моделей |
Оптимальные настройки |
Такт |
|
|||
|
объекта управления |
|
цифровых регуляторов |
квантова- |
|
||
|
|
|
|
|
и компенсаторов |
ния |
|
|
|
1 |
|
|
1 |
сигналов |
|
|
|
WO |
|
|
WP |
|
|
a1 |
1,739410 |
q1 |
0,920878 (т/ч)/% |
T0 4,0 |
|
||
1 |
|
|
0 |
|
|
|
|
|
1 |
0,788627 |
|
1 |
1,653834 (т/ч)/% |
мин |
|
a2 |
q1 |
|
|
||||
b1 |
0,212896 %/(т/ч) |
q1 |
0,768621 (т/ч)/% |
|
|
||
|
|
|
|
2 |
|
|
|
d1 4 |
|
|
|
|
|
||
|
|
W 12 |
|
|
W 21 |
|
|
|
|
O |
|
|
K |
|
|
a12 |
0,961012 |
q |
21 |
0,033432 (т/ч)/% |
|
|
|
1 |
|
|
0 |
|
|
|
|
b12 |
0,044751 %/(т/ч) |
q |
21 |
0,006599 (т/ч)/% |
|
|
|
|
|
|
1 |
|
|
|
|
d12 7 |
q221 0,008321 (т/ч)/% |
|
|
||||
|
|
W 2 |
|
|
W 2 |
|
|
|
|
O |
|
|
P |
|
|
a |
2 |
0,923788 |
q |
2 |
0,603545 (т/ч)/% |
|
|
1 |
|
|
0 |
|
|
|
|
b |
2 |
0,575315 %/(т/ч) |
q |
2 |
0,797134 (т/ч)/% |
|
|
1 |
|
1 |
|
|
|
||
b |
2 |
0,382540 %/(т/ч) |
q |
2 |
0,222656 (т/ч)/% |
|
|
2 |
|
|
2 |
|
|
|
|
d2 2 |
|
|
|
|
|
||
|
|
W 21 |
|
|
W 12 |
|
|
|
|
O |
|
|
K |
|
|
a |
21 |
1,322264 |
q12 |
0,007807 (т/ч)/% |
|
|
|
1 |
|
|
0 |
|
|
|
|
a |
21 |
0108190, |
q12 |
0,001608 (т/ч)/% |
|
|
|
|
2 |
|
1 |
|
|
|
|
a |
21 |
0,451724 |
q12 |
0,000770 (т/ч)/% |
|
|
|
|
3 |
|
|
2 |
|
|
|
b21 0,090779 %/(т/ч) |
|
|
|
|
|
||
d21 10 |
|
|
|
|
|
124

Рис. 3.26. Динамические характеристики оптимальной связанной ЦСУ концентрациями изоамиленов и изопентана на контрольных тарелках колонны ректификации:
125

а) u1 - расход флегмы; u2 - расход теплоносителя в кипятильник ко-
лонны; б) y 31 - задание по концентрации изоамиленов; y1 - концен-
трация изоамиленов; y 32 - задание по концентрации изопентана; y 2 - концентрация изопентана
При расчете связанной ЦСУ по первому способу можно предварительно не выполнять оптимизацию основных цифровых регуляторов, а одновременно оптимизировать настройки основных и компенсирующих регуляторов (компенсаторов). Для этого
формируется общий массив настроек g и выполняется оптими-
зация настроек по критерию (3.84) методом покоординатного спуска. Расчет динамических характеристик ЦСУ проводится по системе уравнений (3.86) при нулевых начальных условиях.
При применении второго способа цифровые компенсаторы определяются из условия автономности (предварительно рассчитываются основные регуляторы в составе одноконтурных схем). Для схемы связанного управления (рис. 3.25) уравнения компенсаторов определяются из следующих условий:
|
|
|
|
|
|
|
|
|
|
|
uk |
12 |
z |
|
|
|||||
|
|
|
|
Wo12 z |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
1 |
|
|
, |
(3.87) |
||||||||
Wk12 z |
|
|
|
|
u |
z |
||||||||||||||
|
|
|
|
|
Wo2 z |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
uk |
|
21 |
|
z |
|
|
|||||
|
|
|
Wo21 z |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
2 |
|
, |
|
|||||||||||
Wk 21 z |
|
|
|
|
|
u |
z |
|
||||||||||||
|
|
|
|
|
Wo1 z |
|
|
|
|
|
|
|
|
|
|
|
|
|
где Wo1 z ,Wo2 z и Wo12 z ,Wo21 z - дискретные передаточные функции основных и перекрестных каналов управления, получаемые по разностным уравнениям из системы (3.34); uk 12 z ,u2 z и uk 21 z ,u1 z - выходы и входы перекрестных компенсаторов.
После проведения аналитических преобразований с помощью оператора сдвига z , получаем уравнения цифровых компенсаторов:
126
|
|
|
|
|
|
uk12 |
|
1 |
uk |
12 |
|
|
2 |
|
uk |
12 |
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
i 2 |
|
|
|
|
|
, |
|
|
(3.88) |
||||||
|
|
|
|
|
|
|
|
u1 |
|
|
|
|
|
|
|
|
|
|
|
|
u1 |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
i 1 d12 d 2 |
|
|
|
|
|
|
i 2 d12 d 2 1 |
|
|
|
|
|||||||||||||||||||
|
|
|
uk 21 |
|
|
1 |
uk |
21 |
|
2 |
uk 21 |
|
|
3 |
uk 21 |
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
i 2 |
|
|
|
|
|
i 3 |
|
, |
||||||||||
|
|
|
u 2 |
|
|
|
|
|
|
|
|
u 2 |
|
|
|
|
|
|
|
u |
2 |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
3 |
|
|
||||||||||||||||||||||||
|
|
|
|
1 |
|
|
i d 21 d1 |
|
|
|
|
|
|
|
|
|
i 1 d 21 d1 |
|
|
|
|
i 2 d 21 d1 |
|
||||||||||||||||||
где |
a12b2 |
|
b2 |
|
; |
|
|
|
a12b2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
1 1 |
|
|
|
|
|
2 |
|
|
|
|
1 |
|
|
2 |
|
; |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
1 |
|
|
|
|
b2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
b2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
b12 |
|
|
|
|
|
|
|
|
a 2b12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
; |
|
|
|
|
|
1 |
|
|
|
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
b2 |
|
|
|
b2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
1 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
a21 ; |
2 |
a |
21 ; |
|
|
3 |
a21 |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
1 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
b21 |
|
; |
|
|
|
|
a11b21 |
|
; |
|
|
a21b21 |
. |
|
|
|
|
|
||||||||||||||||||||||
|
1 |
|
|
|
|
b1 |
|
|
|
2 |
|
|
|
|
|
|
b1 |
|
|
|
|
|
|
|
|
3 |
|
|
b1 |
|
|
|
|
|
|
|
|
После расчета цифровых компенсаторов (табл. 3.12) проводится исследование динамических характеристик связанной ЦСУ
(рис. 3.27).
|
|
|
|
|
Таблица 3.12 |
|
|
Параметры цифровых компенсаторов, |
|
|
|||
|
рассчитанные из условия автономности |
|
|
|||
|
|
|
|
|
||
|
Настройки цифровых компенсаторов |
Такт |
|
|||
|
|
|
|
|
квантова- |
|
|
|
|
|
|
ния |
|
|
|
|
12 |
|
сигналов |
|
|
WK |
|
|
|
||
1 |
0,296089 (т/ч)/% |
|
1 |
0,077785 (т/ч)/% |
T0 4,0 |
|
2 |
0,638999 (т/ч)/% |
|
2 |
-0,071857 (т/ч)/% |
мин |
|
|
W |
21 |
|
|
|
|
|
|
K |
|
|
|
|
1 |
1,322264 (т/ч)/% |
|
1 |
0,426401 (т/ч)/% |
|
|
2 |
0,108190 (т/ч)/% |
|
2 |
0,741686 (т/ч)/% |
|
|
3 |
0,451724 (т/ч)/% |
|
3 |
0,336271 (т/ч)/% |
|
|
|
|
|
|
|
127 |

Рис. 3.27. Динамические характеристики систем управления концентрациями изоамиленов и изопентана на контрольных тарелках
колонны ректификации: y1 * , y 2 * - выходы объекта (последовательная оптимизация регуляторов и компенсаторов (табл. 3.11));
y1 , y 2 - выходы объекта (последовательная оптимизация регуляторов и расчет компенсаторов из условия автономности)
Второй способ обеспечивает лучшее качество управления, но он требует аналитических преобразований, что затрудняет синтез и текущую перенастройку системы связанного управления.
128

3.11. Расчет показателей качества управления по динамическим характеристикам замкнутых систем
После получения дискретного динамического описания объекта управления и расчета оптимальных параметров управляющей части ЦСУ (настроек цифровых регуляторов и компенсаторов выбранных структур) проводится исследование динамических характеристик замкнутой системы управления с целью оценки ее эффективности.
Исследование заключается в расчете показателей качества управления по выходным управляемым параметрам (рис. 3.28).
Рис. 3.28. Динамическая характеристика замкнутой системы при ступенчатом изменении задающего воздействия:
y 3 - задающее воздействие; y - выход объекта
129
Выделяют следующие основные показатели. 1. Интегрально-квадратичная ошибка:
|
|
N |
|
yi 2 , |
|
|
S yi3 |
(3.89) |
|||||
|
|
i mc |
|
|
|
|
где N - общее число тактов; |
mc - номер такта, на котором из- |
|||||
меняется задание регулятору; |
y 3 , y |
i |
- значения задающего воз- |
|||
|
|
|
i |
|
|
|
действия и выхода объекта на i -ом такте. |
|
|||||
2. Время регулирования: |
|
|
|
|
|
|
t |
рег |
i* |
T , |
(3.90) |
||
|
|
|
0 |
|
где i * - номер такта, для которого выполняется проверяемое ус-
ловие при изменении i от |
N до mc с шагом минус единица: |
||||||||||||||||
|
y3 |
y |
i |
|
100 |
% ; 2 |
5 %. |
|
|
|
|
||||||
|
|
|
|
|
|
||||||||||||
|
i |
|
|
|
|
|
|
|
|||||||||
|
|
y3 |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3. Перерегулирование: |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
P |
|
|
A1 |
100% , |
(3.91) |
||||
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
max |
|
y3 |
|
|
|
|
|||
где |
|
A1 |
|
- максимальное |
отклонение |
выхода объекта yi |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
( i mc, N ) от задания y 3 . |
|
|
|
|
|
|
|
|
|
||||||||
Если величина |
A меньше задания y 3 |
(для апериодических про- |
|||||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
цессов), то в этом случае перерегулирование равно нулю. |
|||||||||||||||||
|
|
|
4. Коэффициент затухания: |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
A1 A3 |
100% , |
(3.92) |
||||||
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
A1 |
|
|
|
|
||
где |
A3 - второй максимум выхода объекта yi |
(определяется по- |
|||||||||||||||
сле нахождения минимального значения выхода). |
|||||||||||||||||
При отсутствии перерегулирования ( Pmax |
0 ) |
коэффициент за- |
|||||||||||||||
тухания равен 100 % ( =100 %). |
|
|
|
|
|||||||||||||
|
|
|
5. Статическая ошибка: |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
y |
y3 y |
N |
, |
|
(3.93) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
130 |
|
|
|
|
|
|
|
|
|
|
|
|