Robotica Control Deteccion Vision e Inteligencia
Robotica Control Deteccion Vision e Inteligencia
Robotica Control Deteccion Vision e Inteligencia
.)
() 11
./
'
K. S. FU
R. C. GONZALEZ
C. S. G. LEE
Traduccin
SEBASTIAN DORMIDO BENCOMO
Catedrtico lnfonntica y Automtica
Facultad de Ciencias Fsicas
UNED (Madrid)
Re,isin tcnica
ANTONIO VAQUERO SANCHEZ
Catedrtico Informtica y Automtica
Facultad de Ciencias Fsicas
Universidad Complutense de Madrid
.. RO DE CORTESlJ.
11
McGraw-Hill
MADRID BOGOTA BUENOS AIRES GUATEMALA LISBOA MEXICO
NUEVA YORK PANAMA SAN JUAN SANTIAGO SAO PAULO
AUCKLAND HAMBURGO LONDRES MONTREAL
NUEVA DELHI PARIS SAN FRANCISCO SINGAPUR
ST. LOUIS SIDNEY TOKIO TORONTO
CONTENIDO
Prlogo
l.
Introduccin
l.l.
1.2.
1.3.
1.4.
1.5.
1.6.
1.7.
1.8.
2.
3.
xi
Antecedentes
1
Desarrollo histrico
4
Cinemtica y dinmica del brazo del robot
6
Planificacin de la trayectoria y control del movimiento del manipulador 7
Sensores del robot
9
Lenguajes de programacin de robots
10
Inteligencia del robot
11
11
Referencias
13
2.1.
2.2.
2.3.
2.4.
13
Introduccin
El problema cinemtico directo
El problema cinemtico inverso
Observaciones finales
Referencias
Problemas
Introduccin
Formulacin de Lagrange-Euler
Formulacin de Newton-Euler
Ecuaciones de movimiento generalizadas de d'Alernbert
Observaciones finales
Referencias
Problemas
15
54
78
79
79
85
128
85
146
146
147
87
106
viii
CONTENIDO
Introduccin
Consideraciones generales sobre la planificacin de trayectoria
Trayectorias de articulacin interpoladas
Planificacin de trayectorias de caminos
cartesianos del manipulador
Observaciones finales
Referencias
Problemas
152
152
155
157
178
201
201
203
206
Introduccin
206
5.2.
5.3.
5.4.
5.5.
5.6.
5.7.
5.8.
Control adaptativo
5.9.
Observaciones finales
208
210
228
232
234
239
251
272
Referencias
Problemas
273
274
6.
Deteccin
275
6.1.
Introduccin
275
CONTENIDO
8.4.
8.5.
8.6.
8.7.
Segmentacin y descripcin
Reconocimiento
Interpretacin
Observaciones finales
Referencias
Problemas
de estructuras
tridimensionales
431
440
456
460
460
462
465
465
467
478
487
489
489
491
Introduccin
Bsqueda del espacio de estados
!0.3.
Reduccin del problema
10.4. Uso de la lgica de predicados
!0.5. Anlisis rneans-ends
!0.6.
Resolucin del problema
!0.7.
Aprendizaje del robot
10.8.
Planificacin de tareas de robot
10.9.
Problemas bsicos en la planificacin de tareas
!0.10. Sistemas expertos e ingeniera del conocimiento
!O.! l. Observaciones finales
Referencias
Problemas
540
540
Apndices
542
10.2.
A Vectores y matrices
B Jacobiano del manipulador
Bibliografia
Indice
491
492
502
507
512
515
523
525
528
536
539
542
564
576
591
CAPITULO
UNO
INTRODUCCION
1.1
ANTECEDENTES.
1 ,
una herramienta. Se disea para alcanzar una pieza de trabajo localizada dentro
de su volumen de trabajo. El volumen de trabajo es la esfera de influencia de un
robot cuyo brazo puede colocar el submontaje de la mueca en cualquier punto
dentro de la esfera. El brazo generalmente se puede mover con tres grados de
libertad. La combinacin de los movimientos posiciona a la mueca sobre la
pieza de trabajo. La mueca normalmente consta. de tres movimientos giratorios.
La combinacin de estos movimientos orienta a la pieza de acuerdo a la configuracin del objeto para facilitar su recogida. Estos tres ltimos movimientos se
denominan a menudo elevacin (pitch), desviacin (yaw) y giro (roll). Por tanto,
para un robot con seis articulaciones, el brazo es el mecanismo de posicionamiento, mientras que la mueca es el mecanismo de orientacin. Estos conceptos se
ilustran para el robot Cincinnati Milacron T3 y el robot PUMA de Unimation
que se muestran en la figura 1.1.
Barrido
del brazo
G iro
a)
'
,
17 0i~
\
/
Doblado de mueca
200'
Montaje de la pinza
b)
Figura 1.1
INTRODUCCION
Cilndrico
De revolucin
Cartesiano o xy:
Esfrico
INTRODUCCION
RODOTICA:
CONTROL,
DETECCION,
VISJO:-,. E INTELIGENCIA
Heavy Industries negoci una licencia con Unimation para sus robots. Uno de
los desarrollos ms poco usuales en robots sucedi en 1969, cuando se desarroll
un camin experimental por la General Electric para la Armada Americana. En
el mismo ao se desarroll el brazo Boston y al ao siguiente el brazo Stanford,
que estaba equipado con una cmara y controlado por computadora. Algunos de
los trabajos ms serios en robtica comenzaron cuando estos brazos se utilizaron
como robots manipuladores. Un experimento en el brazo Stanford consista en
apilar automticamente bloques de acuerdo con diversas estrategias. Esto era un
trabajo muy sofisticado para un robot automatizado de esa poca. En 1974,
Cincinnati Milacron introdujo su primer robot industrial controlado por computadora. Lo llam The Tomorrow Tool (la herramienta del maana) o T3, que
poda levantar ms de 100 libras as como seguir a objetos mviles en una linea
de montaje.
Durante los aos setenta se centr un gran esfuerzo de investigacin sobre el
uso de sensores externos para facilitar las operaciones manipulativas. En Stanford, Bolles y Paul [ 1973], utilizando realimentacin tanto visual como de fuerza,
demostraron que un brazo Stanford controlado por computadora, conectado a
una PDP-10, efectuaba el montaje de bombas de agua de automvil. Hacia la
misma poca, Will y Grossman [1975] en IBM desarrollaron un manipulador
controlado por computadora con sensores de contacto y fuerza para realizar
montajes mecnicos en una mquina de escribir de veinte piezas. Inoue [1974],
en el Artificial Intelligence Laboratory del MIT, trabaj sobre los aspectos de
inteligencia artificial de la realimentacin de fuerzas. Se utiliz una tcnica de
bsqueda de aterrizajes. propia ~ la navegacin area, para realizar el posicionado inicial de una tarea de montaje precisa. En el Draper Laboratory, Nevins y
colaboradores [1974] investigaron tcnicas sensoriales basadas en el control
coordinado de fuerza y posicin. Este trabajo desarroll la instrumentacin de un
dispositivo remote center compliance (RCC) (centro remoto de control coordinado de fuerza y posicin) que se uni a la placa de montaje de la ltima articulacin del manipulador para cerrar el conjunto de coincidencias de piezas. Bejczy
[ 1974], en el Jet Propulsion Laboratory, desarroll una tcnica de control de par
basada en computadora sobre su brazo Stanford ampliado para proyectos de
exploracin espacial. Desde entonces han sido propuestos diversos mtodos para
manipuladores mecnicos.
Hoy da vemos la robtica como un campo de trabajo mucho ms amplio
que el que tenamos simplemente hace unos pocos aos, tratando con investigacin y desarrollo en una serie de reas interdisciplinarias, que incluyen cinemtica, dinmica, planificacin de sistemas, control, sensores, lenguajes de programacin e inteligencia de mquina. Estos tpicos, introducidos brevemente en las
secciones que siguen, constituyen el ncleo del material de este libro.
(i
)
,/
"'. ,- La cinemtica del brazo del robot trata con el estudio analtico de la geometra
' 1
~ovimiento de un brazo de robot con respecto a un sistema de coor enaaas
lNTRODuCCIOS
de referencia fijo sin C.Q!!_siderar las fuerzas. o momentos que originan el movimiento. A.fil, la ~nemtica se interesa por la descripcin analtica del desplazamiento espacial del robot comoua funcin del tiempo, en particular de las
relaciones en re la posicin aelasvariables de articulacin y la posiciny
oentac1aeJ-efecror ma -del razo ael rooor.. -Hay dos pro6 emas fundamentales en la cinemtica del robot. El primer.
problema se suele conocer c~mo el problema cinemtico directo, mientras que el
seg_ulao es e --rob!eiacinemtico inverso. Como las variables indepenciien es en
ti"n robot son las variables"' e articulacin, y una tarea se suele dar en trminos
del sistema de coordenadas de referencia, se utiliza de manera ms frecuente el
problema cinemtico inverso. Denavit y Hartenberg [1955] propusieron un enfoque sistemtico y generalizado de utilizar lgebra matricial para describir y
representar la geometra espacial de los elementos del brazo del robot con
respecto a un sistema de referencia fijo. Este mtodo utiliza una matriz de
transformacin homognea 4 x 4 para describir la relacin espacial entre dos
elementos mecnicos rgidos adyacentes y reduce el problema cinemtico directo
a encontrar una matriz de transformacin homognea 4 x 4 que relaciona el
desplazamiento espacial del sistema de coordenadas de la mano al sistema de
coordenadas de referencia. Estas matrices de transformacin homogneas son
tambin tiles en derivar las ecuaciones dinmicas de movimiento del brazo del
robot. En general, el problema cinemtico inverso se puede resolver mediante
algunas tcnicas. Los mtodos utilizados ms comnmente son el algebraico
matricial, iterativo, o geomtrico. En el captulo 2 se da un tratamiento detallado
{ de los problemas cinemticos directo e inverso.
'
- La dinmica del robot, por otra 2arte, trata con la formulacin matemtica
I
l de las ecuaciones del movirmento del brazo. Las ecuaciones c!inmicas de movI='
miento de un manipula or son un conjunto de ecuaciones matemticas que
descrben la conducta dinmica del manipulador. Tales ecuaciones de movimienson ut1les para simulacin en computadora ae'J movimiento del brazo, el diseo
de ecuaciones de control apropiadas para el robot y la evaluacin del diseo y
estructura cinemtica del robot. El modelo dinmico real de un brazo se puede
obtener de leyes fisicas conocidas tales como las leyes de Newton y la mecnica
Iagrangiana. Esto conduce al desarrollo de las ecuaciones dinmicas de movimiento para las distintas articulaciones del manipulador en trminos de los
parmetros geomtricos e inerciales especificados para los distintos elementos. Se
pueden aplicar sistemticamente enfoques convencionales como las formulaciones de Lagrange-Euler y de Newton-Euler para desarrollar las ecuaciones de
movimientos del robot. En el capitulo 3 se presenta una discusin detallada de la
dinmica del brazo del robot. ~
to
1 -1.4
cumplir una tarea deseada controlando al manipulador para que siga un camino
previsto. Antes de mover el brazo, es de inters saber si hay algn obstculo
presente en la trayectoria que el robot tiene que atravesar (ligaduras de obstculos) y si la mano del manipulador necesita viajar a Jo largo de una trayectoria
especificada (ligaduras de trayectoria). El problema de control de un manipulador
se puede dividir convenientemente en dos subproblemas coherentes: el subproblema de planificacin de movimiento (o trayectoria) y el subproblema de control
del movimiento.
La curva espacial que la mano del manipulador sigue desde una localizacin
inicial (posicin y orientacin) hasta una final se llama la trayectoria o camino. La
planificacin de la trayectoria (o planificador de trayectoria) interpola y/o aproxima la trayectoria deseada por una clase de funciones polinomiales y genera una
secuencia de puntos de consignas de control en funcin del tiempo para el
control del manipulador desde la posicin inicial hasta el destino. El captulo 4
discute los distintos sistemas de planificacin de trayectoria para movimientos
libres de obstculo, as como el formalismo para describir el movimiento deseado
del manipulador en trminos de la secuencia de puntos en el espacio a travs de
los cuales debe pasar y la curva espacial que recorre.
CCN~L
En general, el problema de control de movimientos consiste en: 1) obtener los
modelos dinmicos del manipulador, 2) utilizar estos modelos para determinar
leyes o estrategias de control para conseguir la respuesta y el funcionamiento del
sistema deseado. Como la primera parte del problema de control se discute
extensivamente en el captulo 3, el captulo 5 se concentra sobre la segunda parte
del problema de control. Desde el punto de vista de anlisis de control, el
movimiento del brazo de un robot se suele realizar en dos fases de control
distintas. La primera es el control de movimiento de aproximacin en el cual el
brazo se mueve desde una posicin/orientacin inicial hasta la vecindad de la
posicin/orientacin del destino deseado a lo largo de una trayectoria planificada. El segundo es el control del movimiento fino en el cual el efector final del
brazo interacciona dinmicamente con el objeto utilizando informacin obtenida
a travs de la realimentacin sensorial para completar la tarea.
Los enfoques industriales actuales para controlar el brazo del robot tratan
cada articulacin del brazo como un servomecanismo de articulacin simple.
Este planteamiento modela la dinmica de un manipulador de forma inadecuada
porque desprecia el movimiento y la configuracin del mecanismo del brazo de
forma global. Estos cambios en los parmetros del sistema controlado algunas
veces son bastante significativos para hacer ineficaces las estrategias de control
por realimentacin convencionales. El resultado de ello es una velocidad de
respuesta y un amortiguamiento del servo reducido, limitando as la precisin y
velocidad del efector final y hacindolo apropiado solamente para limitadas
tareas de precisin. Los manipuladores controlados de esta forma se mueven a
velocidades lentas con vibraciones innecesarias. Cualquier ganancia significativa
en el rendimiento en esta y otras reas del control del brazo del robot requieren
la consideracin de modelos dinmicos ms eficientes, enfoques de control
sofisticados y el uso de arquitecturas de computadoras dedicadas y tcnicas
de procesamiento en paralelo. El captulo 5 est enfocado a deducir las leyes
INTRODUCCION
i.s
1.6
LENGUAJES DE PROGRAMACION
DE ROBOTS
INTRODUCCION
t.7 INTELIGENCIA
11
DEL ROBOT
1.8
REFERENCIAS
las referencias generales citadas a continuacin son representativas de publicaciones que tratan con temas de inters en robtica y en campos relacionados. Las
referencias que se dan al final de los captulos que siguen se clasifican segn los
ternas especficos que se estudian en el texto. La bibliografa al final del libro se
Organiza en orden alfabtico por autores, y contiene toda la informacin pertinente para cada referencia citada en el texto.
Algunas de las mejores revistas y actas de conferencia que de forma rutinaria
contienen artculos sobre los diversos aspectos de la robtica incluyen:
IEEE Journal o( Robotics and Automation: International Journal of Robotics
Research;
DOS
/
2.1
INTRODUCCION
2.
C
A
P
I
T
U
L
O
Angules de
Cinemtica
directa
las articulaciones
-----.,.
q,(1) .....
q,(1)
---,...
Posicin y orientacin
del efector final
Parmetros
de los
elementos
Angulos de
las articulaciones
q,(I) .....
q,(1)
Figura 2.1
Cinemtica
inversa
~-----...J
CINEMATICA
DEL BRAZO
DEL ROBOT
15
z
p
2.2
otz-"1=
Se utiliza
lgebra"
vectorial y matricial para desarrollar un mtodo
generalizado y sistemtico para describir y representar la localizacin de los
elementos de un brazo con respecto a un sistema de referencia fijo. Como los
elementos de un brazo pueden girar y /o trasladarse con respecto a un sistema
de coordenadas de referencia, se establecer un sistema de coordenadas ligado
al cuerpo a lo largo
del eje de la articulacin para cada elemento. El problema cinemtico directo se
.reduce a_encontrar una matriz de transformcin que relaciona el sistema
de
coordenadas ligado al cuerpo al sistema de coordenadas ae referencia. Se utiliza"
i1a matriz de rotacin 3 x 3 para describir las operaciones rotacionales del
sistema ligado al cuerpo con respecto al sistema de referencia. Se utilizan entonces las coordenadas homogneas para representar vectores de posicin en un
espacio tridimensional, y las matrices de rotacin se ampliarn a matrices de
transformacin homognea 4 x 4 para incluir las operaciones traslacionales del
sistema de coordenadas ligado al cuerpo. Esta representacin matricial de un
elemento mecnico rgido para describir la geometria espacial de un brazo fue
utilizada por primera vez por Denavit y Hartenberg [1955]. La ventaja de
utilizar la representacin de elementos de Denavit-Hartenberg es su universalidad algortmica para derivar las ecuaciones cinemticas de un brazo.
2.2.1
P.i~
Una matriz de rotacin 3 x 3 se puede definir como una matriz de transformacin que opera sobre un vector de posicin en un espacio eucldeo tridimensional
y transforma sus coordenadas expresadas en un sistema de coordenadas rotado
OUVW (sistema ligado al cuerpo) a un sistema de coordenadas de referencia
OXYZ. En la figura 2.2 se dan dos sistemas de coordenadas rectangulares, uno el
Los vectores se representan en letras minsculas en negrita; las matrices en letras maysculas en
negrita.
16
ROBOTICA:
INTELIGENCIA
CONTROL,
DETECCION,
YISIO"'
(2.2-1)
Px,z
Obsrvese que fisicamente el punto Pwvw ha sido girado junto con el sistema de
coordenadas OUVW.
Recordando la definicin de las componentes de un vector tenemos
(2.2-3)
donde Px, p, y p, representan las componentes de palo largo de los ejes OX, O Y
y OZ, respectivamente, o las proyecciones de p sobre los ejes respectivos. As,
utilizando la definicin del producto escalar y la ecuacin (2.2-3),
Px
Py
P,
(2.2-4)
CINEMATICA
[ J'.y" :
]
[
~.
i, .
jy .
k, i"
k,
p,
p,
J,
J,
'. " . k, .
) y
l
,
DEL BRAZO
DEL ROBOT
l [p "
17
(2.2-5)
P v
kz' k ,
p,.
R =
[ i,
i,
)y . l."
k,
I"
i, .
J,
t.
jv
k, 'jv
i,
k,
k,
k, k ,
)y.
(2.2-6)
n ....
Puvw
p.,
p,.
j,.. ~
k,.. 1.,,
Qpxyz
, jyJ,
i.
k,. . jy '
l [;:]
Jv k,k,
i.
k., k,
(2.2-7)
(2.2-8)
Como los productos escalares son conmutativos, se puede ver de las ecuaciones (2.2-6) a (2.2-8) que
(2.2-9)
'Y
QR
(2.2-1 O)
OUVW, tendr coordenadas diferentes (p.,,, P}~ Pzf con respecto al sistema de
referencia OXYZ. La matriz de transformacin necesaria R,,. se llama la matriz
de rotacin respecto al OX con ngulo a. R,,. se puede derivar del concepto de
rnatriz de transformacin anterior, esto es,
Pxyz
R.,,. Puvw
(2.2-11)
con i, - i., y
i, . j.
-s~n
cos
J.: L
k, . j,,
e,.]
(2.2-12)
C1.
</Jl
cos
R=.
<P
-sen
cos
O
= [ se~()
ee
O (2.2-13)
1
Las matrices Rx. ,, R, . .,. y R:. 6 se llaman las matrices de rotacin bsicas.
Se pueden obtener otras matrices de rotacin finitas a partir de estas matrices.
~
'
Ejemplo: Dado dos puntos a., .... = (4, 3, 2)7 y b.v,.. = (6, 2, 4)7 con respecto
al sistema de coordenadas girado OUVW, determinar los puntos correspondientes a_.y,, bxy: con respecto al sistema de coordenadas de referencia si
ha sido rotado 60 respecto del eje OZ.
SOLUCIN:
axy:
xyz
R:.
0,500
[ 0,~66
4(0,5)
60'
a.,...
] [ l
bxy:
- 0,866
0,~00 ~
R,.
+ 3( - 0,866) + 2(0)]
0,500
[
0,~66
60 ~vw
-0.866
0,500 O
O
1
[6]
2
4
0,598
=
=
- l
4,964
2.0
[ 1,268]
6,196
4,0
As, ax,, y bxy, son iguales a (-0,598, 4,964, 2,0)7 y (1,268, 6,196, 4,0)7,
respectivamente, cuando se expresan en trminos del sistema de coordenadas de referencia.
O
Ejemplo: Si a_.y, = (4, 3, 2)7 y bx,: = (6, 2, 4)7 son las coordenadas con
respecto al sistema de coordenadas de referencia, determinar los puntos
correspondientes a.,. .,, b.,.,.. con respecto al sistema de coordenadas rotado
OUVW si ha sido girado 60 respecto del eje OZ.
u
o
X
1
1
1
1
lw
------y
<fJ
= -90
90
19
= [ - 0,866 0,500 0
2,0
0,500 0,866
h.,.... =
- o,~66 o,~oo
;J ln
3(0,866)
2(0)
3
2
= [- ~:~:]
[ 4(0,5)
[-!:I~Jj
(Pi~
Ry.~R,.eR.x,o
[ CO
-~</)
[ coca
SB
-S<j)CB
O SO][CO
I
O
so
C<j)
S S - C<j)S8Ca
CBCct
S<j)S8Ca
+ cos
-SO
CB
] [1
O
I
soca]
-~l
O
O
Six
cosos, +
- cesa
(2.2-14)
CdrCo: - S<j)S8Six
se
R,.R,,R,,
o
= [~ c
s
coco
CaSBC<j) + SaS<j)
SctSBC<j) - Ca.S<j)
-SO
-~[C O]
-SO
C8
s o
o
Ca
COSO
O:x.CB
CctSBS<j) -
SixCB
SixSOS<j)
SixC</)
cecs
] tq s:]
o
-S<j)
C<j)
,, (2.2-15)
CINEMATICA
DEL BRAZO
21
DEL ROBOT
3.
. _ Jr
resultante
Ejemplo:
\;
_ ;..
eje
OU.
SOLUCIN:
;"
R = R, .
'\
-[ c1
-S</>
[ csce
se
-S<f>C()
.,
13 R,..
o s1][ce -se
o se C()
o C</> O o
1
So S - C<f>SeC:1.
COC::r.
S</>SOC::r.
C</>Sx
= R, .
m~
Ru.
.,
R,.. Ru.
-s, ]-
Ccx
csses + s1c,]
-ces
C</>C':1. - S</>SOS::r.
O~~
22
--------------------~
/1
,,
" '---/
l. R,.a
2. R,.-P
3. R,.
4. R,.~
5. K... -a
1
1
1
1
1
1
11
1
1
<,
---..
/'
jI
-,
----------~--1
s -
1
1
1
1
1
1
1
"
,.
'
-2
-~-~-~-------4-
2.
W
respecto de r con ngulo </> y giramos el eje principal del sistema OXYZ para
volver el eje r otra vez a su posicin original. Con referencia a la figura 2.4, el
alineamiento del eje 02 con el eje r se puede hacer girando respecto del eje OX
con ngulo IX (el eje r est en el plano XZ), seguido por una rotacin de ngulo
- f3 respecto del eje O Y (el eje r ahora se alinea con el eje 02). Despus de la
rotacin del ngulo </> respecto del eje OZ o eje r, invertir las secuencias de
rotaciones anteriores con sus respectivos ngulos opuestos. La matriz de rotacin
resultante es
R,.<1>
Rx,
-: Ry.p
o] [ O
[i
C1X
SIX
+S
-S{J
CIX
[cpO o1 -:P] [~
S{J
C{J
e~.
C{J
sp]
o
C1X
S1X
-S</>
S</>
-~]
C</>
CIX
sen
IX
Jr;
r,
+ ,2
%
COS
IX
Jr;
r.
+ ,2
%
,2
;]
sen
f3 ==
cos
f3
Jr;
DEL ROBOT
23
(2.2-16)
'" =
J r;
+ r; + r; =
J3'
=-,
y
J3
l
J3
i3V</> + C</>
l
-S </>
/3V4> -
/) V</> +~:
S</>
i3V</> + C</>
v3
1/3
V</> -
/3
J3
'
R,.~
j3 S</>
/3
V</> +
J3l S</>
V</> +
3V</>-
1/
J31 S</>
fiS</>
i3V4> + C</>
~
2.2.4
24
Tabla 2.1
Secuencia
de
rotaciones
Sistema I ngulos
eulerianos
Sistema II ngulos
de Euler
C<J>Ci/ - S<f>COSI/
S<J>Ci/ + C<f>COSI/
iase
-~o]
[~! ~:"' ~]
ce
z. w
X.U
- C</>Si/1
- S<f>COCI/
1
- S<J>Si/
SOS!/
C<f>COCI/
S8Ci/
S<f>SO
-C<f>SO
ce
(2.2-17)
25
Una rotacin del ngulo <jJ respecto del eje OZ (R, . .,),
U na rotacin del ngulo 8 respecto del eje girado O V (R. 8).
Finalmente una rotacin de ngulo ,/ respecto del eje girado O W (R.... ,;),
sol [e"'
O
C<j)C8C,J - S<j)S,J
S<j)C8C,J + C<j)S,J
[
-S8C,J
SI/
o C8 O
- C<j)COS,J - S<j)C,J
- S<j)C8S,J + C<j)C,J
SOS!/! /
] =
-SI/
C,J
C<j)S8]
S<j)S8
C8
(2.2- 18)
V'
X, U
u
Figura 2.6 Sistema de ngulos eulerianos 11.
[e~s: ] [
c~ce
C<f>
-S~
O
O
1
O
CO
-58
:][~
C8
Ci/
o
Si/
s~si
S</>Ci/ cssoc +
C<f>S8Si/
S</>C8 S</>S8Si/ + C<f>Ci/ S<j)S8Ci/ - C<j)Si/
-58
C8Si/
C8Ci/
~~-]-
(2.2-19)
i:2.5
.: /
Giro
4>
>------Desviacin
y
Elevacin
CINEMATICA
DEL BRAZO
DEL ROBOT
27
t. Cada
coordenadas del sistema de referencia en trminos del sistema de coordenadas rotado OUVW?
Los nuevos vectores unitarios de los ejes de coordenadas se
hacen i,. = (1, O, 0)7, i. = (O, 1, 0)7 y k, = (O, O, 1)7 puesto que se
expresan
en trminos de ellos mismos. Los vectores unitarios originales son entonces
SOLUCIN:
sen exk,..
cos exkw
cos ex
~~
P:
w,
"'
CINEMATICA
DEL BRAZO
DEL ROBOT
29
escalado
x, ~
[I O
=
O cos
O sen
o
T
. , -
sen
o
o
[oos 6
o
C(
C(
-sen
cos
-sen (}
cos e
o
o
C(
C(
o
o
1
~]
~]
T,.
~[ c;s~ o
sen 4>
o
o
cos 4>
-sen iJ>
o
o
~]
(2.2-21)
30
Tt,ans
O O
dy
dz
o o o
(2.2-22)
P,
)'
P,
z
s
s
s
(2.2-25)
CINEMATICA
DEL BRAZO
DEL ROBOT
3)
Por tanto, el cuarto elemento diagonal en la matriz de transformacin homognea tiene el efecto de globalmente reducir las coordenadas si s > 1 y de alargar
tas coordenadas si O < s < 1.
En resumen, una matriz de transformacin homognea 4 x 4 transforma un
vector expresado en coordenadas homogneas con respecto al sistema de coordenadas OUVW en el sistema de coordenadas de referencia OXYZ. Esto es, con
w = 1,
Px)'Z
y
\,
2.2.7
(2.2-26a)
TPuvw
"-.. ;
j.;
["o o
Sx
ny
Sy
nz
Sz
"
P,J
ay
e,
az Pz
= [n
O
o o
~]
(2.2-26b)
En general, una matriz de transformacin homognea para un espacio tridimensional se puede representar como en la ecuacin (2.2-26b). Escojamos un punto p
tjo en el sistema de coordenadas OUVW y expresado en coordenadas homogneas como (O, O, O, 1 )7; esto es, Puvw es el origen del sistema de coordenadas
OVVW. Entonces la submatriz superior derecha 3 x I indica la posicin del
origen del sistema OUVW con respecto al sistema de coordenadas de referencia
OXYZ. Escojamos el punto p como (1, O, O, 1)7; esto es, Puvw = iu. Ms an,
suponemos que los orgenes de ambos sistemas de coordenadas coinciden en un
punto O. Este tiene el efecto de hacer los elementos en la submatriz superior
derecha 3 x I un vector nulo. Entonces la primera columna (o vector n) de la
matriz de transformacin homognea representa las coordenadas del eje OU de
OUVW con respecto al sistema de coordenadas OXYZ. Anlogamente, cogiendo
p como (O, l, O, 1)7 y (O, O, 1, 1)7, se puede identificar que la segunda columna (o
vector s) o la tercera columna (o vector a) de los elementos de la matriz homognea representan, respectivamente, los ejes O V y O W del sistema de coordenadas
OUVW con respecto al sistema de coordenadas de referencia. As, dado un
sistema de referencia OXYZ y una matriz de transformacin homognea T, los
vectores columnas de la submatriz rotacin representan los ejes principales del
sistema de coordenadas OUVW con respecto al sistema de coordenadas de
referencia, y se puede dibujar la orientacin de todos los ejes principales del
sistema de coordenadas OUVW con respecto al sistema de coordenadas de
referencia. El vector cuarta columna de la matriz de transformacin homognea
representa la posicin del origen del sistema de coordenadas OUVW con respecto al sistema de referencia. En otras palabras, una matriz de transformacin
homognea geomtricamente representa la localizacion de un sistema de coorde-
n"
T-1
S:,;
n1 n,
Sy
s,
(2.2-27)
[ a, ay a,
o o o
'--
2.2.8
DEL ROBOT
33
n
-f]
. SOLUCIN:
~
llxy=
o o
o
[~
o 1
o o
hxyz ~
~3(l)1)
1(05)
J [ 9J
3
~ = 2(1) 1(1~(-3) =
[~
o o
1 o
o 1
o o
-:
rn
)J
[' o o
O
= r; .r,. = ~ ~
~
[~
sen ex
cos
C(
b cos
::X
coso
::X
o ex
- sen
sen
cl
cos
[~
b sen ex
b sen ex
(X
f] ~
cos
(X
-sen
e<
c~s ]
1
-l~
-l~
o
cos
- sen a
cos (X
(X
sen a
cos (X
sen a
~]l~ ~]-
h
b sen
(X
(X
o 1
o o
,~so]
o
-sen
cos
o o
(X
Ejemplo: Encontrar una matriz de transformacin homognea T que represente una rotacin de ngulo '.l respecto del eje OX, seguida por una
traslacin de a unidades a lo largo del eje OX, seguida por una traslacin
de d unidades a lo largo del eje OZ, seguida por una rotacin de nguJo O
respecto del eje OZ.
SOLUCIN:
o
o 1
o
-sen (}
cos (}
sen &
-B
l' '
o
o
o
o
o o
1 o
o 1
o o
-l~
eos
sen (}
~H~
cos
cos
(X
(X
sen O
cos (}
] ll
o o
I
O
O
I
o 1
o o
cos
sen a
(X
-sen :x
cos (X
sen (X sen
-sen :x cos O
O
O
~]-
~]- ]
a cos
a sen (}
o
o
sen
(X
cos
(X
d
1
CINEMATICA
DEL BRAZO
DEL ROBOT
35
donde:
(2.2-28)
Tp
Un manipulador mecnico consiste en una secuencia de cuerpos rgidos, llamados elementos, conectados mediante articulaciones prismticas o de revolucin
(vase Fig. 2.8). Cada par articulacin-elemento constituye un grado de libertad.
Articulacin
2
Elemento 1 <,
Articulacl<in !
Elemento 4
Articulac in 6
+,
Figura 2.8
Planar
&7
Cilindrica
Prismtica
Eslerica
De tomiIIo
,.--...._,
'
37
Articulacin i + l
8,+ l
mento i.
En resumen, se asocian cuatro parmetros, a.; IX, d, y 91, con cada elemento de
un manipulador. Si se ha establecido un convenio de signo para cada uno de
estos parmetros, entonces constituyen un conjunto suficiente para determinar
completamente la configuracin cinemtica de cada elemento del brazo del robot.
Obsrvese que estos cuatro parmetros van apareados: los parmetros del elemento (a;, IX;) que determinan la estructura del elemento y los parmetros de la
articulacin (d1, 9) que determinan la posicin relativa de los elementos vecinos.
2.2.10
La representacin de Denavit-Hartenberg
CINEMATICA
DEL BRAZO
39
DEL ROBOT
Para una articulacin giratoria, d.; a; y IX; son los parmetros de articulacin
y permanecen constantes para un robot, mientras que 8 es la variable articulacin que cambia cuando el elemento i se mueve (o gira) con respecto al elemento
; - l. Para una articulacin prismtica, 8, a y IX son los parmetros de la
articulacin y permanecen constantes para un robot, mientras que d es la variable de la articulacin. Para el resto de este libro, la variable de la articulacin se
refiere a 8 (o d), esto es, la cantidad que vara, y los parmetros de articulacin se
refieren a los restantes tres valores geomtricos constantes (d;, a., et) para una
a~ticulacin giratoria o (8;, a.; IX;) para una articulacin prismtica.
Con las tres reglas bsicas anteriores para establecer un sistema de coordenadas ortonormal en cada elemento y la interpretacin geomtrica de los parmetros de la articulacin y del elemento, se presenta en el algoritmo 2.1 un proced-
...E
<.
~~
!'86~
ll3
r::
Y6 (s)
z6
(a)
x6 (n)
1
2
3
4
90
90
o
o
o
431.8 mm
149,09 mm
90
-90
90
-20,32 mm
-45 a 225
433,07 mm
-110 a 170
-100 a 100
-90
56,25 mm
-160 a +160
-225 a 45
-266 a 266
Y5
y6 (Barrido)
Centro de la mano OH
(Alcance) z6~
~z4
v
....._\
.J
(Elevacin) x6
~
Z3
84
d6~
Z5
X4
/,./
,-..:.---
,.,.,.
/---
---
x5
y4
------
los orgenes
coinciden
--
(d4 = ds = O)
x3~
x2
>-
d~~y2
82
/,,(,/~
d
e,
t=z
r \"'
f--=:
Yo--;-Y
=X
~o
a,
81 = -90
82 = -90
-90
-90
90
2
3
4
s
6
84 =
85 = O
116 = O
-90
90
a;
o
o
o
o
o
o
d
d,
d2
d
o
o
dt,
CINEMATICA
DEL BRAZO
DEL ROBOT
41
42
11
=[~
l
o
l
= [~
_
o l o
o o l
[cos
8;
sen O
o o
o
d ~]
[~
cos
(X
sen
tx1
~l
[cos 8;
-sen
sen 01
cos O
o
l
o
-sen
cos
~] =
tx
(X
8;1
- cos IX sen O;
cos (X tx
cos 8
sen
sen :X sen O;
-sen :X cos 8;
cos (X
a; cos
a sen 8
d,
(2.2-29)
1 A]-1
sen 8
cos
cos 8
sen
[
A-1
cos 8;
sen
IX
IX
sen 81
sen 8;
- sen
(X
IX
cos 81
et
COS IX
-d
-d
tx
COS IX
(2.2-30)
CINEMATICA
DEL BRAZO
DEL ROBOT
43
que
- i
[' '
cos
11
cos
- sen
cos
!l
s en
8O;
A -
sen
.o
cos
11
fl
il
(2.2-31)
y su inversa es
;A-1
p-1A;]-1
sen
cos O;
- cos
11
sen 8
cos
11
cos 8
sen 1
cos
il
sen 8
- sen
!l
cos 8;
cos a
-d; ~en
>;l
-d cos
11
(2.2-32)
la matriz i - 1 A se puede relacionar un punto p en reposo en el
elemento i y expresado en coordenadas homogneas con respecto al sistema
de coordenadas i en el sistema de coordenadas i - 1 establecido en el elemento
Utilizando
i - 1 por
(2.2-33)
0~1,
La matriz homognea T que especifica la localizacin del sistema de coordenadas i-simo con respecto al sistema de coordenadas de la base es el producto en
cadena de matrices de transformacin de coordenadas sucesivas ;- 1 A y se expresa como
T
A11A2
...
i-1A
j-1A.
J
para i
1, 2, ... .n
[~
o o
~] =
j=
[~
~]
(2.2-34)
cos O
i-lA;
sen O
cos
o
o
o
o -1o
o o
[ e,S1
A1 =
C3
iA3
S3
o
o
- cos
S3
CL
a3S3
o
o
sen (}
sen
cos O
-sen
11
a cos O
sen O
cos O a sen O
11
o
o
a,c,
-C3
o l
o o
11
[ e,
S4
3~=
o
o
a,c,
aiS2
d
i
o
o
-S4
C4
o
o
-1
,-
._,..,
, .:
-s6
[ e,
c6
o
o
Ss
s
/ 4A -
s-
[ eSs, o
;1
-Cs
o o
-s,
[ e, e,,
T1 s A1
1Ai
iA3
Fipn
-a2S2
Ci3
n
-
d2S,
+ di
a3Si3
C4 Cs C6 - S4S6
-C4CsS6 - S4C6
C4Ss
d6C4S5
S4
-S4CsS6
+ C4C6
S4Ss
d6S4S
s
d6Cs + d4
donde C; cos .; S;
o
o
-Sn
C,Si3
e,
S1 C23
o
T2 3~4Ass~
s~=
e, C6 + C4S6
-SsC6
SsS6
o
:11!
sen .; C;
o
,!!!!
cos(6;
6); S;
Cs
= sen(8;
8),
PUMA.
CINEMATICA
DEL BRAZO
DEL ROBOT
45
donde:
(X, y, zJ
= matriz
de orientacin del sistema de coordenadas i-simo establecido en el elemento i con respecto al sistema de coordenadas
de la base. Es la matriz particionada superior izquierda 3 x 3
de T.
p = vector de posicin que apunta desde el origen del sistema de
coordenadas de la base hasta el origen del sistema de coordenadas i-simo. Es la matriz particionada superior derecha 3 x 1
de T.
Especficamente, para i = 6, obtenemos la matriz T, T = 0T6, que especifica la
posicin y orientacin del punto final del manipulador con respecto al sistema de
coordenadas de la base. Esta matriz T se utiliza tan frecuentemente en la cinemtica del brazo del robot que se llama la matriz del brazo. Considere que la
matriz T sea de la forma:
T
= [~6 o
=
["
p.]
s,. a,
n, s, a, P,
n, s, a, Pz
o o o
[~
a
o
(2.2-35)
= vector
DEL ROBOT
47
[c,c,,
A3
T, =
S,C23
-S23
A11A/A3
-S,
C1
o
o
C1S23
S1S23
C23
(2.2-37)
(2.2-38)
donde Cu
cos (8 + 8) y Si
sen (8 + 8).
La matriz del brazo T para el robot PUMA que se muestra en la figura 2.11
se encuentra que es
(2.2-39)
donde:
nx
C1[CdC4C5C6
n,
=
=
=
=
=
S1[CdC4C5C6
11,
-S2JCC4C5C6
+ C4S6)
- S4S6) - S23S5C6] + C,(S4C5C6 + C4S6)
- S4S6)
- S4S6]
S23S5C6]
- S1(S4C5C6
(2.2-40)
C23S5C6
P,.
=
=
P,
S1[d6(C23C4S5
s"
s,
s,
ax ""
a, =
a,
P: =
C1[ -C23(C4C5S6
o o
o o
O,
-149,09]
O, 85
921,12
20,32
1
1]
OO 10
-1
-1
O 1200
-1
]
10
1
Cul es la posicin del centro del cubo con respecto al sistema de
coordenadas de la base?
b) Suponer que el cubo est dentro del alcance del brazo. Cul es la
matriz de orientacin [n, s, a] si necesita que la pinza (o los dedos) de
la mano se alineen con el eje y del objeto y al mismo tiempo coja el
objeto desde lo alto?
a)
DEL ROBOT
49
';r--x,
Y,/ !2
SOLUCIN:
CllmarTcubo
T,
-[~
o
o
o
o
o
-1
o
"'""T base -- T2 -[ ~O
o
Para encontrar
bmTcubo
-1
o
o
10
o
o
-1
-1100 1
2 0
1
baseT
= [~
-1
~ ~ 1 [~ ~
1]
10
matriz de transformacin resultante:
cubo
o
o
o
-1
O
10
1
O O
O O
Oo
o1
o o
1
1
Por tanto, el cubo est en la posicin ( 11, 10, 1 )7 del sistema de coordenadas de la base. Sus ejes x, y y z son paralelos a los ejes -y, x y z del sistema
de coordenadas de la base, respectivamente.
Para encontrar [n, s, a], hacemos uso de
T6
s
O O
a
O
PJ1
j
O
O
Hl
k
O
-1
o o
~J
-1J
CINEMATICA
DEL BRAZO
DEL ROBOT
51
Utilizando la matriz de rotacin con la representacin de los ngulos euleriano como en la ecuacin (2.2-17), la matriz del brazo 0T6 se puede expresar como:
C</JCt/1 S<jJCeStjJ - C</JSt/1 - S<jJCeCtjJ
p1
T = S<j}Ct/1 + C<jJCeStjJ -S</JSt/1 + C<jJCeCtjJ -C<j)Se
6
Pz
(2.2-44)
ses
[
sec
ssse
ce
Px]
I
csce
T _
S<j)Ce
6
-se
o
C<J>SeStjJ - S<j}Ct/J
S<jJSeStjJ + C<j}Ct/J
CeStjJ
O
C<jJSeCtjJ +
S<j}St/J S<jJSeCtjJ C</JSt/1 CeCtjJ
O
Px]
P,
Pz
1
(2._2-45)
.
.
1
o
donde R6
(RPY).
JJ
(2.2-46)
Coorcknadas clindricas para el subconjunto de posicionamiento. En una representacin de coordenadas cilndricas, la posicin del efector final se puede especificar
por las siguientes traslaciones/rotaciones (vase Fig. 2.15):
l.
2.
J.
La matriz de transformacin homognea que representa las operaciones anteriores se puede expresar como:
He.
I O O
O I d
O O I
O
O
-s o
O I O O = Srx
O O I O
O
O O O I
O
-Srx
Ca
o 1
o o
Ca
,cal
o
o
o 1
o o
rS
~l
(2.2-47)
Como solamente estamos interesados en los vectores de posicin (es decir, la cuarta
columna de Tcillndrical, la matriz de brazo 0T6 se puede obtener utilizando la
ecuacin (2.2-46):
T 6
,ex] [
o
= ['O oI o -s
o o o
y Px
= rCoc, py = rS, P:
(2.2-48)
R6
O
o ~]
d.
2.
3.
DEL ROBOT
53
lea
T,p1r = T,,aRy.{JT:,r =
O OH
-s
Sa
Ca.
O
O O
I O
O I
l~
o o
1o oo
o 1
O
r
_
-
rcp
-Sf3
O
=s
Sa.Cf3
-Sf3
ceO ol Sf3o
o Cf3
o o
,c,spl
CrxSfJ.
o
Ca. S':1.Sflo rSrxSf3
rCf3
o C{J
1
~1
(2.2-49)
T=
6
donde Px
rC'XS{3, py
l
=
o o
o o o
rS':XS{J, P:
rS'rXxSf{3Jl
rC{J
l
(2.2-50)
rC{J.
1/J) y (RPY). El
54
[
Tposicin
Px]
O
O P1
1 P:
~ ~
o o
Orientacin
Cartesiana [n, s, a]
Angulos de Euler (<P, 8, 1/1)
R,P, Y
rol
[[n, a]
R,p,8,.,
s,
~J
2.3
CINEMATICA
DEL BRAZO
DEL ROBOT
55
s, ax. Px.]
s,
a,
P, - o A i A
s,
O
P,
1
2A
JA 4A sA
4
(2.3-1)
6
Cc/>Ci/1 -
Scf>Ci/1
ssces
+ Ccf>C8Sif
S8Sif
- Ccf>Si/1 - ScpC8Cif
- Scf>Si/1 + CcpC8Cif
S8Cif
ScpS8]
- CcpS8
ce
(2.3-2)
Scf>Ci/1
S"
sy
=
=
CcpC8Sif
(2.3-3a)
(2.3-3h)
S8Sif
(2.3-3c)
- C <J>Si/1 - ScpC8Cif
(2.3-3d)
-Scf>Si/1
CcpC8Cif
(2.3-3e)
S8Cif
(2.3-3/)
"
ScpS8
(2.3-3g)
ay
-CcpS8
(2.3-3h)
a, =
ce
(2.3-3,)
sz
DEL ROBOT
57
Utilizando las ecuaciones (2.3-31), (2.3-3f) y (2.3-3h), una solucin a las nueve
ecuaciones anteriores es:
(} =
cos - 1 (a?)
VI =
cos -
1)
(2.3-4)
(2.3-5)
cos - l
( ;~)
(~;y)
(2.3-6)
are tg2 ( y, x)
O ~ 8 ~ 90
i 90 ~ (} ~ 180
-180 ~ (} ~ -90
i
-90 ~ 8 ~
oc
para +x y +
y
para -x y +y
para
-x
(2.3- 7)
y -y
para +x y -y
Utilizando la funcin arco tangente (are tg2) con dos argumentos, tendremos una
visin de la solucin general propuesta por Paul y col. [1981].
En la ecuacin matricial [ecuacin (2.3-2)] se dan los elementos de la matriz en
el lado izquierdo, mientras que los elementos de las tres matrices del lado derecho
son incgnitas y son dependientes de q>, 8, VI Paul y col. [1981] sugieren
premultiplicar la ecuacin matricial anterior por su transformada inversa desconocida
sucesivamente y de los elementos de la ecuacin matricial resultante determinan el
ngulo incgnita. Esto es, movemos una incgnita (por su transformada inversa)
del lado derecho de la ecuacin matricial al lado izquierdo y resolvemos para la
incgnita, a continuacin movemos la siguiente incgnita al lado izquierdo y
repetimos el proceso hasta que se resuelven todas las incgnitas.
Prernultiplicando la ecuz cin matricial anterior por R;: i tenemos una incgnila (q,) en el lado izquierdo y dos incgnitas (8, VI) en el lado derecho de la ecuacin
matricial, as tenemos
Sq> Cq> o
~]
=
SVI
nr
s,
a,
O C(} -S{}
O]} [nxnz
x] [l S(J
O C(}O ] [CVI
s,
s, a,
-SVI
CVI
58
C<J)sx + S<J)sy
+ Sn ;
-S<J)sx + Ces;
- Sen, + C<J)n y
-St/1
C<J)nx
cec
s
sec
n~
-~e]ce
(2.3-8)
tg-1
~:J
(2.3-1 O)
= Con, +
= -C</>sx -
(2.3-1 la)
S</>n,
(2.3-1 lh)
S<J>s,
(St/1)
CiJ
tg -i
= tg
(-C</>sx - S</>s,)
Cdm; + S</>n1
+ S<J>n)
(2.3-12),
se = Soa;
ce = a,
que nos da la solucin para
_1
(se )
-
ce
(2.3-13)
e,
= tg 1(S<l>ax-C<l>a>)
ull =tg
- C<J>a,
] [ o
cOe -seO
I
CO
]
S8
CINEMATICA
DEL BRAZO
DEL ROBOT
59
nxSi/l+sxCi/1
nySi/1 + syCi/1
n,Si/1 + s,Ci/1
"']
[C<j)
ay
S<j)
=
a,
-S<j)C()
C<j)C()
se
S<J>SB]
- C<j)SB (2.3-15)
ce
Otra vez, igualando los elementos (3, 1) de ambas matrices en dicha ecuacin
matricial, tenemos:
(2.3-16)
que da
i/1
tg "
(n)
s.
(2.3-17)
se =
C()
n,Si/1
(2.3-18a)
+ s,Cij
(2.3-18b)
a,
()
tg "
(n,Sij
+ s,Cij, a,)
(2.3-19)
s,Ci/1)
nyCi/1 - sySi/1
(2.3-20a)
(2.3-20b)
que da
<P = tg_1 (nYCiJ - sySi/1)
nxCi/1 - sxSi/1
(2.3-21)
60
PUMA utilizan los smbolos O, A, T para indicar los ngulos de Euler y sus
definiciones se dan a continuacin (con referencia a la figura 2.17):
O (orientacin) es el ngulo formado desde el eje Yo hasta la proyeccin del
eje a de la herramienta sobre el plano XY respecto del eje z0.
A (altitud) es el ngulo formado desd_e.-etplanoXY hasta el eje a de la
herramienta respecto del eje s de 1.',.,herramienta.
T (herramienta) (tool) es el ngulo formado desde el plano XY hasta el ejes
de la herramienta respecto del eje a de la herramienta.
.. '!,.
~{j).....
. ...
,
,
<::: ; :; ,;:~\)
...
..
'D0 ::.
r
,
A una medida rdel ngulo formado
ent e la HERRAMIENTA Z y
..
.....
..
' 'T v,'lJo. ,>..-, . ..,.,0....-~o
;-
HERRAMIENTA
'{
;-.
..
..
....
-,
'!,.
.. ~~{j)
Figura 2.17 Definicin de los ngulos de Euler O, A y T. (Tomado del manual del robot
PUMA 398H.)
CINEMATICA
DEL BRAZO
DEL ROBOT
61
Zo
~
-1
~
O
~]
(2.3-22)
De la definicin de los ngulos OAT y la matriz de alineacin inicial [ecuacin (2.3-22)], la relacin entre la transformacin de la mano y los ngulos OA T
est dada por
l
r
~] u
a,
s,.
-so
so co
= ca
o
["'
>'
=
'1> Sy
n: s,
= R=.O
u]
-1
O1
o
o
[
[CST
-!]
.:
C
AI
R,
,R.'
O SA]
O
CA
-ST
CT
nY sY ay
n, s. =
l[
-ST
CT
] ca
o = so
I
co
O
~] u
]o
-1
CA
-SA
o
1
s;]
CA
62
(2.3-23)
s,CT = O
(2.3-24)
(2.3-25)
lo que la da solucin de T,
T
tg-1
( ~~.)
= - n,CT +
(2.3-26a)
(2.3-266)
s,ST
tg-1
-a,
)
- n,CT + s:ST , '
i=
-n,CT
s,ST)
(2.3-27)
n.S]'
nyST
(2.3-28a)
(2.3-28b)
+ sxCT
+ s)'CT
tg-
are tg2
+ syCT)
+ s.ct
(n,ST + syCT, nxST +
(nyST
nxST
sxCT)
(2.3-29)
La premultiplicacin o la postmultiplicacin anterior de las transformadas inversas desconocidas se pueden tambin aplicar para encontrar la solucin de las
articulaciones de un robot PUMA. Los detalles respecto a la solucin del robot
PUMA se pueden encontrar en Paul y col. [1981].
Aunque la tcnica de la transformada inversa proporciona un mtodo general
para determinar la solucin de las articulaciones de un manipulador, no da una
indicacin clara sobre cmo seleccionar una solucin apropiada de las diversas
soluciones posibles para una configuracin de brazo particular. Esto tiene que
descansar sobre la intuicin geomtrica del usuario. As, un mtodo geomtrico
CINEMATICA
DEL BRAZO
DEL ROBOT
63
T6
=-1~B
-1
ThcrramicnlaH
r
~
s,. a,.
s, a,
o o o
P,
n, s.
l
P,l
(2.3-30)
Pz
Definicin de las diversas configuraciones de brazo. Para un robot PUMA mostrado en la figura 2.11 (y otros robots giratorios), se definen diversas configuraciones
de brazo de acuerdo con la geometra del brazo humano y el sistema de coordenadas de elementos que se estableci utilizando el algoritmo 2.1 como (Fig. 2.19):
BRAZO DERECHO (hombro): 02 positivo mueve la mueca en la direccin z0
positivo mientras la articulacin tres no se activa.
BRAZO IZQUIERDO (hombro): 02 positivo mueve la mueca en la direccin
z0 negativa mientras la articulacin tres no se activa.
BRAZO ARRIBA (codo por encima de la mueca): Posicin de la mueca del
brazo { DERECHO} con respecto al si. stema de coordenadas del hombro
IZQUIERDO
. y2.
negativo} a lo largo del eje
tiene valor de coordenada {
.
.
posiuvo
BRAZO ABAJO (codo por debajo de la mueca): Posicin de la mueca del
.
de coordenadas del hombro
DERECHO} con respecto al sistema
brazo
{ IZQUIERDO
tiene valor de coordenada
positivo}
.
{ negativo
brazo DERECHO
brazo IZQUIERDO
(2.3-31)
CODO
{+1
-1
brazo ARRIBA
brazo ABAJO
(2.3-32)
MUECA
{+1
-1
mueca ABAJO
mueca ARRIBA
(2.3-33)
Brazo izquierdo
65
y arriba
{+
1
_
(2.3-34)
66
v1s10:-.:
E 1:-;TELIGE:-.:CIA
Solucin del brazo para las tres primeras articulaciones. Del diagrama cinemtico del robot PUMA en la figura 2.11, definimos un vector de posicin p que
apunta desde el origen del sistema de coordenadas del hombro (x0, Yo, z0) hasta el
punto donde intersecciona el ltimo de los tres ejes de la articulacin como (vase
figura 2.14):
(2.3-35)
et=<P-l.
r = Jp;
p; - dzz
sen <P = P,
R
dz
sen l.
R
R = Jp;
(2.3-37)
+ p;
cos <P = P- x
R
COS IX
(2.3-38)
(2.3-39)
(2.3-40)
donde los superndices L y R sobre los ngulos de las articulaciones indican las
configuraciones de brazo IZQUIERDO/DERECHO. De las ecuaciones (2.3-37) a
(2.3-40) obtenemos las funciones seno y coseno de 01 para las configuraciones de
brazo IZQUIERDO/DERECHO:
sen
et
cos
et
cos (<P -
ix)
ix
(2.3-41)
(2.3-42)
sen
Of
cos
ef
= sen (n + <P + a)
cos (n
+ </) + a)
(2.3-43)
(2.3-44)
66
v1s10:-;
E ISTELIGENCIA
Solucin del brazo para las tres primeras articulaciones. Del diagrama cinemtico del robot PUMA en la figura 2.11, definimos un vector de posicin p que
apunta desde el origen del sistema de coordenadas del hombro (x0, Yo, z0) hasta el
punto donde intersecciona el ltimo de los tres ejes de la articulacin como (vase
figura 2.14):
(2.3-35)
8=</)-'Y.
r
Jp;
p; -
sen
<P
sen
IX
Py
R
d2
R
(2.3-37)
Jp; + p;
(2.3-38)
<P
Px
R
(2.3-39)
COS IX
r
R
(2.3-40)
d2
2
cos
donde los superndices L y R sobre los ngulos de las articulaciones indican las
configuraciones de brazo IZQUIERDO/DERECHO. De las ecuaciones (2.3-37) a
(2.3-40) obtenemos las funciones seno y coseno de 81 para las configuraciones de
brazo IZQUIERDO/DERECHO:
sen
ef =
cos
ra .
sen
= cos
(<P -
7.)
sen
<P
cos :x - cos
<P
sen 'Y.
(A.
IX)
= cos
,1,
'I'
cos a
+ sen
<P
sen a -
'I'
(2.3-41)
_
Px'
+R2
(2.3-42)
Pyd2
sen
ef
COS
e R
sen (1t
<P
IX)
(2.3-43)
COS
(1t
+ </)
IX)
- Px' + Pyd2
-----=~ -
R 2
(2.3-44)
67
Plano x0y0
OA = d2
A8
1
1
\
08 =
Brazo izquierdo
' ,....._
J p2,
J p2,
+ p2y -d22
+ p2y = R
./
//
....
Yo
OA = di
(p,. p_,.)
,....-
A8 = r =
I
~=r+,t,+a
J P2
X
+ P~-d2
.'f
2
\
\
08 =
lfc"~~~.,...+-'-~~~~-xo
Jp; + p;
-,
Brazo derecho
- BRAZO pyj p;
p;
+ p; + p;
d -
Px'z
(2.3-45)
(2.3-46)
ir..
p ; p,)
O AA = d1 EF = P.,
EG =
BBC = a1
P,. C = a.1
DE =
P,
D= d4
AD = R = J P: + P; + P; - di
E = r =
dl
pl
+ p2 -
,\'
tg
_1
(-BRAZOpyjp; + p; - dJ - Pxdi)
- BRAZO PxJ p; + p; - d'f + pyd2
(2.3-47~
Configuraciones
Brazo
Brazo
Brazo
Brazo
de brazo
IZQUIERDO
IZQUIERDO
DERECHO y
DERECHO y
y ARRIBA
y ABAJO
ARRIBA
ABAJO
82
f3
+ f3
+ f3
- f3
IX :X
IX
IX
BRAZO
CODO
BRAZO CODO
-1
-1
+I
-1
+l
+t
+I
-1
-1
+!
+t
-1
82
+ (BRAZO COD0)/3 = a + K /3
(2.3-48)
r =
R = J p; + p; + p;
sen
P,
!X
Jp; + p; -
df
(2.3-49)
P,
J p; + p; + p; - df
= -- R =
(2.3-50)
2
BRAZO r
cosa
cos
BRAZO Jp;
2 + R2
2
P1
d2
(2.3-51)
/3 =
(2.3-52)
(df + an
2a2R
-
p; + p; + p; + a~ - df - (d + a~)
2a2Jp; + p; + p; - df
sen
/3 = Ji
- cos2
(2.3-53)
sen a cos
=
=
cos
/3) =
sen
!X
cos (K
/3) +
cos
!X
sen (K
sen
/3
/3
(!X
+ K /3)
!X
/3)
(2.3-54)
(2.3-55)
8 2 -_ tg _
(s e n 82)
- cos 82
(2.3-56)
70 ROBOTICA:
Plano
co;-;TROL.
UETEC'C'ION,
YISION
E INTELIGENCIA
,,y,
D
BD=~
AD=
o.,=
ef -
R =.l
/3
~D
0.1
X_
ef, -
/3
90
de (x2, y 2, z2) hasta el punto donde se intersectan los ltimos tres ejes de
articu- lacin.
De la geometra del brazo en la figura 2.22 obtenemos las ecuaciones siguientes para encontrar la solucin a ('3:
(2.3-57)
+ p; - d
R=
cos
<P
Jp;
p;
ai + (d +
2a Jd
+
sen
<P
sen
/3
ai) - R2
+ ai
BRAZO CODO
Jd
d4
+
a~
(2.3-58)
jt
cos2 <j)
/3
laJI
+
cos
jd
a~
(2.3-59)
CINEMATICA
DEL BRAZO
DEL ROBOT
71
(2p4)y
03
) o
tj)-{3
-1
+l
-1
~o
tj)-{3
-l
-1
+l
~o
tj)-{3
+I
+l
+l
)O
tj)-{3
+l
-1
-1
sen (</> -
/3) =
/3) =
sen
cos
</>
</>
cos
cos
/3 /3 +
cos
sen
</>
</>
/3
(2.3-61)
sen /3
(2.3-62)
sen
i
3.
72
a
s
=
=
dado a
(2.3-64)
Z5
dado a
(2.3-65)
Y6
dado s
(2.3-66)
n Ys
si s y5
o
n=
(Z3 X
~Z3 X
(Z3 X
n
llz3 X
I
a)
aI
si
(z , X a)
si S (Z3 X a)
,f.
(2.3-68)
CINEMATICA
1abla 2.5
DEL BRAZO
DEL ROBOT
73
Or entacin de mueca
l=s y5 o n y5
MUECA
ABAJO
ABAJO
ARRIBA
ARRIBA
~o
<o
~o
< o
+I
+1
-1
-1
+I
-1
-1
+1
SI X
1------
z,
1
1
sen 84 = - z,
(
1
~J)
Figura 2.23
Solucin para la articulacin 4.
<
(2.3-71)
Si ocurre el caso degenerado, se puede escoger para 84 cualquier valor conveniente mientras se satisfaga la orientacin de la mueca (ARRIBA/ABAJO). Esto se
puede asegurar siempre fijando 84 igual al valor actual de 84. Adems de esto, el
usuario puede activar el conmutador FLIP para obtener la otra solucin de 84,
esto es, 84 = 84 + 180.
Solucin de la articulacin 5. Para encontrar 85, utilizamos el criterio que alinea
el eje de rotacin de la articulacin 6 con el vector de aproximacin (o a = z5).
Mirando la proyeccin del sistema de coordenadas (x5, y5, z5) sobre el plano
x4y4, se puede demostrar que es cierto lo siguiente (vase Fig. 2.24):
(2.3-72)
donde x4 e y4 son, respectivamente, los vectores columnas x e y de T4 y a es el
vector de aproximatn. As, la solucin a 85 es:
- 85
5)
tg " i (sen 8
cos 85
Si 85
7t ~
5~
7t
(2.3- 73)
Solucin de la articulacin 6. Hasta ahora hemos alineado el eje de la articulacin 6 con el vector de aproximacin. A continuacin necesitamos alinear la
orientacin de la pinza para facilitar la recogida del objeto. El criterio para hacer
X.
CINEMATICA
DEL BRAZO
DEL ROBOT
75
n y5
cos
()6
S .
Y5
(2.3-74)
C1C23S4)nx
C1C23S4}sx
+
+
(C1C4
(C1C4
- S1C23S4)n1
- S1C23S4)s1
+
+
(S4S23)n=J
(S4S23)s,
(2.3- 75)
La derivacin anterior de la solucin cinemtica inversa de un robot PUMA se
basa en la interpretacin geomtrica de la posicin del punto final del elemento 3
y del requisito de orientacin de la mano (o herramienta). Existe un inconveniente en la derivacin anterior para 04, 05 y 06. El criterio para fijar el eje de
movimiento de la articulacin 5 que es igual al producto vectorial de z3 y a puede
no ser vlido cuando sen Os : : : : O, lo que significa que Os : : : : O. En este caso, el
manipulador se hace degenerado con los ejes de movimiento de las articulaciones
4 y 6 alineadas. En este estado, solamente es significativo la suma de 04 y 06. Si
ocurre el caso degenerado, entonces somos libres de escoger cualquier valor para
64, y normalmente se utiliza su valor actual y a continuacin nos gustara tener
64 + 06 igual al ngulo total necesitado para alinear el vector de deslizamiento s
y el vector normal n. Si el conmutador FLIP est activado (es decir, FLIP = 1),
entonces ()4 = ()4 + n, Os = - Os Y 06 = 06 + n.
En resumen, existen ocho soluciones al problema cinemtico inverso de un
robot tipo PUMA de seis articulaciones. La solucin de las tres primeras articu-
X1
laciones (81, 82, 83) posicionan el brazo mientras que las soluciones de las ltimas
tres articulaciones (84, 85, 86) proporcionan la orientacin apropiada de la mano.
Hay cuatro soluciones para las tres primeras articulaciones. Dos para la configuracin del brazo del hombro derecho y dos para la configuracin del brazo
del hombro izquierdo. Para cada configuracin de brazo, las ecuaciones (2.3-47),
(2.3-56),(2.3-63),(2.3-71 ), (2.3-73) y (2.3- 75) dan un conjunto de soluciones (81, 82, 83,
84, 85, 86) y (81, 82, 83, 84 + n, -85, 86 + rr) (con el conmutador FLIP activado)
da otro conjunto de soluciones.
Ecuaciones de decisin para los indicadores de configuracin de brazo. La solucin derivada para un brazo tipo PUMA en la seccin anterior no es nica y
depende de los indicadores de configuracin de brazo especificados por el usuario. Estos indicadores (BRAZO, CODO y MUECA) se pueden determinar
tambin a partir de los ngulos de las articulaciones. En este apartado deducimos
las ecuaciones de decisin respectivas para cada indicador de configuracin de
brazo. El signo de la ecuacin de decisin (positivo. cero o negativo) proporciona
una indicacin de la configuracin de brazo tal como se definieron en las ecuaciones (2.3-31) a (2.3-33).
Para el indicador BRAZO, siguiendo la definicin del brazo DERECHO/IZQUIERDO, se puede encontrar una ecuacin de decisin como:
Z X p'
k
+sen el
Zo. -ijz-,-x-p- Zo.
o [z x p'I
g(8, p)
='j
Px
O
-p>'
sen 81
p., cos 81
(2.3-76)
llz1 x p'I
donde p' = (p.,, py, 0)7 es la proyeccin del vector de posicin p [Ec. (2.3-36)]
sobre el plano x0y0, z1
(sen 1 cos 81, 0)7 del vector tercera columna de 0T1
, y z0 = (O, O, 1)7. Tenemos las posibilidades siguientes:
l.
2.
3.
A sen 81)
(2.3-77)
CINEMATICA
DEL
BRAZO
DEL ROBOT
77
donde la funcin signo se defini en la ecuacin (2.3-70). Sustituyendo las componentes x e y de p de la ecuacin (2.3-36), la ecuacin (2.3-77) se hace:
BRAZO = sign [g(O, p)] = sign [g(O)] = sign ( -d4S23
a3C23
a2C2)
(2.3-78)
+1
{ -1
=>
=>
brazo DERECHO
brazo IZQUIERDO
(2.3-79)
+ 1 si S z4 > 0
{ -1 Si S Z4 <
MUECA
Anulos de articulacin
~------~
Cinemtica directa
Posicin y orientacin
del efector final
+
Ecuaciones de decisin
BRAZO, CODO, MUiil'ECA
Cinemtica inversa
(2.3-81)
sign (n z 4)
(2.3-82)
sin z4 < O
DEL ROBOT
79
CA). Hay ocho soluciones para un robot tipo PUMA -cuatro soluciones para
Jas tres primeras articulaciones y para la configuracin del brazo, dos soluciones
ms para las ltimas tres articulaciones-. Se puede verificar mediante simulacin en computadora la validez de las soluciones cinemticas directa e inversa.
Con modificaciones y ajustes apropiados se puede generalizar el mtodo geomtrico a otros robots industriales simples con articulaciones de tipo giratorio. Los
conceptos cinemticos que se han tratado en este capitulo se utilizarn de forma
amplia en el captulo 3 para deducir las ecuaciones de movimiento que describen
Ja conducta dinmica de un brazo de robot.
REFERENCIAS
Ms informacin sobre matrices se puede encontrar en Bellman [ 1970], Frazer y
colaboradores. (1960] y Gantmacher (1959]. La utilizacin de matrices para
describir la posicin de un elemento mecnico rgido se puede ver en el trabajo de
Denavit y Hartenberg (1955] y en su libro (Hartenberg y Denavit (1964]). Ms
informacin sobre coordenadas homogneas se puede ver en Duda y Hart [ 1973]
y Newman y Sproull (1979]. La discusin sobre cinemtica es una extensin de
un trabajo de Lee [ 1982]. Una ampliacin sobre cinemtica se puede encontrar
en Hartenberg y Denavit (1964] y Suh y RadclifTe(1978]. Aunque la representacin matricial de los elementos presenta un enfoque sistemtico para resolver el
problema cinemtico directo, el mtodo vectorial para el problema cinemtico
presenta una representacin ms concisa de los mismos. Esto se estudia en un
trabajo por Chase (1;963]. Otros libros de robtica que analizan el problema
cinemtico son Paul (1981], Lee, Gonzlez y Fu (1986] y Snyder (1985].
Pieper [ 1968], en su tesis doctoral, utiliz un mtodo geomtrico para resolver el problema cinemtico inverso. El estudio de la tcnica transformada inversa
para encontrar la solucin del brazo se bas en el trabajo de Paul y col. (1981].
El mtodo geomtrico para resolver la cinemtica inversa de un manipulador de
seis elementos con articulaciones giratorias se bas en el trabajo de Lee y Ziegler
(1984]. La solucin del brazo de un robot tipo Stanford se puede ver en un
informe de Lewis [ 1974]. Otras tcnicas para resolver la cinemtica inversa se
pueden ver en los artculos de Denavit (1956], Kohli y Soni (1975], Yang
Y Freudenstein ( 1964], Yang [ 1969], Yuan y Freudenstein [ 1971 ], DufTy y
R.ooney (1975], Uicker y col. (1964]. Finalmente, el libro tutorial editado
por Lee, Gonzlez y Fu (1986] contiene numerosos trabajos recientes sobre
robtica.
PROBLEMAS
('i)
Cul es la matriz de rotacin para una rotacin de 30'' respecto del eje OZ, seguida
una rotacin de 60" respecto del eje OX, seguida por una rotacin de 90 respecto del
CJC OY?
\~
),Cul es la matriz de rotacin para una rotacin de ngulo rjJ respecto del eje OX,
(~i2
-segtida por una rotacin de ngulo i/J respecto del eje O W, seguida por una rotacin de
ngulo O respecto del eje O Y?
..13 Encontrar otra secuencia de rotaciones que sea diferente del problema 2.2, pero que
.,~~lle
2.4 Deducir la frmula para sen (rp + 8) y cos (rp + O) desarrollando simblicamente dos
rotaciones de rjJ y (} utilizando los conceptos de matriz de rotacin estudiados en este
captulo.
(:5) Determinar una matriz T que representa una rotacin de ngulo a. respecto del eje
seguida por una traslacin de b unidades de distancia a lo largo del eje OZ, seguida
por una rotacin de ngulo rjJ respecto del eje O V.
"oz,
'2~Para
Y1
Yo
( ~
~n
Para la figura que _se muestra en la pgina 81, encontrar las matrices de transforma
homognea 4 x 4 - 1 A; y O A; para i = 1, 2, 3, 4.
/~
Se ha preparado una estacin de trabajo de robot con una cmara de TV, tal como
: s~uestra en la seccin 2.2.11. La cmara puede ver el origen del sistema de coordenadas
ele la base donde se fija un robot de seis elementos, y tambin el centro de un cubo que
tiene que ser manipulado por el robot. Si se ha establecido un sistema de coordenada local
en el centro del cubo, entonces este objeto, tal como lo ve la cmara, se puede representar
por una matriz de transformacin homognea T1. Tambin el origen del sistema de
81
.,.,
?\
5 in
1
1
1
1
1
1
1
4in
.j---- --------
,, 1},\
/
~l
'
"/
, ,.,.,.
...
Y1
~ll
1---....,__
---(
3 in
coordenadas de la base tal como lo ve la cmara se puede expresar por una matriz de
transformacin homognea T2, donde
-[
'~]
T,
y
-ri
o
o
-1
o
o
-1
-IO]
20
10
2;0k:stablecer el sistema de coordenadas del elemento ortonormal (x., Y;, z.) para i = 1,
2-,- ... , 6 para el robot PUMA 260 que se muestra en la figura siguiente y completar la
tabla.
Rotacin de cintura 330
ffi
Rotacin de hombro 3 IO
8,0\ ~
Rotacin de codo
/~ci~
:.;:"' ~'(
~:
_>;;;!:'~ndc
13 in
Parmetrosde coordenadasde
los elementos del robot PUMA
Articulacin
11,
~,
d,
a,
2
3
4
5
6
/~
'~i--4= 1, 2, ... , 5 para el robot MINI MOVER que se muestra en la figura siguiente y
completar la tabla.
'
CINEMATICA
DEL BRAZO
DEL ROBOT
83
O,
~,
a,
d,
3
4
5
'i:ii).
~n robot tipo Stanf?rd se ha movid.~ a las posiciones que se_.muestran ,en I~ figu,ra.
= (90, - 120, 22 cm, O", 70, 90 )7.
Establecer el sistema de coordenadas del elemento ortogonal (X, y, z) para i = l, 2, ... , 6
para este brazo y completar la tabla.
3
4
5
6
2.Jl Utilizando las seis matrices i- 1 A (i = l, 2, .... 6) del robot PYMA de la figura 2.13,
encontrar su error de posicin al final del elemento 3 debido al error de medida de los tres
P11meros ngulos de articulacin (601, 602, 603). Es suficiente una solucin aproximada
de primer orden.
2.14 Repetir el problema 2.13 para el robot Stanford que se muestra en la figura 2.12.
~n
la figura se muestra un manipulador con dos grados de libertad. Si la longitud
de cada elemento es de un metro, establecer sus sistemas de coordenadas de elementos y
encontrar O A I y I A2 Encontrar la solucin cinemtica inversa para este manipulador.
/~
-\
'----lW
Para el robot PUMA que se muestra en la figura 2.J 1, suponer que hemos
encontra- do correctamente la solucin de las tres primeras articulaciones (81; 82, 83) y
que se dan
1
; A;, i = 1, 2, ... , 6 y T6. Utilizar la tcnica de transformacin inversa para encontrar la
solucin de los tres ltimos ngulos de las articulaciones (04, 05, 06). Compare
sus soluciones con las ecuaciones dadas (2.3-71 ). (2.3-73) y (2.3-75).
2.17 Para el robot Stanford de la figura 2.12, obtener la solucin de los tres primeros
ngulos de articulacin. Puede utilizar cualquier mtodo con el que se sienta cmodo.
2.18
CAPITULO
TRES
DINAMICA DEL BRAZO DEL ROBOT
Lo inevitable se supera mediante el esfuerzo
Olicer Wendell Holmes
3.l
INTRODUCCION
~}1
La dinmica del robot trata con las formulaciones matemticas de las ecuaciones
de movimiento del brazo. Las ecuaciones de movimiento de un manipulador son
un conjunto de ecuaciones matemticas que describen su conducta dinmica.
Tales ecuaciones son tiles para la simulacin en computadora del movimiento del
robot, el diseo de ecuaciones de control apropiadas para el robot y la evaluacin del diseo y estructura del brazo. En este captulo nos centraremos en la
formulacin, caractersticas y propiedades de las ecuaciones dinmicas de movimiento que son adecuadas con fines de control. El objetivo del control de un
manipulador basado en computadora es mantener la respuesta dinmica del
mismo de acuerdo con algn rendimiento del sistema preespecificado y objetivos
deseados. En general, el rendimiento dinmico de un manipulador depende directamente de la: eficacia de los algoritmos de control y de su modelo dinmico. El
problema de control consiste en obtener modelos dinmicos del brazo del robot
Iisico y a continuacin especificar leyes o estrategias de control correspondientes
para conseguir la respuesta y rendimiento del sistema deseado. Este captulo trata
principalmente con la primera parte del problema de control del manipulador;
esto es, la modelizacin y evaluacin de las propiedades y conducta dinmica de
robots controlados por computadora.
El modelo dinmico de un robot se puede obtener a partir de leyes fisicas
conocidas tales como las leyes de la mecnica newtoniana y lagrangiana, Esto
conduce al desarrollo de las ecuaciones de movimiento dinmico para las diversas articulaciones del manipulador en trminos de los parmetros geomtricos e
inerciales de los elementos. Mtodos convencionales como las formulaciones de
Lagrange-Euler (L-E) y Newton-Euler (N-E) se pueden aplicar entonces sistemticamente para desarrollar las ecuaciones de movimiento del robot. De estas dos
formulaciones se obtienen diferentes formas de describir la dinmica del brazo del
robot, tales como las ecuaciones de Lagrange-Euler de Uicker (Uicker [1965],
Bcjczy [1974]), las ecuaciones recursivas de Lagrange de Hollerbach (Hollerbach ;
[1980]), las ecuaciones de Newton-Euler de Luh (Luh y col. [1980a]) y las
ecuaciones generalizadas d' Alembert y Lee (Lee y col. [ 1983]). Estas ecuaciones
de movimiento son equivalentes unas a otras en el sentido de que describen la
Conducta dinmica del mismo robot fisico. Sin embargo, sus estructuras pueden
DINAMICA
DEL BRAZO
DEL ROBOT
87
3.2 FORMULACION
DE LAGRANGE-EULER
{J~
88
d (L)
dt
o~
oL
l, 2, ... , n
(3.21)
oq
donde
=
=
DINAMICA
DEL BRAZO
DEL ROBOT
89
con las variables articulaciones definidas en cada una de las matrices de transformacin de coordenadas de elementos 4 x 4. As, en el caso de una articulacin
giratoria, q = O;, que es el margen del ngulo de la articulacin; mientras que
para una articulacin prismtica, q = d que es la distancia recorrida por la
articulacin.
La siguiente deduccin de las ecuaciones de movimiento de un manipulador
con n grados de libertad se basa en las matrices de transformacin de coordenadas homogneas desarrolladas en el captulo 2.
3.2.1
~t
, -[n-
(x,, y,,
z,,
!)'
(3.2-
2)
r-
1 A;
_ 1
AI
f
=
e,
cos
sen O;
o
o
a,
cosO,l
a; sen O;
d
1
(3.2-5)
Figura 3.1
[cm 8,
i-1A;
sen u;
-cosa. sen fJ
cos :X cos () i
sen a.;
sen a. sen fJ
- sen a, cos fJ;
cos :X
(3.2-6)
91
i-lA
1A
1.
A;'r
-A
OA
1
1A
+ ...
L .
u:,OA
=i
A;'r
1=1
(3.2- 7)
'r,
-:,-rt
uqj
Q, ~
[f
o
o
o
o
o
o
o
~]
(3.2-8a)
[-~ne
cos 8.'
ai-1A
06
. [I
a,
-1
o
o
[eme,
o
o
- cos 'X sen 8
-a
~ne ]
cos 8;
o
o
sen 'X sen 8;
cose,]
8;, y
o
o
o
== Q/-IA
o
o
o
sen 8;
o
o
- sen a; cos 8
cos '.1.
a;
sen 8;
d
l
92
De aqu, para i
;iOA
_u_
oq1
=
=
1, 2, ... , n,
{A 1 1A
J-2A
J-1
QJ i : lAJ
lA
i-
para}~
para}>
(3.2-10)
A j-1
Q j J-1A
i ~ k ~ j
1-1
i~j~k
(3.2-13)
i<Joi<k
=j =k =
DINAMICA
DEL BRAZO
DEL ROBOT
93
(P~
=
=
(:x + y + i) dm
1/2
traza (vvT) dm = 1/2
1/2
Tr (vvT) dm
(3.2-14)
(3.2-15)
XZdm
f
yf
dm
..------J
== f r. r: dm ==
. dm
/
f y:
{3.2-17)
dm
-r
f -~
'
Tr A ~
L01,
[ f X:
dm fxdm
f y:dm
f : dm
fy,dm
f:dm
f-dm
.
f dm
94
donde \ = (X, y, z;, 1 )7 se define como antes. Si utilizamos el tensor de inercia /11
que se define como
donde los ndices i, j, k indican los ejes principales del sistema de coordenadas
i-simo y i es la delta de Kronecker, entonces J; se puede expresar en un tensor
de inercia como
-(u
+ t., + t.,
t.,
fu
m;i;
ly,
mJ
i; - i:
m.i
t.,
J
/JU
- I,, + i:
2
i:
t;
m.x,
my
l}C}C
2
m.i
m
(3.2-18)
(x;, y1,
- k, 1 + k22 + k33
kf,2
k,3
x,
k23
k, 1 + k22 - k33
kf,2
J;=m
k, 1 - k22 + k33
2
k,3
k23
t,
l
(3.2-19)
donde km es el radio de giro del elemento i respecto del eje yz e i' = (x1, Ji, i,,
J )7 es el vector centro de masa del elemento i desde el sistema de coordenadas del
elemento i-simo y expresado en el sistema de coordenadas del elemento i-si1110,
De aqu que la energa cintica total K de un brazo de robot es
(3.2-20)
DINAMICA
DEL BRAZO
DEL ROBOT
95
que es una cantidad escalar. Obsrvese que las J son dependientes de la distribucin de masa del elemento i y no de su posicin o velocidad de movimiento y se
expresan con respecto al sistema de coordenadas i-simo. Por tanto, la J;
se necesita calcular solamente una vez para obtener la energa cintica de un
robot.
3.2.3
1, 2, ... , n
(3.2-21)
y la energa potencial total del brazo se puede obtener sumando todas las
energas potenciales en cada elemento,
"
(3.2-22)
P =
P;
mg(A; ii')
i= 1
"
i= 1
= K -
P est
dada
por
1/2
m;g(A/r;)
i=I
j=l
k=I
(3.2-23)
i=l
T ~ :
(:~)
;~
==
I,
i=t
s=:
Tr(VikJiVJ)iik
I,
'f. 'f.
j=ik=lm=I
Tr(Vik,,.JiVJ)qkq,,. - I,migV/ri
i=t
(3.2-24)
k=l
DAk
L" L"
k=I
m=I
hikmikim
1, 2, ... , n
(3.2-25)
(3.2-26)
donde
r(t) = n x I vector par generalizado aplicado en las articulaciones i = 1, 2, ... ,
n; esto es,
(3.2-27)
q(t)
un vector n x
expresar como
q(t)
ij(t)
D(q)
o;
h(q. q)
j = nx (i. k)
donde
(3.2-30)
h,
(h1, h2,
..
l = 1 m= 1
hk,,.iki,,.
h.f
1, 2, ... , n
(3.2-32)
hikm
j=mhfi.k,m)
Tr (VjkmJJ Uj)
i, k, m
1, 2, ... , n
(3.2-33)
DINAMICA
c(q)
DEL BRAZO
DEL ROBOT
97
(-migVi/i)
donde
1, 2, ... , n
(3.2-34)
j=i
0(9)
o.,
(3.2-35)
donde
Tr(U11J1Ufi)
+ Tr(U21'J2UL)
+ Tr(U41J4U) + Tr(U51J5Uf)
D11
D12
D21 = Tr(U22J2UL)
+ Tr(U52J5U;)
Du
D14
Du
D,6
+ Tr(U31J3Uft) +
+ Tr(U61J6U)
+ tr(U32J3Uft) + Tr(U42.J4U) +
+ Tr(U62J6U)
D31 = Tr(U33J3Uf)
+ Tr (U63J6 U1)
D41
D51
= Tr(U55J5U;)
Tr(U44J4U)
+ Tr(U43J4U)
+ Tr(U53J5Uf1)
+ Tr(U54J5Ufi)
+ Tr(U64J6U1)
+ Tr(U6sJ6U1)
D61 = Tr (U66J6 U)
D22 = Tr (U22J2U{2)
Tr (U52J5U{2)
D23 = D32
+ Tr (U32J3Uf2)
Tr(U33J3Uf2)
Tr (U63J6 U2)
+ Tr (U42J,U2) +
Tr (U62J6U2)
+ Tr(U,3J,U2) + Tr(U53J5Uf2)
98
D36
=
=
Tr(U55J5U;2)
+Tr(U65J6U2)
D25
Ds2
D26
D62
D33
Tr(U33J3UL) + Tr(U43J4UL)
+ Tr (U63J6 U3)
D34
D43
D3s
D53
D63
D44
=
=
Tr(U66J6U2)
+ Tr(U53J5U;3)
Tr(U44J4UI3)
+ Tr(U54JsU;3)
Tr(U55J5U;3)
+ Tr(U65J6U3)
= Tr (U66J6 U3)
= Tr(U44J4UI4) +
+ Tr (U64J6 U3)
Tr(Us4JsU;4) + Tr(U64J6U4)
D45
D54
= Tr(UssJsU;4) + Tr(U65J6U4)
D46
D64
Tr (U66J6 U4)
o.,
Ds6
D6s
D66
Tr (U66J6 U6)
Tr (U66J6 Us)
H;.v
h23
h4 h24
h4
h;23 h24
h33 h34
h34
h44
h5
h25
hi35
h45
h55
hi56
hil6
hi26
h, 36
hi46
hi56
hi66
h;12
h;12
h22
h13
h3
h5
hil6
h25
hi26
h35
h45
hi36
hi46
= l , 2, ... , 6
(3.2-36)
la velocidad de las
(3.2-37)
7ff.
1,1)
(3.2-38)
se
DINAMICA
DEL BRAZO
DEL ROBOT
99
IJTH3.v
h(O, 9)
h3
h4
9TH4, v
hs
IJTHs,v
h6
{jTH6,v9
C2
= -(m1gU111r1
+ msgUs1 5i\
-(m2gU2/i'2
+ m6gU62 6i'6)
+
+
+
m2gU2/r2
m3gU3/i\
4i'4
m4gU424i'4
+ msgUs/rs +
+
5i's
+
msgU54
Cs
-(msgUs/i's
m6gU6s6r6)
-m6gU66
m4gU4/i'4
C4
+
+
c6
m3gU313i'3
m6gU616i'6)
-(m3gU3/r3
4r4
-(m4gU44
C3
m4gU43
msgUs3
5i's
t m6gU63 6i'6)
m6gU6/r6)
6r6
Los coeficientes C, Dik y hikm en las ecuaciones t).2-31) a ~3.234) son funciones
de ambas, las variables de articulacin y los parmetros inerciales del manipulador y algunas veces se llaman los coeficientes dinmicos del manipulador. El
significado fisico de estos coeficientes dinmicos se puede ver fcilmente de las
ecuacones de movimiento de Lagrange-Euler dadas por las ecuaciones (3.2-26) a
(3.2-34):
l. El. coeficiente c, representa los trminos de carga gravitatoria debido a los
elementos y se define por la ecuacin (3.2-34).
2. El coeficiente
est relacionado con la aceleracin de las variables de
articulacin y se define por la ecuacin (3.2-31). En particular, para i =
k,
D est relacionado con la aceleracin de la articulacin i donde el par
D,.
100
3.
DINAMICA
DEL BRAZO
m:1.
ROHHr
101
Multiplicaciones"
Adiciones
32n(n -
1)
24n(n -
4n(9n -
7)
51n - 45
n--2
f. m gU/rJ
Tr [U.Jt(UtYJ
(128/J)n(n
i'
Tr [UtJt(Ut;)
+ l)(n + 2)
1)
+ t )(n +
(65/2) n(n
2)
1)
h(i,j)
(128/3)n2(n
Tr [U.,1tJ .. (U,.
)7]
. . .i:
t
Tr (U,,.11J,..(U,,.;)7]
l>(q)q
b(q,
q) +
c(q)
+ l)(n + 2)
+ t )(n
(65/2) n2(n
t- 2)
( 1/6)n2(n
(128/3)n4 + (512/3)n3 +
+ (844/3)n2 + (76/3)n
(98/3)n4 + (781/6Jn1 +
+ (637/3)n2 + (107/6)n
- 1 )(n
1)
1
A causa de su estructura matricial, las ecuaciones de L-E son interesantc:s
desde d punto de vista de control en lazo cerrado porque dan un conjunto de
ecuaciones de estado como en la ecuacin (3.2-26). Esta forma permite el diseo
de una ley de control que compensa fcilmente todos los efectos no lineales.
Bastante a menudo en el diseo de un controlador por realimentacin para un
manipulador se utilizan los coeficientes dinmicos para minimizar los efectos no
lineales de las fuerzas de reaccin (Markiewicz [ 1973]).
Es de inters evaluar las complejidades computacionales inherentes a la ohtencin de coeficientes en las ecuaciones (3.2-31) a (3.2-34). La tabla 3. 1 resume las
complejidades computacionales de las ecuaciones de movimiento de L-E en
trminos de las operaciones matemticas requeridas (multiplicaciones y
sumas que se .necesitan para calcular la ecuacin (3.2-26) para cada punto de
consigna en la trayectoria.
Computacionalmente estas ecuaciones de
movimiento son ex- tremadamente ineficientes cuando se las compara con otras
formulaciones. En el siguien.teapartado desarrollaremos
las ecuaciones de
movimiento de un robot Que sern ms eficaces para calcular los pares
nominales.
[}~
[e,
'A, - ~'
A2
A11A2
-S,
OO /C,]
C1
/S1
o
[e "
S1 2
IA2
[e,o -s,
S2
o
-S12
C12
o
o
o
o
l
C2
OO IC,]
/S2
/(C., + C,)]
/(S12
o
1
S1)
DINAMlCA
DEL BRAZO
DEL ROBOT
donde C1 = cos 01; S1 = sen 01; Cii = cos (O; + O); Sii = sen (O;
definicin de la matriz Q;, para una articulacin giratoria, tenemos
o
Q
-1
O
[o
o
o
o
o
o0A1
O
-1
U, 1 =
= Q1
001
[-s,
Ai
o
o
.
[-s.,
0
o0A2
M = Q1 A2
e,
o
o
o
o
o
o
i c ']
o -
IS 1
c..
-IS]
1c
o
o
o O
o o o
o o o
o
o
o
o
-S12
o
o
-1
1
= ~
-C12
C ,2
[ o] [c.,
1
1
o
o
o
o
-S,
Anlogamente,
U21
o
o
o
o
-C,
-Si
e,
o
o
S1
Oi). De la
!]
o
o
o
o
-S12 o
S12
C12
o
o
C,)]
/(S12
+ S1)
=
,)] - l(S., + S
/(C12
+ C1)
o
o
U 22--- - o
A1Q2'A2-oA-2082
[e,
S1
- o
o
[-s.,
C ,2
-s, o
ic,]
e, o IS1
o 1 o
o o 1
-S12
-C12
o o
o
103
[
1
o
o
] [e,
-Si
S2
C2
o o o
o o o
o
o
o
o
-IS.,]
o
o o
o12 o
IC
-1O
o
o
o
o
1
ic,]
IS2
o
1
104
De la ecuacin (3.2-18), suponiendo que todos los productos de inercia son nulos,
podemos deducir la matriz de pseudoinercia J:
[ 'i,~,I'
o
- i/2m1/
-''f"]
o
o
o
o
o
o
o
o
J2
[ 'i,~,I'
o
_ 1/2m2'
m1
o
o
o
o
-'/t]
o
o
o
o
m2
= Tr(U11J1U1)
~Tr [-s,
+Tr \[CSf"
= i3m1/2
+ Tr(U21J2Uf) =
o
o
o
o
-C1
-Si
o
-S12
-r']
-IS] [/,m,1'
o o
/C1
o
o o
o o
oO
- i/2m1/
J(
e 1 2 + e i)
o
o
o
o
TI +
U11
m1
o o
-l(S.,+S,J r,m,I'
-C12
-'t]
- 1/2m2/
o
O
o o
O o
U21 -
m2
+ 4/3m2'2 + m2C2/2
1
Para D 12 tenemos
D12
D21
{[S "
Tr(U22J2Ufi)
-C12
Tr
o
C 12
-S12
o
o
o
o
o
-IS.,l [ 'f,m,I'
IC12
- i/2m2/
o o
o o
o o
o o
-'/r'l
O
m2
}U2T1
= m2/2( - i6
+ 1/2 + 1/2C2)
i3m2l2
+ i2m2l2C2
= Tr
{[s"
-C12
-S12
C12
o
o
o
o
/3m2'2Sf2
-IS.,] ['/,m,l'
o
o
o
o
/3m2'2Cf
IC12
-1/2m2'
'/3m2'2
o
o
o
o
o
o
o
o
-r'] }
o
m2
ur ""'
22
DINAMICA
DEL BRAZO
DEL ROBOT
IOS
h2
L L
= m=
l
hi.mlm
h211f + h21212
+ h221L + h222f =
Por tanto,
e,
-(m1gU11
(0,
-m1
ti\
O, O)
-g,
[-S,
C1
[-s.,
- mz(O, -g, O, O)
1/2m1g/C1
1/2
m2gU212r2)
C12
m2glC12
-C1
-Si
o
o
ol
oo
o
-IS']
o
o
-S12
o
o
+ m2glC1
o
o
-C12
o
o
/C1
-/(Su + s,)]
/(C12
o
o
C1)
--
o
o
106
-m2gU22 22 =
[-s,,
C,2
-C,2
-S,2
o
o
-IS,,]
/C12
-m2(0, -g, O, O)
o
o
o
o
o
o
= -mz(1/2g/C12 - g/C12)
Finalmente, las ecuaciones de movimiento de Lagrange-Euler para el manipulador de dos elementos se encuentra que son
-r(t)
D(8)6.(t)
h(8, )
c(8)
3.3 FORMULACION
DE NEWTON-EULER
D11
D12
D13
D14
D~s
'D12
D22
D23
D24
Dis
D26
116
D13
D23
D33
D34
D3s
D36
D14
D24
D34
D44
D4s
D46
D1s D16
Dis D26
D3s D36
D4s D46
Ds6
Ds6 D66