MECÁNICA
PRÁCTICA 4: Mecánica Hamiltoniana y Sistemas Dinámicos
FCEN - Universidad Nacional de Cuyo - 2023
FORMALISMO HAMILTONIANO
1) Escribir el Hamiltoniano, las ecuaciones de Hamilton y encontrar las ecuaciones de
movimiento para:
a) Un proyectil que puede ser considerado como una partícula, moviéndose bajo la acción
de la gravedad. Despreciar la rotación de la tierra y la fricción del aire.
b) Una masa que pende de un resorte y oscila en un plano.
c) El sistema visto en el ejercicio 8 de la Práctica 1.
d) El sistema visto en el ejercicio 9 de la Práctica 1.
e) Una peonza simétrica con el extremo fijo, en un campo gravitatorio constante.
CORCHETES DE POISSON
2) Verificar las siguientes propiedades de los corchetes de Poisson:
a) [𝑞𝑖 , 𝑞𝑗] = 0 ; [𝑝𝑖 , 𝑝𝑗] = 0 ; [𝑞𝑖 , 𝑝𝑗] = δ𝑖𝑗
b) [𝑢 , 𝑣] = − [𝑣 , 𝑢]
c) [𝑢 , 𝑢] = 0
d) [𝑢 , 𝑐] = 0 con 𝑐 independiente de 𝑞 y 𝑝
e) [𝑢, 𝑣 + 𝑤] = [𝑢, 𝑣] + [𝑢, 𝑤]
f) [𝑢, 𝑣𝑤] = [𝑢, 𝑣]𝑤 + 𝑣[𝑢, 𝑤]
𝑑𝑢
3) Sea 𝑢 = 𝑢(𝑞, 𝑝, 𝑡), probar que 𝑑𝑡
puede escribirse como:
𝑑𝑢 ∂𝑢
𝑑𝑡
= [𝑢, 𝐻] + ∂𝑡
donde 𝐻 es el Hamiltoniano. De esta relación se deduce que si 𝑢 no depende explícitamente
del tiempo y [𝑢, 𝐻] = 0, entonces 𝑢 es una constante de movimiento.
4) Probar que el corchete de Poisson de dos constantes de movimiento es también una
constante de movimiento, aún cuando estas constantes dependan explícitamente del
tiempo.
Ayuda: usar el resultado del ejercicio anterior y la siguiente identidad:
[𝑢, [𝑣, 𝑤]] + [𝑣, [𝑤, 𝑢]] + [𝑤, [𝑢, 𝑣]] = 0 𝐼𝑑𝑒𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝐽𝑎𝑐𝑜𝑏𝑖
Nota: para saciar su curiosidad, en [Link] se puede ver una
demostración corta de la identidad de Jacobi, cuya demostración más directa es
algebraicamente trabajosa.
TRANSFORMACIONES CANÓNICAS
5) Considerando que 𝑞 y 𝑝 son variables canónicas, mostrar que las siguientes
transformaciones son canónicas:
𝑘
a) 𝑄 = 2𝑞 𝑒 cos 𝑝
−𝑘
𝑃 = 2𝑞 𝑒 sin 𝑝
donde 𝑘 es una constante.
2 2
b) 𝑄 = 𝑞 + 𝑘² 𝑝
1
𝑃 = 2 𝑘 arctan ( )𝑘𝑝
𝑞
donde 𝑘 es una constante.
6) Considerando la siguiente transformación canónica
𝑄 = ln(1 + 𝑞 cos 𝑝)
𝑃 = 2 𝑞 (1 + 𝑞 cos 𝑝) sin 𝑝
𝑄 2
probar que 𝐹3(𝑄, 𝑝) = − (𝑒 − 1) tan 𝑝 es la generatriz de dicha transformación.
MÉTODO DE HAMILTON - JACOBI
7) Aplicar completo el método de Hamilton - Jacobi al movimiento de una masa puntual que
desliza sin rozamiento sobre un plano inclinado de ángulo α respecto a la horizontal.
8) Utilizar las ecuaciones de Hamilton - Jacobi para resolver el problema de Kepler, i.e. una
partícula en un potencial central de forma V= - k/r.
9) Variables ángulo-acción: consideremos un sistema mecánico periódico, con trayectorias
cerradas Ci en el espacio de las fases. La integral de línea sobre esa trayectoria cerrada C,
Ji = ∫ pi dqi, se llama la integral de fase o la variable de acción. Podemos encontrar una
función generatriz S=S(qi,Ji) tal que pi=∂S/∂qi; ωι=∂S/∂Ji, donde ωi son las variables de
ángulo. Por lo tanto las nuevas ecuaciones de Hamilton son:
dωi/dt= ∂H’/∂Ji y dJi/dt=- ∂H’/∂ωi ,
donde H’=E depende sólo de las constantes Ji. Entonces ωi=fi t + ci, con fi y ci constantes, y fi
dando las “frecuencias” características del sistema.
a) Si Δωi representa el cambio en ωi luego de un ciclo completo en la coordenada qr, probar
que Δωi=1, si i=r y 0 si i≠r. Asumir que el orden de diferenciación e integración no es
importante.
b) Interpretar físicamente este resultado.
c) Encontrar las frecuencias del oscilador armónico con este método.
SISTEMAS DINÁMICOS
10) Para los siguientes sistemas dinámicos de segundo orden:
i) ẋ(t) = 2x(t) + y(t) ii) ẋ(t) = 4 - 2y(t) iii) ẋ(t) = x(t) - 1
ẏ(t) = x(t) + y(t) ẏ(t) = 12 - 3x2(t) ẏ(t) = x(t) - y2(t)
a) Encontrar los puntos críticos del sistema.
b) Analizar la estabilidad de los mismos.
c) ¿Existe un Hamiltoniano a partir del cual se obtengan estas ecuaciones? En caso de que
exista, hallarlo.
11) Para el sistema dinámico descripto por las ecuaciones de Lotka-Volterra útil para el
problema presa - depredador, expresadas a continuación,
ẋ = a x - b xy
ẏ = - c y + d xy
donde a,b,c,d son constantes positivas y donde x e y son dos variables independientes.
a) Encontrar los puntos críticos del sistema y analizar la estabilidad de los mismos.
b) ¿Existe un Hamiltoniano a partir del cual se obtengan estas ecuaciones de movimiento?
Justificar.
12) Las concentraciones químicas x(t) e y(t) de dos reactivos se relacionan por:
dx/dt = a – (1+b) x + x2 y ; dy/dt = bx - x2 y ; donde a y b son positivos.
Encontrar los puntos críticos del sistema, la matriz Jacobiana en los puntos críticos y describir el
tipo de punto crítico.
Dato interesante: cuando b=1+a2 existe una “bifurcación de Hopf” supercrítica, con un punto
crítico inestable junto a un ciclo limite ([Link]
13) Partiendo de la ecuación de movimiento de un péndulo simple de variable θ, tome como
variables dinámicas 𝑥 = θ y 𝑦 = θ̇, y escriba el problema como un sistema dinámico de
segundo orden. Halle los puntos críticos y el comportamiento del sistema alrededor de los
mismos según la aproximación lineal. Considerando que este sistema es Hamiltoniano, y que H
es una cantidad conservada, muestre que el centro predicho en la aproximación lineal es
verdaderamente un centro. Esquematice el retrato de fases del sistema.
14) Problema desafío:
Analice el siguiente sistema dinámico:
donde a es un parámetro.
a) Muestre que la linealización del sistema predice que el único punto crítico se encuentra en el
origen y es un centro para todo valor del parámetro de control a.
b) Muestre que la predicción del comportamiento alrededor del origen hallada mediante la
linealización del sistema es incorrecta, ya que en realidad el retrato de fases (x,y) muestra una
espiral estable para a < 0 e inestable para a > 0.
Ayuda: Una forma de dar con el comportamiento correcto del sistema es pasar las ecuaciones a
coordenadas polares (𝑟, θ), sugerido por el hecho de la aparición de 𝑥² + 𝑦² en ambas
ecuaciones. Luego es posible combinarlas de forma simple para desacoplarlas, hallando una
ecuación para 𝑟˙ y otra para θ̇. Finalmente, considerando la velocidad radial y la velocidad
tangencial para diferentes valores de 𝑟 en el plano (x,y), esboce los vectores velocidad y
utilícelos como guía para trazar los retratos de fase en (𝑥, 𝑦) correspondientes a los casos
𝑎 < 0 y 𝑎 > 0.
