Dos partículas unidas por una cuerda

Coordenadas polares

La posición de una partícula es (x, y) en coordenadas rectangulares y (r, θ) en coordenadas polares. La relación es

x=r·cosθ, y=r·sinθ

Expresamos la velocidad de la partícula en coordenadas polares

v= dr dt = r ^ dr dt +r d r ^ dt

Calculamos las componentes rectangulares de los vectores unitarios r y θ.

r ^ = i ^ cosθ+ j ^ sinθ θ ^ = i ^ sinθ+ j ^ cosθ

vemos que

d r ^ dt =( i ^ sinθ+ j ^ cosθ ) dθ dt = θ ^ dθ dt d θ ^ dt =( i ^ cosθ j ^ sinθ ) dθ dt = r ^ dθ dt

Las expresión del vector velocidad en coordenadas polares es

v= dr dt r ^ +r dθ dt θ ^

Las expresión del vector aceleración es:

a= dv dt = d 2 r d t 2 r ^ + dr dt d r ^ dt +( r d 2 θ d t 2 + dr dt dθ dt ) θ ^ +r dθ dt d θ ^ dt = d 2 r d t 2 r ^ + dr dt dθ dt θ ^ +( r d 2 θ d t 2 + dr dt dθ dt ) θ ^ r ( dθ dt ) 2 r ^ = ( d 2 r d t 2 r ( dθ dt ) 2 ) r ^ +( r d 2 θ d t 2 +2 dr dt dθ dt ) θ ^

Ecuaciones del movimiento de las partículas

Supongamos dos cuerpos de masas m1 y m2 unidos por una cuerda inextensible y de masa despreciable que pasa a través de un pequeño agujero hecho en un tablero. El primer cuerpo se puede mover en el plano horizontal, mientras que el segundo se puede mover hacia arriba o hacia abajo a lo largo del eje vertical Z, tal como se muestra en la figura. Consideramos únicamente la situación ideal sin rozamiento alguno.

Establecemos un sistema de referencia, de modo que el origen está en el orificio hecho en la tabla, el eje Z apunta hacia abajo, la posición del cuerpo m2 es z y la posición del cuerpo m1 viene determinado por su distancia radial r al orificio y el ángulo θ que hace con el eje X, tal como se muestra en la figura. Si la longitud de la cuerda es l, se cumple que r+z=l.

La ecuación del movimiento de m2 es

m 2 d 2 z d t 2 = m 2 gT

o bien,

m 2 d 2 r d t 2 = m 2 gT

En el primer apartado, hemos obtenido las expresiones de las componentes del vector aceleración en coordenadas polares. Formulamos la segunda ley de Newton para cada una de las componentes, teniendo en cuenta que la única fuerza que actúa sobre la partícula (la tensión de la cuerda) tiene dirección radial y sentido hacia el origen.

m 1 ( d 2 r d t 2 r ( dθ dt ) 2 )=T m 1 ( r d 2 θ d t 2 +2 dθ dt dr dt )=0

La segunda ecuación nos indica que el momento angular L (expresado en coordenadas polares es constante)

d dt ( r 2 dθ dt )=0 r 2 dθ dt = L m 1 =cte

Teniendo en cuenta la constancia del momento angular L y eliminado la tensión T de la cuerda en las ecuaciones del movimiento de las dos partículas, llegamos a la ecuación diferencial que describe el movimiento de la partícula de masa m1 en la dirección radial, r.

m 1 ( d 2 r d t 2 L 2 m 1 2 r 3 )= m 2 g m 2 d 2 r d t 2 ( m 1 + m 2 ) d 2 r d t 2 L 2 m 1 r 3 = m 2 g

Las ecuaciones diferenciales que describen el movimiento de la partícula sobre el tablero en coordenadas polares (r, θ), son por tanto,

d 2 r d t 2 L 2 m 1 ( m 1 + m 2 ) r 3 = m 2 ( m 1 + m 2 ) g dθ dt = L m 1 r 2

Con las condiciones iniciales siguientes, la partícula parte en el instante t=0 de la posición r0, θ=0, con velocidad dr/dt=0, dθ/dt0.

Energía del sistema de partículas

Establecemos el nivel cero de energía potencial en origen z=0. La energía E del sistema de dos partículas es constante

E= 1 2 m 1 ( dr dt ) 2 + 1 2 m 1 r 2 ( dθ dt ) 2 + 1 2 m 2 ( dz dt ) 2 m 2 gz

Teniendo en cuenta la constancia del momento angular L y que z+r=l

E= 1 2 ( m 1 + m 2 ) ( dr dt ) 2 + L 2 2 m 1 r 2 m 2 g( lr )

Derivando esta ecuación con respecto del tiempo y teniendo en cuenta, que la energía E es constante, dE/dt=0, obtenemos la ecuación del movimiento de m1 en la dirección radial.

Casos particulares

Dos situaciones se analizan a nivel elemental

Movimiento en la dirección radial

La situación más sencilla, es aquella en la que m1 en reposo, se libera y se mueve en dirección radial. Las ecuaciones del movimiento en la dirección radial para la partícula de masa m1 y de la partícula de masa m2 a lo largo del eje Z, son

m 1 d 2 r d t 2 =T m 2 d 2 z d t 2 = m 2 gT

Teniendo en cuenta que z+r=l y eliminando la tensión T, vemos que la aceleración de las partículas es constante, m2g/(m1+m2), para m2 y de signo contrario para m1. La distancia radial r de la partícula m1 en función del tiempo t es

r= r 0 1 2 m 2 m 1 + m 2 g t 2

Movimiento circular uniforme

Para que la partícula de masa m1 describa una trayectoria circular de radio r0, bajo la acción de la tensión T=m2g, tiene que llevar una velocidad angular ω0, tal que se cumpla la ecuación de la dinámica del movimiento circular uniforme (fuerza =masa por aceleración normal)

m 2 g= m 1 ω 0 2 r 0 ω 0 = m 2 m 1 g r 0

Por ejemplo, para r0=3, m1=1, m2=0.365. La partícula describe un movimiento circular uniforme si se mueve con velocidad angular constante ω0=1.1030. El periodo o tiempo que tarda en dar una vuelta completa es 2π/ω0=5.6963

Caso general

Para el caso general hay que resolver, el sistema de dos ecuaciones diferenciales por procedimientos numéricos, para obtener la posición de la partícula (r,θ) en cada instante t

m1=1; %masa de la partícula sobre el plano horizontal
m2=0.365; %masa de la partícula que cuelga
w0=0.6;  %velocidad angular inicial en la 
r0=3;		%posición inicial r0, que se mantiene fija
tf=10;	%tiempo final, ajustar

L=m1*r0^2*w0;  %momento angular constante
x0=[r0,0,w0];	%condiciones iniciales
%x(1) es r, x(2) es dr/dt, x(3) es theta
fg=@(t,x) [x(2); L^2/(m1*(m1+m2)*x(1)^3)-m2*10/(m1+m2);
L/(m1*x(1)^2)];
[t,x]=ode45(fg,[0,tf],x0);
plot(x(:,1).*cos(x(:,3)), x(:,1).*sin(x(:,3)))
grid on
xlabel('r')
ylabel('\theta')
axis equal
figure
subplot(2,1,1)
plot(t,x(:,1))
grid on
xlabel('t')
ylabel('r')
subplot(2,1,2)
plot(t,x(:,3))
grid on
xlabel('t')
ylabel('\theta')

Se representa la trayectoria de la partícula sobre el plano horizontal

Se representa sus coordenadas (r, θ) en función del tiempo t

Estos ejemplos reproducen algunas de las figuras del artículo citado en las referencias

m1 m2 dθ/dt
1.0 0.365 0.6
0.3 1.45 1.0
0.333 1.4 0.8
0.157 1.0 1.15
0.335 1.26 0.8
0.3 3.815 1.0

Actividades

Se introduce

Se pulsa el botón titulado Nuevo, se dibuja la trayectoria de la partícula sobre el plano horizontal.

En la parte inferior izquierda, se muestra en color rojo el tanto por ciento de error, la cantidad

| E E 0 | E 0 100

Donde E0 es la energía inicial y E es la energía en el instante t.

Cuando este cociente es mayor que la unidad, 1% se considera que el programa interactivo no resuelve correctamente el sistema de ecuaciones diferenciales y se detiene.

Se pulsa el botón pausa ||, para parar el movimiento y el botón paso a paso >| para acercarnos al momento en el que el móvil completa una trayectoria y vuelve a repetir el movimiento. De este modo, medimos el periodo. En la parte superior izquierda, se muestra el tiempo y en los ejes, se marca las coordenadas del móvil.


Referencias

William Sacks, Alain Mauger. Two balls and a string: from ordered motion to chaos. Eur. J. Phys. 34 (2013) 1487-1506