PROBLEMAS

DE ME C ÁNI CA
Víctor M. Pérez García
Luis Vázquez Martínez y
Antonio Fernández-Rañada

A lin n 7 n F r litn r in l
 V íctor M . P érez G arcía,
Luis V ázq u ez M artínez
y A n to n io F ernández-R añada

100 problemas
de Mecánica

A lia n za E ditorial
 A Luciana, Mary y Conchi

R e serv ad o s to d o s los derechos. E l co n te n id o de esta o b ra está p ro teg id o p o r la L ey, q u e establece

p en as d e prisión y/o m u ltas, ad em ás d e las co rre sp o n d ien tes indem nizaciones p o r d añ o s y perjuicios,
p a ra q u ien es re p ro d u je re n , plagiaren, d istrib u y eren o co m u n icaren p ú b licam en te, e n to d o o e n p arte,
u n a o b ra literaria, artística o científica, o su tran sfo rm ació n , in te rp re ta c ió n o ejecu ció n artística fijada
e n cu alq u ier tip o d e so p o rte o co m u n icad a a través de cu alq u ier m ed io , sin la p recep tiv a autorización.

© V íc to r M . P é re z G a rc ía , L uis V á z q u e z M a rtín e z y A n to n io F e rn á n d e z -R a ñ a d a
© A lian za E d ito ria l, S. A ., M a d rid , 1997
J u a n Ig n acio L u ca d e T e n a , 15; teléf. 393 88 88; 28027 M a d rid
ISBM^fSSSSSSif'.-
DepóW toT égal: M; 37.940-1997 1
F o to co m p o sició n : E F C A , S. A .
P a rq u e In d u stria l «L as M onjas»
28850 T o rre jó n d e A rd o z (M ad rid )
Im p re so e n G ráfic as A N Z O S , S. A . F u e n la b ra d a (M ad rid )
P rin te d in Spain
 índice

Problemas Páginas

I n t r o d u c c i ó n ...................................................................................... XI
1. P a r t íc u l a e n u n a d im e n s ió n ........................................... 1-11 1
i P a r t íc u l a e n d o s y tr e s d im e n sio n e s . P o t e n c ia l e s
c e n t r a l e s ............... ................................................................... 12-21 25
3. C o l i s io n e s .................................................................................. 22-27 53
4. F o r m u l a c ió n l a g r a n g ia n a d e l a m e c á n ic a .......... .28-37 67
5. M o v im ie n t o p l a n e t a r io .................................................... 38-47 91
6. T e o r ía d e l p o t e n c ia l .......................................................... 48-56 107
1. O s c i l a c io n e s ............................................................................ 57-66 125
8. O s c il a c io n e s p e q u e ñ a s e n s is t e m a s c o n v a r io s
g r a d o s d e l ib e r t a d .............................................................. 67-75 147
V. S is t e m a s d e r e f e r e n c ia n o in e r c ia l e s ....................... 76-80 175
10. S ó l id o R íg id o (I): M o v im ie n t o p l a n o ........................ 81-94 185
11. S ó l id o R íg id o (II): M o v im ie n t o g e n e r a l ................ 95-100 223

IX
 introducción

Este libro recoge algunos de los problemas que propusimos a nuestros alum­
nos de la asignatura de «Mecánica y Ondas» de segundo curso de Ciencias Físicas
de la Universidad Complutense a lo largo de más de diez años. Tratan sólo de
Mecánica Newtoniana, sin Ondas, Relatividad ni Fluidos, cuya inclusión alargaría
el texto más allá de las dimensiones de esta colección «100 problemas de...». Al
elegir los que aquí figuran, se ha cuidado mucho que el libro sea adecuado para
la enseñanza de la Mecánica en los planes de estudio actualmente en vigor.
Resolver problemas es imprescindible en el aprendizaje de cualquier ciencia;
sin duda desarrolla en el estudiante la actitud necesaria para enfrentarse creati­
vamente con las dificultades de su oficio futuro. Además, es una tarea enrique-
cedora que ayuda a practicar métodos de razonamiento válidos en todos los ór­
denes de la vida. Hasta el punto que, en sus esfuerzos por entender mejor el
funcionamiento de la mente, los psicólogos estudian los procesos mentales que
llevan a los científicos a resolver problemas.
Comparando la forma de resolver problemas de un estudiante primerizo con
la de otro más veterano o con la de un profesional, se observan enseguida dos
diferencias. El novicio suele lanzarse a hacer los cálculos inmediatamente tras
la primera idea y abandona el trabajo nada más llegar a un resultado. En cam­
bio el que ya está fogueado sabe muy bien que las primeras ideas no suelen ser
las mejores, que es peligroso perder tiempo y energías haciendo cálculos sin
asegurarse de que son necesarios y que ningún resultado debe ser aceptado sin
comprobarlo luego. Por ello, es muy conveniente seguir una táctica dividida en
las cuatro fases siguientes:

(i) Comprender. Al principio, se debe leer cuidadosamente el enunciado

para entender bien qué se pregunta y cuáles son los datos. Hay que presentarse

X!
el problema en términos intuitivos, hacer diagramas y dibujos, intentar una esti­
mación del resultado..., pero sin hacer todavía ningún cálculo.
1
(ii) Formarse un plan. U na vez comprendido el problema en términos intuiti­
Partícula en una dim ensión
vos, hay que elaborar un plan de trabajo, decidiendo qué conceptos, leyes o
ecuaciones se van a usar.
(iii) Calcular y resolver. Se escriben las ecuaciones elegidas en la segunda
fase y se resuelven.
(iv) Comprobar e interpretar. Los resultados obtenidos en la tercera fase de­
ben compararse con las intuiciones de la primera, comprobando que su orden
de magnitud sea adecuado, que las dimensiones sean correctas... Si se descubre
alguna contradicción o simplemente algo que no se entiende bien, hay que pen­
sar más en ello. Es muy importante comparar las intuiciones de la primera fase
con los resultados concretos de la tercera: sin hacerlo, no se puede considerar
terminado un problema.
Problem a 1
Mientras enseñábamos Mecánica y Ondas en la Complutense debatimos a
menudo sobre el mejor modo de hacerlo con varios colegas a quienes agradece­
mos sus aportaciones. Se trata de Roberto González Amado, Ángel Fernández D eterm inar las funcio n e s energía potencial de las
Ortega, Juan Usón, Luis García Gonzalo, Carmen Ballesteros, Amador Álva- que se pueden d e riva r las fuerzas
rez, Francisco Domínguez-Adame y Ángel Sánchez. Nuestro agradecimiento se
extiende a muchos estudiantes con quienes hemos mantenido animadas discu­ 1. F,(x) = a x + b
siones de las que hemos aprendido nosotros también y a José Luis Sánchez Gó­ 2. F¿x) = F0 exp (-¡x|/a), a > 0
mez por invitarnos a publicar esta colección de problemas. Por último, quere­
3. F,(x) = F0 sen (kx)
mos invitar a los lectores a que nos hagan llegar sus sugerencias o comentarios
a la dirección internet: v perez@ ind-cr.uclm.es a \ í +1' Sl s e n (7 rx )> 0
4. FA(x) =
[-1 , si s e n (7 rx )> 0
Bibliografía. Un libro de problemas debe usarse siempre junto con uno de te­
oría. Hay muchos libros de mecánica adecuados para ello; a título de ejemplo
mencionamos los tres siguientes:

A n t o n i o R a ñ a d a , Dinámica clásica, Alianza Universidad Textos, Madrid,

Para obtener el potencial en todos los casos basta in­
1990 (reimpreso en 1994).
tegrar ya que en el caso de movimiento unidimensional
H . G o l d s t e i n , Mecánica clásica, 2.a edición, Reverte, Barcelona, 1987.
P. F r e n c h , Mecánica newtoniana, Reverté, Barcelona, 1974.
F (x )= - ~ ~ = > j d V ( x ) F(x)d x= $ V ( x ) - ~ V ( x 0) = -J F(x)d x
El primero es el que hemos usado los autores al preparar este libro de pro­ C¿JC
blemas.
1.
r üx
Fl(x) =ax + b=>V1(x) =-](a x + b)dx= — -— bx + C
 Observemos que V{2) =- V(0), por el carácter periódi­
co del sistema. Reuniendo los dos resultados obtenemos
que

Consideremos por separado los casos x > 0 y r < 0 e

integremos a partir de cero V(x) = í V(0) + X’ X 6 [° ’ 11
\ v ( 0) + 2 - x , x e ( l,2 )
x > 0, V2(x) - V2(0 ) = -F0 J exp (~xla) dx = F0\a\ (e-*la + X)
Fuera del intervalo [0, 2] el potencial se repite de ma­
nera periódica.
x < 0 , V2(x) - V2(0) - -F 0 J exp (x/a) dx - ~F0\a\ (e*to- 1)

Y uniendo los dos resultados en una sola expresión

encontramos Problem a 2

H allar las frecuencias de las oscilaciones peque­

Vt (x) = -sgn (x) F^a (e~Wa - 1 ) + Vz(0)
ñas de una p a rtícu la de m asa m a lre d e d o r de las
posiciones de e q u ilib rio de los sigu ie n te s potencia­
3. les:
F
F3(x) - F0 sen fáx => V3(x) = eos kx
K
1. V lW - f f j ¿ + 1)
4. Vamos a analizar los dominios de definición
2. V2{z) = Kcoshíz}
sen ttx > 0 =>x e [0,1] U [2,3] U [3,4]... =$x e 0 [2/c, 2k + 1]
k=0 3. V2{z) = /C (z 2 - 4z - 21)
4. l/{z) = /C(z3- 1 0 z 2 -i-21z+15)
sen 77%<0 =>x e (1 ,2)U (3,4)U (5.6) ... =$X <=0 (2k~l,2k)
k=0
donde z = x/a, siendo a > 0 y K constantes con di­
Como la función F es periódica en [0, 2] también lo m ensiones de lo n g itu d y energía respectivam ente.
será el potencial. Basta resolver el problema en este in­
tervalo.

(a) Cuando x < l , S olución

V(x) - V(0) = f dx = x 1. E ste potencial no tiene mínimos con lo cual no
jü
hay puntos de equilibrio estable.
(b) En la región 2 > x > 1, 2. La derivada del potencial V{x) - cosh (x/a) es
K
V'(x) = — senh (x/a). Imponiendo V'(x) = 0 se obtie­
 ne el único punto de equilibrio que es x = 0. Para estu­ Evaluemos la derivada segunda en estos puntos de
diar el carácter del mismo calculemos equilibrio

(P-V K L/x tfV 2K . . in .

= —r cosh --- ~dx2
T T = ~a*
ñ i ^ x ~ 10a)
dx1

y como V"(0) = Kla1 > 0 el punto de equilibrio es un mí­ Entonces V "(xJ = 2 \J Y ¡k ía 2 > 0 y V"(x _) < 0 ..D e
nimo del potencial y por lo tanto un punto de equilibrio donde se deduce que el único mínimo está en x+ y en
estable. En cuanto a la frecuencia de las oscilaciones pe­ torno a él habrá oscilaciones de frecuencia
queñas en torno a este punto de equilibrio, vendrá dada
por las ecuaciones w=lT W . 1 W 2V 37 k
m a V m
V"(*)
m
P roblem a 3
3. De nuevo repetimos el proceso de los casos ante­
riores, pero ahora con V(x) = K[(x/a )2 - 4(x/n) - 21]. su
H a lla r el p e rio d o d e o scila ció n d e una p a rtícu la
derivada es V'(x) - [2x1a2 - 4la] e imponiendo que se
som etida a un potencial de! tip o V (X ) = B \x\h. Par-
anule se encuentra como única solución x = 2. Como
1
V"(2) = 2K/az > 0 el punto de equilibrio es un mínimo. ticu la riza r para el caso B * — k, n = 2.
S olución La frecuencia de las oscilaciones pequeñas será enton­
ces

V"_Q) 2K
S olución
m m
Utilizaremos la conservación de la energía para redu­
cir a cuadraturas el problema
4. V{x) = K [(xlá)3 - 10 (xla )2 + 21 (x/a) + 15]
E = ~ m x2 + B 1*1"
3x2 20x 21
V \x ) = K = ~ [3x2 - 20ax + 21a2] = 0 De aquí
a3 a2 + a

Resolviendo la ecuación de segundo grado se obtie­

nen las raíces

x+ = a (10 + V 37)/3 = 1.3057 a

x_ = a (10 - V37)/3 = 5.3609 a
 Donde se ha hecho uso de que, por la simetría del po­ y entonces el periodo queda
tencial, el tiempo que se tarda en ir y volver al mismo
punto (el periodo) es cuatro veces el tiempo que se tar­
da en ir desde el punto de mayor elongación hasta el
mínimo del potencial (en x - 0). Además cuando la par­
tícula se encuentra en el punto de máxima elongación 2 m „ r í± ,I
debe estar en reposo ya que toda la energía es poten­ l n 2
cial. En este punto E = V (xmax) = B\xmax\n y despejando
se obtiene* ntax = \ /EIB. Particularizando para el caso n = 2 tendremos
n _____
Hagamos el cambio de variables x = u VEIB. Enton­
ces um ax = 1 Jy
r= 4 V — =
V 2 Vb r(i) Yk

Así que en este caso especial el periodo de las oscila­

ciones no depende de la amplitud de las mismas sino
2 V b V e jq V i - u n
sólo de las características del sistema (k y m).

Hemos separado pues la dependencia del periodo en

la amplitud y en la energía. Aunque no supiéramos cal­
cular la integral que nos queda ya tendríamos mucha in­
formación sobre el mismo y la única incógnita sería la
E stu d ia r c u a lita tiv a m e n te el m o v im ie n to de un
función de n que define la integral. La integral puede
péndu lo de lo n g itu d l y masa m. D educir de la ex­
calcularse numéricamente para cada n y tabularse pero
presión del potencial V{a) cuántos tip o s de m o v i­
en nuestro caso puede escribirse usando la función T de
m ie n to s d iferente s existen y representa rlos en un
Euler (ver listado de fórmulas)
diagram a a, á.

r(n)=J'
Va
Sds

de la que usaremos sus propiedades de que T (l/ 2 ) = Vrr

y T(l) - 1. La expresión de la integral usando la función
T es

du Vrr F j

7
 La energía del péndulo es la suma de sus energías ci­ Entonces, cuando a = 2mr ]p Herh¿arja segunda es po­
nética y potencial y se puede escribir en función del án­ sitiva y el equilibrio_es estable. Cuando a = (2n + l)-7rla
gulo formado por el péndulo y el eje Y como .derivada segunda-es~negati-va-y-.eh.equilibxÍQ. inestable..
S olución La forma del potencial es la que se representa en la fi­
gura, apareciendo tres regiones de com portam iento
T =i m (x2+ y2) = -i- m i1(á2 sen2a +cr eos2 a) = ■
- i- m (á/)2 cualitativamente distinto. En primer lugar la línea conti­
2 2 2
nua nos muestra una solución con E > mgl para la cual
V ~ -m gy - -m gl eos a la partícula se va a infinito. La línea de puntos muestra
una solución con E < mgl para las cuales hay oscilacio­
nes. Por último la línea de trazos cortos se llama separa-
E - ~ nía2? - mgl eos a triz, que corresponde a un m ovim iento con energía
igual a la altura del potencial. Entonces la solución ten­
Donde se ha tomado como origen de potenciales el derá a llegar a la máxima amplitud con energía cinética
punto de suspensión del péndulo ya que cuando a = tt!2 y velocidad nula en tiempo infinito. Por otro lado, des­
entonces V = 0. Los puntos de equilibrio del sistema son pejando de la ecuación de la energía se obtiene la forma
exacta de la relación entre c*y cu, que es
aquellos donde las ftierzas son nulas, o sea que satisfa­
cen dV/da = 0.
En este caso . , Í2 E 2g
a = ± \ / — — + — eos a
V mP l
dV
mgl sen a = 0
da E sta ecuación nos permite pintar las curvas que des­
criben las soluciones en el plano (a. á | que es a lo que
cuyas soluciones son a = m r con n = 0, ±1, ±2, ... La es­ se llama el plano o espacio, de-f ases, del sistema. Este
tabilidad de estos puntos de equilibrio nos la dan las de­ tipo de expresiones suele ser complicada de dibujar y
rivadas segundas del potencial norm alm ente se sigue un -procedim iento cualitativo
.pa.ra-p-intaF-fas-airyas^ea.el,.asp.acio„4£„iase^ Veamos
cPV por ejemplo para los tres tipos de movimiento descritos
- mgl eos a
~d¿ anteriormente.
En primer lugar si la energía es menor que la altura
E > mgl ■
del potencial el movimiento consiste (suponemos la ve­
locidad inicial en dirección positiva) en una aproxima­
ción de la partícula a la barrera con velocidad decre­
ciente hasta que llegue a anularse al llegar al punto
donde toda su energía es potencial. En el plano de fases
esto se corresponde con una curva de altura decreciente
y a creciente hasta llegar al punto en el que a es máxi­
mo y á = 0, que se denomina de punto de retroceso. A
partir de este momento la velocidad se hace negativa y
el ángulo a decrece hasta que se llega al otro punto de
retroceso en el que a se hace mínimo y la velocidad

8
vuelve a anularse. Después la velocidad vuelve a crecer, D e s c rib ir c u a lita tiv a m e n te los tip o s de m o v i­
al igual que a, hasta pasar por el centro, donde la velo­ m iento que puede te n e r una partícula de masa m
cidad empieza a decrecer. Cuando se encuentre la con­ en cada uno de los siguien tes potenciales (expresa­
dición inicial el movimiento se repite y corresponde en­ dos en unidades arbitrarias) Problem a 5
tonces a una trayectoria cerrada en el plano de fases. Si
se aumenta la energía la forma de los ciclos cambiará y 1. v^x) = a x + b
su radio aumentará, pero seguirán correspondiendo al
mismo tipo de movimientos cerrados, esto es, el com­ 2. V2{x) = x2 - 4 x ~ 21
portam iento es cualitativamente el mismo. Cuando la 3. l/(x) = x3 - lOx2 + 21 x-h 15
energía coincide con la altura de la barrera el movi­
miento es cerrado pero no periódico ya que —como se
ha indicado anteriormente— se tarda un tiempo infinito
en recorrer la curva. Por último, cuando la energía está S olución
por encima de la de la barrera las órbitas son no acota­
das y corresponden a movimientos con velocidad alter- 1. En este caso sólo hay dos tipos de movimiento. Si
nantemente creciente y decreciente pero siempre positi­ la energía es tal que la velocidad inicial es positiva, el
va (o negativa) y por consiguiente a es siempre movimiento es un acercamiento a la barrera del poten­
creciente (o decreciente). cial hasta que se llega a ella con velocidad nula, después
En la siguiente figura se presenta el aspecto del plano de esto el movimiento es en dirección negativa y unifor­
de fases para este problema. memente acelerado llevando a la partícula a -<». Cuan­
do v0 < 0 la partícula se va directamente a -o» con movi­
miento uniform em ente acelerado. Si escribimos la
ecuación de la energía

E = — m v2 + V(x) = — m v2 + ax + b

y despejando x como función de v se encuentra

Así pues las órbitas en el espacio de fases son parábo­

las.
Todos los movimientos son no acotados. Esto es razo­
nable si entendemos que el potencial corresponde a una
Nótese que dada la periodicidad del potencial es sufi­ fuerza constante y por lo tanto el movimiento va a ser
ciente con analizar el intervalo a g [- 77, 7r]. uniformemente acelerado.
2. El potencial es una parábola con un mínimo en
x - 2. Todas las soluciones van a ser acotadas y corres-

11
 3.
Los puntos de equilibrio para el potencial V3 son las
raíces de

dV
— = 3x2~ 20* + 21 = 0
dx

ponden a movimientos periódicos. La ley de conserva­

ción de la energía es

E —— mx2 + x2 - 4x - 21 = — m f + {x - 2 )2 - 25

ecuación que puede reescribirse como

x2 ( * - 2)2

(VP^)’ ( V £ * a r
que son x ± = --------------. Para analizar su estabilidad

que es la ecuación de una elipse de semiejes^ / — (E + 25) calculemos la derivada segunda, que es V"(x) = 6x - 20.
y /T í Entonces V"(xJ < 0 y V"(xJ > 0. Las formas posibles de
y V is + 25 centrada en x = 2 . movimiento son las siguientes:
(a) E < E2 - V(*+): Sólo hay movimientos posibles
en la zona de * negativos y corresponden a movimientos
de escape al infinito o de aproximación a la barrera, co­
lisión con la misma y escape al infinito.
(b) E 1 = V(* ) > E > E2: Hay dos posibilidades. Si la
partícula está en * > *_ se realizan oscilaciones periódi­
cas en torno a *+, mientras que si está inicialmente con
* < * el comportamiento es el mismo que en el caso (a).
(c) E > Ey Hay escape a
(d) E ~ E l es la separatriz, que corresponde a movi­
mientos que acaban o empiezan en * = *_ con velocidad
nula.
 El espacio de fases queda definido por las curvas definir una energía, que será una cantidad conservada
en la evolución, aunque ahora no tendrá el mismo signi­
x = ± \ / — (E - x 3 + 10x2 - 21x - 15) y puede dibujarse ficado que la energía de una partícula newtoniana. El
y fTT
cualitativamente a partir de la información recogida en potencial del que deriva la fuerza F(x) = - K/x 3 es
el estudio anterior.

Entonces la energía es

2 2x2

y evaluándola para nuestro caso, en el que x( 0 ) = 0 se

obtiene E ~ -K !2x\ < 0.
Esto nos permite calcular de manera exacta la solu­
ción x(t) haciendo

■4¡)
í 40 dx'
ax frx(t) x'dx'
J^ A /** T-» . 7 /O J.
V 2 E + K/x'2 V K V 2 E x '2/K + 1

K
Usando ahora que E = — —y hagamos el cambio de
Problem a 6 ¿xQ
variable u = x'2!x2 entonces du = 2x'lx2dx'.
A l e stu d ia r la p ropaga ción de haces láser en cie r­
to s m e dios se encuentra que la anchura x del haz r <» du _ xl
láser satisface la ecuación de evolución t=- V x ¡ -x 2(t)
2V k V V i —u 2V k , V k

x = - — ^, K> 0 Cuando x(t) = 0 se tiene el tiempo que tarda en ha­

x3 cerse cero la anchura del haz que es
E nco ntrar la fo rm a de la fu n ció n x ( t ) y el tie m p o
que ta rd a en hacerse cero su p o n ie n d o que x{0) = 0.
' ~Vk

Solución y despejando en t se obtiene la forma de la función x(t)

La ecuación de evolución que rige este fenómeno físi­

co es una ecuación de Newton. Este hecho nos permite
 La extensión de este ejemplo a un caso con un poten­ 2. El barco se para cuando v{t) = 0, o sea que
cial más complicado así como con pérdidas puede verse
en [V. M. Pérez García, M. A. Porras, L. Vázquez, Phys.
e~av" + a — t = 1
Lett. A, 202 (1995) 176]. m p

Despejando t
iflp ó b le m a 7 i
m
Un barco con ve lo cid a d in icia l vQ se ve fre n a d o
p o r una fuerza de rozam iento F = - b e xp (a v).
La distancia recorrida hasta pararse será
1. E ncontrar la so lu ció n de su ecuación de m o ­
v im ie n to .
* = f " v{t)dt = - f *— ln erav" + a — t dt —
2. C alcular el tie m p o y la distancia que tarda en P J0 0 Oí m
pararse si ¡nicia lm en te viaja con una ve lo cid a d vQ.
1 + g*»n a _ t dt
=_r ^ ( - av»+in m
S olución
ab
Hagamos el cambio de variables 5 = 1 + eav" — t,
1. La ecuación de movimiento es m
entonces
b
v= eaV
m me-™" f f ' r ,
X = v„r_ ^ J ln Í ds = Vatp 3 ;— [í ln s - s ] p
p 0p a2b h .................°'p a2b
Separando las variables t y v e integrando
Simplificando se encuentra finalmente
^ í = f*' e"xv, dv
- — . ,=—
1 (e~av»
y - e~av)
m Kn a m
x- = ~ ^b +
Despejando el término en v y tomando logaritmos
Es interesante que el lector intente responder la si­
guiente pregunta: ¿qué le ocurre a la velocidad del bar­
In í e~av" + a — 1 1 = ~av co cuando t -»<»? ¿Por qué?

y finalmente

v(r) = — ln
a Un coche de carreras se m ue ve im p u ls a d o p o r
e-“v«+ a — t
m un m o to r que le p ro p o rc io n a un e m p u je (fuerza)

16
P - c te = D. A dem ás de este e m puje sobre el coche E invirtiendo la relación se llega a
actúa una fuerza de fre n a d o p ro p o rc io n a l al cua­
d rado de su vefocidad. E ncontrar la ve lo cid a d del
coche en fu n c ió n de su posición sobre la pista y en v(x) = \ —■V i - e-lkxUn
fu n c ió n de f. ¿Cuál será la velocida d m á xim a alcan­
zada p o r el coche? S upone r que el coche parte in i­
cialm e nte del reposo.
Problem a 9

Una gota de agua cae en una atm ósfera saturada

de va p o r de agua. D urante la caída el va p o r se va
conden sando en la gota, que aum enta así de ta m a ­
La ecuación de movimiento del coche es ño y masa. E studiar el m o v im ie n to aceptando la hi­
pótesis de que el a u m ento de masa p o r unidad de
mv = - kv 2 + D
tie m p o es p ro p o rcio n a l a la su p e rficie de la gota.
C onsiderar que la gota tie n e fo rm a esférica, den si­
Separando las variables e integrando, teniendo en
dad u n ifo rm e y que in ic ia lm e n te se encuentra en
cuenta que v(0) = 0 reposo.

k_ V D ¡k -v \
ln
m V5Tk + 7 /
S o lución
De aquí es fácil obtener v(t), que, tras realizar algu­
nas manipulaciones algebraicas, resulta ser El enunciado del problema nos dice que el aumento
de masa por unidad de tiempo es proporcional a la su­
1 _ g -iV 4 D k /m 2 perficie de la gota, o sea que
v (t)= V D ¡k
dm
= ctAirr1
Cuando t -» se obtiene la velocidad máxima alcan- dt
zable por el coche, que es v = \^ D !k . Para encontrar
v(x) lo mejor es proceder directamente de la ecuación Suponiendo que la gota tiene densidad uniforme la
de movimiento, que reescribiremos en la forma masa y el radio de la gota estarán relacionados por
4 / 3m \ 1/3
la ecuación m - pV = p — rrr3 y de aquí r - I- — , de
, dv 47 rpj
D _ kv2 = mv —— donde
dx

Separando las variables e integrando dm I 3 m \ 212

= a47 t j I - <x(3Ó7t)1/3
dt 1 47t p j I p

19
 Esta ecuación diferencial de primer orden pueden in­
tegrarse de manera sencilla

r m dm r‘ ( 36tt \ 1/3
Problem a 10

y de aquí Un cohete quem a su co m b u stib le a ritm o cons-

dm . _
ta n te = - f. E n c o n tra r su p o s ic ió n al ca b o de
36 tT
m m(t) = t + m.
W í? un tie m p o t despreciando el ca m b io de la fuerza de
gravedad con la altura.
/3 6 tt\ 1/3 , , ,
Definamos las constantes ¡x - a YP = (m 0//x) 1/3
en función de las cuales la expresión para la masa toma Solución
una forma más compacta
La masa del cohete será m{t) = m 0 - ft. La ecuación
m(í) = /¿(r + de movimiento para el cohete es

dv vg dm fv
Llevando este valor de la m asa a la ecuación de
Newton dt m dt ^ m

d{mv) Integrando
= F = mg
dt
/\ r r df ,/ m
v(t) = vQ- g t + fv i — =Vo- g t + ve ln
o mü- f f e
e integrando con v(0 ) = 0 se obtiene
De aquí, volviendo a integrar se encuentra la posición
m(t)v(t) = g f m ( f ) d t = - jm [ ( í + &)* ~ £ 4] en función del tiempo
Jo 4

Explicitando m{t) es posible despejar v(r) obtenién­

dose que

La posición en función de t se obtendrá volviendo a

integrar y el resultado es
 1 v donde m es función de t y v, es la velocidad de escape
= V - ^ g í + y - (m0 - / r ) l n A - A j + /í de los gases.
Eliminando explícitamente el tiempo, agrupando va­
riables iguales e integrando se obtiene
m
(v0 + Ve) (m 0 - m ) + m ve ln

Av = v - v 0 = v ln-^4-
m

1. Cohete de una etapa:

La masa inicial del cohete es M = M l + M2 + m 1 + m 2
Un cohete consta de dos etapas de masas M 1 y M 2 y la masa final por conveniencia la escribiremos como
que contienen adem ás unas cantidades de com bus­ Afj + M v donde M x son los motores de la única etapa
tib le m, y m r Analizar, sin tener en cuenta el efecto del cohete y M 2 el peso muerto. Así pues, el incremento
de la gravedad, cuál será la velocidad final del cohe­ de velocidad alcanzado después de haberse consumido
te y com pararla con la que se obtendría usando un el combustible es
cohete de una sola fase. ¿Qué relación deben satis­
facer las m asas de las fases y el co m b u stib le para
itL . -i -r rrt. T ./r * .
que sea favo ra ble el uso de dos etapas? Av = v ln — í-
Mx+ M2

m 2
2. Cohete de dos etapas:
El cohete parte con una masa inicial M = M x + M 2 +
m x + m 2 que al acabar la primera fase es m = M x + M 2 +
A/, + M 2
m2, el incremento de velocidad obtenido en esta etapa
M\ es

. , M, + M . + m. + m.
Av. = v ln — í ^------1------ 2
1 e AL + AL + wu
Cohete de u n a etapa Cohete d e dos etapas

Al acabarse el combustible de la primera fase se des­

prende el peso muerto M v La masa inicial al encender­
se los cohetes de la segunda fase es m = M 2 + m2 y al
S olución acabarse el combustible la masa es Mv El incremento
de velocidad proporcionado por la segunda fase es
Al no actuar la gravedad la ecuación de movimiento
del cohete es sencilla
. , M . + m.
dv dm
-m - r - = v
dt dt

oo
y entonces el incremento total de la velocidad es
2
Av = AVj + Av2 = ve ln
(Mj + M 2 + m l + m 2) (M2 + m 2) Partícula en dos y tres dim ensiones.
M2(M l + M2 + m2) Potenciales centrales
Si comparamos el cohete de una etapa con el de dos

(M2 + m 2) (M l + M J
R = A V 2 e ta p a s " A v i «apa = In
MAM. + M 2 + m2)

El uso de dos etapas será energéticamente rentable si...

R > 0, condición que es equivalente a exigir

{Ml l m ¿ (M l l M ¿ > i
M 2(M1 + M 2 + m 2)
Problem a 12J
Para facilitar los cálculos hagamos M { = aM, M2 = M,
m 2 - m. Sustituyendo en el miembro derecho de la desi­ E n c o n tra r cuáles de las s ig u ie n te s fuerzas son
gualdad anterior conservativas. E ncontrar la energía potencia l para
las que lo sean ( a y b son constantes y h es un vec­
(1 + a) M (M + m) (1 + a) M + m + am to r constante).
M [ (l + a) M + m\ (1 + a) M + m
1. Fx = a x + b y 2, Fy = a z + 2bxy, F¿= a y + bz2
1 am >1
(1 + a) M + m 2. Fx = ay, Fy = az, F = ax

3. Fr = 2a r sen 6 sen <p, F&= a r eos 9 sen <f>,

Condición que se satisface siempre, de donde se de­ F' « a r eos 0
duce que es más favorable usar un cohete de dos etapas.
Esto es razonable ya que con el cohete de dos etapas el 4. F=h x r
combustible m 2 se emplea únicamente en mover a M 2 y 5. F = hr
no todo el peso del cohete. El resultado que se obtiene
es debido en parte a que no se tiene en cuenta el au­ 6 - = F y = - k yY> F2 = ~ K Z
mento de masa total que supone usar el cohete de dos 1. F = F(r)r
etapas (debido a que requiere una maquinaria más so­
fisticada).
S olución

Una condición necesaria y suficiente para que una

fuerza derive de un potencial es que V x F = 0. Vamos a
 donde m es función de t y ve es la velocidad de escape
1 , v (m 0- / r ) l n í l - A \ + / í de los gases.
Eliminando explícitamente el tiempo, agrupando va­
riables iguales e integrando se obtiene
m
(v0 + Ve) (m 0 - m) + m ve ln
mA
Av =s v - vn = v ln —-
m

1. Cohete de una etapa:

La masa inicial del cohete es M = M x + M 2 + ml + m2
Un cohete consta de dos etapas de masas M 1 y /W2 y la masa final por conveniencia la escribiremos como
que contienen adem ás unas cantidades de com bus­ Afj + M v donde M 1 son los motores de la única etapa
tib le m , y m r A nalizar, sin te n e r en cuenta el efecto del cohete y M 2 el peso muerto. Así pues, el incremento
de la gravedad, cuál será la velocidad final del cohe­ de velocidad alcanzado después de haberse consumido
te y com pararla con la que se obtendría usando un el combustible es
cohete de una sola fase. ¿Qué relación deben satis­
facer las m asas de las fases y el co m b u stib le para
que sea favorab le el uso de dos etapas?
A/j + M2

2. Cohete de dos etapas:

El cohete parte con una masa inicial M = M l + M 2 +
m l + m 2 que al acabar la primera fase es m = M x + M 2 +
Mx+M2 m2, el incremento de velocidad obtenido en esta etapa
es

. , M. + M. + m, + nt>
Av, = v ln — ------ 2------1 L
e A/j + M 2 + m 2
Cohete de una etapa C ohete de dos etapas

Al acabarse el combustible de la primera fase se des­

prende el peso muerto M v La masa inicial al encender­
se los cohetes de la segunda fase es m = M 2 + m2 y al
S o lución acabarse el combustible la masa es M2. El incremento
de velocidad proporcionado por la segunda fase es
Al no actuar la gravedad la ecuación de movimiento
del cohete es sencilla
M2+ m2
Av.2 = ve ln
dv dm M,
-m —— = v ——
dt e dt

oo 23
m

y entonces el incremento total de la velocidad es

2
(Mj +M 2 + m l + m 2) (M2 + m 2) Partícula en dos y tres dim ensiones.
Av = Av.1 + Av,2 = ve ln
M 2(Ml + M2 + m 2) Potenciales centrales
Si comparamos el cohete de una etapa con el de dos

(M2 + m2) {Mx + M 2)

R = Av,2 e ta, p a s - Av,1 e ,ta p a = ve ln
M 2(M x + M 2 + m2)

El uso de dos etapas será energéticamente rentable si.;

R > 0, condición que es equivalente a exigir

(M2 + m2) (M 1 + M2)

>1
M 2(M x + M2 + m 2)

Para facilitar los cálculos hagamos M x = aM, M2 = M,

m 2 = m. Sustituyendo en el miembro derecho de la desi­ E n c o n tra r cuáles de las s ig u ie n te s fu erzas son
gualdad anterior conservativas. E ncontrar la energía potencia l para
las que lo sean ( a y b son constantes y h es un vec­
(1 + a) M (M + m) (1 + a) M + m + am to r constante).
M [ (l + a) M + m] (1 + a) M + m
1. Fx - a x + b y 2, Fy = az + 2 bxy, Fz ~ a y + bz 2
am
1 + >1 2.
(1 + a) M + m Fx = ay> Fy - az, Fz = ax
3. F = 2 a rse n 6 sen 4>, Fe = a rc o s 6 sen
Condición que se satisface siempre, de donde se de­ F' == a rc o s <f>
duce que es más favorable usar un cohete de dos etapas.
Esto es razonable ya que con el cohete de dos etapas ei 4. F=h x r
combustible m 2 se emplea únicamente en mover a M 2 y

u.
X
5.

ii
no todo el peso del cohete. El resultado que se obtiene
es debido en parte a que no se tiene en cuenta el au­ 6. Fx = - k ^ Fy = - k ry, Fz = - k zz
mento de masa total que supone usar el cohete de dos 7. F = F(r)r
etapas (debido a que requiere una maquinaria más so­
fisticada).
S o lución

Una condición necesaria y suficiente para que una

fuerza derive de un potencial es que V x F = 0. Vamos a

OA
usarla para estudiar las distintas fuerzas que se nos Reuniendo los resultados se llega a
plantean en el enunciado del problema.

V(x, y, z) = ~ - y ax2 - by2x - ayz - ~ bz3 + C

1.
V x F = 0,.F - 3 F , 3 F, - 3 3 Fy - 3 F,) =
= (a - a, 0 - 0 , 2 by - 2by) = 0 La constante C que aparece en la expresión es arbi­
traria y corresponde a la libertad que hay para fijar -el
Así pues la fuerza es conservativa. Para calcular explí­ orden de potenciales.
citamente el potencial usaremos la definición F = - W
2.
V x F = (0 -a , 0 -a , 0 -a ) ¿0
V{x,y, z) ^~¡Fxdx+g(y, z) = -— a x t- b fx +g(y, z)
Así pues esta fuerza no deriva de un potencial.
D onde g(y> z) es una función arbitraria de estas va­
riables que podemos determinar a partir de la ecuación 3. Para analizar este caso utilizamos la expresión del
rotacional en coordenadas esféricas

F = - -r—, esto es,

y dy e8 e,¡>
Vx F= 3r 3*d B,<j>
-,U 1L fe(y, z)
az + 2bxy = 2 by x -------^---- F rFn r sen 6 F,
óy

Aplicándolo a nuestra fuerza resulta V x F = 0. Como

y de aquí el rotacional es nulo la fuerza deriva de un potencial.
Para encontrar este potencial tenemos que partir de la
g(y, z) = - / a z d y = -a y z + h (z) expresión de la fuerza en coordenadas esféricas.

Llevando esta expresión a la ecuación de definición

vtí// , dV l dV l dV
de la fuerza a partir del potencial se obtiene que
= (r) = " "07 e ,~~7 ~í}8 e‘ ~ r s m g ~Xji

ay + bzL= a y - -rr— BV t
az F - —~—==>V=~ la r s e n Osen 4>dr=-arl sen 0sen 4>+g(9,4>)
r ór J
Finalmente, integrando la expresión anterior obtene­
mos Para determinar g(9, 4>) utilizaremos las ecuaciones
para las partes angulares de la fuerza

h{z) = - - ~ b z 3 + C
F,=----- \~^>ar eos 6 sen A=— a r eos 9 sen ó —— (9, ó)
6 r b9 r r B0 '
 dg($, ó) 6. Es fácil ver que
D e donde se deduce que —- r - ------= 0, y entonces
OO
1 dV
g(6, Por último, de la ecuación
Vx F = = 0
obtenemos
~KX - - t y ~k¿
1 1 dh
ar eos 0 = ----------- - ar2 sen 0 eos 0 ----------------- -1-7- Para calcular el potencial es posible proceder directa­
r sen 6 r sen 6 o<p mente integrando cada una de las fuerzas por separado
obteniéndose que
que nos permite encontrar que h(<f>) = C. La forma del
potencial es
V = ~ ~2 (k * 2 + kyy l +
V(r, 6, 4>) - -ar1 sen 0 sen 4>+ C
7.
Vx F=0
4.
V x (h x r) - h (V ■r) - (h ■V)r dV
F = - 7V(r) => F(r) - - — => V = V(r)
dr
Examinemos cada uno de los términos
siendo
(h • V)r - h
V ( r ) = ~ j r F{r)dr
H(V ■r) - 3h,

entonces
Una p a rtícu la en m o v im ie n to b id im e n sio n a l ex­
V x (h x r) = 2h * 0 perim enta una fuerza d irig id a según el eje Y y de­
pendiente só lo de la coordenada Y, esto es F{x, y) =
y la fuerza no deriva de un potencial. = F(y)ey , observándose que la tra ye cto ria en el pla­
no X Y es de la fo rm a x y = a2. E ncontrar la fo rm a de
la fuerza que o rig in a dich o m o v im ie n to .
5.
V x F = (Vr) x h + r V x h = j X h í O S olución

La fuerza no deriva de un potencial. Las ecuaciones de movimiento de la partícula son:

mx = 0
my = F(y)

29
 La primera puede integrarse inmediatamente ya que Las ecuaciones de movimiento son
corresponde a un movimiento rectilíneo y uniforme con
x = cte = u. Pero se sabe que y - a2lx entonces mv = mg - k \
S olución
a2 , a2u Descomponiendo en los ejes X e Y

m vx= - k v x \
2a2u 2a2u2 2u2y 2
m vy = -m g - kvy\

La primera ecuación puede integrarse inmediatamen­

Entonces la fuerza es
te y su solución es

v (0 = v (0)e

A este tipo de problemas donde se conoce la trayec­ Para resolver la segunda ecuación separamos las va­
toria y se nos pide inferir la ley de fuerzas se les deno­ riables e integramos
mina y aparecen por ejemplo al es­
tudiar trayectorias’de partículas subatómicas en campos dv k
- — dt
electromagnéticos. mglk + vy m

dv mglk + v
= - ln
-m * = - Jv mglk + vv m glk + v (0)

Un proyectil se dispara con velocidad v0 fo rm a n d o Usando que vy(0) = v0 sen a y despejando se llega a
un á n g u lo a con la h o rizo n ta l. E n co n tra r su p o s i­
ción en fu n c ió n de t y la ecuación de la ó rb ita supo­ vy(t) = -m glk + (v0 sen a + mglk)<rkl/m
n ie n d o que e xiste una fuerza de ro z a m ie n to p ro ­
p orcio na l a la ve lo cid a d (tiro con rozam iento). Para encontrar la posición en función de t integremos
las dos ecuaciones, teniendo en cuenta que Jt(0) = y(0) = 0.

m
x(t) = v0 eos a (1 - e~küm) —
K

m
y (í)= f v dt = -m g/kt + - — (v0 sen a + mglk)erk,lm
y K

= -m g/kt + (1 - (v0 sen a + mglk)

K
 La ecuación de la órbita puede obtenerse a partir de
las ecuaciones paramétricas para x(f) e y(t) sin más que
eliminar t. De la ecuación para x{t) se obtiene

1 - = e~

m kx
t ~~k mv„ eos a

y sustituyendo en la ecuación de y(t) se encuentra

(a) E ncontrar las coordena das del p u n to de im ­
/ \ m 2 8 In
i
kx pacto en el plano in clinad o si se dispara el proyec­
y(x)= — 1 -
til con á ngulo 6 (> a).
(b) ¿Cuál será ei á n gulo de d isparo para el que se
alcanzará la m a y o r distancia p o sible sobre el plano
+ (vn sen a + mglk) inclinado?
(c) ¿Para qué ángulo (s) de in clin a ció n recorrerá
la partícula una cierta distancia s a lo largo de su
trayectoria antes de alcanzar el plano inclinado?

Un p ro ye ctil se m ueve en el plano X Y actuando S olución

sobre él únicam en te la fuerza de la gravedad (tiro
parabólico).
1. S u p o n ie n d o esta p rim e ra situ a ció n s im p lifi­ 1. (a) Las ecuaciones de Newton para el proyectil
son
cada se pregunta
(a) E ncontrar la ecuación de la tra ye cto ria y{x } si
se dispara el p ro ye ctil con velocida d inicial 6 con el
d zx
m — ■■■■■■■= 0
eje X. d i2
(b) E n c o n tra r cu á le s son los v a lo re s p o s ib le s d 2y
para ios que el proyectil pasa p o r el p u n to de co o r­ m = -m g
~di*
denadas X Y prefijadas. ¿Cuántos valores diferentes
del á n gu lo llevan a soluciones?
con la condición inicial ±(0) = U eos 6, y(0) = U sen 0,
2. C onsiderem os ahora el p ro b le m a de d isparar
x(0) = 0, y(0) = 0. Las soluciones son
desde la base de una m ontaña a un p unto situ a d o
m ás arriba de ella. M odelizarem os el problem a por
a:(í) = U t eos $
el tiro desde la base de un plano in clin a d o un án­
g u lo a respecto a la horizontal a un p u nto situado
y{t) = U t sen 6 gv
m ás a rriba del m ism o.

33
 Combinando estas dos ecuaciones se obtiene la ecua­ la ecuación de la trayectoria pase por ese punto se ob­
ción de la trayectoria tiene, tras reagrupar términos

2 U2
y(*) = x t g 0 - - | ^ I sec2 0 x(6) = eos2 0 (tg 0 - tg oc)

(b) Si el proyectil pasa por un cierto punto (*, y) el y(0 )= x(0 ) tg a

ángulo de disparo se puede obtener de la ecuación ante­
rior. Usando la identidad trigonométrica sec2 0 = 1 + tg2 0 (b) La distancia recorrida sobre el plano inclinado es
la ecuación quecla
1 } - x { 0 f + y(6)2 = x(d)2(l + tg2 a) = sec2 a x 2(0)
gx¿
y - * tg 0- tg 2 0 y sustituyendo el valor de x sobre el plano se encuentra
2 U2 2 U2
2 U2
Que es una ecuación cuadrática para tg 0, que en ge­ L ( 0) = eos2 0 (tg 0 - tg a)
neral tiene dos raíces. Resolviendo la ecuación obtene­ §
mos
Esta expresión puede escribirse en una form a más
conveniente
u2
' ± = tg-1 1- í —n
g *\ V Í7 2g eos2 0 (tg 0 - tg a) = (sen 0 eos a - eos 0 sen a) =
eos a }
Cuando el discriminante de la ecuación anterior sea
cero cos ^ t /i \ sen(2 0 ~ a) - sen a
sen( 0 - a ) = *
eos a 2 cos a
U*_
y = entonces

tendremos un único ángulo posible. Los puntos en esta U2

L = ------- — [sen (2 0 - a) - sen a l
curva —denominada en el argot balístico parábola de g eos2 a J
seguridad— sólo pueden ser alcanzados con una posible
inclinación. E l m ayor valor posible de L corresponderá a que
Por último los puntos situados por encima de la re­ sen (2 0 - a) = 1, lo que ocurrirá cuando
gión definida por la parábola de seguridad no pueden
ser alcanzados de ninguna forma ya que el discriminan­ ir a
6 = — + -—
te es negativo. 4 2
2. (a) Para encontrar las coordenadas del impacto
en el plano inclinado nos damos cuenta de que los pun­ que es el ángulo bisector entre el plano y la vertical. En
tos del mismo son de la forma (x, x tg a). Exigiendo que este caso la distancia recorrida será
 (J 2 u2 a I sentr^tg 6) - senh"1(2tg a - tg 0) ~
L = ------- ;— (1 - sen a) =
g eos2 a g (1 + sen a) X V - ( 2 tg a - tg 0)[1 + ( 2 tg a - tg 0 ) 2]

(c) En cuanto a la distancia recorrida por la partícula Que nos da una ecuación que relaciona la distancia a
a lo largo de su trayectoria usamos la fórmula para la lo largo de la curva s con el ángulo de disparo del pro­
distancia recorrida sobre una curva en coordenadas car­ yectil 0.
tesianas, que es

Un cohete se m ueve en un plano vertical bajo la

acción de la gravedad y el im p u lso de sus m otores.
siendo x(0) la coordenada x del punto de caída del cuer­ Si !a ecuación de m o v im ie n to del cohete su p o n ie n ­
po sobre el plano inclinado. Sustituyendo para nuestro do que el m o to r le p ro p o rc io n a un im p u ls o en la
caso concreto dirección de m o v im ie n to es

t 2U 1 eos2 0(tg 9 - tg a)/g A / / ex \ _ dM . dv

^ ) = J0 Sy i + / t g 0 - - | jS ec2 e \ d x F+c —— a M
dt dt

p ro b a r que
siendo x( 0) la coordenada x del punto de caída del cuer­
po sobre el plano inclinado. Sustituyendo para nuestro
caso concreto d c dM
— - !og v = tan <A+ — sec ús
¿«A Mg dt Y
2
2 í / ! C0S: 0 (tg 0 - t g a)h . I PX \
Í ,8y i + t g 0 - - ^ s e c 2 0 )d x donde v e s la velocida d del cohete, M s u m asa, i^-el
á n g u lo fo rm a d o p o r la d ire c c ió n de m o v im ie n to
con la horizo n ta l y c la ve lo cid a d de los gases ex­
Para calcular la integral realizamos el cambio de va- pulsados p o r el cohete. Probar que si la e xpulsió n
£X de gases sigue la ley
riables senh z = tg 0 sec2 0, que lleva a
Ul
c dM
- ~ng
U f senh'1(tg 0) M dt
s(0) = — eos2 0 J cosh2 zdz
g s e n h “'( 2 tg tx - tg 9)
entonces

y de aquí B sec if/

v—
{sec i¡/+ tan ij/)n
U2 U2
s( 0) = — sen 0 + — eos2 0 x
2 g 2g
 Como tanto la velocidad de los gases como la del co­
hete están dirigidos según la tangente a la curva tendre­
mos que v = vt y c = -ct.

Por otro lado

Ay dv dt dv ¡d é t dtp \ dv dtp
rL=_- t+v ~ r= -rt+v H — sen cos ’K l=-r t4~ r n
dt dt dt dt \d t x dt yJ dt dt

Sustituyendo en la ecuación de movimiento encontra­

mos
X

dM . dv dtp
S olución - Mg(t sen tp + n cos tp) - c — ?— t - M —— t + Mv — — n
til til til

Vamos a considerar el movimiento en las direcciones Descomponiendo en las dos direcciones

tangente y normal a las trayectorias. Definamos los vec­
tores unitarios en dichas direcciones a partir de los ejes .. . dM dv
cartesianos - M e sen tp -c —— = M ——
b dt dt

t = ex cos tp + e sen tp 1 ( dtp

- Mg cos ip = M v ——
n= sen tp -e y eos tf.'tJ dt

y la relación inversa es Dividiendo la primera de estas ecuaciones por - Mg

cos ip obtenemos
ex = t eos t p - n sen tp I
= t sen tp + n eos tp j c dM 1dv
tfi tp + -z-------------- sec é = -------- — sec tp
Mg dt g dt
Escribamos ahora la ecuación del movimiento
Y de la segunda de ellas
dM dv
F+c— =M — 1 1 dtp
dt dt —= ------— sec é
v g dt
en componentes. Para ello vamos a ir analizando cada
uno de los términos que aparecen. d ln v \ dtp
■sec tp
En primer lugar la única fuerza que actúa sobre el co­ dv g dt
hete es la fuerza de gravedad
d ln v 1 dv
F = M g ==- M g ey = -M g (t sen tp + n eos tp) = --------— sec tp
dé e dt

38
100 Drobfemast fie M en á n ira
 Con este resultado tenemos identificado el segundo E n c o n tra r la s o lu c ió n de la e cu a ció n de m o v i­
miembro que se obtenía de la primera ecuación y en­ m iento de la partícula en los siguien tes casos
tonces llegamos a
1. Cam po eléctrico u n ifo rm e en la dirección del
d c dM eje X: E = (E, 0, 0)
— - ln v = tg ip + -— - sec <]/
di}/ Mg dt 2. C am po m a g n é tico u n ifo rm e en la dire cció n
del eje Z: B = (0, 0, B)
Para resolver el segundo apartado del problema usa­
mos el dato 3. C am po e lé ctrico en la d ire cció n X y cam po
m agnético en la dire cció n Z

dt S olución

que llevado a la ecuación anterior nos la deja en la for­ 1. En este caso la fuerza de Lorentz es muy sencilla
ma ya que se reduce a F = qE = qEex. Las ecuaciones de
Newton para el sistema son
d ln v
tg i¡/-n sec t/r
di}/ m x = qE
my = 0
Integrando en i}/ encontramos
mi - 0
ln v = J (tg i}/-n sec ifi) di}/-\a sec i}/-n ln (sec if/+ tg ip) + K
cuyas soluciones son inmediatas
Tomando exponenciales en ambos miembros
.. qE
B sec i/r + V°J + X°
(sec i¡/ + tg ip)n
y (0 = V + y0
siendo B = eK. z (0 = V +

2. Ahora la fuerza es de la forma

F = qv X B = q (yB, - xB, 0)

y entonces las ecuaciones de Newton son

La fuerza de Lorentz es la fuerza que actúa sobre
una p artícula cargada en el seno de un cam po elec­ x = c o cy
tro m a g n é tic o y tie ne la fo rm a
y = - <*>J
F = q(E + v x B) miz - 0

40
 Donde ojc = qBím es una cantidad auxiliar con dimen­ para la velocidad debe satisfacer la ecuaciones que te­
siones de inversa de tiempos llamada frecuencia de ci­ níamos antes de derivar, sustituyendo se encuentra
clotrón. Las ecuaciones para la velocidad son
-CjW. sen (úct+C2(¡>c eos coct=(oc(C3 eos (oct+C4 sen (oj) 1
-O j C
sen o)f+C.ü) eos c¡>t=- co(C. eos oit+C.
C * r C C c
Cv ¿ sen wí) J
1 C '
vy = - ojevx
mvz - 0 y entonces C4 = - C: y C2 = Cy Imponiendo ahora las
condiciones iniciales y^O) = 0 , v (0) = v0y se encuentra la
Derivando en la primera de estas ecuaciones y usan­ forma de la velocidad, que es
do la segunda se obtiene
v*(') = % sen
v + o)2v = 0
X e x
= % cos ^
que es la ecuación de movimiento de un oscilador armó­
nico cuyas soluciones son conocidas (sobre estos siste­
mas se discutirá ampliamente en un capítulo posterior) Y de aquí llegamos, tras integrar, a la forma paramé-
trica de la órbita
vx(t) = Cj eos o)ct + C2 sen a>ct
x{t) = x 0 + v0y(l - cos ü)ct)!o>c
siendo C¡ y C2 constantes que determinaremos después. y (t) = y Q+ Voy sen ú>ctlioe
Procediendo de la misma m anera se encuentra para
v.y z{t) = z0+ vt

vyfí)
\ /
= C.3 eos (oct + CA
4
sen <oc t que corresponde a un movimiento en forma de hélice.
La ecuación de la órbita en el plano X Y en forma explí­
Finalmente el movimiento en la dirección Z es libre, cita es inmediata de obtener (aislando en un miembro
entonces las funciones trigonométricas, elevando al cuadrado y
sumando las ecuaciones paramétricas).

2 = v®
^X - X Q % V + (/ y - y ) \2 y
d (ú I (O2
Observemos que — - (v2 + vz) = 0. Así pues las pro- ¡

ai y
yecciones de las órbitas en el plano X Y son circunferen­ que es la ecuación de una circunferencia de radio v0y/a)c
cias con vx2 + vy2 = v2Qx + v2
Qy
. y centro en (x0 + v0y/(oE, y0).
Pongamos el eje Y en la dirección de la velocidad ini­
cial en el plano X Y lo cual puede hacerse sin pérdida de 3. La fuerza es ahora de la forma
generalidad, entonces las condiciones iniciales son del
tipo v = (0, Vq , v0z). La expresión que hemos obtenido F = q(E + v x B) = q(E + yB , - xB , 0)
con lo que las ecuaciones de Newton son Integrando

x - qElm + (ory f1 o
*(í) = x0+ J v d f =X0+ —^sen cot~ (v0},+ EIB) (eos <at-1)
y = - ( ü cx
mi = 0 f' 1
? (*) = ?o + J0 vyM = ^ 0 + — (voy+ £ /S ) s e n a>ct +
c
que reescritas para la velocidad quedan

vx = qElm + (ovy + —— (eos (o t - 1) - EtIB

vy = -ü)e vx
rav = 0 Así pues el efecto de un campo eléctrico y magnético
Z
cruzados apuntando en las direcciones X y Z respectiva­
Procediendo igual que antes se encuentra mente es inducir, aparte de los giros, un movimiento
neto en la dirección Y. Esto queda patente en la ecua­
Í>
.r
+ íe)2 V
ex
= 0 ción de la órbita, que es
i>y + ú)2v
c y
= -qEo) Im

Las soluciones de estas ecuaciones son .r - +

vx(í) = Cj eos ú>cí + c 2 sen o>ct

qE
v (t) = C3 eos ü)J + C4 sen o>t -
moj

Como las soluciones deben satisfacer las ecuaciones Esta ecuación corresponde a un movimiento de tras­
de las que hemos partido para las velocidades sustitui­ lación del centro de la circunferencia con velocidad
mos y encontramos que C4 = - C l y C2 = Cr Entonces, v = (0, EIB, 0) = (E x B ). A este efecto se le llama
imponiendo que en t ~ 0 , y^O) = v0x y y,(0 ) = vQy
deriva de Hall y se observa en metales, en los que apa­
rece una corriente en dirección perpendicular a un par
qE de campos eléctrico y magnético aplicados.
vn0y = C,2 -
ma>_

se encuentran ios valores de C[ y C2 con lo que, sustitu­

yendo la solución queda
raet)C B Una partícula se m ueve en un potencial U{r) = K r2,
K > 0 con m o m e n to a n g u la r J = V 2 m d AK en una
Vj j ) = V0* C0S ^ + ( % + E IB ) Sel1 1 órbita c irc u la r, siendo K y d co nstante s positivas.
vy(r) = (v0>, + EIB) eos (út - v0x sen oa - EIB\ In sta n tá n e a m e n te se le c o m u n ic a una v e lo c id a d
 y entonces
v ~ 2 c f \ J ^ en dire cció n radial. C alcular la energía
E' = E + 8 T = 4K d2
en la nueva ó rbita.

S olución
En sus «Principia M athem atica Philosophiae Natu-
ralis», Isaac N ew ton propuso y resolvió el siguiente
La energía de la partícula puede escribirse como problem a: una partícula de masa m se m ueve en un
potencia) central, de m odo que el ángulo a entre la
velocidad y el vector de posición es constante. Em­
E ~ ~ mr 2 + - í - + U(r)
2 2mr2 v7 pleando argum ento s del m ism o tip o que los usados
en el problem a de los dos cuerpos para hallar el m o­
vim iento de la partícula relativa se pide
como / = V 2m d4fc , el potencial efectivo es
1. Probar que el p ro d u cto v r del m ó d u lo de la
U J r ) = Kri + ^ - velocidad p o r el m ó d u lo del ra d io ve cto r se m antie­
ne constante.
2. H allar la energía to ta l y la fo rm a de la fu n ­
La órbita circular debe corresponder al mínimo del po­ ción de potencia! U(r).
tencial efectivo para el valor de momento angular que
se nos da 3. H allar las funcione s r(t) y <t>(t) que represen­
tan el m o vim ie n to .

düf d4K 4. H allar cuánto tarda la p a rtícula en caer al o ri­

-£ = 2 K r - 2 — - = 0 gen, si los datos iniciales son: r0 = 1.000 m , v0 = 1 0
dr
ms-1 y a = 45°.

y resolviendo la ecuación se encuentra r - d. Así pues la

órbita está en r = d y la energía en ella es

1. Como el potencial_es central sabemos que debe

conservarse el momento angular "

Como el incremento de velocidad se realiza en la di­ J = mvr sen a - Cte.

rección radial deberá ser ortogonal al antiguo valor de
la velocidad, lo cual nos permite hacer y de aquí es evidente que vr es constante

J
vr ------------
2 2 m m sen a
 2. La energía total tam bién debe conservarse y Para encontrar 4>{t) utilizamos la otra ley de conser­
como vación del problem a, que es la conservación del mo­
mento angular

£ = i - m v 2+ U{r) [ — J— — \ \ = U ( r )
(¡>-
mr¿
podemos despejar U(r) que resulta ser
e integrando
J2 1
U(r) = E -
2m sen 2 a r2 dt =
23
m rl - ~ co ' a ( t - Q
Si tomamos el origen de energías en el infinito

£/(*) = E = 0,
=- y tg a log [mr¡ - 23 cot a ( t - ffl)]
y entonces la energía potencial queda
4. Cuando r = 0 se satisfará, despejando de la ecua-
P 1 23
U(r) = ción de movimiento, que r2 = — - cot a ( t - 10), de aquí
2m sen 2 a r 2
puede despejarse el tiempo de caída al origen que es
E l potencial efectivo se obtiene sumando el potencial
centrífugo al potencial físico r \ tg «
t - t Q= m = 70,7 s
23 2 v0 eos a

U« { r ) = U { r ) + 2 ^ = ' l ^ COtía

3. Para reducir el movimiento a cuadraturas usare­

mos la conservación de la energía
Una p a rtícu la se m u e ve en un p o te n cia l U(r) =
kr*, k> 0 .
y m r2 + XJef{r) =0
1. ¿Para qué valores de la energía y el m o m e n ­
to a n g u la r la ó rb ita será una circunferencia de ra­
2 3
r =± \! Uf= cot a dio a con centro en el origen?
m e} mr
2. ¿Cuál será el p e riodo de dicha órbita?
Integrar esta ecuación diferencial es inmediato 3. Si la p artícula se p erturba separándola lige­
ram ente de d ich o m o v im ie n to c irc u la r, ¿cuál será
el p e río d o de las oscilaciones radiales en to rn o a
r = a?
 1. Las órbitas circulares estables corresponden a los 3. D esarrollando el potencial efectivo en serie de
mínimos del potencial efectivo; calculémoslos potencias en torno a r - c i
S olución
U Ár) = + V(r) = + kr* dUef 1 d W f(a)
efK ' 2 m r2 v ’ 2mr2 V ef( a + 8 r ) = U ef( a ) + d r - A +■ (Sr)2 + ...
2 dr2

d V ef P
— —+ 4ka2 = 0 y haciendo uso de que el punto r = a es un extremo del
dr mar potencial efectivo se llega a la ecuación de movimiento

J = 2a3V m k
m
La energía se obtiene imponiendo r = 0 y se encuen­
tra que es un oscilador armónico de periodo

72
E = ~ — r + ka 4 = 3/m4
2ma¿

2. Llamemos A al área encerrada por la órbita, en-

dA Es periodo es —— veces más corto que el de un oscila-
tonces — es la velocidad aerolar. Sabemos que se cum- V6
dt
dor armónico, por lo tanto durante una revolución com­
pie
pleta se completan del orden de V ? oscilaciones en la
dirección radial.
T = —
dA
dt

Cuando el tiempo es t = T el área barrida es la conte­ Una p a rtíc u la se m u e ve en una ó rb ita c irc u la r
nida en toda la órbita A = mz2 y podemos aplicar la ley de ra d io a b ajo la acción de un p o te n cia l central
de las áreas que la atrae hacia un p u n to O. Sean vr, y v2 los va­
lores m á x im o y m ín im o de su v e lo c id a d . P ro b a r
que el p e rio d o de m o v im ie n to en dicha ó rb ita es
" ' 3t / Í íratv, + v2)
dt 2m Vm

obteniéndose el resultado
S olución

Si Vj * v2 el principio de conservación de la energía,

impone que el centro de la circunferencia debe ser dis-

51
100 problemas de Mecánica Partícula en dos y tres dimensiones. Potenciales centrales
tinto del punto O que es el centro de fuerzas. Llamemos
d a la distancia que separa O del centro de la órbita cir­ 3
cular. Colisiones
Como el potencial es central debe conservarse el mo­
mento angular, entonces

J - m (d + a)v2 = m (d - a)vl

El valor de d es desconocido, pero podemos usar las

ecuaciones anteriores para eliminarlo y expresar J en
función de los datos del problema: v2 y v2. Despejando
de la igualdad de la derecha se obtiene que

• v, - v,
d- a
1 ’ r2
P ro b le m a 22
y volviendo a sustituir se llega a

i; i>
J -■ 2ma S upong am os n m asas M colocadas en línea rec­
Vl + V2 ta de m anera que cada una toca a la siguien te. Des­
de la izquierda incide una m asa 2 M con velocidad
Para calcular el periodo podemos usar ahora la ley de v, que choca con la fila de n masas. Probar que es
las áreas para una órbita circular im p o s ib le que, co m o consecuencia de la co lisió n
que se su p o n e e lá stica , sea e x p u ls a d a una sola
A ira2 v +v masa de la derecha o dos m asas con velocidades
T ~ ------= 1- = 77V2 1 ¿
diferentes.
éá. iJ -\ viv2
dt
S olución
como queríamos dem ostrar. En el caso particular de
que Vj = v2 que corresponde al caso de que el centro de
fuerzas coincide con el de la circunferencia, T - Impongamos conservación de energía y m omento
vi para el primer supuesto, como toda la energía y el mo­
que es el resultado esperado. mento pasan a la última masa entonces

Conservación de la energía: 2M v2 = — Mu 2

Conservación de la cantidad de movimiento: 2M v - Mu

53
 1. Las órbitas circulares estables corresponden a los ¡ 3 . D esarrollando el potencial efectivo en serie de
mínimos del potencial efectivo; calculémoslos | potencias en tom o a r = a
S olución n 72 I dU
U J r) = — — + V(r) = -r — r- + kr4 U ef( a + S r ) = Uef( a ) + d r - ^ (* )* +
2m r2 2m r2 dr2
...

dU P y haciendo uso de que el punto r = a es un extremo del

+ 4ka3 = 0
dr potencial efectivo se llega a la ecuación de movimiento

J = 2aiV 'm k
Sr
m
La energía se obtiene imponiendo r = 0 y se encuen- |
tra que es un oscilador armónico de periodo

P TT
£ = + ka 4 = 3 ka4
2m a2

2. Llamemos A al área encerrada por la órbita, en-

dA i Es periodo es —j=r veces más corto que el de un oscila-
tonces —7- es la velocidad aerolar. Sabemos que se cum- ) V6
dt
dor armónico, por lo tanto durante una revolución com­
pie pleta se completan del orden de V 6 oscilaciones en la
dirección radial.
A
T=
dA
dt

Cuando el tiempo es t = T el área barrida es la conte­ Una p a rtíc u la se m u e ve en una ó rb ita c irc u la r
nida en toda la órbita A - ira2 y podemos aplicar la ley de radio a b ajo la acción de un p o te n cia l ce n tra l
de las áreas que la atrae hacia un p u nto O. Sean vy y v2 los va­
lores m á x im o y m ín im o de su v e lo c id a d . P robar
que el p e rio d o de m o v im ie n to en dicha ó rb ita es
7ra(v^ + v2)
dt 2m V tn

obteniéndose el resultado
S olución
17 * ! m
T=
Si Vj * v2 el principio de conservación de la energía
impone que el centro de la circunferencia debe ser dis-

51
100 problemas de Mecánica Partícula en dos y tres dimensiones. Potenciales centrales
tinto del punto O que es el centro de fuerzas. Llamemos
d a la distancia que separa O del centro de la órbita cir­ 3
cular. Colisiones
Como el potencial es central debe conservarse el mo­
mento angular, entonces

J = m (d + á)v2 = m (d - á)vx

El valor de d es desconocido, pero podemos usar las

ecuaciones anteriores para eliminarlo y expresar J en
función de los datos del problema: v, y v2. Despejando
de la igualdad de la derecha se obtiene que

dA =' a —1
V1----
" V2-

y volviendo a sustituir se llega a

V V S upong am os n m asas M colocadas en línea rec­

J=..2ma— l O .
ta de m anera que cada una toca a la siguien te. Des­
de la izquierda incide una masa 2 M con velocidad
Para calcular el periodo podemos usar ahora la ley de v, que choca con la fila de n masas. Probar que es
las áreas para una órbita circular im p o s ib le que, co m o consecuencia de la co lisió n
que se su p o n e e lá stica , sea e x p u ls a d a una sola
„ A tea2 v, + v, masa de la derecha o dos m asas con velocidades
T - = r- = Tea 1 2
dA diferentes.
dt
S olución
como queríamos dem ostrar. En el caso particular de
que vx - v2 que corresponde al caso de que el centro de
fuerzas coincide con el de la circunferencia, T ~ Impongamos conservación de energía y m omento
para el primer supuesto, como toda la energía y el mo­
que es el resultado esperado. mento pasan a la última masa entonces

Conservación de la energía: -i- 2Mv2 = Mu2

Conservación de la cantidad de movimiento: 2M v = Mu

53
 cidad v. C alcular el parám etro de im p a cto para que
La primera ecuación impone que u = V 2v mientras sea m á xim a la co m ponen te perpe n d icu la r a v de la
que de la segunda u = v, que son incompatibles. Hemos velocida d de B después del choque. H allar tam bién
probado que este supuesto es imposible. la dirección del m o v im ie n to de cada bola. S uponer
Analicemos el segundo caso, llamemos u x y u2 a las el c h o q u e c o m p le ta m e n te e lá s tic o y q u e no hay
velocidades con las que salen despedidas las dos par­ n ingún tip o de efectos.
tículas.
1 1 1
Conservación de la energía: — 2Mv2 = — M u\ + — Mu\ S olución

Conservación de la cantidad de movimiento: Llamemos u 'y v 'a las velocidades de las bolas salien­
tes y a y /3 a los ángulos que estas velocidades forman
2 Mv = M u, + Mu2
con la dirección incidente v. Impongamos la conserva­
ción de la energía y la cantidad de movimiento
El sistema a resolver es

Usando la segunda ecuación que relaciona v con u, y

u2 y sustituyendo en la primera se tiene

~ ( u ^ u ^ = ~ u \ + ^ -u \

y reordenando queda

(iq - u2)2 = 0

o sea que las velocidades de las dos partículas deben ser V = V + u

iguales y entonces no pueden salir dos partículas con ve­ |v|2 - KP + |u f
locidades distintas.
D e estas dos ecuaciones se deduce inmediatamente
que las velocidades de las bolas salientes deben ser per­
pendiculares. Para verlo elevamos al cuadrado la prime­
ra de ellas
Dos bolas de b illa r A y B de radio r e iguales m a­
sas se m ueven sin rozam iento sobre el ta p e te de v2 = v'2 + u'2 + 2uV
una mesa. B está en reposo y A se acerca con ve lo ­

55
Colisiones
 Para que se satisfaga la segunda ecuación entonces En el caso concreto del problema
debe ser u V = 0 lo que prueba la ortogonalidad de las
trayectorias salientes. •
A partir de la definición de a es evidente que
'- T
a = v c o T ^ e n to n c e s usando la conservación de la
energía encontramos una relación entre v 'y v
Problem a 24
v' = v sin a

en la que interviene el ángulo a. En 1932, Jam es C hadw ick d e scu b rió el n e u tró n
La componente perpendicular a la velocidad inciden­ p ro d u cie n d o un haz ai b o m b a rd e a r b e rilio con ra­
te de la bola B será diación a. D eterm inó su masa observando las c o li­
siones elásticas de los n e u tro n e s con p ro to n e s y
núcleos de nitró g e n o . O bservó q u e lo s p ro to n e s re­
\u\ = u' sin a'= v cos a sin a = — sin 2 a trocedían con una velocida d m áxim a 3,3 x 107 m s-1
"A... 2
y que los núcleos de Nu lo hacían con una ve lo ci­
7r . . dad m á xim a de 4,7 x 1Q6 ms~1. C alcular con estos
y entonces su valor sera máximo para a = ~ . A partir datos la masa del neutrón y su velocida d inicial.
V 4r2 - b2
de la figura puede verse que cos a - ------- ------ y en-
2r Solución
tonces

Vi V4 r2 - b2 Utilizaremos conservación de energía y cantidad de

movimiento
2 2r

y despejando aquí se encuentra el parámetfo de impac­ 1 2 1 7 1

-2 m/ = -2 « i « Í + Y m/
to que maximiza la velocidad trasversal u’L que es
m v = m bub + mnv'
b = V 2r
Donde el subíndice b se refiere a la partícula blanco,
La dirección de la otra bola se obtiene de la conserva­ sean protones o núcleos de N li. Reordenando
ción de la cantidad de movimiento en dirección trasver­
sal a v. mn (v2- v ' 2) = m bu2b
u' sin a = v' sin (3 mn (v - v') = m bub

v cos a sin a - v sin ¡3 sin a Dividiendo miembro a miembro eliminamos la masa

del neutrón
cos a - sin /3
ut = v + i/

O u in '
a= 2~P v' = ufe- v

56 57
100 problemas de Mecánica Colisiones
 y sustituyendo ahora en la ecuación de conservación de to n en 1923, y que co n stitu ye una prueba m u y im ­
la cantidad de movimiento p o rta n te de la n a tu ra le za c u á n tic a de la luz. Un
fo tó n de lo n g itu d de onda A tiene una energía ciné­
m n (V - Ub + V) = mbUb
tica h c lk y una cantidad de m o v im ie n to h f\, siendo
h la constante de Planck y c la velocidad de la luz.
m n O - «») = C o n sid e re m o s una c o lis ió n e n tre un e le c tró n de
una cierta lo n g itu d de onda A y un electrón en re­
Escribiendo ahora las ecuaciones para cada tipo de poso en la que el p rim e ro se desvía un á n gulo 9 y
partícula blanco y dividiendo miembro a miembro eli­ cam bia su lo n g itu d de onda a A'. R elacionar el cam ­
minamos de nuevo la masa del neutrón y llegamos a bio en la lo n g itu d de onda A - A' con el ángulo' de
desviación 9, su p o n ie n d o que A - A' es una ca n ti­
2 v (mN«N- m u ) = uNup (mN - m ) dad peqeña fre n te a A y A'. Relacionar la dirección
de salida del electrón con el á n gulo de desviación
1 m - m del fo tón.
V = — 1 U„ ll ---------- --------£-
2 N p m j p N - m pu p
(Ey.Py)
Para encontrar la masa basta sustituir en cualquiera (Ey.Py)
de las ecuaciones que involucran a m n

m ^ -r. m »U>___________ = m !JÁN - m pUp

m '1 ' 2 V ~ UP U 0n N Í m N - m n) ^ Up ~ U N
m Nu N - m U p P

Introduciendo valores numéricos se obtiene

mn = 1,96 x 1(P27 kg

v = 3,1 x 107 ms_1 S olución

El resultado para la masa del neutrón se aproxima al

Escribamos las leyes de conservación
valor real mn = 1,67 x 10~27 kg, este error es en parte de­
bido al hecho de que con estas velocidades tan altas de
las partículas vic — 0,1 es necesario considerar correc­ h h .
eos 9 + mv eos
ciones relativistas. A A

0= sin 6 + mv sin (f>

A
he = _he + _1
T mV2

Los fotones pueden c o lisio n a r elásticam ente con

electrones libres, efecto descubie rto p o r A. Com p-

59
m ae rio M p n á n i n a
 Reordenando las dos primeras tenemos ción para las componentes X e Y de la cantidad de mo­
vimiento se obtiene
h h
—----- -7 eos Q= mv eos <$>
A A h .
t t s in ®
h A sin 0„
7 sin 0 = mv sin 4> tan <p
A h h 1 - eos 6
T ~ Á7 c
Elevando al cuadrado y sumando eliminamos <p, que
no nos interesa que es lo último que se nos pedía en el enunciado del
problema.

m V s h Í ^ - +w ~ w cose
Problem a 26
Sustituyendo v2 de la ecuación de conservación de la
Dos esferas, cada una de d iám etro d, tienen m a­
energía llegamos a
sas M y m respectivam ente. Inicialm ente se encuen­
tran en reposo, estando sus centros separados una
2me i 1 1 lV 2 ,
. . . r kM m
h A A' - n ? ( c r a e - 1) distancia a y se atraen con una fuerza -------- — ,
A2
siendo x la distancia entre los centros en un instante
Eliminando denominadores a rb itra rio . E ncontrar la separació n m á xim a de los
centros tras el im pacto, siendo e el coeficien te de
h (A - A')2 restitución.
A -A ' = (eos 0 - 1 )
2me AA' me

S olución
Como el enunciado dice que A - A' es pequeño frente
v ^ podremos
a A, A ^ J • —
despreciar (A-AO
— —2 y entonces la ecua- Llamemos A y B las posiciones iniciales de los cen­
AA' tros, llamemos y B l a las posiciones de los centros en
ción anterior queda
un instante t arbitrario y A 2 y B 2 a sus posiciones en el
instante de la colisión. Llamemos G al centro de grave­
A' - A = (1 - eos 6) dad del sistema.
me

Que nos relaciona el ángulo de desviación del fotón

saliente con el cambio en la longitud de onda. Evidente­ ,— --------------9 0
, - - ■
------------- 1
A A, A2 B 2 Bi B
mente como A' - A > 0 la energía del fotón saliente es
m enor que la del fotón incidente que es lógico ya que
éste cede parte de su energía al electrón. Como la cantidad de movimiento total del sistema es
Por otro lado dividiendo las ecuaciones de conserva­ cero inicialmente y no hay fuerzas externas, esta canti-
 w r-

dad deberá conservarse con el tiempo, de donde se de­ ma

duce que G será un punto fijo. Definamos x = A XG. La Cuando x - la velocidad es nula v = 0 , enton-
M +m
conservación de la cantidad de movimiento es (y = B XG) ces podemos determinar C

M x - my = 0 km 2
C =~
a (M + m)
y entonces, integrando se obtiene y - — -. Como x(0) + 2km 2 m
m ----- v2 =
y(0) - a se encuentra que M +m M +m x a

ma Con este dato podemos calcular la velocidad de M

x justo antes del impacto que será
0 m +M
U<-
Ma 2k 1_ 1
yo u=m
m +M m +M d a

En el momento de la colisión xc + yr = d y entonces y por la conservación de la cantidad de movimiento se

encuentra la de m
md
x. = = 2k J l 1
m +M
Ul wm + M V d a
Md
m +M Llamemos vx y v2 a las velocidades después de la coli­
sión, se satisfará entonces
Llamemos v a la velocidad de M en A v entonces
trasformando la derivada segunda temporal mediante e K + »2) = v, + v2 = v, | l +
d dx dx d dv
x = —.— - = —— v - v — , la ecuación de Newton
dt dt dt dx dx
F = mx queda ya que después de la colisión M vx = mv2 (la cantidad de
movimiento en centro de masas sigue siendo nula); en­
tonces
dv kmM knü’M
-Mv
dx (x + y)1 (M + m )2x 2
me (u. + u0j * j 2k
v. = -i = me
M +m M +m
y de aquí, integrando
2k
km 3 v = ------N 1 L‘-= M p
C 2 M +m \ M +m
2 (M + m )2 x

62______________________ 63
100 oroblemas de M ecánica --------
 La conservación de la energía puede volver a aplicar­ de fo rm a que m , esté en reposo en LAB después
se tras el impacto; pero ahora la constante es diferente y de la colisión.
md
viene determinada por el hecho de que e n r = -------— la
m+
velocidad es v = v]} entonces S olución

km 1 1
C = Como las energías cinéticas son iguales, tendremos
M +m d a

y la ecuación de movimiento tras el impacto es —1 m tu]2 = —1 m2u\2 = —1 m2a2u\

1 km? 1 km 2 1 1
+
2 (M + m )2 x M +m d a de donde se obtiene que a 2 = m í/mr Si m x está en repo­
so después de la colisión la conservación de la energía
Cuando las masas vuelvan a quedar en reposo v = 0, requiere
entonces
y1 m xu 2\ + —1 m2u\2 = ~1 m 2v\
M +m a - é 1 (a-d)

m ad

m xu\2 = y1 m2v|2
Como el momento lineal total es nulo la masa m que­
da en reposo en el mismo instante de tiem po que la
masa M a una distancia de G dada por la expresión y la conservación de la cantidad de movimiento requie­
re que
M +m a - e2 (a - d)
Ü ad m xu, + m2u2 = (m¡ + am2) u x = m 2\ 2

Entonces la separación de los centros es m, + am.

u.
m-,
ad
A 1B 1 = x + y =
a - é2 (a - d)
Introduciendo este resultado en la ecuación de la
energía
Problem a 27

En una c o lisió n elástica de dos partículas de m a­

sas m 1 y m 2 las ve locidades iniciales son u 1 y a u 1
(a > 0). Si las e n e rg ía s c in é tic a s son ig u a le s en
LAB, e n con tra r las condiciones sobre u¿ u2 y m j m 2

64
y como a2 = m 1/m2 4
Form ulación lagrangiana
2 o? - ( a 2 + a2) => a = V 2 - 1 - 0,414 =3 a2 = 0,172
de la Mecánica
Entonces el resultado es

m J m 2 = 0 -1 7 2

u ju z ~ 0,414

Hallar p o r el m é to d o de Lagrange:

1. La ecuación del m o v im ie n to de un péndu lo

sim ple.
2. La de un p é n d u lo sim p le cuyo p u n to de sus­
pensión oscila a rm ó n ic a m e n te en una lín e a h o ri­
zontal con a m p litu d A y frecuencia ñ .

S olución

1. El sistema tiene.im .grado de lihpxtad. Se.a n fh 'n -

^gulo del péndulo^Tornando como origen de potencial la
posición de reposo, las energías cinética y potencial son
T = mPcP-fe y V = mgl{ 1 - eos a), siendo / y m la longi­
tud y la masa del péndulo. El lagrangiano es L = T - V .
^ * j * d dL dL
Cómo la ecuación de Lagrange es — - —----------• = 0, la
dt oá da
del movimiento tiene la forma bien conocida

á. + a)2 sen a = 0 ,

donde é 1 = gil.
 2. También este sistema tiene un solo grado de li­
bertad. Usando coordenadas cartesianas centradas en
la posición inicial del punto de suspensión (siendo x la
horizontal e y la vertical), la posición del punto de
suspensión es (A sen Clt, 0) y la de la lenteja del pén­
dulo, (A sen ü,t + l sen a, - l eos a). Las dos compo­
nentes de la velocidad de la lenteja, (v^, v ) = ( i , y),
son (A fl eos Clt + l eos aá, l sen aa). En consecuencia,
las energías cinética y potencial tienen la forma

T ~ ~ mv2-~¡^m (Pá2+ 2AlClá eos a eos ílt +A 2Í12eos2ílr),

V=
— K.COJT « . — i- cstj2-.

Sustituyendo en la ecuación de Lagrange, como en el

apartado 1, resulta la ecuación del movimiento buscada

A P I2 n 8
los trozos que rodean a las poleas. Sean x. las distancias
a = — -— sen iif eos a ~ ~ sen a en vertical de cada masa a la polea correspondiente, to­
mando como sentido positivo el descendente. Las coor­
denadas verticales de las tres masas son entonces (to­
El primer término del segundo miembro es la fuerza
mando como origen el centro de la polea superior),
de inercia inducida por la oscilación del punto de sus­
pensión.
z x= -T p Z2= L j x x + X 2, z 2 == L j —xj + — X2,

Problem a 29 Como vemos, el problema tiene dos grados de liber­

tad pues, dadas las posiciones de dos de las masas, que­
da determinada la de la tercera.
U tilizando el m é to d o de Lagrange, h a lla r la ace­
Las energías cinéticas son,
leración con la que se m ueven las tre s m asas de
una m áquina de A tw o o d com puesta com o la de la
1 1 1
fig u ra , su p o n ie n d o que las poleas y las cuerdas tie ­
ti = Y mA r 2 = Y m i(Á2- ^i)2’ Tz =■y m 3(-¿2 + ¿i)2-
nen m asas despreciables.

Las energías potenciales son V = -mgz.(el signo me­

S olución nos se debe a que estamos tomando el sentido positivo
de las coordenadas hacia abajo). Con las energías cinéti­
Tomamos el origen de potenciales en el centro de la cas y potenciales, escribimos el lagrangiano L - T - V ,
polea superior. Sean L x y L 2 las longitudes de las dos que resulta valer
cuerdas, desde la masa m x al centro de la polea inferior
y desde la masa m 2 a la m y respectivamente, sin contar
y como a2 = m l/m 2 4
4 Form ulación lagrangiana
2a2 = (a 2 + a2) => a = V 2 - 1 = 0,414 =$>a2 - 0,172
de la Mecánica
Entonces el resultado es

m Jm 2 = 0472

u ju 2 = 0,414

H allar p o r el m é to d o de Lagrange:

1. La ecuación del m o v im ie n to de un péndu lo

sim ple.
2. La de un pé n d u lo s im p le cuyo p unto de sus­
pensión oscila a rm ó n ic a m e n te en una línea h o ri­
zontal con a m p litu d A y frecuencia ft.

S olución

1. El sistema tieneam-giadíldeJihetiad-Sea^LsLán-
■
j>ulo del péndulo..„Tornando como origen de potencial la
posición de reposo, las energías cinética y potencial son
T = mPcd!2 y V = m g l(l - eos a), siendo / y m la longi­
tud y la masa del péndulo. El lagrangiano es L - T - V .
nComo la •' ade rLagrange es - 3d—d—-------
1 ecuación L B L— - 0, la
ut oá da
del movimiento tiene la forma bien conocida

a + cú2 sen a = 0,

donde co2 = gil.

 2. También este sistema tiene un solo grado de li­
bertad. Usando coordenadas cartesianas centradas en
la posición inicial del punto de suspensión (siendo x la
horizontal e y la vertical), la posición del punto de
suspensión es (A sen Oí, 0 ) y la de la lenteja del pén­
dulo, (A sen flí + / sen a, - l cos a). Las dos compo­
nentes de la velocidad de la lenteja, (vt, v^) = (x, y),
son (A íí cos ü t + l cos aóc, l sen a á ). E n consecuencia,
las energías cinética y potencial tienen la forma

T - - - m \P = ~ m (Pá2+ 2AlClá cos a cos Oí + A 2Cl2eos2Oí),

y =?

Sustituyendo en la ecuación de Lagrange, como en el

apartado 1, resulta la ecuación del movimiento buscada

AOP g los trozos que rodean a las poleas. Sean x. las distancias
a = — -— sen Oí cos a ~~j~ sen a en vertical de cada masa a la polea correspondiente, to­
mando como sentido positivo el descendente. Las coor­
denadas verticales de las tres masas son entonces (to­
El primer término del segundo miembro es la fuerza
mando como origen el centro de la polea superior),
de inercia inducida por la oscilación del punto de sus­
pensión.
z¡= xv z2 s Lj - + x2, z3 = L x ~ x x + L 2 ~ x 2.

Problem a 29 Como vemos, el problema tiene dos grados de liber­

tad pues, dadas las posiciones de dos de las masas, que­
da determinada la de la tercera.
U tilizando el m é to d o de Lagrange, ha lla r la ace­
Las energías cinéticas son,
leración con la que se m ueven las tre s masas de
una m áquina de A tw o o d com puesta co m o la de la
fig u ra , su p o n ie n d o que las poleas y las cuerdas tie ­ r r y fflA ^ y ^ W - r 3 = y m3{x2 + x,)2.
nen masas despreciables.

Las energías potenciales son V. = ~mgz.(e 1 signo me­

S olución nos se debe a que estamos tomando el sentido positivo
de las coordenadas hacia abajo). Con las energías cinéti­
Tomamos el origen de potenciales en el centro de la cas y potenciales, escribimos el lagrangiano L - T - V,
polea superior. Sean L x y L 2 las longitudes de las dos que resulta valer
cuerdas, desde la masa m ¡ al centro de la polea inferior
y desde la masa m 2 a la my respectivamente, sin contar ^ = T1 + T2 + T3- V 1- V 1- V 3

68 69
 1. H a lla r la e cu a ció n d e l m o v im ie n to de ta
= - j (m i + m 2 + m 2) x \ + (m 3 - m 2)x xx2 + y (m 2 + m ^x\
masa.
+ (mx - m 2 - m 3)gxl + (m2 - m 2)gx2 + constante 2. H allar la expresión del h a m ilto n ia n o y decir
si es constante y si coin cid e con la energía.
Usando las ecuaciones de Lagrange

d dL dL _ Q S olución
dt dx¡ dx. ’
1. E l sistema tiene un grado de libertad, para el
con ¿ = 1 , 2 , encontramos las ecuaciones del movimiento que tomamos como coordenada el ángulo a que defi­
ne la posición de la masa respecto al punto inferior
(m¡ + m2 + m 2)x1 + (m3- m 2)x2 - g(m l ~ m 2~ m3), del anillo, donde situamos el origen del potencial gra-
(m3 - m 2)xl + (m 3 + m 2)x2 = g(m 2 - m3). vitatorio. La velocidad de la masa es la suma de dos
vectores ortogonales: uno tangente al anillo, de m ódu­
Nótese que el término constante en el lagrangiano no lo aá, y otro perpendicular a su plano, de módulo coa
influye para nada en estas ecuaciones. Las aceleraciones sen a.-Por tanto las energías cinética y-potencial. valen
x. son las soluciones de un sistema de dos ecuaciones li­ T = m(a2á2 + oPa2 sen 2a)/2 , V = mga{ 1 - cos a). Como
neales. Resolviéndolas resulta el lagrangiano es igual a I = T - V y l a ecuación de
Lagrange es •
2m.(m 2 + m3)
(m 2 + m3) M - (m 3 - ra 2) 2 d dL dL q
dt dá da
2m 1(m2 - m 3)g
(m, + m 3)M - (m 3 - m 2)2 la ecuación del movimiento resulta ser

donde M = m 1 + m 2 + m 2 es la masa total. Las acelera­ aP g

ciones de las tres masas se calculan fácilmente a partir á sen 2a + — sen a - 0
2 a
deXj y x2. Valen
2. El sistema no es natural pues la energía cinética
Masa 1: z x - xv
no es una función homogénea de grado 2 de la velo­
Masa 2: z2 = - x x + x2, cidad á. Por el contrario es de la forma T - T2 + T0,
Masa 3: z3 = - x x - xr donde Tk es homogénea de grado k de la velocidad. Si
p a = dLldá es el momento conjugado a a, el hamiltonia­
no vale H = p j x ~ L = T 2- T 0 + V, o sea

Prohiem a 30
H = y m (a2ct- - aPa2 sen2 a) + mga( 1 - cos a)
Un a n illo de ra d io a que tie n e una cu e n ta de
m asa m e n sa rta d a está g ira n d o a lre d e d o r de un La energía en cambio es E = T + V = T2 + T0 + V, es
d iá m e tro vertical con velocida d a n g u la r co. decir

70 71
100 problemas de M ecánica Fnrmiilar.inn lanrannians» Ha ¡a M aránira
 vo hacia abajo, y C la longitud inicial de la cuerda, des­
E = j w i (Vó 2 + w2n2 sen 2 a) + mga (1 - eos a),
contando la parte que rodea a la polea. El lagrangíano
del sistema es
luego H & E. Además, el lagrangíano es independiente
del tiempo, luego H se conserva.
L = ~ m xx \ + i - m 2x \ + (m 1g x l + m 2gx2).
E l ham iltoniano es una constante del movimiento,
pero no coincide con la energía. Esta última no es cons­
tante, debido a que el sistema no está aislado, pues in­ Pero las coordenadas no son independientes pues sa­
tercambia energía con el motor que mantiene constante tisfacen la ligadura
su velocidad angular de giro.
x , = C - x x - j v (t)d t,

Problem a 31 por lo que, como suponemos que la función v(r) está

prefijada, el sistema tiene un solo grado de libertad. Es­
De una polea cuelga p o r un lado una masa m y y cribiendo el lagrangíano en función de x = x x mediante
p o r el o tro un m o n o de m asa m 2(< m T) que trepa la ligadura, queda
con velocidad v{t) respecto de la cuerda. E studiar
e! m o v im ie n to p o r el m é to d o de Lagrange. ¿Es po ­
s ib le qu e el m o n o pueda elevar la m asa m t tre p a n ­
do su ficie nte m e nte deprisa? S upóngase que en el
m o m e n to inicial el sistem a está en reposo. Aplicando las ecuaciones de Lagrange, se llega a que
j,
A. k <■'
la ecuación del movimientq tiene la forma

(m, - m 2) g - m 2v
m¡ + m2

Para que la masa suba, es necesario y suficiente que

x > 0 ya que el sistema parte del reposo. Por lo tanto el
mono sí puede elevar la masa m 2 si su aceleración res­
pecto de la cuerda verifica la condición

- m,
v > — I----- *-g.

S olución
E s tu d ia r p o r el m é to d o de L a g ra n g e el m o v i­
Sean x y x 2 las alturas de la masa y el mono, medidas m ie n to de una p a rtícula m aterial deslizándose sin
desde el centro de la polea y tomando el sentido positi- rozam iento p o r el in te rio r de una cavidad de fo rm a

72

100 problemas de Mecánica

 esférica. A la vista del lagrang iano ¿qué constantes El momento conjugado a 6 no se conserva y su expre­
del m o v im ie n to posee el sistem a? O b te n e r el ha- sión es
m ilto n ia n o y las ecuaciones de H am ilton.
dL
= !TX = mR1<>-
S olución
La expresión del lagrangiano se obtiene de su defini­
Sea m la masa de la partícula y i? el radio de la cavi­ ción H = p e8 + p <j>- L y resulta
dad. Si tomamos coordenadas esféricas con centro en el
de la cavidad, la posición de la partícula queda fijada
por los valores de los dos ángulos 8 y 4>t por lo que el •mgR eos 8 = E.
sistema tiene dos grados de libertad. Para construir el
lagrangiano, notemos que las coordenadas cartesianas
A partir de aquí, sustituyendo en las ecuaciones de
de la partícula valen x = R sen 8 eos <f>,y = R sen 6 sen
Hamilton p k = -dHidqk, qk = dHIdpk, se obtiene la forma
z = R eos 8, de donde explícita de las ecuaciones canónicas del movimiento
x = R8 eos 8 eos 4>- R<j>sen 8 sen <¡>,
p \ eos 8 . pn
y = R$ eos 8 sen <f>+ R<p sen 8 eos <f>, p = mgR sen 8 - * ~ ya
mR2 sen38 ’ mR2 ’
z - -R 8 sen 6.
p =o, ¿ = -------.
Como las energías cinética T = m (x2 -ty2 + i 2) /2 y po­ 0 m R2 sen 2 8
tencial V - - m gz, el lagrangiano es igual a
Finalmente, nótese que este sistema es dinámicamen­
L = T ~ V ~ ~ m R 2(Ú¿ + sen2 dfé) + mgR cog 8. te equivalente al péndulo esférico, pues tiene las mis­
mas ecuaciones v constantes del movimiento.

Se ve que la coordenada 4>es cíclica, es decir quemo

_apar.ece,..en el lagrangiano, por lo que su momento-con­ Problem a
jugado es una constante del movimiento

H allar las ecuaciones del m o v im ie n to de un sis­

p ,- = m R 2 sen28ó = constante.
* d<j> tem a fo rm a d o p o r dos masas iguales que se m ue­
ven en vertical, estando suspendida la prim era de
Esta cantidad-p.0 es eLmomento, angular alrededor del un p u nto fijo O m ed ia n te un m u e lle de constante k
^eje z- Como el lagrangiano no depende del tiempo, se con­ y la segunda de la p rim e ra m e d ia n te o tro m uelle
serva el hamiltoniano y como el sistema es natural (es de­ de la m ism a constante. E ncontrar las posiciones de
cir, la energía cinética es una función homogénea de grado e q u ilib rio del sistem a y escrib ir las ecuaciones de
dos de las velocidades), coincide con la energía. Luego los m o v im ie n to s en to rn o a esos puntos.

H = E = T + V - constante.

74
 wwwwwx Definiendo como nuevas variables y, = x ,- x k;, y2- x 2- x 2e,
quedan las ecuaciones

*1 2 m yx + m% - - k y v
m yx + my2 - - k y 2.

Problem a 34
k
C o m p a ra r las fo rm u la c io n e s n e w to n ia n a y la-
grangiana e studian do el caso de una partícula de
masa m que se desliza sin rozam iento p o r la parte
su perior de un plano inclinad o m ó v il de masa M. El
Solución plano puede deslizar sin rozam iento sobre un suelo
Sean /, y l2 las longitudes de los muelles. Tomemos las horizontal, ta l com o indica la fig u ra y está in ic ia l­
coordenadas como indica la figura: x x y x2 son las distan­ m ente en reposo.
cias entre O y la prim era masa y entre las dos masas, Hallar p o r los dos m étodos las aceleraciones de
respectivamente, con el sentido positivo hacia abajo. la p a rtíc u la y del p la n o , así c o m o la reacción de
Las energías cinética y potencial son pues éste sobre aquella.

T = y m [*2 + ( ^ + x2)2], S olución

1. Formulación newtoniana. Las fuerzas que actúan

1 1 sobre el sistema son la gravedad y la reacción del plano
V = - mgxí - m g(xx +x.2)+ — k (x x - lx)2 + — k (x 2 - 12)2.
sobre la partícula N. Sean (alx, aly), (a^, y ) las acelera­
ciones de la partícula y del plano. La segunda ley de
Una parte del potencial corresponde a la fuerza gravi- Newton aplicada a la primera se expresa
tatoria y la otra, a la elástica. El lagrangiano qs L - T - V
y las ecuaciones de Lagrange del movimiento resultan ser

2m xx + m x2 = 2 mg - k (x x - lx),
m xx + mx2 - mg - k (x 2 - 12).

■En las p osiciones de equilibrio se.anulan las^acele r a ^

piones.y, por tanto, los segundos miembros, por lo que
hay una sola cuyas coordenadas son

2mg y +.
x u ==~ + lv X2

76 77

100 problemas de Mecánica Formulación lagrangiana de la Mecánica

 N sen a = mau que expresa la condición de estar siempre la partícula
en el plano inclinado (C es una constante que no influye
N eos a - mg = m al . en el resultado final). Nótese que no consideramos la
r coordenada y 2 que es constante. Como consecuencia,
Como, el plano no se desplaza en vertical, sus ecuacio­ y = (x2 - Xj) tan », por lo que podemos expresar la coor­
nes del movimiento son
denada y, y su velocidad en función de las otras, toman­
- N sen o: = Ma.2x do el lagrangiano la expresión

0 = a.2y
L - ~ m [(1 + tan 2»)x 2 - 2tan 2ax 1x2] +
Pero hay una ligadura,_pues la masa m está obligada a
mantenerse sobre la superficie incIinadarSu aceleración
relativa al plano es + (M + m tan 2»)x 2 - m g tan a(x2 - Xj),

a rX = a u - a ^ a ry = ...
salvo una constante aditiva irrelevante.
Es muy fácil obtener las ecuaciones de Lagrange
y tiene que cumplirse que a j a =-tan~c£poi lo que
m (l + tan 2a)jf¡ - m tan 2»x 2 - mg tan » = 0
tan£*
-m tan 2axt + (M + m tan 2a)x 2 + mg tan a = 0.
lo que permite eliminar aiy, quedando las ecuaciones
Despejando las aceleraciones, se comprueba que el
mg ~ N eos a = (au - a^)m tan a, resultado es el mismo que obtuvimos en el apartado an­
terior con el método newtoniano, como cabía esperar.
N sen a = ma, ,
Nótese que mx: + M xt - 0, expresión de la constancia
~N sen a = Ma2x. del momento lineal.
Hemos eliminado una coordenada que «sobra» (es
La solución de este sistema, incluyendo la ligadura, es decir y,, aunque podríamos haber eliminado cualquiera
de las otras dos), lo que permite obtener de una manera
Mg sen a eos a (M + m )gsen2a muy sencilla las ecuaciones del movimiento. Pero para
.. - ’ . di —*”
u M + m sen2a ’ iy M + m sen2a calcular las fuerzas de ligadura hay que proceder de
otra m anera (ver A. F. R añada, Dinámica clásica,
-m g sen a eos a , M m s eos a Alianza Editorial, Madrid, 1990, Complemento 4B). Lo
cl, . - — -----------r---- , N = * que se debe hacer es tratar el problema variacional —o
^ M + m sen2» ’ M + m sen2»
sea el principio de Hamilton-— mediante el método de
2. Formulación lagrangiana. El sistema tiene dos los multiplicadores de Lagrange. Es decir, considerando
grados de libertad,, pues entre las tres coordenadas a las tres variables (xp y,, x2) como independientes y to­
xv y v x2, correspondientes las dos primeras a la masa m mando como langrangiano
y la tercera al plano inclinado, se cumple la ligadura
L '= L + A / ^ y " * (x ? + ^ ) + y i M ^ - m g y j + A ^ - Í X j - x ^ t a n a ],
f ( x v y v x2) = y í - (x2 - x j tan » - C = 0,

78
donde A es el multiplicador de Lagrange (prescindimos Una p a rtíc u la de carga q que se m u e ve en el
de una constante aditiva que se puede despreciar por­ seno de un cam po e le ctro m a g n é tico e xperim enta
que no influye en el resultado). Las correspondientes una fuerza, llam ada de Lorentz, igual a
ecuaciones de Lagrange son
Problema 35
r F = qr(E + v x B).
3/
mx. - A —— = m x. - A tan a = 0, D e m o stra r que esa fuerza puede d e d u cirse del
dxl
lagrangiano

my. + mg - A —3/
— - my. + mg - A = 0, 1
L = — m v 2 - qf> + q A ■ v,

3/ donde <f>{x, f) y A (x , t) son los potenciales escalar y

M x. - A ■■■■■ = Mx, + A tan a = 0.
dx2 vector y v, la velocida d. H allar ta m b ié n los m o m e n ­
tos canónicam ente conjuga dos a las coordenadas
Comparando con las obtenidas más arriba, se obtiene de posición y el ham ilto n ia n o .
el valor del multiplicador

mM g eos2 a S olución
A=
M + m sen 2 a
Los campos eléctrico E y magnético B se deducen de
Hay una fuerza de ligadura por cada coordenada y 3A
los potenciales mediante las relaciones E = -V ó - -r— ,
valen dt
B = V x A; en general son funciones del espacio y del
df mM g sen a eos a df tiempo. Sean x = (xp xv x3) las coordenadas cartesianas.
lx ^ dxl . M + m sen 2 a ^ ^ dx2 ,x’ Las ecuaciones de Lagrange son

p ^ df ^ mMg eos2 a d dL dL n , „ „
— -------- = 0, ¿ = 1,2,3.
iy dyy M + m sen2 a dt oxk óxk

Las fuerzas Fu , Fy son las componentes cartesianas El primer término vale

de la reacción del plano sobre la partícula; su módulo
es N, el valor obtenido ya con el método newtoniano. d dL d dA
— =;— = ~r (mx, + qA .) = mx. + ,q----- .
Nótese que Fu + F ^ = 0 es la expresión de la tercera ley dt dxk d t y k H kJ. k H ,dt
de Newton.
Según se hace más complicado un problema, tanto El segundo,
más ventajoso es tratarlo por el método lagrangiano.
dL d(f> dA
dxk ** dxk ' ^ dxk
 Teniendo en cuenta que que coincide con la energía, suma de las cinética y po­
tencial. Como se ve, el campo magnético n o interviene
dAk en la energía, lo que se debe a que no efectúa trabajo,
dt dt 3*. ’’ por ser siempre perpendicular a la velocidad.

pues v. = dx.Jdt y sustituyendo en las ecuaciones de La—

grange, se llega a la del movimiento siguiente Problem a 36

/ 3A . .A, 3A. \ 3<£ 3A El regulador de la figura está fo rm a d o p o r dos bo­

las de masa m ; cuatro v a rilla s a rtic u la d a s de
longitud I las unen a un punto fijo O y a u n anillo que
Usando la expresión de E en términos de los poten­ puede deslizarse a lo largo de un eje ve rtic a l. Las vari­
ciales dada más arriba, esa ecuación se pone en la forma llas pueden girar alrededor de ese eje y u n m uelle de
constante elástica k y longitud natural lQ u n e el punto
/ 3A 3A \ fijo y el anillo. El sistema gira en to rn o al e je con velo­
- s í f v' + v - cidad angular co. Se puede despreciar el rozam iento.

Los dos últimos términos pueden reescribirse en no­ 1. H a lla r las ecuaciones del m o v im ie n to y las
tación vectorial como cantidades conservadas.
2. C alcular la frecuencia de las o scila cio n e s pe­
~(v • V)A + V(v • A), queñas en to rn o al p u nto de e q u ilib rio .

ya que el vector v no depende de xk, sino sólo de t. De­ co

sarrollando, se comprueba que esa expresión es igual a
v x (V x A), con lo que se obtiene

m x = q(E + v x B)

que es la fuerza de Lorentz.

Los momentos conjugados a las xk son

dL
p t = - ^ = mvt +í¡A e

que no coinciden con m v k. E l ham ilto n ian o vale

H = 'Z{pkqk - L = vkdL/dvk - L. Sustituyendo resulta
S olución

H = — mv2 + q4>, El sistema tiene un solo grado de libertad, para el que

tomamos como coordenada el ángulo 6 indicado en la
figura.

82______________________
83
100 problemas de M ecánica Fnrmularión lanrannlana Ha la hAor'án'ma
 1. Tomando ejes centrados en O, las coordenadas 2. Los puntos de equilibrio cumplen 9 - 0 , 9 = 0, por
de las masas son lo que se sigue de la ecuación anterior que hay uno con

x l = l sen 9 eos <ut, - 1sen 9 sen o>t, z i = l eos 9,

9 = 6 = árceos { J f j mS¿. \ .
x2 = l sen 6 eos(<ot + rr), y2 = l sen 9 sen(otf + rr), z2 eos 6. 0 \ 2 k l + m atl I

Para estudiar las oscilaciones pequeñas a su alrede­

La velocidad de la bola 1 es dor, se linealiza la ecuación del movimiento en torno a
él, escribiendo 9 - 90 + e, con e «c 1, lo que indica un
x l = 19 eos 9 eos o>t - Ico sen 9 sen tot, movimiento que se mantiene muy próximo al equilibrio.
Los desarrollos del seno y el coseno son,en primer or­
y x = 19 eos 9 sen coi + lio sen $ sen coi,
den,
= -19 sen 0,
sen 9 ~ sen 90 + e eos 0O,
y una expresión parecida para la bola 2 . Sustituyendo
eos 9 ~ eos 90 - e sen 0.
en la energía cinética, ésta toma la forma
Sustituyendo en la ecuación del movimiento, resulta
T -2 ma}2l2 sen 2 9 + — m&l2j , una ecuación diferencial para e; para obtenerla, hay que
despreciar los términos cuadráticos en e y sustituir el va­
lor antes hallado de 0O; tras algunas operaciones, se lle­
siendo la potencial, suma de un término gravitatorio y ga a
otro elástico,
é+ e sen2 0O= 0 ,
1
V = 2mgl eos 9 + — k(2l eos 9 - l0)2.
2
que es la ecuación de un oscilador armónico de fre­
Por tanto, hay un sólo grado de libertad, expresado cuencia
por el ángulo 9. E l tiem po no aparece de .modo explí-
^cito^__e_Iii,agrangia.no,1,pfi.r.,lp_,^^lxLhamütomana„js,e
fí = sen 60 V ~ " + “ 2-
conserva, pero no coincide con la energía ya que la
energía cinética contiene un térm ino independiente
de la velocidad. Lo que ocurre es que eljsiste.maJU) Luego el sistema oscila armónicamente en torno a la
está aislado, pues intercam bia energía con el m otor posición de equilibrio, con frecuencia O que se puede
que lo mueve, la que hace que n o seco n serv e la ener- escribir como

La ecuación de Lagrange toma la forma

2ml29 = 2mglsen 9 + 2klsen 9(21 eos 0 - / o) + 2raco2/2sen20.

 U na p artícu la de m asa m se m ueve sin rozam ien­ y su velocidad por
to a lo largo de la curva sin u so id e de ecuación
X = x cos cot - (ox sen cot, Y = x sen cot + (ox cos cot,
Problem a 37
ttx ~ . ttx
2 = -a cos — , Z — ttx sen — ,
a

situada en un plano vertical Oxz, en el que Oz es la con lo que el cuadrado de su velocidad en el sistema
vertical ascendente. La partícula está som etida a la inercial será
g ravedad y el plano Oxz g ira en to rn o a Oz con ve­
locidad a ng u la r constante &». V 2= X 2 +Y 2 + Z 2 ^ co2x 2 + ^1 + tt2sen 2 i 2.

1. E n c o n tra r una c o n s ta n te del m o v im ie n to .

¿Coincide con la energía de la partícula? Como la energía potencial es V = mgz, el lagrangiano
vale
2. D e te rm in a r cuál debe ser ía fre c u e n c ia co
para que la p a rtícula pueda m antenerse en la p o si­ TTX TTX
ción x = a/2, z = 0, en reposo respecto a la curva. L =— m co2x 2 + 1 1 + 7T2 sen 2 — I x2 + mga cos
2
3. C alcular la zona en que se m ueve la partícula
Vemos que se trata de un-sistema no natural, porque
si (o = v 'lir g la , su p o n ie n d o que en el tie m p o inicial
su energía cinética no es una función hom ogénea de
está e n x = 0 ,z = -a , con velocidad v0 = V g a {2 - tt/2)
grado dos de la velocidad (contiene el término (Xx2 que
relativa a la curva.
no depende de ella"). Como consecuencia no rninriHpn
4. Expresar m ed ia n te una cuadratura el tie m p o el hamiltonianQJV.hu&n^gía^é.sJa. iilti m a n a ,s(»-mmpfva;
que tarda la p a rtícula en alcanzar la posición más lo que corresponde al hecho de que el sistema no está
alejada del O ze n el m o v im ie n to del apartado ante­ aislado e intercam bia energía con el exterior (con el
rior. motor que mantiene el giro).
no
5. H allar la reacción de la curva sobre la p a rtí­ hai& ilt£^noX QJategfi?L.d.^
cula en el instante in icia l de ese m o vim ie n to .
u -^L l L - rT rx \., 1 TTX
S olución 1 ~ ) T m ^í2C0S— '

Esta cantidad es pues la constante del movimiento

1. Sea O X Y Z un sistema de ejes inerciales cuyo eje pedida.
vertical O Z coincide con Oz. Si en el momento inicial 2. La ecuación de Lagrange del movimiento se ob­
coinciden los ejes Ox y OX, la posición de la partícula tiene fácilmente y es
viene dada por

X = x cos cot, Y = x sen cot, Z - -a cos — ,

a
 Queremos que x - a/2 sea un punto de equilibrio, que
— \ + moPx2 + 2mga eos (7rxfa)
deberá satisfacer x = x = 0 para todo t. Ello implica que
dt V m[ 1 + 7r2 sen 2(7rx/a)J ’
-oPai 2 + TTg sen rrl2 = 0 , por lo que la frecuencia debe
valer ^
expresión que permite expresar el tiempo T buscado en
la forma integral
oj ~ . ! 27r8
1 + tí2 sen2( ttxla)
2hlm + oPx2 + 2ga cos( irxla)
3. Como en x = 0 hay un mínimo, la velocidad es
tangente aleje Ox y x 0 = Vag(2 - ttI2). Sustituyendo
Sustituyendo los valores de h y co y con el cambio
en la integral de Jacobi, resulta que h{0) = -mgairlA y, x = au, resulta
como ésta es constante, se cumple siempre

T - \ l la f I/2W l + ir2 sen 2( 7ru) ^

^ \1 x2
1 + 7r2 sen2o — O - oPx2 -1 eo s ^ = — -—
? 7 - 2ga
'o Y 7t ( 4 m 2 ~ 1 ) + 4COS 7TU
a f a 2

Como el primer término es positivo o nulo, el movi­ 5. En. í = 0, x = 0, la velocidad es x0 - V ag (2 - 7r~/2 );

miento tendrá lugar en la zona en que se cumple sustituyendo en la ecuación del movimiento, resulta que
x0 = 0. Como Z = -a cos(-rrxla), derivando dos veces, se
obtiene
, , ~ irx ^ vag
oPx2 + 2ga eos — > ——
& 2*
„ 7T2 / 7TX \ , . /tT x \ ..
Z = — eos i — 1x 2 + 7r sen | — 1x,
El signo igual en esa relación corresponde al punto de a \a ) \a j
retroceso, en el que se anula la velocidad. Sustituyendo
el valor de a», se tiene con lo que 2 0 = 7r2g (2 - tt/2). Derivando las expresiones
de X , Y en función de x y evaluando en t - 0, resulta
7TX ( 1 X2
eo s > 7T r X Q- 0 , y 0 - 2 x 0ü > = 2g V t 7( 4_ - 7 7 ).
a 1 4 a1
La segunda ley de Newton dice que la reacción R de
La desigualdad se satisface en -ai2 < x <aí2 que es ia la curva sobre la partícula vale
zona accesible'al movimiento (o sea X 2 + Y 2< a2!4). Nó­
tese que x = -ai2 y x = a/2 son dos puntos de retroceso.
R = m (X i + Y¡ + Z k ) + m gk,
4. La posición más alejada del eje z es x - a!2.
Como la integral de Jacobi (o hamiltoniano) es una can­
(pues los. ejes X , Y, Z son inerciales). Sustituyendo los
tidad conservada (h ~ constante), se puede usar para re­
valores de las aceleraciones que acabamos de calcular,
ducir el problem a a cuadraturas de la misma m anera resulta
que se suele hacer con la energía, si ésta se conserva. En
la expresión obtenida en el apartado 1, se puede despe­
/? y = 2 m g V 7 t ( 4 - 7 t), R z ~ m g [ \ + 7t2 ( 2 - 7 t / 2 ) ]
jar la velocidad en la forma
 5
M ovim iento planetario
f

Problem a 38

S uponiendo que la ó rbita de la Luna sea circu la r,

calcular

1. La velocida d de su centro a lre d e d o r del de la

Tierra
2. Las velocidades de los centros de ios dos as­
tro s a lre d e d o r de su ce n tro de masas.

Datos. .Radio de la órbita: f í = 3,844 x 108 m ; M a ­

sas de la T ie rra y la Luna: M r = 5,977 x 1024 kg;
M l = 7,380 x I O 22 kg.

S olución

1. La tercera ley de Kepler dice que

W _ G(M t + M l )
~T* ~~ 4^ ’

siendo T el periodo de la órbita. Despejando, resulta

T = 2,357 x 106 s = 27 d 6 h 40 m. (En realidad, el perio-

91
 do de la Luna es aproximadamente una hora mayor, de­ La tercera ley de Kepler dice:
bido sobre todo a la atracción solar, no incluida aquí,
que tiene un efecto retardador del 0,14%.) La velocidad
fl3 rz c m, + m2 S olución
de la Luna respecto al centro de la Tierra es pires T2 47T2

v = í lü ? = x.025 m/s. siendo a el semieje mayor y T el periodo. Despreciando

T las masas de los satélites respecto a las de los planetas y
dividiendo esa expresión aplicada a los dos casos, se ob­
El centro de masas del sistema Tierra-Luna está en la tiene para el cociente de las masas de la Tierra y Júpiter
>rZ
recta que une los dos centros, estando sus distancias a
ellos en la proporción inversa de las masas. Si a y b son
esas distancias, se tiene

M, AL
a- R, b= R, La masa de una esfera vale M = 47ri?3cr/3 , siendo a la
',vj> M j + Af¿ Mt + Ml densidad; sustituyendo en la expresión anterior, se tiene

o sea a = 4.688 km, b - 379.700 km. Nótese que el cen­

tro de masas está en el interior de la Tierra a 1.670 km
bajo la superficie.
2. Los centros de la Luna y de la Tierra siguen órbi­
tas circulares de radios a y b, por lo que sus velocidades Teniendo en cuenta que T J T G = 39.340/10.300, a(JR / =
S olución = 15, aLlR T = 60, se llega a
serán la de la Luna respecto de la Tierra calculada más
arriba, multiplicada por a/T y b/R, o sea
cr.
= 0,23.
v_ = 12,5 m/s, v = 1.012,5 m/s.

O sea, que Júpiter es entre cuatro y cinco veces me­

nos denso que la Tierra.
P r o b lé tiii

G aním edes, el te rce r satélite de J ú p ite r, y la Luna Problem a 40

siguen ó rb ita s a p ro xim a d a m e n te circu la re s de ra­
dios aG = 15 f í, y = 60 Rr, siendo Rj y Rr los ra­
dio s de J ú p ite r y la Tierra. Sus perio d o s son, res­ S a b ie n d o que la e x c e n tric id a d de la T ie rra es
p ectivam en te, 7 d 3 h 4 0 m y 2 7 d 7 h 4 0 m. Hallar e = 0,01673 y que su sem ieje m a y o r es a = 1,495 x
con estos datos la relación entre las densidades de x 108 km , e n co n tra r las distancias m á xim a y m ín i­
J ú p ite r y la Tierra. ma al Sol y el sem ieje m e n o r b. H allar ta m b ié n las
velocidades de la Tierra en el perih e lio , el afelio y
los extre m o s del ¡atus rectum.

92 93
 De las propiedades de las cónicas, se deduce que las déos p o r tratarse de una ó rbita m u y excéntrica se
distancias al perihelio y al afelio valen r_ = a{ 1 - e), recom ienda to m a r G M = 1,3277 x 1020 uS l usando
r+ - a(l + e); el semieje menor, b = a V l - e2; y el se- al m enos cuatro cifras sig n ifica tiva s en cada paso.)
Solución
milatus rectum, p = a (l - e2).
De la ecuación de la órbita y la relación entre el se­
mieje mayor y la energía E = - GM ml2a, se sigue S o lución

En primer lugar, determinemos el tipo de trayectoria,

v2 = GM para lo que hay que calcular la energía
( 4 4 ) .

(Nótese que se desprecia la masa de la Tierra respecto a

la del Sol.) £ = — - i - v 2- — --3,2295 x l ° 6 J.
Resulta entonces: m 2 r

Distancia al perihelio: r = 1,47 x 108 km. Como E < 0, se trata de una elipse. Para determinar
su excentricidad, hay que calcular su momento angular
Distancia al afelio: r = 1,52 x 108 km.
Semieje menor: b = 1,4945 x 108 km. j7= ~ = vr sen — = 50,9 x 103x 10n —~ = 3,5992 x 10!5m 2/s,
J m 4 VI
Velocidad en el perihelio v en el afelio va y en los
extremos del latus rectum v,Ir :
de donde
GM 1 + e,\M
vp ~ j ) - 30,3 km/s,
a 1 -e
e=V i+ ^ = ° ’976-
GM 1 - e
vn= | -------- r — ~ ) = 29,3 km/s, Se trata de una elipse muy elongada, pues la excentri­
a 1+e cidad es próxima a uno. Su semieje mayor se obtiene de
i a = GAf/2¡£¡, de lo que resulta a = 2,056 x 1012 m, o sea
'GM 1 + e2 a = 2.056 millones de kilómetros. El semilatus rectum es
= 29,8 km/s.
a 1-e 2 p = a( 1 - e2) = 3 zIGM = 9,757 x 1010 m.
La ecuación de la órbita es

~ = 1 + e cos(í> - <pQ). -
Un com eta se m ueve de m o d o que en un cierto
instante está a 108 km del So! y se acerca a é! con Sólo queda calcular el ángulo <pQque da la dirección
la velocid a d de 50,9 km /s fo rm a n d o un á n g u lo de del perihelio. Se cumple
45° con el rad io v e cto r desde el Sol. H allar la ecua­
ción de la ó rbita del com eta.
cos{(p-cp0)
(Al ser este pro b le m a m u y sensible a los redon- er

94 95
100 problemas de Mecánica Mrn/imiontn
 Si tomamos como origen de azimutes la línea que va Al comunicarle el impulso transversal mv, el satélite
desde el Sol a la posición observada del cometa, resulta adquiere un momento angular por unidad de masa y
que para <p= 0, r = Í0 11 m, de donde se obtiene cp0 = 91,42°.
Las distancias mínima y máxima al Sol son ; Rv
y= vh =
1 + e’
Al perihelio: r_ = a( 1 - e) = 4,934 x 10i0 m; al afelio:
r+ = a{1 + é) = 4,063 x 1012 m. Vemos que eLcometa tie­ y su energía cambia de E a E'
ne una órbita muy excéntrica que entra d eítro de la de
la Tierra y sale de lá de Neptuno. f 2 '
Er v = E7 - + —
1
- m v 22 = — m (u¿ + v¿)2 \ G M m
Así queda completamente determinada la órbita. 2 2 R

El semilatus rectum p,
Problem a 42 f R 2v2
P = GM GM{ 1 + e)2’
Un satélite artificial se lanza v e rtica lm e n te hacia
a rrib a con ve lo c id a d in icia l u. En el m o m e n to en La excentricidad y el semieje mayor se deducen fácil­
que llega a pararse se le da una velocidad tra n sve r­ mente.
sa! v. H allar la ó rb ita en fu n ció n de u y v, despre­
ciando la resistencia del aire.
Problem a 43

Solución Dos satélites de la m ism a masa m se m ueven en

dos ó rbita s coplanarias, en el m ism o sentido, a lre ­
d e d o r de la Tierra. Una de ellas es circu la r de radio
La altura h que alcanza en el movimiento vertical /? y la otra elíptica con distancias al apogeo y peri-
(medida desde el centro de la Tierra) se deduce fácil­ geo de 8 f í y fí, respectivam ente. Los dos coinciden
mente de la conservación de la energía en el p u n to de contacto de sus ó rb ita s A y se aco­
plan entre sí para c o n tin u a r su m o v im ie n to unidos.
1 . GMm GMm H allar las características de la nueva órbita.
•mu-
2 R h

Definiendo e como el cociente entre, las energías ini­

ciales cinética y potencial e = - R u 2/2GM, resulta

1+€

(Nótese que e = -1 corresponde a la velocidad de esca­

pe.)

97
96

100 problemas de Mecánica Movimiento planetario

 Los subíndices 1, 2, 3 se refieren a las órbitas circu­ Nótese, finalmente, que en la colisión se pierde la
lar, elíptica y tras la colisión, respectivamente. Los se­ energía
S olución miejes mayores de las dos primeras órbitas son: a1 = R,
a2 = 9RI2. Como
LA EZ7 = rr + E7 -2 - £3
r-
= 1 m(v, - v2)2 = GMm
36R

v = V GM(“ 7 }
Problem a 44
se sigue que, en el instante de la colisión Un satélite de masa m describe una ó rbita circ u ­
lar a una altura de 1.600 km sobre la superficie te ­
GM GM rre stre . U na pieza de d e se ch o espacial de m asa
3 R m /200 se m u e ve so b re la m is m a ó rb ita p e ro en
sentid o co n tra rio , de m o d o que choca con el satéli­
te y queda unida a él. C om o consecuencia se m o d i­
La colisión es inelástica, pues no es posible conservar
fica la ó rb ita . ¿Cuál es la distancia más corta entre
la energía y el momento simultáneamente. Parte de la
la superficie terrestre y el o b je to com puesto p o r el
energía se invierte en deformar los satélites o se trans­ satélite y la pieza?
forma en calor, pero el m om ento se conserva, por lo
que la velocidad del cuerpo es v3 = (Vj + v2)/2 y su ener­
gía

1 2 GMm 23 GMm
3 “ ~2~ 3------- * “ 36 R

(Que corresponde a un cociente de energía por unidad

de masa igual a -23GMI12R.) Por ello el semieje mayor
de la nueva órbita vale

GM2m 36 „ S olución
Ü3~ ~ 2E3 “ 23 '

Las distancias al perigeo y al apogeo valen R y Los dos cuerpos tienen la misma velocidad, pues la de
una órbita no depende de la masa del móvil ni del senti­
do del movimiento. Como se trata de órbitas circulares,
49
r3+^ 2 a 3 - R ^ a3(l + e3) = — R, m v2/(R + h) - G M m/(R + h)2, de donde la velocidad
de los dos cuerpos es v = V G M /(R + h), siendo R el ra­
dio de la Tierra y h = 1.600 km la altura sobre la superfi­
de donde la excentricidad de la órbita final es e3 = 13/36 cie terrestre; notaremos también m, m ' las dos masas y
= 0,36. m'im = e = 1/200. Para determinar la velocidad del obje-

99
100 problemas de Mecánica
 to compuesto tras la colisión V, se emplea la conserva­ Una estación espacial que está en ó rbita circular
ción del momento lineal de radio Rn a lre d e d o r de la Tierra debe ser tra sla d a ­
da a otra ta m b ié n circu la r y de radio R2 > Rr Se de­ Problem a 45
¡71 — ¡71 cide hacerlo sig u ie n d o una órbita de transferencia
(m - m')v = (m + m')V, =» V = -------------v. consistente en una sem ielipse con distancias al pe-
m + m'
rig e o y al a p o g e o fí, y R2 re s p e c tiv a m e n te . Para
Por tanto la energía del objeto compuesto vale ello, se com u n ica en un cierto m o m e n to a la esta­
*
ción un in cre m e n to de velocidad A i' para que pase
7- 1 GM (m + m') GMm l + 6 e + e 2 a la ó rbita de transferencia ; al lle g a r a su apogeo»
E = — (m + m ^V 2------- ---------- - . se com unica o tro increm ento Av2 para colocarlo en
2 ' R +h R +h 2(1 + e)
la circunferencia de radio R2.
Está claro que el punto de choque ^erá uno de los áp­
sides de la nueva órbita pues la velocidad es allí perpen­
1. H a lla r las c a ra c te rís tic a s de la ó rb ita de
dicular al radio vector. Sea t i la altura del otro. Sabe­
transferencia;
mos que la energía y el semieje mayor de la órbita valen
2. H allar Av^ y A v2 en fu n ció n de R, y R2;
_ GM (m + rrí) r +r 2R + h + t i 3. H allar la du ra ció n de la transferencia.
E = --------- , a = --------- ,
la 2 2
(Este tip o de transferencia , conocida com o órbita
donde r_, r+son las distancias al centro de los puntos ap- de H ohm ann, fue propuesta p o r este a s tró n o m o en
sidales (peiigeo y apogeo). Sustituyendo las expresiones 1925.)
anteriores, resulta

y tomando valores numéricos (con R = 6.400 km)

t i = 1.280 km.

Por tanto el punto de choque es el apogeo de la nue­

va órbita. Se entiende que sea así, porque dTchoqueTia
hecho disminuir la velocidad desde v a V = v(l - e)/(l + e)
disminuyendo la fuerza centrífuga, por lo que el objeto
c^eJiaciasuperigeo. ■
/ Como r± = a (l ± e), r+= 8.000 km y f +- 5.120 kiü, re-
/ sulta que a - 6.560 km, e - (r+- r )/2a~=~0',22. También
podría resolverse el problem a calculando prim ero 'el
momento angular J, después la excentricidad en función'
de E y / , y por último r . ■ j
i\

F/¡-
100 101
 1. En la órbita de transferencia, las distancias al pe- S u poniendo que la Tierra sigue una ó rb ita c irc u ­
rigeo y apogeo son evidentemente r_ - <2(1 - e) = R y lar, co nsidere m o s los com etas o m e te o ro s que se
S olución r+ = a( 1 + a) = R v de donde el semieje mayor a y la ex­ m ueven en ó rb ita s parabólicas. P roblem a 46
centricidad e valen
1 . ¿Entre qué lím ite s de ve lo c id a d re la tiva e n ­
Rl+ R2 trarán én la a tm ó sfe ra , d e sprecia ndo la a tra cció n
a= de la Tierra?
^2 + R l
2. ¿Cómo varían esos lím ites si se tiene en cuen­
La energía y m om ento angular de esta órbita son ta la excentricidad de la órbita de nuestro planeta?
iguales a

GMm GMm
E= / = | GA^n 2
2a R i + ^2 R + R2

2. La velocidad en la primera órbita es v = A/GM/R

y en la de transferencia v = V G M (2/r- 1/a), con lo que,
sustituyendo el valor de a y haciendo r = r_, resulta para
el incremento de velocidad

GM 2R,
A v , = -1
R , R í + R2

A nálogam ente, la velocidad en la órbita final es

v = V GMIR2; restándole la de transferencia en el apo­
geo (con r = r+), resulta el segundo incremento de velo­ S olución
cidad
1. La velocidad de la Tierra en su órbita, supo­
GM 2i?, niéndola circular de radio a T = 1,495 x 10 11 m, es
Av2 = 1-
vT = VG M /üj. = 29,9 km/s. Como el cometa o meteorito
viaja en órbita parabólica, su velocidad es vc = \ /2GMIr.
3. La duración de la transferencia rR será la mitad A l c ru z a r la ó r b ita de la T ie rr a lle v a s ie m p re
del periodo de la órbita r, que se calcula según la terce vc = V 2 GMIar = 42,2 km/s, con independencia de la di­
ra ley de Kepler rección que lleve. Por eso las velocidades relativas má­
xima y mínima entre el cometa y la atm ósfera son la
suma y la diferencia de estas dos.
Tr, = = 7r
GM 8GM

vmáx = Y - ^ 7 ( V 2 + ^ = 72,1 km/s>

102
 67tG M 2A7i3r z
— ( V 2 - l ) = 12,3 km/s.

(b) Si se considera la excentricidad de la Tierra, es fá­ más que en la teoría new toniana.
cil comprender que la velocidad máxima se encuentra 3. A p lic a r estos resultados a M ercurio y v e rifi­
cuando el choque de frente tiene lugar en el perihelio y car que el avance acum ulad o de la órbita después
la mínima cuando el choque es por detrás en el afelio. de un siglo es de cerca de 43" de arco.
Sus valores son

GM ( Á I 2 „ / 1 +e S olución
1- e V 1- e 1. Igualando las fuerzas centrífuga y centrípeta so­
bre el planeta, se tiene
GM L 2 J l-e
1+e fv 1 + e „ m v2 GM m ( 6v2
F = ------= ----- =— 1 + ——
y despejando v
Problem a 47
La teoría de la gravitación de Einsteín im plica una pe­ 6GM ~m ( G M \ m
1-
queña corrección a la fuerza newtoniana. Según ella, la rc¿
fuerza gravitatoria producida por un gran cuerpo de
masa M sobre un planeta de masa m que se mueve con Como v T = 27Tr y la tercera ley de Kepler dice que
velocidad ven una órbita circular de radio res Tn = 2777 \'r/G M , resulta tras sustituir
_ GMm í 6 v 2\
F= : 1 + --- , 127rV2 \
T= T, = tA 1-
Tic2
I'
s ien do c la velo cid a d de la luz.
Nótese que 6 GMIc2r es una cantidad pequeña, por lo
1. D em ostrar que si el p e rio d o n e w to n ia n o de que se ha usado la aproximación V 1 + e ~ 1 + e/2. Por
la ó rb ita es T0 el pe rio d o real es ejemplo, en el caso de la Tierra e = 6 x 10-8; en el de
Mercurio e = l,5 x 10-7.
T T / 1277^ 2. Como T < T0, el planeta avanza en el mismo
tiempo T0 un ángulo 27t -1- Ap sobre la misma órbita, con
" 0 i 1“ )' A<p igual a su velocidad angular por la diferencia de pe­
(Sugerencia: Tratar la corrección einsteniana com o riodos
si fuera un pequeño increm ento del valor de G.)
1277-V2
ar = r A- r =
2. D em o strar que, en cada revolución , el plane­ T1 oe 2
ta avanza en ángulo

104 105
100 problemas de Mecánica Movimiento planetario
 por lo que avanza un ángulo
6
vAT 2477V Teoría del potencial
A<P = ------= ~ 2^ 2 '-
3. En el caso del planeta Mercurio: r = rM= 0,39 UA
(UA = unidad astronómica = radio de la órbita terres­
tre). Su periodo se puede calcular en función del de la
Tierra usando la tercera ley de Kepler

- 0,393/2 = 0,244,

o sea T0 = T m = 0,244 años = 7,695 x 106 s, de donde

Atp = 4,77 x lfr 7,
Problem a 48
(en radianes). Cada siglo, Mercurio da un número de
vueltas N = 100/0,244, por lo que el avance total es
¿ ^ =^ = 40" Una em presa de tra n sp o rte s ha c o n stru id o un tú ­
nel rectilíneo que une dos pu n to s de la superficie
que está muy cerca de los 43" que dan las observaciones de la tie rra a través de su in te rio r. Escribir y resol­
astronómicas. ver las ecuaciones de m o v im ie n to para un vehículo
Como la órbita de Mercurio es elíptica, este efecto se que se m ueve en el in te rio r del tú n e l y calcular la
manifiesta en un avance de su perihelio, fenómeno inex­ duración del viaje. ¿Cuál será la velocida d m áxim a
plicable por la física clásica y uno de los grandes éxitos que se alcanzará en el túnel? Sí la aceleración m á­
de la relatividad general de Einstein. U n tratam iento xim a que puede so p o rta r el cuerpo h u m ano sin re­
más exacto da una aproximación aún mejor al valor ob­ sentirse es de unos 3 g, ¿presentará algún p e lig ro
servado. en este se n tid o v ia ja r en el túnel?

106 107

100 problemas de Mecánica

 Vamos a calcular el campo en el interior de la tierra. dio del segmento se obtiene por consideraciones geo­
Por el teorema de Gauss el efecto de una capa cerrada métricas elementales que
S olución exterior de masa es nulo con lo cual, al ir descendiendo
hacia el interior de la tierra sólo hay que considerar la x 0 = R t sen c¿
parte de masa que queda bajo nosotros, entonces
siendo a el ángulo que subtiende e l vector de posición
GM (r) 4
g = - — — = ~ y 7rGPor de uno cualquiera de los puntos c o n la línea que pasa
por el origen y es perpendicular a la cuerda A B que une
ambos puntos (que corresponde a l máximo valor de 8).
Como el valor de la gravedad en la superficie es Si en t = 0 el vehículo se encuentra en la superficie de la
4 tierra <p0 = 0 .
g 0 = — 7rpQG Rr donde R T es el radio de la tierra, basta
Curiosamente la expresión e n c o n trad a para la fre­
sustituir para encontrar que cuencia nos indica que ésta es independiente de la dis­
tancia entre los dos puntos en consideración, esto es, se
h ,
tardaría el mismo tiempo en atravesar la tierra a través
£ So g del ecuador que a través de un agujero de unos pocos
metros cerca de uno de los polos. Calculemos el valor
La componente de la fuerza en la dirección del túnel numérico del tiempo de viaje, que será la mitad del pe­
es g(r) sen 8. Definamos la coordenada de desplaza­ riodo del movimiento,
miento a lo largo del túnel como x = r sen 8. Como el
movimiento en la dirección perpendicular al túnel no es
posible debido a las paredes, el problema es unidimen­ T - —- = 2,53 x 103 s = 42 minutos
10
sional. Las ecuaciones de Newton para el sistema son
Para encontrar la velocidad m áxima derivamos la ex­
m x - -m g 0 - i - presión de la posición
T

o)(t) = x0w sen ojt

Que puede reescribirse en la forma más familiar

x + co2x = 0 El máximo está en í = ttIoj, luego

que es la ecuación de movimiento de un oscilador armó- = V g^sena

A g
nico de frecuencia w = \ l - ^ - . La solución general de
V R j
que sí que depende de la posición inicial. En el caso más
esta ecuación de movimiento es favorable, cerca del ecuador, en el que sen a - 1 se en­
cuentra
x(t) = x0 cos(wt + <p0)

Si tomamos el origen de coordenadas en el punto me­

108
 por lo que avanza un ángulo 6
A vA T 247r3r 2 Teoría del potencial
(P~ r ~ c2T \ '
3. En el caso del planeta Mercurio: r - r M = 0,39 UA
(U A = unidad astronómica = radio de la órbita terres­
tre). Su periodo se puede calcular en función del de la
Tierra usando la tercera ley de Kepler

T / \ 3/2
_K_ ^ / Zal J = o,393/2 = 0,244,

o sea Tq= T m = 0,244 años - 7,695 x 106 s, de donde

A(p = A,11 x 10^7,
Problem a 48
(en radianes). Cada siglo, M ercurio da un número de
vueltas N = 100/0,244, por lo que el avance total es
Una em presa de tra n sp o rte s ha co n stru id o un tú ­
Vsiglo = “ 40"
ne! rectilín e o que une dos puntos de la superficie
que está muy cerca de los 43" que dan las observaciones de la tie rra a través de su in te rio r. Escribir y resol­
astronómicas. ver las ecuaciones de m o v im ie n to para un vehículo
Como la órbita de Mercurio es elíptica, este efecto se que se m ueve en el in te rio r del tú n e l y ca lcular la
manifiesta en un avance de su perihelio, fenómeno inex­ duración del viaje. ¿Cuál será la velocida d m áxim a
plicable por la física clásica y uno de los grandes éxitos que se alcanzará en el túnel? Si la aceleración m á ­
de la relatividad general de Einstein. U n tratam iento xim a que puede so p o rta r el cuerpo h u m ano sin re­
más exacto da una aproximación aún mejor al valor ob­ sentirse es de unos 3 g, ¿presentará algún p e lig ro
servado. , en este se n tid o via ja r en el túnel?

106 107

100 problemas de Mecánica

 Vamos a calcular el campo en el interior de la tierra. dio del segmento se obtiene por consideraciones geo­
Por el teorema de Gauss el efecto de una capa cerrada métricas elementales que
S olución exterior de masa es nulo con lo cual, al ir descendiendo
hacia el interior de la tierra sólo hay que considerar la x Q= R t sen a
parte de masa que queda bajo nosotros, entonces
siendo a el ángulo que subtiende e l vector de posición
GM (r) 4
g = ----- — = - 3 de uno cualquiera de los puntos c o n la línea que pasa
por el origen y es perpendicular a la cuerda A B que une
ambos puntos (que corresponde a l máximo valor de 9).
Como el valor de la gravedad en la superficie es Si en t = 0 el vehículo se encuentra en la superficie de la
4 tierra <p0 = 0 .
= — 7Tp0G RT, donde R T es el radio de la tierra, basta
Curiosam ente la expresión e n c o n trad a para la fre­
sustituir para encontrar que cuencia nos indica que ésta es independiente de la dis­
tancia entre los dos puntos en consideración, esto es, se
*r tardaría el mismo tiempo en atravesar la tierra a través
g = ~8o•»
t\j- del ecuador que a través de un agujero de unos pocos
metros cerca de uno de los polos. Calculemos el valor
La componente de la fuerza en la dirección del túnel numérico del tiempo de viaje, que será la mitad del pe­
es g{r) sen 9. Definamos la coordenada de desplaza­ riodo del movimiento,
miento a lo largo del túnel como x - r sen 9. Como el
movimiento en la dirección perpendicular al túnel no es
posible debido a las paredes, el problema es unidimen­ T = —- = 2,53 x 103 s = 42 minutos
co
sional. Las ecuaciones de Newton para el sistema son
Para encontrar la velocidad m áxima derivamos la ex­
presión de la posición

(ü{t) - x0(ú sen (ot

Que puede reescribirse en la forma más familiar

x+ o j2 x =0 El máximo está en t= irla), luego

que es la ecuación de movimiento de un oscilador armó- Vmí,= w 4 Sen“

/§
nico de frecuencia o> = \ ¡ - ^ - • La solución general de que sí que depende de la posición inicial. En el caso más
V •*'7-
esta ecuación de movimiento es favorable, cerca del ecuador, en el que sen a = 1 se en­
cuentra
x(t) = XQCOS(ú)t + (pQ) *
vm. = 8 km s-1 = 30.000 km/h
Si tomamos el origen de coordenadas en el punto me­

108
 T=
r dty(r)
Al ser la aceleración de la gravedad en el interior de Para calcular el potencial usaremos que g(r) = ---- ~ ~
la Tierra siempre menor que en el exterior, no habría
riesgo biológico, si el autobús está refrigerado.
rr /s Gm
< £(r)= -J d r e ( r ) = - —
V - ?
Problem a 49

Una estrella tie n e p o r densidad Problem a 50

m a2 S u p o n ie n d o que el in te rio r de la T ierra es un flu i­

¡Ar) *
2 irr(r2 + a2)2 do in c o m p re s ib le en e q u ilib rio ca lcu la r la presión
d e n tro de la T ie rra a una distancia r del centro en
d o nd e m es su m asa y a es una constante con d i­ T cm -2.
m ensione s de lo n g itu d . H allar el ca m p o y el p o te n ­
cia! g ra v ita to rio en fu n ció n de la distancia r a l cen­
tro de la estrella.

S olución

Apliquemos el teorema de ¿au ss para una superficie

esférica situada a una distancia r del origen

<J>- 4Trr2¡r(r) = - 4rrGM(r)

Como se nos pide g(r) tenemos que calcular M(r),

pero

M{r) = f f f f(/)d V =
J J JsV) S olución
í*2ir i*w fr m a 2
=J d<j>J de s e n flj r 2— ~— ---- — dr =
Jo Jo Jo z - n r ( r 2 +■ a 2) E n el equilibrio la presión compensa la fuerza gravi-
tatoria en un elemento de volumen cualquiera del flui­
. rrur2 rr
mu r»- rr' ,, mr2
77127r-r Jn (r
jo '2+_i_a„22)2\2
(r'2 a2 + r2
do con lo que se tendrá
dF = g(r)dM - P(r + dr)dS - P{r)dS
Sustituyendo este resultado en la ecuación para el
campo se obtiene Como d M - pdrdS podemos eliminar dS en la ecua­
ción anterior con lo que nos queda
.. GM(r) Gm
g{r) = - dP = pg{r)dr
r2 a2 + r2

111
110
t. I -
 GM (r) GM GQ
y como por el teorem a de Gauss g(r) = se <f>(r) = -------------~t t ~ (3 eos2 0 - 1 )
i2 r 4 r3
obtiene
donde
4 77
f?Grd.r
Q = J J J d r 'p ( r ') ( 2 / 2 - Y 2 - j / 2)

Integrem os la expresión anterior para encontrar la

presión teniendo en cuenta que la presión en la superfi­ S olución
cie de la tierra es nula ya que no hay más masa encima
(P(r = R) = 0)
1. La densidad de masa en la barra es
rP f r 4'jT 4 TT
\ dP = - \ d r ~ ( ? G r = — -? G {R }-P ) M
J o R i 3 A
2L
donde hemos denotado por R al radio de la Tierra.
Para el centro de la Tierra, introduciendo los valo­ Utilizando la definición de potencial
res numéricos (p = 5,5 g/cm3, G = 6,67 x 1011 N m 2kg~2,
R = 6.370 km) se obtiene una presión

P(r = 0) - 1,72 x 104 T cn r 2

para esta geometría concreta tendremos

Problem a 51 GA 'W + L 2- L '

0 (p , 0) = - J -L dz' V/?2 + Z
- GA ln
C onsidére se una m asa M d is trib u id a u n ifo rm e ­ V p2 + U + L
m ente a lo largo de un segm ento co n te n id o en el
eje Z entre los pu n to s z = -L y z = + L
El campo se obtiene de g = -V<£. En coordenadas ci­
lindricas
1. C alcular el potencial y el cam po en el punto
(x, y, 0).
2. E studiar el c o m p o rta m ie n to del potencial le­ , , , , M 2GXL
jo s del s e g m e n to en un p u n to a rb itra rio (x, y, z) ^ P )S S (P )= ~ T P = - 7 W W
(hasta ios té rm in o s cu a d ru p o la re s) realizand o un

desarrollo en serie de potencias de — de la integral 2. A l a lejarn o s de la d istrib u c ió n de m asa

r = V p 2 + z 2 » L. Vamos a calcular el potencial
que define el potencial, siendo V x 2+ y2 + z2» L
3. C alcular d ire cta m e n te el potencia l haciendo <f)(p, z) =? G J dz'
uso del d e sa rro llo m u ltip o la r hasta orden 1/r3 -l V ? T J z - z ') 2

112
 ^ G A j- *L ^ r 1 Usar este resultado para calcular el potencial crea­
~ ~ * -l * V T + {zn - 2 z z ')!F do p o r S a tu rn o en un puntó de su eje p olar consi­
derando que el fa m o so a n illo de S aturno es un dis­
Esta integral puede calcularse de manera exacta pero co de radio in te rio r r, y de radio e x te rio r r2.
nos interesa el desarrollo multipolar. Para encontrarlo
desarrollamos el integrando

1 _ {zn - 2z z ) 3 (zn - 2zz')2

V 1 + {zn ~ 2 z z y r 1 ~ 2F + 8 rA

llevándolo a la integral se tiene

x GA [ 1 r *l
<Kr’ z) = — \ 2 L ~ 2 é L * '( * /2 - 2« ') +

3 r+L
+^ C dz'(z'2~ 2zz')2

y usando que eos 6 - — y 2L k ~ M se obtiene finalmen­

S olución
te el resultado
Para una distribución unidimensional de masa el po­
GM tencial es
(f>(r, Q )= - 1 + L _ (3 cos2 0 - 1 ) + o / L L '

3. En este caso sólo hay que calcular la integral

En nuestro caso la curva y es una circunferencia; así

' ) ¿ z ’2- x '2 -y"2) = Aj *L 2 z 'W =
Q = j j j d r ' p(r> pues el potencial en un punto del eje queda

con lo cual el desarrollo queda / n , AGr . GAr

ó ( 6 ,z ) = - \ d$ — = - 2 -n-— ..
J0 VA + Z2 V r 2+ d
GM GM L2 /
_ 3 _ ( 3 eos 0 —1) para calcular el campo simplemente derivamos

_ . . d<t> . lirG krz

Problem a EÁz) = - — = -
dz (r2 + z2)V2

C alcular el ca m p o gravita cio n a l que crea un a n i­ El potencial creado por el sistema planeta + anillo
llo (aro) con masa u n ifo rm e en un p u nto de su eje. consta de la contribución del planeta, que es la habitual,

114
7 / 7/7 n r n h l o m a e r í a R/l.
 y la contribución del anillo. Para calcular la contribu­ Para calcular la integral que define el momento cua­
ción del anillo hay que integrar el potencial anterior en­ drupolar lo mejor es utilizar coordendas elípticas que se
tre el radio interior y el radio exterior de la corona cir­ definen como
cular debido a la aditividad de los potenciales
.t = ak sen a cos /3
y = ak sen a sen /3
<f?{z) = f ' d r — = )
G, V r- + Z2 2 1 ' z = cÁ. cos a

con 0 < A < 1, 0 < a < ir, o < ¡3 < 2ir. La razón de elegir
donde a es la densidad superficial de masa de la corona
estas coordenadas es que la región de integración se
circular. La integral anterior puede entenderse tanto
simplifica enormemente como veremos enseguida. Ne­
como la suma de los potenciales individuales como la cesitamos calcular el jacobiano de la transformación
definición directa del potencial. Finalmente, el resulta­
do completo es 3 O, y, z )
d(a, (3, A)
4>{z) 7tGí7 ( V ñ + z2- V r i + z 2 )
z ak cos a cos (3 ~ a\ sen a sen (3 a sen a cos (3
ak cos a sen /3 ak sen a cos ¡3 a sen a sen /3
Problem a 53 ~ck sen a 0 c cos a

- a2ck2 sen a
C alcular el m o m e n to c u a d ru p o la r de una d is tri­
bución de m asa con fo rm a de e lipsoid e de re vo lu ­ Tras introducir el jacobiano la integral de volumen
ción de sem iejes a (doble) y c. queda de la forma siguiente:

Q = a2cp í d(3 í A4dA f da sen a (2c2eos2 a - a2sen2 a) =

Jo Jo Jo
S olución 4 77- ? 2
4 - a2 Cp (c2 - a 2) - ~ M{c2 - a2)
Usem os la definición de m om ento cuadrupolar de
una distribución de masa Según que el elipsoide sea achatado por los «polos»
(c < a) o achatado por el ecuador (c > a) el momento
0= fff p (? )(2 z 2- x 2- y 2) d ? cuadrupolar es positivo o negativo debido al signo del
J J J Elipsoide
último término de la ecuación anterior.
En nuestro caso la densidad p es uniforme y el volu­
men de integración es un elipsoide de revolución cuya
ecuación es

x2+ y2 z2 ^ La T ie rra puede c o n sid e ra rse en buena a p ro x i­

~ 02~ + ~d2 m a ció n c o m o un e lip s o id e de re v o lu c ió n achata-

116
100 problemas de Mecánica
 do p o r los p o lo s c u y o m o m e n to c u a d ru p o la r es ]_ £.2
2 se encuentra ( / 2)0 = —------- — . El prim er dato del
Q = - — M (a2 - c 2) siendo a el sem ieje ecuatorial
o problema es el valor de esta magnitud (/,)0, luego
y c el polar. Si la tie rra tuviese densidad uniform e,
G = - 2 Ma2(J2)Qr donde (J2}0 = 1,34 x 10-3. Sin em bar­ I i ?— = 1 , 3 4 x 1 0 ~ 3.
5 a2
go el va lo r observado de J2 es 1,1 x 10-3. Una posi­
Calculemos ahora el momento cuadrupolar de la Tie­
ble explicación de este hecho es co n sid e ra r la es­
rra considerada como formada por dos regiones de den­
tru c tu ra in te rn a de la T ie rra . S u p o n ie n d o que la
sidades distintas. Llamemos Q (x, y ; p) al momento cua­
T ierra está fo rm a d a p o r un corazón con fo rm a de
drupolar creado por un elipsoide de revolución de
elipsoide de revolución y densidad u n ifo rm e pc, ro­
deado p o r una capa de densidad u n ifo rm e pm y que
semiejes x e y y densidad p. Dada la dependencia lineal
del momento cuadrupolar con la densidad podremos es­
la fro n te ra entre am bos tuviera la m ism a fo rm a que
la Tierra, pero con sem iejes de lo n g itu d 4/5 los de cribir
ésta, encon trar qué cociente de las densidades p j p m
explicaría el m o m e n to cua d ru p o la r observado.
Q - Q(a, c ; p j + Q Í ~ a, — c; pc- p \ =

4 1 4 / 4 \ 24 1/42 42

5 c^ > 5 1=

p j c 2 - a2) + (pc - p j — (c2 - a2)

2 4 45
7ra2c(c2 - a2) Pm + l p ( P c - P j
T T

Por otro lado la masa total de la Tierra se escribirá en

función de las densidades de las dos capas como
S olución 43
M = — 7ra2c P,n+ p ( P c ~ P j
3
Como el momento cuadrupolar para una Tierra con
densidad uniforme es
Entonces, escribiendo el momento cuadrupolar en
función de los términos que aparecen en la definición de
q = _ 1 M {a2 - c2) = -2 Ma2(J2\ la masa, se tendrá

118 119
100 problemas de Mecánica Teoría del potencial
 43 estado de e q u ilib rio se sabe que la d istrib u ció n de
Q = - 2 — Trazc(a2 - c2) Pm + - p ( P c - P j
e q u ilib rio tiene sim etría radial y se conoce la te m ­
peratura en su in te rio r T(r). E ncontrar las ecuacio­
1 (1 - 4 5/ 5 5 ) p 45p I55 nes que d e term inan la presión en fu n ció n del p u n ­
xy Pm + 4'(Pe " P j /53 P„ + 43(PC~ P,„)/52 to P {r).

Como Q = -2M a2J2 se obtiejie para J2 ahora el siguien­

te valor S olución

1 c2 - a 2 1 - 4 5/5 5 45/5 5 En un problema anterior vimos que la ecuación de

2= ~ 5 ~ tf~ 1 + 43(a - l)/5 3 l + 4 3( a - l ) / 5 3 equilibrio en el fluido es

donde a = p jp m. Si tenemos en cuenta que dP M M / \ G M {r)

= P(r)g(r) = ~P(r) — ----
1 c2 - a2
5 = ( / 2)0 = 1,34x10-^ si el gas es ideal satisfará

PV= nRT
Y Q116 72 = 1,1 x 10~3. Definiendo y = JJ(J2)Qobtenemos

1 _ 45/55 ■ 45qt'55 siendo R la constante de los gases perfectos. Reescriba-

mos la ecuación en la forma más conveniente
7 “ 1 + 43(a - l)/5 3 + 1 + 4~3(o: - l)/5 3

y despejando a se encuentra P (r)= ^ p (r)T (r)

1- y
D espejando p(r) y sustituyendo en la ecuación de
43y/ 53 - 45/5 5 equilibrio se obtiene

Poniendo el^alor de y =
1,1 = 0,82 se obtiene final­ 1 dp <*A »,/ ^
M{r)
P(r) dr RPT{r)
mente
como la distribución tiene simetría esférica
= 2 ,9
M{r) - 4p í r'2p(r')dr'
Jo
o sea que

1 dP 4 ttGA
C o n s id é re s e un gas p e rfe c to de m o lé c u la s de
r / 2P( /) d /
P{r) dr Rr2T(r) * r,
peso m o le c u la r A. S u p o n g a m o s que en un cie rto

120 121
100 problemas de Mecánica Teoría del potencial
y derivando respecto a r e n los dos miembros se en­ ¿Cuál de los d o s m o d e lo s es m ás re a lis ta , de
cuentra acuerdo con lo expuesto en el enunciado? Consi-
, 1
derar b - ~ .
O
1 * ) =^ G R ^ r) - 4 ^ 1
dr y AP(r) dr J RT(r)

poniendo todas las constantes en un mismo miembro y

definiendo S olución

^ 4 ttA 2G La velocidad orbital en una órbita circular situada a

~~W ~ una distancia r en un campo de fuerzas gravitatorio pue­
de obtenerse a partir de las ecuaciones de Newton
la ecuación queda
GM {r)m v2(r)
p —----- .------= yyi--------
d í r 2T(r) d P \ P(r) A
dr \ P(r) dr J T{r)
y resulta ser
Así pues, conocida T{r) es posible integrar la ecua­
ción diferencial para encontrar P(r). GM{r)
v(r) =

Consideremos el primer caso: en el exterior es eviden­

En el estu dio de la dependencia de la velocidad te que v(r)= , para el interior al ser M(r) = ~ rrpGP
o rb ita l de las estrellas en to rn o al centro de una ga­
laxia suelen encontrarse a n om alía s que hacen ne­ se obtiene que v(r) = irpGr, en conclusión
cesario p o s tu la r la existencia de «halos» de m ate­
ria o sc u ra . En p a rtic u la r se e n c u e n tra que la
v e lo c id a d o rb ita l de las e s tre lla s d e p e n d e m u y GM
r> R
poco de la distancia al centro de la galaxia en un
d e te rm in a d o rango de distancias.
v(r) =
E studiar c u á re s la velocida d o rb ita l en las órbitas — rrpGr, r< R
c irculares para un m odelo de galaxia con densidad
u n ifo rm e en una región de radio f l y nula en el ex­
te rio r y co m pa ra rla con la velocida d o rb ita l que se En el segundo caso en el rango de valores interme­
o b tien e a ñadiendo un «halo» e x te rio r de densidad dios para la distancia al centro de la galaxia se tendrá

4 fr 4
p'(r) = bp , R < r< a fí M (r )= ~ tt R zp + 4‘7rr2p'dr = — 7rR2p+4TrbpR2( r - R ) =
siendo M. la masa de la región más interior de la estre­ O scilaciones
lla. De aquí la velocidad que se encuentra es

v(r) = y 4irGpbR2 + (i __3b), R < r< aR

En particular, cuando b = — obtenemos que la ve­

locidad orbital es casi independiente de la distancia
al centro de la galaxia en el rango comprendido entre
R y aR.
[Adaptado de K. V ander Lugt, «Invisible Galactic
Halos», The Physics Teacher, 31 (1993) 90.]
Problem a 57 j

Una fuerza p o r unidad de masa g0 cos(a>? + <p) ac­

túa so b re un o s c ila d o r a rm ó n ic o con a m o rtig u a ­
m iento su b crítico desde t = 0.

1. ¿C uáles d e b e n ser los d a to s in ic ia le s para

que no haya tra n sito rio ?
2. S i en t = 0 la p a rtíc u la estaba en re p o so
¿cuánto vale el tra n sito rio ?

S olución

1. Para un oscilador am ortiguado regido por una

ecuación diferencial
i
x + 2yx + to p = g0 cos (cot + <p)

la solución general es

(t)=A<r*cos(<ot+8)+ S° T r =r sen(wt+a^l3)
v i((úi-o) 2y + 4 Y (» z\
y derivando respecto a r e n los dos miembros se en­ ¿Cuál de los d o s m o d e lo s es m ás re a lis ta , de
cuentra acuerdo con lo expuesto en el enunciado? Consi-
, 1
derar Jb= — .
1 ( J W * * ) = -4-rrGRr2p(r) - - 4 * 0 d ™ O
dr y AP{r) dr f RT(r)

poniendo todas las constantes en un mismo miembro y

definiendo S olución

^ _ 4r f G La velocidad orbital en una órbita circular situada a

R2 una distancia r en un campo de fuerzas gravitatorio pue­
de obtenerse a partir de las ecuaciones de Newton
la ecuación queda
GAí(r)m v2(r)
/- = ----- .------=m --------
d l *T (r) d P \ _ P{r)
dr y P (r) dr j T{r)
y resulta ser
Así pues, conocida T(r) es posible integrar la ecua­
ción diferencial para encontrar P (r). GM{r)
v(r) =

Consideremos el primer caso: en el exterior es eviden­

En el e stu d io de la dependencia de !a velocidad te que v(r)= \ / ^ - , para el interior al ser M(r) = -í- irpCP
o rb ita l de las estrellas en to rn o al centro de una ga­
laxia suelen en contrarse anom alías que hacen ne­ se obtiene que v(r) = \¡ ~ ^ irpGryen conclusión
cesario p o s tu la r la existencia de «halos» de m ate­
ria o sc u ra . En p a rtic u la r se e n cu e n tra q u e la
v e lo c id a d o rb ita l de las e s tre lla s d e p e n d e m u y GM
r> R
poco de la distancia al centro de la galaxia en un R
d e te rm in a d o ra ngo de distancias. v(r) =
E stu d ia r cuarfes la velocida d o rb ita l en las órbitas TrpGr, r <R
circulares para un m o d e lo de galaxia con densidad
u n ifo rm e en una región de radio R y nula en el ex­
te rio r y co m pa rarla con la velocida d o rbital que se E n el segundo caso en el rango de valores interme­
o b tien e a ñadiendo un «halo» e x te rio r de densidad dios para la distancia al centro de la galaxia se tendrá

p'{r) = bp — , R <r<aR M(r) = — irR2p + J 4/rcr2p'dr = ~ TtR3p + 4TrbpR2(r - R ) -

 M.
1+ 3ih r x 7
siendo Mimla masa de la región más interior de la estre­ O scilaciones
lla. De aquí la velocidad que se encuentra es

v(r) = y 4'wGpbR2 + (1 - 3b), R < r <aR

E n particular, cuando b = — obtenemos que la ve­

locidad orbital es casi independiente de la distancia
al centro de la galaxia en el rango comprendido entre
R yaR .
[A daptado de K. V ander Lugt, «Invisible Galactic
Halos», The Physics Teacher, 31 (1993) 90.]
Problem a 57 j

Una fuerza p o r unidad de masa g0 cos(w f + <p) ac­

túa so b re un o s c ila d o r a rm ó n ic o co n a m o rtig u a ­
m iento s u b crítico desde í = 0.

1. ¿C uáles deben s e r lo s d a to s in ic ia le s para

que no haya tra n sito rio ?
2. .S i en t = 0 la p a rtíc u la estaba en re p o so
¿cuánto vale el tra n sito rio ?

S olución

1. Para un oscilador amortiguado regido por una

ecuación diferencial
4
x + 2yx + upe = g0 cos(a)t + (p)

la solución general es

125
donde A y 8 dependen de los datos iniciales y 0 es A
A = sen(/3 + a 0) Va>2 + (y + w cot(a 0 + 0))2
***y
j.3 = arctan = Vo>2- f
2yco
Problem a 58
Escribamos la solución en t = 0:
C alcular la potencia disipada p o r un oscila d o r ar­
x(0) = •sen(a:0+ 0) + A eos 8 m ónico a m o rtig u a d o y forzado en estado estacio­
V [(^¡~ co 2)2 + 4^a»2] nario d e bido al rozam iento. ¿Para qué frecuencias
& de la fuerza excitadora se hace m áxim a esta pér­
dida?
x (0 )= (ü cos(aQ+0)-Aa)se n 8~yA eos S
V[(ú>2- ctí2)2+ 4y2(o2]
S olución
Denotemos para simplificar A=g^V[(a )2- a)2)2+4y 2ai2].
Si queremos que no haya transitorio deberá cumplirse
A = 0 con lo que la solución está en régimen estaciona­ La energía disipada debido al rozamiento será
rio. Imponiendo esta condición encontramos inmediata­
mente *(0 ) y ¿( 0 ) E= J ^"rozamiento ÚX = 2m? J xdx = 2m y J X2dt
*(0) = A sen (a0 + 0) La energía disipada en un ciclo es
¿( 0) = wA eos (a0 + 0)
E - 2m y j x 2dt
:
2. Si en t - 0 la partícula estaba en reposo x(0) =
y la potencia será el promedio
= x(0) - 0 y entonces podemos calcular A y 8 despejan­
do de las ecuaciones

-Asen(a 0 + 0) = A eos 8
(¡i cos(a 0 + /3)A = A ío sen 8 + yA eos 8 Como el movimiento en el estado estacionario para el
oscilador armónico forzado viene dado por
Dividiendo la segunda de estas ecuaciones por la pri­
mera se obtiene el valor de 8
x(t) = ----- — °._ ■- — sen(wr+ a)
w V [(íü“- (o2)2 + 4 t 2íü2]
tan 5 = - cot (aQ+ 0)
tü
y la potencia queda entonces
Tom ando ahora la prim era ecuación y teniendo en
2y (ú2s 2 rT
cuenta que eos 8 = 1/ V l + tan 2 8 se encuentra el valor P = —=- ------- 2ry° „ J COS2(o>i + a) dt =
m T [(ft>2 - ít»2) 2 + 4 y - ( o 2] Jo
de .A

127
 ysl Cuando la fuerza está en resonancia de energía oí = <ú0
y entonces la ecuación que obtenemos para la amplitud
m - (o +4 f del movimiento es

El máximo de esta potencia con la frecuencia se ob­

tendrá de hacer dP/dco = 0, obteniéndose que oí = o>0. 2 ya> 1
Nótese que esta frecuencia es distinta de la frecuencia
de resonancia de la amplitud, que sí que cambia por y cuando la fuerza oscila con frecuencia w, = 2co0
efecto de la disipación y corresponde al máximo de la
función A(a>) = g jm V(w 2 - col)2 + 4y¿co2. So =a
w2V 9 + 16y2/ai2 2

Dividiendo las dos ecuaciones anteriores se encuen­

P ro b le m a ® :
tra

Un m u elle a rm ó n ic o de frecuencia latera! a> osci- 2yo) A

„ la en el in te rio r de un flu id o viscoso. C uando se le
2V 9 + 16y2/ oí2
Oí■ A.
som ete a una fuerza externa de frecuencia oí = 1 s_1
que c o in c id e con la fre c u e n cia de resonancia en
energía (to0 = oí), el m o v im ie n to se realiza con am- 1 , /3
Introduciendo valores numéricos se obtiene y = — \ ¡—.
/p litu d A x = 10 m m . Si la frecuencia de la fuerza exter­
na se eleva a ^ = 2a>0 m antenie ndo la intensidad de Ahora que conocemos y y podemos calcular g0 toman­
la fuerza im pulsora la am plitud pasa a ser A , = 1 m m . do, por ejemplo, la primera ecuación
Se a um enta la te m p e ra tu ra de! flu id o hasta V para
re d u cir la viscosidad de m odo que la a m p litu d dél
m o v im ie n to en la resonancia pasa a ser A.' - 20 g0 = 2Al7(o = j - V — m m s-
m m . Si en ese m o m e n to (a T') desaparece la fuerza
externa, ¿cuánto tie m p o tardará en decrecer la am ­
Después de elevar la tem peratura y en resonancia se
p litu d de la oscilación del m uelle a m enos de 1 mm?
cumplirá

= A!
S olución 2'y'oí 1

y despejando se encuentra el valor de y', ya que las res­

Recordem os que la solución general para el movi­ tantes cantidades son conocidas
miento del oscilador armónico forzado es

So
x(f)=A£-1"cos(« í+5)h— =______________ sen(aít+a0+fi) 7 80 V2
7 V [ 0 2- « 2) 2+ 4 y ^ 2]

128 129
W n nrnhlfím as Hp Mpr.ánír.p
 Como querem os ver cuando decrece la amplitud a Para calcular la constante de proporcionalidad hare­
menos de 1 mm y el decaimiento de la amplitud es ex­ mos uso de que la probabilidad de encontrar al oscila­
ponencial dor para algún valor de x es 1 , o sea que
A'1e~YT = 1
r A e A C
1= P(x)dx - - ■■■ ■dx
y despejando se obtiene -A J ~A C ú \/Á 2 - X 2

y entonces
7, = ^ 7 ln A / = 3 9 ,ls
y

C J-A coVA2 - X 2

Entonces

C a lcu la r la d e n sid a d de p ro b a b ilid a d de que al

m e d ir la p o s ic ió n de un o s c ila d o r a rm ó n ic o se le P(x) =
rV A ^x2
encuen tre en la posición x.
Puede observarse que la probabilidad de encontrar a
la partícula en tom o a x = 0 es mucho más baja que la
probabilidad de encontrarla cerca de los extremos ya
que en esos puntos la velocidad tiende a anularse y la
partícula invierte mucho tiempo en esas zonas.
La posición en función del tiempo viene dada por

x(t) - A sen(íüt + 5)

La probabilidad de encontrar al oscilador situado en­

tre x y-,x + dx será proporcional al tiempo que tarde en
U n o s c ila d o r a rm ó n ic o de m asa u n id a d se en­
recorrer ese intervalo
cuentra so m e tid o a una fuerza F{x, t) inducid a p o r
una onda viajera de! tip o F{x, t ) = f ( t - x f v ) , siendo v
P(x)dx «: dt .
una constante con dim e n sio n e s de velocida d.

dt 1 1 1
P (x) 1. S u p o n ie n d o q u e = f { ~ = 0 y q u e la
dx x A ío cos((ot + S) A to V 1 - x 2/a2 desviación de e q u ilib rio es una ca n tid a d pequeña
encuéntrese la relación que existe e n tre la fuerza f,
O sea que la c a n tid a d de m o v im ie n to A P y la e n e rg ía A E
tra n sfe rid a s al o scila d o r p o r la on d a cu a n to í - >
p (X) = — - £ = - 2. ¿Qué ocu rre si la fuerza excita d o ra tie n e m e­
coVA2 - x 2 dia nula?
 El problema está motivado por la interacción de un y sustituyendo en la expresión para AP se obtiene
oscilador armónico clásico, que puede representar la
S olución respuesta de un electrón de la materia en una primera r +" AE
aproximación con una onda viajera que podría ser la del AP = f(t)d t + ------
v
campo electromagnético de una onda luminosa, así pues
se trata de un modelo sencillo de interacción luz-mate­ 2. En el caso de que la fuerza excitadora tenga me­
ria. Para más información véase [J. M. Cabrera, F. T. dia temporal nula
López, F. Agulló, Óptica electromagnética, Addison-
Wesley (1993)]. y ~ m d t= o
1. La ecuación del movimiento para la desviación
del oscilador respecto a la posición de equilibrio es la de se tiene en el orden de aproximación analizando
un oscilador forzado

x = (afa = f ( t - x /v ) v

La energía transmitida por la onda al oscilador ven­

drá dada por la fórmula
f +00 í* +00 f +00
AE = J F(x(t), t)dx - J F(x(t), t)d t = j f ( t - x/v)xdt Un o scila d o r a rm ó n ico con a m o rtig u a m ie n to se
ve s o m e tid o a una fuerza tra n s ito ria f(t) = fe-»* de
Como x, que mide la desviación, es una cantidad pe­ fo rm a que su ecuación de m o v im ie n to es
queña, podemos realizar el desarrollo de Taylor d e /e n
tom o a x = 0 lo que nos da X + 2 y x + iú\ x =

A nalizar su m o v im ie n to su p o n ie n d o que y 2 > &>*.

f{ t-x ív ) = /( f ) - /( f ) T

Quedándonos hasta términos de primer orden e n r y x S olución

la energía trasmitida será
La solución general para el movimiento sin fuerza ex­
AE = \ ~ f { t ) x d t terna (ecuación homogénea) es

Para la cantidad de movimiento (impulso) se obtiene x(í) = A ePí‘ + Be&

donde y fi2 son las soluciones de la ecuación caracte­

A P =J F(x, t)dt = J f{t)dt ~ — J f(t)xd t rística. Al ser el amortiguamiento supercrítico y no ha­
ber oscilaciones /31 y ¡32 son constantes reales negativas.
Integrando por partes
A = - 7 + V y r “ <u2
J f{t)xd t - -J f(t)xdt = AE P2 = - r - v 9 ^

132 133
7 / 7/7 n r n h l o m a t : r i o F l A a n á n í / ' o
 Busquem os una solución particular de la ecuación En la configuración inicial el alargamiento del muelle
completa de la forma x {t) - Xer^. El probar esta solu­ cumple
ción viene sugerido por el hecho de que la derivada de S olución
la exponencial es ella misma, lo cual nos hace esperar M g = k{d ~ l0)= $ d = M g ljk + lQ
que xp(t) sea solución si elegimos convenientemente el
coeficiente X . Sustituyendo se obtiene que Denotemos y, x las distancias al punto de suspensión
inicial O, situado en el techo, de las masas M y m res­
pectivamente. Aplicando la segunda ley de Newton se
¡¿ - 2y¡x + oí\
obtiene de m anera inmediata que

La solución general para ¡j? - 2y¡x + to2 * 0 será d 2x

m g + k(y - x - Q = m
dt2

x(t) = A e l-y+v7I-'iao)' + B e{' y- Vy2-*$‘ + —----------------er* d 2y

¡j? - 2 y¡x + íü2 Mg - k(y - x ~ L) = M
~dt2~
Así pues en un caso genérico (con p} - 2y¡x + 0)
el movimiento tiene otra componente de amortiguación Sumando las dos ecuaciones anteriores se encuentra
que será un transitorio más o menos rápido dependien­ que
do de como sea ¡x en relación con -y ± V y 2- co2.
d 2x , d 2y

que puede escribirse como

Un m u elle ideal de m asa despreciable, constante
re cup erad o ra k y lo n g itu d n a tu ra l l0 u n id o en sus d 2 i mx + My \ _
e xtre m o s a dos m asas m y M puntua les se encuen­ dt2 m +M I ^
tra susp e n d id o del te ch o p o r el e x tre m o de masa
m. En t = 0 se suelta el m uelle, se p ide analizar el así pues, el centro de masas, cuya posición viene dada
m o v im ie n to resultante. m x + My
por X -------- —— se mueve en caída libre. Como ini-
m +M
m
cialmente j ( 0) = y( 0) = 0 , *( 0) = 0 , y( 0) = d, entonces
X(0) = ~~rz y X(0) = 0. Integrando dos veces la
M +m
k ecuación de movimiento para el centro de masas se ob­
tiene

2 m+M
M

135
n~~:i
 Análogamente se encuentra operando con las ecua­ la masa m tiene un desplazamiento que es siempre ma­
ciones x e y que yor o igual que su desplazamiento libre mientras que a
la masa M le ocurre justo lo contrario.
d 2y d 2x • [Adaptado de P. Glaister, The Physics Teacher (1993)
239.]

, , , k(m + M)
donde Defmiendo z = y - x ~ l Q, que es
Problem a 64
el alargamiento del muelle respecto a su longitud natu­
ral, la ecuación anterior se escribe E scribir las ecuaciones de m o v im ie n to para los
sistem as de la fig u ra .
d 2z
+ ú)2z = 0
dt2

que es la ecuación de un movimiento armónico simple.

Usando los valores de las condiciones iniciales x(0) = 0,
y(0) = d se encuentra z(0) = d - lQ= M g ljk , y como
z (0 ) = 0 , entonces la solución es

z(f) = —1— eos (üt

rC
A la vista de los resultados o b te n id o s d is c u tir si
observemos como el alargamiento del muelle depende es posible co n stru ir, contand o con dos m uelles de
tanto de la longitud del mismo como de la fuerza de constantes recuperadoras fc, y k2 un d is p o s itiv o que
gravedad. Usando el resultado para el movimiento del recupere elásticam ente con una frecuencia natural
centro de masas y la oscilación en torno al mismo pode­ de o scila ció n m e n o r que la de los c o m p o n e n te s.
mos calcular las posiciones de los extremos S upone r en to d o m o m e n to que las lo n g itu d e s na­
turales de am b o s resortes coinciden.
m x + M y - ( m + M) - + Md
S olución
7 Mgl0
y - x = L + — —1 eos o)t
u k
1. La fuerza total que actúa sobre la masa será la
suma de las fuerzas producidas por los dos muelles, o
encontrándose finalmente que
sea que si el desplazamiento respecto al equilibrio es x
la fuerza recuperadora es
2 k(m + M) F = F1 + F2 ~ - k xx ~ k rx ~ ~(kl + k2)x
y(t) = d + ^ - - p M É s ^ ( i - COSú,t) La ecuación del movimiento es inmediata
2 k(m + M ) '

136 137
 K—
.* + —* + -k2 n
2—x = O
m con la recuperación (ya que w2 = —), tenemos un mon-
m
taje con dos muelles elásticos que da como resultado
y entonces la constante equivalente del sistema es k= k1+k2.
una recuperación efectiva menor que la de cada uno de
sus componentes.
2. En este caso la elongación estará relacionada con
la fuerza total a través de F - -k x , siendo k una constan­
te equivalente al conjunto de los dos muelles. Como la
fuerza que actúa en cada uno de los resortes es la misma Problem a
al igual que la tensión de una cuerda es la misma en to­
dos sus puntos, los desplazamientos de cada uno de los Una masa sujeta p o r un m uelle elástico es sepa­
resortes deberán satisfacer F= kxx- kx2, entonces rada de su p o s ic ió n de e q u ilib rio y llevada a un
p u n to dista n te x0 de la m ism a, en el que se la suel­
ta sin ve lo cid a d inicial. A nalizar el m o v im ie n to su b ­
F F
siguien te de la partícula su p o n ie n d o que la super­
fic ie so b re la que se m u e ve la m is m a es rugosa
y p ro p o rc io n a una fuerza de ro z a m ie n to del tip o
y como x - x í + x2 F= ¡xmg, siendo p. el coeficiente de rozam iento.

F = -k x = -k(x, + x7) = -k F í +-i

\ K K S olución

por tanto la constante equivalente es Como la fuerza de rozamiento es siempre contraria a

la dirección del movimiento, que viene definida por el
L - L J_ sentido de la velocidad, la ecuación de movimiento será
k k , + k.
m x + kx = -fim g sgn(i)
y la ecuación del movimiento
*( 0 ) = x 0

x- + -----
1 — KK
x =0 Resolver directamente esta ecuación diferencial es un
m k x + k2 poco engorroso. Vamos a evitar hacerlo usando en cam­
bio razonamientos físicos.
Supongamos que el muelle de mayor constante recu­ Desde el m om ento en que se suelta la partícula en
peradora es el 1, entonces x = *0 hasta que esta se detiene el movimiento tiene lu­
gar con velocidad negativa (el muelle está recuperado)
■■ , K.....
K = k2
i — —Kk— < k; y entonces la ecuación en este tramo es
k x + k2 kx + k2 1
x + <o2x = p.g
Así pues la constante recuperadora del conjunto es
menor que la de cada uno de los dos muelles individual­ donde w2 = klm. Busquemos una solución de la forma
mente. A l estar la frecuencia relacionada directamente * ( 0 = ■^'homogéneíi + Y teniendo en cuenta que
xp = fxgl &>o = cte. es una solución particular, encontra­
mos
y la velocidad
x(t) = ( x - j cos <Ot +
x(t) = -w 0 ( -X0 + 1 sen <oQ£
donde ya se ha tenido en cuenta que la condición inicial
es *( 0) = xQ.
Esta solución es válida hasta que la velocidad se haga 7r
De nuevo la partícula se vuelve a parar cuanto tp2 = —,
nula, que es cuando cambiará de sentido, y como
y la elongación es

4 ¡xg

el instante de parada es tp¡ = ttIcoq. En este momento la

partícula estará en Así pues, tras cada semiperiodo de oscilación la am­
plitud del movimiento decrece en Ax = 2fxg!oi\. El movi­
miento continuará siendo así hasta que la fuerza elástica
ejercida sobre la partícula en uno de sus puntos extre­
mos (x¡) sea menor, en valor absoluto, que la fuerza de
La acción del rozamiento se ha traducido en que hay rozamiento.
un decaimiento de ía amplitud ya que el oscilador no al­
canza el punto ~-xv Después de llegar a la posición de la ¡x g > c o j jx j
primera parada el oscilador intenta cambiar el sentido
de la velocidad y desplazarse en sentido contrario; pero Cuando esto ocurra la partícula se quedará sin poder
esto sólo lo podrá hacer si la fuerza recuperadora es volver a invertir su movimiento.
mayor que la fuerza de rozamiento, esto es, si se satis­
face

(O * . - 2 >
0 Un b loque de m asa m m o n ta d o sobre una guía
h o riz o n ta l oscila b a jo la acció n de un m u e lle de
Supongamos que se cumple esta condición, entonces constante recuperadora k y lo n g itu d natural /0.
la ecuación de movimiento será
1. A nalizar el m o v im ie n to del b loque partiendo
x + ídJ * = ¡xg de la posición de m áxim a elo n g a ció n su p o n ie n d o
que el a la rg a m ie n to es pequeño y que puede desa­
rrollarse la fuerza elástica en serie de potencias.
con la condición inicial x(0) = - x Q+ — —. La solución de
2. En la m ism a a p ro xim a ció n , e n contrar la de-
esta ecuación diferencial es
 pendencia de! p erio d o en la a m p litu d y com parar La ecuación de movimiento del bloque es entonces
con el caso del o scila d o r a rm ónico.
d 2x k
r - 7
+ ----------- - X3= 0
dt*• 2

Como partimos de la posición de máxima elongación

en t = 0

dx
= 0
X = A ' ~di

Introduzcam os ahora las variables adimensionales:

X - x/A, t = tV k /m y el parámetro adimensional b = A//0.
La ecuación de evolución toma la forma

d2X
+. _b l 3. 0
*- 3=
S olución dr2

1. La deformación del muelle es con condiciones iniciales * ( 0 ) = 1, *( 0 ) = 0 - Para inte­

grar el problema podemos hacer uso de la conservación
A/ = W T / 2 - / 0 de la energía

y la fuerza restauradora sobre el bloque X2 b2 d _ b2

2 8 * 8
F = -k A l sen a
de donde, despejando d x fd res fácil deducir que
donde

x 2 f* ' 1
sen a - t= — — - ¿¿y
\ / x tT1Í

y entonces Haciendo las sustituciones x = (1 ~ ^ 2)1/2 y z = sen se

obtiene
F = -k x 1-
211
7
- T + 1
V2 J ( 2 ,<P)

142 143
100 problemas de Mecánica Oscilaciones
 elípticas. De todas formas, independientemente del va­
donde lor de la integral (que es un número), la expresión ante­
rior nos da la relación funcional entre T y A , pasando a
H k,ip) = j V~— 1 =^ d<p' las variables con dimensiones tendremos
Jo V I - k sen 2 <p

es la integral elíptica de primera clase en forma normal

V2 A V k 0
de Legendre. En nuestro caso k - ^ -. Las propieda­
así que hay una fuerte dependencia del periodo en la
des de las funciones elípticas son bien conocidas [ver amplitud del movimiento inicial, al contrario de lo que
por ejemplo: L. Abellanas, M. Spiegel, Fórmulas y ta­ ocurre en el caso del oscilador armónico en el que el pe­
blas de Matemática Aplicada, Serie Schaum, Me. Graw- riodo no depende de la amplitud de movimiento.
Hill, Madrid (1993)]. Usando la inversa de % que es la [Adaptado de V. Detcheva, V. Spassov, Phys. Educ.
función amplitud am obtenemos 28 (1993) 39].

y como

/ / br \ \ / r \
z = sen <p= sen | am | —-¡=- j j = sn ( ¿>

y para x se obtiene finalmente

A A bT \ b?
x w = V i-” 2 w r n v r

A, / k
x (0 = A c n r - y ^
V o

2. La expresión para el periodo es

— L \ l L . =p díp = 1,8541
V 2 4 Vm 0 V i -A sen2 <p

El valor de la última integral puede ser calculado nu­

méricamente o acudiendo a unas tablas de integrales
 8
O scilaciones pequeñas en sistem as
con varios grados de libertad

Problem a 67

Dos péndulos sim ples de lo n g itu d / y masa m es­

tán unidos m e d ia n te una cuerda elástica de cons­
tante recuperadora k. La cuerda está atada en am ­
bos p é n d u lo s a una d is ta n c ia a de! p u n to de
suspensión. La lo n g itu d natural de la cuerda elásti­
ca es igual a la distancia entre los pu n to s de sus­
pensión de los péndulos. Calcular las frecuencias,
los m odos y las coordenadas norm ales de las os­
cilaciones en to rn o al p u n to de e q u ilib rio del sis­
tem a.

________________ 147
 C o n s id e ra r los p é n d u lo s fo rm a d o s p o r v a rilla s
in deform ab les. = -i-k [V ( a cos (p2 - a cos <px) 2 + (lQ+ a sen <px - a sen <p2) 2 - /()]2

y entonces el lagrangíano exacto del sistema es

S olución 1 1
L ~ — mPcfí h mPi¡P2 - mgl(2 - cos <p, - cos <p ) -
Para escribir las ecuaciones de movimiento necesita­ -2 2
mos construir el lagrangíano y para ello la energía po­
la
tencial, que será suma de la gravitatoria más elásti­ la i.-:
\ __
k[V (a cos <p2- a cos <pjl + (l0+ a sen cp{- a sen <p2)2- /J 2
la
ca almacenada en cuerda elástica, que son

v gr - mSl 0- ~ cos <Pi) + ~ cos <P2) El único punto de equilibrio es <px =<p2 = 0. Vamos a
aproximar el lagrangíano en el entorno de este punto de
equilibrio hasta términos de segundo orden. Para ello
Vd = T k ( d - l af vamos a calcular independientemente las distintas con­
tribuciones
Donde d es la longitud de la cuerda elástica, y /0 su
longitud natural. ^
cos (p = 11 - — sen <p= cp
La distancia d puede determinarse en función de los z*
ángulos por medio de la expresión
Ifisl
mgl( 2 - cos <px - cos <p2) = -— (<p\ +

(‘V ( a cos tp2- a cos tpj)2 + (/0+ n sen <px- a sen<p2)2 - Q 2=

= li + ^4 ~<p0 + (lo + a<Pi - a<p2y -

- 2 lQV (a cos <p2- a cos <p])2+ (l^+a sen <px- a sen <p2)2 -

= P0 + P0 + a2íp\ + a V 2 + 2al0<px - 2alQcp2 - 2a2<p^<p2 -

- 2 PQ~ 2 /0a<p1 + 2alQ(p2 = a2{<p\ + <fP2) - 2a2q>l(Pl

Usando estos resultados, el lagrangíano queda

cP = (a cos <p2 - a cos (p:)2 + (/0 + a sen <px- a sen <p2)2 r 1 12 ■2 i 12 o / ° 2k\

=*Uel = - k { d - l , ) ^ ( mgl a2k \

~i — + - ^ W 2 + k a i Pl<
P2

148 W 149
100 problemas de Mecánica Clsnilaninnes nenueñas
que ya es de la forma L = T - V con Los modos normales son los autovectores correspon­
dientes a cada uno de los autovalores

f - T * *(? i ka2 - ka2

\- k a 2 ka 2
~ 1 . . / mgl + ka2 - ka 2
=y ( ^ ^ ^ka2 mgl + ka2 / ka2 ka 2 ^ f <p,
O).
2 ' ' ka 2 ka1
Este resultado también puede obtenerse directamen­
te haciendo Hemos encontrado los autovalores y autovectores,
que corresponden a las frecuencias normales y los mo­
y= Y J 2 L dos normales de oscilación. El significado físico de estos
d(pd<p¡ modos normales es que son formas de oscilar del siste­
(0,0)
ma que tienen una única frecuencia bien definida por
que es la expresión de la aproximación cuadrática del corresponder a un movimiento al unísono. En este caso
potencial, pero suele ser más sencillo por medio del de­ los dos modos normales corresponden a las dos formas
sarrollo en serie siempre que las funciones sean elemen­ de oscilar que se muestran en la siguiente figura.
tales. Las frecuencias normales satisfacen

det (Tai2 - V) = 0

donde T yV son las matrices de la forma cuadrática que

define las energías cinética y potencial. Explicitando
para nuestro caso la ecuación que queda es

/ mgl + ka2 - oí2m i2 - ka2 \

y - ka2 mgl + ka2 - tdm P j
Con esta estructura de oscilación las frecuencias obte­
Desarrollando el determinante nidas para los modos normales resultan naturales. El
número de n^pdas.normales que se obtienen en un siste-
(mgl + ka 2 - mPür)2 - k 2a4 - 0 => j n a debe coincidir además con el número dé grados de
. libeHaditelsisféma.
=> tó4m2/4 - 2mod[mg + ka1) 4- mgl [2ka2 + mgl] = 0
Por otro lado, las coordenadas normales son, por de­
mPcü2 = mgl + ka 2 ± ka2 finición, aquellas en las que el lagrangiano se escribe

Las raíces positivas de esta ecuación bicuadrada son

las dos frecuencias normales, en este caso

_ \[ s _ \i 8 . o sea que diagonalizan simultáneamente las matrices T

^ V l ^ V ¡ m i2 y V [con esta elección seguimos el criterio más corriente
en la elección de las coordenadas normales, ver por segunda corresponde a un movimiento de los dos en
ejemplo A. F. Rañada, Dinámica Clásica, pp. 524, oposición de fase.
Alianza Editorial]. Los autovectores nos definen un sis­
tema de coordenadas en el que las energías cinética y
potencial son diagonales pero los coeficientes de la dia- , Problema 68
gonal en la energía cinética no son unos, como requie­
ren las coordenadas normales. Para obtenerlas hay que A nalizar las co n trib u cio n e s a la energía potencial
normalizar los autovectores a la energía cinética, esto es, producid as p o r un m u e lle que in icia lm e n te fo rm a
un á n gulo a a rb itra rio con el eje X y cuyos extre­
X 'fX = 1 m os son apartados ligeram en te de su posición de
e q u ilib rio en d ire ccio n e s a rb itra ria s. P articularizar
E n nuestro caso la condición anterior implica que para el caso en que la lo n g itu d de e q u ilib rio c o in ci­
da con la n a tu ra l. U tiliza n d o el re su lta d o de este
=1 p ro b le m a d is c u tir c u a lita tiv a m e n te qué o cu rre en
■>("f i ) (!' el caso de oscilaciones trasversales del sistem a de
la figura.
y de aquí N j = 1/V 2 m i2. Haciendo lo mismo se obtiene
la norma del otro autovector, que resulta ser N 2 ~ N y
Las coordenadas normales serán entonces

= Á 'f q

1 1 1
V 2 mP \1 -1 )

escribiendo explícitamente las matrices S olución

Llamemos r, y r2 a las posiciones iniciales de equili­

I 7] I V 2 m l2 \ 1 brio de los extremos del muelle. Evidentemente la lon­
gitud de equilibrio del muelle es le = ¡| r2 - rj|. Tras aña­
y de aquí dir unos pequeños desplazam ientos obtenem os unas
posiciones para los extremos de los muelles que son
m i2 —- y —» ^

(?! + <P2) r i = ri + ei 1
—$r ~* l
r2 = r2 + *2 )
m i2
(<Pl <P2) Entonces la energía potencial es

La primera corresponde al movimiento acoplado de

los dos péndulos al unísono en la misma dirección y la

153
 brio en la que la longitud coincida con la longitud natu­
ral lQ= le, la energía no depende de la longitud natural.
Es importante hacer hincapié en que esto sólo ocurre
para el caso de lQ= lg. En este caso además la energía
¿Yk 3T depende sólo de la proyección del desplazamiento sobre
la dirección de equilibrio. Si hacemos desplazamientos
/ 'C , sJ , perpendiculares a la dirección de equilibrio la contribu­
ción a la energía es nula.
Con este resultado podemos analizar el sistema de la
figura. Nos piden que estudiemos las oscilaciones tras­
versales, que son las que se producen cuando los despla­
zamientos son perpendiculares a la dirección de equili­
brio. En nuestro caso la longitud de equilibrio coincide
= T fcdi ?2 + ?2- ? , - ? 1ü - / 0)2 = con la natural. Estas condiciones llevan a que no hay
ninguna contribución cuadrática al potencial. Las con­
tribuciones de orden más bajo al potencial serán de or­
den superior y las oscilaciones pequeñas serán anarmó-
nicas.

Definamos ahora r = r2 - r, y e = e2- e,, entonces

1 1
V = ^ h ( \ \ f + t I I - y 2 = y í t ( V ^ + 2 ? -r + ? - V =
A lg u n a s m oléculas tria tó m ic a s co m o la de C 0 2,
están fo rm a d a s p o r dos masas iguales m ligadas a
una tercera M situada en el p unto m edio de la línea
que une las m asas m. C alcular los m od o s y las fre ­
cuencias n o rm a le s de v ib ra c ió n s u p o n ie n d o que
Si desarrollamos por Taylor la raíz obtenemos
las fuerzas que las unen se pueden a p ro x im a r por
\2 fuerzas elásticas y que el m o v im ie n to se realiza a
V = - L * |/ 1+ -L lo largo de! eje de la m olécula.
2 /2

y para el caso en que L = le llegamos al resultado

Y Y
Siendo u ~ — = -ñ^n- • Pues’ cuando realizamos des-
K IMI
plazamientos pequeños desde una posición de equili­
 Las energías cinética y potencial en coordenadas refe­ 1 \ /l\ / i \
ridas a las posiciones de equilibrio son
0 x ,=
2 i X .3 = __ 2M
m_
S olución
1 1 1 -l) U / \ 1 )
T = — m x\ + — M xl + — m x\
2 1 2 2 2 3
M odo 1

V = ~ k(xl - x2y + 4 k(x 2 ~ x2f

A partir de aquí es fácil ver que las matrices T y V M odo 2 M odo 3

son

'm 0 / k -k 0 \
° \
T= 0 M V= ~k 2 k ~k
0 0 J
m
10 lo -k */ y si normalizamos a la energía cinética (X ,'TX¡ = 1) los
modos normales quedan
Las frecuencias normales se obtienen, como es habi­ j 1 \
tual, de la ecuación 1 1 1 \

X ,= ,___ 0 x: - 1
1
V 2 tn * VM + 2m
k -m o ? -k 0 i i /
det | -k 2k ~M (ü 2 ~k
1
i 1 \j
M
0 -k k - mu?
m
2m 1 +
= ( k - mo?) [{k - moj2) (2k - Ma?) - 2k2] - 0 ~2M

El modo X 2 corresponde al movimiento de todas las

De aquí se obtienen las tres frecuencias partículas que recorren la misma distancia en la misma
dirección y se produce con frecuencia cero, lo cual indi­
ca que éste es en realidad un movimiento de traslación
sin vibraciones del conjunto.
[A daptado de H. Goldstein, Dinámica Clásica, Ed.
a>2 = 0 Reverté (1987), pp. 323 y ss.]
2m + M
w3 = \ ¡ k
mM

Los modos normales son los autovectores correspon­ H allar las frecuencias norm ales de oscilación de
dientes a los tres autovalores, que resultan ser una cadena circu la r fo rm a d a por tre s masas e qui­
 distantes ¡guales enlazadas m ediante m uelles idén­
v = - k R 2 [ ( d l ~ e 2y + (d 2 - eyy + (03- e {y \
tic o s de constante recuperadora k.
S up o n e r que las masas sólo pueden m overse so­
bre las guías sin rozam iento.
A la vista de estos resultados las matrices T y V son

0 -1 -i \
i 1 °\ i 2
T - mR2 0 1 V = kR 2 -1 2 -1
0
:: -i- lo 0 1 / 1-1 -1 2/

Las frecuencias normales se obtienen de la ecuación

I 2 k - ma? -k -k \
det I -k 2k - ma? ~k
\ -k ~k 2 k - ma? /

= (2k - ma? ) 2 - 2k 3 -3 (2 * - ma?) k2 =0

Solución que nos lleva a

Como vamos a analizar sólo oscilaciones en la direc­ a) —0 y a? —6(0^0? + 9úJq = 0

ción angular las energías cinética y potencial son
donde a)n = \ l — . Entonces
T= \ m R 2(é\ + Q\ +

3k
Entonces a>2 ~ a)3 = y ^ - con 1ue e* Problema es
degenerado (dos autovalores iguales). Los modos nor­
males son los autovectores correspondientes a cada uno
de los autovalores

( Y - ^ í ) X ¿= 0
3

159
 Para oj, = 0 se tiene * Cualquier vector de la forma anterior es autovector y
por lo tanto satisface la ecuación de modos normales lo
2 k9 \ - kd\ - k d \ = 0 que es consecuencia de la degeneración ya qué todas las
oscilaciones no traslacionales tienen la misma frecuen­
-k e \ + 2k e \ - k d \ = o
cia. Lo m ejor es tom ar los autovectores ortogonales,
-k 0 \ - kQ \ + 2 k d \ =0 por ejemplo el par

Cuya solución es por ejemplo 8\ = 6] = Q\ = Q1 y en­

tonces

X 1 = A '^ n1

es un par ortogonal con la forma apropiada ya que es­

tán contenidos en el subespacio bidimensional X 23.
Evidentemente esta solución corresponde a un modo
de traslación puro sin vibración como ocurre siempre
que tenemos una frecuencia nula.
Para las otras frecuencias propias la única ligadura es

ex + 82 + 82 = 0 Un b loque de m adera de masa m se m ueve sin

rozam iento p o r un p la n o horizonta l. U nida al b lo ­
D e donde se deduce que los autovectores son de la que hay una barra de m asa despreciable que pivo-
forma ta en él, y en el o tro e xtre m o de la barra se coloca
una masa m, que a su vez está unida a un m uelle
fijó la la pared co m o se indica en la fig u ra . La altura
del p u nto de suspensión del m uelle coincide con la
lo n g itu d de la barra.
 k
V- ■k x \ + mgy2 = — (x - 16)2 - mgl 1 - —

Donde se ha hecho uso del desarrollo del potencial

en torno al punto de equilibrio, veámoslo con más de­
talle.
De la figura se desprende que la energía elástica conte­
nida en el muelle es

1. C alcular la posición de e q u ilib rio del sistem a. k[l2( 1 - eos 8)2 + ( x - l sen 0)2].
2. H allar las frecuencias norm ales de vibración.

Solución

1. Resulta evidente que el punto de equilibrio co­

rresponderá al punto en que el muelle no ejerza nin­
guna fuerza, esto es, que su longitud coincida con su
longitud natural (siempre que la altura a la que está sus­
pendido el muelle coincida con la longitud de la barra
del péndulo, como nos dice el enunciado).
2. Empecemos calculando las energías cinética y po­ Entonces si hacemos el desarrollo en torno a x = 8 = 0 y
tencial. Llamemos lQa la longitud natural del muelle, x a 84
la coordenada que describe el movimiento horizontal teniendo en cuenta que (1 - eos 8)2 = — se obtiene para
del bloque a partir de la posición de equilibrio y 8 al án­
gulo que subtiende el péndulo con la vertical. el potencial
Evidentemente tendremos
V ~ — k(x - W)2
x, = x X, = X

y, = o 0
Desarrollando el término de energía potencial y reali­
• x2 = x - 1 sen x2 = xX - 18 eos
X2 = zando la aproximación cuadrática en la energía cinética
.. y2 = - l eos 8 y2 = 18 sen 6 llegamos a que

Así pues
L - -i~~ítt(2x2 + l282 - 21x8) - (x2 + P82 - 2x18) + mgld2
m
T = — (x2j + y\ + x\ + y22) - — (2x2+ Pd2 - 21x8 eos 8)
Escribamos ya las matrices T y V

162 163
700 problemas de Mecánica Oscilaciones oeaueñas
 2m -m i k -kl Tomemos como coordenadas el alargam iento del
T = Y=
-mi m i2 / l-W kl2 + mgl muelle (y) y el ángulo que forma la barra con la vertical
(0). Entonces
S olución
Entonces la ecuación característica queda
x c m - y + L eos 0 x cm = y - L sen
mto2 \ mo? f 2mg \ mg _ =>
y CM= L sen 0 y CM= L d sen d
k ) k l kl ) kl H&
Las energías cinética y potencial son entonces
Las dos raíces positivas de esta ecuación bicuadrada
son las frecuencias normales que son
T = y m (x2CM+ y2CM) + -i IÓ2 =
k —2- 1 ± m 2g2 1
mg -|-----
m kl 2 ~k2F + T
= y m (y2 + L 2é2 - 2yÓL sen 0) + ~ ~ m{2L)262 =

1 4
= 2" m (y2 + — L 202 - 2ySL sen 0)
Una barra hom ogénea de m asa m y lo n g itu d 21
p ivota sin rozam iento a lreded or de uno de sus ex­
tre m o s, m an tenié ndose siem pre en el m is m o pla­ V = -m g(y + L cos 0) + — k(y - 1 )2
no XY. D icho e xtre m o está u n id o a un m uelle ideal
de co n sta n te re cu pe ra d ora k y lo n g itu d n a tu ra l /,
A partir de aquí podemos obtener fácilmente que los
que puede m overse sólo en la dirección vertical. El
términos cuadráticos en torno a la posición de equili­
o tro e xtre m o del m u e lle está u n id o al techo. Calcu­
brio son
lar las frecuencia s norm ales de o scilación del sis­
tem a.

v ~~2 (m8102 + ky2)

1 2
T = — m y2 + — mJJ-02

Es inmediato encontrar las frecuencias normales, ya

que la matriz V - aí2T es diagonal.
1L
4
0
J
0 k - m(i?

165
 1. Las ecuaciones que se presentan tienen la forma
af.2 = ------
k (üí2 = -7 7 - típica de las ecuaciones de Newton, de hecho es fácil ver
1 m 2 4L
que pueden obtenerse del potencial
S olución
Problem a 73
v (v, ’ v v*)= \ +
A l e stu d ia r a p ro x im a d a m e n te el m o v im ie n to de
un co n d e n sa d o de B ose-E instein fo rm a d o p o r un 1 1 1 p.
gas de á to m os neutros a T - 0 K se llega a las si­ 2v2 2v2 2v2 + vxvyvz
g uie nte s ecuaciones de evo lu ció n para las anchu­
ras del c o n ju n to de á tom os condensados en x, y, z 2. Cuando A; = 1 y vx = vy ~ vz = v el problema tiene
simetría esférica y queda por lo tanto unidimensional,
c/2 ,2
V + A V ——
1 p
+ -------------
siendo las ecuaciones
d r2 * x* vi v 2vyvz
d 2v 1 P
d2 12 1 P
V + A ¿V ------ 1---------------
d r 2 Vy + A vVv vi v2 yvxv7
Esta ecuación puede integrarse formalmente usando
d2 1 P la conservación de la energía
v. + AV. = — +
c/r2
' 7 z 7 V 32 v7V yV x

d on de P e s una constante positiva.

■E?= —211------1
7 (dvV + — 21v2+ ------+
\ drj
1 ___p
2
2v2 3v3 ’

1. D is c u tir si las ecuaciones de e vo lu c ió n son obteniéndose

ecuaciones de N e w to n o no, y, en su caso, encon­
tra r el potencial. dv'
2. Supóngase en prim er lu g a r que Ax = Ay = Az = 1 = i:
y v = v = vz = v. Reducir a cuadratura s el proble­ 2 E - —
(2v'2) (3 v'3)
ma y d is c u tir los tip o s de m o v im ie n to posibles.
3. S upóngase ahora que P > 0, Ax = k y = 1 y vx = Los puntos de equilibrio v0 del sistema satisfacen la
v . E scribir las ecuaciones que satisfacen los pun­ ecuación
tos de e q u ilib rio y e n contrar las frecuencias no rm a ­
les de las oscilaciones (trid im e n sio n a le s) en to rn o 1 P
a los m ism o s. C alcular el lím ite de estas fre cu e n ­
o Ko
cias cuando P —> +<*>.
que corresponde a una ecuación algebraica de quinto
grado

166
 Esta ecuación cuando P > 0 tiene una raíz real y dos donde se ha definido P. . = P l(4v'v¿z). Los autovalores
pares de raíces complejas como puede deducirse fácil­ de'V son equivalentes a ios cuadrados de las frecuencias
mente de su gráfica. Esto significa que existe un único normales (las «masas» de estos osciladores son uno
punto de equilibrio, que además es estable ya que el po­ como puede observarse en la ecuación de evolución),
tencial crece en cero y en infinito. Las raíces de esta obteniéndose el resultado
ecuación no son expresables por m edio de radicales
pero pueden obtenerse por métodos numéricos. O)
3. Los puntos de equilibrio ahora satisfacen las
ecuaciones 1/2

“ >2,3 = 2 y (1 + - P13) ± V (1 - A Y P13f - 8P32

P
vo - T X + viv
O K0z Por último, cuando P -» +~ el punto de equilibrio
puede calcularse de forma cerrada ya que las ecuacio­
nes son
A z V 0z 3 + v2 v 2
z Oz 0
donde v0 u = v0u, = v ur
0. P
Para obtener las frecuencias de las oscilaciones pe­ v„ =
viv
queñas en torno a estos puntos de equilibrio podemos 0 Oz

linealizar el potencial o equivalentemente las ecuacio­

nes de evolución haciendo vx(t) - vQ+ Sx(t), v (í) = vQ+ AX = - i
+ 5 (r), vz(t) = v0z + 8z(t), quedando las ecuaciones

P 28 + 8 y entonces la matriz de energía potencial queda

V 0V 0z
v, Oz
/ 3 1
38 P V=
A ± ^ . +i
viv
0 Oz 0
V , V
Oz

P 5 ' + 6 2L+ 25 í de donde se obtienen las frecuencias normales, que son

Í 1 S + a«a = - “ , , 2. . 2
_z

d.7 z z z vi V V,
Oz

cu, = V 2,
De aquí es posible deducir (usando las relaciones que
definen v0 y v0z) que la matriz V tiene la forma
2,3
[4 + 3A2 ± V T 6 l 9A4 - 16A2]1
V2
\
/ 1 + dv i + 2Pu 4.1 3,2
o [Adaptado de V. M. Pérez García, H. Michinel, J. I.
V= P 3,2 Cirac, M. Lewenstein, P. Zoller, Phys. Rev. Lett. 77 5320
4,1
yo (1996).]
P 3.2 1 + + 2P
3.2 ' 2.3
 C onsidérese ei sigu ie n te par de péndulos acopia­ Sumando y restando las ecuaciones de evolución para
dos. x e y se obtienen las de X e Y, que son

X + r X + gX il = F0 eos (ot
Y + r Y + (2k + g!l)Y = F^ eos wt

Usando las expresiones para las soluciones del oscila­

dor armónico forzado podemos escribir las de las ecua­
ciones anteriores siendo las frecuencias asociadas a X
F q COS í0 ! e Y, cox = V g íi ya)y = V 2 k + gil. Para x,y se obtiene el
resultado inmediatamente sin más que sumar o restar
según el caso

+Y |
Sobre el péndulo de la derecha actúa una fuerza
periódica F = F0 eos o>f en dirección horizontal. A de­
y = X -Y j
m ás, am b o s p é n d u lo s su fre n una fuerza de roza­ Cuando r 0 y en estado estacionario se encuentra
m ie nto p ro po rcion al a la velocidad F = - r x . La cons­
tante de rigidez del m uelle in te rm e d io es k. E scribir
y resolver las ecuaciones de m o v im ie n to y p articu­ X = —
Fo
! COS W t
2m
larizar para el estado estaciona rio cuando r — s» 0 y
no se cum ple ning un a condición de resonancia.
F
y = - 2 - eos wt
2m

S olución

Problem a 75
Las ecuaciones de movimiento de las oscilaciones pe­
queñas del sistema pueden escribirse directamente y son
C o n sid é re se la s ig u ie n te cadena de p é n d u lo s
m x - -m gx/l - r x ~ k(x - y) + F0 eos <ot 1 acoolados

my = -m gyll - r y + k{x - y) j

Donde x e y son los desplazamientos del primer y se­

gundo muelle respectivamente respecto a la posición de
equilibrio. Definamos las coordenadas

X * { x + y )l 2
Y = (x -y )i 2

170 171
 1. E n c o n tra r las ecuacion es de e v o lu c ió n efe! El lagrangiano queda entonces
sistem a si la a m p litu d del d e sp la za m ie n to es pe­
queño.
L= mPq, - X m gl(l - cos <p) - ± - kP £ (<P, - <P¡+t)
2. ¿Cuál es el lím ite co n tin u o de las ecuaciones ;=i ^ <=i ¿ ,•=i
discretas que se obtienen?
Aproximando 1 - cos (p. = <p2J2 y usando las ecuacio­
S ugerencia: Para pasar al lím ite co n tin u o consi­ nes de Lagrange se encuentran las ecuaciones de movi­
derar que la distancia entre los p éndu los (lo n g itu d miento, que son
natural de los m uelles) tie n d e a cero, y escrib ir la
masa de los p én d u los co m o m = p/0 y k - T /i0 con g k
y + — (p^ ~ (<Pm ~ 2<P, + *PM) = 0

2. Para pasar al límite continuo hacemos m = p/0,

S olución k = TU: y definimos (oQ= \ /gíl y <p(x = ilQ) = ípr Como

1. La energía potencial es la suma de las contribu­ 1 3^

ciones gravitatoria y elástica

la ecuación continua es
v^ = ^ m Sl^~ co s< p )
i=i
du<p - K (p + ‘t)l (p = 0
v,.„ = ' t ^ K h - Q 2
Í=1 ¿

y como el desplazamiento es a partir de la posición de donde v = \ —

VP
equilibrio
A esta ecuación se le llama ecuación de Klein-Gor-
don y es de gran im portancia en Física Matemática.
Este problema puede repetirse usando una cadena de
osciladores armónicos en vez de péndulos simples y en­
tonces se obtiene la ecuación de ondas. Se recomienda
hacerlo como ejercicio que puede verse detallado en
muchos libros, por ejemplo en H. J. Plain, The Physics
ofVibrations and Waves, Wiley (1993).

i = / sen <p. - 1 sen <pM = - <pM)

1

9
Sistem as de referencia no inerciaies

Problem a 76

Un sistem a de referencia acelerado S c o n base i, j,

k gira respecto a otro S, I, J, K que tiene el m ism o
origen y es inercial. La velocidad angular vale co = 2 f¡
- f2j + (2t + 4}k, siendo t el tie m p o , La posición de
una p a rtíc u la en el sis te m a m ó v il está dada p o r
r = (t2 + 1)i - 6fj + 4f3k. E ncontrar la velocidad y la
aceleración en los dos sistem as en el tie m p o t = 1.

S olución

E n este capítulo usaremos la siguiente notación. Los

sistemas acelerados serán designados por Sxyz y sus ve­
locidades y aceleraciones, por v = r y a = r. Los sistemas
inerciaies, por S¡X YZ y las velocidades y aceleraciones,
por v; = R y a; = R. Los vectores i, j, k serán los vectores
unitarios dirigidos según los ejes de S. Entre las veloci­
dades y las aceleraciones en dos sistemas se dan las si­
guientes relaciones

.a v¡ - v + co x r,
a, = a + 2e> x v + co x (ct> x r) + cb x r,
i

175
 La segunda ecuación expresa el llamado teorema de contiene a su extre m o C. En su in te rio r se encuen­
Coriolis [ver por ejemplo A. F. Rañada, Dinámica clási­ tra una partícula de masa m ; en el m o m e n to inicial,
ca, Alianza Editorial, Madrid, 1990, sección 10.1]. Los su velocida d respecto al tu b o es nula y se encuen­
términos segundo, tercero y cuarto corresponden, res­ tra a una distancia L0 de C. D espreciando el roza­
pectivamente, a la aceleración de Coriolis, a la centrífu­ m ie n to , c a lcu la r la ve lo c id a d con que sale p o r el
ga y a la azimutal. Derivando la expresión de r(í), se ob­ e xtrem o D del tu b o , el tie m p o que tarda en llegar
tiene la velocidad y la aceleración en el sistema a llí y la fuerza que ejerce el tu b o sobre la partícula,
acelerado com paránd ola con la de C oriolis.

v = r = 2t\ - 6j + 12r2k, a = f - 2i + 24ík,

y, evaluando para t - 1, S olución

v = 2i - 6j + 12k, a - 2i + 24k. Consideremos dos sistemas de referencia: uno acele­

rado Cxyz, con el eje Cx según el tubo y el Cz según la
\ Para calcular la velocidad y la aceleración en el siste- vertical y vectores unitarios i, j, k. La posición, la veloci­
v ma inercial necesitamos conocer los tres últimos suman­ dad y la aceleración de la partícula en ese sistema son
dos de la expresión del teorema de Coriolis. Haciendo x = xi, v = i i y a = xi. El vector velocidad angular vale
los productos vectoriales y sustituyendo el valor t = 1, se co = cok. El otro sistema de referencia C X Y Z es inercial;
llega a en el momento inicial, coinciden los ejes primero y ter­
cero de los dos sistemas.
v, = 34¡ - 2j + 2k, a; = 40i + 184j + 36k Se cumple el teorema de Coriolis

a; = a + © x ( c j x r ) + 2 o x v = r i - o^xi + 2aáj.

La ecuación de Newton se puede escribir en la forma

Un tu b o horizontal CD de lo n g itu d L gira u n ifo r­
m em en te con ve lo cid a d a n g u la r constante a> alre­
d e d o r de un eje ve rtica l, perpe n d icu la r al tu b o , que F - m gk = m [ ( x - <ü 2x ) \ + 2co¿j],

donde F es la fuerza perpendicular al tubo ejercida por

éste sobre la partícula (recuérdese que no hay roza­
miento). Proyectando sobre i, se tiene la ecuación en el
sistema del tubo

,x = 0)2X,

cuya solución general es

x (t) = A e 0’1+ BerM.

177
Las condiciones iniciales son x(0) = L0, x(0) = 0, de don­ tante w =V k /m a lre d e d o r del eje vertical que pasa
de A = B = L J 2, por lo que p o r O. ¿Qué ocu rre si ía frecuencia de g iro es d is ­
tinta?
x(r) = L0 coshotf.

Como v = x = L 0 senhwí, se cumple x 2 - v 2I oj2 - L 2, si­

guiéndose así que la velocidad respecto al tubo al llegar S olución
al extremo D vale
La ecuación del movimiento en un sistema inercial
con origen en O y cuyo plano X Y coincide con tt es
v= íoV L 2 - Ll.
m R + k R = 0,
Naturalmente la velocidad en el sistema inercial es
mayor, debido al giro; su valor es v - o>V2L2 - L 2, como
o sea que su aceleración es a; = -cü2R, con (a ='V k /m . La
se puede comprobar de la relación v; = v + co x r. Ei
aceleración en el sistema en rotación vale
tiempo T que tarda en llegar D es
a = a7- o x (co x r) - 2 o x v,
1 L 1 L+ V L ^L 2
T - — arcosh —— = — lo g ----------
(ú L0 o) L0 por lo que se cancelan la fuerza centrífuga y la del mue­
lle. Esto significa que la aceleración aparente en el siste­
En cuanto a la fuerza que ejerce el tubo sobre la par­ m a e n rotación coincide con la de Coriolis que es per­
tícula, una fuerza de ligadura, se deduce de la ecuación pendicular a la velocidad; la componente tangencial de
de Newton escrita más arriba que vale la aceleración es nula y la normal vale |a j = \a\ - v2!p.
Por lo tanto la trayectoria es una circunferencia de ra­
F - 2mcoxj + mgk, dio p = vl2co, recorrida con frecuencia 2co. Como a = -2 co x
x v = 2o)(yi - xj), la ecuación del movimiento en el sis­
o sea que es igual a la fuerza de Coriolis más la de la tema móvil es
gravedad (que vale - mgk), cambiadas de signo; o sea
que es la reacción a la suma de las dos. x = 2o>y, y = -2wx,

cuya solución general es\^

x = R sen(2cot + 8) + a, y = R cos(2fot + 8) + b,

Un p u n to m a te ria l de m asa m se m ueve sobre que es la ecuación de una circunferencia de radio R,

u n p la n o h o riz o n ta l lis o ir, conectada a un p u n to centrada en (a, b). Otro modo de obtener este resultado
fijo O m edia nte un m u e lle ideal de constante recu­ es resolver la ecuación en el sistema inercial. Como se
peradora k. In icia lm en te se encuentra en el punto trata allí de un oscilador armónico bidimensional, la so­
O c o n v e lo c id a d A \/k /m , s ie n d o A una lo n g itu d . lución general es la elipse
A n a liz a r su m o v im ie n to respecto a u n siste m a de
referencia en ro ta ción que g ira con ve lo cid a d cons­ X = R x sen(tot + <£ ), Y = R 2 cos(a¡t + <¿>2).

n o
 Podemos pasar al sistema giratorio mediante el cam­ Tomaremos los ejes del sistema acelerado S, con ori­
bio de coordenadas x - X cos cot + Y sen cot, y = - X sen gen en el punto O de la superficie terrestre, como indica
cot + r e o s cot, que nos conduce a la circunferencia obte­ la figura. El eje Oz es vertical, el Ox apunta al este y el S olución
nida más arriba con Oy al norte. Ese sistema gira en tornó al e je jle la Tierra
con velocidad angular^ = a>(cos Aj + sen~Xk)rsiendo A
2a - sen - R 2 sen <¿>2, 2b = - R ¡ cos 4>: + R 2 cos <f>2, la latitud del lugar. La ecuaciórTdélm ovim iento del
cuerpo en S es, según el teorema de Coriolis [A. F. Ra-
R, sen <b: + i?, sen <j>7 ñada, Dinámica clásica, Alianza Editorial, M adrid,
tan 5 = —é------ -r— zr------ 1990, capítulo 10]
Aj eos + R2 cos <p2
mr = mg - m(o) x (a) x R) + 2co x v),
R 2= ^ { K \ + R l + l R f i , cos (4>~
donde R es el vector con origen en el centro de la Tie­
rra y extremo en O, g es la aceleración de la gravedad y
Si la partícula está en t = 0 en el origen O y tiene velo­ se ha aproximado R + r por R pues r <§; R (r es del or­
cidad A V k/m, dirigida según el eje Ox, que suponemos den de los cien metros, R es el radio de la Tierra). La
coincide en ese instante con el OX, resulta a = 0, aceleración de Coriolis vale (2co(y sen A - z cos A)i -
b = -A¡2, R = A¡2, 6 = 0, o de modo equivalente R : = A , - 2 cox sen Aj + 2 cox cos Ak) y la fuerza centrífuga, (-co2 R
Ri = 0, <¿>t - <¿>2 = 0. cos A sen Aj + co2 R eos2 Ak).
Tomando componentes, la ecuación anterior se escri­
be pues

x = 2<o(y sen A - z cos A),'

D e term ina r la tra y e cto ria y la desviación hacia el y - -2 cox sen A - co2R cos A sen A,
este de un cu e rp o que cae sin velocida d inicial des­ z = -g + 2cox cos A.
de la altura de 200 m en un lugar de la titu d 40° N,
a p ro xim a d a m e n te la de Castellón de la Plana. La componente vertical de la aceleración centrífuga
co2R eos2 A se ha incluido en el valor local de la acelera­
ción de la gravedad g. (Esta componente vertical depen­
de de lá latitud y es responsable de las dos terceras par­
tes de la variación de g entre el ecuador y los polos; el
resto se debe a la forma no esférica de la Tierra).
Para resolver estas ecuaciones, harem os las aproxi­
maciones siguientes, justificadas por el pequeño valor
de lavelocidad angular co « 7,3 x 10"5 s: (i) en la ecua­
ción de z despreciamos el término de la aceleración de
Coriolis, con lo que z(t) = h ~ gt2!2; (ii) también despre­
ciamos las com ponentes horizontales de la velocidad
que son mucho menores que la vertical; por eso tom a­
mos como aceleración de Coriolis = - 2 coz cos Ai = 2cogí

181
eos Ai, que sólo contribuye a la ecuación de x. Esto con­ plano yz de un ángulo pequeño que, en esta latitud, vale
duce al siguiente sistema aproximadamente una décima de grado; ver A. F. R a­
nada, op. dt.). Esta es la razón de que usualmente se ha­
x ■- 2gtto eos A, y = - u)2R sen A eos A, ble sólo de desviación hacia el este.

cuya solución da la ecuación de la trayectoria en fun­

ción del tiempo Problema 80

Un tre n de 2.000 toneladas se desplaza a la v e lo ­

x(t) = ~ g ( o t3 eos A,
cidad de 15 m /s sobre raíles paralelos al m e ridiano
en sentido del sur al norte, en un lugar de latitud 40°.

y(t) - - (o2 R t2 sen A eos A, 1. D e te rm in a r la fuerza latera! que ejerce sobre

los raíles.
2. R epetir la m ism a cuestión si viaja en sentido
z{t) = h - ~ gt2.
del n o rte a) sur.

Eliminando t se encuentra la forma geométrica de la

trayectoria que es S olución

. a ¡o\ n - z y , „ , . n- z Usaremos el mismo sistema de ejes y las mismas no­

x = eo eos A \ ---- -------, y = - a>2 R sen A eos A---------
V 9g g taciones que en el problema anterior. La ecuación del
movimiento del tren es
Con los datos del problem a (h = 200 m, A = 40°,
R - 6.400 km), resulta que el móvil cae en un punto con F = ma = m(g - co x (co x R) - 2o) x v),
coordenadas x = 4,8 cm, y = -34 cm, z = 0 es decir hacia
el sureste. donde F incluye las fuerzas reales y las de inercia. Lo
Pero es conveniente hacer una aclaración. Al definir que nos piden es la componente x de esa fuerza, que es
los ejes, hemos tomado el tercero según la vertical geo­ la que tienen que soportar lateralmente los raíles. Está
gráfica, es decir, según la línea que pasa por el centro de claro que ni la gravedad g ni la fuerza de la máquina ni
la Tierra. Sin embargo, la dirección vertical se suele de­ la centrífuga -meo x (co x R) tienen proyección sobre x.
finir con plomadas, que son sensibles a la fuerza centrí­ Luego hay que calcular la componente x de la fuerza de
fuga pero no a la de Coriolis; por eso, en la práctica, el Coriolis, que vale
tercer eje se define por la dirección de la aceleración de
la gravedad, de modo que la desviación hacia el sur (en F = 2mojv sen A.
el hemisferio norte) desaparece, al estar incluida en la
redefinición de los ejes. (Nótese que los espacios reco­ Con los datos del problema, eso equivale a
rridos según los ejes y y z son proporcionales y que la
aceleración centrífuga contribuye a la ecuación de y con . 1. Senticfo de sur a norte: F = 2.815 N, es decir una
una constante eliminable mediante un giro de ejes en el fuerza hada el este sobre los raíles.
 2. Sentido de norte a sur: Fy = -2.815 N, es decir una
. fuerza hacia el oeste sobre los raíles. 10
Sólido rígido (!): M ovim iento piano
En los dos casos, la fuerza se dirige hacia la derecha
del tren.

D em ostrar que to d o piano de sim e tría es p erpen­

d icu la r a un eje p rin cip a l y que to d o eje de sim e tría
es principal.

S olución

- Todo plano de sim etría es perpendicular a un eje

principal. Llamando X Y al plano de simetría y o* a la
.densidad de masa, eso equivale a decir

, cr{x, y, z) ~ cr(x, y , ~¿)

y por tanto

Ixz = - j o xzd3r = 0 , / f =- J oyzdh - 0

por ser el integrando una función impar integrada sobre

un intervalo simétrico. El tensor de inercia es entonces

184 185
100 nm hlfím as ríp. M aránira
cos Ai, que sólo contribuye a la ecuación de x. Esto con­ plano yz de un ángulo pequeño que, en esta latitud, vale
duce al siguiente sistema aproximadamente una décima de grado; ver A. F. R a­
nada, op. cit.). Esta es la razón de que usualmente se ha­
x ■=2gtco cos A, y = - cú2R sen A cos A, ble sólo de desviación hacia el este.

cuya solución da la ecuación de la trayectoria en fun­

ción del tiempo P roblem a 80

Un tre n de 2.000 to n e la d a s se desplaza a la v e lo ­

*(0 = -jg a tf 3 eos A,
cidad de 15 m /s sobre raíles paralelos al m e ridiano
en sentido del sur al norte, en un lugar de latitud 40°.

y (0 = ~ ~2 0)2 sen ^ cos 1. D e te rm in a r la fuerza lateral que ejerce sobre

los raíles.

z{t) = h ~ ~ g t 2. 2. R epetir la m ism a cuestión si viaja en sentido

del norte al sur.

Eliminando t se encuentra la forma geométrica de la

trayectoria que es S olución

, l8 (h -z y 2 h -z Usaremos el mismo sistema de ejes y las mismas no­

x = oj cos A \ / ---- -------, y = ~(x>¿ R sen A cos A---------
V % § taciones que en el problema anterior. La ecuación del
movimiento del tren es
Con los datos del problem a (h = 200 m, A = 40°,
R = 6.400 km), resulta que el móvil cae en un punto con F = ma = m(g - o x ( o x R) - 2 o x v),
coordenadas x = 4,8 cm, y = -34 era, z = 0 es decir hacia
el sureste. donde F incluye las fuerzas reales y las de inercia. Lo
Pero es conveniente hacer una aclaración. Al definir que nos piden es la componente x de esa fuerza, que es
los ejes, hemos tomado el tercero según la vertical geo­ la que tienen que soportar lateralmente los raíles. Está
gráfica, es decir, según la línea que pasa por el centro de claro que ni la gravedad g ni la fuerza de la máquina ni
la Tierra. Sin embargo, la dirección vertical se suele de­ la centrífuga -m o > x ( o j x R) tienen proyección sobre x.
finir con plomadas, que son sensibles a la fuerza centrí­ Luego hay que calcular la componente x de la fuerza de
fuga pero no a la de Coriolis; por eso, en la práctica, el Coriolis, que vale
tercer eje se define por la dirección de la aceleración de
la gravedad, de modo que la desviación hacia el sur (en Fx = 2mojv sen A.
el hemisferio norte) desaparece, al estar incluida en la
redefinición de los ejes. (Nótese que los espacios reco­ Con los datos del problema, eso equivale a
rridos según los ejes y y z son proporcionales y que la
aceleración centrífuga contribuye a la ecuación de y con . 1. Sentido de sur a norte: F ~ 2.815 N, es decir una
una constante eliminable mediante un giro de ejes en el fuerza hacia el este sobre los raíles.

183
o;.
 2. Sentido de norte a sur: F = -2.815 N, es decir una
. fuerza hacia el oeste sobre los raíles. 10
Sólido rígido (I): M ovim iento plano
En los dos casos, la fuerza se dirige hacia la derecha
del tren.

D em ostrar que to d o plano ele sim e tría es p e rp e n ­

d ic u la r a un eje prin cip a l y que to d o eje de sim e tría
es principal.

S olución

. Todo plano de sim etría es perpendicular a un eje

principal. Llamando X Y al plano de simetría y a a la
-densidad de masa, eso equivale a decir

.. a(x, y, z ) = a(x, y, -z)

y por tanto

lxZ~ ~ J axz& r = 0 , íyz~ ~ \ v y zd 3r ~ 0

por ser el integrando una función impar integrada sobre

un intervalo simétrico. E l tensor de inercia es entonces

184 185
100 problemas de M ecánica
 IIxz Ixy o \ La misma demostración puede hacerse para / y en­
1yx 1 O tonces se prueba que el eje de simetría es principal.
yy

0 0 AJ
entonces el eje Z es principal.
Problema 82
Por otro lado si hay un eje de simetría de orden n la
distribución de masa es invariante bajo giros de ángulo Calcular el te n so r de inercia de un cubo de densi­
a = 2 7rln. En cilindricas en torno al eje dad p, para m asa M y lado b respecto a uno de sus
vértices.
cr(p, c£, z) = cr(p, (f>+ a ,z )

lxz~ J o-pzp eos <j)dz dp d<j>= | p2z dz dp j a eos 0 d(}>

S olución
Calculemos la última integral
Calculemos las integrales que definen el tensor de
í <r(p, <£, z) eos 4>d(f> = f crip, 4>, z)eos 4>dó + ... + inercia con respecto a los ejes naturales del cubo (los
J Jo que son paralelos a las aristas del mismo)
r 2ir
+ cr(p, <£, z)cos <f>d(f>=
llt-CL

- 1 cr(p, 0, z)[eos (f>+ cos(^)H- a) +... + cos(<j>+ (n- 1 ) a)]d(f) Ai = P £ dx^ í0¿ dx2 (x l + XD j* =\ Pb5 = y Mb2

Cb r b fb \ \
Calculemos la suma de cosenos usando que
712 = p j 0 ^ A J 0 * 2 * 2 J„ * 3 = - - ^ p b í = - — Mb1
eos ¡3 = + e~lP).
Por simetría es fácil ver que todos los elementos dia­
S = eos 4>+ cos((¡> + a) + ...+ cos(0 + (n - 1) a) = gonales del tensor de inercia deben ser iguales y lo mis­
mo ocurre con los no diagonales. Entonces si llamamos
a = Mb2 el tensor de inercia queda

/ 2 1 1 \
-w- a —-a -----a
3 4 4
= ~ Re[e'^(l + é a +... + e,Xn~^a) + étí0(1 + e~ia +... + 1 2 1
a —a -----a
4 3 4
Y entonces 1 1 2
a “ ~4a
\ ~ 4 T “/
p"}> j n-1 n-1 \ pi4> / i _ pina 1 _ p-¡na\

187
 C onsidérese un triá n g u lo e q u ilá te ro fo rm a d o p o r entonces
tres va rilla s delgadas de lo n g itu d L y masa m. Cal­
Problema 83 c u la r el m o m e n to de inercia a lre d e d o r de un eje
que pasa de sus vértices y es perpend icular al pla­ I = i m L 2 + m \2 -í- + L2\ = 4~ m L 2
no del triá n g u lo . 4 I 4 4 / 2

S olución P roblem a 84

El m omento de inercia es aditivo así que podemos U na barra u n id im e n sio n a l hom ogén ea de masa
calcular cada una de las contribuciones independiente­ m y lo n g itu d L se m ueve en un plano, estando uno
mente. Conocido el momento de inercia de una varilla de sus e xtre m o s fijo . D e te rm in a r su energía cin é ti­
1 ca y el m o m e n to a ngula r respecto del e xtrem o fijo.
son respecto a su centro de masas, que es Ic = — m L 2,
podemos aplicar el teorema de Steiner para cada una de
las varillas y resulta S o lución
/ = 3IC + m (d 2 + d\ + df) Comencemos calculando el número de grados de li­
bertad del problema. Al tener la barra un punto fijo la
complejidad del problema se reduce en tres grados de
libertad. Por otro lado la barra se mueve en un plano y
por lo tanto los grados de libertad quedan reducidos en
otros 2. Así pues tenemos cinco ligaduras para seis gra­
dos de libertad del caso general y por lo tanto un sólo
grado de libertad.
Con el fin de simplificar los cálculos y facilitar su in­
terpretación, vamos a estudiar los dos sistemas de refe­
rencia (el espacial y el ligado al sólido) de forma que in­
corporen la geometría y ligaduras de la barra:

1. Para ,el sistema de referencia espacial (E) tomare­

mos el origen O coincidiendo con el punto fijo de la ba­
donde d¡ son las distancias de los centros de los lados al rra. Tomaremos como plano O X Y el plano donde la ba­
vértice respecto al cual estamos calculando el momento rra está obligada a moverse y el eje O Z perpendicular al
de inercia que son plano del movimiento.
2. Para el sistema ligado al sólido (S) tomaremos el
Origen O' coincidiendo con el punto fijo de la barra y,
d^ -- d2
a ~ L
- 2 ' ¿3 . I 2 r
por tanto, con el origen del sistema Espacial. Tomaremos
el plano O 'X Y como el plano donde la barra está obli­
gada a moverse, tal que el eje O 'X coincide con el eje

188 . 189
100 problemas de M ecánica
que define la misma barra. El eje O'Z' es, al igual que en Sustituyendo la ecuación de la energía cinética e inte­
el caso anterior, perpendicular al plano del movimiento. grando en r de 0 a L se obtiene
En los cálculos posteriores vamos a utilizar tres dife­
rentes sistemas ligados al sólido cuya diferencia es el T ^ ^ r m L2ó¿í = ~ I á 2
origen: (1) S¡, cuyo origen es el punto fijo; (2) S2 cuyo 6 2
origen está en el centro de la barra (al ser ésta homogé­
nea, coincide con el centro de la masa); y (3) 53 cuyo que también podemos escribir como sigue
origen está en el extremo libre de la barra.

donde I = — m L 2 es el momento de inercia de la barra

respecto de un eje perpendicular a ella y que pasa por
uno de sus extremos. Siendo ó> = (0, 0, ¿) la velocidad
angular de la barra en el sistema de referencia de la ba­
rra y ¿n el inercial.
En este caso, la integración es inmediata debida a la
geometría sencilla del sólido. En general, resulta más
apropiado utilizar la expresión de la energía cinética re­
Una elección apropiada de la coordenada generaliza­ lacionada con el campo de velocidades del sólido rígido.
da es el ángulo a que determina la barra con el eje ver­ 2. Determinación utilizando el sistema S} ligado al
tical OX. sólido y el campo de velocidades.
Obsérvese que el eje instantáneo de rotación pasa En este caso la energía es:
por el punto fijo de la barra y es ortogonal al plano en el
que se mueve la barra.
Vamos a comenzar calculando la energía cinética. En
este caso podemos proceder de varias formas posibles:
ya que VQ/, la velocidad del origen del sistema ligado a
1. Determinación directa. A partir de la misma de­ al barra se anula por coincidir con el extremo fijo de la
finición de energía cinética, se tiene T = v2dm ya misma. Teniendo en cuenta que la densidad de masa es
£ cero en todo el espacio excepto en el intervalo [0, L] del
que se trata de una distribución continua de masa. La eje 0 'X \ tenemos la siguiente expresión del tensor de
contribución de un elemento de masa' (dm = crdr, siendo inercia:
cria densidad lineal) cuyos puntos tienen la misma velo­
cidad es la siguiente

r = (r cos a, r sen a, 0)

v = —- = rá(-sen a, cos a, 0)
di
 A partir del cual se obtiene, de nuevo, la expresión Teniendo en cuenta, como en el caso anterior, que la
anterior de la energía cinética. densidad de masa es cero en todo el espacio excepto en
3. Determinación utilizando el sistema 5, ligado ál un intervalo del eje O 'X ' que ahora es el siguiente
sólido y el campo de velocidades. En este sistema de re­ [-L/2, L/2], tenemos la siguiente expresión del tensor de
ferencia la energía es inercia:

0
/ 0 0 \
t =Y " V + y 0)1' * ' m
0 0
T 2m L l
ya que R = 0, al coincidir el centro de masa de la barra
con el origen del sistema ligado a la misma. El primer 1
0 0 mL2
término de la expresión anterior representa la energía \ 12 /
del centro de masa de la barra. Para su evaluación, te­
nemos que determinar la velocidad del centro de masa La velocidad angular, dado su carácter absoluto, es la
como sigue: misma que la considerada con el sistema de referencia
Sj y, en este caso, la expresión vectorial e también la
misma ya que los dos sistemas de referencia no inercia­
O O ' = eos a, sen a, 0 ies tienen los ejes O 'Z' paralelos.
El segundo término de la energía representa la ener­
¿ O O ' L gía de rotación alrededor de un eje que pasa por el cen­
V = --------- = — a(-sen a , eos a, U) tro de masa (origen del sistema de referencia no iner-
° dt 2
cial) y su expresión es la siguiente:

-5- m L 2á2
24

Sumándo las dos contribuciones, se obtiene de nuevo

la expresión de la energía cinética:

T —-~r m L 2á2
o

4. Determinación utilizando el sistema S3 ligado al

sólido y el campo de velocidades.
A partir de aquí, obtenemos la siguiente contribución En este sistema de referencia la expresión de la ener­
a la energía cinética: gía cinética consta de tres términos:

—- m L 2á2 r = y mV0,2 + mV0, • (o x O O ' + tal ■I ■o

 Y’
Dado que los tres sistemas-de referencia ligados a la
barra son paralelos, los correspondientes tensores de
inercia se relacionan entre sí m ediante una transla­
ción.
La velocidad angular, dado su carácter absoluto, es la
misma que la considerada con el sistema de referencia
Sj y, en este caso, la expresión vectorial es también la
misma ya que los dos sistemas de referencia no inercia-
les tienen los ejes O 'Z' paralelos.
Obtenemos así la contribución del tercer término a la
La contribución total de los tres términos ha de ser energía cinética:
igual a la obtenida con los sistemas de referencia ante­
riores:
~ m L 2á2
6
(a) De forma similar al caso anterior tenemos que
la contribución del primer término es la siguiente: (c) En los dos apartados anteriores hemos obtenido
los vectores necesarios para determinar la contribución
OO' = (L cos a, L sen a, 0) del segundo término de la energía cinética que es la si­
guiente:
V d° ° ' = La(~sen
v 0, = ——— 7 •( a, cos a, n\
0)
- m L 2á2
26
A partir de aquí, obtenemos la siguiente contribución
a la energía cinética: Obsérvese que esta contribución es negativa, lo cual
está relacionado con que el segundo término de la ex­
presión general de la energía cinética no es definido po­
~ m l? c t
2 sitivo. Es de destacar que ello no está en contradicción
con el hecho de que la energía cinética sea definida po­
(b) Teniendo en cuenta, como, en el caso anterior, sitiva que exige que la suma de todas las contribuciones
que la densidad de masa es cero en todo el espacio ex­ no sea negativa.
cepto en un intervalo del eje O 'X', que ahora es [-L, 0], Sumando las tres contribuciones, se obtiene de nuevo
llegamos a la siguiente expresión del tensor de inercia: la expresión de la energía cinética

T = i m L2á2
i 0 0 0 \ o

I= 0 -i- m L 2 0 El problema nos pide determinar también el momen­

to angular respecto del extremo fijo. De nuevo hay va­
rias formas posibles de calcularlo.
o o Y m L2>
 1. Determinación directa. A partir de la misma defi­ En este caso, tenemos que utilizar la expresión gene­
ral del momento angular
nición de momento angular, se tiene T = j ~ r x vdm
ya que se trata de una distribución de masa. De acuerdo L 0= m O O ' x VQ, + O O ' x (co x mR) + m R x V G+ í ■ co
con la expresión anterior de los vectores r y v, se obtie­
ne que el momento angular respecto del punto fijo de la Donde R = (-L/2, 0, 0) es el centro de masa de la ba­
barra es: rra en el sistema Sy Los otros vectores han sido calcula­
dos anteriormente. Es fácil verificar que se obtiene el
mismo valor del momento angular que en los tres pro­
L -- (0, 0, m L2á) - Ico cedimientos anteriores.

)C Al ser la velocidad angular paralela a un eje principal

de inercia de la barra el momento angular es paralelo a Problem a 85
la misma. Obsérvese que en este caso la integración es
muy sencilla debido a la forma del sólido aunque en el Calcular la aceleración de bajada p o r un p la n o In­
caso general resulta más útil la expresión que tiene en clinado un á n gulo a, de una esfera maciza h o m o ­
cuenta el campo de velocidades. génea de radio r y masa m.
2. Determinación utilizando el sistema S¡ ligado al
sólido y el campo de velocidades. En este caso, la expre­
sión del momento angular es S olución

L 0= í ■o
En primer lugar se entiende que la esfera baja por la
línea de máxima pendiente del plano inclinado, que co­
ya que el punto O coincide con el origen del sistema de rresponde a la inclinación definida por el ángulo a. El
referencia y además es un punto fijo del sólido. Basta planteamiento del problema es el mismo para una esfe­
sustituir la expresión del tensor de inercia asociado al ra hueca, así como para un cilindro hueco o macizo
sistema Sx y la velocidad angular para obtener el valor siempre que consideremos movimiento plano. La dife­
del momento angular. rencia es cuantitativa y se debe al diferente momento de
3. Determinación utilizando el sistema S2 ligado al
sólido y el campo de velocidades. En este caso el origen
del sistema ligado al sólido coincide con el centro de
masa del mismo, y la expresión del momento angular es

L 0= m O O ' x V 0, + í •o

D e acuerdo con las expresiones utilizadas para el cál­

culo de la energía cinética, se obtiene, nuevamente, la
expresión anterior del momento angular.
4. Determinación utilizando el sistema S3 ligado al
sólido y el campo de velocidades.

197
inercia con respecto al eje de rotación que pasa por su del plano R. De acuerdo con la ecuación de Newton y
centro de masa. del momento angular, se tiene:
Vamos a estudiar separadamente los casos de roda­
dura con o sin deslizamiento. d 2O O '
R + mg - m - ■
De acuerdo con la figura, podemos ver que el número
de grados de libertad es 2 ya que un punto arbitrario
queda determinado por el desplazamiento a lo largo del
plano inclinado y por una rotación alrededor de un eje -7 ~ L’ = O'A x R
dt 0
ortogonal a la línea de máxima pendiente del plano inli-
nado. En el caso del centro de masa (que coincide con donde R = ( T, N, 0), mg = mg(sen a, - eos a, 0), O O ' =
el centro geométrico al tratarse de una esfera homogé­ = (x, r, 0), O 'A = (0, -r, 0) y L'0 = (0, 0, 7<¿). Siendo
nea) se necesita solamente la coordenada de desplaza­ 2
miento según el plano inclinado. 1= ~ m r2 el momento de inercia de la esfera maciza con
Con el fin de simplificar los cálculos y facilitar su in­
respecto a un eje que pasa por su centro de masa. El
terpretación, se elegirán los dos sistemas de referencia
m om ento angular es paralelo a la velocidad angular
(el espacial y el ligado al sólido) de forma que incorpo­
ren la geometría y ligaduras de la esfera: © = (0 ,0 ,* ).
Las ecuaciones vectoriales anteriores son equivalen­
tes a las siguientes tres ecuaciones escalares con cuatro
1. E l sistema Espacial (E ). El Origen O está en el
incógnitas:
punto más alto del plano inclinado. E l plano inclinado
contiene el plano O X Z de m anera que el eje O X coinci­
d 2x
de con la línea de máxima pendiente del plano inclina­ m = T + m g sen <p
do. El eje O Z es perpendicular al plano del movimiento
que coincide con el plano inclinado.
2. El sistema ligado al sólido (S). El Origen O' es el 0 = N ~ mg eos <f>
centro de masa y geométrico de la esfera maciza. Los
ejes son paralelos a los del sistema espacial. Esto es po­ 7 d 24> _
/ ——- = rT
sible debido a la especial simetría de la esfera: Todo sis­ dt2
tema de ejes ortogonales centrados en el origen consti­
tuyen un sistema de ejes principales y, por tanto, el La cuarta ecuación que necesitamos para resolver el
tensor de inercia es el mismo en cualquiera de dichos sistema anterior se obtiene a partir de si la rodadura es
sistemas de ejes. En este contexto, si considerásemos con o sin deslizamiento:
unos ejes que se mueven solidariamente con la esfera
(no solamente con su centro de masa) el tensor de iner­ 1. M ovimiento de rodadura sin deslizamiento. En
cia es el mismo que en el caso de los ejes que se mueven general, un sólido rígido rueda sin deslizar sobre una su­
paralelos a los ejes del sistema inercial. perficie rígida arbitraria cuando los campos de velocida­
Una elección apropiada de las coordenadas generali­ des de ambos coinciden en los puntos de contacto. Se
zadas es el desplazamiento x del centro de masa y el án­ trata de una ligadura en las velocidades que en algunos
gulo <¡>de giro alrededor de un eje paralelo al eje OZ. casos puede ser integrada y transformarse en una liga­
Pasemos a escribir las ecuaciones de movimiento. So­ dura de las coordenadas, de m anera que se reduce el
bre la esfera actúan las fuerzas del peso mg y la reacción número de grados de libertad.
 En nuestro caso, el plano inclinado está en reposo y, trabajo de la fuerza de rozamiento es nulo y se puede
por tanto, la condición de rodar sin deslizar es simple­ aplicar entonces la formulación lagrangiana.
mente que el campo de velocidades de la esfera maciza 2. Movimiento de rodadura con deslizamiento. E n
se anule en el punto de contacto con el plano inclinado: este caso, la condición de ligadura es en la relación en­
V A = 0. De acuerdo con la expresión del campo de velo­ tre las dos componentes de la fuerza de reacción del
cidades, se tiene plano inclinado: \T\ = fx |iV|, siendo ¡x el coeficiente de
rozamiento por deslizamiento entre las dos superficies y
V A = 0 = V Q, + o x O 'A tal que la componente T se opone al movimiento de la
esfera sobre el plano inclinado, es decir
de donde se obtiene la ecuación escalar
'T ±
— \3C
TT\ tJ-m 8 cos a
x + r<¡) = 0'
Sustituyendo en el sistema de tres ecuaciones anterio­
integrando y eligiendo las condiciones iniciales tales que
la constante de integración sea cero: res, se obtiene las aceleraciones correspondientes, y se
puede estudiar el tipo de movimiento de la esfera en
x = -r<t> función de las condiciones iniciales. ]En este caso no se
puede aplicar la formulación lagrangiana ya que el tra­
bajo realizado por la fuerza de rozamiento no es nulo.
lo cual reduce el número de grados de libertad del siste­
ma a uno. Sustituyendo en las ecuaciones del movimien­
to, se obtiene para la aceleración angular y la de bajada,
la siguiente expresión: Problem a 86

5
x = — g sen a Investigar la diferencia entre el p e río d o de las os­
cilaciones de baja a m p litu d de un pé n d u lo s im p le
puntual y o tro cuya lenteja está co n stitu id a p o r una
éy - -’s sen a esfera pequeña de radio e.
Ir

Por otra parte, las componentes de la reacción del S olución

plano inclinado son
^ El momento de inercia de la esfera alrededor del eje
de suspensión que es un punto fijo puede calcularse uti­
T =- —
T 2 m g sen a
lizando el teorema de Steiner y resulta ser

N = mg cos a Ir0 = —
2 me1
2 + ma2

Observemos que el eje instantáneo de rotación, que

tiene la dirección del eje O Z (O ' Z pasa por el punto donde a es la distancia del punto de suspensión al cen­
instantáneo A , de contacto de la esfera con el plano. El tro del péndulo. La ecuación de movimiento será

201
 T

i
di \ d isco puede rodar sin deslizar sobre un plano hori-
=N=Fxa
dt j zonta!. C a lcu la r la fre c u e n c ia de las o scila cio n e s
> pequeñas e n -to rn o al p u n to de e q u ilib rio . El piano
j que contiene al disco es vertical.

S olución

La posición de equilibrio corresponderá a la situación

en que la masa m se encuentra en su punto más bajo.
Tomemos como coordenada la desviación con respecto
a la dirección de equilibrio 6. La condición de rodar sin
deslizar impone que x CM= -R 9y entonces

T 2 MR2f)2
1 cm = Í2- M v 2cm - J_
Como F x a = -mga sen 0u y —— = I nda la ecuación de
z dt z
movimiento en la aproximación de oscilaciones de pe­
queña amplitud (sen 6 - 9 ) queda

j
é + a2 9 = -g a 6

De aquí se tiene que el periodo es

Por otro lado

El incremento fraccional es
í Xm= XCM+ aSQnd
T ( 6 ) - T ( 0) 6» 1 y = a cos 9
7(0) 5a2
xm~ x cm + a2 cos2 + cos 6
ys m
2 = a2Ó2 sen2 6

La energía cinética de m es entonces

Un disco u n ifo rm e de masa M y radío R está en
el cam po g ra v ita to rio y tiene adherida una masa m
p un tua l a una distancia a del centro del m ism o. El Tm- ~ - m (R 262 - 2aR 62cos 6 + a2Ó2)

202 203
WO nrnhlp.mas d e M erá n ira .i . n i n jt
 Y la energía de rotación del disco es g u la r co. B es una rueda fija de radio b y centro O.
Una tercera rueda C de radio (a - b)¡2 con un exte­
rio r ru g o so toca a A in te rn a m e n te y a B exte rn a ­
m ente y está m o vid a p o r la fric c ió n de los puntos
T Rr - 4y l CrJM 2 = 4A M R2é2 de contacto. Si no hay deslizam ie nto, e n contrar la
velocida d a n g u la r de C(o>') y la de la línea que une
los centros de A y de C.
La energía cinética queda

T
1
= T C M + TR + T m = —
2
— M R2 + m R2 - 2 a R c o s 6 + a2
2

En cuanto a la energía potencial

V - -mga eos 8

El lagrangiano queda

— M R2 + m R2 -2 a R eos 9 +a2 + mga eos 8

- i 2

Haciendo la aproximación de oscilaciones pequeñas

S olución
e2
eos 8 - 1 - —, las ecuaciones de Lagrange quedan
Consideremos el movimiento de C que rueda sobre la
rueda fija B. E l centro instantáneo de rotación es I, que
— M R2 +m R2 -2 a R eos 6+a2 J8 +mga 9 = 0 es el punto de contacto de ambas ruedas. Si co' es la ve­
locidad de C y Q es el punto donde C toca a A la veloci­
que es la ecuación de movimiento de un oscilador armó­ dad de Q es IQ oí o (a - b)oj, perpendicular a IQ, pero
como no hay rozamiento la partícula de la rueda A en Q
nico de frecuencia
se mueve con la misma velocidad que la partícula de la
rueda C en Q, entonces
mga

M R2 + m(a - R )2 (a - b)co' = toa

aoi
(O =
a- b

La velocidad del centro se obtiene de manera análoga

Una rueda A, de fo rm a de a n illo de radio interno
a rota a lre d e d o r de su ce ntro O con velocida d an­ y resulta ser

205
 (a - b)(o' (ú justamente la mitad Ic = 2/5 M R 2. Aplicando dos veces
2 flT el teorem a de Steiner calculemos ahora el momento
de inercia respecto a O, para ello necesitamos conocer
La velocidad angular de la línea OP es entonces las coordenadas del centro de masas

1 1 aoj 1 CR , 1 cR M 5
acó
2 OP a+ b o pzdv = u l — « í * - (*-*)» )*= j *
— ttR 3

Problem a 89 y calculemos entonces l CM

H allar la frecuencia de las oscilaciones pequeñas Ic = ícM + M d(c - C M ) 2

de una sem iesfera u n ifo rm e de masa M que rueda
sin deslizar sobre una superficie horizontal.

z
Con esta información ya es posible calcular el mo­
mento de inercia respecto a O

l o = 1cm + M d( ° ~ CM)2

Por el teorema del coseno, cuando la semiesfera ha

girado un ángulo O

S olución

Para escribir la energía cinética lo mejor es calcularla

respecto al punto O de contacto entre la esfera y el pla­
no horizontal. Necesitamos el momento de inercia co­
rrespondiente al eje perpendicular al plano del movi­
miento y que pasa por O. El momento de inercia que es
conocido es el correspondiente a toda la esfera con res­
pecto a su centro, que es I = 4/5 M R 2. Como estamos
considerando una semiesfera el momento de inercia es

206
y de aquí Una rueda con fo rm a de disco u n ifo rm e de radio
r rueda con velocidad v hasta que encuentra un es­
83 9 3 calón de a ltura 1/4 r. ¿Cuál es el m e n o r v a lo r de v Problem a 90
Io = S ± M R z + ~ M R 2 + M R 2r ^ - M R 2 cos que le p e rm itirá re m ontarlo? ¿Y el m a yo r v a lo r de
v p a ra que después de pasado el escalón no pierda
contacto con el suelo?
11 3
= M R2 — t- eos
5 4

Al ser la velocidad de traslación O nula la energía ci­ S olución

nética es Supondremos que el disco rueda sin deslizar y que no
hay rebote.

T = y /(0 )é 2

y el lagrangiano

1 (7 3 \ . 3
L = — M R2[ y - — cos 0 j02+ y M gP cos 6
Cuando el cilindro llega al escalón se produce en pri­
mer lugar una colisión con el punto P y después un giro
Aplicando las ecuaciones de Lagrange se llega a entorno al punto P hasta que el centro de la rueda que­
da sobre el punto P. En el instante de la colisión el mo­
mento angular de la rueda se conserva. El momento an­
M R2\ ^ - - 4 “ eos 0 10 + Mg P sen 0 + Mi?2Ó2sen 0 = 0 gular en centro de masas es
1 5 4 / o o
1 v i
Jr = — M r2 — = — Mvr
Cuando las oscilaciones son pequeñas despreciamos c 2 r 2
los térm inos de órdenes superiores al prim ero en las
ecuaciones de evolución aproximando sen 0 = 0, cos 0 = 1 y entonces en torno al punto P
y se encuentra
i= —
3 rMv + —
1 Mvr = —
5 Mvr
JP = JC+ M \ioc x y0\
4 2 4

Después de la colisión, y mientras se va levantando,

la rueda gira en torno al punto fijo P, entonces
y de aquí la frecuencia de las oscilaciones pequeñas es
Jp = l p(úp
'15 g
W )tl 26 R donde, por el teorema de los ejes paralelos

209
 la rueda se moverá en una circunferencia de centro P.
El que la rueda no pierda contacto entonces puede ex­
Ip = ~ M r 2 + M r2 = í - M r \ presarse como que la reacción en P sea positiva (R > 0).
2 ¿
Como la ecuación de Newton en dirección radial es
Con estos datos podemos calcular cúp
Mg sen d - R = Mré2
5v
ü)D= R = Mg sen 9 - R = Mré2
6r

La energía cinética que tiene la rueda después del im­ Al crecer 6, sen 6 crece y 62 decrece de forma que si
pacto es R > 0 en el momento del impacto seguirá siéndolo des­
pués y entonces no habrá despegue. Como en el impac­
to sen B - 3/4 y 6 = 5 vl(6r) La condición es

Para que pase el escalón esta energía deberá ser m a­

yor que la diferencia de energías potenciales. O sea que

Uniendo este resultado al anterior se encuentra que

para que la rueda remonte el escalón y continúe rodan­
V l- T 5 gr
do sin perder contacto con el suelo su velocidad deberá
^Por otro lado, si la velocidad iniciál es demasiado estar comprendida en el rango
grande la esfera puede perder contacto con el suelo. E n
la figura se muestra la rueda en posición general cuando v í 2 ,— v n r -
—j — V g r < v < Vgr
/el radio a través de P forma un ángulo arbitrario $ con
/ la horizontal. Mientras haya contacto en P, el centro de

Un disco só lid o de masa m y ra d io r, se sujeta

p o r m e d io de un m u e lle de co n sta n te k y reposa
sobre un plano in clin a d o c o m o se indica en la fig u ­
ra. Si el c ilin d ro puede ro d a r sin deslizar, calcular
las ecuacion es de L a grang e del m o v im ie n to y la
frecuencia de las oscilaciones del cilin d ro .
 Donde el origen de potenciales gravitatorio se ha to­
mado en la posición más alta sobre el plano inclinado
que corresponde a 8 = 0. y /0 es la longitud natural del
resorte. Construyendo el lagrangiano L = T - V y apli­
cando las ecuaciones de Lagrange

— m r2Q^ = - kr28 + mgr sen a + krlQ

2k g 2 k l.
+ —— 8 = -2- sen a +
3m 3r 3m r
S olución
Definiendo una nueva variable 8 ~ d + mg sen al(2kr) +
Tomemos como coordenada generalizada el ángulo 0 + lQ/r la ecuación de movimiento queda
entre el punto de contacto y el punto de contacto que
correspondería a que el centro de masas del disco —que * 2k _
coincide con el punto en que se engancha el muelle— B+ - — 0 = 0
3m
estuviera en la pared. Escribiendo la energía cinética
como suma de la de traslación del centro de masas más
y entonces la frecuencia de las oscilaciones pequeñas en
la de rotación en tom o al mismo obtenemos
torno a la posición de equilibrio es

2k
*>= \ i -=—
3m

Un disco hueco de masa M y radio R puede girar

sin rozam iento a lre d e d o r de su eje, el cual está co­
locado en p o s ic ió n h o rizo n ta l. Un se g u n d o disco
m acizo de m asa m y radio r ( r < R), puede rodar sin
deslizar p o r el in te rio r del p rim e r disco. D eterm inar
las ecuaciones de m o vim ie n to .

Por otro lado la energía potencial es la suma de la

elástica más la gravitatoria que es S o lución

V = y k(r8 - /0)2 - mgrO sen a Definimos las coordenadas como los ángulos girados
por los cilindros con respecto a su posición en la parte

213
inferior. La energía cinética tendrá varias contribucio­ La forma de la energía cinética será
nes. En primer lugar estará la de rotación del cilindro
grande que es
^ T r a s l- C M = \ = \ m (R ~ r f ^ ^

Por último calculemos el ángulo girado por el cuerpo

en tomo a su centro de masas para poder calcular la velo­
cidad angular asociada. Démonos cuenta de que el giro en
tomo al centro de masas está producido tanto por el giro
sobre sí mismo como debido al giro del cilindro grande.
La primera afirmación es evidente. Para comprobar la se­
gunda pensemos en un giro del cilindro grande en tomo a
sí mismo de ángulo <f>mientras que el disco permanece sin
moverse. Entonces, el ángulo girado sobre su centro de
masas es precisamente <f) según se aprecia en la figura.

Por otro lado, para el disco tenemos dos contribucio­

nes:

1. La debida a la traslación de su centro de masas.

2. La debida a la rotación sobre el centro de masas.

Para la traslación del centro de masas consideremos

que el tramo recorrido por el mismo viene dado por

•íc = ( « - r ) ^ + í di5co.¡Men,o
Como los dos ángulos contribuyen cada uno en una direc­
Como la condición de rodar sin deslizar impone ción distinta la velocidad de giro en tomo al centro de masas
R a = rQ, entonces, al ser será íi>c —Ó-<j>, y por consiguiente la energía cinética es

^d isc o -in tern o = (R ~ r) a = (R - r) ~ d r E-cM = Y W c = Y m r2 < - é ~ ' t i 2

de donde Sumando todas las contribuciones

sc = ( R - r ) < } > + ( R ~ r ) ^ 6 T = ^ M R z(f¿ + ~ m r 2(d-<j>)2 + -^ - m ( R -r ) z í <£+ ^ 0

 í

La energía potencial es y la ecuación de Lagrange correspondiente es

V = mg(R - r) 1 - eos I — 6 + $ MR1 + — mr2 + m(R - r)2 <j>+ m -± r> i'9 +

Con esto podemos construir el lagrangiano L - T - V + mg{R - r) sen ^ + — 0j = 0

L ~ — M R + ~ m r2(0 - <f>)2 + - i m(R - r)21<j> + ~ 0 \ - Las dos ecuaciones obtenidas pueden reescribirse de
la forma siguiente

- mg(R - r) 1 - eos 1~ 0 + 4>

4>\MR2 + m (R2 - 2Rr + — r2)) + mrÓ R — —r + —
Á é 2 R

y las ecuaciones de m ovim iento se obtienen como = -~mg(R - r) sen I — 0 + <fi

siempre
mr^>1./?
■Y»_ —— _^ j + m r 2 ^ 1 - 2 _r + ^r¿
_5 rr h+——
dL . 1 r ( r '
w ^ M R 2é + ~ m r 2(é-4>) + m ( R ~ r f - U + - é
{R ~ r)r ( r \
= ~m g sen I 0 + <j>I
dL ( r \ r R R
3 5 -= - " * ( * - ' ) «en *

Calculando la ecuación de Lagrange para 0 se obtiene

j m r2 + m (R - r)2 Probar que la energía de una barra hom ogénea

m (R ~ r)2-~r~~r m r2
R 2 de masa m en m o v im ie n to plano es

í
+ mg{R - r) sen <f>+ ~ 0 j =0 T = ~ m ( u 2 + v 2 + uv)
o

donde u y v son las velocida des de los extrem os.

Por otro lado, para <f>se obtiene

S olución

La energía cinética de la barra es la suma de la ener­

gía cinética de rotación en torno al centro de masas más

217
la de traslación del mismo. Si los extremos son tA y rñ, al La energía cinética de rotación es entonces
ser homogénea la barra el centro de masas quedará si
tuado en el centro de la misma, rc = -i- (r, + rfl), entonces
Tr = — m i2 -ty (v2 + u2 - 2uv) = — m (v2 + u2 - 2uv)

1 1 1
T=— m y ( 4 + ¡a). =~2 m (ti+ v)2=-g- m (u2+v2+2uv) y añadiendo la parte correspondiente a la energía de
traslación del centro de masas
Por otra parte, la energía cinética de rotación en tor­
no al centro de masas es T - 4 - m (v2 + u2 - 2uv)
o

Otra forma de resolver el problema es darse cuenta

de que el campo de velocidades varía de manera lineal
donde / = ~ m i2 es el momento de inercia de la barra entre los dos extremos de la barra de la forma

respecto a su centro y <f>es el ángulo girado. La veloci­

dad de un extremo de la barra respecto al centro de ma­ v(x) - ~ [xu + ( L - x)v]
sas es

1 / x = y1( /v - u ) La energía cinética total de la barra será entonces

VB - V CM = V - y ( u + v)

Pero la velocidad lineal del extremo es la angular por T = ~ j** py (x fd x = Y ~ z i ^ J X + (L - X)Y? dx =

el radio

= ^ j j L [*2»2 + { L - x f v 1 + 2x(L - x)uv]dx =

que junto con la ecuación anterior nos proporciona la

igualdad
= J [x2(u2+ v2- 2uv) + 2Lx(uv - v2) + L2v2] dx-
m
y (T - “ ) = ¿ - i O® ~ O (u2 + v2- 2uv) + 2 (uv - v2) + U x 2
~2V

Elevando al cuadrado la expresión anterior

= ~ m ( u2 + v 2 + uv)
o
v2 + u2 - 2uv = <£2(rfl - r j = 4>2l 2
que es lo que queríamos demostrar.
4>2 = - i (v2 + u2 - 2uv)
 Una p a rtícu la de m asa m se m ueve con v e lo c i­
En resumen, tenemos tres ecuaciones para encontrar
dad v y colisio na con una barra delgada en reposo v', vb, co que son
de masa M y lo n g itu d /. La barra se encuentra s i­
tu a d a p e rp e n d ic u la rm e n te a la tra y e c to ria de la m v' + Mv'. = mv
p artícula, que la golpea a una distancia //4 de uno
de sus extrem os. D espreciando el rozam iento con mvl/4 = mv'l/4 + ¡zboj
el suelo y la rotació n de la barra en to rn o a su p ro ­ v' = v'b + lü)/4
pio eje, h a lla r la ve lo cid a d a ngula r de la barra tras
una co lisión pe rfectam ente plástica.
Usando que Izb = MP podemos resolver el sistema
cuyas soluciones son
S olución
Im M v
v —-------------
Apliquemos las leyes de conservación para obtener 4M + Im
los valores de la velocidad y el momento angular del sis­ , 4mv
tema después de la colisión. La conservación de la can­ Vb= 4M + lm
tidad de movimiento se escribe 12m v
(o ---------------
—

4m + 1M l .
m v' + M v ' = m v

El momento angular respecto al centro de masas an­

tes y después de la colisión será

j-m v li4

J c = mv'l/4 + Iz b. ce/uz

y está dirigido en dirección z por construcción. Igualan­

do se obtiene una relación entre las magnitudes

mvl/4 = mv'l/4 + 1 . 0/
____________________________ ZD

Por otro lado del hecho de que los cuerpos están uni­
dos tras la colisión se deduce que

v' = v; + W 4

220 221
1/1/1 ~ h/!---
 11
Sólido rígido (li): M ovim iento general

H allar la energía cinética de un co n o hom ogén eo

y macizo que rueda sin deslizar, según su genera­
triz, sobre un plano horizonta l y ta l que su vértice
está fijo .e n dich o plano. Su a ltura es h, el ra d io de
la b á se e s R y a e \ se m iá n g u lo del vértice.

S olución
A partir de ia figura podemos ver que, de acuerdo
con las ligaduras geométricas, el número de grados de

223
 libertad es 2 ya que un punto (el vértice) está fijo (-3) y del cono en contacto con.el plano donde se mueve. Ello
el cono está obligado a moverse sobre un plano (-1). quiere decir que dicha generatriz define el eje instantá­
Por otra parte, la condición de rodar sin deslizar implica neo de rotación y por tanto es paralela al vector veloci­
una ligadura en las velocidades generalizadas que en dad angular instantánea cu. Al estar u contenido en el
este caso es integrable según veremos. De todo ello se plano X Y una manera de expresar la propiedad anterior
Cocmc /)t ccCauu- m deriva que el sólido tiene 1 grado de libertad. es la siguiente
Con el fin de simplificar los cálculos y facilitar su in­
¿//a'-'?') & terpretación, se elegirán los dos sistemas de referencia k •o = 0
(el espacial y el ligado al sólido) de forma que incorpo­
ren la geometría y ligaduras del cono:
V '' V1 y sustituyendo la expresión de m:
1. E l sistema Espacial (E). El Origen O coincide con
-/ jjfw c el punto fijo del cono. El plano O X Y es el plano del
Ó + tj/k -i'-O
,y¡ r./tas&fyt-r movimiento. E l eje O Z es perpendicular al plano del
movimiento.
2. El sistema ligado al sólido (£). El Origen O' es el donde k • i' = sen a. Integrando y eligiendo nula la cons­
punto fijo del cono y por tanto coincide con el origen tante de integración, se obtiene la siguiente relación en­
y W « * & h U del sistema ¿spacial. El eje O 'X es el eje de simetría del tre los dos ángulos:
cono. El eje O 'Y' se mueve con el cono de manera que
está contenido en el plano O X Y. Esto es posible debido 8 + fjf sen a = 0
a la especial simetría del cono: Todo sistema de ejes or­
togonales tales que uno de ellos coincida con el eje de Una vez que hemos encontrado una ligadura entre
simetría del cono constituyen un sistema de ejes princi­ los dos grados de libertad del problema pasemos a cal­
pales y, por tanto, el tensor de inercia es el mismo en cular la energía cinética. Dada la elección del sistema li­
cualquiera de dichos sistemas de ejes. gado al sólido, tenem os que esta puede expresarse
Todo ello se indica en la figura. Una elección apro­ como sigue
piada de las coordenadas generalizadas son los ángulos
de rotación alrededor de los ejes OZ: 6 y O 'X : i/í.
Veamos cómo la condición de rodar sin deslizar in­ O)
troduce una relación entre las anteriores coordenadas
de forma que el sistema realmente tiene 1 grado de li­
bertad. La velocidad angular del cono puede escribirse donde la velocidad angular ha de expresarse en el siste­
en la forma ma de referencia del sólido. D e acuerdo con la figura te­
nemos que
o = 6 k + i/fi'
k'-S e-í( i -i coso l'<
k:sen a i ' + eos a k-
donde k e i' son vectores unitarios asociados a los siste­
mas de referencia espacial y del sólido respectivamente con lo cual, podem os expresar la velocidad angular
(k está dirigido según el eje Z e i' según el eje X ) . como sigue
El rodar sin deslizar supone que es nula la velocidad
de los puntos de la generatriz (definida por el vector u) (0 = 8 cos a (-c o t i' + k')

224 225
•' -I - /!» !- " "-
100 problemas de Mecánica
 El tensor de inercia del cono en el sistema elegido de
ejes principales es el siguiente:

0 m (R 2/4 + h2)

En este sistema de ejes resulta

o - i m(RV4 + h2)

J = -i- M R 2(a>, a>, 2a>)

Sustituyendo en la expresión de la energía cinética, se
obtiene:
Z'
rp 3
T = 4Óm e h l F 7 W

siendo tan a = R /h.

Un disco c irc u la r de radio R, m asa m y espesor

d e s p re c ia b le se m u e ve sin ro z a m ie n to a lre d e d o r
de un cen tro que perm anece fijo . E studiar el m o v i­ Al no actuar pares de fuerzas externas deberá satisfa­
m ie n to del disco sab iend o que en el instante inicial cerse en el sistema en movimiento que
gira con una ve lo cid a d o que fo rm a un á n g u lo de
TT
— con el eje del disco. ¿Cuánto vale el m ó d u lo del
— + <a x i = 0
m o m e n to angular? dt

Escribiendo la ecuación en componentes y explicitan-

S o lución do el producto vectorial

Tomamos los ejes móviles coincidentes con los ejes

principales del disco, como se índica en la figura
"z ~ ®
J

22 6 227
 y el módulo del momento angular es
La tercera ecuación nos dice que a> = cte = (o
2 °’

entonces
\ m = V MR2%

V2
°>x= — Y% < °y

Vi
% = — "o" El extre m o A de una barra cilin d rica hom ogénea
de masa m , lo n g itu d 21 y radio R desliza sobre una
recta v e rtic a l; el o tro e x tre m o B p uede m o ve rse
Derivando de nuevo la primera ecuación se llega a deslizando sobre un p la n o h o rizo n ta l. E scribir las
ecuaciones de m o v im ie n to m ediante el fo rm a lis m o
(¿>x + ~ - ( o x = 0 lagrangiano.

Z'
cuya solución es

(o {t) - A eos 1 + 4>

siendo A y <f>constantes que determinaremos después.

Para w (t) se encuentra

,. ^ ( a>0V 2 )
<a(t) = A sen —2-— t + <f>

Vi S olución
Imponiendo las condiciones iniciales a^(0) = — <w0,
Vi 2
toy(0) = 0 se encuentra que <j>~ 0 y A = — o)Qy entonces Escribamos la energía cinética respecto al centro de
masas

T = Trotación en torno a C M + Ttraslación dei C M

donde

^ ro ta c ió n en torno a CM 2 ^ ^

229
i
 T.traslación dei CM = — mví D e aquí se encuentra que las energías cinética y po­
tencial son
T
^ + sen2 0 ) + ~2 cos2 ^ +

+ m l24? sen2 6 + — m l262

2 2

V = mgl cos 9

Apliquemos las ecuaciones de Lagrange

1.
d dL ¿L_
D e la figura se encuentra que dt d<j> d(f>

o) = (9, <f>sen 9, <¡>eos 9) dL

como -=r-- - 0 la coordenada <f>es cíclica.
a<p
respecto a los ejes principales del cilindro x', y', z',
los que el tensor de inercia se escribe
■— [4>sen2 91j + (f>eos2 $I3 - m fél2 sen2 0] = 0
j ± Rz + ± l2 0 0 \
' 4 3 Esta ecuación refleja la conservación de Jz, al ser cí­
I =m clica la coordenada 4>su momento conjugado debe con­
o L servarse.
2.
0 0
(/j + m l2)$ = 4¿ sen 6 cos 0(/j - / 3 + m i2) - mgl sen 9
Por otro lado la posición del centro de masas es

x c := l sen 9 eos <f>

y c ---1 sen 9 sen <f>
H allar la energía cinética de un c ilin d ro de masa
zc - 1 eos 9
M, radio R y altura h, que gira con velocida d angu­
Derivando y elevando al cuadrado lar (o constante a lre d e d o r de un eje que fo rm a 45°
con el eje del c ilin d ro y que pasa p o r su centro de
masas. C alcular el par que debe aplicarse al so p o r­
v 2c = x2c + y2c + z 2c = l2$- + 12 sen2 dtft
te que sostiene el eje de rotación.

230
 En cuanto al par que debe aplicarse, las ecuaciones
S olución de Euler para el movimiento del sólido-son

Ixé i + ( h ~ 12)(0i a}2 = N ,

I2<*)2 + (Ii ~ I 3)(*>1a>3 ^ N 2
h ¿>3+ (I2~ I l)(Ol0>2 = N l
Los momentos de inercia son los momentos principa­
les que para el cilindro ya se vieron en un problema an­
terior, y son

Calculemos la energía cinética en ejes principales

I ^ ~ mR2
1

Además en ejes principales la velocidad angular es

Por la degeneración existente en el plano X Y pode­
mos elegir la dirección de co como la del eje Y con lo G>= (co sen 45, 0, co cos 45)
L V2 V2 \
que queda (o = 10, - y o>, col y evidentemente ó» - 0. Sustituyendo en las ecuaciones
de E uler se encuentra
Calculemos los momentos principales de inercia que
necesitamos
* ,= * 3 = 0, =
/ = f f f d3xr2p = p i 2" d<¡> \ m dz f * drr2r = - 1 M 2
Z J J J V' 0 -h/2 0 2
Este es el resultado en los ejes móviles; para calcular
/ . = 1 1 J fl°x(x2 + z2) = 1 1 1 d3x(rz eos2 <¿>+ z2) = el par en ejes fijos observemos que hemos puesto oj per­
pendicular al eje 2 y entonces N = 0 siempre. En cuanto
i*2ir j*W2 r/? j*2ir í*/i/2 j*/? a las otras componentes N = cos cot, Nx - N 2 sen <ot.
=p eos2 4>\ d z \ drr~ + p | z2d z drr -
JÜ " -h/2 J0 J0 J-hO Jn
I RA „ h2 R2\ ! ^ Í R 2 h2'
= p \ T T h ~ + l p — '2 _ J = M ( t + Í2

Un c o n o circ u la r de m asa M ra d io R y a ltu ra h

Ya podemos calcular la energía cinética tie n e un m o v im ie n to g iro scó p ico . Si su velocida d
de ro ta ció n es O. y su ve lo cid a d de precesión es <j>,
e n c o n tra r la energía cinética del cono si el p u nto
que perm anece fijo es su vértice.

232
 El sistema es el que se presenta en la figura R _ r
Solución X 0, Y0, Z 0: Son los ejes fijos. 1 7
X , Y, Z: Son los ejes móviles, ligados al sólido.
0: ángulo formado entre los ejes Z y Z0. R
r —— 7
4>: ángulo entre la proyección de X y X 0. h Z

entonces

f e 2 ir th fR z / h 3
I = p ) r2(Px = p d<j>\ d z \ r3dr = — M R2 ■
z Jv Jo o o 10

L = I (z2 + x2)dm =1 (z2 + r2 sen2 4>)dm =

y jv Jv
2 ir rh CRzIh f2 n fh fR z /h
ro Í) d<}>1 z2dz I
•'o Jo
rdr + p )
•'o
sen24>d<¡)\ d z \
Jo Jo
r ’d r -

= ^ A Í(iP + 4 ^ )

Con los momentos principales ya podemos calcular la

energía cinética
La energía cinética, expresada en ejes principales será

7' = t (/ X + ^ +/X)

Eligiendo los ejes apropiadamente

<ox = 0 Un astronauta se ve o b lig a d o a efectuar un ate­

(»y = (f>sen 0 rrizaje de em ergencia en un p u n to de 45° de la ti­
tu d , en un asteroide de 100 m de radio, perio d o de
ü)z = O, + 4>eos 0 ro ta ció n ig u a l a 1 hora y d ensida d cinco veces la
dei agua. A l cabo de un rato, el astronauta puede
Necesitamos calcular ahora los momentos.principales despegar, y lo hace en ve rtica l, pero abandonando
de inercia I z e Ix = Iy .
en el astro un pecio de diez tonelad as de m aterial.
C om o consecuencia el m o v im ie n to de rotación se
I z = Jí y r 2dm ~ p Jfy r2d?x v e a lte ra d o . C a lcu la r c u á n to ta rd a rá el a ste ro id e
en g ira r en to rn o a un eje que co rte a su a n tig u o
Como el cono es homogéneo m = pV = p7rr2h/3. Por ecuador.
otro lado la ecuación de la superficie del cono es

234
l f )f ) n m h l a m a c r í a
 donde el segundo término es la contribución del pecio,
cuya masa es m.
El asteroide queda como un trompo simétrico (es de­
cir, con dos ejes principales iguales A - B y otro distin­
to C), en rotación libre, como A = B - (2M/5 + m )R 2,
C = 2MR2!5. Según se deduce de las ecuaciones de Eu-
ler [ver A. F. Ranada, Dinámica clásica, Alianza Uni­
versidad, Madrid, 1990, sección 12.7], m precesiona alre­
dedor del eje 3, de modo que w3 y ¡o | son constantes,
con velocidad angular O dada por

_ C-A (o 5m
£1 = ü),--------- - —-^=r--------------- ,
3 A / V 2 - 2M + 5m
~CG' !'7
o sea que el eje de rotación traza un cono de 45° de se-
S olución
miángulo alrededor del eje 3, que corta a la superficie
en una circunferencia, tal como indica la figura. Al cabo
Antes del naufragio espacial, el astro tiene simetría
de un semiperiodo T = 7r/|ü|, el eje de rotación pasará
.esférica, de modo que su matriz de inercia es diagonal e
por el ecuador.
. igual a (en cualquier sistema de ejes con origen en su
Teniendo en cuenta los datos del problema, M = 2,09 x
centro)
x 1010kg, (o = 1,745 x 10~3 s-1, resulta ñ = -1,473 x 1(Y9
con lo que el eje de rotación pasará por el antiguo ecua­
0 0
1 1 dor al cabo del tiempo
í = 0 I 0
0 / T = 2,13 x 109 s * 68 años terrestres.
, 0

con I ~ 2MRV5 (el momento de inercia de una esfera

respecto a un diám etro cualquiera), siendo M y R la
masa y el radío del astro y suponiendo que la densidad
es uniforme.
' Tras el despegue, se rom pe la simetría esférica. Los
nuevos ejes principales son el que pasa por el centro y el
punto donde queda el pecio, que tom arem os como
eje 3, y otros dos perpendiculares cualquiera. E n esos
nuevos ejes, la matriz de inercia toma la forma

mR2 0 0
0 m R2 0
0 0 0

236
1 0 f í n r n h l a m a o H a A/l.
 /b e fa d o c/^¿(Cd/pCt'¿:^_cío/ts s d a y /a u ^ e . .

d C i 'q a c í o u . ^ f cknc/cuLeo> &+* ¡?c€?adu/<&)^ ( j c^

cj/eUes£ü? ¿r/¿ ¿’i'^ & C ÍU /'Z ^ ■ \^

¿5-y y ¿tV ¿ a 4 ( f c í ^ t- a i c r d & - 0 . \

á -K ct* . ^ ¿L tydif i-Áte dO ~0 $Gcf

rtC u°~ \
¿ fe fe rrn M c u * ^ ^ y , Q-4 q- *#y, d\ ^ d t :
¿Lz* s Ct-T- H ®'

’^d(f(jLS(jac's - d__ Pjs_ — 2JL - c c -v x ^ z ,

crt* ? * > . /

\i =■• * ^ Y ' H / jV c
¿7T ^

ri. í¿ - t í - ^ < ? '?w /a i° ^ ■

cW1 3(5 pff

■'x Z ’d j& . f 'z c o i ¿ y ¿ Ilg íü ^ L x o d ® . c

^ / f t o Lueu.h> ¿ fifí J t

g K Id - j C ^ * /k .
K ¿ * 7 t$ J 1

‘^ f d h u c L O ^ Ú S? OO L ^ d ú £L
¿a* 3 u2s *7 &o :
 ; Universidad de Zaraaoza Bibli<

3497006
i ’ !\ 3201788385

La colección 1 U \J PROBLEMAS ofrece a estudiantes y docen­

tes universitarios un repertorio de libros eminentemente prácticos que
cubren las materias básicas de las licenciaturas científicas y técnicas.
Cada uno de ellos constituye una valiosa herramienta para preparar
una asignatura concreta, puesto que sus autores, profesores con
amplia experiencia didáctica, han recogido en los ejercicios y sus
correspondientes soluciones todos los contenidos básicos que se le
plantean al alumno, desarrollados con claridad y concisión. Para
reforzar el enfoque decididamente funcional de estas obras, se ha
escogido un diseño de página con amplios márgenes que permiten
al estudiante realizar anotaciones.

En 100 PROBLEMAS DE MECÁNICA, los profesores VÍCTOR

M. PÉREZ GARCÍA, LUIS VÁZQUEZ MARTÍNEZ y ANTONIO
FERNÁNDEZ-RAÑADA presentan una selección de ejercicios
resueltos, especialmente adaptados al estudio de la Mecánica
según los planes de estudio actualmente en vigor. Asimismo, los
autores proporcionan una serie de recomendaciones generales,
muy útiles para abordar satisfactoriamente la realización de los
ejercicios.