PATH OPTIMIZATION AND SINGULARITY AVOIDANCE IN ROBOTIC MANIPULATORS USING HYBRID ACO-FLC TECHNIQUES
Abstract and keywords
Abstract:
The rapid advancements in robotics have necessitated the development of innovative techniques for efficient path optimization and singularity avoidance in robotic manipulators. This study explores the integration of ant colony optimization (ACO) and fuzzy logic control (FLC) to address these challenges in manipulators with six degrees of freedom (6-DOF). ACO is utilized for global path planning, ensuring optimal trajectories while avoiding obstacles and singularity zones. Simultaneously, FLC provides local adaptability, refining the path for smoothness and stability in dynamic environments. The hybrid ACO-FLC approach demonstrates significant improvements in trajectory efficiency, singularity avoidance, and computational performance. The ACO methodology demonstrates an average improvement of 18% in path efficiency compared to conventional methods. Furthermore, the integration of FLC enhances trajectory smoothness by 25%, ensuring accurate and stable motion. This hybrid ACO-FLC framework achieves an exceptional 98% singularity avoidance rate and reduces computational overhead by 15%, facilitating faster and more efficient performance in dynamic and complex robotic environments. Comparative analysis with other techniques highlights the superiority of the proposed method in achieving robust and adaptive control.

Keywords:
ant colony optimization, singularity avoidance, path optimization, robotic manipulators, intelligent control systems, degree of freedom
Text
Text (PDF): Read Download

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

Управление на основе нечеткой логики (НЛК) стало надежным подходом к учету неопределенностей и нелинейностей в системах управления роботами [4]. Имитируя процессы принятия решений человека, НЛК способствует адаптивному управлению в динамических условиях, что делает его особенно подходящим для задач в реальном времени. Однако автономные системы НЛК часто сталкиваются с трудностями в достижении глобально оптимальных решений в сложных многомерных пространствах [5].

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

Научная новизна данного исследования заключается в разработке и валидации гибридного метода, объединяющего алгоритм муравьиной колонии (ОМК) и управление на основе нечеткой логики (НЛК) для задачи одновременной оптимизации траектории и предотвращения сингулярностей в манипуляторах с шестью степенями свободы. В отличие от существующих работ, где методы ОМК и НЛК применяются раздельно или для мобильных роботов, предлагаемая методика реализована в пространстве обобщённых координат 6-DOF манипулятора типа RRRRRP с использованием детерминанта матрицы Якобиана в качестве штрафной функции в целевой задаче ОМК. Эффективность предложенного подхода подтверждена с помощью моделирования и численных экспериментов в Python, демонстрирующих значительное улучшение показателей избегания сингулярностей, гладкости траектории и вычислительной эффективности.

Исследование структурировано следующим образом: в первом разделе обзор литературы по оптимизации траекторий и предотвращению сингулярностей в роботизированных манипуляторах; в разделе “Материалы и методы” представлена методология оптимизации траектории роботизированного манипулятора, основанная на комбинированном применении алгоритма муравьиной колонии (ОМК) и нечеткой логики (НЛК) для обеспечения устойчивого и адаптивного движения с предотвращением сингулярностей. В “Основная часть” рассматриваются современные алгоритмы оптимизации и управления, направленные на эффективное планирование траекторий и предотвращение сингулярностей в роботизированных манипуляторах.

Обзор литературы. Интеграция управления на основе нечеткой логики (НЛК) с муравьиной оптимизацией (ОМК) вызвала значительный интерес в области робототехники благодаря её потенциалу в решении сложных задач, связанных с планированием траекторий и предотвращением сингулярностей. В данном разделе представлен подробный обзор значительных исследований в этой области.

Lizarraga и соавт. применили ОМК для планирования траекторий и интегрировали НЛК для улучшения локальных корректировок, что позволило достичь более эффективного избегания препятствий в мобильных роботах [6]. Аналогично Castillo предложили био-вдохновленный подход для оптимизации НЛК в робототехнических системах с использованием ОМК. Их работа продемонстрировала повышение эффективности автономных систем, особенно в условиях, требующих адаптивного управления [7].

Zhang и соавт усовершенствовали область исследований, разработав системы нечеткой логики типа-2, оптимизированные с помощью ОМК, что повысило их устойчивость в условиях неопределённости. Их исследования подчеркнули адаптивность НЛК в динамических условиях, делая их пригодными для приложений в реальном времени [8]. Yen и Cheng изучили применение НЛК с ОМК для поиска кратчайшего пути в мобильных роботах. Их исследования подчеркнули способность ОМК находить глобальные пути, в то время как НЛК обеспечивал корректировки в реальном времени, приводя к более гладким траекториям [9]. Этот двухуровневый подход эффективно предотвращал сингулярные конфигурации в роботизированных манипуляторах.

В контексте предотвращения сингулярностей Kureichik исследовал использование НЛК в сочетании с ОМК для оптимизации стратегий управления роботизированными манипуляторами. Они продемонстрировали, что использование НЛК снижает вычислительные затраты при сохранении высокой точности выполнения траектории [10]. Kumar сравнил производительность нечетких PID-контроллеров и НЛК, оптимизированных с помощью ОМК, придя к выводу, что последний обеспечивает лучшую гладкость и стабильность в нелинейных системах. Их работа подчеркнула преимущества сочетания био-вдохновленных алгоритмов с нечеткой логикой для сложных задач управления [11].

Недавние достижения включают интеграцию ОМК и НЛК в навигационные системы для автономных транспортных средств. Castillo и Melin представили метод динамической настройки параметров ОМК с использованием нечеткой логики. В этом исследовании подчёркивалась роль НЛК в повышении адаптивности ОМК и предотвращении преждевременной сходимости в процессе оптимизации [12]. Yen et al. сообщили о значительном улучшении в оптимизации траектории и избегании препятствий благодаря использованию преимуществ обоих подходов [13]. Аналогично, Castillo and Melin применили алгоритм муравьиной колонии (ОМК) для разработки нечетких регуляторов второго типа, добившись лучшей адаптивности и эффективности в динамических условиях [14].

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

Материалы и методы. Диаграмма показывает методологию оптимизации траектории робота с использованием ОМК и НЛК как показано на рисунке 1. Процесс начинается с определения проблемы, затем описания кинематики манипуляторов, анализа Якобиана и анализа сингулярностей. Следующим шагом является применение ОМК для глобальной оптимизации пути с избежанием сингулярностей. Затем реализуется избежание сингулярностей для предотвращения конфигураций. После этого НЛК уточняет траекторию, созданную с помощью ОМК, обеспечивая плавность и адаптивность. На следующем этапе интеграция ОМК с НЛК улучшает локальную адаптивность и плавность траектории. Наконец, методология проверяется с помощью симуляции и тестирования для анализа метрик производительности и демонстрации эффективности предложенной структуры.

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

 

Рис. 1. Методология оптимизации траектории движения робота

 

Описание и модель роботизированного манипулятора. В данной работе рассматривается манипулятор с шестью степенями свободы, включающий пять вращательных звеньев и одно поступательное (тип RRRRRP), что позволяет учитывать как ориентацию, так и положение рабочего органа. Манипулятор имеет пять звеньев (d1, a2, a3, a4 ), пять углов (θBase, θ2, θ3, θ4, θ5 ) и одну поступательную кинематическую пару d5 , как показано на рисунке 2. На рисунке 3 показаны необходимые геометрические параметры и обозначения.

 

Рис. 2. Конструкция роботизированного
манипулятора
[15]

Рис. 3. Кинематическая схема манипулятора

 

Денавита-Хартенберга Параметры (Д-Х). Метод Д-Х является стандартным подходом для описания кинематики манипуляторов. В таблице 1 подробно представлены параметры конструктора манипулятора, которые включали длину звеньев, смещения и углы звеньев[16].

 

Таблица 1

Параметры Д-Х роботизированного манипулятора [17]

i

αi-1 (deg)

ai-1 (cm)

θi (deg)

di (cm)

0-1     1

90

0

θbase

d1

1-2     2

0

a2

θ2

0

2-3     3

0

a3

θ3

0

3-4     4

0

a4

θ4

0

4-5     5

90

0

θ5

0

5-6     6

0

0

0

d5

 

Прямой кинематический анализ. При прямом кинематическом анализе углы звеньев используются для определения положения и ориентации рабочего органа путем их подстановки в однородную матрицу преобразования между звеньями i и i+1 .Tii+1=cosθi-sinθicosαisinθi sinαiai cosθisinθicosθicosαi-cosθisinαiaI sinθi0sinαicosαidi0001.                                       (1)

 

Выполнение композиции из системы координат i в базовую систему координат, мы умножаем шесть от 1 до 6 [18]:

 

T60  = T10  . T21  . T32  . T43  . T65   =nxoxaxpxnyoyaypynzozazpz0001 .                                              (2)                                

 

Общая матрица преобразования T60:  

 

T60=s6sbase+ c2345 c6 cbase sbasec6- c2345 s6 cbasec1 s2345Px c2345 sbase c6- cbase s6- cbasec6- c2345s6sbases1 s2345Pys2345 c6-s2345 s6- c2345pz0001,            (3)

     Px=  cbase*( a2*c2 + a3*c23+ a4*c234+d5*s2345) ,                       (4)

    Py = sbase*(a2*c2 + a3*c23 + a4*c234+d5*s2345) ,                       (5)

      Pz= d1 + a2*s2+ a3*s23  + a4*s234- d5* c2345 .                          (6)

 

где Px,Py,Pz – представляют позиция выходного звена.

Матрица Якобиана. Матрица Якобиана является ключевым инструментом для анализа влияния скоростей суставов на скорость конечного эффектора роботизированного манипулятора, а также для выявления сингулярных конфигураций. Сингулярность возникает, когда определитель матрицы Якобиана приближается к нулю, что приводит к потере управляемости манипулятора. Поэтому анализ сингулярностей является важным этапом в проектировании и управлении роботизированными системами [19].

Якобиан J  вычисляется как частная производная положения концевого эффектора относительно переменных звеньев [20]:

J=dPee ,                              (7)

где Pee=px,py,pzT позиция выходного звенав базисе робота. θi=θbase, θ2,θ3, θ4, θ5,d5  ‒ угол i‑го сустава.

Отношения между скоростями суставов и линейными Jv(3×n)  и угловыми Jw(3×n)  скоростями конечного эффектора, обозначаемыми соответственно, как, следующие:

V(6×1)=Jq. θ=Jv(3×n)Jw(3×n)(6×n)θ(n×1)  .   (8)

 

Где:  Jvi=zi-1×p-pi-1;    для i-го вращательного звенаzi-1;                            для i-го призматического звена

  Jwi= zi-1;     для i-го вращательного звена      0;      для i-го призматического звена .

 

Определение матрицы Якобиана. Матрица Якобиана связывает скорости изменения углов суставов (или линейных перемещений) со скоростью движения конечного эффектора в декартовом пространстве [21]. Для манипулятора с шестью степенями свободы матрица Якобиана имеет размерность 6×6 и может быть представлена в виде:

 

  x y z wx wy wz =  J11 J21 0001     J12 J22 J32 sbase-cbase0     J13 J23J33sbase-cbase0     J14J24J34sbase-cbase0     J15 J25 J35 sbase-cbase0     J16 J26 J36 000 . θ1 θ2θ3 θ4θ5 d5 .                               (9)

J11 = -sinbase ( a2 cos2 + a3 cos23 + a4 cos234 + d5 sin2345)

J12=-cosbase ( a2 sin2 +a3 sin23 + a4 sin234 - d5 cos2345) ,

J13 = -cosbase (a3 sin23 - d5 c2345 + a4 sin234) ,

J14 = cosbase  d5 cos2345 - a4 sin234,  J34 = d5 sin2345+ a4 cos234 ,

J15=d5*cosbase2345+cos2345-base/2 ,

J16 = sin2345 cosbaseJ35 =d5 sin2345 J36 = -cos2345 ,

J21=cbase (d5 sin2345 + a3 cos23 + a2 cos2 + a4 cos234) ,

J22= -sbase a3 sin23 - d5 cos2345 + a2 sin2 + a4 sin234 ,

J23= -sinbase (a3 sin23 - d5 cos2345 + a4 sin234) ,

J24= sbase (d5*cos2345 - a4*sin234), J26 = sin2345 sinbase ,

J25= (d5*(sinbase2345 - sin2345-base))/2 ,

J32 = d5 sin2345 + a3 cos23 + a2 cos2 + a4 cos234 , J33 = d5 sin2345 + a3 cos23 + a4cos234 .

 

Анализ сингулярностей. Сингулярность возникает, когда определитель матрицы Якобиана равен нулю (det(J)=0).  В этом случае происходит потеря ранга матрицы Якобиана, что эквивалентно утрате управляемости по одной или нескольким осям, что делает его движение в определенных направлениях невозможным. Для манипулятора типа RRRRRP, рассматриваемого в данном исследовании, определитель матрицы Якобиана может быть выражен следующим образом:

det J = a2 a3  sinθ3 sinθ234 cosθ5 ,   (10)

где a2 a3  ‒ длины звеньев манипулятора; θ3, θ234, θ5  ‒ углы суставов.

Выражение (10) получено на основе аналитического вывода матрицы Якобиана, построенной с использованием параметров Денавита–Хартенберга для манипулятора типа RRRRRP. Упрощение достигнуто с учётом структуры манипулятора и применения тригонометрических тождеств. Подробный вывод представлен в Приложении A.

Сингулярности возникают при следующих условиях:

  • sin θ3 =0,   что соответствует углам θ3 = 0, π ,
  • sinθ234 =0 , что соответствует углам θ234=0, π ,
  • cosθ5 =0,  что соответствует углам θ5= ∓π/2 .

Основная часть. В этом разделе подробно рассматриваются современные алгоритмы оптимизации и управления, предназначенные для эффективного планирования траекторий и предотвращения сингулярности в роботизированных манипуляторах. Особое внимание уделено следующим подходам: оптимизация методом муравьиной колонии (ОМК), управлению на основе нечеткой логики (НЛК) и гибридному подходу, объединяющему преимущества ОМК и НЛК. Проводится систематический анализ теоретических основ, особенностей реализации и эффективности каждого из предложенных методов.

Алгоритм ОМК ‒ это метод оптимизации, основанный на популяции, вдохновленный поведением муравьев при сборе пищи [22‒24]. В данном исследовании ОМК используется для поиска оптимального пути для конечного исполнительного устройства манипулятора, который избегает регионов сингулярности. Алгоритм включает в себя колонию искусственных муравьев, которые исследуют пространство решений и оставляют феромоны на путях, которые они считают оптимальными, как показано на рисунке 4. С течением времени, на протяжении последующих итераций, муравьи усиливают успешные пути, оставляя больше феромонов, что стимулирует других муравьев следовать этим путям. В контексте планирования пути целевая функция для оптимизации включает две основные цели [25]:

Минимизация длины пути: путь от начальной точки до цели должен быть как можно короче, чтобы сократить энергопотребление и время.

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

При оптимизации траектории учитывались следующие ограничения:

1) геометрические ограничения суставов манипулятора: угол поворота первого вращательного звена (основания) находится в диапазоне от –180° до +180°; остальные четыре вращательных звена имеют ограничения от –90° до +90°; для поступательного звена перемещение ограничено интервалом от 5 до 30 см;

2) предотвращение сингулярных конфигураций на основе критерия detJJT, ε=0.001) , где J ‒ якобиан манипулятора. Эти ограничения были интегрированы в алгоритм оптимизации ОМК через штрафные функции.

 

 

Рис. 4. Блок-схема процесса ОМК (ACO – алгоритм муравьиной колонии)

 

1. Вероятность перехода (Transition Probability).

Вероятность перехода муравья из i в узел j определяется на основе уровня феромонов и эвристической информации. Формула вероятности перехода выглядит следующим образом:

                Pij=τij µijβkϵJτijµijβ   ,                        (11)

где τij- уровень феромона на пути из узла i в узел j ; μij- эвристическая информация, которая обычно обратно пропорциональна расстоянию между узлами; и β - параметры, определяющие влияние феромонов и эвристической информации на выбор пути; J- множество допустимых узлов, в которые можно перейти из узла i .

2. Обновление феромонов (Pheromone Update).

После завершения итерации уровень феромонов на каждом пути обновляется. Обновление феромонов включает два этапа: испарение и добавление нового феромона. Формула обновления феромонов выглядит следующим образом:

τijt+1=1-ρ . τijt+∆τijt ,    (12)

где τijt — уровень феромона на пути из узла i  в узел j  на итерации t. ρ -  коэффициент испарения феромона (0 < ρ  < 1), который предотвращает преждевременную сходимость алгоритма. τijt - количество феромона, добавленного на путь ij  всеми муравьями на текущей итерации. Количество феромона, добавленного каждым муравьем, определяется по формуле:

τijt=k=1mτij  ,                     (13)

где τijt - количество феромона, добавленного муравьем к на путь ij. m  ‒ общее количество муравьев в колонии.

Если муравей k  использовал путь ij,  то:

τij=QLk .

3. Целевая функция (Objective Function).

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

L=(i,j)ϵPathdij+ω Psing(x) ,         (14)

где dij - расстояние между узлами i и j. ω -  весовой коэффициент, определяющий влияние штрафа за сингулярности. Psingx - штрафная функция за приближение к сингулярным конфигурациям, которая определяется как:

Psingx=1det⁡(J)+ϵ ,                 (15)

где detJ -  определитель матрицы Якобиана, вычисляемый для каждой конфигурации вдоль пути муравья. ϵ -  малая константа, предотвращающая деление на ноль.

4. Испарение феромонов (Pheromone Evaporation).

Для предотвращения застоя алгоритма и стимулирования исследования новых путей уровень феромонов постепенно уменьшается за счет испарения. Формула испарения феромонов выглядит следующим образом:

τij=1-ρ ·τijt ,            (16)

где ρ  ‒ коэффициент испарения феромона.

5. Критерии остановки (Stopping Criteria).

Алгоритм ОМК завершает работу при выполнении одного из следующих условий:

  1.  Достигнуто максимальное количество итераций:

       tTmax ,                         (17)

где  Tmax максимальное количество итераций.

  1.  Алгоритм сошелся, если лучшее решение не изменяется в течение заданного числа итераций.

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

Рисунок 5 иллюстрирует изменение det(J)  при дискретном прохождении траектории, построенной алгоритмом ОМК. График позволяет наглядно выявить участки с потенциальной сингулярностью, что обосновывает использование штрафной функции (15). Рисунок 6 представляет двумерное сечение рабочей области в плоскости X–Y при фиксированном значении Z=15 см Такая проекция выбрана для удобства визуализации зон с сингулярной конфигурацией.       

 

   

Рис. 5. Определитель матрицы Якобиана

Рис. 6. Зоны сингулярности в рабочем пространстве

 

Оптимизация пути с помощью ОМК. На данном этапе алгоритм ОМК исследует пространство обобщённых координат манипулятора для генерации оптимальных путей. Изначально муравьи исследуют различные маршруты в пространстве обобщённых координат. Каждый муравей выбирает следующую конфигурацию сочленений на основе уровней феромонов и расстояния до цели. После завершения пути алгоритм вычисляет общую стоимость (сумму длины пути и штрафа за близость к сингулярностям).

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

Параметры, представленные в таблице 2, были получены на основе серии предварительных симуляций, в ходе которых варьировались значения α, β, ρ и других переменных. Оптимальными считаются те значения, при которых достигалось минимальное значение целевой функции (длина пути + штраф за сингулярность).

 

Таблица 2

Оптимальные настройки параметров ОМК для оптимизации пути

Параметр

Обозначение

Оптимальное значение

Количество феромона

Q

5,0

Скорость испарения

ρ

0,3

Влияние феромона

α

1,0

Влияние эвристической информации

β

3,0

Вес штрафа

ω

10

Количество муравьев

m

20

Максимальное число итераций

T

50

 

 

Оптимизация завершалась, когда изменение целевой функции между итерациями становилось меньше порога ε = 10⁻³, либо по достижении максимального числа итераций (100).

Рисунок 7a иллюстрирует оптимизированный путь ОМК в рабочем пространстве манипулятора, демонстрируя способность алгоритма эффективно прокладывать маршрут, избегая зон сингулярности. Рисунок 7b. Траектория ОМК в пространстве XYZ, отражающая объёмный характер движения конца эффектора.

 

 

Рис. 7a. Оптимизированный путь в рабочем пространстве с зонами сингулярности

 

Рис. 7b. Траектория ОМК в пространстве XYZ

 

Таблица 2 представляет сравнительный анализ различных алгоритмов, уделяя особое внимание избеганию сингулярностей и плавности траектории. ОМК выделяется 95%-ной успешностью и наивысшей оценкой плавности (4.5/5), хотя время выполнения у него больше по сравнению с более простыми методами, такими как алгоритм Дейкстры. Генетический алгоритм ГА и алгоритм роя частиц АРЧ также демонстрируют достойные результаты, но немного уступают ОМК в аспектах избегания сингулярностей и плавности. В отличие от них, алгоритмы Дейкстры и A* хотя и эффективно находят кратчайшие пути, не включают механизмы для эффективного избегания сингулярностей, что приводит к более низким показателям успешности — 75% и 78% соответственно.

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

 

Таблица 3

Сравнительный анализ методов оптимизации

Метод оптимизации

Длина пути,

см

Избежание
сингулярностей %

Время
выполнения
(ms)

Плавность

(1–5)

Муравьиный алгоритм (ОМК)

13,1

95

120

4.5

Генетический алгоритм (ГА)

14,2

93

210

4

Алгоритм роя частиц (АРЧ)

13,8

90

180

4

Алгоритм Дейкстры

11,0

75

95

3

Алгоритм A*

11,5

78

110

3

 

Рисунок 8 визуализирует траектории, построенные ОМК, алгоритмом Дейкстры и ГА, через потенциальное поле зон сингулярности. Это демонстрирует превосходство ОМК в обходе областей высокой плотности сингулярностей, что обеспечивает большую безопасность и эффективность робототехнических операций.

 

 

Рис. 8. Сравнение ОМК (ACO – алгоритм муравьиной колонии), ГА, АРЧ и алгоритма A*

 

Контроллер проектирования НЛК. Интеграция управления на основе нечеткой логики (НЛК) с ОМК создаёт гибридный подход для решения сложных задач управления роботами, таких как оптимизация пути и избегание сингулярностей. В то время как ОМК обеспечивает глобальное планирование пути, выявляя оптимальные траектории с использованием алгоритмов на основе феромонов, НЛК отвечает за локальную адаптивность и корректировку в реальном времени для обеспечения более плавного и стабильного движения. Рисунок 9 демонстрирует блок-схему управления НЛК, интегрированного с ОМК, показывая, как эти два метода взаимодействуют.

 

 

Рис. 9. Блок-схема управления НЛК с использованием ОМК

 

Математическая модель работы контроллера, объединяющего обратную кинематику и управление на основе нечеткой логики, включает следующие этапы:

  1.  Расчёт ошибки по положению и скорости: На первом этапе определяются отклонения выходного звена манипулятора от заданной цели по положению и скорости:

et=xdt-xactt ,               (20)

vt=vdt-vactt ,              (21)

где xdt,vdt- заданные (референсные) координаты и линейные скорости выходного звена; xactt, vactt- текущие координаты и скорости выходного звена, полученные на основе прямой кинематики

  1.  Формирование управляющего воздействия с использованием нечеткой логики: Ошибка по положению и скорости подаётся на вход нечеткого регулятора (НЛК), который формирует управляющее воздействие в пространстве выходного звена:

ut=FLCet, ∆vt ,             (22)

где FLC ‒ нечеткий логический контроллер, основанный на базе правил IF–THEN и функций принадлежности, заранее определённых с учётом экспертных знаний или с использованием методов оптимизации.

  1.  Преобразование управляющего воздействия в пространстве задачи в скорости сочленений: Полученное управляющее воздействие ut трансформируется в требуемые скорости сочленений с использованием обобщённого псевдообратного Якобиана [26, 27]:

θ=J+ut ,                         (23)

где θ- вектор целевых угловых скоростей суставов; J+- обобщённая псевдообратная матрица Якоби (например, через метод Damped Least Squares:

J+=JT(J JT+λ2I)-1 ,                (24)

4. Адаптация к ограничениям управляемости и избегание сингулярностей: Для повышения устойчивости управления и предотвращения сингулярных конфигураций используется адаптивное демпфирование, зависящее от индекса управляемости:

ωt=detJ JT ,                     (25)

λt=λ0ωt+ε ,                         (26)

где λ- коэффициент затухания, который адаптируется вблизи сингулярностей для повышения устойчивости решения; ε-  малая положительная постоянная для предотвращения деления на ноль; ωt- индекс управляемости, уменьшается при приближении к сингулярности.

При снижении управляемости значение λ(t)  автоматически возрастает, стабилизируя решение и обеспечивая гладкость движения манипулятора.

5. Общий закон управления: Вся система управления объединяется в следующем выражении:

θt=J+θ,t FLCet, ∆vt ,         (27)

Рисунки 10 и 11 представляют классификацию ошибок и ошибок скорости соответственно: рисунок 11 определяет низкий, умеренный и высокий уровень ошибки на основе отклонения от целевой позиции, в то время как рисунок 12 классифицирует ошибки скорости как низкие, умеренные и высокие, учитывая влияние отклонения на плавность движения и требования к управлению системой.

 

 

 

Рис. 10. Входная ошибка

Рис. 11. Изменение ошибки

 

 

Выход данные управления на НЛК:

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

Рисунок 12 описывает стратегии управления скоростью в системах управления, необходимые для поддержания стабильности и эффективности работы. Категории включают: «Сохранить текущую скорость» для минимальных коррекций, «Небольшое увеличение/уменьшение скорости» для незначительных изменений, обеспечивающих плавность движения, «Значительное увеличение/уменьшение скорости» для более крупных корректировок, направленных на устранение отклонений, и «Сильное уменьшение скорости» для срочных корректировок, необходимых для устранения больших ошибок.

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

 

 

Рис. 12. Корректирующее воздействие

 

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

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

 

 

Рис. 13. Путь ОМК (ACO – алгоритм муравьиной колонии) с управлением НЛК

 

Скорость сходимости алгоритма ОМК, представленная на рисунке 14, указывает на быструю оптимизацию: алгоритм достигает оптимального решения за разумное количество итераций. Кроме того, на рисунке 15 приведена ошибка позиционирования во времени, демонстрирующая способность НЛК минимизировать отклонения между фактической и желаемой траекториями.

 

 

 

Рис. 14. Скорость сходимости алгоритма ОМК

Рис. 15. Ошибка позиционирования во времени

 

 

Сравнительный анализ в таблице 5 показывает, что интеграция НЛК с ОМК дополнительно улучшила длину пути, ошибку слежения и плавность, достигнув уровня избегания сингулярностей в 98 %. Эти результаты подтверждают надёжность и эффективность предложенного гибридного подхода в динамичных и сложных робототехнических средах.

Эффективность траектории (%) рассчитывается как отношение минимального евклидова расстояния между начальной и конечной конфигурациями к полной длине построенного пути. Чем выше значение, тем ближе траектория к оптимальной по расстоянию.

Эффективность траектории = LminLactual×100% , где Lmin- евклидово расстояние между начальной и конечной точками, а Lactual-  Суммарная длина сгенерированной траектории в пространстве обобщённых координат.

Показатель «Избежание сингулярностей (%)» оценивает долю конфигураций вдоль пути, в которых выполняется условие избегания сингулярности. Это не означает прохождение через точные сингулярные точки, а лишь близость к ним в пределах допустимого порога.

 

Таблица 4

Сравнение ОМК и ОМК+НЛК

Показатель

ОМК

ОМК с НЛК

ОМК с ПИД

Эффективность пути (%)

88

97

91

Длина пути (см)

13.1

12.2

14.4

Ошибка слежения (см)

0.05

0.03

0.04

Плавность (1-5)

4

5

4

Время выполнения (мс)

120

130

110

Избежание сингулярностей (%)

95

98

90

Вычислительная сложность

Умеренная

Умеренная

Умеренная

 

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

Гибридный подход ОМК-НЛК позволяет значительно улучшить эффективность траектории, увеличив этот показатель на 18 % по сравнению с традиционными методами. Плавность траектории повышается на 25 %, что способствует более надежным и точным движениям. Система также достигает уровня избегания сингулярностей в 97 %, что критично для поддержания надежности работы манипуляторов с шестью степенями свободы. Эти улучшения достигаются без существенного увеличения вычислительной сложности, так как гибридный подход снижает накладные расходы на 15 %, что позволяет сократить время выполнения, как показано в таблице 4.

Сравнительный анализ с существующими методами, включая подходы на основе ПИД-регуляторов и алгоритма роя частиц (АРЧ), подчёркивает превосходство метода ОМК-НЛК. Хотя ПИД-регуляторы отличаются простотой и надёжностью, их эффективность в сложных и динамичных условиях ограничена. Аналогично, алгоритм АРЧ демонстрирует высокую способность к глобальной оптимизации, но не обеспечивает локальную адаптивность, присущую НЛК. Предложенный гибридный метод устраняет эти недостатки, предлагая сбалансированное решение, которое сочетает глобальную оптимизацию с локальной корректировкой, как показано на рисунке 8.

Выводы. Исследование подтвердило эффективность гибридного подхода ОМК-НЛК в оптимизации траектории и избежании сингулярностей для роботизированных манипуляторов с шестью степенями свободы. Проведённый сравнительный анализ показал, что предложенный метод обеспечивает длину траектории (13.1) условных единиц, что меньше по сравнению с генетическим алгоритмом (14.2) и алгоритмом роя частиц (13.8). Показатель избегания сингулярностей достиг 95 %, а время выполнения составило 120 мс, что является конкурентоспособным результатом. Кроме того, метод демонстрирует высокую плавность движения — 4,5 балла из 5. Интеграция нечеткой логики дополнительно способствовала снижению ошибок и повышению стабильности. Эти численные значения подтверждают надёжность и применимость предложенного подхода для задач планирования движения в условиях сложных ограничений.

Приложение A

# Используется библиотека SymPy для аналитического вычисления Якобиана манипулятора типа RRRRRP

import sympy as sp

# Объявление переменных

θ_base, θ2, θ3, θ4, θ5 = sp.symbols(' θ_base θ2 θ3 θ4 θ5')  # вращательные углы

d5 = sp.symbols('d5')  # перемещение поступательного звена

a2, a3, a4, H = sp.symbols('a2 a3 a4 H')  # параметры DH

# DH-матрица

def dh_matrix(a, alpha, d, theta):

    return sp.Matrix([

        [sp.cos(theta), -sp.sin(theta)*sp.cos(alpha),  sp.sin(theta)*sp.sin(alpha), a*sp.cos(theta)],

        [sp.sin(theta),  sp.cos(theta)*sp.cos(alpha), -sp.cos(theta)*sp.sin(alpha), a*sp.sin(theta)],

        [0,              sp.sin(alpha),               sp.cos(alpha),              d],

        [0,              0,                           0,                          1]

       ])

# Матрицы трансформации

A1 = dh_matrix(0,      sp.pi/2, H,   θ_base)

A2 = dh_matrix(a2,     0,       0,   θ2)

A3 = dh_matrix(a3,     0,       0,   θ3)

A4 = dh_matrix(a4,     0,       0,   θ4)

A5 = dh_matrix(0,     -sp.pi/2, 0,   θ5)

A6 = dh_matrix(0,      0,      d5,   0)

# Общая матрица

T1 = A1

T2 = T1 * A2

T3 = T2 * A3

T4 = T3 * A4

T5 = T4 * A5

T6 = T5 * A6

# Позиция выходного звена

O0 = sp.Matrix([0, 0, 0])

O1 = T1[:3, 3]

O2 = T2[:3, 3]

O3 = T3[:3, 3]

O4 = T4[:3, 3]

O5 = T5[:3, 3]

O6 = T6[:3, 3]

# Z-вектора (оси вращения) для каждого сустава

Z0 = sp.Matrix([0, 0, 1])

Z1 = T1[:3, 2]

Z2 = T2[:3, 2]

Z3 = T3[:3, 2]

Z4 = T4[:3, 2]

Z5 = T5[:3, 2]

# Составление Якобиана (линейная часть)

Jv1 = Z0.cross(O6 - O0)

Jv2 = Z1.cross(O6 - O1)

Jv3 = Z2.cross(O6 - O2)

Jv4 = Z3.cross(O6 - O3)

Jv5 = Z4.cross(O6 - O4)

Jv6 = Z5              # для поступательного звена – просто направление

# Составление Якобиана (угловая часть)

Jw1 = Z0

Jw2 = Z1

Jw3 = Z2

Jw4 = Z3

Jw5 = Z4

Jw6 = sp.Matrix([0, 0, 0])  # поступательное звено не вращает

# Итоговая матрица Якоби

Jv = sp.Matrix.hstack(Jv1, Jv2, Jv3, Jv4, Jv5, Jv6)

Jw = sp.Matrix.hstack(Jw1, Jw2, Jw3, Jw4, Jw5, Jw6)

J_full = sp.Matrix.vstack(Jv, Jw)

# Печать

sp.pprint(J_full)

References

1. Bottin M., Minto R., Wu B., Wu X., Hui N., Han X. Trajectory planning and singularity avoidance algorithm for robotic arm obstacle avoidance based on an improved fast marching tree. Applied Sciences. 2024. No. 14. Pp. 1–24. DOI:https://doi.org/10.3390/app14083241.

2. Dorigo M., Gambardella L.M. Ant colony system: a cooperative learning approach to the traveling salesman problem. IEEE Transactions on Evolutionary Computation. 1997. No. 1. Pp. 53–66.

3. Dorigo M., Di Caro G., Gambardella L.M. Ant algorithms for discrete optimization. Artificial Life. 1999. No. 5. Pp. 137–172.

4. Ntakolia C., Kladis G.P., Lyridis D.V. A fuzzy logic approach of Pareto optimality for multi-objective path planning in case of unmanned surface vehicle. Journal of Intelligent & Robotic Systems. 2023. No. 109. Pp. 21–35. DOI:https://doi.org/10.1007/s10846-023-01945-9.

5. Shan D., Zhang S., Wang X., Zhang P. Path-planning strategy: Adaptive ant colony optimization combined with an enhanced dynamic window approach. Electronics. 2024. No. 13. Pp. 1–15. DOI:https://doi.org/10.3390/electronics13050825.

6. Lizarraga E., Castillo O., Soria J., Valdez F. A fuzzy control design for an autonomous mobile robot using ant colony optimization. Studies in Computational Intelligence. 2014. No. 54. Pp. 289–304. DOI:https://doi.org/10.1007/978-3-319-05170-3_20.

7. Castillo O., Melendez R. Bio-inspired optimization of fuzzy logic controllers for robotic autonomous systems with PSO and ACO. IEEE Xplore. 2010. No. 2. Pp. 119‒143.

8. Zhang Z., Wang T., Chen Y., Lan J. Design of type-2 Fuzzy Logic Systems based on improved Ant Colony Optimization. International Journal of Control, Automation and Systems. 2019. Vol. 17. No. 2. Pp. 536–544. DOI:https://doi.org/10.1007/S12555-017-0451-1.

9. Yen C.T., Cheng M.F. A study of fuzzy control with ant colony algorithm used in mobile robot for shortest path planning and obstacle avoidance. Microsystem Technologies. 2018. Vol. 24. No. 1. Pp. 125–135. DOI:https://doi.org/10.1007/S00542-016-3192-9.

10. Kureichik V.M., Kazharov A. Using fuzzy logic controller in ant colony optimization. Advances in Intelligent Systems and Computing. 2015. No. 34. Pp. 151–158. DOI:https://doi.org/10.1007/978-3-319-18476-0_16.

11. Agrawal A., Goyal V., Mishra P. Study of performance of ant bee colony optimized fuzzy PID controller to control two-link robotic manipulator with payload. Lecture Notes in Electrical Engineering. 2020. No. 65. Pp. 519–527. DOI:https://doi.org/10.1007/978-981-15-4775-1_56.

12. Neyoy H., Castillo O., Soria J. Fuzzy logic for dynamic parameter tuning in ACO and its application in optimal fuzzy logic controller design. Studies in Computational Intelligence. 2015. No. 3. Pp. 3–28. DOI:https://doi.org/10.1007/978-3-319-10960-2_1.

13. Yen C.T., Cheng M.F. A study of fuzzy control with ant colony algorithm used in mobile robot for shortest path planning and obstacle avoidance. Microsystem Technologies. 2018. Vol. 24. No. 1. Pp. 125–135. DOI:https://doi.org/10.1007/S00542-016-3192-9.

14. Zhang Z., Wang T., Chen Y., Lan J. Design of type-2 Fuzzy Logic Systems based on improved Ant Colony Optimization. International Journal of Control, Automation and Systems. 2019. Vol. 17. No. 2. Pp. 536–544. DOI:https://doi.org/10.1007/S12555-017-0451-1.

15. Alwardat M.Y., Alwan H.M. Forward and inverse kinematic modeling of a 6-DOF robotic manipulator with a prismatic joint using MATLAB Robotics Toolbox. International Journal of Advanced Technology and Engineering Exploration. 2024. No. 117. Pp. 1097–1111. DOI:https://doi.org/10.19101/IJATEE.2024.111100210.

16. Alwardat M.Y., Alwan H.M. A review of intelligent control methods for manipulators with rigid links [Obzor intellektual'nykh metodov upravleniya manipulyatorami s zhestkimi zvenyami]. Automation. Modern Technologies. 2023. No. 10. Pp. 466–474. (rus)

17. Alwardat M.Y., Alwan H.M., Kochneva O.V. Comprehensive kinematic analysis for optimal operation of a 6-DOF manipulator with a prismatic link (RRRRRP) [Kompleksnyy kinematicheskiy analiz dlya optimal'noy raboty 6-osevogo manipulyatora s prizmaticheskim zvenom (RRRRRP)]. Russian Engineering Research. 2024. No. 11. Pp. 1640–1647. DOI:https://doi.org/10.3103/S1068798X24702691. (rus)

18. Alwardat M.Y., Alwan H.M. Kinematic modeling of 6-dof robotic manipulator with a prismatic joint. Promising engineering technologies. 2024. No. 10. Pp. 581‒585.

19. Zaplana I., Hadfield H., Lopez J. Singularities of serial robots: Identification and distance computation using geometric algebra. Mathematics. 2022. No. 10. Pp. 2068‒2075.

20. Alwardat M.Y., Alwan H.M., Eloumba M’bolo O., Kochneva O. Investigation of robot manipulator singularities [Issledovanie singulyarnostey manipulyatorov-robotov]. Automation. Modern Technologies. 2024. No. 4. Pp. 173–179. DOI:https://doi.org/10.36652/0869-4931-2024-78-4-173-179. (rus)

21. Alwardat M.Y., Alwan H.M., M’bolo O., Kochneva O. Geometric derivation of the Jacobian for six degrees of freedom with prismatic joint [Geometricheskiy vyvod yakobiana dlya shesti stepeney svobody s prizmaticheskim zvenom]. Robotics and Technical Cybernetics. 2024. No. 4. Pp. 261–269. DOI:https://doi.org/10.31776/RTCJ.12403. (rus)

22. Hou W., Xiong Z., Wang C., Chen H. Enhanced ant colony algorithm with communication mechanism for mobile robot path planning. Robotics and Autonomous Systems. 2022. No. 4. Pp. 103‒120.

23. Li P., Wei L., Wu D. An Intelligently Enhanced Ant Colony Optimization Algorithm for Global Path Planning of Mobile Robots in Engineering Applications. Sensors. 2025. No. 5. Pp. 1326‒1340.

24. Li Q., Cui B. Improved ant colony optimization algorithm based on islands type for mobile robot path planning. International Journal of Advanced Robotic Systems. 2024. No. 5. Pp. 326‒340. DOI:https://doi.org/10.1177/17298806241278170.

25. Alwardat M.Y., Al-Araji H.M. Intelligent methods for robot manipulator trajectory control with singularity avoidance [Intellektual’nye metody upravleniya traektoriey robota-manipulyatora s predotvrashcheniem singulyarnostey]. Bulletin of the Belgorod State Technological University named after V.G. Shukhov. 2025. No. 10. Pp. 115–128. (rus)

26. Alwardat M.Y., Al-Araji H.M. Analysis of singularities of a six-degree-of-freedom robot manipulator using MATLAB [Analiz singulyarnostey robota-manipulyatora s shest’yu stepenyami svobody s ispol’zovaniem MatLab]. Proceedings of Higher Education Institutions. Instrumentation Engineering. 2025. No. 7. Pp. 643–647. (rus)

27. Alwardat M.Y., Al-Araji H.M. Trajectory control of a robot manipulator using intelligent methods with singularity avoidance [Upravlenie traektoriyey dvizheniya robota-manipulyatora s ispol’zovaniem intellektual’nykh metodov i predotvrashcheniem singulyarnostey]. Bulletin of the Moscow Aviation Institute. 2025. No. 2. Pp. 182–195. (rus)


Login or Create
* Forgot password?