Bonjour a tous
J ai un petit probleme de differences finies. Je souhaite modeliser le comportement d'une plaque d'Aluminium, lorsqu on applique une force dynamique.
J ai donc ouvert matlab, fait mon petit programme en rentrant l equation de Kirchoff, et hop, ca diverge. C est pourquoi je voulais savoir ce que vous me conceille de prendre comme schema pour une equation en laplacien carre.
Plus precisement, voila mon equation
w_new = (2*w_new-w) ... -interval^2*(f(t)*heaviside(-abs(x)-abs(y)+0.1)/mass ... +1/2*(16*D/mass/10^-6*del2(del2(w_new))+D/mass10^-6*del2(del2(w))));
la premiere ligne est m*d^2w/dt^2
la deuxieme est une force appliquee sur le point (0,0)
et la troisieme est le Laplacien carre.
Sachant que pour les conditions aux limites, j impose a chaque pas, une bande de 2 indice de large qui est egale a zero (ce qui me permet d avoir un bloquage)
En vous remerciant
-----