Exercice 2 : On se propose de résoudre l’équation différentielle suivante par la méthode des
éléments finis :
𝑑2 𝑢
− 𝑑𝑥 2 + 𝑢 = 𝑥 𝑝𝑜𝑢𝑟 𝑥 ∈ ]0,1[ (1)
où u est l’inconnue du problème. Les conditions aux limites associées à l’équation (1) sont définies
par :
𝑢(0) = 0
{ (2)
𝑢(1) = 0
La solution exacte de ce problème est donnée par :
sinh (𝑥)
𝑢(𝑥) = 𝑥 −
sinh (1)
1°/ Montrer que la formulation variationnelle du problème (1) s’écrit sous la forme :
1 1
𝑑𝑢 𝑑𝛿𝑢
∫ ( + 𝑢𝛿𝑢) 𝑑𝑥 = ∫ 𝑥𝛿𝑢𝑑𝑥
0 𝑑𝑥 𝑑𝑥 0
Réponse :
En multipliant l’équation (1) par 𝛿𝑢 et en intégrant sur l’intervalle ]0,1[, on trouve :
1 1
𝑑2𝑢
∫ (− + 𝑢) 𝛿𝑢𝑑𝑥 = ∫ 𝑥𝛿𝑢𝑑𝑥
0 𝑑𝑥 2 0
Ensuite, en appliquant la méthode d’intégration par partie au premier terme et en tenant compte
des conditions aux limites, on obtient :
1 1
𝑑𝑢 𝑑𝛿𝑢
∫ ( + 𝑢𝛿𝑢) 𝑑𝑥 = ∫ 𝑥𝛿𝑢𝑑𝑥
0 𝑑𝑥 𝑑𝑥 0
2/ On fait la même démarche que dans la question 2/ de l’exercice 1, écrire la formulation
variationnelle faible dans un élément 𝑖 et donner la matrice élémentaire [𝐾 𝑖 ] et le vecteur
élémentaire second membre {𝐹 𝑖 }.
Réponse :
Soit un élément fini 𝑖 limité par les nœuds 𝑖 et 𝑖 + 1 de coordonnées respectivement 𝑥𝑖 et 𝑥𝑖+1.
Alors la formulation variationnelle faible dans l’élément 𝑖 est donnée par :
𝑥𝑖+1 𝑥𝑖+1
𝑖
𝑑𝑢𝑖 𝑑𝛿𝑢𝑖
𝑊 =∫ ( + 𝑢 𝛿𝑢 ) 𝑑𝑥 − ∫ 𝑥𝛿𝑢𝑖 𝑑𝑥
𝑖 𝑖
𝑥𝑖 𝑑𝑥 𝑑𝑥 𝑥𝑖
𝑢𝑖
En introduisant l’approximation 𝑢𝑖 (𝑥) = 𝑁1 (𝑥)𝑢𝑖 + 𝑁2 (𝑥)𝑢𝑖+1 =< 𝑁1 (𝑥), 𝑁2 (𝑥) > {𝑢 } dans 𝑊 𝑖 , on
𝑖+1
obtient :
𝑊 𝑖 =< 𝛿𝑈 𝑖 > ([𝐾 𝑖 ]{𝑈 𝑖 } − {𝐹 𝑖 })
où [𝐾 𝑖 ] et {𝐹 𝑖 } sont donnés par :
Page 1 sur 6
𝑑𝑁1 2 𝑑𝑁1 𝑑𝑁2 3 + (𝑥𝑖+1 − 𝑥𝑖 )2 −6 + (𝑥𝑖+1 − 𝑥𝑖 )2
𝑥𝑖+1 ( ) 2
𝑖
[𝐾 ] = ∫ 𝑑𝑥 𝑑𝑥 𝑑𝑥 + [(𝑁1 ) 𝑁1 𝑁2
] 𝑑𝑥 =
3(𝑥𝑖+1 − 𝑥𝑖 ) 6(𝑥𝑖+1 − 𝑥𝑖 )
𝑥𝑖 𝑑𝑁1 𝑑𝑁2 𝑑𝑁2 2 𝑁1 𝑁2 (𝑁2 ) 2
−6 + (𝑥𝑖+1 − 𝑥𝑖 )2 3 + (𝑥𝑖+1 − 𝑥𝑖 )2
( ) ]
([ 𝑑𝑥 𝑑𝑥 𝑑𝑥 ) [ 6(𝑥𝑖+1 − 𝑥𝑖 ) 3(𝑥𝑖+1 − 𝑥𝑖 ) ]
1 2 1 2 1
𝑥𝑖+1 − 𝑥 + 𝑥 + 𝑥𝑥
{𝐹 𝑖 } = ∫
𝑁
𝑥 { 1 } 𝑑𝑥 = { 3
𝑖
6 𝑖+1
6 𝑖 𝑖+1 }
𝑁2 1 1 2 1
𝑥𝑖 − 𝑥𝑖2 + 𝑥𝑖+2 − 𝑥𝑖 𝑥𝑖+1
6 3 6
3/ D’après la technique d’assemblage, calculer la solution nodale dans le cas de deux et trois
éléments finis.
Réponse :
Discrétisation en deux éléments finis :
1
Pour l’élément 𝑖 = 1, 𝑥1 = 0 et 𝑥2 = 2, la matrice élémentaire [𝐾 1 ] et le second membre élémentaire
{𝐹1 } sont donnés par :
13 23
−
[𝐾 1 ] = [ 6 12]
23 13
−
12 6
1
{𝐹1 } = {24}
1
12
1
Pour l’élément 𝑖 = 2, 𝑥2 = 2 et 𝑥3 = 1, la matrice élémentaire [𝐾 2 ] et le second membre élémentaire
{𝐹 2 } sont donnés par :
13 23
−
[𝐾 2 ] = [ 6 12]
23 13
−
12 6
1
{𝐹 2 } = { 6 }
5
24
Les matrices et les seconds membres augmentés correspondants sont :
Page 2 sur 6
13 23
− 0
6 12
[𝐾 1𝑔 ] = 23 13
− 0
12 6
[ 0 0 0 ]
1
24
{𝐹1𝑔 } = 1
12
{0}
0 0
0 13 23
0 −
[𝐾 2𝑔 ] = 6 12
23 13
0 −
[ 12 6 ]
0
1
2𝑔
{𝐹 } = 6
5
{24}
Alors la matrice globale et le second membre global sont donnés par :
13 23
− 0
6 12
23 13 23
[𝐾] = [𝐾 1𝑎 ] + [𝐾 2𝑎 ] = − −
12 3 12
23 13
[ 0 −
12 6 ]
1
24
1
{𝐹} = {𝐹1𝑎 } + {𝐹 2𝑎 } =
4
5
{24}
Le système vérifié par le vecteur nodal {𝑈} qui tient compte des conditions aux limites est :
1 0 0 𝑈1 0
23 13 23 1
[− − ] {𝑈2 } = { }
12 3 12 𝑈 4
3
0 0 1 0
Les solutions sont : 𝑈1 = 0, 𝑈2 = 0.05769230769 et 𝑈3 = 0.
Discrétisation en trois éléments finis :
1
Pour l’élément 𝑖 = 1, 𝑥1 = 0 et 𝑥2 = 3, la matrice élémentaire [𝐾 1 ] et le second membre élémentaire
{𝐹1 } sont donnés par :
Page 3 sur 6
28 53
−
[𝐾 1 ] = [ 9 18]
53 28
−
18 9
1
{𝐹1 } = {54}
1
27
1 2
Pour l’élément 𝑖 = 2, 𝑥2 = 3 et 𝑥3 = 3, la matrice élémentaire [𝐾 2 ] et le second membre élémentaire
{𝐹 2 } sont donnés par :
28 53
−
[𝐾 2 ] = [ 9 18]
53 28
−
18 9
2
{𝐹 2 } = {27}
5
54
2
Pour l’élément 𝑖 = 3, 𝑥3 = 3 et 𝑥4 = 1, la matrice élémentaire [𝐾 3 ] et le second membre élémentaire
{𝐹 3 } sont donnés par :
28 53
−
[𝐾 3 ] = [ 9 18]
53 28
−
18 9
7
{𝐹 3 } = {54}
4
27
Les matrices et les seconds membres augmentés correspondants sont :
28 53
− 0 0
9 18
[𝐾 ] = − 53
1𝑔 28
0 0
18 9
0 0 0 0
[ 0 0 0 0 ]
Page 4 sur 6
1
54
{𝐹 } = 1
1𝑔
27
0
{0}
0 0 0
0 28 53
0 − 0
[𝐾 2𝑔 ] = 9 18
53 28
0 − 0
18 9 0
[0 0 0 ]
0
2
{𝐹 2𝑔 } = 27
5
54
{0}
0 0 0 0
0 0 0 0
28 53
[𝐾 3𝑔 ] = −
0 0 9 18
0 0 53 28
[ −
18 9 ]
0
0
7
{𝐹 3𝑔 } =
54
4
{27}
Alors la matrice globale et le second membre global sont donnés par :
28 53
− 0 0
9 18
56 53
−53/18 − 0
[𝐾] = [𝐾 1𝑎 ] + [𝐾 2𝑎 ] + [𝐾 3𝑎 ] = 9 18
56 53
53 −
0 − 9 18
0 18 53 28
0 −
[ 18 9 ]
Page 5 sur 6
1
54
1
{𝐹} = {𝐹1𝑎 } + {𝐹 2𝑎 } + {𝐹 2𝑎 } = 9
2
9
4
{27}
Le système vérifié par le vecteur nodal {𝑈} qui tient compte des conditions aux limites est :
1 0 0 0 0
53 56 53 𝑈1 1
− − 0 𝑈
18 9 18 { 2} = 9
53 56 53 𝑈3 2
0 − −
18 9 18 𝑈4 9
[ 0 0 0 1 ] {0}
Les solutions sont : 𝑈1 = 0, 𝑈2 = 0.04478685157, 𝑈3 = 0.05690806369 et 𝑈4 = 0.
4/ Comparer la solution obtenue par la méthode des éléments finis à celle analytique.
La solution analytique de l’équation (1) est donnée par :
1
𝑢(𝑥) = (𝑒 𝑥+1 − 𝑒 1−𝑥 ) + 𝑥(1 − 𝑒 2 )
1 − 𝑒2
avec 𝑒 = 𝑒𝑥𝑝(1).
Pour le cas de deux éléments finis :
𝑢(0) = 0
{ 𝑢(0.5) = 0.0565905579
𝑢(1) = 0
On remarque que les solutions obtenues par la méthode des éléments finis sont proches de celles
analytiques
Pour le cas de trois éléments finis :
𝑢(0) = 0
1
𝑢 ( ) = 0.0444121179
3
2
𝑢 ( ) = 0.0564235314
3
{ 𝑢(1) = 0
On remarque que les solutions obtenues par la méthode des éléments finis sont plus proches de
celles analytiques que dans le cas de deux éléments finis.
Page 6 sur 6
