УДК 681.5.11
Синтез следящего электропривода азимутальной оси телескопа с эталонной моделью в контуре положения
Демидова Г.Л., ассист., Ловлин С.Ю., Цветкова М.Х., аспиранты
Предлагается структура алгоритма управления следящим электроприводом с эталонной моделью в контуре положения. Приведены результаты структурно-параметрической оптимизации, направленной на устранение «шагающего» режима.
Ключевые слова: система управления, регулятор, эталонная модель, шагающий режим, сухое трение.
Synthesis of Follow-up Electric Drive of telescope’s azimuth axis with reference model in position contour
G.L. Demidova, Assistant, S.Yu. Lovlin, M.Kh. Tsvetkova, Post Graduate Students The article presents the structure of control algorithm for the follow-up electric drive with reference model
in the position contour. The authors present the results с mode eliminating.
Keywords: control system, regulator, reference model,
Введение. Создание высокоточных, надежных и долговечных систем автоматического управления, работающих в жестких условиях эксплуатации, в большинстве случаев возможно только на основе реализации цифровых безре-дукторных следящих систем (ЦБСС) [1-3]. При этом наилучшая гибкость управления механическим движением достигается за счет использования высокомоментного синхронного двигателя с самокоммутацией (вентильного двигателя). Такой электропривод обычно называется безредукторным или прямым электроприводом с вентильным двигателем [3].
По своим точностным и динамическим характеристикам, а также по показателям надежности и долговечности работы приборы с безредукторными системами существенно превосходят приборы с редукторными приводами, однако преимущества последних по массогабаритным показателям и энергопотреблению в отдельных случаях делают их применение единственно возможным [1-3]. Так, мощности отдельных образцов прямого электропривода достигают сотен киловатт, а типовые значения погрешностей углового и линейного позиционирования - единиц секунд и микрометров, что с полным правом позволяет назвать такой привод прецизионным (высокоточным). Сфера применения прецизионного прямого электропривода включает машиностроение, робототехнику, электронную индустрию, измерительную технику и специальные технические устройства [3].
К числу актуальных проблем современной электромеханики и преобразовательной техники относится проблема управления прецизионными следящими электроприводами для оптико-механических систем и, в частности, для комплексов позиционирования и слежения.
При управлении прецизионными следящими электроприводами, такими как, например,
f structural and parametric optimization, aimed at step
step mode, dry friction.
электромеханическая система телескопа, на инфранизких скоростях (скорости, сравнимые со скоростью вращения Земли) в условиях не-жесткости оси точность слежения за сигналом задания уменьшается [5-7].
Наличие сухого трения тоже отрицательно влияет на точность реализации заданных программных траекторий систем управления. Это связано с тем, что момент, создаваемый сухим трением, описывается функцией, разрывной относительно обобщенной скорости (сила трения изменяется скачкообразно при изменении направления движения) [4]. Как правило, при движении с инфранизкой скоростью этот факт приводит к «шагающему» режиму, колебаниям системы, вызванным реакцией на частые скачки возмущающего воздействия. Условие отсутствия «шагающего» режима - движение электропривода со скоростью, средняя величина которой превышает амплитуду ее колебаний. Таким образом, для повышения точности таких систем вводится нижняя граница скоростей.
Ниже предлагается алгоритм управления, позволяющий расширить нижнюю границу скоростей.
Постановка задачи. Наличие «шагающего» режима свойственно электроприводам, движущимся на инфранизких скоростях, вплоть до нулевых скоростей (режим позиционирования). Стандартные методы повышения точности путем увеличения порядка астатизма не могут бороться с нелинейными возмущениями типа сухого трения. Кроме того, именно наличие астатизма по возмущающему воздействию способствует возникновению «шагающего» режима, а в статических системах он отсутствует.
Одно из наиболее простых решений -создание таблицы моментов трения. Значения моментов трения в зависимости от угла поворота записываются в таблицу и затем исполь-
зуются в качестве добавки к управлению с соответствующим знаком. Однако это требует огромных вычислительных затрат (необходимо где-то хранить эту таблицу). Кроме того, возможно изменение трения при изменении температуры с течением времени из-за механических воздействий и изменения свойств смазки.
Метод адаптивной компенсации сухого и вязкого трения [4] позволяет эффективно решить проблему только в случае отсутствия активного трения на валу.
Метод, рассматриваемый в данной работе, ориентирован на структурно-параметрическую оптимизацию системы управления по возмущающему воздействию. Кроме стандартных контуров подчиненной системы регулирования [5] в систему управления вводится дополнительный контур положения с эталонной моделью (рис. 1).
Система управления с эталонной моделью в контуре положения обладает одним привлекательным свойством: введение дополнительного контура никак не влияет на передаточную функцию системы по задающему воздействию при условии того, что эталонная модель будет соответствовать передаточной функции контура положения. Действительно, при отсутствии возмущающего воздействия на входе регулятора контура с эталонной моделью WРЭм(p) будет всегда ноль (рис. 1). Таким образом, можно добиться различных передаточных функций контура положения по возмущению при одинаковых передаточных функциях по заданию.
Итак, задача состоит в том, чтобы найти такую настройку контура с эталонной моделью, что влияние сухого трения будет минимальным.
Настройка основных контуров. Отсутствие датчиков на второй массе в прецизионных следящих электроприводах для оптикомеханических систем с нежесткими осями не
позволяет применять традиционные методы настройки в связи с возможностью возникновения механического резонанса [5]. Поэтому оптимизация системы осуществляется настройкой контура положения на «симметричный оптимум» на основании выражений, представленных в [5].
В качестве примера объекта управления был выбран вентильный двигатель (ВД) ДБМ 1200-1200 с двухмассовой нагрузкой с моментами инерции первой массы ^ = 35,28 кг м2, второй массы и2 = 94,34 кг м2, с жесткостью Сі2 = 1,874 Н м/рад. На двигатель действует момент сопротивления на первой массе МС1 = 40 Нм и момент сопротивления на второй массе МС2 = 40 Нм. В качестве модели трения выбрана модель с эффектом Штрибека. При настройке основных контуров двухмассовая модель (рис. 1) представляется в виде эквивалентной одномассовой [5], поэтому моменты сопротивления масс далее при настройке контура с эталонной моделью будут рассматриваться как эквивалентное возмущающее воздействие. В таблице - сопротивление фазной обмотки ВД; Тф - электрическая постоянная времени фазной обмотки ВД; СЕ - коэффициент по про-тиво-ЭДС; КСТ - статический коэффициент передачи широтно-импульсного преобразователя (ШИП); ТЗ - некомпенсируемая постоянная времени запаздывания ШИП; КдТ, КдС, КдП - коэффициенты передачи датчиков тока, скорости и положения соответственно; ю0 - полоса пропускания частот в скоростной подсистеме, вычисляемая по следующей формуле:
С|2(^1 + ^2)
^2
Рис. 1. Структурная схема системы с эталонной моделью в контуре положения: азад - сигнал задания; а1 и а2 - углы первой и второй масс; а1 и а2 - скорости первой и второй масс; СМ - коэффициент по моменту
Параметры регуляторов ^ основных контуров
ПИ-регулятор контура тока Коэффициент передачи П-регулятора К КфТф ПТ 2Г3КстКдт • 1,5
Время изодрома И-регулятора ТИТ _ ТФ
Контур скорости Коэффициент передачи П-регулятора внутреннего контура ,У К ДТ “о(^1 + ^2) КПС к с к дссе
Время изодрома И-регулятора внешнего контура Т _2 1 ИС _ “о
ПИ-регулятор контура положения Коэффициент передачи П-регулятора £ 1 _ ЕЕ
Время изодрома И-регулятора Т _ 8 1 ИП “о
На основании выбранной настройки эталонная модель будет иметь следующую передаточную функцию:
4Тцрр +1
8Тш3 Р3 + 8ТР р2 + 4Т№Р +1
^ (р) =
8|цр р -8'№ где Тцр = 2Тцз = 2/ ^0■
Синтез и настройка регулятора контура с эталонной моделью. Синтез передаточной функции по возмущению■ Передаточная функция контура положения по возмущающему воздействию (рис. 1) в общем случае будет иметь следующий вид:
Мво
а (Р) =
ТИСКДТ КПСКДС
рТитР +1)
(2Гц52р + 2Гцзр + 1)р + \WPЭM + Кпп ТпР +1) 1 Кш.
ДС
(1)
Тогда при подстановке значений из таблицы в (1) получим
64СЕ „2(т р , 1)
--4еР (ТфР + 1)
а ( Р)-
64 „4 , 64 „3 , 32 р2 + 8 +1 + ^ 4Кдп
Тр +-----+ 1 + ИРЭМ-------
р +-
;Р +-
“0
“оК ДС
где + ^2 ■
Реакцию скорости на возмущение можем узнать, добавив дифференциальное звено в (1):
64СЕ 3
—4Е Р (ТфР +1)
Ивозм.“( р)
4 I
“о
64 4 64 3 32 2 8 . ...
—ГР +—з Р +—2;Р + — Р +1 + Ин “о “о “о “о
4К
дп
I/
“оК ДС
(2)
Рассмотрим возможные варианты передаточной функции (2) в зависимости от 1ЛРЭМ. Анализировать их свойства будем по реакции на скачок возмущающего воздействия при наличии момента сухого трения. Результаты приведены на рис. 2, в качестве выхода был выбран угол, а не скорость, так как именно он является конечной настраиваемой величиной.
Графикам 1-4 (рис. 2) соответствует следующая передаточная функция контура положения с эталонной моделью:
И
А р) _-
КмР (ТфР +1)
(Т 12р2 + 2^р + 1)(Т 2 Р + 1)(Т 3 р +1)
где % = 1, Т2 = Тф - для графика 1; % < 1, Т2 = Тф -для графика 2; % < 1, Т2 > Тф, Т1 < Тф - для графика 3; % = 1, Т2 > Тф, Т1 > Тф - для графика 4. Графику 5 соответствует передаточная функция
А р) _■
КмР (ТфР +1)
(Т1Р + 1)3(Т 2 р + 1)
Очевидно, что наиболее приемлемыми вариантами являются передаточные функции, соответствующие графикам 1 и 4, т.е. передаточная функция (2) должна обладать следующими свойствами: отсутствием комплексных корней в характеристическом уравнении и минимальным порядком астатизма:
А Р) = -
КмР (ТфР +1)
—. (3)
(Т|Р + 1)(Т 2 Р + 1)(Т3 Р + 1)(Т4 Р +1) Синтез РегулятоРа эталонной модели. Для настройки системы управления с передаточной функцией (2) необходим следующий регулятор:
ИРЭМ ( р) = К 6р 4 + К 5р3 + К 4р2 + К3р + К2 +-~.
Р
Передаточная функция (2) примет вид
где
0( Р) =
“о ^ЕК1
К6 5 64К5 4 64К4 3 32К3 2 8К
—6 р5 +—-г-5- р4 +—р3 +—^ р2 +-----------------2 +1
К “о4к1 “о3К1 “о2К1 “оК
(4)
К1= 4КДПКі +1; к2= КДПК2 +1;
“оКДС Кгт ггЮп К
2К
ДС
Кгт пЮп Кл
к пп“г\Къ г к пКа
К3= дп о 3 +1; К4 = ДП о 4 +1;
16К
ДС
8Кдс
к5 —кДП ю°3к5 +1.
16кдс
Коэффициент К6 обеспечит относительную степень передаточной функции (4), равную единице, придаст ЛАЧХ передаточной функции наклон -2° Дб/дек в области высоких частот и позволит регулировать коэффициент усиления по возмущающему воздействию. Однако четыре последовательно-соединенных дифференциальных звена в цифровых прецизионных системах создают существенные шумы дискретизации.
Такая же проблема и с коэффициентами К5 и К4, поэтому К6 = К5 = К4 = 0, а коэффициент усиления по возмущающему воздействию в области высоких частот будет постоянным и равным 0ЕТФ/и1. С учетом выше сказанного, а также формул (3) и (4) получим
1111 «ч
— + — + — + — — ®°. (5)
Т Т2 Т3 Т4 ° ^
Одна из постоянных времени Т^ Т2, Т3, Т4 должна быть меньше ТФ. Это одно из необходимых условий минимальности коэффициента усиления по возмущающему воздействию на всем диапазоне частот. Согласно (5), настраиваемые постоянные времени ограничены конечной величиной, и данное условие неосуществимо при ТФ< ю0"1. Тогда одна из постоянных времени (пусть Т4) должна быть минимальной, а максимальный коэффициент усиления будет
на частоте 1/Т4 равным 0ЕТ4/и1, т.е.
Т4<<(Т1,Т2,Т3). В таком случае
111 ®°К3
— + — + — — ° 3 , (6)
Т Т2 Т3 2
Т1Т2Т3 _
64
“К'
(7)
Согласно вышеописанному условию и уравнениям (6) и (7), получаем Т = Т2 = Т3. Очевидно, что К’1 = 1 - наиболее оптимальный вариант для выполнения Т4<<(Ъ,Т2,Т3), т.е. К = о.
Итак, для настройки контура с эталонной моделью необходим ПД-регулятор со следующей передаточной функцией:
ИРЭМ (р) _ кд р + КП .
НастРойка РегулятоРа эталонной модели■ Коэффициенты регулятора рассчитываются исходя из вышеописанных условий и формулы (4):
К Г Кп _■
ДС
К
дп
КД _■
4К
ДС
Кдп ®°
Результаты эксперимента. Результаты моделирования реакции на линейное задающее воздействие с инфранизкой скоростью 1''/с изображены на рис. 3. В системе подчиненного регулирования без эталонной модели отчетливо видны колебания, превышающие среднюю скорость в несколько раз. В системе с эталонной моделью эти колебания минимальны и обусловлены только лишь конечной разрядностью датчиков. Реверсивный режим работы электропривода не рассматривался, так как он возможен только в режиме переброса [2], а в этом режиме высокие точности движения электропривода не требуются.
1 1
Мі
1 і \ М2 і
I и и: и 1 1
\ Мі у
1 і І
о 15 1 1,5 2 2.5 3 0 0.5 1 1.5 2 25 3
а) б)
Рис. 3. Отработка сигнала задания с постоянной скоростью 17с традиционной системой подчиненного регулирования (а) и системой с эталонной моделью в контуре угла (б)
Заключение Список литературы
Предложенная структура с эталонной моделью в контуре положения настроена на уменьшение автоколебаний и устранение «шагающего» режима при максимальном быстродействии системы. Уникальность предложенного метода настройки систем управления заключается в том, что можно получить желаемую передаточную функцию не только по задающему воздействию, но и по возмущающему. Благодаря данной настройке была существенно снижена ошибка слежения и расширена нижняя граница скоростей. Данная система управления используется в прецизионных следящих электроприводах с нежесткими осями, в опорно-поворотных устройствах телескопов траекторных измерений.
1. Электропривод с вентильными двигателями: Докл. науч.-практич. семинара. - М.: Изд. дом МЭИ, 2оо7.
2. Томасов В.С., Овчинников И.Е., Егоров А.В. Энергоподсистема большого алтайского телескопа траекторных измерений // Известия тульского государственного университета. - 2оіо. - Вып. 3. - Ч. 3. - С. 216-222.
3. Балковой А.П., Цаценкин В.К. Прецизионный электропривод с вентильными двигателями. - М.: Изд. дом МЭИ, 2оіо.
4. Стещенко Е.В. Адаптивная компенсация сухого и вязкого трения в системах автоматического управления // Вестник ЧитГУ. - 2оо8. - Т. 48. - № 3. - С. 123-128.
5. Толмачев В.А. Синтез следящего электропривода оси опорно-поворотного устройства // Изв. вузов. Приборостроение. - 2оо8. - Т. 51. - № 6. - С. 68-72.
6. Ключев В.И. Теория электропривода: учеб. для вузов. - М.: Энергоатомиздат, 1997.
7. Денисов К.М., Жданов И.Н. Оптимизация каналов скорости и положения прецизионных следящих электроприводов // Изв. вузов. Приборостроение. - 2оо8. -Т. 51. - № 6. - С. 45-51.
Работа выполнена ПРи финансовой поддеРЖке ФедеРального агентства по обРазованию по НИР НК-430П ГК №П2479 от 19.11.2009 г. (меРопРиятие №1.2.2) по напРавлению «Многофункциональное пРибоРостРоение для пРомышленных систем упРавления» по пРоблеме «РазРаботка систем упРавления движением пРецизионных электРопРиводов пРибоРных комплексов нового поколения» в Рамках ФЦП «Научные и научно-педагогические кадРы инновационной России на 2009-2013 г».
Демидова Галина Львовна,
Санкт-Петербургский государственный университет информационных технологий, механики и оптики, ассистент кафедры электротехники и прецизионных электромеханических систем, телефон (812)233-83-36, e-mail: demidofffa@mail.ru
Ловлин Сергей Юрьевич,
Санкт-Петербургский государственный университет информационных технологий, механики и оптики, аспирант кафедры электротехники и прецизионных электромеханических систем, e-mail: s.lovlin@ets.ifmo.ru
Цветкова Мадина Хасановна,
Санкт-Петербургский государственный университет информационных технологий, механики и оптики, аспирант кафедры электротехники и прецизионных электромеханических систем, телефон (812)233-83-36, e-mail: madina1986@bk.ru