Formalismo lagrangiano y hamiltoniano para péndulo amortiguado

Supongamos que tengo un péndulo simplificado (una cuerda sin masa de longitud yo , bola de masa metro y algo de volumen distinto de cero). Quiero derivar las ecuaciones de movimiento, pero quiero tener en cuenta la resistencia del aire. Sé cómo hacerlo aplicando la segunda ley de Newton, pero ¿hay alguna forma de hacerlo usando lagrangiana o hamiltoniana? ¿Cuál sería la forma de lagrangiano/hamiltoniano en tal caso?

Respuestas (3)

Para una fuerza general de amortiguamiento no existe un Principio de Acción Mínima. Sin embargo, puede obtener las ecuaciones de movimiento del principio de d'Alembert. El procedimiento se explica aquí y las ecuaciones de los movimientos son

d d t L q ˙ i L q i = q i pag ,
dónde L = T V es la función de Lagrange y
q i pag α F α r α q i ,
es la llamada fuerza generalizada que depende de la fuerza de amortiguamiento F α actuando de partícula α .

Dado que para el amortiguamiento general no podemos formular un principio variacional, L en sí mismo no es suficiente para determinar la evolución dinámica del sistema. Por lo tanto no hay hamiltoniano H satisfactorio

d S = d t 1 t 2 [ i pag i q ˙ i H ( q , pag , t ) ] d t = 0.
Recuerde que las ecuaciones canónicas también se pueden obtener del principio variacional anterior, por lo que no hay función H satisfaciendo las ecuaciones canónicas y las ecuaciones de movimiento del sistema (que tendrán en cuenta el amortiguamiento) simultáneamente.

En el caso de la amortiguación lineal, puede seguir la idea de Bateman que se menciona en la respuesta de ZeroTheHero o usar la función de disipación de Rayleigh

Depende de la forma del amortiguamiento, pero si escribes

L = mi α t ( 1 2 metro q ˙ 2 V ( q ) )
entonces obtienes una ecuación de movimiento de la forma
mi α t ( metro q ¨ + metro α q ˙ + d V d q ) = 0
lo que puede explicar la fricción dependiendo de su modelo de fricción.

El hamiltoniano se calcularía de la forma habitual pero claramente L depende explícitamente de t por lo que la energía no se conserva en su sistema.


Editar: cambiado α α después de un comentario correcto.

Por qué mi a t ? No veo de dónde vino. También olvidé mencionar que me gustaría que la amortiguación dependiera del cuadrado de la velocidad. Bueno, la razón por la que hago esa pregunta es que, en última instancia, quiero simular un péndulo doble con resistencia del aire y derivar ecuaciones del lagrangiano parece ser más fácil que el método de la fuerza.
Creo que quieres el signo opuesto para α .
@JG De hecho. Acabo de hacer la corrección.
@G.Fil Hasta donde yo sé, esto mi α t factor es fenomenológico, es decir, se observa que produce el término de amortiguamiento adecuado. No conozco ningún truco que capture una amortiguación cuadrática en q ˙ : es posible que deba usar las ecuaciones de movimiento generalizadas, donde el término que no proviene de un potencial se agrega "a mano".
No es fenomenológico; hay una razón profunda para la elección. Esto es un poco exagerado, pero el mi α t El factor es análogo al | gramo | factor en la relatividad general si consideramos el tiempo como un 1 Variedad riemanniana bidimensional. Entonces q es un campo en esa variedad, con la ecuación de movimiento q ˙ término proporcional al término de Christoffel Γ 00 0 = 1 2 gramo 00 0 gramo 00 = 0 en gramo , de donde el gramo = mi 2 α t sigue la elección. (Efectivamente, es una variedad de De Sitter).
@JGwow! ¡Acabo de aprender algo! Pero, ¿por qué querría pensar en el tiempo de esta manera? ¿Proporciona esto información adicional que ayudaría a manejar casos adicionales (más complicados)?

En cierto sentido sí; deberá escribir la fuerza de la resistencia del aire como si pudiera derivarse de un potencial dependiente de la velocidad.

Un ejemplo común de esto es una partícula en un campo magnético . La fuerza del campo magnético sobre la partícula depende de su velocidad, por lo que el potencial se escribe como A v .