Vibraciones longitudinales de una varilla

Varilla en posición horizontal

La ecuación diferencial del Movimiento Ondulatorio es

2 ψ t 2 = c 2 2 ψ x 2

Donde c es la velocidad de propagación de las ondas longitudinales en la varilla.

c= Y ρ

Condiciones iniciales y de contorno

Para resolver la ecuación diferencial en derivadas parciales que es de segundo orden en x y en t, se precisa especificar dos condiciones iniciales: la deformación inicial de la varilla y la velocidad inicial de cada uno de sus puntos

t=0{ ψ(x,0)= ψ 0 (x) ψ(x,t) t | t=0 = ψ ˙ 0 (x)

y otras dos condiciones de contorno:

Modos de vibración

Para resolver la ecuación diferencial dividimos la solución en producto de dos funciones una dependiente de x y la otra de t.

ψ(x,t)=Φ(x)·T(t) 1 T(t) d 2 T d t 2 = c 2 Φ(x) d 2 Φ d x 2 = ω 2

Como el miembro izquierdo depende solamente de t y el derecho solamente de x, entonces, igualamos ambos a una constante que denominaremos -ω2

d 2 T d t 2 + ω 2 ·T=0 d 2 Φ d x 2 + ω 2 c 2 Φ=0

Hemos sustituimos una ecuación diferencial en derivadas parciales por dos ecuaciones diferenciales. Ambas ecuaciones tienen soluciones conocidas, ecuación de las oscilaciones libres

T(t)=Acos(ωt)+Bsin(ωt) Φ(x)=Ccos( ω c x )+Dsin( ω c x )

Donde A, B, C y D se determinan a partir de las dos condiciones iniciales y las otras dos de contorno.

Si la varilla está sujeta por el extremo izquierdo x=0, entonces Φ(0)=0 y C=0.

Como el extremo x=l es libre

dΦ dx | x=l =D ω c cos( ω c l )=0 ω n c l= 2n+1 2 π,n=0,1,2,3...

ωn son las frecuencias de los modos normales de vibración de la varilla

Φ n (x)=sin( ω n c x )=sin( 2n+1 2 π x l )

La solución de la ecuación diferencial correspondiente al modo n de vibración es

ψ n (x,t)= Φ n (x) T n (t)= sin( 2n+1 2 π x l ){ A n cos( 2n+1 2 π ct l )+ B n sin( 2n+1 2 π ct l ) }

donde los coeficientes An y Bn se determinarán a partir de las condiciones iniciales

Superposición

La solución general que satisface las condiciones de contorno, es la superposición

ψ(x,t)= n=0 ψ n (x,t) = n=0 Φ n (x)( A n cos( ω n t)+ B n sin( ω n t) ) ψ(x,t)= n=0 sin( 2n+1 2 π x l ){ A n cos( 2n+1 2 π ct l )+ B n sin( 2n+1 2 π ct l ) }

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

ψ(x,0)= ψ 0 (x) n=0 A n sin( 2n+1 2 π x l ) = ψ 0 (x) ψ(x,t) t | t=0 = ψ ˙ 0 (x) n=0 2n+1 2 π c l B n sin( 2n+1 2 π x l ) = ψ ˙ 0 (x)

Teniendo en cuenta el resultado de las integrales

0 l sin 2 ( 2n+1 2 π x l )dx = l 2 0 l sin( 2n+1 2 π x l )sin( 2m+1 2 π x l )dx =0mn

>> syms n x L;
>> assume(n,'integer');
>> int(sin((2*n+1)*pi*x/(2*L))^2,x,0,L)
ans =L/2 - (L*sin(pi*(2*n + 1)))/(2*pi*(2*n + 1))
>> simplify(ans)
ans =L/2

Los valores de An y Bn se determinan de modo análogo a los coeficientes de un desarrollo en serie de Fourier.

A n = 2 l 0 l ψ 0 (x)sin( 2n+1 2 π x l )dx B n = 4 (2n+1)πc 0 l ψ ˙ 0 (x)sin( 2n+1 2 π x l )dx

Perturbación inicial

Supongamos que aplicamos una fuerza al extremo de la varilla, que la estira y luego, la soltamos

ψ(x,0)=kx,0xl ψ(x,t) x | t=0 =0

Calculamos los coeficientes An integrando por partes. Los coeficientes Bn=0, son nulos.

A n = 2 l 0 l kxsin( 2n+1 2 π x l )dx A n = 8kl (1) n ( 2n+1 ) 2 π 2

Aproximamos la deformación inicial de la cuerda Ψ(x,0) con los primeros cinco términos del desarrollo en serie

L=1; %longitud de la varilla
k=0.1; %constante
x=linspace(0,L,200);
y=k*x;
yy=zeros(1,length(x));
for n=0:4 %aproxima con cinco términos
    yy=yy+sin((2*n+1)*pi*x/(2*L))*(-1)^n/(2*n+1)^2;
end
yy=yy*(8*k*L/pi^2);
plot(x,y,x,yy)
xlabel('x')
ylabel('\Psi(x,0)')
title('Deformación inicial de la varilla')
grid on

Solución

La solución de la ecuación diferencial en derivadas parciales con las condiciones iniciales y de contorno especificadas es

ψ(x,t)= 8kl π 2 n=0 (1) n ( 2n+1 ) 2 sin( 2n+1 2 π x l )cos( 2n+1 2 π ct l )

El desplazamiento del extremo x=l de la varilla es

ψ(l,t)= 8kl π 2 n=0 (1) n ( 2n+1 ) 2 sin( 2n+1 2 π )cos( 2n+1 2 π ct l ) = 8kl π 2 n=0 1 ( 2n+1 ) 2 cos( 2n+1 2 π ct l )

Representamos el desplazamiento ψ(x,t) en función del tiempo, de los puntos situados en x=l/4, x=l/2, y x=l (extremo de la varilla), en este último caso, utilizamos la expresión simplificada ψ(l,t)

L=1; %longitud de la varilla
c=1; %velocidad de propagación
k=0.1; %constante
t=linspace(0,4*L/c,200);
y=zeros(1,length(t));
hold on
x=L/4; %posición
for n=0:20 %aproxima con 20 términos
    y=y+(-1)^n*cos((2*n+1)*pi*c*t/(2*L))*sin((2*n+1)*pi*x/(2*L))/(2*n+1)^2;
end
y=y*(8*k*L/pi^2);
plot(t,y,'b')
x=L/2; %posición
y=zeros(1,length(t));
for n=0:20 
    y=y+(-1)^n*cos((2*n+1)*pi*c*t/(2*L))*sin((2*n+1)*pi*x/(2*L))/(2*n+1)^2;
end
y=y*(8*k*L/pi^2);
plot(t,y,'r')
%para x=L; posición
y=zeros(1,length(t));
for n=0:20 
    y=y+cos((2*n+1)*pi*c*t/(2*L))/(2*n+1)^2;
end
y=y*(8*k*L/pi^2);
plot(t,y,'k')
hold off
legend('L/4','L/2','L','Location','north')
xlabel('t')
ylabel('\Psi(x,t)')
title('Deformación de la varilla')
grid on

Los puntos interiores permanecen en reposo un tiempo

El periodo de la oscilación del extremo de la varilla es 4l/c

Actividades

Se introduce

Se pulsa el botón titulado Nuevo

En la parte inferior, aparece la varilla en posición horizontal sujeta por su extremo izquierdo, el extremo derecho se ha desplazado kl y un punto de la varilla situado en x se ha desplazado kx. Observamos el desplazamiento los puntos situados en x=l/2 en color azul

Cuando se pulsa el botón titulado , la varilla deformada se suelta en el instante t=0

Observamos la vibración de los puntos de la varilla. En la parte superior, se representa a otra escala ψ(x,t)

En la parte superior derecha, se proporcionan los datos del tiempo t y de la posición l+ψ(l,t) del extremo derecho de la varilla


Varilla en posición vertical

Consideremos el caso de una varilla de longitud l sin deformar que se mantiene en posición vertical sujeta por su extremo superior. En el instante t=0, se suelta el extremo inferior, la varilla se deforma.

Deducimos la ecuación de propagación de las ondas longitudinales en la varilla cuando está en posición vertical

La parte superior de la varilla ejerce una fuerza F sobre el elemento de varilla de anchura dx, la parte inferior de la barra ejerce una fuerza F' sobre dicho elemento. La diferencia dF=F'-F es

F'F=dF=SY 2 ψ x 2 dx

A la que hay que sumar el peso de dicho elemento, (ρS·dx)g. La segunda ley de Newton afirma que la fuerza sobre dicho elemento es igual al producto de la masa (densidad por volumen) por la aceleración (derivada segunda del desplazamiento)

(ρSdx)g+dF=(ρSdx) 2 ψ t 2

La ecuación diferencial de un movimiento ondulatorio cuando la varilla está en posición vertical es

2 ψ t 2 = c 2 2 ψ x 2 +g,c= Y ρ

Condiciones iniciales y de contorno

Resolvemos esta ecuación diferencial con las condiciones de contorno: El extremo superior x=0, está fijo. El extremo inferior x=l está libre para t≥0

ψ(0,t)=0 ψ(x,t) x | x=l =0

Las conciones iniciales son

En el instante t=0, la varilla está en posición vertical sin deformar

ψ(x,0)=0 ψ(x,t) t | t=0 =0

En la sección anterior, hemos obtenido la solución de la ecuación diferencial homogénea bajo las mismas condiciones de contorno, tiene la forma

ψ(x,t)= n=0 u n (t)sin( 2n+1 2 π x l )

Donde un(t) era una función proporcional a cos(ωnt). Para este caso, vamos a determinar la forma de un(t) sustituyendo la solución en la ecuación diferencial del movimiento ondulatorio cuando la varilla está en posición vertical.

2 ψ t 2 = c 2 2 ψ x 2 +g n=0 ( d 2 u n d t 2 + c 2 ( 2n+1 2l π ) 2 u n )sin( 2n+1 2 π x l )=g

Multiplicamos por sin((2m+1)πx/(2l)) e integramos entre 0 y l. Teniendo en cuenta el resultado de las integrales

0 l sin 2 ( 2n+1 2 π x l )dx = l 2 0 l sin( 2n+1 2 π x l )sin( 2m+1 2 π x l )dx =0mn

Transformamos una ecuación diferencial en derivadas parciales en infinitas ecuaciones diferenciales de segundo orden, fáciles de resolver

n=0 ( d 2 u n d t 2 + c 2 ( 2n+1 2l π ) 2 u n ) 0 l sin( 2n+1 2 π x l )sin( 2m+1 2 π x l )dx =g 0 l sin( 2m+1 2 π x l )dx l 2 { d 2 u n d t 2 + c 2 ( 2n+1 2l π ) 2 u n }=g 2l ( 2n+1 )π d 2 u n d t 2 + c 2 ( 2n+1 2l π ) 2 u n =g 4 ( 2n+1 )π

La solución de la homogénea es Ancos(ωnt)+Bnsin(ωnt)

La solución particular de la ecuación diferencial es

C= g 4 ( 2n+1 )π c 2 ( 2n+1 2l π ) 2 = 16g l 2 c 2 ( 2n+1 ) 3 π 3

La forma funcional de un(t) es

u n (t)= A n cos( 2n+1 2 π ct l )+ B n sin( 2n+1 2 π ct l )+ 16g l 2 c 2 ( 2n+1 ) 3 π 3

El desplazamiento de cualquier punto x de la barra en el instante t es

ψ(x,t)= n=0 { A n cos( 2n+1 2 π ct l )+ B n sin( 2n+1 2 π ct l )+ 16g l 2 c 2 ( 2n+1 ) 3 π 3 }sin( 2n+1 2 π x l )

Donde An y Bn se calculan a partir de las condiciones iniciales

A n = 16g l 2 c 2 ( 2n+1 ) 3 π 3 , B n =0 ψ(x,t)= 16g l 2 c 2 π 3 n=0 1 ( 2n+1 ) 3 ( 1cos( 2n+1 2 π ct l ) )sin( 2n+1 2 π x l ) ψ(x,t)= 32g l 2 c 2 π 3 n=0 1 ( 2n+1 ) 3 sin 2 ( 2n+1 4 π ct l )sin( 2n+1 2 π x l )

Calculamos el desplazamiento del extremo de la varilla x=l en cada instante

ψ(l,t)= 32g l 2 c 2 π 3 n=0 (1) n ( 2n+1 ) 3 sin 2 ( 2n+1 4 π ct l )

Su valor máximo se alcanza cuando ct/(2l)=1,3,5.... El periodo de la oscilación, por ejemplo, del extremo inferior es 4l/c

ψ máx = 32g l 2 c 2 π 3 n=0 (1) n ( 2n+1 ) 3 = 32g l 2 c 2 π 3 π 3 32 = g l 2 c 2

Comprobamos que la suma de la serie vale π3/32

>> syms n;
>> symsum((-1)^n/(2*n+1)^3,n,0,inf)
ans =pi^3/32   

La longitud de la varilla varía entre l y l+gl2/c2

Representamos el desplazamiento ψ(x,t) en función del tiempo, de los puntos situados en x=l/4, x=l/2, y x=l (extremo de la varilla), en este último caso, utilizamos la expresión simplificada ψ(l,t)

L=1; %longitud de la varilla
c=2; %velcidad de propagación
t=linspace(0,4*L/c,200);
hold on
y=zeros(1,length(t));
x=L/4; %posición
for n=0:20 %aproxima con20 términos
    y=y+sin((2*n+1)*pi*x/(2*L))*sin((2*n+1)*pi*c*t/(4*L)).^2/(2*n+1)^3;
end
y=y*(32*9.8*L^2/(c^2*pi^3));
plot(t,y,'b')
x=L/2; %posición
y=zeros(1,length(t));
for n=0:20 %aproxima con20 términos
    y=y+sin((2*n+1)*pi*x/(2*L))*sin((2*n+1)*pi*c*t/(4*L)).^2/(2*n+1)^3;
end
y=y*(32*9.8*L^2/(c^2*pi^3));
plot(t,y,'r')
x=L; %posición
y=zeros(1,length(x));
for n=0:20 %aproxima con 20 términos
    y=y+(-1)^n*sin((2*n+1)*pi*c*t/(4*L)).^2/(2*n+1)^3;
end
y=y*(32*9.8*L^2/(c^2*pi^3));
plot(t,y,'k')
hold off
legend('L/4','L/2','L')
xlabel('t')
ylabel('\Psi(x,t)')
title('Deformación de la varilla')
grid on

El periodo de la oscilación del extremo de la varilla es 4l/c=2, y la máxima deformación del extremo x=l de la varilla es gl2/c2=2.45

Actividades

Se introduce

Se pulsa el botón titulado Nuevo

Aparece la varilla en posición vertical sujeta por su extremo superior, se sostiene el extremo inferior en la posición x=l, varilla sin deformar

Cuando se pulsa el botón titulado , el extremo inferior de la varilla se suelta en el instante t=0

Observamos la vibración de los puntos de la varilla y en particular, los puntos situados en x=l/2 en color azul

En la parte superior derecha, se proporcionan los datos del tiempo t y de la posición l+ψ(l,t) del extremo inferior de la varilla

El periodo de oscilación es


Referencias

N. S. Kosholyakov, M. M. Smirnov, E. B. Gliner. Differential equations of Mathematical Physics. North-Holland Publishing Company (1964) pp. 123-126, 139-141