Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Живоглядов, В. П. Адаптация в автоматизированных системах управления технологическими процессами

.pdf
Скачиваний:
18
Добавлен:
19.10.2023
Размер:
6.85 Mб
Скачать

4. 2. 4. Алгоритмы коррекции задания локальной системы регулирования

Вследствие статистической связи косвенных показателей q с выходом w стабилизация сигнала Q устройства распре­ деленного контроля, особенно в условиях изменяющихся про­ изводственных ситуаций и колебания размалываемости клин­ кера, не гарантирует высокой точности поддержания ш на за­ данном уровне w*. Повышение эффективности работы систе­ мы может быть достигнуто при использовании текущей ин­ формации о w и коррекции с помощью УВМ задания Q* и па­ раметров а / локальной САР. Для коррекции задания можно применить метод стохастической аппроксимации:

Q*[s] = Q*[s— 1]—тЫ(-даЫ—ш*).

Алгоритм гарантирует асимптотическую сходимость w к иу* при неизвестной динамике объекта, если объект устойчив и возмущение — случайная величина. Скорость сходимости невелика и существенно зависит от выбора коэффициента f. Применим для синтеза алгоритма коррекции задания теорию дуального управления. Поскольку реализация алгоритма предполагается на цифровой управляющей вычислительной машине, все переменные будем рассматривать- в дискретном времени. Заметим, что и контроль тонкости помола цемента возможен в настоящее время лишь в дискретные моменты.

Динамику мельницы по каналу «Задание Q* регулятору загрузки клинкера — тонкость помола ш» приближенно опи­ шем уравнениями:

о;Ы= вЫ + н-Ы,

(4. 12)

0 Ы= (1 —d)0 [s— 1 ]+ Q[s—т] Kd,

где

 

 

A t

 

 

d--

&-

T06

Тоб, тн К — постоянная времени, чистое запаздывание и ко­ эффициент усиления. Корреляционную функцию эквивалент­ ного возмущения ц аппроксимируем экспонентой (3. 36). По­ грешности измерения тонкости помола считаем центрирован­ ными, независимыми, гауссовыми. Функция потерь квад­ ратичная.

Введем в управляющее устройство оператор, компенси­ рующий влияние инерционности объекта:

©*[s]=

-U [s] - m[s-

1] ( -i- —1

,(4. 13)

 

d

\ d

 

При одинаковых

начальных условиях 0 и и имеем

©Ы —

/(«[s—т]. То же

справедливо в стационарном режиме (s —>оо)

при любых начальных условиях. Закон оптимального управ­

ления получен в разделе 1 и определяется формулой

(1. 131).

Перепишем ее в обозначениях данного раздела:

 

и[s]= (ш*[s+ т]—bs-\W* [s+ 1 — 11+&s- i и[s—1:])—

y[s— 1 ]

(4- 14)

— w [ s — т — 1 ])

К

 

Совместное использование (4. 13) и (4. 14) дает искомый алгоритм управления заданием 0* с накоплением информа­ ции. Переменные коэффициенты bs, ds определяются выраже­ ниями (1. 132), (1. 133). Они могут быть заранее вычислены и введены в память УВМ, либо заранее выполняется только часть расчетов по формулам (1. 132) и (1. 133).

Для хранения в памяти машины текущих значений и и у

требуется (-д-^+З) ячейки оперативного запоминающего уст­

ройства.

bs и й5при s-s-'oo..

Вычислим стационарные значения

На основании изложенного в разделе 1

имеем

 

Д- - U r n

r

- 2lAl

 

где

СО

/—»СО

Pjl

 

 

 

 

 

 

 

3i= (l + А х+Р2 A ОхНЛ 1+ AH-PsAi f - 4р2 Ах2,

limps _ Х_ 1Л

1 4~ А1

г -2р A j

.2 1 + Ал

2р2 Aj2

S—>со

2°л2

 

К

2оу

 

6=

 

г

2р Ах

*'(4. 15)

lim bs=

 

й ’

 

5—>оо

■'оо

 

°1

 

 

 

 

 

 

 

 

 

T+l

 

 

 

: ton g s = ____ P

2'p2 A !2

 

(4. 16)

(1 + Д xjSx— 2p2 A i*

S—уoo

1 + Ax

 

 

 

 

 

 

 

 

Полученный алгоритм был испытан в автоматизированной

системе

управления

цементной

мельницей

(размером

2,6Х

X ГЗ м)

Чимкентского цементного завода. Числовые характе­

ристики

объекта определены

экспериментально

и

равны

тн = 15 мин, сг.[а2 = 4

(%)2, й| =0,03 (мин)~1. Дисперсия

слу­

чайной

погрешности

в приборе

контроля

тонкости

помола

равна 0,3. Интервал квантования был принят At = 5 мин.

 

По экспериментальным данным и приведенным выше фор­ мулам найдены значения коэффициентов. Алгоритм был при­

веден к виду

 

 

Qls] = 10(w[s]—0,5a[s—11),

(4.

17)

«Ы =до*[5+т] —0,2bu*[s-H —Ц-f 0,2«[s—1 ] —0,44(y[s— 1] —

—a k - x - l] ) ,

(4.

18)

где y[s], эд*[э] — измеренное и заданное значение тонкости помола в % (остаток на сите 008).

Для реализации алгоритма на ЦВМ необходимо (г+8) ячеек памяти, 5 умножений и 5 сложений, т. е. необходимое количество запоминающих ячеек и вычислительных операции очень незначительно.

При синтезе алгоритма коррекции можно применить и бо­ лее простую модель по каналу «задание Q* — тонкость типа (4. 6). В. этом-случае в стационарном режиме

0'*Ы— -rr-(pS4s-\-i]—bw*[s-\- -—ll)-)-£©*[s — 1] —

Ко

S

y[s— 1 ] — Q*[s — т— 1]

(4. 19)

 

К0

 

Былоустановлено, что алгоритм дуального управления мало чувствителен к изменению модели объекта, и алгоритмы типа (1. 130) и (4. 19) могут быть применены и для инерционных объектов с запаздыванием. Обратим внимание на то, что при

Т

оптимизации по. критерию типаWs= e2[s]-(--^L (e[sl—e[s—I])2

2G2

в канал управления вводится инерционное звено (см. подраз­ дел 1.3. 1). Можно предположить, что процессы в послед­ ней системе п в системе управления инерционным объектом по алгоритму (4. 19) близки.

Несколько отлична от рассмотренной выше методика ин­ ститута ВИАСМ [4. 4]. В развитие алгоритма (1. 136) и (1. 24) для компенсации инерционности мельницы в УУ вве­ дено форсирующее звено (1 + 7 0бР), а для устранения мед­ ленного дрейфа и постоянного смещения выхода — интегра­ тор, после чего по методу z-форм был осуществлен переход от

.непрерывного времени к дискретному. Алгоритм реализован в результате работ института ВИАСМ на УВМ «Тбилиси-1». Распределенный контроль в этой системе отсутствует.

4. 2. 5. Алгоритмы прямого цифрового управления (ПЦУ) загрузкой мельницы

АСУ с прямым цифровым управлением загрузкой реали­ зуется по схеме, представленной на рис. 4. 26.

Управление положением исполнительного механизма ИМ (или отсекающего ножа тарельчатого питателя) можно осу­ ществлять непосредственно по выходу га» [s] без внутреннего контура стабилизации в соответствии с алгоритмами типа (4. 13), (4. 14) или (4. 19). Однако информация о состоянии мельницы, содержащаяся в косвенных показателях q(xi) , при этом не используется. Здесь мы рассмотрим применение ПЦУ во внутреннем контуре стабилизации. Для описания объекта по каналам «загрузка материала — косвенные пока­ затели qK=q(xK)= fqlqA(xK).qB(xK)\>>воспользуемся линеаризован­ ными моделями типа (4. 5) или (1. 148), приведенными в не­ сколько отличных обозначениях. Если предположить, что имеется лишь одно возмущение ц° с экспоненциальной корре­

ляционной функцией типа (3. 36),

то закон управления за­

грузкой задается формулой (1. 156)

и алгоритм ПЦУ в обоз­

начениях данного раздела имеет вид

 

 

s-1 /

 

/=1к=0

-Л -и*и[/-тя])],

(4. 20)

где ук,/]— измеренное значение показателя qK в /-ft момент

203

времени. Если у связан с q некоторым известным взаимно­

однозначным

преобразованием

ybc,j]—f(qlK]+h[K,j]),

то>

вместо y[K,j] в (4. 20) следует

подставить / у-1(у[и,Л),

где

Ду-1— функция,

обратная f у-

 

для

Весовые коэффициенты ак введены в алгоритм (4. 20)

настройки при изменении характеристик объекта и возму­ щающих воздействий. Если некоторые q(x, к) не измеряются,, достаточно положить 11^=0 в формуле (1. 157) для вычисле­ ния qKj . Блок-схема алгоритма ПЦУ для произвольного чи­ сла дискретных датчиков контроля состояния мельницы по­ казана на рис. 4. 5. Практически достаточно ограничиться двумя-тремя индукционными датчиками у первой и второй ка­

мер мельницы. Обозначено: а\,ак,Ки1,Кик,К ----множитель-

'Аи

ные звенья на соответствующие коэффициенты; Ф — дискрет­

ный линейный фильтр с весами qKj \

ИМ — исполнительный

механизм

(позиционер); т1, ... хк

-- задержки

на время

хк=хк—тк_!.

Приближенно Ф можно заменить фильтром пер­

вого порядка, на вход которого поступает сумма У

М У[ас,/] —

- К икЦ - х к}).

Эффективность алгоритма подтверждается промышлен­ ными испытаниями, проведенными с использованием УВМ ти­ па «Днепр-1» на мельнице размером 3X4 як

Приведем еще один алгоритм прямого цифрового управ­ ления загрузкой, который получим из (3. 35) путем дискрети­ зации с интервалом At:

M[s]= aju[s—l]-[-aae[s]—a3e[s— 1],

(4.

21)

<?[s]= Q*[s]—Q[s] = Q*[s]—^jixKy[K,s],

(4.

22)

где

к

 

 

^из(^“ЬКр7Г0С)

 

 

a _

 

 

1

TH3( l + K pK 0C) + K 0CAt ''

 

 

 

тпзкр+ дt

 

 

fl2= ТНЗ(\+ К РКОС)+КОСА ( -

 

 

_

_________7U3 ^ ________

 

 

3

^из(1+^р^ос) + ^ос Д ^

 

 

204

Сигналы дотиикоб

Рис. 4.5

■ш.

Простота реализации на ЦВМ обоих ПЦУ очевидна. Дей­ ствительно, для управления по алгоритму (4. 21) при двух датчиках заполнения мельницы (^i и q2) требуется 11 ячеек памяти, 5 умножений и 4 сложения. При этом мы не учитыва­ ем загрузку машины программами контроля, опроса датчи­ ков, предварительной обработки сигналов, выполнения раз­ личных логических операций, программой-диспетчером и другими служебными программами, необходимыми для нор­ мального функционирования УВМ в системе управления комплексом технологических агрегатов.

4. 2. 6. Алгоритмы адаптации

Общий подход к синтезу алгоритмов адаптации для си­ стем с распределенным контролем изложен в разделе 3. Здесь мы приведем несколько модифицированный алгоритм наст­ ройки вектора а коэффициентов локальной системы управле­ ния процессом измельчения. При квадратичной функции по­ терь он имеет вид

a[s]=a[s — 1]— r i[s]G[s](co[s]— u>*)—Г2[s]G[s —l](m[s—1]—до*).

(4. 23)

Вектор G[s] получаем в результате преобразования дискрет­ ным динамическим фильтром с запаздыванием (моделью чувствительности) вектора сигналов косвенных показателен в схеме рис. 3, 6, либо сигналов, поступающих на входы бло­ ков ак в схеме рис. 4. 5. Матрицы переменных коэффициентов обозначены /’lfs] и r 2[s].

На рис. 4. 6 показан процесс настройки двух коэффициен­ тов ai и а2 при моделировании АСУ процессом помола на ЦВМ. Модель объекта состояла из ряда звеньев с запазды­ ванием. Возмущение аппроксимировали двумерным марков­

ским процессом. Полученные при

моделировании оценки

R6 (т) второго начального момента

A4(e[s]e[s—т]}

ошибки е

приведены

на рис. 4. 7. Кривая 1

соответствует

режиму

a 1= a2= 0)

и = 0, а кривая 2 характеризует качество управле­

ния при ai = a2=0,4, лежащих в окрестности точки оптимума, т. е. когда процесс настройки в основном завершен. Из ри­ сунка видно, что управляющие воздействия u[s] компенсиру­ ют низкочастотные составляющие возмущений, приложенные ж объекту.

Изложим некоторые результаты производственных испы­ таний алгоритма адаптации и управления цементной мельни*

-206

/>о

to

о Рис. 4.6

ч

Рис. 4.7

цей с использованием УВМ типа «Днепр». Модели чувст­ вительности были аппроксимированы звеньями с чистым за­ паздыванием на время 30 мин. Интервал дискретности изме­ рения выхода принят At=5 мин, а матрицы Л Ы , A [s] в алгоритме (4. 23) выбраны вида

Гi[s]=diag{0,001; 0,001}, r a[s]= diag{0,0005;t0,0005b

•Испытания проходили при значениях задания w*, соответ­ ствующих режимам работы агрегата с различной производи­ тельностью: номинальной, выше и ниже номинальной. Сред­

нее значение квадрата ошибки csl находилось в диапазоне

Рис. 4.8

208

2247 14

ГО

о

о

Соседние файлы в папке книги из ГПНТБ