Например, Бобцов

АЛГОРИТМ УПРАВЛЕНИЯ ПОДАТЛИВЫМ ДВИЖЕНИЕМ МАНИПУЛЯТОРА С УПРУГИМИ ЭЛЕМЕНТАМИ В ШАРНИРНЫХ СОЧЛЕНЕНИЯХ

Алгоритм управления податливым движением манипулятора

29
УДК 007.52

В. В. ТИТОВ
АЛГОРИТМ УПРАВЛЕНИЯ ПОДАТЛИВЫМ ДВИЖЕНИЕМ МАНИПУЛЯТОРА С УПРУГИМИ ЭЛЕМЕНТАМИ В ШАРНИРНЫХ СОЧЛЕНЕНИЯХ
Рассматриваются подходы к управлению податливым движением манипуляторов с упругими элементами в шарнирных сочленениях. Приведен краткий анализ подходов и представлен вариант улучшения динамических показателей существующих регуляторов.
Ключевые слова: манипулятор, робот, силомоментное управление, шарнирные сочленения.
Введение. Одним из основных требований к манипуляционным операциям в робототехнике является успешное осуществление контактных взаимодействий с внешней средой. Необходимость выполнения подобных операций дала начало направлению в робототехнике, связанному с очувствлением по моментам и усилиям взаимодействия (или силомоментному управлению). Это направление активно развивалось последние 30—40 лет. В результате в робототехнике сформировались несколько подходов к выполнению силомоментных операций (см., например, [1]).
Для осуществления силомоментного управления требуется использовать датчики сил и моментов, расположенные в шарнирах манипулятора либо в запястье схвата (имеются и другие точки зрения: см., например, [2]). В основе измерения усилий, возникающих в датчиках момента, лежит измерение величины деформаций механических элементов конструкции датчика [3].
Использование датчиков момента в шарнирах манипуляторов, а также волновых редукторов обусловливает существенную упругость/податливость механических передач сочленений робота. Упругость механических передач приводит к значительному количеству негативных эффектов, сказывающихся на динамике манипулятора и усложнении процесса управления им. Наиболее податливыми элементами конструкции манипулятора являются шарнирные сочленения, кистевой датчик сил и моментов, а также звенья манипулятора.
Основными целями при управлении манипуляторами с упругими элементами в шарнирных сочленениях являются:
— минимизация дополнительных динамических эффектов, вызванных наличием упругого элемента между управляемым выходным валом сочленения и управляющим валом двигателя, развивающего момент;
— увеличение механической податливости манипулятора при сохранении точности отработки манипуляционных операций.
Отметим, что в области управления манипуляторами с упругими элементами в шарнирных сочленениях за последние 30 лет были достигнуты значительные успехи, некоторые из которых упомянуты в обзоре [4].
Модель объекта управления. Примем, для упрощения дальнейших выкладок, что все сочленения манипулятора являются вращательными парами (степенями подвижности). Кроме того, будем считать, что датчики моментов расположены в сочленениях манипулятора.
Уравнения динамики движения таких манипуляторов, построенные исходя из уравнения Лагранжа второго рода (см., например, [5]), получаются сильносвязными и трудно анализируемыми. В этой связи в большинстве работ применяется упрощенная модель Спонга, получаемая введением ряда допущений и имеющая следующий вид [6]:

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 11

30 В. В. Титов

B(q)q + C(q, q)q + g(q) =
Jmα + T + DK −1T = Tm ; T = −Kθ = K (α − q),

T

+

DK

−1T

+

Text

;⎫ ⎪ ⎬ ⎪ ⎭

(1)

где α — вектор углов поворота выходных валов двигателя до достижения валом упругого

элемента,

α=

d2α dt 2

;

q

— вектор

углов поворота

выходных

валов,

расположенных

после

уп-

ругого

элемента,

q

=

dq dt

,

q

=

d 2q dt 2

,

qi

= θi + αi ;

θ



вектор

угловых

деформаций

упругих

элементов; B(q) — n×n-матрица инерции манипулятора; C(q, q) — моменты кориолисовых

и центробежных сил; g(q) — моменты силы тяжести; T — момент деформации упругих эле-

ментов,

T

=

dT dt

;

Tm —

развиваемый

двигателем

момент;

Text



момент

внешних

сил;

D



диагональная матрица коэффициентов внутреннего вязкого трения в упругих элементах; K —

диагональная матрица коэффициентов жесткости упругих элементов; Jm — диагональная матрица инерции роторов двигателей (приведенных к выходным валам).

Запись уравнений (1) относительно переменных состояния (q, q,T ,T ) , являющихся, как

правило, доступными для измерения, имеет следующий вид:

B(q)q + C(q, q)q + g(q) = T Jm K −1T + DK −1T + T = Tm

+ DK −1T + Text , − Jm B−1(q)(T + DK −1T

+ Text



C(q, q)q



⎫⎪ g (q)).⎬⎪⎭

(2)

Алгоритмы управления. В задачу системы управления силомоментно-очувствленными роботами входит как регулирование положения характерной точки манипулятора (например, схвата), так и регулирование усилия, развиваемого в этой точке.
Следуя методу, изложенному в работе [7], построим систему управления манипулято-
ром с помощью точной линеаризации. Выражение для Tm , позволяющее точно линеаризовать
системы уравнений (1), (2), имеет вид

Tm

=

J m−1 D −1 ( B(q)u

+

B(q)q

+ n(q, q)) +

Jm−1q + T

+

(DK −1



J

−1 m

D−1

)T

,

(3)

где n(q, q) = C(q, q)q + g(q) , и — сигнал управления.
При замыкании системы (1) с использованием выражения (3) для момента привода выполняется уравнение q = u .
Точная линеаризации не является самым простым способом управления, так как требует значительных вычислительных затрат, а именно: вычисления линеаризующей обратной связи (3) и вычисления двух производных по времени от q , а также первой производной от момента Т.
В работе [8] управление роботом с упругими элементами в шарнирных сочленениях исследовалось с позиций теории сингулярных возмущений, согласно которой первые два уравнения системы (1) могут рассматриваться независимо, так как переходные процессы, описываемые вторым выражением, происходят намного быстрее, чем в первом. Поэтому предлагается использовать раздельные контуры управления: контур управления моментами упругих элементов (быстрый контур, с локальной линеаризацией) и контур управления положением по типу импедансного (медленный контур). Основным недостатком предложенного подхода является необходимость использования в быстром контуре управления данных, поставляемых медленным контуром управления. Для линеаризации в быстром контуре также требуется вы-

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 11

Алгоритм управления податливым движением манипулятора

31

числение или изменение величин Text либо q , что вызывает дополнительные трудности и не

всегда представляется возможным.

В более поздних исследованиях были предложены варианты управления манипулято-

ром на основе свойства пассивности (passivity based control) динамических систем. Так, в ра-

боте [9] используется обратная связь по α и α , за счет чего обеспечивается свойство пассив-

ности самого контроллера. В этой работе рассматривается модель, аналогичная модели Спон-

га (1), а в качестве контура управления моментом предложен регулятор, описываемый выра-

жением

Tm

=

Jm Jm−1u + T

+

DK −1T



J

m

J

−1 m

(T

+ Ds K −1T ) ,

(4)

где Jm — диагональная матрица требуемых моментов инерции роторов двигателей (приве-
денных к выходным валам), Ds — положительно-определенная диагональная матрица демпфирования.
Регулятор положения выходных валов, используемый совместно с регулятором (4), характеризуется выражением [9]

u = −Kα (α − αd ) − Dαα + g(qd ), αd = qd + K −1g(qd ),

(5)

где αd и qd — векторы желаемых углов поворота выходных валов до и после достижения
упругого элемента; Kα , Dα — положительно-определенные матрицы жесткости и демпфирования.
Уравнения динамики робота после замыкания имеют вид

B(q)q Jmα +

+ T

C(q, q)q + + Ds K −1T

g(q) = = −Kα

−Kα (α −

(α αd

− )

α −

d ) − Dαα + g(qd Dαα + g(qd ).

)



J



+

Text

,⎫⎪ ⎬ ⎪⎭

(6)

Регулятор (5) реализует аналог импедансного управления. Однако следует отметить, что
установившаяся ошибка отработки конечного положения выходных валов Kα (α − αd ) = Text
при q = q = 0 не гарантирует (и не может гарантировать) высокую точность отработки при
движении. Основным недостатком регулятора (5) является значительная погрешность отработки положения при выполнении динамичных манипуляционных операций с сохранением требуемых импедансных свойств. Интересным представляется закон управления выходным моментом (4), позволяющий виртуально уменьшить значение момента инерции ротора двигателя.
Наиболее предпочтительным способом построения системы управления манипулятором является, как и в рассмотренных выше случаях, разделение процессов управления моментами упругих элементов и управления положением валов. Данное разделение целесообразно проводить также и на физическом уровне, принимая во внимание пропускную способность каналов связи. Кроме того, желательна как можно более точная отработка траектории при сохранении свойств податливости манипулятора.
Воспользуемся моделью динамики манипулятора, описываемого системой уравнений (2), совместно с регулятором момента (4). Тогда выражения для динамики замкнутой системы примут следующий вид:

B(q)q + C(q, q)q + g(q) = u + (D − Jm K −1T + T + Ds K −1T = u − Jmq.

Ds

)K

−1T



Jm

q



Jm

K

−1T

+

Text

,⎫⎪ ⎬ ⎪⎭

(7)

Рассмотрим, далее, подход к линеаризации движения манипулятора при отсутствии уп-

ругости конструкции.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 11

32 В. В. Титов

Следуя методике, изложенной в работе [1], введем закон управления

u = [B(q)qd + KD (qd − q) + KP (qd − q)] + C(q, q)q + g(q) ,

(8)

где qd , qd , qd — желаемые программные ускорения, скорости и положения выходных валов;

KD , KP — положительно-определенные матрицы, определяющие динамику изменения

погрешности отработки в отсутствие внешнего момента (Text = 0) . При Text = 0 движение

манипулятора можно описать уравнением

B(q)(qd − q) + KD (qd − q) + KP (qd − q) = 0 ,

а при наличии внешнего момента Text — уравнением

B(q)(qd − q) + KD (qd − q) + KP (qd − q) = −Text . Заметим, что полностью исключить взаимозависимость обобщенных координат мани-

пулятора таким образом не удается, так как остается составляющая B(q) . Однако удается до-

биться взаимной независимости статических погрешностей отработки положений выходных

валов:

KP (qd − q) = −Text .

Результаты моделирования. Компьютерное моделирование предложенного алгоритма

(см. формулы (4), (8)) проводилось на модели двухстепенного пространственного робота-

манипулятора. В качестве задающего воздействия на каждый из шарниров подавался косину-

соидальный сигнал по положению единичной амплитуды (1 рад). В качестве контура управ-

ления моментом использовался регулятор (4) с одинаковой настройкой для всех эксперимен-

тов; настройки параметров сравниваемых регуляторов (5) и (8) были идентичны ( Kα = KP ,

Dα = Kd ). На рис. 1 и 2 соответственно приведены результаты моделирования работы регу-

лятора (5) и регулятора (8) с неполным задающим сигналом (т.е. при qd = 0 , qd = 0 ,

qd = qd (t) ). Сравнительный анализ результатов показал, что регулятор (8) проявил себя лучше, имея практически только фазовое отставание от задающего сигнала.

q, рад 1

0,5

0

–0,5 –1
–1,5 0
q, рад

0,5

1 1,5

положение 1-го шарнира положение 2-го шарнира

задающий сигнал

2 2,5

3 3,5 4

4,5

Рис. 1

t, c

1

0,5

0

–0,5 положение 1-го шарнира –1 положение 2-го шарнира

–1,5 0

0,5

1 1,5

задающий сигнал

2 2,5

3 3,5 4

4,5 t, c

Рис. 2

Применение полного задающего сигнала ( qd = qd (t) , qd = qd (t) , qd = qd (t) ) намного

улучшает результат, как видно из рис. 3.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 11

Алгоритм управления податливым движением манипулятора

33

q, рад

1

0,5

0

–0,5 положение 1-го шарнира –1 положение 2-го шарнира

–1,5 0

0,5

1 1,5

задающий сигнал

2 2,5

3 3,5 4

4,5 t, c

Рис. 3

Предлагаемый регулятор (8) имеет, тем не менее, ряд недостатков: значительную зави-

симость качества управления от точности задания массоинерционных параметров модели ро-

бота, используемой в регуляторе; затрудненный анализ устойчивости вследствие большей

связности по сравнению с регулятором (5).

СПИСОК ЛИТЕРАТУРЫ

1. Robotics: Modelling, Planning and Control / B. Siciliano, L. Sciavicco, L. Villani, G. Oriolo. London: Springer, 2009.

2. Кулаков Ф. М. Активное силомоментное управление роботами без использования запястных силомоментных
сенсоров // Изв. РАН. Сер. Теория и системы управления. 2012. № 1. C. 149—170.

3. Гориневский Д. М., Формальский А. М., Шнейдер А. Ю. Управление манипуляционными системами на основании информации об усилиях / Под ред. В. С. Гурфинкеля и Е. А. Девятина. М.: Физматлит, 1994.

4. Ozgoli S., Taghirad H. D. A survey on the control of flexible joint robots // Asian J. of Control. 2006. Vol. 8, N 4.
P. 1—15.

5. Козлов В. В., Макарычев В. П., Тимофеев А. В., Юревич Е. И. Динамика управления роботами. М.: Наука, 1984.

6. Spong M. Modeling and control of elastic joint robots // IEEE J. of Robotics and Automation. 1987. P. 291—300.

7. Дыда А. А. Адаптивное и нейросетевое управление сложными динамическими объектами. Владивосток: Дальнаука, 2006.

8. Ott C., Albu-Schaffer A., Kugi A., Hirzinger G. Decoupling based cartesian impedance control of flexible joint
robots // Proc. IEEE Intern. Conf. on Robotics and Automation. 2003. Vol. 3. P. 3101—3107.

9. Albu-Schaffer A., Ott C., Hirzinger G. A unified passivity-based control framework for position, torque and
impedance control of flexible // The Intern. J. of Robotics Research. 2007. Vol. 26. P. 23—39.

Виктор Викторович Титов

Сведения об авторе — аспирант; СПИИРАН, лаборатория информационных технологий в управ-
лении и робототехнике; E-mail: victortitov2005@mail.ru

Рекомендована СПИИРАН

Поступила в редакцию 10.06.12 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 11