Péndulo con resorte, lagrangiano y hamiltoniano
Simulación de Matlab
Sede Veraguas, Veraguas, Segundo Semestre 2024
Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Dinámica de Sistemas Electromecánicos
Simulación en Matlab de péndulo con resorte
Estudiantes:
Sebastián De León 9-765-1563
Josue Gaitán 9-764-1748
Energía Potencial en el Resorte
En este caso de un péndulo se nos presenta con un resorte, el
radio no será constante, y esta variación conlleva a por que se
tenga más de una coordenada generalizada.
Podemos describir la energía potencial asociada al estiramiento
del resorte como:
En el cual podemos observar que expresa la contracción del
resorte, expresado de la siguiente forma:
Dándonos como resultado:
Energía Potencial Gravitatoria
Para obtener la energía potencial se tomará como referencia el
origen del péndulo debido a la variación del radio, lo que nos
termina dando:
Donde es la distancia de la masa con respecto al origen, si
embargo esta tiene una dependencia de r y va con dirección al eje
negativo de las ordenadas.
por lo tanto:
Energía Cinética
Para que podamos obtener la energía cinética necesitamos la velocidad de la masa.
Donde es la velocidad del sistema en el eje x y es la velocidad del sistema en el eje de las y, reescribiremos X y Y en
coordenadas polares para así obtener nuestras coordenadas generalizadas
Al derivar respecto al tiempo para obtener la velocidad tendremos:
Posteriormente al elevarlas al cuadrado tendremos:
Después de obtener los valores requeridos se procede a reemplazar:
Al remplazar dentro de T, tendremos:
Lagrangiano
Ahora para calcular el Lagrangiano remplazando T y V en
la ecuación tendremos tendremos que:
Ecuación de Lagrange para r
Para la coordenada generalizada r tendremos:
Al derivar parcialmente respecto a r tendremos:
Al derivar parcialmente respecto a tendremos:
Luego derivamos respecto al tiempo la derivada parcial de :
Finalmente remplazando en nuestra ecuación de Lagrange para r
tendremos:
Ecuación de Lagrange para
Para la coordenada generalizada tendremos:
Al derivar parcialmente respecto a tendremos:
Al derivar parcialmente respecto a tendremos:
Luego derivamos respecto al tiempo la derivada parcial de :
Finalmente remplazando en nuestra ecuación de Lagrange para
tendremos:
Momentos Canónicos y Hamiltonaino
Momentos Canónicos
El Hamiltoniano
Ecuaciones de Movimientos de Hamilton
Las ecuaciones de movimiento en el formalismo de Hamilton son:
, ;
Se sustituye el Hamiltoniano Se sustituye el Hamiltoniano
obtenido anteriormente, para r obtenido anteriormente, para
obtenemos: obtenemos:
CODIGO EN MATLAB PARA SIMULAR EL PENDULO CON
RESORTE
• Como primeros pasos para la
creación del código pedimos la
entrada de los perímetros
• Que las longitudes sean validas
• nombramos las constantes como
la gravedad, la constante del
resorte y condiciones iniciales
• Como segundo paso resolvemos
las ecuaciones de movimiento
utilizando ode45
• Animación de la simulación
• Dibujar los elementos y el bucle
de animación
INTERFAZ DEL PROGRAMA
Simulación