УДК 62-523.8:519.71 DOI:10.30724/1998-9903-2021-23-6-119-133
АЛГОРИТМ ПАРАМЕТРИЧЕСКОЙ ИДЕНТИФИКАЦИИ ЭЛЕКТРОПРИВОДА ПОСТОЯННОГО ТОКА С ПРИМЕНЕНИЕМ ИНВЕРСНОЙ МОДЕЛИ
Малёв1 Н.А., Погодицкий1 О.В., Чиляева1 М.Р., Имамиев2 А.Р.
1Казанский государственный энергетический университет, г. Казань, Россия 2ФЛ АО «Сетевая компания» - «Дирекция строящихся объектов», г. Казань, Россия
maleeev@mail. ru
Резюме: ЦЕЛЬ. В современных системах управления различных промышленных агрегатов основой для автоматизации технологических процессов является электрический привод. Высокие требования к качеству процесса управления, надёжности и энергетической эффективности технологических процессов и промышленных агрегатов определяют постановку и решение научно-практических задач по разработке и созданию новых перспективных систем управления электроприводов (СУЭП), позволяющих сохранить требуемое качество функционирования в условиях влияния дестабилизирующих факторов. Данные факторы могут воздействовать на систему управления в виде внешних координатных возмущений и суммироваться с соответствующими сигналами замкнутой СУЭП, а также выражаться в форме параметрических возмущений, приводящих к нестабильности параметров СУЭП в процессе эксплуатации. Возникает необходимость создания такой системы управления, которая позволит осуществить идентификацию параметров СУЭП за счет введения в структуру системы соответствующих алгоритмов идентификации. МЕТОДЫ. При решении поставленной задачи применялся беспоисковый градиентный метод адаптивной идентификации, реализованный средствами программной среды MatLab. РЕЗУЛЬТАТЫ. В работе рассматривается задача параметрической идентификации электропривода с двигателем постоянного тока на основе определения функций чувствительности. Показатель качества процесса идентификации выбран из условия обеспечения близости дополнительного движения, формируемого моделью, и фактического дополнительного движения при вариациях нестабильных параметров. Для построения алгоритма параметрической идентификации используется инверсная модель исследуемой СУЭП, показателем качества является квадрат невязки, а идентифицируемым параметром - общий коэффициент передачи электропривода. ЗАКЛЮЧЕНИЕ. Моделирование в программной среде MatLab показало высокую робастность разработанного алгоритма идентификации к параметрическим возмущениям, которые не оказывают влияния на установившееся значение идентифицируемого параметра. Предложенный в работе способ компенсации статического момента сопротивления электропривода позволил обеспечить также малую чувствительность алгоритма и к внешним возмущениям. Идентификация общего коэффициента передачи электропривода постоянного тока осуществляется с погрешностью, не превышающей 0,5% в режиме реального времени в условиях влияния возмущений различной физической природы.
Ключевые слова: алгоритм; нестабильные параметры; электропривод; инверсная модель; модель чувствительности; дополнительное движение; невязка.
Для цитирования: Малёв Н.А., Погодицкий О.В, Чиляева М.Р., Имамиев А.Р. Алгоритм параметрической идентификации электропривода постоянного тока с применением инверсной модели // Известия высших учебных заведений. ПРОБЛЕМЫ ЭНЕРГЕТИКИ. 2021. Т. 23. № 6. С. 119-133. doi:10.30724/1998-9903-2021-23-6-119-133.
ALGORITHM FOR PARAMETRIC IDENTIFICATION OF THE DC ELECTRIC DRIVE
USING THE INVERSE MODEL
NA. Malev1, OV. Pogoditsky1, MR. Chilyaeva1, AR. Imamiev2
1Kazan State Power Engineering University, Kazan, Russia 2Branch of JSC «Grid Company» - «Directorate of facilities under construction»,
Kazan, Russia
maleeev@mail.ru
Abstract: THE PURPOSE. In modern control systems of various industrial units, the basis for the automation of technological processes is an electric drive. High requirements for the quality of the control process determine the formulation and solution of scientific and practical tasks for the development and creation of new promising control systems for electric drives (EDCS), allowing to maintain the required quality of functioning under the influence of destabilizing factors. These factors can affect the control system in the form of external disturbances and be summed up with the corresponding signals of the closed EDCS, and also be expressed in the form of parametric disturbances. The need arises to create such a control system that will allow for the identification of the parameters of the EDCS by introducing the appropriate identification algorithms into its structure. METHODS. When solving the problem, a search-free gradient method of adaptive identification was used, implemented by means of the MatLab software environment. RESULTS. The paper considers the problem ofparametric identification of an electric drive with a DC motor based on the definition of sensitivity functions. Wherein to construct an algorithm for parametric identification, an inverse model of the studied EDCS is used, the quality indicator is the squared discrepancy, and the identifiable parameter is the overall transmission coefficient of the electric drive. CONCLUSION. Modeling in the MatLab software environment showed a high robustness of the developed identification algorithm to parametric disturbances that do not affect the steady-state value of the identified parameter. The proposed method for compensating the moment of resistance electric drive also made it possible to provide a low sensitivity algorithm to external disturbances. The identification of the overall transmission coefficient of a DC electric drive is carried out with an error not exceeding 0.5% in real time under the conditions of the influence of disturbances of various physical nature.
Keywords: algorithm; unstable parameters; electric drive; inverse model; sensitivity model; additional motion; discrepancy.
For citation: Malev NA, Pogoditsky OV, Chilyaeva MR, Imamiev AR. Algorithm for parametric identification of the dc electric drive using the inverse model. Power engineering: research, equipment, technology. 2021;23(6):119-133. doi:10.30724/1998-9903-2021-23-6-119-133.
Введение
Процесс эксплуатации систем управления электроприводов сопровождается, как правило, влиянием дестабилизирующих факторов, под действием которых изменяются параметры системы и, в свою очередь, показатели качества процесса управления [1 - 4]. В этой связи в современных перспективных СУЭП актуальным является применение автоматических систем идентификации, позволяющих вычислять оценки нестабильных параметров на этапе функционирования в режиме реального времени [5 - 7]. Способность системы изменять свои свойства при изменении совокупности параметров определяют её параметрическую чувствительность [8, 9].
Условие нормального функционирования СУЭП можно сформулировать следующим образом: множество значений фазовых координат х, зависящих от вектора параметров СУЭП х и вектора внешних воздействий f которые реально могут иметь место в процессе функционирования СУЭП должно быть строго включено в множество значений х, при которых состояние СУЭП удовлетворяет предъявляемым требованиям. Для выполнения данного условия необходимо предъявить к системе требование нулевой параметрической
чувствительности или параметрической инвариантности \xy,f\= 0, что практически
неосуществимо [10, 11]. Поэтому требование параметрической инвариантности выражается
в форме е-чувствительности |x"x/| < s , которое означает, что принятый показатель качества
функционирования СУЭП J находится в заданных пределах Jm;n < J < Jmax. Одним из методов обеспечения ограничений на показатель качества является выбор номинальных значений компонент вектора параметров СУЭП х с учётом технологических отклонений в процессе серийного производства [12]. В настоящей работе определение экстремального значения показателя качества J осуществляется аналитическим путем в сочетании с градиентным способом идентификации.
Методы
Показатели качества функционирования СУЭП могут быть различными. В инженерной практике широко применяются показатели, основанные на невязке [13]
s(t) = Уо )"Ут (t) , С1)
формируемой как разность между действительным у0 (7) и эталонным ут (7)
значениями выходной координаты СУЭП. Относительная простота реализации алгоритма идентификации параметров СУЭП обеспечивается при показателе качества функционирования[14, 15]
3 = е2 (7) (2)
или в векторно-матричной форме
3 = (Уо-Ут) (Уо-Ут) - (3)
Рассмотрим математическую постановку задачи идентификации вектора параметров х на основе выражения (3). Пусть собственное движение электропривода описывается выражением
х = Р{х,и). (4)
Тогда дополнительное движение от вариаций вектора параметров СУЭП запишется с помощью соотношения
Ах = Г (Дх). (5)
Для получения оценки вектора х решается обратная задача, которой соответствует формальной отношение
Дх = Г* (Дх). (6)
Показатель качества (3) формируется следующим образом. Текущее значение вектора состояния СУЭП
хо (X) = хт+ххДХ , (7)
где хт - требуемое значение вектора состояния модели СУЭП.
С учётом (7) показатель качества (3) для вектора промежуточных фазовых координат электропривода принимает вид
3 = (Дх - ххДх)Т (Дх - ххДх), (8)
где Дх= х0-хт.
Показатель качества (8) характеризует степень близости дополнительного движения, формируемого моделью, и фактического дополнительного движения [16, 17].
Определим оценку вектора х, исходя из условия минимума критерия J. Преобразуем выражение показателя качества (8)
з = ДхГ Дх - ДхГхГ Дх - Дхг ххДх+ДхГхГххДх и продифференцируем полученное выражение по вектору А/:
Т '■у '■у '■у
^(Д -) = -хх Дх - Дх хх+ 2Дх хх хх . (9)
Произведем в (9) матричные преобразования и приравняем данное выражение к
нулю:
2ДхТхх х{ = х^Дх + ДхТх{. (10)
Введём обозначения
т т т
Р1 = 2хТхх и Р2 = хТ Дх + ДхТ хх. (11)
С учётом введённых обозначений, из выражения (10) следует, что
Р ДхТ = Р2 , (12)
откуда вариация вектора параметров х определяется по формуле
ДхТ = Р1-1 Р2. (13)
При выполнении условия (13) оценка вектора параметров СУЭП находится из выражения
х = хт + Дх . (14)
Матрицы Р1 и Р2 зависят от матрицы чувствительности хх. Поэтому для
реализации выражения (13) в алгоритме идентификации параметров СУЭП необходимо предусмотреть вычислитель функций чувствительности [18, 19].
Объектом настоящего исследования является электропривод постоянного тока, получающий питание от управляемого выпрямителя и охваченный отрицательной обратной
121
связью по скорости. В качестве идентифицируемого параметра принят общий коэффициент передачи электропривода К.
Закон изменения идентифицируемого параметра
Щ = -№/, (15)
д
где = — - символ градиента; Xк - коэффициент, характеризующий скорость
дК
идентификации параметра.
С учетом принятого показателя качества функционирования СУЭП (2) вычислим градиент
к дК дК w дК
Передаточную функцию электропривода в разомкнутой структуре можно представить дробно-рациональным выражением в операторной форме
КX aks к im + am_lsml +... + axs + a0 )
W ( s ) = —^=°----,
Y^b^k bnSn + bn_\sn +... + bs + bo
k=0
где K - общий коэффициент передачи электропривода (добротность по скорости). Передаточная функция замкнутой системы относительно сигнала ошибки
ÁU(s) 1 B(s)
(17)
®Á (s) =
Uin (s) 1 + W (s) B(s) + KA(s)
(18)
Принимая за основу выражение (18), запишем передаточную функцию инверсной
модели
^m{,_Uim (s) _ B(s) + KA(s) _ | | KA(s) Á ( ) ÁU(s) B(s) B(s) .
Дифференцируя (19) по общему коэффициенту передачи системы K, получим
(19)
соответствующую функцию чувствительности передаточной функции инверсной модели
ЗФд (s) A(s)
(20)
дК В(5)
которая не зависит от настраиваемого параметра, то есть функция чувствительности инвариантна к изменению К.
По выражениям (15) и (16) определим алгоритм беспоисковой градиентной идентификации
ав(0
ед = -№/ = -а*2в(0
ЗК (t)
Переходя к преобразованию Лапласа, получим
К (s) = _ ^ [Uin (s) _ Uim (s)]
3Uim (s) дК
!%k
[Uin (s) _ Uim (s)]
ЗФ? (S)
дК
ÁU (s)
и тогда, с учётом формулы (20), добротность по скорости найдём из выражения
К(s) = _^[Uin (s) _ Uim (s)]
A( s) B(s)
ÁU (s)
(21)
.(22)
(23)
где (5) - входное задающее воздействие; ищ (5) - выходной сигнал инверсной модели; Ди(5) - сигнал рассогласования.
Следует отметить, что в рассматриваемом случае роль эталонной модели выполняет инверсная модель, коэффициенты полиномов которой зависят от номинальных значений параметров электропривода за исключением добротности по скорости К, а невязка, соответственно, вычисляется по формуле
е(0= ит к)-ищ (О. (24)
На основании полученных выражений структурная схема динамической модели инверсного идентификатора общего коэффициента передачи электропривода примет вид, показанный на рис. 1.
m
Рис. 1. Структурная схема динамической модели Fig. 1. Block diagram of the dynamic model of the
инвеРсного идентификатора общего inverse identifier of the total transmission coefficient
коэффициента передачи электропривода of the electric drive
Результаты и обсуждение
В качестве электромеханического преобразователя в исследуемом электроприводе применен серводвигатель постоянного тока типа PIVT-6-25/3A с возбуждением от постоянных магнитов и встроенным тахогенератором с крутизной выходного напряжения 0,1 Вмин/об. Технические характеристики машины сведены в таблицу 1 [20].
Технические характеристики двигателя PIVT-6-25/3 А.
Таблица 1
Номинальное напряжение Ua, В Коэффициент передачи по моменту Km, Нм/А Электромагни тная постоянная времени Та, с Электромехан ическая постоянная времени Тт, с Частота вращения при максимально м моменте n, об/мин Максимальны й момент Mmax, Н-м
30 0,072 0,002 0,025 3000 0,225
Передаточная функция двигателя постоянного тока
jjr , ч 1/С
Wm (s) =-2-,
TaTms + Tms +1
(25)
где c = Km = 0,072 В-с/рад - коэффициент противо-ЭДС.
Двигатель получает питание от однофазного двухполупериодного мостового стабилизированного выпрямителя типа ВСА-5К с передаточной функцией
K
Wr = , (26)
Гrs +1
где —г = 15 - коэффициент усиления выпрямителя; Гг = 0,005 с - постоянная времени выпрямителя [20].
На выходе тахогенератора установлен ЛС-фильтр для сглаживания пульсаций выходного напряжения с постоянной времени Г/ = 0,001 с. Передаточная функция
тахогенератора с фильтром
Щ„ (*) =-—^, (27)
/ + 1
где —ю = 0,0104 В-с/рад - коэффициент передачи тахогенератора, численно равный крутизне выходного напряжения [20].
Передаточная функция замкнутого электропривода относительно сигнала рассогласования
ди^) 1
Фд (s) =
Uin (s) 1 + Wr (s)Wem (s)Wm (s)
(t>+i)(ra
Tffîf S + Tffîf s +1
ITs +1
)(/ +1)
(28)
|T„T„s2 + Tms + 1)(T/S +1) + KrKemKw
(Гrs + 1)(ГаГт" 1 гт" где —ет = 1/ с = 13,89 рад/В-с - коэффициент передачи двигателя.
Запишем добротность по скорости в виде произведения коэффициентов передачи всех элементов электропривода
К — КгКетК&
(29)
и определим передаточную функцию инверсной модели:
Фгт ( s)—
Uim (s )
(T> +1)(2
„s + 1)[TaTms + Tms + 1UT/-s +1) + К
1)(T/s +1)"
AU ( s) (Trs +1)( — 1 +
s + 1)(TaTmS2 + TmS + Hl^s +1)
К
1)(T/s +1)
(30)
(Trs + 1)(TaTms2 + Tms + 1|Ts +1)
1)(T/s +1)'
Продифференцируем выражение (30) по общему коэффициенту передачи электропривода
дФ'т ( s)
1
5К (Trs + 1)(TaTms2 + Tms + 1)(T/s +1)
(31)
и убедимся, что результат дифференцирования не зависит от идентифицируемого параметра, т.е. функция чувствительности инвариантна к изменению К (см. выражение (20)).
Структурная схема динамической модели идентификации общего коэффициента передачи исследуемого электропривода постоянного тока показана на рис. 2.
Um (s) AU (s) Кг 1/с
V Trs +1 T a^ms + Tms + 1
К
T/s + 1
Q(s)
(Trs + 1XTgTms2 + Tms + 1)(T/s +1)
(Trs + 1)(TaTms2 + Tms + 1)(T/s + 1)
->
Uim (s)
К?)
2Хь
s
Инверсный идентификатор
Рис. 2. Структурная схема динамической модели инверсного идентификатора общего
коэффициента передачи электропривода постоянного тока
______________J
Fig. 2. Block diagram of the dynamic model of the inverse identifier of the total transmission coefficient of a DC electric drive
На основании полученных выше выражений и схемы, показанной на рисункке 2, реализуем &ши11пк-модель инверсного идентификатора общего коэффициента передачи электропривода постоянного тока с применением программной среды ЫаЛаЬ. Соответствующая схема представлена на рисунке 3.
Рис. 3. Simulink-модель инверсного Fig. 3. Simulink-model of the inverse identifier of the идентификатора общего коэффициента передачи total transmission coefficient of a DC electric drive электропривода постоянного тока
На рисунке 3 блоки, выделенные жёлтым цветом, соответствуют исследуемой системе, блоки, выделенные голубым цветом, реализуют алгоритм идентификации общего коэффициента передачи электропривода согласно выражению (23), а блоки Scope K и Display K регистрируют и отображают идентифицируемый параметр.
В командное окно MatLab вводим значения параметров электропривода, соответствующие выражениям (25) - (27), характеризующие номинальный режим работы системы. Коэффициент, характеризующий скорость идентификации, принимаем X к = 500.
Результаты моделирования при номинальных значениях параметров представлены на рисунке 4.
K (t) K (t) ,o.e.
1
»
t, c
О 0.02 0.04 006 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Рис. 4. Результаты моделирования процесса идентификации общего коэффициента передачи электропривода К: 1 - номинальное значение К(Г); 2 - оценка К (/)
Fig. 4. Results of modeling the process of identification of the total transmission coefficient of the electric drive K: 1 - nominal value K(t); 2 - evaluation
В соответствии с выражением (29) номинальное значение общего коэффициента передачи электропривода Кпот = КгКетКю = 15-13,89 • 0,0104 « 2,167, при этом полученная оценка практически не отличается от номинального значения.
Дальнейшие исследования проведем при детерминированном изменении коэффициента К, а также при вариациях постоянных времени элементов электропривода. Результаты моделирования показаны на рисунках 5 и 6.
Проблемы энергетики, 2021, том 23, № 6 K (t), K (t) ,o.e.
4.5 Л 35 3 2.5 2 1.5 1
0.5 0 -0.5
О 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.16 0.2
Рис. 5. Результаты моделирования процесса Fig. 5. Results of modeling the process of идентификации общего кюффщиетта жредачи identification of the total transmission coefficient of электропривода K\ 1 - значение K = 2Knom; 2 - the electric drive K: 1 - value; 2 - estimate ; 3 -оценка K (t) = 2Knom; 3 - значение к = 0,5Knom; value; 4 - estimate 4 - оценка
Из анализа графиков, приведенных на рисунке 5 следует, что вычисление оценок общего коэффициента передачи электропривода при заданных изменениях K осуществляется с минимальным значением погрешности, не превышающем 0,001%, т.е. K « K. При этом время, затрачиваемое на идентификацию, не превышает 0,02 с.
Произведем подобное исследование при изменении электромеханической постоянной времени машины постоянного тока Tm и постоянной времени фильтра Tf и
K = Knom (см. рис. 6).
1
2
3
4
t, c
K (t), K (t), 2 o.e.
1
3
t, c
О 0.02 0.04 0.06 0.06 0.1 0.12 0.14 0.16 0.16 0.2
Рис. 6. Результаты моделирования процесса Fig. 6. Results of modeling the process of
идентификации общего коэффициента передачи identification of the total transmission coefficient of
электропривода K: 1 - номинальное значение the electric drive K: 1 - the nominal value of K(t); 2 -
K(t); 2 - оценка K (t) при уменьшении Tm и Tf an estimate with a decrease and by 20%; 3 - an
на 20%; 3 - оценка K(t) при увеличении Tm и estimate with an increase and by 20% Tf на 20%
Вариации указанных постоянных времени в заданных пределах приводят к изменениям переходных процессов вычисления оценок общего коэффициента передачи электропривода, которые характеризуются снижением быстродействия и временем установления около 0,1 с как для случая с наличием перерегулирования (график 2), так и при монотонном процессе (график 3). При этом установившиеся значения вычисляемых оценок для обоих случаев соответствует номинальному значению коэффициента передачи.
Проведенные исследования показали, что рассматриваемый в работе алгоритм идентификации является инвариантным к параметрическим возмущениям, обеспечивая
высокую точность процесса вычисления оценок общего коэффициента передачи исследуемого электропривода при сравнительно малых временных затратах.
Следует отметить, что функционирование электропривода в реальных условиях сопровождается внешними координатными возмущениями, к которым можно отнести момент сопротивления нагрузки Мс, а также шумы и наводки в цепях управления f. Структурная схема динамической модели исследуемого электропривода с учётом внешних возмущений представлена на рисунок 7 [21, 22].
F (s)
Um (s) S±\ AU(s)
Uout( s)
Kr
E ( s)
Fs +1
U„
Ua (s) 1/R la (s)
Tas +1
Me (s).
К„
Ie (s)
T/s +1
Q(s)
Рис. 7. Структурная схема динамической модели Fig. 7. Block diagram of a dynamic model of a DC электропривода постоянного тока с учетом electric drive taking into account disturbing возмущающих воздействий. influences.
На основании показанной на рисунке 7 схемы сформируем Simulink-модель электропривода. Неизвестное сопротивление якорной цепи R было измерено с помощью измерительного моста постоянного тока типа ММВ с погрешностью, не превышающей +5% [20]. Результаты измерения составили R = 5,15 Ом. Полученная Simulink-модель представлена на рисунке 8. Устанавливаем значение момента нагрузки Мс = 0,8Mmx = 0,8 • 0,225 = 0,18 Н- м, а сигнал шума зададим с помощью блока Signal Generator, выбирая во вкладке Wave form наименование сигнала Random и устанавливая амплитуду сигнала 0,3 В, частоту - 1 кГц.
Рис. 8. Simulink-модель инверсного Fig. 8. Simulink-model of the inverse identifier of the идентификатора общего коэффициента передачи total transmission coefficient of a DC electric drive электропривода постоянного тока с учётом taking into account external disturbances внешних возмущений
Моделирование производим при набросе нагрузки через 0,1 с после пуска. Результаты вычисления оценки K представлены на рисунке 9.
с
R
eTms
K (t), K (t), o.e.
1
2
1
t, c
О 0.05 0.1 0.1 Б 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Рис. 9. Результаты моделирования процесса Fig. 9. Results of modeling the process of
идентификации общего коэффициента передачи identification of the total transmission coefficient of
электропривода K: 1 - номинальное значение the electric drive K: 1 - nominal value K(t); 2 -
K(t); 2 - оценка K (t) при влиянии внешних evaluation under the influence of extemal
disturbances andf
возмущений Мс иf
Показанные на рисунке 9 зависимости иллюстрируют существенное влияние внешних возмущений на результат идентификации. Полученное в процессе вычисления значение К (/) сопровождается флуктуациями с амплитудой и частотой наведенной помехи
/, что затрудняет определение установившегося значения К.
Влияние момента сопротивления Мс проявляется в значительном отклонении
установившегося значения оценки Куст ) от заданной величины К, вызванном перепадом скорости вращения электропривода под действием нагрузки. Данное влияние приводит к изменению сигнала рассогласования Ди(5) и, в свою очередь, невязки Е(5) (24), в
результате чего возникает погрешность результата идентификации общего коэффициента передачи К согласно формулы (23). Из рисунка 9 следует, что относительная погрешность вычисления К = Кпот = 2,167 при полученном значении оценки Куст (/) = 1,352
составляет около 37%, что говорит о неудовлетворительном решении. Максимальная погрешность идентификации обусловлена именно моментной составляющей ошибки, поскольку возмущающее воздействие / несущественно влияет на Куст ), а к ошибке по
управлению рассматриваемая система инвариантна.
В целях снижения чувствительности алгоритма параметрической идентификации с применением инверсной модели к статическому моменту сопротивления нагрузки необходимо разработать способ определения Мс с последующим формированием корректирующего воздействия, компенсирующего влияние Мс и сводящего ошибку идентификации 1С (/) к минимуму. Запишем уравнение механического равновесия электропривода, выразив моменты через токи (см. схему на рис. 7):
'а(0-'с(0='(0, (32)
где 1а (/) - ток якоря; 1с (/) - ток, соответствующий моменту сопротивления;
, ч сТт dю ^)
1 и) = ——--— - ток, соответствующий динамическому моменту.
К dt
Из выражения (32) найдём
(33)
т.е. на основании информации о токе 1а ^) с использованием датчика тока и об ускорении
d ю и )
-— с использованием датчика скорости получим искомый сигнал, пропорциональный
dt
моменту сопротивления Мс .
Корректирующее воздействие, компенсирующее влияние Мс , формируется в виде сигнала напряжения
(34)
,(t ) — Ve (t )
где кс = 0,744 Ом - коэффициент пропорциональности, подобранный путем компьютерного эксперимента. Невязка (24) становится функцией трех переменных ) = / (ит, иш, ис ). 81шиИпк-модель инверсного идентификатора общего коэффициента передачи электропривода постоянного тока с компенсацией Мс представлена (рис. 10).
Рис. 10. Simulink-модель инверсного Fig. 10. Simulink-model of the inverse identifier of
идентификатора общего коэффициента передачи the total transmission coefficient of a DC electric
электропривода постоянного тока с drive with compensation компенсацией Мс
На рисунке 10 блоки, выделенные сиреневым цветом, иллюстрируют процедуру формирования компенсирующего воздействия ис (^).
Компенсацию возмущающего воздействия / осуществим с помощью фильтра нижних частот первого порядка с единичным коэффициентом усиления и постоянной времени, равной 0,012 с. Фильтрация осуществляется по отношению к сигналу рассогласования Ди(5). Результаты моделирования процесса идентификации общего коэффициента передачи электропривода показаны на рисунке 11.
4 3.5
3 2.5
2 15
О 0.05 0.1 0.15 0.2 0.25 03 0.35 04 0.45 05
К (t ), К (t ), o.e.
l\
\
- V -4- i
1 t, c
Рис. 11. Результаты моделирования процесса Fig. 11. Results of modeling the process of
идентификации общего коэффициента передачи электропривода К с компенсацией возмущений Мс и/ 1 - номинальное значение Щ); 2 - оценка
К ) при влиянии внешних возмущений Мс и/
identification of the total transmission coefficient of the electric drive K with compensation of disturbances and f: 1 - nominal value of K(t); 2 -evaluation under the influence of external disturbances and f
График 2 на рисунке 11 характеризует динамические выбросы вначале процесса идентификации с перерегулированием 36% на 0,0225 с, а также при набросе нагрузки через
0,1 с от начала моделирования с амплитудой Kmax около 4 о.е. Процесс идентификации устанавливается примерно через 0,2 с и установившееся значение общего коэффициента передачи электропривода ¿уст » Knom = 2,167. Флуктуации от возмущающего воздействия f
благодаря фильтрации становятся незначительными и вносят пренебрежимо малую погрешность около 0,5 %. В этой связи можно считать, что предложенные выше решения позволили получить инвариантный как к параметрическим, так и к координатным возмущениям алгоритм идентификации общего коэффициента передачи электропривода постоянного тока.
Выводы
В настоящей статье проведены разработка и исследование алгоритма идентификации общего коэффициента передачи электропривода постоянного тока с применением инверсной модели, не зависящей от идентифицируемого параметра. Алгоритм реализован с применением градиентного метода вычисления интеграла от квадрата невязки с использованием в качестве эталонной модели объекта исследования инверсной модели, а также модели чувствительности по общему коэффициенту передачи электропривода K. Имитационное моделирование в программной среде MatLab показало, что разработанный алгоритм инвариантен к параметрическим возмущениям и изменение параметров исследуемого объекта, выраженных через постоянные времени элементов электропривода, не оказывает влияния на результат идентификации K, тогда как вариации коэффициентов передачи элементов электропривода, являющихся сомножителями K, идентифицируются в виде оценки ¿уст (t) с высокой точностью. При исследовании влияния координатных
возмущений была выявлена высокая чувствительность алгоритма идентификации к данному виду возмущающих воздействий, причем максимальная погрешность вычисления оценки к (t) обусловлена моментной составляющей ошибки, которая алгебраически
суммируется с сигналом рассогласования Ам (t) и, соответственно, приводит к
неконтролируемому изменению невязки e(t) и, в итоге, вычисляемой оценки ¿уст (t). Для
минимизации негативного влияния момента сопротивления Мс предложен способ
формирования компенсирующего воздействия uc (t), полученного с применением датчиков
тока якоря ia (t) и скорости вращения двигателя ra(t) . В результате невязка становится
функцией трех переменных s(t)= f (uin, ujm, uc) и сигнальная компенсация момента
сопротивления обеспечивает определение установившегося значения общего коэффициента передачи электропривода к:уст (t) с погрешностью, не превышающей 0,5%. Влияние шумов
и наводок не вносит существенной ошибки в результат идентификации и при необходимости может быть снижено путем фильтрации сигнала рассогласования Ам (t) за
счет применения фильтра нижних частот. Таким образом, разработанный алгоритм идентификации позволяет с высокой точностью за малые промежутки времени произвести оценку общего коэффициента передачи электропривода постоянного тока в условиях текущей эксплуатации. Алгоритм не требует сложного измерительного оборудования и больших вычислительных мощностей и реализуется на основе стандартно применяемых в электроприводе датчиков. Процедура идентификации может быть автоматизирована путем получения соответствующего цифрового алгоритма и его реализации в виде рабочей программы микропроцессора. Результаты идентификации могут быть использованы как для анализа и контроля функционирования электропривода, так и для построения адаптивных параметрически инвариантных систем управления электроприводов с применением беспоисковых градиентных алгоритмов.
Литература
1. Calliess J., Roberts S.J., Rasmussen C.E., Maciejoswki J. Adapted Constant Kinky Inference for nonparametric regression and model-reference adaptive control. Automatica, 2020.
2. Malev N.A., Mukhametshin A.I. Pogoditsky O.V. and Mwaku W.M. Method of analysis and monitoring of the electromechanical converters parameters based on a linear integral criterion using sensitivity models. International Scientific and Technical Conference Smart Energy Systems 2019 (SES-2019). V. 124. Доступно по: https://doi.org/10.1051/e3sconf/201912402005.
3. Girish Chowdhary, Hassan A. Kingravi, Jonathan How, and Patricio A. Vela. Nonparametric adaptive control of time-varying systems using Gaussian processes. In American Control Conference (ACC), 2013.
4. Saushev A.V. Solution of problems of parametric optimization and control of electric drives state based on information about operability area boundary // IOP Conference Series: Materials Science and Engineering. 2018. №5. V.327.
5. Deisenroth M.P., Fox D., and Rasmussen C.E. Gaussian processes for data-efficient learning in robotics and control. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2015.
6. Малёв Н.А., Погодицкий О.В., Малацион А.С. Метод формирования g-таблиц для автоматизированного контроля параметров электромеханических преобразователей с применением линейного интегрального критерия. Известия высших учебных заведений. Проблемы энергетики. 2020. № 22(2). С. 86-97.
7. Furtat I., Fradkov A., Tsykunov A. Robust synchronization of linear dynamical systems with compensation of disturbances // Int. J. Robust and Nonlinear Control. 2014. V. 24. No. 17. -P. 2774-2784.
8. Ostroverkhov M., Pyzhov V., Korol S. Control of the Electric Drive under Conditions of Parametric Uncertainty and Coordinates' Interrelation // International Conference on Modern Electrical and Energy Systems (MEES). 2017.
9. Кучер Е.С., Комазенко М.А., Ромащенко А.И. Синтез систем векторного управления малочувствительных к изменениям параметров асинхронного электропривода // Доклады академии наук высшей школы российской федерации. 2017. № 2(35). С.61-72.
10. Kovela I., Ivaniuk O., Vlakh-Vyhrynovska G., Stakhiv R. Параметрический синтез комбинированных автоматических систем регулирования с цифровыми пид-регуляторами // Восточно-европейский журнал передовых технологий. 2017. № 2(90). С.37-34.
11. Томчина О.П., Кондратова М.Е., Горохов М.М., Журов К.В. Адаптивный регулятор с нелинейной эталонной моделью // Научные достижения современной науки: новация, история, действительность, перспективы и практика реализации. СПб.: Изд-во «КультИнформПресс», 2017. C. 206-208.
12. Анисимов А. А. Параметрическая оптимизация электромеханических систем с регуляторами и наблюдателями состояния // Вестник Ивановского государственного энергетического университета. 2016. № 2. С. 21-26.
13. Рутковский В.Ю., Глумов В.М. Особенности динамики адаптивной системы управления с нелинейной эталонной моделью. Автомат. и телемех., 2017. № 4. С.92-105.
14. Anisimov A.A., Kotov D.G., Tararykin S.V. Analysis of Parametric Sensitivity and Structural Optimization of Modal Control Systems with State Controllers // Journal of Computer and System Sciences International. 2011. V. 50. 5. pp. 698-719.
15. Сольницев Р.И., Каримов А.И., Каримов Т.И. Синтез цифровых регуляторов гироскопических командных приборов // Гироскопия и навигация. 2017. Т. 25. № 1 (96).
16. Саушев А.В., Бова Е.В. Показатели качества и критерии оптимальности при структурно-параметрическом синтезе автоматизированных электроприводов // Вестник Государственного университета морского и речного флота имени адмирала С.О. Макарова. 2019. Т.11. № 2. С. 380-395.
17. Pshikhopov V., Medvedev M., Kupovikh G., et al. Position Control of Mobile Robots With Multi-Contour Adaptation, 2018 IEEE International Conference on Autonomous Robot Systems and Competitions (Icarsc), eds.Costelha H., Calado J., Bento L., Lopes N., Oliveira P., IEEE, 2018, 54-59.
18. Афанасьев А.Ю., Макаров В.Г., Загирова В.Н. Идентификация параметров и частоты вращения ротора асинхронного двигателя на основе функций чувствительности // Математические методы в технике и технологиях-ММТТ. 2020.Т.6. С.128-131.
19. Малёв Н.А., Погодицкий О.В., Цветкович А.М. Особенности применения теории чувствительности для анализа влияния параметрических возмущений на динамические свойства электромеханических преобразователей. Известия высших учебных заведений. Проблемы энергетики. 2019;21(6):101-110.
20. Малёв Н.А., Мухаметшин А.И., Погодицкий О.В., Чичков П.В. Сравнительный анализ аппаратно-программного обеспечения метода контроля функционирования электромеханических преобразователей постоянного тока. Известия высших учебных заведений. Проблемы энергетики. 2020;22(5):142-154.
21. Сташинов Ю.П. К вопросу о настройке системы управления электропривода постоянного тока на модульный оптимум. Ч. 1 // Электротехника. 2016. № 1. С. 2-7.
22. Гарькина И.А., Данилов А.М., Тюкалов Д.Е. Сложные системы: идентификация динамических характеристик, возмущений и помех // Современные проблемы науки и образования. 2015. № 1. Ч. 1. С. 88.
Авторы публикации
Малёв Николай Анатольевич - ст. преп. кафедры Приборостроение и мехатроника, Казанский государственный энергетический университет.
Погодицкий Олег Владиславович - канд. техн. наук, доцент кафедры Приборостроение и мехатроника, Казанский государственный энергетический университет.
Чиляева Маргарита Романовна - студент, Казанский государственный энергетический университет.
Имамиев Айрат Рустемович - инженер группы сопровождения проектов КРиС ФЛ АО «Сетевая компания» - «Дирекция строящихся объектов».
References
1. J. Calliess SJ, Roberts CE, Rasmussen J, Maciejoswki. Lazily Adapted Constant Kinky Inference for nonparametric regression and model-reference adaptive control. Automatica, 2020.
2. Malev NA, Mukhametshin AI, Pogoditsky OV and Mwaku WM. Method of analysis and monitoring of the electromechanical converters parameters based on a linear integral criterion using sensitivity models. International Scientific and Technical Conference Smart Energy Systems 2019 (SES-2019).V. 124. Available At: https://doi.org/10.1051/e3sconf/201912402005.
3. Girish Chowdhary, Hassan A. Kingravi, Jonathan How, and Patricio A. Vela. Nonparametric adaptive control of time-varying systems using Gaussian processes. In American Control Conference (ACC), 2013.
4. Saushev AV. Solution of problems of parametric optimization and control of electric drives state based on information about operability area boundary. IOP Conference Series: Materials Science and Engineering. 2018;5:327.
5. M. P. Deisenroth, D. Fox, and C. E. Rasmussen. Gaussian processes for data-efficient learning in robotics and control. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2015.
6. Malev NA, Pogoditsky OV, Malacion AS. Q-tables formation method for automated monitoring of electromechanical converters parameters with application of linear integral criterion. Power engineering: research, equipment, technology. 2020;22(2):86-97.
7. Furtat I, Fradkov A. Tsykunov A. Robust synchronization of linear dynamical systems with compensation of disturbances. Int. J. Robust and Nonlinear Control..2014;24:17:2774-2784.
8. Ostroverkhov M, Pyzhov V, Korol S. Control of the Electric Drive under Conditions of Parametric Uncertainty and Coordinates' Interrelation. International Conference on Modern Electrical and Energy Systems (MEES). 2017.
9. Kucher ES, Komazenko MA, Romashchenko AI. Synthesis of systems of vector control of induction motor drive parameters low-sensitive to changes. Reports of the Academy of Sciences of the Higher School of the Russian Federation. 2017;2(35):61-72.
10. Kovela I, Ivaniuk O, Vlakh-Vyhrynovska G, Stakhiv R. Parametric synthesis of combined automatic regulating systems with digital pid-controllers. Eastern european journal of advanced technologies. 2017;2(90):С.37-34.
11. Tomchina OP, Kondratova ME, Gorokhov MM, Zhurov KV. An adaptive controller with a nonlinear reference model - in the collection: Scientific achievements of modern science: novation, history, reality, prospects and implementation practice: Collection of scientific articles based on the results of an international scientific and practical conference. St. Petersburg: Publishing House «KultInform-Press», 2017. P. 206-208.
12. Anisimov AA. Parametrical optimization of regulators and state observers in electromechanical systems. VestnikIGEU. 2016;2:21-26.
13. Rutkovskii VYu. , M. Glumov V. Dynamics peculiarities of an adaptive control system with nonlinear reference model. Autom. Remote Control. 2017;78:4:654-665.
14. Anisimov AA, Kotov DG, Tararykin SV. Analysis of Parametric Sensitivity and Structural Optimization of Modal Control Systems with State Controllers. Journal of Computer and System Sciences International. 2011;50(5):698 - 719.
15. Solnitsev RI., Karimov AI., Karimov TI.. Synthesis of digital regulators of gyroscopic command appliance. Gyroscopy and navigation. 2017;25(1):(96).
132
16. Saushev AV, Bova EV. The quality indicators and optimality criterion at the structural-parametric synthesis of automated electric drives. Bulletin of the State University of Marine and River Fleet named after admiral S.O. Makarov. 2019;2:380-395.
17. Pshikhopov V, Medvedev M, Kupovikh G, Shibanov V. Position Control of Mobile Robots With Multi-Contour Adaptation, 2018 IEEE International Conference on Autonomous Robot Systems and Competitions (Icarsc), eds. Costelha H., Calado J., Bento L., Lopes N., Oliveira P., IEEE, 2018;54-59.
18. Afanasiev AY, Makarov VG, Zagirova VN. Identification of parameters and rotor speed of the induction motor with the sensitivity functions. Mathematical methods in engineering and technology-MMTT. 2020;6:128-131.
19. Malev NA, Pogoditsky OV, Cvetkovich AM. Features of application of sensitivity theory for analysis of influence of parametric disturbances on dynamic properties electromechanical converters. Power engineering: research, equipment, technology. 2019;21(6):101-110.
20. Malev NA, Mukhametshin AI, Pogoditsky OV, Chichkov PV. Comparative analysis of hardware and software of the monitoring method of fun ctioning electromechanical DC converters. Power engineering: research, equipment, technology. 2020;22(5):142-154.
21. Stashinov YP. On the issue of control system adjustment of a direct current drive on the modular optimum. Pt . Russian Electrical Engineering. 2016;1:2-7.
22. Garkina IA, Danilov AM, Tyukalov DE. Complex systems: identification of dynamic characteristics, disturbances and interferences. Modern problems of science and education. 2015;1. Pt 1:88.
Authors of the publication
Nikolai A. Malev - Kazan State Power Engineering University. Oleg V. Pogoditsky - Kazan State Power Engineering University. Margarita R. Chilyaeva - Kazan State Power Engineering University
Ayrat R. Imamiev - engineer of the projects support group for KRiS branch JSC «Grid Company» - «Directorate of facilities under construction».
Получено 07.04.2021 г.
Отредактировано 06.10.2021 г.
Принято 01.12.2021 г.