Cuerda que gira alrededor de un eje
Consideremos un elemento diferencial ds de cuerda. Dibujamos las fuerzas sobre dicho elemento.
- La fuerza aplicada dF(x)
- La tensión T(x), que ejerce la parte izquierda de la cuerda sobre dicho elemento
- La tensión T(x+dx), que ejerce la parte derecha de la cuerda sobre dicho elemento
Ecuaciones de equilibrio
Equilibrio a lo largo del eje X
Equlibrio a lo largo del eje Y
T(x+dx)·cosα(x+dx)-T(x)·cosα(x)=0
T(x+dx)·sinα(x+dx)-T(x)·cosα(x)-dF(x)=0
La primera ecuación es equivalente a d(T(x)·cosα(x))=0, o bien,
siendo T0 es la tensión de la cuerda para α=0
Dividiendo la segunda entre la primera
La pendiente de la recta tangente a la curva y=y(x) en x es
La ecuaciones de equilibrio se escriben
Dependiendo de la expresión de la fuerza externa dF/dx, obtenemos distintas ecuaciones diferenciales y así distintas soluciones para la forma de la curva
La catenaria
El peso de un elemento ds de cuerda es dF=λg·ds, siendo λ la densidad lineal y , la longitud del elemento diferencial de cuerda
La ecuación diferencial de la catenaria que se ha estudiado en la página titulada La catenaria es
Forma de la cuerda que gira
Supongamos que los extremos de la cuerda están atados a un eje (el X) que gira con velocidad angular constante ω. En el sistema de referencia en rotación los elementos de la cuerda experimentan una fuerza centrífuga proporcional al radio -y (en la figura y es negativo) de la circunferencia que describen. La fuerza centrífuga es similar al peso, es proporcional a la masa. La aceleración de dicho elemento es -ω2y. En la ecuación diferencial de la catenaria reemplazamos g por -ω2y
Llamamos 1/L2 a la constante de proporcionalidad, que tiene dimensión de longitud.
Resolvemos la ecuación diferencial con las condiciones en los extremos y(0)=0 e y(a)=0.
Solución aproximada
Si la deformación de la cuerda no es muy grande, se cumple que la pendiente dy/dx en todos los puntos de la cuerda es pequeña
Por lo que la solución de la ecuación diferencial es sencilla, recuerda a la de un Movimiento Armónico Simple
ym es la amplitud que tiene que ser pequeña, la frecuencia angular es 1/L y φ la fase inicial. La distancia entre los extremos de la cuerda sujetos al eje es medio periodo, Δx=π/(1/L)=πL, tal como se ve en la figura
L=1/pi; g=@(x) sin(x/L-pi/2); fplot(g,[-0.5,0.5]) grid on xlabel('x') ylabel('y') title('Cuerda que gira')
Como vemos la cuerda está sujeta en los puntos x=-0.5 y x=0.5. Falta multiplicar la función seno por una amplitud ym pequeña, para simular la deformación de la cuerda que gira con velocidad angular ω pequeña.
Solución analítica
Escribimos esta ecuación de la forma equivalente
donde ym es una constante a determinar. Elevamos al cuadrado y ordenamos los términos para que nos salga una ecuación diferencial del tipo
cuya solución es z=sn(u,p), véase Wikipedia, Jacobi elliptic functions o Integrales elípticas
En este problema
con
La solución z=sn(u, p) se transforma en
Como se ha comprobado en la página titulada Integrales elípticas, medio periodo de la función sn(u,p) es Δu=2K(p) lo que corresponde a
Siendo K(p) la integral elíptica completa de primera especie. Δx es la distancia entre los extremos de la cuerda e ym es la amplitud
Sea
L=1/pi; ym=1; p=sqrt(ym^2/(4*L^2+ym^2)); K=ellipke(p^2); Dx=2*K*L/(sqrt(1+ym^2/(4*L^2))); %distancia entre los extremos alfa=sqrt(1+ym^2/(4*L^2))/L; f=@(x) -ellipj(alfa*x+K,p^2); hold on fplot(f,[-Dx/2,Dx/2]) g=@(x) ym*sin(x/L-pi/2); fplot(g,[-pi*L/2,pi*L]) hold off ylim([-1,0]) legend('exacta','aproximada', 'location','best') grid on xlabel('x') ylabel('y') title('Cuerda que gira')
La solución exacta y aproximada difieren
Cuando L=1, la solución exacta y aproximada casi coinciden
Procedimiento numérico
Este problema nos da la oportunidad de resolver una ecuación diferencial con condiciones en los extremos y comparar la solución numérica aplicando el procedimiento
Los pasos para aplicar el procedimiento
Una funcion que define la ecuación diferencial, lo mismo que en el procedimiento
ode45
f=@(~,x) [x(2), -x(1)*sqrt(1+x(2)^2)/L^2];
Una función que especifica las condiciones en los extremos
Para -Δx/2, y=0 (extremo izquierdo de la cuerda), para Δx/2, y=0 (extremo derecho de la cuerda)
g=@(xa,xb) [xa(1), xb(1)];
La aproximación inicial, que va a consistir en la función
xIni=linspace(-Dx/2, Dx/2,100); aprox=@(x) [ym*sin(x/L-pi/2), ym*cos(x/L-pi/2)/L]; solinit = bvpinit(xIni, aprox); %primera aproximación
Representamos la función exacta y la función calculada mediante el procedimiento
L=1/pi; ym=1; p=sqrt(ym^2/(4*L^2+ym^2)); K=ellipke(p^2); Dx=2*K*L/(sqrt(1+ym^2/(4*L^2))); %distancia entre los extremos alfa=sqrt(1+ym^2/(4*L^2))/L; fe=@(x) -ellipj(alfa*x+K,p^2); %exacta %aproximada f=@(~,x) [x(2), -x(1)*sqrt(1+x(2)^2)/L^2]; g=@(xa,xb) [xa(1), xb(1)]; xIni=linspace(-Dx/2, Dx/2,100); aprox=@(x) [ym*sin(x/L-pi/2), ym*cos(x/L-pi/2)/L]; solinit = bvpinit(xIni, aprox); %primera aproximación sol1 = bvp4c(f, g, solinit); y = deval(sol1,xIni); hold on fplot(fe,[-Dx/2, Dx/2]) plot(xIni,y(1,:)); hold off grid on xlabel('x') ylabel('y') title('Cuerda que gira'
Cambiamos el parámetro L=1
Referencias
Arne Nordmark, Hanno Essén. The skipping rope curve. Eur. J. Phys. 28 (2007) 241-247