Modos de vibración de una cuerda con cuentas

Sea una cuerda de longitud 2L y densidad lineal μ0=m/L, sujeta por ambos extremos, cuya tensión es T en la que se colocado una masa puntual (una cuenta) M en su punto medio.

Fuerza sobre un elemento diferencial de cuerda

Consideremos una porción de cuerda comprendida entre x y x+dx

La componente Y de las fuerzas, dFy, que actúan sobre dicho elemento diferencial es

d F y =T(x+dx)·sin(θ+dθ)T(x)sinθ ( T(x)+ dT dx dx )·(sinθ+cosθ·dθ)T(x)sinθ= d dz ( sinθ·T(x) )·dx d dz ( tanθ·T(x) )·dx= d dx ( dy dx ·T(x) )·dx

Para ángulos pequeños, hacemos la aproximación sinθ≈tanθ

Esta componente dFy de fuerza sobre el elemento diferencial de masa μ(x)dx produce una aceleración a lo largo del eje Y. Teniendo en cuenta, y es una función continua de x y t escribimos, la segunda ley de Newton, fuerza=masa×aceleración

x ( y x T(x) )dx=( μ·dx ) 2 y t 2 x ( y x T(x) )=μ 2 y t 2

Donde μ(x) es la densidad de una cuerda no homogénea y T(x) es la tensión en cada punto de la cuerda

y(x,t) es una función continua de x y su derivada dy/dx también lo es.

Fuerza sobre una masa puntual situada en la cuerda

Consideremos una masa puntual Mi sujeta a la cuerda en la posición xi

Como vemos en la figura, en un instante t, la función y(x) es continua pero no su derivada dy/dx en xi

La componente Y de las fuerzas, Fy, que actúan sobre dicha partícula es

F y =Tsin θ 2 T'sin θ 1 Ttan θ 2 T'tan θ 1 =( ( T y dx ) x i + ( T y dx ) x i )

T es la fuerza que ejerce la parte izquierda de la cuerda sobre la partícula y T' es la fuerza que ejerce la parte derecha de la cuerda sobre la partícula. θ1 y θ2 son las direcciones de dichas fuerzas, tangentes a la cuerda en el punto xi

La segunda ley de Newton a lo largo del eje Y se escribe

( ( T y x ) x i + ( T y x ) x i )= M i 2 y t 2

Modos de vibración de una cuerda con una cuenta

Escribimos la solución de la ecuación diferencial en derivadas parciales como producto, y(x,t)=Φ(x)·cos(ωt)

d dx ( dΦ dx T(x) )+μ ω 2 Φ=0 ( ( T dΦ dx ) L + ( T dΦ dx ) L )+ M ω 2 Φ=0

Consideramos que la tensión T(x)=T de la cuerda es constante, su densidad μ0=m/L es constante y la posición de la cuenta de masa M es xi=L.

d 2 Φ d x 2 + k 2 Φ=0 ω=k T μ 0

El término, bajo la raíz cuadrada es la velocidad de propagación de las ondas transversales en una cuerda. La solución de esta ecuación diferencial tiene la forma conocida, Asin(kx)+Bcos(kx)

La masa puntual M situada en la posición x=L, divide a la cuerda en dos partes. Probamos la solución

Φ(x)={ Asin(kx)+Bcos(kx)0x<L Csin( k(xL) )+Dcos( k(xL) )Lx<2L

En la figura, se muestran los modos de vibración de una cuerda de longitud 2L. Los modos de vibración de orden par n=2,4,6,... presentan un nodo en x=L

Por tanto, los valores de knL se obtienen de dos ecuaciones

Representamos la función f(x)= M m xsin(x)2cos( x ) , para el valor M/m=1. Se sugiere cambiar este valor para apreciar el efecto de la masa M de la partícula

L=1; %longitud
M_m=1; %cociente masas

f=@(x) M_m*x.*sin(x)-2*cos(x);
fplot(f,[0,10])
grid on
xlabel('kL')
ylabel('y');
title('Cuerda con una cuenta')

Las raíces son los puntos de corte con el eje X. Utilizamos la función raices definida en el apartado titulado 'Raíces múltiples' de la página Funciones MATLAB para calcular las raíces de una ecuación, para calcularlas. Asignamos al cociente M/m=1. Donde M es la masa de la partícula y μ0=m/L es la densidad lineal de la cuerda

M_m=1; %cociente M/m
f=@(x) M_m*x.*sin(x)-2*cos(x);
x=linspace(0,20,20);
k=raices(f,x); % raíces
disp(k);
    1.0769    3.6436    6.5783    9.6296   12.7223   15.8336   18.9547

Modos de vibración

Representamos los tres primeros modos de vibración para M/m=1

L=1; %longitud
M_m=1; %relación entre masas
f=@(x) M_m*x.*sin(x)-2*cos(x);

subplot(3,1,1)
k=fzero(f,1);
area=L-sin(2*k*L)/(2*k);
hold on
fplot(@(x) sin(k*x)/sqrt(area),[0,L])
fplot(@(x) sin(k*(2*L-x))/sqrt(area),[L,2*L])
plot(L,sin(k*L)/sqrt(area),'ro','markersize',4,'markeredgecolor','r'
,'markerfacecolor','r')
hold off
grid on
xlabel('x')
ylabel('\Phi')
title(sprintf('k=%1.2f',k))

k=pi/L;
subplot(3,1,2)
hold on
fplot(@(x) sin(k*x),[0,2*L])
plot(L,0,'ro','markersize',4,'markeredgecolor','r','markerfacecolor','r')
hold off
grid on
xlabel('x')
ylabel('\Phi')
title(sprintf('k=%1.2f',pi/L))

k=fzero(f,4);
area=L-sin(2*k*L)/(2*k);
subplot(3,1,3)
hold on
fplot(@(x) sin(k*x)/sqrt(area),[0,L])
fplot(@(x) sin(k*(2*L-x))/sqrt(area),[L,2*L])
plot(L,sin(k*L)/sqrt(area),'ro','markersize',4,'markeredgecolor','r',
'markerfacecolor','r')
hold off
grid on
xlabel('x')
ylabel('\Phi')
title(sprintf('k=%1.2f',k))

Observamos que para los modos impares, la función Φ(x) es continua y también su derivada. La partícula pernanece en reposo, no vibra, en un nodo. Para los modos pares, la función Φ(x) es continua pero no lo es su derivada, la partícula oscila.

Actividades

Se introduce

Se pulsa el botón titulado Nuevo

Se pulsa el botón titulado >>, para observar el siguiente modo de vibración.

Se pulsa el botón titulado <<, para observar el modo de vibración anterior

En la parte superior, se muestra el modo de vibración y su frecuencia f en Hz.

Condiciones iniciales

La solución general y(x,t) que satisface la ecuación diferencial del movimiento ondulatorio

T μ 0 2 y x 2 = 2 y t 2

con las condiciones de contorno y(0,t)=0 e y(2L,t)=0 es

y(x,t)= n=1 Φ n (x)( A n cos( ω n t )+ B n sin( ω n t ) ) ω n = k n T μ 0

Esta ecuación describe todos los posibles modos de vibración de la cuerda. La vibración particular que experimenta la cuerda está únicamente determinada por las condiciones iniciales, que a su vez determinan los valores de las constantes An y Bn.

y(x,0)= y 0 (x) n=1 A n Φ n (x) = y 0 (x) y(x,t) t | t=0 = y ˙ 0 (x) n=1 ω n B n Φ n (x) = y ˙ 0 (x)

Supongamos que la cuerda se estira desde su punto medio donde se encuentra la cuenta, hasta adoptar la forma inicial la siguiente.

y(x,0)={ h L x0x<L h L x+2hLx<2L

Como velocidad inicial de la cuerda es nula, los coeficientes Bn=0.

Las relaciones de ortogonalidad de Φn(x) para una cuerda no homogénea de densidad μ(x)

0 2L μ(x) Φ m (x) Φ n (x)dx=0mn

Para una cuerda homogénea de densidad μ0 con una cuenta de masa M situada en x=L se escribe

0 L μ 0 Φ m (x) Φ n (x)dx+ L L+ μ(x) Φ m (x) Φ n (x)dx+ L+ 2L μ 0 Φ m (x) Φ n (x)dx= μ 0 0 L Φ m (x) Φ n (x)dx+M Φ n (L) Φ m (L)+ μ 0 L+ 2L Φ m (x) Φ n (x)dx=0 μ 0 0 2L Φ m (x) Φ n (x)dx+M Φ n (L) Φ m (L)=0mn


Nota: Formalmente, la densidad de la cuerda uniforme con una cuenta en x=L se escribe

μ(x)= μ 0 +Mδ(xL) δ(xL)f(x)dx=f(L)

Véase la función delta de Dirac δ(x) al final de la página titulada Transformada de Fourier


Utilizamos la relación de ortogonalidad de las funciones Φn(x) y Φm(x) para calcular los coeficientes An

0 2L μ(x) y 0 (x) Φ n (x)dx= m=1 A m 0 2L μ(x) Φ m (x) Φ n (x)dx 0 2L μ(x) y 0 (x) Φ n (x)dx= A n 0 2L μ(x) Φ n 2 (x)dx μ 0 0 2L y 0 (x) Φ n (x)dx +M y 0 (L) Φ n (L)= A n { μ 0 0 2L Φ n 2 (x)dx+M Φ n 2 (L) }

La forma de la cuerda es en el instante t es

y(0x<L,t)= n=1 2h ( k n L ) 2 ( 1+ M 2m sin( k n L ) ) sin( k n x)cos( ω n t ) y(Lx<2L,t)= n=1 2h ( k n L ) 2 ( 1+ M 2m sin( k n L ) ) sin( k n (x2L) )cos( ω n t ) ω n = k n T μ 0

knL son las raíces de la ecuación transcendente.

Elaboramos un script, para representar la forma de la cuerda con una cuenta en su punto medio, para varios instantes. La longitud de la cuerda 2L= 2, la deformación máxima de la cuerda es h=0.1 y la relación de masas M/m=1. (M es la masa de la partícula y μ0=m/L la densidad lineal de la cuerda). Se sugiere cambiar este parámetro y observar sus efectos por ejemplo, para M/m=0.5 y M/m=2.

El periodo o tiempo que tarda en viajar una perturbación desde un extremo a otro de la cuerda y regresar al punto de partida es τ=4L/v. Donde v=1 es la velocidad de propagación. Representamos la forma de la cuerda en los instantes t/τ=0 (forma inicial), 0.5, 1, 2, 3, 3.5

M_m=1;%relación entre masas M/m
L=1; %longitud de la cuerda
h=0.1; %deformación de la cuerda

f=@(x) M_m*x.*sin(x)-2*cos(x);
x=linspace(0.1,60,50);
k=raices(f,x); %20 raíces
An=2*h*sin(k*L)./(((k*L).^2).*(1+M_m*sin(k*L).^2/2)); %coeficientes

hold on
x=linspace(0,2*L,200);
%forma inicial
y=(h*x/L).*(heaviside(x)-heaviside(x-L))+(2*h*(1-x/(2*L))).*
(heaviside(x-L)-heaviside(x-2*L));
plot(x,y, 'r','displayName',num2str(0))

for t=[0.5,1,2,3,3.5]
    yy=zeros(1,length(x)); 
    for n=1:length(k)
        yy=yy+An(n)*(sin(k(n)*(x<=L).*x)-sin(k(n)*(x>L).*(x-2*L)))*cos(k(n)*t/L);
    end
    plot(x,yy, 'displayName',num2str(t))
end
hold off
grid on
legend('-DynamicLegend','location','northeast')
xlabel('x')
ylabel('y(x,t)')
title('Cuerda con una cuenta')

Modificamos el script para representar la forma de la cuerda sin cuenta M/m=0, en los instantes t/τ=0 (forma inicial), 0.5, 1, 1.5, 2. Comprobamos que vuelve a la forma inicial en el instante t/τ=4

Actividades

Se introduce

Se pulsa el botón titulado Nuevo

Dos partículas sujetas a la cuerda

Consideremos una cuerda homogénea de longitud 3L y densidad uniforme m/L con dos masas puntuales M situadas en las posiciones x=L y x=2L tal como se muestra en la figura.

Las masas puntuales M situadas en la posiciones x=L y x=2L, dividen a la cuerda en tres partes. Probamos la solución

Φ(x)={ Asin(kx)+Bcos(kx)0x<L Csin( k(xL) )+Dcos( k(xL) )Lx<2L Esin( k(x2L) )+Fcos( k(x2L) )2Lx<3L

Tenemos un sistema de cinco ecuaciones homogéneas con cinco incógnitas

{ Esin(kL)+Fcos(kL)=0 Asin(kL)D=0 Acos(kL)+C+ M m (kL)D=0 Csin(kL)+Dcos(kL)F=0 Ccos(kL)+Dsin(kL)+E+ M m (kL)F=0 ( 0 0 0 sin(kL) cos(kL) sin(kL) 0 1 0 0 cos(kL) 1 M m (kL) 0 0 0 sin(kL) cos(kL) 0 1 0 cos(kL) sin(kL) 1 M m (kL) )( A C D E F )=0

El determinante de la matriz de los coeficientes deberá ser cero

Llamamos x a kL y q al cociente M/m

>> syms q x;
>> M=[0,0,0,sin(x),cos(x); sin(x),0,-1,0,0;-cos(x),1,q*x,0,0;
0,sin(x),cos(x),0,-1;0,-cos(x),sin(x),1,q*x];
>> eq=det(M);
>> factor(eq)
ans =[ sin(x), 3*cos(x)^2 - sin(x)^2 + q^2*x^2*sin(x)^2 - 4*q*x*cos(x)*sin(x)]

La segunda expresión entre corchetes, se escribe, 4cos(x)2 - 1 + q2x2sin(x)2 - 4qxcos(x)sin(x). Comprobamos que es el producto de (qxsin(x)-2cos(x)+1)(qxsin(x)-2cos(x)-1).

>> eq=4*cos(x)^2 - 1 + q^2*x^2*sin(x)^2 - 4*q*x*cos(x)*sin(x);
>> factor(eq)
ans =[ q*x*sin(x) - 2*cos(x) + 1, q*x*sin(x) - 2*cos(x) - 1]

El determinante es el producto de tres factores

sin(kL)( M m (kL)sin(kL)2cos(kL)+1 )( M m (kL)sin(kL)2cos(kL)1 )=0

Los valores de knL se obtienen de tres ecuaciones

Calculamos cuatro raíces knL de cada una de las tres ecuaciones. Las ordenamos mediante la función sort de MATLAB

L=1;
M_m=1; %cociente masas
%primera función
f=@(x) M_m*x.*sin(x)-2*cos(x)+1;
x=linspace(0,10,50);
r1=raices(f,x);
disp(r1);

%segunda función
g=@(x) M_m*x.*sin(x)-2*cos(x)-1;
x=linspace(0,10,50);
r2=raices(g,x);
disp(r2);

r3=(1:4)*pi/L;
disp(r3)

r=[r1,r2,r3];
sort(r)
    0.7314    3.8528    6.4356    9.7284
    1.3876    3.4159    6.7158    9.5288
    3.1416    6.2832    9.4248   12.5664
ans =    0.7314    1.3876    3.1416    3.4159    3.8528    6.2832    6.4356
    6.7158    9.4248    9.5288    9.7284   12.5664

Despejamos los coeficientes C, D, E y F del sistema de ecuaciones

{ D=Asin(kL) C=Acos(kL) M m (kL)D F=Csin(kL)+Dcos(kL) E=Ccos(kL)Dsin(kL) M m (kL)F

Cuando knL=nπ, sin(kL)=0, los coeficientes, D=0, C=Acos(kL), F=0 y E=Acos2(kL)=A

Φ(x)={ Asin(kx)0x<L Csin( k(xL) )Lx<2L Esin( k(x2L) )2Lx<3L Φ(x)={ Asin(kx)0x<L Acos(kL)( sin(kx)cos(kL)cos(kx)sin(kL) )Lx<2L A( sin(kx)cos(2kL)cos(kx)sin(2kL) )2Lx<3L Φ(x)={ Asin(kx)0x<3L

Modos de vibración

Φ(x)={ Asin(kx)0x<L Csin( k(xL) )+Dcos( k(xL) )Lx<2L Esin( k(x2L) )+Fcos( k(x2L) )2Lx<3L

Los coeficientes C, D, E, F se expresan en términos de A

{ D=Asin(kL) C=Acos(kL) M m (kL)D F=Csin(kL)+Dcos(kL) E=Ccos(kL)Dsin(kL) M m (kL)F

Determinamos el valor del coeficiente A de modo que

0 3L Φ n 2 (x)·dx=1

El resultado de calcular las integrales es

A 2 0 L sin 2 (kx)dx + L 2L ( Csin( k(xL) )+Dcos( k(xL) ) ) 2 dx+ L 2L ( Esin( k(x2L) )+Fcos( k(x2L) ) ) 2 dx=1 A 2 0 L sin 2 (kx)dx + C 2 L 2L sin 2 ( k(xL) )dx + D 2 L 2L cos 2 ( k(xL) )dx +CD L 2L sin( 2k(xL) )dx + + E 2 L 2L sin 2 ( k(x2L) )dx + F 2 L 2L cos 2 ( k(x2L) )dx +EF L 2L sin( 2k(xL) )dx = 1 2 ( A 2 + C 2 + D 2 + E 2 + F 2 )+ 1 4k sin( 2kL )( A 2 C 2 + D 2 E 2 + F 2 )+ 1 2k ( 1cos( 2kL ) )(CD+EF)

Representamos los tres primeros modos de vibración cuyos valores k1=0.7314, k2=1.3876, k3=3.1416, hemos calculado

L=1; %longitud
M_m=1; %relación entre masas
A=1;

i=1;
for k=[0.7314, 1.3876, 3.1416]
    subplot(3,1,i)
    D=A*sin(k*L);
    C=A*cos(k*L)-M_m*(k*L)*D;
    F=C*sin(k*L)+D*cos(k*L);
    E=C*cos(k*L)-D*sin(k*L)-M_m*(k*L)*F;
    area=(A^2+C^2+D^2+E^2+F^2)/2+(-A^2-C^2+D^2-E^2+F^2)*sin(2*k*L)/(4*k)+
(C*D+E*F)*(1-cos(2*k*L))/(2*k);
    f1=@(x) sin(k*x)/sqrt(area);
    f2=@(x) (C*sin(k*(x-L))+D*cos(k*(x-L)))/sqrt(area);
    f3=@(x) (E*sin(k*(x-2*L))+ F*cos(k*(x-2*L)))/sqrt(area);
    hold on
    fplot(f1,[0,L])
    fplot(f2,[L,2*L])
    fplot(f3,[2*L,3*L])
    plot(L,f1(L),'ro','markersize',4,'markeredgecolor','r','markerfacecolor','r')
    plot(2*L,f2(2*L),'ro','markersize',4,'markeredgecolor','r','markerfacecolor','r')
    hold off
    grid on
    xlabel('x')
    ylabel('\Phi')
    title(sprintf('k=%1.2f',k))
    i=i+1;
end

La tercera figura representa el tercer modo de vibración de una cuerda de longitud 3L. Para los modos múltiplos de tres, la función Φ(x) es continua y también su derivada. La partícula pernanece en reposo, no vibra, en un nodo. Para los otros modos, la función Φ(x) es continua pero no lo es su derivada, la partícula oscila.

Actividades

Se introduce

Se pulsa el botón titulado Nuevo

Se pulsa el botón titulado >>, para observar el siguiente modo de vibración.

Se pulsa el botón titulado <<, para observar el modo de vibración anterior

En la parte superior, se muestra el modo de vibración y su frecuencia f en Hz.

Referencias

Gómez B. J. Repetto C. E., Stia C. R., Welti R. Oscillations of a string with concentrated masses. Eur. J. Phys. 28 (2007) pp. 961-975