Sistema de Refrigeración
Sistema de Refrigeración
Sistema de Refrigeración
El doctorando El director
Impreso en Espana
Primera Edicion, noviembre 2015
A mis padres y mi hermano, que con su compresion y apoyo me
motivaron a culminar con exito este paso en mi proyecto de vida.
Agradecimientos
Muchas gracias,
Vctor Hugo Jaramillo Velasquez, noviembre 2015.
Indice general
Indice de figuras xi
Indice de tablas xv
1 Introduccion 1
1.1 Motivacion de la tesis . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 Objetivos de la tesis . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3 Organizacion de la tesis . . . . . . . . . . . . . . . . . . . . . . . 5
v
INDICE GENERAL
vi
INDICE GENERAL
vii
INDICE GENERAL
viii
INDICE GENERAL
ix
INDICE GENERAL
8 Apendices 195
8.1 Especificaciones de los componentes mecanicos . . . . . . . . . . 195
8.2 Especificaciones de los componentes electricos . . . . . . . . . . 199
8.3 Generacion de trayectorias . . . . . . . . . . . . . . . . . . . . . 205
Referencias 207
x
Indice de figuras
xi
INDICE DE FIGURAS
xii
INDICE DE FIGURAS
xiii
INDICE DE FIGURAS
xiv
Indice de tablas
xv
INDICE DE TABLAS
xvi
CAPITULO
1
Introduccion
Desde que el primer robot industrial fuera instalado en 1960, un gran numero
de aplicaciones han surgido en esta rama de la ingeniera. Un ejemplo de esto es
el ambito industrial, donde se pueden apreciar las tareas de soldadura, corte, pin-
tura, manipulacion entre otras. Aparte de estas tareas, se puede tomar en cuenta el
uso de los manipuladores roboticos en la industria aeroespacial, como el Canadarm
y el Canadarm2 de la Figura 1.1 y el Dextre [3], presentes en la Estacion Espa-
cial Internacional (en ingles ISS). Dichos robots son usados en tareas de ensamble,
transporte de equipos y mantenimiento durante las caminatas espaciales de la es-
tacion espacial internacional [4]. Estos robots al ser utilizados para aplicaciones
aeroespaciales poseen el peso mnimo necesario. Este factor, sumado a su parti-
cular esbeltez, hace que las deflexiones experimentadas en el robot sean bastante
apreciables. Desde este punto de vista se puede entonces catalogar a este tipo de
robots como Manipuladores Flexibles o Robot Flexibles. Esta flexibilidad se
refiere a flexibilidad de la estructura, y se puede considerar como concentrada en
las articulaciones o distribuida a lo largo de la estructura del robot (en particular en
los eslabones). En el desarrollo de esta tesis, se hara enfasis en los robot que poseen
eslabones flexibles.
El exito en el uso de los brazos roboticos en la industria radica en su velocidad
y precision. Un brazo robotico puede llevar a cabo una tarea 24 horas al da sin
1
1. INTRODUCCION
2
1.1 Motivacion de la tesis
3
1. INTRODUCCION
4
1.3 Organizacion de la tesis
5
1. INTRODUCCION
caciones obtenidas y se discute sobre los trabajos futuros que resultara interesante
continuar desarrollando.
6
CAPITULO
2
Descripcion de las
Plataformas Experimentales
7
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
grado de libertad (1GDL) y la mesa de aire, se puede apreciar en la Figura 2.1, por
otro lado los componentes mecanicos del robot de un grado de libertad se enumeran
a continuacion:
1 Mesa de aire
2 Carga util
3 Brazo flexible
4 Acoplamiento Barra-Motor
Figura 2.1: Apariencia fsica del Robot de 1GDL y su entorno de trabajo (tomado de
[2])
8
2.1 Descripcion del Robot de un grado de libertad
9
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
de un grado de libertad. Dichas cargas utiles no son mas que discos de poliamida
(nailon) cuyo peso ha sido calibrado y ajustado para poder ser utilizadas en la mesa
de aire. Dichos discos poseen un eje metalico que permite que sean conectados al
brazo, y poder as transmitirles el movimiento.
La Figura 2.3 muestra el aspecto de una de las cargas utiles tpicamente uti-
lizadas para los experimentos con el robot de un grado de libertad. En la Figura
tambien se puede apreciar como el brazo flexible se conecta con la carga util por
medio de una articulacion de rotula. Esto permite que el par de torsion ocasionado
por la inercia de rotacion de la carga util, no se transmita al brazo, de esta forma,
el brazo vela carga util (o masa del extremo) como una masa concentrada cuya
inercia de rotacion es cero.
10
2.1 Descripcion del Robot de un grado de libertad
estan fabricados en fibra de carbono [27]. Esto permite al robot levantarse a s mis-
mo (alrededor de 17 kg). En el caso del robot de un grado de libertad disponible en
el Laboratorio de Automatica y Robotica de la Universidad de Castilla-La Mancha
(UCLM), su eslabon tambien esta construido con materiales compuestos.
El robot de un grado de libertad tiene dos opciones de brazo flexible. La primera
opcion es una barra de fibra de vidrio con un metro de largo y un diametro de 3 mm.
La segunda opcion es una barra de fibra de carbono de 0,7m de largo y un diametro
de 2,8mm.
Para los experimentos desarrollados durante el trabajo de investigacion de esta
tesis se utilizo el brazo de fibra de carbono de 0,7m de largo. Las propiedades
mecanicas de dicho brazo se listan la Tabla 8.1.
11
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
amplificador del motor y el puente de galgas necesarios para el control del motor y
el sistema de medicion, respectivamente.
12
2.1 Descripcion del Robot de un grado de libertad
Los datos de velocidad y par estan dados con respecto al eje de salida de la reduc-
tora.
2.1.2.2 Servo-amplificador
Para poder controlar efectivamente la posicion del motor es necesario contar con
la electronica de potencia adecuada. Se usan servo-amplificadores especficamente
disenados para motores DC sin escobillas. La electronica de potencia del servo-
amplificador se encarga de enviar el tren de pulsos respectivo a los IGBT para
alternar la alimentacion a las fases del motor y poder as generar el movimiento. El
servo-amplificador utilizado es de la marca Maxon, modelo 4 Q DC Servo-
amplifier LSC30/2. La Tabla 8.9 del apendice indica los parametros electricos del
servo-amplificador.
13
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
2.1.3.1 Encoder
Para medir la posicion del motor se utiliza un encoder optico incremental que viene
incorporado en el motor. El encoder genera un tren de pulsos electricos que, depen-
diendo de la resolucion del encoder, permiten determinar cuanto ha rotado el eje
del motor. Al ser incremental, el encoder no mide la cantidad de giro respecto a un
punto de referencia fijo, sino que la posicion cero dependera de la posicion inicial
en la que comience a girar el motor. La Tabla 8.10 muestra las especificaciones del
encoder.
14
2.1 Descripcion del Robot de un grado de libertad
15
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
16
2.1 Descripcion del Robot de un grado de libertad
17
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
Amplificador de seal
Seal de las KYOWA
galgas
Servo-amplificador PC
Seal de corriente al del Motor
motor
Contador de Pulsos
del encoder
Seal del encoder
18
2.1 Descripcion del Robot de un grado de libertad
desarrollo de sistemas de control sin la necesidad de escribir codigo como tal, sino
a traves de diagramas de bloques. Y segundo, el uso del RTW permite la conver-
sion (y posterior compilacion) del modelo en Simulink
R y/o codigo de Matlab
R
19
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
Carga util
20
2.2 Descripcion del Robot de tres grado de libertad
21
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
22
2.2 Descripcion del Robot de tres grado de libertad
23
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
24
2.2 Descripcion del Robot de tres grado de libertad
Nodo 3
Barra 2 Barra 3
Nodo 2
Nodo 4
Barra 5
Nodo 1
motor que permite cambiar la elevacion (motor 2). La articulaciones del nodo 1 y
un par de sensores de carrera se pueden apreciar en la Figura 2.20.
La Tabla 8.5 resume las principales caractersticas de las barras y la Tabla 8.6
resume las principales caractersticas de las articulaciones pasivas del robot.
25
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
26
2.2 Descripcion del Robot de tres grado de libertad
27
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
2.2.2.2 Servo-amplificadores
Para poder controlar la posicion de cada uno de los motores, es necesario utilizar
los servo-amplificadores que permiten recibir la consigna de control y generar una
senal de control adecuada (corrientes) para ser enviada al motor. Ya que los tres
motores utilizados son de la misma referencia, los servo-amplificadores utilizados
tambien son del mismo modelo. Estos tambien son de la marca Baldor, referencia
BT S 10R200 2,5 (1993). La Tabla 8.15 contiene las caractersticas tecnicas
de los servo-amplificadores.
2.2.3.1 Encoders
Para garantizar el correcto control de movimiento del robot es necesario incorpo-
rar un encoder en cada motor. El encoder permite realimentar la posicion de cada
motor y de esta forma se puede cerrar el lazo de control de posicion, no solo de los
motores, sino tambien de la posicion del extremo del robot. Los 3 encoders utiliza-
dos son de la marca Heidenhain, referencia ROC 426 3600 03 (1994). La
Tabla 8.16 ilustra las caractersticas tecnicas de dichos encoders.
28
2.2 Descripcion del Robot de tres grado de libertad
29
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
Los sensores de final de carrera, aunque son elementos pasivos dentro del siste-
ma de control, garantizan la integridad fsica del robot en caso de un movimiento
no deseado. Estos sensores, son sensores inductivos ubicados en la base del robot.
Los sensores estan conectados a un sistema de compuertas logicas que, en cuanto
reciben senal, cortan la alimentacion del motor que se corresponde con el sensor
que envio la senal. Este es un mecanismo de proteccion muy pertinente durante
cualquier tipo de operaciones de prueba o funcionamiento normal del robot. Los
sensores de final de carrera utilizados fueron los sensores inductivos Telemecani-
que, referencia XS1N 08P B349.
30
2.2 Descripcion del Robot de tres grado de libertad
31
2. DESCRIPCION DE LAS PLATAFORMAS EXPERIMENTALES
Amplificador de las
seales
Seales de las de las galgas
galgas
Seales de los PC
encoders
32
CAPITULO
3
Modelado de Manipuladores
Flexibles
3.1 Introduccion
Para comprender el comportamiento dinamico de los robots flexibles, as como para
el desarrollo de sistemas de control para los mismos, es necesario el desarrollo de
modelos matematicos que describan el sistema de la mejor forma posible. Este
requisito fue identificado por autores como Usoro [28] o Cannon[29] durante sus
primeras investigaciones con manipuladores flexibles, en las cuales se percataron
de la importancia de un modelo lo mas exacto posible para garantizar un adecuado
control.
La correcta estimacion de los estados en un manipulador flexible depende en
gran medida de la precision del modelo utilizado, ya que este tiene que garanti-
zar una completa caracterizacion de los estados internos del sistema flexible. Para
los manipuladores roboticos tradicionales (o robots industriales) los estados inter-
nos asociados con las deflexiones de los eslabones pueden ser despreciados en la
mayora de las aplicaciones, ya que las velocidades son bajas y la rigidez del ma-
nipulador es alta.
33
3. MODELADO DE MANIPULADORES FLEXIBLES
34
3.2 Modelado Dinamico de Robots Flexibles
35
3. MODELADO DE MANIPULADORES FLEXIBLES
metodo de Expansion Modal (tambien conocido como metodo de los modos asu-
midos - en ingles Assumed Modes Method o AMM), el metodo FEM y el metodo de
los parametros concentrados (que es una particularizacion del FEM, considerando
elementos tipo viga) - ver la Figura 3.1. Esta popularizacion de algunos metodos
esta evidenciada en el trabajo de Dwivedy y Eberhard [40], un artculo dedicado a
resenar el estado del arte en el modelado de manipuladores flexibles desde su apari-
cion en el ambito cientfico hasta el 2006. En [40], los autores agrupan los trabajos
de diferentes autores dependiendo de la tecnica utilizada para el modelado, que-
dando entonces tres grupos significativos de autores, los enfocados en AMM, los
enfocados en FEM y los enfocados en el metodo de los parametros concentrados.
36
3.2 Modelado Dinamico de Robots Flexibles
37
3. MODELADO DE MANIPULADORES FLEXIBLES
numero extenso de elementos esto puede ser un proceso tedioso. Sin embargo, al
final es mas facil incorporar las no-linealidades sin complicar el modelo basico
[30].
38
3.2 Modelado Dinamico de Robots Flexibles
39
3. MODELADO DE MANIPULADORES FLEXIBLES
40
3.2 Modelado Dinamico de Robots Flexibles
permite conocer todas las frecuencias y formas modales del manipulador en cual-
quier configuracion del mismo. Esta informacion es posteriormente usada para su
control con un controlador LQR. Dong et al. desarrollaron en [84] un modelo para
un manipulador flexible de 3 GDL usando formulacion Lagrangiana. El robot en
cuestion posea eslabones y articulaciones flexibles. En ultima instancia, los auto-
res determinaron que la baja frecuencia de vibracion causada por la flexibilidad en
las articulaciones era la que mas impacto tena en el comportamiento dinamico del
manipulador. Qian et al. [85] parten del trabajo de Book [70] y formulan un metodo
recursivo para el modelado dinamico de manipuladores flexibles de multiples gra-
dos de libertad espaciales que, a diferencia del metodo inicialmente propuesto por
Book, incluye las deformaciones torsionales y considera la flexibilidad de las arti-
culaciones al igual que el aumento de la rigidez debido a las aceleraciones. Estas
consideraciones adicionales le brindan mayor fidelidad al modelo.
Para efectos practicos, en la mayora de las aplicaciones, el metodo propuesto
por Book en [70] para modelado de manipuladores flexibles de multiples grados
de libertad es valido, sin restarle merito al trabajo desarrollado por Qian et al. [85].
Es por esto que, a pesar de los diversos metodos propuestos para el modelado de
robots flexibles, el trabajo de Book en [70] ha sido uno de los mas difundidos para
el modelado de manipuladores flexibles. Dicho metodo se caracteriza por varios
aspectos: 1) es eficiente computacionalmente, 2) es una formulacion completa y
3) es un metodo simple y directo. Este metodo utiliza matrices de transformacion
4 4 (similares a las utilizadas en robotica) para representar el movimiento, tanto
el de las articulaciones como el debido a las deflexiones de los eslabones. Gracias a
la formulacion Lagrangiana y a la naturaleza recursiva del metodo, las ecuaciones
resultantes no ignoran la interaccion entre las velocidades angulares y las deflexio-
nes y permiten llevar a cabo los calculos de forma eficiente. Adicionalmente, las
ecuaciones asumen pequenas deformaciones de los eslabones, las cuales son des-
critas como la suma de las formas modales considerando un modelo elastico lineal.
Las principales desventajas de este metodo es que para el modelado solo se pue-
den considerar articulaciones rotacionales y unicamente es util para determinar la
dinamica de manipuladores flexibles cuya configuracion sea una cadena cinemati-
ca abierta. Un ejemplo claro de la aplicacion de esta tecnica de modelado se puede
apreciar en el trabajo de Yoshikawa et al. [86].
41
3. MODELADO DE MANIPULADORES FLEXIBLES
42
3.3 Modelos Matematicos de las Plataformas Utilizadas
mb la masa de la barra,
43
3. MODELADO DE MANIPULADORES FLEXIBLES
Jt, mt
Eb, Ib, Lb r
mb
t
c
m, , Jm
m
X
Figura 3.2: Robot Flexible 1GDL
Una forma de obtener el modelo del robot, es a traves del analisis Lagrangiano.
Este metodo parte del analisis de las energas cinetica y potencial, y permite obtener
el modelo del sistema de forma rapida. Primero, se define el Lagrangiano L como
L=T V (3.6)
44
3.3 Modelos Matematicos de las Plataformas Utilizadas
45
3. MODELADO DE MANIPULADORES FLEXIBLES
L2b 2
1 2 2
2 2
T = mb c + mt Lb t + Jt t 2t r + r (3.14)
2 4
Por otro lado, la energa potencial elastica se define como
1 L2b
1 2 1 16 2
V = kb (m t ) + kb (m c ) + kb (t r )2 (3.15)
2 2 5 2 3
3Eb Ib
kb = . (3.16)
Lb
Una vez se conocen los terminos del Lagrangiano se definen las coordenadas
generalizadas que se utilizaran para determinar la dinamica. Para el caso del Robot
de 1GDL, las coordenadas generalizadas son m , t , c y r . Se tiene entonces que
d L L
m =
dt m m
d L L
0=
dt c c
d L L
0=
dt r r
d L L
0= (3.17)
dt t t
46
3.3 Modelos Matematicos de las Plataformas Utilizadas
La ecuacion (3.18) permite ver como la dinamica del brazo afecta el par motor m .
Por otro lado, las ecuaciones (3.19)-(3.21) permiten obtener la relacion dinamica
entre el angulo del motor m y el angulo del extremo t . Aunque este modelo con-
sidera la masa de la barra y la inercia de rotacion de la carga util, estos valores
son generalmente muy pequenos en comparacion con la carga util. Si se considera
Jt 0 y mt mb , la funcion de transferencia que se obtiene entre el par motor
m y el angulo del extremo a partir de (3.19)-(3.21) es:
mt L2b
m = t + t (3.22)
kb
Esta dinamica se asemeja a la de un sistema masa resorte
q de un grado de libertad,
kb
en el que la frecuencia natural b se define como b = mt L2b
.
El modelo obtenido es equivalente al modelo que se habra obtenido utilizando
la tecnica de los Modos Asumidos (AMM) y posteriormente truncando el sistema a
un modo de vibracion. En situaciones en las que la masa de la carga util del robot no
es significativamente mayor a la masa de la barra, el modelo adquirira complejidad
y no sera tan facil el despreciar los efectos de modos de vibracion de orden superior.
El modelado del robot de 1GDL tambien se puede llevar a cabo con un analisis
Newtoniano. Sin embargo, para ser congruentes con el proceso llevado a cabo para
el robot de 3GDL (que se vera en la siguiente seccion) se decidio modelar el de
1GDL por medio de la mecanica Lagrangiana.
La validacion de este modelo no era parte esencial del desarrollo de esta tesis
doctoral, esto sumado al hecho de que este tipo de modelos ya han demostrado ser
lo suficientemente validos por otros autores, fue considerado como razon suficiente
para no incorporar un analisis de la validacion del modelo en este documento.
L1 , L2 , ..., L5 : son las longitudes de la barra 1, barra 2, ..., barra 5, del robot.
47
3. MODELADO DE MANIPULADORES FLEXIBLES
nodo 3
L2 L3
nodo 4
nodo 2
L5 F
Z
L1 3 L4 F
nodo 5
2
R
nodo 1 F
Y R zR
R
xR X
yR
Figura 3.3: Robot flexible
F , F , F : son las fuerzas que actuan en el extremo del robot en las direc-
ciones de las coordenadas esfericas.
En coordenadas esfericas, la posicion del extremo del robot (como si fuera rgido)
se define como:
Pesf = R r (3.23)
48
3.3 Modelos Matematicos de las Plataformas Utilizadas
donde
r = cos(R ) cos(R )i + sin(R ) cos(R )j + sin(R )k (3.24)
y
p
R = x2 + y 2 + z 2
y
R = arctan (3.25)
x !
z
R = arctan p
x2 + y 2
+ cos()
Pesf = r + (3.27)
donde
La energa cinetica de una partcula que esta localizada en el extremo del robot, se
define como:
1
T = mPesf PT
esf (3.30)
2
49
3. MODELADO DE MANIPULADORES FLEXIBLES
VG = mg sin() (3.32)
50
3.3 Modelos Matematicos de las Plataformas Utilizadas
donde
1 (L5 )3
= (3.38)
K5f lex 3EI
51
3. MODELADO DE MANIPULADORES FLEXIBLES
52
3.3 Modelos Matematicos de las Plataformas Utilizadas
5 = 5 cos()
5 = 5 sin() (3.47)
Esto permite expresar las deflexiones del extremo del robot en coordenadas esferi-
cas.
Para obtener las deflexiones del nodo 4 se procede de la siguiente forma. Se
desplazan las fuerzas F y F del nodo 5 al nodo 4. De igual forma que para la
barra 5, las fuerzas F y F se pueden descomponer en un par de componentes,
una paralela a la barra 4 (o fuerza Fx4 ), y la otra perpendicular a la barra 4 (o fuerza
Fy4 ). Dichas fuerzas se definen como:
F4y
4 = (3.49)
K4f lex
53
3. MODELADO DE MANIPULADORES FLEXIBLES
donde
1 (L4 )3
= (3.50)
K4f lex 3EI
La deformacion causada por la fuerza F4x es una deformacion axial y su mag-
nitud se considera despreciable en comparacion con la deformacion 4 .
4 = 4 cos()
4 = 4 sin() (3.51)
tip = 4 + 5
tip = 4 + 5 (3.52)
Ahora, para obtener los desplazamientos absolutos en funcion de las fuerzas del
extremo hay que sustituir las ecuaciones (3.51), (3.47), (3.48) y (3.35) en (3.52),
quedando as:
54
3.3 Modelos Matematicos de las Plataformas Utilizadas
sin sin
tip = [F cos() F sin()] + [F sin() + F cos()] (3.53)
K4f lex K5y
cos cos
tip = [F cos() F sin()] f lex + [F sin() + F cos()] y (3.54)
K4 K5
55
3. MODELADO DE MANIPULADORES FLEXIBLES
La magnitud de dichos pares esta en funcion del angulo . La deflexion del nodo 5
en la direccion dependera de la deflexion torsional de la barra 4 ocasionada por el
par x4 y de la deflexion a flexion de la barra 5 ocasionada por la fuerza F . As se
tiene que:
L35 L5
5 = F + F L5 cos( /2) cos( /2) (3.59)
3EI JP G
1 1
5 = F + F (3.60)
K5f lex K5T orsBarra4
siendo
1 L4
= L25 cos2 ( /2) (3.61)
K5T orsBarra4 JP G
56
3.3 Modelos Matematicos de las Plataformas Utilizadas
1
5 = F (3.62)
K5
1 K5T orsBarra4 + K5f lex
=
K5 K5T orsBarra4 K5f lex
F F
4 = + (3.63)
K4f lex K4f
donde
1 L5 L24 sin( /2)
= (3.64)
K4f 2EI
Que finalmente se puede simplificar a:
F
4 = (3.65)
K4
siendo
1 1
tip = F
+ F (3.66)
K4 K5
57
3. MODELADO DE MANIPULADORES FLEXIBLES
K4 + K5
C33 = (3.68)
K4 K5
donde:
58
3.3 Modelos Matematicos de las Plataformas Utilizadas
d L L
= (3.73)
dt R R
donde, F , y son las fuerzas y pares aplicados en el robot y que generan
movimiento en las coordenadas R , R y R respectivamente, como se muestra en
la Figura 3.3. Desarrollando las anteriores expresiones se obtiene:
1 K11
F = 2K11 (R ) + (R )2 + 2K12 (R )
2 R
K12 K22 2
+2 (R ) (R ) + (R )2
R R
K 33
+ 2 cos2 () (R )2 (3.74)
R
59
3. MODELADO DE MANIPULADORES FLEXIBLES
J1 1 = 1 1 + m1 L1
J2 2 = 2 2 + m2 L2 (3.85)
J3 3 = 3 3 + m3 L3
donde Ji representa el momento de inercia del eje del motor i, i la friccion visco-
sa del motor i y mi el par motor del motor i. Considerando las ecuaciones de la
60
3.3 Modelos Matematicos de las Plataformas Utilizadas
dinamica de los motores y las ecuaciones dinamicas del robot, es posible encontrar
la relacion entre las entradas del sistema (pares) y las salidas (angulos de los moto-
res y posiciones del extremo del robot). Combinando las ecuaciones (3.74), (3.75),
(3.76), (3.80), (3.81), (3.82), (3.84) y (3.85) se define la dinamica total, quedando
un sistema no lineal expresado como x = f (x, u). Donde el vector de estados x
esta constituido por
x= (3.86)
1
1
2
2
3
3
61
3. MODELADO DE MANIPULADORES FLEXIBLES
donde
2 2 2
K11 K21
f1 (x) = cos () + + (R ) (R )(R )
m m
K21 K22
+ (R ) (R )2 cos2 ()K33 (R )2
m m m
c
g sin() (3.89)
m
2 K21 K22
f2 (x) = cos() sin()2 + (R ) + (R )
m m
K33 g c
+ cos() sin()(R )2 cos() (3.90)
m m2
2 K33 c
f3 (x) = 2 tan() + (R ) (3.91)
m m2 cos2 ()
(3.92)
Por simplicidad, en las anteriores expresiones las coordenadas rgidas del robot
R , R y R se dejan expresadas como tal. Sin embargo dichas coordenadas se
definen en funcion de los angulos de los motores por medio de la transformada
cinematica directa como
p
R = 2L2 (1 cos ()) (3.93)
R = 2 3 + (3.94)
2 2
R = 1 (3.95)
L = L4 = L5 (3.96)
p
L6 = 2L2 (1 cos (3 )) (3.97)
62
3.3 Modelos Matematicos de las Plataformas Utilizadas
que:
x1
x2
x3
x= = (3.98)
x4
x5
x6
R u1
u = R = u2 (3.99)
R u3
K11 K21
f1 (x, u) = x2 cos2 (x4 )x25 + x23 + (u1 x2 ) (u2 x4 )(u1 x2 )
m m
K21 K22 x2
+ x2 (u2 x4 ) x2 (u2 x4 )2 cos2 (x4 )K33 (u3 x6 )2
m m m
cx2
g sin(x4 ) x1 (3.101)
m
2 K21 K22
f2 (x, u) = x1 x3 cos(x4 ) sin(x4 )x25 + (u1 x2 ) + (u2 x4 )
x2 mx2 m
K33 g cx
+ cos(x4 ) sin(x4 )(u3 x6 )2 cos(x4 ) 42 x3 (3.102)
m x2 mx2
2x1 x5 K33 cx6 x 5
f3 (x, u) = 2 tan(x4 )x5 x3 + (u3 x6 ) 2
(3.103)
x2 m mx2 cos2 (x4 )
(3.104)
De igual forma se pueden escribir las ecuaciones de las salidas del sistema en fun-
cion de los estados. Teniendo en cuenta las galgas extensiometricas instaladas en
la estructura del robot, se consideraran las senales de estas como las salidas del
63
3. MODELADO DE MANIPULADORES FLEXIBLES
sistema 4 , 5 , 4 , 5 , 1 , as:
1
4 = [mf2 (x)x2 cos((u1 )) mf1 (x) sin((u1 ))] (3.105)
K4f lex
1 1
5 = + [mf2 (x)x2 cos((u1 )) + mf1 (x) sin((u1 )] (3.106)
K4f lex KR (u1 )
1 1
4 = + mf3 (x)x2 cos(x4 ) (3.107)
K4f lex K4f (u1 )
1 1
5 = + mf3 (x)x2 cos(x4 ) (3.108)
K5f lex K5T orsBarra4 (u1 )
1 = K1f (u1 ) [mf2 (x)x2 cos((u1 )) + mf1 (x) sin((u1 )] (3.109)
donde
64
3.3 Modelos Matematicos de las Plataformas Utilizadas
1/G(s) coup
Posicin
ref + mot real robot
Dinmica
C(s) G(s)
del brazo
65
3. MODELADO DE MANIPULADORES FLEXIBLES
Xi = 0,5222 [m]
Yi = 0 [m]
Zi = 0,3077 [m]
donde los ejes cartesianos estan definidos segun la Figura 3.3. Para el caso de la
trayectoria circular, el punto final coincide con el punto inicial. Para el caso de la
trayectoria rectilnea, el punto final tiene la coordenada cartesiana:
Xf = 0,3573 [m]
Yf = 0,4669 [m]
Zf = 0 [m]
66
3.3 Modelos Matematicos de las Plataformas Utilizadas
Tipo de trayectoria
Rectilnea Circular
Xi = 0.5222 Xi = 0.5222
Posicion Inicial [m] Yi = 0.0 Yi = 0.0
Zi = 0.3076 Zi = 0.3076
Xf = 0.3573 Xf = 0.5222
Posicion final [m] Yf = 0.4669 Yi = 0.0
Zf = 0 Zi = 0.3077
c = 0,07
c = 0,12
c = 0,09
Para poder ejecutar las simulaciones era necesario conocer las condiciones ini-
ciales del sistema, estas condiciones iniciales incluyen:
67
3. MODELADO DE MANIPULADORES FLEXIBLES
Para las simulaciones (y para los experimentos) se considera que el robot parte
del reposo, en el instante inicial (t0 ) la unica fuerza que actua sobre el extremo del
robot es el peso de la carga util en la direccion del eje z, y ya que R = = 0 la
ecuacion (3.117) se hace igual a cero y las expresiones (3.115) y (3.116) se reducen
a
Para estimar las posiciones flexibles iniciales del robot, se resuelve el conjunto
de ecuaciones (3.118) y (3.119). Para resolverlas es necesario conocer los valores
de R y R que son los valores de la posicion del extremo del robot considerando
que es rgido. Estos valores dependeran de la trayectoria a describir por el robot.
A continuacion se muestran las graficas comparativas entre los resultados de las
simulaciones y los resultados de los experimentos realizados con el robot flexible.
Las Figuras 3.9 y 3.10 muestran los resultados para las trayectorias circulares, y las
Figuras 3.11 y 3.12 muestran los resultados para las trayectorias rectilneas.
68
3.3 Modelos Matematicos de las Plataformas Utilizadas
0.6 0.5
Simulation Simulacion Simulation
Optotrack Optotrack Optotrack
0.54 0.45
0.5
0.4
0.53
0.4 0.35
0.52
0.3
X (m)
Y (m)
Z (m)
0.3
0.51 0.25
0.2
0.2
0.5
0.15
0.1
0.49
0.1
0
0.48
0.05
0 1 2 3 4 5 1 2 3 4 5 0 1 2 3 4 5
t(s) t (s) t (s)
0.53
Simulation Simulacion Simulation
Optotrack Optotrack 0.5 Optotrack
0.52 0.6
0.51
0.5
0.4
0.5
0.4
0.49 0.3
X (m)
Y (m)
Z (m)
0.3
0.48
0.2 0.2
0.47
0.1
0.46
0.1
0.45 0
0.44 0.1 0
0 1 2 3 4 5 0 1 2 3 4 5 0 1 2 3 4 5
t(s) t (s) t (s)
0.52 0.4
Simulation Simulacion Simulation
Optotrack Optotrack Optotrack
0.51
0.35 0.25
0.5
0.3
0.2
0.49
0.25
0.48 0.15
X (m)
Y (m)
Z (m)
0.2
0.47 0.1
0.46 0.15
0.05
0.45 0.1
0
0.44
0.05
0.43 0.05
0
0 1 2 3 4 5 0 1 2 3 4 5 0 1 2 3 4 5
t(s) t (s) t (s)
69
3. MODELADO DE MANIPULADORES FLEXIBLES
0.52 0.4
Simulation Simulacion 0.3 Simulation
Optotrack Optotrack Optotrack
0.51 0.35
0.25
0.5
0.3
0.2
0.49
0.25
0.15
0.48
X (m)
Y (m)
Z (m)
0.2
0.1
0.47
0.15
0.05
0.46
0.1 0
0.45
0.43 0 0.1
0 1 2 3 4 5 0 1 2 3 4 5 0 1 2 3 4 5
t(s) t (s) t (s)
Crculo Crculo
Tipo de Error
con 1,2 [kg] con 1,6 [kg]
X Y Z Total X Y Z Total
Error max. [m] 0,0083 0,0814 0,0327 0,0881 0,0099 0,0911 0,038 0,0992
Error abs.
0,0014 0,0104 0,0052 0,0117 0,0009 0,0109 0,0058 0,0123
promedio [m]
Error relativo
0,3 16,3905 1,1383 16,4327 0,1969 12,7678 1,1341 12,8196
promedio %
70
3.4 Conclusiones
Lnea Lnea
Tipo de Error
con 1,2 [kg] con 1,6 [kg]
X Y Z Total X Y Z Total
Error max. [m] 0,0047 0,0106 0,0117 0,0881 0,005 0,0814 0,0327 0,0163
Error abs.
0,0010 0,0017 0,0016 0,0024 0,0010 0,0016 0,0052 0,0024
promedio [m]
Error relativo
0,1497 1,4508 0,4882 1,538 0,1424 1,2411 1,1383 1,3427
promedio %
3.4 Conclusiones
Se han desarrollado los modelos analticos para los manipuladores de 1GDL y 3
GDL mediante la formulacion de Euler-Lagrange. El modelo del manipulador de
71
3. MODELADO DE MANIPULADORES FLEXIBLES
3GDL esta simplificado en el sentido que se ha considerado que la masa del mani-
pulador esta concentrada en el extremo del mismo. De esta forma el sistema solo
posee 3 modos de vibracion y 3 frecuencias naturales. De esta forma se pretende
obtener un modelo de facil implementacion y de bajo costo computacional, que a
la vez capture la mayor parte de la dinamica del sistema en cuestion. La validez del
modelo del robot de 3GDL se verifico experimentalmente demostrando un nivel de
ajuste aceptable.
72
CAPITULO
4
Mejoras en el Diseno y
Control de Actuadores para
Manipuladores Flexibles
73
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
r ua u0 u y
C(s) G(s)
Observador
de
Perturbaciones
74
4.1 Introduccion a los observadores de perturbaciones
x = Ax + bu + ed (4.1)
y = cx (4.2)
A = A0 + A (4.4)
b = b0 + b. (4.5)
75
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
De esta forma las ecuaciones de espacio de estados (4.1) se pueden expresar como:
x = A0 x + b0 u + ed (4.6)
x = A0 x + b0 u
y = c0 x (4.7)
h iT
donde x = d,
d,
xT , d, . . . , dp1 R(n+p)1 ,
A0 e 0n(p1)
A0 = R(n+p)+(n+p) ; (4.8)
0pn 0p1 0pp1
T
R(n+p)1 ;
b0 = bT0 01p (4.9)
c0 = c0 01p R1(n+p) .
(4.10)
76
4.2 Introduccion a los reguladores fraccionarios
n Z t
1 d f (x)
a Dt f (t) = dx, n 1 < < n, (4.14)
(n ) dt a (t x)n+1
77
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
78
4.3 Control robusto de actuadores ante rozamiento y grandes cambios en la
carga
propuesto por Oustaloup con el objetivo de alcanzar lo que el llamo Robustez frac-
cionaria [110, 111], una propiedad que establece que el amortiguamiento de un
sistema esta exclusivamente ligado al orden fraccionario de la ecuacion diferen-
cial que define dicho sistema [112]. El controlador CRONE ha demostrado ser
una herramienta util para el control robusto y ha encontrado multiples aplicaciones
reales, como el control de la suspension de coches [111], transmisiones flexibles
[110] y actuadores hidraulicos [113]. Incluso se ha desarrollado una Toolbox de
Matlab [114] para su empleo. Por ultimo, las tpicas redes de adelanto/retraso tam-
bien se han replanteado desde la perspectiva del calculo fraccionario. Incluso se
han desarrollado tecnicas de autosintonizacion de estas redes de adelanto/retraso
fraccionarios [115].
Entre las diversas aplicaciones que ha encontrado el control fraccionario en
problemas de ingeniera, uno de gran interes es el area de control de movimiento.
Gracias a las caractersticas de robustez de los controladores fraccionarios, estos
son idoneos en tareas de control de movimiento, tal y como se discute en [116, 117].
79
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
Km V = JL + f (m , V ) + m (4.17)
80
4.3 Control robusto de actuadores ante rozamiento y grandes cambios en la
carga
F(s)
d
r + ua u0
+ u y
C(s) G(s)
Observador
de
Perturbaciones
d
JL0/Km0
81
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
la dinamica del motor y ii) las perturbaciones ocasionadas por las variaciones en
los parametros nominales del motor. Estas las ultimas pueden ser definidas como:
JL0 = Km0 V + JL0 d, (4.21)
1
d = (f + J Km0 V ) (4.22)
JL0
Considerando que la perturbacion se puede asumir como una variable constante
a tramos, se tendra que la primera derivada de d es cero y, haciendo uso del conjunto
de ecuaciones (4.7), se puede reescribir el sistema en espacio de estados aumentado
como:
0 1 0 0
d Km0
= 0 0 1 + Jm0
V (4.23)
dt
d 0 0 0 d 0
y = 1 0 0 (4.24)
d
Cuando se puede medir al menos uno de los estados del sistema, es posible
utilizar un observador de orden reducido en vez de un observador de completo,
algo que permite reducir la carga computacional del observador. Para el caso del
sistema analizado, se puede medir la posicion por medio de un encoder. El diagra-
ma de bloques de la Figura 4.3 ilustra la forma del observador de orden reducido
implementado. En la Figura 4.3 se definen:
f11 f12 g11
F= , g= (4.25)
f21 f22 g21
t11 1 0
T= , C = c0 , I = Tb0 (4.26)
t21 0 1
82
4.3 Control robusto de actuadores ante rozamiento y grandes cambios en la
carga
V z z
1
M ++
+ . [ ]
C
T
F d
As se tiene que:
0 1 0
f11 f12 t11 1 0 t11 1 0 0 0 1 g11
= 1 0 0
f21 f22 t21 0 1 t21 0 1 g21
0 0 0
(4.28)
donde:
2
f11 = t11 , g11 = (f11 + f21 )
f12 = 1
f22 = 0 (4.29)
83
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
(s + )2 = s2 + 2s + 2 (4.35)
Con la ecuacion (4.36) se ubicaran los dos polos del observador en el mismo
lugar del eje real del semi-plano negativo. Utilizando las variables de estado del
84
4.3 Control robusto de actuadores ante rozamiento y grandes cambios en la
carga
JL0
V = V0 d (4.38)
Km0
Km0
= (4.39)
V0 JL0 s2
4.3.2.2 Pre-alimentacion
Como se aprecia en la Figura 4.2, existe un termino F (s) que representa la funcion
de transferencia de la pre-alimentacion. Para el caso del esquema de control pro-
puesto, F (s) es equivalente a la dinamica inversa del sistema compensado. De esta
forma se puede hacer uso de la ecuacion (4.39) para obtener:
Jm0 s2
F (s) = (4.40)
Km0
Este termino adicional permitira al sistema seguir con mayor precision las trayec-
torias al disminuir el tiempo de asentamiento [120].
85
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
del sistema por medio de una fase constante, es decir, a pesar de las posible incer-
tidumbres en el modelo del sistema, la fase del controlador se mantiene constante
en un ancho de banda determinado [121]. La ecuacion de este controlador es:
C(s) = Ks (4.41)
En [122, 123] se han establecido criterios para la relacion entre los parametros
temporales y frecuenciales de forma que, de ser necesario, los controladores frac-
cionarios se pueden sintonizar considerando especificaciones temporales o frecuen-
ciales. Para la sintonizacion en frecuencia del regulador CRONE de 1ra generacion,
se necesitan dos parametros:
2m
= (4.42)
c2 JL0
K = (4.43)
Km0
Este tipo de controlador es util cuando la planta a controlar posee fase constante,
al menos en un rango de frecuencias alrededor de la frecuencia de corte. En ese
caso, el lazo de control sera robusto a incertidumbres en la planta, y aun cuando
el margen de ganancia cambie, el margen de fase de la planta y el del controlador
permaneceran constantes.
86
4.4 Ejemplo de aplicacion
JL0 c2 cos(m )
KP = (4.44)
Kt0
JL0 c sin(m )
KD = (4.45)
Kt0
87
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
Margen de Fase (m ) 80 [o ]
Frecuencia de cruce (c ) 50 [rad/s]
Ganancia (K) 1.42 Ecuacion (4.43)
Orden fraccionario () 0.8889 Ecuacion (4.42)
Ganancia Proporcional (KP ) 15.2034 Ecuacion (4.44)
Ganancia Derivativa (KD ) 1.7245 Ecuacion (4.45)
Tiempo de muestreo (Ts ) 0.001 [s] -
Ubicacion de los polos del observador -200 Ecuacion (4.46)
88
4.4 Ejemplo de aplicacion
1.2
0.8
Posicin angular [rad]
0.6
0.4
0.2
Ks + Obsv.
0 PD + Obsv.
Referencia
-0.2
0 0.5 1 1.5 2 2.5 3
Tiempo [s]
se puede decir de la senal de control, como se ve en la Figura 4.5. Las Figuras 4.6 y
4.7 muestran la estimacion de la perturbacion para ambos sistemas. Se puede apre-
ciar que en ambos casos el observador de perturbaciones estima correctamente el
valor de la perturbacion durante la trayectoria, lo que permite alcanzar la referen-
cia de posicion sin problemas. En estas figuras se puede apreciar una oscilacion de
alta frecuencia despues de un segundo de trayectoria, estas oscilaciones se deben
al fenonemo llamado en ingles hunting y que resulta de los cambios de friccion
estatica a friccion de Coulomb y viceversa [125].
89
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
PD + Obsv.
1.5 Ks + Obsv.
1
Seal de control [V]
0.5
-0.5
-1
0 0.5 1 1.5 2 2.5 3
Tiempo [s]
4 4
Perturbacin [m/s2]
Perturbacin [m/s2]
2 2
0 0
-2 -2
Error absoluto
Error absoluto
4 4
2 2
0 0
0 0.5 1 1.5 2 2.5 3 0 0.5 1 1.5 2 2.5 3
Tiempo [s] Tiempo [s]
90
4.4 Ejemplo de aplicacion
1 1
Posicin Angular [rad]
0.8 0.8
0.6 0.6
0.4 0.4
1 1
Posicin Angular [rad]
0.8 0.8
0.6 0.6
0.4 0.4
Figura 4.8: Posicion angular del motor para diferentes valores de inercia de rotacion
de la carga util
91
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
1.2
1.0
0.8
Posicin angular [rad]
0.6
0.4
0.2
Ks + Obsv.
0.0 PD + Obsv.
Referencia
-0.2
0 0.5 1 1.5 2 2.5
Tiempo [s]
92
4.4 Ejemplo de aplicacion
1.0 1.0
0.8 0.8
Posicin angular [rad]
0.4 0.4
J=15J025J0 J=15J025J0
0.2 Ks+Obsv. 0.2 Ks+Obsv.
PD+Obsv. PD+Obsv.
Referencia Referencia
0 0
0 0.5 1.0 1.5 2.0 2.5 3.0 0 0.5 1.0 1.5 2.0 2.5 3.0
Tiempo [s] Tiempo [s]
Cambio en t = 2.0 [s]
1.0
0.8
Posicin angular [rad]
0.6
0.4
J=15J025J0
0.2 Ks+Obsv.
PD+Obsv.
Referencia
0
0 1.0 2.0 3.0 4.0 5.0 6.0
Tiempo [s]
Figura 4.10: Respuesta del sistema para cambios repentinos en la carga inercial (de
JL = 15JL0 to JL = 25JL0 )
extremo del eje del motor. En este escenario se somete el sistema a dos pares ex-
ternos a la dinamica del motor, el par de friccion y el par debido al peso del brazo
rgido movido por el motor. La Figura 4.11 ilustra el brazo rgido y la fuerza de
gravedad actuando en su centro de masa. El par que perturba el sistema d debido
a los efectos de la gravedad se puede definir con la siguiente expresion:
L
cos(). d = mg (4.47)
2
Ahora, suponiendo una trayectoria entre 0 y rad y una inercia igual a la iner-
cia nominal, la respuesta del sistema es la que se ilustra en la Figura 4.12. En esta
figura se puede apreciar una ligera sobreoscilacion en la respuesta del sistema con-
trolado con el regulador Ks . A pesar de que dicha respuesta esta dentro del ancho
93
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
L/2
mg
d
Figura 4.11: Motor con brazo rgido
de banda del error admisible, su convergencia al error en estado estable cero es len-
ta. Para compensar estos efectos, se puede resintonizar el controlador Ks con la
tecnica descrita en [123]. Esta tecnica de sintonizacion permite que el controlador
PD y el fraccionario posean las mismas especificaciones en el dominio del tiempo
(antes tenan las mismas especificaciones en el dominio de la frecuencia).
Las especificaciones de los controladores utilizados en este escenario de simu-
lacion estan contenidas en la Tabla 4.2, el control Ks (resintonizado), se refiere al
control cuyas especificaciones en el dominio del tiempo son las mismas que las del
control PD.
Una vez que ambos controladores tenan una respuesta similar en el dominio
del tiempo, se llevaron a cabo simulaciones adicionales para identificar la robustez
ante cambios repentinos en la inercia de la carga util ademas del par externo no
modelado. La variacion de la inercia considerada durante estas simulaciones fue
94
4.4 Ejemplo de aplicacion
3.5
2.5 3.1
3.0
Posicin angular [rad]
2
2.9
1.5 2.8
2.7
1
2.5
1.0 1.4 1.8 2.0 2.6 3.0
0.5 PD + Observador
Ks + Observador
0 Ks (resintonizado) + Obsv.
Referencia
-0.5
0 1 2 3 4 5 6 7 8 9 10
Tiempo [s]
95
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
3.0 3.0
2.5 2.5
Posicin angular [rad]
3.0
2.5
Posicin angular [rad]
2.0
1.5
1.0
0.5
J=15J025J0 Ks resintonizado + Obsv.
0 PD + Observador
Referencia
-0.5
0 1 2 3 4 5 6 7 8 9 10
Tiempo [s]
Figura 4.13: Respuesta del sistema ante variaciones en la inercia de la carga util (de
JL = 15JL0 a JL = 25JL0 )
96
4.5 Diseno optimo de la relacion de reduccion
coup = kb (m t ). (4.49)
mL2b
Jm n ....
m = 2 t + Jm n + t . (4.50)
b n
97
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
mL2b
m = Jm n + t (4.51)
n
En este caso, para garantizar una senal de control finita, la trayectoria debera ser
como mnimo de 2do orden, lo que garantiza que la aceleracion siempre estara de-
finida. En este sentido, la principal restriccion en el diseno de la trayectoria para
un robot de 1GDL rgido se determina con la maxima aceleracion a la que el brazo
puede someterse antes de que se alcance su maximo esfuerzo admisible. Para el
caso de un robot flexible de 1GDL, se hace necesario tener en cuenta la dinamica
definida por la ecuacion (4.50). Del termino de 4to orden presente en la ecuacion
(4.50) se desprende que para poder obtener una senal de control finita, las trayec-
torias deseadas deberan ser al menos derivables 4 veces, es decir, las trayectorias
deberan ser al menos de 4to orden [127].
Ahora, en primera instancia, la trayectoria del robot debera ser de 4to orden. Sin
embargo, para garantizar que la optimizacion de la reductora es valida para cual-
quier trayectoria, es conveniente considerar una trayectoria que represente el caso
mas exigente posible para el sistema. Esta trayectoria debera permitirle al robot
desplazarse la mayor distancia posible, en el menor tiempo, partiendo del reposo y
culminando en reposo. Con estos aspectos considerados se podra entonces seleccio-
....
nar la aceleracion maxima (t ) y el snap maximo ( t ) a los que se puede someter
el robot, y en base a dichos valores se podra obtener la relacion de reduccion que
minimiza el par suministrado por el motor.
98
4.5 Diseno optimo de la relacion de reduccion
en el sentido de que permite obtener el mnimo par motor para una aceleracion t
determinada con un par motor mnimo. Para el caso de un manipulador rgido, el
proceso de minimizacion del par se puede realizar optimizando la ecuacion (4.50)
en funcion de la relacion de reduccion. De esta forma se tiene:
m
=0 (4.52)
n
cuya solucion optima es: s
mL2b
n = (4.53)
Jm
....
| t (t)| t (t) t (4.54)
mL2b
Jm n
max (t1 ) = max1 = 2 + Jm n + (4.55)
b n 2
mL2b
max (t2 + ) = max2 = Jm n + (4.56)
n
Para poder optimizar la relacion de reduccion para el caso del robot flexible es
necesario determinar cual de los dos pares definidos por las ecuaciones (4.55) y
(4.56) es el par maximo durante la trayectoria.
El par maximo es max1 si se cumple la condicion
mL2b
Jm n
> Jm n + (4.57)
b2 n 2
99
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
t(t)
t(t)
t(t) t2
t1
t
t(t)
t(t)
t
-
Figura 4.14: Trayectorias de 4to orden
100
4.5 Diseno optimo de la relacion de reduccion
= (4.60)
b2
De esta forma, las ecuaciones (4.58) y (4.59) se pueden reescribir de la siguiente
forma:
2
n
2 1 > (4.61)
n
2
n
2 1 < (4.62)
n
Para poder determinar cual de los candidatos a par maximo es en verdad el par
maximo durante la trayectoria se establecen las siguientes dos hipotesis:
mL2b
max1 Jm
= 0 = 2 + Jm 2 (4.63)
n b n 2
101
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
no 1
= (4.65)
n 1 + 2
2
n
= 1 + 2 > 2 1 (4.66)
no
Esto implica que la ecuacion (4.61) deja de ser valida, pero a la vez se verifica
la ecuacion (4.62).
mL2
max2
=0= Jm 2 (4.67)
n n
s
mL2
no = (4.68)
Jm
no
=1 (4.69)
n
no
si < 1 =1 (4.70)
n
102
4.5 Diseno optimo de la relacion de reduccion
Para > 1 se puede establecer que, sino se cumplen las condiciones de las
ecuaciones (4.61) y (4.62), entonces la siguiente ecuacion sera valida 2 1 =
no 2
n
, por lo que si > 1
no 1
= (4.71)
n 2 1
Lo que significa que el criterio que permite seleccionar la relacion de transmision
optima se encuentra en la frontera entre las dos hipotesis.
Las ecuaciones (4.70) y (4.71) permitiran seleccionar la relacion de transmision
optima de acuerdo al valor de .
no
La figura 4.15 muestra el comportamiento de n
a medida que crece. Se ob-
serva que para < 1 la relacion es igual a 1, lo que coincide con la ecuacion (4.70)
y permite obtener la relacion de transmision optima para el caso en el que el brazo
es considerado rgido. Ahora para obtener los valores de par maximo se hara uso
de las ecuaciones (4.71) y (4.55) o (4.56), con lo que se obtiene:
si > 1 omax = 2Jm n (4.72)
2 1
En el caso que < 1 se utilizaran las ecuaciones (4.70), (4.55) y (4.56) para
calcular los valores maximos de par as:
Del par de ecuaciones arriba mostradas se puede decir que el maximo valor
(si < 1) es el dado por la ecuacion (4.74). La Figura 4.16 muestra como el
maximo valor de par aumenta en la medida que aumenta. Este comportamiento
es obvio debido a la dependencia del par maximo con el valor de .
103
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
1.0
0.9
0.8
0.7
n o
n 0.6
0.5
0.4
0.3
0.2
0 1 2 3 4 5 6 7 8 9 10
n0
Figura 4.15: n vs.
2.5
2.0
0max / 2Jn*
1.5
1.0
0 1 2 3 4 5 6 7 8 9 10
0max
Figura 4.16: 2Jn vs.
104
4.5 Diseno optimo de la relacion de reduccion
max Ib
coup = (4.75)
d
Donde max es el esfuerzo flector maximo al que se pretende someter el brazo
flexible y d es la distancia perpendicular al eje neutro de la seccion transversal del
brazo (usualmente igual al radio para una seccion transversal circular). Teniendo en
cuenta que el par de acoplamiento tambien se puede definir como coup = t mL2b ,
se puede reemplazar coup en (4.75) para obtener la siguiente expresion para la
maxima aceleracion en funcion del esfuerzo flector maximo:
max
max = I . (4.76)
mL2b d
El valor maximo del esfuerzo flector, dependera en gran medida del material del
brazo. En el caso de metales, debera ser el valor maximo antes de que el material
experimente deformacion plastica, y en el caso de plasticos y materiales compues-
tos debera limitarse a la zona de comportamiento lineal del material.
105
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
2 2 + = 0 (4.77)
donde: !2
sat
= (4.78)
2Jm n
sat
La solucion real para la ecuacion (4.77) solo existe si 2Jm n
1. En ese
caso, puede ser definido as:
p
= b2 + ( 1) (4.79)
106
4.6 Conclusiones
Criterio
Parametros Rgido Flexible
Maxima Acceleracion () [m/s2 ] 15,3
Snap () [m/s4 ] - 12359,4
Lambda () - 7,5
Mu () - 4
Relacion trasmision optima (no ) 2 0,5
Par maximo (max1 ) [N.m] 0,958 0,452
Par de saturacion (sat ) [N.m] 0,452
Reduccion par maximo 52,8 %
4.6 Conclusiones
En este captulo se han estudiado varios aspectos para mejorar el control de los
actuadores de los manipuladores flexibles.
Por un lado se diseno un esquema de control de movimiento para actuadores
que es robusto a perturbaciones y se valido su desempeno por medio de simulacio-
nes. Este sistema de control esta compuesto por dos elementos clave, un observador
de perturbaciones y un controlador de orden fraccionario del tipo CRONE de 1ra
generacion. Los resultados obtenidos por medio de simulaciones demuestran el
buen desempeno del control propuesto en el caso de grandes cambios de la carga
inercial. La robustez se aumento de forma significativa (en el rango de los valores
simulados) por medio del uso del controlador fraccionario, el cual fue compara-
do con un regulador PD. El control fraccionario se sintonizo en el dominio de la
frecuencia y en el dominio del tiempo, dependiendo del escenario de simulacion,
demostrando la versatilidad del metodo y su robustez superior en comparacion con
107
4. MEJORAS EN EL DISENO Y CONTROL DE ACTUADORES PARA
MANIPULADORES FLEXIBLES
1.0
Par rgido
0.8 Par "flexible"
0.6
0.4
Par motor (N.m) m
0.2
-0.2
-0.4
-0.6
-0.8
-1.0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Tiempo (s)
Figura 4.17: Diferencia entre el par optimo y par para el sistema rgido
108
4.6 Conclusiones
109
CAPITULO
5
Control Robusto del Robot
Flexible de un grado de
libertad
111
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
5.1 Introduccion
Uno de los artculos seminales en el tema de control de manipuladores flexibles fue
el de Robert Cannon y Eric Schmitz [29]. En este artculo se demuestra la efec-
tividad del control de vibraciones en el control de la posicion del extremo de un
robot flexible, en el cual se mide la posicion del extremo del brazo, pero se efectua
la accion de control (par motor) en la base del mismo, demostrando entonces las
dificultades que trae consigo el control de sistemas no colocados, como fue co-
mentado inicialmente por Gevarter en [129], pero a la vez describe las tecnicas de
control utilizadas de forma exitosa, la identificacion del modelo del sistema y los
resultados experimentales.
Posterior al trabajo de Cannon y Schmitz, se han escrito un gran numero de
publicaciones que han apuntado a resolver los problemas que plantea el control de
manipuladores flexible. Debido a esto resulta engorroso citarlos todos. Sin embar-
go, es posible encontrar algunos estados del arte bastante detallados, por ejemplo el
escrito por Benosman y Le Vey [26] o Kiang et al. [6]. Igualmente se pueden encon-
trar libros que contienen estados del arte (aunque no tan detallados), por ejemplo
el de Tokhi y Azad [30] o el de Wang y Gao [130].
Entre las tecnicas de control propuestas por diversos autores para el control de
manipuladores flexibles, es posible encontrar tecnicas robustas ante variaciones en
la carga, siendo esto algo evidente y necesario ya que un robot puede estar someti-
do, de forma intencional o no, a variaciones en la carga util, lo que puede poner en
riesgo su desempeno. Entre las tecnicas mas eficientes para el control de manipu-
ladores flexible ante cargas variables, se pueden considerar las tecnicas de control
adaptativo, control basado en redes neuronales, control deslizante y control robus-
to. Las tecnicas de control adaptativo requieren una estimacion de parametros en
lnea para poder modificar los parametros del controlador en tiempo real. Como
ejemplo se tiene el control no-lineal adaptable usado por Yang et al. en [131] para
el control de un manipulador flexible de dos eslabones en el caso de incertidum-
bre en los parametros del robot. Por otro lado, el problema especfico de cambios
112
5.1 Introduccion
significativos en la carga util se puede ver en [132] y [133]. Del control adaptativo
se pueden identificar las siguientes desventajas: i) la dificultad para garantizar la
estabilidad en lazo cerrado bajo todas las circunstancias; ii) la necesidad de una
excitacion persistente del sistema para poder estimar de forma precisa los parame-
tros del robot y iii) el deterioro de la respuesta transitoria del extremo del robot
en el periodo de tiempo en que se lleva a cabo la estimacion de parametros y se
efectua la accion de control, lo que puede llevar al sistema a presentar errores no
aceptables. Por otro lado, se pueden ver aplicaciones del control robusto a variacio-
nes en la carga de brazos flexibles basados en redes neuronales en [134] y [135].
Estas tecnicas incluyen leyes de control que cambian el comportamiento del con-
trolador en tiempo real, dependiendo de ciertas variables medidas. Al igual que
para las tecnicas de control adaptativo, es difcil garantizar la estabilidad del sis-
tema de control con esta tecnica, sin contar con que se requiere un trabajo previo
(usualmente tedioso) de entrenamiento de la red neuronal. Las tecnicas de control
deslizante tambien se han utilizado para el control ante variaciones de la carga en
manipuladores flexibles, como se puede apreciar en [136] y [137]. Los autores pro-
ponen tecnicas similares, en las cuales la superficie deslizante del controlador se
adapta dependiendo de las condiciones iniciales de funcionamiento de un mani-
pulador flexible. Las tecnicas de control deslizante necesitan satisfacer la llamada
en ingles Matching Condition para poder ser utilizados, y satisfacer esta condicion
requiere que las incertidumbres del sistema permanezcan en el rango espacial de la
senal de entrada, lo que permite a su vez garantizar la propiedad de invarianza del
comportamiento del sistema mientras este en el modo deslizante. Adicionalmente,
el correcto desempeno de la dinamica del brazo no se puede garantizar en la fase
inicial mientras se alcanza la superficie deslizante.
Otro conjunto de tecnicas de control robusto se puede apreciar en [138, 139,
140, 141, 142]. En todas estas metodologas de control robusto, la estabilidad del
sistema controlado se garantiza bajo ciertas condiciones de carga o considerando
efectos despreciables del spill-over, pero la calidad en el seguimiento de trayec-
torias, en el caso de mover la carga util nominal, no es tan buena como la que
se lograra con un controlador especficamente disenado para dicha carga nominal
(asumiendo que este tipo de controlador posee la misma complejidad que el con-
trolador robusto), p. ej. se aumenta el margen de estabilidad a expensas de hacer el
113
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
114
5.2 Generalidades del control Proporcional-Integral Generalizado
115
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
116
5.3 Esquema de Control GPI con Observador de Perturbaciones para el
Control de Vibraciones de un Manipulador Flexible de 1GDL
Lazo externo
Lazo interno
*m u*c Motor DC
-
Externo + +- +- Interno
+ uc Compensacin Flexible
del par de
acoplamiento
kb
-
^
~
1/ 2
b
d Observador de
perturbaciones +
t
1, (o sea mi = m ). De esta forma el lazo externo de control (el del brazo flexi-
ble) se puede disenar sin considerar dinamicas adicionales que puedan complicar
su diseno. Entre las ventajas de este esquema de control se pueden enumerar las
siguientes:
117
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
Si se desea controlar la posicion del extremo del brazo flexible, t , para que siga
una trayectoria determinada por t , con m actuando como senal de control, debe
118
5.3 Esquema de Control GPI con Observador de Perturbaciones para el
Control de Vibraciones de un Manipulador Flexible de 1GDL
existir una senal de control en lazo abierto m que permita hacer el seguimiento de
t , definida como:
1
m = + t
2 t
(5.3)
b
Substrayendo (5.3) de (3.22), surge una expresion en terminos de los errores de
seguimiento de la posicion del motor y de la posicion del extremo del brazo flexible:
et = b2 (em et ) (5.4)
donde em = m m y et = t t . Si se pudiese medir el error en la posicion
del extremo t , se podra proponer el siguiente controlador PID,
Z t
1
em = et + 2 k2t et k1 et k0 et ( )d (5.5)
b 0
Los parametros de diseno del controlador {k2t , k1t , k0t }, se seleccionan de forma
que el polinomio caracterstico del sistema en lazo cerrado se convierta en un poli-
nomio de Hurwitz, cuyas races se encuentran en la posicion deseada. Sin embargo,
para evitar la necesidad de medir la velocidad del error de seguimiento, se puede
utilizar un reconstructor integral para la velocidad del error de la posicion del extre-
mo del robot et . Para obtener dicho reconstructor integral se integra la expresion
(5.4), pero ignorando las condiciones iniciales resultantes despues del proceso de
integracion. De esta forma, la estimacion de velocidad del error [et ]e se puede
calcular de la siguiente forma:
Z t
[et ]e = et (t) et (0) = b2 (em ( ) et ( )) d (5.7)
0
119
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
1 s + 0
(m m ) = (t t ) (5.8)
s + 2
donde
Aunque la posicion del extremo del brazo flexible t no puede ser medida de
forma directa, se puede estimar a partir de la salida del sistema coup , la cual pue-
de medirse mediante el uso de galgas extensiometricas y la salida m , usando la
ecuacion (5.2):
1
t = m coup (5.9)
kb
Para seleccionar los parametros {k2t , k1t , k0t }, es necesario primero seleccionar
en que posicion en el lado izquierdo del plano complejo se desean ubicar los polos
(o races del polinomio) en lazo cerrado. Los tres polos pueden ser ubicados en el
mismo punto del eje real, sea s = a, siendo a estrictamente positivo, en cuyo
caso el polinomio caracterstico del sistema tiene la siguiente forma:
k2t = 3a
k1t = 3a2
k0t = a3 (5.11)
120
5.3 Esquema de Control GPI con Observador de Perturbaciones para el
Control de Vibraciones de un Manipulador Flexible de 1GDL
de esta senal de control proviene del lazo de control del brazo flexible. Dicha refe-
rencia se define como mi . En el captulo 4 se definio la dinamica de un motor que
mova una carga inercial JL y no se consideraron perturbaciones externas. Para este
caso, se considerara que la carga inercial sera igual a la inercia de rotacion del eje
del motor Jm . Adicionalmente se incorporara en la ecuacion el par de acoplamiento
generado por el movimiento del brazo flexible coup . De esta forma se obtiene:
coup
Km V = Jm m + f (m , V ) + (5.12)
n
donde n es la relacion de reduccion de la reductora del motor.
Considerando el modelo dinamico del motor en la ecuacion (5.12), el par de
acoplamiento se puede eliminar anadiendo una compensacion de dicho par, tal co-
mo se ve en la Figura 5.1 (igualmente se describe esto en [147]). En este caso el
voltaje utilizado como senal de control se convierte en
coup
V = Vc + (5.13)
nKm
donde Vc es la senal de control antes de introducir el termino de compensacion. De
esta forma, la dinamica del motor descrita por la ecuacion (5.12) se convierte en:
Km Vc = Jm m + f (m , V ). (5.14)
Jm
Vc = Vc + u+ e
Km Km m
Z t Z tZ 1
u = k3m em k2m em k1m em ( )d k0m em (2 )d2 d1
0 0 0
Z t
K
em = em ( )d e
Jm 0 Jm m
Jm
Vc = m + m (5.15)
Km Km
121
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
122
5.3 Esquema de Control GPI con Observador de Perturbaciones para el
Control de Vibraciones de un Manipulador Flexible de 1GDL
= 0 + (5.22)
donde = b2 y 0 = b0
2
, son el cuadrado de la frecuencia natural actual y la
frecuencia natural nominal, respectivamente. Reemplazando la ecuacion (5.22) en
(5.1), se obtiene:
t = 0 t + 0 m + d (5.23)
donde
d = d + (m t ) . (5.24)
Ya que el angulo del extremo del brazo se puede estimar a partir de la ecuacion
(5.2) y midiendo el par de acoplamiento, coup , y la posicion del motor, m , se
123
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
puede disenar un observador de orden reducido para estimar los estados restantes,
t y d:
2
d z1 f11 1 z1 (f11 + f21 + 0 ) 0
= + t + m (5.27)
dt z2 f21 0 z2 (f11 f21 ) 0
donde z1 y z2 son metaestados que se utilizan para estimar los estados de la si-
guiente forma:
t
1 0 0 t
t = f11 1 0 z1 (5.28)
d f21 0 1 z2
De esta forma, si se desean ubicar los dos polos del observador de forma arbitra-
ria en un lugar , la ecuacion caracterstica sera (s + )2 = 0. Igualando los
terminos de este polinomio y de la expresion (5.29) se obtienen los coeficientes del
observador (siendo > 0):
f11 = 2 (5.30)
f21 = 2 (5.31)
d d
mi = mr = mr 2 . (5.32)
0 b0
2
t (s) b0
(s)
= 2
(5.33)
mr s2 + b0
que corresponde con una dinamica donde el sistema es equivalente a un sistema
cuyo comportamiento es nominal.
124
5.4 Ejemplo de aplicacion
125
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
126
5.4 Ejemplo de aplicacion
Como parte del estudio de la robustez del controlador, se procuro simular los
efectos que pueden ser producidos por dinamicas no modeladas. Las simulaciones
llevadas a cabo consideraron los siguientes efectos: (a) el brazo flexible es contro-
lado para seguir una trayectoria t pero la masa del extremo es diferente a la masa
nominal mt 6= mt0 ; (b) el brazo flexible es controlado para seguir una trayectoria
dada por t con una masa real que es diferente a la masa nominal mt 6= mt0 , y
durante la trayectoria la masa del extremo cambia su valor a un nuevo valor mtf
diferente del nominal, de esta forma mt 6= mt0 , mtf 6= mt0 .
127
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
libre
OBSV. PERT.
(t)
libre
(t)
OBSV. PERT.
libre
(t)
OBSV. PERT.
libre
(t)
OBSV. PERT.
Figura 5.2: Escenario 1: evolucion del seguimiento de trayectoria para la posicion del
extremo simulada, cuando: a) mt = mt0 ; b) mt = 2mt0 ; c) mt = 3mt0 ; d) mt = 4mt0
En los dos casos en que se utiliza el lazo de control GPI se logra estabilizar el
sistema. Sin embargo, el error de seguimiento se logra disminuir de forma signifi-
cativa en el caso en que se incluye el observador de perturbaciones. En la Figura
5.2 es evidente que controlando solo la posicion del motor no es posible anular las
vibraciones debidas a la flexibilidad del brazo.
La Figura 5.3 muestra la senal de control aplicada al sistema bajo las diferentes
condiciones de cambio de la masa real respecto a la nominal. Se puede ver como la
magnitud de la senal de control aumenta cuando el valor real de la masa de la carga
util, mt , aumenta, evidenciando el mayor esfuerzo que tiene que realizar el control
para estabilizar el sistema y seguir la trayectoria de referencia.
En la Figura 5.4 se muestra como la estimacion de la perturbacion cambia de
acuerdo con la diferencia entre el valor real de la masa de la carga util mt y el
valor nominal mt0 , evidenciando un valor mayor de la perturbacion para mayores
variaciones de la masa.
128
5.4 Ejemplo de aplicacion
OBSV. PERT.
OBSV. PERT.
OBSV. PERT.
OBSV. PERT.
Figura 5.3: Escenario 1: senal de control u(t), para el manipulador flexible simulado,
cuando a) mt = mt0 ; b) mt = 2mt0 ; c) mt = 3mt0 ; d) mt = 4mt0
Figura 5.4: Escenario 1: perturbacion estimada d para el manipulador flexible simu-
lado, cuando a) mt = mt0 ; b) mt = 2mt0 ; c) mt = 3mt0 ; d) mt = 4mt0
129
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
1. La masa inicial de la carga util mt era dos veces la masa nominal para la cual
se disenaba el controlador, mt = 2mt0 ,
130
5.4 Ejemplo de aplicacion
Polea
Mesa
Motor DC
de aire
Carga til
Robot
flexible
Figura 5.5: Diagrama esquematico: Dispositivo para cambio de masa durante la tra-
yectoria
131
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
a b
libre
OBS. PERT.
c d
OBS. PERT.
Z
IAE = |e(t)|dt (5.35)
0
Z
IT AE = t|e(t)|dt (5.36)
132
5.4 Ejemplo de aplicacion
GPI
GPI+OBS. PERT.
Jt / (mtL2)
mb / mt
Figura 5.7: Regiones de desempeno satisfactorio para el controlador cuando mb y Jt
cambian, usando el criterio IAE
Las Figuras 5.7 y 5.8 ilustran las regiones en las cuales los controladores (GPI y
GPI + observador) demuestran un desempeno satisfactorio, bajo los criterios IAE
e IT AE. Se puede observar que los valores de la masa de la barra mb y de la inercia
de rotacion de la carga util Jt , pueden tomar valores mayores antes de abandonar
la region de desempeno satisfactorio para el caso del control GPI+Observador de
perturbaciones que para el control GPI estandar. Esto demuestra que el algorit-
mo de control propuesto mejora de forma significativa la robustez del controlador
sin compensacion de perturbaciones, para el caso de dinamicas de alto orden no
modeladas (spill over).
133
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
Jt / (mtL2)
mb / mt
Figura 5.8: Regiones de desempeno satisfactorio para el controlador cuando mb y Jt
cambian, usando el criterio IT AE
La carga util utilizada para los experimentos es un disco plastico con un diame-
tro de 0,09[m] y una masa nominal de mt0 = 0,05[kg]. Para los experimentos
en los que se cambiaba la masa nominal, se utilizaron discos cuyas masas fueron
mt = 4mt0 , mt = 3mt0 y mt = 2mt0 . Estos valores se suponan desconocidos
para el controlador.
134
5.4 Ejemplo de aplicacion
b
GPI+OBS. PERT.
c
GPI+OBS. PERT.
d
GPI+OBS. PERT.
Figura 5.9: Escenario 1: Posicion del extremo del manipulador flexible durante los
experimentos, cuando a) mt = mt0 ; b) mt = 2mt0 ; c) mt = 3mt0 ; d) mt = 4mt0 .
135
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
embargo, el sistema logra alcanzar mas rapidamente el estado estable, y con menor
sobreoscilacion, como se aprecia en la Figura 5.9.
GPI+OBS. PERT
GPI+OBS. PERT
GPI+OBS. PERT
GPI+OBS. PERT
136
5.4 Ejemplo de aplicacion
La Figura 5.12 ilustra la trayectoria seguida por el brazo robotico al usar ambos
controladores (GPI estandar y GPI + Observador de perturbaciones) cuando la ma-
sa del extremo cambia durante el movimiento. En ambos casos, el sistema alcanza
el estado estable sin problemas. Sin embargo, en las Figuras 5.12(a)-(b) se puede
ver que el GPI con el observador de perturbaciones presenta un desempeno mejor
que el GPI estandar ya que se tiene un error de seguimiento menor, se reduce el
error en estado estable y el sistema alcanza el valor de consigna mas rapidamen-
te a pesar de experimentar dinamicas no modeladas. En las Figuras 5.12(c)-(d) se
puede apreciar como la senal de control y la perturbacion estimada experimentan
variaciones bruscas alrededor de t = 1 [s], momento de la trayectoria en el cual se
cambio la masa.
137
5. CONTROL ROBUSTO DEL ROBOT FLEXIBLE DE UN GRADO DE
LIBERTAD
a b
OBS. PERT.
OBS. PERT.
c d
OBS. PERT.
Figura 5.12: Escenario 2: cambio de la masa del extremo durante la trayectoria, para
la plataforma experimental. a) Posicion del extremo; b) Error de seguimiento de la
trayectoria, et = t t , de la posicion del extremo; c) Senal de control, u(t); d)
Estimacion de la perturbacion d.
138
5.5 Conclusiones
5.5 Conclusiones
Del comportamiento del controlador robusto a variaciones en la masa, observado
tanto en simulaciones como en experimentos, puede concluirse que:
139
CAPITULO
6
Control Robusto de un Robot
Flexible de tres grados de
libertad
141
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
142
6.1 Introduccion al Control de Robots Flexibles de multiples grados de
libertad
143
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
R(s)G(s)-1
b) R(s) G(s)-1 G(s) R(s)
Figura 6.1: a) Control en lazo abierto, b) Control en lazo abierto con dinamica inversa
dinamica inversa en lazo abierto la salida del sistema sera identica (teoricamente)
a la senal de referencia.
Otro reto que plantea el considerar la dinamica inversa para el control en lazo
abierto de manipuladores flexibles es que estos manipuladores, ademas de ser siste-
mas de parametros distribuidos, son sistemas de fase no-mnima. Matematicamente
hablando, los sistemas de fase no-mnima son aquellos que tienen al menos un cero
positivo en la funcion de transferencia, lo que implica que al invertir la dinamica
del sistema aparecera un polo positivo, haciendo as que el compensador G1 (s) se
inestable. Diversos autores se percataron de este aspecto y propusieron alternativas
para solucionar este problema. Entre ellos esta Asada y col. [162] [163], De Luca
y Lucibello [164] y de Bayo [165], obteniendo resultados satisfactorios.
Otra alternativa para el control en lazo abierto es el uso de los llamados mol-
deadores de entrada [166] (Input Shapers en ingles). Los Input Shapers son filtros
que toman la senal de referencia y la convolucionan con una serie de impulsos para
generar una nueva entrada que anula frecuencias especficas, determinadas inicial-
mente al disenar el Input Shaper, y que se corresponden con las frecuencias del
sistema que se desean anular.
2 2
1
* 1
0 0
0 1 2 0 0 1 2
Figura 6.2: Funcionamiento del Input Shaper
144
6.1 Introduccion al Control de Robots Flexibles de multiples grados de
libertad
La Figura 6.2 ilustra de forma grafica de que forma se hace uso del Input Sha-
per. A la izquierda de la Figura 6.2 se muestra la respuesta del sistema ante una
entrada escalon (indicada en lnea punteada). Esta entrada escalon se convoluciona
con un par de impulsos para generar una entrada nueva, como se ve en la figura
del lado derecho en lnea punteada. Esta entrada nueva lleva el sistema al punto de
operacion deseado y evitando sobreoscilaciones.
Los metodos de diseno para Input Shapers fueron propuestos inicialmente en
los anos 50 [167, 168], pero debido a la inexistencia de ordenadores digitales su
implementacion resultaba practicamente imposible. Las tecnicas convencionales
de diseno de Input Shapers se basan en el conocimiento de la frecuencia natural
y el coeficiente de amortiguamiento del sistema en cuestion. Sin embargo, cuando
hay incertidumbres en los parametros del sistema o variacion de los mismos, el
Input Shaper disenado puede no reducir las vibraciones como se espera y por ende
arrojar resultados no satisfactorios. Para compensar esta falta de robustez de las
tecnicas convencionales de diseno de Input Shapers han surgido diversas tecnicas
como: diseno de Input Shapers por aprendizaje [169] [170] [171], Input Shapers
adaptativos [172] [173] [174] e Input Shapers robustos [175]. Estas tecnicas de
diseno de Input Shapers mas sofisticados tienen el proposito de disminuir la sensi-
bilidad de la senal de referencia a incertidumbres y no linealidades en el modelo,
a costa de aumentar la duracion del tiempo de subida (directamente proporcional a
la duracion del Input Shapers) [175].
En ultima instancia, se puede decir que el uso de Input Shapers robustos para el
control en lazo abierto es una opcion bastante simple y conveniente, ya que es un
tipo de control de diseno sencillo y con un nivel relativo de robustez, que permite
compensar incertidumbres en la dinamica del sistema y o cambios en los parame-
tros del mismo. Sin embargo, como se comenta en [176] este metodo de control no
es robusto ante perturbaciones externas, lo que implica el uso de controladores en
lazo cerrado para el control de las vibraciones causadas.
145
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
146
6.1 Introduccion al Control de Robots Flexibles de multiples grados de
libertad
Sensor PZT
Brazo Flexible
V Detalle de montaje
de los PZT Actuador PZT
1 2 3 4 t
Figura 6.4: Respuesta Temporal de Sistemas de Fase Mnima y Fase no-Mnima. a=1,
a=0.5; Fase mnima, a=-1; Fase no-mnima, b=c=1
147
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
Para evitar las complicaciones que trae consigo el controlar un sistema de fase
no-mnima, diversos autores han optado por medir o realimentar otras variables (o
estados), como por ejemplo el esfuerzo de la base del brazo [159], la tasa de cambio
de la deformacion del extremo [184], o incluso se realimentan estados virtuales
[44].
En la mayora de las aplicaciones de robotica flexible, los sensores no estan
ubicados en el mismo punto que los actuadores, lo que conlleva a un comporta-
miento de fase no mnima debido a la flexiblidad del sistema [185]. En este sen-
tido, multiples autores han demostrado que el control de las vibraciones del brazo
se puede llevar a cabo mas adecuadamente con una configuracion de actuador y
sensor ubicados en el mismo punto (o muy cerca), convirtiendo al sistema en un
sistema colocado [186], como por ejemplo Park y Asada [187] o Shin y col. [181].
Gracias a los desarrollos recientes en materia de piezo-actuadores (PZT), se han
venido desarrollando mas a fondo tecnicas de control de vibraciones colocadas pa-
ra manipuladores flexibles. Sin embargo, cabe anotar que los robots en los que es
mas facilmente aplicable esta tecnologa de actuadores/sensores son los robots pla-
nares, ya que usualmente sus eslabones son laminas delgadas (ver Figura 6.3), lo
que facilita la instalacion de tales elementos.
Desde el surgimiento del interes en la robotica flexible, multiples tecnicas de
control han sido utilizadas con el proposito de cancelar las vibraciones del mani-
pulador durante su movimiento. Las tecnicas utilizadas abarcan desde el control
clasico como un simple control Proporcional-Derivativo, hasta tecnicas de control
moderno como la linealizacion por realimentacion de estados, control deslizante o
control predictivo basado en modelo.
Debido a la inherente naturaleza de los manipuladores flexibles, resulta tenta-
dor el poner a prueba tecnicas de control avanzadas con plataformas de este tipo.
Por un lado, con el fin de evaluar la validez y eficacia de la tecnica y, por el otro, re-
solver (tal vez de manera mas adecuada) el problema de las vibraciones durante el
movimiento del manipulador flexible. Motivados por el reto que representa el con-
trol de las vibraciones de un manipulador flexible, multiples autores han propuesto
un diverso numero de tecnicas de control (aparte de las tecnicas de control clasi-
co), desde el mismo comienzo de la robotica flexible en el ambito cientfico hace
cerca de 30 anos. Estados del arte en el control de manipuladores flexibles como el
148
6.1 Introduccion al Control de Robots Flexibles de multiples grados de
libertad
149
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
150
6.1 Introduccion al Control de Robots Flexibles de multiples grados de
libertad
151
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
152
6.1 Introduccion al Control de Robots Flexibles de multiples grados de
libertad
153
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
un componente que representa el par computado del robot como si fuera rgido y
finalmente un componente encargado de cancelar las vibraciones del manipulador
flexible, ya que esta asociado a las coordenadas flexibles del sistema. Este ultimo
elemento de la ley de control es el que le otorga al metodo su caracter adaptativo, ya
que esta accion de control posee ganancias variables que varan en funcion de los
parametros del sistema. Yang y col. demuestran la validez del metodo por medio
de un exhaustivo analisis de estabilidad mediante el criterio de Lyapunov, al igual
que con los exitosos resultados experimentales llevados a cabo en un robot flexible
de 2GDL, cuyos eslabones no estan sujetos a los efectos de la gravedad.
El control por espacio de estados es una tecnica tambien utilizada para el con-
trol de manipuladores flexibles, tal es el caso del trabajo de Bai y col. [201]. En esta
publicacion los autores controlan un robot de 2GDL por medio de realimentacion
de estados y la tecnica de sintonizacion LQR. El modelo del sistema se lineali-
za entre dos puntos de operacion, pero dicho modelo es un modelo extendido que
posee dos estados adicionales (no controlables pero observables) que representan
las dinamicas no modeladas (p. ej. la friccion y las dinamicas de alto orden). Para
estimar/observar este vector de estados expandido, los autores proponen el uso de
un filtro llamado STF - Strong Tracking Filter, cuyo funcionamiento es similar al
del filtro de Kalman extendido (EKF). Los resultados experimentales indican erro-
res en estado estable diferentes a cero tanto para el caso de control punto a punto
como para el control de seguimiento. En este aspecto, los autores comentan que
estos errores se deben seguramente a la histeresis que se presenta en la reductora
de los motores. Sin embargo, a grandes rasgos, el metodo propuesto demuestra un
buen comportamiento ante las variaciones de la carga (hasta cerca del 30 % mas
de la carga nominal). Otros autores que hacen uso del EKF son Atashzar y col. ,
en [202], que lo utilizan para estimar perturbaciones en el control de manipulado-
res flexibles. Como tecnica de control usan una linealizacion por realimentacion
de estados que redisenan por medio del criterio de Lyapunov para garantizar su
estabilidad ante perturbaciones acotadas desconocidas. Los autores ponen a prueba
la tecnica propuesta por medio de simulaciones, demostrando resultados satisfac-
torios para la estimacion de perturbaciones (en particular, fuerzas externas). Como
se comentaba anteriormente, el sistema de actuadores y sensores utilizado para
el control del manipulador puede influenciar en gran medida el buen desempeno
154
6.1 Introduccion al Control de Robots Flexibles de multiples grados de
libertad
de una tecnica particular para el control. En torno a estos aspectos, Shin y Choi
[181] proponen lo que ellos llaman un control hbridopara el control de posicion
de un manipulador flexible de 2GDL. Los autores llaman a la tecnica de control
hbrida, en el sentido de que el sistema de control consta de dos conjuntos de
actuadores de diferente naturaleza, dos servo-motores en las articulaciones, y dos
piezo-ceramicos (PZT) adheridos a las superficies de los eslabones flexibles. La
posicion de los servo-motores es controlada por medio de controladores de modos
deslizantes, mientras que el control preciso de la posicion del extremo del bra-
zo y la cancelacion de las vibraciones se lleva a cabo con un control de amplitud
constante por medio de los actuadores piezo-ceramicos. En los resultados experi-
mentales estos autores demuestran que el metodo funciona adecuadamente para el
control punto a punto y para el control de seguimiento. Tambien se han utilizado
recientemente tecnicas modernas de control como el Control Predictivo basado en
Modelo, (en ingles Model Predictive Control, o MPC) para el control de manipu-
ladores flexibles. En [14] Dubay y col. proponen el uso de un controlador MPC en
combinacion con un modelo del manipulador desarrollado por medio de FEM. El
aspecto de modelar el manipulador flexible mediante FEM se considera innovador
en comparacion con otros artculos que han utilizado el MPC como metodo de con-
trol. Las vibraciones se cancelan por medio de tres actuadores PZT distribuidos a lo
largo del brazo flexible y con el control del motor DC que mueve el brazo flexible.
Los autores demuestran por medio de simulaciones y experimentos, que el uso de
un MPC en combinacion con un modelo de FEM brinda mejores resultados en el
control de vibraciones en comparacion con el uso de un MPC estandar (es decir,
en combinacion con un modelo no basado en FEM). Pradhan y Subudhi, desarro-
llan otra aplicacion [15] donde utilizan un control predictivo-adaptativo basado en
un modelo no lineal (NAMPC) para el control de posicion del extremo de un ma-
nipulador flexible de 2GDL planar. Los autores desarrollan el modelo del sistema
usando la formulacion de Lagrange y, por medio de la tecnica NARMAX (en ingles
Non-Linear Autoregressive Moving Average with Exogenous Input), identifican los
parametros del modelo, para luego utilizarlos en el control NAMPC. La tecnica
propuesta demuestra un buen desempeno para el control ante una carga variable,
tanto en simulacion como experimentalmente. El control MPC inteligente, o IMPC
fue tambien utilizado por Tao Fan en [203]. Dicha tecnica combina dos esquemas
155
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
de control, por un lado el control MPC y por el otro tecnicas de logica difusa, lo
que permite sintonizar el controlador e identificar el sistema por medio de metodos
determinsticos y refinar los resultados con tecnicas heursticas.
De forma similar a otros trabajos citados previamente, Halim y col. propo-
nen un control de vibraciones que incluye el uso de actuadores / sensores piezo-
ceramicos (PZT) en los eslabones del robot [13]. El artculo en particular posee
dos aspectos innovadores, por un lado se propone el uso de un modelo no-lineal
del manipulador de multiples GDL basado en el Metodo de Elementos Finitos co-
rotacionales. Segun los autores, este metodo permite tratar de forma mas eficiente
las no-linealidades geometricas que se presentan en el manipulador. El otro aspec-
to a resaltar en el trabajo de Halim y col. es el uso de lo que ellos llaman control
descentralizado. Para implementar este tipo de controlador, se adhieren un par de
PZT a ambos lados de cada uno de los eslabones del manipulador. De esta forma,
uno de los PZT funciona como actuador y el otro como sensor. El control de vi-
braciones de cada eslabon esta desacoplado y no esta ligado a los parametros del
sistema, lo que le brinda cierto nivel de robustez al metodo. El metodo es valida-
do mediante simulaciones con un manipulador flexible de 2GDL, demostrando su
buen desempeno para la cancelacion de vibraciones.
En las tecnicas de control encontradas en el estado del arte se puede evidenciar
como los avances tecnologicos, tanto en el sentido del acceso a ordenadores con
mayor capacidad de calculo as como el desarrollo de nuevos dispositivos de sen-
sado y actuacion (como los piezoelectricos), han posibilitado la evolucion de las
tecnicas de control para estructuras flexibles, no solo desde una perspectiva de la
teora de control y sus avances, sino tambien desde la sensorica y la combinacion de
nuevas tecnologas de actuadores, todo en alas de facilitar y mejorar el rendimiento
de los sistemas de control. Es por esto que se puede afirmar que, aunque el control
de vibraciones en manipuladores flexibles se considera un problema resuelto [26],
siempre estara sujeto a mejoras, ya sea desde el punto de vista teorico de la Inge-
niera de control o desde un punto de vista practico como el de la instrumentacion
industrial.
156
6.2 Esquema de Control Multivariable para un robot de 3 Grados de
Libertad bajo efectos de la gravedad
6.2.1 Cinematica
La posicion del extremo considerando que el robot es rgido se puede describir tanto
T
en coordenadas cartesianas Pr = xR yR zR como en coordenadas esfericas
T
r = R R R , donde ambos sistemas miden la posicion del extremo res-
pecto a la base del robot. Los angulos de las articulaciones (o coordenadas articula-
T
res), se expresan por medio del vector = 1 2 3 . La cinematica directa
157
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
L2 L3 Robot rgido
G3 Robot flexible
L5 F
Z Z
L1 3 L4
nodo 5
2 G2
R F Pr P
G1 P
F
Y R zR Y
R
xR X X
yR
6.2.2 Dinamica
La dinamica utilizada para el diseno del controlador propuesto en este captulo se
deriva directamente de la dinamica obtenida en la subseccion 3.3.2. En particu-
lar, se parte del conjunto de ecuaciones (3.74), (3.75) y (3.76), pero considerando
despreciables los terminos cuadraticos, obteniendose:
F = K P (6.1)
F K11 K12 0
F = K21 K22 0 (6.2)
F 0 0 K33
donde F es el vector de fuerzas (en coordenadas esfericas) que generan las defor-
maciones elasticas del robot, y K es la matriz de rigidez en coordenadas esferi-
cas. Ki,j son sus respectivos elementos como fueron definidos en las ecuaciones
158
6.2 Esquema de Control Multivariable para un robot de 3 Grados de
Libertad bajo efectos de la gravedad
(3.70) en funcion del angulo , pero que con las ecuaciones (3.58) se pueden de-
finir en funcion del angulo de apertura del robot , o en funcion del angulo del
motor 3 (3 ), de forma que, en terminos generales, K ' K (). P es el vec-
tor de deformaciones del extremo en coordenadas esfericas, siendo = R ,
= (R ), = cos()(R ). Las expresiones (6.1)-(6.2) tambien se
pueden expresar en funcion de coordenadas cartesianas como:
h i
T
T()F = T()K ()T() T()P
F = K()P; (6.3)
F = F (6.5)
Esta fuerza F induce pares de reaccion o pares de acoplamiento como los des-
critos por la ecuaciones (3.83) y (3.84). Estos se agruparan en un vector de pares
de acoplamiento:
L1
= L2 (6.6)
L3
que, igualmente, se puede expresar por medio del Jacobiano J() en funcion de
las fuerzas soportadas por el extremo del robot:
= JT ()F (6.7)
159
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
m = QV = m + m + coul + N1 (6.8)
= N1 m (6.9)
La Tabla 6.1 lista los elementos de la ecuacion (6.8). Cabe anotar que las matri-
ces 33 indicadas en la tabla, son matrices diagonales, cuyos valores corresponden
con la descripcion indicada en la misma tabla.
Tamano Descripcion
160
6.2 Esquema de Control Multivariable para un robot de 3 Grados de
Libertad bajo efectos de la gravedad
motores y el lazo de control externo controla la posicion del extremo y cancela las
vibraciones. Ambos lazos de control se explicaran en las siguientes secciones.
Modelo del brazo flexible con el lazo interno
R(s)=(s)=I3x3
+ - + Fg D 2P DP
P* + + Pu u * + P F + P
E(D) fp-1 R(s) (s) fp C-1 1/m
- -
QV0 = m + m (6.11)
161
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
diagonales, y se sintonizan de forma que la dinamica interna del lazo de control sea
mucho mas rapida que la dinamica estructural del manipulador. Dichos valores se
deben hacer tan grandes como sea posible, siempre y cuando no se induzca a los
motores a entrar en sus zonas de saturacion durante el movimiento y se pueda ob-
tener un error de estado estable cercano a cero. Una vez se cierra el lazo de control
expresado en (6.12), la relacion entre m y m se puede expresar como:
Kp m = m + (v + Kv )m + Kp m . (6.13)
donde (s) = N1 m (s) es el vector con los angulos de referencia para las arti-
culaciones del lazo interno. De esta forma se tiene un conjunto de actuadores con
dinamica lineal que no son influenciados por la dinamica de la estructura flexible.
Si la dinamica descrita por (s) es mucho mas rapida que la de la estructura fle-
xible se puede asumir entonces que las senales de referencia para los angulos de
las articulaciones sean usadas como las senales de control en vez de las senales
de voltaje. Este desacoplo se llevo a cabo de la misma forma que el desacoplo
desarrollado para el control del robot de 1GDL mostrado en el captulo 5.
162
6.2 Esquema de Control Multivariable para un robot de 3 Grados de
Libertad bajo efectos de la gravedad
donde = 1 2 3 son las senales entregadas por las galgas y KGM es la
matriz constante que depende de la geometra de los eslabones y de la ubicacion de
las galgas. Conociendo entonces el valor de Fe , se podra estimar la deflexion del
extremo del robot,
Pe = C()Fe (6.16)
P = fP () + P (6.17)
P = fP () + C()F. (6.18)
Pe = fP () + Pe
Pe = fP () + C()Fe . (6.19)
163
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
anotar que la validez de este compensador esta limitada por los niveles de satura-
cion de los actuadores. El esquema de control propuesto para el lazo externo es
entonces:
Pu = P C()Fe + E(D)(P Pe ) (6.21)
donde Ei son las matrices diagonales que representan las ganancias derivativas del
lazo de control externo, y Di es el operador derivada. Una vez se obtiene el vector
para el control de posicion Pu , este se puede convertir en senales de referencia para
los actuadores usando la transformada cinematica inversa en tiempo real
u = fP 1 (Pu ). (6.23)
164
6.3 Analisis de Estabilidad y Robustez
P = Pu + u + P = Pu + u + C()F, (6.24)
C(u )F = P Pu u . (6.25)
E(D)(P P ) = u , (6.26)
165
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
166
6.3 Analisis de Estabilidad y Robustez
= R(D)u (6.27)
Para poder analizar las implicaciones que trae consigo una compensacion im-
perfecta, las ecuaciones (6.27) y (6.14) se pueden expresar como:
Pu = J(f )u . (6.31)
167
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
168
6.3 Analisis de Estabilidad y Robustez
169
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
donde
(1 + i s)2
pi (s) = , i , i , i > 0 (6.44)
(1 + i s)(1 + i s)
siendo i , i , y i el iesimo elemento de la diagonal de las matrices , A, y B
respectivamente, y
donde
e1 s + (e2 mcd,i (f ))s2
gj (s, f ) = (6.46)
1 + mcd,i (f )s2
y cd,i (f ) es el iesimo elemento de la diagonal de Cd (f ). Adicionalmente, ya que
C (f ) es una matriz simetrica de la forma
c,11 (f )
0 0
C (f ) = 0 c,22 (f ) c,23 (f ) , (6.47)
0 c,23 (f ) c,33 (f )
I + H(s)G(s) = 0, (6.49)
donde
h1,1 (s) 0 0
H= 0 h2,2 (s) h2,3 (s) (6.50)
0 h3,2 (s) h3,3 (s)
170
6.3 Analisis de Estabilidad y Robustez
siendo
Esto implica tener que encontrar las races de un polinomio de 12do (el determinante
de la matriz indicada en (6.49)), lo que puede resultar numericamente complicado.
Por otro lado, esta expresion no permite obtener condiciones analticas simples para
conservar la estabilidad del sistema en lazo cerrado (condiciones de robustez).
171
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
y por (p2 (s) y p3 (s)). Siendo as, la condicion (6.53) debe ser verificada para
todo el rango de angulos de apertura y para todo el rango de valores esperados
para las matrices A y B (incluyendo posibles valores complejos).
La condicion establecida por la ecuacion (6.53) puede verificarse numerica-
mente para el rango de configuraciones del robot f [, ] y para el rango de
variacion de los parametros de los motores. Los parametros de los motores que son
mas propensos a cambiar son la friccion dinamica y la inercia de rotacion, mientras
que los cambios en la ganancia del motor son poco probables. De esta forma, se
asume un rango de variacion dado por el maximo porcentaje de los valores nomi-
nales de la friccion dinamica del motor y la inercia de rotacion. En la seccion de
simulaciones se mostrara que la condicion de diagonal estrictamente dominante se
satisface en el robot flexible de 3GDL para todos los valores del angulo de apertura
y para el rango dado de variaciones en los parametros.
2
1 c,i e2 > 0; (6.54)
(ii) que la funcion de transferencia hi,i (s) es una red de adelanto de fase en todo el
rango de frecuencia, de forma que
172
6.3 Analisis de Estabilidad y Robustez
I =c,i
IV
-1 Eje
real
=c,i III
VI
II
173
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
nerar un encerramiento no deseado del punto (1, 0). De esta forma, la necesidad
de la condicion 2 queda demostrada. Finalmente, se puede observar igualmente en
la Figura 6.7 que, para poder garantizar que la ultima parte del diagrama de Ny-
quist (cuando ) permanezca al lado derecho del punto (1, 0), la condicion
2
lm hi,i (j)(c,i e2 1) > 1 debe ser verificada, lo que conlleva a la condi-
cion 3.
2
(1+i s)
Lema 2. Considere pi (s) = (1+ i s)(1+i s)
, entonces una condicion suficiente que
garantiza la estabilidad del sistema resultante al cerrar el lazo de control alrededor
de gi (s)pi (s) con una realimentacion negativa unitaria es:
s !
i + i i i
< i < mn 2(i + i ), 2
, (6.59)
2 1 e2 c,i
i i i , i i i , 1 i 3 (6.60)
, (6.61)
m m m, (6.62)
174
6.3 Analisis de Estabilidad y Robustez
e2 m mn cd,i (6.63)
1i3
v i + Kv,i v + Kv,i
< i < 2 i ,1 i 3 (6.64)
2Kp,i Kp,i
s
1 m cd1
1 < (6.65)
Kp1 (m cd1 e2 )
Kp,2 2 Kp,3 mcd,2 ()
cos (())32 + sin2 (())32 < (6.66)
2 3 mcd,2 () e2
Kp,2 Kp,3 2 mcd,3 ()
sin2 (())22 + cos (())32 < (6.67)
2 3 mcd,3 () e2
donde las condiciones de (6.66) y (6.67) deben satisfacerse para todos los valo-
res [, ] ya que y cd,i dependen de . Las ganancias Kp,i y Kv,i son los
parametros del controlador P D del iesimo motor.
(1 + i s)2
pj (s) = i 2
, (6.68)
1 + vi K
+Kvi
pi
s + Kp,i
s
175
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
6.4 Simulaciones
Dos conjuntos de estudios se llevaron a cabo:
176
6.4 Simulaciones
treo para las simulaciones de seguimiento de trayectoria fue de 1 ms. Las Tablas 6.2
y 6.3 contienen los parametros del robot considerados para todas las simulaciones.
Tabla 6.2: Rangos para los angulos de las articulaciones y de apertura, y longitudes
de eslabones
Como se comento previamente, los parametros del controlador del lazo de con-
trol interno, Kp y Kv , se han seleccionado de forma que se pueda alcanzar un lazo
de control lo mas rapido posible, pero sin que esto implique la saturacion de los
actuadores (2,3 V para todos los motores). La funcion de transferencia equiva-
1
lente para los lazos de control esta dada por i (s) = (1+ai s)2
(1 i 3), donde el
177
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
valor para ai se ha seleccionado como 1/30 para todos los motores. Los parametros
obtenidos para cada uno de los controladores se agrupan en la Tabla 6.4. El valor
nominal de la carga util es m0 = 1,975kg. La variacion de la carga util va desde su
valor mnimo, m = 0,925 kg, hasta su valor maximo, m = 2,975 kg.
Motor Kp Kv
1 154.175 9.1082
2 88.616 5.1849
3 103.478 6.2047
Tabla 6.4: Rangos para los angulos de las articulaciones y de apertura, y longitudes
de eslabones
178
6.4 Simulaciones
2 3
0.02 0.02
Eje imaginario
Eje imaginario
0.01 0.01
0 0
0.01 =0.008627 0.01 3=0.009038
2
0.02 0.02
0.03 0.03
0 0.01 0.02 0.03 0.04 0.05 0 0.01 0.02 0.03 0.04 0.05
Eje real Eje real
2 3
0.02 0.02
Eje imaginario
179
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
0.12
0.1 (,)
0.08
0.06 Plano = 0
0.04
0.02
0 0
100
80
0.005
60
(grados) 40
20 0.01
En primer lugar, cabe anotar que la condicion (6.63) establecio el valor maximo
permitido para e2 . En esta simulacion, los valores de la matriz de flexibilidad y su
matriz diagonalizada, C() y Cd () se han obtenido para el rango de variacion
del angulo de apertura . El valor mnimo obtenido de cd,i , cd,i , se incluye en la
expresion (6.63), de forma que:
180
6.4 Simulaciones
e =1.8088e005
2
0.2
0.18
0.16
0.14
0.12
3
0.1 2 y 3 seleccionados
0.08
0.06
zona de
0.04 estabilidad
0.02
0
0 0.05 0.1 0.15 0.2
2
181
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
182
6.4 Simulaciones
0.6 0.6
0.55 0.55
z
z
0.5 0.5
0.45 0.45
0.4 0.4
0.3 0.2 0.3 0.2
0.4 0.1 0.4 0.1
0.5 0.5
0.6 0 0.6 0
0.1 y 0.1 y
x 5 x
x 10
0.01 1
seguimiento et (m)
seguimiento et (m)
0.008 0.8
0.006 0.6
Error de
Error de
0.004 0.4
0.002 0.2
0 0
0 1 2 3 4 5 6 0 1 2 3 4 5 6
t(s) t(s)
lazo interno que controla los motores siempre esta cerrado), y considerando que la
carga util vara de m a m. Esta figura tambien representa el error de seguimiento,
et , definido como:
p
et = (x x)2 + (y y)2 + (z z)2 , (6.79)
muestran que el error de seguimiento crece a medida que la carga util crece. Cuan-
do el lazo externo de control esta abierto, el error de seguimiento esta en el rango
et [103 , 102 ], y cuando este se cierra, es alrededor de mil veces menor.
La Figura 6.12 representa las senales de control para el sistema de control pro-
puesto. Estas senales se corresponden con los voltajes aplicados a los motores para
llevar a cabo el seguimiento de trayectoria con las diferentes cargas.
183
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
2 m=0.925 Kg.
Seal control
motor1 (V)
m=1.975 Kg.
0 m=2.975 Kg.
2
0 1 2 3 4 5 6
t(s)
2
Seal control
2
0 1 2 3 4 5 6
t(s)
2
Seal control
motor3 (V)
2
0 1 2 3 4 5 6
t(s)
184
6.4 Simulaciones
3 3 3
x 10 x 10 x 10
2.5
2 1.5 1.5
(t)
(t)
(t)
1.5
1
3
1 1
1
0.5 0.5
0.5
0 2 4 6 0 2 4 6 0 2 4
t(s) t(s) t(s)
4 4 4
x 10 x 10 x 10
1.6 2
2.3
1.9
1.5
2.2
1.8
(t)
(t)
(t)
1.4
1
3
2
2.1
1.7
2 1.3
1.6
1.9
1.2 1.5
0 2 4 6 0 2 4 6 0 2 4 6
t(s) t(s) t(s)
185
6. CONTROL ROBUSTO DE UN ROBOT FLEXIBLE DE TRES GRADOS
DE LIBERTAD
5
x 10
3
referencia
salida
0.7 2.5
0.65
e (m)
2
0.6
seguimiento t
0.55 1.5
Error de
z
0.5
1
0.45
0.4 0.5
0.3 0.2
0.4 0.1
0.5
0.6 0 0
0.7 0.1 0 1 2 3 4 5 6
y
x t(s)
0.015
Referencia con lazo de control externo
e (m)
0.65 0.005
0.6
0
0 1 2 3 4 5 6
0.55
z
3
0.45
seguimiento t
2
0.4
Error de
0.3 0.2 1
0.4 0.1
0.5 0
0.6 0
0.7 0.1
y 0 1 2 3 4 5 6
x t(s)
186
6.5 Conclusiones
6.5 Conclusiones
Se ha disenado una algoritmo de control novedoso que permite controlar de forma
satisfactoria un robot flexible de 3GDL. Debido a las simplificaciones en el mo-
delo dinamico y a las consideraciones de masa concentrada, solo tres modos de
vibracion estan presentes en el robot.
La tecnica de control propuesta es bastante simple y es capaz de controlar un
sistema de orden superior, como lo es el robot, el cual tiene 12 estados. El algoritmo
propuesto esta compuesto por dos lazos de control, el lazo interno que controla la
posicion de los actuadores, y el lazo externo que controla la posicion del extremo
y cancela las vibraciones.
Se desarrollaron analisis de su estabilidad y robustez, proporcionando las con-
diciones de estabilidad para el caso en que hay incertidumbre en los parametros del
motor. En la seccion de simulaciones se validaron las condiciones de estabilidad
para el robot y se demostro que el algoritmo de control propuesto funciona de for-
ma satisfactoria para el seguimiento de trayectorias, garantizando la estabilidad y
la no saturacion de los actuadores. Por otra parte, los cambios en la carga nominal
del robot no fueron obstaculo para que el control mantuviera un buen desempeno
en el control de trayectorias, garantizando la estabilidad y el error en estado estable
cercano a cero.
187
CAPITULO
7
Conclusiones, Aportaciones y
Trabajos Futuros
7.1 Conclusiones
Se ha desarrollado un metodo de control de movimiento de actuadores (con posi-
bles usos a manipuladores flexibles), basado en reguladores fraccionarios y obser-
vadores de perturbaciones. Este metodo de control es robusto ante perturbaciones
externas como la friccion o pares externos no modelados, as como a grandes va-
riaciones de la carga inercial que mueve el motor. El regulador se puede sintonizar
tanto en el dominio del tiempo como en el dominio de la frecuencia, siendo el
primer metodo de sintonizacion mas conveniente para poder garantizar especifica-
ciones temporales de la respuesta del sistema.
Se ha establecido un criterio de seleccion de la relacion de reduccion optima
para manipuladores flexibles de 1GDL. La relacion de reduccion obtenida permite
obtener senales de control de menor intensidad, en comparacion con la relacion
189
7. CONCLUSIONES, APORTACIONES Y TRABAJOS FUTUROS
190
7.2 Aportaciones
7.2 Aportaciones
A continuacion se enumeran las aportaciones cientficas de esta tesis.
7.2.1 Congresos
Partida, G., Jaramillo, V. H., and Feliu, V. (2010, Noviembre). Mass varia-
tions robust controller, for a 3-DOF flexible robot. In IECON 2010-36th An-
nual Conference on IEEE Industrial Electronics Society (pp. 1558-1564).
IEEE.
Castillo, F. J., Feliu, V., Sanchez L., Rivas, R., Jaramillo, V. H., (2010, Oc-
tubre) Robust PD Time-Domain Tuning Rules for Controlling DC-Motors.
4th IFAC Workshop on Fractional Differentiation and Its Applications. IFAC.
191
7. CONCLUSIONES, APORTACIONES Y TRABAJOS FUTUROS
Gomez L.M., Jaramillo, V. H., Castillo, F. J., and Feliu, V. (2009, Noviem-
bre). Time Domain Tuning of a Class of a Smith Predictor Based Fractional
Control. In Fractional Signals and Systems, 2009 Symposium on. M. Orti-
gueira et al. (eds.).
7.2.2 Revistas
Feliu, V., Castillo, F., Jaramillo, V., and Partida, G. (2013). A Robust Contro-
ller for A 3-DOF Flexible Robot with a Time Variant Payload. Asian Journal
of Control, 15(4), 971-987.
Morales, R., Feliu, V., and Jaramillo, V. (2012). Position control of very light-
weight single-link flexible arms with large payload variations by using dis-
turbance observers. Robotics and Autonomous Systems, 60(4), 532-547.
192
7.3 Trabajos futuros
193
APENDICE
8
Apendices
Propiedad Valor
Longitud [m] 0,7
2
Area [m ] 1 108
Momento de Inercia del area [m4 ] 3,017 1012
Modulo de Young (E) [P a] 122,7 109
195
8. APENDICES
Propiedad Valor
Diametro Interior 13[mm]
Diametro exterior 16[mm]
Momento de inercia del area 1,815 109 [m4 ]
Area seccion transversal 6,8330 10 5[m2 ]
Modulo de elasticidad 6,6847 1010 [N/m2 ]
Modulo de elasticidad a torsion 4,9811 1010 [N/m2 ]
Carga util 1
Masa 1,2[kg]
Diametro [mm]
Espesor [mm]
Momentos de Inercia Ixx e Iyy [kg.m2 ]
Momento de inercia Izz [kg.m2 ]
Carga util 2
Masa 1,6[kg]
Diametro [mm]
Espesor [mm]
Momentos de Inercia Ixx e Iyy [kg.m2 ]
Momento de inercia Izz [kg.m2 ]
196
8.1 Especificaciones de los componentes mecanicos
Parametro Valor
Reductora motor 1 (RFS-25-100-CC)
Relacion de reduccion 100 : 1
Par nominal 67[N.m]
Par de pico 157[N.m]
Inercia vista desde el rotor 1,71 104 [kg.m2 ]
Rigidez torsional 2,5 104 [N.m/rad]
Reductoras motor 2 y 3 (RFS-25-160-CC)
Relacion de reduccion 160 : 1
Par nominal 98[N.m]
Par pico 157[N.m]
Inercia vista desde el rotor 1,71 104 [kg.m2 ]
Rigidez torsional 2,5 104 [N.m/rad]
197
8. APENDICES
Caractersticas de la barra 1
Nodos que une 1 con 2
Longitud 0,625[m]
Masa 0,098[kg]
Caractersticas de la barra 2
Articulaciones pasivas
Nodos que une 2 con 3
Nodo 2
Longitud 0,7004[m]
Masa 0,1162[kg] Eslabones que une 1 con 2
Caractersticas de la barra 3 Masa de la articulacion 0,1078[kg]
Nodos que une 3 con 4 Nodo 3
Longitud 0,3125[m]
Eslabones que une 2 con 3
Masa 0,0465[kg]
Masa de la articulacion 0,1078[kg]
Caractersticas de la barra 4
Nodo 4
Nodos que une 4 con 5
Eslabones que une 3 y 4 con 5
Longitud 0,625[m]
Masa de la articulacion 0,3020[kg]
Masa 0,1023[kg]
Caractersticas de la barra 5 Tabla 8.6: Caractersticas de las ar-
ticulaciones pasivas
Nodos que une 1 con 4
Longitud 0,625[m]
Masa 0,0903[kg]
198
8.2 Especificaciones de los componentes electricos
Parametro Valor
Potencial nominal [W ] 8,6
Voltaje nominal [V ] 24
Corriente nominal [A] 1
Corriente de pico [A] 1,6
Par nominal [N.m] 1,4
Par pico [N.m] 2,7
Velocidad nominal [rpm] 60
Constante de par [N.m/A] 2,1
Constante de voltaje [V /rpm] 0,22
Relacion de reduccion 50
Resistencia del estator [] 10
Inductancia del estator [mH] 2,2
Parametro Valor
Inercia del eje del motor (Jm ) 6,873 105 [kg.m2 ]
Friccion viscosa () 1,041 103 [N.m.s]
Constante de par del motor (Km ) 0,21[N.m/V ]
Relacion de reduccion (nr ) 50
199
8. APENDICES
Parametro Valor
Voltaje de alimentacion 12 30[V DC]
Voltaje maximo de salida 25[V ]
Corriente maxima de salida 2[A]
Potencia maxima de salida 50[W ]
Especificaciones de las entradas
Senal de control Configurable +/- 10 o +/- 3.9 [V]
Entrada para tacometro DC 2 50[V DC]
Entrada para encoder Frequencia maxima 100kHz (tipo TTL)
Especificaciones de las salidas
Salida al motor 25[V DC]
Voltaje de alimentacion del encoder 5[V DC]
Parametro Valor
Voltaje de alimentacion 5 [V DC]
Tipo de senal TTL
Tipo de encoder Incremental
Resolucion 500[pulsos/rev]
Frecuencia maxima 100[kHz]
Parametro Valor
Voltaje de alimentacion 6 16[V ]
Corriente de alimentacion 5[mA]
Voltaje de alimentacion al encoder 5[V ]
Maxima frequencia de entrada 1,2[M Hz]
Impedancia de entrada 10[k]
200
8.2 Especificaciones de los componentes electricos
Especificaciones
Familia NI 6024E
Bus PCI, PCMCIA
Entradas analogicas
Numero de entradas 16
Resolucion 12 bits
Tiempo de muestreo maximo 200[kS/s]
Rango de entrada 10[V ]
Salidas analogicas
Numero de salidas 2
Resolucion de salida 12 bits
Tiempo de muestreo a la salida 10kS/s
Especificaciones
Numero de canales 24 I/O
Salida nivel alto 3,0[V ]@ 2,5[mA]
Salida nivel bajo 0,4[V ]@2,5[mA]
Entrada nivel alto 2,8[V ] min, 5,5[V ] max. abs.
Entrada nivel bajo 0,8[V ] max, 0,5[V ] min. abs.
Consumo 170[mA] nominal, 270[mA] max. @5[V ]
201
8. APENDICES
Parametro Valor
Velocidad Nominal 6000[rpm]
Par nominal 0,77[N.m]
Par de pico 3,10[N.m]
Constante de par 0,38[N.m/A]
Constante de velocidad 0,22[V /rad/s]
Resistencia por fase 8,4[]
Inductancia por fase 12,8[mH]
Inercia del rotor 0,29[kg.cm2 ]
Parametro Valor
Voltaje de alimentacion 200[V DC]
Voltaje por fase 130[V AC]
Corriente nominal por fase 2,5[A]
Corriente pico 5[A]
202
8.2 Especificaciones de los componentes electricos
Parametro Valor
Voltaje de entrada 2[mV /V ]
Resistencia del puente 87 5000[]
Voltaje de excitacion 5[V DC]
Ancho de banda 250[Hz]
Voltaje de salida 5[V DC]
Voltaje de alimentacion 10,5 18[V DC]
Consumo de corriente 30 + corriente del puente [mA]
Temperatura de funcionamiento 10 + 65[o C]
Especificaciones
Familia NI-6229
Bus PCIe
Entradas analogicas
Numero de entradas 32
Resolucion 16 bits
Tiempo de muestreo maximo 250[kS/s]
Rango de entrada 10[V DC]
Salidas analogicas
Numero de salidas 4
Resolucion de salida 16 bits
Tiempo de muestreo a la salida 625kS/s (por salida)
Rango de salida 10[V DC]
203
8. APENDICES
Especificaciones
Familia NI-6602
Bus PCI
Contadores/Temporizadores
Numero de contadores 8
Resolucion 32 bits
Tiempo de muestreo maximo 80[M Hz]
Compatibilidad 5[V ] TTL/CMOS
Especificaciones
Maxima frecuencia de muestreo 240 [Hz]
Rango de medida 0,2 70 [m]
Maximo numero de marcas 150 a 60 [Hz]
Numero de Leds 250
Rayo de luz Infrarrojo a 880 [nm]
Resolucion del sensor de imagen 658 500 pixeles
Velocidad de transferencia MCU al computador 2 [M bps]
204
8.3 Generacion de trayectorias
% ENTRADAS
% p: Orden del polinomio
% t0: Tiempo inicial de la trayectoria
% tf: Tiempo final de la trayectoria
% d: Numero de pasos (steps)
% SALIDAS
% theta: Vector de posicion
% t: Vector de tiempo
delt=(tf-t0)/d;
205
8. APENDICES
for i=1:p+1
for j=1:p+1
A(i,j)=(-1)(j+1)*prod(1:p+j-1)/prod(1:p+j-i)*(1/(tf-t0)(i-1));
end
end
P=inv(A);
for r=1:p+1
for k=1:d+1
t1=t0+(k-1)*delt;
sum(r)=P(1,1)*prod(1:p)/prod(1:p-r+1);
for l=2:p+1
fac(r)=prod(1:p+l-1)/prod(1:p+l-r);
psi(l,k,r)=(-1)(l+1)*P(l,1)*(abs(t1-...
t0)/(tf-t0))(l-1)*fac(r);
sum(r)=sum(r)+psi(l,k,r);
end
bez(r,k)=1-((abs(t1-t0)/(tf-t0))(p+...
1-r))/(tf-t0)(r-1)*sum(r);
end
end
t=[t0:delt:tf];
phi=bez(1,:);
206
Referencias
[1] Canadian Space Agency. Canadarm2, the Canadian Robotic Arm of the
International Space Station. [Online]. Disponible en: http://www.
asc-csa.gc.ca/eng/iss/canadarm2/, 2015. [Visitado el 24 de
noviembre 2015].
[3] Canadian Space Agency. Dextre, the International Space Stations Robotic
Handyman. [Online]. Disponible en: http://www.asc-csa.gc.ca/
eng/iss/dextre/default.asp, 2015. [Visitado el 01 Julio 2015].
[5] H.N. Rahimi and M. Nazemizadeh. Dynamic analysis and intelligent con-
trol techniques for flexible manipulators: a review. Advanced Robotics,
28(2):6376, 2014.
[6] ChangTai Kiang, Andrew Spowage, and ChanKuan Yoong. Review of con-
trol and sensor system of flexible manipulator. Journal of Intelligent and
Robotic Systems, 77(1):187213, 2015.
207
REFERENCIAS
[8] B. Siciliano and L. Villani. Parallel force and position control of flexi-
ble manipulators. Control Theory and Applications, IEE Proceedings -,
147(6):605612, Nov 2000.
[9] Zhao-Hui Jiang. Impedance control of flexible robot arms with parametric
uncertainties. Journal of Intelligent and Robotic Systems, 42(2):113133,
2005.
[10] Alin Albu-Schaffer, Christian Ott, and Gerd Hirzinger. A unified passivity-
based control framework for position, torque and impedance control of fle-
xible joint robots. International Journal of Robotics Research, 26(1):2339,
2007.
[12] Said G. Khan, Guido Herrmann, Mubarak Al Grafi, Tony Pipe, and Chris
Melhuish. Compliance control and humanrobot interaction: Part 1 - survey.
International Journal of Humanoid Robotics, 11(3), 2014.
[13] Dunant Halim, Xi Luo, and Pavel M. Trivailo. Decentralized vibration con-
trol of a multi-link flexible robotic manipulator using smart piezoelectric
transducers. Acta Astronautica, 104(1):186 196, 2014.
[14] Rickey Dubay, Marwan Hassan, Chunying Li, and Meaghan Charest. Finite
element based model predictive control for active vibration suppression of a
one-link flexible manipulator. ISA Transactions, 53(5):16091619, 2014.
[15] Santanu Kumar Pradhan and Bidyadhar Subudhi. Nonlinear adaptive model
predictive controller for a flexible manipulator: An experimental study. IEEE
Transactions on Control Systems Technology, 22(5):1754 1768, 2014.
208
REFERENCIAS
[16] M.S. de Queiroz, D.M. Dawson, M. Agarwal, and F. Zhang. Adaptive nonli-
near boundary control of a flexible link robot arm. Robotics and Automation,
IEEE Transactions on, 15(4):779787, Aug 1999.
[19] Matteo Laffranchi, Lisha Chen, Navvab Kashiri, Jinoh Lee, Nikos G. Tsaga-
rakis, and Darwin G. Caldwell. Development and control of a series elastic
actuator equipped with a semi active friction damper for human friendly ro-
bots. Robotics and Autonomous Systems, 62(12):1827 1836, 2014.
[23] Yuri A Ivanov and Rodney Brooks. Robotic placement and manipulation
with enhanced accuracy, October 10 2014. US Patent App. 14/511,718.
[24] Bruce Blumberg, Rodney Brooks, Christopher J Buehler, Noelle Dye, Gerry
Ens, Natan Linder, Michael Siracusa, Michael Sussman, and Matthew M
Williamson. Training and operating industrial robots, March 31 2015. US
Patent 8,996,175.
209
REFERENCIAS
[25] Rethink Robotics. BaxterTM Safety and Compliance Overview - How this
unique collaborative robot safely manages operational risks. [Online].
Disponible en: http://www.active8robots.com/wp-content/
uploads/Baxter-Safety-and-Compliance-Overview.pdf,
2013. [Visitado el 24 Junio 2015].
[28] Usoro P. B., Nadira R., and Mahil S. S. A finite element/lagrange approach
to modeling lightweight flexible manipulators. Journal of Dynamic Systems,
Measurement, and Control, 108(3):198205, 1986.
[29] Robert H. Cannon and Eric Schmitz. Initial experiments on the end-point
control of a flexible one-link robot. The International Journal of Robotics
Research, 3(3):6275, 1984.
[30] M.O. Tokhi and A.K.M. Azad. Flexible Robot Manipulators - Modeling,
Simulation and Control. The institution of Engineering and Technology,
2008.
[33] G.H. Handelman, W.E. Boyce, and H. Cohen. Vibrations of a uniform rota-
ting beam with tip mass. In Proceedings of the Third U.S. National Congress
of Applied Mechanics, pages 175180, Providence, Rhode Island, 1958.
210
REFERENCIAS
[34] W.E. Boyce and G.H. Handelman. Vibrations of rotating beam with a tip
mass. Zeitschrift fur angewandte Mathematik un Physik, 12:369392, 1961.
[35] R. Kumar. Vibrations of space booms under centrifugal force field. Tran-
sactions of the C.A.S.I, 7:15, 1974.
[36] S. Putter and H. Manor. Natural frequencies of radial rotating beams. Jour-
nal of Sound and Vibration, 56:175185, 1978.
[37] Book W. Study of design and control of remote manipulators part ii: Vi-
bration considerations in manipulator design. Journal of Dynamic Systems
Measurement and Control, February 1973.
[38] Whitney D., Book W., and Lynch P.M. Design and control considerations
for industrial and space manipulators. Technical report, Department of Me-
chanical Engineering, MIT, 1974.
[39] H. Kanoh, S. Tzafestas, H.G. Lee, and J. Kalat. Modelling and control of
flexible robot arms. In Proceedings of the 25th Conference on Decision and
Control, pages 18661870, Athens, Greece, 1986.
[40] Santosha Kumar Dwivedy and Peter Eberhard. Dynamic analysis of fle-
xible manipulators, a literature review. Mechanism and Machine Theory,
41(7):749 777, 2006.
[41] Moheimani S.O.R. and Clark R.L. Minimizing the truncation error in assu-
med modes models of structures. In Proceedings of the American Control
Conference, volume 4, page 23982402, 2000.
[44] M. Vidyasagar D. Wang. Transfer function for a single flexible link. The
International Journal of Robotics Research, 10(5):540549, 1991.
211
REFERENCIAS
[45] F. Belleza, L. Lanari, and G. Ulivi. Exact modeling of the flexible slewing
link. In Proceedings of the IEEE International Conference on Robotics and
Automation, pages 734739, 1990.
[46] E. Bayo. Timoshenko versus bernoulli beam theories for the control of fle-
xible robots. In Proceeding of IASTED International Symposium on Applied
Control and Identification, pages 178182, 1986.
[48] H. Kanoh and H.G. Lee. Vibration control of one-link flexible arm. In
Proceedings of the 24th IEEE Conference on Decision and Control, pages
11721177, Fort Lauderdale, FL, USA, 1985.
[49] W.J. Book and M. Majette. Controller design for flexible, distributed pa-
rameter mechanical arms via combined state space and frequency domain
techniques. ASME Journal of Dynamic Systems, Measurement, and Control,
105:245254, 1983.
[53] AAM Al-khafaji and IZM Darus. Finite element method to dynamic mode-
lling of an underwater flexible single-link manipulator. Journal of Vibroen-
gineering, 16(7):36203636, 2014.
212
REFERENCIAS
[54] Z. Mohamed and M.O. Tokhi. Command shaping techniques for vibration
control of a flexible robot manipulator. Mechatronics, 14(1):6990, 2004.
[56] Moulin H. and Bayo E. On the accuracy of end-point trajectory tracking for
flexible arms by noncausal inverse dynamic solutions. Journal of Dynamic
Systems Measurement and Control, 113(2):320324, 1991.
[57] J.D. Lee and B.L. Wang. Optimal control of a flexible robot arm. Computers
and Structures, 29(3):459467, 1988.
[58] M. Dado and A. Soni. A generalized approach for forward and inverse dy-
namics of elastic manipulators. In Robotics and Automation. Proceedings.
1986 IEEE International Conference on, volume 3, pages 359364, Apr
1986.
[59] Feliu V., Rattan K.S., and Brown Jr. H.B. Modeling and control of single-
link flexible arms with lumped masses. Journal of Dynamic Systems Mea-
surement and Control, 114(1):5969, 1992.
[62] Ge S.S., Lee T.H., and Zhu G. A new lumping method of a flexible manipu-
lator. In Proceedings of the American Control Conference, volume 3, page
14121416, 1997.
[63] Zhu G. Ge S.S. and Lee T.H. Simulation studies of tip tracking control of a
single-link flexible robot based on a lumped model. Robotica, 17(1):7178,
1999.
213
REFERENCIAS
[65] J.C. Simo and L. Vu-Quoc. On the dynamics of flexible beams under large
overall motions the plane case: Part i. Journal of Applied Mechanics,
53(4):849854, 1986.
[66] J.C. Simo and L. Vu-Quoc. On the dynamics of flexible beams under large
overall motions the plane case: Part ii. Journal of Applied Mechanics,
53(4):855863, 1986.
[69] W. Book. Analysis of massless elastic chains with servo controlled joints.
ASME Journal of Dynamic Systems Measurement and Control, 101, Sep-
tember 1979.
214
REFERENCIAS
[75] E. Bayo. Timoshenko versus bernoulli-euler beam theories for inverse dyna-
mics of flexible robots. International Journal of Robotics and Automation,
4(1):5356, 1989.
[77] J.A. Somilinos, V. Feliu, and L. Sanchez. Design, dynamic modeling and
experimental validation of a new three-degree-of-freedom flexible arm. Me-
chatronics, 12(7):919948, 2002.
[78] R.J. Theodore and A. Ghosal. Comparison of the assumed modes and finite
element models for flexible multi-link manipulators. International Journal
of Robotics Research, 14(2):91111, 1995.
[81] Turcic D.A. and Midha A. Dynamic analysis of elastic mechanism systems.
part i: applications. ASME Journal of Dynamic Systems, Measurement and
Control, 106(4):249254, 1984.
215
REFERENCIAS
[83] Carmelo di Castri and Arcangelo Messina. Exact modeling for control of
flexible manipulators. Journal of Vibration and Control, 18(10):15261551,
2012.
[84] Pan Dong, Zhao Yang, Zhang Yue, Wei Cheng, and E Wei. Dynamic mode-
ling and analysis of space manipulator considering the flexible of joint and
link. Advanced Materials Research, 823:270275, 2013.
[85] Zhenjie Qian, Dingguo Zhang, and Jun Liu. Recursive formulation for dy-
namic modeling and simulation of multilink spatial flexible robotic manipu-
lators. Advances in Mechanical Engineering, 5, 2013.
[88] Francisco Ramos de La Flor. Vibration Control Strategies for a Very Light-
weight One Degree-Of-Freedom Flexible Arm Built with Composite Mate-
rial. PhD thesis, Escuela Tecnica Superior de Ingenieros Industriales de la
Universidad de Castilla-La Mancha, 2009.
[89] Yuan-Qing Xia, Meng-Yin Fu, Zhi-Hong Deng, and Xue-Mei Ren. Recent
developments in sliding mode control and active disturbance rejection con-
trol. Control Theory & Applications, 30(2):137147, 2013.
[91] Zhiqiang Gao, Yi Huang, and Jingqing Han. An alternative paradigm for
control system design. In Decision and Control, 2001. Proceedings of the
40th IEEE Conference on, volume 5, pages 45784585. IEEE, 2001.
216
REFERENCIAS
[94] Lili Dong and David Avanesian. Drive-mode control for vibrational mems
gyroscopes. Industrial Electronics, IEEE Transactions on, 56(4):956963,
2009.
[96] K. Ohnishi, M. Shibata, and Murakami T. Motion control for advanced me-
chatronics. Mechatronics, IEEE/ASME Transactions on, 1(1):5667, March
1996.
[100] Xiaoxia Yang and Yi Huang. Capabilities of extended state observer for
estimating uncertainties. In American Control Conference, 2009. ACC09.,
pages 37003705. IEEE, 2009.
217
REFERENCIAS
[104] Shigehiko Yamamoto and Iori Hashimoto. Present status and future needs:
the view from japanese industry. In Chemical Process Control-CPCIV. Pro-
ceedings of the Fourth International Conference on Chemical Process Con-
trol, Texas, 1991.
[105] W.L. Bialkowski. Dreams versus reality: a view from both sides of the gap:
manufacturing excellence with come only through engineering excellence.
Pulp & Paper Canada, 94(11):1927, 1993.
[107] Concepcion A Monje, Blas M Vinagre, Vicente Feliu, and YangQuan Chen.
Tuning and auto-tuning of fractional order controllers for industry applica-
tions. Control Engineering Practice, 16(7):798812, 2008.
[108] YangQuan Chen, Ivo Petras, and Dingyu Xue. Fractional order control-
a tutorial. In American Control Conference, 2009. ACC09., pages 1397
1411. IEEE, 2009.
[110] Alain Oustaloup, Benot Mathieu, and Patrick Lanusse. The crone control
of resonant plants: application to a flexible transmission. European Journal
of control, 1(2):113121, 1995.
218
REFERENCIAS
[112] Alain Oustaloup, Jocelyn Sabatier, and Xavier Moreau. From fractal robust-
ness to the crone approach. In Proc., ESAIM, volume 5, pages 177192,
1998.
[113] Patrick Lanusse, Valerie Pommier, and Alain Oustaloup. Fractional control
system design for a hydraulic actuator. In Proc. of 1st IFAC Conference on
Mechatronics Systems, Mechatronic 2000, 2000.
[114] A Oustaloup, P Melchior, P Lanusse, O Cois, and F Dancla. The crone tool-
box for matlab. In Computer-Aided Control System Design, 2000. CACSD
2000. IEEE International Symposium on, pages 190195. IEEE, 2000.
[116] Ying Luo and YangQuan Chen. Fractional order motion controls. John
Wiley & Sons, 2012.
[117] Chengbin Ma and Yoichi Hori. Fractional-order control: theory and applica-
tions in motion control. IEEE Industrial Electronics Magazine, 1(4):616,
2007.
[119] Alain Oustaloup and Muriel Bansard. First generation crone control. In
Systems, Man and Cybernetics, 1993.Systems Engineering in the Service
of Humans, Conference Proceedings., International Conference on, pages
130135. IEEE, 1993.
219
REFERENCIAS
[120] Jocelyn Sabatier, Patrick Lanusse, Pierre Melchior, and Alain Oustaloup.
Fractional order pid and first generation crone control system design. In
Jocelyn Sabatier, Patrick Lanusse, and Alain Oustaloup, editors, Fractio-
nal Order Differentiation and Robust Control Design. Springer Netherlands,
2015.
[124] S.S. Haykin. A unified treatment of recursive digital filtering. IEEE Tran-
sactions Automatic Control, AC-17:113116, 1972.
[126] Dar-Zen Chen and Lung-Wen Tsai. The generalized principle of inertia
match for geared robotic mechanisms. In Proceedings of the 1991 IEEE
International Conference on Robotics and Automation, pages 12821287,
Sacramento, California, 1991.
[127] Ramos F., Feliu V., and Payo I. Design of trajectories with physical cons-
traints for very lightweight single link flexible arms. Journal of Vibration
and Control, 14(8):10911110, August 2008.
220
REFERENCIAS
[129] William B Gevarter. Basic relations for control of flexible vehicles. AIAA
journal, 8(4):666672, 1970.
[130] Fei-Yue Wang and Yanqing Gao. Advanced Studies of Flexible Robotic Ma-
nipulators: Modeling, Design, Control, and Applications, volume 4. World
Scientific, 2003.
[131] Jung Hua Yang, Feng-Li Lian, and Li Chen Fu. Nonlinear adaptive control
for flexible-link manipulators. Robotics and Automation, IEEE Transactions
on, 13(1):140148, Feb 1997.
[133] Jorge J Feliu, Vicente Feliu, and Carlos Cerrada. Load adaptive control of
single-link flexible arms based on a new modeling technique. Robotics and
Automation, IEEE Transactions on, 15(5):793804, 1999.
[135] Masahiro Isogai, Fumihito Arai, and Toshio Fukuda. Modeling and vibration
control with neural network for flexible multi-link structures. In Robotics
and Automation, 1999. Proceedings. 1999 IEEE International Conference
on, volume 2, pages 10961101. IEEE, 1999.
[136] An-Chyau Haung and Kuo-Kai Liao. Fat-based adaptive sliding control for
flexible arms: Theory and experiments. Journal of Sound and Vibration,
298(1):194205, 2006.
[137] S-B Choi, C-C Cheong, and H-C Shin. Sliding mode control of vibration in
a single-link flexible arm with parameter variations. Journal of Sound and
Vibration, 179(5):737748, 1995.
221
REFERENCIAS
[138] V Feliu, JA Somolinos, Carlos Cerrada, and Jose Antonio Cerrada. A new
control scheme of single-link flexible manipulators robust to payload chan-
ges. Journal of Intelligent and Robotic Systems, 20(2-4):349373, 1997.
[139] S.S. Ge, T.H. Lee, and G. Zhu. Energy-based robust controller design for
multi-link flexible robots. Mechatronics, 6(7):779 798, 1996.
[140] Ioan Dore Landau, Jochen Langer, Daniel Rey, and Jean Barnier. Robust
control of a 360 flexible arm using the combined pole placement/sensitivity
function shaping method. Control Systems Technology, IEEE Transactions
on, 4(4):369383, 1996.
[141] Jamal Daafouz, Germain Garcia, and Jacques Bernussou. Robust control of
a flexible robot arm using the quadratic d-stability approach. IEEE Transac-
tions on Control Systems Technology, 6(4):524533, 1998.
[145] Emiliano Pereira, Sumeet S Aphale, Vicente Feliu, and SO Reza Moheimani.
Integral resonant control for vibration damping and precise tip-positioning of
a single-link flexible manipulator. Mechatronics, IEEE/ASME Transactions
on, 16(2):232240, 2011.
[146] Abdullah Al-Mamun, Ehsan Keikha, Charanjit Singh Bhatia, and Tong Heng
Lee. Integral resonant control for suppression of resonance in piezoelectric
222
REFERENCIAS
[147] Feliu V. and Ramos F. Strain gauge based control of single-link flexible very
lightweight robots robust to payload changes. Mechatronics, 15:547571,
2005.
[148] Vicente Feliu, Fernando J. Castillo, Francisco Ramos, and Jose A. Somoli-
nos. Robust tip trajectory tracking of a very lightweight single-link flexible
arm in presence of large payload changes. Mechatronics, 22(5):594 613,
2012.
[149] Michel Fliess, Richard Marquez, Emmanuel Delaleau, and Hebertt Sira-
Ramrez. Correcteurs proportionnels-integraux generalises. ESAIM: Con-
trol, Optimisation and Calculus of Variations, 7:2341, 2002.
[150] Adrian-Rene Ramrez Lopez. Analisis del desempeno del controlador pro-
porcional integral generalizado en un rectificador activo multinivel mo-
nofasico. Msc., Universidad Autonoma de San Luis de Potos, Diciembre
2010.
[152] Alfonso Parreno Torres, Pedro Roncero-Sanchez, Xavier Del Toro Garca,
and Vicente Feliu Batlle. Generalized proportional-integral control for
voltage-sag compensation in dynamic voltage restorers. In Compatibi-
lity and Power Electronics (CPE), 2011 7th International Conference-
Workshop, pages 9297. IEEE, 2011.
223
REFERENCIAS
[155] Jonathan Becedas, Juan Ramon Trapero, Hebertt Sira-Ramrez, and V Feliu-
Battle. Fast identification method to control a flexible manipulator with pa-
rameter uncertainties. In Robotics and Automation, 2007 IEEE International
Conference on, pages 34453450. IEEE, 2007.
[157] Vicente Feliu, Kuldip S Rattan, and H Benjamin Brown. Control of a two-
degree-of-freedom lightweight flexible arm with friction in the joints. Jour-
nal of robotic systems, 12(1):117, 1995.
[159] V. Feliu, J.A. Somolinos, and A. Garcia. Inverse dynamics based control sys-
tem for a three-degree-of-freedom flexible arm. Robotics and Automation,
IEEE Transactions on, 19(6):10071014, Dec 2003.
[160] Eduardo Bayo, Philip Papadopoulos, James Stubbe, and Miguel Angel Ser-
na. Inverse dynamics and kinematics of multi- link elastic robots: An ite-
rative frequency domain approach. The International Journal of Robotics
Research, 8(6):4962, 1989.
[161] Takahiko Tsujisawa and Wayne J Book. A reduced order model deriva-
tion for lightweight arms with a parallel mechanism. In Robotics and Auto-
mation, 1989. Proceedings., 1989 IEEE International Conference on, pages
728735. IEEE, 1989.
224
REFERENCIAS
[162] Haruhiko Asada, Zeng-Dong Ma, and Hidekatsu Tokumaru. Inverse dyna-
mics of flexible robot arms for trajectory control. In ASME Winter Annual
Meeting, Proceedings of the, pages 329336, 1987.
[163] Asada H., Ma Z.-D., and Tokumaru H. Inverse dynamics of flexible robot
arms: Modeling and computation for trajectory control. Journal of Dynamic
Systems Measurement and Control, 112(2):177185, 1990.
[164] De Luca A., Lucibello P., and Ulivi G. Inversion techniques for trajectory
control of flexible robot arms. Journal of Robotic Systems, 6(4):325344,
1989.
[165] E. Bayo. Computed torque for the position control of open-chain flexible
robots. In Robotics and Automation, Proceedings of the 1988 IEEE Interna-
tional Conference on, pages 316321 vol.1, Apr 1988.
[167] John F. Calvert and Donald J. Gimpel. Method and apparatus for control
of system output in response to system input, July 30 1957. US Patent
2,801,351.
[169] C.F. Cutforth and L.Y. Pao. Adaptive input shaping for manoeuvring flexible
structures. Automatica, 4(40):685693, 2004.
[170] Park J. and Chang P. H. Learning input shaping technique for non-lti
systems. ASME Journal of Dynamic Systems, Measurement, and Control,
2:288293, 2001.
[171] J. Park, P. H. Chang, H.S. Park, and E Lee. Design of learning input shaping
technique for residual vibration suppression in an industrial robot. IEEE/AS-
ME Transactions on Mechatronics, 11(1):5565, 2006.
225
REFERENCIAS
[172] Jianda Han, Zhiqiang Zhu, Yuqing He, and Juntong Qi. A novel input
shaping method based on system output. Journal of Sound and Vibration,
335(0):338 349, 2015.
[173] E. Pereira, J.R. Trapero, I.M. Daz, and V. Feliu. Adaptive input shaping for
manoeuvring flexible structures using an algebraic identification technique.
Automatica, 45(4):1046 1051, 2009.
[174] M. Bodson. An adaptive algorithm for the tuning of two input shaping met-
hods. Automatica, 34(6):771776, 1998.
[175] Joshua Vaughan, Aika Yano, and William Singhose. Comparison of robust
input shapers. Journal of Sound and Vibration, 315(45):797 815, 2008.
[176] Pereira E., Becedas J., Payo I., Ramos F., and Feliu V. Control of flexible
manipulators. In Agustin Jimenez and Basil M. Al Hadithi, editors, Robot
Manipulators Trends and Development. InTech, 2004.
[177] Eric Schmitz. Experiments on the End-Point Position Control of a Very Fle-
xible One-Link Manipulator. phd, Stanford University, Stanford, CA 94305,
June 1985.
[178] Paul T Kotnik, Stephen Yurkovich, and Umit Ozguner. Acceleration feed-
back for control of a flexible manipulator arm. Journal of Robotic Systems,
5(3):181196, 1988.
[179] Jorn Malzahn and Torsten Bertram. Collision detection and reaction for a
multi-elastic-link robot arm?. In IFAC World Congress, Cape Town (South
Africa), 2014.
[180] Jorn Malzahn, Anh Son Phung, and Torsten Bertram. Scene adaptive rgb-d
based oscillation sensing for a multi flexible link robot arm in unstructu-
red dynamic environments. In Intelligent Robots and Systems (IROS), 2012
IEEE/RSJ International Conference on, pages 10171022. IEEE, 2012.
[181] Ho-Cheol Shin and Seung-Bok Choi. Position control of a two-link flexible
manipulator featuring piezoelectric actuators and sensors. Mechatronics,
11(6):707 729, 2001.
226
REFERENCIAS
[182] Jun-qiang Lou, Yan-ding Wei, Yi-ling Yang, and Feng-ran Xie. Hybrid pd
and effective multi-mode positive position feedback control for slewing and
vibration suppression of a smart flexible manipulator. Smart Materials and
Structures, 24(3), 2015.
[183] Dunant Halim, Xi Luo, and Pavel M Trivailo. Decentralized vibration con-
trol of a multi-link flexible robotic manipulator using smart piezoelectric
transducers. Acta Astronautica, 104(1):186196, 2014.
[185] David Magee and Wayne Book. Control and control theory for flexible ro-
bots. Journal of the Society of Instrument and Control Engineers, 32(4),
1993.
[186] Eric Schmitz. Experiments on the End-Point Position Control of a Very Fle-
xible One-Link Manipulator. phd, Stanford University, Stanford, CA 94305,
June 1985. Also published as SUDAAR 548.
[187] Jahng-Hyon Park and Haruhiko Asada. Design and analysis of flexible arms
for minimum-phase endpoint control. In American Control Conference,
1990, pages 12201225, 1990.
[188] Wang W.J., Lu S.S., and Hsu C.F. Experiments on the position control of a
one-link flexible robot arm. IEEE Transactions on Robotics and Automation,
5(3):373377, 1989.
[189] Tokhi M.O. and Azad A.K.M. Control of flexible manipulator systems. Pro-
ceedings of the Institution of Mechanical Engineers, Part I: Journal of Sys-
tems and Control Engineering, 210(2):113130, 1996.
[190] Lanari L. and Wen J. Asymtotically stable set point control laws for flexible
robots. Systems and Control Letters, 19:119129, 1992.
[191] M.-T. Ho and Y.-W. Tu. Position control of a single-link flexible manipulator
using h/spl infin/-based pid control. Control Theory and Applications, IEE
Proceedings, 153(5):615622, Sept 2006.
227
REFERENCIAS
[193] Levent Gumusel and Nurhan Gursel Ozmen. Modelling and control of ma-
nipulators with flexible links working on land and underwater environments.
Robotica, 29:461470, 5 2011.
[196] B Siciliano, JVR Prasad, and AJ Calise. Output feedback two-time scale
control of multilink flexible arms. Journal of Dynamic Systems, Measure-
ment, and Control, 114(1):7077, 1992.
[197] D.S. Kwon and W.J. Book. A time-domain inverse dynamic tracking control
of a single-link flexible manipulator. Journal of Dynamic Systems, Measu-
rement and Control, 116(2):193200, 1994.
[199] V. Moudgal, W.A. Kwong, and K. Passino. Fuzzy learning control for a
flexible-link robot. IEEE Transactions on Fuzzy Systems, 3(2):199210,
1995.
[200] T.H. Lee, S.S. Ge, and Z.P. Wang. Adaptive robust controller design for
multi-link flexible robots. Mechatronics, 11(8):951 967, 2001.
[201] M. Bai, D.H. Zhou, and H. Schwarz. Adaptive augmented state feedback
control for an experimental planar two-link flexible manipulator. Robotics
and Automation, IEEE Transactions on, 14(6):940950, 1998.
228
REFERENCIAS
[203] Tao Fan. Intelligent Model Predictive Control of Flexible Link Robotic Ma-
nipulators. phd, The University of British Columbia, Mayo 2007.
229