Модальное управление температурным полем распределенного объекта Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

СПИСОК ЛИТЕРАТУРЫ

1. Филиппов М.М., Бабушкин Ю.В., Грибенюков А.И., Гин-сар В.Е. Система управления многозонной термической установкой для выращивания кристаллов по методу Бриджмена // Известия Томского политехнического университета. - 2010. -Т. 316. - №5. - C. 146-151.

2. Марков А.В. Выращивание монокристаллов арсенида галлия с высоким структурным совершенством методом вертикально направленной кристаллизации // Известия вузов. Материалы электронной техники. - 2006. - № 6. - С. 16-19.

3. Бутковский А.Г. Характеристики систем с распределенными параметрами. - М.: Наука, 1979. - 224 с.

4. Рапопорт Э.Я. Структурное моделирование объектов и систем управления с распределенными параметрами. - М.: Высшая школа, 2003. - 229 с.

5. Пантелеев А.В., Бортаковский А.С. Теория управления в примерах и задачах. - М.: Высшая школа, 2003. - 583 с.

Поступила 12.04.2011 г.

УДК 62-533.65+681.52

МОДАЛЬНОЕ УПРАВЛЕНИЕ ТЕМПЕРАТУРНЫМ ПОЛЕМ РАСПРЕДЕЛЕННОГО ОБЪЕКТА

М.М. Филиппов, А.И. Грибенюков, Ю.В. Бабушкин*, В.А. Кочегуров*

Институт мониторинга климатических и экологических систем СО РАН, г. Томск *Томский политехнический университет E-mail: imces@yandex.ru

Рассмотрена задача стабилизации температурного поля объекта с распределенными параметрами с помощью системы модального управления, включающей программируемый контроллер, который производит анализ гармоник и выработку управляющих воздействий. На примере длинного стержня, с встроенными в него нагревателями показано, что точность стабилизации температурного поля, достигаемая рассмотренным методом, соизмерима с точностью системы управления на основе ПИД-ре-гуляторов.

Ключевые слова:

Система с распределенными параметрами, модальное управление, стабилизация, температурное поле.

Key words:

Distributed parameters system, modal control, stabilization, thermal field.

В современной теории автоматического регулирования все большее внимание уделяется проблемам управления объектами с распределенными параметрами (ОРП). К настоящему времени получены теоретические результаты, позволяющие решать конкретные задачи оптимального управления в различных областях [1].

Интерес представляет использование разработанных теорий и для таких объектов, как многозонные термические установки для выращивания кристаллов [2], параметры которых распределены в пространстве. Однако, подбирать алгоритмы управления непосредственно на установке или ее модельном представлении не всегда возможно, поэтому в качестве модели объекта управления рассматривается упрощенный объект, например, стержень длиною Ь, характеризующийся одномерным распределением теплофизических свойств вдоль независимой пространственной переменной х и зависящим от времени распределением температуры Т(х,0 [3].

Целью работы является определение структуры и параметров регулятора для системы модального управления температурным полем стержня и сравнение возможностей стабилизации температурного поля с системой управления, включающей в контур управления независимые ПИД-регуляторы.

Модальное представление объекта управления в виде длинного стержня, математическая модель которого приведена в [3], записывается следующим образом

йТп (Я ’t] = -(аЯ + к )Тп (Рп, t) +

Ш

+4А (Я ,) + кТ°с + К( Я X (1)

Т(Яп,0) = Т0(Рп), п = 1,

_ I

где Т(я, 0 = |т(х, 0ф(я, х)г(х)оХ - конечное ин-

0

тегральное преобразование температуры; (р(л,х) -ядро конечного интегрального преобразования, которое определяется в результате решения соответствующей задачи Штурма-Лиувилля; /л - параметр преобразования; г(х) - весовая функция, зависящая от структуры математической модели объекта управления (для модели стержня г (х)=аг1); а, к, д„, Я(л„) - коэффициенты, учитывающие температуропроводность стержня, теплоотдачу с поверхности стержня в окружающую среду, тепло, выделяемое источниками тепловыделения и граничные условия; Т“, -„(л,,,0 - конечные интегральные преобразования температуры окружающей среды и управляющего воздействия.

Полагая, что температура стержня в начальный момент времени равна температуре окружающей среды, уравнения (1) можно записать в форме отклонений мод

Л ЛТп (Яп, 0 =

= -(ал1 + к)ЛТп (Яп, t) + Чу Лип (Яп,п =1 ОТ

где

__ I

АТп (Лп , t) = |(Т* (х) - Т(X t))(Рп ( Лп , х)г(х)^

0

I

АК(Лп, О = \(й(X) -и(х, t))(pn (Лп, X)г(X)Лх,

0

Т*(х) - заданное распределение температуры, ~(х) - распределение управляющих воздействий на нагреватели, необходимое для реализации Т*(х).

Модальное представление ОРП отличается от типовых моделей систем с сосредоточенными параметрами только бесконечным порядком указанной системы, учитывающим пространственную распределенность функции состояния. Распределенное управление Ли(х,/) при модальном описании объекта заменяется его эквивалентным представлением в виде бесконечного числа независимых сосредоточенных воздействий Л-„(д,,0, образующих Ли(х,() в форме ряда

от

Ли ( х г1) = £Лип (Лп, (Лп, х)■

п=1

Функционал качества, заданный как взвешенная с постоянными коэффициентами сумма интегральных квадратичных ошибок приближения управляемой функции состояния Т(х,/) на всем протяжении процесса управления к невозмущенному процессу Т*(х), а также энергетических затрат, обычно оцениваемых величиной интеграла от квадрата управляющего воздействия по области его определения, при модальном представлении объекта имеет вид [1]

(от ______ ^

Р1 ^(Лп, t) +

'

+Р2ЧІ 1Аи1(Лп, 0)

(2)

Щтт + [-(аЛ2„ + к)ЛТп (Лп , t) + ЧуЛип (Я , t)] +

дt дЛТп

+ р1ЛТп (Лп , О + р2ЧуЛип (Лп , О = 0 п = 1 от,

где щ - функция Беллмана.

Оптимальное управление находится из условия достижения максимума

-(ол1 + к )ЛТп (Лп, 0 +

+Чу Лип (Лп, )

дУп + дУп_

ді дАТп . , „„ чг. ,

п _ 'Іу пУГ“п 5

+ РіАТп (Лп , І) + Р2<йАК (Лп , О

дАип (Лп > О

= 0.

Откуда оптимальное управление и-й моды примет вид

АГ =-

К.

Р2 Чу

где

Кп =Р2

-(ар1 + к х + 4 ¡(арп + к )2 +"РL

Р2

- решение алгебраического уравнения Рикатти, тогда

Аип =-----------

Чу

-(ар + к) +

Р2

А Тп (Л , О-

где р1, р2 - весовые коэффициенты, определяющие вклад отклонений температур и энергетических затрат, выбор которых представляет собой самостоятельную проблему.

При возможности автономного регулирования всех мод ЛТ„(л„,0 с помощью несвязанных друг с другом мод управляющих воздействий Л-„(л„,0 решение задачи аналитического конструирования регуляторов при минимизации функционала распадается на независимые задачи минимизации всех его отдельных гармоник.

Для поиска оптимального управления по каждой моде используется метод динамического программирования. Согласно этому методу необходимо записать для каждой моды уравнение Беллмана:

Полученное управление реализуется в замкнутом контуре регулирования каждой из N мод выхода объекта с пропорциональным регулятором Ки и воздействует только на контролируемые N мод АТ„(л„,0, и=1,Ыуправляемой величины Т(х,/).

Согласно модальному описанию объекта системой независимых друг от друга дифференциальных уравнений пер-ого порядка относительно N учитываемых мод АТ„(л„,0, и^^ величины Т(х,і) с автономными управлениями Ам„(р„,/), и= 1 , N по модам распределенного воздействия и(х,і), система модального управления объектом с распределенным воздействием должна быть построена в виде совокупности N несвязанных контуров регулирования учитываемых мод.

Управление и(х,і) на входе объекта рассчитывается по его модам Ам„(р„,/) в форме суммы Nслага-емых ряда

N

и(х, Ґ) = й( х) -^Мп(Лп, №п(Рп >х)■

п=1

Значения АТ„(л„,^) вычисляются анализатором гармоник НАТ

I

АТп (л„, 0 = |(Т *(х) - Т (х 0Ж (л ,х)г( х)йх =

0

= Т (Лп Х - Т(Лп > (Х п =1 N■

Вариант реализации модального управления представлен на рисунке.

Рисунок. Система модального управления ОРП. Обозначения: НАТ - анализатор гармоник; БУ - блок умножения; БО - блок ограничения; Р - регуляторы мод

Рассмотрим применение системы модального управления для формирования температурного поля, используемого при выращивании кристаллов. По условиям технологического процесса необходимо реализовать распределение температуры Т*(х), состоящее из трех температурных зон

Т*(х) =

Т *, 0 < х < х:,

Т* + (Тг - Т*)(х- х{)/(х2 - х1),

х1 < х< х2,

х2 < х < Ь,

где Р, Рг - температуры низко- и высокотемпературной зон; X;, х2 - координаты начала переходной и высокотемпературной зон.

В табл. приведены результаты расчета заданного стационарного распределения температуры стержня, достигаемого с использованием системы модального управления, при различных весовых коэффициентах функционала (2) и системы управления на основе независимых ПИД-регуляторов.

Результаты расчетов показывают, что система модального управления обеспечивает точность стабилизации заданного распределения температуры (средняя ошибка ~0,07 К), соизмеримую с точностью системы управления на основе независимых ПИД-регуляторов (средняя ошибка ~0,08 К). Изменение весовых коэффициентов функционала (2) влияет на точность стабилизации температуры незначительно.

Таблица. Уставки регуляторов и показания температур при реализации заданного распределения температуры, К, р=1

Номер секции 1 2 3 4 5 6 7

71* 353,00 353,00 353,00 363,00 373,00 373,00 373,00

Р2=1 353,08 353,11 353,11 363,05 373,05 373,07 373,02

Р2=0,1 353,09 353,04 353,04 363,00 372,96 372,98 373,03

р2=0,01 352,98 352,99 353,11 363,04 372,83 373,08 372,86

ПИД 352,91 352,98 353,07 363,08 373,10 373,10 373,10

Выводы

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

2. Точность стабилизации температурного поля при соответствующем выборе числа гармоник и весовых коэффициентов квадратичного критерия качества достигает 0,05 %, что соизмеримо с точностью системы управления на основе ПИД-регуляторов.

СПИСОК ЛИТЕРАТУРЫ

1. Рапопорт Э.Я. Оптимальное управление системами с распределенными параметрами. - М.: Высшая школа, 2009. - 677 с.

2. Филиппов М.М., Бабушкин Ю.В., Грибенюков А.И., Гин-сар В.Е. Система управления многозонной термической установкой для выращивания кристаллов по методу Бриджмена // Известия Томского политехнического университета. - 2010. -Т. 316. - №5. - С. 146-151.

3. Филиппов М.М., Грибенюков А.И., Бабушкин Ю.В., Кочегу-ров В.А. Оптимальная стабилизация температурного поля распределенного объекта // Известия Томского политехнического университета. - 2011. - Т. 319. - № 4. - С. 21-26.

Поступила 20.04.2011 г.