Выбор настроек и анализ замкнутой системы регулирования с ПИ-регулятором

Размещено на http: //www. allbest. ru/

1. Выбор настроек и анализ замкнутой системы регулирования с ПИ-регулятором

simulink автоматический регулятор математический

1. Цель работы

Целью работы является моделирование на ПЭВМ с помощью инструментария Simulink замкнутой системы автоматического регулирования (САР) по заданному математическому описанию и подбор настроек ПИ-регулятора.

2. Краткая теория

2.1. Замечания по структуре САР

Замкнутая CAP состоит из объекта регулирования и регулятора.

Объект может иметь различные передаточные функции по возмущающему каналу и по регулирующему каналу

Wвк Wрк.

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

Если возмущающий и регулирующий каналы совпадают

Wвк = Wрк = Wо, то схема CAP упрощается, рис. 2.

Такая схема называется схемой стабилизирующего регулирования. Во многих случаях объект имеет транспортное запаздывание по каналу регулирования. В этом случае для моделирования объекта необходимо дополнительно использовать блок постоянного запаздывания (БПЗ), рис. З.

При определении настроек регулятора постоянное запаздывание по возмущающему каналу не имеет значения и его можно принять равным нулю.

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

4. Подбор настроек ПИ-регулятора

4.1. Моделирование ПИ-регулятора.

Уравнение ПИ-регулятора записывают следующим образом:

(8)

где

yрег - выходная величина регулятора;

x = ( x - x3 ) -

входная величина регулятора, - рассогласование, - разница между текущим значением регулируемой величины x и заданием x3;

S1, S0 - настройки регулятора соответственно пропорциональная и интегральная части.

На рис. 1 показан ПИ-регулятор, собранный на элементах, взятых из библиотек пакета Simulink. Сумматор Sum1 является компаратором, т. е. устройством, сравнивающим текущую величину x(t) с заданием x3. Следует иметь в виду, что знак суммирования величины x3 должен быть противоположен знаку суммирования x(t).



Рис. 1

Пример 3.

Пусть регулятору следует поддерживать величину + 76.4 оС. В этом случае величина задания будет х3 = 76.4, поданное на вход сумматора со знаком " - ". На рис. 1 величина задания, реализованная на элементе

Константа (Constant), x3 = 100, поэтому следует изменить ее на величину x3 = 76.4. Тогда действительно, если текущая величина x ( t) = 76.4, т. е. равна величине задания, то на выходе сумматора Sum1 рассогласование x = 0. Если x(t) > 76.4, то x(t) > 0, если x(t) < 76.4, то x(t) < 0.

Пропорциональное и интегральное воздействия регулятора складываются на Cумматоре2 (Sum2), на выходе которого получается выходная величина регулятора y рег.

Все коэффициенты в схеме на рис. 1, кроме S1 и S0, постоянные и равны 1. Если S1 = 0, то получается

И-регулятор, если, наоборот, S0 = 0, - П-регулятор.

Задание. Установить (через рабочую область MatLab) S1 = 1 и S0 = 0.1. параметры скачка элемента Step установить следующие: время скачка - 3 с, начальное значение 76.4, конечное значение 60.0.

Запустить модель и наблюдать на осциллографе переменную yрег(t). После окончания моделирования, поменять произвольным образом величины настроек S1 и S0, снова запустить модель и наблюдать y рег(t).

4.2. Включение регулятора в обратную связь объекта регулирования

Пусть к примеру необходимо стабилизировать

ПИ-регулятором объект с передаточной функцией

 (9)

Составляем структурную схему моделирования в среде Simulink этого объекта без запаздывания, рис. 1. Для этого следует передаточную функцию переписать в виде дифференциального уравнения, а затем составить структурную схему (работа 1 ).

Рис. 2

Так как объект имеет запаздывание, то следует на выходе модели объекта включить блок запаздывания из библиотеки пакета Simulink. Выход объекта x ( t) подключаем ко входу регулятора "переменная". Управляющее воздействие регулятора yрег, подаваемое на вход объекта, должно суммироваться с возмущающим воздействием u(t) и компенсировать его (знак минус на Sum4 по входу yрег ).

Таким образом,

во-первых, для yрег должен быть предусмотрен отдельный вход с коэффициентом, равным 1;

во-вторых, регулятор должен быть включен в отрицательную обратную связь. Для примера придадим возмущению u(t) знак плюс. Если регулятор включен в отрицательную обратную связь, то возмущение, прошедшее через замкнутую систему "объект - регулятор", вернется в виде компенсирующего воздействия yрег со знаком " - " (минус). Если же получится знак " + " (плюс), то значит, регулятор включен в положительную обратную связь и замкнутая схема регулирования будет неустойчивой.

Если объект регулирования линейный, то полученные настройки регулятора и переходный процесс не зависят от величины задания. Поэтому для простоты принимают величину задания x3 = 0. Если x3 = 0, то вход для него не нужен и, в принципе, не нужен и сумматор Sum1.

4.3. Критерии качества регулирования

Зададимся целью подобрать такие настройки S1 и S0, чтобы переходный процесс в замкнутой САР был оптимальным.

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

Также надо точно оговорить стандартный вид возмущающего воздействия u(t), чтобы была возможность сравнивать переходные процессы. Скачкообразное возмущение (функция Хэвиссайда) u(t)=const является наиболее неблагоприятным, и вследствие, этого, принято стандартным для изучения САР, рис. 3 Оно реализуется элементом Step раздела Source.

Рис. 3

Когда возмущение является более сложной функцией времени t, то необходимо специальное устройство генерирующее u (t). Элементы-генераторы функций независимой переменной t легко переносятся в блок-схему модели из библиотеки ИСТОЧНИКИ (Source) пакета Simulink.

Критерии оптимальности J переходного процесса x(t) могут быть различными:

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

При аналитическом расчете настроек ПИ-регулятора по расширенной АЧХ переходной процесс должен отвечать двум критериям качества:

1. Степень затухания равна расчетной, например ш = 0.75 (ей cответствует степень колебательности m = 0.221). Величина ш показывает на интенсивность затухания и определяется как

при этом A i+2 > A i+3,

(Для ш = 0.75, Ai= 4*A i+2)

2. Выполняется интегральный критерий (10), т.е. минимум амплитуды переходного процесса и минимум времени регулирования.

На рис. показана плоскость параметров настройки с линией равного затухания, а на рис. переходные процессы, соответствующие различным парам настроек регулятора S1, S0 на линии равного затухания.

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

Таким образом, можно подобрать бесчисленное множество пар настроек S1 и S0, которые будут обеспечивать заданную степень затухания; однако необходимо найти такие, которые давали бы одновременно минимум интегрального критерия.

Рис. 4

Рис. 5



Рис. 6

Рис. 7

Рис. 8

4.4. Определение оптимальных настроек

Пусть необходимо подобрать оптимальные настройки для объекта с передаточной функцией (9). Степень затухания, к примеру, ш= 0.75.

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

Проверяется правильность реализации в среде Simulink объекта без запаздывания с отключенным регулятором. Для этого на вход u(t) подают скачкообразное возмущение, рис.. Начальные условия должны быть нулевыми.

1. Если кривая разгона совпадает с предполагаемой (время переходного процесса, его форма и коэффициент усиления в статике), то настраивают и подключают БПЗ. По осциллографу (Scope) наблюдают одновременно переходный процесс до БПЗ и после него. Для этого надо использовать перед входом осциллографа элемент МUХ. Если видно, что время запаздывания равно или больше времени переходного процесса в объекте, то можно не продолжать отыскание настроек, так как ПИ-регулятор дает неустойчивую замкнутую систему регулирования почти при всех настройках. Проверьте еще раз правильность полученного постоянного запаздывания.

2. Если БПЗ работает правильно, то подключают

3. ПИ-регулятор. Выход х(t) до БПЗ уже не нужен, и его можно отключить от осциллографа. Ставят настройки в нулевые состояния, т.е. S1 = 0 и S0 = 0. Это равносильно разрыву цепи обратной связи. Потом проверяют, получается ли первоначальная кривая разгона. Если кривая разгона объекта не получается, то необходимо проверить коэффициенты S1 и S0 в схеме регулятора, а также обязательно начальные условия на интеграторах, - они должны быть нулевыми.

4. Если при S1 = 0 и S0 = 0 кривая разгона объекта получилась правильно (рис.), то приступают к подбору настроек регулятора.

Simulink - модель должна находиться в режиме "исходное положение". Затем несколько увеличивают настройку S1. После чего переводят Simulink - модель в режим "решение" ("интегрирование") и наблюдают, как изменился переходный процесс. Затем снова переводят в режим "исходное положение", изменяют коэффициент, пускают машину в режим "интегрирование" и наблюдают новый переходный процесс и т. д.

При увеличении настройки S1 уменьшается амплитуда, кривые 1 и 2 на рис. При дальнейшем увеличении переходный процесс станет колебательным кривые 4 и 5. Для удобства наблюдения можно увеличить вертикальное усиление осциллографа и постараться подобрать такую настройку S1 , чтобы получить степень затухания, близкую к заданной ш= 0.75, кривая 1 на рис. Для того, чтобы оценить степень колебательности, следует мысленно провести линию установившегося состояния, рис. проверить величины амплитуд переходного процесса. Каждая i+2 амплитуда должна быть приблизительно в 4 раза меньше i-й для ш= 0,75. При дальнейшем увеличении настройки S1 степень затухания уменьшается, кривая 2 на рис..

5. Если, приблизительно, настройкой S1 П-регулятора достигнута заданная степень затухания, то следует перейти к подбору настройки S0, не изменяя S1.

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

При дальнейшем увеличении S0 можно получить переходный процесс, отвечающий ш = 0.75 и условию минимума интегрального критерия, кривая 2 на рис. При этом, может быть, понадобится немного изменить настройку S0.

Таким образом, сначала мы двигались в плоскости параметров настроек по оси S1 , а затем поднимались вверх, увеличивая S0, (рис., отмечено пунктиром).

Следовательно, вручную, мы осуществляли процедуру оптимизации типа " покоординатного спуска ".

6. Если переходный процесс получился оптимальным, то следует записать полученные коэффициенты, т.е. полученные настройки S0 и S1.

5. Интерпретация полученных значений настроек

Пусть известно, что у моделируемого объекта были следующие значения размерностей переменных:

И пусть также полученные в среде Simulink настройки имеют следующие значения: s1 = 2.6 и s0 = 0.85. Каковы - же их размерности ?! Зная размерности входа, выхода и размерность независимой переменной t (14-16), легко получаем

где ВРЕМЯ - это или секунды [c], или минуты [мин], или [часы] , в зависимости от того, какая размерность постоянной времени Т перед первой производной имела место в математической модели объекта. При этом коэффициент при x(t) должен быть единицей, например, для уравнения (3) величена Т равна

6. Пересчет настроек из машинных в реальные (для справки)

Предположим, что моделирование проводилось после предварительного масштабирования дифференциального уравнения, описывающего объект регулирования. Предположим также, что были при этом применены следующие масштабы:

где [с] реальное время ;

[см] модельное время в Simulink-модели.

Маш.ед. - модельные безразмерные единицы измерения величин u(t), и(t) или x(t) в среде Simulink, которая наблюдается на осях осциллографа (Scope).

Теперь требуется пересчитать полученные модельные настройки s0 и s1 в реальные S0 и S1.

Запишем уравнение регулятора

Так как выходная величина регулятора yрег является входным воздействием для объекта, таким же как и u(t), то можно записать

Для реальных переменных также можно записать



Переменные Simulink-модели x, y, t м связаны с реальными переменными G, и, t через масштабы

Подставим в (21) машинные переменные, выраженные через реальные и масштабы

(23)

Сгруппируем масштабы с модельными коэффициентами s0 и s1

(24)

Сравнивая теперь выражения (22) и (24), можно записать

Размещено на http: //www. allbest. ru/

После подстановки численных значений нашего примера получим



Степень обратной связи теперь легко определить

В нашем примере

Полезно бывает подсчитать также время изодрома

7. Оформление отчета

Отчет по лабораторной работе выполняется в виде связного (читаемого) текста в файле формата Microsoft Word (шрифт основного текста Times New Roman, 12 пунктов, через 1,5 интервала, выравнивание по ширине). Он должен включать

· название предмета, номер и название лабораторной работы

· фамилию и инициалы авторов, номер группы

· фамилию и инициалы преподавателя

· номер варианта

· краткое описание исследуемой системы

· результаты выполнения всех пунктов инструкции: результаты вычислений, графики, ответы на вопросы.

При составлении отчета рекомендуется копировать необходимую информацию через буфер обмена из рабочего окна среды визуального моделирования (Scilab, Matlab и др.). Для этих данных используйте шрифт Courier New, в котором ширина всех символов одинакова.

8. Задание

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

Время запаздывания взять из табл. 1.

Таблица 1

№ вар.	1	2	3	4	5	6	7	8	9	10	11	12
фр	0.4 мин	5 мин	0.16 с	11.8 с	2.6 с	1.,2 мин	42 мин	3.1 мин	4.2 мин	1.0 мин	0.003 с	0.65 с

Литература

1. Бесекерский, В.А. Теория систем автоматического управления : [Учебное издание] / В.А.Бесекерский, Е.П.Попов.-- 4-е изд.,перераб.и доп. -- СПб. : Профессия, 2004.-- 752с. :

2. Ким, Д.П. Сборник задач по теории автоматического управления. Линейные системы : учеб. пособие для вузов / Д. П. Ким, Н. Д. Дмитриева.-- М. : ФИЗМАТЛИТ, 2007.-- 168 с.

3. Андриевский, Б.Р. Элементы математического моделирования в программных средах МАТLАВ 5 и Scilab / Б.Р.Андриевский,А.Л.Фрадков.-- СПб. : Наука, 2001.-- 286с.

4. Цыпкин Я.З. Основы теории автоматических систем. - Наука, 1977.-559 с.

5. Пупков, К.А. Методы классической и современной теории автоматического управления. Т.3, Синтез регуляторов систем автоматического управления/К.А.Пупков [и др.];под ред.К.А.Пупкова,Н.Д.Егупова : учебник для вузов:в 5 т. -- 2-е изд.,перераб.и доп. -- М. : МГТУ им.Баумана, 2004.-- 616с.

Размещено на Аllbest.ru