Решение прямой задачи кинематики KUKA youBot
Мобильный манипулятор KUKA youBot на рис. 16 состоит из установленных на платформе пяти последовательно соединённых звеньев, связанных соединениями вращательного типа с платформой и между собой.
Рис. 16 Схема мобильного манипулятора KUKA youBot
Будем считать, что звено 1 – твёрдое тело, вращающееся вокруг вертикальной оси С1Z1, точка С1 принадлежит основанию платформы, звенья 2, 3, 4 – стержни, которые вращаются вокруг горизонтальных осей, проходящих через точки C2, C3, C4 соответственно. Звено 5 – твердое тело, ладонь, вращается вокруг оси С5z, направленной вдоль звена 4. На рис. 16 ладонь представлена отрезком С5С'5 оси С5z. Со звеном 5 в точке С'5 соединён схват манипулятора в виде двух пальцев, которые могут перемещаться поступательно относительно звена 5, точка E на оси С5z – крайняя точка схвата при сомкнутых пальцах. Платформа – твёрдое тело с горизонтальным верхним основанием, перемещающееся по неподвижной горизонтальной плоскости с помощью четырёх меканум колёс. Звенья 1–4 и ось вращения звена 5 расположены в одной вертикальной плоскости.
Введём системы координат:
OXYZ – неподвижная, абсолютная, ось OZ вертикальна, горизонтальная плоскость OXY расположена в плоскости верхнего основания платформы;
СXPYPZP – подвижная, жёстко связанная с платформой, C – фиксированная точка верхнего основания платформы, ось СZP вертикальна;
Exyz – жёстко связана со схватом, вращается вокруг оси С5z;
С1X1Y1Z1 – вращается вокруг вертикальной оси С1Z1 вместе со звеном 1, ось С1X1 направлена от оси С1Z1 к проекции точки C2 на плоскость OXY.
Обозначим:
h – расстояние между точками С и С1; a1 – расстояние между точкой C1 и проекцией точки C2 на горизонтальную плоскость; d1 – расстояние между точкой C1 и проекцией точки C2 на вертикальную ось С1Z1, a2 = |C2C3|, a3 = |C3C4|, a4 = |C4C5|, a5 = |C5C'5|, a6 = |C'5E| .
Положение платформы и руки манипулятора определим с помощью обобщённых координат XC,YC, θ, θi, i = 1 − 5 , XC, YC - координаты точки C относительно осей OX, OY, θ - угол между осями OX и СXP, θ1 - угол между осями СXP и С1X1, θ2 − угол от оси С1Z1 до вектора C2C3, θ3 – угол от вектора C2C3 до вектора C3C4, θ4 – угол от вектора C3C4 до вектора C4C5 , θ5 – угол поворота звена 5 вокруг оси С5z.
За положительное направление отсчёта углов в координатных плоскостях OXY, С1X1Z1, Exy примем направление против часовой стрелки, если смотреть из положительного направления осей OZ, С1Y1, С5z соответственно [26]. Числовые значения углов должны удовлетворять [22] техническим ограничениям:
-1690 ≤ θ1 ≤ 1690, − 650 ≤ θ2 ≤ 900 , − 1510 ≤ θ3 ≤ 1460, − 102,50 ≤ θ4 ≤ 102,50, − 167,50 ≤ θ5 ≤167,50.
При решении прямой задачи кинематики считаются известными обобщённые координаты XC,YC, θ, θi, i = 1 − 5. Требуется определить параметры локации рабочего органа, за которые примем координаты XE,YE, ZE точки E относительно системы координат OXYZ и углы Эйлера ψ, θE, φ [26], задающие ориентацию базиса системы координат Exyz относительно базиса системы OXYZ. Интервалы изменения этих углов определяются неравенствами
0 ≤ ψ < 2 π , 0 ≤ θE ≤ π , 0 ≤ φ < 2 π.
Из рис. 16 видно
XE= XC + h cosθ + cos( θ + θ1 ) ( a1 + a2 sinθ2 + a3 sin( θ2 + θ3 ) + ( a4 + a5 + a6 ) sin( θ2 + θ3 + θ4 ) ) ,
YE = YC + h sinθ + sin( θ + θ1 ) ( a1 + a2 sinθ2 + a3 sin( θ2 + θ3 ) + ( a4 + a5 + a6 ) sin( θ2 + θ3 + θ4 ) ) ,
ZE = d1 + a2 cosθ2 + a3 cos( θ2 + θ3 ) + ( a4 + a5 + a6 ) cos( θ2 + θ3 + θ4 )
Углы Эйлера и обобщённые координаты связаны соотношениями:
если 0 < θ2 + θ3 + θ4 < π , то
ψ = θ + θ1 + π /2 ,
θE = θ2 + θ3 + θ4,
φ = θ5 − π /2 ;
если 0 < θ2 + θ3 + θ4 >π , то
ψ =θ + θ1 − π /2 ,
θE = 2 π − ( θ2 + θ3 + θ4 ) ,
φ = θ5 + π /2 ;
если -π < θ2 + θ3 + θ4 < 0 , то
ψ =θ + θ1 − π /2 ,
θE = − ( θ2 + θ3 + θ4 ) ,
φ = θ5 + π /2 ;
если θ2 + θ3 + θ4 < -π , то
ψ =θ + θ1 + π /2 ,
θE = 2 π + θ2 + θ3 + θ4 ,
φ = θ5 − π /2 ;
если θ2 + θ3 + θ4 = 0 , то
ψ + φ = θ + θ1 + θ5 ,
θE = 0;
если θ2 + θ3 + θ4 =±π, то
ψ + φ = θ + θ1 − θ5 + π,
θE = π;
В несингулярных случаях, оси Ez и OZ непараллельны, т. е. θ2 + θ3 + θ4 ≠ 0 и θ2 + θ3 + θ4 ≠ ±π, каждый из углов ψ, θE, φ может быть вычислен по представленным выше формулам. При этом выполняются условия:
ψ ∈ [ − 336,5 π /180 , 3 π + 336,5 π /180 ] , θE ∈ ( 0, π ) , φ ∈ [ − 257,5 π / 180 , 257,5 π /180 ].
Если же оси Ez и OZ параллельны (сингулярные случаи), т. е. θ2 + θ3 + θ4 = 0 или θ2 + θ3 + θ4 = ±π, то возможно нахождение только суммы
ψ + φ ∈ [ − 156,5 π /180 , 876,5 π /180 ].
Реализовать вычисления по приведённым соотношениям можно, обратившись к представленной на следующей странице функции FK на языке системы Mathematica [27].