Control de un balancín aerodinámico
Proyecto Final - Métodos Numéricos
Sergio Agudelo$^1$
Jueves 27 de Junio de 2013
Universidad Nacional de Colombia
$^1$Código: 285675, e-mail: seagudelobe@unal.edu.co
Introducción
In [2]:
from IPython.core.display import Image
Image(filename='Sketch.png')
Out[2]:
El balancín consta de dos contrapesos, uno de los cuales es capaz de ejercer una fuerza estabilizadora, gracias a la hélice y motor que lleva montados. Este sistema es básicamente un péndulo, con la adición de que puede mantenerse en una posición diferente a la horizontal.
Para medir la posición en la que se encuentra la barra transversal, es usado un potenciómetro (más adelante se caracterizará).
Modelo Matemático
La ecuación dinámica del sistema es:
En donde:
$l_1:=\ $ Distancia desde la masa izquierda hasta el centro
$l_2:=\ $ Distancia desde la masa derecha hasta el centro
$W_m:=\ $ Fuerza gravitatoria relativa al conjunto motor-hélice
$F_e:=\ $ Fuerza de empuje generada por la hélice
$W_c:=\ $ Fuerza gravitatoria relativa al contrapeso
$\beta_1:=\ $ Coeficiente de rozamiento relativo al conjunto motor-hélice
$\beta_2:=\ $ Coeficiente de rozamiento relativo al contrapeso
$I:=\ $ Momento de inercia de la sección móvil del montaje
$\alpha:=\ $ Aceleración angular
El sistema se limitará a una zona de operación en que se cumple $sin\:\theta \approx\theta$. Entonces, la siguiente aproximación es correcta:
Reescribiendo los términos como diferenciales:
Para que sea una EDO se debe hacer un cambio de variables, con el fin de que todas las ecuaciones diferenciales resultantes sean de primer orden.
Ahora, si se tiene en cuenta el valor numérico de las constantes usadas:
El sistema resultante es:
[ \left{ \begin{array}{l l} x = \frac{d\theta}{dt}\\ \\ \frac{dx}{dt} = \frac{l_1(y - W_m\:\theta - \beta_1x) + l_2(W_c\:\theta-\beta_2x)}{I} \end{array} \right.]
In [ ]: