Physics · Libro 5 · Bachelor Year 3

Física universitaria — tercer año

Física universitaria — tercer año · Bachelor Year 3

1Mecánica lagrangiana

Inténtese escribir la ley de Newton para un péndulo doble: dos fuerzas de tensión, ninguna conocida de antemano, ambas cambiando de dirección a cada instante, y cuatro ecuaciones escalares que desenredar solo para averiguar cómo evolucionan dos ángulos. Hágase lo mismo ahora con una cuenta ensartada en un aro que gira, con una cadena de resortes acoplados, con el brazo de un robot. Las fuerzas que mantienen unido un sistema — varillas, carriles, bisagras — no realizan trabajo y, sin embargo, el método de Newton nos obliga a arrastrarlas por todo el cálculo. Este capítulo presenta la reformulación que Lagrange publicó en 1788: describir el sistema mediante las pocas coordenadas que realmente pueden variar, escribir una única función L=EkEpL = E_k - E_p y aplicar una sola receta que produce las ecuaciones del movimiento en cualesquiera coordenadas, con todas las varillas y todos los carriles ya eliminados. Mejor aún: la nueva formulación hace visible lo que la de Newton oculta, pues cada simetría de LL proporciona una magnitud conservada — la cantidad de movimiento nace de la uniformidad del espacio, la energía de la uniformidad del tiempo —, un resultado de Emmy Noether que se ha convertido en el principio organizador de la física mucho más allá de la mecánica.

1.1 De las fuerzas a las coordenadas

Definición 1.1 (Ligaduras, grados de libertad, coordenadas generalizadas)

Una ligadura es una condición geométrica impuesta sobre las posiciones de un sistema — una cuenta permanece en su alambre, la varilla de un péndulo conserva una longitud fija, dos ruedas de un mismo eje giran solidariamente. Una ligadura que puede expresarse como una ecuación f(r1,,rN,t)=0f(\vect r_1, \dots, \vect r_N, t) = 0 entre las coordenadas (y eventualmente el tiempo) se llama holónoma. El número de maneras independientes en que la configuración todavía puede variar es el número de grados de libertad nn; todo conjunto de nn magnitudes independientes q1,,qnq_1, \dots, q_n que fije por completo la configuración es un conjunto de coordenadas generalizadas — ángulos, longitudes o cualquier mezcla cómoda. Sus derivadas temporales q˙1,,q˙n\dot q_1, \dots, \dot q_n son las velocidades generalizadas.

Ejemplo 1.2 (Recuento de grados de libertad)

Un punto sobre una mesa: n=2n = 2. Un péndulo plano de longitud fija: un ángulo, n=1n = 1. Un péndulo doble: dos ángulos, n=1+1=2n = 1 + 1 = 2. Una cuenta sobre un aro rígido: un ángulo, n=1n = 1 — incluso si se obliga al aro a girar, pues la rotación impuesta no añade libertad alguna. Un sólido rígido libre en el espacio: tres coordenadas de su centro más tres ángulos, n=6n = 6. Un gas de NN moléculas libres (puntuales): n=3Nn = 3N. Cada ligadura holónoma suprime un grado de libertad: dos puntos unidos por una varilla tienen 3+31=53 + 3 - 1 = 5.

Coordenadas generalizadas: cada sistema queda descrito por los ángulos que realmente pueden variar, no por las coordenadas cartesianas de sus masas. Las tensiones de las varillas y la reacción normal del aro no aparecen jamás.
Coordenadas generalizadas: cada sistema queda descrito por los ángulos que realmente pueden variar, no por las coordenadas cartesianas de sus masas. Las tensiones de las varillas y la reacción normal del aro no aparecen jamás.

1.2 El principio de mínima acción

Definición 1.3 (Lagrangiano y acción)

El lagrangiano de un sistema mecánico cuyas fuerzas derivan de una energía potencial EpE_p es la función de las coordenadas, de las velocidades y, eventualmente, del tiempo

L(q,q˙,t)=EkEp,L(q, \dot q, t) = E_k - E_p ,

energía cinética menos energía potencial, ambas expresadas en las coordenadas generalizadas. La acción de un movimiento concebible q(t)q(t) entre extremos fijos q(t1)q(t_1) y q(t2)q(t_2) es el número

S[q]=t1t2L(q(t),q˙(t),t) ⁣dt.S[q] = \int_{t_1}^{t_2} L\big(q(t), \dot q(t), t\big)\,\dd t .

SS es un funcional: se traga un camino entero y devuelve un solo número, en julios-segundo — las unidades de la constante de Planck.

Teorema 1.4 (Principio de Hamilton y ecuaciones de Euler–Lagrange)

Entre todos los movimientos concebibles que unen los mismos dos extremos en el mismo tiempo, el movimiento real es aquel que hace la acción estacionaria (principio de Hamilton, o principio de mínima acción). Equivalentemente, el movimiento obedece a las ecuaciones de Euler–Lagrange

 ⁣d ⁣dtLq˙iLqi=0,i=1,,n:\frac{\dd}{\dd t}\frac{\partial L}{\partial\dot q_i} - \frac{\partial L}{\partial q_i} = 0 , \qquad i = 1, \dots, n :

una ecuación de segundo orden por grado de libertad, en las coordenadas que se hayan elegido.

Demostración. Deformemos el camino: qi(t)qi(t)+δqi(t)q_i(t) \to q_i(t) + \delta q_i(t) con δqi(t1)=δqi(t2)=0\delta q_i(t_1) = \delta q_i(t_2) = 0. A primer orden,

δS=t1t2i(Lqiδqi+Lq˙iδq˙i) ⁣dt=t1t2i(Lqi ⁣d ⁣dtLq˙i)δqi ⁣dt,\delta S = \int_{t_1}^{t_2}\sum_i\Big( \frac{\partial L}{\partial q_i}\,\delta q_i + \frac{\partial L}{\partial\dot q_i}\,\delta\dot q_i\Big)\dd t = \int_{t_1}^{t_2}\sum_i\Big( \frac{\partial L}{\partial q_i} - \frac{\dd}{\dd t}\frac{\partial L}{\partial\dot q_i}\Big)\delta q_i\,\dd t ,

tras integrar por partes el segundo término (δq˙i= ⁣d(δqi)/ ⁣dt\delta\dot q_i = \dd(\delta q_i)/\dd t) y suprimir el término de borde, que se anula en los extremos fijos. Si δS=0\delta S = 0 para toda deformación, el corchete debe anularse en cada instante y para cada ii: si fuera positivo en algún punto, una protuberancia δqi\delta q_i concentrada allí daría δS0\delta S \neq 0. El recíproco se lee en la misma línea.

Principio de Hamilton: entre todos los caminos con los mismos extremos y la misma duración, el real hace estacionaria S = ∈t L\, t — las deformaciones q de primer orden cambian S solo a segundo orden.
Principio de Hamilton: entre todos los caminos con los mismos extremos y la misma duración, el real hace estacionaria S=L ⁣dtS = \int L\,\dd t — las deformaciones δq\delta q de primer orden cambian SS solo a segundo orden.

Proposición 1.5 (Newton recuperado)

Para una partícula en coordenadas cartesianas, L=12m(x˙2+y˙2+z˙2)Ep(x,y,z)L = \tfrac12 m(\dot x^2 + \dot y^2 + \dot z^2) - E_p(x,y,z), y las ecuaciones de Euler–Lagrange se escriben mx¨=Ep/xm\ddot x = -\partial E_p/\partial x y análogamente para yy, zz: exactamente ma=Fm\vect a = \vect F. Las mecánicas lagrangiana y newtoniana coinciden dondequiera que ambas se apliquen; la nueva forma sobrevive sin más a un cambio de coordenadas, cosa que las ecuaciones en componentes de Newton no hacen.

Demostración. L/x˙=mx˙\partial L/\partial\dot x = m\dot x, L/x=Ep/x\partial L/\partial x = -\partial E_p/\partial x; la ecuación de Euler–Lagrange es  ⁣d(mx˙)/ ⁣dt+Ep/x=0\dd(m\dot x)/\dd t + \partial E_p/\partial x = 0.

Observación 1.6 (Por qué desaparecen las fuerzas de ligadura)

La tensión de una varilla o la reacción normal de un carril actúan perpendicularmente a todo desplazamiento que la ligadura permite: no realizan trabajo en ningún movimiento compatible con ella. Como la acción se construye con energías evaluadas solo sobre tales movimientos, esas fuerzas no entran nunca en LL; ese es el milagro práctico del método. La justificación general (el principio de los trabajos virtuales de d’Alembert) se admite aquí; para todos los sistemas de este libro la receta siguiente puede comprobarse directamente frente a Newton, como empezó a hacer la Proposición 1.5. Si lo que se quiere es la propia fuerza de ligadura — ¿se romperá la varilla? —, se vuelve a Newton solo para esa fuerza, con el movimiento ya conocido.

Método 1.7 (La receta lagrangiana)

(1) Contar los grados de libertad y elegir las coordenadas qiq_i — ángulos para las rotaciones, abscisas a lo largo de los carriles. (2) Expresar las posiciones r(qi,t)\vect r(q_i, t), derivarlas para obtener las velocidades y escribir EkE_k; escribir EpE_p. (3) L=EkEpL = E_k - E_p, suprimiendo cualquier constante aditiva. (4) Una ecuación de Euler–Lagrange por coordenada. (5) Antes de resolver, cosechar las magnitudes conservadas: las coordenadas cíclicas (Definición 1.10) y la función energía (Proposición 1.14). (6) Comprobar los límites: ángulos pequeños, rotación desconectada, casos particulares conocidos.

Ejemplo 1.8 (El péndulo, en tres líneas)

Una sola coordenada θ\theta; v=θ˙eθ\vect v = \ell\dot\theta\,\vect e_\theta, de modo que Ek=12m2θ˙2E_k = \tfrac12 m\ell^2\dot\theta^2 y Ep=mgcosθE_p = -mg\ell\cos\theta. Entonces L=12m2θ˙2+mgcosθL = \tfrac12 m\ell^2\dot\theta^2 + mg\ell\cos\theta, y  ⁣d(m2θ˙)/ ⁣dt=mgsinθ\dd(m\ell^2\dot\theta)/\dd t = -mg\ell\sin\theta:

θ¨=gsinθ,\ddot\theta = -\frac{g}{\ell}\sin\theta ,

la ecuación que el volumen del Año 1 obtuvo a partir del momento del peso — sin mencionar jamás la tensión.

Ejemplo 1.9 (Cuenta sobre un aro que gira)

Una cuenta de masa mm desliza sobre un aro circular vertical de radio RR obligado a girar en torno a su diámetro vertical con ω\omega constante (Ejemplo 1.2). Una sola coordenada: el ángulo polar θ\theta medido desde el punto más bajo. La velocidad de la cuenta tiene una componente Rθ˙R\dot\theta a lo largo del aro y otra ωRsinθ\omega R\sin\theta debida a la rotación impuesta, perpendicular a aquella:

L=12mR2θ˙2+12mω2R2sin2θ+mgRcosθ,mR2θ¨=mω2R2sinθcosθmgRsinθ.L = \tfrac12 mR^2\dot\theta^2 + \tfrac12 m\omega^2R^2\sin^2\theta + mgR\cos\theta , \qquad mR^2\ddot\theta = m\omega^2R^2\sin\theta\cos\theta - mgR\sin\theta .

Equilibrios donde θ¨=0\ddot\theta = 0: el punto más bajo θ=0\theta = 0 (y el más alto, siempre inestable) y, cuando ω2>g/R\omega^2 > g/R, una nueva pareja cosθeq=g/ω2R\cos\theta_{\text{eq}} = g/\omega^2R. Escribir la ecuación en la forma mR2θ¨= ⁣dUef/ ⁣dθmR^2\ddot\theta = -\dd U_{\text{ef}}/\dd\theta con la energía potencial efectiva

Uef(θ)=mgRcosθ12mω2R2sin2θU_{\text{ef}}(\theta) = -mgR\cos\theta - \tfrac12 m\omega^2R^2\sin^2\theta

hace visible la geometría: por debajo del régimen crítico ωc=g/R\omega_{\text{c}} = \sqrt{g/R} el fondo es un pozo; por encima, el fondo se convierte en una cima y la cuenta se instala en la ladera, tanto más arriba cuanto más deprisa gira el aro — el principio del regulador centrífugo que gobernaba las máquinas de vapor.

El potencial efectivo de la cuenta sobre el aro que gira. Cuando  supera √g/R, el pozo único del fondo se desdobla en dos pozos simétricos: el ángulo de equilibrio crece con la velocidad de giro.
El potencial efectivo de la cuenta sobre el aro que gira. Cuando ω\omega supera g/R\sqrt{g/R}, el pozo único del fondo se desdobla en dos pozos simétricos: el ángulo de equilibrio crece con la velocidad de giro.

1.3 Simetrías y leyes de conservación

Definición 1.10 (Momento generalizado, coordenada cíclica)

El momento generalizado conjugado de la coordenada qiq_i es

pi=Lq˙i.p_i = \frac{\partial L}{\partial\dot q_i} .

Para una coordenada cartesiana es la cantidad de movimiento ordinaria mx˙m\dot x; para un ángulo es un momento angular. Una coordenada que no aparece en LL (aunque sí su velocidad) se llama cíclica.

Proposición 1.11 (Las coordenadas cíclicas dan leyes de conservación)

Si qiq_i es cíclica, su momento conjugado se conserva: L/qi=0\partial L/\partial q_i = 0 implica  ⁣dpi/ ⁣dt=0\dd p_i/\dd t = 0. Elegir las coordenadas de modo que sean cíclicas cuantas más mejor es el paso individual más eficaz para resolver un problema de mecánica.

Demostración. Inmediato a partir de la ecuación de Euler–Lagrange para qiq_i.

Ejemplo 1.12 (Fuerza central, resuelta por inspección)

Una partícula en un potencial central Ep(r)E_p(r), en coordenadas polares en su plano de movimiento:

L=12m(r˙2+r2φ˙2)Ep(r).L = \tfrac12 m(\dot r^2 + r^2\dot\varphi^2) - E_p(r) .

φ\varphi es cíclica, luego pφ=mr2φ˙p_\varphi = mr^2\dot\varphi — el momento angular — se conserva: la ley de las áreas de Kepler, que el volumen del Año 1 dedujo de la ecuación de momentos, aparece aquí antes de resolver ecuación alguna. La ecuación radial restante es mr¨=mrφ˙2Ep(r)m\ddot r = mr\dot\varphi^2 - E_p'(r), es decir, el movimiento unidimensional en el potencial efectivo Ep(r)+pφ2/2mr2E_p(r) + p_\varphi^2/2mr^2 de aquel volumen.

Teorema 1.13 (Teorema de Noether)

A toda simetría continua del lagrangiano le corresponde una magnitud conservada. Con precisión: si el desplazamiento qiqi+εKi(q)q_i \to q_i + \varepsilon\,K_i(q) deja LL invariante a primer orden en ε\varepsilon para todos los movimientos, entonces

Q=ipiKi(q)Q = \sum_i p_i\,K_i(q)

es constante a lo largo de todo movimiento real. La uniformidad del espacio (invariancia por traslación) da la cantidad de movimiento; la isotropía del espacio (invariancia por rotación) da el momento angular; la uniformidad del tiempo da la energía (Proposición 1.14).

Demostración. La invariancia a primer orden significa que 0=δL=i(qiLKi+q˙iLK˙i)ε0 = \delta L = \sum_i\big(\partial_{q_i}L\,K_i + \partial_{\dot q_i}L\,\dot K_i\big)\varepsilon. Sobre un movimiento real, qiL=p˙i\partial_{q_i}L = \dot p_i por Euler–Lagrange, así que el corchete es i(p˙iKi+piK˙i)= ⁣dQ/ ⁣dt\sum_i(\dot p_iK_i + p_i\dot K_i) = \dd Q/\dd t. Una coordenada cíclica es el caso particular Ki=δijK_i = \delta_{ij}. El caso del desplazamiento temporal exige el cálculo aparte de la Proposición 1.14; el teorema completo, para transformaciones que cambian también tt o alteran LL en una derivada total, se demuestra en los cursos de mecánica analítica y se admite aquí en esa generalidad.

Proposición 1.14 (La función energía)

A lo largo de todo movimiento, la función energía

h=iq˙iLq˙iLcumple ⁣dh ⁣dt=Lt.h = \sum_i \dot q_i\,\frac{\partial L}{\partial\dot q_i} - L \qquad\text{cumple}\qquad \frac{\dd h}{\dd t} = -\frac{\partial L}{\partial t} .

Si LL no depende explícitamente del tiempo, hh se conserva. Si además las relaciones r(q)\vect r(q) entre posiciones y coordenadas no involucran el tiempo — ninguna rotación impuesta, ningún soporte móvil —, entonces EkE_k es una forma cuadrática en las q˙i\dot q_i y h=Ek+Eph = E_k + E_p: la energía mecánica. Con una ligadura dependiente del tiempo, hh sigue conservándose cuando L/t=0\partial L/\partial t = 0, pero no es la energía: el motor que impone la ligadura intercambia trabajo con el sistema.

Demostración.  ⁣dh/ ⁣dt=i(q¨ipi+q˙ip˙i)i(qiLq˙i+q˙iLq¨i)tL\dd h/\dd t = \sum_i(\ddot q_ip_i + \dot q_i\dot p_i) - \sum_i(\partial_{q_i}L\,\dot q_i + \partial_{\dot q_i}L\,\ddot q_i) - \partial_tL. Los términos en q¨i\ddot q_i se cancelan; las ecuaciones de Euler–Lagrange convierten p˙i\dot p_i en qiL\partial_{q_i}L, lo que cancela la pareja siguiente; queda tL-\partial_tL. Si Ek=12ajk(q)q˙jq˙kE_k = \tfrac12\sum a_{jk}(q)\dot q_j\dot q_k, la identidad de Euler para las formas cuadráticas da q˙iEk/q˙i=2Ek\sum\dot q_i\,\partial E_k/\partial\dot q_i = 2E_k, luego h=2Ek(EkEp)=Ek+Eph = 2E_k - (E_k - E_p) = E_k + E_p.

Ejemplo 1.15 (El aro que gira conserva hh, no EE)

Para la cuenta del Ejemplo 1.9, LL no contiene tt explícitamente, así que h=12mR2θ˙2+Uef(θ)h = \tfrac12 mR^2\dot\theta^2 + U_{\text{ef}} (\theta) se conserva — pero la energía mecánica E=12mR2θ˙2+12mω2R2sin2θmgRcosθE = \tfrac12 mR^2\dot\theta^2 + \tfrac12 m\omega^2R^2\sin^2\theta - mgR\cos\theta no: E=h+mω2R2sin2θE = h + m\omega^2R^2\sin^2\theta varía a medida que la cuenta desliza. La diferencia es el trabajo del motor que mantiene ω\omega constante mientras cambia la distancia de la cuenta al eje.

1.4 Partículas cargadas y pequeñas oscilaciones

Proposición 1.16 (Lagrangiano de una partícula cargada)

En un campo electromagnético descrito por los potenciales VV y A\vect A (con E=VtA\vect E = -\vect\nabla V - \partial_t\vect A y B=curlA\vect B = \operatorname{\vect{curl}}\vect A, como en el volumen del Año 2), el lagrangiano

L=12mv2qV+qvAL = \tfrac12 m\vect v^{\,2} - qV + q\,\vect v\cdot\vect A

produce, a través de las ecuaciones de Euler–Lagrange, exactamente la fuerza de Lorentz mv˙=q(E+vB)m\dot{\vect v} = q(\vect E + \vect v\wedge\vect B). La fuerza magnética, que no realiza trabajo y no deriva de ninguna energía potencial ordinaria, entra por un término lineal en la velocidad; el momento conjugado pasa a ser p=mv+qA\vect p = m\vect v + q\vect A, y ya no mvm\vect v solo.

Demostración. Para la componente xx: px=mx˙+qAxp_x = m\dot x + qA_x y xL=qxV+qvxA\partial_xL = -q\,\partial_xV + q\,\vect v\cdot\partial_x\vect A. La ecuación de Euler–Lagrange da mx¨=qxVq ⁣dAx/ ⁣dt+qvxAm\ddot x = -q\,\partial_xV - q\,\dd A_x/\dd t + q\,\vect v\cdot\partial_x\vect A. A lo largo del movimiento  ⁣dAx/ ⁣dt=tAx+(v)Ax\dd A_x/\dd t = \partial_tA_x + (\vect v\cdot\vect\nabla)A_x, luego mx¨=qEx+q[vxA(v)Ax]m\ddot x = qE_x + q\big[\vect v\cdot\partial_x\vect A - (\vect v\cdot \vect\nabla)A_x\big], y el corchete es la componente xx de v(A)=vB\vect v\wedge(\vect\nabla\wedge\vect A) = \vect v\wedge\vect B — desarróllense ambos para comprobarlo.

Proposición 1.17 (Pequeñas oscilaciones y modos normales)

Cerca de un equilibrio estable qeqq^{\text{eq}}, desarróllese LL a segundo orden en los desplazamientos ui=qiqiequ_i = q_i - q_i^{\text{eq}}: L12miju˙iu˙j12kijuiujL \approx \tfrac12\sum m_{ij}\dot u_i\dot u_j - \tfrac12\sum k_{ij}u_iu_j con matrices simétricas constantes. Las ecuaciones del movimiento son lineales y todo movimiento es superposición de modos normales: oscilaciones colectivas ui(t)=aicos(Ωt+ϕ)u_i(t) = a_i\cos(\Omega t + \phi) en las que todas las coordenadas vibran a una misma frecuencia común; las amplitudes y las frecuencias resuelven j(kijΩ2mij)aj=0\sum_j(k_{ij} - \Omega^2m_{ij})\,a_j = 0 — un problema matricial de valores propios, con nn modos para nn grados de libertad.

Demostración. Los términos lineales se anulan en un equilibrio; las ecuaciones de Euler–Lagrange del LL cuadrático dan jmiju¨j=jkijuj\sum_jm_{ij}\ddot u_j = -\sum_jk_{ij}u_j. Al insertar la oscilación de prueba se obtiene el sistema lineal enunciado, que tiene vector de amplitudes no nulo solo cuando det(kijΩ2mij)=0\det(k_{ij} - \Omega^2m_{ij}) = 0: nn valores de Ω2\Omega^2, todos positivos en un equilibrio estable. Que el movimiento general sea superposición de los modos es la diagonalización simultánea de dos formas cuadráticas, un resultado del álgebra lineal del volumen de matemáticas del Año 2; la completitud se admite.

Ejemplo 1.18 (Dos péndulos acoplados por un resorte)

Dos péndulos iguales (masa mm, longitud \ell) cuyas lentejas están unidas por un resorte de rigidez kk que queda relajado cuando ambos cuelgan verticales. Para ángulos pequeños,

L=12m2(θ˙12+θ˙22)12mg(θ12+θ22)12k2(θ2θ1)2.L = \tfrac12 m\ell^2(\dot\theta_1^2 + \dot\theta_2^2) - \tfrac12 mg\ell(\theta_1^2 + \theta_2^2) - \tfrac12 k\ell^2(\theta_2 - \theta_1)^2 .

La simetría sugiere las combinaciones s=θ1+θ2s = \theta_1 + \theta_2 y d=θ1θ2d = \theta_1 - \theta_2, que desacoplan las ecuaciones: el modo en fase (θ1=θ2\theta_1 = \theta_2, con el resorte inerte) a Ω1=g/\Omega_1 = \sqrt{g/\ell}, y el modo en oposición (θ1=θ2\theta_1 = -\theta_2, con el resorte estirado el doble) a Ω2=g/+2k/m\Omega_2 = \sqrt{g/\ell + 2k/m}. Láncese un solo péndulo — una mezcla a partes iguales de los dos modos — y la energía migrará por completo de un péndulo al otro y de vuelta con la frecuencia de batido (Ω2Ω1)/2π(\Omega_2 - \Omega_1)/2\pi: la demostración de los péndulos acoplados y el mecanismo que hay detrás de toda transferencia resonante de energía, de los circuitos sintonizados a las vibraciones moleculares.

Los dos modos normales de los péndulos acoplados. En fase, el resorte no se estira nunca y la frecuencia es la del péndulo libre; en oposición, cada lenteja siente el resorte duplicado.
Los dos modos normales de los péndulos acoplados. En fase, el resorte no se estira nunca y la frecuencia es la del péndulo libre; en oposición, cada lenteja siente el resorte duplicado.
Un péndulo doble trazado con una larga exposición: dos coordenadas, un lagrangiano y un movimiento que ninguna fórmula predice por mucho tiempo — la maquinaria de mínima acción de este capítulo escribe las ecuaciones; el caos mantiene humildes a sus soluciones.
Un péndulo doble trazado con una larga exposición: dos coordenadas, un lagrangiano y un movimiento que ninguna fórmula predice por mucho tiempo — la maquinaria de mínima acción de este capítulo escribe las ecuaciones; el caos mantiene humildes a sus soluciones.

1.5 Ejercicios

Ejercicio 1.1

Cuéntense los grados de libertad y propónganse coordenadas generalizadas: (a) una partícula en el interior de un cuenco fijo; (b) un cilindro que rueda sin deslizar por un plano inclinado fijo; (c) un péndulo doble cuyo pivote superior desliza por un carril horizontal; (d) una mancuerna (dos masas y una varilla rígida) en el espacio; (e) dos cuentas sobre el mismo alambre circular fijo. ¿Cuál de estas ligaduras relaciona velocidades en lugar de posiciones y por qué es, pese a ello, integrable a una ligadura holónoma?

Solución

Solución de Ejercicio 1.1.

(a) 2 (dos ángulos sobre la superficie del cuenco). (b) 1: la abscisa xx a lo largo de la pendiente, pues el ángulo de rotación queda ligado a ella por la rodadura, x=Rϕx = R\phi. (c) 3: XX, θ1\theta_1, θ2\theta_2. (d) 5: tres para el centro y dos para la dirección de la varilla. (e) 2: un ángulo cada una. La condición de rodadura es una relación entre velocidades, x˙=Rϕ˙\dot x = R\dot\phi; en este problema plano se integra de inmediato a x=Rϕ+constx = R\phi + \text{const}, luego es holónoma. (Para una bola que rueda sobre un plano no se integra, y la mecánica lagrangiana necesita una extensión).

Ejercicio 1.2

Para el péndulo plano (Ejemplo 1.8): (a) verifíquense las dimensiones de LL y de pθ=L/θ˙p_\theta = \partial L/\partial\dot\theta; (b) identifíquese físicamente pθp_\theta; (c) dedúzcase el periodo para ángulos pequeños; (d) calcúlese la acción SS de una oscilación pequeña completa de amplitud θ0\theta_0 — y explíquese la respuesta antes de calcular (¿cuál es el promedio temporal de EkEpE_k - E_p en una oscilación armónica?).

Solución

Solución de Ejercicio 1.2.

(a) [L]=J[L] = \mathrm{J}; [pθ]=kgm2/s[p_\theta] = \mathrm{kg}\,\mathrm{m}^{2}/\mathrm{s} — un momento angular, ya que θ\theta es adimensional. (b) pθ=m2θ˙p_\theta = m\ell^2\dot\theta es el momento angular de la lenteja respecto del pivote. (c) θ¨=(g/)θ\ddot\theta = -(g/\ell)\theta da T=2π/gT = 2\pi\sqrt{\ell/g}. (d) Con EpE_p medida desde el equilibrio, una oscilación armónica tiene promedios temporales iguales de energía cinética y potencial, de modo que L=0\langle L\rangle = 0 y S=0S = 0 en un periodo completo — el cálculo lo confirma: S=0T12m2θ02[ω2sin2ωt(g/)cos2ωt] ⁣dt=0S = \int_0^T\tfrac12 m\ell^2\theta_0^2[\omega^2\sin^2 \omega t - (g/\ell)\cos^2\omega t]\,\dd t = 0 puesto que ω2=g/\omega^2 = g/\ell.

Ejercicio 1.3

Una máquina de Atwood: las masas m1m_1 y m2m_2 cuelgan de un hilo ideal que pasa por una polea sin masa. (a) Elíjase una coordenada y escríbase LL. (b) Hállese la aceleración. (c) El tratamiento de Newton necesitaba la tensión: ¿adónde ha ido a parar? (d) ¿Cómo recuperaría la tensión una vez conocido el movimiento?

Solución

Solución de Ejercicio 1.3.

(a) Sea xx el descenso de m1m_1 (de modo que m2m_2 sube xx): L=12(m1+m2)x˙2+(m1m2)gxL = \tfrac12(m_1 + m_2)\dot x^2 + (m_1 - m_2)gx. (b) (m1+m2)x¨=(m1m2)g(m_1 + m_2)\ddot x = (m_1 - m_2)g: a=(m1m2)g/(m1+m2)a = (m_1 - m_2)g/(m_1 + m_2). (c) La tensión actúa en los dos extremos de un hilo inextensible: en cualquier desplazamiento permitido sus trabajos se cancelan, así que no entra nunca en LL. (d) Newton sobre m1m_1 solo: T=m1(ga)=2m1m2g/(m1+m2)T = m_1(g - a) = 2m_1m_2g/(m_1 + m_2).

Ejercicio 1.4

Una partícula libre en coordenadas cilíndricas: L=12m(r˙2+r2φ˙2+z˙2)L = \tfrac12 m(\dot r^2 + r^2\dot\varphi^2 + \dot z^2). (a) ¿Qué coordenadas son cíclicas y cuáles son los momentos conservados? (b) ¿Por qué no se conserva prp_r pese a que no actúa ninguna fuerza? (c) Escríbase la ecuación de Euler–Lagrange para rr e interprétese el término mrφ˙2mr\dot\varphi^2. (d) Compruébese que una recta recorrida a velocidad constante la resuelve.

Solución

Solución de Ejercicio 1.4.

(a) φ\varphi y zz: se conservan pφ=mr2φ˙p_\varphi = mr^2\dot\varphi (momento angular respecto del eje) y pz=mz˙p_z = m\dot z. (b) rr aparece en LL a través de r2φ˙2r^2\dot\varphi^2, luego no es cíclica: p˙r=mrφ˙20\dot p_r = mr\dot\varphi^2 \neq 0. Nada falla — pr=mr˙p_r = m\dot r es la componente radial de un vector constante p\vect p, y una componente a lo largo de una dirección que gira no tiene por qué ser constante. (c) mr¨=mrφ˙2m\ddot r = mr\dot\varphi^2: el término centrífugo, el precio de usar coordenadas ligadas a direcciones que rotan. (d) Para una recta a distancia bb recorrida a velocidad vv: r=b2+v2t2r = \sqrt{b^2 + v^2t^2}, r2φ˙=bvr^2\dot\varphi = bv; entonces r¨=b2v2/r3=r(bv/r2)2=rφ˙2\ddot r = b^2v^2/r^3 = r(bv/r^2)^2 = r\dot\varphi^2.

Ejercicio 1.5 ★★

Un bloque de masa mm desliza por la cara sin rozamiento (de ángulo α\alpha) de una cuña de masa MM, libre a su vez de deslizar sobre un suelo sin rozamiento. (a) Elíjanse dos coordenadas: la abscisa XX de la cuña y la distancia ss recorrida por el bloque a lo largo de la cara; escríbase LL. (b) ¿Qué coordenada es cíclica y qué ley de conservación expresa? (c) Hállense las dos aceleraciones. (d) Compruébense los límites MM \to \infty y α90\alpha \to 90^\circ.

Solución

Solución de Ejercicio 1.5.

(a) Posición del bloque (X+scosα, ssinα)(X + s\cos\alpha,\ -s\sin\alpha), luego

L=12MX˙2+12m(X˙2+2X˙s˙cosα+s˙2)+mgssinα.L = \tfrac12 M\dot X^2 + \tfrac12 m\big(\dot X^2 + 2\dot X\dot s\cos\alpha + \dot s^2\big) + mgs\sin\alpha .

(b) XX es cíclica: se conserva P=(M+m)X˙+ms˙cosαP = (M + m)\dot X + m\dot s\cos\alpha — la cantidad de movimiento horizontal total, ya que no actúa ninguna fuerza horizontal exterior. (c) La ecuación de ss es s¨+X¨cosα=gsinα\ddot s + \ddot X\cos\alpha = g\sin\alpha; con X¨=ms¨cosα/(M+m)\ddot X = -m\ddot s\cos\alpha/(M + m) obtenido en (b),

s¨=(M+m)gsinαM+msin2α,X¨=mgsinαcosαM+msin2α.\ddot s = \frac{(M + m)g\sin\alpha}{M + m\sin^2\alpha} , \qquad \ddot X = -\frac{mg\sin\alpha\cos\alpha}{M + m\sin^2\alpha} .

(d) MM \to \infty: s¨gsinα\ddot s \to g\sin\alpha, el plano inclinado fijo; α90\alpha \to 90^\circ: s¨g\ddot s \to g, X¨0\ddot X \to 0 — caída libre a lo largo de una cara vertical, sin empujar la cuña.

Ejercicio 1.6 ★★

El péndulo esférico: una lenteja en el extremo de una varilla de longitud \ell, libre en ambos ángulos (θ\theta desde la vertical descendente y φ\varphi alrededor de ella). (a) Escríbase LL. (b) Identifíquense la coordenada cíclica y el pφp_\varphi conservado. (c) Redúzcase el movimiento de θ\theta a un potencial efectivo y esbócese. (d) Para el movimiento cónico θ=θ0\theta = \theta_0, recupérese la relación del péndulo cónico cosθ0=g/ω2\cos\theta_0 = g/\ell\omega^2 del volumen del Año 1.

Solución

Solución de Ejercicio 1.6.

(a) L=12m2(θ˙2+sin2θφ˙2)+mgcosθL = \tfrac12 m\ell^2(\dot\theta^2 + \sin^2\theta\,\dot\varphi^2) + mg\ell\cos\theta. (b) φ\varphi cíclica: pφ=m2sin2θφ˙p_\varphi = m\ell^2\sin^2\theta\,\dot\varphi, el momento angular vertical. (c) Al eliminar φ˙\dot\varphi, se conserva 12m2θ˙2+Uef(θ)\tfrac12 m\ell^2\dot\theta^2 + U_{\text{ef}}(\theta) con

Uef(θ)=pφ22m2sin2θmgcosθ:U_{\text{ef}}(\theta) = \frac{p_\varphi^2}{2m\ell^2\sin^2\theta} - mg\ell\cos\theta :

una barrera en θ=0\theta = 0 y en θ=π\theta = \pi (si pφ0p_\varphi \neq 0) con un único mínimo intermedio — la lenteja nuta entre dos circunferencias. (d) En el mínimo, θ˙=0\dot\theta = 0 con θ=θ0\theta = \theta_0 constante: pφ2cosθ0/m2sin3θ0=mgsinθ0p_\varphi^2\cos\theta_0/m\ell^2\sin^3\theta_0 = mg\ell\sin\theta_0; al insertar pφ=m2sin2θ0ωp_\varphi = m\ell^2\sin^2\theta_0\,\omega resulta ω2cosθ0=g\ell\omega^2\cos\theta_0 = g.

Ejercicio 1.7 ★★

Un péndulo (masa mm, longitud \ell) cuelga de un carro de masa MM libre de rodar sobre un carril horizontal. (a) Con las coordenadas XX (carro) y θ\theta, escríbase LL. (b) ¿Qué se conserva y por qué, físicamente? (c) Linealícese para θ\theta pequeño y muéstrese que la frecuencia de oscilación es Ω=(1+m/M)g/\Omega = \sqrt{(1 + m/M)\,g/\ell}. (d) Explíquense los límites MM \to \infty y M0M \to 0: ¿por qué un carro ligero eleva la frecuencia?

Solución

Solución de Ejercicio 1.7.

(a) Lenteja en (X+sinθ, cosθ)(X + \ell\sin\theta,\ -\ell\cos\theta):

L=12(M+m)X˙2+mcosθX˙θ˙+12m2θ˙2+mgcosθ.L = \tfrac12(M + m)\dot X^2 + m\ell\cos\theta\,\dot X\dot\theta + \tfrac12 m\ell^2\dot\theta^2 + mg\ell\cos\theta .

(b) XX cíclica: se conserva (M+m)X˙+mcosθθ˙(M + m)\dot X + m\ell\cos\theta\,\dot\theta — la cantidad de movimiento horizontal del conjunto (el carril solo empuja verticalmente). (c) Ángulos pequeños: (M+m)X¨+mθ¨=0(M + m)\ddot X + m\ell\ddot\theta = 0 y θ¨+X¨+gθ=0\ell\ddot\theta + \ddot X + g\theta = 0; al eliminar X¨\ddot X, θ¨[1m/(M+m)]=gθ\ell\ddot\theta\,[1 - m/(M + m)] = -g\theta, luego Ω2=(M+m)g/M=(1+m/M)g/\Omega^2 = (M + m)g/M\ell = (1 + m/M)\,g/\ell. (d) MM \to \infty: el pivote fijo, Ω2=g/\Omega^2 = g/\ell. Con MM pequeña, el carro retrocede en sentido opuesto a la lenteja y la oscilación se produce en torno a un punto intermedio (el centro de masas, que está fijo), lo que acorta el péndulo efectivo — de ahí la frecuencia más alta, que diverge cuando M0M \to 0.

Ejercicio 1.8 ★★

Una partícula de carga qq en un campo uniforme B=Bez\vect B = B\vect e_z, descrito por A=12Br\vect A = \tfrac12\vect B\wedge\vect r. (a) Escríbase LL en coordenadas cartesianas. (b) Dedúzcanse las ecuaciones del movimiento y compruébese que describen la circunferencia ciclotrónica a ωc=qB/m\omega_{\text{c}} = qB/m. (c) Calcúlense los momentos conjugados pxp_x, pyp_y: ¿son mx˙m\dot x, my˙m\dot y? ¿Se conservan? (d) Muéstrese que LL escrito en coordenadas cilíndricas tiene φ\varphi cíclica e identifíquese el pφp_\varphi conservado para una circunferencia centrada en el eje.

Solución

Solución de Ejercicio 1.8.

(a) A=12B(y,x,0)\vect A = \tfrac12 B(-y, x, 0): L=12m(x˙2+y˙2+z˙2)+12qB(xy˙yx˙)L = \tfrac12 m(\dot x^2 + \dot y^2 + \dot z^2) + \tfrac12 qB(x\dot y - y\dot x). (b) La ecuación de xx:  ⁣d(mx˙12qBy)/ ⁣dt=12qBy˙\dd(m\dot x - \tfrac12 qBy)/\dd t = \tfrac12 qB\dot y, es decir, mx¨=qBy˙m\ddot x = qB\dot y; análogamente my¨=qBx˙m\ddot y = -qB\dot x: movimiento circular a ωc=qB/m\omega_{\text{c}} = qB/m, con z˙\dot z constante. (c) px=mx˙12qBymx˙p_x = m\dot x - \tfrac12 qBy \neq m\dot x; ni xx ni yy son cíclicas (L/x=12qBy˙\partial L/\partial x = \tfrac12 qB\dot y), así que no se conserva ninguno de los dos momentos — solo combinaciones como mx˙qBym\dot x - qBy (compruébese que su derivada se anula). (d) En coordenadas cilíndricas Aφ=12BrA_\varphi = \tfrac12 Br: L=12m(r˙2+r2φ˙2+z˙2)+12qBr2φ˙L = \tfrac12 m(\dot r^2 + r^2\dot\varphi^2 + \dot z^2) + \tfrac12 qBr^2\dot\varphi; φ\varphi es cíclica y pφ=mr2φ˙+12qBr2p_\varphi = mr^2\dot\varphi + \tfrac12 qBr^2. Sobre una circunferencia de radio RR centrada en el eje, φ˙=ωc\dot\varphi = -\omega_{\text{c}}, luego pφ=qBR2+12qBR2=12qBR2p_\varphi = -qBR^2 + \tfrac12 qBR^2 = -\tfrac12 qBR^2.

Ejercicio 1.9 ★★

Para la cuenta sobre el aro que gira (Ejemplo 1.9): (a) calcúlese la función energía hh y verifíquese que la ecuación del movimiento la conserva; (b) calcúlese la energía mecánica EE y muéstrese que Eh=mω2R2sin2θE - h = m\omega^2R^2\sin^2\theta; (c) hállese la potencia entregada por el motor en función de θ\theta y θ˙\dot\theta; (d) hállese la frecuencia de las pequeñas oscilaciones en torno al equilibrio inclinado cuando ω2>g/R\omega^2 > g/R y muéstrese que se anula cuando ω2g/R\omega^2 \to g/R — la ralentización que anuncia el desdoblamiento del pozo.

Solución

Solución de Ejercicio 1.9.

(a) h=θ˙L/θ˙L=12mR2θ˙212mω2R2sin2θmgRcosθ=12mR2θ˙2+Uefh = \dot\theta\,\partial L/\partial\dot\theta - L = \tfrac12 mR^2\dot\theta^2 - \tfrac12 m\omega^2R^2\sin^2\theta - mgR\cos\theta = \tfrac12 mR^2\dot\theta^2 + U_{\text{ef}};  ⁣dh/ ⁣dt=θ˙[mR2θ¨+Uef(θ)]=0\dd h/\dd t = \dot\theta\,[mR^2\ddot\theta + U_{\text{ef}}'(\theta)] = 0 por la ecuación del movimiento. (b) E=12mR2θ˙2+12mω2R2sin2θmgRcosθ=h+mω2R2sin2θE = \tfrac12 mR^2\dot\theta^2 + \tfrac12 m\omega^2R^2\sin^2\theta - mgR\cos\theta = h + m\omega^2R^2\sin^2\theta. (c) P= ⁣dE/ ⁣dt=mω2R2sin2θθ˙P = \dd E/\dd t = m\omega^2R^2\sin2\theta\,\dot\theta: positiva mientras la cuenta se aleja del eje (el motor trabaja contra la inercia de la cuenta) y negativa a la vuelta. (d) Uef=mgRcosθmω2R2cos2θU_{\text{ef}}'' = mgR\cos\theta - m\omega^2R^2\cos2\theta; en cosθeq=g/ω2R\cos\theta_{\text{eq}} = g/\omega^2R vale mω2R2sin2θeqm\omega^2R^2\sin^2\theta_{\text{eq}}, de modo que ωosc=ωsinθeq=ω1(g/ω2R)20\omega_{\text{osc}} = \omega\sin\theta_{\text{eq}} = \omega\sqrt{1 - (g/\omega^2R)^2} \to 0 cuando ω2g/R\omega^2 \to g/R: la fuerza recuperadora se aplana justo cuando los pozos se funden.

Ejercicio 1.10 ★★★

El péndulo doble de brazos iguales (m1=m2=mm_1 = m_2 = m, 1=2=\ell_1 = \ell_2 = \ell). (a) Muéstrese que para ángulos pequeños L=12m2(2θ˙12+2θ˙1θ˙2+θ˙22)12mg(2θ12+θ22)L = \tfrac12 m\ell^2(2\dot\theta_1^2 + 2\dot\theta_1\dot\theta_2 + \dot\theta_2^2) - \tfrac12 mg\ell(2 \theta_1^2 + \theta_2^2). (b) Escríbanse las dos ecuaciones del movimiento. (c) Hállense las frecuencias normales Ω±2=(22)g/\Omega_\pm^2 = (2 \mp \sqrt2)\,g/\ell y la forma de cada modo (θ2=±2θ1\theta_2 = \pm\sqrt2\,\theta_1). (d) A gran amplitud este sistema es un ejemplo clásico de caos: explíquese en unas líneas qué rompe el análisis de ángulos pequeños y por qué no se pierde ninguna ley de conservación.

Solución

Solución de Ejercicio 1.10.

(a) Posiciones x2=(sinθ1+sinθ2)x_2 = \ell(\sin\theta_1 + \sin\theta_2), etc.; conservando los términos cuadráticos, el término cruzado de velocidades es m2θ˙1θ˙2m\ell^2\dot\theta_1\dot\theta_2, lo que da el LL enunciado. (b) 2θ¨1+θ¨2=2ω02θ12\ddot\theta_1 + \ddot\theta_2 = -2\omega_0^2\theta_1 y θ¨1+θ¨2=ω02θ2\ddot\theta_1 + \ddot\theta_2 = -\omega_0^2\theta_2, con ω02=g/\omega_0^2 = g/\ell. (c) Al insertar θi=aicosΩt\theta_i = a_i\cos\Omega t: det(2ω022Ω2Ω2Ω2ω02Ω2)=Ω44ω02Ω2+2ω04=0\det\begin{pmatrix} 2\omega_0^2 - 2\Omega^2 & -\Omega^2\\ -\Omega^2 & \omega_0^2 - \Omega^2\end{pmatrix} = \Omega^4 - 4\omega_0^2\Omega^2 + 2\omega_0^4 = 0, luego Ω±2=(22)ω02\Omega_\pm^2 = (2 \mp \sqrt2)\omega_0^2; la segunda línea da a2/a1=Ω2/(ω02Ω2)=±2a_2/a_1 = \Omega^2/(\omega_0^2 - \Omega^2) = \pm\sqrt2: lentejas juntas (modo lento), lentejas opuestas (modo rápido). (d) A gran amplitud, los acoplamientos en sin\sin y cos\cos vuelven no lineales las ecuaciones; las soluciones ya no se superponen y condiciones iniciales vecinas se separan exponencialmente (caos). La energía sigue conservándose exactamente — LL no depende explícitamente del tiempo —: el caos es cuestión de predictibilidad, no de conservación.

Ejercicio 1.11 ★★★

La braquistócrona. Una cuenta desliza sin rozamiento partiendo del reposo en el origen por una curva y(x)y(x) (yy hacia abajo) hasta un punto (a,b)(a, b). (a) Usando la conservación de la energía, muéstrese que el tiempo de descenso es T=0a(1+y2)/2gy ⁣dxT = \int_0^a\sqrt{(1 + y'^2)/2gy}\,\dd x — un funcional en el que xx desempeña el papel del tiempo. (b) El integrando F(y,y)F(y, y') no depende explícitamente de xx: muéstrese que h=yF/yFh = y'\,\partial F/\partial y' - F es constante a lo largo de la curva óptima (el mismo cálculo que en la Proposición 1.14). (c) Dedúzcase y(1+y2)=2ry(1 + y'^2) = 2r para cierta constante rr y verifíquese que la cicloide x=r(ϕsinϕ)x = r(\phi - \sin\phi), y=r(1cosϕ)y = r(1 - \cos\phi) la satisface. (d) Muéstrese que el descenso de la cuenta hasta el fondo de un arco dura πr/g\pi\sqrt{r/g} y compárese con la rampa recta hasta el mismo punto.

Solución

Solución de Ejercicio 1.11.

(a) v=2gyv = \sqrt{2gy} y  ⁣ds=1+y2 ⁣dx\dd s = \sqrt{1 + y'^2}\,\dd x dan el funcional. (b) El cálculo de la Proposición 1.14 con xx en el papel del tiempo:  ⁣dh/ ⁣dx=F/x=0\dd h/\dd x = -\partial F/\partial x = 0. (c) h=1/2gy(1+y2)h = -1/\sqrt{2gy(1 + y'^2)}, luego y(1+y2)=2ry(1 + y'^2) = 2r. Para la cicloide, y=sinϕ/(1cosϕ)y' = \sin\phi/(1 - \cos\phi) y 1+y2=2/(1cosϕ)1 + y'^2 = 2/(1 - \cos\phi), de donde y(1+y2)=2ry(1 + y'^2) = 2r. (d)  ⁣dt= ⁣ds/v=r/g ⁣dϕ\dd t = \dd s/v = \sqrt{r/g}\,\dd\phi (toda la dependencia en ϕ\phi se cancela), así que el fondo (ϕ=π\phi = \pi) se alcanza en πr/g\pi\sqrt{r/g} — sea cual sea el punto de partida: la cicloide es además la tautócrona. La rampa recta hasta (πr,2r)(\pi r, 2r) tarda π2+4r/g3.72r/g\sqrt{\pi^2 + 4}\,\sqrt{r/g} \approx 3.72\sqrt{r/g}, alrededor de un 18%18\% más que πr/g\pi\sqrt{r/g}.

Ejercicio 1.12 ★★★

Fermat como mínima acción. La luz en un medio de índice n(y)n(y) viaja entre dos puntos en el tiempo mínimo (volumen del Año 2). (a) Muéstrese que el tiempo de recorrido a lo largo de y(x)y(x) es T=1cn(y)1+y2 ⁣dxT = \tfrac1c\int n(y)\sqrt{1 + y'^2}\,\dd x. (b) Como el integrando no depende explícitamente de xx, úsese la magnitud conservada hh del ejercicio anterior para probar que n(y)/1+y2=constn(y)\big/\sqrt{1 + y'^2} = \text{const}, y compruébese que es la ley de Snell nsini=constn\sin i = \text{const} para un rayo medido desde la vertical. (c) Sobre una carretera caliente el índice crece con la altura según n(y)n0(1+βy)n(y) \approx n_0(1 + \beta y); muéstrese que un rayo casi horizontal se curva con radio de curvatura R1/βR \approx 1/\beta. (d) Con β=1.2×105m1\beta = 1.2 \times 10^{-5}\,\mathrm{m}^{-1}, ¿desde qué distancia ve un conductor cuyos ojos están a 1.2m1.2\,\mathrm{m} sobre la carretera el espejismo de “agua” sobre ella?

Solución

Solución de Ejercicio 1.12.

(a)  ⁣dt= ⁣ds/(c/n)\dd t = \dd s/(c/n). (b) Se conserva h=n(y)/1+y2h = -n(y)/\sqrt{1 + y'^2}; 1/1+y2=cosθ1/\sqrt{1 + y'^2} = \cos\theta (θ\theta es el ángulo de la pendiente) =sini= \sin i con ii medido desde la vertical: nsini=constn\sin i = \text{const} — la ley de Snell, aplicada de forma continua. (c) Para un rayo casi horizontal, ncosθconstn\cos\theta \approx \text{const} con θ\theta pequeño da θ ⁣dθ= ⁣dn/nβ ⁣dy\theta\,\dd\theta = \dd n/n \approx \beta\,\dd y; como  ⁣dy=θ ⁣dx\dd y = \theta\,\dd x, la curvatura es  ⁣dθ/ ⁣dx=β\dd\theta/\dd x = \beta, es decir, R=1/β83kmR = 1/\beta \approx 83\,\mathrm{km}, y el rayo se curva hacia arriba (hacia nn mayor). (d) Un rayo que sale del ojo y llega tangente a la carretera tras curvarse con radio RR la toca a d=2hR=2×1.2×8.3e4450md = \sqrt{2hR} = \sqrt{2 \times 1.2 \times 8.3e4} \approx 450\,\mathrm{m}: más allá de esa distancia ya no se ve la carretera — se ve cielo refractado hacia arriba, el “agua” temblorosa.

El péndulo de Foucault en el Panteón de París. Coordenadas generalizadas, una ligadura y un sistema de referencia que gira lentamente: la rotación de la sala aparece en las ecuaciones del movimiento — y la deriva del plano de oscilación permite a un sótano medir el giro de la Tierra. Fotografía: Olga Khomitsevich, CC BY 2.0.
El péndulo de Foucault en el Panteón de París. Coordenadas generalizadas, una ligadura y un sistema de referencia que gira lentamente: la rotación de la sala aparece en las ecuaciones del movimiento — y la deriva del plano de oscilación permite a un sótano medir el giro de la Tierra. Fotografía: Olga Khomitsevich, CC BY 2.0.

1.6 Problema: la escoba que se sostiene cabeza abajo

Problema 1.1

Problema de fin de semana — el péndulo invertido de Kapitza

Un péndulo rígido puede mantenerse establemente por encima de su pivote si se agita este arriba y abajo con suficiente rapidez — un descubrimiento analizado por Kapitza en 1951, lo bastante llamativo como para parecer un truco de magia, y el principio por el que los campos eléctricos oscilantes atrapan iones aislados. Modelizamos el péndulo como una masa puntual mm en el extremo de una varilla rígida sin masa de longitud =40cm\ell = 40\,\mathrm{cm}; θ\theta es el ángulo medido desde la vertical descendente; g=9.81m/s2g = 9.81\,\mathrm{m}/\mathrm{s}^{2}.

Parte I — El péndulo rígido. El pivote se mantiene fijo al principio.

  1. Escríbanse el lagrangiano y la ecuación del movimiento.
  2. Dese la frecuencia f0f_0 de las pequeñas oscilaciones en torno a θ=0\theta = 0 y su valor.
  3. Muéstrese que aquí h=Eh = E y úsese su conservación para hallar la velocidad mínima de lanzamiento de la lenteja, en el punto más bajo, que la lleva hasta arriba.
  4. Linealícese la ecuación cerca de θ=π\theta = \pi (hágase θ=π+ϵ\theta = \pi + \epsilon) y muéstrese que ϵ¨=+(g/)ϵ\ddot\epsilon = +(g/\ell)\epsilon: ¿con qué rapidez crece una inclinación inicial de una milésima de grado? Dese el tiempo que tarda en crecer un factor ee.
  5. Una escoba en equilibrio sobre la yema de un dedo cae en cosa de un segundo y, sin embargo, nunca se sostiene sola: enúnciese en una frase lo que la ecuación linealizada dice del equilibrio θ=π\theta = \pi.

Parte II — Agitando el pivote. El pivote oscila ahora verticalmente, con altura ys(t)=acosΩty_{\text{s}}(t) = a\cos\Omega t, donde a=2.0cma = 2.0\,\mathrm{cm} y Ω\Omega es ajustable.

  1. Escríbanse las coordenadas de la lenteja y muéstrese que v2=2θ˙22aΩsin(Ωt)sinθθ˙+a2Ω2sin2Ωt\vect v^{\,2} = \ell^2\dot\theta^2 - 2a\Omega\ell\sin(\Omega t)\sin\theta\,\dot\theta + a^2\Omega^2\sin^2\Omega t.
  2. Muéstrese que dos lagrangianos que difieren en una derivada temporal total  ⁣dF(q,t)/ ⁣dt\dd F(q,t)/\dd t dan las mismas ecuaciones de Euler–Lagrange.
  3. Usando esa libertad, redúzcase el lagrangiano a L=12m2θ˙2maΩ2cos(Ωt)cosθ+mgcosθL = \tfrac12 m \ell^2\dot\theta^2 - ma\Omega^2\ell\cos(\Omega t)\cos\theta + mg\ell\cos\theta. (Indicación: sin(Ωt)sinθθ˙\sin(\Omega t)\sin\theta\,\dot\theta se combina con un término cos(Ωt)cosθ\cos(\Omega t)\cos\theta para dar una derivada total; los términos que solo dependen de tt pueden suprimirse).
  4. Dedúzcase la ecuación del movimiento

    θ¨=gsinθ+aΩ2cos(Ωt)sinθ.\ddot\theta = -\frac{g}{\ell}\sin\theta + \frac{a\Omega^2}{\ell}\cos(\Omega t)\sin\theta .

    Interprétese el segundo término como la sustitución del peso por g+y¨sg + \ddot y_{\text{s}}: el péndulo vive en un ascensor.

  5. ¿Se conserva ahora la función energía hh? ¿Y la energía? ¿Qué bombea energía hacia dentro y hacia fuera?
  6. En el régimen de Kapitza aa \ll \ell y Ωω0=g/\Omega \gg \omega_0 = \sqrt{g/\ell}: compruébense para a=2cma = 2\,\mathrm{cm}, =40cm\ell = 40\,\mathrm{cm}, Ω/2π=40Hz\Omega/2\pi = 40\,\mathrm{Hz}, y explíquese físicamente por qué la lenteja no puede entonces seguir la excitación.

Parte III — Separar lo rápido de lo lento. Búsquese el movimiento en la forma θ(t)=Θ(t)+ξ(t)\theta(t) = \Theta(t) + \xi(t): una deriva lenta Θ\Theta más un pequeño rizado ξ\xi a la frecuencia de la excitación.

  1. Conservando solo el término dominante de cada miembro, muéstrese que el rizado cumple ξ¨+(aΩ2/)cos(Ωt)sinΘ\ddot\xi \approx +(a\Omega^2/\ell)\cos(\Omega t) \sin\Theta, con Θ\Theta congelada en la escala de tiempos de la excitación.
  2. Dedúzcase ξ(t)=(a/)cos(Ωt)sinΘ\xi(t) = -(a/\ell)\cos(\Omega t)\sin\Theta y compruébese que su amplitud es pequeña, del orden de a/a/\ell.
  3. Desarróllese sinθ=sin(Θ+ξ)\sin\theta = \sin(\Theta + \xi) a primer orden en ξ\xi e insértese en la ecuación del movimiento.
  4. Promédiese sobre un periodo de la excitación, con Θ\Theta fija: usando cosΩt=0\langle\cos\Omega t\rangle = 0 y cos2Ωt=12\langle\cos^2\Omega t\rangle = \tfrac12, muéstrese que

    Θ¨=gsinΘa2Ω222sinΘcosΘ.\ddot\Theta = -\frac{g}{\ell}\sin\Theta - \frac{a^2\Omega^2}{2\ell^2}\sin\Theta\cos\Theta .
  5. Muéstrese que se trata de un movimiento en el potencial efectivo

    Uef(Θ)=mg(cosΘ+a2Ω24gsin2Θ).U_{\text{ef}}(\Theta) = mg\ell\Big({-\cos\Theta} + \frac{a^2\Omega^2}{4g\ell}\sin^2\Theta\Big) .
  6. Compárese con la cuenta sobre el aro que gira (Ejemplo 1.9): la misma matemática, con el signo opuesto en el término nuevo. ¿Qué le hace la agitación al equilibrio inferior que la rotación no hacía?
  7. Esbócese UefU_{\text{ef}} para excitación lenta y para excitación rápida y descríbanse todos los equilibrios y su estabilidad en cada caso.

Parte IV — La escoba se pone de pie.

  1. Desarrollando UefU_{\text{ef}} cerca de Θ=π\Theta = \pi, muéstrese que la posición invertida es estable exactamente cuando

    a2Ω2>2g.a^2\Omega^2 > 2g\ell .
  2. Calcúlese la frecuencia crítica de excitación fc=Ωc/2πf_{\text{c}} = \Omega_{\text{c}}/2\pi para nuestro péndulo, así como la velocidad máxima del pivote aΩca\Omega_{\text{c}} y la aceleración aΩc2a\Omega_{\text{c}}^2 (en unidades de gg) que exige.
  3. Para Ω=2Ωc\Omega = 2\Omega_{\text{c}}, hállese la frecuencia del lento balanceo del péndulo erguido en torno a Θ=π\Theta = \pi, y su valor; compruébese que es en efecto lenta comparada con la excitación.
  4. Todavía con Ω=2Ωc\Omega = 2\Omega_{\text{c}}, ¿cuál es la amplitud del rizado ξ\xi para Θ\Theta ligeramente apartada de π\pi, en grados, si a/=0.05a/\ell = 0.05? ¿Delataría una fotografía el truco?
  5. ¿Cuánto puede inclinarse la escoba respecto de la vertical y volver todavía? Muéstrese que el pozo de la posición erguida se extiende sobre Θπ<arccos(2g/a2Ω2)|\Theta - \pi| < \arccos(2g\ell/a^2\Omega^2), y evalúese en Ω=2Ωc\Omega = 2\Omega_{\text{c}}.
  6. Un malabarista que mantiene una escoba en equilibrio sobre la yema quieta de un dedo también la sostiene: ¿mediante qué mecanismo completamente distinto? Nómbrese la característica del péndulo de Kapitza que no necesita realimentación.
  7. Resúmase el resultado con nombre propio: un pivote agitado con amplitud 2cm2\,\mathrm{cm} a 45Hz45\,\mathrm{Hz} — el doble de la frecuencia crítica — sostiene cabeza abajo un péndulo de 40cm40\,\mathrm{cm}, en un pozo que alcanza unos 7575^\circ desde la vertical y que se balancea suavemente a unos 1.4Hz1.4\,\mathrm{Hz}. ¿Dónde se usa en física el mismo atrapamiento promediado para retener una única partícula cargada?
Solución

Solución de Problema 1.1.

1. L=12m2θ˙2+mgcosθL = \tfrac12 m\ell^2\dot\theta^2 + mg\ell\cos\theta; θ¨=(g/)sinθ\ddot\theta = -(g/\ell)\sin\theta. 2. f0=g//2π=0.79Hzf_0 = \sqrt{g/\ell}/2\pi = 0.79\,\mathrm{Hz}. 3. La suspensión está fija, luego h=E=12m2θ˙2mgcosθh = E = \tfrac12 m\ell^2\dot\theta^2 - mg\ell\cos\theta; del fondo a la cima, 12mv2=2mg\tfrac12 mv^2 = 2mg\ell: v=2g=4.0m/sv = 2\sqrt{g\ell} = 4.0\,\mathrm{m}/\mathrm{s}. 4. ϵ¨=+(g/)ϵ\ddot\epsilon = +(g/\ell)\epsilon: ϵet/τ\epsilon \propto \eu^{t/\tau} con τ=/g=0.20s\tau = \sqrt{\ell/g} = 0.20\,\mathrm{s}. 5. El equilibrio invertido existe, pero es exponencialmente inestable: cualquier inclinación, por pequeña que sea, se multiplica por ee cada 0.2s0.2\,\mathrm{s}. 6. Lenteja en (sinθ, acosΩtcosθ)(\ell\sin\theta,\ a\cos\Omega t - \ell\cos\theta); dedúzcase y eleve al cuadrado: el v2\vect v^{\,2} enunciado, con el término cruzado 2aΩsin(Ωt)sinθθ˙-2a\Omega\ell\sin(\Omega t)\sin\theta\,\dot\theta. 7. Si L=L+ ⁣dF(q,t)/ ⁣dtL' = L + \dd F(q,t)/\dd t, la acción cambia en F(q2,t2)F(q1,t1)F(q_2, t_2) - F(q_1, t_1), una constante frente a variaciones con extremos fijos: los mismos caminos estacionarios, las mismas ecuaciones. 8. El término cruzado es maΩsin(Ωt)sinθθ˙= ⁣d[maΩsin(Ωt)cosθ]/ ⁣dtmaΩ2cos(Ωt)cosθ-ma\Omega\ell\sin(\Omega t)\sin\theta\, \dot\theta = \dd[ma\Omega\ell\sin(\Omega t)\cos\theta]/\dd t - ma\Omega^2\ell\cos(\Omega t)\cos\theta; al suprimir la derivada total y los términos que solo dependen de tt (12ma2Ω2sin2Ωt\tfrac12 ma^2\Omega^2\sin^2\Omega t y mgacosΩt-mga\cos\Omega t) queda el LL enunciado. 9. m2θ¨=mgsinθ+maΩ2cos(Ωt)sinθm\ell^2\ddot\theta = -mg\ell\sin\theta + ma\Omega^2\ell\cos(\Omega t)\sin\theta. Como y¨s=aΩ2cosΩt\ddot y_{\text{s}} = -a\Omega^2\cos\Omega t, esto es θ¨=[(g+y¨s)/]sinθ\ddot\theta = -[(g + \ddot y_{\text{s}})/\ell]\sin\theta: en el sistema del pivote la gravedad aparente oscila. 10. Ahora LL depende explícitamente de tt: no se conserva hh, y tampoco EE — el agitador mete y saca energía por el pivote. 11. a/=0.05a/\ell = 0.05; Ω=251rad/s\Omega = 251\,\mathrm{rad}/\mathrm{s} frente a ω0=4.9rad/s\omega_0 = 4.9\,\mathrm{rad}/\mathrm{s}: razón 5151. En un periodo de la excitación (25ms25\,\mathrm{ms}) la gravedad apenas cambia θ˙\dot\theta: la lenteja es demasiado perezosa para seguirla y solo tiembla. 12. ξ¨\ddot\xi es la derivada mayor (Ω2\propto\Omega^2) y la excitación el mayor término de fuerza: ξ¨=(aΩ2/)cos(Ωt)sinΘ\ddot\xi = (a\Omega^2/\ell)\cos(\Omega t)\sin\Theta. 13. Integrando dos veces con Θ\Theta fija: ξ=(a/)cos(Ωt)sinΘ\xi = -(a/\ell)\cos(\Omega t)\sin\Theta, de amplitud como mucho a/=0.05a/\ell = 0.05: un temblor de dos grados. 14. sinθsinΘ+ξcosΘ\sin\theta \approx \sin\Theta + \xi\cos\Theta, luego

Θ¨+ξ¨=g(sinΘ+ξcosΘ)+aΩ2cosΩt(sinΘ+ξcosΘ).\ddot\Theta + \ddot\xi = -\frac{g}{\ell}(\sin\Theta + \xi\cos\Theta) + \frac{a\Omega^2}{\ell}\cos\Omega t\,(\sin\Theta + \xi\cos\Theta) .

15. El promedio mata ξ¨\ddot\xi, cosΩt\langle\cos\Omega t\rangle y ξ\langle\xi\rangle; el término cruzado que sobrevive es (aΩ2/)cosΘξcosΩt=(a2Ω2/22)sinΘcosΘ(a\Omega^2/\ell)\cos\Theta\,\langle\xi\cos\Omega t\rangle = -(a^2\Omega^2/2\ell^2)\sin\Theta\cos\Theta, lo que da la ecuación enunciada para Θ\Theta. 16. Θ¨=(1/m2)Uef(Θ)\ddot\Theta = -(1/m\ell^2)\,U_{\text{ef}}'(\Theta) con Uef=mg[cosΘ+(a2Ω2/4g)sin2Θ]U_{\text{ef}} = mg\ell[-\cos\Theta + (a^2\Omega^2/4g\ell) \sin^2\Theta] — dedúzcase para comprobarlo. 17. El mismo término en sin2\sin^2 que el del aro, pero con el signo opuesto: la rotación cavaba pozos en los flancos y solo podía aplanar el fondo; la agitación vertical endurece el pozo del fondo y cava un pozo nuevo arriba. 18. Excitación lenta (a2Ω2<2ga^2\Omega^2 < 2g\ell): mínimo en Θ=0\Theta = 0 y máximo en π\pi — nada nuevo. Excitación rápida: mínimos en 00 y en π\pi, separados por máximos en cosΘ=2g/a2Ω2\cos\Theta^* = -2g\ell/a^2\Omega^2; tanto el péndulo colgante como el erguido oscilan establemente. 19. Cerca de π\pi, con Θ=π+ϵ\Theta = \pi + \epsilon: ϵ¨=[g/a2Ω2/22]ϵ\ddot\epsilon = [g/\ell - a^2\Omega^2/2\ell^2]\,\epsilon; la estabilidad exige que el corchete sea negativo: a2Ω2>2ga^2\Omega^2 > 2g\ell. 20. Ωc=2g/a=140rad/s\Omega_{\text{c}} = \sqrt{2g\ell}/a = 140\,\mathrm{rad}/\mathrm{s}: fc=22Hzf_{\text{c}} = 22\,\mathrm{Hz}; velocidad máxima aΩc=2.8m/sa\Omega_{\text{c}} = 2.8\,\mathrm{m}/\mathrm{s}, aceleración máxima aΩc2=392m/s240ga\Omega_{\text{c}}^2 = 392\,\mathrm{m}/\mathrm{s}^{2} \approx 40g. 21. ωlento=a2Ω2/22g/\omega_{\text{lento}} = \sqrt{a^2\Omega^2/2\ell^2 - g/\ell}; en Ω=2Ωc\Omega = 2\Omega_{\text{c}}, a2Ω2=8ga^2\Omega^2 = 8g\ell, luego ωlento=3g/=8.6rad/s\omega_{\text{lento}} = \sqrt{3g/\ell} = 8.6\,\mathrm{rad}/\mathrm{s}: 1.4Hz1.4\,\mathrm{Hz}, treinta veces más lento que la excitación de 45Hz45\,\mathrm{Hz}. 22. ξmax=(a/)sinΘ\xi_{\max} = (a/\ell)\sin\Theta: inclinado 1010^\circ respecto de la vertical, ξmax=0.05sin170=8.7mrad0.5\xi_{\max} = 0.05\sin170^\circ = 8.7\,\mathrm{mrad} \approx 0.5^\circ — una fotografía ordinaria muestra una escoba inmóvil. 23. El pozo de la posición erguida llega hasta los máximos que lo flanquean: Θπ<arccos(2g/a2Ω2)|\Theta - \pi| < \arccos(2g\ell/a^2\Omega^2); en Ω=2Ωc\Omega = 2\Omega_{\text{c}}, arccos14=75\arccos\tfrac14 = 75^\circ — un pozo notablemente indulgente. 24. El malabarista usa realimentación: los ojos miden la inclinación y la mano acelera lateralmente para cancelarla. La estabilización de Kapitza es en lazo abierto — la excitación nunca sabe dónde está el péndulo. 25. Agitado a 45Hz45\,\mathrm{Hz} con 2cm2\,\mathrm{cm} de amplitud, el péndulo de 40cm40\,\mathrm{cm} se mantiene invertido en un pozo de 7575^\circ y se balancea a 1.4Hz1.4\,\mathrm{Hz}. El mismo potencial efectivo promediado, creado con un campo eléctrico cuadrupolar oscilante en lugar de un pivote agitado, confina iones aislados en la trampa de Paul — el caballo de batalla de los relojes atómicos y de la computación cuántica con iones atrapados.

Términos definidos en este capítulo

Ver los 431 términos del glosario