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:

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 , 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-R
cosθ

Las componentes de la velocidad partícula se obtienen derivando, respecto del tiempo t, las coordenadas x e y.

v x = dx dt =R(1cosθ) dθ dt v y = dy dt =Rsinθ dθ dt

Principio de conservación de la energía

La energía cinética es

E k = 1 2 m( v x 2 + v y 2 )=2m R 2 ( dθ dt ) 2 sin 2 θ 2

la energía potencial es

E p =mgy=2mgR sin 2 θ 2

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.

2m R 2 ( dθ dt ) 2 sin 2 θ 2 +2mgR sin 2 θ 2 =2mgR sin 2 θ 0 2

Separando las variables t y θ, la ecuación del movimiento se escribe

-sin(θ/2)dθ sin 2 ( θ 0 /2) sin 2 (θ/2) = g R dt

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

-sin(θ/2)dθ -cos 2 ( θ 0 /2)+ cos 2 (θ/2) = g R dt

Haciendo el cambio de variable

u= cos(θ/2) cos( θ 0 /2) u 0 u 2du u 2 1 =Ωt

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

ln| u+ u 2 1 |ln u 0 = Ωt 2

Como 0≤θ/2≤θ0/2≤π/2, se cumple que  u≥1

u+ u 2 1 =exp( Ωt 2 )

Despejamos la raíz cuadrada, elevamos ambos miembros al cuadrado, simplificamos y despejamos u.

u= exp(Ωt/2)+exp(Ωt/2) 2 =cosh( Ωt 2 )

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

exp( ΩP 8 )= 1+sin( θ 0 /2) cos( θ 0 /2) P=8 R g ln( 1+sin( θ 0 /2) cos( θ 0 /2) )

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

{ x c =RθRsinθ y c =RRcosθ { d x c dt =R( 1cosθ ) dθ dt d x c dt =Rsinθ dθ dt { d 2 x c d t 2 =Rsinθ ( dθ dt ) 2 +R( 1cosθ ) d 2 θ d t 2 d 2 y c d t 2 =Rcosθ ( dθ dt ) 2 +Rsinθ d 2 θ d t 2

Las fuerzas sobre el sistema formado por el aro y la partícula son:

Las ecuaciones del movimiento son

m d 2 x c d t 2 = F r m d 2 y c d t 2 =Nmg F r y c N( Rθ x c )= I c d 2 θ d t 2

El momento de inercia Ic=0, ya que la posición de la partícula coincide con el centro de masas

( Rsinθ ( dθ dt ) 2 +R( 1cosθ ) d 2 θ d t 2 )( RRcosθ )( g+Rcosθ ( dθ dt ) 2 +Rsinθ d 2 θ d t 2 )Rsinθ=0 2sin θ 2 d 2 θ d t 2 +( g R + ( dθ dt ) 2 )cos θ 2 =0

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

L= E k E p =2m R 2 ( dθ dt ) 2 sin 2 θ 2 2mgR sin 2 θ 2

La ecuación del movimiento es

d dt ( L θ ˙ ) L θ =0 d dt ( 4m R 2 ( dθ dt ) sin 2 θ 2 )2m R 2 ( dθ dt ) 2 sin θ 2 cos θ 2 +2mgRsin θ 2 cos θ 2 =0 4m R 2 ( d 2 θ d t 2 ) sin 2 θ 2 +4m R 2 ( dθ dt ) 2 sin θ 2 cos θ 2 2m R 2 ( dθ dt ) 2 sin θ 2 cos θ 2 +2mgRsin θ 2 cos θ 2 =0 2( d 2 θ d t 2 )sin θ 2 + ( dθ dt ) 2 cos θ 2 + g R cos θ 2 =0

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

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 , ha girado un ángulo θ. La posición del centro del disco y del aro son, respectivamente

{ x m =Rθrsinθ y m =Rrcosθ { x M =Rθ y M =R

La posición, velocidad y aceleración del centro de masas es

{ x c = m x m +M x M m+M =Rθ m m+M rsinθ y c = m y m +M y M m+M =R m m+M rcosθ { d x c dt =( R m m+M rcosθ ) dθ dt d x c dt = m m+M rsinθ dθ dt { d 2 x c d t 2 = m m+M rsinθ ( dθ dt ) 2 +( R m m+M rcosθ ) d 2 θ d t 2 d 2 y c d t 2 = m m+M rcosθ ( dθ dt ) 2 + m m+M rsinθ d 2 θ d t 2

Ecuación del movimiento

Las fuerzas sobre el sistema formado por el aro y el disco son:

Las ecuaciones del movimiento son

(m+M) d 2 x c d t 2 = F r (m+M) d 2 y c d t 2 =N(m+M)g F r y c N( Rθ x c )= I c d 2 θ d t 2

El momento de inercia del sistema formado por el aro y el disco, respecto del centro de masas es

I c = I m +m( ( x m x c ) 2 + ( y m y c ) 2 )+ I M +M( ( x M x c ) 2 + ( y M y c ) 2 )= I m + I M +m( ( M m+M rsinθ ) 2 + ( M m+M rcosθ ) 2 )+M( ( m m+M rsinθ ) 2 + ( m m+M rsinθ ) 2 ) I c = I m + I M + mM m+M r 2

La fuerza de rozamiento Fr vale.

F r =(m+M) d 2 x c d t 2 =mrsinθ ( dθ dt ) 2 ( (m+M)Rmrcosθ ) d 2 θ d t 2

La reacción del plano N

N=(m+M)g+mrcosθ ( dθ dt ) 2 +mrsinθ d 2 θ d t 2

La ecuación del movimiento

{ mrsinθ ( dθ dt ) 2 ( (m+M)Rmrcosθ ) d 2 θ d t 2 }( R m m+M rcosθ ) { (m+M)g+mrcosθ ( dθ dt ) 2 +mrsinθ d 2 θ d t 2 }( m m+M rsinθ )=( I m + I M + mM m+M r 2 ) d 2 θ d t 2 mgrsinθmRrsinθ ( dθ dt ) 2 ( m 2 m+M r 2 +(m+M) R 2 2mrRcosθ ) d 2 θ d t 2 =( I m + I M + mM m+M r 2 ) d 2 θ d t 2 ( I m + I M +m r 2 +(m+M) R 2 2mrRcosθ ) d 2 θ d t 2 +mRrsinθ ( dθ dt ) 2 +mgrsinθ=0

La ecuación diferencial del movimiento es

I(θ) d 2 θ d t 2 +mRrsinθ ( dθ dt ) 2 +mgrsinθ=0

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

( 2m R 2 2m R 2 cosθ ) d 2 θ d t 2 +m R 2 sinθ ( dθ dt ) 2 +mgRsinθ=0 4 sin 2 θ 2 d 2 θ d t 2 +2sin θ 2 cos θ 2 ( dθ dt ) 2 + g R 2sin θ 2 cos θ 2 =0 2sin θ 2 d 2 θ d t 2 +( g R + ( dθ dt ) 2 )cos θ 2 =0

La misma ecuación del movimiento obtenida en el primer apartado

Energía

Energía cinética

E k = 1 2 I c ( dθ dt ) 2 + 1 2 (m+M)( ( d x c dt ) 2 + ( d y c dt ) 2 )= 1 2 ( I m + I M + mM m+M r 2 ) ( dθ dt ) 2 + 1 2 (m+M)( ( ( R m m+M rcosθ ) ) 2 + ( m m+M rsinθ ) 2 ) ( dθ dt ) 2 = E k = 1 2 ( I m + I M +m r 2 +(m+M) R 2 2mRrcosθ ) ( dθ dt ) 2

Energía potencial. La energía potencial del aro no cambia

E p =mgr(1cosθ)

Ecuaciones de Lagrange

La lagrangiana

L= E k E p = 1 2 ( I m + I M +m r 2 +(m+M) R 2 2mRrcosθ ) ( dθ dt ) 2 mgr(1cosθ) = 1 2 I( θ ) ( dθ dt ) 2 mgr(1cosθ)

La ecuación del movimiento es

d dt ( L θ ˙ ) L θ =0 d dt ( I(θ) dθ dt ) 1 2 dI(θ) dθ ( dθ dt ) 2 +mgrsinθ=0 I(θ) d 2 θ d t 2 +2mRrsinθ ( dθ dt ) 2 mRrsinθ ( dθ dt ) 2 +mgrsinθ=0 I(θ) d 2 θ d t 2 +mRrsinθ ( dθ dt ) 2 +mgrsinθ=0

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

1 2 ( I m + I M +m r 2 +(m+M) R 2 2mRrcosθ ) ( dθ dt ) 2 +mgr(1cosθ)=mgr(1cos θ 0 ) ( dθ dt ) 2 = 2mgr(cosθcos θ 0 ) I m + I M +m r 2 +(m+M) R 2 2mRrcosθ ( dθ dt ) 2 = 2mgr(cosθcos θ 0 ) I+2mrR(1cosθ) ,I= I m + I M +M R 2 +m ( Rr ) 2 ( dθ dt ) 2 = mgr I 2(cosθcos θ 0 ) 1+ k 2 (1cosθ) ,k= 4mrR I

El periodo es cuatro veces el tiempo que tarda en desplazarse desde 0 a θ0

0 P dt =4 I mgr 0 θ 0 1+ k 2 (1cosθ) 2(cosθcos θ 0 ) P=4 I mgr 0 θ 0 1+k sin 2 θ 2 4( sin 2 θ 0 2 sin 2 θ 2 )

Hacemos el cambio de variable

sin θ 2 =sin(x)sin θ 0 2 1 2 cos θ 2 ·dθ=cos(x)sin θ 0 2 ·dx

La integral se transforma

P=4 I mgr 0 π/2 1+k sin 2 (x) sin 2 θ 0 2 4( sin 2 θ 0 2 sin 2 x· sin 2 θ 0 2 ) 2cosx·sin θ 0 2 1 sin 2 x· sin 2 θ 0 2 dx P=4 I mgr 0 π/2 1+k sin 2 (x) sin 2 θ 0 2 1 sin 2 x· sin 2 θ 0 2 dx P= P 0 2 π 0 π/2 1+k sin 2 (x) sin 2 θ 0 2 1 sin 2 x· sin 2 θ 0 2 dx, P 0 =2π I mgr

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

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 θ

N=(m+M)g+mrcosθ ( dθ dt ) 2 +mrsinθ d 2 θ d t 2

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

F r =mrsinθ ( dθ dt ) 2 ( (m+M)Rmrcosθ ) d 2 θ d t 2

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