Параметрический синтез закона управления движением мобильного робота в городских условиях Текст научной статьи по специальности «Механика и машиностроение»

УДК 681.51:007.52:631.324

Све Лин Хту Аунг

ПАРАМЕТРИЧЕСКИЙ СИНТЕЗ ЗАКОНА УПРАВЛЕНИЯ ДВИЖЕНИЕМ МОБИЛЬНОГО РОБОТА В ГОРОДСКИХ УСЛОВИЯХ

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

E-mail: swelinhtooaung@gmail.com

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

Разработка методов управления автономным движением мобильных роботов (МР), в последнее время привлекает внимание все большего числа исследователей. Управление естественным образом строится на основе информации, поступающей с бортовых сенсорных систем, а также внешних источников информации (GPS, Глонасс). В качестве бортовых систем наиболее часто используются системы зрения, выполненные на базе СТЗ и сканирующих лазерных дальномеров [1,2]. Хорошим доказательством актуальности проблемы развития методов автономного управления является проведение соревнований МР [3], в которых участвуют ведущие университеты развитых стран. Так, в России в течение ряда лет проходит Девянинский фестиваль, в рамках которого проводятся соревнования автономных МР, выполняющих ряд заданий на специально оборудованном полигоне. В 2010 г. состоялись первые в России соревнования МР, движущихся по трассе, проложенной по пересеченной местности (берег озера Селигер) [4]. Аналогичные соревнования проводились в США в 2004 и 2005 гг. [3], а в 2007 г. там же проходили соревнования роботов в городских условиях. В процессе этих соревнований роботы, сконструированные на базе серийно выпускаемых автомобилей, должны были достичь предварительно указанной на карте точки, соблюдая в процессе движения соответствующие дорожные правила [3, 4].

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

Рис. 1. Движение робота в городских условиях

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

Постановка задачи. Пусть колесный МР движется по плоской горизонтальной поверхности в городских условиях. При этом целью управления является обеспечение движения МР вдоль тротуара с использованием имеющейся на борту системы очувствления (рис. 1), в качестве которой выступает сканирующий лазерный 3Э дальномер. При помощи лазерного дальномера составляется цифровое описание рельефа рабочей зоны (картография), определяется положение автомобиля на дороге (локализация), осуществляется слежение за соседней полосой, контролируется дистанция от тротуара. В работе [5] для решения 3Э 8ЬЛМ-проблемы был предложен метод, позволяющий постоянно отслеживать координаты положения робота и интерпретировать информацию об окружающей среде по измерениям дальномера.

Рассматриваемая в настоящей работе задача состоит в определении и выборе параметров, входящих в закон управления, обеспечивающих устойчивость и приемлемое качество траектории.

Предположим, что абсолютная система координат расположена так, что ось Х0 параллельна тротуару, а сам тротуар представляет собой (в плане) последовательность отрезков, лежащих на одной прямой, расположенной на оси Х0. Таким образом, не рассматриваются более сложные движения робота, например повороты.

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

Кинематическая и динамическая модели колесного МР. Рассмотрим результаты, полученные в работах [6-8]. На рис. 2 приведена

Рис. 2. Схема движения робота вдоль непрерывного тротуара

схема движения колесного МР с соответствующими системами координат. Учитывая динамические свойства робота, а также наличие исполнительной системы [8], получим уравнения движения МР. Считаем, что исполнительная система МР представляет собой два электромеханических привода, выполненных на базе двигателей постоянного тока, каждый из которых приводит в движение ведущие колеса правого и левого бортов шасси транспортного средства МР. Тогда уравнения, связывающие компоненты фазового вектора с управляющими напряжениями, можно записать следующим образом:

x = и cos ф; У = и sin ф,

¥ = и;

1

и = -Ъи2 +-(rL + tr );

pm

• Ътии W . ч

и = -:- + (Tr - Tl);

2piz

R

TL =

TR =

M

2L

RM

TL

2L

Tr

¿2 км кт Г ky

il км кт г ky

Lp

и

и +

wWN

wWs

ip kM kyTT

+--;— Ul ;

+

L

ip kM ky

Ur ,

(1)

где Ь — координата центра масс платформы, лежащего на оси симметрии; т — масса робота; W — расстояние между центрами ведущих колес; гг — момент инерции платформы относительно вертикальной оси ZR; р — радиус колеса; Ям — сопротивление обмотки якоря исполнительного двигателя; гР — передаточное отношение редуктора; км — коэффициент, связывающий ток и развиваемый двигателем момент; ктг — коэффициент передачи тахогенератора; ку — коэффициент усилителя мощности; Ь — индуктивность обмотки якоря; UR,Ць — напряжения, подаваемые на обмотку якоря.

Соотношения (1) представляют собой систему нелинейных дифференциальных уравнений седьмого порядка с фазовым вектором

z

(х, у, ф, и, ш, тя, гь)т и вектором управления (ия, Цъ)т. При этом угловые скорости правого (шя) и левого колес, угловая ш и линейная и скорости робота связаны следующими соотношениями (линейной скоростью робота мы называем скорость начала связанной системы координат):

1 / ш^

шя = - и +

Р V 2

1 f uW

шь = - и -

р V 2

Примем следующий способ формирования управляющих напряжений, подаваемых на приводы ведущих колес МР:

ur = Ud +

(2)

ul = Ud - Uv,

где Ud = const — напряжение, обеспечивающее заданную линейную скорость движения МР; U^ — напряжение, управляющее угловой скоростью движения МР.

Запишем следующий закон управления:

U^ = -kip + k2(-y + y*), (3)

где k1,k2 — некоторые параметры; p — угол между продольной осью робота и направлением вдоль тротуара; у* и y — желаемая и текущая координаты тротуара в собственной системе координат. Параметры у и p определяются по сканам, получаемым от 3D лазерного дальномера.

Сделаем теперь следующие предположения. Считаем, что при движении МР на плоскости вдоль гладкой программной траектории управляющие моменты ведущих колес, развиваемые электродвигателями, изменяются незначительно (это допущение часто используют при исследовании электромеханических систем), т.е. tr = tl = 0. Это допущение позволяет понизить до пятого порядок системы.

С учетом сказанного и соотношений (2) и (3) дифференциальные уравнения для моментов правого (tr) и левого (tl) ведущих колес МР системы (1) можно записать следующим образом:

• rm ivkM ктгку f uW\

TL = -"2MTL--LP" +

+^Mky (Ud + kip + k2y - k2y*) = 0;

• rm ipkM Kky

TR = -~2L TR - [V + +

+ipkMky (Ud - kip - k2y + k2y*) = 0.

Введем обозначения

mктгky

Rm p

и

^ _ 2ip км ky (5)

RM

учитывая которые, получаем выражения

TR + TL _ -2^V + 25Ud;

tr - tl _ -^uW + 2S(-kip - k2y + k2y*). Тогда уравнения управляемого объекта примут следующий вид:

x _ и cos p;

У _ и sin p;

Р _ u; (6)

и _ -bu2 - Av + BUd; u _ Cuu - Du - Eklp - Ek2y + Ek2y*,

где

A = B = C = pm pm iz

„ W2v ^ WS

D =--, E =

2ргг ргг

Соотношения (6) представляют собой систему нелинейных дифференциальных уравнений пятого порядка, описывающих движение МР на плоскости под действием напряжений Ць, подаваемых на якорные обмотки двигателей правого и левого бортов и формируемых в соответствии с соотношениями (2) и (3).

Задача состоит в том, чтобы определить и выбрать параметры к1,к2, обеспечивающие приемлемое качество траектории движения МР вдоль тротуара.

Выбор параметров к1, к2, входящих в закон управления. Полученные уравнения движения МР представляют собой систему нелинейных дифференциальных уравнений пятого порядка, что значительно осложняет задачу дальнейшего анализа и выбора необходимых коэффициентов к1, к2 для реализации алгоритма управляющего напряжения Ц^. Поэтому для выбора коэффициента к1,к2 проведем линеаризацию уравнений движения так, чтобы можно было использовать классические методы анализа линейных стационарных систем.

В качестве программной траектории выберем х (£) = ис£, у (£) = у*,

* * *

ф (£) = 0, и (£) = ис, ш (¿) = 0, т.е. равномерное и прямолинейное движения робота параллельно тротуару.

Тогда уравнения в отклонениях относительно программной траектории будут иметь вид

Д х = Ди;

Д у = исДф;

Д ф = Дш; (7)

Д и = -АДи;

Д ш = -Ек2Ду - ЕкгДф - БДш.

Видно, что систему (7) можно разбить на две независимые подсистемы: одна относительно вектора (Дх, Ди)т, другая — относительно вектора (Ду, Дф, Дш)т. Отметим вторую подсистему, в которую входят параметры кг,к2.

Соответствующий характеристический многочлен /(Л) имеет следующий вид:

/ (Л) = Л3 + БЛ2 + ЕкгЛ + Ек2ис. (8)

Проанализируем зависимость расположения корней от коэффициентов кг,к2, воспользовавшись критерием Гурвица, в соответствии с которым необходимым и достаточным условиями расположения корней характеристического уравнения в левой полуплоскости являются положительные значения главных миноров матрицы Гурвица:

Б Ек2ис 0

Н = 1 Екг 0 0 Б Ек2ис

Дг = Б > 0;

Д2 = БЕкг - Ек2ис > 0;

Дз = Ек^сД) > 0.

Полученные соотношения приводят к неравенству

0 < ^ < к\. (9)

Отметим, что для типового мобильного робота [1] со значениями параметров т = 30 кг, {г = 4,2 кг^м2, р = 56,7, £ = 4,4, Ш = 500 мм, р = 87мм, Ь = 0 имеем А = 43,4483; В = 3,3716; С = 0; Б = 19,3966; Е = 6,02.

Исследуем теперь случай, когда переходный процесс имеет апериодический характер, т.е. корни характеристического уравнения (8) не только лежат в левой полуполоскости, но еще и действительны.

Запишем соотношение коэффициентов кубического уравнения (8), при выполнении которого все три корня действительны [9]:

где

p3 q2

— + — < 0 27 4

(10)

Ь2 с 3а2 а

- d ^ 27а3 3а2 а'

а = 1, Ь = Д с = Е&1, d = Ек2ис.

Тогда одновременное выполнение условий (9), (10) позволяет выделить область значений параметров к1 ,к2, обеспечивающих апериодический переходный процесс.

На рис.3 показана плоскость параметров к1 ,к2. В области 1 система (в приращениях) неустойчива, в области 2 — переходный процесс имеет колебательный характер, а в области 3 — апериодический. Безусловно, полученный результат имеет приближенной характер, поскольку мы рассматривали систему в малых отклонениях от программной траектории.

Определение параметров у, входящих в закон управления.

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

Рис. 3. Области значений параметров к±, к2, при которых траектория движения МР неустойчива (1), имеет колебательный (2) и апериодический (3) характер

изображения — сечение рельефа в плоскости дороги, тротуара и парапета. Необходимо определить угол ^ и расстояние y между роботом и тротуаром по данным лазерного дальномера.

Угол и расстояние между роботом и тротуаром вычисляются в результате определения линии пересечения плоскостей дороги, тротуара и парапета каждым сканом дальномера относительно координат дальномера в процессе движения робота. Приведем уравнения линий пересечения дороги Lld, тротуара Llt и парапета Llh, связанные с функцией сечения дальномера:

где в — угол отклонения плоскости скана от горизонтальной плоскости; к — расстояние от центра координат робота до дороги; к — высота тротуара; ^ — угол между продольной осью робота и направлением вдоль тротуара; у — получаемое расстояние между роботом и тротуаром. На рис. 4 показано движение МР вдоль тротуара.

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

На рис. 5 показаны схемы движения робота вдоль непрерывного тротуара. Угол наклона отрезка АВ, а также координаты его концевых точек позволяют легко вычислить параметры у, входящие в

Lid = sinßXi - h = 0;

Lit = sin ßXi - h + ht = 0;

(11)

Lih = sin ^ cos ßX/ + cos pYl + y = 0,

z

h

W

Рис. 4. Схема движения МР вдоль тротуара

хг

h sinß h hT sinß А \lh V/t

В \

xv ¿/Д

h sinß h-hT sinß А Llh Lh

В

У\

а в

У\ У\

У\

-У=У б

Рис. 5. Движение робота вдоль разрывного тротуара

выражение для закона управления (3):

1

^ =

cos ß

tg м; y =

cos фъYA - kaYB )

ka - къ

(12)

где ka = h ctg ß, къ = (h - h) ctg ß.

Рассмотрим алгоритм управления роботом в случае движения вдоль разрывного тротуара (рис. 6). В каждый момент времени от 3D дальномера на борту получаем набор сканов, имеющих разные углы в вертикальной плоскости и, следовательно, в разных точках пересекающих тротуар. В том случае, когда робот должен миновать разрыв в тротуаре (например, перекресток), необходимо выбрать ближайший скан, которому соответствуют приемлемые значения параметров ^ и У из (12).

Экспериментальные исследования. Для реализации закона управления угловой скоростью движения автономного МР в форме (3) необходимо выбрать значения параметров управления ki, k2, для чего выполним компьютерное моделирование движения МР на динамической модели, описываемой системой линейных дифференциальных уравнений (7). Проверку полученного результата проведем путем мо-

Рис. 6. Движение МР под некоторым углом и параллельно вдоль разрывного тротуара в системе координатах дальномера х\, у1,

делирования движения МР на нелинейной динамической модели, описываемой системой дифференциальных уравнений (6).

Были выбраны следующие начальные условия: x(0) = 0; y (0) = 1,5; ^(0) = 0,345 рад; и(0) = 0; и(0) = 0; Ul(0) = Ur(о) = 0, т.е. робот неподвижен, находится на расстоянии 1,5 м от тротуара с ориентацией 20о. В качестве программной траектории запишем

x (t) = uct, y (t) = 1, ^ (t) = 0, и (t) = uc, и (t) = 0.

Рассмотрим движение системы при ki = 15, k2 = 2. Результаты сравнения параметров движения при использовании нелинейной и линеаризованной моделей приведены на рис. 7, из которых следует,

Рис. 7. Изменение угловых параметров движения (а) и и (£) (б), а также координат ж(£) и у (£) (в) при использовании линеаризованной (сплошная кривая) и нелинейной (штриховая) моделей МР

что при выбранных значениях параметров k1,k2, лежащих в области 2 (см. рис. 3), переходные процессы имеют апериодический характер и для нелинейной (штриховая кривая) и линеаризованной (сплошная кривая) моделей отличаются незначительно.

Выводы. Предложен способ управления колесным МР при его движении в городских условиях вдоль прямолинейного тротуара. В качестве сенсорной системы использован 3D лазерный сканирующий дальномер, поставляющий в систему управления робота информацию о дальности до объектов внешней среды, находящихся в передней полусфере. Рассмотрен метод выбора значений параметров закона управления, обеспечивающих апериодический характер переходного процесса при движении робота вдоль тротуара. Учтены динамические свойства объекта управления, а также параметры приводов. Выделены три области значений параметров, при которых траектория движения вдоль тротуара неустойчива и имеет апериодический или колебательный характер. Полученные результаты позволяют найти первые приближения значений параметров закона управления. Приведены результаты вычислительного эксперимента.

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

1. Зенкевич С. Л., Космачев П. В. Управление движением мобильного робота в неподвижную точку // Мехатроника, автоматизация, управление. -2010. -№ 3. - С. 55-60.

2. Н о с к о в В. П., Р у б ц о в И. В. Опыт решения задачи автономного управления движением мобильных роботов // Мехатроника, автоматизация, управление. -2005. -№ 12. - С. 16-24.

3. Официальная страница соревнований «Grand Challenge» [электронный ресурс] http://www.darpa.mil/GRANDCHALLENGE/

4. Официальная страница соревнований «Робокросс Селигер-2010» [электронный ресурс] http://www.robosport.ru

5. Зенкевич С. Л., Минин А. А., Све Лин Хту Аунг. Еще одно решение задачи построения 3D-карты // Мехатроника, автоматизация, управление. -2010. -№ 6. - С. 55-61.

6. Д е в я н и н Е. А. О движении колесных роботов // Труды конф. "Мобильные роботы и мехатронные системы". - М., 1998. - С. 169-200.

7. Кобрин А. И., Марты ненко Ю. Г. Неголономная динамика мобильных роботов и ее моделирование в реальном времени // Труды конф. "Мобильные роботы и мехатронные системы". - М., 1998. - С. 107-123.

8. З е н к е в и ч С. Л., Н а з а р о в а А. В. Система управления мобильного робота // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. - 2006. - № 3. -С. 31-51.

9. Формула Кардано [электронный ресурс] http://ru.wikipedia.org/wiki/Фор-мула_Кардано

Статья поступила в редакцию 24.09.2010 Све Лин Хту Аунг родился в 1980 г. (Союз Мьянма). Аспирант кафедры "Робото-технические системы" МГТУ им. Н.Э. Баумана. Степень магистра защитил в МГТУ им. Н.Э. Баумана в 2007 г. Автор трех научных работ в области систем автоматизированного управления робототехническими системами.

Swe Lin Htoo Aung (b. 1980), Union of Myanmar. Master's degree was defended in the Bauman Moscow State Technical University in 2007. Post-graduate of "Robotic Systems" department of the Bauman Moscow State Technical University. Author of 3 publications in the field of systems of automated control of robotic systems.