el Péndulo
Nombre: Nikki braguero
Número de Estudiantes: 09369481
Abstracto:
En este experimento, modelamos el movimiento de un péndulo simple, un oscilador amortiguado, y
un oscilador amortiguado accionado, utilizando aproximaciones lineales y no lineales de las
ecuaciones de movimiento. Los datos obtenidos se representan a continuación en gnuplot. Se
concluyó que la aproximación lineal es siempre lo suficientemente precisa para pequeños ángulos
iniciales, pero para ángulos iniciales más grandes se necesitaban las aproximaciones no lineales.
También se concluyó que, si bien utilizando la regla trapezoidal para obtener las ecuaciones de
movimiento era aceptable para el péndulo simple, que esto no fue el caso para el péndulo
amortiguado, para lo cual el método de Runge-Kutta era mucho más precisa. Finalmente, el período
de duplicacióny comportamiento caótico que se produce como la amplitud de la fuerza
motriz de una amortiguada del péndulo impulsado se incrementa, se observó a través de
retratos de fase.
Objetivos:
Los objetivos de este experimento fueron;
Para investigar la ecuación de movimiento de un péndulo simple
Para investigar la ecuación de movimiento de un péndulo amortiguado
Para investigar el comportamiento de un péndulo amortiguado accionado
Para investigar las soluciones numéricas de las ecuaciones diferenciales usando la regla
trapezoidal y el método de Runge-Kutta
Para aprender cómo escribir los datos de salida en un archivo de programación C.
Para aprender a utilizar gnuplot para crear gráficos a partir de los ficheros de datos
Introducción y antecedentes Teoría:
El péndulo simple
De acuerdo con la segunda ley de Newton del movimiento, la ecuación para el péndulo es
es el componente de la aceleración debida a la
Donde L es la longitud del péndulo y la
gravedad en la dirección hacia abajo.
Esta es una ecuación no lineal, pero puede ser linealizado mediante el uso de la aproximación de
que, para la pequeña
anglos, . Entonces obtenemos la siguiente,
Esta ecuación se puede resolver para la solución habitual para un oscilador armónico,
Donde A es la amplitud, y es la constante de fase (inicial).
Ecuaciones diferenciales de segundo orden
Para resolver ecuaciones diferenciales de segundo orden numérico, se debe introducir una nueva
variable y transformar la ecuación de segundo orden en dos ecuaciones diferenciales de primer
orden. Esto se puede realizar en los dos casos lineales y no lineales, utilizando la velocidad angular
de la sacudida,, que se define como
El método de Euler simple
Los métodos de Euler para resolver las ecuaciones diferenciales simples péndulo consiste en elegir
los valores iniciales de θ ω y, a continuación, utilizando una expansión de Taylor para encontrar los
valores de estos después de algún tiempo ha pasado. Sin tener en cuenta los términos de segundo
orden o superior, nos quedamos con las siguientes expansiones;
Hay términos de error de orden y más alto.
La regla del trapezoide
La regla trapezoidal hace mejores aproximaciones para la actualización de la frecuencia angular y
el ángulo, mediante la aproximación del área bajo la curva como un trapezoide. Esto es exacto el
segundo orden para un paso de tiempo h, pero esto no es lo suficientemente precisa para un
péndulo amortiguado accionado.
El-Kutta Runge Método
La cuarta orden algoritmo de Runge-Kutta es una aproximación lineal dada por
En estas ecuaciones, la k de son los valores de la pendiente de la función en puntos a través del
intervalo.
Este algoritmo es exacta a la 4º ordenar términos, por lo que se llama el cuarto orden de Runge-
Kutta.
El oscilador amortiguado Driven
El modelo anterior es sólo exacta para un péndulo simple, sin embargo esto puede ser extendido
para tener en cuenta
para la amortiguación y la conducción de fuerzas, dada por - y respectivamente. Esto transforma la ecuación de movimiento para el siguiente;
Donde k es un coeficiente de amortiguación, A es la amplitud de la fuerza motriz, y es la
frecuencia angular de esta fuerza.
Esta fórmula, una vez más puede ser linealizado utilizando la aproximación para ángulos pequeños
para obtener
Método experimental:
Un programa en C se utilizó para simular el sistema del péndulo, y para escribir los datos
en un archivo. A continuación, los datos se representan gráficamente. Esto se realizó para
un número de los casos;
i. El péndulo simple, tanto para el lineales y no lineales ecuaciones de movimiento
utilizando la regla trapezoidal
ii. El péndulo simple, tanto para el ecuaciones no lineales de movimiento lineal y el
º
uso de la 4 ordenar método de Runge-Kutta
iii. El oscilador amortiguado tanto para el ecuaciones no lineales de movimiento lineal
º
y el uso de la 4 ordenar método de Runge-Kutta
iv. El oscilador amortiguado impulsado tanto para el ecuaciones no lineales de
º
movimiento lineal y el uso de la 4 ordenar método de Runge-Kutta
En cada caso, las condiciones de partida que se utilizaron eran y
. Una variedad de gráficos han sido entonces trazada usando gnuplot para cada caso.
resultados:
Caso i) El péndulo simple que utiliza la regla del trapecio
Aproximación lineal:
Fig 1. Los gráficos de y para varios valores de partida
La aproximación lineal da valores adecuados para los pequeños ángulos iniciales, pero para
grandes ángulos de partida los resultados siguen siendo una sinusoide perfecta, que es
inexacta.
Aproximación no lineal:
Fig. 2 Los gráficos de y para varios valores de partida
Para ángulos de partida pequeños y grandes ángulos de partida se observa un
comportamiento muy diferente. La aproximación no lineal es mucho más preciso para
grandes ángulos iniciales.
Comparación de Linear y aproximaciones no lineales:
La Fig. 3 Comparación de lineal [En rojo] y no lineales [en verde]
En estos diagramas se puede observar que a medida que aumenta el ángulo inicial, la
aproximación lineal se vuelve cada vez más incorrecta. La aproximación lineal solamente
produce una imagen sinusoidal debido a la aproximación para ángulos pequeños, mientras
que la aproximación no lineal no lo hace. Esto puede verse, por ejemplo, en el último caso,
que está modelando el comportamiento cuando un péndulo simple es liberado de casi
vertical.
º
Caso ii) El péndulo simple que utiliza la 4 Orden de Runge-Kutta Método
Fig. 4 Gráficos de y para varios valores de partida
Cuando los datos obtenidos a partir de la aproximación no lineal Runge-Kutta se
representó gráficamente, y en comparación con la aproximación no lineal usando la regla
trapezoidal se observó que no había ninguna diferencia en la precisión para un péndulo
simple, como se muestra por la similitud entre las figuras arriba y aquellos para la
aproximación no lineal en el caso i.
º
Caso iii) El oscilador amortiguado utilizando la 4 Orden de Runge-Kutta Método
Fig. 5 Comparación de lineal [En rojo] y no lineales [En verde] con amortiguación
De la comparación de aproximaciones lineales y no lineales para el oscilador amortiguado,
podemos observar una vez más que para el ángulo de partida pequeño, la aproximación
lineal es suficientemente exacto, pero para ángulos iniciales más grandes esto se vuelve
menos exacto. Sin embargo, ambas aproximaciones finalmente convergen al mismo valor
final cerca de cero de equilibrio estático se alcanza, como se muestra a continuación.
La Fig. 6 Comparación de comportamiento lineal y no lineal sobre una escala de
tiempo más largo
º
iv caso) La amortiguada oscilador forzado utilizando la 4 Orden de Runge-Kutta Método
La Fig. 7 retratos de fase que muestra el comportamiento de un oscilador
amortiguado accionado
Las cifras anteriores son retratos de fase que describen el comportamiento de un
amortiguado impulsado
oscilador, con θ representa frente ω, por las siguientes amplitudes,. Se puede observar a
partir de los diagramas anteriores que cuando la amplitud es suficientemente grande, se
produce ese período de duplicación. Cuando la amplitud alcanza un valor de 1,5,
comportamiento caótico se puede observar durante un intervalo de tiempo más largo, como
se ve a continuación.
Fig. 8 Comportamiento caótico observa cuando A = 1,5
Discusiones y Conclusiones:
Se puede concluir a partir de las simulaciones anteriores que tanto para el péndulo simple
y el péndulo amortiguado, y para los pequeños ángulos iniciales, la aproximación lineal es
adecuada. Sin embargo, a medida que aumenta el ángulo, se necesitan modelos tales
como la aproximación no lineal para mantener la precisión, aunque para el oscilador
amortiguado tanto en el lineal y las aproximaciones no lineales convergido en el mismo
valor con el tiempo.
Se demostró que para ángulos pequeños, la regla trapezoidal y el método de Runge-Kutta
fueron igualmente exacto, pero para ángulos más grandes del método de Runge-Kutta era
un algoritmo muy superior. Utilizando el método de Runge-Kutta, era posible crear retratos
de fase para mostrar el periodo de duplicación y comportamiento caótico que se produce
como la amplitud de la fuerza de accionamiento de un péndulo amortiguado accionado está
aumentado.
