4AN01 1617 Ch3

4AN01 : Méthodes numériques avancées
Ch. III :
La méthode des éléments finis
1. Principe de la méthode
2. Exemple 1D
3. Cas 2D ‐ Eléments finis triangulaires P1
1 florence.ossart@upmc.fr
1. Principe de la méthode
⊂
Exemple de problème traité :
ouvert connexe
 ( : trouver telle que : ∆
∆
∪
par morceaux
A priori, on cherche dans Ω , e.v. de dimension infinie !!!

Approximation MEF : on cherche dans , s.e.v. de Ω , de dimension finie N,
avec ∑ .
où , est une famille de fonctions formant une base de , construites à
partir d’un maillage du domaine Ω. On les appelle « fonctions de forme ».
Les inconnues sont les coefficients , coordonnées de dans cette base.
Exemple : densité de courant dans une plaque conductrice
soumise à une différence de potentiel
0, ⟹ . grad 0
V=10 volts
Δ 0
dV/dn=0
(isolation)
V=0
(masse)
Maillage de et approximation de :
avec , ∑ . ,
Inconnues du problème : ,
Maillage :
 Approximation du domaine d’étude Ω par un découpage polygonal.
La MEF s’applique à des problèmes 1D, 2D ou 3D. Différentes formes de mailles
existent pour paver le domaine :
 1D : segment (éventuellement courbe)
 2D : triangle ou rectangle
 3D : tétraèdre, pyramide, prisme, cube
Maillage :
 L’utilisation d’éléments triangulaires (en 2D) ou tétraédriques (en 3D) :
 permet de suivre des formes courbes quelconques,
 permet de raffiner le maillage là ou c’est utile.
Fonctions de forme (d’interpolation) :
 Caractéristiques générales :
 Fonction polynomiale par morceau, continue (morceau = maille)
 1 fonction par nœud du maillage (nœud = sommet des mailles)
 1 et 0 ‐ Autre écriture :
Conséquence : est nulle partout, sauf sur les mailles dont est un
sommet
 Exemples :
P1 (affine par morceau) en 1D P1 (affine par morceau) en 2D :
 Conséquences :
 ∑ . est une fonction polynomiale par morceau, continue

Les coefficients sont directement la valeur de aux nœuds.
Problème discrétisé :
La MEF s’appuie sur la formulation variationnelle du problème.
 (PV) : Trouver ∈ telle que ,

∀ ∈ ∈ H Ω et 0
 Problème approximé :
∀ ∈ ,
Trouver ∑ . ∈ telle que
∑ .
Vrai ∀ ssi : ∀ 1, ,
⟹ ∀ 1, ∑ . ,
⟹ ∀ 1, ∑ . ,
Système d’équations linéaires : ∑ . , dont les inconnues sont ,
avec : , et
2. Exemple 1D : traité au tableau, prendre des notes
 Equation de l’électrostatique en 1D : u est le potentiel électrique

⟹
en tout point de Ω , .
et sont des fonctions d’espace connues.
Conditions aux limites :
 Dirichlet sur le bord gauche : 0
 Neuman homogène sur le bord droit : ′
3. Cas 2D – Interpolation par éléments finis triangulaires P1
 Problème plan, recherche de , – Maillage par des triangles

 Triangulation du domaine = partition formée de triangles ,
 , sommets des triangles
= nœuds du maillage
 Règles de construction :
 le sommet d’un triangle ne doit jamais
appartenir au côté d’un autre triangle
 Aucun triangle ne doit être plat (éviter
les triangles aplatis)
Interpolation P1 par morceaux
 Forme de la solution :
Interpolation P1 : affine par morceau, continue
 Sur chaque triangle de sommets , et , on cherche une approximation
de sous la forme : ,
 Soient , et , les valeurs de aux sommets , et
, et sont solutions du système : .
Si le triangle n’est pas plat, il existe une solution unique sur chaque triangle.
 Une arête est commune à 2 triangles, mais le long de cette arête,
, ne dépend que de et , pas de au 3ème sommet. Il s’ensuit que la
fonction , est continue lors du passage d’un triangle à un autre.
Base canonique sur un triangle
 Triangle quelconque, avec numérotation locale des nœuds
 On introduit les fonctions affines , , , et , définies par :
, 0 0
, 0 0
, 0 0
 , , est une base canonique de P1:
 , , , 0⇒ 0
 Soit un polynôme de degré 1. Il est entièrement défini par ses valeurs
aux nœuds , pour 1,3
, . , . , . ,
 , , est une famille libre, génératrice de l’espace des polynômes
de degré 1, donc c’est une base, dite canonique relativement au triangle
et 1 3.
Base canonique sur tout entier
 Numérotation globale des nœuds de Ω ⟹ nœuds (sommets des triangles)
 : espace vectoriel des fonctions affines par triangle et continues sur Ω
 Les fonctions de sont déterminées de manière unique par leurs valeurs aux
sommets des triangles (c'est‐à‐dire aux nœuds), donc
 Soit , , une famille de fonctions de définies par :
∀ 1, (ce sont les fonctions chapeaux)
 Considérons une fonction de la forme : ∑ .

∀ 1, ∑ . ∑ . , donc .
En particulier, si est la fonction nulle : 0.

Conséquence : , est une famille de fonctions linéairement
indépendantes, c’est donc une base de .
On l’appelle base canonique de
∀ ∈ , ∑ . , ∑ . ,
Approximation du problème variationnel
Problème discrétisé :
La MEF s’appuie sur la formulation variationnelle du problème.
 (PV) : Trouver ∈ telle que ,

∀ ∈ ∈ H Ω et 0
 Problème approximé :
∀ ∈ ,
Trouver ∑ . ∈ telle que
∑ .
MAIS… 0 ⟹ 0 pour les fonctions attachées à ces nœuds,

donc dim
Par un choix adéquat de la numérotation des nœuds, on peut se ramener à :
∑ . et ∑ .
Résolution du problème discrétisé
 Pour avoir , ∀ ∈ , il suffit que ce soit vérifié pour les fonctions
de la base de .
 ∀ 1, ∑ . ,
 ∀ 1, ∑ . ,
 Convention de notation d’Einstein : . , ‐ sommation sur

l’indice répété
 Notation matricielle : . ,
où : , avec , est la matrice de « rigidité » du système
et est le terme « source ».

 NB : 0 uniquement si ou si et sont adjacents, donc la matrice A
est très creuse, et à diagonale dominante si on numérote astucieusement les
nœuds.
Exemple d’application
∆
(PC)
MEF :
On cherche ∑ .
Système d’équations linéaires :

∑ .
avec : .
et . . . .
1
Nœud intérieur au domaine :
3 2 calculs développés au tableau
4 0
0 1 1
4
.
6
5
⟹ 4. .
Nœud sur la frontière :

calculs développés au tableau
2 0
3 2 . .
1 1 1
0 ⟹ 2.
4 2 2
. .
Matrice élémentaire, associée à (012) :
Le triangle (012) contribue aux coefficients
, , , mais aussi , , , ,
, ,
2
calculs développés au tableau
1
1
0 012 0
0
012 constitue la contribution de (012)

dans la matrice globale A, laquelle est
construite par assemblage des matrices
élémentaires de chaque élément.
On fait de même pour le second membre.
Eléments finis de Lagrange
 Définition : triplet , ,
 : polyèdre compact (triangle, pavé, …) courbe ou droit
 : ensemble de points appelés nœuds
 : espace vectoriel de fonctions d’interpolation, de dimension égale au
nombre de nœuds
 Dans l’élément : doit être : les valeurs d’une fonction
de P aux nœuds de Σ doivent la définir de manière unique.
 Sur chaque face : doit être . Σ et sont les
restrictions de Σ et de à une face de .
 Exemple : élément finis triangulaire P1
Convergence de la méthode des éléments finis :
 Soit , la solution de (PV) : On appelle interpolée de sur , que l’on note

, la fonction de définie par : ∀ .
On peut montrer que, pour une discrétisation de plus en plus fine :

‖ ‖ est du même ordre de grandeur que ‖ ‖ .
 Lemme de Bramble‐Hilbert :
Soit , une mesure de la taille du plus grand des éléments finis du maillage, qui est
supposé ne pas comporter d’éléments aplatis.
Si ⊂ et si u est assez régulière, alors l’approximation est d’ordre :
‖ ‖ et donc ‖ ‖
Pour les éléments triangulaires , nous avons construit une approximation

d’ordre 1 (en fait, ordre 2 sur et ordre 1 sur ).
Eléments finis P2 du second ordre :
Pour augmenter l’ordre, il faut augmenter le degré des polynômes d’interpolation,
et donc augmenter le nombre de nœuds.
 Eléments finis droits :
o : triangle ( )
o Σ : nœuds , , , , et
(sommets et milieux des 3 côtés)
o : espace des polynômes de degré 2

, 2
Toute fonction de est déterminée de façon unique par 6 coefficients, donc dim 6.
o Base canonique : , , , , , avec :
o . 2 1 pour 1,3
o 4 . pour 12, 13, 23
 Elément de référence : le plan « chapeau »
Expression des fonctions de forme :

; ;
1
d’où :
. 2 1
. 2 1
1 . 1 2 2
4 .
4 . 1
4 . 1
Eléments finis Q1, quadrilatéral du premier ordre :
 Dans le plan chapeau :
o : carré ( )
o Σ : nœuds , , ,
o : espace des polynômes de degré 1 par

rapport à et
,
Toute fonction de est déterminée de façon unique

par 4 coefficients, donc dim 4.
o Base canonique :
o 1 1 o
o 1 o 1
4. Pour conclure…
 Ce cours n’est qu’un aperçu de la méthode des éléments finis. Il existe bien
d’autres types d’éléments…
 Vous serez sans doute amené un jour au l’autre, et sans doute rapidement, à
utiliser un logiciel de simulation par la MEF. Vous aurez alors une petite idée de
ce qu’il y a sous le capot…

4AN01 1617 Ch3

Transféré par

Droits d'auteur :

Formats disponibles

4AN01 1617 Ch3

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

4AN01 1617 Ch3

Transféré par

Droits d'auteur :

Formats disponibles

4AN01 : Méthodes numériques avancées

A priori, on cherche dans Ω , e.v. de dimension infinie !!!

 ∑ . est une fonction polynomiale par morceau, continue

 (PV) : Trouver ∈ telle que ,

Système d’équations linéaires : ∑ . , dont les inconnues sont ,

 Equation de l’électrostatique en 1D : u est le potentiel électrique

 Problème plan, recherche de , – Maillage par des triangles

 , sommets des triangles

, et sont solutions du système : .

∀ 1, (ce sont les fonctions chapeaux)

 Considérons une fonction de la forme : ∑ .

En particulier, si est la fonction nulle : 0.

 (PV) : Trouver ∈ telle que ,

MAIS… 0 ⟹ 0 pour les fonctions attachées à ces nœuds,

Par un choix adéquat de la numérotation des nœuds, on peut se ramener à :

 Convention de notation d’Einstein : . , ‐ sommation sur

et est le terme « source ».

Système d’équations linéaires :

Nœud sur la frontière :

012 constitue la contribution de (012)

 Exemple : élément finis triangulaire P1

 Soit , la solution de (PV) : On appelle interpolée de sur , que l’on note

On peut montrer que, pour une discrétisation de plus en plus fine :

Pour les éléments triangulaires , nous avons construit une approximation

o : espace des polynômes de degré 2

Expression des fonctions de forme :

o : espace des polynômes de degré 1 par

Toute fonction de est déterminée de façon unique

Vous aimerez peut-être aussi