ELÉMENTS FINIS 2D / ℙ

Séance 6

MEC4122
Prof. Anas RACHID
Département de mécanique
ENSAM – Casablanca | Université Hassan II de Casablanca
 Problème modèle

Mettre en évidence les techniques de calcul par éléments finis sur le problème modèle
suivant:

∆ , + , + , = ( , ), ∀ , ∈Ω
, = 0, ∀ , ∈ Ω

où Ω est un domaine assez régulier de ℝ .

Dans un premier temps, on prends = = 0.

Formulation variationnelle
Chercher ∈ (Ω) telle que:
( , ) , = , , , ∀ ∈ Ω (1)

Anas RACHID | ENSAM – Casablanca | Université Hassan II de Casablanca 2

 Procédure éléments finis

Création du Matrices Technique Système

maillage élémentaires d’assemblage linéaire

La procédure de l’approximation par la méthode des éléments finis

Anas RACHID | ENSAM – Casablanca | Université Hassan II de Casablanca 3

 Procédure éléments finis

Création de maillage

Définition (maillage admissible)

Un maillage (triangulation) de Ω est une famille finie = d’ouvert borné telle
que:
• ∀ ∈ 1, … , , est non-dégénéré,
• Ω=⋃
l ensemble vide
• ∀ , ∈ 1, … , , ≠ , ∩ = un sommet commun
un coté commun

Anas RACHID | ENSAM – Casablanca | Université Hassan II de Casablanca 4

 Définition (maillage admissible)
Un maillage (triangulation) de Ω est une famille finie = d’ouvert borné telle
que:
Création de Maillage

• ∀ ∈ 1, … , , est non-dégénéré,
• Ω=⋃
l ensemble vide
• ∀ , ∈ 1, … , , ≠ , ∩ = un sommet commun
un coté commun

On pose
ℎ = max ℎ ℎ est appelé la taille du maillage.
∈
Avec ℎ = = max − .
, ∈

Anas RACHID | ENSAM – Casablanca | Université Hassan II de Casablanca 5

 Exemple de maillage non
conforme (non admissible).
Création de Maillage

Maillage triangulaire

Maillage quadrilatère

6
Création de Maillage Exemple de maillage triangulaire

7
Procédure éléments finis

Calcul des matrices élémentaires

Les fonctions de forme

L’élément de référence

La matrice élémentaire

8
 En deux dimensions d’espace les fonctions polynômiales peuvent êtres présentées de deux
manières: ℙ ou ℚ

Les polynômes de type ℙ Les polynômes de type ℚ

ℙ = :Ω → ℝ | , = ℚ = :Ω → ℝ | , =
Fonctions de forme

Où: Où:
• Polynôme de degré 1 • Polynôme de degré 1
∈ℙ ⟺ , = + + ∈ℚ ⟺ , = + + +

• Polynôme de degré 2 • Polynôme de degré 2

∈ℙ ⟺ , = + + + + + ∈ℚ ⟺ ,
= + + + + + + +ℎ
+

9
 Espace d’approximation
Soit Ω = une maillage de Ω par des triangles ou des quadrangles. Sur , on
définit l’espace d’approximation dans le quel on approche la solution de 1 :

→ Si est un maillage triangulaire → Si est un maillage quadrangulaire

Fonctions de forme

= ∈ Ω | | ∈ ℙ ,∀ ∈ Ω = ∈ Ω | | ∈ ℚ ,∀ ∈ Ω
= =

Où, dans chaque triangle la fonction de forme Où, dans chaque quadrangle la fonction de forme
vérifie: vérifie:
1. , = + + 1. , = + + +
2. , = , 1≤ ≤3 2. , = , 1≤ ≤3
, ,
,

, ,
,
,
10
 Espace d’approximation: cas des triangles
Soit Ω = une maillage de Ω par des triangles. Sur , on définit l’espace
d’approximation dans le quel on approche la solution de 1 :

→ Si est un maillage triangulaire

Fonctions de forme

= ∈ Ω | | ∈ ℙ ,∀ ∈ Ω
=

Où, dans chaque triangle la fonction de forme

vérifie:
1. , = + +
2. , = , 1≤ ≤3
,

,
,
11
 Espace d’approximation: cas des triangles

Dans l’élément nous avons 3 fonctions de forme :

Nœud ❶ Nœud ❷ Nœud ❸

Fonctions de forme

, =1 , =0 , =0
, =0 , =1 , =0
, =0 , =0 , =1

❸ ❸ ❸

❶ ❶ ❶

❷ ❷ ❷
12
 Espace d’approximation: cas des triangles ,
❸

❶ ❷
, ,
Dans chaque élément les 3 fonctions de forme
Vérifient:
Fonctions de forme

Nœud ❶

, = + + =1
, = + + =0
, = + + =0

Nœud ❷ Nœud ❸

, = + + =0 , = + + =0
, = + + =1 , = + + =0
, = + + =0 , = + + =1

13
 Espace d’approximation: cas des triangles
Dans chaque élément les 3 fonctions de forme vérifient:
Nœud ❶ Nœud ❷ Nœud ❸
, = + + =1 , = + + =0 , = + + =0
, = + + =0 , = + + =1 , = + + =0
, = + + =0 , = + + =0 , = + + =1
Fonctions de forme

D’où: D’où: D’où

1 1 1 0 1 0
1 = 0 1 = 1 1 = 0
1 0 1 0 1 1

1 1 0 0
1 = 0 1 0
1 0 0 1

=
( )
14
 Espace d’approximation: cas des triangles
Dans chaque élément les 3 fonctions de forme vérifient:
Nœud ❶ Nœud ❷ Nœud ❸
, = + + =1 , = + + =0 , = + + =0
, = + + =0 , = + + =1 , = + + =0
, = + + =0 , = + + =0 , = + + =1
Fonctions de forme

= ( )

https://fr.wikipedia.org/wiki/Aire_d%27un_triangle

det ( |=2 ≔| |

15
 Espace d’approximation: cas des quadrangles
15 16 17 18
9 Dans chaque élément les 4 fonctions de forme
12 Vérifient:
11 13 Nœud ❶

6 14 , = + + + =1
10 , = + + + =0
7 8
, = + + + =0
Fonctions de forme

4 , = + + + =0
1 2 3 5
❸ , Nœud ❷
, ❹
, = + + + =0
, = + + + =1
❷ , = + + + =0
❶ , , = + + + =0
,
Nœud ❸ Nœud ❹
, = + + + =0 , = + + + =0
, = + + + =0 , = + + + =0
, = + + + =1 , = + + + =0
, = + + + =0 , = + + + =1
16
 Espace d’approximation: cas des quadrangles
Dans chaque élément les 4 fonctions de forme vérifient:

Nœud ❶: est donnée par Nœud ❷: est donnée par Nœud ❸ : est donnée par

1 1 1 0 1 0
1 1 1
= 0 = 1 = 0
1 0 1 0 1 1
Fonctions de forme

1 0 1 0 1 0

Nœud ❹: est donnée par

1 1 0 0 0
1 0 1
1
= 0 = 0 1 0 0
1 1 1 0 0 1 0
1 0 1 0 0 0 1

1
1
=
1
1
17
 Quadratures de Gauss en 2D
→ Cas des triangles → Cas des quadrangles
t
t
−1,1 1,1
Élément de référence

0,1

s −1, −1 1, −1
0,0 1,0

= , = , = , = ,

≈ , ≈ ,

18
 Quadratures de Gauss en 2D
Exacte pour les polynomes de degré 1
→ Cas des triangles

t ≈
1 1 1
,
2 3 3
Élément de référence

0,1

Exacte pour les polynomes de degré 2

0,0
s 1 1 1 1 1 1 1
1,0 ≈ , + ,0 + 0,
6 2 2 6 2 6 2

= , = ,
Exacte pour les polynomes de degré 3

≈ ,
−27 1 1 25
≈ , + 0.2,0.6
96 3 3 96
25 25
+ 0.2,0.2 + 0.6,0.2
96 96
19
 Quadratures de Gauss en 2D
→ Cas des quadrangles
t

−1,1 1,1
Élément de référence

−1, −1 1, −1

= , = ,

≈ ,

20
Élément de référence Élément de référence : cas des triangles

21
 Élément de référence : cas des quadrangles
t

−1,1 1,1

s
Élément de référence

−1, −1 1, −1

22
 Fonctions de forme dans l’élément de référence
→ Cas des triangles → Cas des quadrangles t
t
−1,1 1,1
0,1
Élément de référence

0,0 1,0
s −1, −1 1, −1

23
Matrices élémentaires Matrices élémentaires

24
Procédure éléments finis

Technique d’assemblage

Séance de TP

25
Assemblage Assemblage

26
Procédure éléments finis

Système linéaire

Séance de TP

27
Extensions
• L’extension du cours aux cas des éléments
finis de degrés élevés est facile.
• La condition aux limites de type Neumann
est une intégrale qui s’ajoute au second
membre.

28