Caída de una varilla inclinada

Movimiento de la partícula

Ya hemos estudiado el movimiento de caída de los cuerpos en el capítulo de Cinemática. La partícula se deja caer desde una altura y0=γLcos θ0  partiendo del reposo

y= y 0 1 2 g t 2

El tiempo Pp que tarda en llegar al suelo y=0 es

P p = 2 y 0 g = 2γLcos θ 0 g

Movimiento de rotación la varilla

La varilla es un sólido rígido en rotación alrededor de un eje fijo perpendicular a la varilla y que pasa por su extremo O.

Las fuerzas que actúan sobre la varilla son:

La ecuación de la dinámica de la rotación alrededor de un eje fijo, es

I 0 α=M I 0 d 2 θ d t 2 =M

Para obtener el ángulo θ que hace la varilla con el suelo en función del tiempo t, se integra la ecuación diferencial mediante procedimientos numéricos.

d 2 θ d t 2 3g 2L sinθ=0

con las condiciones iniciales siguientes: en el instante t=0, la varilla está inclinada un ángulo θ=θ0 (10°) y parte del reposo, ω=dθ/dt=0

L=2; %longitud de la varilla
k=3*9.8/(2*L);
f=@(t,x) [x(2);k*sin(x(1))]; 
ang=10*pi/180; %ángulo inicial
opts=odeset('events',@varilla_ode45);
[t,x,te]=ode45(f,[0,4],[ang,0], opts);
plot(t,x(:,1))
set(gca,'YTick',0:pi/12:pi/2)
set(gca,'YTickLabel',{'0','\pi/12','\pi/6','\pi/4','pi/3','5\pi/12','\pi/2'})
text(0.1,pi/2-0.2,num2str(te)) %tiempo de caída
grid on
xlabel('t')
ylabel('\theta');
title('Posición angular en función del tiempo')

Definimos una función para que el proceso de integración se detenga cuando la posición de la varilla sea π/2, el instante t en el que ocurre se guarda en la variable te

function [detect,stopin,direction]=varilla_ode45(~,x)
    detect=x(1)-pi/2; 
    stopin=1;
    direction=1; 
end

Estudio energético

La energía potencial de la varilla es la energía potencial de una partícula de masa m situada en el c.m. de la varilla E=mg(L/2)·cosθ0. La energía potencial se convierte en energía cinética de rotación. El principio de conservación de la energía se escribe.

E= 1 2 ( 1 3 m L 2 ) ω 2 +mg L 2 cosθ

Comprobamos que la energía se mantiene aproximadamente constante e igual a la inicial, mgLcos(θ0)/2

>> L*9.8*cos(ang)/2
ans =    9.6511
>> E=L^2*x(:,2).^2/6+9.8*L*cos(x(:,1))/2
    9.6511
    9.6511
    ....
    9.6603
    9.6611
    9.6622

La ecuación del movimiento nos permite calcular la posición angular θ y la velocidad angular ω de la varilla en función del tiempo t. Sin embargo, el principio de conservación de la energía nos permite calcular la velocidad angular ω en función de la posición θ.

ω= dθ dt = 3g L (cos θ 0 cosθ)

Supongamos que la varilla parte de la posición θ0≈0, (la posición θ0=0 es de equilibrio inestable)

dθ dt = 6g L sin( θ 2 )

Integramos, véase la la tabla de integrales

θ 0 0 θ dθ sin( θ 2 ) = 6g L 0 t dt 2log( tan( θ 4 ) )2α= 6g L t tan( θ 4 )=exp( 3g 2L t+α ) θ=4arctan{ exp( 3g 2L t+α ) }

La varilla llega al suelo θ=π/2, en el instante tf

t f = tan( π 8 )α 3g 2L

Comparamos la solución numérica y la analítica aproximada cuando la varilla se deja caer desde la posición angular casi vertical, θ0=1° (π/180)

function varilla
    L=2; %longitud de la varilla
    ang=1*pi/180;  %posición de partida
    k=3*9.8/(2*L);
    f=@(t,x) [x(2);k*sin(x(1))]; 
    opts=odeset('events',@varilla_ode45);
    [t,x, te]=ode45(f,[0,4],[ang,0], opts);
    hold on
    plot(t,x(:,1))
    set(gca,'YTick',0:pi/12:pi/2)
    set(gca,'YTickLabel',{'0','\pi/12','\pi/6','\pi/4','pi/3','5\pi/12','\pi/2'})
  
    alfa=log(tan(ang/4));
    tf=(log(pi/8)-alfa)/sqrt(9.8*3/(2*L));
    th=@(t) 4*exp(sqrt(9.8*3/(2*L))*t+alfa);
    fplot(th,[0,tf])
    hold off
    disp([te,tf])
    grid on
    legend('numérica', 'analítica','location','best')
    xlabel('t')
    ylabel('\theta');
    title('Posición angular en función del tiempo')

    function [detect,stopin,direction]=varilla_ode45(~,x)
        detect=x(1)-pi/2; 
        stopin=1;
        direction=1; 
    end
end

Los tiempos calculados son, respectivamente

    1.9352    1.6598

El primero, 1.9352 s, es el que más se aproxima al tiempo exacto, que vamos a calcular en el próximo apartado

Tiempo que tarda la varilla en caer al suelo

Integrando entre θ0 y π/2, obtenemos el tiempo Pv que tarda la varilla en llegar al suelo.

P v = L 3g θ 0 π/2 dθ cos θ 0 cosθ

Empleando la fórmula del coseno del ángulo doble cos2A=cos2A-sin2A y la relación sen2A+cos2A=1

P v = L 6g θ 0 π/2 dθ cos 2 ( θ 0 2 ) cos 2 ( θ 2 )

Haciendo el cambio de variable

sinφ= cos( θ 2 ) k k=cos( θ 0 2 ) dθ= 2kcosφ·dφ sin( θ 2 ) = 2kcosφ·dφ 1 k 2 sin 2 φ

Llegamos finalmente, a la expresión deseada para el tiempo Pv que tarda la varilla en llegar al suelo cuando se suelta en la posición angular θ0. Una expresión similar a la obtenida para el periodo de un péndulo para cualquier amplitud.

P v = 2L 3g φ 0 π/2 dφ 1 k 2 sin 2 φ sin φ 0 = 1 2 k

La integral se puede calcular numéricamente o se puede encontrar en tablas, si la escribimos como diferencia entre la integral elíptica completa cuyos límites son 0 y π/2 y la integral elíptica de primera especie de límites 0 a φ0.

P v = 2L 3g ( 0 π/2 dφ 1 k 2 sin 2 φ 0 φ 0 dφ 1 k 2 sin 2 φ )sin φ 0 = 1 2 k

El código MATLAB calcula la diferencia entre las dos integrales elípticas, cuando introducimos la posición inicial θ0 de la varilla en grados.

angulo=10;  %ángulo en grados
L=2; %longitud de la varilla
k=cos(angulo*pi/360);
phi=asin(1/(sqrt(2)*k));
t=(ellipke(k^2)-ellipticF(phi,k^2))*sqrt(2*L/(3*9.8));
fprintf('El tiempo de caída es %1.4f\n',t);
El tiempo de caída es 1.0866

Comparación de los dos movimientos

Comparación de aceleraciones

Comparamos la aceleración de la partícula cuando se encuentra en la posición P a una distancia γL del eje de rotación, con la aceleración inicial del punto P de la varilla. En la figura se muestran ambas aceleraciones

El punto P de la varilla describe una circunferencia de radio γL. Hallamos su aceleración tangencial at.

a t =γLα= 3g 2 γ·sinθ

La aceleración de la gravedad en dicha dirección es sinθ

Si la primera es mayor que la segunda, se cumplirá que la varilla cae más deprisa que la partícula.

3g 2 γ·sinθ>g·sinθγ> 2 3

Como vamos a comprobar, esta es una condición necesaria pero no suficiente para que la varilla llegue antes al suelo que la partícula.

Comparación de tiempos de vuelo

Comparemos los tiempos que requiere cada cuerpo para alcanzar el suelo

La partícula sale de la posición y0=γLcosθ0 y llega al suelo en el instante Pp

P p = 2γLcos θ 0 g

La varilla parte de la posición angular inicial θ0 y tarda un tiempo Pv

P v = 2L 3g I e

Donde Ie es la diferencia entre los valores de las dos integrales elípticas de primera especie del apartado anterior.

La varilla llegará antes que la partícula si se cumple  que Pp>Pv, o bien

γ< I e 2 3cos θ 0

La posición inicial de la partícula γL depende de la posición angular inicial de la varilla θ0, para que la varilla llegue antes que la partícula. En la siguiente tabla:

θ0 Ie γc
30 1.8454 1.31
40 1.5467 1.04
50 1.3009 0.878
60 1.0783 0.775
70 0.8547 0.712
80 0.5941 0.6775
89.99 0.0187 0.667≈2/3

Los números en las columnas de la tabla nos sugieren las siguientes conclusiones:

Para cada posición angular inicial θ0, se tiene que

Como γ es menor que la unidad, para las posiciones angulares θ0=30º, 40º, la partícula siempre llega antes que la varilla. Como podemos comprobar, mediante el programa interactivo, el ángulo mínimo para el cual γc=1 es θ0≈42º

Para un valor dado γ<0.666…=2/3, la varilla llegará al suelo después de la partícula cualquiera que sea el ángulo inicial θ0 de partida tal como puede apreciarse en la tabla.

Actividades

Se introduce

Se pulsa en el botón titulado Nuevo

Se observa el movimiento de rotación de la varilla y el movimiento de la partícula. Investigar para que posiciones y ángulos de partida la varilla llega al suelo antes que la partícula.


Fuerzas Fx y Fy que se ejercen sobre la varilla en su eje de rotación.

El centro de masas de la varilla describe un movimiento circular de radio L/2. Por tanto, la aceleración del c.m. tiene dos componentes:

En la parte izquierda de la figura, se muestra las direcciones de la aceleración tangencial y normal cuando la varilla hace un ángulo θ con la vertical. En la parte derecha de la figura, se muestra las componentes rectangulares de dichas aceleraciones.

Aplicando la segunda ley de Newton escribimos

m(at·cosθ-an·sinθ)=Fx
m
(an·cosθ+at·sinθ)=mg-Fy

a t =α L 2 = 3g 4 sinθ a n = ω 2 L 2 = 3 2 g(cos θ 0 cosθ)

Dado el ángulo θ, despejamos Fx y Fy del sistema de ecuaciones.

F x = 3mg 4 sinθ(3cosθ2cos θ 0 ) F y =mg 3mg 4 (1+2cosθcos θ 0 3 cos 2 θ)

Cuando la varilla parte de una posición casi vertical, θ0≈0, la conponente vertical vale, aproximadamente

F y mg( 1 4 3 2 cosθ+ 9 4 cos 2 θ )= mg 4 ( 3cosθ1 ) 2

Se anula para el ángulo tal que cosθ=1/3, θ=70.5°

La componente horizontal vale, aproximadamente, F x = 2 2 mg

Ejemplo:

Calcular Fx y Fy cuando θ=60º

Primero calculamos las componentes de la aceleración

at=6.37 m/s2
an
=3.04 m/s2

Luego, calculamos las componentes de la fuerza que se ejerce sobre la varilla en el eje

Fx=0.546 N
Fy
=2.765 N

Fuerzas y momento internos

En Internet encontramos vídeos que muestra la caída de una chimenea y observamos que se rompe a una cierta altura antes de llegar al suelo. Para predecir la altura es necesario conocer las fuerzas y momento internos.

Una chimenea es una estructura hecha de ladrillo o cemento de sección variable que no es fácil de analizar. En esta sección, nos limitaremos a calcular las fuerzas internas Pr y Pθ y el momento Nb que se ejercen sobre la porción r inferior de una varilla homogénea que cae desde la posición vertical

Cuando la varilla de masa m y longitud L, forma un ángulo θ con la dirección vertical, las fuerzas que actúan sobre la varilla son:

El centro de masas describe un movimiento circular de radio L/2.

La última ecuación nos proporciona la aceleración angular

d 2 θ d t 2 = 3 2 g L sinθ

Integramos esta ecuación diferencial, sabiendo que la velocidad angular w=dθ/dt=0 para θ=0, parte del reposo en la posición vertical

d 2 θ d t 2 = dω dt = dω dθ dθ dt =ω dω dθ 1 2 d ω 2 dθ = 3 2 g L sinθ ω 2 0=3 g L 0 θ sinθ·dθ ( dθ dt ) 2 =3 g L ( 1cosθ )

Alternativamente, aplicamos el principio de conservación de la energía al movimiento de rotación

mg L 2 = 1 2 ( 1 3 m L 2 ) ( dθ dt ) 2 +mg L 2 cosθ ( dθ dt ) 2 =3 g L ( 1cosθ )

Despejamos las componentes de la fuerza que ejerce el eje de rotación sobre la varilla.

F θ = 1 4 mgsinθ F r = 5 2 mg( cosθ 3 5 )

Fuerzas y momento interno

Vamos a estudiar el movimiento de una porción de varilla de longitud r y de masa m'=mr/L, el centro de masas está situado a r/2 del eje O

Las fuerzas que actúan sobre esta porción de varilla son externas e internas (que ejerce la porción L-r de varilla)

El centro de masas de la porción de varilla describe un movimiento circular de radio r/2.

Conocemos la velocidad y la aceleración angular de la varilla y las componentes de la fuerza en el eje de rotación, despejamos las fuerzas internas Pr, Pθ y el momento Nb

P θ = 3 4 mgsinθ( r 2 L 2 4 3 r L + 1 3 ) P r = 1 2 mg( 1 r L )( ( 5+3 r L )cosθ3( 1+ r L ) ) N b = 1 4 mgr ( 1 r L ) 2 sinθ

Comprobamos que Nb y Pθ están relacionados

P θ = N b r

Las posiciones importantes de la varilla son r/L=1/3 y 2/3 medidas desde el eje de rotación

Referencias

Para el apartado "Tiempo que tarda la varilla en caer al suelo"

Theron W. The "faster than gravity" demostration revisited. Am. J. Phys. 56 (8) August 1988, pp. 736-739

Peter F Hinrichsen. The physics of falling chimney stacks. Phys. Educ. 56 (2021) 055023

Gabriele Varieschi, Kaoru Kamiya. Toy models for the falling chimney. Am. J. Phys. 71 (10), October 2003, pp. 1025-1031