CC BY
66
ИДЕНТИФИКАЦИЯ И НАСТРОЙКА СИСТЕМЫ УПРАВЛЕНИЯ ЭЛЕКТРОПРИВОДА ...
УДК 681.5.015
ИДЕНТИФИКАЦИЯ И НАСТРОЙКА СИСТЕМЫ УПРАВЛЕНИЯ ЭЛЕКТРОПРИВОДА АЗИМУТАЛЬНОЙ ОСИ ТЕЛЕСКОПА
С.А. Тушев
Описывается способ получения математической модели электропривода азимутальной оси телескопа с помощью идентификации его параметров. Полученная модель используется для настройки прецизионной системы автоматического управления объектом, работа которой проверяется на реальном устройстве. Теоретические расчеты в значительной степени совпадают с практическими результатами, что показывает достоверность данного метода построения системы управления.
Работа представляет интерес для студентов, аспирантов, работников электротехнических специальностей, а также специальностей и направлений подготовки, связанных с теорией автоматического управления. Ключевые слова: идентификация, электропривод, настройка, прецизионный электропривод, телескоп.
Введение
В настоящее время проблемы, связанные с созданием математических моделей технических объектов, технологических процессов, экономики и живой природы, формируют одно из основных направлений науки и техники - моделирование. Это объясняется тем, что математические модели объектов широко применяются на всем периоде работы с ними - от разработки до эксплуатации.
Система управления - неотъемлемая часть современного высокотехнологичного комплекса. Ее создание подразделяется на основные этапы: получение математической модели объекта, настройка системы автоматического управления (САУ) для этой модели и проверка работоспособности на реальном объекте. Самым сложным является первый этап, так как он подразумевает структурную и параметрическую идентификацию.
Существует множество методов решения задачи идентификации. В настоящей работе будет рассмотрена параметрическая идентификация модели электропривода азимутальной оси телескопа траек-торных измерений типа «Сажень» с использованием пакета программ МЛТЬЛБ и создание системы управления этим объектом. Выбор данного способа обусловливается тем, что он наиболее точен, так как все коэффициенты передаточной функции объекта определяются отдельно друг от друга. Это позволяет искать передаточную функцию объекта в любом виде - от простейшего апериодического звена до сложных звеньев высоких порядков, в зависимости от требований к точности определения модели [1].
Идентификация параметров объекта
Известно, что реальная установка содержит вентильный двигатель (ВД), однако анализ переходных характеристик [2, 3] показал, что модель ВД с достаточной степенью точности соответствует модели двигателя постоянного тока (ДПТ), который описывается системой уравнений (1), учитывающей двух-массовый характер механической нагрузки. Двухмассовость обусловлена тем, что азимутальная ось телескопа в реальности не является абсолютно жесткой. По этой причине она условно разделяется на две массы, между которыми присутствует нежесткость. са аг
= ^I „ +т.а2 _к • а,
аг 1 1 1 г 1 г 2 г 2 г
^ = а (1) аг а
1 . Ce- Q + 1 —u,
T 1f L 1 L
N -m Q1 _ J1 c12 —— a. + J1 m Q2 + J1 C12 — a 2 - J1 Mr1 J1
= «1,
L= m Q J2 c12 +—— a J 2 - m q2 J2 2 - ^ a 2 J2 2 Mr 2 J2 '
= «2-
at
da2
at
Здесь K - коэффициент, учитывающий величину момента нагрузки от кабельного перехода; Mr1 и Mr2 - моменты трения на 1-ой и 2-ой массе соответственно; Ce и Cm - постоянные электродвигателя; J1 и J2 - приведенные моменты инерции; m - коэффициент демпфирования; с12 - коэффициент нежесткости связи между массами; Tf- электрическая постоянная времени.
На исследуемый объект воздействуют внешние возмущения - моменты трения и момент сопротивления, создаваемый кабельным переходом. Они определяются следующими выражениями:
Mc (b) = b • MJ, b = sign (Q12);
M (ai revive ) = J1 • K • «1 revive -
Кабельный переход, являющийся конструктивной особенностью опорно-поворотного устройства телескопа, создает момент нагрузки, пропорциональный относительному углу поворота оси o.i-relative и действует только на первую массу.
С помощью пакета программ MATLAB/System Identification Toolbox [4] и MATLAB/Simulink осуществляется процедура идентификации параметров модели. Вначале регистрируется реакция разомкнутого объекта на определенное входное воздействие. Затем задается математическое описание системы в матрицах состояния (2).
X(t) = Ax(t) + Bu (t) + Ke(t), • y(t) = Cx(t) + Du(t) + e(t), (2)
x(0) = x0.
Начальные значения коэффициентов могут выбираться произвольно, однако рекомендуется принимать их близкими к реальным, так как это может сильно повлиять на результат идентификации. В данной структуре неизвестны следующие параметры: K, Mrl, Mr2, Ce, Cm, m. Остальные значения были определены опытным путем или были указаны в техническом задании. В качестве тестового сигнала выбран двусторонний меандр, так как при таком воздействии меняется знак скорости и могут быть учтены особенности, связанные с направлением вращения оси. оз, рад/с
Выход объекта
t, с
Рис. 1. Результат параметрической идентификации
После завершения выполнения задачи программа МАТЬАБ возвращает модель системы в уравнениях состояния с найденными коэффициентами. График сравнения реальной системы и полученной модели представлен на рис. 1.
Выходом модели и объекта является скорость вращения азимутальной оси. Неполное совпадение (около 87,5%) обусловливается тем, что структурная схема модели не до конца соответствует объекту управления, а также погрешностью определения параметров. Несмотря на эти неточности, проведенная идентификация позволяет синтезировать систему САУ для модели с ДПТ.
Настройка САУ электропривода
Создание системы управления предполагает множество решений при помощи различных типов регуляторов. Распространенным решением является система подчиненного регулирования координат объекта, таких как вращающий момент или ток двигателя, скорость, положение.
Из-за особенностей объекта управления контур скорости замыкается по первой массе, так как на вторую массу невозможно установить датчик положения в рассматриваемом объекте. Свое влияние оказывает и тот факт, что система имеет высокий коэффициент жесткости, что позволяет представить ее в виде одномассовой для упрощения процедуры настройки САУ. При расчетах используется суммарный момент инерции, равный сумме приведенных моментов инерции обеих масс.
Контур скорости настраивается на технический оптимум (ТО) с помощью пропорционально-интегрально-дифференциального регулятора с замедлением. Объект управления представляется как колебательное звено, а числитель регулятора компенсирует его. Следовательно, постоянная времени настроенной системы зависит только от постоянной времени фильтра Тц, входящего в структуру регулятора. Графики переходных процессов при скачке задающего воздействия представлены на рис. 2.
Реакции систем - эталонной и настроенной с ДПТ - на задание скорости имеют основные характеристики, представленные в табл. 1.
ИДЕНТИФИКАЦИЯ И НАСТРОЙКА СИСТЕМЫ УПРАВЛЕНИЯ ЭЛЕКТРОПРИВОДА
Рис. 2. Реакции модели с ДПТ, настроенной на ТО, и эталонной модели
Характеристика Эталон ДПТ Реальная установка
Время первого согласования, Ти 4,7 4,5 4,2
Перерегулирование, % 4,3 4,6 18
Таблица 1. Сравнение параметров реакций моделей и реальной установки на скачок задания скорости
Контур положения настраивается на симметричный оптимум (СО) с помощью пропорционально-интегрального регулятора для обеспечения второго порядка астатизма по управляющему воздействию. Графики переходных процессов представлены на рис. 3.
Реакции систем - эталонной и настроенной с ДПТ - на задание углового положения имеют основные характеристики, представленные в табл. 2. На рис. 4 представлены реакции настроенной реальной системы на скачок задающего воздействия по скорости и углу поворота соответственно.
Анализ данных, сведенных в табл. 1, 2, показывает, что значения параметров достаточно близки, что позволяет сделать вывод о целесообразности выбора такого способа синтеза САУ. В результате реальный объект удается настроить на постоянную времени скоростного контура, равную 5 мс, т.е. время переходного процесса в такой системе - порядка 30 мс, а в контуре регулирования положения - 60 мс.
0,0 0,1 0,2 0,3 и с
Рис. 3. Реакции систем с ДПТ и ВД, настроенных на СО, и эталонной модели
Характеристика Эталон ДПТ Реальная установка
Время первого согласования, Tu 3,1 2,9 3
Перерегулирование, % 43 50 34,5
Таблица 2. Сравнение параметров реакций моделей и реальной установки на скачок задания угла оо, рад/с а' рад
250 ty мс
100 200 300 400 500 600 700 U мс б
Рис. 4. Реакция реального объекта на скачок задания: скорости (а); угла поворота (б)
Заключение
Основным преимуществом использования параметрической идентификации в среде MAT-LAB/System Identification Toolbox является то, что есть возможность задавать любую структурную схему объекта управления - от простейшей до самой точной, учитывающей все возможные постоянные времени и внешние воздействия. Программный комплекс подбирает значения параметров модели объекта так, чтобы максимально приблизить смоделированную реакцию системы к реальной. Наглядно анализируя точность совпадения, можно решить, стоит ли дальше усложнять структуру модели или в данных условиях можно принять текущую.
Получив модель объекта с найденными параметрами и величинами возмущающих воздействий, можно построить систему управления в желаемом виде и проверить ее работу в MATLAB/Simulink. Завершающим этапом настройки системы является проверка ее работоспособности на реальном объекте. Если статические и динамические характеристики настроенной системы не удовлетворяют требованиям, процедура идентификации повторяется.
Предложенный метод настройки системы автоматического управления электропривода позволяет получить в короткие сроки достаточно высокие показатели качества переходных процессов, что дает возможность использовать его на начальном этапе настройки системы управления цифрового электропривода.
Литература
1. Семенов А.Д., Артамонов Д.В., Брюхачев А.В. Идентификация объектов управления: Учебное пособие. - Пенза: Изд-во Пенз. гос. ун-та, 2003. - 211 с.
2. Аракелян А.К., Афанасьев А.А. Вентильные электрические машины. - М.: Энергоатомиздат, 1997. -508 с.
3. Балковой А.П., Цаценкин В.К. Прецизионный электропривод с вентильными двигателями. - М.: Издательский дом МЭИ, 2010. - 328 с.
4. System Identification Toolbox User's Guide. - The MathWorks, Inc., 2008. - 528 p.
Тушев Сергей Александрович
Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант, [email protected]
а
