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

АЛГОРИТМ УПРАВЛЕНИЯ АВТОНОМНЫМ ДВУХКОЛЕСНЫМ МОБИЛЬНЫМ РОБОТОМ «МОТОБОТ»

А.А. Бобцов, А.С. Боргуль, К.А. Зименко, А.А. Пыркин

Литература
1. Lauwers T.B., Kantor G.A., Hollis R.L. A dynamically stable single-wheeled mobile robot with inverse mouse-ball drive // IEEE International Conference on Robotics and Automation. – 2006. – 2884 р.
2. Nagarajan U., Mampetta A., Kantor G., Hollis R. State transition, balancing, station keeping, and yaw control for a dynamically stable single spherical wheel mobile robot // IEEE International Conference on Robotics and Automation. – 2009. – P. 3161–3166.
3. Колюбин С.А., Пыркин А.А. Управление нетривиальными маятниковыми системами в условиях параметрической и функциональной неопределенностей // Научно-технический вестник СПбГУ ИТМО. – 2010. – № 69. – С. 34–39.
4. Anderson B.D.O., Moore J.B. Optimal control: linear quadratic methods. – Prentice-Hall, 1989. – 394 p. 5. Ha Y.-S., Yuta S. Trajectory tracking control for navigation of self-contained mobile inverse pendulum //
Proc. IEEE/RSJ Int’l. Conf. on Intelligent Robots and Systems. – 1994. – Р. 1875–1882. 6. Linear-Quadratic-Regulator (LQR) design – MATLAB [Электронный ресурс]. – Режим доступа:
http://www.mathworks.com/help/toolbox/control/ref/lqr.html, св. Яз. англ. (дата обращения 18.03.2011). 7. Быстров С.В., Григорьев В.В., Рабыш Е.Ю., Черевко Н.А. Экспоненциальная устойчивость непре-
рывных динамических систем // Научно-технический вестник СПбГУ ИТМО. – 2011. – № 73. – С. 44–47.

Боргуль Александр Сергеевич

– Санкт-Петербургский государственный университет информационных

технологий, механики и оптики, студент, borgulalexandr@gmail.com

Громов Владислав Сергеевич

– Санкт-Петербургский государственный университет информационных

технологий, механики и оптики, студент, object253@yandex.ru

Зименко Константин Александ- – Санкт-Петербургский государственный университет информационных

рович

технологий, механики и оптики, студент, kostyazimenko@gmail.com

Маклашевич Сергей Юрьевич

Санкт-Петербургский государственный университет информационных

технологий, механики и оптики, студент, s.maklashevich@gmail.com

УДК 681.513.1, 681.513.3
АЛГОРИТМ УПРАВЛЕНИЯ АВТОНОМНЫМ ДВУХКОЛЕСНЫМ МОБИЛЬНЫМ РОБОТОМ «МОТОБОТ»
А.А. Бобцов, А.С. Боргуль, К.А. Зименко, А.А. Пыркин
Рассмотрена задача управления мобильным роботом типа мотоцикл, построенном на базе Lego Mindstorms NXT. Для макета робота разработан регулятор, обеспечивающий устойчивое вертикальное положение при перемещении по горизонтальной плоскости с возможностью дистанционного управления траекторией движения. Ключевые слова: мобильные роботы, автономный мотоцикл, неустойчивые системы, управление по выходу.
Введение
Роботы сегодня входят в нашу жизнь в различных областях. Они летают в космос, исследуют другие планеты, помогают в военных целях – разминируют бомбы и разведывают обстановку с воздуха. В промышленности многие отрасли уже немыслимы без роботов: они собирают автомобили, помогают находить новые лекарства. Многие устройства, принимающие решения на основе полученных от сенсоров данных, тоже можно считать роботами, например, лифты, стиральные машины, системы антиблокировочного торможения, помогающие избежать аварий.
Одну из наиболее важных проблем, которую можно решить с помощью робототехники – это вопрос, связанный с транспортом. В условиях современной загруженности дорог (пробки в городе и, наоборот, бездорожье в сельской местности) особую ценность представляет одноколейный транспорт, об-

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 5 (75)

63

АЛГОРИТМ УПРАВЛЕНИЯ АВТОНОМНЫМ ДВУХКОЛЕСНЫМ МОБИЛЬНЫМ …
ладающий маневренностью, малыми габаритами и пониженным энергопотреблением. Так, в частности, с помощью мотоцикла можно проехать в пробке, найти тропинку на проселочной дороге и решить многие другие задачи, связанные с современным темпом жизни. С развитием робототехники, теории и систем автоматического управления стало возможным решение задачи автономного движения мотоцикла.
Задача стабилизации одноколейного двухколесного транспорта в вертикальном положении была поставлена еще в середине прошлого века и относится к классическим задачам механики и теории управления. Известны работы [1, 2], в которых уделено внимание задаче управления «роботомвелосипедистом», включая математические модели и возможные алгоритмы управления. Достижение устойчивости на одноколейном транспорте является гораздо более сложной проблемой, чем на любой другой платформе, поэтому никто не преуспел в создании автономных мотоциклов, несмотря на все их преимущества.
С развитием высокоточного приборостроения и робототехники поменялся подход к созданию аппаратов данного типа. Самой известной является разработка фирмы DARPA GrandChallenge [3].
Рис. 1. Автономно-управляемый мотоцикл DARPA GrandChallenge
В этой разработке стабилизация осуществляется за счет скорости поворота руля. Также модель оснащена GPS навигацией для определения карты местности. Система технического зрения по трем каналам вместе с GPS позволяют создать подробную карту местности. Преимущества двухколесной системы перед трех-, четырех-, шестиколесными конструкциями – в большей управляемости и уменьшении жесткости привязки к земной поверхности. Следует отметить, что для одноместного мотоцикла требуется постоянная настройка, которая чрезвычайно трудно достигается (в частности, в данной разработке более 140 тысяч строк программного кода для успешной стабилизации).
Устойчивость мотоцикла зависит от баланса сил, которыми управляет водитель. При длительном или напряженном темпе движения проявляются такие факторы, как усталость, нервное напряжение – причина большинства дорожно-транспортных происшествий. Другой проблемой при удержании равновесия мотоцикла является зависимость гироскопического стабилизирующего момента от линейной скорости движения. При низких скоростях исчезает управляемость в системе.
В работе описывается макет автономного мобильного робота «Мотобот», являющегося моделью двухколесного мотоцикла. Как и в [4], для построения макета используется конструктор Lego Mindstorms NXT. Анализируется математическая модель движения двухколесного одноколейного робота на плоскости, на основе которой синтезируются алгоритмы автоматического управления. Задача может быть решена путем управления рулевым колесом и скоростью ведущего колеса. Решены локальные задачи стабилизации вертикального положения равновесия, движения вдоль прямой с заданной скоростью, а также реализовано дистанционное управление траекторией движения с помощью джойстика и радиоканала связи.
Математическая модель робота «Мотобот» Математическая модель мотоцикла является нелинейной и включает несколько важных параметров, таких как масса мотоцикла, коэффициент трения, расстояние между проекциями колес, моменты инерции составных частей относительно центра масс и др. Математическая модель состоит из четырех основных частей: рамы, двух колес и рулевой вилки. Приемлемая модель имеет 7 степеней свободы, из них 3 – степени свободы угловых скоростей. Модель может быть соответственно параметризованной с помощью следующих переменных: угол крена вдоль поперечной оси ; угол поворота с помощью рулевого колеса  , изменение уровня вращения заднего колеса  . За основу была взята модель сбалансированного одноколейного транспорта с симметрией масс (рис. 2, а) [5]. Рассмотрим декартову систему координат с центром в точке O и осями x , y и z (рис. 2, б). Заднее колесо R, соединенное с рамой, при наличии водителя имеет центр масс в точке B.
64 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 5 (75)

А.А. Бобцов, А.С. Боргуль, К.А. Зименко, А.А. Пыркин

Центр масс B размещен в следующей области: xB  0 , yB  0 , zB  0 . Момент инерции рамы относительно центра масс представлен матрицей моментов инерции размерности (3 3) , где все массы симмет-

ричны относительно проекции x z . Центр масс передней части рамы с вилкой и рулем H находится в сле-

дующей области xH  0 , yH  0 , zH  0 . Матрицы моментов инерции передней и задней сборок равны

 IBxx 0 IBxz 

IHxx 0 IHxz 

IB



 

0

I Byy

0

 

,

IH



 

0

I Hyy

0

 

.

 IBxz 0 IBzz 

 IHxz 0 IHzz 

аб

Рис. 2. Прототип мобильного робота «Мотобот» (а) и схематическое изображение углов положения составных частей (б)

Представим линеаризованные уравнения движения:
P  Fi  i  T  T  TR  R ,
где T , T , TR – обобщенные силы каждой линейной комбинации компонентов различных приложен-
ных сил F . Принятие полной линеаризации и движения без завала, по прямой, может быть достигнуто при
любой скорости, удовлетворяющей условию    RrR . Боковая симметрия системы, совмещенная с
линейностью уравнений, позволяет получить зависимость между прямолинейным движением и углами падения и поворота. Исходя из этого, линеаризованное уравнение движения для первой степени свободы (движения вперед) просто получить, исследуя систему по осям x, z :

[rR2mT  IRyy  (rR / rF )2 IFyy ]R  TR .
Линеаризованные уравнения для двух других степеней свободы, угла падения и поворота рулевого колеса – это пара совмещенных дифференциальных уравнений второго порядка. Уравнение в канонической форме, где первое уравнение отражает крен, а второе поворот руля, имеет вид M q  Cq  K q  f .

Для модели мотоцикла запишем его в форме M q   C1 q  [g K 0  2 K 2 ]q  f , где изменяемые во вре-

мени переменные q  [, ]T и f  [T ,T ]T ; M – симметричная матрица масс, которая отображает кине-

тическую энергию мотоцикла. Матрица затухания C   C1 линейной скорости  и охватывает косо-

симметрические гироскопические моменты, возникающие из-за уровней крена мотоцикла и рулевого

управления. C1 также содержит инерциальные реакции, возникающие во время рулевого управления. Матрица жесткости K состоит из двух частей: независимой от скорости симметричной части g K 0 , про-

порциональной гравитационному ускорению, которая может быть использована при расчете изменения

потенциальной энергии, и части 2 K2 , которая пропорциональна квадрату скорости и зависит от гиро-

скопического и центробежного эффектов. Таким образом, имеем характеристическое уравнение:

det(M 2   C1   g K 0  v2 K 2 )  0 .

(1)

Далее следует рассматривать решения [4] при различных скоростях. На низких и высоких скоро-

стях асимптотическая устойчивость не достижима, в силу ряда факторов. На начальных этапах рассмат-

ривается диапазон средних скоростей, где решается задача асимптотической устойчивости.

В общем виде модель движения выглядит следующим образом:

q  S(q,t)u , u  [M(q,t)]1Q(q,u,t) ,

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 5 (75)

65

АЛГОРИТМ УПРАВЛЕНИЯ АВТОНОМНЫМ ДВУХКОЛЕСНЫМ МОБИЛЬНЫМ …

где q – обобщенная координата; u – обобщенная скорость; S – кинематическая матрица, которая ото-
бражает отношение уровней обобщенных координат к обобщенным скоростям; M – система матриц масс; Q содержит все силы и зависимости от скорости инерции.

Нахождение параметров линеаризованной модели движения робота

Будем использовать следующие индексы: R – для заднего колеса; B – для задней части рамы совместно с мотоциклистом; H – для передней части с рулем; F – для переднего колеса; T – для системы в целом; A – для сборки передней части рамы с рулем.
Найдем общую массу и расположение центра масс (относительно точки контакта заднего колеса с землей P):
mT  mR  mB  mH  mF , xT  (xBmB  xH mH  wmF ) / mT ,
zT  (rRmR  zBmB  zH mH  rF mF ) / mT . Для вывода момента инерции всего мотоцикла вдоль оси z используем соответствующие моменты масс и значения инерций по отношению к контактной точке P вдоль общих осей, с учетом зависимости осесимметричных моментов инерции переднего и заднего колес:
ITzz  IRzz  IBzz  IHzz  IFzz  mB xB2  mH xH2  mF w2 .
Таким же способом описывается передняя сборка A. Рассмотрим вектор   sin , 0, cos T , на-

правленный вниз вдоль рулевой оси, где  – угол на плоскости x z . Центр масс передней сборки распо-

ложен выше руля, перпендикулярно u A  (xA  w  c) cos   z A sin  . Для передней сборки требуются три
специальные величины инерции: момент инерции относительно оси руля и значения моментов инерции относительно смещенных и косых осей, взятых около точки, где они пересекаются. Из этих соотношений берется момент около одной оси, возникающий из-за углового ускорения около других.
Отношение механического трейла (расстояние между продолжением рулевой вилки и проекцией
передней оси на поверхность) к базе колес   (c / w) cos  .

Угловые моменты заднего и переднего колес вдоль оси y, вычисленные с помощью линейной скорости, вместе с их суммой из гиростатических коэффициентов, равны
SR  I Ryy / rR , SF  I Fyy / rF , ST  SR  SF .
Сформируем параметры линеаризованных уравнений движения.

M



M 

 

M



M  M 

  

,

K0



 K0

 

K

0

K0 K0

  

,

K2



 K2

 

K

2

K2 

K 2 

 

,

C1



C1 C1

C1 C1

  

,

где M   ITxx , M   I Ax  ITxz , M   M  , M   I A  2I Az  2 ITzz – элементы симметричной матрицы масс M ; K0  mT zT , K0  SA , K0  K0 , K0  SA sin  – гравитационные коэффициенты жесткости (умноженные на g ); K2  0 , K2  0 , K2  ((ST  mT zT ) / w) cos  , K2  ((SA  SF sin ) / w) cos  – скоростные коэффициенты жесткости (умноженные на v2 ); C1  0 , C1  ST  SF cos   (ITxz / w)cos   mT zT , C1  (IAz / w)cos  (SA  (ITxz / w)cos) , C1  (ST  SF cos) – коэффициенты матрицы «затухания» C .

Расчет параметров регулятора
Запишем уравнения баланса угловых моментов. Баланс углового момента крена относительно продольной оси:
mT yP zT  ITxx  ITxz  I Ax   ST  SF cos   TB  gmT zT   gSA .
Баланс углового момента рыскания относительно точки P: mT yP xT  ITxz  ITzz  I Az   ST  SF sin   wFFy .
Баланс углового момента руля для передней сборки: mA yPuA  I Az  I A  SF ( cos    sin )  TH  cFFy cos   g(   sin )SA .
Используя выражение, характеризующее изменение угла крена относительно изменения угла рыскания рамы, получим итоговое выражение баланса мотоцикла:
yP  ((2  c) / w) cos  . Уравнения задания положения мотоцикла относительно задней контактной точки: xP   cos , yP   sin  .

66 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 5 (75)

А.А. Бобцов, А.С. Боргуль, К.А. Зименко, А.А. Пыркин

Для модели вход-состояние-выход, где x   , , , T – вектор состояния, найдем закон управле-
ния u  T . Расчет регулятора для случая адаптивной системы будем проводить с помощью алгоритма «полоска». Для неадаптивной системы воспользуемся линейно-квадратичным регулятором с критерием

качества J  (xT Qu x  uT Ruu)dt , где Qu и Ru – диагональные матрицы штрафов по состоянию и 0
управлению. Выберем Qu  diag{0,0,1,0} , так как итоговая цель балансирования мотоцикла, удержание вертикального положения, может быть достигнута минимизацией угла крена. Для контроля усилий возьмем Ru  1 .
Алгоритм управления рулевым колесом по измерению угла крена выберем в виде пропорционально-интегрально-дифференциального (ПИД) регулятора, коэффициенты которого подберем методом Циглера–Никольса: k p  8,23 – для пропорционального, ki  1,48 – для интегрального и kd  2,76 – для
дифференциального составляющих сигнала управления.

Заключение

На базе макета мобильного робота «Мотобот» были проведены экспериментальные исследования системы управления. На рис. 3 приведены временные диаграммы реальных процессов в системе во время эксперимента.


50 0 –50
–100 0

2000

4000 а

6000 t, с

20 0 –20 –40
0

2000

4000 б

6000 t, с

Рис. 3. Временные диаграммы процессов в роботе «Мотобот»:
угол крена мотоцикла  [град/c] (а); угол поворота рулевого колеса  [град] (б)

Как видно из рис. 3 синтезированная система обладает достаточным запасом устойчивости не только для постоянного поддержания мотоцикла в вертикальном положении, но и для возможности управления направлением и скоростью движения.
В отличие от ближайшего аналога [5], в этой работе были сделаны ключевые изменения в конструкции, позволившие обеспечить стабильное передвижение робота по горизонтальной плоскости в течение длительного периода времени. Была уменьшена длина рамы и трейла, изменена конструкция руля. Расположение гироскопа ближе к центру сделало конструкцию реалистичнее. В отличие от [5], вместо пропорционального регулятора был использован ПИД-регулятор, что позволило компенсировать статическую ошибку слежения и демпфировать колебания. Было реализовано дистанционное управление с джойстика для управления траекторией движения «Мотобота». На рис. 4 представлены кадры, полученные во время проведения эксперимента.

Рис. 4. Движение мобильного робота «Мотобот» с коррекцией траектории
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 5 (75)

67

АЛГОРИТМ УПРАВЛЕНИЯ АВТОНОМНЫМ ДВУХКОЛЕСНЫМ МОБИЛЬНЫМ …

Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009–2013 годы (государственный контракт № П498 от 08.08.2010 г.; проекты 2011-1.2.1-402014/057 и 2011-1.2.2-215-009/023).

Литература
1. Фомин В.Н., Фрадков А.Л., Якубович В.А. Адаптивное управление динамическими объектами. – М.: Наука. Гл. ред. физ.-мат. лит., 1981. – 448 с.
2. Селиванов А.А., Никитин Д.А. Адаптивное управление роботом-велосипедистом // XII конференция молодых ученых «Навигация и управление движением». – СПб, 2010. – 16 с.
3. Архив агентства по перспективным оборонным научно-исследовательским разработкам США [Электронный ресурс]. – Режим доступа: http://archive.darpa.mil/, свободный. Яз. англ.
4. Бобцов А.А., Капитанюк Ю.А., Капитонов А.А., Колюбин С.А., Пыркин А.А., Чепинский С.А., Шаветов С.В. Технология LEGO MINDSTORMS NXT в обучении студентов основам адаптивного управления // Научно-технический вестник СПбГУ ИТМО. – 2011. – № 1. – С. 103–108.
5. Meijaard J.P., Papadopoulos J.M., Ruina A., Schwab A.L. Linearized dynamics equations for the balance and steer of a bicycle: a benchmark and review // Proceedings of the Royal society A: mathematical, physical and engineering sciences. – 2007. – V. 463. – № 2084. – Р. 1955–1982.

Бобцов Алексей Алексеевич

– Санкт-Петербургский государственный университет информационных

технологий, механики и оптики, доктор технических наук, профессор,

декан, bobtsov@mail.ifmo.ru

Боргуль Александр Сергеевич

– Санкт-Петербургский государственный университет информационных

технологий, механики и оптики, студент, borgulalexandr@gmail.com

Зименко Константин Александ- – Санкт-Петербургский государственный университет информационных

рович

технологий, механики и оптики, студент, kostyazimenko@gmail.com

Пыркин Антон Александрович

– Санкт-Петербургский государственный университет информационных

технологий, механики и оптики, кандидат технических наук, ассистент,

a.pyrkin@gmail.com

68 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 5 (75)