- •Московский государственный технический университет им.Н.Э.Баумана
- •Манипуляционные роботы
- •Ису (су манипулятора)
- •Раздел 1. Кинематика им мр
- •Координаты и параметры им
- •Описание им мр
- •Системы координат
- •Координаты шарниров
- •Параметры им
- •Обобщенные координаты им
- •Матрицы преобразований поворота
- •1.6.1. Матрицы поворота вокруг одной координатной оси (матрицы элементарных поворотов).
- •1.6.2. Матрицы композиции элементарных поворотов
- •Этот факт иллюстрирует рисунок 1.10. Зеленый цвет – последовательность поворотов ск. Синий – последовательность преобразований векторов, заданных координатами в подвижной ск, в неподвижную.
- •1.6.3. Матрицы сложных поворотов
- •1.6.4. Матрицы поворота ск смежных звеньев им
- •Очевидно,
- •1.6.5. Матрица преобразования поворота ск произвольных звеньев
- •Положение им в рабочем пространстве. Линейные и угловые координаты звеньев
- •2.1. Линейные координаты звеньев им относительно базовой ск
- •Прямая позиционная кинематическая задача
- •Понятие обратной позиционной кинематической задачи
- •2.4.1 Сингулярности окз
- •2.4.2 Обратная кинематическая задача в управлении манипулятором
- •Однородные координаты и однородные преобразования в робототехнике
- •2.6 Блочные матрицы
- •2.7 Специальные системы координат
- •Угловая и линейная скорости звеньев им
- •3.1. Угловые скорости звеньев им
- •Введем обозначения
- •3.2. Свойства векторов сj(0) .
- •3.3. Матрицы-якобианы угловых скоростей им.
- •3.4. Производные по времени матриц преобразований поворота.
- •3.4.1. Запись производной по времени матрицы I-1i
- •3.4.2. Свойства матрицы векторного произведения (.)
- •3.4.3. Вычисление `I-1i .
- •3.4.4. Производная по времени матрицы 0i.
- •3.5. Линейные скорости звеньев им
- •(Различные индексы суммирования во втором и третьем слагаемых приняты для удобства дальнейших преобразований).
- •Развернем полученное выражение, подставив в него
- •Получим: I I
- •3.6. Свойства векторов Dij(0)
- •3.7. Матрицы-якобианы линейных скоростей им.
- •3.8. Матрица-якобиан последнего звена им.
- •3.9. Прямая кинематическая задача о скорости им
- •3.10. Понятие обратной задачи о скорости им
- •3.11. Примеры решения прямой и обратной задач о сокорости звеньев им
- •3.11.1. Прямая и обратная кинематические задачи об угловой скорости звеньев им
- •3.11.1. Прямая и обратная кинематические задачи о линейной скорости звеньев им
- •Угловые и линейные ускорения звеньев
- •4.1. Угловые ускорения
- •4.2. Линейные ускорения
- •Раздел 2. Динамика им мр
- •Силы и моменты.
- •Силы, действующие на звенья им
- •Моменты, действующие на звенья им.
- •Сумма двух последних слагаемых дает
- •С учетом этого выражение для I (0) примет вид:
- •Выражение для Mиj(0) имеет вид:
- •Силы, действующие вдоль осей шарниров им.
- •1.4. Моменты, действующие вдоль осей шарниров.
- •1.5. Силы и моменты, развиваемые приводами манипулятора. Обратная задача динамики им.
- •1.6. Уравнения статического равновесия им.
- •Симметрия уравнений озд и уравнений для определения скоростей звеньев им.
- •Управление по вектору силы
- •Уравнения движения им мр
- •Вывод уравнений движения
- •С учетом выполненных преобразований и введенных обозначений первое выражение п.2.1 примет вид:
- •2.4. Запись уравнений движения им с использованием матриц преобразования между ск звеньев.
- •2.3. Пример расчета матрицы инерционных коэффициентов.
Прямая позиционная кинематическая задача
Определение линейных координат точек на ЗУМ и ориентации СК ЗУМ относительно базовой СК (определение углов Эйлера) при заданных значениях координат шарниров и заданных параметрах звеньев ИМ носит название прямой позиционной кинематической задачи.
Решение прямой позиционной задачи выражается формулами
n
tn(0)(0) = 0n tn(n)(n) + 0j (j (1-j)qj + lj) ,
j =1
где
n
0n = П k-1 k
k=1
Прямая позиционная кинематическая задача всегда имееет решение и притом единственное. Место прямой задачи в управлении МР - РИСУНОК, ИЛЛЮСТРИРУЮЩИЙ КОНТУРНОЕ УПРАВЛЕНИЕ
Понятие обратной позиционной кинематической задачи
Определение координат сочленений по заданным линейным координатам ЗУМ и углам его ориентации относительно базовой СК при известных параметрах звеньев ИМ носит название обратной кинематической задачи о положении ИМ или обратной позиционной кинематической задачи (сокращенно, ОКЗ).
В общей постановке для многих ИМ обратная позиционная кинематическая задача может иметь несколько решений, бечсисленное множество решений или не иметь ни одного решения.
Первый случай связан с тем, что одному и тому же положению и ориентации ЗУМ относительно базовой СК может соответствовать несколько разных наборов значений координта шарниров.
Рисунок иллюстрирует один из таких случаев. Одному и тому же положению и ориентации ЗУМ изображенного на этом рисунке ИМ соответствуют 2 набора значений координат шарниров:
q1 = 0, q2 = -30 град., q3 = 0 град. (на рисунке слева)
и
q1 = 180 град, q2 = +30 град., q3 = 180 град. (на рисунке справа)
Второй случай соответствует определённым конфигурациям кинематической цепи. Для изображенногго на вышеприведенном рисунке ИМ это случай соответствует случаю, когда все звенья как бы вытянуты с линию.
Третий случай связан с ограничениями рабочей зоны вследствие конструктивных ограничений возможных значений координат сочленений и геометрических параметров звеньев ИМ.
Для обеспечения единственности решения ОКЗ применяют специальные меры.
Траектории перемещения ЗУМ предварительно планируют (ТСУ) так, чтобы исключить выход шарниров на ограничения. При этом также исключают попадание в зоны бесчисленного множества решений.
Для обеспечения единственности решения ОКЗ (в «неоднозначном» случае) в рассмотрение вводятся дополнительные параметры (индексы конфигурации кинематической цепи), позволяющие однозначно определить значения координат шарниров. Выбор этих параметров зависит от вида кинематической схемы ИМ и алгоритма решения ОКЗ.
Для изображенного на рисунке ИМ в качестве такового удобно использовать значение величины Ks, заданной нижеследующим выражением
Ks = sign (Z1 · ( Z3 х Z0)),
где
Z1 , Z3 , Z0 - единичные векторы в направлении осей Z СК1, СК3 и СК0, соответственно,
· - символ скалярного, а х - символ векторного произвеления векторов.
Для рассматриваемого ИМ в положении, изображенном на рисунке 1 слева Ks = 1. а в положении справа, Ks = -1.
Введение в рассмотрение индекса конфигурации Ks удобно проиллюстрировать графически. Так, в случае Ks = 1 при увеличении значения кординаты q2 ЗУМ приближается к линии, проведенной вдоль вертикальной оси (плечо находится как бы "впереди" по ходу движения ЗУМ). В случае Ks = -1 при увеличении значения кординаты q2 ЗУМ удаляется от вертикальной линии (плечо находится как бы "сзади" по ходу движения ЗУМ).
ПРИМЕР решения ОКЗ.
Рассмотрим 3-звенный ИМ на рис.1.
Будем считать известными углы ориентации (Pitch, Yaw, Roll) осей этой СК3 относительно базовой СК (СК0).
Pitch = -30 град,
Yaw = Roll = 0.
Будем считать также заданным значение индекса конфигурации Ks:
Ks = 1.
По формулам п. получим значения матрицы τ 03 . В нашем случае роль СК {xyz} играет СК0, а СК {uvw} – СК3.
τ03 = τ xyz, uvw (ось u, угол Pitch) τ xyz, uvw (ось v`, угол Yaw) τ xyz, uvw (ось w``, угол Roll)
τ 03 = 1 0 0
0 √3/2 ½
0-1/2 √3/2
Последний столбец этой матрицы (обозначим его Z3(0)) есть вектор Z3, заданный в проекциях на оси СК0, т.е. Z3(0).
Z3(0)= (0 1/2 √3/2)T.
Зная Z3, вычислим вектор вектор Z1:
Z1 = Ks (Z3 х Z0)/|(Z3 х Z0)|.
В координатной форме последнее выражение записывается так
Z1(0) = Ks λ(Z3(0)) Z0(0) /|(Z3 х Z0)|= - Ks λ(Z0(0)) Z3(0) /|(Z3 х Z0)|.
где
λ(Z3) - матрица векторного произведения, составленная из элементов вектора Z3 = (x y z)T:
λ(Z3) = 0 -z y
z 0 –x
-y x 0
Z0(0) = (0 0 1)T.
Z1(0) = - (+1) 0 -1 0
1 0 0
0 0 0 (0 ½ √3/2)T = - (- ½ 0 0 )T / (1/2) = (1 0 0)T.
Вычислим угол q1:
q1 = Arctg (sin q1 / cos q1)
cos q1 = X0(0) · Z1(0) = 1,
(X0(0) - единичный вектор в направлении оси X СК0 в проекциях на оси этой же СК; X0(0) = (1 0 0 )T.
sin q1 = Z0(0) · λ (X0(0)) Z1(0) = (0 0 1) 0 0 0
0 0-1
01 0 ( 0 0 1)T = (0 0 1) (0-1 0)T = 0.
Отсюда
q1 = 0.
Вычислим угол q3:
q3 = Arctg (sin q3 / cos q3)
cos q3 = X3(0) · Z1(0).
X3(0) - единичный вектор в направлении оси X СК3 в проекциях на оси СК0; этот вектор - первый столбец матрицы τ 03 :
X3(0) = (1 0 0 )T.
сos q3 = X3(0) · Z1(0) = (1 0 0 ) (1 0 0)T = 1.
sin q3 = Z3(0) · (λ (Z1(0)) х X3(0)) = (0 1/2 √3/2) 0 -1 0
1 0 0
0 0 0 (100)T = 0.
Отсюда
q3 = 0.
Вычислим угол q2.
q2 = Arctg (sin q2 / cos q2)
cos q3 = Z0(0) · Z3(0) = (001) (0 ½ √3/2 )T = √3/2.
sin q3 = Z1(0) · (λ (Z0(0)) Z3(0)) = (1 0 0) 0-1 0
1 0 0
0 0 0 (01/2√3/2)T = (1 0 0) (-1/2 0 0)T = -1/2.
q2 = -30 град.
Отметим, что в общем случае решение ОКЗ требует задания помимо углов ориентации, координат точки t3(0)(0) начала СК3. Однако в данном случае в силу ограниченности числа степеней свободы ИМ, достаточно задать только параметры ориентации. В данном примере решение ОКЗ можно получить, задав только линейные кординаты ЗУМ. Последовательность вычислений при этом останется прежней с той лишь разницей, что вектор Z3(0) следует вычислить как единичный вектор в направлении линии,связывающей начала 1-й и 3-й СК.
