Cours 5
MEF en 1D
• Maillage: connectivité
• Élément de référence
• Technique d’assemblage
• Résolution et postraitement
1
Maillage avec plusieurs éléments
◼ Exemple : 3 éléments chacun à 2 noeuds
4 noeuds
◼ Un maillage est bien décrit par 2 tables:
… coordonnées : vcor = x1 x2 x3 x4
1 2
… connectivités : conec = 2 3
3 4
2
Remarques
Remarque 1: la numérotation des éléments peut-être quelconque
1 2
conec = 2 3
3 4
2 4
conec = 4 1
1 3
4 1
conec = 1 3
3 2
Remarque 2 : les éléments peuvent-être de tailles différentes
3
Évaluation de la forme intégrale
L
d ( x ) dT ( x ) L
W = + ( x) k dx − ( x ) f dx + ( L ) h (T ( L ) − Text ) − ( 0 ) q ( 0 ) = 0
0 dx dx 0
W
Wint CL
La décomposition du domaine et le choix de fonctions intégrables permet d’utiliser:
L N −1 xi +1
(...) dx = x (...) dx
0 i =1
with x1 = 0 and xN = L
i
element contribution
Ainsi :
nelt
W = Wint + WCL = Winte + WCL = 0, avec nelt = nombre d'élements.
e =1
4
Élément de référence
◼ Les fonctions de forme semblent dépenadantes des coordonnées
Besoin de les recalculer pour chaque élément
◼ En fait, avec la transformation géométrique vers l’élément de référence,
les fonctions de forme sont unique en fonction des coordonnées de
référence.
5
Reference element
( xi +1 − x)
◼ Changement de coordonnées x xi , xi +1 → 0, 1 , =
Le
avec Le = xi +1 − xi
xi +1 1
D’où : dx dx
x (...) dx = 0 (...) d d with
d
= Le
i
Nodes xi +1 nelt 1
(...) dx = (...) d
i =1 e =1 0
x i
◼ Les fonctions de forme linéaire en 1D sont:
N1 ( ) = 1 − , N 2 ( ) =
Vérifier que : N1 ( 0 ) = 1 N 2 ( 0 ) = 0
and
N1 (1) = 0 N 2 (1) = 1
6
Élément de réference
◼ On peut aussi choisir l’élément de référence: −1, 1
x xi , xi +1 → −1, 1 , x = N1 ( ) xi + N 2 ( ) xi +1
1− 1+
avec N1 ( ) = et N 2 ( ) =
2 2
ainsi: xi +1 1
dx dx Le
(...) dx = (...) d avec =
xi −1
d d 2
Les intégrales peuvent être aussi calculées numériquement avec
la quadrature de Gauss-Legendre:
Nodes −1 xi +1 nelt 1
Le
nelt NPG
Le
( f ( x) ) dx = ( f ( x( ) ) d = w j f ( x( j ))
2
i =1
x i
e =1 −1 e =1 j =1 2
7
Solution sur le domaine global : reconstruction
élément-par-élément
8
Calcul des intégrales sur chaque élément: exemple
Problème de diffusion avec K et f constants:
e
T1
T ( ) = N1 ( ) N 2 ( )
◼ Approximation par élément :
T2
e
1
( ) = N1 ( ) N 2 ( )
2
◼ Calcul des contributions de chaque élément :
e
e k 1 −1 T Le 1
Wint = 1 2
e
− 1 2
1
f
Le −1 1 T2 2 1
e
T
= 1 2
e
K e 1 − 1 2
e
F
e
T2
◼ [Ke ] : matrice élémentaire de rigidité
◼ {Fe } : vecteur élémentaire de charge
9
Assemblage
Comment mettre sous forme matricielle la contribution de chaque élément?
e
e 1
nelt T nelt
W = 1 2
e
K − 1 2
e
F + W
e
CL =0
e =1 T2 e =1
3. Assembler les matrices élémentaire dans une
matrice globale [K]
Idem pour le vecteurs de charge dans {F}
T1
W = 1 N K − F − R = 0
N N T
N N 1
Deux techniques sont possibles!
10
Assemblage avec un surdimensionnement (à
éviter!)
◼ Augmenter la dimension de chaque matrice élémentaire et remplir par
zéro
◼ Exemple :
1 −1 0 0 T1 1
k
−1 1 0 0 T2 L(1)
1
Wint = Wint1 + Wint2 + Wint3 = 1 2 3 4 (1) − f +
L 0 0 0 0 T3 2 0
0 0 0 0 T4 0
0 0 0 0 T1 0 0 0 0 0 T1 0
(3)
k 0 1 −1 0 T2 L 1 k 0 0 0 0 T2
(2)
L 0
− f + − f
L(2) 0 −1 1 0 T3 2 1 L(3) 0 0 1 −1 T3 2 1
0 0 0 0 T4 0 0 0 −1 1 T
4 1
11
Assemblage par un algorithme de localisation
Utikiser la table des connectivities
1 2
conec = 2 3 N° lgine = numéro de l’élément
3 4
Colonnes = numéros des nœuds de l’élément courant
12
Assemblage
◼ Boucle sur les éléments :
1. Calculer [Ke] et {Fe}
2. Extraire la table des connectivités de l’élément
3. Trouver la correspondance entre le numéro local et le numéro global
1. Ajouter [Ke] à sa place dans [K]
2. Ajouter {Fe} à sa place dans {F}
fin
Introduction des conditions aux limites de Dirichlet!
13
exemple:
1 2 3 4
1 −1 0 0 T1 1
1
Assemblage de l’élément 1:
k −1 1 0 0 T2 e
1
= f
◼ L
2
Conec(1,[1 2])=[1 2] Le 0 0 0 0 T3 2 0
3
0 0 0 0 T4 0
4
N° de l’élément Liste des noeuds
1 2 3 4
1 −1 0 0 T1 1
1
k −1 1 + 1 −1 0 T2 Le
1 + 1
élément 2 : = f
2
◼
Le 0 −1 1 0 T3 2 1
Conec(2,[1 2])=[2 3]
3
0 0 0 0 T4 0
4
1 2 3 4
1 −1 0 0 T1 1
1
k −1 2 −1 0 T2 Le 2
2
◼ élément 3 : = f
Le
0 −1 1 + 1 −1 T
3 2 1 + 1
Conec(3,[1 2])=[3 4]
3
0 0 −1 1 T4 1
4
Remarque : pour simplifier on a supposé L(1) = L (2) = L (3) = Le
14
Cas d’un maillage non structuré
◼ Exemple : conec(e, [1 2]) = [ 1 3]
k k
Le −
Le
K =
e
1 2 3 4
− k k
Le Le 1
2
K =
3
4
Technique : utilisation des the connectivités!
15
◼ Example : conec(e, [1 2]) = [ 4 2]
k k
Le −
Le
K =
e
1 2 3 4
− k k
Le Le 1
2
K =
3
4
16
C.L. de type Dirichlet
◼ exemple : T ( x = 0) = T = 30
1
Remplacer la ligne correspondante
1 0 0 0 T 30
1
k k k e
− e 2 − 0 T fL
L Le Le 2
=
k k k
0 − e 2 e − e T fL
e
L L L 3
k k L e
0 0 − e T f
L Le 4 2
Remarque : le vecteur des réaction disparait!
17
Autre approche
Ajouter un grand terme à la diagonale de la ligne correspondante
k k Le
Le + Big − e
L
0 0 T1 Big 30 + f
2
−k k k
2 e − 0 T2 fLe
Le L Le
=
k k k
0 − e 2 e − e T3 fLe
L L L
k k Le
0 0 − e T
e 4 f
L L
2
avec Big = 1012 x max(K)
18
C.L. de type Cauchy
◼ Sont déjà introduites dans la forme faible :
Wcl = ( L ) h (T ( L ) − Text ) = 4 hT4 − 4 hText
attention au
signe!
1 0 0 0 T1 30
k k k
− e 2 − 0 T fLe
L Le Le 2
=
k k k
0 − e 2 e − e T fL e
L L L 3
L
e
k k
0 0 − e + h T f + h Text
L Le 4 2
19
Résultats: vérification et validation
◼ Vérification: debuguer le code
◼ Vérifier les c.l. Drichlet
Vérifier les conditions de Cauch
compute the gradient
Calculer les réactions
vérifier:
▪ L’équilibre global
▪ Conservation globale de l’énergie
◼ Convergence avec le maillage:
vérifier l’independance
20
Post-traitement: calculer le the gradient
dT
flux : q = −k
dx
du
réaction : N = EA
dx
e e
e e T1 T1
Méthode : dT 1 dT 1 1 1 T2 − T1
◼ = e = e N1' ( ) N 2' ( ) = − e =
dx L d L T L Le Le
2 T2
Question : comment définir le gradient au nœud « i » ?
21
Calculer du gradient au noeud
0
◼C approximation entraine :
Continuité de la solution;
Discontinuité des dérivées!
◼ Comment lisser les discontinuités:
Utiliser une approximation par EF de type C1
❑ Ou par une moyenne pondérée: by Averaging :
i
q (e)
Ae
q i = e=i −1 i with i = 2,.., N − 1.
A
e =i −1
e
❑ Moyenne par les moindres carrés (Voir le cours 1)
22
Post-traitement : calcul des réactions
◼ Méthode : R =
K T solution − F
avant Dirichlet !
Réaction
On doit garantir :
F = R + P = 0
weight
Model réel Model EF
23
Convergence avec le maillage
◼ Objectif :
vérifier [Link] de convergence (qualité et efficacit): est-ce qu’il y a une
taille du maillage pour laquelle la solution numérique est satisfaifante?
24
Cas spécial : super convergence!
Solution de température convergée
Température
Flux
Flux non convergé !
La solution numérique T est exacte aux nœuds!
25