Материал: Моделирование систем автоматического регулирования теплоэнергетических установок (120

РАБОТА 8. ПАРАМЕТРИЧЕСКАЯ ОПТИМИЗАЦИЯ САР

В работе проводятся параметрическая оптимизация и анализ влияния параметров регулятора на выбранные критерии качества системы стабилизации давления (см. работу № 6). Критериями качества были выбраны:

– квадратичный функционал

J t1 tp ρк2dt;

p 0

– максимальное отклонение относительного давления в камере

Δρк max в переходном процессе (рис. 3);

– время переходного процесса tр;

– комплексный аддитивный критерий (сумма частных критериев с весовыми коэффициентами).

Рис. 3. Переходный процесс САР

Оптимизируемыми параметрами системы являются коэффициент передачи исполнительного устройства kс и постоянная времени чувствительного элемента Т1.

16

Порядок выполнения работы

I. Подготовка схемы моделирования к оптимизации.

Схему моделирования исходной САР (см. работу № 6) необходимо дополнить блоками, формирующими критерии оптимальности. Удалите из схемы блок «Фазовый портрет». Перенесите в схемное окно блоки, необходимые для формирования критериев оптимальности: квадратичный функционал – из библиотеки «Динамические», максимальное значение – из библиотеки «Нелинейные». На выходах критериев установите блоки «В память» из библиотеки «Субструктуры». Проведите необходимые линии связи.

В библиотеках типовых блоков ПК «МВТУ 3.7» нет блока измерения времени переходного процесса. Создайте такой блок в виде субмодели. Из библиотеки «Субструктуры» перенесите блок «Субмодель» в поле Схемного окна и установите его около выходной координаты САР. Инициализируйте этот блок «мышью» и нажмите клавишу PgDn, при этом откроется Схемное окно субмодели. Сформируйте в окне субмодели схему измерения времени переходного процесса (рис. 4). Типовые блоки схемы расположены в следующих библиотеках: блок «Порт входа» (входной сигнал) и блок «В память»

– в библиотеке «Субструктуры», модуль – в библиотеке «Операции», ключ 3 – в библиотеке «Ключи», задержка на шаг – в библиотеке «Дискретные», часы (модельное время) – в библиотеке «Источники».

Рис. 4. Субмодель – время переходного процесса

Схема работает следующим образом. На средний входной порт ключа подается выходной сигнал САР. Если этот сигнал выходит за

17

пределы допустимой зоны, на выход ключа передается сигнал с третьего нижнего входного порта, т. е. текущее модельное время. Если выходной сигнал САР находится в пределах допустимой зоны, на выход ключа передается сигнал с первого верхнего входного порта, т. е. тот же сигнал, но задержанный на один шаг интегрирования. Двукратным щелчком «мыши» откройте Диалоговое окно ключа и установите ширину допустимой зоны нестабильности 0,05. Закройте окно субмодели, нажав на клавишу PgUp. На выходе блока измерения времени переходного процесса блок «В память» включать в схему не надо, так как он расположен внутри субмодели.

II. Установка параметров оптимизации.

Оптимизируемыми параметрами в рассматриваемой САР являются коэффициент передачи ИУ (интегрирующее звено) k и постоянная времени ЧУ (апериодическое звено первого порядка) Т. Оптимизируемые параметры САР задаются как глобальные параметры. Щелкнув «мышью» в поле Схемного окна, нажмите клавишу F8, при этом откроется окно «Глобальные параметры», реализуемое через окно редактора «Интерпретатор математических функций». Введите в первой строке k = 0,077, а во второй Т = 1. Закройте окно редактора, нажав на кнопку «ОK».

Двойным щелчком «мыши» откройте Диалоговое окно блока ИУ и вместо коэффициента передачи введите символ k. Аналогичным образом в блоке ЧЭ вместо постоянной времени введите символ Т.

Нажмите на кнопку «Оптимизация» и в открывшемся меню выберите опцию «Параметры», при этом откроется Диалоговое окно «Оптимизация» с инициализированной закладкой «Параметры». Нажмите «мышью» на специальную кнопку «<» в правой части Диалогового окна – откроется дополнительное окно «Глобальные параметры», в котором введены варьируемые при оптимизации параметры k и Т. При нажатой клавише Shift сделайте щелчки «мышью» на строках с параметрами k и Т. Нажмите на кнопку «ОK» – произойдет возврат в Диалоговое окно параметров оптимизации. Введите минимальные и максимальные значения варьируемых параметров: для k – 0,07 и 2, а для Т – 0,1 и 2.

III. Подготовка и проведение оптимизационного расчета.

В схеме моделирования установлены блоки для трех критериев оптимальности. Оптимизационные расчеты выполняют в следующей

18

последовательности: вначале проводится оптимизация по каждому из критериев в отдельности (частные критерии), а затем формируется комплексный критерий как сумма частных критериев с весовыми коэффициентами (аддитивный метод свертки частных критериев).

Для проведения оптимизационного расчета по первому частному критерию оптимальности – квадратичному функционалу – необходимо выключить из схемы моделирования блоки, которые формируют другие частные критерии. Это осуществляется инициализацией выключаемого блока однократным щелчком «мыши» и нажатием «мышью» на командную кнопку «Исключить блок», при этом фон блока становится черным. После исключения из схемы всех частных критериев (включая соответствующие блоки памяти), кроме блока «Квадратичный функционал», нажмите на кнопку «Оптимизация», выберите опцию «Параметры» и в открывшемся Диалоговом окне нажмите на закладку «Критерии», в результате чего откроется незаполненное Диалоговое окно для задания локальных критериев оптимальности. Нажмите на кнопку «<» в правой части окна, при этом откроется дополнительное окно «Глобальные переменные» со списком переменных Y1, Y2 и Y3 соответствующих блоков «В память». Выделите переменную выбранного критерия (например, Y1) следующим образом: переместите курсор на строку с именем переменной, нажмите клавишу Shift и сделайте щелчок «мышью». Закройте данное окно, нажав на кнопку «ОK» – произойдет возврат в Диалоговое окно «Критерии» с выбранной переменной в качестве критерия оптимальности. Задайте минимальное (0) и максимальное (0,005) значения этого критерия. Нажмите на кнопку «Метод» – в открывшемся окне будут представлены методы оптимизации и другая информация. Выберите метод оптимизации Поиск-2, который отличается одним из простейших алгоритмов прямого поиска минимума функционала, обладает высоким быстродействием и эффективностью. Задайте максимальное число циклов моделирования 200. Две последние диалоговые строки заполняются только при использовании стохастического метода оптимизации. Закройте диалоговое окно «Параметрическая оптимизация», нажав на кнопку «Да».

Нажмите на кнопку «Оптимизация» и выберите опцию «Расчет» – начнется расчет, результаты которого отображаются на экра-

19

не. В окне «Результаты оптимизации» выводятся значения оптимизируемых параметров и выбранного критерия оптимальности. В окне «График» строится соответствующий переходный процесс. Для удобства наблюдения переходных процессов щелкните «мышью» в поле окна «График» – это окно займет положение поверх других окон. На экране также появится окно «Возникшие ошибки» с предупреждениями, что заданная точность не обеспечивается. На эти сообщения можно не обращать внимания. По завершении процесса оптимизации на экране появится окно с вопросом: «Обновить параметры оптимизации?». Переместив экранные окна таким образом, чтобы были видны результаты оптимизации, которые заносятся в таблицу и отображаются в виде графиков, внесите в соответствующую строку таблицы оптимальные значения параметров и распечатайте график переходного процесса. Откажитесь от обновления параметров оптимизации, нажав кнопку «Нет» в соответствующем окне. Закройте окно «Возникшие ошибки». После этого можно перейти к оптимизационному расчету по другому критерию.

С помощью кнопок «Исключить блок» и «Включить блок» исключите из схемы моделирования блок «Квадратичный функционал» и соответствующий ему блок «В память» и включите блок «Максимальное значение» с блоком «В память». Нажмите на кнопку «Оптимизация» и выберите опцию «Параметры» – откроется Диалоговое окно параметров оптимизации. Нажмите на закладку «Критерии», в открывшемся окне нажмите на кнопку «>» – строка критерия оптимальности станет пустой. Нажмите на кнопку «<» и в соответствующем окне «Глобальные переменные» выберите в списке переменных ту, которая соответствует блоку памяти, например, Y2. Выбор осуществляется щелчком «мыши» на нужной строке списка. Нажмите на кнопку «Да», при этом окно «Глобальные переменные» закроется и произойдет возврат в окно «Критерии» с заполненной строкой нового критерия оптимальности – максимального значения регулируемого параметра. Введите максимальное значение критерия, равное 0,3, и закройте окно «Параметры оптимизации», нажав на кнопку «Да». Нажмите на кнопку «Оптимизация» и, выбрав опцию «Расчет», начните оптимизационный расчет по этому критерию. По окончании оптимизации заполните соответ-

20