Ce programme a pour but de calculer la fréquence en tenant compte de la dispersion. Il calcule cette fréquence en utilisant deux méthodes différentes la première consistant à calculer oméga et la deuxième en déterminant la pulsation.
Voici le code réalisé sous Matlab, qui nous a permit de faire ce programme.
figure
h=.352/2;
Vt=4670;
Vl=9134;
a=5.2e-2;
b=3.6e-2;
dx=0.01e-2;
x=0:dx:a;
y=0:dx:b;
e=450e9;%modul d'young
%h=352e-6;%epaisseur suposé a changer!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
nu=0.42;%coef poisson
rho=2330;%masse volumique
%i=1
for m=1:5;
for n=1:5;
%i=1+i;
wx=cos(m*pi*x/a);
wy=cos(n*pi*y/b);
w=wx'*wy;
subplot(2,1,1)
%%surf(w)
%%shading interp
%%colorbar
imagesc(x,y,w);
%savename = (['m=',num2str(m),'n=',num2str(n),'image']);
%saveas(gcf,savename,'jpg');
k=sqrt((n*pi/a)^2+(m*pi/b)^2);
Vp=2*Vt*(1-(Vt^2/Vl^2));
ohmega=1.15*(Vp./sqrt(3)).*k.^2*h*1e-3;
freqreel=ohmega/(2*pi);
subplot(2,1,2)
masse=rho*h;
D=e*h^3/(12*(1-nu^2));
%puls=pi^2*sqrt(D/masse)*((m^2/a^2)+(n^2/b^2));
%fre=puls/(2*pi);
fre=freqreel;
fre=round(fre);
vectfre=zeros(1,20000);
vectfre(fre)=10;
plot(vectfre)
hold on
drawnow,
pause,
close
end
end
