UNIDADIV
UNIDADIV
opavtor91@[Link]
Semestre 2021-B
Consideremos sólo sistemas para los cuales el trabajo virtual total de las fuerzas de
ligadura es cero. Por lo tanto, tenemos
X (a)
Fi · δri = 0
i
A esta relación se le conoce como principio de trabajo virtual. Como los coeficientes
(a)
de δri no pueden seguir siendo iguales a cero; i.e., en general Fi 6= 0, ya que δri
no son completamente independientes pero están conectadas por ligaduras. Con el
fin de igualar los coeficientes a cero, debemos transformar el principio en una forma
que involucre los desplazamientos virtuales de qi , las cuales son independientes.
Nuestro caso satisface que no existe una dependencia de fi , pero sólo está rela-
cionado con la estática; queremos un caso general del movimiento del sistema.
Para hacerlo debemos considerar el siguiente procedimiento originalmente ideado
por James Bernoulli y desarrollado por D’Alembert. Partimos de
Fi = ṗi
o
Fi − ṗi = 0
Lo cual quiere decir que el sistema está en equilibrio respecto a una fuerza igual a
la fuerza actual, pero en sentido contrario
X
(Fi − ṗi ) · δri = 0
i
ri = ri (q1 , q2 , ..., qn , t)
La velocidad vi puede ser expresada de la siguiente manera
dri X ∂ri ∂ri
vi = = q̇k +
dt ∂qk ∂t
k
Los desplazamientos virtuales arbitrarios δri pueden ser conectados con los desplaza-
mientos virtuales δqi por
X ∂ri
δri = δqj
∂qj
j
X ∂ri X d ∂ri
d ∂ri
mi r̈i · = mi ṙi · − mi ṙi ·
∂qj dt ∂qj dt ∂qj
i i
Tenemos
( # )
X d ∂T ∂T
− − Qj δqj = 0
dt ∂ q̇j ∂qj
Entonces tenemos que
d ∂T ∂T
− = Qj
dt ∂ q̇j ∂qj
1
Las fuerzas generalizadas pueden escribirse como
X ∂ri X ∂ri ∂V
Qj = Fi · =− ∇i V · =−
∂qj ∂qj ∂qj
i i
−V (r1 , r2 , ..., rN , t)
respecto a qj
Entonces tenemos que
d ∂T ∂(T − V )
− =0
dt ∂ q̇j ∂qj
1 Goldstein, Poole and Safko, Classical dynamics, tercera edición (en inglés) página 21,
L=T −V
entonces podemos escribir la ecuación de la siguiente manera
d ∂L ∂L
− =0
dt ∂ q̇j ∂qj
dL X d ∂L
= q̇i
dt dt ∂ q̇i
i
Para ello usamos el teorema de Euler de las funciones homogéneas, el cual nos dice
que:
Una función f (x, y , z) es una función homogénea de grado k si
Reemplanzando el valor en
X ∂L
E= q̇i −L
∂ q̇i
i
obtenemos
∂(T − U)
= mẋi = pi = constante
∂ ẋi
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 19 / 72
Conservación del momento angular
Si el sistema rota alrededor de cierto eje un ángulo infinitesimal δθ, el radio vector
r a un punto dado cambia una cantidad r + δr, donde
δr = δθ × r
Los vectores de velocidad también cambian bajo la rotación del sistema, y ya que la
ecuación de transformación para todos los vectores es la misma, entocnes tenemos
δṙ = δθ × ṙ
Consideremos una sola partı́cula y expresamos el lagrangiano en coordenadas rect-
angulares. El cambio en L causado por la rotación infinitesimal es
X ∂L ∂L
δL = δxi + δ ẋi = 0
∂xi ∂ ẋi
i
Sabemos que
∂L
pi =
∂ ẋi
por lo tanto las ecuaciones de Lagrange pueden ser expresadas como
∂L
ṗi =
∂xi
Por lo tanto tenemos que
X
δL = ṗi δxi + pi δ ẋi = 0
i
o
ṗ · δr + p · δṙ = 0
δθ · (r × ṗ) + δθ · (ṙ × p) = 0
δθ · [(r × ṗ) + (ṙ × p)] = 0
Podemos ver fácilmente que
d
δθ ·
(r × p) = 0
dt
como δθ es arbitrario entonces tenemos que
r × p = constante
donde el producto vectorial representa al momento angular.
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 23 / 72
Ejemplo 1
Ejemplo 1. Una masa se mueven en una dimensión de tal manera que su la-
grangiano tiene la forma
m2 ẋ 4
+ mẋ 2 V (x) − V 2 (x)
L=
12
donde V es una función diferenciable de x. Encontrar la ecuación de movimiento
para x(t) y describir la naturaleza fı́sica del sistema sobre la base de esta ecuación.
Solución
∂L dV dV
= mẋ 2 − 2V (x)
∂x dx dx
∂L m2 ẋ 3
= + 2mẋV (x)
∂ ẋ 3
d ∂L d
= m2 ẋ 2 ẍ + 2mẍV (x) + 2mẋ V (x)
dt ∂ ẋ dt
ẋ 2
dV
mẍ + (m + V (x)) = 0
dx 2
dV
donde podemos identifiar los siguientes términos F = − y T = mẋ 2 /2, podemos
dx
escribir lo anterior como
(F − mẍ)(T + V ) = 0
La que nos quiere decir la expresión anterior, es que en cualquier tiempo, ya sea la
diferencia F y ma es cero o la suma de la energı́a cinética más la energı́a potencial
es cero.
Ejemplo 2. Un pendulo doble está atado a un carro de masa 2m que se mueve sin
fricción sobre una superficie horizontal. Cada pendulo tiene una longitud b y masa
m. Encontrar las ecuaciones de movimiento.
x1 = bcosθ1 + x; x2 = x1 + bcosθ2
y1 = bsenθ1 ; y2 = y1 + bsenθ2
La energı́a cinética está dada por
2
2m 2 1 1
T = ẋ + m b 2 θ̇12 + ẋ 2 + 2b θ̇1 ẋcosθ1 + m ẋ + b θ̇1 cosθ1 + b θ̇2 cosθ2
2 2 2
2
1
+ m b θ̇1 sinθ1 + b θ̇2 sinθ2
2
La energı́a potencial es
El Lagrangiano es
1
L = T − U = 2mẋ 2 + mb 2 θ̇12 + 2mb θ̇1 ẋcosθ1 + mb 2 θ̇22
2
∂L d ∂L
=
∂θ1 dt ∂ θ̇1
∂L d ∂L
=
∂θ2 dt ∂ θ̇2
Solución
Para m1 :
x1 = acosθ
y1 = rcosΩt = asenθcosΩt
z1 = −rsenΩt = −asenθsenΩt
Para m2 :
x2 = acosθ
y2 = −rcosΩt = −asenθcosΩt
z2 = rsenΩt = asenθsenΩt
donde r = asenθ
xM = 2acosθ
yM = 0
zM = 0
La velocidad para m1 es:
ẋ1 = −aθ̇senθ
ẏ1 = aθ̇cosθcosΩt − aΩsenθsenΩt
ż1 = −(aθ̇cosθsenΩt + aΩsenθcosΩt)
La velocidad para m2 es:
ẋ2 = −aθ̇senθ
ẏ2 = −aθ̇cosθcosΩt + aΩsenθsenΩt
ż2 = aθ̇cosθsenΩt + aΩsenθcosΩt
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 33 / 72
Ejemplo 3
ẋM = −2aθ̇senθ
ẏM = żM = 0
La expresión para el Lagrangiano es L = T − U, la energı́a cinética T es
1 1 1
T = m1 (ẋ12 + ẏ12 + ż12 ) + m2 (ẋ22 + ẏ22 + ż22 ) + M(ẋM
2 2
+ ẏM 2
+ żM )
2 2 2
ya que m1 = m2 obtenemos
1 1
T = m(ẋ12 + ẏ12 + ż12 + ẋ22 + ẏ22 + ż22 ) + M(ẋM
2 2
+ ẏM 2
+ żM )
2 2
x1 = x + Lsenφ
y = −Lcosφ
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 36 / 72
Ejemplo 4
L = Lm + LM
1 2 1 2 1
L= mẋ − kx + M[L2 φ̇2 + 2Lẋcosφφ̇ + ẋ 2 ] + MgLcosφ
2 2 2
1 1 1
L= (m + M)ẋ 2 + M[L2 φ̇2 + 2Lẋcosφφ̇] − kx 2 + MgLcosφ
2 2 2
senφ ≈ φ
cosφ ≈ 1
Reemplazando en las ecuaciones para x y φ, obtenidas con anterioridad, tenemos:
U = U(r) = U(x, y , z)
El lagrangiano es
L=T −U
∂L ∂U
=− = Fx
∂x ∂x
y
∂L ∂T
= = mẋ = px
∂ ẋ ∂ ẋ
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 40 / 72
Ecuaciones de Lagrange para un movimiento sin
rectricciones
Derivando de nuevo la segunda ecuación respecto al tiempo y recordando la segunda
ley de Newton Fx = ṗx vemos que
∂L d ∂L
=
∂x dt ∂ ẋ
Podemos obtener un resultado similar para y y z. Por lo tanto, podemos demostrar
que la segunda ley de Newton implica que:
∂L d ∂L ∂L d ∂L ∂L d ∂L
= ; = ; =
∂x dt ∂ ẋ ∂y dt ∂ ẏ ∂z dt ∂ ż
Para una sola partı́cula en coordenadas rectangulares la segunda ley de Newton es
equivalente a las tres ecuaciones de Lagrange.
La trayectoria de la partı́cula determinada por la segunda ley de Newton es la misma
que la determinada por las tres ecuaciones de Lagrange. Podemos reconocer que
las tres ecuaciones de Lagrange tienen la misma forma de las ecuaciones de Euler-
Lagrange. Lo cual implica que la integral S es estacionaria para la trayectoria que
sigue la partı́cula.
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 41 / 72
Equivalencia de las formulaciones de Lagrange y Newton
∂(T − U) d ∂(T − U)
− = 0, i = 1, 2, 3
∂xi dt ∂ ẋi
En coordenadas rectangulares y para un sistema conservativo tenemos que T =
T (ẋi ) y U = U(xi ), entonces tenemos que las ecuaciones de Lagrange son:
∂U d ∂T
− =
∂xi dt ∂ ẋi
entonces por
∂U d ∂T
− =
∂xi dt ∂ ẋi
tenemos que
Fi = ṗi
En general podemos pasar de un sistema expresado en coordenadas rectangulares
a un sistema en coordenadas generalizadas qj .
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 43 / 72
Equivalencia de las formulaciones de Lagrange y Newton
Entonces tenemos
xi = xi (qj , t)
entonces
X ∂xi ∂xi
ẋi = q̇j +
∂qj ∂t
j
entonces
∂ ẋi ∂xi
=
∂ q̇j ∂qj
entonces el momento generalizado pj asociado con qj puede ser fácilmente deter-
minado por
∂T
pj =
∂ q̇j
Encontremos ahora las ecuaciones de Lagrange para una partı́cula que se mueve en
dos dimeniones en coordenadas polares.
x = rcosθ
y = rsenθ
Podemos encontrar los componentes de la velocidad, como lo hicimos anteriormente
entonces tenemos que
v 2 = ẋ 2 + ẏ 2 = ṙ 2 + r 2 φ̇2
Por lo tanto, la energı́a cinética es:
1
T = m(ṙ 2 + r 2 φ̇2 )
2
F = m(r̈ − r φ̇2 )
fk (qj , t) = 0, k = 1, 2, 3, ...
Si las ligaduras pueden ser expresadas como se mencionó en la condición 2, estas
se llaman ligaduras holónomicas u holónomas . Si las ecuacioens no contienen
explicitamente al tiempo, las ligaduras se dice que son fijas o esclerónomas; las
ligaduras que se mueven son llamadas reónomas2 .
Las ligaduras que pueden ser expresadas como relaciones algebraicas entre las coor-
denadas son ligaduras holónomicas. Si un sistema sólo está sujeto a tales ligaduras,
podemos econtrar siempre un conjunto apropiado de coordenadas generalizadas a
modo que las ecuaciones de Lagrange sean libres de una referencia explicita a estas
ligaduras.
Cualquier ligadura que deba ser expresada en términos de las velocidades de las
partı́culas en el sistema son de la forma:
f (qj , q̇j , t) = 0
y constituyen ligaduras no holónomicas a menos que la ecuación pueda ser integrada
para darnos las relaciones entre las coordenadas3 .
lo cual es
df
=0
dt
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 53 / 72
Ecuaciones de Lagrange con multiplicadores
indeterminados
f (qj , t) − constante = 0
por lo cuál tenemos que la ligadura es holónomica.
Podemos concluir que las ligaduras expresables en la forma diferencial
X ∂f ∂f
q̇j + dt = 0
∂qj ∂t
i
Ejemplo 1. Una partı́cula de masa m parte del reposo desde la cima de un hemis-
ferio liso fijo de radio a. Encuentre la fuerza de ligadura, y determine el ángulo en
el cual la partı́cula deja el hemisferio.
Solución
Para nuestro caso las coordenadas generalizadas son r y θ. La ecuación de la
restricción (o ligadura) es
f (r , θ) = r − a = 0
El lagrangiano está determinado por las energı́a cinética y potencial:
m 2
(ṙ + r 2 θ̇2 )
T =
2
U = mgrcosθ
m
L = T − U = (ṙ 2 + r 2 θ̇2 ) − mgrcosθ
2
La energı́a potencial es cero cuando está en el fondo del hemisferio.
Reemplazando r = a, ṙ = r̈ = 0
maθ̇2 − mgcosθ + λ = 0
magasinθ − ma2 θ̈ = 0
de la segunda ecuación obtenemos que
g
θ̈ = sinθ
a
Realizando el siguiente ansatz
d dθ d θ̇ d θ̇ dθ d θ̇
θ̈ = = = = θ̇
dt dt dt dθ dt dθ
Z Z
g
θ̇d θ̇ = sinθdθ
a
entonces
θ̇2 −g g
= cosθ +
2 a a
donde la constante de integración es g /a, ya que θ̇ = 0 en t = 0, después de
resolver para λ. Sustituyendo en
maθ̇2 − mgcosθ + λ = 0
obtenemos
λ = mg (3cosθ − 2)
la cual es la fuerza de ligadura.
λ = mg (3cosθ − 2)
0 = mg (3cosθ0 − 2)
despejando tenemos
2
θ0 = cos −1
3
r = (rcosθ, rsinθ, z)
y la velocidad es
f (r , θ, z) = −r 2 − αz = 0
−2rdr + αdz = 0
Las ecuaciones de Lagrange son
d ∂L ∂L
− = Qi
dt ∂ q̇i ∂qi
donde Qi son las fuerzas generalizadas de restricción.
mr 2 θ̇ = constante = J
r2
usando la ecuación de restricción z = obtenemos
α
2r ṙ 2rr̈ 2ṙ 2
ż = , z̈ = +
α α α
Podemos escribir la expresión para la energı́a total como
1
E= m(ṙ 2 + r 2 θ̇2 + ż 2 ) + mgz
2
−1
J2 2gr 2 4r 2
2E
ṙ 2 = − 2 2− 1+
m m r α α2
entonces de la ecuación mz̈ + mg = αλ obtenemos
m
(2rr̈ + 2ṙ 2 ) + mg = αλ
α
y usando las ecuaciones restantes obtenemos
4r 2 2mṙ 2 2J 2
αλ 1 + 2 = + mg +
α α mαr 2
Fácilmente podemos obtener la siguiente expresión
−2 −2
8J 2 4r 2 4r 2
4E mg
λ= 2
+ + 1+ = constante · 1 +
α mα4 α α α
r − 3
4r 2 4r 2
2
f = αλ 1 + 2 = constante · 1 +
α α
Esta fuerza se encuentra en el plano rz y es perpendicular al interior de la superficie
del paraboloide. (La cual forma un ángulo arctan(−α/2r ) con el eje r mientra que
la pendiente de la parábola es 2r /α.
Ejemplo 3. Una masa m1 está situada sobre una mesa horizontal sin rozamiento.
Enganchado a ella hay una cuerda que corre horizontalmente hasta el borde de la
mesa, donde pasa por una polea pequeña sin rozamiento. Al final de la cuerda cuelga
una masa m2 . Usar como coordenadas x e y las distancias de m1 y m2 desde la
polea. Estas cumplen la ecuación de ligadura f (x, y ) = x + y = constante. Escribir
las dos ecuaciones de Lagrange modificadas y resolverlas (junto con la ecuación de
∂f
ligadura) para ẍ, ÿ , y el multiplicador de Lagrange λ. Usar la ecuación λ = Fxlig
∂x
(y la ecuación en y correspondiente) para encontrar las tensiones sobre las dos
masas. Verificar la respuesta resolviendo el problema mediante la aproximación
Newtoniana elemental.
Solución
Para la masa 1 (m1 ) tenemos:
1
T = m1 ẋ 2
2
λ = m1 ẍ
−m2 g + λ = m2 ÿ
despejando λ = m2 ÿ + m2 g e igualando ambas expresiones de λ tenemos:
m1 ẍ = m2 ÿ + m2 g
m1 ẍ − m2 ÿ = m2 g
Podemos notar que:
x +y =l
ẍ = −ÿ
entonces
m1 ẍ + m2 ẍ = m2 g
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 69 / 72
Ejemplo 3
m2 g
ẍ =
m1 + m2
m2 g
a=
m1 + m2
Usando la formulación de Newton
Para m1 tenemos:
m1 ax = T
N − m1 g = 0
Dr. en C. Omar Pavón Torres (UAEM) Ecuaciones de Lagrange Semestre 2021-B 70 / 72
Ejemplo 3
Para m2 tenemos:
T − m2 g = may
Despejando T de alguna de las expresiones y sustituyéndola en la otra obtenemos:
m1 ax − m2 g = may
m1 ax − m2 ay = m2 g
Sabemos que ax = −ay , entonces tenemos que
m1 ax + m2 ax = m2 g
Despejando ax = a, tenemos
m2 g
a=
m1 + m2