Oscilaciones de una partícula situada en el borde de un aro que rueda
La cicloide se produce cuando se hace rodar sin deslizar un disco sobre una superficie horizontal. Un punto del borde del disco describe una curva que se denomina cicloide (palabra griega que significa circular). A un giro del disco le corresponde un arco de la cicloide
Si el punto P en el instante inicial está en la parte superior del disco, al cabo de un cierto tiempo t las coordenadas del punto P serán, tal como se muestra en la figura

x=vc·t+R·sinθ
y=R-R·cosθ
donde R es el radio del círculo y θ el ángulo girado en el tiempo t, θ=ω·t.
La relación entre la velocidad de traslación del centro de masas vc y de rotación ω alrededor de un eje que pasa por el c.m. es vc=ω·R.
La ecuaciones paramétricas de la cicloide son
x=R(θ +sinθ)
y=R(1-cosθ)
Oscilaciones de una partícula situada en el borde de un aro que rueda
Supongamos que una partícula de masa m se sitúa en el bode de un aro de radio R, que supondremos de masa despreciable, y que rueda sin deslizar sobre una superficie horizontal. En la figura, se muestra la situación inicial:
- la partícula forma un ángulo θ0 con la dirección vertical
- el centro del aro está a una distancia Rθ0 del origen O, posición de equilibrio estable.
Trataremos de determinar la posición de la masa puntual a cabo de un cierto tiempo t, así como el periodo P de las oscilaciones que describe este sistema mecánico.
La condición de rodar sin deslizar, implica que el aro gira un ángulo θ, a la vez que su centro se traslada una distancia Rθ, tal como se muestra en la figura.
Las coordenadas de la partícula respecto al sistema de referencia, situado en la posición de equilibrio estable son
x=Rθ-Rsinθ
y=R-Rcosθ
Las componentes de la velocidad partícula se obtienen derivando, respecto del tiempo t, las coordenadas x e y.
Principio de conservación de la energía
La energía cinética es
la energía potencial es
La energía total, es la suma de ambas contribuciones, e igual a la energía potencial cuando la partícula se encuentra en el extremo de su trayectoria, θ=θ0, dθ/dt=0.
Separando las variables t y θ, la ecuación del movimiento se escribe
El signo menos aparece porque el ángulo θ disminuye con el tiempo. En el instante inicial t=0, la partícula se encuentra en la posición de máximo desplazamiento, θ=θ0 y se dirige hacia el origen O.
Para integrar la ecuación diferencial, transformamos el denominador en otra forma equivalente mediante la relación sin2θ+cos2θ=1
Haciendo el cambio de variable
donde u0=1 es el valor de u para θ=θ0, que es la posición inicial de partida tal como se muestra en la primera figura. Ω es la frecuencia angular de las oscilaciones de un péndulo de longitud igual al radio R de aro.
Esta integral es inmediata y puede consultarse en la tabla de integrales
Como 0≤θ/2≤θ0/2≤π/2, se cumple que u≥1
Despejamos la raíz cuadrada, elevamos ambos miembros al cuadrado, simplificamos y despejamos u.
Finalmente, deshacemos el cambio de variable, para obtener la ecuación del movimiento.
cos(θ/2)=cos(θ0/2)cosh(Ωt/2)
Periodo de las oscilaciones
El tiempo que tarda en alcanzar la posición θ=0, mide un cuarto del periodo P de este oscilador mecánico
1=cos(θ0/2)cosh(ΩP/8),
De esta ecuación, despejamos el periodo P. Para ello, utilizamos la relación cosh2θ-sinh2θ=1
El periodo P es función del máximo desplzamiento angular, θ0, tiende a cero cuando la amplitud se hace pequeña, a infinito cuando θ0 tiende a 180º (posición de equilibrio inestable). En el péndulo simple, el periodo es función de la amplitud, pero es aproximadamente independiente de la amplitud cuando se separa ligeramente de la posición de equilibrio estable.
R=1; %radio que genera la cilcoide P=@(x) (8*sqrt(R/9.8))*log((1+sin(x/2))./cos(x/2)); fplot(P, [30,170]*pi/180) set(gca,'XTick',0:pi/6:pi) set(gca,'XTickLabel',{'0','\pi/6','\pi/3','\pi/2','2\pi/3','5\pi/6','\pi'}) grid on xlabel('\theta_0') ylabel('P(s)') title('Periodo en función de la amplitud')
Ecuaciones del movimiento

Como el aro no tiene masa, el centro de masa coincide con la posición de la partícula. La posición, velocidad y aceleración del centro de masas son
Las fuerzas sobre el sistema formado por el aro y la partícula son:
- El peso, mg
- La reacción del plano horizontal, N
- La fuerza de rozamiento, Fr necesaria para que el aro ruede sin deslizar
Las ecuaciones del movimiento son
- Movimiento en la dirección horizontal
- Movimiento en la dirección vertical
- Rotación alrededor de un eje perpendicular al plano del aro y que pasa por el centro de masas
El momento de inercia Ic=0, ya que la posición de la partícula coincide con el centro de masas
Esta es la ecuación difeencial del movimiento
Ecuaciones de Lagrange
Hemos calculado la energía cinética y potencial del sistema formado por el aro y la partícula. La lagrangiana L es
La ecuación del movimiento es
Obtenemos la misma ecuación diferencial del movimiento
Resolvemos la ecuación diferencial por procedimientos numéricos, con las siguientes condiciones iniciales: en el instante t=0, θ=θ0, dθ/dt=0. Interrumpimos el procedimiento poco antes de llegar a θ=0, ya que en esa posición la aceleración angular se hace infinita y el procedimiento numérico deja de funcionar
function cicloide_mas R=1; %radio alfa=pi/2; %ángulo inicial f=@(t,x) [x(2); -(9.8/R+x(2)^2)*cos(x(1)/2)/(2*sin(x(1)/2))]; opts=odeset('events',@stop_cicloide); [t,x]=ode45(f,[0,3],[alfa,0], opts); plot(t,x(:,1)) grid on xlabel('t') ylabel('\theta'); title('Oscilador') function [value,isterminal,direction]=stop_cicloide(t,x) value=x(1)-pi/360; %se detiene muy cerca del origen (a 0.5 grados) isterminal=1; direction=-1; % 1 crece, -1 decrece, 0 no importa end end
El tiempo de un cuarto de periodo es 0.5631 s, el periodo es 2.2524 s, el mismo valor que se obtendría por la fórmula del periodo
Actividades
Se introduce
- El ángulo inicial θ0, en el control titulado Ángulo inicial
- El radio del aro se ha fijado en R=1 m
Se pulsa el botón titulado Nuevo
Se observa las oscilaciones de la partícula situada en el borde del aro que rueda sin deslizar.
Oscilaciones de un aro con un pequeño disco fijado en su interior

La situación descrita en el primer apartado no es realista, por lo que vamos a estudiar un sistema análogo formado por un aro de masa M, radio R y momento de inercia IM y un disco sujeto al interior del aro cuyo centro dista r del centro del aro, su masa es m y su momento de inercia Im
El aro rueda sin deslizar a lo largo del plano horizontal, cuando el aro se ha desplazado Rθ, ha girado un ángulo θ. La posición del centro del disco y del aro son, respectivamente
La posición, velocidad y aceleración del centro de masas es
Ecuación del movimiento

Las fuerzas sobre el sistema formado por el aro y el disco son:
- El peso, mg en el centro del disco y Mg en el centro del aro
- La reacción del plano horizontal, N en el punto de contacto
- La fuerza de rozamiento, Fr necesaria para que el aro ruede sin deslizar
Las ecuaciones del movimiento son
- Movimiento en la dirección horizontal
- Movimiento en la dirección vertical
- Rotación alrededor de un eje perpendicular al plano del aro y que pasa por el centro de masas
El momento de inercia del sistema formado por el aro y el disco, respecto del centro de masas es
La fuerza de rozamiento Fr vale.
La reacción del plano N
La ecuación del movimiento
La ecuación diferencial del movimiento es
Consideremos el sistema formado por el aro y la partícula, el aro tiene masa M=0 y momento de inercia IM=0 y la partícula de masa m y momento de inercia Im=0 situada en la interior del aro r=R
La misma ecuación del movimiento obtenida en el primer apartado
Energía
Energía cinética
Energía potencial. La energía potencial del aro no cambia
Ecuaciones de Lagrange
La lagrangiana
La ecuación del movimiento es
Obtenemos la misma ecuación diferencial del movimiento
Periodo
El máximo desplazamiento se produce cuando la velocidad angular dθ/dt=0, el ángulo girado es θ0, la energía cinética es nula y la potencial E=mgr(1-cosθ0) es igual a la total. Despejamos la velocidad angular
El periodo es cuatro veces el tiempo que tarda en desplazarse desde 0 a θ0
Hacemos el cambio de variable
La integral se transforma
El máximo desplazamiento θ0 ha de ser menor que π
Representamos el cociente P/P0 en función del máximo desplazamiento θ0, para varios valores del parámetro k. θ0 tiene que ser algo menor que π para no producir error
th0=linspace(0,3,50); P=zeros(1,length(th0)); hold on for k=[0,0.5,1,2,8] i=1; for th_0=th0 f=@(x) sqrt((1+k*sin(x).^2*sin(th_0/2)^2)./(1-sin(x).^2*sin(th_0/2)^2)); P(i)=integral(f,0,pi/2)*2/pi; i=i+1; end plot(th0,P,'displayName',num2str(k)) end set(gca,'XTick',0:pi/6:pi) set(gca,'XTickLabel',{'0','\pi/6','\pi/3','\pi/2','2\pi/3','5\pi/6','\pi'}) xlabel('\theta_0') ylabel('P/P_0') legend('-DynamicLegend','location','best') grid on title('Periodo')
Cuando el momento de inercia I→0 o bien, k→∞, se obtiene el periodo para el sistema formado por el aro y la partícula estudiado en el primer apartado
Movimiento del sistema
Resolvemos la ecuación diferencial del movimiento para el siguiente sistema
- masa del aro, M=25 g
- radio del aro, R= 5 cm
- momento de inercia, IM=MR2
- masa del disco, m=30 g
- distacia del centro del disco al centro del aro, r=R-0.8 cm
- radio del disco, 0.8 cm
- momento de inercia, Im=m·0.82/2 g·cm2
Representamos el desplazamiento angular θ en función del tiepo t durante un periodo P. Las condiciones iniciales son: en el instante t=0, el aro parte de la posición de máximo desplzamiento θ0 con velocidad angular dθ/dt=0
M=25/1000; %masa aro R=5/100; %radio aro m=30/1000; %masa disco r=R-8/1000; %posición disco IM=M*R^2; Im=m*(8/1000)^2/2; %periodo I=Im+IM+M*R^2+m*(R-r)^2; k=4*m*r*R/I; th_0=pi/2; g=@(x) sqrt((1+k*sin(x).^2*sin(th_0/2)^2)./(1-sin(x).^2*sin(th_0/2)^2)); P=4*sqrt(I/(m*9.8*r))*integral(g,0,pi/2); %movimiento I0=Im+IM+m*r^2+(m+M)*R^2; f=@(t,x) [x(2); -m*r*(9.8+R*x(2)^2)*sin(x(1))/(I0-2*m*r*R*cos(x(1)))]; [t,x]=ode45(f,[0,P],[th_0,0]); plot(t,x(:,1)) grid on xlabel('t') ylabel('\theta'); title('Oscilador')
Verificamos que la energía se conserva
>> E0=m*9.8*r*(1-cos(th_0)) %inicial E0 = 0.0123 >> E=(I0-2*m*r*R*cos(x(:,1))).*x(:,2).^2/2+m*9.8*r*(1-cos(x(:,1))) E = 0.0123 0.0123 ... 0.0123 0.0123
Reacción N
Representamos la reacción del plano horizontal N/((M+m)g) en función del desplazamiento angular θ
De la conservación de la energía, obtenemos el cuadrado de la velocidad angular (dθ/dt)2. De la ecuación diferencial del movimiento, obtenemos la aceleración angular d2θ/dt2,
M=25/1000; %masa aro R=5/100; %radio aro m=30/1000; %masa disco r=R-8/1000; %posición disco IM=M*R^2; Im=m*(8/1000)^2/2; I0=Im+IM+m*r^2+(m+M)*R^2; hold on for th_0=[pi/4,pi/2, 3*pi/4] %cuadrado de la velocidad vel2=@(x) 2*m*9.8*r*(cos(x)-cos(th_0))./(I0-2*m*r*R*cos(x)); %aceleración acel=@(x) -m*r*(9.8+R*vel2(x)).*sin(x)./(I0-2*m*r*R*cos(x)); %reacción N N=@(x) ((m+M)*9.8+m*r*cos(x).*vel2(x)+m*r*sin(x).*acel(x))/((M+m)*9.8); fplot(N,[-th_0,th_0]) end hold off set(gca,'XTick',-pi:pi/4:pi) set(gca,'XTickLabel',{'-\pi','-3\pi/4','-\pi/2','-\pi/4','0','\pi/4', '\pi/2','3\pi/4','\pi'}) legend('\pi/4', '\pi/2','3\pi/4','location','best') grid on xlabel('\theta_0') ylabel('N/((M+m)g)'); title('Reacción N')
La característica más importante, es el máximo pronunciado de la reacción N cuando el disco pasa por el origen, θ=0, posición de equilibrio estable
Fuerza de rozamiento Fr
M=25/1000; %masa aro R=5/100; %radio aro m=30/1000; %masa disco r=R-8/1000; %posición disco IM=M*R^2; Im=m*(8/1000)^2/2; I0=Im+IM+m*r^2+(m+M)*R^2; hold on for th_0=[pi/4,pi/2, 3*pi/4] %cuadrado de la velocidad vel2=@(x) 2*m*9.8*r*(cos(x)-cos(th_0))./(I0-2*m*r*R*cos(x)); %aceleración acel=@(x) -m*r*(9.8+R*vel2(x)).*sin(x)./(I0-2*m*r*R*cos(x)); %Fuerza de rozamiento Fr=@(x) (m*r*sin(x).*vel2(x)+((M+m)*R-m*r*cos(x)).*acel(x))/((M+m)*9.8); fplot(Fr,[-th_0,th_0]) end hold off set(gca,'XTick',-pi:pi/4:pi) set(gca,'XTickLabel',{'-\pi','-3\pi/4','-\pi/2','-\pi/4','0','\pi/4', '\pi/2','3\pi/4','\pi'}) legend('\pi/4', '\pi/2','3\pi/4','location','best') grid on xlabel('\theta_0') ylabel('F_r/((M+m)g)'); title('Fuerza de rozamiento F_r')
La fuerza de rozamiento suele ser negativa para θ>0, tal como se dibuja en la figura de las fuerzas, y positiva para θ<0. Se opone al movimiento de traslación del centro de masas. Para pequeños desplazamientos angulares iniciales θ0, tiene un comportamiento lineal (línea de color azul). Para desplazamientos angulares iniciales grandes, el comportamiento es diferente
El análisis de la fuerza de rozamiento Fr y la reacción N es importante, para determinar el ángulo a partir del cual, el sistema formado aro y el disco empieza a deslizar, cuando la fuerza de rozamiento es igual o mayor que μsN
Referencias
Levy-Leblond J.M. Rock and roll: Non-isochronous small oscillations (an example). Am. J. Phys. 46 (1) January 1978, pp. 106-107
A. Filipponi. A non-isochronous rocking oscillator. Am. J. Phys. 82 (12), December 2014, pp. 1142-1148