Introduction à la méthode des éléments finis

On a vu dans l’article [AF 503] qu’une équation aux dérivées partielles ellip-tique pouvait être exprimée sous diverses formulations équivalentes, en ce sens que toute solution d’une formulation est solution d’une autre formulation et réciproquement. La formulation forte du problème présente un intérêt dans la mesure où l’utilisation de la méthode des différences finies est envisagée pour effectuer une approximation du problème. La formulation équivalente du problème basée sur la formulation d’un problème d’optimisation associé à la fonctionnelle $v\to J\left(v\right)$ , avec $J\left(v\right)$ définie par : $J\left(v\right)=\frac{1}{2}a\left(v,\mathrm{\hspace{0.17em}}v\right)-L\left(v\right)$

nécessite que la forme bilinéaire $a\left(u,\mathrm{\hspace{0.17em}}v\right)$ soit symétrique, ce qui en soit est restrictif dans la mesure où certains phénomènes sont modélisés à partir d’opérateurs non autoadjoints. Cependant, lorsque a(.,.) est symétrique, cette formulation du problème conduit à la méthode de Ritz ; numériquement, l’idée est de chercher à minimiser J(.) non plus sur l’ensemble E tout entier, mais sur un sous-espace de E construit à partir de fonctions facilement calculables ; la fonction inconnue qui réalise le minimum est représentée comme combinaison linéaire de fonction de forme (ou de tout autre famille physiquement admissible) et les coefficients de cette combinaison linéaire sont les inconnues du problème. J(.) est alors transformée en une fonctionnelle quadratique et déterminer le minimum de cette nouvelle fonctionnelle revient alors à annuler les dérivées partielles de celle-ci par rapport à ces inconnues, ce qui conduit classiquement à la résolution d’un système linéaire. Nous ne développerons pas cette méthode .

L’autre formulation équivalente mise en évidence est la formulation faible ou...