Исследование динамики исполнительного механизма экзоскелета нижних конечностей с учётом реакций опорной поверхности

Тип работы:
Реферат
Предмет:
ТЕХНИЧЕСКИЕ НАУКИ


Узнать стоимость

Детальная информация о работе

Выдержка из работы

Наука к Образование
МГТУ им. Н.Э. Баумана
Сетевое научное издание
Наука и Образование. МГТУ им. Н. Э. Баумана. Электрон. журн. 2014. № 12. С. 128−136.
ISSN 1994−0408
Б01: 10. 7463/0815. 9 328 000
Представлена в редакцию: Исправлена:
© МГТУ им. Н.Э. Баумана
УДК 621. 865. 8
Исследование динамики исполнительного механизма экзоскелета нижних конечностей с учётом реакций опорной поверхности
*
Верейкин А. А., Ковальчук А. К., Каргинов Л. А.
##. ##. 2014 ##. ##. 2014
1МГТУ им. Н. Э. Баумана, Москва, Россия
Представлены результаты исследования динамики исполнительного механизма экзоскелета нижних конечностей. Уравнения динамики получены с использованием принципа Д'-Аламбера. Для описания кинематической структуры исполнительного механизма применялись системы координат Денавита-Хартенберга и элементы теории графов. На базе разработанной в программном комплексе SolidWorks 3−0 модели исполнительного механизма экзоскелета найдены массо-инерционные параметры его звеньев. Определены законы изменения обобщённых координат в сочленениях экзоскелета от времени, соответствующие наиболее нагруженным случаям типовых движений человека. Получены значения моментов и мощностей в степенях подвижности исполнительного механизма как без учёта, так и с учётом реакций опорной поверхности шагания.
Ключевые слова: робототехника, древовидная кинематическая структура, расчёт привода, экзоскелет, человек-машина, уравнение динамики робота, динамический анализ, математическое моделирование
Введение
Робототехника на сегодняшний день является одной из самых динамично развивающихся отраслей научного знания, что находит своё объяснение в широте и разнообразии областей её возможного применения. При этом значительная часть исследований, имеющих отношение к робототехнике, посвящена шагающим роботам (ШР), среди несомненных преимуществ которых стоит упомянуть более высокие в сравнении с колёсными или гусеничными машинами проходимость и энергоэффективность в условиях передвижения по пересечённой местности. Ещё одним достоинством ШР является возможность частичного или полного замещения человеческого труда роботами в условиях, изначально созданных для человека (здания, кабины различной техники и т. п.).
Одним из наиболее перспективных видов шагающих роботов в плане частичной замены человеческого труда являются роботизированные экзоскелеты [1−5], разгружающие человека-оператора и повторяющие его биомеханику. Широкое распространение получили экзоскелеты, оснащённые электрогидравлическими приводами в силу высокой энергоёмкости и приемлемых динамических свойств последних. Однако практическое применение экзоскелетов ограничено в силу отсутствия на данный момент подходящего бортового источника энергии, способного поддерживать их автономность на приемлемом уровне. Исследования [6−10] посвящены вопросам создания силовых установок мобильных роботов.
Ещё одной трудностью в процессе разработки экзоскелетов является проблема управления ими человеком-оператором в силу того, что человеческий организм представляет собой сложную нелинейную систему с распределёнными параметрами. Данный вопрос рассматривается в работах [11−13].
По результатам обзора, приведённого в [14], можно сделать вывод, что в РФ наблюдается значительное отставание в этой области техники. Актуальной ведущейся разработкой, насколько можно судить по [15], является проект Ехоа1! е1 НИИ Механики МГУ им. М. В. Ломоносова. Таким образом, разработка исполнительного механизма (ИМ) экзоске-лета и исследование его параметров и характеристик представляют собой актуальные задачи. В настоящем исследовании используется метод построения математических моделей ИМ, прошедший апробацию при расчёте параметров шагающих роботов как с древовидной [16−20], так и с линейной [21] кинематическими структурами. В качестве приводов выбраны электрогидравлические следящие приводы (ЭГСП) в силу указанных выше причин.
В работах [22−24] представлены результаты синтеза кинематической схемы экзоскеле-та нижних конечностей, проведённого по методике, изложенной в [25], определены диапазоны изменения обобщённых координат сочленений экзоскелета, оснащённых ЭГСП, назначены системы координат (СК) для которых определены параметры Денавита-Хартенберга [26], получены граф и матрица достижимости.
В первом разделе записано уравнение динамики ИМ экзоскелета, на базе которого проводится математическое моделирование. Во втором разделе приведены исходные данные, необходимые для проведения моделирования: параметры Денавита-Хартенберга основных и вспомогательных СК ИМ экзоскелета, массо-инерционные параметры его звеньев, точки приложения сил реакции опоры, законы изменения обобщённых координат в степенях подвижности ИМ экзоскелета, оснащённых приводами. В третьем разделе пред-
ставлены результаты моделирования уравнения динамики ИМ экзоскелета, включающие в себя значения моментов и мощностей в его сочленениях как без учёта, так и с учётом реакций опорной поверхности.
1. Уравнение динамики исполнительного механизма экзоскелета
Уравнение динамики ИМ экзоскелета с учётом воздействия внешних сил и моментов, а также внешних наложенных связей запишем в следующем виде [27, 28]:
Аъ) — 4я (я)
vJt (я) 0
А я ^ (В (я, я)^
У
, К
V f У
+
Р (я)

0 г Л
0
ГтЛ
V 0 у
(1)
где q — вектор обобщённых координат ИМ-
А (я), В (я, я) — матричные коэффициенты-
Jt (я) — Якобиева матрица, определяющая кинематические ограничения, накладываемые на ИМ связями-
JvR (я) — транспонированная Якобиева матрица, определяющая кинематические ограничения, накладываемые кинематикой ИМ-
0 Ку — блочная матрица сил и моментов реакций связи-
Р (я) — вектор-столбец, определяемый из кинематических выражений-
Дя) — Якобиева матрица размерностью N X (2 ¦ Ы), объединяющая матричные коэффициенты внешних сил и внешних моментов, наложенных на ИМ-
0 — блочный вектор внешних сил и моментов, приложенных к звеньям ИМ-
т — матрица моментов, развиваемых приводами ИМ.
Значения компонентов уравнения (1) определяются в соответствии с [29, 30].
2. Исходные данные для расчёта динамики ИМ экзоскелета
Исходными данными для расчета компонентов уравнения (1) являются:
¦ [м], аг [м], а.1 [рад] - параметры Денавита-Хартенберга основных и вспомогательных СК экзоскелета, параметры /(г) и т (г) основных систем координат, где /(г) — номер звена-отца [31] звена г и т (1) — параметр, показывающий, каким по счёту сыном [31] является звено г для звена /(г), а также параметры Ыг и Ык вспомогательных систем координат, где Ыг — номер основной СК звена, Ык — номер вспомогательной СК-
¦ диапазоны изменения обобщённых координат сочленений ИМ-
¦ параметры звеньев ИМ: масса, тензоры инерции, координаты центров масс звеньев, получаемые по результатам 30-моделирования конструкции-
¦ координаты точек приложения внешних сил-
¦ координаты точек приложения внешних связей.
2.1. Параметры Денавита-Хартенберга
Синтезированная в [22, 23] кинематическая схема и соответствующая ей 3-D модель ИМ экзоскелета представлены на рис. 1 и 2 соответственно. Значения параметров Денавита-Хартенберга для основных и вспомогательных СК ИМ экзоскелета указаны в таблицах 1 и 2.
Рис. 1. Кинематическая схема ИМ экзоскелета
Рис. 2. 3-D модель ИМ экзоскелета
Таблица 1. Значения параметров Денавита-Хартенберга основных СК ИМ экзоскелета
№ СК 0, [рад] d, [м] a, [м] а, [рад] f (i) ns (i)
1 0 0 0 0 0 1
2 0 0 0 0 1 1
3 0 0 0 0 2 1
4 0 0 0 0 3 1
5 0 0 0 0 4 1
6 — п/2 0 0,238 — п/2 5 1
7 0 0,198 0,167 п/2 6 1
8 п/2 — 0,107 0 — п/2 7 1
9 п/2 0 0,526 0 8 1
10 0 0 0,674 0 9 1
11 0 0 0,094 п/2 10 1
12 п/2 0 0 п/2 11 1
13 0 0,101 0 — п/2 12 1
14 п 0,198 — 0,167 п/2 6
15 п/2 — 0,107 0 — п/2 13 1
16 п/2 0 0,526 0 14 1
17 0 0 0,674 0 15 1
18 0 0 0,094 п/2 16 1
19 п/2 0 0 п/2 17 1
20 0 0,101 0 — п/2 18 1
Таблица 2. Значения параметров Денавита-Хартенберга вспомогательных СК ИМ экзоскелета
№ СК 0, [рад] и, [м] а, [м] а, [рад] N1
6,2 п 0 0,476 0 6 2
2.2. Массо-инерционные параметры звеньев ИМ
Массо-инерционные параметры звеньев ИМ экзоскелета были получены по результатам 3−0 моделирования в программном комплексе SolidWorks. На рис. 3−10 показаны модели звеньев ИМ и в таблицах 3−9 приведены соответствующие им массо-инерционные параметры. Отметим, что грузоподъёмность разрабатываемого ИМ экзоскелета составляет 50 кг в предположении, что груз располагается в геометрическом центре звена 6. Масса силовой установки составляет 15 кг.
Таблица 3. Массо-инерционные параметры звеньев 6, 7 и 14
Звено 6 Звено 7 Звено 14
Масса, [кг] М = 50 М = 2. 40 525 758 М = 2. 40 525 758
Координаты центра масс, [м] X = -0. 23 510 159 X = -0. 4 822 968 X = 0. 4 822 968
V = -0. 22 231 140 V = -0. 10 362 463 V = -0. 10 362 463
г = -0. 12 147 542 г = -0. 6 311 183 г = -0. 6 311 183
Тензоры инерции, [кг ¦ м2] 4 = 7. 12 970 399 1х = 0. 2 850 829 1х = 0. 2 850 829
1у = 5. 96 302 003 1у = 0. 1 966 455 1у = 0. 1 966 455
4 = 9. 36 554 832 4 = 0. 3 135 668 4 = 0. 3 135 668
1ху = -0. 1 278 172 1ху = 0. 1 087 169 1ху = -0. 1 087 169
1Х2 = -0. 1 018 372 1Х2 = 0. 622 898 1Х2 = -0. 622 898
1уг = 0. 2 476 764 1уг = 0. 1 208 787 1у2 = 0. 1 208 787
Таблица 4. Массо-инерционные параметры звеньев 8 и 15
Звено 8 Звено 15
Масса, [кг] М = 0. 83 994 875 М = 0. 83 994 875
Координаты центра масс, [м] X = 0. 473 085 X = 0. 473 085
V = -0. 6 957 595 V = -0. 6 957 595
г = 0. 16 564 г = -0. 16 564
Тензоры инерции, [кг ¦ м2] 1х = 0. 213 518 4 = 0. 213 518
1у = 0. 52 669 1у = 0. 52 669
4 = 0. 240 389 4 = 0. 240 389
1ху = 0. 24 127 1ху = 0. 24 127
1хг = 0. 187 1хг = -0. 187
1уг = 0. 240 1уг = -0. 240
т
Рис. 5. 3-Б модель звена 8
Таблица 5. Массо-инерционные параметры звеньев 9 и 16
Звено 9 Звено 16
Масса, [кг] М = 8. 91 380 755 М = 8. 91 380 755
Координаты центра масс, [м] X = -0. 31 796 049 X = -0. 31 796 049
V = -0. 5 298 032 V = -0. 5 298 032
г = 0. 285 154 г = -0. 285 154
4 = 0. 2 483 047 4 = 0. 2 483 047
1У = 0. 14 443 780 1У = 0. 14 443 780
Тензоры инерции, [кг ¦ м2] 4 = 0. 16 161 733 4 = 0. 16 161 733
1ху = 0. 257 638 1Ху — 0. 257 638
хг = 0. 322 018хг -0. 322 018
1уг = 0. 114 324 1уг ~ -0. 114 324
Рис. 6. 3−0 модель звена 16 Рис. 7. 3−0 модель звена 17 Рис. 8. 3−0 модель звена
11
Таблица 6. Массо-инерционные параметры звеньев 10 и 17
Звено 10 Звено 17
Масса, [кг] М = 12. 78 978 439 М = 12. 78 978 439
Координаты центра масс, [м] X = -0. 43 625 373 X = -0. 43 625 373
V = -0. 415 647 V = -0. 415 647
г = 0. 161 494 г = -0. 161 494
Тензоры инерции, [кг ¦ м2] 1х = 0. 7 880 506 /ж = 0. 7 880 506
1у = 1. 6 478 425 1у = 1. 6 478 425
4 = 1. 13 450 491 4 = 1. 13 450 491
1ху = -0. 23 544 191 1ху = -0. 23 544 191
1хг = 0. 595 373 1хг = -0. 595 373
1у2 = 0. 114 324 1у2 = -0. 114 324
Таблица 7. Массо-инерционные параметры звеньев 11 и 18
Звено 11 Звено 18
Масса, [кг] М = 0. 75 006 227 М = 0. 75 006 227
Координаты центра масс, [м] X = -0. 7 648 757 X = -0. 7 648 757
V = 0. 13 199 V = -0. 13 199
г = 0. 929 483 г = 0. 929 483
Тензоры инерции, [кг ¦ м2] 4 = 0. 62 512 1х = 0. 62 512
1у = 0. 134 389 1у = 0. 134 389
4 = 0. 105 048 4 = 0. 105 048
1ху = 0. 8 1ху = -0. 8
1хг = -0. 14 631 1хг = -0. 14 631
1уг = -0. 86 1уг = 0. 86
Рис. 9. 3−0 модель звена 12 Рис. 10. 3−0 модель звена 13
Таблица 8. Массо-инерционные параметры звеньев 12 и 19
Звено 12 Звено 19
Масса, [кг] М = 0. 84 265 474 М = 0. 84 265 474
Координаты центра масс, [м] X = 0. 0 X = 0. 0
V = -0. 5 627 V = -0. 5 627
г = 0. 4 704 563 г = 0. 4 704 563
Тензоры инерции, [кг ¦ м2] 4 = 0. 108 224 4 = 0. 108 224
1у = 0. 102 586 1у = 0. 102 586
4 = 0. 38 008 4 = 0. 38 008
1ху = -0. 1 1ху = -0. 1
1хг = 0. 13 1хг = 0. 13
1уг = -0. 6 1уг = -0. 6
Таблица 9. Массо-инерционные параметры звеньев 13 и 20
Звено 13 Звено 20
Масса, [кг] М = 2. 40 532 807 М = 2. 40 532 807
Координаты центра масс, [м] X = 0. 2 536 075 X = -0. 2 536 075
V = 0. 2 160 510 V = 0. 2 160 510
Z = 0. 3 175 562 Z = 0. 3 175 562
Тензоры инерции, [кг ¦ м 2] 4 = 0. 1 470 942 Ix = 0. 1 470 942
1у = 0. 1 737 613 ly = 0. 1 737 613
4 = 0. 680 497 4 = 0. 680 497
1ху = -0. 151 431 Ixy = 0. 151 431
Ixz = 0. 383 638 Ixz = -0. 383 638
lyz = -0. 189 032 Iyz = -0. 189 032
Необходимо заметить, что силы реакции опоры приложены в двух точках на левой и правой стопах, которые находятся на осевой линии стопы на мысках, что соответствует наименее благоприятному случаю нагружения.
2.3. Законы изменения обобщённых координат в степенях подвижности ИМ экзоскелета, оснащённых приводами
С целью определения законов изменения обобщённых координат в степенях подвижности ИМ экзоскелета, оснащённых приводами, были рассмотрены характерные типовые движения человека-оператора, совершаемые за счёт изменения обобщённой координаты (угла & lt-р н, а гр (t)) во времени в соответствующем сочленении в определённом ранее диапазоне (таблица 10) согласно гармоническому закону вида:
Фнагр (0 = & lt-Ра нагр & quot- 5 (шt) + & lt-РH нагр,
2)
где & lt-р, а н, а гр — амплитуда изменения угла в сочленении, [град]- h н, а гр — начальное значение угла в сочленении [град]-
ш — скорость изменения угла в сочленении [град/c].
Таблица 10. Диапазоны изменения обобщённых координат в сочленениях экзоскелета, оснащённых приводами
Сочленения Значение диапазона
Чт & lt-?14 [135°- 210°]
& lt-?9-<-?16 [120°- 184°]
& lt-7l0- & lt-?17 [90°- 180°]
Чи'- 4i в [140°- 230°]
Чтобы определить закон изменения обобщённой координаты для спинного сочленения ИМ, соответствующий наиболее нагруженному случаю, были рассмотрены следующие типовые движения:
• изменение угла крена корпуса, при котором человек-оператор остаётся на месте-
• поочерёдные движения ногами, при которых человек-оператор остаётся на месте и положение корпуса с течением времени не изменяется. В этом случае, во-первых, несущественна инерционная составляющая нагрузки, а во-вторых, масса ноги исполнительного механизма мала в сравнении с суммарной массой корпуса и груза, что делает исследование данного типа движения необязательным-
• поочерёдные движения ногами, совершаемые в процессе бокового перемещения, при этом вертикальная координата корпуса меняется во времени.
Результаты расчёта показали, что наиболее тяжёлому случаю нагружения соответствует изменение крена корпуса ИМ, описываемое следующим законом изменения угла в спинном сочленении:
ЫО = 914(0 = 37. 5° ¦ 5 т (128. 572 ¦ 0 + 172. 5°, 135° & lt- д7(0 = д14(0 & lt- 210°.
С целью определения закона изменения обобщённой координаты в бедренном сочленении ИМ, соответствующего наихудшему случаю нагружения, были рассмотрены следующие типовые движения:
• одновременное изменение угла для двух ног, при котором человек-оператор остаётся на месте, сопровождаемое изменением тангажа корпуса ИМ, т. е. наклоны корпуса вперёд-назад-
• поочерёдные махи ногами, при которых человек-оператор остаётся на месте и положение корпуса с течением времени не изменяется. В этом случае, во-первых, мала инерционная составляющая нагрузки, а во-вторых, масса ноги исполнительного механизма мала в сравнении с суммарной массой корпуса и груза, что делает исследование данного типа движения необязательным-
• поочерёдные махи ногами, совершаемые в процессе ходьбы, при этом вертикальная координата корпуса ИМ меняется во времени.
Результаты расчёта показали, что наименее благоприятному случаю нагружения соответствуют поочерёдные махи ногами, совершаемые в процессе ходьбы, при которых закон изменения угла в бедренном сочленении ИМ задаётся соотношением: = Ч1б (0 = 32° ¦ 5 т (156. 532 ¦ 0 + 152°,
1 2 0° & lt- д9 (0 = д1б (0 & lt- 1 84°.
Для коленного сочленения в качестве типового движения рассматривалось приседание на одной ноге с полным контактом поверхности стопы ИМ и опорной поверхности. Закон изменения угла для коленного сочленения примет вид: д1 о (0 = ?1 7 (0 = 45° ¦ 5 т (143. 984 ¦ О + 1 3 5°, 9 0 ° & lt- д 1 о (0 = д 1 7 (0 & lt- 1 8 0°.
Для щиколоточного сочленения рассматривались движения типа поднятие — опускание человека на мысках. Итак, закон изменения угла для щиколоточного сочленения ИМ запишется следующим образом:
д11 (0 = д1 8 (0 = 45° ¦ 5 т (1 3 3. 3 2 7 ¦ 0 + 9 5°- 140° & lt- д1 1 (0 = д1 8 (0 & lt- 2 30°.
3. Результаты
В результате моделирования уравнения динамики (1) в программе моделирования древовидных ИМ ТТТР с учётом внешних наложенных связей [29] были получены численные значения его элементов. Моменты и мощности в сочленениях, оснащённых ЭГСП, представлены на рис. 11−18. При этом моменты и мощности, рассчитанные без учёта реакции опоры обозначены как Мг и Ыг, а моменты и мощности, рассчитанные с учётом реакции опоры, — Мгр& gt- и соответственно, где г — номер сочленения.
1×103
? МЖ
— -4×10-
— 5×103-----
-бхЮ3'-------
0 0.1 0.2 0.3 0. 4
т
[с]
Рис. 11. Зависимость момента в сочленении 7 от времени без учёта и с учётом реакции опоры
[с]
Рис. 12. Зависимость момента в сочленении 9 от времени без учёта и с учётом реакции опоры
[с]
Рис. 13. Зависимость момента в сочленении 10 от времени без учёта и с учётом реакции опоры
[с]
Рис. 14. Зависимость момента в сочленении 11 от времени без учёта и с учётом реакции опоры
[с]
Рис. 15. Зависимость мощности в сочленении 7 от времени без учёта и с учётом реакции опоры
Рис. 16. Зависимость мощности в сочленении 9 от времени без учёта и с учётом реакции опоры
Рис. 17. Зависимость мощности в сочленении 10 от времени без учёта и с учётом реакции опоры
Рис. 18. Зависимость мощности в сочленении 11 от времени без учёта и с учётом реакции опоры
На рис. 19 приведён график суммарной мощности, развиваемой в сочленениях ИМ экзоскелета как без учёта, так и с учётом реакции опоры. Критерий минимальной суммарной (затрачиваемой) мощности может служить одним из базовых критериев при оценке конструкции.
Рис. 19. Зависимость суммарной мощности в сочленении ИМ экзоскелета без учёта и с учётом реакции
опоры
Как видно из приведённых рисунков, часто учёт реакций опоры приводит к уменьшению значений моментов и мощностей в сочленениях, что может быть объяснено разнона-правленностью сил тяжести и реакции опоры.
Заключение
Одним из критических параметров экзоскелета является потребляемая мощность, которая, в свою очередь, определяет степень его автономности. Стремление к уменьшению потребляемой мощности заставляет разработчиков проектировать компактные силовые установки, исследовать различные виды приводов и способы энергоэффективного управления ими. Очевидна необходимость применения методов, позволяющих оперативно определять энергетические параметры ИМ как функцию от множества факторов (кинематической схемы, конструкции, грузоподъёмности, законов управления и т. д.).
Предложенный метод расчёта энергетических параметров экзоскелета, используемый совместно с описанным в работах [22, 23] методом синтеза кинематической схемы ИМ экзоскелета, позволяет на протяжении всего процесса разработки изделия, начиная с этапа эскизного проектирования, когда возможно оперативное изменение принятых решений, с минимальными трудозатратами прогнозировать энергетические характеристики приводов, отслеживать влияние изменений, вносимых в ИМ, и оценивать различные варианты конструкции.
Представленные в настоящей статье результаты являются основой дальнейших исследований, направленных на создание экзоскелета.
Список литературы
1. Vukobratovic M. Active exoskeletal systems and beginning of the development of human-oid robotics // In: Monograph of ANS: Academy of Nonlinear Sciences. Advances in Nonlinear Sciences II — Sciences and Applications. Vol. 2. Belgrade, 2008. P. 329−348.
2. Hanlon M. Raytheon XOS 2: second generation exoskeleton. Second-Generation Robotics Suit, United States of America // Gizmag: website. Режим доступа: www. gizmag. com/raytheon-significantly-progresses-exoskeleton-design/16 479 (дата обращения 20. 10. 2014).
3. Kazerooni H., Steger R. The Berkeley Lower Extremity Exoskeleton // ASME Journal of Dynamics Systems, Measurements, and Control. 2005. Vol. 128, is. 1. P. 14−25. DOI: 10. 1115/1. 2 168 164
4. Kazerooni H., Steger R., Huang L. Hybrid Control of the Berkeley Lower Extremity Exo-skeleton (BLEEX) // The International Journal of Robotics Research. 2006. Vol. 25, no. 5−6. P. 561−573. DOI: 10. 1177/278 364 906 065 505
5. Sankai Y. HAL: Hybrid Assistive Limb based on Cybernics // Sankai Laboratory. Cybernics Group: website. Режим доступа: http: //sanlab. kz. tsukuba. ac. jp/sonota/ISSR_Sankai. pdf (дата обращения 20. 10. 14).
6. Raade J.M., Amundson K.R., Kazerooni H. Development of Hydraulic-Electric Power Units for Mobile Robots // Proceedings of IMECE 2005. ASME International Mechanical Engineering Congress and Exposition (5−11 November, 2005, Orlando, Florida USA). ASME, 2005. Paper no. IMECE2005−80 138. P. 27−34. DOI: 10. 1115/IMECE2005−80 138
7. Raade J.M., Kazerooni H., McGee T.G. Analysis and Design of a Novel Power Supply for Mobile Robots // Proceedings of IEEE 2004. International Conference on Robotics & amp- Automation (April, 2005, New Orleans, LA USA). Vol. 5. IEEE Publ., 2004. P. 4911−4917. DOI: 10. 11Q9/ROBOT. 2004. 1 302 496
8. Barth E. J, Riofrio J.A. Design of a Free Piston Pneumatic Compressor as a Mobile Robot Power Supply // Proceedings of the 2005 IEEE International Conference on Robotics and Automation (18−22 April, 2005, Barcelona, Spain). IEEE Publ., 2005. P. 235−240. DOI: 10. 11Q9/ROBOT. 2005. 1 570 125
9. Riofrio J.A., Barth E.J. A Free Piston Compressor as a Pneumatic Mobile Robot Power Supply: Design, Characterization and Experimental Operation // International Journal of Fluid Power. 2007. Vol. 8, no. 1. P. 17−28.
10. Amundson K.R., Raade J.M., Harding N., Kazerooni H. Hybrid Hydraulic-Electric Power Unit for Field and Service Robots // 2005 IEEE/RSJ International Conference on Intelligent Robots and Systems (2−6 August, 2005, Shaw Convention Center, Edmonton, AB, Canada). IEEE Publ., 2005. P. 3453−3458. DOI: 1Q. 11Q9/IROS. 2005. 1 545 452
11. Воробьев Е. И., Скворчевский А. К., Сергеев А. М. Проблемы создания алгоритмов управления антропоморфными биопротезами рук и ног // Медицина и высокие технологии. 2012. № 1. С. 7−12.
12. Скворчевский А. К., Воробьев Е. И., Скворчевский К. А. Концепция развития технических средств реабилитации (ТСР) людей с ограниченными физическими возможностями на основе антропоморфных биопротезов рук // Медицина и высокие технологии. 2012. № 2. С. 7−22.
13. Скворчевский А. К., Воробьев Е. И., Скворчевский К. А., Сергеев А. М., Акентьев А. А., Ковалев Н. С. Технические средства реабилитации (ТСР) на основе антропоморфных биопротезов рук и ног // Информатизация и связь. 2013. № 4. С. 43−58.
14. Верейкин А. А. Расчёт исполнительных гидроцилиндров экзоскелета // Молодёжный научно-технический вестник МГТУ им. Н. Э. Баумана. Электронный журнал. 2013. № 5. Режим доступа: http: //sntbul. bmstu. ru/doc/569 290. html (дата обращения 20. 10. 2014).
15. Лавровский Э. К., Письменная Е. В. О регулярной ходьбе экзоскелетона нижних конечностей при дефиците управляющих воздействий // Российский журнал биомеханики. 2014. Т. 18, № 2. С. 208−225.
16. Ковальчук А. К., Семенов С. Е., Каргинов Л. А., Кулаков Д. Б., Верейкин А. А., Яроц
B.В. Выбор кинематической структуры и исследование динамики древовидного исполнительного механизма робота-треножника // Инженерный вестник МГТУ им. Н. Э. Баумана. Электрон. журн. 2013. № 10. Режим доступа: http: //engbul. bmstu. ru/doc/654 481. html (дата обращения (20. 10. 2014).
17. Ковальчук А. К., Верейкин А. А., Каргинов Л. А., Цепкова А. А. Проектирование и расчёт мощности приводов медицинского робота с учётом динамики его исполнительного механизма // XII Международная научно-практическая конференция «Техника и технология: новые перспективы развития» (25. 02. 2014): матер. М.: & quot-Спутник +& quot-, 2014.
C. 28−41.
18. Ковальчук А. К., Ромашко А. М., Верейкин А. А., Каргинов Л. А., Малякина Е. А., Ста-тива В. А. Кинематический и динамический анализ исполнительного механизма шес-тиногого шагающего робота // XIV Международная научно-практическая конференция & quot-Современное состояние естественных и технических наук& quot- (14. 03. 2014): матер. М.: & quot-Спутник +& quot-, 2014. С. 60−79.
19. Ковальчук А. К. Выбор кинематической структуры и исследование древовидного исполнительного механизма робота-собаки // Известия вузов. Машиностроение. 2011. № 8. С. 65−73.
20. Ковальчук А. К. Проектирование и исследование динамики древовидного исполнительного механизма шестиногого шагающего робота // Естественные и технические науки. 2014. № 4 (72). С. 104−110.
21. Kovalchuk A.K. Designing drives of a medical robot actuator // Life Science Journal. 2014. № 11 (11s). P. 337−340.
22. Верейкин А. А., Ковальчук А. К., Кулаков Д. Б., Семёнов С. Е., Каргинов Л. А., Кулаков Б. Б, Яроц В. В. Синтез кинематической схемы исполнительного механизма экзоскеле-та // XIII Международная научно-практическая конференция & quot-Актуальные вопросы науки& quot- (25. 04. 2014): матер. М.: & quot-Спутник +& quot-, 2014. С. 68−76.
23. Верейкин А. А., Ковальчук А. К., Кулаков Д. Б., Семенов С. Е. Анализ и выбор кинематической структуры исполнительного механизма экзоскелета // Наука и образование. МГТУ им. Н. Э. Баумана. Электронный журнал. 2014. № 7. С. 72−93. DOI: 10. 7463/0714. 717 676
24. Верейкин А. А., Ковальчук А. К., Кулаков Д. Б., Семёнов С. Е., Каргинов Л. А., Кулаков Б. Б, Яроц В. В. Динамика исполнительного механизма экзоскелета // XIII Международная научно-практическая конференция & quot-Техника и технология: новые перспективы развития& quot- (26. 05. 2014): матер. М.: & quot-Спутник +& quot-, 2014. C. 5−16.
25. Ковальчук А. К. Метод синтеза древовидных исполнительных механизмов шагающих роботов // Естественные и технические науки. 2014. № 3 (71). С. 127−130.
26. Denavit J., Hartenberg R.S. A kinematic notation for Lower-Pair Mechanisms Based on Matrices // J. Appl. Mech. 1955. Vol. 77. P. 215−221.
27. Ковальчук А. К., Кулаков Б. Б., Кулаков Д. Б., Семенов С. Е., Яроц В. В. Основы теории исполнительных механизмов шагающих роботов. М.: Изд-во Рудомино, 2010. 170 с.
28. Ковальчук А. К. Метод математического описания кинематики и динамики древовидных исполнительных механизмов шагающих роботов // Естественные и технические науки. 2014. № 5 (73). С. 87−90.
29. Ковальчук А. К., Каргинов Л. А., Кулаков Б. Б., Кулаков Д. Б., Семенов С. Е., Яроц В. В., Верейкин А. А. Моделирование древовидных исполнительных механизмов шагающих роботов с учётом внешних наложенных связей: свидетельство о гос. регистрации программ для ЭВМ № 2 014 612 547 от 28. 02. 2014.
30. Ковальчук А. К., Каргинов Л. А., Кулаков Б. Б., Кулаков Д. Б., Семенов С. Е., Яроц В. В. Программа моделирования древовидных исполнительных механизмов шагающих роботов: свидетельство о гос. регистрации программ для ЭВМ № 2 012 610 398. 10. 01. 2012.
31. Алексеев В. Е., Таланов В. А. Графы и алгоритмы. Структуры данных. Модели вычислений. М.: Изд-во Бином, 2006. 319 с.
Science and Education of the Bauman MSTU, 2014, no. 12, pp. 128−136.
DOI: 10. 7463/0815. 9 328 000
Received: Revised:
##. ##. 2014 ##. ##. 2014
Science^Education
of the Bauman MSTU
ISS N 1994−0408 © Bauman Moscow State Technical Unversity
The Lower Extremities Exoskeleton Actuator Dynamics Research Taking into Account Support Reaction
*
A. A. Vereikin, A.K. Kovalchuk, L.A. Karginov
: Bauman Moscow State Technical University, Moscow, Russia
Keywords: robotics, tree-like kinematic structure, drive design, exoskeleton, human-machine, dynamic equation of robot, dynamics analysis, mathematical simulation
The article shows high relevance of research, aimed at the robotic exoskeleton creation. It defines some problems related to the process of their designing- including a lack of power supply to provide enough autonomy, and difficulties of man-machine complex control. There is a review of literature on the walking robots with tree-like kinematic structure development. This work reflects the continuing investigations, currently conducted by the authors, and relies heavily on the results of previous works, devoted to this subject.
The article presents the exoskeleton dynamics equation, taking into account the impact of external forces and torques, as well as external relations imposed. Using a model of lower extremities exoskeleton developed in SolidWorks software complex, baricentric parameters of the actuator links were found. The different types of movements, committed due to harmonic changes of generalized coordinates in exoskeleton degrees of mobility, equipped with electrohydraulic actuators, were analyzed. The laws of generalized coordinate changes in time, corresponding to the worst case loading, were found. All the necessary input data for the exoskeleton dynamics equation solution were determined.
The numerical values of all components of the dynamics equation were obtained as result of the dynamics equation simulation. In this case, the exoskeleton actuator load capacity was assumed to be 50 kg. The article shows dependences of torque and power in the actuator degrees of mobility on the time, as well as a curve of total capacity of all drives both, ignoring and taking into consideration the support surface reactions. Obtained dependences are the initial data for the calculation of the drive system.
The proposed method for determination of exoskeleton energy parameters allows developer to perform a prompt evaluation of various options for the actuator design in accordance with the selected criteria. As a main evaluation criterion, a minimum consumed power is offered to choose. This, under otherwise equal conditions, provides a maximum autonomy for the projecta-ble exoskeleton.
References
1. Vukobratovic M. Active exoskeletal systems and beginning of the development of humanoid robotic. In: Monograph of ANS: Academy of Nonlinear Sciences. Advances in Nonlinear Sciences II — Sciences and Applications. Vol. 2. Belgrade, 2008, pp. 329−348.
2. Hanlon M. Raytheon XOS 2: second generation exoskeleton. Second-Generation Robotics Suit, United States of America. Gizmag: website. Available at: www. gizmag. com/raytheon-significantly-progresses-exoskeleton-design/16 479, accessed 20. 10. 2014.
3. Kazerooni H., Steger R. The Berkeley Lower Extremity Exoskeleton. ASME Journal of Dynamics Systems, Measurements, and Control, 2005, vol. 128, is. 1, pp. 14−25. DOI: 10. 1115/1. 2 168 164
4. Kazerooni H., Steger R., Huang L. Hybrid Control of the Berkeley Lower Extremity Exoskel-eton (BLEEX). The International Journal of Robotics Research, 2006, vol. 25, no. 5−6, pp. 561−573. DOI: 10. 1177/278 364 906 065 505
5. Sankai Y. HAL: Hybrid Assistive Limb based on Cybernics. Sankai Laboratory. Cybernics Group: website. Available at: http: //sanlab. kz. tsukuba. ac. jp/sonota/ISSR_Sankai. pdf, accessed 20. 10. 14.
6. Raade J.M., Amundson K.R., Kazerooni H. Development of Hydraulic-Electric Power Units for Mobile Robots. Proceedings of IMECE 2005. ASME International Mechanical Engineering Congress and Exposition, 5−11 November, 2005, Orlando, Florida USA). ASME, 2005, paper no. IMECE2005−80 138, pp. 27−34. DOI: 10. 1115/IMECE2005−80 138
7. Raade J.M., Kazerooni H., McGee T.G. Analysis and Design of a Novel Power Supply for Mobile Robots. Proceedings of IEEE 2004 International Conference on Robotics and Automation, April, 2005, New Orleans, LA USA. Vol. 5. IEEE Publ., 2004, pp. 4911−4917. DOI: 10. 1109/ROBOT. 2004. 1 302 496
8. Barth E. J, Riofrio J.A. Design of a Free Piston Pneumatic Compressor as a Mobile Robot Power Supply. Proceedings of the 2005 IEEE International Conference on Robotics and Automation, 18−22 April, 2005, Barcelona, Spain. IEEE Publ., 2005, pp. 235−240. DOI: 10. 1109/ROBOT. 2005. 1 570 125
9. Riofrio J.A., Barth E.J. A Free Piston Compressor as a Pneumatic Mobile Robot Power Supply: Design, Characterization and Experimental Operation. International Journal of Fluid Power, 2007, vol. 8, no. 1, pp. 17−28.
10. Amundson K.R., Raade J.M., Harding N., Kazerooni H. Hybrid Hydraulic-Electric Power Unit for Field and Service Robots. 2005 IEEE/RSJ International Conference on Intelligent Robots and Systems, 2−6 August, 2005, Shaw Convention Center, Edmonton, AB, Canada. IEEE Publ., 2005, pp. 3453−3458. DOI: 1Q. 11Q9/IROS. 2005. 1 545 452
11. Vorob'-ev E.I., Skvorchevskiy A.K., Sergeev A.M. The development of dynamic models and control algorithms bioprostheses anthropomorphic hands and legs. Meditsina i vysokie tekhnologii, 2012, no. 1, pp. 7−12. (in Russian).
12. Skvorchevskiy A.K., Vorob'-ev E.I., Skvorchevskiy K.A. The concept of the development of technical means of rehabilitation people with disabilities on the basis of anthropomorphic hands bioprostheses. Meditsina i vysokie tekhnologii, 2012, no. 2, pp. 7−22. (in Russian).
13. Skvorchevskiy A.K., Vorob'-ev E.I., Skvorchevskiy K.A., Sergeev A.M., Akent'-ev A.A., Kovalev N.S. Technical devices of rehabilitation (TDR) on the basis of anthropomorphic prosthetic hands and feet. Informatizatsiya i svyaz'-, 2013, no. 4, pp. 43−58. (in Russian).
14. Vereikin A.A. Calculation of actuating hydraulic cylinders of exoskeleton. Molodezhnyi nauchno-tekhnicheskii vestnik MGTU im. N.E. Baumana = Youth Science and Technology Herald of the Bauman MSTU, 2013, no. 5. Available at: http: //sntbul. bmstu. ru/doc/569 290. html, accessed 20. 10. 2014. (in Russian).
15. Lavrovskiy E.K., Pis'-mennaya E.V. On the lower limb exoskeleton regular locomotion under input control deficit. Rossiyskiy zhurnal biomekhaniki = Russian Journal of Biome-chanics, 2014, vol. 18, no. 2, pp. 208−225. (in Russian).
16. Koval'-chuk A.K., Semenov S.E., Karginov L.A., Kulakov D.B., Vereikin A.A., Yarots V.V. Choice of kinematic structure and study of dynamics of tree-like actuator of robot-tripod. Inzhenernyi vestnik MGTU im. N.E. Baumana = Engineering Herald of the Bauman MSTU, 2013, no. 10. Available at: http: //engbul. bmstu. ru/doc/654 481. html, accessed 20. 10. 2014. (in Russian).
17. Koval'-chuk A.K., Vereykin A.A., Karginov L.A., Tsepkova A.A. Design and calculation of power of drives of medical robot taking into account the dynamics of its actuator. 12 Mezhdunarodnaya nauchno-prakticheskaya konferentsiya & quot-Tekhnika i tekhnologiya: novye perspektivy razvitiya& quot- (25. 02. 2014): mater. [Proc. of the 12 International Scientific-Practical Conference & quot-Engineering and technology: new prospects of development& quot- (25. 02. 2014)]. Moscow, & quot-Sputnik+"- Publ., 2014, pp. 28−41. (in Russian).
18. Koval'-chuk A.K., Romashko A.M., Vereykin A.A., Karginov L.A., Malyakina E.A., Stativa V.A. Kinematic and dynamic analysis of actuator of six-legged walking robot. 14 Mezhdunarodnaya nauchno-prakticheskaya konferentsiya & quot-Sovremennoe sostoyanie estestvennykh i tekhnicheskikh nauk& quot- (14. 03. 2014): mater. [Proc. of the 14 International Scientific-Practical Conference & quot-The current state of natural sciences and engineering& quot- (14. 03. 2014)]. Moscow, & quot-Sputnik+"- Publ., 2014, pp. 60−79. (in Russian).
19. Koval'-chuk A.K. Choice of a kinematic structure and reseach of a & quot-Robot-dog"- tree-like actuator. Izvestiia vysshikh uchebnykh zavedenii. Mashinostroenie = Proceedings of Higher Educational Institutions. Machine Building, 2011, no. 8, pp. 65−73. (in Russian).
20. Koval'-chuk A.K. Design and investigation of the dynamics of tree-like actuator of six-legged walking robot. Estestvennye i tekhnicheskie nauki, 2014, no. 4 (72), pp. 104−110. (in Russian).
21. Kovalchuk A.K. Designing drives of a medical robot actuator. Life Science Journal, 2014, no. 11 (11s), pp. 337−340.
22. Vereikin A.A., Koval'-chuk A.K., Kulakov D.B., Semenov S.E., Karginov L.A., Kulakov B. B, Yarots V.V. Synthesis kinematic scheme of the actuating mechanism of the exoskeleton. 13 Mezhdunarodnaya nauchno-prakticheskaya konferentsiya & quot-Aktual'-nye voprosy nauki& quot- (25. 04. 2014): materialy [Proc. of the 13 International Scientific-Practical Conference & quot-Actual Problems of Science& quot- (25. 04. 2014)]. Moscow, & quot-Sputnik+"- Publ., 2014, pp. 68−76. (in Russian).
23. Vereikin A.A., Koval'-chuk A.K., Kulakov D.B., Semenov S.E. Analisys and Choice of the Exoskeleton'-s Actuator Kinematic Structure. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the Bauman MSTU, 2014, no. 7, pp. 72−93. DOI: 10. 7463/0714. 717 676 (in Russian).
24. Vereykin A.A., Koval'-chuk A.K., Kulakov D.B., Semenov S.E., Karginov L.A., Kulakov B. B, Yarots V.V. Dynamics of actuator of exoskeleton. 13 Mezhdunarodnaya nauchno-prakticheskaya konferentsiya & quot-Tekhnika i tekhnologiya: novye perspektivy razvitiya& quot- (26. 05. 2014): mater. [Proc. of the 13 International Scientific-Practical Conference & quot-Engineering and technology: new prospects of development& quot- (26. 05. 2014)]. Moscow, & quot-Sputnik+"- Publ., 2014, pp. 5−16. (in Russian).
25. Koval'-chuk A.K. Method for the synthesis of tree-like actuators of walking robots. Estestvennye i tekhnicheskie nauki, 2014, no. 3 (71), pp. 127−130. (in Russian).
26. Denavit J., Hartenberg R.S. A kinematic notation for Lower-Pair Mechanisms Based on Matrices. J. Appl. Mech, 1955, vol. 77, pp. 215−221.
27. Koval'-chuk A.K., Kulakov B.B., Kulakov D.B., Semenov S.E., Yarots V.V. Osnovy teorii ispolnitel'-nykh mekhanizmov shagayushchikh robotov [Fundamentals of the theory of actuators of walking robots]. Moscow, Rudomino Publ., 2010. 170 p. (in Russian).
28. Koval'-chuk A.K. Method of mathematical description of the kinematics and dynamics of tree-like actuators of walking robots. Estestvennye i tekhnicheskie nauki, 2014, no. 5 (73), pp. 87−90. (in Russian).
29. Koval'-chuk A.K., Karginov L.A., Kulakov B.B., Kulakov D.B., Semenov S.E., Yarots V.V., Vereykin A.A. Modelirovanie drevovidnykh ispolnitel'-nykh mekhanizmov shagayushchikh robotov s uchetom vneshnikh nalozhennykh svyazey [Simulation of treelike actuators of walking robots, taking into account external imposed external relations]. Certificate of state registration of computer program, no. 2 014 612 547 RF, 2014. (in Russian).
30. Koval'-chuk A.K., Karginov L.A., Kulakov B.B., Kulakov D.B., Semenov S.E., Yarots V.V. Programma modelirovaniya drevovidnykh ispolnitel'-nykh mekhanizmov shagayushchikh robotov [Simulation program of tree-like actuators of walking robots]. Certificate of state registration of computer program, no. 2 012 610 398 RF, 2012. (in Russian).
31. Alekseev V.E., Talanov V.A. Grafy i algoritmy. Struktury dannykh. Modeli vychisleniy [Graphs and algorithms. Data structures. Models of computation]. Moscow, Binom Publ., 2006. 319 p. (in Russian).

ПоказатьСвернуть
Заполнить форму текущей работой