Modos de Vibración de Una Membrana Circular
Modos de Vibración de Una Membrana Circular
Modos de Vibración de Una Membrana Circular
membrana circular
En esta página, se estudian los modos de vibración de una membrana circular de
radio R y de un anillo de radio exterior a e interior b
Consideremos una membrana circular de radio R, que está fija por su borde. La
ecuación del movimiento ondulatorio en coordenadas polares se escribe.
∂2Ψ∂t2=c2(∂2Ψ∂r2+1r∂Ψ∂r+1r2∂2Ψ∂θ2)∂2Ψ∂t2=c2(∂2Ψ∂r2+1r∂Ψ∂
r+1r2∂2Ψ∂θ2)
Simetría radial
Estudiamos primero, el caso más sencillo, las soluciones Ψ(r,t) independientes
del ángulo θ.
∂2Ψ∂t2=c2(∂2Ψ∂r2+1r∂Ψ∂r)∂2Ψ∂t2=c2(∂2Ψ∂r2+1r∂Ψ∂r)
Ψ(R,0)=0t≥0Ψ(R,0)=0 t≥0
t=0⎧⎨⎩Ψ(r,0)=Ψ0(r)∂Ψ(r,t)∂t∣∣t=0=˙Ψ0(r)t=0
{Ψ(r,0)=Ψ0(r)∂Ψ(r,t)∂t|t=0=Ψ˙0(r)
Modos de vibración
Aplicamos el método de separación de variables
Ψ(r,t)=Φ(r)⋅T(t)Ψ(r,t)=Φ(r)·T(t)
1T(t)⋅d2Tdt2=c2Φ(r)
(d2Φdr2+1rdΦdr)=−ω2⎧⎨⎩d2Φdr2+1rdΦdr+ωc22Φ=0d2Tdt2+ω2T(t)=01T(t)·d
2Tdt2=c2Φ(r)
(d2Φdr2+1rdΦdr)=−ω2{d2Φdr2+1rdΦdr+ωc22Φ=0d2Tdt2+ω2
T(t)=0
d2Φds2+1sdΦds+Φ=0s2d2Φds2+sdΦds+s2Φ=0d2Φds2+1sdΦds+Φ=0s2
d2Φds2+sdΦds+s2Φ=0
Φ=C⋅J0(s)+D⋅Y0(s)Φ=C·J0(s)+D·Y0(s)
Φ(r)=J0(ωcr)Φ(r)=J0(ωcr)
>> alfa=linspace(0,20,200);
>> plot(alfa,besselj(0,alfa),'b')
>> xlabel('\alpha')
>> ylabel('y')
>> title('Raíces de J_0')
>> grid on
Con la herramienta Data cursor, determinamos los valores aproximados de las
raíces: 2.4, 5.5, 8.6, 11.8, 14.9...
ans = 2.4048
ans = 5.5201
ans = 8.6537
ans = 11.7915
ans = 14.9309
La diferencia entre dos raíces consecutivas se aproxima a π.
alpha=[2.4048,5.5201,8.6537,11.7915,14.9309];
R=1; %radio de la membrana circular
r=0:0.01:1;
ang=0:pi/50:2*pi;
[r,ang]=meshgrid(r,ang);
x=r.*cos(ang);
y=r.*sin(ang);
z=besselj(0,alpha(3)*r/R); %tercer modo
surf(x,y,z)
xlabel('x')
ylabel('y')
zlabel('z')
title('Vibración de una membrana circular')
d2Tdt2+ω2⋅T=0Tn(t)=Ancos(ωnt)
+Bnsin(ωnt)d2Tdt2+ω2·T=0Tn(t)=Ancos(ωnt)+Bnsin(ωnt)
Ψn(r,t)=Φn(r)Tn(t)=J0(αnRr)(Ancos(ωnt)
+Bnsin(ωnt))Ψn(r,t)=Φn(r)Tn(t)=J0(αnRr)(Ancos(ωnt)
+Bnsin(ωnt))
Solución completa
La solución general que satisface las condiciones de contorno es la
superposición
Ψ(r,t)=∞∑n=1Ψn(r,t)=∞∑n=1J0(αnRr)(Ancos(ωnt)
+Bnsin(ωnt))Ψ(r,t)=∑n=1∞Ψn(r,t)=∑n=1∞J0(αnRr)(Ancos(ωnt)
+Bnsin(ωnt))
Esta ecuación describe todos los posibles modos de vibración de la membrana
circular. La vibración particular que experimenta la membrana está únicamente
determinada por las condiciones iniciales, que a su vez determinan los valores
de las constantes An y Bn.
Ψ0(r)=∞∑n=1J0(αnRr)AnΨ0(r)=∑n=1∞J0(αnRr)An
∫ ∫
R 0rΨ0(r)J0(αnRr)dr=∞∑n=1AnR∫0rJ0(αnRr)J0(αmRr)drR 0rΨ0(r)J0(αnRr)dr=
R2∞∑n=1An1∫0uJ0(αnu)J0(αmu)duAn=2R2J21(αn)R∫0rΨ0(r)J0(αnRr)drn=1,2,3.
..∫0RrΨ0(r)J0(αnRr)dr=∑n=1∞An∫0RrJ0(αnRr)J0(αmRr)dr∫0Rr
Ψ0(r)J0(αnRr)dr=R2∑n=1∞An∫01uJ0(αnu)J0(αmu)duAn=2R2J1
2(αn)∫0RrΨ0(r)J0(αnRr)dr n=1,2,3...
∫
∂Ψ(r,t)∂t∣∣t=0=∞∑n=1J0(αnRr)ωnBnBn=2R2J21(αn)ωnR 0r(∂Ψ(r,t)∂t∣∣t=0)J0(αnRr)drn
=1,2,3...∂Ψ(r,t)∂t|
t=0=∑n=1∞J0(αnRr)ωnBnBn=2R2J12(αn)ωn∫0Rr(∂Ψ(r,t)∂t|
t=0)J0(αnRr)dr n=1,2,3...
∂2Ψ∂t2=c2(∂2Ψ∂r2+1r∂Ψ∂r+1r2∂2Ψ∂θ2)∂2Ψ∂t2=c2(∂2Ψ∂r2+1r∂Ψ∂
r+1r2∂2Ψ∂θ2)
Ψ(r,t)=Φ(r,θ)⋅T(t)1T(t)⋅d2Tdt2=c2Φ(r,θ)
(∂2Φ∂r2+1r∂Φ∂r+1r2∂2Φ∂θ2)=−ω2Ψ(r,t)=Φ(r,θ)·T(t)1T(t)·d2Tdt2=c2Φ(r,
θ)(∂2Φ∂r2+1r∂Φ∂r+1r2∂2Φ∂θ2)=−ω2
⎧⎨⎩∂2Φ∂r2+1r∂Φ∂r+1r2∂2Φ∂θ2+ω2c2Φ=0d2Tdt2+ω2T(t)=0{∂2Φ∂r2+1r∂Φ∂r
+1r2∂2Φ∂θ2+ω2c2Φ=0d2Tdt2+ω2T(t)=0
Aplicamos el método de separación de variables a la primera ecuación
diferencial
Φ(r,θ)=R(r)⋅Θ(θ)Θ(θ)(r2d2Rdr2+rdRdr)+R(r)d2Θdθ2+ω2c2r2R(r)Θ(θ)=0−1R(r)
(r2d2Rdr2+rdRdr+ω2c2r2R(r))=1Θ(θ)d2Θdθ2=−m2Φ(r,θ)=R(r)·Θ(θ)Θ(θ)
(r2d2Rdr2+rdRdr)+R(r)d2Θdθ2+ω2c2r2R(r)Θ(θ)=0−1R(r)
(r2d2Rdr2+rdRdr+ω2c2r2R(r))=1Θ(θ)d2Θdθ2=−m2
⎧⎪⎨⎪⎩d2Rdr2+1rdRdr+(ω2c2−m2r2)R=0d2Θdθ2+m2Θ=0{d2Rdr2+1rdRdr+
(ω2c2−m2r2)R=0d2Θdθ2+m2Θ=0
{R(r)=CrJm(ωcr)+DrYm(ωcr)Θ(θ)=Cθcos(mθ)+Dθsin(mθ)
{R(r)=CrJm(ωcr)+DrYm(ωcr)Θ(θ)=Cθcos(mθ)+Dθsin(mθ)
x=linspace(0,20,20);
for m=0:3;
alfa=raices(@(x) besselj(m,x),x);
disp([m,alfa])
end
0 2.4048 5.5201 8.6537 11.7915 14.9309
18.0711
1.0000 3.8317 7.0156 10.1735 13.3237 16.4706
19.6159
2.0000 5.1356 8.4172 11.6198 14.7960 17.9598
3.0000 6.3802 9.7610 13.0152 16.2235 19.4094
Nota: la función raices se describe en la página titulada Funciones de
Bessel donde se calculan las raíces y se representan las funciones Jm(x),
para m=0,1,2..,
Ψm,n(r,θ,t)=Jm(αm,nrR)⋅(Cncos(mθ)+Dnsin(mθ))⋅(Ancos(ωm,nt)
+Bnsin(ωm,nt))m=0,1,2,...n=1,2,3...Ψm,n(r,θ,t)=Jm(αm,nrR)·(Cn
cos(mθ)+Dnsin(mθ))·(Ancos(ωm,nt)
+Bnsin(ωm,nt))m=0,1,2,... n=1,2,3...
Los primeros modos de vibración de una membrana circular. Las líneas en color
rojo señalan los nodos, los puntos de la membrana que no vibran para un modo
particular. Más abajo, veremos la representación en 3D de los distintos modos
(m,n) de vibración
Jm(αm,nrR)=0Jm(αm,nrR)=0
Consideremos el caso n=2. Supongamos que hemos calculado las raíces (véase
la tabla más arriba) α01=5.5201 (segunda raíz de J0(x)), α11=7.0156 (segunda
raíz de J1(x)), α21=8.4172 (segunda raíz de J2(x))
alfa=[5.5201,7.0156,8.4172];
for m=0:2
f=@(x) besselj(m,alfa(m+1)*x);
r=fzero(f, 0.5);
disp(r)
end
Radios r/R de los nodos
0.4356
0.5462
0.6101
Consideremos el caso n=3. Supongamos que hemos calculado las raíces (véase
la tabla más arriba) α02=8.6537 (tercera raíz de J0(x)), α12=10.1735 (tercera
raíz de J1(x)), α22=11.6198 (tercera raíz de J2(x))
alfa=[8.6537,10.1735,11.6198];
for m=0:2
f=@(x) besselj(m,alfa(m+1)*x);
r1=fzero(f, 0.2);
r2=fzero(f, 0.8);
disp([r1,r2])
end
Radios r/R de los nodos
0.2779 0.6379
0.3766 0.6896
0.4420 0.7244
Representamos el modo de vibración de la membrana circular, para m=1 y
para n=3, tercera raíz de J1(x)
Ψm,n(r,θ,t)=Jm(αm,nrR)⋅cos(mθ)⋅cos(ωm,nt)Ψm,n(r,θ,t)=Jm(αm,nrR
)·cos(mθ)·cos(ωm,nt)
frame=getframe;
im = frame2im(frame);
[imind,cm] = rgb2ind(im,256);
if t==0
imwrite(imind,cm,fichero,'gif','DelayTime',0,'loopcount',inf);
else
imwrite(imind,cm,fichero,'gif','DelayTime',0,'writemode','append');
end
end
{R(r)=CrJm(ωcr)+DrYm(ωcr)Θ(θ)=Cθcos(mθ)+Dθsin(mθ)
{R(r)=CrJm(ωcr)+DrYm(ωcr)Θ(θ)=Cθcos(mθ)+Dθsin(mθ)
La función R(r) se tiene que hacer cero para r=a y para r=b. Obtenemos un
sistema homogéneo de dos ecuaciones con dos incógnitas C y D.
{CJm(ωca)+DYm(ωca)=0CJm(ωcb)+DYm(ωcb)=0{CJm(ωca)
+DYm(ωca)=0CJm(ωcb)+DYm(ωcb)=0
Jm(ωca)Ym(ωcb)−Ym(ωca)Jm(ωcb)=0Jm(ωca)Ym(ωcb)
−Ym(ωca)Jm(ωcb)=0
a=4; %radios
b=2;
x=linspace(0,10,20);
for m=0:3;
f=@(x) bessely(m,x*b).*besselj(m,x*a)-
besselj(m,x*b).*bessely(m,x*a);
k=raices(f,x);
disp([m,k])
end
0 1.5615 3.1367 4.7091 6.2807 7.8520
9.4231
1.0000 1.5983 3.1562 4.7222 6.2906 7.8599
9.4297
2.0000 1.7035 3.2139 4.7614 6.3202 7.8837
9.4496
3.0000 1.8644 3.3080 4.8261 6.3692 7.9231
9.4825
La solución R(r) es
Rmn(r)=Ym(ωmncb)Jm(ωmncr)
−Jm(ωmncb)Ym(ωcr)Rmn(r)=Ym(ωmncb)Jm(ωmncr)
−Jm(ωmncb)Ym(ωcr)
Ψm,n(r,θ,t)={Ym(ωmncb)Jm(ωmncr)−Jm(ωmncb)Ym(ωcr)}⋅(Cncos(mθ)
+Dnsin(mθ))⋅(Ancos(ωm,nt)
+Bnsin(ωm,nt))m=0,1,2,...n=1,2,3...Ψm,n(r,θ,t)={Ym(ωmncb)Jm
(ωmncr)−Jm(ωmncb)Ym(ωcr)}·(Cncos(mθ)
+Dnsin(mθ))·(Ancos(ωm,nt)+Bnsin(ωm,nt))m=0,1,2,...
n=1,2,3...
FUENTE: http://www.sc.ehu.es/sbweb/fisica3/ondas/membrana/membrana.html