Bonjours à tous,

Je me permet de vous demander de l'aide pour un projet que j'effectue dans le cadre de ma scolarité.
Il s'agit de vanter les mérites de la méthode par élements finis avec l'aide de scilab. Je développe en exemple un problème d'acoustique en 1D.
après obtention de la formulation faible et discrétisation, nous obtenons :


equ1.jpg


p(x,t) la surpression dans notre conduit;
phi est la fonction de forme.

Les fonction de formes choisis sont : N1 = (1-s)/2 et N2=(1+s)/2 , s est compris entre -1 et 1

J'ai des problèmes pour le passage en forme matriciel (a gauche):
On a l'integrale sur 1 element du maillage :

equ2.jpg


Visiblement, il y a une erreur quelque part... Le n'apparait nul part...
Comment intervient le changement de variable? Un jacobien apparait il?
D'autre part, phi n'a pas l'air de vouloir disparaitre, une petite piste pour que je puisse m'en débarasser?


La méthode des élements finis, que je trouve très interessant, me pose problème malgré tout... Je suis au bord de la déprime


Je vous serai fort reconnaissant pour votre aide...