Movimiento de una partícula atada a una cuerda que se enrolla en un cilindro horizontal.

Se ata una partícula de masa m al extremo de una cuerda de longitud L, inextensible y de masa despreciable. El otro extremo se fija a la parte superior de un cilindro de radio R. La cuerda se mantiene inicialmente tensa en posición horizontal tal como se muestra en la figura.

Se suelta y la cuerda se enrolla completamente en el cilindro o bien, la cuerda se enrolla parcialmente hasta que su tensión se hace cero. En el segundo caso, supondremos que la partícula se desprende de la cuerda describiendo una trayectoria parabólica.

Cuando la partícula cae, en el instante t la parte de la cuerda que está enrollada en el cilindro es la longitud del arco de una circunferencia de radio R. Describiremos el movimiento de la partícula en términos del ángulo φ.

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 ) θ ^

Para describir el movimiento de la partícula adoptamos el sistema de referencia que se muestra en la figura

Ecuaciones de la dinámica

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

Las ecuaciones del movimiento en la dirección radial r y en la dirección θ son, respectivamente

m( d 2 r d t 2 r ( dθ dt ) 2 )=mgsinθTcosδ m( r d 2 θ d t 2 +2 dr dt dθ dt )=mgcosθTsinδ

Donde sinδ=R/r y cosδ=(L-Rφ)/r

Despejamos la tensión T de la segunda ecuación y la sustituimos en la primera

T m = r R { gcosθ( r d 2 θ d t 2 +2 dr dt dθ dt ) } ( d 2 r d t 2 r ( dθ dt ) 2 )=gsinθ r R { gcosθ( r d 2 θ d t 2 +2 dr dt dθ dt ) } (LRϕ) r

Esta es la ecuación diferencial del movimiento de la partícula en términos de r y θ  y el ángulo φ que como veremos se simplifica notablemente al expresarla únicamente en términos de la variable φ.

El arduo trabajo algebraico consiste en expresar r y θ y sus derivadas en términos del ángulo φ y sus derivadas.

Como vemos en la figura

{ rsinθ=(LRϕ)sinϕRcosϕ rcosθ=(LRϕ)cosϕ+Rsinϕ tanθ= (LRϕ)tanϕR (LRϕ)+Rtanϕ r 2 = ( LRϕ ) 2 + R 2

Expresamos las derivadas con respecto del tiempo t de r y θ en términos de φ y sus derivadas.

r dr dt =R(LRϕ) dϕ dt r 2 dθ dt = ( LRϕ ) 2 dϕ dt

A continuación, las derivadas segundas

r d 2 r d t 2 =R(LRϕ) d 2 ϕ d t 2 + R 4 r 2 ( dϕ dt ) 2 r 2 d 2 θ d t 2 = ( LRϕ ) 2 d 2 ϕ d t 2 2 R 3 r 2 ( LRϕ ) ( dϕ dt ) 2

Las componentes de la aceleración

( d 2 r d t 2 r ( dθ dt ) 2 )= R r (LRϕ) d 2 ϕ d t 2 + 1 r ( R 2 ( LRϕ ) 2 ) ( dϕ dt ) 2 ( r d 2 θ d t 2 +2 dr dt dθ dt )= ( LRϕ ) 2 r d 2 ϕ d t 2 2R r (LRϕ) ( dϕ dt ) 2

Finalmente, expresamos la tensión de la cuerda y la ecuación del movimiento en términos del ángulo φ y sus derivadas.

T m =gsinϕ+(LRϕ) ( dϕ dt ) 2 d 2 ϕ d t 2 = R (LRϕ) ( g R cosϕ+ ( dϕ dt ) 2 )

Resolvemos la ecuación diferencial por procedimientos numéricos con las siguientes condiciones iniciales: en el instante t=0, φ=0, dφ/dt=0.

Una vez que conocemos el ángulo φ en función del tiempo, las coordenadas polares r y θ de la partícula se calculan del siguiente modo.

θ=arctan( (LRϕ)sinϕRcosϕ (LRϕ)cosϕ+Rsinϕ )r= ( LRϕ ) 2 + R 2

A continuación, la abscisa x y la ordenada y de la partícula

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

Conservación de la energía

Ponemos el nivel cero de energía potencial en el origen.

La energía inicial de la partícula es Ei=mgR

La energía final de la partícula es

E f = 1 2 m( v r 2 + v θ 2 )mgrsinθ E f = 1 2 m( ( dr dt ) 2 + r 2 ( dθ dt ) 2 )mgrsinθ

Expresamos r, θ y sus derivadas en términos de φ y sus derivadas

gR= ( LRϕ ) 2 2 ( dϕ dt ) 2 g( (LRϕ)sinϕRcosϕ ) ( dϕ dt ) 2 = 2g( R+(LRϕ)sinϕRcosϕ ) ( LRϕ ) 2

La tensión T de la cuerda es una función del ángulo φ.

T m = g (LRϕ) ( 3(LRϕ)sinϕ+2R2Rcosϕ )

La tensión de la cuerda se anula cuando

3(LRϕ)sinϕ+2R2Rcosϕ=0 3 2 ( ϕ L R )=tan( ϕ 2 )

Se resuelve por procedimientos numéricos esta ecuación trascendente.

La cuerda alcanzará la tensión nula siempre que no se enrolle completamente, es decir, Rφ<L. El primer término de la igualdad ha de ser negativo, por lo que se tiene que igualar a la rama negativa de tan(φ/2), es decir, el ángulo 180<φ<360, tal como se aprecia en la figura.

En la figura representamos la recta y=(3/2)(φ-L/R) y la función y=tan(φ/2) para L=6 m y R=1 m.

Denominamos longitud crítica de la cuerda Lc a aquella que se enrolla completamente en el cilindro y a la vez, alcanza la tensión cero en el punto de contacto de la partícula.

Como vemos en la figura, la recta y=(3/2)(φ-L/R) y la función y=tan(φ/2) se cortan en el punto de tangencia. En el punto de corte la pendiente de la función es 3/2.

3 2 = 1 2 1 cos 2 (ϕ/2) ϕ c =250.5º=4.372rad

Introducimos el valor del ángulo φ en la ecuación trascendente y despejamos el cociente L/R.

L c R = ϕ c 2 3 tan( ϕ c 2 ) L c R =5.315

Tiro parabólico

Cuando la cuerda alcanza la tensión cero, supondremos que se corta y la partícula se mueve bajo la acción del peso, describiendo un tiro parabólico.

Se calcula la posición inicial x0, y0 de la partícula en el momento en el que la tensión es cero.

x0=r·cosθ
y0=-r·sinθ

Se calculan las componentes de la velocidad inicial

v 0x = v r cosθ v θ sinθ=( dr dt )cosθr( dθ dt )sinθ v 0x = v r sinθcosθ=( dr dt )sinθr( dθ dt )cosθ

Expresamos v0x y v0y en términos del ángulo φ y sus derivadas

v 0x =(LRϕ) dϕ dt sinϕ v 0y =(LRϕ) dϕ dt cosϕ

Las ecuaciones del movimiento de la partícula bajo la acción de su propio peso son:

x= x 0 + v 0x t y= y 0 + v 0y t+ 1 2 (g) t 2

Actividades

Se introduce

Observamos el movimiento de la partícula:

stokesApplet aparecerá en un explorador compatible con JDK 1.1.