Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Механика промышленных роботов и манипуляторов с электроприводом

..pdf
Скачиваний:
19
Добавлен:
12.11.2023
Размер:
9.74 Mб
Скачать

Следовательно, их можно заменить одним поворотом на угол

+

*Тогда

 

т 3 “

1 1 Z 2 '

 

 

 

Размеры звеньев в связи с указанным выбором связанны* систем задаются векторами р

"О, р 2 "

[/2»0. 0]Т, р з - [/3 ,0 , О]7 . Выберем на схвате некоторую точку 5 (ею, напримео,

может быть центр масс объекта манипулирования), определенную вектором ^

3 ) - (лс^З),

У'§'.

, и исследуем ее движение. Зависимости <4.55) —(4.58) позволяют

составить

алгоритм для звеньев 7—3.

 

 

 

Операторы угловой скорости и углового ускорения можно найти согласно формулам (4.57). Однако вычисления будут проще, если поступить следующим образом. Отыскать векторы

абсолютной угловой скорости S/ и абсолютного углового ускорения ? /(7- 1 , 2 ,3) по формулам (4,58) и из их компонент образовать операторы А(0 /) и А(У t ):

ш1 ^1 ’ ш2 "

ш\ + ^ 2 ’ ^3

" Ш2 + ^3 ’

где

 

 

S5J “ T| SÎ J , Q2 —^ 2 ^ 2

* ^ з * т 3 Q3 ,

z , ^ 2 **^2 ^ 2 ^ ^ 3 3 ^ 3 ^

Аналогичным образом получаются формулы для абсолютных угловых ускорений звеньев:

е 2 ■ £ | + Т?2 + А( 0 | ) П 2 *

Г 3 “ f 2 + ^3 + ^

> « 3 ’

 

где ^

^

^

 

 

ÈJ - т J

, Е2- т 2 ^ 2 . 1 з -

т 3 È3 ,

ЛГ7?, Я2- ?2XTz%£ 3- р 3* 1 г.

Для скоростей точек Ai, А2 , A3 , А$ имеют место формулы, вытекающие из выражения (4.55):

7Л1

"

7Аг

" ° ’

ГА3

~

Я ( " 2 ) т 2 ? 2 '

K s

-

%

+ *<S3>*3'5S

Ускорения тех же точек подсчитываются по формулам

РА\ - ?Л2 - 0 ,

7A3 -

а

7A S

А ( г 2

Од

+

+ т\ ) т 2 ^ 2

Д (Г 3 + ® 3 ) T 3^S

4.9.Обратная задача кинематики манипулятора

Приведенные выше рекуррентные соотношения (4.46), (4.55) — (4.58) позволяют определить положение, скорость и ускорения любых точек и звеньев манипулятора, если заданы обобщенные-координаты и их произ­ водные по времени. Рассмотрим обратную задачу: определим обобщенные координаты, их скорость и ускорение, если задано движение исполнитель­ ного звена манипулятора — кисти, несущей захватное устройство или инструмент.

В общем случае положение исполнительного звена задается шестью

эйлеровыми координатами: положением полюса, в качестве которого мо­ жет быть задан центр масс захватного устройства вместе с объектом мани­ пулирования, и тремя углами Эйлера, определяющими ориентацию захватного устройства связанной системы в неподвижном пространстве. Возникает шестимерная обратная задача по определению шести обобщен­ ных координат. Целесообразно максимально снизить размерность задачи. Для этого следует выяснить движение точек, максимально близко располо­ женных к началу кинематической цепи. Такая возможность обычно суще­ ствует благодаря тому, что присоединение кисти робота к несущей руке можно смоделировать сферической кинематической парой. Тогда, рассмат­ ривая кисть как твердое тело, осуществляющее предписанное пространст­ венное движение, по известным формулам кинематики определим движение центра сферической пары. Задача стала трехмерной: требуется определить три обобщенные координаты для трехзвенной кинематической цепи позиционирующих движений по заданным координатам исполни­ тельной точки. Для этого необходимо решить векторное уравнение

f

( Q

Q 2 * ^ 3 ) ” 7 Л

r A( V\ » ^2 ’

^ = 0 »

(4.59 )

где

 

 

 

 

 

/

-

[ / , < 5 ) , / 2<?>,

/ 3<?>],

 

 

гл 0 — заданный радиус-вектор точки; f A — вектор-функция положения точки, найденная из решения прямой задачи.

Для численного решения уравнения (4.59) можно применить метод Ньютона. Предположим, что найденор-е приближение искомых координат q Pss [01(р\ ? 2(р)»<7з(р) I7Тогда следующее приближение можно представить

в виде

 

*<Р+1> = д(Р) + ё,

 

гдее — поправка.

 

Поправка с находится из линейной системы уравнений

 

We, = - r A(g (P h ,

(4.60)

где W — матрица Якоби системы скалярных функций/j (^), / 2<<7>*/3(^)1

df\ ! d q j W = d f 2 / d q x d f 3 / d q x

d f {/ d q 2

d f x/ d q 3

d f 2 / d q 2

d f 2 / d q 3

d f 3 l d q 2

d f 3l d q 3

Матрица Якоби вычисляется при значениях qv соответствующих р-му приближению. При исследовании движения на границах рабочей зоны мо­ гут возникнуть затруднения в связи с вырождением матрицы W, а значит, и невозможностью решения линейного уравнения (4.60).

Для ортогональных манипуляторов существует возможность составле­ ния алгоритма решения обратной задачи, который более эффективен, чем

юз

метод последовательных приближений. Покажем это на примере шарнир­ ного манипулятора (рис. 4.23). Вданном случае система координат отлича­ ется от используемой для манипулятора, схема которого представлена на рис. 4.22.

Считаем, что известны параметры движения исполнительной точки 5, а также задан закон изменения ориентации рабочего органа в осях х, у, z. Это имеет место, например, для манипулятора, ведущего сварку по простран­ ственной траектории. Задача состоит в обеспечении перемещения свароч­ ной горелки по заданной траектории с заданной скоростью и в сохранении требуемого наклона горелки к поверхности сварки. В других случаях ори­ ентация рабочего органа может оставаться неизменной, но для этого ею все равно необходимо управлять.

Итак, определены движение точки S и сферическое движение связанной системыжоординат рабочего органа с началом в точке S. Воспользовавшись формулами кинематики для общего случая движения твердого тела, приведенными в§ 4.4, найдем движение точки Л4-центра сферического шарнира, соединяющего руку манипулятора с кистью:

= г с

T 4 r A S '

Что касается определения обобщенных координат механизма ориента­ ции, эта задача решается на основании приведенных в § 4.3 соотношений между направляющими косинусами и эйлеровыми углами.

Решение задачи ориентации обеспечивает переход к задаче позициони­ рования точки Л4: по заданному движению этой точки следует определить угловые кооординаты х, ф2, $, их скорости и ускорения.

Рассмотрим проекцию манипулятора на плоскость ху (рис.4.24). Из

Рис. 4.23. К обратной задаче кинематики шарнирного манипулятора

Рис. 4.24. Схема для определения угла j (правосторонняя схема манипулятора)

простых геометрических соображений следует зависимость

Ч>\

= <Р0 + /3

=

a rc t g( y / x ) + a r c t g ( / i 3/ p 3) ,

(4.61)

где y,

x обозначают координаты точки А по осям у, х\

углы (р0 и /3, а

также отрезки Л3 и р3 указаны н£ рисунке.

 

Очевидно, что

 

 

 

р 3 = \/JC2 +

у 2

Л2

 

Дифференцируя зависимость (4.61) по времени, находим выражение

для определения скорости по координате<р{:

 

<р{ = ( х 2 + y 2 ) ~ l (xy

ух)

h 3( x 2 + у2

h 3) "3 / 2 (JCX + уу).

(4.62) Рассматриваемую схему манипулятора можно назвать правосторонней.

Ее зеркальным отображением относительно главной плоскости манипу­ лятора, проходящей через оси jt(1 \ z(1), является левосторонняя схема. Для левосторонней схемы формула (4. 61) отличается лишь знаком /3:

1 - <Р0 &

Для центральной схемы (А3 = 0) формулы (4.61) и (4.62) значительно упрощаются.

Последующим дифференцированием формулы (4.62) можно опреде­ лить угловое ускорение по координате(р^.

Для определения движения по координатам и следует рассмотреть движение эквивалентной кинематической цепи АВС в плоскости

На основании уравнений преобразования координат для плоскости и условий ортогональности (4.8) составляется система уравнений:

X

II

н

+

/ 2

=

2 .

‘ 2

 

ХВ

+

1 Ъа 1 1 ' ZA ZB + 1 3й

ZB > 1 = а \ х + а \ , .

Здесь верхний индекс, указывающий на то, что проекции берутся на осях Л(1) и z(1), опущен. Параметры /3и /2определяют длины звеньев 3 и 2, кроме

того, введены обозначения aj j= cosa, Û2I= sina, причем a = ^ 3 + ^>2 ~ 2ir. Система (4.63) нелинейная, соответствующими преобразованиями ее

можно свести к квадратному уравнению, решение которого имеет вид

•ч

а 2 1 = I t z A ± ^ А

( х А + ZA ) i { 2

ф 1 / < ^ + г 1 > »

< = ( x ^ + z ^ + / 3

/ | ) / ( 2 / 3).

 

(4.64) Два знака перед радикалом указывают на наличие двух решений, соот­

ветствующих двум возможным конфигурациям кинематической цепи, обычно не переходящим одна в другую (их переход возможен только при особом положении манипулятора, когда звенья 2 и 3 вытянуты в одну линию). Для конфигурации, показанной на рис. 4.23, в уравнении (4.64) следует принять знак минус перед радикалом.

Для определения скоростей продифференцируем систему уравнений (4.63) по времени, получим линейную систему, в векторно-матричной за­ писи имеющую вид

V х = Avyi

(4.65)

где vx = [хл , i A, 0,0 F , vy= [ûibû2b*z?. 2в F ,

 

h

0

1

0 '

0

0

1

0

0

х в

г в

J

а и .

а 2 1 0

0

Из уравнения (4.65) следует:

 

v y = A~1V x9

(4.66)

где A~l — обратная матрица.

 

Вычисление обратной матрицы производится по правилам, приведен­ ным в курсах линейной алгебры. Для данного случая матрица А'1 имеет вид

/3(а11ZB

а 2\ ХВ^

Л

 

 

 

~а2 \ х В

~а2 1

а21

/Зг.

а М х В

а \ \ 2В

-а11

13х

1за п г в

l 3a2

\ z B

~13а21

 

- 13а \ 1х в

~1За 2

\ ХВ

13а 11

Ф.

Дифференцируя уравнение (4.66) еще один раз, находим ускорения

звеньев. Удобнее, однако, продифференцировать исходную систему (4.65), а затем произвести ее обращение. При дифференцировании получим

V x = у + А 9у .

Несложными преобразованиями это уравнение приводится к виду

где

**х

=

ZA*

+

» " ^ п + à \ \ ) ] T

а у

~

^ 2 1 » *2?»

 

 

Полученные соотношения позволяют найти параметры дп , а2х, хв,гв и их производные. Исходя из того, что aj j =cosa, tg<P2 =ZB ^XB ^и принимая во внимание геометрические соотношения между углами, нетрудно опреде­ лить <р2 и <р$ и их производные.

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

Ях =

h x = z,

 

я 2

Л2 “ У'

 

<?з

=

Л з = *

V

с цилиндрической системой

<7i =<Pi = arctg(*/z),

Q 2 = h 2 =y,________

qz=hz =Jx 2+ z2 ;

со сферической системой

^ i = ^ P l = a r c t g ( * / z ) ,

Q2 = $2 = a r c t g ( y / / * 2 + z 2) ,

Q3 = h 3 = V x 2 + y 2 + z 2

Для отыскания связи между скоростями изменения координат <7р <?2» и параметрами движения исполнительной точки 5 в пространстве xyz сле­ дует продифференцировать приведенные выше соотношения.

Представленные на рис. 1.8 и 1.9 схемы манипуляторов являются цент­ ральными. В том случае, если анализируются манипуляторы с левосторон­ ними или правосторонними схемами, при определении угла <Pi и скорости его изменения следует использовать формулы (4.61) и (4.62).

4.10. Кинематика дифференциальных передаточных механизмов

Манипуляторы с электроприводом, размещенным на основании, имеют сложные схемы передачи движения. Во многих случаях они представляют системы дифференциальных и планетарных механизмов с коническими и цилиндрическими зубчатыми колесами, соединенными коаксиальными трансмиссионными валами. Дифференциальные механизмы используются практически во всех механизмах передачи движения на кисть робота. Ус­ тановим основные особенности кинематического расчета таких передач.

Допустим, что рассматривается дифференциальный механизм с тремя входами и одним выходом. Угол поворота выходного вала представляет

некоторую функцию углов поворота входных валов р

^ 3:

* о ы х = / < * 1 ’ V2, v3)

(4 .6 7 )

Угловая скорость на выходе найдется дифференцированием выражения

(4.67) по времени:

 

 

 

 

3 /

d*Pj

£ /

dtp2

а /

d<f з

д(р j

d t

дф2

dt

д<р3

dt

 

или

 

 

 

 

 

 

=

к хфх +

к 2ф2 +

к 3ф2

 

(4.68)

где к1, к2,

— частные производные от функции (4. 67).

Для определения, например, кj можно рассмотреть частный случай, когда ф2= Фз = 0- Тогда, очевидно,

к,

= ф

/ ф .

1

г в ы х

^ 1

Аналогичным образом выражаются и остальные частные производные. Таким образом, коэффициенты кх, к2, къ в формуле (4.68) представляют передаточные отношения от выходного вала к входному в механизме, в котором остальные входные валы неподвижны. Будем называть их переда­ точными коэффициентами. Степень подвижности таких частных механиз­ мов равна единице. По структуре это либо замкнутые дифференциалы, планетарные механизмы, либо рядовые передачи. Расчет таких механиз­ мов производится известными методами.

В механизмах привода кисти робота из-за габаритных ограничений раз­ местить двигатели непосредственно в кинематических парах невозможно. На рис. 4.25 представлена схема привода механизма ориентации с провод­ кой движения от двигателей Ml, М2, М3, расположенных на предплечьи манипулятора. Механизм обеспечивает движение по трем эйлеровым уг­ лам.

Угловые скорости коаксиальных валов обозначим ы2, w3>скорости относительных вращений в шарнирах кисти — qx, q2, à$. Установим связь между этими скоростями. Очевидно, что qx = о) j . Для определения q2 рас­

смотрим соответствующий двухстепен­ ной дифференциал, включающий води­ ло Я 7 и шестерни 7 и 2. Воспользуемся формулой (2. 1), тогда при соответствую­ щих обозначениях будем иметь

Яг

=

(W2

- ,Wl ) Z l / z 2

 

 

 

Для определения <у3 следует рассмот­

 

 

реть трехстепенной дифференциал с тре­

 

 

мя входами от Ml, М2, М3. Для этого

 

 

механизма

 

 

 

 

 

 

ЯЪ

 

^ \ы \

+ ^2W2 +

^3^3*

 

 

 

Чтобы

найти

рассмотрим меха­

 

 

низм, полученный из исходного диффе­

 

 

ренциала при о)2 ~ 0,

= 0. Это

 

 

 

комбинированная система, в состав ко­

 

 

торой входят два коаксиальных плане­

 

 

тарных

механизма

с коническими

 

 

 

колесами (один образован шестернями 1

 

 

и 2 и водилом Я7, второй — шестернями

 

 

 

3 и 4 и водилом Н1) и дифференциаль­

 

 

 

ный механизм с цилиндрическими коле­

 

 

 

сами (водило Н2 и шестерни 5 и б). Из

 

 

 

расчета планетарных механизмов имеем

 

U

 

Ы2Н1

=

 

 

 

 

 

 

 

 

 

( Г

W4 ///

= (Z3/ Z l K

 

 

гдеи 2нI и о)4//у — угловые скорости звена

 

2 и 4 относительно водила Н1.

Рис. 4.25. Механизм ориентации кисти с

Проведя соответствующие вычисле­

коническими дифференциалами

ния, получим

 

 

 

 

 

 

к

. Î A _ £ т£5 (£ з

£ i)

 

 

 

 

Ш1

^ 6

Z 4

Z 2

 

 

 

Для

определения к2 предполагаем,

что

= о>3 = 0.

В этом случае

приведенный дифференциал можно рассматривать как двухступенчатый механизм, первая ступень которого содержит рядовую коническую переда­ чу 7—2, а вторая — планетарный механизм, состоящий из водила Я2, неподвижного колеса 5 и сателлита 6. Выполнив кинематический расчет, получим

к

= Ь

= Z 7Z 5Z l

2

U2

Z bZ bZ 2 '

Для определения Л3 полагаем, чтоu j = 0, и>2 т0. Тогда из рассмотрения системы рядовых зубчатых передач следует:

,

_ Я з _ Z3Z5Z7

Л3

ш 3 - “z4 z6 z8 *

Итак, для механизма, показанного на рис. 4.25, установлена зависи­ мость

 

 

Z7 Z5

+

Z7 Z5 Z2 ы

Z3 Z5 Z7

h

=

 

------------

 

Z 8 Z6

 

Z 8 Z 6 Z 1

Z4 Z6 Z8

 

 

 

Скорость выходного вала есть линейная форма с постоянными коэффици­ ентами от скоростей входных звеньев. Эту связь можно представить с по­ мощью некоторой квадратной матрицы U передаточных коэффициентов:

? =

Uû.

 

 

(4.69)

Для рассматриваемого механизма

 

Ъ = [ ? , . Я2 > Я3 ]Т о = [ « , , ш2 , w3 ]r .

 

1

 

 

о

 

~г \ !

2" i

2 1 / Z2

 

U=

Z7 Z5 ( Z3

z l }

Z7 Z5 Z2

Z7 Z5 Z3

 

Z 8 Z 6 z 4

Z 2

Z 8 Z6 Z I

Z8 Z6 Z4

Если принять z, = z2, z3 = z4, z5 = z6, z7 = z8, матрица передаточных коэффициентов приобретает частный вид:

 

 

1

0

0

1

U

=

-1

1

 

 

0

1

- ? ]

 

 

т.е.

 

w 1, <?2=

 

1 +w2»^3 = а)2 " w3*

Используя уравнение (4.69), можно найти скорости двигателей, обеспе­

чивающие заданное поворотное движение кисти:

 

ü = U ' xq.

(4.70)

Подобным же образом производится анализ дифференциальных схем передачи движения звеньям цепи позиционирования для манипуляторов типа, представленного на рис. 3.3. Характерной особенностью всех схем передачи движения посредством зубчатых колес является постоянство мат­ рицы передаточных коэффициентов, ее независимость от конфигурации манипулятора. Вторая особенность — кинематическая взаимосвязь движе­ ний отдельных звеньев, что усложняет управление манипулятором и явля­ ется существенным недостатком таких схем.

Если приводы звеньев разместить непосредственно в кинематических парах, то

по

Соседние файлы в папке книги