Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
CORRIGE SERIE N°03 : Techniques d’interpolations nodales (suite)
Solution de l’exercice 01 :
On a une barre maillée avec deux éléments. Le premier élément est un élément 1D quadratique, par
contre le deuxième élément est un élément 1D linéaire.
Elément 1 Elément 2
1 2 3 4
X=0 X=3
1/ Détermination de l’interpolation (pour chaque élément) 𝑢(𝑥) en passant par l’élément de référence
𝑈(1) (𝑥) et 𝑈(2) (𝑥) ?
Commençons par l’élément quadratique 1D à 3 nœuds 𝑈(1) (𝑥)
3
𝑥(𝜉) = ∑ 𝑁𝑖 (𝜉) 𝑥𝑖 = 𝑁1 (𝜉) 𝑥1 + 𝑁2 (𝜉) 𝑥2 + 𝑁3 (𝜉)𝑥3
𝑖=1
𝜉 𝜉
𝑥(𝜉) = − (1 − 𝜉)0 + (1 − 𝜉 2 )1 + (1 + 𝜉)2
2 2
𝑥(𝜉) = 1 − 𝜉 2 + 𝜉 + 𝜉 2 = 1 + 𝜉
⇒ 𝑥(𝜉) = 1 + 𝜉 ⇒ 𝜉 =𝑥−1
𝜉 =𝑥−1
Et
3
𝑢(1) (𝑥) = ∑ 𝑁𝑖 (𝜉) 𝑢𝑖 = 𝑁1 (𝜉) 𝑢1 + 𝑁2 (𝜉)𝑢2 + 𝑁3 (𝜉)𝑢3
𝑖=1
𝜉 𝜉
𝑢(1) (𝑥) = − (1 − 𝜉)0 + (1 − 𝜉 2 )2 + (1 + 𝜉)4
2 2
𝑢(1) (𝑥) = 2 − 2𝜉 2 + 2𝜉 + 2𝜉 2 ⇒ 𝑢(1) (𝑥) = 2 + 2𝜉 = 2(1 + 𝜉)
1
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
𝑢(1) (𝑥) = 2(1 + 𝜉) or 𝜉 =𝑥−1
Donc
𝑢(1) (𝑥) = 2(1 + (𝑥 − 1))
𝑢(1) (𝑥) = 2𝑥
Pour le deuxième élément linéaire 1D à 2 nœuds 𝑈(2) (𝑥)
2
𝑥(2) (𝜉) = ∑ 𝑁𝑖 (𝜉) 𝑥𝑖 = 𝑁1 (𝜉) 𝑥1 + 𝑁2 (𝜉) 𝑥2
𝑖=1
Avec 𝑥1 = 2 et 𝑥2 = 3
(1 − 𝜉) (1 + 𝜉) 3 3𝜉 5 𝜉
𝑥(2) (𝜉) = 2+ 3=1−𝜉+ + = +
2 2 2 2 2 2
1
𝑥(2) (𝜉) = 2 (5 + 𝜉) ⇒ 𝜉 = 2𝑥 − 5
Et
2
𝑢(2) (𝜉) = ∑ 𝑁𝑖 (𝜉) 𝑢𝑖 = 𝑁1 (𝜉) 𝑢3 + 𝑁2 (𝜉)𝑢4
𝑖=1
1−𝜉 1+𝜉
𝑢(2) (𝑥) = 𝑢3 + 𝑢4
2 2
Or 𝑢3 = 4 et 𝑢4 = 0
Cela permet d’écrire :
1−𝜉
𝑢(2) (𝜉) = (4) = 2(1 − 𝜉)
2
Mais
𝜉 = 2𝑥 − 5
𝑢(2) (𝜉(𝑥)) = 𝑢(2) (𝑥) = 2(1 − (2𝑥 − 5)) = 2(6 − 2𝑥)
⇒ 𝑢(2) (𝑥) = 12 − 4𝑥
Pour 𝑥 = 1.5
𝑑𝑢1 (𝑥) 𝑑(2𝑥)
𝜉𝑥 = = =2
𝑑𝑥 𝑑𝑥
2
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
𝑑𝑢 𝑑𝑢
2/ Donnons l’expression de en fonction de (pour chaque élément) pour l’élément de référence
𝑑𝑥 𝑑𝜉
associé à ces types d’éléments.
Donc pour le premier élément :
𝑑𝑢 𝑑𝜉 𝑑𝑢 𝑑(𝑥 − 1) 𝑑(2 + 2𝜉)
=( ) =
𝑑𝑥 𝑑𝑥 𝑑𝜉 𝑑𝑥 𝑑𝜉
𝑑𝑢 𝑑(2𝑥) 𝑑(𝑥 − 1)
( = = 2) = 2
𝑑𝑥 𝑑𝑥 𝑑𝑥
𝑑𝑢 𝑑𝑢
=1
𝑑𝑥 𝑑𝜉
Pour le deuxième élément :
𝑑𝑢 𝑑𝜉 𝑑𝑢 𝑑(2𝑥 − 5) 𝑑𝑢
=( ) =
𝑑𝑥 𝑑𝑥 𝑑𝜉 𝑑𝑥 𝑑𝜉
𝑑𝑢 𝑑𝑢
⇒ =2
𝑑𝑥 𝑑𝜉
Solution de l’exercice 02 :
1/ Déterminons ‶la matrice Jacobienne″ associée à cet élément :
On a une forme particulière de la matrice Jacobienne en EF, cas du TRI 3 avec nœud 1 (𝑥1 , 𝑦1 ),
nœud 2 (𝑥2 , 𝑦2 ) et nœud 3 (𝑥3 , 𝑦3 ).
𝜕𝑥 𝜕𝑦
𝜕𝜉 𝜕𝜉
𝐽 = [𝜕𝑥 𝜕𝑦
]
𝜕𝜂 𝜕𝜂
Avec
3
𝑥(𝜉, 𝜂) = ∑ 𝑁𝑖 (𝜉, 𝜂) 𝑥𝑖
𝑖=1
3
𝑦(𝜉, 𝜂) = ∑ 𝑁𝑖 (𝜉, 𝜂) 𝑦𝑖
𝑖=1
3
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
Remplaçons 𝑥(𝜉, 𝜂) et 𝑦(𝜉, 𝜂) dans 𝐽
On aura :
𝜕𝑁1 𝜕𝑁2 𝜕𝑁3
𝜕𝜉 𝜕𝜉 𝜕𝜉
𝐽= [{𝑥𝑖 }{𝑦𝑖 }]
𝜕𝑁1 𝜕𝑁2 𝜕𝑁3
[ 𝜕𝜂 𝜕𝜂 𝜕𝜂 ]
𝜕𝑁1 𝜕𝑁2 𝜕𝑁3
𝑥 𝑦1
𝜕𝜉 𝜕𝜉 𝜕𝜉 𝑥1
𝐽= [ 2 𝑦2 ]
𝜕𝑁1 𝜕𝑁2 𝜕𝑁3 𝑥 𝑦3
3
[ 𝜕𝜂 𝜕𝜂 𝜕𝜂 ]
𝜕𝑁1 𝜕𝑁2 𝜕𝑁3 𝜕𝑁1 𝜕𝑁2 𝜕𝑁3
( 𝑥1 + 𝑥2 + 𝑥3 ) ( 𝑦1 + 𝑦2 + 𝑦 )
𝜕𝜉 𝜕𝜉 𝜕𝜉 𝜕𝜉 𝜕𝜉 𝜕𝜉 3
𝐽=
𝜕𝑁1 𝜕𝑁2 𝜕𝑁3 𝜕𝑁1 𝜕𝑁2 𝜕𝑁3
( 𝑥1 + 𝑥2 + 𝑥3 ) ( 𝑦1 + 𝑦2 + 𝑦 )
[ 𝜕𝜂 𝜕𝜂 𝜕𝜂 𝜕𝜂 𝜕𝜂 𝜕𝜂 3 ]
Avec
𝑁1 (𝜉, 𝜂) = 1 − 𝜉 − 𝜂 , 𝑁2 (𝜉, 𝜂) = 𝜉 , 𝑁3 (𝜉, 𝜂) = 𝜂
Et
𝜕𝑁1 𝜕𝑁2 𝜕𝑁3
= −1, = 1, =0
𝜕𝜉 𝜕𝜉 𝜕𝜉
𝜕𝑁1 𝜕𝑁2 𝜕𝑁3
= −1, = 0, =1
𝜕𝜂 𝜕𝜂 𝜕𝜂
Donc
−𝑥 + 𝑥 −𝑦1 + 𝑦2
𝐽 = [−𝑥1 + 𝑥2 −𝑦1 + 𝑦3 ]
1 3
𝑥2 − 𝑥1 𝑦2 − 𝑦1
𝐽 = [𝑥 − 𝑥 𝑦3 − 𝑦1 ]
3 1
Calculons le déterminant de 𝐽
det 𝐽 = (𝑥2 − 𝑥1 )(𝑦3 − 𝑦1 ) − (𝑥3 − 𝑥1 )(𝑦2 − 𝑦1 )
1 𝑦3 − 𝑦1 𝑥1 − 𝑥3 𝑇
𝑗 = [𝐽]−1 = ( 𝑥2 − 𝑥1 )
det 𝐽 𝑦1 − 𝑦2
4
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
1 𝑦3 − 𝑦1 𝑦1 − 𝑦2
𝑗= ( 𝑥2 − 𝑥1 )
det 𝐽 𝑥1 − 𝑥3
Solution de l’exercice 03 :
On a :
Nœud 1 (0, 0), nœud 2 (2, 0), nœud 3 (2, 1), nœud 4 (0, 1)
Et
𝑇1 = 20, 𝑇2 = 21, 𝑇3 = 18, 𝑇4 = 15
Revenons à l’élément de référence
4
𝑥(𝜉, 𝜂) = ∑ 𝑁𝑖 (𝜉, 𝜂) 𝑥𝑖
𝑖=1
Et
4
𝑦(𝜉, 𝜂) = ∑ 𝑁𝑖 (𝜉, 𝜂) 𝑦𝑖
𝑖=1
Avec les fonctions d’interpolation d’un élément quadrangle linéaire suivantes :
1
𝑁1 (𝜉, 𝜂) = 4 (1 − 𝜉)(1 − 𝜂)
1
𝑁2 (𝜉, 𝜂) = 4 (1 + 𝜉)(1 − 𝜂)
1
𝑁3 (𝜉, 𝜂) = 4 (1 + 𝜉)(1 + 𝜂)
1
𝑁4 (𝜉, 𝜂) = (1 − 𝜉)(1 + 𝜂)
4
𝑥(𝜉, 𝜂) = 𝑁1 (𝜉, 𝜂)0 + 𝑁2 (𝜉, 𝜂)2 + 𝑁3 (𝜉, 𝜂)2 + 𝑁4 (𝜉, 𝜂)0
2
𝑥(𝜉, 𝜂) = 2(𝑁2 (𝜉, 𝜂) + 𝑁3 (𝜉, 𝜂)) = (1 + 𝜉)((1 − 𝜂) + (1 + 𝜂))
4
𝑥(𝜉, 𝜂) = (1 + 𝜉) ⇒ 𝜉 =𝑥−1
Et
𝑦(𝜉, 𝜂) = 𝑁1 (𝜉, 𝜂)0 + 𝑁2 (𝜉, 𝜂)0 + 𝑁3 (𝜉, 𝜂)1 + 𝑁4 (𝜉, 𝜂)1
1
𝑦(𝜉, 𝜂) = 𝑁3 (𝜉, 𝜂) + 𝑁4 (𝜉, 𝜂) = (1 + 𝜂)((1 + 𝜉) + (1 − 𝜉))
4
5
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
1
⇒ 𝑦(𝜉, 𝜂) = (1 + 𝜂) ⇒ 𝜂 = 2𝑦 − 1
2
Or dans l’élément de référence
𝑇(𝜉, 𝜂) = 𝑁1 (𝜉, 𝜂)𝑇1 + 𝑁2 (𝜉, 𝜂)𝑇2 + 𝑁3 (𝜉, 𝜂)𝑇3 + 𝑁4 (𝜉, 𝜂)𝑇4
𝑇(𝜉, 𝜂) = 𝑁1 (𝜉, 𝜂)20 + 𝑁2 (𝜉, 𝜂)21 + 𝑁3 (𝜉, 𝜂)18 + 𝑁4 (𝜉, 𝜂)15
20 21 18 15
𝑇(𝜉, 𝜂) = (1 − 𝜉)(1 − 𝜂) + (1 + 𝜉)(1 − 𝜂) + (1 + 𝜉)(1 + 𝜂) + (1 − 𝜉)(1 + 𝜂)
4 4 4 4
21 9 15
𝑇(𝜉, 𝜂) = 5(1 − 𝜉)(1 − 𝜂) + (1 + 𝜉)(1 − 𝜂) + (1 + 𝜉)(1 + 𝜂) + (1 − 𝜉)(1 + 𝜂)
4 2 4
21 9 15
𝑇(𝜉, 𝜂) = (1 − 𝜂) [5 (1 − 𝜉) + (1 + 𝜉)] + (1 + 𝜂) [ (1 + 𝜉) + (1 − 𝜉)]
4 2 4
20 20 21 21 18 18 15 15
𝑇(𝜉, 𝜂) = (1 − 𝜂) [ − 𝜉+ + 𝜉] + (1 + 𝜂) [ + 𝜉 + − 𝜉]
4 4 4 4 4 4 4 4
41 1 33 3
𝑇(𝜉, 𝜂) = (1 − 𝜂) [ + 𝜉] + (1 + 𝜂) [ + 𝜉]
4 4 4 4
41 1 41 1 33 3 33 3
𝑇(𝜉, 𝜂) = + 𝜉− 𝜂− 𝜉𝜂+ + 𝜉+ 𝜂+ 𝜉𝜂
4 4 4 4 4 4 4 4
74 4 8 2
𝑇(𝜉, 𝜂) = + 𝜉− 𝜂+ 𝜉𝜂
4 4 4 4
37 1
𝑇(𝜉, 𝜂) = + 𝜉 − 2𝜂 + 𝜉 𝜂
2 2
Or pour
𝜉 =𝑥−1 et 𝜂 = 2𝑦 − 1
37 1
𝑇(𝜉(𝑥, 𝑦), 𝜂(𝑥, 𝑦)) = 𝑇(𝑥, 𝑦) = − 2(2𝑦 − 1) + (2𝑥𝑦 − 2𝑦 − 𝑥 + 1) + (𝑥 − 1)
2 2
37 1 1
𝑇(𝑥, 𝑦) = − 4𝑦 + 2 + 𝑥 − 1 + 𝑦𝑥 − 𝑥 − 𝑦 +
2 2 2
1
𝑇(𝑥, 𝑦) = 𝑥 − 5𝑦 + 𝑦𝑥 + 20
2
1
𝑇(𝑥, 𝑦) = 𝑥 − 5𝑦 + 𝑦𝑥 + 20
2
2/ Le vecteur gradient de la température aura pour composantes :
6
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
𝜕𝑇(𝑥, 𝑦)
1
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑇(𝑥, 𝑦) = 𝜕𝑥 +𝑦 )
𝑔𝑟𝑎𝑑 𝜕𝑇(𝑥, 𝑦) = ( 2
−5 + 𝑥
( 𝜕𝑦 )
3/Calculons et dessinons la courbe isovaleur correspondante.
A 𝑇 = 20 C’est-à-dire, trouvons les 𝑇(𝑥, 𝑦) telle que 𝑇(𝑥, 𝑦) = 20
1 1
𝑥 − 5𝑦 + 𝑦𝑥 + 20 = 20 ⇒ 𝑥 − 5𝑦 + 𝑦𝑥 = 0
2 2
𝑥
⇒ + (𝑥 − 5)𝑦 = 0
2
𝑥
⇒ = (5 − 𝑥)𝑦
2
𝑥
𝑦=
2(5 − 𝑥)
0 X
Solution de l’exercice 04 :
On a deux triangles linéaires :
L’élément 1 avec les nœuds 1, 2 et 4.
L’élément 2 avec les nœuds 2, 3 et 4.
7
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
(0, 1) (2, 1)
4 3
2
1
(0, 0) (2, 0) X
1/ Donc calculons en passant par l’élément de référence l’expression de l’interpolation 𝑇(1) (𝑥, 𝑦) et
𝑇(2) (𝑥, 𝑦).
Pour l’élément (1) :
3
𝑥(𝜉, 𝜂) = ∑ 𝑁𝑖 (𝜉, 𝜂) 𝑥𝑖
𝑖=1
Et
3
𝑦(𝜉, 𝜂) = ∑ 𝑁𝑖 (𝜉, 𝜂) 𝑦𝑖
𝑖=1
Avec les fonctions d’interpolation d’un élément de type triangle linéaire suivantes :
𝑁1 (𝜉, 𝜂) = 1 − 𝜉 − 𝜂
𝑁2 (𝜉, 𝜂) = 𝜉
𝑁3 (𝜉, 𝜂) = 𝜂
𝑥(𝜉, 𝜂) = (1 − 𝜉 − 𝜂) × 0 + 2 × 𝜉 + 𝜂 × 0
1
⇒ 𝑥(𝜉, 𝜂) = 2 𝜉 ⇒ 𝜉 = 2𝑥
Et pour 𝑦(𝜉, 𝜂) :
𝑦(𝜉, 𝜂) = (1 − 𝜉 − 𝜂) × 0 + 0 × 𝜉 + 𝜂 × 1
⇒ 𝑦(𝜉, 𝜂) = 𝜂 ⇒ 𝜂=𝑦
Or
𝑇(1) (𝜉, 𝜂) = (1 − 𝜉 − 𝜂) 𝑇1 + 𝜉 𝑇2 + 𝜂 𝑇4 = (1 − 𝜉 − 𝜂) 20 + 𝜉 21 + 𝜂 15
8
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
⇒ 𝑇(1) (𝜉, 𝜂) = 20 − 20 𝜉 − 20 𝜂 + 21 𝜉 + 15 𝜂 = 20 + 𝜉 − 5 𝜂
𝑥
⇒ 𝑇(1) (𝜉, 𝜂) = 20 + 𝜉 − 5 𝜂 où 𝜉=2 et 𝜂=𝑦
Donc
𝑥
𝑇(1) (𝜉(𝑥, 𝑦), 𝜂(𝑥, 𝑦)) = 20 + − 5𝑦
2
𝑥
⇒ 𝑇(1) (𝑥, 𝑦) = 20 + − 5𝑦
2
Pour l’élément (2) :
3
𝑥(𝜉, 𝜂) = ∑ 𝑁𝑖 (𝜉, 𝜂) 𝑥𝑖
𝑖=1
Et
3
𝑦(𝜉, 𝜂) = ∑ 𝑁𝑖 (𝜉, 𝜂) 𝑦𝑖
𝑖=1
Donc
𝑥(𝜉, 𝜂) = (1 − 𝜉 − 𝜂) × 2 + 2 × 𝜉 + 𝜂 × 0
⇒ 𝑥(𝜉, 𝜂) = 2 − 2 𝜉 − 2 𝜂 + 2 𝜉
𝑥
⇒ 𝑥(𝜉, 𝜂) = 2(1 − 𝜂) ⇒ 𝜂 = −2 +1
Et pour 𝑦(𝜉, 𝜂) :
𝑦(𝜉, 𝜂) = (1 − 𝜉 − 𝜂) × 0 + 𝜉 + 𝜂
𝑥 𝑥
𝑦(𝜉, 𝜂) = 𝜉 + 𝜂 mais on a déjà 𝜂 =1−2 donc 𝜉 =𝑦+2−1
Or
𝑇(2) (𝜉, 𝜂) = (1 − 𝜉 − 𝜂) 𝑇2 + 𝜉 𝑇3 + 𝜂 𝑇4 = (1 − 𝜉 − 𝜂) 21 + 𝜉 18 + 𝜂 15
⇒ 𝑇(2) (𝜉, 𝜂) = 21 − 21 𝜉 − 21 𝜂 + 18 𝜉 + 15 𝜂
⇒ 𝑇(2) (𝜉, 𝜂) = 21 − 3 𝜉 − 6 𝜂
On remplace 𝜉 et 𝜂 par leurs expressions respectives
𝑥 𝑥
𝑇(2) (𝜉(𝑥, 𝑦), 𝜂(𝑥, 𝑦)) = 𝑇(2) (𝑥, 𝑦) = 21 − 3 (𝑦 + − 1) − 6(1 − )
2 2
3
𝑇(2) (𝑥, 𝑦) = 21 − 3𝑦 − 𝑥 − 3 − 6 − 3𝑥
2
9
Université A. Mira Bejaia. Niveau : Première année Master.
Faculté de Technologie. Option :
Département de Génie mécanique. Matière : Méthode des Eléments Finis.
3
𝑇(2) (𝑥, 𝑦) = 18 + 𝑥 − 3𝑦
2
2/ Calculons le gradient de température.
𝜕𝑇(1)
1
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝜕𝑥
(𝑥,
𝑔𝑟𝑎𝑑 𝑇(1) 𝑦) = 𝜕𝑇 =( 2 )
(1)
−5
( 𝜕𝑦 )
𝜕𝑇(2)
3
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝜕𝑥
𝑔𝑟𝑎𝑑 𝑇(2) (𝑥, 𝑦) = 𝜕𝑇 =( 2 )
(2)
−3
( 𝜕𝑦 )
On peut remarquer que le gradient n’est pas situé à la frontière entre les deux éléments.
3/ Concernant l’équation de l’isovaleur 𝑇 = 20 c’est-à-dire l’ensemble des points (x, y) telles que :
Sur l’équation (1) 𝑇(1) (𝑥, 𝑦) = 20 :
𝑥 𝑥
𝑇(1) (𝑥, 𝑦) = 20 + − 5𝑦 = 20 ⇒ − 5𝑦 = 0
2 2
𝑥
⇒ 𝑦(1) = 10
Sur l’équation (2) 𝑇(2) (𝑥, 𝑦) = 20 :
3
𝑇(2) (𝑥, 𝑦) = 18 + 𝑥 − 3𝑦 = 20
2
3
𝑇(2) (𝑥, 𝑦) = 𝑥 − 2 = 3𝑦
2
2 1
⇒ 𝑦(2) = − 3 + 2 𝑥
On peut remarquer que le champ de température est bien situé à la frontière.
10