УДК 621.34: 62-50
Исследование применения нечетких регуляторов в системах управления технологическими процессами
Куленко М.С., канд. техн. наук, Буренин С.В., ст. преп.
Рассмотрены вопросы использования регуляторов на основе нечеткой логики в системах автоматического управления различной конфигурации. Приведены результаты исследования систем с нечеткими регуляторами. Эффективность использования нечетких регуляторов показана на примере построения системы нечеткого регулирования натяжения транспортируемого гибкого материала.
Ключевые слова: нечеткая логика, термы, фаззификация, пропорционально-интегральный регулятор.
Researching Application of Fuzzy Controllers in Technological Processes Control Systems
M.S. Kulenko, Candidate of Engineering, S.V. Burenin, Senior Teacher
The author considers questions of using the controllers based on fuzzy logic in automatic control systems of different configurations. The author presents the research results of fuzzy controllers systems. The efficiency of using fuzzy controllers is shown by the example of system integration of fuzzy regulation of pulling transported flexible material.
Key words: fuzzy logic, terms, fuzzification, proportional and floating controller.
Свойства нечеткой логики обрабатывать неполную информацию, моделировать человеческие знания и выдавать обоснованные решения предполагают ее интенсивное использование для наблюдения в реальном времени за технологическими процессами, а также решение проблем, связанных с практической реализацией систем управления технологическими процессами [1]. Использование нечеткой логики в системах управления позволяет уменьшить вмешательство оператора в процесс управления и, следовательно, позволяет разработать новые методики управления, более адаптированные к промышленной среде.
Регуляторы, построенные на базе нечеткой логики, в ряде случаев способны обеспечить более высокие показатели качества переходных процессов по сравнению с классическими регуляторами. Кроме того, используя методы синтеза нечетких алгоритмов управления, можно выполнить оптимизацию сложных контуров регулирования без проведения всесторонних математических исследований.
Исследуем поведение систем на основе элементарных динамических звеньев. Выберем в качестве простейшего объекта регулирования инерционное звено первого порядка.
Так как объект не содержит интегрирующих звеньев, то можно провести оптимизацию модуля передаточной функции контура регулирования.
Рассмотрим применение пропорционально-интегрального (ПИ) регулятора. Структурная схема модели объекта с одним инерционным звеном 1-го порядка с классическим и нечетким ПИ-регулятором, выполненная в про-
граммном комплексе МаНаЬ/Бітиііпк, представлена на рис. 1.
Нечеткая система выполнена по типу Мамдани [2] с двумя входами: для пропорциональной х1 и интегральной х2 составляющих. Для лингвистического описания каждой входной переменной выбраны семь треугольных термов (ЫВ, ЫМ, УБ, 2Б, РБ, РМ, РВ) и симметричные диапазоны изменения.
Рис. 1. Структурная схема объекта, представленного инерционным звеном 1-го порядка с ПИ-регулятором и с ПИ-подобным (І^у-регулятором
Лингвистические правила для ПИ-подобного Ішгу-регулятора приведены в табл. 1.
Таблица 1
Х2 N Z P
N NB NM ZE
Z NM ZE PM
P ZE PM PB
Правила формируются по типу ЕСЛИ ... И ..., ТО.... Переходные процессы при ступенчатом изменении задания на входе и реакции на возмущающее воздействие на объект представлены на рис. 2.
' ■ ■
1
^ а): :
г
Рис. 2. Переходный процесс для системы с инерционным звеном 1-го порядка: а - с ПИ-регулятором; б - с ПИ-подобным (І^у-регулятором
Исследуем применение пропорционально-интегрально-дифференциального (ПИД)
регулятора. Структурная схема модели объекта с одним инерционным звеном 1-го порядка с ПИД-подобным Иггу-регулятором представлена на рис. 3. Здесь в структурную схему модели добавлена дифференциальная составляющая сигнала рассогласования.
Рис. 3. Структурная схема объекта, представленного одним инерционным звеном 1-го порядка с ПИД-регулятором и ПИД-подобным (шгу-регулятором
Список логико-лингвистических правил дополнен следующими тремя правилами:
1) IF (x3=N) THEN (y=N);
2) IF (x3=ZE) THEN (y=ZE) ;
3) IF (x3=P) THEN (y=P).
Действительно, если скорость изменения
ошибки положительна, то выходной сигнал регулятора должен увеличиться, чтобы предупредить увеличение сигнала рассогласования, и наоборот. Более точная настройка диапазона изменения Д-составляющей осуществляется методом подбора значений.
Переходные процессы для модели с аналоговым ПИД-регулятором и ПИД-подобным fuzzy-регулятором представлены на рис. 4.
"I----------------------------------------------------------------1-------------------------------------------------------------1------------------------------------------------------------1------------------------------------------------------------1-------------------------------------------------------------1------------------------------------------------------------1-------------------------------------------------------------г
: : !/«;
: : : Г7
; : ;
О 10 20 30 40 50 60 70 80 90 100
Рис. 4. Переходный процесс для системы, представленной одним инерционным звеном 1-го порядка: а - с ПИД-регулятором; б - с ПИД-подобным ^гу-регулятором
Наличие Д-составляющей уменьшило динамическую ошибку системы, однако время регулирования в случае использования ПИД-подобного (иггу-регулятора увеличилось.
Предположим теперь, что объект регулирования содержит два инерционных звена первого порядка с одинаковыми постоянными времени. Исследуем применение ПИ-
регулятора.
Лингвистические правила в этом случае приведены в табл.2.
Таблица 2
X2 N Z P
N NB PM NB
Z NM ZE PM
P ZE PM PB
Переходные процессы для модели с аналоговым регулятором и модели с ІІггу-регулятором при ступенчатом изменении задания на входе и отработке единичного возмущающего воздействия на объект представлены на рис. 5.
Рис. 5. Переходный процесс для системы, представленной двумя инерционным звеньями 1-го порядка: а - с ПИ-регулятором; б - с ПИ-подобным (шгу-регулятором
В случае применения ПИД-подобного (иггу-регулятора для объекта с двумя инерционными звеньями 1-го порядка получим графики переходных процессов, представленные на рис. 6.
6)
1 ЙЯ
.. г ..../ у.
\ ... ,| ....
А*
й2 :/ : и : И 7
10 20 ЗЭ 40 S ' і , » , 0
Рис. 6. Переходный процесс для системы, представленной двумя инерционными звеньями 1-го порядка: а - с ПИД-регулятором; б - с ПИД-подобным (икгу-регулятором
Очевидно, время регулирования и число колебаний больше для системы с (иггу-регулятором, однако динамическая ошибка значительно уменьшилась.
Предположим, что объект регулирования содержит не только инерционное звено первого порядка, но и интегрирующее звено. Оптимизацию контура регулирования с объектом, содержащим интегрирующее звено, осуществим по симметричному оптимуму.
Создадим нечеткую систему, реализующую типовой аналоговый ПИ-регулятор.
Лингвистические правила для такого ПИ-подобного Ішгу-регулятора приведены в табл. 3.
Таблица З
N Z P
N NB PM NB
Z NM ZE PM
P ZE PM PB
Переходные процессы для модели с аналоговым ПИ-регулятором и модели с ПИ-подобным fuzzy-регулятором представлены на рис. 7.
'О 10 20 30 40 50 60 70 S0 90 100
Рис. 7. Переходный процесс для системы, представленной одним инерционным звеном 1-го порядка и интегрирующим звеном: а - с ПИ-регулятором; б - с ПИ-подобным fuzzy-регулятором
Очевидно, колебательность несколько выше в системе с fuzzy-регулятором, однако динамическая ошибка значительно меньше.
Результаты реализации ПИД-подобного fuzzy-регулятора представлены на рис. 8.
1.Ь 1 0.5 0 ■0.5 •1 1 !
Г~а)
10 20 30 40 50 ео 70 SO 90 100
Рис. В. Переходный процесс для системы, представленной одним инерционным звеном 1-го порядка и интегрирующим звеном: а - с ПИД-регулятором; б - с ПИД-подобным fuzzy-регулятором
По аналогичному алгоритму выполнен анализ использования нечетких регуляторов для объектов, описываемых комбинацией элементарных звеньев: интегральных и апериодических с большими и малыми постоянными времени.
В качестве примера реализации нечеткого регулятора рассмотрим использование fuzzy-регулятора в системах управления натяжением гибкого транспортируемого материала.
В большинстве случаев при проектировании систем управления натяжением [З] требуется модель, отражающая лишь наиболее общие характеристики системы и отличающаяся достаточной для анализа точностью при ее максимальной простоте. В практике проектирования автоматических систем регулирования натяжения, построенных на базе петлеоб-разователей, традиционно используются ско-
ростная и метрическая [4] модели зоны деформации. В основу их положено известное дифференциальное уравнение деформации материала в петлеобразователе:
. 1 + 8п О.
-~л ЗТ = Vр2 - Ур1 + Ур2е0 - ^р181 ± кг^к ,
I + &1 ОТ
где - - длина материала в зоне деформации; е1, е0 - относительное удлинение в зоне деформации и на ее входе соответственно; ур1, ур2 - линейные скорости материала на входе и выходе зоны; ук - линейная скорость перемещения измерительного ролика; кг - коэффициент, учитывающий геометрию охвата измерительного ролика ветвями материала.
В скоростной модели относительное удлинение вычисляется непосредственно по скоростям ур1 и vр2, в метрической - на основе длин нерастянутого материала на входе и выходе зоны деформации. Возьмем за основу скоростную модель петлеобразователя (рис. 9) как наиболее отвечающую цели управления с точки зрения формирования управляющего воздействия по каналу регулирования скорости ведомого двигателя.
Рис. 9. Структурная схема петлеобразователя
Определение входных и выходных переменных системы, так же как и правил нечеткого регулятора, основано на практических знаниях о системе и на интуиции. Поскольку регулируемой величиной является положение измерительного ролика петлеобразователя, то целесообразно выбрать одной из входных переменных ошибку по положению ролика:
_ Лзад hр,
где Лр - текущее положение ролика; Лзад - заданное положение. Причем ошибка может принимать как положительные, так и отрицательные значения.
Регулирование положения ролика петлеобразователя происходит путем изменения рассогласования линейных скоростей движения материала в зоне за счет изменения скорости вращения электродвигателя ведомого канала, причем величина рассогласования зависит от положения ролика: чтобы уменьшить натяжение материала, рассогласование следует снизить, а с другой стороны, если натяжение находится на заданном уровне, рассогласование также должно соответствовать заданному значению.
Таким образом, выходной переменной регулятора имеет смысл назначить сигнал управления скоростью ведомого электропривода, который и определяет рассогласование линейных скоростей приводных точек, а еще
З
одной дополнительной входной величиной будем считать скорость перемещения измерительного ролика или, соответственно, скорость изменения ошибки регулирования положения, что эквивалентно использованию пропорционально-интегрального регулятора.
Структурная схема системы управления с Иггу-регулятором положения ролика представлена на рис. 10.
Рис. 10. Структурная схема системы управления с ^гу-регулятором положения
Петлеобразователь представлен редуцированной передаточной функцией [4]:
4(S) =
kn
Канал ведомого электропривода представляем, также пренебрегая старшими производными в оптимизированной структуре:
Нэп($ ) =:
k
ЭП
ТэпЯ +1
Рассмотрим систему с Иггу-регулятором положения и систему, использующую аналоговый ПИ-регулятор при работе на объект, содержащий два апериодических звена с существенно различными постоянными времени.
При разработке Иггу-регулятора количественное задание параметров регулятора на его входе преобразуем в коррекцию диапазона изменения его входных координат: диапазона изменения ошибки положения ролика х1_рг (соответствующей пропорциональной составляющей регулятора) и диапазона возможных значений скорости изменения ошибки положения ролика х2_іП (соответствующей интегральной составляющей).
Формирование базы правил ІІггу-регулятора ведется на основе уравнения для линейного непрерывного ПИ-регулятора
y(t) = ^ - s(t) + -I -Je(t)dt,
которое можно заменить близким по логике управления (иггу-регулятором, если в качестве его выходной переменной рассматривать приращение управляющего воздействия Ду. Тогда закон регулирования ПИ-регулятора можно представить в дифференциальной форме Оу (Т)
dt
. d s(t) 1 ...
= kn—— + —s(t)
dt
Ти
или в разностной форме
Ду(k) = y(k) - y(k -1) = ^ Дєik) -
Дt
Таким образом, для входных переменных е(к) и Де(к) и выходной Ду(к) может быть синтезирован (иггу-регулятор, реализующий нелинейный закон
ду (к) = Р [Де(к), е(к)],
эквивалентный в определенном смысле аналоговому ПИ-регулятору.
Для нашего случая х1_рг = ДЛ соответствует сигналу рассогласования е(к), х2_П = ДЛ/ДТ соответствует приращению сигнала рассогласования Де(к), а у соответствует ДЛ(к). Лингвистические правила для такого ПИ-подобного (иггу-регулятора приведены в табл. 4.
Таблица 4
Дії NB NM NS ZE PS PM PB
^/At
NB NB NB NB NB NM NS ZE
NM NB NB NB NM NS ZE PS
NS NB NB NM NS ZE PS PM
ZE NB NM NS ZE PS PM PB
PS NM NS ZE PS PM PB PB
PM NS ZE PS PM PB PB PB
PB ZE PS PM PB PB PB PB
Структурная схема математической модели системы в комплексе МаНаЬ представлена на рис. 11.
t
ЧЕЬ
-m—щ—ДЬпі
Рис. 11. Структурная схема математической модели системы управления натяжением с нечетким регулятором
На рис. 12 представлена реакция системы с fuzzy-регулятором и аналоговым ПИ-регулятором на единичное задание положения. В момент времени t = б с происходит отработка возмущающего воздействия в виде увеличения скорости ведущего электропривода.
Следует отметить, что перерегулирования для обоих типов регуляторов примерно одинаковы, а по быстродействию нечеткий регулятор более эффективен. При отработке возмущения по каналу скорости ведущего электропривода очевидна более быстрая и качественная реакция нечеткого регулятора.
На рис.1З представлена реакция системы с fuzzy-регулятором и аналоговым ПИ-регулятором на единичное задание положения при уменьшенном на 50 % коэффициенте усиления объекта регулирования.
s(k).
и
Рис. 12. Реакция системы на единичное задание положения и возмущение по каналу скорости ведущего двигателя: 1 - ПИ-регулятор; 2 - ПИ-подобный (гкгу-регулятор
ния при двукратной вариации коэффициента усиления объекта регулирования: 1 - ПИ-регулятор; 2 - ПИ-подобный (шгу-регулятор
В результате проведенных исследований применения нечетких регуляторов при управлении объектами с различными динамическими свойствами получены рекомендации по формированию базы нечетких правил в зависимости от конфигурации систем с ПИ- и ПИД-подобными регуляторами для обеспечения характеристик, не уступающих характеристикам классических регуляторов.
Применение fuzzy-регулятора положения измерительного ролика датчика натяжения, эффективность использования которого показана на примере системы взаимосвязанного электропривода с нечетким регулятором натяжения для транспортирования гибкого длинномерного материала, позволило стабилизировать регулируемую координату в широком диапазоне вариаций параметров объекта.
Список литературы
1. Заде Л.А. Понятие лингвистической переменной и его применение к принятию приближенных решений. -М.: Мир, 197б.
2. Takagi T., Sugeno M. Fuzzy identification of systems and its applications to modeling and control // IEEE Trans. Systems Man Cybernet. - 1985. - Vol. 15. - №. 11б. - P. 11б-1З2.
3. Глазунов В.Ф., Тарарыкин С.В., Спичков Ю.П. О рациональном построении датчика натяжения ткани в поточных линиях текстильной промышленности // Изв. вузов. Технология текстильной промышленности. - 1981. -№ 1. - С. 78-82.
4. Глазунов В.Ф., Куленко М.С. Редуцирование математических моделей петлеобразователей в системе управления натяжением ткани // Изв. вузов. Технология текстильной промышленности. - 2008. - № б. - С.101-104.
Заключение
Куленко Михаил Сергеевич,
Ивановский государственный энергетический университет,
кандидат технических наук, доцент кафедры электропривода и автоматизации промышленных установок, телефон 8-920-670-87-91, e-mail: [email protected]
Буренин Сергей Владимирович,
Ивановский институт ГПС МЧС России,
старший преподаватель кафедры высшей математики и информатики, e-mail: [email protected]
б