Escalera apoyada en una pared móvil

En la página titulada, Escalera apoyada en dos paredes perpendiculares, se estudia el movimiento de una escalera cuyos extremos se apoyan en dos superficies lisas. El movimiento de la escalera desde la posición inicial hasta que impacta con el plano horizontal, se divide en dos etapas

El ángulo crítico, que forma la escalera con la vertical, es

θ l =arccos( 2 3 cos θ 0 )

Donde θ0 es el ángulo inicial de partida en reposo

Consideremos una pared móvil de masa M situada en el origen y una escalera de masa m y longitud 2l apoyada en en suelo y en la pared tal como se muestra en la figura, cuando el ángulo que forma la escalera con la vertical es θ0. Supondremos que las superficies sobre las que deslizan los extremos de la escalera son lisas (sin rozamiento). Se suelta la escalera, sus extremos empiezan a deslizar hasta que llega al suelo. Vamos a estudiar el movimiento de los dos cuerpos que forman este sistema aislado.

La escalera está apoyada en la pared móvil

En el instante t, la posición de la pared es x, la posición del centro de masas de la escalera es

r { x c =x+lsinθ y c =lcosθ v { d x c dt = dx dt +lcosθ dθ dt d y c dt =lsinθ dθ dt

Donde θ es el ángulo que forma la escalera con la vertical

La energía cinética de la pared es

E k1 = 1 2 M ( dx dt ) 2

La energía cinética de la escalera es la suma de la energía cinética de traslación del c.m. (centro de masas) y de la energía cinética de rotación alrededor de un eje perpendicular que pasa por el c.m.

E k2 = 1 2 m{ ( d x c dt ) 2 + ( d y c dt ) 2 }+ 1 2 ( 1 3 m l 2 ) ( dθ dt ) 2 = 1 2 m{ ( dx dt ) 2 +2lcosθ dx dt dθ dt }+ 2 3 m l 2 ( dθ dt ) 2

El momento de inercia de la escalera de masa m y longitud 2l respecto de un eje perpendicular que pasa por el c. m. es m(2l)2/12=ml2/3

La pared se mueve en el plano horizontal, su energía potencial no cambia. La energía potencial del c.m. de la escalera es Ep=mglcosθ.

La lagrangiana L=Ek1+Ek2-Ep es

L= 1 2 (m+M) ( dx dt ) 2 +mlcosθ dx dt dθ dt + 2 3 m l 2 ( dθ dt ) 2 mglcosθ

La lagrangiana no depende de x, por lo que hay una constante del movimiento que denominamos P

L x ˙ =cte (m+M) dx dt +mlcosθ dθ dt =P

Integramos, para obtener la posición de la pared móvil x en función del tiempo t y de la posición angular θ

x= P m+M t ml m+M θ 0 θ cosθ·dθ x= P m+M t ml m+M ( sinθsin θ 0 )

Como es un sistema aislado, la energía se mantiene constante

E= 1 2 (m+M) ( dx dt ) 2 +mlcosθ dx dt dθ dt + 2 3 m l 2 ( dθ dt ) 2 +mglcosθ E= 1 2 (m+M) { 1 m+M ( Pmlcosθ dθ dt ) } 2 +mlcosθ{ 1 m+M ( Pmlcosθ dθ dt ) } dθ dt + 2 3 m l 2 ( dθ dt ) 2 +mglcosθ E= P 2 2(m+M) +( 2 3 m cos 2 θ 2(m+M) )m l 2 ( dθ dt ) 2 +mglcosθ

Expresamos la lagrangiana en términos del ángulo θ, del mismo modo que la energía

L= P 2 2(m+M) +( 2 3 m cos 2 θ 2(m+M) )m l 2 ( dθ dt ) 2 mglcosθ

La ecuación del movimiento es

d dt ( L θ ˙ ) L θ =0 ( 4 3 m m+M cos 2 θ )l d 2 θ d t 2 + m m+M lsinθcosθ ( dθ dt ) 2 gsinθ=0

Se resuelve esta ecuación diferencial por procedimientos numéricos con las siguientes condiciones iniciales: en el instante t=0, la escalera está inclinada un ángulo θ0 y se suelta, (dθ/dt)0=0. La pared móvil está situada en el origen en reposo, x=0, dx/dt=0.

Estas condiciones iniciales hacen que la constante P=0 y la energía inicial E=mglcosθ0

Fuerza de interacción

La escalera ejerce una fuerza Fx sobre la pared móvil que hace que se mueva hacia la izquierda

F x =M d 2 x d t 2 = Mm m+M l( sinθ ( dθ dt ) 2 cosθ d 2 θ d t 2 )

Situación límite

La escalera va cayendo, deslizando su extemo superior por la pared móvil hasta que dejan de estar en contacto, en ese momento tl, la fuerza Fx=0, la aceleración de la pared es nula. A partir de ese instante, la pared móvil y el c.m. de la escalera se mueven con velocidad constante

d 2 θ d t 2 = sinθ cosθ ( dθ dt ) 2

Sustituyendo en la ecuación del movimiento, obtenemos la velocidad angular en el instante tl

( 4 3 m m+M cos 2 θ )l sinθ cosθ ( dθ dt ) 2 + m m+M lsinθcosθ ( dθ dt ) 2 gsinθ=0 ( dθ dt ) 2 = 3gcosθ 4l

Sustituyendo en la expresión de la energía, obtendremos el ángulo límite θl que forma la escalera con la vertical en el instante tl

E= P 2 2(m+M) +( 2 3 m cos 2 θ 2(m+M) )m l 2 3gcosθ 4l +mglcosθ E= P 2 2(m+M) + 3 2 mglcosθ 3 8 m 2 gl m+M cos 3 θ

Las condiciones iniciales especificadas hacen que la constante P=0 y la energía inicial E=mglcosθ0. Igualamos la energía inicial y la energía en el instante límite tl

cos θ 0 = 3 2 cosθ 3 8 m m+M cos 3 θ

Obtenemos una ecuación cúbica en cosθ. La raíz real buscada cosθl estará comprendida entre 0 y 1

En el instante límite tl, la posición xl y velocidad vl de la pared móvil serán

x l = P m+M t l ml m+M ( sin θ l sin θ 0 ) dx dt | l = P m+M ml m+M cos θ l ( dθ dt ) l = P m+M ml m+M cos θ l 3gcos θ l 4l

Las condiciones iniciales especificadas hacen que la constante P=0

x l = ml m+M ( sin θ l sin θ 0 ) v l = ml m+M cos θ l 3gcos θ l 4l

En la situación incial, la pared móvil de masa M se encuentra en la posición x=0 y el centro de masas de la escalera de masa m en lsinθ0. En el instante tl, la pared móvil se encuentra en xl y el centro de masa de la escalera en lsinθl+xl. En un sistema aislado, si el centro de masas del conjunto de los dos cuerpos estaba inicialmente en resposo, sigue en reposo en la misma posición

M·0+mlsin θ 0 m+M = M x l +m( mlsin θ l + x l ) m+M x l = ml m+M ( sin θ l sin θ 0 )

Como el centro de masas del conjunto de los dos cuerpos está en reposo, se cumple que M·vl+m·vcm=0

Siendo vcm la velocidad horizontal del centro de masas de la escalera

Primera etapa del movimiento

Con estas condiciones iniciales la constante P=0

Calculamos el ángulo límite θl, cuando la escalera deja de estar en contacto con la pared móvil. Calculamos las raíces de la ecuación cúbica y elegimos la raíz real entre 0 y 1. Utilizamos la función raices_3 o alternativamente, el comando roots

Resolvemos la ecuación diferencial en θ, mediante ode45 e interrumpimos el proceso de integración cuando θ=θl definiendo la función stop_escalera

function escalera_bloque
    L=0.5; %longitud de la escalera, 2L
    M=2; %masa de la pared móvil
    m=1; %masa de la escalera
    th_0=pi/6; %ángulo inicial

%ángulo límite final
    raiz=raices_3([3*m/(8*(m+M)),0,-3/2,cos(th_0)]);
    for i=1:length(raiz)
        if raiz(i)>0 && raiz(i)<1
            th_L=acos(raiz(i));
            break;
        end
    end
    f=@(t,x) [x(2); (9.8/L-m*cos(x(1))*x(2)^2/(m+M))*sin(x(1))/
(4/3-m*cos(x(1))^2/(m+M))]; 
    opts=odeset('events',@(t,x) stop_escalera(t,x, th_L));
    [t,x]=ode45(f,[0,10],[th_0,0],opts);
    plot(t,x(:,1));
    
    %posición final e instante
    disp([th_L*180/pi,t(end)])
    %velocidad y posición final de la pared móvil
    vPared=-m*cos(th_L)*sqrt(3*9.8*L*cos(th_L))/(2*(m+M)); 
    xPared=-m*L*(sin(th_L)-sin(th_0))/(m*M);
    disp([xPared,vPared])

    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'})
    grid on
    xlabel('t')
    ylabel('\theta')
    title('Escalera, pared móvil')

%conservación de la energía
    E0=m*9.8*L*cos(th_0); %energía inicial
    %energía final
    E=m*L^2*(2/3-m*cos(x(end,1)).^2/(2*(m+M))).*x(end,2).^2+m*9.8*L*cos(x(end,1));
    disp([E0,E])

%detiene el proceso de integración ode45 cuando se alcanza el ángulo límite
    function [value,isterminal,direction]=stop_escalera(~,x, th_L)
        value=x(1)-th_L;  
        isterminal=1;
        direction=1; 
    end

%raíces de una ecuación cúbica
    function x = raices_3(p)
        for i=4:-1:1
            p(i)=p(i)/p(1);
        end
        Q=(p(2)*p(2)-3*p(3))/9;
        R=(2*p(2)^3-9*p(2)*p(3)+27*p(4))/54;
        x=zeros(3,1); %reserva memoria para un vector de tres elementos
        if (R*R)<(Q^3)
            tetha=acos(R/sqrt(Q^3));
            x(1)=-2*sqrt(Q)*cos(tetha/3)-p(2)/3;
            x(2)=-2*sqrt(Q)*cos((tetha+2*pi)/3)-p(2)/3;
            x(3)=-2*sqrt(Q)*cos((tetha-2*pi)/3)-p(2)/3;
        else
            A=-sign(R)*nthroot(abs(R)+sqrt(R*R-Q^3),3);
            if A==0
                B=0;
            else
                B=Q/A;
            end
            x(1)=(A+B)-p(2)/3;
            x(2)=-(A+B)/2-p(2)/3+(sqrt(3)*(A-B)/2)*1i; 
            x(3)=-(A+B)/2-p(2)/3-(sqrt(3)*(A-B)/2)*1i;
        end
    end
end
53.4949    0.2917
-0.0760   -0.2932
4.2435    4.2438

Los resultados son:

Representamos la posición angular θ de la escalera en función del tiempo t hasta el instante tl

Sea θ0=π/6, la posición angular inicial de la escalera. Si hacemos que la masa de la pared móvil sea grande, por ejemplo M=20, obtenemos un ángulo θl=54.6°, muy próximo al que se obtiene considerando una pared fija (véase la fórmula al principio de esta página)

>> acos(2*cos(pi/6)/3)*180/pi
ans =   54.7356

Segunda etapa del movimiento

La pared móvil de masa M se mueve hacia la izquierda con velocidad constante vl

La constancia del momento lineal en un sistema aislado, implica que el centro de masa de la escalera se mueve con velocidad constante vcm hacia la derecha, tal que, m·vcm+Mvl=0

En el sistema de referencia que se mueve con el centro de masa de la escalera, su posición y velocidad es

y c =lcosθ d y c dt =lsinθ dθ dt

La energía cinética de la escalera es la suma de la energía cinética de traslación del c.m. y de la energía cinética de rotación alrededor de un eje perpendicular que pasa por el c.m.

E k = 1 2 m ( dy dt ) 2 + 1 2 ( 1 3 m l 2 ) ( dθ dt ) 2 = 1 2 m l 2 ( 1 3 + sin 2 θ ) ( dθ dt ) 2

La energía potencial del c.m. de la escalera es Ep=mglcosθ. La lagrangiana es

L= 1 2 m l 2 ( 1 3 + sin 2 θ ) ( dθ dt ) 2 mglcosθ

La ecuación del movimiento

d dt ( L θ ˙ ) L θ =0 ( 1 3 + sin 2 θ ) d 2 θ d t 2 +sinθcosθ ( dθ dt ) 2 g l sinθ=0

La misma que hemos obtenido en la página titulada Escalera apoyada en dos paredes perpendiculares, en la sección, 'El extremo B de la escalera deja de estar en contacto con la pared vertical'

Resolvemos esta ecuación diferencial por procedimientos numéricos con las siguientes condiciones iniciales: en el instante tl, el ángulo que forma la escalera con la vertical es θl y su velocidad angular de rotación es

( dθ dt ) l = 1 2 3gcos θ l l

Cuando la escalera impacta con el suelo, θ=π/2, emplea un tiempo tf en alcanzarlo. La posición de la pared móvil y del c.m. de la escalera serán

x f = x l + v l ( t f t l ) x cm =lsin θ l + v cm ( t f t l )=lsin θ l M m v l ( t f t l )

ya que durante este intervalo de tiempo se mueven con velocidad constante vl y vcm, respectivamente

Movimiento completo

Resolvemos la ecuación diferencial en θ, mediante ode45 e interrumpimos el proceso de integración cuando θ=π/2 definiendo la función stop_escalera_1

function escalera_bloque
L=0.5; %longitud de la escalera 2L
M=2; %masa de la pared móvil
m=1; %masa de la escalera
th_0=pi/6; %ángulo inicial

%ángulo límite
raiz=raices_3([3*m/(8*(m+M)),0,-3/2,cos(th_0)]);
for i=1:length(raiz)
    if raiz(i)>0 && raiz(i)<1
        th_L=acos(raiz(i));
        break;
    end
end

%primera etapa
f=@(t,x) [x(2); (9.8/L-m*cos(x(1))/(m+M))*sin(x(1))/(4/3-m*cos(x(1))^2/(m+M))]; 
opts=odeset('events',@(t,x) stop_escalera(t,x, th_L));
[t,x]=ode45(f,[0,10],[th_0,0],opts);
hold on
plot(t,x(:,1));
tL=t(end); %final de la primera etapa
%velocidad final de la pared móvil
vPared=-m*cos(th_L)*sqrt(3*9.8*L*cos(th_L))/(2*(m+M)); 
line([tL,tL],[0,th_L],'color','k','lineStyle','--')
line([0,tL],[th_L,th_L],'color','k','lineStyle','--')

%segunda etapa
opts=odeset('events',@stop_escalera_1);
f=@(t,x) [x(2);(9.8/L-cos(x(1))*x(2)^2)*sin(x(1))/(1/3+sin(x(1))^2)]; 
[t,x]=ode45(f,[0,10],[th_L,sqrt(3*9.8*cos(th_L)/L)/2],opts);
plot(t+tL,x(:,1))
disp(t(end)+tL) %tiempo total
hold off
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'})
grid on
xlabel('t')
ylabel('\theta')
title('Escalera, pared móvil')

%conservación de la energía
E0=m*9.8*L*cos(th_0);
E=M*(m+M)*vPared^2/(2*m)+m*L^2*(1/3+sin(x(end,1))^2)*x(end,2).^2
/2+m*9.8*L*cos(x(end,1));
disp([E0,E])
%velocidad angular final
wf=3*9.8*(cos(th_0)-3*M*cos(th_L)^3/(8*(m+M)))/(2*L);
disp([sqrt(wf),x(end,2)])

%detiene el proceso de integración ode45 cuando se alcanza el ángulo límite
function [value,isterminal,direction]=stop_escalera(~,x, th_L)
    value=x(1)-th_L;  
    isterminal=1;
    direction=1; 
end

%detiene el proceso de integración ode45 cuando se alcanza el ángulo pi/2
function [value,isterminal,direction]=stop_escalera_1(~,x)
    value=x(1)-pi/2;  
    isterminal=1;
    direction=1; 
end

%raíces de una ecuación cúbica
function x = raices_3(p)
    for i=4:-1:2
        p(i)=p(i)/p(1);
    end
    %p(4)=p(4)/p(1);
    %p(2)=p(2)/p(1);
    %p(3)=p(3)/p(1);
    p(1)=1;
    Q=(p(2)*p(2)-3*p(3))/9;
    R=(2*p(2)^3-9*p(2)*p(3)+27*p(4))/54;
    x=zeros(3,1); %reserva memoria para un vector de tres elementos
    if (R*R)<(Q^3)
        tetha=acos(R/sqrt(Q^3));
        x(1)=-2*sqrt(Q)*cos(tetha/3)-p(2)/3;
        x(2)=-2*sqrt(Q)*cos((tetha+2*pi)/3)-p(2)/3;
        x(3)=-2*sqrt(Q)*cos((tetha-2*pi)/3)-p(2)/3;
    else
        A=-sign(R)*nthroot(abs(R)+sqrt(R*R-Q^3),3);
        if A==0
            B=0;
        else
            B=Q/A;
        end
        x(1)=(A+B)-p(2)/3;
        x(2)=-(A+B)/2-p(2)/3+(sqrt(3)*(A-B)/2)*1i; 
        x(3)=-(A+B)/2-p(2)/3-(sqrt(3)*(A-B)/2)*1i;
    end
end
end

Representamos la posición angular θ de la escalera en función del tiempo t desde el instante t=0, hasta el instante tl y desde el instante tl hasta el instante tf cuando la escalera ha impactado en el suelo θ=π/2

La energía en la segunda etapa del movimiento es

E= 1 2 m v cm 2 + 1 2 M v l 2 + 1 2 m l 2 ( 1 3 + sin 2 θ ) ( dθ dt ) 2 +mglcosθ= 1 2 M( 1+ M m ) v l 2 + 1 2 m l 2 ( 1 3 + sin 2 θ ) ( dθ dt ) 2 +mglcosθ

Igualamos la energía inicial mglcosθ0 a la final, cuando θ=π/2, despejamos la velocidad angular dθ/dt en la posición final

mglcos θ 0 = 1 2 M( 1+ M m ) v l 2 + 2 3 m l 2 ( dθ dt ) 2 ( dθ dt ) 2 = 3g 2l ( cos θ 0 3 8 M ( M+m ) cos 3 θ l )

Comprobamos que la energía incial es similar a la final y que la velocidad angular dθ/dt final es similar a la que nos proporciona el procedimiento numérico x(end,2)

    4.2435    4.2417 %energía
    4.8902    4.8890 %velocidad angular final

Actividades

Se introduce

El programa interactivo ha fijado los parámetros

Se pulsa el botón titulado Nuevo

Observamos que la escalera empieza a caer, a la vez que la pared se mueve hacia la izquierda del origen

Una flecha de color azul muestra la fuerza Fx que ejerce la escalera sobre la pared móvil, mientras la escalera está en contacto con dicha pared. Medimos el ángulo critico θl cuando dicha fuerza se hace cero, el mensaje 'En contacto' desaparece

En la parte derecha, se proporcionan los datos de

Después de que la escalera impacte con el suelo, la pared móvil y la escalera siguen moviéndose con velocidad constante, aunque en la simulación parece que se detengan

Referencias

Daniel Arovas. Lecture Notes on Classical Mechanics (A Work in Progress). University of California, San Diego. May 8, 2013, pp. 393-396. https://courses.physics.ucsd.edu/2010/Fall/physics200a/LECTURES/200_COURSE.pdf