Файл: Живоглядов, В. П. Адаптация в автоматизированных системах управления технологическими процессами.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 62

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Все автоматические системы в зависимости от объема ап­ риорной информации об объекте и цели управления могут быть разделены на два больших класса систем: с полной ин­ формацией и с неполной [В .15—В .17].

Теория оптимальных систем управления до недавнего времени развивалась в основном как теория систем с полной информацией [В. 15, В. 18—В. 20]. Несколько лет назад чис­ ло публикаций, равно как и научных результатов в области управления с накоплением информации, было сравнительно невелико. Однако на практике в большинстве систем инфор­ мацию об объекте нельзя считать полной. В распоряжении разработчиков алгоритмов, как правило, имеются лишь не­ которые приближенные математические модели процессов, зависящие от набора параметров, которые определяют по результатам промышленного эксперимента. Ошибки экспери­ мента, наличие различного рода случайных помех, дрейф ха­ рактеристик объектов в процессе эксплуатации существенно усложняют задачу автоматического управления.

Класс систем с неполной информацией довольно широк, причем в нем можно выделить системы без накопления ин­ формации, с пассивным и активным накоплением ее. В каче­ стве математического аппарата исследования нашли приме­ нение корреляционные теории [В. 21, В. 22], теория управ­ ляемых случайных процессов [В. 23—В. 24] и условных мар­ ковских процессов [В. 25—В. 26], стохастическая аппрокси­ мация [В. 17, В. 27, В. 28], теория статистических решений [В. 22] и теория дуального управления [В. 15]. В данной ра­ боте рассматривается класс систем с неполной информацией и накоплением ее в процессе работы, т. е. класс адаптивных систем. Основное внимание уделяется малоизученным зада­ чам стохастического и адаптивного управления объектами с чистым запаздыванием и распределенными параметрами. Теория систем с распределенными параметрами интенсивно стала развиваться в последние годы. В литературе можно найти сведения об исследовании оптимального детерминиро­ ванного управления [В. 19, В. 29—В. 31], устойчивости [В. 32], идентификации [В. 33—В. 37], вероятностном анализе [В. 38—

В. 40] распределенных

систем, фильтрации

[В. 41—В. 44],

прогнозе случайных

полей [В. 45, В. 46]

и управлении

[В. 46—В. 47].

 

 

Ниже приведены математические модели объектов с за­ паздыванием и распределенными параметрами, дана харак­ теристика их динамических свойств и сформулированы в об­

10


щих чертах задачи синтеза алгоритмов управления в усло­ виях неполной информации о параметрах, возмущениях и со­ стоянии управляемых объектов.

В. 2. МАТЕМАТИЧЕСКИЕ МОДЕЛИ ОБЪЕКТОВ С ЗАПАЗДЫВАНИЕМ И РАСПРЕДЕЛЕННЫМИ ПАРАМЕТРАМИ

Известно, что уровни математического описания производ­ ственных объектов в зависимости от назначения математи­ ческой модели могут быть различными. Так, для выбора на­ строек линейных регуляторов в системах стабилизации зача­ стую можно ограничиться приближенной моделью, получен­ ной из условий теплового и материального балансов, либо путем аппроксимации экспериментально снятых кривых раз­ гона. При оптимизации технологических режимов, разработке систем оптимального или близкого к оптимальному управ­ ления, систем управления с УВМ, решении задач прогноза вы­ хода возникает необходимость получения более точных ма­ тематических моделей, учитывающих химизм протекающих: реакций, распределенность процессов в пространстве, случай­ ный характер возмущений и дрейфа параметров и т. д. При других подходах достаточно статистического описания объек­ тов. Применение адаптивных алгоритмов снижает требования к точности моделей объектов управления.

Многие промышленные процессы в цементном производст­ ве, металлургии, обогащении полезных ископаемых, нефте­ химической промышленности характеризуются распределен­ ностью процессов в пространстве. К объектам такого типа относятся, например, вращающиеся обжиговые печи, бара­ банные шаровые мельницы, трубчатые реакторы, классифи­ каторы, флотомашины и т. д. Качественные показатели их работы, как правило, определяются не только значениями технологических переменных (температуры, давления, кон­ центрации веществ) в одной точке, например, на выходе объекта, но и их распределениями по длине, высоте агрега­ та или какой-либо другой координате.

С распределенностью объектов связан эффект транспорт­ ного запаздывания в них, поскольку необходимо определен­ ное время на перемещение обрабатываемого материала из одной пространственной точки технологического аппарата в., другую. Математическая модель такого рода объектов может быть задана в виде уравнений в частных производных или в «частных» разностях, в виде системы иррациональных пере­

11

даточных функций, либо интегральными соотношениями. Сто­ хастические объекты адекватно и наиболее полно могут быть описаны путем задания условных вероятностных распреде­ лений (или условных плотностей вероятности) выходных сиг­ налов объектов при фиксированных входных контролируе­ мых воздействиях и задания вероятностных характеристик возмущений.

В. 2. 1. Линейные непрерывные модели

Простейшим Примером распределенной системы является

•одномерный объект с чистым запаздыванием, описываемый уравнением в частных производных первого порядка

M + V ? x = ° {t>0’ 1*>Х>У-

(В. 1)

 

9(0,г)=/С[а(/)+|*°(/)]

снекоторыми начальными условиями.

Здесь q— q(x,

i) — функция состояния объекта,

зависящая

от пространственной координаты х и времени t.

Параметр и

характеризует

скорость перемещения среды

в аппарате,

u(t), \i(t) — управляющее и возмущающее воздействия, К — коэффициент усиления. Решение уравнения (ВЛ) при

q(x,t) = Ku°(t ^ )= K u [ t - z x },~x = Л

(В. 2)

и передаточная функция имеет вид

х

(В. 3)

В частном случае, когда х=1н и считая q(lH^ ) —qo(0 выхо­ дом объекта, имеем

(В. 4)

тде хн— время чистого или транспортного запаздывания.

-12


Примерами таких объектов могут служить ленточныетранспортеры, ленты с магнитной записью, когда запись и счи­ тывание сигналов производятся в различных пространствен­ ных точках [В. 47, В. 48], клеть прокатнего стана, рассмат­ риваемая как объект регулирования толщины проката [В. 49] . В последнем случае управляющим воздействием служит раст­ вор валков стана. Запаздывание хн создается за счет того, что толщиномер по техническим условиям устанавливается на некотором удалении от оси валков. К устройствам, обладаю­ щим чистым запаздыванием, относятся также гидравлические или пневматические каналы связи между различными час­ тями автоматических систем.

В книге [В. 50] для описания трубчатого химического ре­ актора непрерывного действия приведено следующее урав­ нение:

+ v -^ -= К < 72, (/«>/*>0, / > 0),

(В. 5)

<7i(0,0 = 0.

В реакторе протекает необратимая реакция qz->qi первого-

порядка, где k — константа скорости реакции, сек~х; qx и q2 — концентрации веществ.

Уравнением первого порядка описывается также [В. 19] процесс нагрева тонкого слоя, движущегося в положительном, направлении оси х материала

 

(В.6)

<7(0,0= ц°(*).

/>0),

где q = q (х, t) — температура материала; и(х, t) — температура преющей среды.;-

р,°(t) — температура материала на входе в агрегат. Нагрев осуществляется за счет теплопередачи между грею­ щей средой и материалом. Коэффициент а(х, t) определяет­ ся теплофизическими параметрами материала.

Динамика широкого класса промышленных объектов мо­ жет быть описана линейными дифференциальными уравнения­ ми в частных производных второго порядка следующего вида:

d2q

d2q .

d2q , d

q . d q

.

s i r

+ “■- w + “ ■m r

+ а ‘ Гх + a4 t + a 'q~

(B. 7)

 

 

 

 

 

с некоторыми начальными и граничными условиями,

причем

4aia2— а2з =

6<!0. При 6=

0 уравнение (В. 7) является урав­

нением параболического

типа,

а

при 6 < 0 — гиперболиче­

ского.

 

 

 

 

 

Рассмотрим некоторые частные случаи.

 

Характерная особенность многих

процессов обогащения

полезных ископаемых — наличие в технологической схеме ап­ паратов, представляющих собой гидравлические емкости (флотомашины, классификаторы) [В. 51, В. 52]. Через них непрерывно протекает пульпа, в которой происходят физико­ химические превращения частиц обогащаемого минерально­ го сырья и реагентов.

Динамика участков технологической схемы, включающих указанные аппараты, определяется главным образом макро­ процессами — переносом и перемешиванием компонент пуль­ пы и переходом их из одной фазы в другую. В качестве ос­ новной математической модели, отражающей характерные ди­

намические особенности таких участков,

в работе [В. 51]

принято уравнение турбулентной диффузии

в движущейся

среде

 

 

° - 0 —

“ | г = 0 ’

(В. 8)

где q = q(x, t) — объемная концентрация

вещества (напол­

нение);

 

 

v — средняя скорость потока;

 

D — коэффициент турбулентной диффузии.

Граничные условия определяются потоком и характеристи­ ками загружаемого в аппарат материала и потоком разгруз­

ки. Они имеют, например,

такой вид:

 

 

4 ! -

= м ° ( о ,

(в. 9)

дд_

=bl [ul(t)—q{lHJ)]

(В. 10)

дх

Х = 1 ц

 

 

 

 

 

Воздействия u°(t) и u!(t) могут быть управляющими или воз­ мущающими.

14


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

Мокрое измельчение является непрерывным технологиче­ ским процессом. Исходный материал, проходя через мельни­ цу, измельчается до требуемой тонкости. Активное перемеши­ вание, возникающее при вращении мельницы, способствует образованию потоков материала из сечений с большей кон­ центрацией в сечения с меньшей концентрацией. Отсутствие специальных механизмов для транспортировки обусловлива­ ет возможность перемещения материала за счет разности кон­ центраций в соседних сечениях и вытеснения материала. Для описания динамических свойств мельниц можно применить уравнение (В. 8), понимая под D коэффициент переме­ шивания.

Некоторые авторы (например, [В. 53]) указывают на воз­ можность использования в отдельных случаях более простой модели, а именно:

§ = Д

(/«> *> 0 . t > 0)

(В. 11)

с граничными условиями второго и третьего рода.

 

Коэффициенты D,

v, k зависят от ряда факторов и в про­

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

Управление такими объектами, как правило, производится по границам либо является параметрическим. Например, уп­ равление может осуществляться путем изменения скорости прохождения материала через агрегат.

Звено, описываемое уравнением (В. 11), иногда называют звеном с «распределенным» запаздыванием [В. 54] или полуинерционным [В. 47, В. 55]. Выражение (В. 11) определяет процесс нагрева массивных тел, например, слитков металла, и называется в этом случае уравнением теплопроводности, а коэффициент D — коэффициентом температуропроводности.

Выше мы рассмотрели сравнительно простые модели звеньев с распределенными параметрами. Зачастую для более точного описания промышленных объектов желательно ис­ пользовать не одно уравнение, а систему двух и более диф­ ференциальных уравнений в частных производных первого и второго порядков [В. 56].

15


В. 2. 2. Передаточные функции звеньев с распределенными параметрами

Особенность передаточных функций звеньев с распреде­ ленными параметрами состоит в том, что они включают транс­ цендентные и иррациональные по отношению к оператору р

выражения

[В. 47, В. 54].

Пусть объект описывается уравне­

нием

(В.

7)

в области

 

(х>0,^>-0) при

р(х, t ) = 0,6 = 0,

ai>0,

о6< Д

Начальные условия нулевые,

а граничные име­

ют вид

 

 

 

 

 

 

 

 

c04 (0,t)+ c^

= и°(0<°° ,

 

 

 

 

 

*=0

(В. 12)

 

 

 

q(x,t)<co

 

 

 

 

 

при х—»°о и /->оо.

Считая u°(t) — входным воздействием, a q(x, t) — выходным и применяя к уравнению (В. 7) преобразование Лапласа, по­ лучим следующую передаточную функцию:

ф, , _ д(х,р) _

ехр{—х(тяр +

/га.,+/т|Д + отг)

а°(р)

c0—ci(m3p + m i+ / m lp + m.,)

где

 

а3а4—2агав

ОТ‘

2аЧ

1

 

 

т2 =

а24— 4fltae

 

 

 

4a2i

 

т3= 2ai

/л4-

_ а 4_

2

(В 13)

(В- 14)

Рассмотрим некоторые частные случаи.

a^= —v,

a5—1,

^

Пусть в уравнении

(В. 7)

ax= D,

^а6=

^ а 2= а з = а 4= : |х = 0 и в условии

(В. 12)ci = 0, т.

е. зве­

но оййсывается уравнением

(В. 8).

 

протекающих в

Приходим к соотношению для описания

обогатительных аппаратах процессов диффузионного перено­

са вещества без массопередачи [В. 57] ^ - j6j.