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

УДК 621.01

DOI: 10.25206/2588-0373-2018-2-3-71-76

ОПРЕДЕЛЕНИЕ СФЕРИЧЕСКОЙ КРИВОЙ, ЗАДАЮЩЕЙ УГОЛ СЕРВИСА РУКИ АНДРОИДНОГО РОБОТА МЕТОДОМ СИНТЕЗА МАЛЫХ ДВИЖЕНИЙ

Ф. Н. Притыкин, В. И. Небритов

Омский государственный технический университет, Россия, 644050, г. Омск, пр. Мира, 11

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

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

I ■

л

О

1Я 1> N1

ОИ О О Е н Т х

>О 2 А

■ К > О

1 о

О

< К ОО

Введение

При моделировании процессов перемещения объектов манипулирования с использованием роботов необходимо на начальных этапах выяснять двигательные возможности и собственные свойства исполнительных механизмов [1—4]. Двигательные возможности механизмов манипуляторов так же необходимо учитывать при интеллектуальном управлении движением робототехнических систем [5 — 8]. Одной из задач на этом этапе является определение угла сервиса. Решение данной задачи на виртуальном уровне возможно с помощью моделирования малых движений механизмов манипуляторов по вектору скоростей выходного звена с обеспечением заданного фиксированного положения центра выходного звена. Заметим, что обход запретных зон при моделировании движений возможен только при наличии двигательной избыточности. Эту избыточность определяют разностью р = п —г, где параметр п соответственно задает размерность вектора обобщенных скоростей О механизма манипулятора (или число обобщенных координат), г — размерность вектора скоростей выходного звена V В работе [9] приведен способ определения угла сервиса на основе использования векторного уравнения замкнутости эквивалентного механизма. Однако при сложной структуре кинематических цепей механизмов манипуляторов, заданных различных предельных значений обобщенных координат и двигательной избыточности, равной двум и более, использование указанного способа затруднено. В работе [10] предложен способ определения телесного угла на основе исследования линейной зависимости векторов V и О [11, 12]. Однако синтез движений оси схватоносителя в различных направлениях при неподвижном центре захвата в указанной работе не исследовался, а определялись лишь отдельные положения оси схватоно-сителя, соответствующие допустимым значениям вектора обобщенных скоростей для заданной конфигурации. Изображения указанного угла серви-

са, полученного в работе [10], не в полной мере графически отражают двигательную избыточность и манипулятивность манипулятора в различных точках конфигурационного пространства с учётом заданных предельных значений обобщённых координат. В работе [13] предложен обобщенный способ определения угла сервиса с обеспечением синтеза движений оси схватоносителя по часовой и против часовой стрелки при неподвижном центре захвата для плоских механизмов манипуляторов. В данных работах исследовались плоские семизвенные механизмы манипуляторов, имеющие вращательные и поступательные кинематические пары. Однако в этих работах предложенный метод определения угла сервиса может быть использован только для плоских механизмов манипуляторов, имеющих различное число обобщенных координат и различную структуру кинематических цепей. В настоящей работе продолжены геометрические исследования, связанные с определением телесного угла сервиса для механизма руки андроидного робота с использованием синтеза движений оси схватоносителя по различным направлениям при неподвижном центре захвата.

Теория

Сущность предлагаемого метода основана на нахождении такой конфигурации, при которой угол сервиса в заданном секторе развёртки сферы единичного радиуса увеличивается [10]. Данная конфигурация находится реализацией дискретных значений векторов О из области допустимых значений. Область допустимых значений вектора О указанного механизма для различных конфигураций исследована в работах [14, 15]. Проведенные в работе исследования позволяют геометрическим методом оценить такие собственные свойства манипуляторов, как манипулятивность и маневренность для заданных предельных значений обобщённых координат.

Пусть задан пространственный семизвенный механизм руки андроидного робота ЛК600Б (рис. 1а).

а) б)

Рис. 1. Определение параметров, задающих положение точки L на развёртке сферы: а) кинематическая схема механизма руки андроидного робота и положение углов а и р в инерциальной системе координат; б) положение точки L на развёртке сферы Fig. 1. Determination of the parameters defining the position of the point L on the development of the sphere: a) the kinematic scheme of the arm's mechanism of the android robot and the position of the angles а and p in the inertial coordinate system; b) position of the point L on the development of the sphere

Рис. 2. Определение углов а и р в инерциальной системе координат O0, задающих положение точки L на развёртке сферы Fig. 2. The determination of the angles а and p in the inertial coordinate system O0 defining the position of the point L on the development of the sphere

На рисунке показано положение системы координат O x y z , связанной с характерной точкой кисти

n nJ n n' 1 L

руки, движение которой осуществляется по заданной траектории.

Точка O12 = On на рисунке определяет центр выходного звена, который не должен отклоняться от заданной точки траектории на заданную величину 5, называемую точностью позиционирования. При проведении исследований значение параметра принято 5 = 1 мм.

Реализация векторов Q, из области допустимых значений [10] обеспечивает новое положение кинематической цепи и оси схватоносителя. Под реализацией понимается новое положение кинематической цепи, определяемое значениями обобщённых координат, которые получают путём сложения обобщённых координат старой конфигурации и рассчитанных значений компонент вектора Q по уравнениям (при принятом допущении Aq е qi):

qi =qi + qi' qi = qi + q i,

где q — компоненты вектора Q.

(1)

С целью графического представления телесного угла и для пространственного шестизвенного механизма руки (рис. 1а) изобразим условную развертку поверхности сферы Ф единичного радиуса. Между точками этой условной развертки и точками сферы Ф существует взаимно однозначное соответствие. Любая точка Б поверхности сферы Ф радиуса ОпБ задает единственную точку Ь на этой развертке, определяемую полярными координатами ар и рр (рис. 1б). Отрезок ОпБ при этом равен единице (рис. 1). Развёртка, как видно из рис. 1б, имеет форму круга. В этом случае совокупность близлежащих крайних точек Ь будет определять условную развертку некоторого участка поверхности сферы Ф с центром в точке Ор.

Множеству отдельных положений оси От,п схватоносителя (рис. 1) при реализации значений вектора О и 8 = 1 мм будут соответствовать положения точки Ь на развёртке сферы единичного радиуса, принадлежащие различным секторам развёртки. Зададим сектор развёртки углами в рах и в рр111. В данном секторе угол в изменяется в интервале в < в р < в рах.

Точка Ь на развёртке сферы определяется параметрами ар и рр. Параметр ар на этой развертке задает расстояние от центра круга Ор до точки Ь развёртки и определяется углом, образованным между начальным положением продольной оси схватоносителя zn и осью схватоносителя , полученной реализацией вектора О (рис 1а). Единичные направляющие векторы этих осей определяют векторы ^ и Ш2 (рис. 2). Угол задает направление при нахождении точки Ь на условной развертке с помощью полярных координат (рис. 1б). Этот угол определяется углом р между проекцией вектора оси схватоносителя на плоскость П±О z и вектором

пп

задающим направление оси хп (рис. 1а и рис. 2) Ш1 е П. Вектор ^ находится проекцией вектора W5 на плоскость П. При этом вектор Ш5 = Ш3 - Ш'. Для нахождения точки О'!? (рис. 2) используется вектор Ш3П = Ш' + Ш7. При этом вектор является проекцией вектора Ш5 по направлению вектора Шд.

Положение вектора определяется точками Сп и О'П и векторами W8 и где Ш8 = Ш4 - Ш3', Ш = Ш — W. Точка Оявляется проекцией точ-

а) б)

Рис. 3. Проекции угла а и положение точки L на развёртке сферы при синтезе движения в одном из секторов развёртки для конфигурации q = (25°, 20°, -65°, 25°, -65°) при: а) -75° < q . < 75°; б) -120° < q. < 120° Fig. 3. Projections of the angle and the position of the point L on the development of the sphere during the synthesis of motion in one of the sectors for the configuration q. = (25°, 20°, -65°, 25°, -65°) at : a) -75° < q . < 75°; b) -120° < q . < 120°

о о s s О К s И ûl

О А

О Н

го z

S J= О И О О E Н T I >0 о >

û Т

I!

Рис. 4. Схема алгоритма вычисления положения крайних точек L, определяющих сферическую кривую lp на развёртке сферы Fig. 4. The scheme of the algorithm for calculating the position of the extreme points L defining the spherical curve lp on the development of the sphere

о т "О

» о

ки С>'п на плоскость П. Оси и хп принадлежат системе координат, связанной с выходным звеном первоначально заданной конфигурации (рис. 1а). Косинусы углов а и в вычисляют скалярными произведениями векторов (рис. 2):

W2 • W2

cos а = 2 2

W21 • W21 '

6

_ W • W6

cos ß = !-1 . , .

WilW1

(2)

ся допущение O12 « O{2. Точки O^ и O'2 определяют реализацией мгновенного состояния механизма руки. Начальное положение конфигурации руки андроидного робота, определяющей положение точки 012, зададим значениями обобщенных координат q.(25°, 20°, —65°, 25°, —65°). На каждом шаге расчетов определяют максимальное значение угла amax из значений a и a', где a и a' — значение угла соответственно на k-ой и k+1 итерации.

(max a, a')

^ a"

(3)

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

После определения конфигурации k +1, которая удовлетворяет условиям:

ak+1 > aш

ß г > ß > ß m

(4)

находят следующую конфигурацию к + 2 с сохранением списка обобщенных координат д. = д..

а)

б)

в)

Рис. 5. Положение сферической кривой l на развёртке сферы при -120° < q. < 120°для конфигураций: а) q. = (25°, 20°, 10°, 60°, 10°); б) q"= (25°, 20°, 35°, 35°, 35°); в) q. = (25°, 20°, 60°, 10°, 60°); Fig. 5. The position of the spherical curve l on the development of the sphere at -120 ° < q < 120 ° for configurations: a) q. = (25°, 20°, 10°, 60°, 10°p b) q. = (25°, 20°, 35°, 35°, 35°); c) q. = (25°, 20°, 60°, 10°, 60°)

Hp Hp

поверхности сферы. На рис. 3а указанный сектор задают соответственно значения ß™m = 90° и ßm<lx = 120°. При этом за исходно заданную конфигурацию принимается конфигурация k +1, для которой выполняется условие (4).

г, min

Задание различных значений углов pp и вртах позволяет последовательно выполнять синтез движений оси схватоносителя при нахождении точки О внутри сферы радиуса 1 мм в заданных секторах развёртки.

Алгоритм вычисления положений точек L на развёртке сферы представлен на рис. 4.

Обозначения, принятые на схеме: 1 — ввод зна-

ттлтттттт л- irr rr rr ^ st max гт min „ max ^ та nmin nmax

чений qt(qv q2 ..., q5^ qi , qi , ap = 0, v, вр , вp , 5 . = 1 мм, q. = q.; 2 — q. = q. ; 3 — вычисление

mm ' ^ ip ^r ^ i ^ ip'

матриц M0k, определяющих положение звеньев и значений векторов — W4; 4 — расчёт матриц M0k, векторов W/, W2, W3 и W5 — W9 и реализацией вектора Q; 5 — вычисление погрешностей линеаризации 5p; 6 — 5p < 5; 7 — qmin < qmax; 8 — вычисление а и ß; 9 — определение amax (3); 10 — а > amax, в™* > в > emax; 11 — qip = qi; 12 — изменение параметров ki, задающих значение вектора Q при наличии двигательной избыточности [10, 14] k. = k. + 1; 13 — k. > kmax [14]; 14 — расчёт вектора Q; 15 — вывод координат последней точки L, принадлежащей сферической кривой l и сектору ß , где вm'n < вР > epmax (рис. 3); 16 —p в^ > 360°; П -задание следующего сектора развёртки в^ и вmax' 18 — конец.

p

Результаты экспериментов

На рис. 3а представлены результаты исследований, связанных с определением проекций угла а, образованного между крайними положениями оси схватоносителя, на трех плоскостях проекций и на условной развертке сферы Ф. При данном синтезе движений оси схватоносителя точка Ь развёртки смещается в области заданного сектора развёртки. Предельные значения обобщённых координат заданы в интервале —75° < q. < 75°. На рис. 3б данный угол соответственно построен для интервала —120° < qí < 120°. Углы Ц, и2 и из на (рис. 3а, б) определяют проекции угла а. На рис. 3 отображена линейчатая поверхность От с помощью её горизонтальной, фронтальной и профильной проекций Б, Вт2 и Бт3. Линейчатая

поверхность Бт образована движением оси схватоносителя (заданной отрезком 012Ои (рис. 1)) при неподвижном центре захвата 012. Данная точка при синтезе движений не выходит за пределы окружности 5 = 1 мм. Линейчатая поверхность в полной мере характеризует маневренность манипулятора. Изображения проекций Бт1, Вт2 и От3 области Бт дают лишь некоторое геометрическое представление об указанной маневренности манипулятора в заданном секторе развёртки сферы при заданных значениях q¡max и q¡min.

На рис. 3а, б отражено движение оси схватоно-сителя при заданной точности позиционирования 5 = 1 мм и модуле вектора |У(Ух, V, V) — 1 мм Л. Вектор V при этом параллелен оси х неподвижной системы координат. На рис. 3а, б показаны положения точки Ь в заданном секторе. Крайние точки данной траектории задают положение сферической кривой 1р. Площадь замкнутого контура заданного этой кривой определяет телесный угол сервиса руки андроидного робота.

На рис. 5 показано положение сферической кривой для трех конфигураций, заданных значениями обобщённых координат q. (25°, 20°, 10°, 60°, 10°), qi (25°, 20°, 35°, 35°, 35°) и qi(25°, 20°, 60°, 10°, 60°). Предельные значения обобщённых координат приняты -120° < qi < 120°.

Основные результаты и выводы

Результаты проведенных исследований позволяют сделать следующие выводы:

1. Форма и положение сферической кривой 1р в значительной степени зависят от положения кинематической цепи манипулятора.

2. При вычислении вектора О, следует использовать максимальные значения параметров ктах [14]. При этом значительно сокращается время расчетов при нахождении сферической кривой.

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

нения ориентации ВЗ может быть использован в интеллектуальных системах управления роботов с целью смены положения конфигурации по отношению к заданным запретным зонам.

Список источников

1. Lopatin P. Investigation of a target reachability by a manipulator in an unknown environment // IEEE International Conference on Mechatronics and Automation, 2016 August 7 — 10. Harbin, China. 2016. P. 37-42.

2. Егоров А. С., Лопатин П. К. Использование алгоритма полиномиальной аппроксимации в задаче управления манипулятором в среде с неизвестными препятствиями // Мехатро-ника, автоматизация, управление. 2013. № 3. С. 24-29.

3. Карташев В. А., Богуславский А. А., Карташев В. В. [и др.]. Задачи управления манипуляционным роботом для обеспечения безопасности перемещений // Мехатроника, автоматизация, управление. 2015. № 1. С. 24-28.

4. Крахмалев О. Н., Петрешин Д. И., Крахмалев Г. Н. Математические модели систем управления для калибровки ориентации инструмента промышленных роботов // Мехатрони-ка, автоматизация, управление. 2017. Т. 18, № 10. С. 664-668. DOI: 10.17587/mau.18.664-669.

5. Лохин В. М., Манько С. В., Александрова Р. И. [и др.]. Механизмы интеллектуальных обратных связей, обработки знаний и самообучения в системах управления автономными роботами и мультиагентными робототехническими группировками // Мехатроника, автоматизация, управление. 2015. Т. 10, № 8. С. 545-555. DOI: 10.17587/mau.16.545-555.

6. Лохин В. М., Манько С. В., Романов М. П. Повышение адаптивных свойств автономных роботов на базе интеллектуальных технологий // Экстремальная робототехника. 2015. № 1 (1). С. 59-67.

7. Юсупов Р. М., Тимофеев А. В. Интеллектуализация процессов управления и навигации робототехнических систем // Экстремальная робототехника. 2014. Т. 1, № 1. С. 16-21.

8. Mitchell S. Fuzzy logic decision making for autonomous robotic applications // IEEE 6th International Conference on Awareness Science and Technology (iCAST). Paris, 2014. P. 1-6. DOI: 10.1109/ICAwST.2014.6981843.

9. Лебедев П. А. Аналитический метод определения коэффициента сервиса манипулятора // Проблемы машиностроения и надежности машин. 1991. № 5. С. 93-98.

10. Притыкин Ф. Н. Графическое представление телесного угла и окружающего пространства руки при реализации мгновенных состояний манипуляторов // Проблемы машиностроения и надежности машин. 2002. № 3. С. 93-101.

11. Кобринский А. А., Кобринский А. Е. Манипуляцион-ные системы роботов. М.: Наука, 1985. 344 с.

12. Корендясев А. И., Саламандра Б. П., Тывес П. Н. Определение числа степеней свободы исполнительного органа промышленного робота // Машиноведение. 1985. № 6. С. 44 — 53.

13. Pritykin F., Gordeev O. Defining a service angle for planar mechanisms of manipulators based on the instantaneous states analysis // MEACS2015. IOP Conf. Series Materials Science and Engineering. 2016. № 124. P. 1-7. DOI:10.1088/1757-899X/124/1/012025.

14. Притыкин Ф. Н., Небритов В. И. Исследование размеров и формы области в многомерном пространстве обобщённых скоростей, задающей допустимые мгновенные состояния механизма андроидного робота // Омский научный вестник. 2016. № 5 (149). С. 29-34.

15. Pritykin F. N., Nebritov V. I. Studying tolerance range of generalized velocities vector under android motion synthesis // Dynamics of Systems, Mechanisms and Machines (Dynamics). 2016, Nov. 15-17, 2016. 7819064. D0I:10.1109/Dynamics.2016.7819064.

ПРИТЫКИН Федор Николаевич, доктор технических наук, доцент (Россия), профессор кафедры «Инженерная геометрия и САПР». SPIN-код: 7628-8023 ORCID: 0000-0001-8081-6840 AuthorID (SCOPUS): 6507269253 Адрес для переписки: [email protected] НЕБРИТОВ Валерий Иванович, аспирант кафедры «Инженерная геометрия и САПР». SPIN-код: 8196-1026 AuthorID (SCOPUS): 57193406556 Адрес для переписки: [email protected]

Для цитирования

Притыкин Ф. Н., Небритов В. И. Определение сферической кривой, задающей угол сервиса руки андроидного робота методом синтеза малых движений // Омский научный вестник. Сер. Авиационно-ракетное и энергетическое машиностроение. 2018. Т. 2, № 3. С. 71-76. DOI: 10.25206/2588-03732018-2-3-71-76.

Статья поступила в редакцию 27.06.2018 г. © Ф. Н. Притыкин, В. И. Небритов