УДК 681.5.03
Р.С. ВОЛЯНСКИЙ, канд. техн. наук, доц., докторант, Днепродзержинский государственный технический университет, Днепродзержинск,
А.В. САДОВОЙ, д-р техн. наук, проф., проректор по научной работе, Днепродзержинский государственный технический университет, Днепродзержинск,
Н.В. СЛИПЧЕНКО, ст.преп., Днепродзержинский государственный технический университет, Днепродзержинск
ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ ЗАМКНУТОЙ ЭЛЕКТРОМЕХАНИЧЕСКОЙ СИСТЕМЫ С ПИ-РЕГУЛЯТОРОМ КОМПЛЕКСНОГО ПОРЯДКА
Составлено матричное уравнение движения обобщенной замкнутой динамической системы (ОЗДС) с интегродифференциалом комплексного порядка. Найдена частотная интерпретация интегродифференциала комплексного порядка. Составлена матричная частотная передаточная функция (МЧПФ) исследуемой ОЗДС. Рассмотрен пример определения МЧПФ для электропривода постоянного тока. Ил.: 1. Библиограф.: 8 назв.
Ключевые слова: динамическая система, интегродифференциал комплексного порядка, матричная частотная передаточная функция, электропривод.
Постановка проблемы. Последние достижения микропроцессорной и микроконтроллерной техники в разработке недорогих мощных микроконтроллеров создали предпосылки для усовершенствования существующих и разработки новых систем управления электромеханическими объектами и технологическими процессами. Оптимизация таких систем по различным критериям позволяет улучшить характеристики систем управления и, как следствие, повысить качество выпускаемой продукции при одновременном снижении ее себестоимости. Одним из направлений оптимизации является использование интегродифференциалов дробной размерности координат объекта управления [1, 2].
Анализ литературы. Системам управления, регуляторы которых работают на дробномерных линиях равновесия, посвящено большое количество работ как отечественных [3 - 5], так и зарубежных авторов [6]. Приведенные в них результаты показывают, что использование в составе алгоритма управления производных и интегралов дробной
© Р.С. Волянский, А.В. Садовой, Н.В. Слипченко, 2013
размерности позволяет улучшить качество процессов управления, создавая предпосылки для обобщения операций
интегродифференцирования на случай комплексной размерности [7]. Последнее обуславливает возникновение актуальной научной задачи, посвященной исследованию свойств замкнутых динамических систем с интегродифференциалами комплексного порядка. Одним из подходов к решению этой задачи является анализ динамических систем в частотной области, базирующийся на определении и исследовании частотной передаточной функции (ЧПФ) [8].
Цель статьи - разработка метода определения ЧПФ замкнутой динамической системы с интегродифференциалом комплексного порядка и исследование ее частотных характеристик.
1. Теоретические основы метода. Рассмотрим обобщенную электромеханическую систему, уравнения движения которой в относительных единицах можно представить следующим образом
рУ = ЛУ + ми, (1)
где р - оператор Лапласа; У - вектор переменных состояния; Л -матрица коэффициентов, определяющих свободное движение объекта управления; м - вектор коэффициентов, которые определяют управляемое движение объекта; и - управляющее воздействие.
Будем считать, что движение динамического объекта (1) происходит под действием управляющего воздействия
и = (кр + р)х(у* - у), (2)
где Кр, К- матрицы-строки коэффициентов усиления
пропорциональной и интегральной частей регулятора; У , У - матрицы-столбцы желаемых и текущих значений переменных состояний.
Регулятор, реализующий алгоритм (2), является обобщением ПИа-регулятора с дробномерным интегратором [6]. Такой регулятор будем называть ПИс-регулятором комплексного порядка.
Подставив управляющее воздействие (2) в уравнение движения объекта (1), представим его следующим образом
[рЕ - Л + М(Кр + р-“^Кг )]У = М(Кр + р-“^Кг )У*. (3)
Заменив в уравнении (3) оператор Лапласа р на оператор Фурье ую, получим следующее выражение для амплитудно-фазовой характеристики (АФХ) динамической системы (3)
[/юЕ - А + М(Кр + /ю-“-/РКг )]¥ = М(Кр + (/ю)-“-/Р К )У*. (4)
Из уравнения (4) может быть получена матричная частотная передаточная функция системы (3)
У = [/юЕ - А + М(Кр + /ю-“-/ЄКг)]-1 М(Кр + (/ю)-“-/Р Кг) . (5)
Выражения (4) и (5) содержат комплексный оператор (/га)"
“-jP
главное значение которого может быть определено через элементарные функции следующим образом
= fRe (â) + jfIm (â) =
_ ep"/2cos(|3lnю + œrc /2) . ep"/2 sin((3lnю + œrc /2) (6)
" â“ + j â“ '
Подстановка выражения (6) в зависимость (5) позволяет записать МЧПФ следующим образом
:Yjâ\ = [jâE - A + M(Kp + fRe (â)Ki + fm (Ю)К, J-1 x
Y (jâ) (7)
x м(Кр + /Re (ffl)K, + j/m (ffl)K, )
или
Y(jra) adj[-A + M(Kp + /Re(ю)К,)+ j(raE + fIm(ffl)MK,)J ^
Y* (jro) det (- A + M(Kp + /Re (Ю)К, )+ j(raE + /m (ю)МК, )J) X (8)
Lp + J R^^/KU + ЛШЕ + J Im\
X lKp + (/Re (ю) + j/Im (ю))Кі J>
где adj[.] - функция, возвращающая матрицу алгебраических
дополнений аргумента.
Выражение (8) является МЧПФ системы с ПИс-регулятором комплексного порядка и позволяет определять частотные характеристики широкого класса замкнутых систем не только с П, ПИ, ПИа, ПИс -регуляторами, но и, принимая во внимание, что тип операции интегродифференцирования определяется знаком степени оператора дифференцирования, систем с ПД-, ПДа-, ПДс-регуляторами. В случае
исследования систем с ПИс1Дс2-регуляторами или более сложными, выражение (8) дополняется соответствующими составляющими, не ограничивая при этом область применения разработанного метода.
2. Пример использования метода. В качестве примера определим частотную передаточную функцию для контура регулирования скорости электропривода постоянного тока с ПИс-регулятором, движение которого описывается следующими уравнениями динамики:
руі = аиУ 2;
РУ2 = а2іУі + а22У2 + т2и;
(9)
и =
к
к + к2
Р “+;РУ
(у* -У1 )•
Составим для системы (9) матрицы, входящие в выражение (8)
0 (0 1
А = ; М =
V а21 а22 У Vт2 у
Кр =(кі 0); (10)
К =(к2 0) •
Подставив матрицы (10) в выражение (8), получим МЧПФ системы
(9)
У _ adj[W ]
( 0 1
У* det (^)
[(к1 0)+ /не (ю) + ІЇіт (ю)(к2 0)]
(11)
Vт2 У
где
( 0 а121 (0 1
W = - +
V а21 а22 У Vт2 У
[(к1 0)+ /Не (юХк2 0)]
+ І
(1 01 (01
ю + /1т (ю) (к2 0)
,0 1У чт2 У
+
После выполнения алгебраических преобразований МЧПФ (11) примет вид
( (
У
У*
12т2
к^“ + к2Є2
рл
сое) ріп(ю)+ “Л 1 + і'віпі Р1п(ю)+ “Л
-Ю2'“ - Іа22ю1+“ +ю“(к1а12т2 - а12а21) + В
Jm2Ю
к1юа + к2е2
сов| р1п(ю)+“л 1+ріп(ю)+“л
— ю 2+а - Іа22ю1+“ +ю“(к1а12т2 - а12а21) + В
(13)
рл
где Б = к2а12т2е 2
I ™ ✓ ч “Л 1 • • ( , ч “л
С0б| ріп(ю) + — I + ^т! ріп(ю) + —
Размерность МЧПФ (13) задается числом независимых координат состояния, при этом количество ненулевых столбцов определяется числом задающих воздействий, подаваемых на систему. Выделив из МЧПФ (13) первый элемент первого столбца, получим искомую частотную передаточную функцию
( рл
^га“ + к2е 2
Ф(/ю) =
I , ч “л 1 . ( _ ч “л
COS| ріп(ю) + — I + JSinl ріп(ю) + —
-га Іа22ю +Ю (к1а12т2 а12а21)+ В
(14)
Приняв “ = 1, р = 0 и воспользовавшись зависимостью (14), определим ЧПФ замкнутой системы с классическим ПИ-регулятором
Ф(ію) =
а12т2 (к1Ю + Ік2 )
-га3 -Іа22ю2 +ю(к1а12т2 - а12а21)+ Ік2а12т2
(15)
Умножение числителя и знаменателя ЧПФ (15) на } приводит к выражению, которое может быть получено классическим методом. Считая параметры исследуемой системы единичными, построим на основании ЧПФ (14) и (15) соответствующие АФХ (рис.).
0
0
12"*2
Рис. АФХ исследуемой системы при различных параметрах интегродифференциала комплексного порядка.
Анализ полученных АФХ при а = в = 0 и а = 1, в = 0, описываемых выражением (14) и найденные классическим способом совпадают. Приведенные на рис. АФХ позволяют сделать вывод, что при частоте ю равной нулю единичным коэффициентом усиления обладает только система с классическим ПИ-регулятором. Указанный недостаток систем с интегродифференциалом комплексного порядка может быть устранен путем увеличения порядка интегродифференциала.
Выводы. Использование при анализе частотных свойств систем матричного исчисления совместно с теорией функции комплексного переменного позволяет получить простые выражения для МЧПФ динамических систем с интегродифференциалом комплексного порядка.
Анализ зависимостей (13) и (14) показывает, что ЧПФ замкнутой системы с интегродифференциалом комплексного порядка независимой переменной, по которой находится передаточная функция, содержит гармонические составляющие, определяющиеся порядком интегродифференциала и частотой входного воздействия. Принимая во внимание, что интегродифференциал комплексного порядка является обобщением операций дифференцирования и интегрирования, можно сделать вывод, что ЧПФ с гармоническими составляющими являются обобщением известных ЧПФ динамических систем.
Анализ полученных АФХ показывает, что наличие в ЧПФ гармонических составляющих в определенном частотном диапазоне
позволяет повысить коэффициент усиления системы, обеспечивая тем самым форсировку замкнутой системы в динамических режимах.
Список литературы: 1. Васильев В.В. Дробное исчисление и аппроксимационные методы в моделировании динамических систем / В.В. Васильев, Л.А. Симак. - К.: НАН Украины, 2008. - 256 с. 2. Нахушев А.М. Дробное исчисление и его применение / А.М. Нахушев. - М: ФИЗМАТЛИТ, 2003. - 272 с. 3. Бушер В.В. Динамические свойства систем управления с дробным порядком астатизма / В.В. Бушер // Электротехнические и компьютерные системы.
- 2010. - № 1. - С.13-16. 4. Васильев А.В. Математические модели ПИД-контроллеров динамических систем целого и дробного порядков на основе ^-преобразования /А.В. Васильев // Информационные и телекоммуникационные технологии. - 2013. - Т. 322.
- № 2. - С. 21-26. 5. Учайкин В.В. Метод дробных производных / В.В. Учайкин. - Ульяновск: "Артишок", 2008. - 512 с. 6. Podlubny I. Fractional differential equation / I. Podlubny. - San Diego: Academic Press, 1999. - 340 p. 7. Самко С.Г. Интегралы и дифференциалы дробного порядка и некоторые их приложения / С.Г. Самко, А.А. Килбас, О.И. Маричев. - Минск: Наука и техника, 1987. - 688 с. 8. Авсиевич А.В. Частотные критерии устойчивости Михайлова и Найквиста для моделей систем автоматического управления с дробным порядком / А.В. Авсиевич, В.В. Авсиевич // Вестник транспорта Поволжья. - 2010. - № 1. -С. 35-41.
Поступила в редакцию 07.06.2013
УДК 681.5.03
Частотні характеристики замкнутої електромеханічної системи з ПІ-регулятором комплексного порядку / Волянський Р.С., Садовой О.В., Сліпченко Н.В. // Вісник НТУ ''ХПГ'. Серія: Інформатика та моделювання. - Харків: НТУ 'ХПГ". - 2013. - № 19 (992). - С. 3 - 9.
Складено матричне рівняння руху узагальненої замкнутої динамічної системи (ОЗДС) з інтегродифференціалом комплексного порядку. Знайдена частотна інтерпретація інтегродифференціала комплексного порядку. Складена матрична частотна передатна функція (МЧПФ) досліджуваної ОЗДС. Розглянуто приклад визначення МЧПФ для електроприводу постійного струму. Гл.: 1. Бібліограф.: 8 назв.
Ключові слова: динамічна система, інтегродифференціал комплексного порядку, матрична частотна передатна функція, електропривод.
UDC УДК 681.5.03
Frequency response of closed-loop electromechanical system with PI controller of complex order / Volyanskiy R.S., Sadovoy A.V., Slipchenko N.V. // Herald of National Technical University "KhPI". Subject issue: Information Science and Modeling. - Kharkov: NTU "KPI". - 2013. - № 19 (992). - Р. 3 - 9.
Matrix equation of generalized closed-loop complex order dynamic system (GCDS) motion is given. Frequency interpretation of complex order integrator is found. Frequency transfer function matrix (MCHPF) investigated GCDS is shown. MCHPF for DC electric drive is determined. Figs. 1. Ref: 8 titles.
Keywords: dynamic system, integrator of complex order, matrix frequency transfer function, electric drive.