Articulo de referencia

El péndulo de Kapitza

Dibujo que muestra cómo se puede construir un péndulo de Kapitza: un motor hace girar una manivela a alta velocidad, la manivela hace vibrar un brazo de palanca hacia arriba y h...

Dibujo que muestra cómo se puede construir un péndulo de Kapitza: un motor hace girar una manivela a alta velocidad, la manivela hace vibrar un brazo de palanca hacia arriba y hacia abajo, al cual está unido el péndulo mediante un pivote.

El péndulo de Kapitza es un péndulo rígido cuyo punto de pivote vibra verticalmente, hacia arriba y hacia abajo. Recibe su nombre del físico ruso y premio Nobel Pyotr Kapitza , quien en 1951 desarrolló una teoría que explica con éxito algunas de sus propiedades inusuales. [ 1 ] La característica única del péndulo de Kapitza es que la suspensión vibrante puede hacer que se equilibre de forma estable en una posición invertida , con la masa por encima del punto de suspensión. En el péndulo convencional con suspensión fija, la única posición de equilibrio estable es con la masa colgando por debajo del punto de suspensión; la posición invertida es un punto de equilibrio inestable , y la más mínima perturbación desplaza el péndulo fuera del equilibrio. En la teoría de control no lineal, el péndulo de Kapitza se utiliza como ejemplo de oscilador paramétrico que demuestra el concepto de "estabilización dinámica".

El péndulo fue descrito por primera vez por Andrew Stephenson en 1908, quien descubrió que la posición vertical superior del péndulo podía ser estable cuando la frecuencia de excitación era alta. [ 2 ] Sin embargo, hasta la década de 1950 no existía una explicación para este fenómeno tan inusual y contraintuitivo. Pyotr Kapitza fue el primero en analizarlo en 1951. [ 1 ] Llevó a cabo varios estudios experimentales y también proporcionó una comprensión analítica de las razones de la estabilidad al dividir el movimiento en variables "rápidas" y "lentas" e introducir un potencial efectivo. Este trabajo innovador creó una nueva disciplina en física: la mecánica vibracional . El método de Kapitza se utiliza para describir procesos periódicos en física atómica , física de plasmas y física cibernética . El potencial efectivo que describe el componente "lento" del movimiento se describe en el volumen "Mecánica" (§30) del Curso de Física Teórica de Landau . [ 3 ]

Otra característica interesante del sistema de péndulo de Kapitza es que la posición de equilibrio inferior, con el péndulo colgando por debajo del pivote, ya no es estable. Cualquier pequeña desviación de la vertical aumenta su amplitud con el tiempo. [ 4 ] También puede producirse resonancia paramétrica en esta posición, y pueden darse regímenes caóticos en el sistema cuando hay atractores extraños en la sección de Poincaré . [ 5 ]

Notación

El plan del péndulo de Kapitza

Denotemos el eje vertical comoy{\displaystyle y}y el eje horizontal comoincógnita{\displaystyle x}de modo que el movimiento del péndulo ocurra en el (incógnita{\displaystyle x}-y{\displaystyle y}) plano. Se utilizará la siguiente notación.

  • gramo{\displaystyle g}— aceleración de caída libre,
  • l{\displaystyle l}— longitud del péndulo rígido y ligero,
  • metro{\displaystyle m}— masa del péndulo (idealizado, matemático),
  • ν{\displaystyle \nu }—frecuencia de las oscilaciones verticales de la suspensión,
  • a{\displaystyle a}— amplitud de las oscilaciones de la suspensión,
  • ω0=gramo/l{\displaystyle \omega _{0}={\sqrt {g/l}}}— frecuencia propia del péndulo matemático.

Denotando el ángulo entre el péndulo y la dirección hacia abajo comoφ{\displaystyle \varphi }La dependencia temporal de la posición del péndulo se escribe como

{incógnita=lpecadoφy=lporqueφaporqueνt{\displaystyle {\begin{casos}x&=l\sin \varphi \\y&=-l\cos \varphi -a\cos \nu t\end{casos}}}

Energía

La energía potencial del péndulo se debe a la gravedad y se define, en términos de la posición vertical, como

miPAGOT=metrogramoy=metrogramo(lporqueφ+aporqueνt).{\displaystyle E_{\mathrm {POT} }=mgy=-mg(l\cos \varphi +a\cos \nu t).\,}

La energía cinética incluye, además del término estándarmiKInorte=metrol2φ˙2/2{\displaystyle E_{\mathrm {KIN} }=ml^{2}{\dot {\varphi }}^{2}/2}para la velocidad de un péndulo matemático, una contribución debida a las vibraciones de la suspensión y un término cruzado

miKInorte=metrol22φ˙2+metroalν pecado(νt)pecado(φ) φ˙+metroa2ν22pecado2(νt).{\displaystyle E_{\mathrm {KIN} }={\frac {ml^{2}}{2}}{\dot {\varphi }}^{2}+mal\nu ~\sin(\nu t)\sin(\varphi )~{\dot {\varphi }}+{\frac {ma^{2}\nu ^{2}}{2}}\sin ^{2}(\nu t)\;.}

La energía total viene dada por la suma de las energías cinética y potencial.mi=miKInorte+miPAGOT{\displaystyle E=E_{\mathrm {KIN} }+E_{\mathrm {POT} }}y el lagrangiano por su diferenciaL=miKInortemiPAGOT{\displaystyle L=E_{\mathrm {KIN} }-E_{\mathrm {POT} }}.

La energía total se conserva en un péndulo matemático, por lo que el tiempot{\displaystyle t}dependencia del potencialmiPAGOT{\displaystyle E_{\mathrm {POT} }}y cinéticamiKInorte{\displaystyle E_{\mathrm {KIN} }}La energía es simétrica con respecto a la línea horizontal. Según el teorema del virial, las energías cinética y potencial medias en un oscilador armónico son iguales. Esto significa que la línea de simetría corresponde a la mitad de la energía total.

En el caso de una suspensión vibratoria, el sistema ya no es cerrado y la energía total ya no se conserva. La energía cinética es más sensible a la vibración que la potencial. La energía potencial está limitada tanto por debajo como por encima.metrogramo(l+a)miPAGOTmetrogramo(l+a){\displaystyle -mg(l+a)\leq E_{\mathrm {POT} }\leq mg(l+a)}mientras que la energía cinética está limitada solo desde abajomiKInorte0{\displaystyle E_{\mathrm {KIN} }\geq 0}Para vibraciones de alta frecuenciaν{\displaystyle \nu }La energía cinética puede ser grande en comparación con la energía potencial.

Ecuaciones de movimiento

El movimiento del péndulo satisface las ecuaciones de Euler-Lagrange . La dependencia de la faseφ{\displaystyle \varphi }del péndulo en su posición satisface la ecuación: [ 6 ]

ddtLφ˙=Lφ,{\displaystyle {\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {\varphi }}}}={\frac {\partial L}{\partial \varphi }},}

donde el LagrangianoL{\displaystyle L}lecturas

L=metrol22φ˙2+metrol(gramo+a ν2porqueνt)porqueφ,{\displaystyle L={\frac {ml^{2}}{2}}{\dot {\varphi }}^{2}+ml(g+a~\nu ^{2}\cos \nu t)\cos \varphi ,}

hasta términos de derivadas temporales totales irrelevantes . La ecuación diferencial

φ¨=(gramo+a ν2porqueνt)pecadoφl,{\displaystyle {\ddot {\varphi }}=-(g+a~\nu ^{2}\cos \nu t){\frac {\sin \varphi }{l}},}

que describe el movimiento del péndulo es no lineal debido a lapecadoφ{\displaystyle \sin \varphi }factor.

Posiciones de equilibrio

El modelo de péndulo de Kapitza es más general que el péndulo simple . El modelo de Kapitza se reduce a este último en el límite.a=0{\displaystyle a=0}En ese límite, la punta del péndulo describe un círculo:incógnita2+y2=l2=constante{\displaystyle x^{2}+y^{2}=l^{2}={\text{constante}}}Si la energía en el momento inicial es mayor que el máximo de la energía potencialmi>metrogramol{\displaystyle E>mgl}entonces la trayectoria será cerrada y cíclica. Si la energía inicial es menormi<metrogramol{\displaystyle E<mgl}Entonces el péndulo oscilará cerca del único punto estable.(incógnita,y)=(0,l){\displaystyle (x,y)=(0,-l)}.

Cuando la suspensión vibra con una amplitud pequeñaal{\displaystyle a\ll l}y con frecuenciaνω0{\displaystyle \nu \gg \omega _ {0}}mucho más alta que la frecuencia adecuadaω0{\displaystyle \omega _{0}}, el ánguloφ{\displaystyle \varphi }puede considerarse como una superposiciónφ=φ0+ξ{\displaystyle \varphi =\varphi _ {0}+\xi}de un componente "lento"φ0{\displaystyle \varphi _{0}}y una oscilación rápidaξ{\displaystyle \xi }con pequeña amplitud debido a las pequeñas pero rápidas vibraciones de la suspensión. Técnicamente, realizamos una expansión perturbativa en las " constantes de acoplamiento ".(a/l),(ω0/ν)1{\displaystyle (a/l),(\omega _ {0}/\nu )\ll 1}mientras se trata la proporción(a/l)(ν/ω0){\displaystyle (a/l)(\nu /\omega _ {0})}como fijo. El tratamiento perturbativo se vuelve exacto en el límite de doble escala.a0,ν{\displaystyle a\to 0,\nu \to \infty }. Más precisamente, la rápida oscilaciónξ{\displaystyle \xi }se define como

ξ=alpecadoφ0 porqueνt.{\displaystyle \xi ={\frac {a}{l}}\sin \varphi _{0}~\cos \nu t.\,}

La ecuación de movimiento para el componente "lento".φ0{\displaystyle \varphi _{0}}se convierte

φ¨0=φ¨ξ¨=(gramo+a ν2porqueνt)pecadoφlal(φ¨0porqueφ0 porqueνtφ˙02pecadoφ0 porqueνt2νφ˙0porqueφ0 pecadoνtν2pecadoφ0 porqueνt)=gramolpecadoφ0(gramo+a ν2porqueνt)1l(ξporqueφ0+O(ξ2))al(φ¨0porqueφ0 porqueνtφ˙02pecadoφ0 porqueνt2νφ˙0porqueφ0 pecadoνt).{\displaystyle {\begin{aligned}{\ddot {\varphi }}_{0}={\ddot {\varphi }}-{\ddot {\xi }}&=-(g+a~\nu ^{2}\cos \nu t){\frac {\sin \varphi }{l}}\\&{}\quad {}-{\frac {a}{l}}\left({\ddot {\varphi }}_{0}\cos \varphi _{0}~\cos \nu t-{\dot {\varphi }}_{0}^{2}\sin \varphi _{0}~\cos \nu t-2\nu {\dot {\varphi }}_{0}\cos \varphi _{0}~\sin \nu t-\nu ^{2}\sin \varphi _{0}~\cos \nu t\right)\\[8pt]&=-{\frac {g}{l}}\sin \varphi _{0}-(g+a~\nu ^{2}\cos \nu t){\frac {1}{l}}\left(\xi \cos \varphi _{0}+O(\xi ^{2})\right)\\&{}\quad {}-{\frac {a}{l}}\left({\ddot {\varphi }}_{0}\cos \varphi _{0}~\cos \nu t-{\dot {\varphi }}_{0}^{2}\sin \varphi _{0}~\cos \nu t-2\nu {\dot {\varphi }}_{0}\cos \varphi _{0}~\sin \nu t\right).\end{aligned}}}

Promedio temporal sobre el rápidoν{\displaystyle \nu }-la oscilación cede al orden principal

φ¨0=gramolpecadoφ012(aνl)2pecadoφ0porqueφ0.{\displaystyle {\ddot {\varphi }}_{0}=-{\frac {g}{l}}\sin \varphi _{0}-{\frac {1}{2}}\left({\frac {a\nu }{l}}\right)^{2}\sin \varphi _{0}\cos \varphi _{0}.}

La ecuación de movimiento "lenta" se convierte en

metrol2φ¨0=VmiFFφ0,{\displaystyle ml^{2}{\ddot {\varphi }}_{0}=-{\frac {\partial V_{\mathrm {eff} }}{\partial \varphi _{0}}}\;,}

al introducir un potencial efectivo

VmiFF=metrogramolporqueφ0+metro(aν2pecadoφ0)2.{\displaystyle V_{\mathrm {eff} }=-mgl\cos \varphi _{0}+m\left({\frac {a\nu }{2}}\sin \varphi _{0}\right)^{2}.}

Resulta que [ 1 ] el potencial efectivoVmiFF{\displaystyle V_{\mathrm {eff} }}tiene dos mínimos si(aν)2>2gramol{\displaystyle (a\nu )^{2}>2gl}, o equivalentemente,(a/l)(ν/ω0)>2{\displaystyle (a/l)(\nu /\omega _{0})>{\sqrt {2}}}El primer mínimo está en la misma posición.(incógnita,y)=(0,l){\displaystyle (x,y)=(0,-l)}como el péndulo matemático y el otro mínimo está en la posición vertical superior(incógnita,y)=(0,l){\displaystyle (x,y)=(0,l)}Como resultado, la posición vertical superior, que es inestable en un péndulo matemático, puede volverse estable en el péndulo de Kapitza.

Potencial efectivo
Potencial efectivo

Bifurcación y caos

Para una frecuencia de forzamiento fijaν{\displaystyle \nu }y aumentando la amplitud de la fuerzaa{\displaystyle a}, el punto invertidoφ=π{\displaystyle \varphi =\pi }se convierte en un equilibrio estable, y luego volvería a ser inestable. El punto experimenta una bifurcación de Hopf : el equilibrio estable se expande en una oscilación estable alrededor deφ=π{\displaystyle \varphi =\pi }, a frecuenciaν/2{\displaystyle \nu /2}y amplitudaado{\displaystyle \propto {\sqrt {a-a_{c}}}}dóndeado{\displaystyle a_{c}}es la amplitud crítica. [ 5 ]

Para amplitudes aún mayores, la oscilación estable se vuelve inestable y el péndulo comienza a girar. Con amplitudes aún mayores, el modo de rotación también se destruye y aparece un atractor extraño por cascada de duplicación de período . [ 7 ] [ 8 ]

Soluciones rotativas

Las soluciones rotatorias del péndulo de Kapitza se producen cuando el péndulo gira alrededor del punto de pivote con la misma frecuencia con la que se impulsa dicho punto. Existen dos soluciones rotatorias, una para cada dirección de rotación. Cambiamos al sistema de referencia giratorio utilizandoφφ±νt{\displaystyle \varphi \rightarrow \varphi ^{\prime }\pm \nu t}y la ecuación paraφ{\displaystyle \varphi }se convierte en:

φ¨=1l[12aν2pecado(φ)+gramopecado(φ±νt)+12aν2pecado(φ±2νt)].{\displaystyle {\ddot {\varphi }}^{\prime }=-{\frac {1}{l}}\left[{\frac {1}{2}}a\nu ^{2}\sin(\varphi ^{\prime })+g\sin(\varphi ^{\prime }\pm \nu t)+{\frac {1}{2}}a\nu ^{2}\sin(\varphi ^{\prime }\pm 2\nu t)\right]\;.}

Considerando nuevamente el límite en el queν{\displaystyle \nu }es mucho más alta que la frecuencia adecuadaω0{\displaystyle \omega _{0}}, encontramos que el rápido-ν{\displaystyle \nu }lento-φ0{\displaystyle \varphi _{0}^{\prime }}El límite conduce a la ecuación:

φ¨0=12laν2pecadoφ0.{\displaystyle {\ddot {\varphi }}_{0}^{\prime }=-{\frac {1}{2l}}a\nu ^{2}\sin \varphi _{0}^{\prime }\;.}

El potencial efectivo es simplemente el de una ecuación de péndulo simple. Existe un equilibrio estable enφ0=0{\displaystyle \varphi _{0}^{\prime }=0}y un equilibrio inestable enφ0=π{\displaystyle \varphi _{0}^{\prime }=\pi }.

Retrato de fase

Se pueden obtener retratos de fase interesantes en regímenes que no son accesibles dentro de las descripciones analíticas, por ejemplo en el caso de una gran amplitud de la suspensión.al{\displaystyle a\approx l}. [ 9 ] [ 10 ] Aumentar la amplitud de las oscilaciones impulsoras a la mitad de la longitud del pénduloa=l/2{\displaystyle a=l/2}conduce al retrato de fase que se muestra en la figura.

Aumento adicional de la amplitud aal{\displaystyle a\approx l}conduce al llenado completo de los puntos internos del espacio de fases: si antes algunos puntos del espacio de fases no eran accesibles, ahora el sistema puede alcanzar cualquiera de los puntos internos. Esta situación también se cumple para valores mayores dea{\displaystyle a}.

Datos interesantes

  • Kapitza señaló que un reloj de péndulo con suspensión de péndulo vibratorio siempre va más rápido que un reloj con suspensión fija. [ 11 ]
  • La marcha se define por un movimiento de péndulo invertido en el que el cuerpo se eleva sobre la extremidad o extremidades rígidas con cada paso. Una mayor estabilidad durante la marcha podría estar relacionada con la estabilidad del péndulo de Kapitza. Esto se aplica independientemente del número de extremidades, incluso en artrópodos con seis, ocho o más extremidades. [ 12 ]

Literatura

  1. 1 2 3 Kapitza PL (1951). "Estabilidad dinámica de un péndulo cuando su punto de suspensión vibra". Soviet Phys. JETP . 21 : 588– 597.; Kapitza PL (1951). «Péndulo con suspensión vibratoria» . Usp. Fiz. Nauk . 44 : 7– 15. doi : 10.3367/UFNr.0044.195105b.0007 .
  2. Stephenson Andrew (1908). "XX. Sobre la estabilidad inducida" . Philosophical Magazine . 6. 15 (86): 233– 236. doi : 10.1080/14786440809463763 .
  3. LD Landau , EM Lifshitz (1960). Mecánica . Vol. 1 (1.ª ed.). Pergamon Press. ASIN B0006AWV88 .   
  4. Бутиков Е. И. «Маятник с осциллирующим подвесом (к 60-летию маятника Капицы»), учебное пособие .
  5. 1 2 Blackburn, James A.; Smith, HJT; Gro/nbech-Jensen, N. (octubre de 1992). "Estabilidad y bifurcaciones de Hopf en un péndulo invertido" . American Journal of Physics . 60 (10): 903– 908. Bibcode : 1992AmJPh..60..903B . doi : 10.1119/1.17011 . ISSN 0002-9505 . 
  6. VP Krainov (2002). Métodos matemáticos selectos en física teórica (1.ª ed.). Taylor & Francis . ISBN  978-0-415-27234-6.
  7. McLaughlin, John B. (febrero de 1981). "Bifurcaciones de duplicación de período y movimiento caótico para un péndulo paramétricamente forzado" . Journal of Statistical Physics . 24 (2): 375– 388. Bibcode : 1981JSP....24..375M . doi : 10.1007/BF01013307 . ISSN 0022-4715 . S2CID 55706875 .  
  8. Koch, BP; Leven, RW; Pompe, B.; Wilke, C. (1983-07-04). "Evidencia experimental del comportamiento caótico de un péndulo paramétricamente forzado" . Physics Letters A. 96 ( 5): 219– 224. Bibcode : 1983PhLA...96..219K . doi : 10.1016/0375-9601(83)90336-5 . ISSN 0375-9601 . 
  9. GE Astrakharchik, NA Astrakharchik «Estudio numérico del péndulo Kapitza» arXiv : 1103.5981 (2011)
  10. El movimiento temporal del péndulo de Kapitza se puede modelar en applets de Java en línea en los siguientes sitios:
    • "Simulación física de péndulo vibrante invertido" . Archivado del original el 1 de octubre de 2011. Consultado el 8 de abril de 2011 .
    • Se pueden utilizar parámetros arbitrarios del sistema y se pueden insertar manualmente.
  11. ^ Butikov, Eugene I. "Péndulo de Kapitza: una explicación simple físicamente transparente" (PDF) . pag. 8 . Consultado el 1 de septiembre de 2020 . 
  12. Quintanilla, José; Pérez, Moisés; Balderas, Rodolfo; González, Alejandro; Cárdenas, Antonio; Maya, Mauro; Piovesan, Davide (mayo de 2021). "Estabilización inercial de la postura erguida al caminar". 2021 Décima Conferencia Internacional IEEE/EMBS sobre Ingeniería Neural (NER) . págs. 849– 852. doi : 10.1109/NER49283.2021.9441114 . ISBN  978-1-7281-4337-8. S2CID 235307555 .