автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.01, диссертация на тему:Разработка и исследование модернизированного многоканального ПД-регулятора для стабилизации режимов работы теплоэнергетического котла

На правах рукописи

ПОНОМАРЕВ АЛЕКСАНДР АЛЕКСАНДРОВИЧ

РАЗРАБОТКА И ИССЛЕДОВАНИЕ МОДЕРНИЗИРОВАННОГО МНОГОКАНАЛЬНОГО ПД - РЕГУЛЯТОРА ДЛЯ СТАБИЛИЗАЦИИ РЕЖИМОВ РАБОТЫ ТЕПЛОЭНЕРГЕТИЧЕСКОГО КОТЛА

Специальность: 05.13.01 - Системный анализ, управление и обработка информации (промышленность)

АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата технических наук

1В МАЙ 2013

005058638

Новосибирск - 2013

005058638

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Новосибирский государственный технический университет»

Научный руководитель: заслуженный деятель науки РФ,

доктор технических наук, профессор Востриков Анатолий Сергеевич

Официальные оппоненты: Мухопад Юрий Федорович,

заслуженный деятель науки и образования РФ, доктор технических наук, профессор, Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Иркутский государственный университет путей сообщения», кафедра управления техническими системами, профессор Толстиков Александр Сергеевич, доктор технических наук, доцент, Федеральное государственное унитарное предприятие «Сибирский государственный ордена Трудового Красного Знамени научно-исследовательский институт метрологии», отдел государственной службы времени, частоты и определения параметров вращения Земли, начальник

Ведущая организация: Федеральное государственное бюджетное учреждение науки Институт автоматики и электрометрии Сибирского отделения Российской академии наук (ИАиЭ СО РАН)

Защита диссертации состоится «16» апреля 2013 года в 12.30 часов на заседании диссертационного совета Д 212.173.05 при Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Новосибирский государственный технический университет» по адресу: 630073, г. Новосибирск, пр. К. Маркса, 20

С диссертацией можно ознакомиться в библиотеке Новосибирского государственного технического университета

Автореферат разослан » марта 2013 г. Ученый секретарь

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Для качественного ведения технологического процесса выработки пара необходимо автоматическое регулирование параметров. Улучшение поддержания заданных параметров осуществляется благодаря синтезу систем регулирования с усовершенствованными законами регулирования. Модернизированные законы регулирования позволяют устанавливать высокие требования к показателям качества.

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

Метод локализации предлагает теоретически обоснованную, регулярную методику синтеза систем управления для нелинейных и нестационарных объектов. Исследование является перспективным для развития метода локализации. Область применения метода локализации распространяется на теплоэнергетические объекты. Совершенствуется его практическое применение и увеличивается прикладное значение.

Объект исследований. Теплоэнергетическое оборудование, входящее в состав современных тепловых электрических станций необходимо оснащать системами автоматического регулирования высокого качества и надежности. Отдельно можно выделить газоводушный тракт, являющийся неотъемлемой частью любого котла. Согласно исследованиям, проведенным в США, при увеличении эффективности работы газовоздушного тракта может происходить энергосбережение на уровне 2-20%. С точки зрения теории регулирования объект представляет особый интерес ввиду его многоканальное™ и нелинейных и нестационарных свойств.

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

Задачи исследования:

- разработать математическую модель на основе физических законов протекающих процессов в топке котла,

- провести синтез регуляторов параметров объекта на основе метода локализации,

- усовершенствовать математическую модель путем добавления дополнительных инерционностей, рассмотрения нелинейных и нестационарных свойств объекта,

- применить регулятор для усовершенствованной модели, проверить его работоспособность,

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

Методы исследований. Математические модели строятся аналитическим методом с использованием уравнений, описывающих физические свойства объекта. Синтез системы регулирования проводится на основе метода локализации с использованием различных схем включения регуляторов. Эксперименты над системой и ее последующая корректировка проводится в пакете прикладных программ БтиНпк МаОаЬ. Разработка метода настройки регулятора, не требующего перенастройки в ходе работы, включает применение БипиНпк МайаЬ для идентификации объекта на основе экспериментальных данных. Синтез адаптивного регулятора и проверка его работы для экспериментальных моделей проводится с помощью ЯСАОА-системы.

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

Научная новизна. Разработаны новые математические модели газовоздушного тракта аналитическим методом. К ним относятся линейная, нелинейная и нелинейная с электроприводом. Каждая из моделей может применятся для синтеза систем регулирования в зависимости от необходимой точности описания объекта. Рассмотрено применение метода локализации к теплоэнергетическим объектам, разработаны системы автоматического регулирования на его основе, вследствие чего расширена область прикладного применения метода. Новые модели, регуляторы и результаты экспериментов позволили создать метод настройки регуляторов газовоздушного тракта на основе принципа локализации. В диссертационной работе описано создание нового способа автоматической настройки регулятора (АСНР) на основе принципа локализации, представляющего собой компьютерную программу. АСНР используется для быстрой настройки регуляторов. В рамках исследований предложен новый алгоритм настройки регуляторов для различных объектов с использованием электронно-вычислительной техники.

Основные положения, выносимые на защиту:

- разработана новая линейная модель газовоздушного тракта (ГВТ), для которой проведен синтез ПД-регуляторов по методу локализации,

- разработаны новые нелинейные модели ГВТ и ГВТ с электроприводом, для которых проведен синтез ПД-регуляторов по методу локализации,

- разработан метод настройки регулятора газовоздушного тракта на основе принципа локализации,

- разработана программа автоматической настройки регулятора,

- разработан метод автоматизированной настройки регулятора с использованием ЭВМ.

Практическая ценность. Результаты исследований данной работы могут быть применены для реальных технологических объектов. Математическими

моделями могут описываться множества типовых газовоздушных трактов с различными техническими характеристиками. Для этого в модели должны быть перерасчитаны числовые коэффициенты. Применение регулятора, синтез которого проведен по методу локализации, также возможно для множества газовоздушных трактов с проверкой работоспособности и корректировкой коэффициентов в случае необходимости. Метод настройки регулятора газовоздушного тракта опробован для теплоэнергетического котла ТПЕ-214, установленного на Новосибирской ТЭЦ-5. Автоматизированный метод настройки с использованием ЭВМ применим для технологических объектов, которые требуют регулирования параметров. Метод используется в работе фирмы ЗАО «СИНТЭП».

Внедрение результатов работы. Результаты диссертационной работы использованы в производственной деятельности ЗАО «СИНТЭП»: при разработке программного обеспечения, наладке оборудования. Использование указанных результатов позволяет повысить эффективность работы технологического оборудования.

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

Апробация работы. Основные положения диссертационной работы и отдельные ее результаты докладывались на: VI Международной молодежной научной конференции «Тинчуринские чтения» (г. Казань, КГТУ, 2011 г.), 2-nd Indo-Russian Joint Workshop on Computing Intelligence and Modern Heuristics in Automation and Robotics (Novosibirsk, NSTU, 2011), Конференции «Системный анализ и информационные технологии» (г. Новосибирск, НГТУ, 2011 г.), 1Х-ой Международной научно-практической конференции «Современные инструментальные системы, информационные технологии и инновации» (г. Курск, Юго-Зап. гос. ун-т, 2012 г.), Международной научно-технической конференции студентов, аспирантов и молодых ученых "Прогрессивные направления развития машино-приборостроительных отраслей и транспорта" (г. Севастополь, СевНТУ, 2012 г.).

Публикации. Основные положения и результаты диссертационной работы опубликованы в 9 работах, 3 из которых входят в перечень рецензируемых научных журналов и изданий для опубликования основных научных результатов диссертаций, 2 статьи в сборниках научных трудов и 4 статьи в сборниках материалов международных конференций.

Структура и объем работы. Диссертация состоит из введения, трех глав, заключения, списка использованной литературы из 114 наименований и приложения. Общий объем работы составляет 122 страницы, из которых основное содержание работы - 119 страниц, включая 11 таблиц и 53 рисунка.

СОДЕРЖАНИЕ РАБОТЫ

В первой главе описывается разработка линейной модели газовоздушного тракта (ГВТ) и синтез ПД-регуляторов по методу локализации.

Построение модели ГВТ проводится с последовательным приближением к наиболее точному математическому описанию. Такой способ позволяет проводить проверку соответствия модели реальному объекту по всем его

свойствам. Основные методы проверки это построение переходных процессов и сравнение с разгонными характеристиками технологического объекта, а также анализ и синтез систем автоматического регулирования (САР).

Линейная модель это первый этап математического описания ГВТ. Модель объекта составляется на основании законов физики.

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

Входные переменные (управляющие величины): угловая скорость вращения дутьевого вентилятора - (ОдВ, рад/с, угловая скорость вращения

дымососа - (ОдС, рад/с. Выходные переменные (регулируемые величины): расход общего воздуха - £)вх, м3/с, разрежение в топке - рт, Па. Промежуточная величина - расход дымовых газов - ()вых, м3/с.

Объект имеет несколько взаимовлияющих параметров управления, что позволяет считать его многоканальным и многосвязным. Соответственно управление ГВТ также целесообразно осуществлять в виде многоканальной САР.

Изменение массы вещества т, кг, переносимого ГВТ, описывает закон сохранения массы: = (2вх~йеых- При рассмотрении данного уравнения,

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

Уравнение состояния, как известно из термодинамики, имеет вид: ртУ = тИТ/\х, где рт, Па - давление (разрежение) в топке, V, м3 - объем топки, т, кг - масса газа в топке, ¡1 = 44, кг/Кмоль - молярная масса газа для С02, Я = 8,31-103, Дж/°С'Кмоль - универсальная газовая постоянная, Т,°С -температура газа в топке.

Уравнение позволяет получить соотношение, связывающее массу и давление: т = ксрт, где кс = р^/КГ - кг/Па - коэффициент состояния.

Подставив соотношение в закон сохранения массы, можно получить формулу изменения давления: кс <1рт= (¿^ - <2вЬ1Х.

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

11вАйвх1$в. & = Рдв ~Рт~ Рсопр.6. К^бвыхЛгУ' = Рдс + Рт~ Рсопр.г.'

где Ів, м - длина воздушного участка, , м2 - площадь поперечного сечения воздушного участка, 1г, м - длина газового участка, , м2 - площадь поперечного сечения газового участка, рдв, Па - давление, нагнетаемое дутьевым вентилятором, Рсопр.в. > Па - аэродинамическое сопротивление воздушного участка, рдс, Па - давление, нагнетаемое дымососом, рсопрг , Па -

аэродинамическое сопротивление газового участка.

Величины аэродинамического сопротивления воздушного и газового участков связаны с расходами воздуха и газов соответственно. На начальном этапе для простоты принимается допущение, что расходы воздуха и газов ламинарные. Соотношения описываются линейными гидродинамическими уравнениями: рсопр.в.=^в.Овх,рсопр.г.=^г.Оеых, где ^ , Па-с/м3 - коэффициент аэродинамического сопротивления воздушного участка, \г, Па с/м3 — коэффициент аэродинамического сопротивления газового участка.

В баланс давлений входят величины давлений (напоров) дутьевого вентилятора и дымососа, связанные с угловыми скоростями вращения. На начальном этапе для простоты допускается, что давления дутьевого вентилятора и дымососа зависят линейно от угловых скоростей вращения электродвигателей ДВ и ДС. Соотношения описываются линейными уравнениями дутьевого вентилятора и дымососа: рДВ = к дВ(йдВ, рдс = кдс^дс >

где кдВ, Па-с - коэффициент дутьевого вентилятора, кдС, Па-с - коэффициент

дымососа.

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

йвх=~ £>в.ЗвЯвх/1в. + кдв$вшдв/1«. ~ $а.Р,п1К.

• йеых = - ^гЯеых/к + ЬдсЗгЫдс/к + $г.Рт/к ■ (2)

Рт ~ О-вхі^сост. ~ @выхI^сост.

Система уравнений ГВТ может быть записана в качестве дифференциальных уравнений состояния с числовыми значениями (3)

і] = -6,56*] - 0,2х3 + 3,34«! ■ х2= -1,11*2 + 0Дх3 +1,91 и2 , (3)

і3 = 5000х, - 5000х2

где применены следующие обозначения: ()вх = хх, <2вых = х2, рт=х2, ЮдВ = м1,

содс=и2. Используются параметры теплоэнергетического котла ТПЕ-214,

установленного на Новосибирской ТЭЦ-5. Данный пример используется во всей диссертационной работе.

Так как модель линейная, дифференциальные уравнения состояния ГВТ записываются в векторно-матричном виде:

-6,56 0 -0,2" "3,34 0 "

х = Ах + Ви, А = 0 — 1,11 ОД , в = 0 0

5000 5000 0 0 1,91

Для изучения свойств объекта проводится построение переходных процессов с использованием пакета прикладных программ БтиНпк. Кривые получены как реакция объекта на ступенчатое возмущение по входу. То есть в точке времени 0,2 с изменяется угловая скорость вращения электродвигателей тягодутьевых машин. Такое возмущение характерно при эксплуатации теплоэнергетических установок.

Переходный процесс переменной х3 носит колебательный характер из-за наличия больших коэффициентов в системе (5000 - коэффициент при переменной). В технике такое явление известно как пульсации давления, возникающие вследствие изменения расходов воздуха или уходящих газов в ГВТ. В данной работе математическая модель отражает такое свойство объекта, как пульсации давления. Таким образом, решение проблемы может быть найдено при синтезе регулятора ГВТ.

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

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

Итак, требуется обеспечить выполнение условия Нтх(I) = у, где V -входное воздействие. Условие в статике формулируется так: А(с°)=у(оо)-у<А, где г) - ошибка регулирования, а Д - заданная ошибка

(исходя из опыта практики А = 5%).

Требования к динамике системы, а именно: г„ - время переходного процесса ист- перерегулирование, выражается в виде желаемой передаточной функции, получаемой в соответствии с модальным методом: \У( р ) = С~1( р ).

В представленной системе задаются значения гп = 1 с (с учетом того, что в модели на данном этапе отсутствуют нелинейности и возмущения, а в реальных системах постоянная времени равна 2 - 4 с) и перерегулирование ст = 5%, как допустимое отклонение при эксплуатации САР в теплоэнергетике.

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

разрываются. Каждый полученный подобъект должен иметь свою передаточную функцию, которая записывается как W0 = к0В( р )А'1( р).

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

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

Передаточная функция объекта первого канала, согласно дифференциальным уравнениям, но без учета связей со вторым каналом, выглядит следующим образом: И^ = кък\1( к^1р +1 )~1.

Полиномы передаточной функции, соответственно, записываются: В(р) = 1, А(р) = к?р + \, к0 = к3к[\

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

Порядок относительной производной выхода п = 1. Используя структурный метод, вид желаемого полинома принимает вид С( р ) = Тр +1, или с числовыми значениями для примера: С( р) = р +1.

Согласно способу большого коэффициента и производной от ошибки передаточная функция регулятора будет иметь следующий вид: р )=С( р )к, или с числовыми значениями для примера ГВТ: \Г1(р) = к(Тр + 1).

Характеристическое уравнение замкнутой системы записывается как А( р) + к0кС( р ) = 0. Раскрывая характеристическое уравнение замкнутой системы можно получить следующее выражение:

(к{]1 р + и + к0к(Тр +1) = (1 + к0кТ )р+(к[1 +кк0) = 0. (5) Из выражения (5) следует, что при достаточно большом коэффициенте усиления к, то есть при к = 10 система не зависит от параметров объекта и можно переписать характеристическое уравнение как к0к( Тр +1) = 0.

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

Уравнение, связывающее вход и выход системы выглядит следующим образом: х = С( р )кк0( А( р) + кк0С( р))~1у.

Для данного выражения, при к = 10, выходная величина х = V. Таким образом, система отрабатывает всевозможные входы V с ограниченной старшей производной. Такая система может использоваться в различных задачах слежения за входным сигналом. Для отработки скачкообразных входов необходимо производить сглаживание специальным дифференцирующим фильтром.

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

Рисунок 1. Расчетная схема системы регулирования с дифференцирующим фильтром

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

Дифференцирующий фильтр является вспомогательной следящей системой, его передаточная функция записывается следующим образом: \Уф = Р[( 0(\хр ))~1. Полином, обеспечивающий нужные фильтрующие свойства,

можно записать в общем виде: Б(\1р) = \1Г рг + с1г_1\х.г~1 рг~1 +... +1, где значение коэффициента на порядок меньше времени переходного процесса Гя объекта, а значения коэффициентов выбираются из условий устойчивости быстрых движений. Порядок дифференцирующего фильтра должен быть выше или равен порядку полинома С( р). С учетом этого порядок фильтра следует

выбирать равным п > 1.

В соответствии с требованиями полином регулятора первого канала записывается в виде: цр ) = ц.р + е?0, или с числовыми значениями примера ГВТ Щ\}р) = 0,1р + 1.

Синтез регулятора второго канала проводится по тому же алгоритму, что и для первого канала. Для изолированного второго канала полиномы передаточной функции \У() объекта записываются в виде выражений: В( р) = 1,

А(р) = к6~1р2 + к4к6~1р + к2, к0 = ~к5'1.

Полином В( р) также как и для первого канала имеет значение, при котором выполняется условие разрешимости задачи синтеза.

Порядок относительной производной выхода этой передаточной функции п = 2. Согласно структурному методу желаемый полином получается в виде: С( р ) = Т2р2+ 2с1Тр + 1, или с предварительными числовыми значениями: С(р)=р2 + р +1.

Характеристическое уравнение замкнутой системы выглядит следующим образом:

( к6~1р2 + к4 к6~1р + к2)+ к0к( Т2р2 + 2Тр + 1) =

(о)

= (к6~1 + к0кТ2 )р2+(к4к6~1 + к0к2Т)р+(к2 + к0к) = 0 и при достаточно большом к (для второго канала равного 100) сводится к характеристическому уравнению первого канала, что обеспечивает и для второго канала независимость свойств системы от параметров объекта и способности отработки системой ступенчатого входного сигнала V.

Для регулятора второго канала, как и для первого, необходимо определить дифференцирующий фильтр. Поскольку порядок относительной производной выхода п = 2, то фильтр должен иметь порядок п> 2. Соответствующий полином выглядит следующим образом: \хр ) = ц2/;2 + й?,|хр + с1{), или с числовыми значениями для используемого примера Э( цр) = 0,0\р2 + 0,5/? +1.

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

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

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

С целью уточнения коэффициентов настройки регуляторов для САР ГВТ теплоэнергетического барабанного котла ТПЕ-214 проводится числовая имитация.

Для моделирования в 51пш1шк с учётом того, что большой разброс коэффициентов объекта не соответствует возможностям программы, необходимо замедлить объект пропорционально величине корней системы дифференциальных уравнений. Проводится построение переходных процессов разрежения в топке х3 для различных значений коэффициентов = 6,56; -2; -10 и к2 = ~0,2', 3; 10 в системе дифференциальных уравнений. Данные коэффициенты изменяются вследствие нестационарности объекта. Заданный регулятором вид переходного процесса, одинаков для различных режимов работы и параметров объекта даже при сильном изменении коэффициентов объекта.

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

Допущением линейной модели является то, что потоки воздуха и газов ламинарные. Однако более точно поток газовоздушной среды характеризуется уравнением Дарси-Вейсбаха (соотношением между скоростью движения газа и потерями давления на трение для газопровода конечной длины): р =Лу2р//2^, где р,Па - потери давления в газопроводе на трение, X = 0,06 - безразмерный коэффициент гидравлического трения (принимается одинаковым для воздушного и газового участков при турбулентном движении в зоне шероховатых труб), й,м - диаметр газопровода, v,м|c - скорость движения

газа, р = \,Ъкг/м3 - плотность общего воздуха и уходящих газов (принимаем одинаковой), 1,м - длина газопровода.

С использованием уравнения Дарси-Вейсбаха записываются нелинейные соотношения между расходами общего воздуха, уходящих газов и аэродинамическими сопротивлениями участков: РСОпр.в.=^0вхЬ./п2Р^». >

Рсопр.г. = ^'■0вых1г./п2р^г.' гДе = 2,76л» - диаметр газопровода воздушного участка, с1г = 3,51м - диаметр газопровода газового участка.

Уравнения упрощаются с использованием коэффициентов Рсопрл. = ксоИр*&> Рсопр.г =ксопр.г.0.вых, ГДе ксопр = Шв .М^в.= 0,005 кг/Л( • С2 - коэффициент аэродинамического сопротивления воздушного участка, Копр г. - Я2Р = 0,002 кг/м-с1 - коэффициент аэродинамического

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

На основе линейной модели ГВТ и с использованием уравнений коэффициентов на основе уравнения Дарси-Вейсбаха записывается нелинейная математическая модель:

йвх = - kconp.eSe.Ql ¡К. + $в.Рт/1в. + *аАЯаА -

•<2вых=- кСопр.г^г.01„/к ~ «г.Рт/Ь. + кдс$г^с1к, - (7)

гт

В основе математической модели ГВТ (7) используются уравнения, описывающие баланс между входящей и исходящей средой, соотношения между давлением и расходом воздуха или удаляемыми газами. Эти процессы можно назвать аэродинамическими. Однако, физика ГВТ подразумевает также процессы тягодутьевых машин, приводящих в движение среду. Для приведения в движение газовоздушной среды необходимо заставить вращаться лопасти вентилятора и дымососа. Эту задачу выполняет электродвигатель. Уравнения электромеханического взаимодействия между электродвигателем, лопастями тягодутьевой машины и движимой средой являются электродинамической составляющей модели ГВТ.

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

тягодутьевых установок. Изучается влияние электродинамических инерционностей на переходные процессы.

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

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

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

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

Движение двухмассовой системы описывается уравнениями:

^с1щ/сИ = М-с12((р1-(р2)-0,05Мвном ^

/2 <ко2М = с12( Ф, -Ф2)- Шв + (Мв пом + АМв )(м2/м„ „ш, )2 ' где кг-м2 - момент инерции двигателя и жестко с ним связанной первой полумуфты, С]2, Н-м - жесткость механической связи, 32, кг-м2 - момент инерции вентилятора и жестко с ним связанной второй полумуфты, М, Н-м -электромагнитный момент привода, щ, <р2 - угловые перемещения масс, й)Ь рад/с - угловая скорость вращения ротора двигателя, Н-м - момент

механических потерь вентилятора, Мв ном = Рв ном,/юв ном.= 30Рвном,/я«в ном , Н-м - номинальный момент вентилятора, Рвном., кВт - номинальная мощность вентилятора, пв ном., об/мин - номинальная скорость вращения вентилятора, со2, рад/с - угловая скорость вращения вентилятора.

Переход от двухмассовой упругой системы к одномассовой возможен при абсолютно жесткой механической связи между массами С12=оо. Для такой схемы имеет место ключевое равенство: со - со 1 = со2> рад/с - угловая скорость вращения ВА.

С учетом обозначений: М = Мс\ = Мс2, Н-м - момент сопротивления ВА,./ = J\ - У2, кг-м2 - момент инерции ВА, одномассовое уравнение движения электропривода имеет вид: М -Мс = 7й?со/Л .

Уравнение механической характеристики электродвигателя, на линейном участке, можно записать в операторном виде: (\ + Т3р)М = |3(ю0-со,), где: р -оператор дифференцирования, Тэ = 1/ш0эл.еюм.лк >с " постоянная времени, соозл.„ом., рад/с - номинальная электрическая угловая скорость, 5К - критическое скольжение электропривода, Р = 2Мк/щиом хк, Н-м-с - статическая жесткость, Мк, Н-м - критический момент электропривода, сооном.. рад/с - номинальная угловая скорость, ш0 = 2л//рп, рад/с - угловая скорость идеального холостого хода, /, Гц - частота тока электродвигателя, ра - число пар магнитных полюсов обмотки статора.

Структурная схема ГВТ с учетом электропривода изображена на схеме (Рисунок 2), где: /дв, Гц - частота тока электродвигателя (ЭД) дутьевого вентилятора (ДВ), Мдв, Н-м - электромагнитный момент привода ДВ, содв, рад/с - угловая скорость вращения ДВ, <2ВХ, кг/с - расход общего воздуха, рт, Па - разрежение в топке, бВЬ1Х, кг/с - расход уходящих газов, шдс, рад/с -угловая скорость вращения дымососа (ДС), Мдс, Н-м - электромагнитный момент привода ДС, /дс, Гц - частота тока ЭД ДС. Точками отмечены места измерения параметров (моменты приводов могут быть измерены косвенным путем).

Топка

Рисунок 2. Схема ГВТ

Объект является двухканальным. К входным (управляющим) величинам относятся - /дв, /дс, к выходным (управляемым) - <2ВХ, рт, промежуточные

величины: Мдв, Мдс, мдв, содс. На основе разработанной нелинейной модели ГВТ и уравнений электропривода можно составить систему уравнений ГВТ:

бвх -_^сопр.вАв.бвх/'в. + ^"двшдв + Sb.Pt/L.

шдв =-М^в//дв +МдВ/7дВ

^дв = -Л/двАдэв - РдвЮдвДдв + 2пРдв/дв/РдаТдв

' ¿вых = _^сопр.г.^г.2вых/'г. + £дсюдс — ^г.Рт/'г.

Юдс =-Mcдc/Jдc+MRC/JRC

Л/дС =-Мдс/гдэс -рдссодс/7дэс +2яРдС/дс/^сГдЭс

Рт=Овх/кс. - ' ^вых /К. где ^Сопр.в. > кг/м-с2 - коэффициент аэродинамического сопротивления воздушного участка, 5В, м2 - площадь поперечного сечения воздушного участка, 1В , м - длина воздушного участка, &дв, Пас - коэффициент зависимости давления ДВ от угловой скорости вращения (принимается по аэродинамической характеристике), МдВ=/(ЮдВ), Н-м - момент сопротивления ДВ (пропорционален квадрату угловой скорости вращения), УдВ, кг-м - момент инерции ДВ, 7дВ, с - электромагнитная постоянная

времени ДВ, Здв, Н-м-с - статическая жесткость ДВ, рдВ - число пар магнитных полюсов обмотки статора ДВ, £сопр г, кг/м-с2- коэффициент аэродинамического сопротивления газового участка, 5Г , м2 - площадь поперечного сечения газового участка, 1Т , м - длина газового участка, £дс, Па-с

- коэффициент зависимости давления ДС от угловой скорости вращения (принимается по аэродинамической характеристике), /дс, кг-м2 - момент

инерции ДС, Л/дС = f( шдс ), Н-м - момент сопротивления ДС (пропорционален квадрату угловой скорости вращения), ТдС, с - электромагнитная постоянная времени ДС, Рдс, Н-м-с - статическая жесткость ДС, /?дС - число пар магнитных полюсов обмотки статора ДС, кс, м-с2- коэффициент состояния газовоздушной среды (отношение массы и давления).

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

' . 2

. _ 2

х2 --а22х2 +а23х3 ^

Х3 = ^зз-^з — Л32Х2 1^1

У1=х1

где Яц =^сопр.в5в.Д. . а\2=к№' а22-У]№ ' а23 ~ У^ДВ ' Я33 ~ У^ДВ '

«32 = РдвДдв> ^31 = 2яРдв/рдв7дв ■

Можно переписать систему, преобразовав уравнения относительно выходной величины У1 к уравнению вида: у\= ^(г,х) + Ьх(1,х)иу, где: х = {у1,у[,у^}={х1,х2,х3}.

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

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

Переходные процессы в указанной области имеют апериодический характер. Корни характеристического уравнения указанной области определяют вид желаемого уравнения: 6у"+ 6+ >1 = V; •

Закон управления со старшей производной в обратной связи записывается на основе характеристического уравнения щ = кх(-6 у"- 6>>, -У\+Ух-ух). Статическую точность требуется обеспечивать на уровне Д= ( 0,05-0,0 и/С г,у).

Следовательно, коэффициент усиления регулятора к\ следует выбирать из диапазона 20-г100.

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

Щ \ихр) = Ц, р3 + А^р2 +ВЩР +1.

Для данного случая малый параметр |Д, выбирается равным 0,1, коэффициенты демпфирования (А, В) равны двум. Передаточная функция дифференцирующего фильтра выглядит следующим образом:

\¥хф( щ>) = 0,001р3 + 0,02р2 +0,2/? + 1.

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

У/^Ъ)=к1(ръ+6р2+6р + \ )/0,001 рг + 0,04 р2 + 0,2 р +1.

Возможны различные схемные реализации регуляторов. При построении САР ГВТ используется два типа схем. Одна из них это включение ПД-регулятора в обратном контуре замкнутой системы. Вторая схема это модернизация ПД-регулятора до ПИД-регулятора путем параллельного присоединения интегрирующего звена. Такая схема включения ПИД регулятора называется двухконтурной.

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

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

Проверка работы регулятора осуществляется путем моделирования САР в 8тиНпк. Коэффициенты дифференциальных уравнений рассчитываются автоматически с помощью программного приложения в МаЙаЬ.

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

ГВТ.

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

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

Таким образом, метод настройки регулятора газовоздушного тракта предполагает следующую последовательность:

1. Разработка модели ГВТ на основе физических уравнений.

1.1. Разработка линейной модели.

1.2. Усовершенствование модели с рассмотрением нелинейных и нестационарных свойств модели.

1.3. Усовершенствование модели с рассмотрением малых инерционностей. Увеличение порядка модели.

2. Синтез регулятора по методу локализации.

2.1. Синтез ПД-регулятора для линейной модели.

2.2. Проверка работы ПД-регулятора для нелинейной модели.

2.3. Синтез ПЛ-регуляторов (порядок дифференцирующего звена больше 1) для усовершенствованных моделей с малыми инерционностями.

3. Составление структурной схемы САР на основе дифференциальных уравнений модели.

4. Автоматический расчет коэффициентов дифференциальных уравнений с помощью программы.

5. Снятие переходных процессов, проверка работы регулятора.

6. Применение регулятора для множества однотипных технологических объектов (повторение пунктов 4 и 5).

Проводится описание разработки автоматической системы настройки регулятора (АСНР). Определены типы технологических объектов, с которыми может работать система. Для каждого конкретного случая рассматривается функционально-алгоритмическая схема и принципы настройки регулятора.

Особое внимание уделено произвольному объекту: рассмотрена блок-схема выбора типа регулятора и пример работы программы. Проведено сравнение АСНР с адаптивным регулятором.

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

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

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

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

Информацию удобно представлять в виде функционально-алгоримитеской схемы АСНР. В ходе исследований были получены две математические модели газовоздушного тракта: аэродинамическая и аэродинамическая с электродинамической составляющей. Для известных объектов и произвольного объекта функционально-алгоритмическая схема имеет различия.

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

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

Для аэродинамической модели ГВТ алгоритм работы АСНР следующий. Выбирая тип технологического объекта - ГВТ, программа анализирует нелинейную математическую модель ГВТ. Результатом анализа модели можно представить следующее: объект является нелинейным и нестационарным, тип

входного воздействия - ступенчатый; на объект воздействуют ступенчатые возмущения. Требования к качеству регулирования следующие: статическая ошибка А = 5%, t„ - время переходного процесса менее 10с и перерегулирование о = 5%. На основании проведенных ранее исследований по синтезу регулятора для нелинейной модели ГВТ выбирается двухконтурный ПИД-регулятор.

Алгоритм работы АСНР для аэродинамической модели ГВТ с электродинамической составляющей содержит следующие действия. В случае выбора типа технологического объекта - ГВТ с электроприводом, в программе анализируется соответствующая объекту нелинейная модель. Результатом анализа модели можно представить следующее: объект является нелинейным и нестационарным, тип входного воздействия - ступенчатый; на объект воздействуют ступенчатые возмущения. Требования к качеству регулирования следующие: статическая ошибка А = 5%, tn - время переходного процесса менее 20с и перерегулирование о = 5%. На основании проведенных ранее исследований по синтезу регулятора для модели ГВТ с электроприводом выбирается двухконтурный Р120-регулятор.

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

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

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

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

Дальнейшее развитие АСНР предполагает добавление в программу нелинейных моделей других объектов. Также могут быть проведены попытки автоматической идентификации объектов нелинейными моделями. На данном этапе АСНР может применяться при настройке регуляторов для ГВТ теплоэнергетических котлов, а также других объектов при наличии математических моделей. АСНР позволяет сравнивать системы регулирования с регуляторами, синтез которых проведен по методу локализации с системами с обычным ПИД-регулятором.

Рассматривается автоматизированный способ настройки регулятора с использованием ЭВМ.

Ручной процесс настройки ПИД-регулятора имеет следующий порядок действий:

- Установка всех коэффициентов регулятора в нулевое значение.

- Плавное увеличение пропорционального коэффициента кпр, и фиксация

значения при статической ошибке около 10%.

- Плавное увеличение времени интегрирования Ти в случае необходимости уменьшения статической ошибки.

Плавное увеличение времени дифференцирования Тд в случае

склонности системы к излишним колебаниям.

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

Использование современных технологий регулирования позволяет избежать вышеперечисленных проблем. Виртуальные модели с адаптивными регуляторами в БСАБА-системе используются для расчета коэффициентов реальных регуляторов на объекте.

Алгоритм настройки регулятора имеет следующую последовательность действий.

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

Для сглаживания функции проводится аппроксимация полиномами в приложении МаНаЬ. Таким образом, представлены возможности аппроксимации исходной функции удобными способами с помощью программы МаЙаЬ.

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

Таблица. Зависимость разрежения в топке от времени

р, Па t, с

0 0

17,5 5

22,5 10

24 15

25 20

Аппроксимация полиномом четвертой степени производится с помощью следующего кода: >>х=[0,5,10,15,20]; » у=[0,17.5,22.5,24,25]; » p=polyfit(x,y,4)

р = -0.0004 0.0240 -0.5400 5.6500 -0.0000 Полученный полином:

у = -0,0004л:4 + 0,024л;3 - 0,54л:2 + 5,65х. Вычисление значений функции от х осуществляется записью: » f=polyval(p,x)

f = -0.0000 17.5000 22.5000 24.0000 25.0000 (совпали с исходными полностью).

Программное приложение Matlab позволяет идентифицировать математическую функцию в виде передаточной функции объекта. Данные для приложения Matlab/Toolboxes/System identification формируются в приложении Simulink в блоке IDDATA SINK. На вход Input блока IDDATA SINK подается константа, которая формирует кривую, соответствующую математической функции объекта, на вход Output сама кривая непосредственно. Задание параметров производится в блоке data. Кривая задается блоком Polival, реализующим функцию полинома от линейно-нарастающей величины, реализованной блоком Integrator. В Matlab полином записывается как р= -0.0004 0.0240 -0.5400 5.6500 -0.0000.

Данные, сгенерированные в Simulink импортируются в приложение System Identification Tool: Import Data - Data Object - ввод имени «data» в поле Workspace variable Object - Import.

На основе импортированной модели данных осуществляется идентификация модели объекта в виде передаточной функции Estimate -Process models.

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

На основе полученной модели объекта с помощью SCADA-системы, например Trace Mode определяются коэффициенты настройки

ПИД-регулятора. Для этой цели проводится построение контура регулирования с добавлением блоков адаптивного ПИД-регулятора.

ЗАКЛЮЧЕНИЕ

В диссертационной работе получены следующие результаты, при разработке и исследовании модернизированного многоканального пд - регулятора для стабилизации режимов работы теплоэнергетического котла:

1) разработана линейная математическая модель ГВТ на основе физических законов протекающих процессов,

2) проведен синтез ПД-регуляторов по методу локализации для линейной модели,

3) в результате анализа нелинейных и нестационарных свойств объекта получена нелинейная математическая модель; путем добавления дополнительных инерционностей получена нелинейная математическая модель с электроприводом,

4) применен ПД-регулятор для усовершенствованных моделей, проведен синтез ПД-регуляторов более высокого порядка для усовершенствованных моделей,

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

6) разработан метод автоматизированной настройки коэффициентов ПИД-регулятора.

Математическая модель отражает такое физическое свойство объекта, как пульсации разрежения, являющиеся известной проблемой в практике.

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

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

Рассмотрен метод настройки регулятора ГВТ, объединяющий разработку модели и синтез регулятора. Важность метода в том, что объект широко распространен в теплоэнергетике и требует высокого качества регулирования.

Существует возможность применения методики для сходных технологических объектов в других отраслях, например, горнодобывающей промышленности. Полученный регулятор учитывает особенности ГВТ, такие как нестационарность и нелинейность.

Разработанная программа для автоматического расчета коэффициентов АСНР позволяет работать со всем множеством моделей ГВТ. Полезным свойством программы является автоматический перенос рассчитанных коэффициентов на заранее построенную структурную схему САР в Simulink. Регуляторы АСНР имеют возможность отрабатывать «нестационарности» объекта без изменения настройки регулятора.

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

Публикации по теме работы Статьи в журналах, рекомендованных ВАК для публикации результатов диссертации на соискание ученой степени доктора и кандидата наук

1. Пономарев A.A. Исследование двухконтурного ПИД-регулятора газовоздушного тракта теплоэнергетического котла на основе метода локализации // Научный вестник НГТУ - 2012. № 4(49). - С. 191-196.

2. Вострнков A.C., Пономарев A.A. Метод настройки регулятора газовоздушного тракта на основе принципа локализации // Автометрия - 2012. Т. 46, № 2. - С. 58-66.

3. Востриков A.C., Пономарев A.A. Синтез двухканального регулятора газовоздушного тракта теплоэнергетического котла // Доклады академии наук ВШ РФ - 2011. № 1 (16). - С. 95-105.

Статьи в научных сборниках и периодических научных изданиях

4. Пономарев A.A. Метод настройки регулятора газовоздушного тракта теплоэнергетического котла // Сборник научных трудов НГТУ - 2012. -№ 1(67). - С. 57-68.

5. Пономарев A.A. Модель газовоздушного тракта теплоэнергетического котла как объекта регулирования // Сборник научных трудов НГТУ - 2010. -№3(61).-С. 19-29.

Публикации в материалах научных мероприятий (международные конференции)

6. Пономарев A.A. Разработка автоматической системы настройки регулятора // Материалы международной научно-технической конференции студентов, аспирантов и молодых ученых "Прогрессивные направления развития машино-приборостроительных отраслей и транспорта". -Севастополь: СевНТУ, 2012. - с. 29-31.

7. Пономарев A.A. Разработка автоматической системы настройки регулятора // Современные инструментальные системы, информационные технологии и инновации. Материалы IX-ой Международной научно-практической конференции. - Курск: Юго-Зап. гос. ун-т, 2012. - с. 172-180.

8. Пономарев A.A. Разработка и исследование двухканальной системы автоматического регулирования газовоздушного тракта теплоэнергетического

котла // Материалы VI Международной молодежной научной конференции «Тинчуринские чтения». - Казань: Казан, гос. энерг. ун-т, 2011. - С. 14-15. 9. Vostrikov A.S. Ponomarev А.А. Research and synthesis of control for seven degree air-gas circuit model // Proceedings of the 2-nd Indo-Russian Joint Workshop on Computing Intelligence and Modern Heuristics in Automation and Robotics. -Novosibirsk: NSTU, 2011. - P. 42-45. [Исследование и синтез регулятора для модели газовоздушного тракта седьмого порядка]

Отпечатано в типографии Новосибирского государственного технического университета 630073, г. Новосибирск, пр. К. Маркса, 20, тел./факс (383) 346-08-57

формат 60x84/16, объем 1.75 п. л., тираж 100 экз., заказ № 431, подписано в печать 12.03.2013 г.

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ

УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ «НОВОСИБИРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ

УНИВЕРСИТЕТ»

РАЗРАБОТКА И ИССЛЕДОВАНИЕ МОДЕРНИЗИРОВАННОГО МНОГОКАНАЛЬНОГО ПД - РЕГУЛЯТОРА ДЛЯ СТАБИЛИЗАЦИИ РЕЖИМОВ РАБОТЫ

!

ТЕПЛОЭНЕРГЕТИЧЕСКОГО КОТЛА

Специальность: 05.13.01 - Системный анализ, управление и обработка

информации (промышленность)

На правах рукописи

ПОНОМАРЕВ Александр Александрович

Диссертация на соискание учёной степени кандидата технических наук

Научный руководитель д.т.н., проф. Востриков А.С.

Новосибирск - 2012

СОДЕРЖАНИЕ

СПИСОК СОКРАЩЕНИЙ.....................................................................................3

ВВЕДЕНИЕ.................................................................................................................4

1 РАЗРАБОТКА ЛИНЕЙНОЙ МОДЕЛИ ГАЗОВОЗДУШНОГО ТРАКТА (ГВТ) И СИНТЕЗ ПД - РЕГУЛЯТОРОВ ПО МЕТОДУ ЛОКАЛИЗАЦИИ 12

1.1 Линейная модель ГВТ.......................................................................12

1.2 Синтез ПД-регуляторов для модели................................................24

1.3 Постановка задач для исследования................................................36

2 УСОВЕРШЕНСТВОВАНИЕ МОДЕЛИ ГВТ.............................................37

2.1 Нелинейная модель ГВТ и синтез ПД-регуляторов......................37

2.2 Модель ГВТ с электроприводом и синтез двухконтурных ПИД-регуляторов ............................................................................................................50

2.3 Выводы по главе................................................................................69

3 ИССЛЕДОВАНИЕ МЕТОДОВ НАСТРОЙКИ РЕГУЛЯТОРОВ ГВТ.. 70

3.1 Методика настройки регулятора ГВТ на основе принципа локализации ............................................................................................................70

3.2 Автоматическая система настройки регулятора (АСНР)..............80

3.3 Автоматизированный способ настройки регулятора с использованием ЭВМ................................................................................................92

3.4 Выводы по главе..............................................................................106

ЗАКЛЮЧЕНИЕ.....................................................................................................107

СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ.........................................109

ПРИЛОЖЕНИЕ Свидетельство о государственной регистрации программы для эвм и акт о внедрении.............................................................120

СПИСОК СОКРАЩЕНИЙ

ПИД-регулятор - пропорционально-интегрально-дифференциальный регулятор.

ПД-регулятор - пропорционально-дифференциальный регулятор.

ПЛ-регулятор - регулятор, синтез которого проведен по методу локализации.

ГВТ - газовоздушный тракт.

ДВ - дутьевой вентилятор.

ДС - дымосос.

ЭД - электродвигатель.

КА - котельный агрегат.

ВА - вентиляторный агрегат.

ТДУ — тягодутьевая установка.

САР - система автоматического регулирования.

АСНР - автоматическая система настройки регулятора.

ВВЕДЕНИЕ

Актуальность проблемы. Для качественного ведения технологического процесса выработки пара необходимо автоматическое регулирование параметров. Улучшение поддержания заданных параметров осуществляется благодаря синтезу систем регулирования с усовершенствованными законами регулирования. Модернизированные законы регулирования позволяют устанавливать высокие требования к показателям качества.

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

Метод локализации предлагает теоретически обоснованную, регулярную методику синтеза систем управления для нелинейных и нестационарных объектов. Исследование является перспективным для развития метода локализации. Область применения метода локализации распространяется на теплоэнергетические объекты. Совершенствуется его практическое применение и увеличивается прикладное значение.

Степень изученности. В настоящее время проблема синтеза регуляторов газовоздушных трактов является многогранной и рассматривается с различных сторон. Предпринимаются попытки адекватно идентифицировать технологический объект математической моделью с целью построения точных структурных схем систем регулирования и оптимального управления [62]. Имеет место практическая идентификация объекта с помощью прикладного программного обеспечения. Другим направлением исследований является изучение применения современных технических средств, таких как контроллеры и частотные преобразователи [8]. Изобретательские и

нестандартные подходы включают в себя комбинирование различных методов регулирования параметров. Например, изучается одновременное применение регулирующих заслонок и частотных преобразователей для регулирования параметров газовоздушного тракта [63].

Объект исследований. Теплоэнергетическое оборудование, входящее в состав современных тепловых электрических станций необходимо оснащать системами автоматического регулирования высокого качества и надежности. Отдельно можно выделить газоводушный тракт, являющийся неотъемлемой частью любого котла. Согласно исследованиям, проведенным в США, при увеличении эффективности работы газовоздушного тракта может происходить энергосбережение на уровне 2-20%. С точки зрения теории регулирования объект представляет особый интерес ввиду его многоканальности и нелинейных и нестационарных свойств.

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

- разработать математическую модель на основе физических законов протекающих процессов в топке котла,

- провести синтез регуляторов параметров объекта на основе метода локализации,

- усовершенствовать математическую модель путем добавления дополнительных инерционностей, рассмотрения нелинейных и нестационарных свойств объекта,

- применить регулятор для усовершенствованной модели, проверить его работоспособность,

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

Методы исследований. Математические модели строятся аналитическим методом с использованием уравнений, описывающих физические свойства объекта. Синтез системы регулирования проводится на основе метода локализации с использованием различных схем включения регуляторов. Эксперименты над системой и ее последующая корректировка проводится в пакете прикладных программ 81гпиНпк МаНаЬ. Разработка метода настройки регулятора, не требующего перенастройки в ходе работы, включает применение БшидИпк МайаЬ для идентификации объекта на основе экспериментальных данных. Синтез адаптивного регулятора и проверка его работы для экспериментальных моделей проводится с помощью 8САОА-системы.

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

Научная новизна. Разработаны новые математические модели газовоздушного тракта аналитическим методом. К ним относятся линейная, нелинейная и нелинейная с электроприводом. Каждая из моделей может применятся для синтеза систем регулирования в зависимости от необходимой точности описания объекта. Рассмотрено применение метода локализации к теплоэнергетическим объектам, разработаны системы автоматического регулирования на его основе, вследствие чего расширена область прикладного применения метода. Новые модели, регуляторы и результаты экспериментов позволили создать метод настройки регуляторов газовоздушного тракта на основе принципа локализации. В диссертационной работе описано создание нового способа автоматической настройки регулятора (АСНР) на основе принципа локализации, представляющего собой компьютерную программу. АСНР используется для быстрой настройки регуляторов. В рамках исследований предложен новый алгоритм настройки регуляторов для различных объектов с использованием электронно-вычислительной техники.

Основные положения, выносимые на защиту:

- разработана новая линейная модель газовоздушного тракта (ГВТ), для которой проведен синтез ПД-регуляторов по методу локализации,

- разработаны новые нелинейные модели ГВТ и ГВТ с электроприводом, для которых проведен синтез ПД-регуляторов по методу локализации,

- разработан метод настройки регулятора газовоздушного тракта на основе принципа локализации,

- разработана программа автоматической настройки регулятора,

- разработан метод автоматизированной настройки регулятора с использованием ЭВМ.

Практическая ценность. Результаты исследований данной работы могут быть применены для реальных технологических объектов. Математическими моделями могут описываться множества типовых газовоздушных трактов с различными техническими характеристиками. Для этого в модели должны быть перерасчитаны числовые коэффициенты. Применение регулятора, синтез которого проведен по методу локализации, также возможно для множества газовоздушных трактов с проверкой работоспособности и корректировкой коэффициентов в случае необходимости. Метод настройки регулятора газовоздушного тракта опробован для теплоэнергетического котла ТПЕ-214, установленного на Новосибирской ТЭЦ-5. Автоматизированный метод настройки с использованием ЭВМ применим для технологических объектов, которые требуют регулирования параметров. Метод используется в работе фирмы ЗАО «СИНТЭП».

Внедрение результатов работы. Результаты диссертационной работы использованы в производственной деятельности ЗАО «СИНТЭП»: при разработке программного обеспечения, наладке оборудования. Использование указанных результатов позволяет повысить эффективность работы технологического оборудования.

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

Апробация работы. Основные положения диссертационной работы и отдельные ее результаты докладывались на: VI Международной молодежной научной конференции «Тинчуринские чтения» (г. Казань, КГТУ, 2011 г.), 2-nd Indo-Russian Joint Workshop on Computing Intelligence and Modern Heuristics in Automation and Robotics (Novosibirsk, NSTU, 2011), Конференции «Системный анализ и информационные технологии» (г. Новосибирск, НГТУ, 2011 г.), 1Х-ой Международной научно-практической конференции «Современные инструментальные системы, информационные технологии и инновации» (г. Курск, Юго-Зап. гос. ун-т, 2012 г.), Международной научно-технической конференции студентов, аспирантов и молодых ученых "Прогрессивные направления развития машино-приборостроительных отраслей и транспорта" (г. Севастополь, СевНТУ, 2012 г.).

Публикации. Основные положения и результаты диссертационной работы опубликованы в 9 работах, 3 из которых входят в перечень рецензируемых научных журналов и изданий для опубликования основных научных результатов диссертаций, 2 статьи в сборниках научных трудов и 4 статьи в сборниках материалов международных конференций.

Структура и объем работы. Диссертация состоит из введения, трех глав, заключения, списка использованной литературы из 114 наименований и приложения. Общий объем работы составляет 122 страницы, из которых основное содержание работы - 119 страниц, включая 11 таблиц и 53 рисунка.

Благодарность. Автор выражает благодарность научному руководителю - заслуженному деятелю науки РФ, доктору технических наук, профессору Вострикову Анатолию Сергеевичу за вклад в виде руководства работой, практических советов и консультаций по теме диссертации.

В первой главе описывается разработка линейной модели ГВТ и синтез ПД-регуляторов по методу локализации.

Построение модели ГВТ проводится с последовательным приближением к наиболее точному математическому описанию. Такой способ позволяет проводить проверку соответствия модели реальному объекту по всем его

свойствам. Основные методы проверки это построение переходных процессов и сравнение с разгонными характеристиками технологического объекта, а также анализ и синтез САР.

Линейная модель это первый этап математического описания ГВТ. Модель объекта составляется на основании законов физики.

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

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

Допущением линейной модели является то, что потоки воздуха и газов ламинарные. Однако более точно поток газовоздушной среды характеризуется уравнением Дарси-Вейсбаха (соотношением между скоростью движения газа и потерями давления на трение для газопровода конечной длины.

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

Для модели проводится процедура синтеза регулятора третьего порядка по методу локализации. Проводятся необходимые эксперименты, делаются

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

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

гвт.

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

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

Проводится описание разработки автоматической системы настройки регулятора (АСНР). Определены типы технологических объектов, с которыми может работать система. Для каждого конкретного случая рассматривается функционально-алгоритмическая схема и принципы настройки регулятора. Особое внимание уделено произвольному объекту: рассмотрена блок-схема выбора типа регулятора и пример работы программы. Проведено сравнение АСНР с адаптивным регулятором.

Рассматривается автоматизированный способ настройки регулятора с использованием ЭВМ.

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

исследовании модернизированного многоканального пд - регулятора для стабилизации режимов работы теплоэнергетического котла.

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

1 РАЗРАБОТКА ЛИНЕЙНОЙ МОДЕЛИ ГАЗОВОЗДУШНОГО ТРАКТА (ГВТ) И СИНТЕЗ ПД - РЕГУЛЯТОРОВ ПО МЕТОДУ ЛОКАЛИЗАЦИИ

1.1 Линейная модель ГВТ

Термин математическая модель [70] определяется как набор уравнений, описывающих поведение физической системы. Поскольку невозможно учесть все реальные свойства объекта, модель не может являться абсолютно точной. При составлении модели необходимо находить баланс между адеква�