Movimiento sobre una cúpula semiesférica

Sin rozamiento

La partícula se encuentra inicialmente en reposo sobre el vértice de la cúpula, en una posición de equilibrio inestable. Cuando se desvía ligeramente de esta posición, la partícula desliza sin rozamiento, incrementando su velocidad hasta que llega un momento en el que deja de tener contacto con la cúpula. En este apartado, calcularemos la posición θc para la cual la reacción N de superficie semiesférica es nula.

Conservación de la energía

La energía de la partícula en la posición inicial θ=0, es

Ei=mgR

La energía de la partícula en la posición θ es

E f = 1 2 m v 2 +mgRcosθ

Aplicando el principio de conservación de la energía Ei=Ef, calculamos la velocidad del móvil v en la posición θ

v 2 =2gR(1cosθ)v=2 gR sin( θ 2 )

Dinámica

Las fuerzas que actúan sobre la partícula son

La partícula describe un movimiento circular con aceleración tangencial at y aceleración normal an. Estas aceleraciones se determinan aplicando la segunda ley de Newton a un movimiento circular de radio R

mg·sinθ=mat

La segunda ecuación, junto al principio de conservación de la energía, nos permite calcular la reacción del plano N, en la posición θ

N=mgcosθm v 2 R N=3mgcosθ2mg

La partícula deja de tener contacto con la cúpula cuando la reacción N se anule. Para el ángulo θc tal que

cosθ= 2 3

Aproximadamente, 48º medidos desde la vertical. Como vemos el ángulo límite es independiente del radio de la cúpula, de la masa de la partícula y de la aceleración de la gravedad g.

La velocidad de la partícula cuando alcanza en esta posición es

v 0 2 = 2 3 gR

La primera ecuación nos permite calcular la posición angular θ en función del tiempo t.

R d 2 θ d t 2 =gsinθ

Que se resuelve por procedimientos numéricos. Si se establece las condiciones iniciales t=0, θ=0, dθ/dt=0. La partícula permanece todo el tiempo en esta posición de equilibrio inestable. Las condiciones iniciales serán t=0, θ=θ0, dθ/dt=ω0.Donde θ0 es un pequeño ángulo y ω0 es la velocidad angular de la partícula correspondiente a este ángulo, que se calcula aplicando el principio de conservación de la energía. En la figura, se representa la posición θ de la partícula en función del tiempo t

R=1; %radio
x0=[0.01,sqrt(9.8/R)*2*sin(0.01/2)];
f=@(t,x) [x(2);9.8*sin(x(1))/R]; 
opts=odeset('events',@(t,x) stop_cupula(t,x,acos(2/3)));
tspan=[0,100]; %hasta un tiempo de 3
[t,x]=ode45(f,tspan,x0,opts);
tFin=t(end);

plot(t,x(:,1)*180/pi)
grid on
xlabel('t')
ylabel('\theta')
title('Cúpula semiesférica sin rozamiento')

El procedimiento de cálculo se termina cuando se alcanza la posición final θ=arccos(2/3), para ello se define la función

function [value,isterminal,direction]=stop_cupula(~,x,xFin)
    value=x(1)-xFin; 
    isterminal=1;
    direction=0; 
end

Movimiento bajo la aceleración constante de la gravedad

Una vez que la partícula deja de tener contacto con la cúpula, se mueve bajo la acción de su propio peso, es decir, describe una trayectoria parabólica desde el punto de coordenadas

x0=R·sinθ
y0=R·cosθ

Con velocidad inicial

v 0x = v 0 cosθ v 0y = v 0 sinθ

Las ecuaciones del movimiento son

x= x 0 + v 0 cosθt y= y 0 v 0 sinθt 1 2 g t 2

El punto de impacto sobre el suelo se calcula poniendo y=0 en la segunda ecuación, despejando el tiempo t y sustituyéndolo en la primera.

Ejemplo:

Sea el radio de la cúpula es R=1 m. En el momento en el que la partícula deja de tener contacto con la cúpula N=0, su posición angular es cosθ= 2/3 y su velocidad es,

v 0 = 2 3 9.8·1 =2.56m/s

Cuando llega al suelo y=0. Se resuelve la ecuación de segundo grado t=0.222 s

Se calcula el alcance medido desde el centro de la cúpula x=1.12 m

R=1;
angulo=acos(2/3);
v0=sqrt(2*R*9.8/3);
%tiempo que tarda en llegar al suelo, y=0
t=(sqrt((v0*sin(angulo))^2+2*9.8*R*cos(angulo))-v0*sin(angulo))/9.8;
%alcance
x=R*sin(angulo)+v0*cos(angulo)*t
x =1.1246

Actividades

El radio R de la cúpula se ha fijado en 1 m

Se pulsa el botón Nuevo.

El círculo situado en la parte superior izquierda representa la energía total de la partícula, la porción de color rojo representa la energía cinética, y la porción azul, la energía potencial. Observamos que la energía potencial se va transformando en energía cinética, pero la suma de los valores de ambas clases de energía se mantiene constante a lo largo de la trayectoria de la partícula.

Mediante una línea roja a trazos se señala el ángulo límite arccos(2/3), para el cual la reacción N=0