Ce programme a pour but de calculer la fréquence des modes sans tenir compte de la dispersion. Apres l avoir calculé il affiche la fréquence trouvée.


Voici le code réalisé sous Matlab, qui nous a permit de faire ce programme avec le resultat obtenue.

figure
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=480e-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');

   
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=round(fre);
vectfre=zeros(1,20000);
vectfre(fre)=10;
plot(vectfre)
hold on
drawnow,
pause,
close
 end
 end