Enonce tp4

4M054 - TP4
Pierre Marchand : [email protected]
Février 2017
Objectifs de la séance :
• Débuter la mise en œuvre de la méthode des éléments finis.
• Calculez les matrices élémentaires pour les matrices de masse et de rigidité.
• Implémentez les routines d’assemblage de ces dernières.
Méthodes des éléments finis - assemblage

1 Problème considéré
Soit Ω ⊂ R2 , n le champs de vecteur normal à ∂Ω dirigé vers l’extérieur de Ω et f ∈ L 2 (Ω). Considérons

le problème suivant : trouver u ∈ H 1 (Ω) telle que
 −∆u + u = f dans Ω

∂u (1)
 = 0 sur ∂Ω.
∂n
Exercice 1 : Formulation variationnelle
→ Proposez une formulation variationnelle pour le problème (1) et montrez que celle-ci est bien
posée.
Soit Th une triangulation du domaine Ω. Notons (Tl )l =1,L les triangles de Th , (S i )i =1,N les nœuds du
maillage et (ϕi )i =1,N les fonctions de base associées à Th dont on note l’espace vectoriel engendré Vh .
Exercice 2 : Formulation variationnelle discrétisée
→ Proposez une formulation variationnelle discrète pour le problème (1) et montrez que celle-ci est
bien posée.
Par construction, la solution approchée u h a la forme suivante : pour tout x ∈ Ω
N
u hi ϕi (x).
X
u h (x) =
i =1
Exercice 3 : Système linéaire

1. Mettez la formulation variationnelle discrète sous la forme d’un système linéaire :
(M + K)u = f (2)
où M et K sont des matrices carrées de taille N tandis que u et f sont des vecteurs de taille N .
2. Donnez des propriétés intéressantes de M et K. Ces propriétés pourront notamment être utilisées
pour vérifier l’assemblage des matrices.
1
2 Assemblage des matrices et du second membre
2.1 Matrice de masse
Soit un triangle Tl de sommets S 1 (x 1 , y 1 ),S 2 (x 2 , y 2 ) et S 3 (x 3 , y 3 ). Les coordonnées barycentriques λ1 λ2

et λ3 du triangle ont pour expression :
1
λ1 (x, y) = (y 23 (x − x 3 ) − x 23 (y − y 3 ))
D
1
λ2 (x, y) = (y 31 (x − x 1 ) − x 31 (y − y 1 ))
D
1
λ3 (x, y) = (y 12 (x − x 2 ) − x 12 (y − y 2 ))
D
où x i , j = x i − x j , y i , j = y i − y j pour i , j entre 1 et 3, et D = x 23 y 31 − x 31 y 23 . Remarquez que D est égale à

deux fois l’aire du triangle au signe près.
Exercice 4 : Matrice élémentaire - matrice de Masse Soit T un triangle de la triangulation Th .
1. Calculez les termes suivants pour tout 1 ≤ i , j ≤ N
Z
ϕi (x, y)ϕ j (x, y)d xd y
T
2. Définissez la fonction MassElem qui prend en argument trois points et qui renvoie la matrice de
Masse élémentaire associé au triangle T formé de ces trois points :
µZ ¶
MT = ϕi ϕ j
T
1
où ϕi est la fonction de base P associée au nœud i .
3. Vérifiez que vous trouvez bien les bonnes valeurs pour un triangle de référence.
Exercice 5 : Assemblage - matrice de Masse
1. Implémentez la boucle d’assemblage de la matrice M.
2. Vérifiez que vous retrouvez bien les propriétés remarquables de cette matrice.
3 Matrice de rigidité
Exercice 6 : Matrice élémentaire - matrice de rigidité

1. Calculez les termes suivants pour tout 1 ≤ i , j ≤ N
Z
∇ϕi (x, y) · ∇ϕ j (x, y)d xd y
Tl
2. Définissez la fonction RigElem qui prend en argument trois points et qui renvoie la matrice de
rigidité élémentaire associé au triangle T formé de ces trois points :
µZ ¶
KT = ∇ϕi · ∇ϕ j
T
1
où ϕi est la fonction de base P associée au nœud i .
3. Vérifiez que vous trouvez bien les bonnes valeurs pour un triangle de référence.
Exercice 7 : Assemblage - matrice de rigidité
1. Implémentez la boucle d’assemblage de la matrice K.
2. Vérifiez que vous retrouvez bien les propriétés remarquables de cette matrice.
2
4 Second membre et convergence
Nous définissons l’opérateur d’interpolation I h sur le maillage Th de la façon suivante :
I h : C 0 (Ω) → Vh
N
X
f 7→ I h f = f (S i )ϕi .
i =1
Nous admettons alors le résultat suivant :

Théorème 1 (Erreur d’interpolation). Il existe des constantes C 0 et C 1 indépendantes de la finesse du
maillage h, telle que pour toute fonction f ∈ H 2 (Ω)
k f − I h f kL 2 (Ω) ≤ C 0 h 2 | f |H 2 (Ω) et | f − I h f |H 1 (Ω) ≤ C 1 h| f |H 2 (Ω) .
Exercice 8 : Convergence
→ En utilisant le théorème précédent et le lemme de Céa, démontrez une convergence d’ordre 1 avec
la semi-norme H1 de la solution approximée vers la solution exacte, si celle-ci est assez régulière,
c’est-à-dire
|u − u h |H1 (Ω) ≤ C h|u|H 2 (Ω) ,
où C est une constante indépendante de h.

Exercice 9 : Second membre
→ On assimilera la fonction f à son interpolée P 1 , déduisez une expression de f faisant intervenir la
matrice de masse.
5 Résolution du problème
Pour tester l’assemblage de nos matrices, nous pouvons résoudre un problème dont nous connaissons
la solution afin de vérifier que notre solution approximée s’en approche bien.
Exercice 10 : Résolution numérique
1. Soit Ω un domaine rectangulaire, choisissez une fonction u qui vérifie les conditions aux limites
de (1). (on pourra prendre un produit de cosinus bien choisi)
2. Déduisez une fonction f telle que la fonction u de la question précédente soit solution de (1).
3. Résolvez numériquement (1) en utilisant le système linéaire (2).
4. Vérifiez graphiquement que la solution numérique u h est une approximation correcte de u en la
représentant graphiquement sur le maillage ainsi que la solution exacte u.
Exercice 11 : Étude de l’erreur
1. Donnez une expression de la semi-norme H 1 de l’erreur, |u − u h |H 1 (Ω) = k∇u − ∇u h kL 2 (Ω) , faisant
intervenir la matrice de rigidité K. Tracez log h 7→ log(|u − u h |H 1 /|u h |H 1 ) pour différentes valeurs
de h. Qu’observez vous ? Expliquez. (on pourra assimiler u à son interpolée)
2. Faites de même avec la norme L 2 de l’erreur, ku − u h kL 2 (Ω) . Que remarquez vous cette fois-ci ?

Enonce tp4

Transféré par

Droits d'auteur :

Formats disponibles

Enonce tp4

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

Enonce tp4

Transféré par

Droits d'auteur :

Formats disponibles

4M054 - TP4

Pierre Marchand : [email protected]

Méthodes des éléments finis - assemblage

Soit Ω ⊂ R2 , n le champs de vecteur normal à ∂Ω dirigé vers l’extérieur de Ω et f ∈ L 2 (Ω). Considérons

Par construction, la solution approchée u h a la forme suivante : pour tout x ∈ Ω

Exercice 3 : Système linéaire

2.1 Matrice de masse

Soit un triangle Tl de sommets S 1 (x 1 , y 1 ),S 2 (x 2 , y 2 ) et S 3 (x 3 , y 3 ). Les coordonnées barycentriques λ1 λ2

où x i , j = x i − x j , y i , j = y i − y j pour i , j entre 1 et 3, et D = x 23 y 31 − x 31 y 23 . Remarquez que D est égale à

Exercice 6 : Matrice élémentaire - matrice de rigidité

Nous définissons l’opérateur d’interpolation I h sur le maillage Th de la façon suivante :

Nous admettons alors le résultat suivant :

k f − I h f kL 2 (Ω) ≤ C 0 h 2 | f |H 2 (Ω) et | f − I h f |H 1 (Ω) ≤ C 1 h| f |H 2 (Ω) .

|u − u h |H1 (Ω) ≤ C h|u|H 2 (Ω) ,

où C est une constante indépendante de h.

Vous aimerez peut-être aussi