Le document décrit la procédure de calcul par éléments finis sur un problème modèle en deux dimensions. La procédure inclut la création d'un maillage, le calcul des matrices élémentaires et l'assemblage du système linéaire.
Le document décrit la procédure de calcul par éléments finis sur un problème modèle en deux dimensions. La procédure inclut la création d'un maillage, le calcul des matrices élémentaires et l'assemblage du système linéaire.
Le document décrit la procédure de calcul par éléments finis sur un problème modèle en deux dimensions. La procédure inclut la création d'un maillage, le calcul des matrices élémentaires et l'assemblage du système linéaire.
Le document décrit la procédure de calcul par éléments finis sur un problème modèle en deux dimensions. La procédure inclut la création d'un maillage, le calcul des matrices élémentaires et l'assemblage du système linéaire.
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
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
≈ , −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.
