Oscilaciones de un cuerpo que rueda sobre superficie cóncava.
Superficie cilíndrica cóncava
Como hemos demostrado, para que un cuerpo en forma de aro, cilindro o esfera de radio r, ruede sin deslizar en un plano horizontal o inclinado, se tiene que cumplir que vc=ωr, a fin de que la velocidad del punto P de contacto del cuerpo con el plano sea nula. Donde
- vc es la velocidad del centro de masas
- ω es la velocidad angular de rotación alrededor de un eje que pasa por el centro
En el caso de que el cuerpo ruede sobre una superficie cilíndrica, la velocidad del centro de masas es
La velocidad angular de rotación es
Para que la velocidad del punto P de contacto sea nula, se tiene que cumplir que
La relación entre el desplazamiento angular del centro del cuerpo θ y el ángulo φ girado por el cuerpo es (R-r)θ=rφ
Vamos a comprobar esta relación de forma geométrica
Un cuerpo de radio r rueda sin deslizar sobre una superficie cilíndrica de radio R. El centro del cuerpo describe un arco s=(R-r)θ (en color verde), mientras el cuerpo rueda una longitud u=Rθ sobre la pista (en color negro) e igual al arco de color rojo pintado sobre el cuerpo.
En la figura, se señala la posición inicial del punto P del cuerpo y su posición angular en el instante t, φ=u/r-θ=Rθ/r-θ o bien, (R-r)θ=rφ
Parte de la figura ha sido elaborada con el siguiente código
th=pi/3; %ángulo r=0.25; %radio del cuerpo que rueda hold on %pista fplot(@(t) cos(t), @(t) sin(t), [-pi/2,0], 'linewidth',1.5) %cuerpo ang=(1:360)*pi/180; x1=r*cos(ang); y1=-(1-r)+r*sin(ang); fill(x1,y1,[0.7 0.7 0.7]) line([0,0],[0,-1],'linestyle','--') line([0,1],[0,0],'linestyle','--') x1=(1-r)*sin(th)+r*cos(ang); y1=-(1-r)*cos(th)+r*sin(ang); fill(x1,y1,[0.7 0.7 0.7]) line([0,sin(th)],[0,-cos(th)],'linestyle','--','color','k') fplot(@(t) cos(t), @(t) sin(t), [-pi/2,-(pi/2-th)], 'linewidth', 1.5,'color','k') fplot(@(t) (1-r)*sin(th)+r*cos(t), @(t) -(1-r)*cos(th)+r*sin(t), [-th/r-(pi/2-th),-(pi/2-th)], 'linewidth',1.5, 'color','r') fplot(@(t) (1-r)*cos(t), @(t) (1-r)*sin(t), [-pi/2,-(pi/2-th)], 'linewidth',1,'color','g') plot(0,-1,'bo','linewidth',1,'markersize',4,'markeredgecolor','b', 'markerfacecolor','b') plot((1-r)*sin(th)+r*cos(-th/r-(pi/2-th)),-(1-r)*cos(th)+ r*sin(-th/r-(pi/2-th)),'bo','linewidth',1,'markersize',4,'markeredgecolor', 'b','markerfacecolor','b') hold off axis equal axis off
Ecuaciones del movimiento
Cuando el cuerpo asciende disminuye su velocidad vc y ω

Las fuerzas sobre el cuerpo:
- El peso mg
- La reacción N de la superficie cilíndrica en el punto de contacto P
- La fuerza de rozamiento Fr, cuyo momento hace que disminuya la velocidad angular ω
Formulamos las ecuaciones del movimiento:
Movimiento de traslación del centro de masas con aceleración ac=d2s/dt2
Movimiento de rotación alrededor de un eje que pasa por el centro de masas con aceleración angular α=d2φ/dt2
Para un cuerpo redondo (aro, cilindro o esfera) de masa m y de radio r y cuyo momento de inercia Ic=k·mr2 siendo k=1 para un cuerpo en forma de un aro, k=1/2 para un cilindro o disco y k=2/5 para una esfera.
La aceleración angular es negativa, disminuye la velociad angular y también la aceleración ac. Teniendo en cuenta la relación entre aceleraciones, ac=αr, sumamos las dos ecuaciones para eliminar Fr
Para una superficie cilíndrica de radio R, la relación entre el arco s que describe el c.m. y la posición angular θ del c.m. del cuerpo que rueda es, s=(R-r)θ. La ecuación diferencial del movimiento es
Ecuaciones de Lagrange
La energía cinética del cuerpo que rueda es la suma de la energía cinética de traslación del centro de masa, con velocidad vc y la energía cinética de rotación alrededor de un eje que pasa por el c.m. con velocidad angular dφ/dt
La energía potencial, tomando como nivel cero, el eje de la superficie cilíndica es, -mg(R-r)cosθ
La lagrangiana L=Ek-Ep es
La ecuación del movimiento se obtiene
El resultado es
Cuando las oscilaciones son de amplitud pequeña, hacemos la aproximación, sinθ≈θ
El cuerpo que rueda sobre la superficie cilíndrica, describe aproximadamente un Movimiento Armónico Simple de frecuencia angular ω o periodo P=2π/ω
Superficie cóncava cualesquiera

En la figura, se muestra la curva C’ sobre la que rueda el cuerpo redondo de radio r y la curva C que describe el centro de dicho cuerpo. La variable u mide el desplazamiento OP a lo largo de la curva C’. Mientras el c.m. se desplaza s a lo largo de la curva C, el cuerpo gira un ángulo φ.
El ángulo θ mide la pendiente de la recta tangente a la curva C’ en P
Movimiento de rotación y traslación del c.m.

En la figura, vemos que en un intervalo de tiempo comprendido entre los instantes t y t+dt,
-
el punto de contacto P se ha movido la distancia PP1= du a lo largo de la curva C’,
-
el c.m. se ha movido la distancia QQ1=ds a lo largo de la curva C.
Si trazamos las tangentes a la curva C’ en los puntos P y P1 y las normales a dichas direcciones, vemos que se cortan en un punto O denominado centro de curvatura, la distancia entre dicho punto y P o P1 se denomina radio de curvatura ρ.
Movimiento de traslación del cuerpo.
El ángulo que forman las rectas tangentes dθ es igual al ángulo que forman las direcciones normales.
El arco PP1 es igual al producto del radio OP por el ángulo comprendido dθ.
du=ρ· dθ
El arco QQ1 es igual al producto del radio OQ por el ángulo comprendido dθ.
ds=(ρ-r)· dθ
Movimiento de rotación del cuerpo.
Ahora, nos fijamos en el movimiento de rotación en el intervalo de tiempo comprendido entre los instantes t y t+dt
Como el cuerpo rueda sin deslizar, la longitud del arco PP1=P1P’=du. La pendiente de la curva se ha incrementado en dθ, el ángulo girado por el cuerpo en dicho intervalo dt de tiempo es
Eliminando dθ, tenemos las relaciones
Teniendo en cuenta esta última relación dφ=ds/r, la ecuación del movimiento del cuerpo que rueda se obtiene eliminando Fr, es decir, sumando la ecuación del movimiento de traslación del c.m. y la de rotación alrededor de un eje que pasa por el c.m.
La trayectoria del centro de masas
Vamos a demostrar que cuando la trayectoria del centro de masas es una cicloide, el cuerpo que rueda sin deslizar describe un MAS.
El c.m. describe un Movimiento Armónico Simple (MAS) si la aceleración d2s/dt2 es proporcional al desplazamiento s y de signo contrario a éste. Si escribimos
donde R es una constante con dimensiones de longitud que no es el radio r del cuerpo como veremos más adelante. El c.m. describe un MAS
de frecuencia angular y periodo
Para una esfera k=2/5
Ecuación de la trayectoria C del centro de masas
s=4R·sinθ
ds=4R·cosθ·dθ
dx=ds·cosθ=4Rcos2θ·dθ=2R(1+cos2θ)·dθ
dy=ds·sinθ=4R·sinθ·cosθ·dθ=2R·sin2θ·dθ
Integrando x e y con las condiciones iniciales θ=0, x=0, y=r.
x=R(2θ+sin2θ)
y=r+ R(1-cos2θ)
Estas son las ecuaciones paramétricas de una cicloide
Ecuación de la curva C' sobre la que que rueda el cuerpo
ds=4R·cosθ·dθ. De la relación entre ds y du y como du=ρ·dθ
Obtenemos el valor del radio de curvatura ρ
ρ=r+4R·cosθ
La relación entre el desplazamiento infinitesimal du del punto P y la pendiente θ es
du=(r+4R·cosθ)·dθ
dx’=du·cosθ=(r·cosθ+4R·cos2θ)·dθ=(r·cosθ+2R·(1+cos2θ))·dθ
dy’=du·sinθ=(r·sinθ+4R· sinθ·cosθ)·dθ=(r·sinθ+2R· sin2θ)·dθ
Integrando x’ e y’ con las condiciones iniciales θ=0, x=0, y=0.
x’=r·sinθ+R(2θ+sin2θ)
y’=r(1-cosθ)+
R(1-cos2θ)
Esta no es la ecuación de una cicloide salvo que r→0
Representamos la pista y la trayectoria del centro de masas, una cicloide, utilizando código MATLAB
R=1; r=0.25; %radio del cuerpo que rueda hold on fplot(@(t) r*sin(t)+R*(2*t+sin(2*t)), @(t) r*(1-cos(t))+R*(1-cos(2*t)), [-pi/2,pi/2], 'linewidth',1.5, 'color','k') %pista fplot(@(t) R*(2*t+sin(2*t)), @(t) r+R*(1-cos(2*t)),[-pi/2,pi/2]) %trayectoria c.m. hold off axis off axis equal
Comprobación
En la figura, vemos que al punto P(x’, y’) es un punto de la curva C’ le corresponde el punto Q(x, y) de la curva C. La relación entre x’ y x, y entre y e y’ es
x=x’-r·sinθ
y=y’+ r·cosθ
La curva C que describe el c.m. es una cicloide invertida (-π/2<θ< π/2) cuyo mínimo está en (0, r) y cuyo disco generador tiene un radio R. Sin embargo, la curva C’ sobre la que rueda el cuerpo no es una cicloide salvo que el radio del cuerpo sea r=0.
Movimiento Armónico Simple
Cuando un cuerpo describe un MAS su desplazamiento (en este caso el arco s) se expresa en función del tiempo de acuerdo con la ecuación.
s=A·sin(ω t+φ )
vc=ω s0·cos(ω t+φ)
Donde la amplitud A y la fase inicial φ se determina a partir de las condiciones iniciales, en nuestro caso t=0, s=s0, vc=0. La rueda parte de la posición s0 con velocidad inicial cero.
s=s0·sin(ω t+π/2)=s0·cos(ω ·t)
La velocidad vc del c.m. de la rueda se obtendrá derivando s respecto del tiempo
vc=-ω s0·sin(ω t)
Balance energético
El balance energético es similar al que efectuamos al estudiar el movimiento de un cuerpo que baja rodando por un plano inclinado.
Cuando la rueda se encuentra en la posición dada por el arco s a lo largo de la cicloide, o a una altura y sobre el origen. La energía potencial es
La energía cinética es la suma de la energía cinética de traslación del c.m. y de rotación alrededor de un eje que pasa por c.m.
donde se ha tenido en cuenta la relación entre ambas velocidades vc=ω ·r para que ruede sin deslizar.
La suma de la energía cinética Ek y potencial Ep es constante e igual a
Actividades
Se introduce
- La posición inicial s0 (longitud inicial del arco de la cicloide), en el control titulado Arco inicial
- El cuerpo que rueda, en el control titulado Cuerpo.
- El radio del disco que genera la cicloide se ha fijado en R=1 m, con ello queda fijada la forma de la superficie cóncava
- El radio del cuerpo redondo se ha fijado en r=0.4 m
Se pulsa el botón titulado Nuevo.
Observamos como el cuerpo rueda sin deslizar por un camino en forma de cicloide. La ecuación del movimiento es
s=s0·cos(ω ·t)
Donde s es la longitud del arco de cicloide y ω la frecuencia angular del MAS
En la parte superior, se muestra el balance energético en un diagrama en forma de tarta. La energía potencial Ep y la energía cinética Ek como suma de la energía cinética de traslación del c.m. y de rotación alrededor de un eje que pasa por el c.m.
Braquistócrona en el campo gravitatorio terreste
En la página titulda Braquistócrona en el campo gravitatorio terreste, tratamos de encontrar la curva que hace que una partícula que cae bajo la acción de la gravedad, emplee un tiempo mínimo en recorrerla. Utilizaremos el cálculo de variaciones para obtener la ecuación de la curva que como vimos es una cicloide
En este apartado, se estudia un ejemplo similar. El centro de masas de una esfera que rueda sin deslizar bajo la acción de su peso, describe una cicloide.
Trayectoria del centro de la esfera
En la figura vemos una esfera de masa m y radio r que rueda sin deslizar a lo largo una pista de color azul paralela a la trayectoria del centro de la esfera (en color rojo). Vamos a demostrar que esta trayectoria del c.m. es una cicloide, pero la pista no lo es.

El centro de la esfera parte del reposo cuando coincide con el origen. Cuando el centro ha descendido una altura y, la velocidad de traslación del c.m. es v y la velocidad de rotación alrededor del un eje que pasa por el c.m. es v/r. La conservación de la energía se escribe
El tiempo que tarda el centro de la esfera en desplazarse de A a B es
Aplicando la ecuación de Euler-Lagrange a la función f que no depende de x
se obtiene las ecuaciones paramétricas de la curva denominada cicloide
La pista sobre la que rueda la esfera

En la figura, vemos que en un intervalo de tiempo comprendido entre los instantes t y t+dt,
-
el punto de contacto P se ha movido la distancia PP1= du a lo largo de la pista C'
-
el centro de la esfera se ha movido la distancia QQ1=ds a lo largo de la curva C.
Si trazamos las tangentes a la curva C’ en los puntos P y P1 y las normales a dichas direcciones, vemos que se cortan en un punto O denominado centro de curvatura, la distancia entre dicho punto y P o P1 se denomina radio de curvatura ρ.
El ángulo que forman las rectas tangentes dθ es igual al ángulo que forman las direcciones normales.
El arco PP1 es igual al producto del radio OP por el ángulo comprendido dθ.
du=ρ·dθ
El arco QQ1 es igual al producto del radio OQ por el ángulo comprendido dθ.
ds=(ρ-r)·dθ
La longitud del arco ds entre Q y Q1 es

La longitud del arco du entre P y P1 es
Sus componentes son
Integramos respecto del parámetro φ
Se determinan las constantes c1 y c2 sabiendo que para φ=0, x'=-r e y'=0
Las ecuaciones paramétricas de la pista sobre la que rueda la esfera son
Vamos a determinar el valor del parámetro R para que el centro de la esfera pase por el punto (x1, y1). Para ello, hay que calcular la raíz φ1 de la ecuación trascendente
Conocidos y1 y φ1 calculamos R,
Calculamos el parámetro R de la trayectoria del centro de la esfera, sabiendo que pasa por el punto (x1=2, y1=1). Representamos la trayectoria del c.m. y la pista sobre la que rueda sin deslizar una esfera de radio r=0.25
x1=2; %punto y1=1; r=0.25; %radio hold on f=@(x) (1-cos(x))/(x-sin(x))-y1/x1; phi=fzero(f,pi); R=y1/(1-cos(phi)); %parámetro R ang=(1:360)*pi/180; %esfera fill(x1+r*cos(ang),y1+r*sin(ang),'c') %la Tierra %cicloide fplot(@(th) R*(th-sin(th)),@(th) R*(1-cos(th)),[0,2*pi]) %pista fplot(@(th) R*(th-sin(th))-r*cos(th/2),@(th) R-R*cos(th)+r*sin(th/2),[0,2*pi]) plot(x1,y1,'o','markersize',3,'markeredgecolor','r','markerfacecolor','r') set (gca,'Ydir','reverse') hold off axis equal grid on xlabel('x') ylabel('y') title('Esfera que rueda')
Tiempo que tarda la esfera en discribir un arco de cicloide
Supongamos que la esfera parte del reposo desde el punto A (x0, y0), vamos a calcular el tiempo que tarda en llegar al punto B
El principio de conservación de la energía se escribe
Calculamos el tiempo T que tarda en desplazarse entre dos puntos A y B
Haciendo el cambio z=cos(φ/2)/cos(φ0/2)
Si A es el origen φ0=0 y si B es el punto (x1, y1), calculamos la raíz φ1 de la ecuación trascendente y después, el parámetro R, tal como hemos hecho anteriormente
El tiempo T que tarda en llegar la partícula al punto B describiendo un arco de cicloide, partiendo del origen en reposo es
x1=2; %punto y1=1; r=0.25; %radio f=@(x) (1-cos(x))/(x-sin(x))-y1/x1; phi=fzero(f,pi); R=y1/(1-cos(phi)); %parámetro R T=sqrt(7*R/(5*9.8))*phi; %tiempo disp(T);
0.9536
El periodo de la oscilación es 4 veces el tiempo que tarda en ir de A, φ0=0 a B, φ1=π (punto más bajo de la trayectoria)
Que es el resultado obtenido en el primer apartado
Referencias
Gillespie D. T. Simple harmonic motion of a round body rolling on a concave curve. Am. J. Phys. 52 (2) February 1984, pp. 180-182.
Eric Rodgers. Brachistochrone and Tautochrone Curves for Rolling Bodies. Am. J. Phys. 1946, pp. 249-252