Simulación del péndulo simple y amortiguado

Usando ecuaciones diferenciales

Contenido

  1. Explorar el modelo matemático para el péndulo
    • Simple y Amortiguado
  2. Usar GeoGebra para modelar el movimento del péndulo.

Modelación

  • Barra rígida, ingrávida y no experimenta fricción.
  • Todo el peso se concentra en su masa sujeta al extremo.

Condiciones:

Modelación

Parámetros:

  • Masa mm.
  • Longitud de la barra LL.
  • Posición inicial.
  • Ángulo θ\theta.
  • Arco ss.
  • Tiempo tt.

Variables:

Constantes:

  • Gravedad gg.

Modelación

  • Aceleración angular
  • Aceleración gravitacional

Fuerza tangencial:

F=maF = m\cdot a
F = m\cdot a

Segunda Ley del Movimiento de Newton

  • Aceleración angular
θ=sL\theta = \dfrac{s}{L}
\theta = \dfrac{s}{L}
s=Lθs = L \cdot \theta
s = L \cdot \theta
s=Lθs'' = L \cdot \theta ''
s'' = L \cdot \theta ''
a=Lθa = L \cdot \theta ''
a = L \cdot \theta ''
ma=mLθm\cdot a = m\cdot L \cdot \theta ''
m\cdot a = m\cdot L \cdot \theta ''
F=ma=mLθF = m\cdot a = m\cdot L \cdot \theta ''
F = m\cdot a = m\cdot L \cdot \theta ''
F=mLθF = m\cdot L \cdot \theta ''
F = m\cdot L \cdot \theta ''
s=Lθs' = L \cdot \theta '
s' = L \cdot \theta '

Modelación

  • Aceleración angular
F=mLθF = m\cdot L \cdot \theta ''
F = m\cdot L \cdot \theta ''

Modelación

  • Aceleración gravitacional
F=mgsen θF = - m\cdot g \cdot \text{sen } \theta
F = - m\cdot g \cdot \text{sen } \theta

Modelación

  • Aceleración gravitacional
F=mgsen θF = - m\cdot g \cdot \text{sen } \theta
F = - m\cdot g \cdot \text{sen } \theta

Modelación

Aceleración angular y Aceleración gravitacional

F=mgsen θF = - m\cdot g \cdot \text{sen } \theta
F = - m\cdot g \cdot \text{sen } \theta
F=mLθF = m\cdot L \cdot \theta ''
F = m\cdot L \cdot \theta ''
mLθ=mgsen θm\cdot L \cdot \theta '' =- m\cdot g \cdot \text{sen } \theta
m\cdot L \cdot \theta '' =- m\cdot g \cdot \text{sen } \theta
θ=gLsen θ\theta '' =- \dfrac{g}{L} \text{sen } \theta
\theta '' =- \dfrac{g}{L} \text{sen } \theta
θ+gLsen θ=0\theta '' + \dfrac{g}{L} \text{sen } \theta =0
\theta '' + \dfrac{g}{L} \text{sen } \theta =0

y

Modelación

θ+gLsen θ=0\theta '' + \dfrac{g}{L} \text{sen } \theta =0
\theta '' + \dfrac{g}{L} \text{sen } \theta =0

Modelación

Aceleración angular y Aceleración gravitacional

θ+gLsen θ=0\theta '' + \dfrac{g}{L} \text{sen } \theta =0
\theta '' + \dfrac{g}{L} \text{sen } \theta =0

Posición del péndulo

xp=Lsen θyyp=Lcosθx_p = L \,\text{sen }\theta \quad \text{y} \quad y_p = -L \cos \theta
x_p = L \,\text{sen }\theta \quad \text{y} \quad y_p = -L \cos \theta

Movimiento del péndulo simple

Simulación

Resuelve numéricamente un sistema de ecuaciones diferenciales.

θ+gLsen θ=0\theta '' + \dfrac{g}{L} \text{sen } \theta =0
\theta '' + \dfrac{g}{L} \text{sen } \theta =0
ω=θ\omega = \theta '
\omega = \theta '
ω=θ\omega' = \theta ''
\omega' = \theta ''
    \implies
\implies
ResuelveNEDO()
{θ=ωω=gLsen θ\left\{ \begin{array}{rll} \theta ' &= &\omega \\ \omega' &=& - \dfrac{g}{L}\, \text{sen }\theta \end{array} \right.
\left\{ \begin{array}{rll} \theta ' &= &\omega \\ \omega' &=& - \dfrac{g}{L}\, \text{sen }\theta \end{array} \right.
θ+gLsen θ=0\theta '' + \dfrac{g}{L} \text{sen } \theta =0
\theta '' + \dfrac{g}{L} \text{sen } \theta =0
    \implies
\implies

Simulación

ResuelveNEDO()

Resuelve numéricamente un sistema de ecuaciones diferenciales.

Simulación

ResuelveNEDO(<List-der>, <Val-ini-x>, <List-val-ini-y>, <Val-fin-x>)
  1. Lista de derivadas: En nuesto caso un sistema de ecuaciones diferenciales.
  2. Valor inicial de xx: Aquí xx representa el tiempo t=0.t=0.
  3. Lista de valores iniciales de yy: Es decir θ0\theta_0 y ω0.\omega_0.
  4. Valor final de xx: Esto es el tiempo final t.t.

Simulación - Script

# Gravedad
g = 9.81
# Longitud de la barra
L = 2
# Posición inicial θ0 y velocidad inicial ω0
θ0 = 3 pi/4
ω0 = 0
# Sistema de ecuaciones diferenciales
θ'(t, θ, ω) = ω
ω'(t, θ, ω) = -g / L sen(θ)
# Resuelve el sistema
ResuelveNEDO[{θ', ω'}, 0, {θ0, ω0}, 17.3]
# Lo siguiente obtiene los valores de la posición
len = Longitud[IntegralNumérica1]
c = Deslizador[0, 1, 1 / len, 1, 100]
# Posición
xp = L sen(y(Punto[IntegralNumérica1, c]))
yp = -L cos(y(Punto[IntegralNumérica1, c]))
# Dibuja la masa y la barra que la sostiene
A = (xp, yp)
Segmento[(0, 0), A]
# Anima
IniciaAnimación()

Péndulo amortiguado

{θ=ωω=γmωgLsen θ\left\{ \begin{array}{rll} \theta ' &= &\omega \\ \omega' &=& -\dfrac{\gamma}{m}\omega-\dfrac{g}{L}\, \text{sen }\theta \end{array} \right.
\left\{ \begin{array}{rll} \theta ' &= &\omega \\ \omega' &=& -\dfrac{\gamma}{m}\omega-\dfrac{g}{L}\, \text{sen }\theta \end{array} \right.
γ=\gamma =
\gamma =

Constante de amortiguamiento

Péndulo amortiguado

{θ=ωω=γmωgLsen θ\left\{ \begin{array}{rll} \theta ' &= &\omega \\ \omega' &=& -\dfrac{\gamma}{m}\omega-\dfrac{g}{L}\, \text{sen }\theta \end{array} \right.
\left\{ \begin{array}{rll} \theta ' &= &\omega \\ \omega' &=& -\dfrac{\gamma}{m}\omega-\dfrac{g}{L}\, \text{sen }\theta \end{array} \right.

🌈 Ondas del péndulo

Referencias

Gracias por su atención...

😃

ABC of Mathematics:

An interactive experience

🔗 eBook para iPad

Complex Analysis:

A visual and interactive introduction

🔗 complex-analysis.com

Visita 🔗 jcponce.com 

para ver otros proyectos

Topología General

🔗 topologia-general.github.io

Simulación del péndulo

By Juan Carlos Ponce Campuzano

Simulación del péndulo

Presentación para la V Jornada de GeoGebra de la Comunitat Valenciana.

  • 1,606