CC BY
41
УДК 621.867.2-83
Е.С. Ребенков
СИНТЕЗ НЕЧЕТКОГО РЕГУЛЯТОРА ЧАСТОТНО-УПРАВЛЯЕМОГО ЭЛЕКТРОПРИВОДА КОНВЕЙЕРА С УПРУГИМИ СВЯЗЯМИ
Получена в векторно-матричной форме передаточная функция электромеханической системы частотно-управляемого двухдви-гательного электропривода конвейерной линии. Разработана математическая модель частотно-управляемого электропривода, учитывающая распределенную нагрузку на конвейерной линии. Показано возникновение в механической части конвейера колебательных явлений, обусловленных наличием упругости рабочей части конвейера. Получены переходные процессы в электроприводе при управляющих и возмущающих воздействиях на электропривод, выполненный по схеме «частотно-управляемый асинхронный двигатель - преобразователь частоты». Синтезирован нечеткий регулятор скорости перемещения конвейера, демпфирующий колебания скоростей масс и деформаций в упругих связях. Получены переходные процессы координат привода с нечетким регулятором. Проведено сравнение переходных процессов в электроприводе конвейера с нечетким регулятором и без него. Показана эффективность демпфирования упругих колебаний скоростей элементарных масс конвейера и деформаций его упругих связей. Ключевые слова: нечеткий регулятор, многомассовая модель привода конвейера, упругие связи, демпфирование колебаний, фиксатор нулевого порядка, функции принадлежности, универсальное множество, база правил, лингвистические оценки.
Кинематическая схема привода конвейера представлена в виде замкнутой последовательности сосредоточенных масс, соединенных упругими связями [1]. Перемещение конвейерной линии производится частотно-регулируемыми асинхронными двигателями (АД) на концевых барабанах. При управляющих воздействиях u1 и u2 на преобразователях частоты (ПЧ) в механической части привода возникают динамические колебания скоростей га. элементарных масс и деформаций 5.. упругих
ISSN 0236-1493. Горный информационно-аналитический бюллетень. 2016. № 4. С. 68-75. © 2016. Е.С. Ребенков.
связей, снижающие ресурс и надежность электропривода конвейера [2], [3].
Математическая модель механической части привода конвейера представлена в виде уравнения передаточной функции (ПФ) в векторно-матричной форме [4], [5]:
Ш(5) = М(з) •
(1)
где
Ш (5) = [юДз), ю2(з),..., юп (5)] — вектор приведенных скоростей элементарных масс;
М(з) = [ М1 (5), М2 (5),..., Мп (5)] — вектор моментов сил и моментов сопротивлений на массах;
5 (5), Х(з) — матрицы ПФ п-го порядка, где п —количест
элементарных масс.
Матрицы ПФ конвейера имеют вид:
' т11 а 12 0 .. 0 0 0 ... 0 ап1
а21 т22 а23 .. 0 0 0 ... 0 0
Х(5) = 0 0 0 .. ак,к-1 тии а к,к+1 ... 0 0
0 0 0 .. 0 0 0 ... тп-1,п-1 ап-1,п
ч а1п 0 0 .. 0 0 0 ... ап,п-1 тпп
' т11 а12 0 ... 0 0 0 ... 0 ап1
а21 т22 0 ... 0 0 0 ... 0 0
5 (5) = 0 0 5 ... ак,к-1 ткк ак,к+1 ... 0 0
0 0 0 ... 0 0 0 ... тп-1,п-1 ап-\,п
ч а1п 0 0 ... 0 0 0 ... ап,п-1 т пп
Элементы матриц определяются параметрами механической части конвейера:
-(са + V) , ткк = Зк52 - £ ац, (2)
ац =-1
где с. — жесткость упругих связей, Нм; Ь — вязкость упругих связей, Нс/м; Jk — приведенный момент инерции элементарной сосредоточенной массы, кгм2.
Двухдвигательный регулируемый электропривод конвейера выполнен по схеме «асинхронный двигатель — преобразователь частоты» (АД — ПЧ), в которой напряжение управления ПЧ иу, скорости ротора АД ю с учетом редуктора и вращающие моменты М1(к) двигателей связаны ПФ вида [6]:
иу(з) • и(з) = ют(з) • F(s) + Мт(з) • N(s). (3)
В операторных соотношениях ПФ (3) заданы параметры системы АД — ПЧ:
и(s) = вдкп, Г(s) = вд(Тз +1), Ы(s) = (Т +1)(Тs + 1),
где р„ — динамическая жесткость механической характеристики АД, Нмс; кП — коэффициент усиления по каналу напряжения и частоты ПЧ, 1/Вс; Т1 и Т2 — постоянные времени соответственно АД и ПЧ, с.
При к = 5 и п = 8 совместное решение ПФ механической части (1) и электропривода АД — ПЧ (3) дает ПФ скоростей масс в зависимости от сигнала управления ПЧ:
^ = иу ^ (4)
где полиномы ПФ (4) выражены через параметры системы:
ТО = N(5) • Х(з)[F(s) • Я11(з) + F(s) • Я55(з) + Ы(з) • Х(з)] + +Г2(з) [Я11(з) • Я55(з) - Я15(з) • Я51(з)] .
(5)
Г(з)
[Яц(з) • Яа(з) + Яй(з) • Я^(з) -
Г (з)
е1и) = ии)\ [—Я55(з) • яй(з) -яе(з) • яп(з) ] [. (6)
-N(5) • Х(з)[Яй(з) + Яа(з)]
ПФ деформаций упругой связи между /-ым и у-тым элементарными массами определяется как интеграл от разности их
скоростей:
(3) = иу (3) ^
ю
(з) -ю- (з)
з • К(з) '
Исследования на модели показали, что максимальные значения динамических колебаний скоростей элементарных масс и деформаций упругих связей возникают в местах соединений
-10
\ \
/ /
10
2В
30
■с
1,С
Рис. 1. Переходная функция скоро- Рис. 2. Переходная функция дефор-сти массы 12 мации 512
с приводными барабанами, поэтому обратные связи в регулятор введены по скоростям га2 и гак1. Переходные функции скорости (рис. 1) и деформации (рис. 2) на грузовых участках конвейера, примыкающих к приводным барабанам построены по рассчитанным ПФ для параметров механической части конвейера и электропривода:
/1 = /5 = 1,163 кг-м2, /2 = /3 = /4 = 0,567 кг-м2, /6 = /7 = /8 = 0,113 кг-м2, жесткости и вязкости на всех участках приняты равными соответственно
ci¡ = 0,126 Нм, Ьу = 0,132 Нс/м, РД = 63,7 Нмс, ' ^ = 15,7 1/Вс, Т1 = 0,02 с, Т2 = 0,0167 с, Приведенные угловые скорости элементарных масс на графиках показаны в виде линейных скоростей.
Характеристический полином электромеханической системы получен в виде:
¥(5) = 1,103 • 10-25 + 9,26 • 10-23 • 5 + 9,303 • 10-9 • 52 + 4,83 • 10-6 • 53 + +6,57 • 10-5 • 54 + 5,79 • 10-4 • 55 + 3,68 • 10-3 • 56 + 0,0 1 86 • 57 + 0,0773 • 58 + +0,272 • 59 + 0,819 • 510 + 2,14 • 511 + 4,897 • 512 + 9,84 • 513 + +17,43 • 514 + 27,3 • 515 + 37,7 • 516 + 46,1 • 517 + 49,7 • 518 + 47,02 • 519 + +38,9 • 520 + 27,9 • 521 +17,2 • 522 + 8,99 • 523 + 3,92 • 524 +1,39 • 525 + +0,387 • 526 + 0,081 • 527 + 0,0118 • 528 + 0,00107 • 529 + 5,21 • 10-5 • 530 +
+1,56 • 10-6 • 531 + 3,11 • 10-8 • 532 + 4,09 • 10-10 +1,428 • 10-14 • 535 .
533 + 3,39 • 10-
Рис. 3. Структурная схема одного канала нечеткого регулятора
При расчете регулятора использованы только доминирующие корни характеристического уравнения, рассчитанные по разработанной нами методике.
Для демпфирования колебаний деформаций упругих элементов конвейера разработан двухканальный нечеткий регулятор (НР) на основе алгоритма «минимаксного» нечеткого вывода Мамдани (Max — Min Inferense). Для входных и выходных параметров НР использованы одни и те же функции принадлежности треугольного вида [7], [8]. Структурная схема НР в системе MATLAB представлена на рис. 3.
НР состоит из аналого-цифрового преобразователя (АЦП) — фиксатора нулевого порядка Zero-Order Hold 1 с шагом квантования h, блоков оценки первой e и второй e производной ошибки e системы, блоков нормировки входных (normin) и выходных (normout) сигналов, центрального блока НР Fuzzy Logic Controller и выходного цифроаналогового преобразователя — фиксатора Zero-Order Hold 2. Функции принадлежности (ФП) (membership functions) и база правил (rules) задаются в центральном блоке НР [9], [10].
Диапазоны изменения входных параметров НР e, e, e и выходных параметров и1(2) отображены на единое универсальное множество М = [0, 1], на котором задаются три нечетких подмножества, функции принадлежности которых имеют треугольную форму. Для каждой лингвистической величины получены аналитические выражения вида [11]:
^ (x) = 1 - x, x е [0, 1];
x) = x, x е [0, 1];
[2x, x е [0, 0,5] x [2(1 - x), x е [0,5, 1]
Расчет четких величин x*, x*, x* определяется по формулам:
Х* = (ё — ^min ) / (emax — 6min X
Х2 = (ё - ^min ) / (^max - ^min );
Х3 = (ё - emin ) / (emax - emin );
Лингвистическое правило управления НР имеет вид: Если (e* = а/) и (e* = «2) и (ё* = «3) ,
то (м*(2) = fl') , j = 1,3 , (8)
где а/, «2, «3 — лингвистические оценки входных параметров, а aj — лингвистические оценки управляющих воздействий на объект. Таким образом, все лингвистические переменные разделены на отрицательные (j = 1), положительные (j = 2) или близкие к нулю (j = 3) [12].
Для определения конкретного значения управляющего воздействия формируется «результирующая фигура», ограниченная результирующей ФП. Производится поиск абсциссы xc «центра тяжести» результирующей фигуры и полученное значение по формуле
U3(2) = Umin + (Umax - ^min ) ' ХС (9)
преобразуется в четкое значение управляющего воздействия на объект управления [13].
По результатам математического моделирования в системе MATLAB с использованием пакета нечеткой логики (Fuzzy Logic Toolbox) получены переходные функции частотного привода конвейера с НР — скорости ю. и ю2 (рис. 4 ).
м!С 10
«ц
/ / у < >
/ /
/У
О 1С 20 30 Ф 50 i. С Рис. 4. Переходные функции скоростей ю., ю2 конвейера с НР
Сравнение с соответствующими переходными функциями в электроприводе конвейера без регулятора показывает практически полное демпфирование колебаний скоростей элементарных масс конвейера и снижение деформаций упругого звена до 70% от значений деформаций в приводе конвейера без НР.
СПИСОК ЛИТЕРАТУРЫ
1. Ребенков Е. С. Оптимальное управление электроприводом ленточного конвейера / Математические методы в технике и технологиях: сборник трудов XXI Международной научной конференции. Т. 6. -Саратов: СГТУ, 2008. - С. 27-29.
2. Ребенков Е. С. Демпфирование колебаний в приводе ленточного конвейера // Горный информационно-аналитический бюллетень. -2009. - № 8. - С. 183-188.
3. Поляков А. Е., Чесноков А. В., Филимонов Е. М. Электрические машины, электропривод и системы интеллектуального управления электромеханическими комплексами. - М.: Форум: Инфра-М, 2015. -224 с.
4. Ребенков Е. С. Цифровой регулятор демпфирования упругих колебаний в многомассовом электроприводе // Горный информационно-аналитический бюллетень. - 2014. - № 2. - С. 330-335.
5. Чикуров Н. Г. Моделирование систем и процессов. - М.: РИОР: Инфра-М, 2013. - 398 с.
6. Шишмарев В. Ю. Теория автоматического управления. - М.: Издательский центр «Академия», 2012. - 352 с.
7. Конышева Л.К., Назаров Д.М. Основы теории нечетких множеств. - СПб.: Питер, 2011. - 192 с.
8. Mamdani E. H. Application of fuzzy logic to approximate reasoning using linguistic synthesis // IEEE Transactions on Computers. 1977. Vol. 26, no. 12. pp. 1182 - 1191.
9. Емельянов С. Г., Титов В. С., Бобырь М. В. Автоматизированные нечетко-логические системы управления. - М.: Научная мысль. Инфра-М, 2012. - 176 с.
10. Ross T. J. Fuzzy logic with engineering applications. McGraw-Hill: 1995. 600 p.
11. Zadeh L. A. Fuzzy logic // IEEE Transactions on Computers. 1988, Vol. 21, no. 4. pp. 83-93.
12. Takai T, Sugeno M. Fuzzy identification of systems and its applications to modeling and control // IEEE Transactions on Systems, Man and Cybernetics. 1985, Vol. 15, no. 1, pp. 116-132.
13. Пегат А. Нечеткое моделирование и управление; пер. с англ. -М.: БИНОМ. Лаборатория знаний, 2009. - 798 с. ЕШ
КОРОТКО ОБ АВТОРЕ
Ребенков Евгений Степанович - кандидат технических наук, доцент, e-mail: diolallin@mail.ru, Новомосковский институт (филиал) Российского химико-технологического университета имени Д.И. Менделеева. 301650, Новомосковск.
Gornyy informatsionno-analiticheskiy byulleten'. 2016. No. 4, pp. 68-75. E.S. Rebenkov
THE INDISTINCT REGULATOR OF THE FREQUENCY OPERATED ELECTRIC DRIVE OF THE CONVEYOR WITH ELASTIC TIES
Transfer function of electromechanical system of the frequency operated two-engine electric drive of the conveyor line is developed in a vector-matrix form. The indistinct regulator of speed of movement of the conveyor, damping fluctuations of speeds of masses and deformations in elastic communications is synthesized. Transients of coordinates of the drive with the indistinct regulator are given.
Key words: indistinct regulator, multimass model of the drive of the conveyor, elastic ties, damping of fluctuations.
AUTHOR
Rebenkov E.S., Candidate of Technical Sciences, Assistant Professor, e-mail: diolallin@mail.ru,
Novomoskovsk Institute (branch) of D. Mendeleev University of Chemical Technology of Russia, 301650, Novomoskovsk, Russia.
REFERENCES
1. Rebenkov E. S. Matematicheskie metody v tekhnike i tekhnologiyakh: sbornik trudov XXI Mezhdunarodnoy nauchnoy konferentsii. T. 6. (Mathematical methods in techniques and technologies: 21st International Conference Proceedings, vol. 6), Saratov, SGTU, 2008, pp. 27-29.
2. Rebenkov E. S. Gornyy informatsionno-analiticheskiy byulleten'. 2009, no 8, pp. 183-188.
3. Polyakov A. E., Chesnokov A. V., Filimonov E. M. Elektricheskie mashiny, elektrop-rivod i sistemy intellektual'nogo upravleniya elektromekhanicheskimi kompleksami (Electrical machines, electrical drives and intellectual control of electromechanical assemblies), Moscow, Forum: Infra-M, 2015, 224 p.
4. Rebenkov E. S. Gornyy informatsionno-analiticheskiy byulleten'. 2014, no 2, pp. 330-335.
5. Chikurov N. G. Modelirovaniesistem iprotsessov (Modeling systems and processes), Moscow, RIOR: Infra-M, 2013, 398 p.
6. Shishmarev V. Yu. Teoriya avtomaticheskogo upravleniya (Theory of automated control), Moscow, Izdatel'skiy tsentr «Akademiya», 2012, 352 p.
7. Konysheva L. K., Nazarov D. M. Osnovy teorii nechetkikh mnozhestv (Elements of theory of fuzzy sets), Saint-Petersburg, Piter, 2011, 192 p.
8. Mamdani E. H. Application of fuzzy logic to approximate reasoning using linguistic synthesis // IEEE Transactions on Computers. 1977. Vol. 26, no. 12. pp. 1182 1191.
9. Emel'yanov S. G., Titov V. S., Bobyr' M. V. Avtomatizirovannye nechetko-logicheskie sistemy upravleniya (Automated fuzzy-logic control systems), Moscow, Nauchnaya mysl'. Infra-M, 2012, 176 p.
10. Ross T. J. Fuzzy logic with engineering applications. McGraw-Hill: 1995. 600 p.
11. Zadeh L. A. Fuzzy logic. IEEE Transactions on Computers. 1988, Vol. 21, no. 4. pp. 83-93.
12. Takai T., Sugeno M. Fuzzy identification of systems and its applications to modeling and control. IEEE Transactions on Systems, Man and Cybernetics. 1985, Vol. 15, no. 1, pp. 116-132.
13. Pegat A. Nechetkoe modelirovanie i upravlenie; per. s angl. (Fuzzy modeling and control. English-Russian translation), Moscow, BINOM. Laboratoriya znaniy, 2009, 798 p.
UDC 621.867.2-83
