CC BY
22
Доклады БГУИР
2003
апрель-июнь
Том 1, № 2
ТЕХНОЛОГИИ
УДК 621.9.048
ПРОГРАММНАЯ КОМПЕНСАЦИЯ ПОГРЕШНОСТИ ПОЗИЦИОНИРОВАНИЯ
ЛИНЕЙНОГО ШАГОВОГО ПРИВОДА
В работе показано, что для повышения точности позиционирования ЛШД при разомкнутом управлении необязательно использовать внешние датчики. Система управления может быть построена на основе внутреннего датчика и цифровой обработки результатов изменения индуктивности в рабочем зазоре индуктор-статор. Приведены уравнения исследуемой системы. Реализован алгоритм коррекции погрешности позиционирования, позволяющий без использования внешних датчиков повысить точность позиционирования в 5-10 раз.
Ключевые слова: линейный шаговый привод, погрешность позиционирования, программная компенсация.
Применение линейных шаговых двигателей (ЛШД) при разомкнутом микрошаговом управлении в системах позиционирования, где требуется точное и равномерное движение, в настоящее время ограничено из-за их погрешности позиционирования [1-4]. Эту проблему можно существенно уменьшить путем управления ЛШД несинусоидальными токами.
Погрешность позиционирования как периодическая функция координаты вызывает две проблемы при движении ЛШД:
отклонение положения ЛШД от желаемого при позиционном управлении с амплитудой примерно от 5 до 50 мкм;
колебания индуктора с частотой, пропорциональной скорости движения, и амплитудой, резко возрастающей вблизи резонансной частоты ЛШД.
Ошибка позиционирования ЛШД Ах вызывается как физической природой используемых материалов, так и конструктивными и технологическими причинами. В общем случае погрешность позиционирования А х отражает все различия между реальным ЛШД и идеальной моделью с тяговым усилием, описываемым уравнением
где в стационарном состоянии Ах=0.
Эти отклонения реально возникают и зависят от положения индуктора, его скорости и амплитуд токов в фазах.
Причины погрешности позиционирования можно разделить на две группы [1, 5, 6]:
С.Е. КАРПОВИЧ., Ю.С. МЕЖИНСКИЙ, В В. ЖАРСКИЙ
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 25 мая 2003
Введение
(1)
дополнительные силы, действующие на индуктор, кроме тягового усилия (1), такие, как силы трения, нагрузки и остаточные силы из-за индукционных и вихревых токов, а также намагничивание статора;
физические взаимосвязи, которыми пренебрегли при проектировании и моделировании ЛШД, такие, как отклонение изменения магнитной проводимости в воздушном зазоре от синусоидальной формы, нелинейная зависимость магнитного потока от тока в обмотке вследствие насыщения магнитных материалов и рассеяния магнитного поля в воздушном зазоре, и, наконец, гистерезисное явление при перемагничивании полюсов индуктора и статора.
Принцип компенсации погрешности позиционирования
В упрощенной идеальной математической модели ЛШД [2, 3] тяговое усилие Р может быть рассчитано по выражению (1) для синусоидальных фазовых токов ^ и ¡2 :
11 () = 1тахС00 1, () = ^ах^0 1 , (2)
2п
где о =-.
Т
Полагая, что Р — это вся сила, действующая на индуктор, можно записать функцию положения индуктора от времени х(() в стационарном состоянии в виде
2пх
--о г = 0
Р
или
) = Т. (3)
В рассмотренной идеальной модели (3) индуктор, управляемый синусоидальными токами, движется в стационарном состоянии равномерно со скоростью V = т/Т .
Однако реальный ЛШД имеет погрешность позиционирования, которая зависит от положения х почти периодически с периодом зубцовой структуры статора т . В общем случае форма зависимости и амплитуда Ах являются также функцией амплитуды фазовых токов
1 тах .
Погрешность позиционирования ЛШД Ах в стационарном состоянии может быть значительно компенсирована изменением управления в фазах. Для описания измененных фазовых токов введем составляющую Ах в уравнение (3):
2п(х-Ах )-ог = о (4)
т
или
2п ( 2пААх Л л
--1 о 1 +-1 = 0. (5)
Погрешность позиционирования Ах будет компенсирована, когда реальный ЛШД управляется токами, определяемыми по выражениям:
2пАх Л
11 () = 1шсх:01 -■
т
¡2(() = *- . (6)
При этом выражения (5) и (3) становятся идентичными. Для того чтобы компенсировать
погрешность позиционирования А— , ток управления должен быть подан с опережением по фа-
2п Ах зе на-.
г
Компенсация погрешности позиционирования с помощью внешних измерений
Погрешность позиционирования можно компенсировать с использованием внешних устройств измерения положения, например, лазерного интерферометра. Измерение погрешности проводится для каждого микрошага к, массив погрешностей запоминается в памяти контроллера, и затем компенсация погрешности Ахк для каждого микрошага к по выражениям:
¡1 (к ) = 1тах
2П к А
I п т
¡2 (к ) = 1п
2П к-
I п т
(7)
где п — количество дроблений одного периода управления.
На рис. 1 показана погрешность позиционирования в пределах одного периода зубцо-вой структуры до и после компенсации.
Рис. 1. Компенсация погрешности позиционирования на к-м шаге: 1 — идеальная траектория; 2 — реальная траектория (с погрешностью)
Описанный способ позволяет уменьшить погрешность позиционирования примерно наполовину, так как несинусоидальные токи ЛТТТД рассчитывались по упрощенной модели.
Измерив реальную траекторию х- и зная длину микрошага хз, получим массив погрешностей позиционирования Ах-, где каждому микрошагу к соответствуют пара токов управления ^(к) и ¡2 (к);
начальное положение
кх,.
В этом случае реальная траектория определится по выражению
хк = Ахк + кхх.
(8)
Как показано на рис. 1, алгоритм поиска для микрошага к находит индекс "точного" положения г и помещает г -ю пару токов на к -е место. Это происходит последовательно для всех к при условии минимизации Ахг:
(1>*2 )к )г , Ах, = (9)
Результатом является массив индексов, элементы которого показывают последовательность фазовых токов для управления ЛШД с минимальной погрешностью.
Серьезным недостатком данного подхода является необходимость повторного измерения погрешности позиционирования при любом изменении внешних условий (амплитуд токов в фазах, температуры, силы нагрузки), что зачастую оказывается или невозможным, или попросту непрактичным.
Алгоритм калибровки и квазизамкнутого управления
Согласно конструкции и принципу действия идеального шагового двигателя, подача синусоидальных управляющих сигналов вызывает движение ЛШД без погрешности позиционирования. Стремление повысить тяговое усилие увеличением уровня токов в обмотках управления приводит к насыщению магнитных материалов ЛШД, и, как следствие, реальный ЛШД, имеющий идеальную геометрию, движется с погрешностью. Теоретически идеальное движение без погрешности возможно, однако управляющие токи должны для этого быть несинусоидальными. Повышение точности позиционирования ЛШД за счет изменения адаптации сигналов управления и обеспечивает разработанный алгоритм.
Целью данного алгоритма является уменьшение погрешности позиционирования ЛШД посредством минимизации его остаточного динамического тягового усилия за счет изменения формы токов управления с учетом результатов измерения индуктивностей обмоток. Алгоритм позволяет учесть насыщение магнитопроводов ЛШД несинусоидальный характер изменения магнитного сопротивления зубцовой структуры полюсных наконечников, модуляцию магнитного потока постоянного магнита в полюсных наконечниках.
Расчет остаточного тягового усилия ЛШД, необходимый для проведения коррекции формы управляющих токов в фазах, осуществляется в такой последовательности:
Магнитные потоки постоянного магнита в полюсах:
Фо = = Ф-ЛУ2Ь° , ФУ= = = , Ф, = Ф_А3, ь,, фт = Ф=ЛТ/ Ь,,
/4~ /4 /4 /4~ /4 /4
где Лк(х) = —^г — магнитное сопротивление в полюсах зубцовых структур, к = 0, Т4, Т2 , Ь(х) / 4 / 2
т/ /4.
Суммарные магнитные потоки в полюсах:
Ф = Ф + кг Ф , = Ф , — кг
^о ^о = о 'о' т/2 т/2= '
Ф/ = Ф/_ + 4/ К/ , Фт = Фт — 4/ т:
/4 /4~ /4 /4 /4 /4~ /4 /4
где Ь { — индуктивности в полюсах зубцовых структур; гк — токи в полюсных обмотках,
к = о т/ т/ т/ к = о, /4, /2, /4 •
Электромагнитные силы в полюсных наконечниках:
Р = -Ф2 дЛ I = от/ т/ 3т/ Р ф дх ' ' /4/2' /4 •
Остаточное динамическое тяговое усилие ШД:
Р = Р0 + Ру4 + Ру + Р3у4 .
Расчет погрешности позиционирования ШД:
= Р Т
6 = Рт 2п
Коррекция формы тока:
Т* т 51 I = I - е— .
дх
Блок-схема алгоритма калибровки и квазизамкнутого управления показана на рис. 2. Исходные данные для реализации алгоритма: Т — период нарезки зубцовой структуры; 1т — номинальная амплитуда тока в обмотке; Ф= — поток постоянного магнита; Ьп — номинальное значение индуктивности; Рт — максимальное значение статического синхронизирующего тягового усилия. Первоначальная форма тока — синусоидальная.
Предложенный алгоритм основан на измерении индуктивности в рабочем зазоре ЛТТТД, по результатам которой осуществляется расчет магнитного сопротивления, электромагнитных сил в полюсных наконечниках и остаточного тягового усилия, по которому производится необходимая коррекция тока в обмотках управления.
Измерение индуктивности ЛШД
Индуктивность шагового двигателя существенно влияет на его динамические характеристики. Вследствие присутствия насыщаемых магнитных материалов магнитный поток не пропорционален току, как это было бы в линейной магнитной системе. Когда магнитопровод входит в режим насыщения, изменение магнитного потока, вызванное изменением управляющего тока, становится меньше. В индуктивной системе любое изменение тока вызывает изменение магнитного поля, однако этот процесс инерционный. Так как магнитное поле имеет инерционный эффект на процесс задания управляющего тока, формирование тока происходит с задержкой. Таким образом, изменение тока и, следовательно, изменение магнитного потока, является зависимым во времени от задающего напряжения.
В достаточно общем случае, который учитывает насыщение, но пренебрегает потерями в магнитной цепи, ЛШД описывается электрическими уравнениями:
и, = ч, , (10)
где , = 1..т, т — число фаз;
магнитным уравнением:
т
^ = Е1 л, (11)
к=1
Задание исходных данных Задание синусоидальной формы тока Измерение индуктивности при малом токе Задание равномерной координатной сетки
-Да-
Инициализация?
Расчет дифференциального магнитного
сопротивления магнитопровода Расчет идеального дифференциального магнитного сопротивления зубцовой структуры Коррекция дифференциального магнитного сопротивления зубцовой структуры при малом токе
-Нет-
Коррекция координатной сетки
Коррекция дифференциального магнитного сопротивления зубцовой структуры при номинальном токе Расчет магнитного сопротивления зубцовой структуры Расчет остаточного динамического тягового усилия ШД Расчет погрешности позиционирования ШД
Да
Расчет дифференциального магнитного сопротивления магнитопровода
Измерение индуктивности при
номинальном токе
Рис. 2. Блок-схема алгоритма коррекции погрешности позиционирования электромеханическим уравнением:
I =
т г1 т т г ЛТ"
1=10 дх 1=11=10 дх
(12)
1=1 о
механическим уравнением:
й2х йх т—- + в— + ¥, = Р. йг2 Ж
Продифференцировав уравнения (11) по времени, получим:
йг
=2
к=1
( йЬ
йг
Ьзк
Ак I-йг
=2
к=1
Ь1к + ч
дЬ,к } А
дг.
к
йг
к. + г дЬ,к йх
дх йг
(13)
(14)
Таким образом, необходима характеристика магнитного потокосцепления, чтобы корректно моделировать двигатель. Это потокосцепление является функцией (т +1) переменных, что делает ее исследование достаточно затруднительным. Однако на практике эта функция может быть упрощена до зависимости только от двух переменных, так как взаимным влиянием фаз двигателя можно пренебречь.
Уравнения (11) и (12) оказывают, что характеристика магнитного потокосцепления, а следовательно, и тягового усилия, определяется функцией индуктивности. Соответственно очень важно иметь, с одной стороны, достаточно точные измерения и, с другой стороны, модель, позволяющую быстрый расчет согласно уравнению (14).
В электромеханической системе с насыщением не существует прямого метода измерения средней индуктивности Ь . Однако измерение дифференциальной индуктивности возможно. Индуктивность рассматривается как функция только двух переменных — управляющего тока и положения индуктора относительно статора.
На рис. 3 показана зависимость потокосцепления Ш от функции тока г. Средняя индуктивность определяется как следующий коэффициент:
ш
Ь(о ) = —. (15)
го
Дифференциальная индуктивность определяется углом наклона касательной к функции
Ш = / (г): д Ш
I = дф-. (16)
д г
Из соотношений (15) и (16) получим:
(17)
1 (го ) = Ь(г0 ) + го
дг
Таким образом, можно сказать, что если дифференциальная индуктивность известна, то потокосцепление может быть рассчитано как интеграл:
го
Ш(г0)=|Ш. (18)
0
Средняя индуктивность Ь получается делением потокосцепления на го:
1 го
ь(г0 ) = — | ыг. (19)
го о
Для нулевого тока го = 0, Ь(о) = I (о).
В системе с переменным магнитным сопротивлением соотношения (15)-(19) определены для каждого положения индуктора.
Рассмотрим электромеханическую систему (ЛШД) с насыщением, где подвижная часть (индуктор) фиксирована в определенном положении. Прикладывая скачек напряжения амплитуды и1 в одной обмотке двигателя, описываемого уравнением (10), где отдельные обмотки разомкнуты, ток в этой фазе задается соотношением (рис. 4):
и1 = Я, +
г ^ л лг
Т ■ дЬ,, Ьи + г1 д1
\ д г1 J
~Т. (20)
Л
о
1=10
Рис. 3. Зависимость потокосцепления как функции тока
Рис. 4. Переходный процесс в системе с насыщением из-за единичного изменения напряжения
Дифференциальная самоиндуктивность, соответствующая току г10, получается как
1 ( ) = Т + г д)Ьп = и1 - г 10
дг1 йг^йг
(21)
- и,
Рис. 5. Напряжение и ток во время измерения дифференциальной индуктивности В другой фазе, где ток равен нулю, получим напряжение и2 как
и й (т ) и2 =-ЦТ = л'=
(
Ь12 + г
. дЬ,
Л
дг
1 У
йг1
йг
Дифференциальная взаимоиндуктивность будет равна:
112 ((10 ) Ь12 + г10
дЬ1
и2
дг1 йг11йг
Таким образом, если измерить и 1, Я1 и функции
% = I (г) и и 2 = I (1),
(22)
(23)
то можно
получить дифференциальные само- и взаимоиндукции по выражениям (21) и (23).
Возможность реализации системы измерения индуктивности
Принципиальная трудность данной измерительной системы — это дифференцирующий элемент, который должен дать напряжение, пропорциональное производной по времени. Очевидно, что дифференцирование электрического сигнала обычно сопровождается рядом проблем. Среди них принципиальными являются чувствительность к возмущениям, возможная нестабильность и время реакции. Однако с помощью правильного масштабирования сигналов спроектировать подходящее устройство возможно.
Для того чтобы система была действительно качественной, она должна быть специально рассчитана для каждого конкретного применения из-за масштабирования сигналов и подбора постоянных времени. Время реакции такого дифференцирующего элемента ненулевое, и для предотвращения случайного изменения di/dt в момент приложения напряжения в фазу подается скачек напряжения от — U\ до + U\ (рис. 5).
Переходный процесс заканчивается за время между моментом коммутации tc прохождения значения тока через нуль, и измерение выполняется для i = 0..Imax .
В качестве примера реализации системы измерения индуктивности фаз ЛТТТД рассмотрим контроллер SoftStep [3], созданный нами совместно с центром лазерной техники LST Laser & Strahl Technik (г. Вена, Австрия). Функциональная блок-схема измерительной части контроллера показана на рис. 6 для двухфазного ЛШД.
Цифроаналоговые преобразователи (TP1 и TP4) формируют управляющие синусоидальные токи фаз ЛШД, а также наложенные на них пробные прямоугольные сигналы для измерения индуктивностей. Частота пробных прямоугольных сигналов должна быть значительно выше рабочей частоты управления ЛШД. После интегрирующего звена C/R1 (TP2 и TP5) пробные прямоугольные импульсы преобразуются в сигналы треугольной формы, в то время как управляющие напряжения ЛШД остаются неизменными. Следующий за интегратором каскад служит для токового управления фазовых обмоток ЛШД (L1 и L2) пропорционально напряжению на входе интегратора.
Измерение напряжения на обмотке осуществляется с помощью дифференциального усилителя, причем следующий аналогово-цифровой преобразователь с выборкой и хранением позволяет измерять напряжение как с дифференциального усилителя (TP3 и TP6), так и с выхода интегратора. Ключи S1 и S2 служат также для компенсации погрешностей изготовления электрических компонентов интегратора путем проведения цикла калибровки измерительной электрической цепи контроллера.
Описанная измерительная схема позволяет проводить измерение дифференциальной индуктивности, не прерывая нормального управления ЛШД.
Расчет управляющих токов, прямоугольных импульсов и обработка измерительных данных производится процессором цифровой обработки сигналов ADSP-2100, который хранит также и скорректированные токи.
Заключение
1. Установлено, что погрешность позиционирования ЛШД — это периодическая функция координаты с амплитудой от 5 до 50 мкм и частотой, пропорциональной скорости движения.
2. Погрешность позиционирования ЛШД может быть существенно уменьшена как с помощью алгоритмов управления, построенных на использовании результатов внешних измерений, так и на основе программной коррекции по результатам измерения индуктивности в рабочем зазоре. При прочих равных условиях программная коррекция является предпочтительной и экономически обоснованной, поскольку не требует дорогостоящей аппаратуры внешних измерений.
3. Предложенные математическая модель и алгоритм программной коррекции основаны на измерениях индуктивности ЛШД в процессе движения и калибровки токов управления по результатам расчета остаточного динамического тягового усилия.
Компьютерное моделирование системы управления ЛТТТД на основе разработанного алгоритма калибровки и квазизамкнутого управления и проведенные нами экспериментальные исследования позволили разработать и создать контроллер управления многокоординатным приводом с программной коррекцией погрешности позиционирования.
движения положение
Рис. 6 Функциональная блок-схема системы измерений контроллера 8ой81ер
PROGRAM COMPENSATION OF LINEAR STEPPING DRIVE POSITIONING ERROR
S.E. KARPOVICH, YU.S. MEZHINSKY, V.V. ZHARSKY
Abstract
External sensors are not necessary for increasing LSM positioning accuracy when open-loop control is used. The control system can be built on the base of internal sensor and digital processing of inductance measurement in the air gap "inductor-stator". Equations of considered system are obtained. The positioning error correction algorithm that allows increasing positioning accuracy in 5-10 times without external sensors is realized.
Литература
1. Карпович С.Е., Русецкий А.М., Ляшук Ю.Ф. Теория построения прецизионных механизмов оборудования производства электронной техники. Мн., 1999.
2. Кенио Т. Шаговые двигатели и их микропроцессорные системы управления. М., 1987.
3. Карпович С. Е., Жарский В. В., ЛяшукЮ. Ф., Межинский Ю. М. Прецизионные координатные системы на основе электропривода прямого действия. Мн., 2001.
4. Русецкий А.М. Координатные позиционеры гибких производственных систем для электронного машиностроения. Мн., 1998.
5. Исии Т., Симояна И., ИноуэХ. и др. Мехатроника. М., 1988.
6. Онегин Е.Е. Точное машиностроение для микроэлектроники. М., 1986.
