Выпуск #3/2021
Е. ГЕБЕЛЬ, Л. ГАВРИЛИНА, В. ГЛАЗУНОВ, А. ЛАСТОЧКИН, А. РОМАНОВ, А. ДУХОВ
К АНАЛИЗУ СИНГУЛЯРНЫХ ЗОН МЕХАНИЗМОВ ПАРАЛЛЕЛЬНОЙ СТРУКТУРЫ С ЛИНЕЙНЫМИ ДВИГАТЕЛЯМИ
К АНАЛИЗУ СИНГУЛЯРНЫХ ЗОН МЕХАНИЗМОВ ПАРАЛЛЕЛЬНОЙ СТРУКТУРЫ С ЛИНЕЙНЫМИ ДВИГАТЕЛЯМИ
Просмотры: 1546
DOI: 10.22184/2499-9407.2021.24.3.92.98
Рассмотрены различные особые положения (сингулярности) механизмов параллельной структуры с линейными двигателями, а также особое (сингулярное) положение, при котором был смещен центр выходного звена по сравнению с исходной конфигурацией. Проведено сравнение кинематических и взаимных винтов, соответствующих различным особым положениям (сингулярностям).
Рассмотрены различные особые положения (сингулярности) механизмов параллельной структуры с линейными двигателями, а также особое (сингулярное) положение, при котором был смещен центр выходного звена по сравнению с исходной конфигурацией. Проведено сравнение кинематических и взаимных винтов, соответствующих различным особым положениям (сингулярностям).
Теги: linear motors parallel mechanism reciprocal twist singularities special positions взаимный кинематический винт линейные двигатели механизм параллельной структуры особые положения сингулярности
К АНАЛИЗУ СИНГУЛЯРНЫХ ЗОН МЕХАНИЗМОВ ПАРАЛЛЕЛЬНОЙ СТРУКТУРЫ С ЛИНЕЙНЫМИ ДВИГАТЕЛЯМИ
Елена ГЕБЕЛЬ, Любовь ГАВРИЛИНА, Виктор ГЛАЗУНОВ,
Алексей ЛАСТОЧКИН, Андрей РОМАНОВ, Андрей ДУХОВ
Рассмотрены различные особые положения (сингулярности) механизмов параллельной структуры с линейными двигателями, а также особое (сингулярное) положение, при котором был смещен центр выходного звена по сравнению с исходной конфигурацией. Проведено сравнение кинематических и взаимных винтов, соответствующих различным особым положениям (сингулярностям).
Введение
Для решения задач развития технических систем, применяемых в различных областях промышленности, добывающих отраслях, медицине и робототехнике [1‒4], эффективным средством являются механизмы параллельной структуры [5‒8], при исследовании которых целесообразно использовать винтовое исчисление [9].
Следует отметить, что данное направление исследований активно развивается. В частности, речь идет о механизмах с тросовыми приводами [10, 11], о применении этих устройств в медицинских системах и антропоморфных роботах [12, 13]. Углубленно исследуются различные свойства этих объектов и, в частности, возможности их балансировки и снижения энергозатрат [14, 15], способы устранения возможных особых положений (сингулярностей) [16, 17].
При анализе указанных технических систем важно учитывать конкретные условия, для которых предназначен данный объект [18‒22]. Кроме того, необходимо принимать во внимание свойства технических объектов, в частности возможность наличия особых положений [7, 8, 23, 24]. В этих положениях резко снижается нагрузочная способность, может присутствовать неуправляемая подвижность. Наличие таких положений ухудшает функциональные возможности рассматриваемых устройств.
В качестве примера приведем перспективный технологический робот METROM, имеющий пять степеней свободы и пять кинематических цепей, каждая из которых включает линейный двигатель [25, 26]. Этот робот (рис. 1) имеет достаточно большую зону обслуживания в силу того, что его особые положения располагаются вне рабочего объема.
Вместе с тем, для большинства роботов параллельной структуры приходится ограничивать рабочую зону конфигурациями, далекими от особых, это сокращает рабочий объем.
Существует подход, позволяющий избежать этого недостатка – он заключается в переключении приводной кинематической пары при подходе к особым положениям, в частности, такие работы выполнены для трипода – механизма с тремя степенями свободы [27]. Однако, важно знать такие возможные движения робота, которые наиболее быстро выводят его из особого положения, а также такие движения, которые переводят в соседние особые положения. Дело в том, что дополнительный привод должен включаться кратковременно.
В связи с изложенным, в данной статье ставится задача определить кинематические винты – движения робота, которые соответствуют неуправляемой подвижности в особых положениях, а также наиболее быстрому выходу из особых положений. Эта задача рассматривается на относительно простом примере: когда все оси линейных двигателей пересекают одну прямую линию. В дальнейшем этот анализ будет распространен на более сложные случаи.
Определение винта-градиента, выводящего механизм из особого положения
Перед рассмотрением вопросов вывода механизма из особых положений укажем на некоторые известные определения [9]. Кинематический винт – это дуальный вектор, включающий в себя вектор угловой скорости и вектор линейной скорости твердого тела.
Силовой винт – это дуальный вектор, включающий вектор силы и вектор момента, действующие на твердое тело. Механизм параллельной структуры – это механизм, включающий в себя основание и выходное звено, соединенное несколькими кинематическими цепями. Особое положение (сингулярность) механизма параллельной структуры – это такое положение, при котором либо уменьшается число степеней свободы, либо имеет место неуправляемая подвижность, не связанная с движениями в приводах.
Плюккеровы координаты силового и кинематического винтов – это шесть координат, три из которых представляют вектор силы или угловой скорости, а три других представляют координаты вектора момента или линейной скорости.
Рассмотрим механизм параллельной структуры с линейными двигателями, который находится в состоянии сингулярности, поскольку все силовые винты, передаваемые со стороны кинематических цепей на выходное звено, пересекают одну прямую линию νν (рис. 2).
Отметим, что силовые винты нулевого параметра — это скользящие векторы. Скользящий вектор, пересекающий оси всех силовых винтов, взаимен всем указанным винтам. Точки А1–А6 и В1–В6 являются центрами сферических шарниров, расположенных соответственно на выходном звене (подвижной платформе) и основании (базе). Векторы Е1–Е6 ‒ единичные винты, направленные вдоль осей линейных двигателей, представляющих собой силовые винты. Линия νν пересекает оси всех силовых винтов (отсюда следует, что все винты Е1–Е6 – линейно зависимые).
Нужно найти винт-градиент с плюккеровыми координатами (ξ, η, ζ, ξ0, η0, ζ0), наиболее быстро выводящий из состояния сингулярности, а также соотнести его с кинематическим винтом, взаимным пяти независимым силовым винтам. Кроме того, покажем, что при расположении всех точек пересечения силовых винтов с линией νν на основании любое движение не будет выводить механизм из сингулярности.
Не нарушая общности, предположим, что линяя νν совпадает с осью Ох (рис. 3).
Точки А1–А6 и В1–В6 имеют следующие координаты:
Используем разработанный ранее подход [8, 23, 24] к поиску кинематического винта-градиента, который ведет из особого положения по «наикратчайшему пути». Основу данного подхода представляет определение кинематического винта рабочего органа, а также связанных с ним плюккеровых координат векторов единичной длины силовых винтов, лежащих на осях двигателей линейного перемещения.
Рассмотрение проведем для конкретного механизма. В результате решения задачи о положениях для заданных координат центров сферических шарниров на выходном звене и основании механизма с линейными двигателями (см. рис. 3) получены обобщенные координаты (м):
Q = (2,449 2,449 3,742 3,000 4,359 4,243).
Матрица плюккеровых координат ортов силовых винтов формируется согласно соотношению:
где первые три строки – направляющие косинусы единичного винта, а три оставшиеся – моментная часть винта, XBi, YBi, ZBi – координаты точки Bi; XAi, YAi, ZAi‒ координаты точки Ai; Li – длина отрезка AiBi, где i = 1, …, 6.
Матрица плюккеровых координат для заданного положения механизма:
. (1)
Ранг матрицы (1) равен пяти, а определитель ‒ нулю. Запишем частные производные по плюккеровым координатам кинематического винта.
Для каждой строки матрицы, описывающей частные производные от плюккеровых координат ортов силовых винтов по соответствующей плюккеровой координате кинематического винта выходного звена, имеют место соотношения, зависящие от положения механизма [8, 23, 24].
Для нахождения координат винта-градиента далее в матрице плюккеровых координат последовательно нужно заменить столбцы на соответствующие столбцы из матрицы частных производных. В итоге, последовательно суммируя определители полученных матриц, можно с точностью до множителя найти координаты винта-градиента.
Далее для иллюстрации результатов представим полученный винт-градиент в виде двух векторов. Обозначим векторную часть винта через V (ξ, η, ζ), а моментную – через M (ξ0, η0, ζ0). Модуль вектора V определяется как корень квадратный из суммы квадратов его координат. Принимаем модуль вектора, равный двум, изменяя при этом ξ, η, ζ. Параметр винта p определяется как скалярное произведение вектора на момент, деленное на квадрат модуля вектора.
Радиус-вектор точки В, принадлежащей оси винта, равен:
.
Радиус-вектор точки А, соответствующей концу векторной части:
ρА = ρВ + V.
Радиус-вектор точки С, соответствующей концу моментной части:
ρС = ρВ + ρV.
Матрицы частных производных по координате x имеют вид:
Аналогично вычислены значения частных производных для остальных плюккеровых координат. В соответствии с найденными параметрами получены значения плюккеровых координат винта-градиента:
(ξ, η, ζ, ξ0, η0, ζ0) = (0, 0, 0, 0, 0, 0).
Полученный результат объясняется тем, что при любых движениях выходного звена сохраняются условия сингулярности, поскольку все оси силовых винтов всегда пересекают ось Oх.
Определение винта-градиента для смещенного положения выходного
звена
Для другого механизма параллельной структуры с линейными двигателями (рис. 4) точки В1–В6 смещены вдоль оси Oz на расстояние: ‒1 м.
Соответствующие значения координат точек Ai и Bi выходного звена и основания исследуемого механизма запишутся следующим образом:
При этом матрица плюккеровых координат остается неизменной. Интерес представляет взаимное соответствие кинематического винта, взаимного ортам силовых винтов, и кинематического винта-градиента.
В результате решения обратной задачи о положениях получены следующие значения обобщенных координат механизма (м):
Q = (4,899 3,674 4,989 4,500 5,812 8,485).
Матрица плюккеровых координат для заданного положения центров сферических шарниров выходного звена и основания механизма с линейными двигателями не изменилась. Следовательно, заданное положение исследуемого механизма также является сингулярным.
Запишем частные производные по плюккеровым координатам кинематического винта. Частная производная по ξ:
.
Частные производные по остальным плюккеровым координатам вычисляются аналогичным образом.
Винт-градиент, наиболее быстро выводящий из состояния сингулярности, имеет плюккеровы координаты:
(ξ, η, ζ, ξ0, η0, ζ0) = (1,603 –1,520 0,081 0,000 –0,241 0,294).
В рассматриваемом положении механизма параллельной структуры кинематический винт-градиент, наиболее быстро выводящий из особого положения, не совпадает с единственным кинематическим винтом, взаимным ортам осей пяти независимых силовых винтов, имеющим координаты (1, 0, 0, 0, 0, 0).
Зная винт-градиент, находим плюккеровы координаты пяти кинематических винтов, переводящих в бесконечно близкие соседние особые положения:
Для рассматриваемого случая координаты точек B, A, C, характеризующих винт-градиент:
.
Приведем другие виды изображений исследуемого механизма (рис. 5).
При сканировании положения выходного звена с постоянной его ориентацией можно получить «срез» сингулярной поверхности [24]. Рассмотрим одно из положений, соответствующее этой поверхности, при котором центр выходного звена сместился по координатам (0,5; 0,2; 0,133).
Координаты точек выходного звена примут вид:
.
Координаты единичного вектора силового винта для каждой кинематической цепи найдем аналогично описанному выше случаю, шесть кинематических цепей соответственно дают матрицу PL:
.
Ранг полученной матрицы 5 – новое положение механизма также является сингулярным (особым). Винт, взаимный ортам осей силовых винтов, не изменился.
Найдем винт-градиент с плюккеровыми координатами (ξ, η, ζ, ξ0, η0, ζ0).
Частная производная по ξ:
.
Остальные частные производные аналогичны рассмотренному случаю.
В результате получаем винт-градиент
(1,403 –1,27 0,00253 0 –0,163 0,292).
Координаты точек B, A, C, характеризующих винт-градиент:
.
Полученное решение представим с учетом расположения звеньев механизма и винта-градиента (рис. 6).
Таким образом, винт-градиент изменился.
Построение сингулярных «поверхностей» при изменении ориентации выходного звена
Далее рассмотрим сингулярные «поверхности», получаемые при вращениях выходного звена. Матрица, описывающая ориентацию выходного звена, получается путем перемножения матриц, соответствующих вращениям вокруг каждой из осей.
Вращение вокруг оси Ох:
Матрицы, описывающие вращения вокруг осей Оy и Оz, имеют аналогичный вид. В результате перемножения этих матриц получим матрицу M, описывающую поворот вокруг трех осей. Смещенное положение точек выходного звена будет таковым:
.
Зададим максимальные углы поворота от 0 до 60° для поворотов относительно осей Оx и Оy, для оси Оz максимальный поворот примем от 0 до 180°.
При перемещении из начального положения в конечные будем последовательно увеличивать углы на фиксированные значения, равные одной двадцатой от максимального отклонения. В каждом из полученных положений вычисляем ранг матрицы, сформированной из силовых винтов, действующих на выходное звено со стороны кинематических цепей. При приближении к особому положению определитель матрицы плюккеровых координат стремится к нулю, здесь принято в качестве признака для особого положения, что определитель матрицы силовых винтов по модулю должен быть меньше 0,005. В результате можно получить «поверхности» в координатах α, β, γ, в которых механизм находится в особых положениях.
Заключение
В данной статье представлен анализ особых положений механизмов параллельной структуры с линейными двигателями. Данные механизмы имеют широкое поле применения – это двигательные системы, технологические роботы, измерительные и обучающие устройства. Особые положения существенно снижают функциональные возможности, ограничивая рабочую зону. Для того чтобы исключить этот недостаток, можно использовать дополнительные двигатели, однако требуется знать, какие движения переводят в соседние особые положения и какие выводят из них.
На примере относительно простого случая, когда оси всех линейных двигателей пересекают одну прямую линию, установлено, что кинематический винт, описывающий неуправляемое движение, не соответствует кинематическому винту, описывающему движение, наиболее быстро выводящему из особого положения. При этом при изменении точек крепления приводов имеет место изменение кинематического винта, соответствующего наиболее быстрому выводу из особого положения.
Данное исследование в дальнейшем будет распространено на более сложные случаи особых положений, соответствующих конкретным механизмам параллельной структуры с линейными двигателями, применяемыми в различных областях техники, в частности в станкостроении.
Литература
1. Ганиев Р. Ф. Проблемы механики машин и технологий. Перспективы развития Института машиноведения им. А. А. Благонравова РАН. Ч. 1 // Проблемы машиностроения и надежности машин. 2010. № 1. С. 3–20.
2. Ганиев Р. Ф. Проблемы механики машин и технологий. Перспективы развития Института машиноведения им. А. А. Благонравова РАН. Ч. 2 // Проблемы машиностроения и надежности машин. 2010. № 3. С. 3–17.
3. Каляев И. А., Капустян С. Г. Групповое управление роботами: проблемы, решения // Известия высших учебных заведений. Машиностроение. 2011. № 12. С. 7–12.
4. Лопота А. В., Спасский Б. А. Robotics. A View of the Future. // Мехатроника, автоматизация, управление. 2017. Т. 18. № 9. С. 612–615.
5. Davidson J. K., Hunt K. H. Robot and Screw Theory: Applications of Kinematics and Statics of Robotics // Oxford University Press, 2004. 467 р.
6. Воробьев Е. И. Осуществление заданного относительного движения двух твердых тел двуруким роботом // Механика твердого тела. 2018. № 2. С. 122–128.
7. Bohigas O., Zlatanov D., Ros L., Manubens M., Porta J.M. Numerical computation of manipulator singularities // Proceedings of the IEEE International Conference on Robotics and Automation, 2012, article number 6225083, pp. 1351–1358.
8. Алешин А. К., Глазунов В. А., Рашоян Г. В., Оффер Шаи Анализ кинематических винтов, определяющих топологию сингулярных зон роботов параллельной структуры // Проблемы машиностроения и надежности машин. 2016. № 4. С. 3–8.
9. Диментберг Ф. М. Теория винтов и ее приложения. М.: Наука, 1978. 336 с.
10. Sugahara Y. et al. A Suspended Cable-Driven Parallel Robot for Human-Cooperative Object Transportation // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. PP. 109–117.
11. Tempel P., Alfeld M., Wijk V. Design and Analysis of Cable-Driven Parallel Robot CaRISA: A Cable Robot for Inspecting and Scanning Artwork // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. PP. 136–144.
12. Russo M., Ceccarelli M. Dynamics of a Humanoid Robot with Parallel Architectures // Advances in Mechanism and Machine Science. Proceedings of the 15th IFToMM World Congress on Mechanism and Machine Science. Ed. T. Uhl. Springer. 2019. P. 1799–1808.
13. Li R., Fan X., Li X., Bai S., Zhang J. Kinematic Design of a 2-SPS/PU&R 4-DOF Hybrid Robot for Ankle Rehabilitation // Advances in Mechanism and Machine Science. Proceedings of the 15th IFToMM World Congress on Mechanism and Machine Science. Ed. T. Uhl. Springer. 2019. P. 1849–1858.
14. Balderas Hill R., Briot S., Chriette A., Martinet P. Minimizing the Energy Consumption of a Delta Robot by Exploiting the Natural Dynamics // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. P. 213–224.
15. Nguyen V. L., Kuo C. H., Lin C. Y. Gravity Compensation of Delta Parallel Robot Using a Gear-Spring Mechanism // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. P. 280–287.
16. Geng J., Arakelian V. Balancing of Planar 5R Symmetrical Parallel Manipulators Taking into Account the Varying Payload // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. P. 372–379.
17. Harada T., Kunishige Y. Singularity Free Mode Changes of a Redundantly Driven Two Limbs Six-Dof Parallel Robot // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. P. 405–413.
18. Balchanowski J., Szrek J., Wudarczyk S. Analysis of constraint equations of the parallel mechanisms with 3 DoF in singular configurations // Advances in Mechanism and Machine Science. Proceedings of the 15th IFToMM World Congress on Mechanism and Machine Science. Ed. T. Uhl. Springer. 2019. P. 607–616.
19. Fomin A., Antonov A., Glazunov V. Forward Kinematic Analysis of a Rotary Hexapod // CISM International Centre for Mechanical Sciences 2021. G. Venture et al. (Eds.): ROMANSY 2020, CISM 601, 2021. PP. 486–494.
20. Велиев Е. И., Ганиев Р. Ф., Глазунов В. А., Скворцов С. А., Чернецов Р.А. Разработка и исследование механизмов с постоянной точкой ввода инструмента в рабочую область, предназначенных для хирургических операций и исследования свойств плазмы // Проблемы машиностроения и надежности машин. 2020. № 6. С. 3–15.
21. Хейло С. В., Ларюшкин П. А. Определение рабочей зоны манипуляторов параллельной структуры // Справочник. Инженерный журнал. 2013. № 2. С. 27–31.
22. Пащенко В. Н., Романов А. В., Артемьев А. В., Меньшова Е. В., Логинов Н. А. Решение прямой задачи о положении шестистепенного манипулятора параллельной структуры на базе кривошипно-шатунного механизма // Электронные информационные системы. 2017. № 4(15). С. 91–101.
23. Гебель Е. С., Глазунов В. А. Определение взаимного соответствия кинематического винта выходного звена и винта-градиента в особом положении механизма параллельной структуры // Известия высших учебных заведений. Машиностроение. 2017. № 7 (688). С. 3–9.
24. Гебель Е. С., Гаврилина Л. В., Глазунов В. А., Демидов С. М., Швец П. А. Построение зоны сингулярности для одного вида механизмов параллельной структуры // Проблемы машиностроения и автоматизации. 2020. № 3. С. 4–10.
25. Патент Германии № WO 03/037562 Al. Device for the displacement and/or positioning of an object in five axes / МПК В23Q 1/54, опубл. 08.05.2003.
26. Чунихин А. Ю., Глазунов В. А. Разработка механизмов параллельной структуры с пятью степенями свободы, предназначенных для технологических роботов // Проблемы машиностроения и надежности машин. 2017. № 4. С. 3–11.
27. Arakelian V., Briot S., Glazunov V. Improvement of functional performance of spatial parallel manipulators using mechanisms of variable structure // Proceedings of the Twelfth World Congress in Mechanism and Machine Science IFToMM, Besancon, France. 2007. V. 5. P. 159–164.
ГЕБЕЛЬ Елена Сергеевна –
кандидат технических наук, доцент, заведующая кафедрой «Автоматизация и робототехника» Омского государственного технического университета (ОмГТУ)
ГАВРИЛИНА Любовь Васильевна –
научный сотрудник Федерального государственного бюджетного учреждения науки Института машиноведения им. А. А. Благонравова РАН
ГЛАЗУНОВ Виктор Аркадьевич –
доктор технических наук, профессор, директор Федерального государственного бюджетного учреждения науки Института машиноведения им. А. А. Благонравова РАН
ЛАСТОЧКИН Алексей Борисович –
кандидат технических наук, старший научный сотрудник Федерального государственного бюджетного учреждения науки Института машиноведения им. А. А. Благонравова РАН
РОМАНОВ Андрей Александрович –
инженер Федерального государственного бюджетного учреждения науки Института машиноведения им. А.А. Благонравова РАН
ДУХОВ Андрей Валерьевич –
начальник группы АО НИИМА «Прогресс»
Елена ГЕБЕЛЬ, Любовь ГАВРИЛИНА, Виктор ГЛАЗУНОВ,
Алексей ЛАСТОЧКИН, Андрей РОМАНОВ, Андрей ДУХОВ
Рассмотрены различные особые положения (сингулярности) механизмов параллельной структуры с линейными двигателями, а также особое (сингулярное) положение, при котором был смещен центр выходного звена по сравнению с исходной конфигурацией. Проведено сравнение кинематических и взаимных винтов, соответствующих различным особым положениям (сингулярностям).
Введение
Для решения задач развития технических систем, применяемых в различных областях промышленности, добывающих отраслях, медицине и робототехнике [1‒4], эффективным средством являются механизмы параллельной структуры [5‒8], при исследовании которых целесообразно использовать винтовое исчисление [9].
Следует отметить, что данное направление исследований активно развивается. В частности, речь идет о механизмах с тросовыми приводами [10, 11], о применении этих устройств в медицинских системах и антропоморфных роботах [12, 13]. Углубленно исследуются различные свойства этих объектов и, в частности, возможности их балансировки и снижения энергозатрат [14, 15], способы устранения возможных особых положений (сингулярностей) [16, 17].
При анализе указанных технических систем важно учитывать конкретные условия, для которых предназначен данный объект [18‒22]. Кроме того, необходимо принимать во внимание свойства технических объектов, в частности возможность наличия особых положений [7, 8, 23, 24]. В этих положениях резко снижается нагрузочная способность, может присутствовать неуправляемая подвижность. Наличие таких положений ухудшает функциональные возможности рассматриваемых устройств.
В качестве примера приведем перспективный технологический робот METROM, имеющий пять степеней свободы и пять кинематических цепей, каждая из которых включает линейный двигатель [25, 26]. Этот робот (рис. 1) имеет достаточно большую зону обслуживания в силу того, что его особые положения располагаются вне рабочего объема.
Вместе с тем, для большинства роботов параллельной структуры приходится ограничивать рабочую зону конфигурациями, далекими от особых, это сокращает рабочий объем.
Существует подход, позволяющий избежать этого недостатка – он заключается в переключении приводной кинематической пары при подходе к особым положениям, в частности, такие работы выполнены для трипода – механизма с тремя степенями свободы [27]. Однако, важно знать такие возможные движения робота, которые наиболее быстро выводят его из особого положения, а также такие движения, которые переводят в соседние особые положения. Дело в том, что дополнительный привод должен включаться кратковременно.
В связи с изложенным, в данной статье ставится задача определить кинематические винты – движения робота, которые соответствуют неуправляемой подвижности в особых положениях, а также наиболее быстрому выходу из особых положений. Эта задача рассматривается на относительно простом примере: когда все оси линейных двигателей пересекают одну прямую линию. В дальнейшем этот анализ будет распространен на более сложные случаи.
Определение винта-градиента, выводящего механизм из особого положения
Перед рассмотрением вопросов вывода механизма из особых положений укажем на некоторые известные определения [9]. Кинематический винт – это дуальный вектор, включающий в себя вектор угловой скорости и вектор линейной скорости твердого тела.
Силовой винт – это дуальный вектор, включающий вектор силы и вектор момента, действующие на твердое тело. Механизм параллельной структуры – это механизм, включающий в себя основание и выходное звено, соединенное несколькими кинематическими цепями. Особое положение (сингулярность) механизма параллельной структуры – это такое положение, при котором либо уменьшается число степеней свободы, либо имеет место неуправляемая подвижность, не связанная с движениями в приводах.
Плюккеровы координаты силового и кинематического винтов – это шесть координат, три из которых представляют вектор силы или угловой скорости, а три других представляют координаты вектора момента или линейной скорости.
Рассмотрим механизм параллельной структуры с линейными двигателями, который находится в состоянии сингулярности, поскольку все силовые винты, передаваемые со стороны кинематических цепей на выходное звено, пересекают одну прямую линию νν (рис. 2).
Отметим, что силовые винты нулевого параметра — это скользящие векторы. Скользящий вектор, пересекающий оси всех силовых винтов, взаимен всем указанным винтам. Точки А1–А6 и В1–В6 являются центрами сферических шарниров, расположенных соответственно на выходном звене (подвижной платформе) и основании (базе). Векторы Е1–Е6 ‒ единичные винты, направленные вдоль осей линейных двигателей, представляющих собой силовые винты. Линия νν пересекает оси всех силовых винтов (отсюда следует, что все винты Е1–Е6 – линейно зависимые).
Нужно найти винт-градиент с плюккеровыми координатами (ξ, η, ζ, ξ0, η0, ζ0), наиболее быстро выводящий из состояния сингулярности, а также соотнести его с кинематическим винтом, взаимным пяти независимым силовым винтам. Кроме того, покажем, что при расположении всех точек пересечения силовых винтов с линией νν на основании любое движение не будет выводить механизм из сингулярности.
Не нарушая общности, предположим, что линяя νν совпадает с осью Ох (рис. 3).
Точки А1–А6 и В1–В6 имеют следующие координаты:
Используем разработанный ранее подход [8, 23, 24] к поиску кинематического винта-градиента, который ведет из особого положения по «наикратчайшему пути». Основу данного подхода представляет определение кинематического винта рабочего органа, а также связанных с ним плюккеровых координат векторов единичной длины силовых винтов, лежащих на осях двигателей линейного перемещения.
Рассмотрение проведем для конкретного механизма. В результате решения задачи о положениях для заданных координат центров сферических шарниров на выходном звене и основании механизма с линейными двигателями (см. рис. 3) получены обобщенные координаты (м):
Q = (2,449 2,449 3,742 3,000 4,359 4,243).
Матрица плюккеровых координат ортов силовых винтов формируется согласно соотношению:
где первые три строки – направляющие косинусы единичного винта, а три оставшиеся – моментная часть винта, XBi, YBi, ZBi – координаты точки Bi; XAi, YAi, ZAi‒ координаты точки Ai; Li – длина отрезка AiBi, где i = 1, …, 6.
Матрица плюккеровых координат для заданного положения механизма:
. (1)
Ранг матрицы (1) равен пяти, а определитель ‒ нулю. Запишем частные производные по плюккеровым координатам кинематического винта.
Для каждой строки матрицы, описывающей частные производные от плюккеровых координат ортов силовых винтов по соответствующей плюккеровой координате кинематического винта выходного звена, имеют место соотношения, зависящие от положения механизма [8, 23, 24].
Для нахождения координат винта-градиента далее в матрице плюккеровых координат последовательно нужно заменить столбцы на соответствующие столбцы из матрицы частных производных. В итоге, последовательно суммируя определители полученных матриц, можно с точностью до множителя найти координаты винта-градиента.
Далее для иллюстрации результатов представим полученный винт-градиент в виде двух векторов. Обозначим векторную часть винта через V (ξ, η, ζ), а моментную – через M (ξ0, η0, ζ0). Модуль вектора V определяется как корень квадратный из суммы квадратов его координат. Принимаем модуль вектора, равный двум, изменяя при этом ξ, η, ζ. Параметр винта p определяется как скалярное произведение вектора на момент, деленное на квадрат модуля вектора.
Радиус-вектор точки В, принадлежащей оси винта, равен:
.
Радиус-вектор точки А, соответствующей концу векторной части:
ρА = ρВ + V.
Радиус-вектор точки С, соответствующей концу моментной части:
ρС = ρВ + ρV.
Матрицы частных производных по координате x имеют вид:
Аналогично вычислены значения частных производных для остальных плюккеровых координат. В соответствии с найденными параметрами получены значения плюккеровых координат винта-градиента:
(ξ, η, ζ, ξ0, η0, ζ0) = (0, 0, 0, 0, 0, 0).
Полученный результат объясняется тем, что при любых движениях выходного звена сохраняются условия сингулярности, поскольку все оси силовых винтов всегда пересекают ось Oх.
Определение винта-градиента для смещенного положения выходного
звена
Для другого механизма параллельной структуры с линейными двигателями (рис. 4) точки В1–В6 смещены вдоль оси Oz на расстояние: ‒1 м.
Соответствующие значения координат точек Ai и Bi выходного звена и основания исследуемого механизма запишутся следующим образом:
При этом матрица плюккеровых координат остается неизменной. Интерес представляет взаимное соответствие кинематического винта, взаимного ортам силовых винтов, и кинематического винта-градиента.
В результате решения обратной задачи о положениях получены следующие значения обобщенных координат механизма (м):
Q = (4,899 3,674 4,989 4,500 5,812 8,485).
Матрица плюккеровых координат для заданного положения центров сферических шарниров выходного звена и основания механизма с линейными двигателями не изменилась. Следовательно, заданное положение исследуемого механизма также является сингулярным.
Запишем частные производные по плюккеровым координатам кинематического винта. Частная производная по ξ:
.
Частные производные по остальным плюккеровым координатам вычисляются аналогичным образом.
Винт-градиент, наиболее быстро выводящий из состояния сингулярности, имеет плюккеровы координаты:
(ξ, η, ζ, ξ0, η0, ζ0) = (1,603 –1,520 0,081 0,000 –0,241 0,294).
В рассматриваемом положении механизма параллельной структуры кинематический винт-градиент, наиболее быстро выводящий из особого положения, не совпадает с единственным кинематическим винтом, взаимным ортам осей пяти независимых силовых винтов, имеющим координаты (1, 0, 0, 0, 0, 0).
Зная винт-градиент, находим плюккеровы координаты пяти кинематических винтов, переводящих в бесконечно близкие соседние особые положения:
- (1,520 1,603 0,000 0,000 0,000 0,000);
- (–0,081 0,000 1,603 0,000 0,000 0,000);
- (0,000 0,081 1,520 0,000 0,000 0,000);
- (0,000 0,000 0,000 0,000 0,294 0,241);
- (0,000 0,000 0,000 0,241 0,000 0,000).
Для рассматриваемого случая координаты точек B, A, C, характеризующих винт-градиент:
.
Приведем другие виды изображений исследуемого механизма (рис. 5).
При сканировании положения выходного звена с постоянной его ориентацией можно получить «срез» сингулярной поверхности [24]. Рассмотрим одно из положений, соответствующее этой поверхности, при котором центр выходного звена сместился по координатам (0,5; 0,2; 0,133).
Координаты точек выходного звена примут вид:
.
Координаты единичного вектора силового винта для каждой кинематической цепи найдем аналогично описанному выше случаю, шесть кинематических цепей соответственно дают матрицу PL:
.
Ранг полученной матрицы 5 – новое положение механизма также является сингулярным (особым). Винт, взаимный ортам осей силовых винтов, не изменился.
Найдем винт-градиент с плюккеровыми координатами (ξ, η, ζ, ξ0, η0, ζ0).
Частная производная по ξ:
.
Остальные частные производные аналогичны рассмотренному случаю.
В результате получаем винт-градиент
(1,403 –1,27 0,00253 0 –0,163 0,292).
Координаты точек B, A, C, характеризующих винт-градиент:
.
Полученное решение представим с учетом расположения звеньев механизма и винта-градиента (рис. 6).
Таким образом, винт-градиент изменился.
Построение сингулярных «поверхностей» при изменении ориентации выходного звена
Далее рассмотрим сингулярные «поверхности», получаемые при вращениях выходного звена. Матрица, описывающая ориентацию выходного звена, получается путем перемножения матриц, соответствующих вращениям вокруг каждой из осей.
Вращение вокруг оси Ох:
Матрицы, описывающие вращения вокруг осей Оy и Оz, имеют аналогичный вид. В результате перемножения этих матриц получим матрицу M, описывающую поворот вокруг трех осей. Смещенное положение точек выходного звена будет таковым:
.
Зададим максимальные углы поворота от 0 до 60° для поворотов относительно осей Оx и Оy, для оси Оz максимальный поворот примем от 0 до 180°.
При перемещении из начального положения в конечные будем последовательно увеличивать углы на фиксированные значения, равные одной двадцатой от максимального отклонения. В каждом из полученных положений вычисляем ранг матрицы, сформированной из силовых винтов, действующих на выходное звено со стороны кинематических цепей. При приближении к особому положению определитель матрицы плюккеровых координат стремится к нулю, здесь принято в качестве признака для особого положения, что определитель матрицы силовых винтов по модулю должен быть меньше 0,005. В результате можно получить «поверхности» в координатах α, β, γ, в которых механизм находится в особых положениях.
Заключение
В данной статье представлен анализ особых положений механизмов параллельной структуры с линейными двигателями. Данные механизмы имеют широкое поле применения – это двигательные системы, технологические роботы, измерительные и обучающие устройства. Особые положения существенно снижают функциональные возможности, ограничивая рабочую зону. Для того чтобы исключить этот недостаток, можно использовать дополнительные двигатели, однако требуется знать, какие движения переводят в соседние особые положения и какие выводят из них.
На примере относительно простого случая, когда оси всех линейных двигателей пересекают одну прямую линию, установлено, что кинематический винт, описывающий неуправляемое движение, не соответствует кинематическому винту, описывающему движение, наиболее быстро выводящему из особого положения. При этом при изменении точек крепления приводов имеет место изменение кинематического винта, соответствующего наиболее быстрому выводу из особого положения.
Данное исследование в дальнейшем будет распространено на более сложные случаи особых положений, соответствующих конкретным механизмам параллельной структуры с линейными двигателями, применяемыми в различных областях техники, в частности в станкостроении.
Литература
1. Ганиев Р. Ф. Проблемы механики машин и технологий. Перспективы развития Института машиноведения им. А. А. Благонравова РАН. Ч. 1 // Проблемы машиностроения и надежности машин. 2010. № 1. С. 3–20.
2. Ганиев Р. Ф. Проблемы механики машин и технологий. Перспективы развития Института машиноведения им. А. А. Благонравова РАН. Ч. 2 // Проблемы машиностроения и надежности машин. 2010. № 3. С. 3–17.
3. Каляев И. А., Капустян С. Г. Групповое управление роботами: проблемы, решения // Известия высших учебных заведений. Машиностроение. 2011. № 12. С. 7–12.
4. Лопота А. В., Спасский Б. А. Robotics. A View of the Future. // Мехатроника, автоматизация, управление. 2017. Т. 18. № 9. С. 612–615.
5. Davidson J. K., Hunt K. H. Robot and Screw Theory: Applications of Kinematics and Statics of Robotics // Oxford University Press, 2004. 467 р.
6. Воробьев Е. И. Осуществление заданного относительного движения двух твердых тел двуруким роботом // Механика твердого тела. 2018. № 2. С. 122–128.
7. Bohigas O., Zlatanov D., Ros L., Manubens M., Porta J.M. Numerical computation of manipulator singularities // Proceedings of the IEEE International Conference on Robotics and Automation, 2012, article number 6225083, pp. 1351–1358.
8. Алешин А. К., Глазунов В. А., Рашоян Г. В., Оффер Шаи Анализ кинематических винтов, определяющих топологию сингулярных зон роботов параллельной структуры // Проблемы машиностроения и надежности машин. 2016. № 4. С. 3–8.
9. Диментберг Ф. М. Теория винтов и ее приложения. М.: Наука, 1978. 336 с.
10. Sugahara Y. et al. A Suspended Cable-Driven Parallel Robot for Human-Cooperative Object Transportation // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. PP. 109–117.
11. Tempel P., Alfeld M., Wijk V. Design and Analysis of Cable-Driven Parallel Robot CaRISA: A Cable Robot for Inspecting and Scanning Artwork // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. PP. 136–144.
12. Russo M., Ceccarelli M. Dynamics of a Humanoid Robot with Parallel Architectures // Advances in Mechanism and Machine Science. Proceedings of the 15th IFToMM World Congress on Mechanism and Machine Science. Ed. T. Uhl. Springer. 2019. P. 1799–1808.
13. Li R., Fan X., Li X., Bai S., Zhang J. Kinematic Design of a 2-SPS/PU&R 4-DOF Hybrid Robot for Ankle Rehabilitation // Advances in Mechanism and Machine Science. Proceedings of the 15th IFToMM World Congress on Mechanism and Machine Science. Ed. T. Uhl. Springer. 2019. P. 1849–1858.
14. Balderas Hill R., Briot S., Chriette A., Martinet P. Minimizing the Energy Consumption of a Delta Robot by Exploiting the Natural Dynamics // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. P. 213–224.
15. Nguyen V. L., Kuo C. H., Lin C. Y. Gravity Compensation of Delta Parallel Robot Using a Gear-Spring Mechanism // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. P. 280–287.
16. Geng J., Arakelian V. Balancing of Planar 5R Symmetrical Parallel Manipulators Taking into Account the Varying Payload // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. P. 372–379.
17. Harada T., Kunishige Y. Singularity Free Mode Changes of a Redundantly Driven Two Limbs Six-Dof Parallel Robot // ROMANSY 23 – Robot Design, Dynamics and Control. 2021. P. 405–413.
18. Balchanowski J., Szrek J., Wudarczyk S. Analysis of constraint equations of the parallel mechanisms with 3 DoF in singular configurations // Advances in Mechanism and Machine Science. Proceedings of the 15th IFToMM World Congress on Mechanism and Machine Science. Ed. T. Uhl. Springer. 2019. P. 607–616.
19. Fomin A., Antonov A., Glazunov V. Forward Kinematic Analysis of a Rotary Hexapod // CISM International Centre for Mechanical Sciences 2021. G. Venture et al. (Eds.): ROMANSY 2020, CISM 601, 2021. PP. 486–494.
20. Велиев Е. И., Ганиев Р. Ф., Глазунов В. А., Скворцов С. А., Чернецов Р.А. Разработка и исследование механизмов с постоянной точкой ввода инструмента в рабочую область, предназначенных для хирургических операций и исследования свойств плазмы // Проблемы машиностроения и надежности машин. 2020. № 6. С. 3–15.
21. Хейло С. В., Ларюшкин П. А. Определение рабочей зоны манипуляторов параллельной структуры // Справочник. Инженерный журнал. 2013. № 2. С. 27–31.
22. Пащенко В. Н., Романов А. В., Артемьев А. В., Меньшова Е. В., Логинов Н. А. Решение прямой задачи о положении шестистепенного манипулятора параллельной структуры на базе кривошипно-шатунного механизма // Электронные информационные системы. 2017. № 4(15). С. 91–101.
23. Гебель Е. С., Глазунов В. А. Определение взаимного соответствия кинематического винта выходного звена и винта-градиента в особом положении механизма параллельной структуры // Известия высших учебных заведений. Машиностроение. 2017. № 7 (688). С. 3–9.
24. Гебель Е. С., Гаврилина Л. В., Глазунов В. А., Демидов С. М., Швец П. А. Построение зоны сингулярности для одного вида механизмов параллельной структуры // Проблемы машиностроения и автоматизации. 2020. № 3. С. 4–10.
25. Патент Германии № WO 03/037562 Al. Device for the displacement and/or positioning of an object in five axes / МПК В23Q 1/54, опубл. 08.05.2003.
26. Чунихин А. Ю., Глазунов В. А. Разработка механизмов параллельной структуры с пятью степенями свободы, предназначенных для технологических роботов // Проблемы машиностроения и надежности машин. 2017. № 4. С. 3–11.
27. Arakelian V., Briot S., Glazunov V. Improvement of functional performance of spatial parallel manipulators using mechanisms of variable structure // Proceedings of the Twelfth World Congress in Mechanism and Machine Science IFToMM, Besancon, France. 2007. V. 5. P. 159–164.
ГЕБЕЛЬ Елена Сергеевна –
кандидат технических наук, доцент, заведующая кафедрой «Автоматизация и робототехника» Омского государственного технического университета (ОмГТУ)
ГАВРИЛИНА Любовь Васильевна –
научный сотрудник Федерального государственного бюджетного учреждения науки Института машиноведения им. А. А. Благонравова РАН
ГЛАЗУНОВ Виктор Аркадьевич –
доктор технических наук, профессор, директор Федерального государственного бюджетного учреждения науки Института машиноведения им. А. А. Благонравова РАН
ЛАСТОЧКИН Алексей Борисович –
кандидат технических наук, старший научный сотрудник Федерального государственного бюджетного учреждения науки Института машиноведения им. А. А. Благонравова РАН
РОМАНОВ Андрей Александрович –
инженер Федерального государственного бюджетного учреждения науки Института машиноведения им. А.А. Благонравова РАН
ДУХОВ Андрей Валерьевич –
начальник группы АО НИИМА «Прогресс»
Отзывы читателей