student from 01.01.2018 to 01.01.2021
Russian Federation
Russian Federation
employee
doctoral candidate
Belgorod, Belgorod, Russian Federation
student
Belgorod, Russian Federation
UDK 004.94 Компьютерное моделирование
The article presents a model, tools, methodology and results of identifying special provisions when planning the movement of a parallel manipulator such as the Gough-Stewart platform. Special positions are understood as positions having geometric uncertainty, singularity. Singular positions are estimated by analyzing the Jacobian matrix using its singular value decomposition. The calculation of the Jacobi matrix and its singular values is performed using special libraries of the high-level Python programming language. The presence of special positions in the workspace is analyzed using the dexterity index, which characterizes the ratios of the singular values of the singular value decomposition of the Jacobian matrix. A series of computational experiments are carried out to calculate the dexterity index, identify and analyze special positions in the working area. The results of a computational experiment for various design versions of the platform, which differ in the ratio of the radii of the movable platform and the fixed base, are presented. The analysis of the characteristic features of the formation of singularity zones in their working zones is carried out and design options are identified that have the maximum dexterity index of the working area. The presented approach is applicable in the framework of the development, modernization, optimization of structural elements of parallel structure manipulators.
parallel structure manipulators, special positions, dexterity index, singularity, motion planning, Gough-Stewart platform
Введение. Параллельные манипуляторы, имеющие несколько кинематических звеньев, параллельно передающих движение подвижной платформе относительно неподвижного основания, обладают рядом преимуществ по сравнению с механизмами с последовательной конфигурацией звеньев: простота элементов механизма и их сборки, возможность простого масштабирования схемы устройства, высокая жесткость конструкции при относительно малом весе и низкой металлоёмкости, способность манипулировать большими нагрузками при низкой энергоемкости, высокая точность перемещения и позиционирования, хорошие динамические характеристики, большие скорости и ускорения. Однако, конструктивная особенность в виде параллельных звеньев является также причиной серьезных недостатков. Нелинейная зависимость кинематики и динамики механизма в зависимости от точки рабочей области приводит к анизотропии и неоднородности динамических, упругих и скоростных свойств, а также возможности интерференции (пересечения) отдельных кинематических цепей. Это в свою очередь является причиной особых положений в виде геометрической неопределенности, в которых может произойти потеря устойчивости. Геометрическая неопределенность ставит под угрозу работу манипулятора, при прохождении таких зон существенно снижается производительность, а в худшем случае может произойти полный отказ при заклинивании механизма. Поэтому усилия ряда ученых направлены на разработку надежных алгоритмов планирования движения [1-7, 12-20].
Планирование движения параллельного манипулятора связано с решением прямой задачи кинематики, заключающейся в определении координат рабочего органа в зависимости от длин выдвижения штоков, и обратной задачи кинематики — определения длин выдвижения штоков по заданным координатам рабочего органа. Ряд исследователей применяют аналитический подход, при этом прямая задача кинематики не имеет четко формализованного решения и требует применения более сложных математических методов по сравнению с решением обратной задачи, решение которой в общем случае сводится к решению шести нелинейных уравнений для шестистепенных манипуляторов типа платформы Гью-Стюарта 6-6. Математические подходы, применяемые для решения прямой задачи кинематики, можно условно разделить на две группы: использование аппарата векторной алгебры и аналитической геометрии. При использовании аппарата векторной алгебры координаты рабочего органа выражаются через орты подвижной системы координат посредством векторов, связывающих основание и подвижную платформу, или через составляющие матрицы поворота. При использовании аналитической геометрии конструкция представлена в виде геометрических соотношений пространственной конфигурации. Решение задачи прямой кинематики для параллельного манипулятора считается сложным в вычислительном отношении и трудоемким процессом из-за связанных нелинейных уравнений, которые сложно решать аналитически. В данной работе предлагается применить искусственные нейронные сети (ИНС) для решения прямой кинематики параллельного манипулятора. ИНС известны своей способностью представлять нелинейные отношения между входными и выходными данным, они обладают большими способностями к обучению и способны учиться на примерах.
В настоящее время существуют две стратегии решения проблемы геометрической неопределенности манипулятора при прохождении зоны сингулярности: методы, основанные на безопасном прохождении зоны сингулярности, и алгоритмы планирования движения в обход зон сингулярности. В рамках первого варианта стратегии предложены методы: ограничения крутящего момента двигателей при прохождении зоны сингулярности и установлении критериев устойчивости в виде определения производных высших порядков кинематических величин или создания оптимальных динамических условий [6, 7]; использования активных масс посредством отключения основных приводов исполнительных механизмов и включения массовых приводов в близости зон сингулярности [8]; применения физических критериев для прохождения зон сингулярности [12]. Однако наиболее надежными по мнению ряда исследователей являются алгоритмы планирования движения с целью исключения и обхода зон сингулярности.
Виейра и др. в своей статье [5] представили следующую методику планирования движения параллельного манипулятора: для определения положений платформы, в которых возникает неопределенность, используется моделирование методом Монте-Карло, на основе этого моделирования вычисляются вероятности отказа для нескольких конфигураций рабочего пространства манипулятора и полученные результаты используются для обучения искусственной нейронной сети. Нейронная сеть применена в этой методике с целью преодоления вычислительной неэффективности метода Монте-Карло. Полученная модель оценивается с помощью итеративной стратегии, которая использует генетические операторы для вычисления оптимальных траекторий, избегая областей, на которые в значительной степени влияют неопределенности. В качестве примера применения и отработки предложенной методики используется манипулятор 3RRR.
Мерле [6] пишет о концепциях матрицы Якоби (Jacobian matrix), манипулируемости (manipulability) и числах обусловленности (condition number) как о показателях точности параллельных роботов с учетом оптимальности конструкции, заключающейся в исключении таких положений платформы, при которых могут возникать ошибки позиционирования платформы. Основным критерием вероятности возникновения неопределенности в текущей работе принят индекс управляемости (dexterity index). Шанг [7] производит сравнение двух параллельных манипуляторов по критерию точности, базируясь в том числе на индексе управляемости.
В данной работе применена идея Виейра, заключающаяся в использовании результатов моделирования в виде положений платформы и вероятностей возникновения в них неопределенностей, для последующего обучения нейронной сети с целью быстрой работы по планированию оптимальной траектории платформы с избеганием положений, которые могут привести к отказу манипулятора. Представлена методика, инструментарий и результаты первого этапа: моделирования особых положений манипулятора типа платформы Гью-Стюарта с использованием расчета индекса управляемости по методике, описанной в работе Шанга.
Основная часть. Целью работы является создание интеллектуальной системы (ИС) для планирования движения роботизированной платформы Гью-Стюарта, а первым этапом реализации поставленной цели – разработка модели платформы, с использованием которой можно выполнять исследование особых положений. Модель платформы должна осуществлять возможность решения прямой и обратной задачи кинематики, а также проведения статического анализ сил. Поскольку эта задача является достаточно трудоемкой, было принято решение использовать специализированную библиотеку Basic robotics [8], реализующую перечисленные функции на высокоуровневом языке программирования Python. Также были использованы другие необходимые библиотеки: Numpy [9], Pandas [10], Matplotlib [11].
Используя библиотеку Basic robotics [8], представляется возможным создать объект платформы по заданным параметрам, после чего вызывать функции для сгенерированного объекта платформы, например, изменять положение подвижного основания (рис. 1).
Рис. 1. Визуализация изменения положения верхней платформы:
а – начальное положение, б – промежуточное положение, в – конечное положение
Анализ сингулярных положений платформы, то есть положений, в которых происходит геометрическая неопределённость, выполняют обычно по матрице Якоби. Матрица Якоби связывает скорости обобщенных координат и рабочего органа. Содержание и структура матрицы Якоби зависит от конструктивных параметров робота и изменяется в каждом новом положении при осуществлении движения по заданной траектории. Сингулярные положения могут быть вызваны различными причинами, например бесконечным множеством решений в уравнениях обратной задачи кинематики или располагаться в окрестности границы рабочей области. В сингулярных конфигурациях ограниченные значения скорости, силы и вращающего момента рабочего органа могут соответствовать неограниченным значениям обобщенных скоростей и вращающих моментов. Нулевое значение определителя матрицы Якоби говорит о том, что соответствующая конфигурация является сингулярной.
Матрица Якоби для платформы Гью-Стюарта в общем виде выглядит следующим образом:
где 
– матрицы частных производных от неявной функции по
и от неявной функции по обобщенным координатам 
– абсолютные скорости центра выходного звена; 
– обобщенные скорости во вращательных шарнирах 
Неявная функция между обобщенными и абсолютными координатами имеет вид:
Для решения задачи необходимо учесть положения входных звеньев, описываемых углами 
а также значения координат 
положения центра платформы 
и углов Эйлера 
поворота платформы.
В общем случае может возникнуть три вида сингулярности: сингулярной является матрица B (первый тип), сингулярной является матрица A (второй тип), сингулярными являются обе матрицы A и B (третий тип). В любом из вариантов сингулярность является препятствием движению манипулятора и ее необходимо избегать при планировании движения манипулятора. Существуют различные способы анализа сингулярности матрицы Якоби. В данной работе для этой цели применяется сингулярное разложение, декомпозиция вещественной матрицы и приведение ее к каноническому виду.
Предлагается генерировать множество точек пространства, после чего для каждой точки вычислять и анализировать соответствующую этому положению матрицу Якоби при условии нахождения центра верхней платформы в сгенерированной точке.
От количества точек будет зависеть точность результатов: чем больше будет точек и чем ближе они будут друг к другу, тем большее количество возможных положений платформы будет проанализировано. Однако с ростом количества точек растет и трудоемкость вычислений. Была написана функция для генерации пространства точек в виде куба. Предполагается, что все точки данного пространства находятся в рабочей области. Такое допущение связано с тем, что рабочая область платформы Стюарта представляет собой сложное пространство и для упрощения анализа механизма робота используется не все пространство, а только некоторая кубическая область, каждая точка которой входит в рабочую зону. На рис. 2 представлены сформированные области с разным количеством точек.
Рис. 2. Сгенерированные области точек: а – 1000 точек, б – 27000 точек
Далее основная идея заключается в том, чтобы для каждой точки сгенерированного пространства рассчитать индекс управляемости. Существуют разные способы расчета данного индекса, однако в данной работе был использован следующий метод:
где 
и 
x – это минимальное и максимальное сингулярные значения матрицы Якоби.
Для получения матрицы Якоби была использована соответствующая встроенная функция библиотеки Basic robotics, а для расчета сингулярных значений матрицы Якоби – функция svd из библиотеки Numpy.
Для примера приведем указанные параметры для положений платформы, представленных на рис.1. Координаты центра платформы соответственно в начальном, промежуточном и конечном положениях – [0, 0, 0, 0, 0, 0], [-0.3, 0.3, 1, 0.2, 0, 0], [-0.5, 0.4, 1, 0.2, 0.15, -0.05].
Матрицы Якоби соответственно для начального, промежуточного и конечного положений:
Сингулярные значения матрицы Якоби соответственно для начального, промежуточного и конечного положений:
Индекс управляемости соответственно для начального, промежуточного и конечного положений: 0.254, 0.145, 0.093. Из представленных результатов явно прослеживается, что индекс управляемости снижается при движении от центра к периферии рабочей области, что соответствует традиционным представлениям.
В общем случае индекс управляемости может принимать значения в диапазоне [0, 1]. Чем более индекс управляемости приближен к нулю, тем сильнее возникает геометрическая неопределенность. Таким образом, нулевое значение индекса управляемости представляет собой сингулярность, а единичное значение изотропное положение.
Результаты. По представленной выше методике выполнен вычислительный эксперимент и на рис. 3 приведены его результаты: визуализация рассчитанных значений индексов управляемости в кубическом пространстве для 27000 точек. Зоны высокой сингулярности отмечены темно синим цветом и расположены симметрично относительно исследуемой рабочей зоны, что связано с геометрической симметричностью конструкции.
Рис. 3. Рассчитанные значения индекса управляемости для 27 000 точек
Для более детального понимания полученного результата была написана подпрограмма для возможности “рассечения” куба, результат работы которой представлен на рис. 4.
Рис. 4. Визуализация значений индексов управляемости для 27000 точек с разрезом
Результаты по полученным индексам управляемости для 27000 точек обобщены в таблице 1, где cw – это выраженное в процентах соотношение области точек, в которых платформа имеет допустимый индекс управляемости, т. е. выполняется условие di>cdi.
Таблица 1
Результаты вычисления индексов управляемости
|
dimin |
dimax |
diavg |
cw, % |
rb/rt |
|
0,19 |
6,4510-6 |
0,07 |
22 |
1,8 |
В качестве критического значения индекса управляемости принято значение 0,1 [12]. Исходя из полученных результатов, можно сделать вывод о том, что только 22% точек находятся в допустимом значении индекса управляемости для исследуемого конструктивного исполнения.
Как отмечалось ранее, сингулярные положения обусловлены конструктивными параметрами платформы, в данном случае важно обратить внимание на соотношение радиусов основания к подвижной платформе. Было сделано предположение, что от этого соотношения зависит объем пространства с допустимым индексом управляемости. Оставляя прежними все другие параметры платформы, а также рассчитывая аналогичным образом пространство с точками в количестве 27000, была проведена серия экспериментов, в которых изменялся радиус подвижной платформы. Результаты эксперимента представлены в таблице 2.
Таблица 2
Результаты эксперимента с изменением радиуса подвижной платформы
|
rb/rt |
0,75 |
0,80 |
0,85 |
0,90 |
0,95 |
1,00 |
1,05 |
1,10 |
1,15 |
1,20 |
1,25 |
|
cw, % |
49,70 |
53,44 |
56,73 |
58,99 |
60,60 |
61,19 |
59,83 |
58,57 |
55,80 |
53,07 |
50,13 |
|
diavg |
0,12 |
0,13 |
0,15 |
0,17 |
0,19 |
0,2 |
0,19 |
0,17 |
0,16 |
0,14 |
0,13 |
Для более наглядной демонстрации результатов серий экспериментов были построены графики, представленные на рис. 5.
Рис. 5. Результаты вычислительного эксперимента по исследованию конструктивных исполнений соотношений радиусов подвижного и неподвижного оснований платформы: а – график зависимости diavg от rb/rt, б – график зависимости cw от rb/rt
Анализируя полученные зависимости, можно сделать вывод о том, что при заданных геометрических параметрах платформы, наиболее удовлетворительными показателями с точки зрения среднего индекса управляемости и процентного соотношения количества точек, в которых платформа имеет допустимый индекс управляемости, обладает конструкция платформы с равными радиусами подвижной платформы и основания. Далее в работе будет использована именно такая конфигурация.
Кроме того, представляет интерес визуализация значений индексов управляемости для 27000 точек для случая равенства нижней и верхней платформ (рис. 6).
Рис. 6. Визуализация значений индексов управляемости для 27000 точек с разрезом для конструктивного исполнения при равенстве радиусов подвижной платформы и неподвижного основания
Выводы. Представленная модель манипулятора типа платформы Гью-Стюарта и методика определения особых положений посредством вычисления индекса управляемости показывает хорошие результаты и может использоваться при планировании движения параллельного манипулятора с целью обхода зон высокой сингулярности, что повысит надежность его функционирования. Методика позволяет прорабатывать конструктивные варианты параллельного манипулятора, оценивать влияние отдельных конструктивных элементов на особенности формирования рабочей зоны, выявлять проблемные зоны с высокой сингулярностью, выполнять оптимизацию конструкции по критерию максимального индекса управляемости рабочей зоны. Представленный подход применим в рамках разработки, модернизации, оптимизации конструктивных элементов манипуляторов параллельной структуры.
1. Ermilov G.S. Development of simulator management systems based on the Stuart platform [Razrabotka sistem upravleniya trenazhera na osnove platformy Styuarta]. Colloquium-Journal. 2020. No. 12-1 (64). Pp. 11-21. DOIhttps://doi.org/10.24411/2520-6990-2020-11786. (rus)
2. Lopatin A.A. Parsing of mechanisms with six degrees of freedom and practical application on the example of the Gew-Stewart platform [Razbor mekhanizmov s shest'yu stepenyami svobody i prakticheskogo primeneniya na primere platformy G'yu-Styuarta]. Modern problems of machine theory. 2020. No. 9. Pp. 33-36. DOIhttps://doi.org/10.26160/2307-342X-2020-9-33-36. (rus)
3. Greshnyakov P.I. Stewart platform pneumatic control system development and characteristic research [Sovershenstvovanie pnevmaticheskoj sistemy upravleniya platformy Styuarta]. Bulletin of the Samara State Aerospace University named after Academician S.P. Korolev (National Research University). 2016. Vol. 15. No. 1. Pp. 151-162. DOIhttps://doi.org/10.18287/2412-7329-2016-15-1-151-162 (rus)
4. Koltsov A.G., Blokhin D.A., Khabarov A.V. Influence of the kinematic characteristics of the Stuart platform on the accuracy of movement of the measuring mechanism [Vliyanie kinematicheskih harakteristik platformy Styuarta na tochnost' peremeshcheniya izmeritel'nogo mekhanizma]. Dynamics of systems, mechanisms and machines. 2014. No 2. Pp. 296-299. (rus)
5. Sopin P.K., Gaynullina Ya.N. Design booth-simulation based on the stewart-gough platform [Razrabotka konstrukcii stenda-simulyatora na osnove platformy Styuarta-Gaufa]. Fundamental and applied problems of engineering and technology. 2015. No. 5(313). Pp. 142-145. (rus)
6. Ider S.K. Inverse dynamics of parallel manipulators in the presence of drive singularities. Mech. Mach. Theory. 2005. Vol. 40. Pp. 33-44. DOI:https://doi.org/10.1016/j.mechmachtheory.2004.05.007.
7. Briot S., Arakelian V. Optimal force generation in parallel manipulators for passing through the singular positions. Int. J. Robot. Res. 2008. Vol. 27. Pp. 967-983. DOIhttps://doi.org/10.1177/0278364908094403.
8. Parsa S.S., Boudreau R., Carretero J.A. Reconfigurable mass parameter s to cross direct kinematic singularities in parallel manipulators. Mech. Mach. Theory. 2015. Vol. 85. Pp. 53-63. DOIhttps://doi.org/10.1016/j.mechmachtheory.2014.10.008
9. Briot S., Pagis G., Bouton N., Martinet P. Degeneracy conditions of the dynamic model of parallel robots. Multibody Syst. Dyn. 2006. Vol. 37. Pp. 371-412. DOIhttps://doi.org/10.1007/s11044-015-9480-9.
10. Vieira H.L., da Silva M.M, Beck A.T. Reliable motion planning for parallel manipulators. Mech. Mach. Theory. 2019. Vol. 140. Pp. 553-566. DOIhttps://doi.org/10.1016/j.mechmachtheory.2019.06.022.
11. Duyun I.A., Kornilov A.V., Duyun T.A. Evaluation of kinematic and dynamic parameters of the Hough-Stewart platform using a virtual prototype [Primenenie interfejsa Python-Adams dlya modelirovaniya raboty platformy G'yu-Styuarta]. Bulletin of BSTU named after V.G. Shukhov. 2022. Vol. 7. No. 10. Pp. 96-106. DOIhttps://doi.org/10.34031/2071-7318-2022-7-10-96-106. (rus)
12. Duyun I.A., Gorlov A.S., Duyun T.A. Co-simulation parallel manipulator movements using Adams-Matlab [Sovmestnoe modelirovanie dvizheniya parallel'nogo manipulyatora c ispol'zovaniem Adams-Matlab]. Bulletin of the Belgorod state technological university named after V.G. Shukhova. 2022. Vol. 7. No 11. Pp. 108-119. DOI:https://doi.org/10.34031/2071-7318-2022-7-11-108-119. (rus)
13. Duyun I.A., Gorlov A.S., Duyun T.A. Modeling the trajectory of a robotic platform using a virtual prototype and the Python-Adams interface [Modelirovanie traektorii dvizheniya robotizirovannoj platformy s ispol'zovaniem virtual'nogo prototipa i interfejsa Python-Adams]. STIN. 2022. No. 12. Pp. 34-37. (rus)
14. Merlet J. P. Jacobian, Manipulability, Condition Number, and Accuracy of Parallel Robots. Journal of Mechanical Design. 2006. Vol. 128. Pp. 199-206. DOI:https://doi.org/10.1115/1.2121740.
15. Shang W., Cong S. Dexterity and Adaptive Control of Planar Parallel Manipulators With and Without Redundant Actuation. Journal of Computational and Nonlinear Dynamics 2015. Vol. 10. Pp. 12-23. DOI:https://doi.org/10.1115/1.4027581.
16. Vieira H.L., de Carvalho Fontes J.V., Beck A.T., da Silva M.M. Robust Critical Inverse Condition Number for a 3RRR Robot Using Failure Maps, Springer International Publishing/ Pp. 285-294. DOIhttps://doi.org/10.1007/978-3-319-67567-1_27.
17. Library Basic robotics. URL: https://github.com/64-B1T/basic_robotics (date of treatment: 09.01.2023).
18. Library Numpy. URL: https://numpy.org/ (date of treatment: 09.01.2023).
19. Library Pandas. URL: https://pandas.pydata.org/ (date of treatment: 09.01.2023).
20. Library Matplotlib. URL: https://matplotlib.org/ (date of treatment: 09.01.2023).
