Péndulo Invertido
Péndulo Invertido
Péndulo Invertido
on y de Estudios Avanzados
del IPN
Departamento de Ingenier
a El
ectrica
Secci
on de Mecatr
onica
Maestr
a en Ciencias opci
on Mecatr
onica
Modelado de Sistemas Din
amicos
Asignaci
on:
Modelado, Control y Simulaci
on
de un P
_
R
6
(1)
Denici on 2 (Energa Cinetica) La energa cinetica del cuerpo A con masa constante
m se dene y denota como:
k(q, q) =
1
2
n
i=1
m
i
v
2
i
(2)
Denicion 3 (Energa Potencial) la energa potencial es energa que mide la capa-
cidad que tiene dicho sistema para realizar un trabajo en funcion exclusivamente de su
posicion o conguracion. Puede pensarse como la energa almacenada en el sistema, o
como una medida del trabajo que un sistema puede entregar.
U = U(q) (3)
Dicha relacion entre las ecuaciones 2 y 3 est a dada por el Lagrangiano, el cual est a denido
como:
Denici on 4 (Lagrangiano) El lagrangiano de un sistema mecanico esta dado como la
diferencia de su energa cinetica y energa potencia, es decir,
L = L(q, q) = k(q, q) U(q) (4)
La relaci on descrita por la ecuaci on 4 es utilizada para encontrar las ecuaciones de movi-
miento de cualquier sistema mecanico (en este el pendulo invertido de la gura 1) mediante
la aplicacion de las ecuaciones las llamadas Ecuaciones de Euler Lagrange.
Denici on 5 (Ecuaciones de Euler Lagrange[3]) Las ecuaciones que describen el mo-
vimiento de un sistema mecanico de n grados de libertad
2
vienen dadas por:
d
dt
L
q
k
L
q
k
=
k
(5)
donde
k
son las fuerzas y pares ejercidos externamente y k = 1, ..., n.
Aplicando lo lo descrito al pendulo invertido, se tiene lo siguiente. Primero que nada se
establece un sistema de referencia como el mostrado en la gura 1. El sistema requiere
de un par ametro para caracterizar la posicion angular de la masa respecto al eje de giro
del pendulo sobre el carro, en el mismo orden de ideas, se requiere un parametro para
describir la posicion del carro sobre el plano en una sola direccion, por lo que se toma
como coordenada generalizada a q = [x ]
T
.
2
El n umero de grados de libertad es la cardinalidad del espacio de conguracion Q
2
Cinvestav-IPN Maestra en Ciencias opcion Mecatronica
Para la masa M.
La posici on esta descrita por x
M
= x, su velocidad es x
M
= x, por lo que la energa
cinetica es:
k
M
=
1
2
M x
2
(6)
Su energa potencial respecto al sistema de referenciay considerando la entrada es:
U
M
= Fx (7)
donde F es la fuerza aplicada al carro para que se mueva.
Para la masa m.
Su posicion:
x
m
= x + l sin x
m
= x + l
cos
y
m
= l cos y
m
= l
sin
Por lo anterior se tiene la velocidad al cuadr ado dada por:
x
2
m
+ y
2
m
= x
2
+ 2l x
cos + l
2
2
(8)
La energa cinetica para m es:
k
m
=
1
2
m( x
2
+ 2l x
cos + l
2
2
) (9)
La energa potencial viene dada por:
U
m
= mgl cos (10)
Energa cinetica y Potencia total
La cinetica total esta dada por la suma de las expresiones 6 u 9:
k =
1
2
(M + m) x
2
+ ml x
cos +
1
2
ml
2
2
(11)
La energa Potencial total esta dada por la suma de 7 y 10:
U = mgl cos Fx (12)
El Lagrangiano del sistema es la diferencia de la ecuaci on 11 y 12:
L =
1
2
(M + m) x
2
+ ml x
cos +
1
2
ml
2
2
mgl cos + Fx (13)
Como uno de los objetivos es simular el sistema, entonces el tener diferentes unidades de
medici on implica m as cuidado y gasto, por lo quese normaliza el sistema de tal forma que
se obtenga un sistema adimensional. Para lograr lo anterior propuesto se procede de la
siguiente manera.
3
Cinvestav-IPN Maestra en Ciencias opcion Mecatronica
Sea el tiempo normalizado (adimensional):
=
_
g
l
t d =
_
g
l
dt dt =
l
g
d (14)
de tal forma que
dx
dt
=
_
g
l
dx
d
y
d
dt
=
_
g
l
d
d
(15)
Sustituyendo las expresiones 15 en el lagrangiano dado por la ecuaci on 13, se tiene:
L =
1
2
(M+m)
__
g
l
dx
d
_
2
+ml
__
g
l
dx
d
___
g
l
d
d
_
cos +
1
2
ml
2
__
g
l
d
d
_
2
mgl cos +Fx
Simplicando
L =
1
2
(M +m)
g
l
_
dx
d
_
2
+ml
g
l
_
dx
d
__
d
d
_
cos +
1
2
ml
2
g
l
_
d
d
_
2
mgl cos +Fx (16)
Dividiendo todo por mgl, tenemos:
L =
1
2
_
M
m
+ 1
_
1
l
2
_
dx
d
_
2
+
1
l
dx
d
d
d
cos +
1
2
_
d
d
_
2
cos +
Fx
mgl
(17)
Considerando los cambios de variable
=
M
m
, =
x
l
, =
F
mg
Finalmente se tiene el Lagrangiano normalizado
L =
1
2
_
+ 1
_
2
+
cos +
1
2
2
cos + u (18)
Para la ecuacion 18, las coordenadas generalizadas normalizadas son y .
Ahora se determinar an las ecuaciones de movimiento a partir del Lagrangiano norma-
lizado.
4
Cinvestav-IPN Maestra en Ciencias opcion Mecatronica
Ecuaci on de Movimiento para coordenada , considerando la fuerza externa dentro
de la energa potencial.
d
dt
L
L
= 0 (19)
L
=
_
+ 1
_
+
cos
d
dt
L
=
_
+ 1
_
+
cos
2
sin
L
= u
Finalmente, al sustituir lo anterior en 19 se tiene la ecuaci on de movimiento para la
coordenada generalizada :
_
+ 1
_
+
cos
2
sin = u (20)
Ecuaci on de Movimiento para coordenada , analogamente a .
d
dt
L
= 0 (21)
L
= cos +
d
dt
L
= cos
sin +
sin + sin
Luego, al sustituir lo anterior en 21 se tiene la ecuaci on de movimiento para la coordenada
generalizada :
2
sin = u (23)
=
...
F
= F
...
F
= F
(4)
=
F F
(4)
(28)
Sustituyendo la parametrizaci on en las expresiones 28 en las ecuaciones linealizadas 26:
( + 1)
F F
(4)
= u (29)
0 = 0 (30)
3
Planitud diferencial es una propiedad de algunos sistemas dinamicos, que permite trivializar las tareas
de planicacion de trayectoria. Tambien permite simplicar el control por retroalimentacion estado a un
conjunto de sistemas desacoplados.
6
Cinvestav-IPN Maestra en Ciencias opcion Mecatronica
El sistema de ecuaciones anterior es isomorfo al sistema de ecuaciones 26 y 25 normali-
zadas linealizado del pendulo invertido. Con lo anterior se puede ver que la dinamica del
pendulo alrededor de su punto de equilibrio y la parametrizaci on 28, esta completamente
descrita por la ecuaci on 29.
Ahora, supongamos que quiere que la variable F del sistema descrito por la ecuaci on
29 siga a una trayectoria deseada F
sea lo
m as peque na posible. A dicha diferencia se le llamara error denido como:
F = F
F (31)
Despejando F
(4)
de la ecuaci on 29, se tiene:
F
(4)
=
_
1 +
1
_
F
1
u (32)
Se propone la ley de control siguiente, basada en la tecnica de control por par calculado[3].
u =
__
1 +
1
_
F F
(4)
k
3
(
...
F
...
F
) + k
2
(
F
F) + k
1
(
F
F) + k
0
(F
F)
_
(33)
Donde k
i
> 0, i = 0, 1, 2, 3. Al tomar las derivadas necesarias de la expresi on 31 y
sustituyendo el resultado en la expresion del control 33, se tiene:
u =
__
1 +
1
_
F F
(4)
k
3
...
F k
2
F k
1
F k
0
F
_
(34)
Luego, colocando el control dado por la expresion 34 en el modelo 32 y al hacer algunas
simplicaciones, se tiene:
F
(4)
= F
(4)
+ k
3
...
F + k
2
F + k
1
F + k
0
F (35)
Pasando el termino del miembro izquierdo de la ecuaci on 35 y considerando que
F
(4)
=
F
(4)
F
(4)
, se tiene un sistema equivalente, cuyas variables son los errores de lso estados
del sistema de la salida plana.
F
(4)
+ k
3
...
F + k
2
F + k
1
F + k
0
F = 0 (36)
El sistema anterior es asint oticamente estable
4
, y basta hacer una sintonizaci on de las
ganancias, lo cual se hace a partir de su polinomio caracterstico:
s
4
+ k
3
s
3
+ k
2
s
2
+ k
1
s + k
0
= 0 (37)
4
Es posible establecer las condiciones de las ganancias de tal forma que el sistema sea asintoticamente
estable a cero. Se puede hacer desde el enfoque de funcion de transferencia mediante el criterio de Routh-
Hurwitz o el enfoque de variables de estado usando el segundo metodo de Lyapunov.
7
Cinvestav-IPN Maestra en Ciencias opcion Mecatronica
Con el n de poder relacionar el sistema con factor de amortiguamiento relativo
y frecuencia natural w
n
y as poder manipular la ubicacion de los polos del sistema, se
establece la siguiente relaci on entre la ecuaci on caracterstica de la funci on de transferencia
del modelo general de un sistema de segundo orden y la ecuaci on caracterstica 37 como
sigue:
s
4
+ k
3
s
3
+ k
2
s
2
+ k
1
s + k
0
= (s
2
+ 2w
n
s + w
2
n
)
2
(38)
De tal forma que se establece la relacion entre el factor de amortiguamiento relativo y la
frecuencia natural del sistema con las constantes de dise no.
k
3
= 4w
n
k
2
= 2w
2
n
+ 4
2
w
2
n
k
1
= 4w
3
n
k
0
= w
4
n
4. Simulaci on
Para tener una aproximaci on de lo que ocurrira en la realidad al implementar el
pendulo invertido con el esquema de control 34 se realiza una simulacion numerica. La
simulaci on se har a en Simnon .
Como se mencion o en la parte de modelado se ha considerado una dinamica normalizada
con nes de simplicacion de operaciones y facilidad de an alisis. El esquema de control
dise nado previamente es dise nado para un sistema lineal, pero como tambien se men-
cion o en la parte de control, el rango de operaci on y u es una vecindad del equilibrio,
entonces la ley de control debe de funcionar para el modelo no lineal. Dada la forma de
dar instrucciones en el simulador, se debe escribir el modelo en formato de variables de
estado.
z
1
= z
2
z
2
=
1
_
u + z
2
4
sin(z
3
) sin(z
3
) cos(z
3
)
z
3
= z
4
z
4
=
1
_
( + 1) sin(z
3
) ucos(z
3
) z
2
4
sin(z
3
) cos(z
3
)
12
Cinvestav-IPN Maestra en Ciencias opcion Mecatronica
Finalmente en la gura 10 se muestra la salida plana F y la salida plana deseada F