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



УДК 681.5

Формирование упрощенной траектории движения двуногого шагающего робота

А.К. Ковальчук, Д.Б.Кулаков, С.Е.Семенов

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

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

The article presents a method of forming a motion path of an actuator of a bipedal walking robot in the dynamic walking mode. The path is formed as periodic cycles of motion, described by piecewise harmonic functions with a small number of variable parameters.

Keywords: bipedal walking robot, motion path, harmonic functions, motion control synthesis.

Для проведения научно-исследовательских работ, направленных на изучение динамики двуногих шагающих роботов (ДШР), в МГТУ им. Н.Э. Баумана на кафедре «Гидромеханика, гидромашины и гидропневмоавтоматика» (Э10) совместно с Межотраслевым институтом повышения квалификации (при участии авторов) создан лабораторный исследовательский стенд ДШР [1]. Исполнительный механизм робота (рис. 1) состоит из корпуса и двух ног, имеет 12 степеней подвижности, оснащенных электрогидравлическими следящими приводами (ЭГСП). Высота исполнительного механизма ДШР составляет 2,2 м, вес 220 кг.

Применение ЭГСП обусловлено их преимуществами, по сравнению с другими типами приводов, по удельной мощности и динамическим характеристикам, что является особенно важным при создании ДШР больших размеров и масс. В корпусе ДШР размещены насосная станция, аппаратура системы управления, аппаратура системы ориентации. В стопы ног встроены силомоментные датчики. Снизу на стопах для увеличения податливости закреплены резиновые накладки. Энергопитание осуществляется по кабелю от внешнего источника трехфазного напряжения (200 В, 400 Гц). Программное обеспечение системы управления выполняется на трех совместно работающих ЭВМ, объединенных в локальную сеть Ethernet.

Задача синтеза динамически устойчивого движения ДШР чрезвычайно сложна. Насколько можно судить, к настоящему времени в общем виде она не решена. Однако управление динамической ходьбой

КОВАЛЬЧУК Александр Кондратьевич

кандидат технических наук, доцент, директор МИПК (МГТУ им. Н.Э. Баумана)

КУЛАКОВ Дмитрий Борисович

кандидат технических наук, доцент кафедры «Гидромеханика, гидромашины и гидропневмоавтоматика» (МГТУ им. Н.Э. Баумана)

СЕМЕНОВ Станислав Евгеньевич

кандидат технических наук, доцент, зав. кафедрой «Гидромеханика, гидромашины и гидропневмоавтоматика» (МГТУ им. Н.Э. Баумана)

Рис. 1. Лабораторный ДШР (вид спереди)

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

На рис. 2 приведена кинематическая схема исполнительного механизма ДШР, а также показано расположение абсолютной и связанных систем координат. Ось I абсолютной системы координат направлена вверх, робот движется вдоль оси х. Поскольку исполнительный механизм ДШР не закреплен на неподвижном основании, для описания положения исполнительного механизма ДШР в окружающем пространстве были введены фиктивные звенья с тремя поступательными и тремя вращательными сочленениями, соединяющие неподвижную абсолютную систему координат с корпусом робота.

Рис. 2. Кинематическая схема ДШР и назначение связанных систем координат

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

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

РсТ' 0,2'

0,15

0,1

0,05

-0,05

—1-

Рс1У/ |п Рстх \

-г/ V

// \ 1 к Рщ

54 54,5 55 55,5 56 56,6 57 57,5 58

58,5 59 /, с

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

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

Положение ДШР в декартовом пространстве в каждый момент времени определяется 11-мерным вектором положения S:

^ (х ст , У ст , £ ст , а ст , Р ст , У ст , Х цм , У цм , а к , Р к , У к ) ,

(1)

где IX

ст , У ст , £ с

)Т

координаты начала систе-

мы координат переносимой стопы в абсолютной системе координат; (аст, Рст, уст) — углы ориентации переносимой стопы в абсолютной системе координат; (хцм, уцм) — проекция центра масс ДШР на плоскость ОХУ абсолютной системы координат; (ак, Рк, ук ) — углы, определяющие ориентацию корпуса робота.

Движение ДШР в декартовом пространстве задается как программное изменение во времени координат вектора положения S (1). Процесс формирования траектории движения представляется в виде определения последовательности кадров положения. Каждый кадр положения определяется значением вектора положения S робота в декартовом пространстве в данный момент времени. Движение ДШР представляется как прохождение роботом положений в декартовом пространстве, определенных в последовательности кадров положений (рис. 5). Частота следования кадров задается в соответствии с требуемой точностью позиционирования робота в пространстве и динамическими свойствами исполнительного механизма.

Изменения значений компонент вектора положения (КВП) S (1) (с частотой дискретизации по времени, равной частоте следования кадров положений) определяются соответствующими функциями изменений переменных положения (ФИПП). Каждая ФИПП в общем случае может состоять из последовательности элементов ФИПП, представляющих собой простые функции с малым числом определяющих их параметров (рис. 6). Начало и продолжительность изменения /-й КВП, определяемого у-м элементом /-й ФИПП, определяется временем tb и Т у соответственно.

и ^

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

Рис. 5. Представление движения ДШР в виде последовательности кадров положений

Рис. 6. Временная диаграмма изменений компонент вектора положения X , Xст, Zст, определяемых элементами функций изменения переменных положения

векторных ФИПП S (У) (с постоянным составом элементов ФИПП). Такую векторную ФИПП S (У), определяющую траекторию движения робота в пределах шага, определим как модуль движения ДШР.

Модуль движения состоит из функций изменений следующих компонент вектора положения S (У):

— перенос центра масс исполнительного механизма в продольном направлении (хцм (У));

— перенос центра масс исполнительного механизма в поперечном направлении с переносимой стопы на опорную (уцм(У));

— подъем переносимой стопы (гст (У), а ст (У),

Рст(у), г ст(У));

— дополнительное угловое отклонение корпуса ДШР от первоначально заданного положения (аК(У), РК(У), уК(У));

— перенос переносимой стопы (хст (У),

у ст (У ),а ст (У ) Рст (У ) У ст (У));

— опускание переносимой стопы (гст (У),

аст(У), Рст (У), уст(У));

— возвращение углового положения корпуса ДШР к первоначально заданному значению (ак (У), РК(У), ук (У)).

Таким образом в нашем случае модуль движения состоит из 21 элемента ФИПП, определяющих изменения 11 КВП в пределах одного шага ДШР.

Из анализа траектории движения ДШР, полученной с помощью упрощенной аналитической модели (рис. 3 и рис. 4), видно, что изменения переменных положения ДШР S (1) в декартовом пространстве хорошо аппрокси-

мируются кусочными гармоническими функциями. В нашей работе все элементы ФИПП, аппроксимирующие изменения переменных положения ДШР, имеют вид

As(тs ) = As (1- cos(œsтs )).

ФИПП для каждой КВП в пределах модуля движения определяется выражением

N / \

Si (t) = (tнц ) + 2Js,(1 - cos(œssTsj )) x

j=1

x

' S i, |T S j

x

(t,- тsl) + \Tij - тS,

2\Tij - тs,

(2)

где s{ (t) — функция изменения i-й КВП (i-я ФИПП); Nt — количество элементов, состав-

Hj

ляющих i-ю ФИПП; As =-— амплитуда

sij 2

гармонического закона изменения i-й КВП на j-м элементе соответствующей ФИПП; H tj — требуемое изменение i-й КВП на j-м элементе

п

соответствующей ФИПП; ю? =— — частота

sij t ij

гармонического закона изменения i-й КВП на j-м элементе соответствующей ФИПП; Tj — продолжительность j-го элемента i-й ФИПП;

т^ = У — нц + Уь ) — относительное время у-го

элемента г'-й ФИПП: 0 < т^ <Тп; У — текущее время движения ДШР; Унц — время начала модуля движения; Уь — время начала у-го элемента 1-й ФИПП относительно времени начала модуля движения (рис. 6).

Формирование траектории ходьбы ДШР сводится к определению параметров Уь , Ту, Ну элементов ФИПП на каждом шаге робота (для каждого модуля движения). Эти параметры определяются при исследовании движения робота на имитационной (наиболее подробной) математической модели.

На рис. 7 представлены изменения КВП и параметры, характеризующие элементы ФИПП в пределах одного модуля движения ДШР. Условные обозначения, приведенные на рис. 7:

Xцм, ^цм — координаты центра масс робота;

х 12, х 18, у12, у18, г 12, г 18 — координаты начал систем координат стоп в абсолютной системе координат;

а — угол наклона корпуса во фронтальной плоск ости;

Уьцму и Тцму — момент начала поперечного движения центра масс робота и длительность движения центра масс робота;

Рис. 7. Параметры ФИПП, определяющих модуль движения ДШР

Т^ и Т^ — время начала и продолжительность фазы подъема левой стопы;

Т

bSdown1

и Т^^ — время начала и продолжительность фазы быстрого опускания правой стопы;

Т

bSdown2

и Т8Аот2 — время начала и продолжительность фазы замедленного опускания правой стопы;

tbакxfiom и tа* — время начала и продолЖИтельность фазы корректирующего наклона корпуса;

baкxreturn и Taкxretum — время начала и прОДолжительность фазы возврата корпуса к вертикальному положению;

Н^ — величина подъема переносимой стопы.

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

Коррекция траектории движения исполнительного механизма ДШР производится по результатам исследования подробной имитационной математической модели. Критерии качества могут быть самые разнообразные. Применительно к лабораторному ДШР, разработанному авторами, в качестве критерия качества синтезируемой траектории движения выбраны функции моментов на стопах робота, которые по абсолютной величине должны быть минимальными, а в идеале нулевыми на всем протяжении ходьбы. Таким способом обеспечивается «запас на стабилизацию» по моментам на стопах, необходимый для удер-

жания робота на выбранной траектории в процессе движения.

Коррекция движения ДШР с использованием инерционных свойств исполнительного механизма

При ходьбе робота неучтенные в упрощенной математической модели динамические эффекты оказывают сильное влияние на характер его движения. Так, при подъеме переносимой стопы силы инерции «утягивают» робот в сторону переносимой стопы, опускание переносимой стопы приводит к обратному результату: робот «выносит» наружу. Одним из способов преодоления этих возмущений является придание роботу завышенных скоростей и ускорений в поперечном направлении, но при этом происходит сильное раскачивание корпуса и трудно получить устойчивую ходьбу с требуемыми скоростями перемещения в продольном направлении. Вторым способом компенсации сил инерции, возникающих при перемещении переносимой стопы по вертикали, является встречное угловое перемещение корпуса ДШР. При повороте корпуса возникает инерционный момент, направленный вокруг оси поворота, в направлении, противоположном направлению углового ускорения корпуса. Величина момента определяется инерционными свойствами исполнительного механизма и величиной придаваемого углового ускорения корпуса. Варьируя параметры соответствующего элемента движения, можно в значительной степени компенсировать возмущения при подъеме и опускании переносимой стопы. Этот прием позволил уменьшить взаимосвязь подъема и опускания переносимой стопы с переносом центра масс робота в поперечном направлении.

Для определения требуемого углового положения корпуса вводится вспомогательная система координат направления движения (СКНД) (рис. 8). Это подвижная система координат, назначаемая по следующим правилам:

— начало СКНД находится в середине проекции на горизонтальную плоскость вектора, соединяющего центры стоп;

— ось 02 СКНД направлена вертикально;

Рис. 8. Формирование заданного углового положения корпуса

— ось ОХ СКНД совпадает с направлением движения ДШР в данный момент времени, которое задается генератором траектории.

Положение СКНД робота характеризуется матрицей Твд:

Т =

нд

/0 X,,

0 д

0 Р нд А

0 У

нд

0г

0 Р нд А

Первоначальное угловое положение корпуса вд Я'6 = ( вд X6 вд у6 нд £6) определяется относительно СКНД направлением от начала СКНД до центра масс. Плоскость 0X2 системы координат корпуса совпадает с плоскостью 0У2 СКНД, и орты осей системы координат корпуса, определяющей его первоначально задаваемую ориентацию, определяются в СКНД выражениями:

нд =/

¿6

нд г'

цмуг

нд г'

цмуг

(У6 = -I

Вд — вд V' V нд X 6 = У 6 Х г 6,

цмУг = ТжН ■ • °ГЦМ; 0ГЦМ — требуемое

где г

положение центра масс, определяемое гене

ратором траектории движения ДШР, выраженное в абсолютной системе координат; ТУН = Ша£(0 КН 1 1) — матрица, определяющая проекцию центра масс робота на плоскость 0У2 системы координат направления движения робота, учитывающая недоворот корпуса. (Здесь и далее при преобразовании векторов с помощью матрицы Т4 х4 используются обобщенные векторы.)

Недоворот корпуса относительно определенного начального углового положения вносится для уменьшения углового раскачивания робота при ходьбе. Определяется коэффициентом недоворота КН, который выступает в качестве варьируемого параметра ходьбы и может меняться при синтезе походки.

Для нахождения итогового требуемого углового положения корпуса, определяемого мат-

0д

рицей

6

осуществляем его отклонение от

первоначально заданной ориентации, определяемой матрицей 44Я'6, на углы доворота а^), Р^), у^), задаваемые генератором траектории движения ДШР, относительно осей 44 X6

44 У6,

г6 соответственно.

0 » — 0 » Д6 = Дцд

\

ео8(у) — 8т(у) 0 8т(у) ео8(у) 0

0

0

1

X

X

X

ео8(Р) 0 8т(Р) 0 1 0 — 8Ш(Р) 0 ео8(Р)

А

X

1 0 0 ео8(а) 0 81п(а)

0

— 8т(а) ео8(а)

А

1Д'.

Введенные углы доворота a(t), P(t), y(t) входят в качестве компонент в вектор положения ДШР л (1).

Авторами был проведен синтез траектории движения лабораторного ДШР при динамической ходьбе с разгоном на начальном участке траектории. В ходе синтеза были определены параметры элементов ФИПП, обеспечивающие ходьбу ДШР с характеристиками, приведенными в табл. 1.

Таблица 1

Характеристики динамической ходьбы

Средняя скорость движения ДШР 0,08 м/с

Величина подъема переносимой стопы 0,04 м

Величина продольного переноса переносимой стопы за один шаг 0,4 м

Расстояние между центрами стоп в поперечном направлении 0,35 м

Величина переноса ЦМ ДШР в поперечном направлении 0,29 м

Максимальный угол дополнительного отклонения корпуса ДШР от заданного углового положения во фронтальной плоскости 0,05 рад

На рис. 9 показаны ФИПП, определяющие траекторию движения ДШР. Параметры элементов ФИПП (2) приведены в табл. 2.

Таблица 2

Параметры функций изменения компонент вектора положения ДШР

^ЦМХ 0 с T цмх 2,48 c ASx 0,2 м

^-Цму 0,5 с T цму 2 c Асму 0,73 м

T bSup 2,25 с T -'Sup 0,5 с ASup 0,02 м

^¿down 3.25 с T ¿down1 0,75 с AÄdown1 0,016 м

tb5down2 4 с T ¿down2 0,5 с AÄdown2 0,004 м

ba кх ге(иГП 2,25 с ^a^xreturn 0,5 с ^aXform 0,025 рад

^bareturn 3,25 с T areturn 1,25 с AaXreturn 0,025 рад

Значения КВП из кадров траектории движения передаются в программные модули решения обратной задачи кинематики, где траектория движения ДШР в декартовом пространстве преобразуется в траекторию движения в пространстве обобщенных координат. При этом дополнительно учитываются довороты стоп относительно первоначально заданных угловых положений, обеспечивающие требуемые моменты для стабилизации движения ДШР при ходьбе.

Выводы

Предложенная методика синтеза алгоритмов управления движением исполнительного механизма робота, основанная на представлении изменения переменных его состояния

L, и а, рад 0,6f—г

34 36 38 40 42 44 46 48

U с

Рис. 9. Траектория движения ДШР в декартовом пространстве при ходьбе: Хцм, Уцм — координаты центра масс робота; х12, x18, у12, у 18, z12, z18 — координаты начал систем координат стоп; акх — угол наклона корпуса во фронтальной плоскости

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

Для компенсации инерционного воздействия неучтенных в упрощенной модели динамических эффектов необходимо управлять дополнительным изменением пространственной конфигурации исполнительного механизма ДШР.

Литература

1. Исследование принципов управления и стабилизации двуногих шагающих роботов в режиме динамической ходьбы: Отчет по проекту РНП 2.1.2.730. Федеральное агентство по образованию / Руководитель А.К. Ковальчук. Исполнители: Медведев В.С., Шахтарин Б.И., Семенов С.Е. и др. РК № 01200703234. М., 2007. 134 с.

2. Ковальчук А.К., Кулаков Д.Б., Семенов С.Е. Управление исполнительными системами двуногих шагающих роботов. Теория и алгоритмы / Под ред. А.К. Ковальчука. М.: Изд-во МГОУ, 2007. 160 с.

3. Ковальчук А.К., Кулаков Д.Б., Семенов С.Е., Ломакин В.О. Синтез функций управления ходьбой двуногого шагающего робота при помощи упрощенной математической модели // Наука и образование: Электронное научно-техническое издание. М., 2010. № 1.

4. Loffler K., Gienger M., Pfeiffer F. Sensor and Control Design of a Dynamically Stable Biped Robot // Robotics and Automation. 2003: Proceeding. ICRA apos; 03. IEEE International Conference. W.t. 2003. Vol. 1. Issue 14—19. P. 484—490.

Статья поступила в редакцию 16.02.2011 г.