Имитационное моделирование САУ с ПИД регулятором при задающем воздействии на фоне помех

 

Федеральное агентство по образованию и науке Российской Федерации

ГОУ ВПО «Норильский индустриальный институт»

Кафедра электропривода и автоматизации технологических процессов и производств

Курсовая работа

по дисциплине «Моделирование объектов и процессов в цветной металлургии»

Тема: «Имитационное моделирование САУ с ПИД регулятором при задающем воздействии на фоне помех»

Выполнил: Арламов А.С.

Группа: АПм - 06

Проверил: Писарев А. И.

Ботвиньев К. М.

Норильск 2009

Содержание

1.Идентификация объекта управления

.1 Краткое описание технологии объекта управления

.2 Производственный пассивный эксперимент над ОУ

.3 Передаточная функция и частотные характеристики ОУ

1.4 Статистические характеристики случайного процесса

2.Определение передаточной функции формирующего фильтра сигнала помехи

3.Имитационное моделирование САУ при действии сигнала помехи

3.1 Структурная схема имитационной модели САУ

.2 Оптимальные параметры настройки ПИД-регулятора

.3 Соотношение «сигнал/шум» на входе САУ

.4 Соотношение «сигнал/шум» на выходе САУ

.Фильтрация сигнала помехи на входе регулятора

.1 Модуль спектральной плотности задающего сигнала

.2 Передаточная функция фильтра задающего сигнала

.3 Включение фильтра задающего сигнала в модель САУ

.4 Имитационное моделирование САУ с фильтром задающего сигнала. Определение соотношения «Сигнал/шум» на выходе фильтра и выходе САУ

5. Список использованной литературы

1.Идентификация объекта управления

1.1Краткое описание технологического объекта управления (ОУ)

имитационное моделирование сигнал помеха

В качестве объекта управления была выбрана флотомашина ФПМ-16 2-4 цепи аппаратов II секции никель-пирротиновой флотации НОФ.

Транспортно-технологическая схема переработки ЛПК состоит из его гидродобычи из пирротинохранилища №1 с помощью земснаряда, перекачки в виде пульпы в отстойник «Южный», затем материал из «Южного» подается в цех сгущения ИФЦ НОФ на 7 (5) сгуститель. Сгущенный до 40-45% тв. ЛПК насосами 7-2, 7-2а подается на доизмельчение и классификацию (мельница 4-5, 11) с целью «обдирки» окисленных поверхностей ЛПК.

Далее материал самотеком поступает в пульподелитель перед основной флотацией. Перед основной флотацией проводится агитация во флотомашине ФПМ-16 № 2-6 (камеры 1,2) не менее 10 минут с целью улучшения контакта сульфидных минералов с реагентами для последующей стабилизации процесса флотации. Основная никель-пирротиновая флотация проводится во флотомашине ФПМ-16 №№ 2-6 (камеры 3,4), 8, 10, 12. Пенный продукт основной флотации насосами №№ 2-41, 42 направляется насосами в пульподелитель перед перечистной флотацией. Перечистная флотация проводится во флотомашинах ФПМ-16 №№ 2-2,4). Пенный продукт перечистной флотации является готовым никель-пирротиновым концентратом и насосами №№ 2-45,46 перекачивается на никелевый узел. Камерный продукт основной и перечистной флотаций объединяется и является отвальным продуктом от обогащения ЛПК. Отвальный продукт направляется на гравиообогащение текущих отвальных хвостов.

Учитывая переменный состав ЛПК по содержанию никеля (от 1,0 до 2,0%), определены зависимости содержания никеля в хвостах флотации ЛПК от содержания никеля в исходном питании.

В данной работе математическая модель объекта управления (ОУ) будет определяться по каналу «расход воздуха - уровень пульпы».

На рисунке 1.1.1 приведена структурная схема объекта управления.

X(t)Y(t)

Рис.1.1.1 Структурная схема ОУ

В качестве входного сигнала X(t) был взят расход воздуха на выходе насоса, измеряемый в м3/ч; а в качестве выходного Y(t) - уровень пульпы во флотомашине, измеряемый в мм.

1.2Производственный пассивный эксперимент над ОУ

Над одномерным ОУ был проведен производственный пассивный эксперимент, результатом которого стали реализации сигналов на входе X(t) и выходе Y(t) ОУ в установившемся режиме работы.

Объём выборки составляет 200 точек; дискретность - 1 минута.

Данные пассивного эксперимента приведены в таблицах 1.2.1 и 1.2.2.

Табл. 1.2.1. Массив входных данных

X(t), м3/ч1 2 3 4 1897,97083900,70105898,70831897,881232900,23749902,17188898,44269899,155213900,50836899,62396897,68024902,707284897,94269901,57501896,95624903,018745898,85626900,70831897,04999904,709356900,2323901,46667897,06976907,44697898,97815901,10315895,38019907,721868899,39789902,48541893,98645906,546889900,98956904,2146894,12396907,2093510896,91565904,57916896,09479905,0041511898,89478905,06458898,4729904,3614512900,6781907,74585900,49585903,7458513897,28125905,80627896,99792902,2177114897,61249906,28748896,66669903,303115896,2052907,63123896,72498904,4927416898,15314905,54272897,95935904,8822617898,06146907,18439896,2323902,3864718897,17291907,92603896,83752904,8791519896,94794905,60938897,27502904,6677220894,19269907,92706896,45831904,2614721897,78644908,06976897,50623905,8947822898,67499908,61249896,3125906,3635323894,65936908,21039896,37292906,1437424896,12604907,65002900,74481906,2781425897,74689907,75311904,88751906,6989726897,96667908,97919903,38751905,1541727898,76563902,60211901,22815904,1333628896,94586897,05103899,54376904,5072629899,23334898,27081901,33856905,2249830900898,65625901,59479907,1500231896,49896896,58228900,8479905,7489632898,78333895,38544900,50311899,5468833899,72919898,27399898,93854897,2354134900,81561898,4823897,34894895,2125235898,51666899,09998894,63751894,6947636895,90625895,52814896,53229894,2645937898,98645898,67706896,93231894,0760538896,75415897,98022896,33856894,972939900,95935898,73126896,99377897,3531540898,05939897,41248896,4552896,4260341896,78021899,40002896,76874896,3468642900,54895898,13959894,13226894,5104443899,00208895,6073896,29999899,0312544898,55731897,41876897,08228900,2874845902,38226896,83539897,53748898,9447646899,62811895,77081894,25898,3864747899,8844897,69269896,37396898,9760148900,44165897,76044897,53229902,044849902,0719896,26666899,89063899,2000150900,8656897,22919900,50104897,13959

Табл. 1.2.2. Массив выходных данных

Y(t), мм 1 2 3 4 190,37584,81562587,5937567,59375290,9588,25937585,19062559,659375389,57812587,262592,92574,196875488,1187588,30937584,87582,33125591,4062588,30312585,51562573,84375687,59687588,287586,1687589,196875786,52187590,2687585,0586,16875890,59687586,7312582,86562585,275991,712588,69687581,25312589,2968751086,39687588,85312580,2437585,9968751185,5812580,57584,2589,551290,37187588,56562584,7937590,9906251387,39062588,562581,38437587,83751486,7812587,25937579,04062592,81251585,462589,1937578,10937588,6156251685,77812589,3812581,34062590,1843751784,0187587,62581,33437589,08751886,6187590,23437580,0312589,2751986,66562588,5937579,6589,518752078,04062588,987578,40312589,33752181,45937588,49062577,7937589,552285,00312588,2187577,212586,7968752381,18437585,34062575,17589,252480,0312592,24062574,22812593,0906252582,70937590,287570,49062587,468752682,585,90312572,737584,756252782,20937588,037575,4812587,5468752881,89687591,7062564,987586,806252990,96562587,06562570,35312583,118753083,6187592,92570,25937598,6343753182,97812590,12187569,90937588,0906253284,17583,3562569,581,6656253385,72187591,64687569,1592,093753486,9062589,3937566,8583,353587,17187589,77812565,99062587,5968753679,5562586,262565,0812590,768753789,33437592,78437564,1437579,993753880,05937586,1563,212583,7253990,3589,462562,27187589,2218754086,59062588,30937567,3062585,94181,3187592,70937567,90937592,2031254289,87187589,44687564,2562584,5906254384,94687589,212564,35312591,7593754484,6812585,59062567,40312586,58754589,83437584,82812570,96562593,318754681,44062588,42566,55937587,5218754787,00312578,6937568,47187585,893754887,17187587,96562569,4062587,8906254988,15312583,83437569,37812590,555088,5312581,85937567,2687588,26875

1.3Передаточная функция и частотные характеристики ОУ

По полученным экспериментальным данным с помощью программы «MATLAB System Identification Toolbox» рассчитаем передаточную функцию объекта управления. Окно работающей программы представлено на рисунке 1.3.1.

По результатам работы программы получена передаточная функция ОУ вида:

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

Рис. 1.3.1 Определение передаточной функции ОУ.

Погрешность модели при 8 итерациях составляет 0.000415 %, что допустимо.

Рассчитанный коэффициент передачи: K = 0.095274

Рассчитанная постоянная времени : T = 5.8438

Таким образом, в безразмерных величинах передаточная функция примет вид:

С помощью программы «MATLAB Simulink» найдём частотные характеристики модели ОУ. На рисунке 1.3.2 представлена амплитудно-частотная характеристика (АЧХ) ОУ.

Рис.1.3.2 Амплитудно-частотная характеристика ОУ.

На рисунке 1.3.3 представлена фазо-частотная характеристика (ФЧХ) ОУ.

Рис.1.3.3 Фазо-частотная характеристика ОУ.

На рисунке 1.3.4 представлена логарифмическая амплитудно-частотная характеристика (ЛАЧХ) ОУ.

Рис.1.3.4 Логарифмическая амплитудно-частотная характеристика ОУ

Ниже приведены временные характеристики ОУ: переходная функция на рисунке 1.3.5 и импульсная переходная функция на рисунке 1.3.6.

Рис. 1.3.5 Переходная функция h(t) ОУ

Рис. 1.3.6 Импульсная переходная функция ?(t) ОУ.

1.4 Статистические характеристики случайного процесса на входе ОУ (дисперсия, корреляционная функция, спектральная плотность)

На рисунке 1.4.1 представлена реализация случайного процесса Х(t) на входе ОУ (расход воздуха).

Рис. 1.4.1 Реализация случайного процесса X(t) на входе ОУ.

С помощью программы «MathLab» рассчитаны математическое ожидание, дисперсия и среднеквадратическое отклонение процесса X(t).

На рисунке 1.4.2 приведен график центрированного случайного процесса.

Рис. 1.4.2 Центрированный случайный процесс X*(t).

Приведём автокорреляционную функцию к ненормированному виду (домножим все её значения на дисперсию) и аппроксимируем её экспонентой.

Аппроксимирующая функция будет иметь вид:

На рисунке 1.4.3 приведена автокорреляционная функция и её аппроксимация.

Рис. 1.4.3 Автокорреляционная функция процесса X(t) и её аппроксимация.

В результате работы программы «Advanced Grapher» был определён коэффициент затухания аппроксимирующей функции:

Таким образом, принимаем автокорреляционную функцию случайного процесса X(t) равной:

Для нахождения спектральной плотности случайного процесса X(t) воспользуемся формулой Винера-Хинчина:

Подставим вместо KX(?) полученное выражение корреляции:

Подставляя численные значения, получим:

График спектральной плотности S(?) случайного процесса X(t) представлен на рисунке 1.4.4.

Рис. 1.4.4 Спектральная плотность S(?) случайного процесса X(t).

2.Определение передаточной функции формирующего фильтра сигнала помехи

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

В нашем случае на вход формирующего фильтра поступает сигнал помехи «белый шум» ( SБШ(?) = S0 = const ), на выходе же требуется получить сигнал со спектральной плотностью, рассчитанной в п. 1.4.

Отсюда находим комплексную частотную передаточную функцию формирующего фильтра:

Переходя от мнимой частоты j? к комплексной p, получаем передаточную функцию формирующего фильтра сигнала помехи:

Структурная схема имитационной модели формирующего фильтра сигнала помехи показана на рис. 2.1.

Рис. 2.1 Модель формирующего фильтра сигнала помехи.

Графики белого шума единичной мощности и отфильтрованного сигнала приведены на рисунках 2.2 и 2.3 соответственно.

Рис. 2.2 График сигнала белого шума на входе формирующего фильтра.

Рис. 2.3 График отфильтрованного сигнала на выходе формирующего фильтра.

Рассчитаем амплитудно-частотную характеристику полученного формирующего фильтра:

Амплитудно-частотная характеристика (АЧХ) и логарифмическая амплитудно-частотная характеристика (ЛАЧХ) формирующего фильтра сигнала помехи приведены на рисунках 2.4 и 2.5 соответственно.

Рис. 2.4 АЧХ формирующего фильтра сигнала помехи.

3.Имитационное моделирование САУ при действии сигнала помехи

.1 Разработка структурной схемы имитационной модели САУ

В качестве регулирующего воздействия на систему был взят прямоугольный импульс амплитудой А = 80 (А) и длительностью tИ = 10 минут.

Модель входного импульса реализуем с помощью источника «Step».

Воздействие, поступающее на вход САУ, представляет собой сумму входного импульса, являющегося полезным сигналом, и сигнала помехи.

Имитационная модель воздействия поступающего на вход САУ представлена на рис. 3.1.1. График полезного сигнала с наложенной на него помехой представлен на рис. 3.1.2.

Рис. 3.1.1 Структурная схема модели входного воздействия.

Рис. 3.1.2 График входного воздействия (полезный сигнал + помеха).

Теперь, получив входной сигнал, разработаем структурную схему САУ. Она будет состоять из сигнала задания (с наложенной на него помехой), ПИД-регулятора и собственно объекта управления.

На рисунке 3.1.3 приведена структурная схема разрабатываемой САУ.

Рис. 3.1.3 Структурная схема САУ.

На рисунке WРЕГ (p) - передаточная функция ПИД-регулятора, имеющая вид:

WОР (p) - передаточная функция объекта управления, равная:

.2 Расчет оптимальных параметров настроек ПИД-регулятора по заданным показателям качества регулирования при отключенном источнике шума с помощью программы «MATLAB Simulink Response Optimization»

Для получения настроек ПИД-регулятора составлена схема имитационной модели, представленная на рисунке 3.2.1. В качестве задания взята «ступенька».

Рис. 3.2.1 Схема получения настроек ПИД-регулятора

Настройку произведём исходя из следующих требований:

·Перерегулирование не более 15%

·Отсутствие колебательности на выходе САУ.

·Минимальное время переходного процесса.

На рисунке 3.2.2 приведено окно работы программы «MATLAB Simulink Response Optimization», рассчитывающей настройки ПИД-регулятора.

Рис. 3.2.2 Получение настроек ПИД-регулятора методом итерации.

Ниже представлен результат работы программы:

max Directional First-orderS-count f(x) constraint Step-size derivative optimality Procedure

1 0 0.02752

14 0 0.004149 25.4 0 97

21 0 0.000142 5.62 0 1.49termination.a feasible or optimal solution within the specified tolerances.

Таким образом, имеем передаточную функцию ПИД-регулятора:

Теперь, когда входной сигнал и ПИД-регулятор рассчитаны, проведём имитационное моделирование работы САУ согласно схеме на рисунке 3.2.3.

Рис. 3.2.3 Структурная схема модели САУ при воздействии помехи

На рисунке 3.2.4 представлен график реакции САУ на воздействие полезного сигнала с наложившейся на него помехой.

Рис.3.2.4 Реакция САУ на полезный сигнал с помехой

Для определения качества фильтрации полезного сигнала рассчитываются соотношения « сигнал / помеха » на входе и выходе САУ с использованием фильтра помехи и без него.

3.3 Определение сигнала ошибки на входе ПИД-регулятора путём имитации работы САУ. Расчёт соотношения « сигнал / помеха » на входе системы

Отношение « сигнал / помеха » во всех случаях будем рассчитывать как отношение энергий соответствующих сигналов.

Средняя энергия сигнала s(t) рассчитывается следующим образом:

Примем время интегрирования от 0 до 10 минут, так как полезный сигнал обладает ненулевой энергией только на этом промежутке.

На рисунке 3.3.1 представлена структурная схема определения средней энергии полезного сигнала, подаваемого на САУ.

Рис. 3.3.1 Структурная схема определения средней энергии полезного сигнала

Энергия полезного сигнала составляет 3,84?106 (А2).

Найдём аналогичным способом среднюю энергию сигнала помехи, подаваемого на вход САУ. На рисунке 3.3.2 представлена структурная схема определения средней энергии сигнала помехи.

Рис. 3.3.2 Структурная схема определения средней энергии сигнала помехи

Энергия сигнала помехи составляет 3,617?104 (А2).

Рассчитаем соотношение « сигнал / помеха » на входе САУ:

Рассчитаем дисперсию входного сигнала с помощью программы «Microsoft Excel».

3.4 Определение реакции САУ путём имитационного моделирования. Расчёт соотношения « сигнал / помеха » на выходе системы

Определим реакцию САУ на полезный сигнал и найдём среднюю энергию полезного сигнала на выходе САУ по схеме, приведённой на рисунке 3.4.1.

Рис. 3.4.1 Структурная схема определения реакции САУ на полезный сигнал и средней энергии полезного сигнала на выходе САУ.

На рисунке 3.4.2 приведена реакция САУ на полезную составляющую сигнала задания.

Рис. 3.4.2 График реакции САУ на полезную составляющую сигнала задания.

Определим реакцию САУ на сигнал помехи и найдём среднюю энергию сигнала помехи на выходе САУ по схеме, приведённой на рисунке 3.4.3.

Рис. 3.4.3 Структурная схема определения реакции САУ на сигнал помехи и средней энергии сигнала помехи на выходе САУ.

На рисунке 3.4.4 приведена реакция САУ на сигнал помехи.

Рис. 3.4.4 График реакции САУ на сигнал помехи.

Рассчитаем соотношение « сигнал / помеха » на выходе САУ

Как видно, соотношение « сигнал / помеха » на входе и выходе САУ примерно одинаково.

Для увеличения этого соотношения на выходе САУ, необходимо отфильтровать сигнал задания на входе регулятора.

.Фильтрация сигнала помехи на входе регулятора

.1 Расчёт модуля спектральной плотности задающего сигнала

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

Рассчитаем модуль спектральной плотности входного импульса, используя прямое преобразование Фурье:

На рисунке 4.1.1 приведён график модуля спектральной плотности полезного сигнала.

Рис. 4.1.1 Модуль спектральной плотности полезного сигнала.

4.2 Определение АЧХ фильтра сигнала по графику модуля спектральной плотности задания. Выбор передаточной функции фильтра

В качестве фильтра входного сигнала выберем фильтр Баттерворта второго порядка. Его передаточная функция имеет вид:

Выразим из передаточной функции фильтра его АЧХ:

где Т - постоянная времени фильтра.

?ср - частота среза, равная: ?ср = 1/Т.

По графику на рисунке 4.1.1 определим частоту среза исходя из того, что большая часть площади под графиком модуля спектральной плотности полезного сигнала приходится на интервал от 0 до 0,011 рад/с.

сек

На рисунке 4.2.1 приведена АЧХ полученного фильтра.

Рис. 4.2.1 Амплитудно-частотная характеристика фильтра Баттерворта.

Таким образом, передаточная функция фильтра будет равна:

.3 Включение блока фильтра задающего сигнала на вход ПИД-регулятора в структуру модели САУ

Включим полученный фильтр на вход регулятора по схеме, приведённой на рисунке 4.3.1.

Рис. 4.3.1 Структурная схема САУ с фильтром на входе регулятора

На выходе САУ был снят отфильтрованный сигнал, изображённый на рисунке 4.3.2.

Рис. 4.3.2 Реакция САУ на отфильтрованный сигнал

Как видно, форма выходного сигнала искажается вследствие подавления фильтром высокочастотных гармоник полезного сигнала.

Однако фильтр устраняет влияние помехи. В спектр помехи теперь входят лишь низкочастотные гармоники.

.4 Определение соотношения « сигнал / помеха » на выходе фильтра (входе регулятора) и выходе САУ

Соотношение « сигнал / помеха » на входе и выходе САУ с фильтром определим тем же методом, что применялся ранее.

Во всех случаях расчёта средней энергии сигналов примем интервал времени от 0 до 1400 секунд.

На рисунке 4.4.1 приведена схема определения средней энергии полезного сигнала на выходе фильтра и выходе САУ.

Рис. 4.4.1 Структурная схема определения средней энергии полезного сигнала на выходе фильтра и выходе САУ

На рисунке 4.4.2 приведены графики сигналов на выходе фильтра и выходе САУ при отключенном сигнале помехи.

Рис. 4.4.2 Графики сигналов на входе регулятора и выходе САУ при отключенном сигнале помехи

На рисунке 4.4.3 приведена схема определения средней энергии сигнала помехи на выходе фильтра и выходе САУ.

Рис. 4.4.3 Структурная схема определения средней энергии сигнала помехи на выходе фильтра и выходе САУ

На рисунке 4.4.4 приведены графики сигналов на выходе фильтра и выходе САУ при отключенном полезном сигнале.

Рис. 4.4.4 Графики сигналов на входе регулятора и выходе САУ при отключенном полезном сигнале

Рассчитаем соотношение « сигнал / помеха » на выходе САУ с фильтром:

Сравним этот результат с соотношением « сигнал / помеха » на выходе САУ без фильтра (Ксигнал/помеха_1 = 111,53) :

Из этого следует, что фильтр подавляет энергию помехи на выходе САУ примерно в 45 раз.

Такой показатель фильтрации нас вполне устраивает, однако полезный сигнал вместе с этим существенно искажается. Поэтому такая фильтрация приемлема не везде.

Наиболее подходящее применение такая фильтрация может найти в тех системах, где требуется преобразование входного прямоугольного импульса в импульс трапецеидальной формы (см. рис. 4.4.2).

5.Список использованной литературы

1.Тарасик В. П. Математическое моделирование технических систем: учебник для вузов. - М.: ДизайнПРО, 2004. - 640 с.: ил.

2.Зарубин В. С. Математическое моделирование в технике: Учебник для вызов Под ред. В.С.Зарубина, А.П.Крищенко. - М.: Изд-во МГТУ им. М.Э.Баумана, 2001. - 496 с.

.Дьяконов В. П. Simulink 4. Специальный справочник. - СПб: Питер, 2002. - 528 с.: ил.

.Дьяконов В. П. MatLab 6.5 SP1/7+Simulik 5/6 в математике и моделировании. Серия «Библиотека профессионала». - М.: СОЛОН-Пресс, 2005. - 576 с.

.Лукас В. А., Теория автоматического управления. - М.: Недра, 1990. - 416 с.: ил.

.Гмурман В. Е., Теория вероятностей и математическая статистика. - М: Высшая школа, 2003.

.Технологическая инструкция НОФ. - Норильск, 2009.

Федеральное агентство по образованию и науке Российской Федерации ГОУ ВПО «Норильский индустриальный институт» Кафедра электропривода и автомати

Больше работ по теме: