Алгоритм управления автономным двухколесным мобильным роботом «Мотобот»

Из полученных графиков видно, что использование качественно-экспоненциального регулятора приводит к уменьшению колебательности процесса стабилизации болбота и более точному позиционированию конструкции.

Заключение

В результате проделанной работы был получен опытный образец и синтезирован алгоритм, позволяющий стабилизировать в вертикальном положении объект на сферическом катке, обладающий достаточным запасом качества для отработки дополнительных управляющих воздействий. Был произведен синтез линейно-квадратичного регулятора и регулятора на основе качественной экспоненциальной устойчивости. Анализ полученных данных показал, что регулятор на основе качественной экспоненциальной устойчивости обеспечивает лучшие точностные и качественные показатели, а, следовательно, является более предпочтительным в условиях поставленной задачи.

Литература

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. Для макета робота разработан регулятор, обеспечивающий устойчивое вертикальное положение при перемещении по горизонтальной плоскости с возможностью дистанционного управления траекторией движения. Ключевые слова: мобильные роботы, автономный мотоцикл, неустойчивые системы, управление по выходу.

Введение

Роботы сегодня входят в нашу жизнь в различных областях. Они летают в космос, исследуют другие планеты, помогают в военных целях - разминируют бомбы и разведывают обстановку с воздуха. В промышленности многие отрасли уже немыслимы без роботов: они собирают автомобили, помогают находить новые лекарства. Многие устройства, принимающие решения на основе полученных от сенсоров данных, тоже можно считать роботами, например, лифты, стиральные машины, системы антиблокировочного торможения, помогающие избежать аварий.

Одну из наиболее важных проблем, которую можно решить с помощью робототехники - это вопрос, связанный с транспортом. В условиях современной загруженности дорог (пробки в городе и, наоборот, бездорожье в сельской местности) особую ценность представляет одноколейный транспорт, об-

ладающий маневренностью, малыми габаритами и пониженным энергопотреблением. Так, в частности, с помощью мотоцикла можно проехать в пробке, найти тропинку на проселочной дороге и решить многие другие задачи, связанные с современным темпом жизни. С развитием робототехники, теории и систем автоматического управления стало возможным решение задачи автономного движения мотоцикла.

Задача стабилизации одноколейного двухколесного транспорта в вертикальном положении была поставлена еще в середине прошлого века и относится к классическим задачам механики и теории управления. Известны работы [1, 2], в которых уделено внимание задаче управления «роботом-велосипедистом», включая математические модели и возможные алгоритмы управления. Достижение устойчивости на одноколейном транспорте является гораздо более сложной проблемой, чем на любой другой платформе, поэтому никто не преуспел в создании автономных мотоциклов, несмотря на все их преимущества.

С развитием высокоточного приборостроения и робототехники поменялся подход к созданию аппаратов данного типа. Самой известной является разработка фирмы БЛЯРЛ вгаМСИаИе^е [3].

Рис. 1. Автономно-управляемый мотоцикл DARPA GrandChallenge

В этой разработке стабилизация осуществляется за счет скорости поворота руля. Также модель оснащена GPS навигацией для определения карты местности. Система технического зрения по трем каналам вместе с GPS позволяют создать подробную карту местности. Преимущества двухколесной системы перед трех-, четырех-, шестиколесными конструкциями - в большей управляемости и уменьшении жесткости привязки к земной поверхности. Следует отметить, что для одноместного мотоцикла требуется постоянная настройка, которая чрезвычайно трудно достигается (в частности, в данной разработке более 140 тысяч строк программного кода для успешной стабилизации).

Устойчивость мотоцикла зависит от баланса сил, которыми управляет водитель. При длительном или напряженном темпе движения проявляются такие факторы, как усталость, нервное напряжение -причина большинства дорожно-транспортных происшествий. Другой проблемой при удержании равновесия мотоцикла является зависимость гироскопического стабилизирующего момента от линейной скорости движения. При низких скоростях исчезает управляемость в системе.

В работе описывается макет автономного мобильного робота «Мотобот», являющегося моделью двухколесного мотоцикла. Как и в [4], для построения макета используется конструктор Lego Mindstorms NXT. Анализируется математическая модель движения двухколесного одноколейного робота на плоскости, на основе которой синтезируются алгоритмы автоматического управления. Задача может быть решена путем управления рулевым колесом и скоростью ведущего колеса. Решены локальные задачи стабилизации вертикального положения равновесия, движения вдоль прямой с заданной скоростью, а также реализовано дистанционное управление траекторией движения с помощью джойстика и радиоканала связи.

Математическая модель робота «Мотобот»

Математическая модель мотоцикла является нелинейной и включает несколько важных параметров, таких как масса мотоцикла, коэффициент трения, расстояние между проекциями колес, моменты инерции составных частей относительно центра масс и др. Математическая модель состоит из четырех основных частей: рамы, двух колес и рулевой вилки. Приемлемая модель имеет 7 степеней свободы, из них 3 - степени свободы угловых скоростей. Модель может быть соответственно параметризованной с помощью следующих переменных: угол крена вдоль поперечной оси ф; угол поворота с помощью рулевого колеса 5 , изменение уровня вращения заднего колеса 6 .

За основу была взята модель сбалансированного одноколейного транспорта с симметрией масс (рис. 2, а) [5]. Рассмотрим декартову систему координат с центром в точке O и осями х, y и z (рис. 2, б). Заднее колесо R, соединенное с рамой, при наличии водителя имеет центр масс в точке B.

Центр масс В размещен в следующей области: хв = 0, ув = 0, гв < 0. Момент инерции рамы относительно центра масс представлен матрицей моментов инерции размерности (3 х 3) , где все массы симметричны относительно проекции х2 . Центр масс передней части рамы с вилкой и рулем Н находится в следующей области Хн = 0 , ун = 0, 2н < 0. Матрицы моментов инерции передней и задней сборок равны

1 Вхх 0 ^ Вх 1Нхх 0

IВ = 0 1Вуу 0 , I н = 0 1Нуу 0

_ ^Вхг 0 1' Вгг _ } Нх2 0 1Нгг

а б

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

составных частей

Представим линеаризованные уравнения движения: Р = X Р  Аv,. = ГфАФ + Т5 А8 + Гвв Д0,, где Тф , Т§ , - обобщенные силы каждой линейной комбинации компонентов различных приложенных сил Р.

Принятие полной линеаризации и движения без завала, по прямой, может быть достигнуто при любой скорости, удовлетворяющей условию V = -(3. Боковая симметрия системы, совмещенная с линейностью уравнений, позволяет получить зависимость между прямолинейным движением и углами падения и поворота. Исходя из этого, линеаризованное уравнение движения для первой степени свободы (движения вперед) просто получить, исследуя систему по осям х, 2 :

2 2 " [гктт + 1Куу + (ГК /Гр) 1Руу ]0Д = Т0К .

Линеаризованные уравнения для двух других степеней свободы, угла падения и поворота рулевого колеса - это пара совмещенных дифференциальных уравнений второго порядка. Уравнение в канонической форме, где первое уравнение отражает крен, а второе поворот руля, имеет вид Мс( + Сс| + Kq = f .

Для модели мотоцикла запишем его в форме Мq + vС q + [gК0 +v К2] q = f, где изменяемые во времени переменные q = [ф, 8]т и f = [Тф, Тъ]т ; М - симметричная матрица масс, которая отображает кинетическую энергию мотоцикла. Матрица затухания С = v С! линейной скорости V и охватывает косо-симметрические гироскопические моменты, возникающие из-за уровней крена мотоцикла и рулевого управления. С1 также содержит инерциальные реакции, возникающие во время рулевого управления. Матрица жесткости К состоит из двух частей: независимой от скорости симметричной части g К0 , пропорциональной гравитационному ускорению, которая может быть использована при расчете изменения потенциальной энергии, и части V2 К2 , которая пропорциональна квадрату скорости и зависит от гироскопического и центробежного эффектов. Таким образом, имеем характеристическое уравнение:

<1е1(М X2 + V Сх Х + g К 0 + V2 К 2) = 0.

Далее следует рассматривать решения [4] при различных скоростях. На низких и высоких скоростях асимптотическая устойчивость не достижима, в силу ряда факторов. На начальных этапах рассматривается диапазон средних скоростей, где решается задача асимптотической устойчивости. В общем виде модель движения выглядит следующим образом:

q = 8^,г)и , и = [М^,и,0,

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

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

Будем использовать следующие индексы: R - для заднего колеса; B - для задней части рамы совместно с мотоциклистом; H - для передней части с рулем; F - для переднего колеса; T - для системы в целом; A - для сборки передней части рамы с рулем.

Найдем общую массу и расположение центра масс (относительно точки контакта заднего колеса с землей P):

mT = mR + mB + mH + mF , xt = (xBmB + xHmH + wmF)/mT ,

zT = (~rRmR + zBmB + zHmH -rFmF)/mT .

Для вывода момента инерции всего мотоцикла вдоль оси z используем соответствующие моменты масс и значения инерций по отношению к контактной точке P вдоль общих осей, с учетом зависимости осесимметричных моментов инерции переднего и заднего колес:

2 2 2

ITzz = IRzz + IBzz + IHzz + IFzz + mBxB + mH xH + mFw

Таким же способом описывается передняя сборка A. Рассмотрим вектор X = (sin X, 0, cos Х), направленный вниз вдоль рулевой оси, где X - угол на плоскости xz . Центр масс передней сборки расположен выше руля, перпендикулярно ua = (xa - w - с) cos X - za sin X. Для передней сборки требуются три специальные величины инерции: момент инерции относительно оси руля и значения моментов инерции относительно смещенных и косых осей, взятых около точки, где они пересекаются. Из этих соотношений берется момент около одной оси, возникающий из-за углового ускорения около других.

Отношение механического трейла (расстояние между продолжением рулевой вилки и проекцией передней оси на поверхность) к базе колес | = (с / w) cos X .

Угловые моменты заднего и переднего колес вдоль оси y, вычисленные с помощью линейной скорости, вместе с их суммой из гиростатических коэффициентов, равны

SR = IRyy / rR , SF = 1 Fyy / TF , ST = SR + SF -

Сформируем параметры линеаризованных уравнений движения.

M =

~M фф M ф," , K 0 = К л0фф К 0ф5 , K 2 = К 2фф К2ф5 , C1 = ^фф с ' 1ф5

_ M5ф M55 _ К 05ф К055 ' 2 _ К25ф К 255 ' 1 с ^ 15ф с 155

где Мфф = ITxx , Мф5 = IAXx

M 5ф= Mф,

M55 = IAXX + 2|IAXz + ц ITz - элементы симметричной мат-

рицы масс M ; К0фф = mTzT , К0ф5 = -SA , К05ф = К0ф5, К055 = -SA sin X - гравитационные коэффициенты жесткости (умноженные на g); К2фф = 0, К25ф = 0, К2ф5 = ((ST - mTzT) / w) cosX ,

K255 = ((SA + SF sin X)/ w)cos X - скоростные коэффициенты жесткости (умноженные на v2); С1фф = 0, С1ф5 = |ST + SF cos X+(ITxz / w) cos X - imTzT, C155 = (IAXz / w) cos X+|SA + (ITxz / w) cos X), С15ф = -(|ST + SF cos X) -коэффициенты матрицы «затухания» C .

Расчет параметров регулятора

Запишем уравнения баланса угловых моментов. Баланс углового момента крена относительно продольной оси:

-mTy pzt + ITxxф + IT,z Ф + iaXx5 + 4>vst +5VSf cos X = ^ - gmTzT ф+ gSA 5 .

Баланс углового момента рыскания относительно точки P:

mTypXT + ITxzф + ITzzФ + iaXz5 - Фvsт - 5vSf sin X = wFFy .

Баланс углового момента руля для передней сборки:

mAypUA +1axzФ + Iaxx 5 + vSf (-ф cos X + Ф sin X) = TH5 - cFfv cos X + g(ф + 5 sin X)Sa . Используя выражение, характеризующее изменение угла крена относительно изменения угла рыскания рамы, получим итоговое выражение баланса мотоцикла:

yP = ((v25 + vc5) / w) cos X .

Уравнения задания положения мотоцикла относительно задней контактной точки: Хр =v cos Ф , Ур = v sin Ф .

Для модели вход-состояние-выход, где х = ^ср, 5, ф, - вектор состояния, найдем закон управления и = Тд. Расчет регулятора для случая адаптивной системы будем проводить с помощью алгоритма «полоска». Для неадаптивной системы воспользуемся линейно-квадратичным регулятором с критерием

да

Г т т

качества 3 = I (х Qux + и ^ии)Ж, где Qu и Яи - диагональные матрицы штрафов по состоянию и

0

управлению. Выберем Qu = diag{0,0,1,0}, так как итоговая цель балансирования мотоцикла, удержание вертикального положения, может быть достигнута минимизацией угла крена. Для контроля усилий возьмем Яи = 1.

Алгоритм управления рулевым колесом по измерению угла крена выберем в виде пропорционально-интегрально-дифференциального (ПИД) регулятора, коэффициенты которого подберем методом Циг-лера-Никольса: кр = 8,23 - для пропорционального, к^ = 1,48 - для интегрального и к^ = 2,76 - для

дифференциального составляющих сигнала управления.

Заключение

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

Ф

50 0

-50

-100

5 20

0

-20

-40

2000

4000 а

6000

г, с

2000

4000 б

6000

г, с

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

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

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

Рис. 4. Движение мобильного робота «Мотобот» с коррекцией траектории

0

0

Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной

России» на 2009-2013 годы (государственный контракт № П498 от 08.08.20i0 г.; проекты 20ii-i.2.i-402-

014/057 и 20ii-i.2.2-2i5-009/023).

Литература

1. Фомин В.Н., Фрадков А.Л., Якубович В.А. Адаптивное управление динамическими объектами. - М.: Наука. Гл. ред. физ.-мат. лит., i98i. - 448 с.

2. Селиванов А.А., Никитин Д.А. Адаптивное управление роботом-велосипедистом // XII конференция молодых ученых «Навигация и управление движением». - СПб, 20i0. - i6 с.

3. Архив агентства по перспективным оборонным научно-исследовательским разработкам США [Электронный ресурс]. - Режим доступа: http://archive.darpa.mil/, свободный. Яз. англ.

4. Бобцов А.А., Капитанюк Ю.А., Капитонов А.А., Колюбин С.А., Пыркин А.А., Чепинский С.А., Шаветов С.В. Технология LEGO MINDSTORMS NXT в обучении студентов основам адаптивного управления // Научно-технический вестник СПбГУ ИТМО. - 20ii. - № i. - С. i03-i08.

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. - Р. i955-i982.

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

технологий, механики и оптики, доктор технических наук, профессор, декан, bobtsov@mail.ifmo.ru

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

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

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

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

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

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

ев