méthode des éléments finis

Imaginez une pièce mécanique dont le contour comporte des trous et des courbes. Chercher d’un seul coup sa déformation en chaque point serait difficile. La méthode des éléments finis découpe alors la pièce en petites zones simples, reliées par des nœuds. Sur chaque zone, elle décrit la déformation par une fonction élémentaire.

Les valeurs inconnues aux nœuds sont calculées ensemble. Une fois ces valeurs trouvées, les fonctions locales se raccordent pour donner une approximation sur toute la pièce. Un maillage plus fin peut mieux suivre les variations, mais demande davantage de calculs.

La méthode des éléments finis est une famille de méthodes numériques pour approcher la solution d’une équation aux dérivées partielles sur un domaine continu. Le domaine est partagé en éléments simples : intervalles en une dimension, souvent triangles ou quadrilatères en deux dimensions, tétraèdres ou hexaèdres en trois dimensions. Cet ensemble d’éléments forme le maillage.

La formulation variationnelle remplace l’équation ponctuelle par des égalités intégrales testées contre des fonctions admissibles. On choisit ensuite un espace de dimension finie engendré par des fonctions de base locales. Si φ_i désigne la fonction associée au nœud i et U_i son coefficient inconnu, l’approximation u_h s’écrit $u_h={\sum }_{i=1}^N{U}_i{\phi }_{i}u\_h=\backslash sum\_\{i=1\}^\{N\}U\_i\backslash varphi\_i$. Le support local signifie qu’une fonction de base n’agit que sur quelques éléments voisins.

En reportant cette combinaison dans la formulation variationnelle, on obtient un système algébrique pour les coefficients U_i. Les conditions aux limites déterminent quelles valeurs sont imposées et quelles valeurs restent inconnues. La qualité du résultat dépend du maillage, du degré des fonctions polynomiales et de la régularité de la solution recherchée.

Pour appliquer la méthode, on établit d’abord une formulation variationnelle avec ses conditions aux limites. On maille le domaine, puis on choisit des fonctions de base locales compatibles avec ce cadre. L’approximation est injectée dans les égalités intégrales.

Si K désigne la matrice d’assemblage, U le vecteur des coefficients et F le vecteur des données, l’assemblage aboutit au système $KU=FKU=F$. Lorsque des valeurs essentielles U_b sont imposées au bord, on sépare les coefficients inconnus U_u et on résout le système réduit $K_{uu}{U}_u=F_u-K_{ub}{U}_{b}K\_\{uu\}U\_u=F\_u-K\_\{ub\}U\_b$. L’ensemble des coefficients ainsi obtenus détermine entièrement la solution approchée u_h.

La méthode s’applique lorsque le problème possède une formulation variationnelle bien posée et que le domaine peut être maillé. Les données comprennent l’équation, les coefficients physiques, les sollicitations et les conditions aux limites. L’espace discret doit respecter les régularités et les conditions essentielles exigées par cette formulation.

Trois contrôles sont concrets : les éléments couvrent le domaine sans lacune ni recouvrement indésirable ; les fonctions voisines respectent les raccordements imposés par la formulation — la continuité requise dans un espace conforme, ou les conditions de couplage dans une méthode discontinue ; le système assemblé possède une solution unique après imposition des conditions au bord. Si un problème de déformation ne bloque aucun déplacement rigide, la matrice reste singulière et le calcul ne fournit pas une solution unique. Il faut alors ajouter des conditions d’appui appropriées ou reformuler le modèle.

On cherche une approximation sur l’intervalle [0 ; 1] de la fonction u soumise à l’équation −u″ = 1, avec u(0) = 0 et u(1) = 0. Le maillage comporte deux éléments de longueur 0,5 et trois nœuds placés en 0, 0,5 et 1. On utilise des fonctions affines par morceaux.

Données. Les deux valeurs de bord sont imposées. Il reste une seule inconnue, notée U, au nœud central. La figure matérialise les deux éléments et la fonction affine par morceaux obtenue.

1. Sur un élément de longueur 0,5, la matrice locale vaut :

2. Les deux éléments contribuent au nœud central : le coefficient assemblé vaut 2 + 2 = 4. La charge constante apporte 0,25 par élément, soit 0,5 au total.

3. Le système réduit est 4U = 0,5, donc U = 0,125. L’approximation relie alors les valeurs 0, 0,125 et 0 par deux segments.

Contrôle. La solution exacte est u(x) = x(1 − x)/2 ; elle vaut aussi 0,125 en x = 0,5. L’égalité au nœud central ne signifie pas que les deux courbes coïncident entre les nœuds.

En mécanique des structures, on maille une pièce pour estimer ses déplacements et ses contraintes sous des charges données. Une formule analytique reste préférable pour une forme très simple ; les éléments finis deviennent utiles lorsque la géométrie, les matériaux ou les appuis compliquent le problème.

En thermique, le maillage sert à approcher la température dans un objet dont les bords sont chauffés, refroidis ou isolés. Le praticien raffine surtout les zones où la température varie rapidement, puis compare plusieurs tailles de maille.

En électromagnétisme ou en mécanique des fluides, le même principe transforme des champs continus en coefficients calculables. Le choix des éléments et de la formulation dépend toutefois de la continuité physique à préserver ; une discrétisation standard n’est pas interchangeable avec toutes les autres.

Éléments finis et différences finies. Les éléments finis partent d’une formulation variationnelle et de fonctions de base portées par un maillage. Les différences finies remplacent directement les dérivées par des relations entre valeurs voisines. Sur un rectangle régulièrement quadrillé, les deux calculs peuvent sembler proches ; une géométrie triangulée révèle leur construction différente.

Éléments finis et volumes finis. Les volumes finis imposent des bilans de flux sur chaque volume de contrôle. Les éléments finis imposent les égalités de la formulation variationnelle sur un espace de fonctions tests. Dans leurs formulations standard, les deux approches se distinguent donc par leur construction ; certaines formulations d’éléments finis, notamment mixtes ou discontinues, peuvent toutefois satisfaire elles aussi des bilans locaux.

Maillage et solution. Le maillage est le découpage géométrique ; la solution approchée est la combinaison des fonctions de base. Deux calculs peuvent donc employer le même maillage mais des degrés polynomiaux différents et produire des approximations différentes.

Maille trop grossière. Dans l’exemple, deux éléments donnent bien u(0,5) = 0,125, mais l’approximation reste formée de segments alors que la solution exacte est courbe. Une valeur nodale juste ne prouve donc pas que l’erreur est nulle entre les nœuds. Il faut mesurer l’erreur ou recommencer avec un maillage raffiné.

Mauvaise qualité des éléments. Des éléments très aplatis ou déformés peuvent dégrader la précision et le conditionnement du système. Le symptôme est une forte sensibilité du résultat ou une résolution numérique difficile. Il faut contrôler la qualité géométrique et remesher les zones concernées.

Modèle mal posé. Raffiner le maillage ne corrige ni une condition au bord manquante ni une loi physique inadéquate. Une matrice singulière ou des solutions qui changent sans se stabiliser signalent ce problème. Le modèle continu et sa formulation variationnelle doivent être vérifiés avant la discrétisation.

Convergence non automatique. Réduire la taille des éléments n’assure une amélioration que si les espaces choisis approchent correctement la solution et respectent la formulation. Les singularités, discontinuités ou couches très fines demandent souvent un raffinement local, un autre degré polynomial ou des éléments adaptés.

La méthode de Galerkin fournit le cadre qui relie fonctions tests, fonctions de base et système algébrique. L’étude de la convergence compare ensuite la solution approchée à la solution continue dans une norme adaptée au problème.

Deux voies améliorent le calcul : le raffinement h réduit la taille de certains éléments, tandis que le raffinement p augmente le degré des polynômes locaux. Les estimateurs d’erreur servent à choisir où concentrer les inconnues plutôt qu’à raffiner uniformément tout le domaine.