Formalismo de Lagrange
Formalismo de Lagrange
Formalismo de Lagrange
1
Capı́tulo 1
Tenemos la impresión de nuestro curso previo que todos los problemas de mecáni-
ca clásica de reducen a resolver ecuaciones diferenciales del siguiente tipo:
.. (e)
X
mi r i = F i + F ji , (1.1)
j
1. Cuerpo rı́gido, donde la constricción consiste en que la distancia entre las dife-
rentes partı́culas rij no debe de cambiar.
2. Las cuentas de un abaco que estan restringidas a moverse en una sola dimen-
sión.
4. Una partı́cula limitada a moverse sobre la superficie externa de una esfera soli-
da.
2
a) Cuerpo rı́gido b) Abaco.
c) Moléculas contenidas en un recipiente. d) Partı́cula limitada a moverse sobre una superficie esférica.
3
por la superficie de una esfera, puede escribirse como una desigualdad
r2 − a2 ≥ 0, (1.4)
Otro ejemplo de una fuerza de restricción es la que hace que la una partı́cula se
mantenga siempre sobre la superficie de una esfera. Notemos que la fuerza normal
4
sobre la partı́cula es variable ya que depende de la posición y de su velocidad, por lo
tanto la restricción dependerá también de la posición y velocidad (figura 2).
Figura 1.2: Fuerza de restricción sobre una partı́cula que se mueve sobre una super-
ficie esferica.
5
varilla de masa despreciable e inextensible y una barra similar que fija a la primer
partı́cula), es claro que las coordenadas generalizadas son los ángulos θ1 , θ2 . (figura 3).
La coordenadas generalizadas, tiene un sentido distinto en términos de coordenadas
a las Cartesianas, estas son útiles cuando no hay constricción
Sin embargo, si las cosntricción es noholónimica, las ecuaciones que expresan las
constricciones pueden no ser usadas para eliminar las coordenadas dependientes.
Un ejemplo de esta situación, es un disco que rueda sobre un plano horizontal xy
limitado a moverse ta tal forma que siempre el plano del disco se encuentre vertical.
Las coordenadas utilizadas para describir el movimiento podrı́a ser las coordenadas
x, y del centro del disco, el ángulo de rotación φ sobre el eje del disco, y un ángulo θ
entre los ejes del disco por decir el eje x (figura 4). Como resultados de la constricción
.
la velocidad del centro del disco, v , tiene una magnitud proporcional a φ,
.
v=aφ (1.5)
.
x = vsenθ (1.6)
.
y = −vcosθ. (1.7)
6
Figura 1.4: Disco que rueda sobre un plano horizontal.
constricciones:
dx − asenθdφ = 0 (1.8)
dy + acosθdφ = 0. (1.9)
Es claro que las ecuaciones para las constricciones no puede reducirse a la forma
de la ecuación (2) y por lo tanto son noholonomicas. El usos de coordendas gene-
ralizadas tiene como resultado que las ecuaciones de movimiento una forma más
simplificada. Desafortunadamente no existe una regla para seleccionar el conjunto de
coordendas generalizadas las más adecuadas para un problema dado. Además como
hemos visto las coordenadas generalizadas, pueden estar definidas por un conjun-
. . .
to de cantidades que pueden incluir derivadas temporales q 1 , q 2 , ..., q i las cuales son
llamadas velocidades generalizada. En general, las componentes Cartesianas de la ve-
locidades dependen de las coordenadas generalizadas, las velocidades generalizadas
y el tiempo:
. . .
r i =r j (qj , qj , t) (1.10)
qj = qj (rr i , t) (1.11)
. . .
q j =q j (rr i , r , t) (1.12)
Ejemplo.
Encontrar un conjunto de coordenadas generalizadas para una partı́cula puntual
7
moviéndose sobre la superficie de una semiesfera de radio R la cual esta centrada en
el origen. Debido a que el movimiento siempre ocurre sobre la superficie, tenemos
x2 + y 2 + z 2 − R = 0, z ≥ 0 (1.13)
Escojamos nuestras coordenadas generalizadas como los cosenos del ángulo entre los
ejes x, y y z y la linea que conecta la partı́cula con el origen. Por lo tanto
x
cosθ1 = q1 = ,
R
y
cosθ2 = q2 = ,
R
z
cosθ3 = q3 = (1.14)
R
estos son los llamados cosenos directores y sabemos que la suma del cuadrado de
estos es igual a la unidad. Esto es,
8
el sistema esta en equilibrio; esto implica que el total de las fuerzas sobre la cada
partı́cula se anula, F i = 0. Por lo tanto el trabajo realizado por la fuerza F i en el
desplazamiento δrr i, también se anula. La suma de todos los trabajos sobre todas las
partı́culas deberı́an de ser cero también:
X
F i · δrr i = 0. (1.18)
i
(a)
F i = F i + f i, (1.19)
(a)
X X
F i · δrr i + f i · δrr i = 0. (1.20)
i i
Nos hemos restringido a sistemas para los cuales el trabajo total virtual de las
fuerzas de constricción es cero. Por lo tanto tenemos como la condición de equilibrio
de un sistema sobre el cual el trabajo virtual de las fuerzas aplicadas se anula:
(a)
X
F i · δrr i = 0. (1.21)
i
.
F i =p i , (1.22)
9
puede escribirse como
.
F i − p i = 0, (1.23)
la cual establece que las parı́ulas en un sistema que esta en equilibrio bajo una una
.
fuerza igual a la fuerza actual más un ”fuerza efectiva reversible”− p i . Sustituyendo
en la ecuación 18 tenemos
.
F (a)
X
(F i − p i ) · δr
r i = 0, (1.24)
i
el cual es llamado el principio de D’Alembert. El cambio de del lenguaje de r i a qj
comienza por las ecuaciones de tranfromación (6),
donde las Qj son llamadas las componentes de las fuerzas generalizadas, definidas
como
X ∂rr i
Qj = Fi · . (1.30)
∂qj
i
Notemos que las q no necesariamente tendrán dimensiones de distancia, y por lo tanto
las Q no necesariamente tendrán unidades de fuerza, pero el producto Qj δqj siempre
deberán tener dimensiones de trabajo. Consideremos ahora el según termino de la
ecuación 24, el cual puede escribirse como
X . X ..
p ·δrr i = mi r i ·δrr i . (1.31)
i i
10
Utilizando la ecuación 27, tenemos que
X .. ∂rr i
mi ±ri · δqj . (1.32)
∂qj
i,j
∂vv i ∂.
hP
∂rr i q. ri
i
. = k ∂qk k
+ ∂t (1.37)
∂ qj ∂ qj
.
P ∂rr i ∂ q k ∂ 2r i
= .
k ∂qk ∂ q j
+ . (1.38)
∂ q j ∂t
∂rr i
= ∂qj (1.39)
entonces el segundo termino del lado izquierdo de la ecuación (24) resulta ser
( " !# ! )
X d ∂ X1 ∂ X1
. mi vi2 − mi vi2 − Qj δqj . (1.43)
dt ∂ q j 2 ∂qj 2
j i i
11
P 1 2
Podemos identificar la energı́a cinética del sistema T como i 2 m i vi , el principio de
D’Alambert es: (" ! # )
X d ∂T ∂T
. − − Qj δqj = 0. (1.44)
dt ∂ qj ∂qj
j
por lo tanto !
d ∂T ∂T
. − = Qj . (1.46)
dt ∂ qj ∂qj
En total tendremos n ecuaciones. Cuando las fuerzas provienen de un potencial es-
calar V ,
F i = −∇i V. (1.47)
∂V
Qj = − . (1.50)
∂qj
12
es una función explicita del tiempo el sistema es conservativo. Por lo tanto podemos
.
incluir un termino que se la derivada parcial de V con respecto de q j :
!
d ∂(T − V ) ∂(T − V )
. − =0 (1.52)
dt ∂ qj ∂qj
L = T − V, (1.53)
13
función explicita solamente de las coordenadas de posición, entonces el sistema es
monogénico y también es conservativo.
Para sistemas sistemas monogénicos el principio de Hamilton puede escribirse
como: El movimiento del sistema del tiempo t1 al tiempo t2 es tal que la linea integral
(llamada acción o integral de acción)
Z t2
Ldt, (1.55)
t1
14
principio se considera como un postulado básico mejor que las leyes de Newton, ya
que tiene las sigientes ventajas:
∂T
=0y
∂xi
∂V
. =0
∂ xi
∂V d ∂T
− = . (1.59)
∂xi dx ∂ xi
∂U
− = Fi (1.60)
∂xi
15
y
3
d ∂T X 1 . 2 d . .
. m xj = (m xi ) =pi = Fi . (1.61)
dt ∂ xi 2 dt
j=1
Por lo tanto, las ecuaciones de Lagrange y Newton son idénticas si las coordenadas
generalizadas son las coordenadas rectangulares.
x = r cos θ
y = r sin θ (1.66)
16
Las componentes de la velocidad estan dadas por
. . .
x=r cos θ − r θ sin θ,
. . .
y =r sin θ + r θ cos θ. (1.67)
La energı́a cinética es
2 2
. . . .
1 r cos θ − r θ sin θ
T = 2m + r sin θ + r θ cos θ
h.
2 . . .2 .2 . . .2 i
1 r cos2 θ − 2 r r θ cos θ sin θ + r2 θ sin2 θ+ r sin2 θ + 2 r r θ sin θ cos θ + r2 θ cos2 θ
T = 2m
.
2 θ + sin2 θ + r 2 2 cos2 θ + sin2 θ
h. i
1 2
T = 2 m r cos θ
. i
2 2 .
h.
1 2
T = 2 m r +r θ (1.68)
∂rr
Qr = F · = F · r̂r = Fr ,
∂r
∂θθ
Qθ = F · = F · θ̂θ = rFθ , (1.70)
∂r
d ∂T ∂T
. − = Qr
dt ∂r ∂r
. .2
2 2
. .2
1 2
d ∂ 2 m r +r θ ∂ 21 m r +r2 θ
. − = Fr
dt ∂r ∂r
d . .2
(m r) − mr θ = Fr
dt
.. .2
m r −mr θ = Fr (1.71)
17
ción es la asociada a la coordenada θ, tenemos que
d ∂T ∂T
. − = Qθ
dt ∂ θ ∂θ
. .
2 2 2 2
. .2
1 2 1
d 2∂ m r +r θ ∂ 2 m r +r θ
. − = Fθ
dt ∂θ ∂θ
d .
(mr2 θ) = Fθ
dt
..
2 2 .
mr θ +2mr r= rFθ (1.72)
diente x, la posición dela otra masa esta determinado por la constricción de que la
longitud de la cuerda es siempre la misma (no se estira). El potencial es
La energı́a cinética es
1 .2
T = (M1 + M2 ) x . (1.74)
2
1 .2
L = T − V = (M1 + M2 ) x +M1 gx + M2 g(l − x). (1.75)
2
18
Esto implica que solo hay una ecuación de movimiento,
d ∂L ∂L
. − =0
dt ∂ x ∂x
d .
(M1 + M2 ) x +(M2 − M1 )g = 0
dt
..
(M1 + M2 ) x= (M1 − M2 )g (1.76)
o
.. M1 − M2
x= g (1.77)
M1 + M2
3. Una cuenta que se desliza sobre una alambre que gira uniformemente libre de fuer-
za (figura 7). Este ejemplo es una ilustración simple de constricciones dependientes
del tiempo, con la rotación entorno al eje z y el alambre se encuentra en el plano xy.
Las ecuaciones de tranformación contienen una dependencia explicita del tiempo
1 .2
T = r +r2 ω 2 (1.80)
2
19
La ecuación de movimiento es entonces
d ∂T ∂T
. − =0
. dx
.∂ r ∂r
1 2 2
d ∂ 2 m r +r ω
2 2 ∂ 12 m r +r2 ω 2
. − =0
dx ∂r ∂r
d .
(m r) − mrω 2 = 0
dx
..
m r −mrω 2 = 0
..
r= rω 2 (1.81)
esta ecuación es parecida a la del oscilador armónico pero con el signo cambiado. La
solución es r = r0 eωr muestra que la cuenta se mueve exponencialmente hacia afuera
debido a la aceleración centrı́peta. El momento angular, L = mr02 ω 2 eωt , proviene de la
fuerza F = mr0 ω 2 eωt , que actúa perpendicularmente al alambre y al eje de rotación.
20