Al comienzo de la teórica 2 reciclamos el código Python que habíamos hecho en la teórica 1 para un péndulo ideal, para el caso de un péndulo con rozamiento sin la aproximación de ángulos pequeños. Hicimos la simulación numérica con un valor del ángulo inicial igual a 3 radianes (casi coincidente con la vertical) y con un empujón que generaba una velocidad inicial de 10 1/s. En estas condiciones y sin rozamiento, el péndulo se quedaría dando vueltas eternamente. Pero debido a la fuerza de rozamiento, sólo da N vueltas y llega un momento en que el comportamiento del péndulo comienza a ser idéntico al predicho por la aproximación de ángulos pequeños.
En la simulación que vimos en clase, el número de vueltas era N=6. Ahora les dejo el resultado de una simulación con los mismos valores, excepto que bajé el valor del coeficiente de rozamiento. En consecuencia, el número de vueltas aumenta y ahora toma el valor N=11. El gráfico de posición en función del tiempo es el siguiente
Y para que vean que son 11 vueltas, acá dejo la animación (click en imagen)
El desafío consiste en encontrar las condiciones iniciales para que el péndulo realice N vueltas, con N arbitrario.