№
п.п
|
=20
сек
|
|
Т1,
сек
|
Т2,
сек
|
Т3,
сек
|
кр
|
кр/Ти
|
6
|
20
|
200
|
0
|
0,51
|
0,028
|
1. Теоретическая
часть
.1 Методы моделирования
.1.1 Моделирование систем
автоматического регулирования
Характеристики промышленных объектов сведены к
типовым характеристикам, поэтому многочисленные законы функционирования
регуляторов, работающих с промышленными объектами тоже можно свести к типовым
законам, а именно:
Пропорциональный закон
Интегральный закон
Пропорционально - интегральный закон
Пропорционально - дифференциальный закон
Пропорционально - интегрально - дифференциальный
закон
Пропорциональный регулятор
Описывается в динамике следующим уравнением:
(t) = K1*Dy(t)
Любой регулятор по отношению к объекту,
являющейся неизменной частью системы, можно считать последовательным
корректирующим звеном, поэтому для определения влияния регулятора на систему
найдем его передаточную функцию и определим, к каким типовым динамическим
звеньям его нужно отнести:
(p) = x(p)/ Dy(p) = K1
Пропорциональный регулятор относится к
безинерционным звеньям и обеспечивает хорошие динамические свойства системы.
Интегральный регулятор
Описывается в динамике следующим уравнением:
¶x(t)/¶t
= K1*Dy(t)(p)
= x(p)/ Dy(p)
= K1/p
Интегральный регулятор относится к интегрирующим
звеньям и обеспечивает хорошие статические свойства системы, но ухудшает при
этом динамику систему.
Пропорционально - интегральный регулятор
Описывается в динамике следующим уравнением:
x(t) = K1*Dy(t)
+ K0* (p)
= K1 + K0/p
Пропорционально - интегральный регулятор можно
представить как последовательное соединение безинерционного, интегрирующего и
пропорционально - дифференцирующего звеньев. Улучшает и статику и динамику
системы.
Пропорционально - дифференциальный регулятор
Описывается в динамике следующим уравнением:
x(t) = K1*Dy(t)
+ K2* ¶Dy(t)/¶t(p)
= K1 + K2*p
Пропорционально - дифференциальный регулятор
можно представить как последовательное соединение безинерционного и
пропорционально - дифференцирующего звеньев. Улучшает динамику системы, внося
при этом статическую ошибку в систему.
Пропорционально - интегрально - дифференциальный
регулятор
Описывается в динамике следующим уравнением:
x(t) = K1*Dy(t)
+ K2* +
K3* ¶Dy(t)/¶t
W(p) = K1 + K2/p + K3*p
Пропорционально - интегрально - дифференциальный
регулятор можно представить как последовательное соединение безинерционного
звена, интегрирующего звена и пропорционально - дифференцирующего звена второго
порядка. Улучшает как статику, так и динамику системы.
1.1.2 Моделирование объекта
управления
Под моделированием какого - либо объекта,
понимается воспроизведение и исследование другого объекта подобного оригиналу в
форме, удобной для исследования, и перенос полученных результатов на
моделируемый объект. При этом объекты считаются подобными, если характеристики
процессов, протекающих в каком - либо из них отличаются от соответствующих
характеристик другого объекта вполне определенными и постоянными коэффициентами
в течение одного процесса.
Существуют различные методы моделирования:
геометрическое моделирование
физическое моделирование
моделирование методом прямых аналогий
математическое моделирование
полунатурное моделирование
Каждый из методов имеет свои плюсы и минусы.
Применение конкретного метода обусловлено исследуемой системой и условиями ее
работы.
В данной курсовой работе рассматривается только
математическое моделирование. При математическом моделировании в качестве
объекта моделирования выступают исходные уравнения, представляющие
математическую модель объекта, а в качестве модели системы выступают те
процессы, которые воспроизводятся на ЭВМ и решают по заданным уравнениям
поведение отдельных компонент и их взаимосвязь в системе.
Моделирование на вычислительных машинах часто
называют аналоговым и цифровым моделированием. Имитационное моделирование
рассматривают также как управляемый эксперимент, проводимый не на реальном
объекте, а на модели подобной оригиналу. В этом случае оценивается поведение
отдельных компонент, входящих в систему в условиях направленного эксперимента:
заданных входных координатах возмущений, параметрических и внешних случайных
воздействий.
Моделирование звеньев
Инерционное звено первого порядка.
В динамике описывается дифференциальным
уравнением первого порядка, которое может быть приведено к виду:
где T - постоянная времени звена; -
статический коэффициент передачи звена;
В операторной форме уравнение имеет
вид:
Т py(p) + y(p) = kx(p)
А передаточная функция находится
как:
Переходная характеристика в операторской форме:
Значение выходного сигнала
инерционного звена 1-го порядка находится по формуле Yn+1 = Yn + Y¢HT, где НТ -
приращение времени и Y¢ = .
у
¢HT
yn+1t
Рис. 2
Инерциальное звено второго порядка.
Колебательное звено.
Дифференциальное уравнение
инерционного звена второго порядка:
в операторной форме:
Т22p2yвых(p) + T1pyвых(p) + yвых(p)
= kxвх(p)
Передаточная функция:
Звено запаздывания (безинерционное).
Передаточная функция:
Для моделирования запаздывания, как
правило, создается вспомогательный массив длины N (N=t/h, h - шаг моделирования), в
который в соответствии с заданным временем запаздывания t осуществляется выборка и
загрузка координат входящих переменных на время запаздывания.
Выходным значением звена является
последний элемент очереди, при поступлении на вход нового значения все элементы
очереди сдвигаются вправо, последний теряется, а в начало записывается новый.
.1.3 Описание метода Рунге-Кутта
В вычислительной практике наиболее
часто используется метод Рунге Кутта четвертого порядка, имеющий погрешность R
~ (h5),. Метод Рунге-Кутта для дифференциального уравнения 2-го порядка вида:
''==F(x,y,y')
В этом методе величины yi+1
вычисляются по следующим формулам:
yi+1 = yi + Δyi Δyi=h(k1
+ 2k2 + 2k3 + k4)/6, i = 0, 1, ... k1 = f(xi, yi),
k2 = f(xi+h/2, yi+hk1/2), k3 = f(xi+h/2, yi+hk2/2), k4 = f(xi+h, yi+hk3).
1.1.4 Параметрическая оптимизация
Интегральные показатели качества служат для
анализа качества процесса регулирования. Оценка по интегральным критериям осуществляется
следующим образом:
IAE
При анализе и синтезе систем
используют обобщенные критерии с учетом сигналов ошибки и их производных,
взятых с весовыми коэффициентами. При перемещении симплекса в процедуре
оптимизации при поиске локального экстремума осуществляется оценка близости
нахождения к локальному экстремуму. Данная ситуация оценивается путем анализа
происхождения исследуемых точек и закрутки симплекса. Процедура закрутки
симплекса считается итерацией и после ее возникновения осуществляется изменение
размеров симплекса (уменьшения) после чего процедура повторяется.
1.1.5 Оценка динамических свойств
системы
Качество САР
Устойчивость является необходимым, но не
достаточным показателем САР. При исследовании систем автоматического
регулирования приходится решать задачу обеспечения требуемых показателей
качества переходного процесса: быстродействия, колебательности,
перерегулирования, характеризующих точность и плавность протекания процесса.
Показатели качества принято определять по кривой
переходного процесса и называть прямыми. Кривая переходного процесса может быть
получена теоретически (как решение дифференциального уравнения системы, когда
правая часть уравнения [входной сигнал] единичная ступенька) или экспериментально.
Пусть кривая переходного процесса системы имеет
вид:
Рис. 3
Максимальное динамическое отклонение -
максимальная разность между заданными и действительными значениями регулируемой
величины в переходном режиме.
∆max дин = hуст
2. Максимальное перерегулирование - максимальное
отклонение переходной характеристики от установившегося значения переходной
величины, выраженное в относительных единицах. Обычно σmax
≤
20÷30%.
.Колебательность процесса:
(определяется как отношение разности
двух соседних амплитуд, направленных в одну сторону, к большей из них в
относительных единицах)
Для работоспособных систем ψ ≥ 75÷90%
4.Время регулирования - tрегул -
минимальное время от начала нанесения возмущения до момента, когда регулируемая
величина будет оставаться близкой к установившемуся значению с заданной
точностью; т.е.
|h(t) - hуст| ≤ ∆, где
∆ - постоянная величина,
значение которой нужно оговаривать (обычно ∆=2÷5% hуст).
.2 Описание пакета ИМОДС
Автоматическое проектирование систем
управления является в настоящее время весьма актуальной задачей. Оно
основывается на имитационном моделировании, которое стало возможным с
появлением ЭВМ. Имитационное моделирование является продолжением существующих
ранее физического и математического моделирования, и в большинстве случаев
вытесняет их.
Физическое моделирование требует
реальной установки, стоимость которой может быть высокой. Кроме того,
физическое моделирование некоторых процессов практически невозможно.
Математическое моделирование требует большого количества расчетов,
высококлассных специалистов в области математики. Не все существующие в природе
объекты
Можно достаточно полно описать
математически. Существует большой класс объектов, математическое описание
которых очень сложно.
Система имитационного моделирования
ИМОДС предназначена для решения задач, связанных с машинным моделированием,
построения и исследования систем управления.
.2.1 Назначение и функции
Система имитационного моделирования
ИМОДС предназначена для автоматизированного решения задач цифрового
моделирования, исследования и оптимизации сложных динамических систем.
Система моделирования выполняет
следующие основные функции:
автоматизированное построение (сборка)
имитационной модели (ИМ),
сохранение ИМ на диске,
параметрическая оптимизация ИМ,
получение динамических характеристик
ИМ,
графическое представление
результатов моделирования.
Система имитационного моделирования
ИМОДС реализована при помощи транслятора Turbo Pascal 5.0.
1.2.2 Система файлов
Система имитационного моделирования ИМОДС
включает в себя комплекс основных программ, ряд вспомогательных файлов,
библиотеку типовых элементов структуры (ТЭС) и библиотеку имитационных моделей
(ИМ), которую пользователь может пополнять.
.2.3 Структура меню пользователя
Выбор режима работы можно осуществить
посредством перемещения светового маркера при помощи клавиш управления курсором
и нажатия на клавишу "ENTER" для активизации выбранной команды.
Команды также могут быть активизированы посредством нажатия на функциональные
клавиши, связанные с нужными командами. Например, активизация табличного
редактора может быть осуществлена посредством нажатия на клавишу F4.
Основное меню содержит следующие команды:
"Таблица F4"- осуществляет запуск
табличного редактора для создания и корректировки имитационных моделей (ИМ);
"Оптимизация F6" - осуществляет
задание параметров оптимизации для параметрического синтеза ИМ;
"Время F7" - осуществляет задание
времени начала и конца переходного процесса, а также шаг моделирования;
"Опции" - осуществляет активизацию
подменю "Опции" , имеющее следующие команды:
"Запуск F9" - осуществляет процесс
моделирования работы ИМ;
"Возврат F5" - возврат к графическому
экрану с результатами предыдущего моделирования ИМ;
"Оптимизация F8" - запуск
параметрической оптимизации ИМ в работу;
"Выход Alt-X" - выход в среду MS DOS/.
1.2.4 Табличный редактор
Табличный редактор представляет собой средство
для редактирования и создания ИМ. При работе он использует табличную форму
представления данных. Перелистывание страниц осуществляется посредством нажатия
на клавиши "PgUp", "PgDn". Номер страницы отображается в
правом верхнем углу.
Нижняя строка содержит сведения о командах,
которые могут быть выполнены посредством нажатия на горячие клавиши. Такими
командами являются:
"Запись" (горячая клавиша F2) -
осуществляет сохранение на диске построенной ИМ. Нужно указывать полные путь и
имя файла для записи. Если путь не указан, то файл помещается в текущий
каталог;
"Чтение" (горячая клавиша F3) -
осуществляет считывание ранее созданной ИМ с диска для редактирования или
работы с ней;
"Новая" (горячие клавиши Alt+F3) -
осуществляет создание новой ИМ. Применяется для очистки таблицы при создании
новой ИМ.
Имя активного файла ИМ находится в левом нижнем
углу экрана.
Колонки таблицы имеют следующее информационное
содержание:
первая - входы типового элемента структуры (ТЭС)
соответственно структурной схеме;
вторая - коды ТЭС соответственно структурной
схеме и библиотеки ТЭС;
третья - выходы ТЭС соответственно структурной
схеме;
четвертая - имена выходов ТЭС соответственно
структурной схеме; имя выхода ставится только в случае необходимости
использования параметра с данным именем при графическом представлении
результатов моделирования. Максимальное число выходов, которым соответствует
имя, не должно превышать 16;
пятая - "параметры настройки" имеет
восемь подколонок. Здесь задаются параметры ТЭС. При перемещении светового
маркера по вертикали в "шапке" пятой колонки появляются названия
параметров настройки, соответствующие коду ТЭС, который проставлен во второй
колонке.
При редактировании таблицы используются
следующие клавиши:- Y - для удаления строк;- для вставки пустых строк;
Клавиши управления курсором - для перемещения светового
маркера и активизации различных полей таблицы;, PgUp - для постраничного
листания таблицы вперед и назад.
.2.5 Оптимизация
Оптимизация реализована методом
последовательного симплексного планирования (симплексным методом). Режим
оптимизации включает в себя два следующих подрежима:
задание начальных условий и границ оптимизации;
запуск оптимизации в работу.
При выборе в основном меню режима
"Оптимизация" на экране появляется диалоговое окно "Начальные
условия и границы оптимизации". При помощи этого окна можно модифицировать
следующие параметры:
код оптимизируемого ТЭС согласно библиотеке ТЭС
(Приложение 5) и структурной схеме;
порядковый номер строки, в которой находится ТЭС
в таблице конкретной ИМ;
граничные и начальные условия оптимизируемых
параметров.
Также указывается:
нужно или не нужно выводить промежуточные
значения экспериментов при каждом перемещении симплекса (результаты оптимизации
на момент осуществления каждой итерации выводятся во всех случаях);
нужно или не нужно выводить результаты оптимизации
на принтер.
Оптимизация ведется по интегральному
квадратичному критерию качества:
где х - ошибка регулирования.
Код ТЭС, реализующий этот критерий - 328. В
случае оптимизации этот ТЭС должен быть помещен в имитационную модель. Его
входом, как правило, является ошибка регулирования. Выход должен обязательно
иметь номер 100.
После указания начальных условий и границ
оптимизации производится запуск оптимизации в работу посредством активизации
команды "Оптимизация" из подменю "Опции".
.2.6 Система графического
представления результатов
Основная часть экрана используется для
построения графиков. Внизу экрана расположена строка - подсказка, содержащая
назначение функциональных клавиш. Рассмотрим команды, которые могут быть
активизированы посредством нажатия на них.
"F1 - помощь" - осуществляет вывод
краткой справки на экран.
"F3 - повтор" - при активизации этой
команды на экране появляется следующий по порядку график. Порядковый номер
графика находится под графиком. Максимальное количество графиков - 10. Просмотр
графика можно также осуществить путем нажатия клавиш "1, 2, 3…".
"F4 - окно" - выдает окно настройки
графической подсистемы, содержащее 3 поля. Первое и третье поля являются
информационными. На первом поле изображены допустимые размеры окон от 1 до 10.
На третьем находятся имена выходов ТЭС для просмотра на графике (они задавались
в четвертой колонке табличного редактора). Перемещение по именам осуществляется
посредством нажатия клавиши "+" (для перемещения вверх) и
"-" (для перемещения вниз). На среднем поле задаются параметры
графиков. В первой его колонке задается номер окна, соответствующее изображению
в левом поле. Во второй колонке задается аргумент функции. Выбирается один из
трех входов ТЭС, находящихся на правом поле. В оставшихся колонках задаются
непосредственно функции. В них также выбирается один из выходов ТЭС из правого
поля. Каждая система координат может содержать до трех функций одного
аргумента.
При работе с окном настройки используются
следующие клавиши:
клавиши управления курсором - используются для
перемещения светового маркера по среднему полю;
+, - - используются для выбора номеров окон и
имен выходов ТЭС; клавиша "+" служит для увеличения параметра, а
"-" - для уменьшения;- используется для выхода из окна настройки;-
используется для сохранения параметров настройки подсистемы графического
представления данных.
Продолжим рассмотрение команд, которые могут
быть активизированы посредством нажатия функциональных клавиш.
"F5 - очистить" - используется для
очистки экрана от изображенных на нем графиков.
"F7 - линия" - используется для
изменения толщины линии графика с тонкой на жирную и наоборот.
"F8 - сетка" - используется для
включения и выключения режима отображения масштабной сетки на экране.
"F9 - печать" - осуществляет формирование
графика в удобной для вывода на форматный лист бумаги форме. При помощи этой
команды можно указать надпись под графиком. Максимальная длина надписи - 50
символов. Распечатка производится путем посылки на принтер графической копии
экрана.
"F10 - меню" - осуществляет выход из
подсистемы графического представления данных в основное меню.
В главе приведено описание различных методов
моделирования систем автоматического управления, промышленных регуляторов,
объектов управления. Изложен метод Рунге-Кутта для решения дифференциальных
уравнений второго порядка, способы вычисления динамических показателей качества
для системы автоматического регулирования, а также представлены интегральные
критерии качества, применяемые для анализа и синтеза систем при имитационном
моделировании. На основе вышеперечисленных методов и будет базироваться
программная имитация системы.
Во второй части главы 1 приводится описание
пакета ИМОДС, его назначение и функции, описание работы в этом пакете, система
графического представления результатов и оптимизация, что в дальнейшем поможет
эффективнее использовать возможности пакета.
2. Практическая часть
.1 Программная реализация
В качестве среды разработки предпочтительна
среда C++ Builder 6. Выбор обусловлен удобством среды в плане простоты
написания и отладки кода. Простота языка в консоли позволяет сконцентрировать
внимание на самих методах моделирования объекта управления, а не на создании
"красивого интерфейса".
Цели, которые ставятся перед этой средой,
ограничиваются написанием имитационной модели САУ, приведенной в задании.
Задачи. На месте задачи стоит вывод результатов
моделирования в файл, в соответствии с требованиями задания.
моделирование программа параметр
файл
2.1.1 Требования к программе
К программной реализации объекта управления
предъявляются следующие требования:
реализация программным методом;
с помощью пакета ИМОДС симплекс методом;
программа на C++ Builder 6;
запись результатов в файл.
.1.2 Описание программы и
моделируемых объектов
Структурная схема уже была приведена в задании,
внесем обозначения в соответствии с программной реализацией:
Рис. 4
Получает
входной сигнал ys и, путем простого сложения с величиной обратной связи
y_tau, вырабатывает входной сигнал e для регулятора.
|
На
входы пропорциональной и интегральной составляющих приходит выход с первого
сумматора (e). Пропорциональная составляющая моделируется путем умножения
входного сигнала на коэффициент k1.
|
Первый
сумматор
=
ys - y_tau
ПИ-регулятор
_reg1
= e * k1;
y_reg2
= y_reg2 + k2*e*step;_reg = y_reg1 + y_reg2;
//
Выход регулятора
Второй
сумматор
Получает
входной сигнал y_reg и, путем простого сложения с величиной параметрического
возмущения в канале f, вырабатывает входной сигнал y_f для объекта
регулирования.
|
Сам
объект регулирования является инерционным звеном второго порядка.
|
_f
= y_reg + f
Объект
регулирования
По полученной формуле производим решение
дифференциального уравнения второго порядка методом Рунге-Кутта. Значения
производных хранятся в двух переменных: dy1 и dy2, которые считают
соответственно первую и вторую производные.
= (y_f - y - dy1*T1)/T2; // Вторая производная
rk1 = dy2;= dy2+rk1/2;= dy2+rk2/2;=
dy2+rk3/2;
dy1 = dy1 + step*(rk1 + 2*rk2 + 2*rk3 + rk4)/6;
// Первая производная
rk1 = dy1;= dy1+rk1/2;= dy1+rk2/2;=
dy1+rk3/2;= y + step*(rk1 + 2*rk2 + 2*rk3 + rk4)/6; // Выход
Моделирование запаздывания
Рис. 5
Для реализации запаздывания создан массив Z,
размерность которого n равна n = tau/step. В начале, массив необходимо
обнулить.
Реализация выполнена следующим образом:
y_tau = Z[i]; // Запаздывание[i]
= y;
i++;(i == n)= 0;
Код программы представлен в Приложении 1.
2.1.3 Результаты моделирования
Результаты работы выводятся в файл Lika.xls в 4
колонки. В первой колонке записано время, во второй колонке - y_tau, в третьей
колонке - сигнал ошибки e, в четвёртой колонке - скорость изменения сигнала
ошибки de. Файл можно открыть программой Microsoft Excel.
При помощи программного графического пакета
Microsoft Excel строим график переходного процесса системы и график фазового
портрета САУ. графики представлены для двух заданных случаев - входное
воздействие либо 0, либо 1. Динамические свойства системы оцениваем с помощью
фазового портрета САУ для двух вышеуказанных случаев. Фазовый портрет системы
представляет собой зависимость производной (скорости изменения) сигнала ошибки
от сигнала ошибки.
В ходе программного моделирования оптимизация не
производится.
Результаты, полученные при программном
моделировании представлены в Приложении 2.
.2 Моделирование в системе ИМОДС
.2.1 Прогон модели в ИМОДС
Для построения имитационной модели рассмотрим
подробнее, что из себя представляет данная система. Система автоматического
регулирования включает в себя ПИ - регулятор, сумматор на два входа,
сравнивающее устройство и объект управления, представленный инерционным звеном
II порядка с запаздыванием. Для начала примем задающее воздействие равное нулю,
на вход объекта управления подадим возмущение в виде "единичной
ступеньки". (Затем, наоборот, возмущение будет равно 0, а задающее воздействие
1).
Рис. 6
Структурная схема системы управления,
реализованная в ИМОДС
Перед заполнением таблицы по структурной схеме
САУ представляются порядковые номера связей. Они нумеруются в порядке
возрастания от 2 до 999 (0 и 1 являются зарезервированными числами, служащими
для задания соответственно нулевого и единичного возмущения).
Для реализации структурной схемы модели на
системе имитационного моделирования представим ее в качестве таблицы:
Вход
ТЭС
|
Код
ТЭС
|
Выход
ТЭС
|
Имя
Выхода
|
Параметры
Настройки
|
|
|
|
|
k1
|
k0
|
T2
|
0
|
150
|
10
|
e
|
|
|
|
50
|
150
|
10
|
|
|
|
|
10
|
200
|
20
|
|
0.51
|
0.028
|
|
1
|
152
|
30
|
|
|
|
|
20
|
152
|
30
|
|
|
|
|
30
|
104
|
40
|
|
1
|
200
|
20
|
40
|
130
|
50
|
y
|
20
|
|
|
10
|
120
|
60
|
de
|
|
|
|
Для второго случая ys = 1; f = 0 необходимо
будет лишь поменять местами в таблице входы ТЭС первой и четвертой строчки.
В случае, если мы хотим увидеть графическое
представление значения какого-нибудь выходного параметра, необходимо присвоить
этому параметру имя (произвольный набор ASCII - символов). Каждому ТЭС должен
быть присвоен код, зарезервированный в библиотеке ТЭС (см. Приложение 5).
Заполнение таблицы удобнее вести построчно,
начиная с первой строки. Структурная схема поэлементно переносится в таблицу.
Первые две строки в этой таблице соответствуют сумматору на два входа, имеющему
входы "+" и
"-"; об этом говорит код ТЭС во второй
колонке - 150. В первой колонке стоят два входа ТЭС: 0 и 50 (выход y с объекта
поступает по цепи обратной связи на сравнивающий элемент); в третьей колонке -
выход 10.
В четвертой колонке - имя выхода s1. Третья строка
соответствует ПИ - регулятору, код ТЭС - 200, вход - 10, выход- 20, имя выхода
- r. При нахождении светового маркера на этой строке, в "шапке" пятой
колонки появляются названия параметров ТЭС, которые нужно задавать в
соответствующих подколонках. Подобным образом заполняются все строки. В
последней строке стоит дифференциатор, код ТЭС - 120, на его вход подается
ошибка - 10, на выходе - производная; он нужен для построения фазовых
траекторий.
После заполнения таблицы производится запуск ИМ
в работу; для этого выполняются следующие операции:
выход в главное меню (посредством нажатия на
клавишу F10);
задание начала, конца и шага переходного
процесса (посредством активизации подменю "Время");
запуск ИМ в работу (посредством активизации
команды "Работа" из подменю "Опции").
Графические результаты прогона модели
(переходные характеристики и фазовые портреты) представлены в Приложении 3 для
двух заданных случаев.
2.2.2 Оценка параметров системы
Для случая ys = 0 & f = 1
А1 = 0,92; А2 = 0,36; А3 = 0,2;
Максимальное динамическое отклонение: hуст = 0;
Максимальное перерегулирование: smax
= 92%
(для работоспособных систем σmax
≤
20÷30%);
Колебательность: y = 73,3%
(для работоспособных систем ψ ≥ 75÷90%)
Время регулирования: tpег = 327 с.
Для случая ys = 1 & f = 0
А1 = 0,374; А2 = 0,22; А3 = 0,1;
Максимальное динамическое
отклонение: hуст = 1;
Максимальное перерегулирование: smax = 37,4%;
Колебательность: y = 73,3%
Время регулирования: tpег = 340 с.
Из вышеперечисленных расчетов можно сделать
вывод о неработоспособности и неэффективности данной системы в обоих случаях.
Показатели (smax, y) не удовлетворяют
основным показателям качества, поэтому требуется коррекция (оптимизация)
параметров системы.
.2.3 Оптимизация САУ в ИМОДС
Для оптимизации необходимо добавить ещё одну
строчку в нашу таблицу:
Вход
ТЭС
|
Код
ТЭС
|
Выход
ТЭС
|
Имя
Выхода
|
Параметры
Настройки
|
|
|
|
|
k1
|
T2
|
0
|
150
|
10
|
e
|
|
|
|
50
|
150
|
10
|
|
|
|
|
10
|
200
|
20
|
|
0.51
|
0.028
|
|
1
|
152
|
30
|
|
|
|
|
20
|
152
|
30
|
|
|
|
|
30
|
104
|
40
|
|
1
|
200
|
20
|
40
|
130
|
50
|
y
|
20
|
|
|
10
|
120
|
60
|
de
|
|
|
|
10
|
328
|
100
|
|
|
|
|
). Установка параметров оптимизации для ys =0
& f = 1
Рис. 7
Результат оптимизации
Рис. 8
Результаты для разных вариантов количества
итераций
Число
итераций
|
J
|
k1
|
k0
|
6
|
43,735
|
0,1906
|
0,0221
|
10
|
43,204
|
0,1966
|
0,0203
|
16
|
43,192
|
0,1968
|
0,0207
|
По результатам работы системы ИМОДС квадратичный
критерий качества (J), а также оптимизируемые параметры, начиная с 16 итерации
не меняются, т.е. можно сказать, что достигнут оптимум.
). Установка параметров оптимизации для ys =1
& f = 0
Рис. 9
Результат оптимизации
Рис. 10
Результаты для разных вариантов количества
итераций
Число
итераций
|
J
|
k1
|
k0
|
6
|
40,406
|
0,5528
|
0,0221
|
10
|
40,406
|
0,5547
|
0,0222
|
16
|
40,406
|
0,5551
|
0,0222
|
Интегральный критерий качества и параметры k1,
k2 практически не изменились с шестой итерации, таким образом, можно сказать,
что для данного случая достаточно и шести итераций, так как нет существенных
изменений.
Подставляя параметры ПИ - регулятора, полученные
в результате оптимизации, значительно улучшились качественные показатели САУ.
Для того чтобы оценить качество оптимизации,
составим сводную таблицу по критериям качества САУ до оптимизации и после
оптимизации:
Критерий
Качества
|
До
оптимизации
|
После
оптимизации
|
|
ys=0,
f=1
|
ys=1,
f=0
|
ys=0,
f=1
|
ys=1,
f=0
|
Время
регулирования
|
323
c
|
340
c
|
275
с
|
217
с
|
Максимальное
перерегулирование
|
92%
|
37,4%
|
96%
|
15,8%
|
Колебательность
|
73,3%
|
73,3%
|
93,8%
|
89,9%
|
Количество
оборотов фазового портрета
|
4
|
3
|
2
|
3
|
Как видно из сводной таблицы, качество системы
заметно улучшилось, и во втором случае система даже стала работоспособной и
эффективной.
В ситуации ys=0 & f = 1 значение
максимального перерегулирования (s=96%) было слишком
высоко, чтоб удовлетворять стандартам; вероятно, это обуславливается исходными
параметрами объекта (инерционного звена второго порядка).
Графические результаты, полученные с помощью
системы имитационного моделирования «ИМОДС», после проведения оптимизации
представлены в Приложении 4.
Выводы по главе 2
Сравним результаты, полученные из созданной
программы и пакета ИМОДС:
Динамические показатели качества для ys = 0
& f = 1
А1 = 0,900607; А2 = 0,3059; А3 = 0,159506
Реализация
|
s,%
|
y,%
|
tpег,
с
|
Программа
|
90
|
82,3
|
327
|
ИМОДС
|
92
|
73,3
|
327
|
Динамические показатели качества для ys = 1
& f = 0
А1 = 0,323686; А2 = 0,283526; А3 = 0,079816
Реализация
|
s,%
|
y,%
|
tpег,
с
|
Программа
|
32,4
|
75,3
|
330
|
ИМОДС
|
37,4
|
73,3
|
340
|
В главе описывается выбор и обоснование среды
программирования. На основании приведенных в главе 1 методов моделирования САУ,
произведена программная реализация системы.
В программной интерпретации был использован
математический метод решения дифференциальных уравнений (метод Рунге-Кутта) при
моделировании объекта управления.
По окончании работы программы получен текстовый
файл с результатами САУ. С помощью графического пакета Microsoft Excel были
построены графики переходного процесса системы, которые представляются для двух
случаев (ys = 0 & f = 1; ys = 1 & f = 0).
С помощью пакета ИМОДС были построены графики
переходного процесса и фазового портрета САУ, также для двух случаев входного
сигнала.
На основании непосредственного сравнения
графиков полученных с помощью разработанной программы и пакета ИМОДС, можно
сделать вывод об их сходстве, что говорит о правильности выполненного
моделирования.
Оценка параметров системы показала, что наша
система является неработоспособной, поэтому с помощью пакета ИМОДС была
проведена оптимизация системы. Добившись наилучших значений параметров системы
при оптимизации и еще раз прогнав систему, получили хорошие показатели
качества, соответствующие работоспособной системе.
Заключение
В курсовой работе приведено описание различных
методов моделирования систем автоматического управления, промышленных
регуляторов, объектов управления. Изложен метод Рунге-Кутта для решения
дифференциальных уравнений второго порядка. Приводится описание пакета ИМОДС,
его назначение и функции, описание работы в этом пакете, систему графического
представления результатов и оптимизацию.
Были выполнены поставленные перед нами задачи:
В данной курсовой работе была выполнена
программная реализация на языке С++ в среде Borland C++ Builder имитационной
модели системы управления, состоящей из ПИ-регулятора и инерционного объекта
2-го порядка с запаздыванием.
Моделирование работы системы управления
выполнено для двух случаев: ys=0 & f = 1 и ys =1 & f = 0. Данные,
полученные в результате моделирования, использованы для построения графиков
переходных процессов и фазовых портретов в MS Excel. Определены динамические
показатели качества САР.
Изучены принципы работы с системой моделирования
IMODS, и с ее помощью выполнен прогон имитационной модели для двух выше
указанных случаев. Также были определены характеристики переходных процессов,
полученных при помощи IMODS.
Выполнена оптимизация параметров (k0, k1)
ПИ-регулятора для обоих случаев средствами IMODS. В качестве критерия
оптимизации использовался интегральный квадратичный критерий качества.
Оптимизация проводилась симплекс - методом. Сравнение показателей качества
переходных процессов, построенных при заданных настройках регулятора и
оптимальных, показывает их значительное улучшение при использовании оптимальных
настроек.
Список литературы
1. Клиначёв Н.В. Теория
систем автоматического регулирования и управления: Учебно-методический
комплекс. - Offline версия 3.6. - Челябинск, 2005. - 652 файла.
2. Кельтон В., Лоу А.
Имитационное моделирование. - СПб.: Питер; Киев: Издательская группа ВНV, 2004.
- 847 с.
. Советов Б.Я.
Моделирование систем: Учеб. для вузов/ - М: Высш. Шк., 2005.- 343 с.
. Автоматизированное
проектирование систем управления. Под ред. М. Джамшиде, Ч. Дж. Херчета. - М.:
Машиностроение, 1989. - 340 с.
. Максимей И.В.
Имитационное моделирование на ЭВМ. -М.: Радио и связь, 1988. - 231 с.
. Труды Московского
ордена Ленина Энергетического института. Вып. XVII. Под ред. Г.И. Круга. - М.:
МЭИ, 1986. - 212 с.
. Дейтел Х, Дейтел П.
Как программировать на С++. М.: BHV, 2001.
Приложение 1
Код программы на языке Borland C/C++
#include <conio.h>
#include <stdio.h>
#include <math.h>main()
{*p;str1[9],str2[9],str3[9],str4[9];tab
= '\t',enter = '\n';tn = 0, tk = 1000;k1 = 0.51, k2 = 0.028; // Параметры
регулятора
float T1 = 20, T2 = 200; // Параметры объектаtau
= 20;
float step = 0.1, e = 0; // Шаг
и
ошибкаrk1,
rk2, rk3, rk4;
float ys = 0, f = 1; // Задание и возмущениеep =
0; // Предыдущее значение ошибки по заданиюde = 0; // Производная ошибки по
заданиюy_reg1 = 0, y_reg2 = 0,y_reg = 0; // Выход регулятораy_f = 0; // Выход
регулятора с возмущениемy = 0, y_tau = 0; // Выход системыdy1 = 0, dy2 = 0; //
Производные 1 и 2 порядка на выходе объекта
float Z[200], n; // Массив
запаздыванияi
= 0, print = 0;= fopen("Lika.xls","wt");(p == NULL)
{("File not open");1;
}= tau/step;(i = 0; i<n; i++)[i]
= 0;= 0;(tn < tk)
{= ys - y_tau; // Значение ошибки= (e -
ep)/step; // Производная ошибки= e; // Запоминаем предыдущее_reg1 = e * k1;
y_reg2 = y_reg2 + k2*e*step;_reg =
y_reg1 + y_reg2; // Выход регулятора_f
= y_reg + f;= (y_f - y - dy1*T1)/T2;= dy2;= dy2+rk1/2;= dy2+rk2/2;= dy2+rk3/2;=
dy1 + step*(rk1 + 2*rk2 + 2*rk3 + rk4)/6;= dy1;= dy1+rk1/2;= dy1+rk2/2;=
dy1+rk3/2;= y + step*(rk1 + 2*rk2 + 2*rk3 + rk4)/6; // Выход
объекта_tau
= Z[i]; // Запаздывание[i]
= y;++;(i == n)= 0;++;(print%100==0)
{(str1,"%f",tn);(str2,"%f",y_tau);(str3,"%f",e);(str4,"%f",de);(&str1,sizeof(str1),1,p);(&tab,sizeof(tab),1,p);(&str2,sizeof(str2),1,p);(&tab,sizeof(tab),1,p);(&str3,sizeof(str3),1,p);(&tab,sizeof(tab),1,p);(&str4,sizeof(str4),1,p);(&enter,sizeof(enter),1,p);
}= tn + step;
}(p);
return 0;
}
Приложение 2
Результаты работы программы на С++
Builder 6
). Графический результат для ys = 0 & f = 1
Переходная характеристика
Фазовый портрет
2). Графический результат для ys = 1 & f = 0
Переходная характеристика
Фазовый портрет
Приложение 3
Результаты, полученные с помощью
пакета ИМОДС до оптимизации системы
). Графический результат для ys =0 & f = 1
Переходная характеристика
Фазовый портрет
). Графический результат для ys = 1 & f = 0
Переходная характеристика
Фазовый портрет
Приложение 4
Результаты, полученные после
оптимизации системы
). Графический результат для ys = 0 & f = 1
Переходная характеристика
Фазовый портрет
1) Графический результат для ys = 1 & f = 0
Переходная характеристика
Фазовый портрет
Приложение 5
Базовый набор библиотеки ТЭС
Типовые динамические звенья:
А) Статические звенья:
- безинерционное звено (пропорциональное звено);
- инерционное звено 1-го порядка;
- инерционное звено 2-го порядка;
- инерционное звено 3-го порядка;
- колебательное звено;
- консервативное.
Б) Интегрирующие звенья:
- идеальное интегрирующее звено;
- реальное интегрирующее звено (с замедлением);
- изодромное звено.
В) Дифференцирующие звенья:
- идеальное дифференцирующее звено;
- реальное дифференцирующее звено;
- звено чистого запаздывания;
Алгоритмы управления:
- ПИ - алгоритм;
- ПД - алгоритм;
- ПИД - алгоритм;
- Аm + I - алгоритм;
- Аm - алгоритм;
- Аm - алгоритм с
переключением;
- Аm2 - алгоритм;
- Аm2 + Im
- алгоритм;
- А - алгоритм;
¸ 221 - Аm
- ПИ - алгоритм.
Вспомогательные элементы:
- сумматор на 2 входа;
- сумматор на 2 входа;
- сумматор на 3 входа;
- сумматор на 3 входа;
- сложное ступенчатое воздействие;
- синусоидальное воздействие;
- генератор шума (случайных воздействий).