2004
Доклады Б ГУ ИР
ОКТЯБРЬ-ДЕКАБРЬ
№ 4
УДК 681.5
ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ ДВИЖЕНИЕМ ИНСТРУМЕНТА В РОБОТОТЕХНИЧЕСКОМ КОМПЛЕКСЕ РЕЗКИ С МАНИПУЛЯТОРОМ ИЗДЕЛИЯ
А.П. ПАШКЕВИЧ, О .А. ЧУМАКОВ
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 15 июня 2004
Предложены методы, позволяющие автоматизировать процесс генерации управляющей программы для РТК резки трехмерных объектов сложной формы. При этом оптимально используется кинематическая избыточность РТК с шестистепенным роботом и двухстепенным позиционером. В результате генерируется плавная траектория движения инструмента манипулятора, которая может быть воспроизведена в реальном масштабе времени с заданной точностью.
Ключевые слова: робот, позиционер, off-line программирование, автоматизация проектирования.
Введение
Автоматизация современного производства связана с использованием промышленных роботов, которые выполняют большинство основных технологических операций: сварку, сборку, окраску, и др. В настоящее время программирование многих технологических операций осуществляется при помощи САПР робототехнических комплексов (РТК), которые позволяют создавать управляющую программу на основе графической модели, включающей робот-манипулятор инструмента, манипулятор изделия (позиционер), систему управления, элементы ограждения и транспортировки деталей [1]. Такой подход не требует полной остановки РТК и его исключения из производственного процесса. Однако существующие САПР РТК [2] ориентированы в основном на технологию сварки и неприменимы для автоматизации программирования процесса резки трехмерных объектов сложной формы. Для операции резки таких объектов достаточно пять степеней подвижности, хотя в настоящее время наибольшее распространение получили манипуляторы с шестью степенями подвижности, которые применительно к данной технологической операции обладают кинематической избыточностью (это обусловлено возможностью произвольного вращения инструмента вокруг режущего луча без влияния на качество процесса). Кроме того, дополнительную избыточность образуют и степени подвижности позиционера. Эта избыточность существенно усложняет программирование, однако ее можно использовать при разработке управляющих программ для оптимизации движений суставов в пространстве обобщенных координат. Поэтому в данной работе рассматривается проблема автоматизации программирования РТК резки для системы с шестистепенным манипулятором и двухстепенным позиционером.
Формирование множества допустимых локаций инструмента
В рассматриваемых РТК промышленный робот используется в качестве "манипулятора инструмента", т.е. для перемещения инструмента вдоль траектории, а позиционер — в качестве "манипулятора изделия", т.е. для изменения положения детали в пространстве. Предположим, что расположение системы координат робота "Base" и позиционера "Positioner" относительно
глобальной системы "World" описывается 4х4-матрицами однородных преобразований и
positioner соответственно (рис. 1), а каждое звено робота и позиционера описываются матричными функциями Т¡(qi) и Т,(фу) от обобщенных координат звеньев <у, и ф,. Тогда кинематическую модель рассматриваемого РТК можно описать матричным уравнением
Tjfo)Т2(д2)..Л5(д5)Т6(д6)Т7(Ф,)Т8(ф2) Шап*еТТоо1, (1)
где Flangc Т]оЫ — матрица, описывающая положение и ориентацию системы координат инструмента "Tool" по отношению к системе координат крепежного фланца "Flange".
Для формализации рассматриваемой задачи необходимо представить траекторию резки в виде совокупности координатных фреймов. При этом для каждого сегмента траектории примем следующие соглашения: ось Xt направлена вдоль траектории резки; ось Z, указывает направление нормали к поверхности изделия; ось Уг дополняет оси Xh Z, до правой тройки. Тогда соответствующий фрейм сегмента может быть описан матричным уравнением
О >у _0 >у Positioner у Wpiece у
Segm ^ Positioner Wpiece Segm ' V /
О т
где однородная матрица iPositloner описывает расположение системы координат позиционера
« Wpiece т* Positioner гр
относительно мировои системы координат, а матрицы и * wpiece — располо-
жение координатного фрейма 7-го сегмента траектории относительно детали и фрейма детали относительно системы координат позиционера соответственно.
Траектория резки, используемая в современных САПР РТК, создается на основе каркасной или поверхностной модели изделия и представляется в виде набора узловых точек ш' [3], каждая из которых включает вектор pt (декартовы координаты рабочей точки инструмента) и единичный вектор ориентации инструмента пъ перпендикулярный к обрабатываемой поверхности. Для каждой узловой точки траектории ш' определим систему координат и обозначим как "Node". Начало координат этой системы разместим в точке р', ось X направим вдоль траектории резки, ось Z - в направлении режущего луча а\, а ось Y — так, чтобы X, Y, Z образовывали правую тройку векторов. Соответствующую матрицу однородных преобразований обозначим °TSegm . Тогда манипуляционную задачу робота можно представить как совмещение координатных осей фреймов и фреймов ° T7W (q) таким образом, ЧТО ОСЬ Xfool
совпадала с осью ZSegm, а взаимное расположение оставшихся осей определяется с точностью до произвольного скалярного параметра у, описывающего поворот рассматриваемых фреймов вокруг оси резки:
С К фланца робота
"Flange"
Flange гр
С К детали "Workpiece" и позиционера
"FlangePos"
1
позиционер^ "Positioner" \
0-у ' Pontioner
Мировая СК
"World' j
С К основания робота "Base"
Рис. 1. Системы координат манипулятора инструмента (робота) и манипулятора изделия (позиционера)
¡■цт, ^.-W ^-'fool
(3)
где постоянная матрица Н = КЛ.(71)К1 (к/2) учитывает порядок совмещения осей. аЯх Д.-
операторы вращения вокруг X. У. 2 соответственно. В результате, применяя соответствующий алгоритм решения обратной задачи кинематики
<ШФ1,ср2) = 1пуКп1"Т^ТлЛ^(ф1,ф2)К.(у)и ']. (4)
исходная совокупность сегментов траектории резки О' преобразуется в трехпараметрическое семейство обобщенных координат О,(у. фь ф). На практике целесообразно ограничиться дискретными значениями параметров у. фь ф и множества {0,(у, фь фг): уе |(),7г |; ф1е [ф1 ппп, ф1 пигх ]; ф2 е [фшш- ф2шш;]} описать многомерной матрицей, столбцы которой соответствуют одной и той же физической точке траектории резки, а элементами являются векторы обобщенных координат ^ либо соответствующие им локации инструмента Ь. При этом для каждого элемента этой матрицы необходимо проверить существование решения обратной задачи кинематики, близость к точкам сингулярности и отсутствие столкновений робота и режущего инструмента с другими элементами РТК. Применяя указанные операции ко всем узлам траектории резки П'. получаем, что множество допустимых локаций инструмента и соответствующих обобщенных координат представляется в виде "разряженной" четырехмерной матрицы {Луи, ц^ыУ, каждый элемент которой представляет собой совокупность фрейма Ь,^, шестимерного вектора (¡^ и вектора ориентации изделия фи.
Синтез оптимальной траектории движения инструмента
Набор матриц локаций описывается многослойным направленным графом с вершинами V ={Ьу} и ребрами Е = {(Ьу,Ьи); 1-к 1; \Zi\j\kJ\ (рис. 2), где каждый слой соответствует определенному положению позиционера для заданного индекса конфигурации робота М. Тогда функция у(0 и две координаты позиционера изделия фД/), к= 1,2, определяют шесть тра-
екторий с} ставов робота qM. к= 1...6. каждая из которых может быть оценена скалярным критерием который характеризует "объем движений" суставов, амплитуду приращений
координат (или их взвешенную сумму) либо диапазон изменения координат.
В результате рассматриваемая задача планирования движения робота сводится к следующей оптимизационной задаче на графе: для заданных множеств вершин V и ребер Е найти
"наилучший" путь блины п П(уп,...уп) = —>L1; —>...Ьи; ^ с начальным состоянием
V0e{Loj} и конечным состоянием V„e{L,v}. который обеспечивает Парето-оптималъностъ векторного критерия качества (5). Для упрощения описания алгоритмов обобщенные координаты, соответствующие положению суставов робота и позиционера Ьщ, обозначим как qi(i.j.k.l). а траектории, соответствующие вектору решения Г, — как qi^ij k, /,). Таким образом, используя эти обозначения, исходная задача формулируется в терминах теории комбинаторной оптимизации как поиск "наилучшего" пути, который минимизирует желаемый показатель каче-
Для минимизации объема движения для каждой обобщенной координаты
j[k (П=Zk ,Д-, in
i
применим метод динамического программирования. Предположим, что для всех конечных точек. принадлежащих множеству допустимых значений Qp{у, cpb срг), получены оптимальные решения для всех траекторий длины р. Тогда оптимальные решения F{.) для траекторий длины р+1 (с конечными точками Qp \ (у. фь <р )) могут быть получены на основе следующей рекуррентной процедуры:
f 1(у,ф1,ф2)= min{f;(x,e,,e2),+|^(/7+i, у х,е,,е2)|}
Таким образом, начиная с р= 1 и последовательно увеличивая длину последовательности Т(р, у, ф), для каждого конечного состояния могут быть найдены и оптимальный путь, и соответствующее значение показателя качества. Как в предыдущем случае, на последнем шаге происходит выбор лучшего конечного состояния из последовательности у, феГ„. Подобная рекурсия также дает оптимальное решение для взвешенной суммы и "наихудшего" компонента соответственного векторного показателя качества.
Рис. 2. Графовая модель пространства поиска (СО) — отсутствие решений обратной задачи кинематики; — нарушение теста столкновений; — приближение к точкам сингулярности)
Таким образом, задача синтеза траектории резки сводится к нахождению некоторой оптимальной кривой в четырехмерном пространстве с запретными областями и линейными граничными условиями. Для рассматриваемой задачи, требующей только пять степеней свободы, избыточный параметр у и две степени подвижности позиционера используются для сглаживания траектории в пространстве обобщенных координат, чтобы избежать резких поворотов режущего инструмента.
Рассматриваемые технические требования невозможно полностью описать с помощью скалярной целевой функции, хотя в идеальном случае все рассмотренные выше критерии качества должны стремится к нулю, на практике минимизация одной компоненты критерия может увеличить значение другой. Поэтому разработчик должен выбрать одну из методик, которые обычно используются в векторной оптимизации. В этой статье не отдается предпочтение какому-либо одному критерию качества, а предлагается перенести окончательное решение на этап проектирования, когда оно может быть выбрано конструктором путем назначения приоритетов или построения глобального критерия.
Заключение
Разработанные методы позволяют автоматизировать процесс генерации управляющей программы для РТК резки трехмерных объектов сложной формы. Предложенные алгоритмы реализованы в подсистеме Laser САПР ROBOMAX и были использованы при решении реальных производственных задач по автономному программированию РТК резки. Как следует из опыта эксплуатации, они позволяют эффективно синтезировать и оптимизировать программы роботизированной резки трехмерных деталей сложной формы и достичь существенного повышения экономической эффективности при мелкосерийном производстве.
OPTIMAL CONTROL OF INSTRUMENT MOVEMENT IN ROBOT-TECHNICAL CUTTING SYSTEM WITH ARTICLE MANIPULATOR
A.P. PASHKEVICH, O A. CHUMAKOV Abstract
Means for generation process automation of resident programme for robot-technical three-dimensional complex shape objects cutting system are proposed. At that kinematic redundancy of robot-technical system with sixdegree robot and bidegree positioner is optimally used. As a result smooth mechanical trajectory of manipulator instrument which can be reproduced in real time with prescribed accuracy is generated.
Литература
1. Geiger M., Otto A. Laser in der Elektronikproduktion & Feinwerktechnik, Tagungsband des 3. Erlanger Seminars LEF 2000. Bamberg: Meisenbach, 2000.
2. Пашкевич А.П. И Справочник. Инженерный журнал. 1997. № 5. С. 27-33