Université A.

Mira, Béjaïa, Année Universitaire 2020-2021

Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI

D ELEMENTS FINIS
(SERIE N°2 (ELEMENT BARRE))

Exercice 1
Le système représenté sur la figure ci-après montre une barre de longueur 𝐿 et de section (a x b),
la barre est accrochée au point B par une articulation et elle est soumise à son poids propre.

Figure 1. Barre suspendue soumise à son poids propre

La figure 2, montre trois maillages différents de la barre respectivement avec un, deux et
trois éléments finis.
x=0

x=L

Maillage a Maillage b Maillage c

(1 élément fini) (2 éléments finis) (3 éléments finis)
Figure 2. Modélisation de la barre avec 1, 2 et 3 éléments finis
1. Si on considère uniquement un seul degré de liberté par nœud qui est le déplacement
axial vers le bas, calculer manuellement en utilisant la méthode des éléments finis
pour chaque type de maillage (a, b et c) le déplacement de chaque nœud.
2. Pour chaque type de maillage, calculer le déplacement de la barre en utilisant la
formule analytique et comparer les résultats déterminés par la méthode des éléments
finis.
3. Tracer les déplacements calculés par le maillage (a, b et c) et donner votre conclusion.
On donne :
La solution analytique : avec W represente le poids total de la barre.
𝐸𝐴 1 −1
Matrice de rigidité d’une barre avec un DDL par Nœud : 𝐾𝑒 = [
𝐿 −1
]
1
𝜌 = 7800 𝐾𝑔/𝑚3 𝐸 = 2,1 105 𝑀𝑝𝑎, A = 0.04 m², L = 1 m , g = 10 m/s².

1
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI

Exercice 2
Le système représenté sur la figure ci-après est constitué de trois barres assemblées avec des
articulations. La section des barres est A  2500 mm² et elles ont un module d’élasticité E 
210 000 MPa. Le système est posé sur deux appuis (l’un est simple, l’autre est double) et
supporte une charge horizontale F  65 KN.
y
3
F

Elè 02 4m
Elè 03

1 Elè 01 2 x

2m 1m

1) Calculer la matrice de rigidité de chacune des trois barres.

2) Assembler la matrice de rigidité globale.
3) Calculer le vecteur force global.
4) Appliquer les conditions aux limites et écrire le système d’équations réduit.
5) Résoudre le système réduit et donner les déplacements aux nœuds.
6) Calculer l’effort interne dans labarreN° 03 et donner la réaction au nœud 1.

On donne :
𝐶2 𝐶𝑆 −𝐶 2 −𝐶𝑆
𝐴𝐸 𝑆2 −𝐶𝑆 −𝑆 2
𝐾𝑏 = [ ]
𝐿 𝐶2 𝐶𝑆
2
𝑆𝑦𝑚 𝑆

2
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI

Solution Exercice 01
On a une barre suspendue soumise à son poids propre avec 𝜌 = 7800 𝐾𝑔/𝑚3 𝐸 = 2,1 105 𝑀𝑝𝑎, A =
0.04 m², L = 1 m , g = 10 m/s².

Figure 3. Barre suspendue soumise à son poids propre

1. Calcul des déplacements de la barre en utilisant les trois types de maillages (a, b et c) :

 Pour le calcul des déplacements on doit résoudre le système suivant : K U = F (K

représente la matrice de rigidité de l’élément, F le vecteur force et U le vecteur
déplacement).
 Si on considère uniquement un seul degré de liberté par nœud, alors, la matrice de
𝐸𝐴 1 −1
rigidité d’un élément élémentaire s’écrit : 𝐾𝑒 = 𝐿 [ ]
−1 1
a. Maillage a (un élément)

Le système à résoudre est le suivant : K U = F avec :

𝐸𝐴 1 −1 210000000∗0.04 1 −1 8400000 −8400000

𝐾 = 𝐾𝑒 = [
𝐿 −1
]= [ ]=[ ] 𝐾𝑁/𝑚
1 1 −1 1 −8400000 8400000

1/2
𝐹 = 𝐹𝑒 = 𝜌𝑔𝐴𝐿𝑒 [ ] 𝐿𝑒 : Représente la longueur élémentaire d’un élément fini.
1/2
1
1.56
𝐹 = 7.8 ∗ 10 ∗ 0.04 ∗ 1 [2] = [ ] 𝐾𝑁
1 1.56
2
Ce qui donne le système suivant :
8400000 −8400000 𝑈1 1.56
[ ][ ] = [ ]
−8400000 8400000 𝑈2 1.56
Et comme le nœud 1 est articulé alors U1 = 0, donc après application des conditions aux limites on
trouve le système réduit suivant :
8400000 ∗ 𝑈2 = 1.56
Donc : 𝑈2 = 1.8571e-07 m avec 𝑈1 = 0 m

3
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI

b. Maillage b (deux éléments)

Elément 01 :
Matrice de rigidité élémentaire de l’élément 1 :
Elé 1
210000000 ∗ 0.04 1 −1 16800000 −16800000
𝐾𝑒1 = [ ]=[ ] 𝐾𝑁/𝑚
𝐿𝑒 = 0.5 −1 1 −16800000 16800000
vecteur force élémentaire de l’élément 1 :
Elé 1
1
0.7800
𝐹𝑒1 = 7.8 ∗ 10 ∗ 0.04 ∗ 0.5 [2] = [ ] 𝐾𝑁
1 0.7800
2
Elément 02 :
Matrice de rigidité élémentaire de l’élément 2 :
210000000 ∗ 0.04 1 −1 16800000 −16800000
𝐾𝑒2 = [ ]=[ ] 𝐾𝑁/𝑚
𝐿𝑒 = 0.5 −1 1 −16800000 16800000
vecteur force élémentaire de l’élément 2 :
1
0.7800
𝐹𝑒2 = 7.8 ∗ 10 ∗ 0.04 ∗ 0.5 [2] = [ ] 𝐾𝑁
1 0.7800
2
Assemblage des deux éléments 1 et 2
16800000 −16800000 0 𝐾𝑁 0.7800
𝐾 =[ 33600000 −16800000] 𝐹 = [1.5600] 𝐾𝑁
𝑆𝑦𝑚 16800000 𝑚
0.7800

Le système global s’écrit :

16800000 −16800000 0 𝑈1 0.7800
[ 33600000 −16800000] [𝑈2] = [1.5600]
𝑆𝑦𝑚 16800000 𝑈3 0.7800
Apres application des CAL on trouve :

33600000 −16800000 𝑈2 1.5600

[ 𝑆𝑦𝑚 ][ ] = [
16800000 𝑈3
]
0.7800

Et après résolution du système :

𝑈2 = 1.393e-07 m , 𝑈3 = 1.8571e-07 m avec 𝑈1 = 0 m

c. Maillage 3 (trois éléments)

Elément 01 :
Elé 1
Matrice de rigidité élémentaire de l’élément 1 :
210000000 ∗ 0.04 1 −1 25200000 −25200000
𝐾𝑒1 = [ ]=[ ] 𝐾𝑁/𝑚
𝐿𝑒 = 0.333 −1 1 −25200000 25200000 Elé 2

vecteur force élémentaire de l’élément 1 :

1 Elé 3
0.5200
𝐹𝑒1 = 7.8 ∗ 10 ∗ 0.04 ∗ 0.333 [2] = [ ] 𝐾𝑁
1 0.5200
2

4
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI

Elément 02 et 3 :
On retrouve les mêmes matrices élémentaires Ke et vecteur Fe car on a les mêmes longueurs
élémentaires Le = 0.333 m.
Apres Assemblage, Le système global s’écrit :
2.5200 −2.5200 0 0 𝑈1 0.5200
5.0400 −2.5200 0 𝑈2 1.0400
1.0e + 07 ∗ [ ][ ] = [ ]
5.0400 −2.5200 𝑈3 1.0400
𝑆𝑦𝑚 2.5200 𝑈4 0.5200
0.5200Apres application des CAL on trouve :
5.0400 −2.5200 0 𝑈1 1.0400
1.0e + 07 ∗ [ 5.0400 −2.5200] [𝑈2] = [1.0400]
𝑆𝑦𝑚 2.5200 𝑈3 0.5200
Et après résolution du système :
𝑈2 = 1.032e-07 m, 𝑈3 = 1.651e-07 m, 𝑈3 = 1.857e-07 m avec 𝑈1 = 0 m

2. Calcul des déplacements en utilisant la formule analytique :

La solution analytique est donnée par :

avec W represente le poids total de la barre. ,𝑊 = 𝜌𝑔𝐴𝐿 = 7.8 ∗ 10 ∗ 0.04 ∗ 1 = 1.0400 𝐾𝑁
Les résultats de u(x) pour les différents maillages sont représentés dans le tableau suivant :
Maillage 01 Maillage 02 Maillage 03
x solution solution x solution solution x solution Solution
(m) MEF analytique (m) MEF analytique (m) MEF analytique
0 0 0 0 0 0 0 0 0
L 0.1857 e-06 0.1857 e-06 L/2 0.1393 e-06 0.1393 e-06 L/3 0.1032 e-6 0.1032 e-6
L 0.1857 e-06 0.1857 e-06 2L/3 0.1651 e-6 0.1651 e-6
L 0.1857 e-6 0.1857 e-6

3. Représentation graphique des résultats :

On remarque que les résultats MEF et Analytique sont identiques.

5
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI
Solution Exercice 02
1. Calcul de déplacement du nœud 3
Pour le calcul des déplacements on doit résoudre le système suivant : K U = F avec :
 K : matrice de rigidité globale du système.
 F : vecteur force global.
 U : vecteur déplacement.
Calcul de la matrice de rigidité globale K
Pour le cas d’une barre inclinée avec un angle θ par rapport à l’horizontale, l’expression de la
matrice de rigidité élémentaire de la barre s’écrit :
𝐶2 𝐶𝑆 −𝐶 2 −𝐶𝑆
𝐴𝐸 𝑆2 −𝐶𝑆 −𝑆 2
𝐾𝑏 = [ ] avec c  cos(θ) et s  sin(θ).
𝐿 𝐶2 𝐶𝑆
2
𝑆𝑦𝑚 𝑆

Remarque :
D’une manière générale pour le calcul de cos(θ) et sin(θ), il est préférable d’utiliser les expressions
suivantes :

a. Matrice de rigidité élémentaire de la barre N°1

E = 210000 MPa A = 2500 mm² L = 3m
3−0 0−0
𝑐= = 1 𝑒𝑡 𝑠 = =0
3 3
On remplace les valeurs de E, A, L, c et s dans la matrice élémentaire de l’élément barre :
Ke1 = DDL
175000 0 -175000 0 U1
0 0 0 0 V1
-175000 0 175000 0 KN/m U2
0 0 0 0 V2

b. Matrice de rigidité élémentaire de la barre N°2

E = 210000 MPa A = 2500 mm² L = 1²+4² = 4.1231 m
2−3 4−0
𝑐= = −0.2425 𝑒𝑡 𝑠 = = 0.9701
4.1231 4.1231

6
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI

Ke2 = 1.0e+05 * DDL

0.0749 -0.2996 -0.0749 0.2996 U2
-0.2996 1.1984 0.2996 -1.1984 V2
-0.0749 0.2996 0.0749 -0.2996 KN/m U3
0.2996 -1.1984 -0.2996 1.1984 V3

c. Matrice de rigidité élémentaire de la barre N°3

E = 210000 MPa A = 2500 mm² L = 2²+4² = 4.4721 m
2−0 4−0
𝑐= = 0.4472 𝑒𝑡 𝑠 = = 0.8944
4.4721 4.4721

Ke3 = 1.0e+04 * DDL

2.3479 4.6957 -2.3479 -4.6957 U1
4.6957 9.3915 -4.6957 -9.3915 V1
-2.3479 -4.6957 2.3479 4.6957 KN/m U3
-4.6957 -9.3915 4.6957 9.3915 V3

2. Matrice de rigidité globale des trois éléments 1, 2 et 3

Matrice de rigidité globale vide :
K = DDL
0 0 0 0 0 0 U1
0 0 0 0 0 0 V1
0 0 0 0 0 0 U2
0 0 0 0 0 0 V2
0 0 0 0 0 0 U3
0 0 0 0 0 0 V3
Emplacement de Ke1 dans la matrice globale
K = DDL
175000 0 -175000 0 0 0 U1
0 0 0 0 0 0 V1
-175000 0 175000 0 0 0 U2
0 0 0 0 0 0 V2
0 0 0 0 0 0 U3
0 0 0 0 0 0 V3

Emplacement de Ke2 dans la matrice globale

K = 1.0e+05 *
DDL
0 0 0 0 0 0 U1
0 0 0 0 0 0 V1
0 0 0.0749 -0.2996 -0.0749 0.2996 U2
0 0 -0.2996 1.1984 0.2996 -1.1984 V2
0 0 -0.0749 0.2996 0.0749 -0.2996 U3
0 0 0.2996 -1.1984 -0.2996 1.1984 V3

7
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI

Emplacement de Ke3 dans la matrice globale

K = 1.0e+04 *
DDL
2.3479 4.6957 0 0 -2.3479 -4.6957 U1
4.6957 9.3915 0 0 -4.6957 -9.3915 V1
0 0 0 0 0 0 U2
0 0 0 0 0 0 V2
-2.3479 -4.6957 0 0 2.3479 4.6957 U3
-4.6957 -9.3915 0 0 4.6957 9.3915 V3

Assemblage des trois matrices Ke1, Ke2 et Ke3 dans la matrice globale

K = 1.0e+05 *
DDL
1.9848 0.4696 -1.7500 0 -0.2348 -0.4696 U1
0.4696 0.9391 0 0 -0.4696 -0.9391 V1
-1.7500 0 1.8249 -0.2996 -0.0749 0.2996 U2
0 0 -0.2996 1.1984 0.2996 -1.1984 V2
-0.2348 -0.4696 -0.0749 0.2996 0.3097 0.1700 U3
-0.4696 -0.9391 0.2996 -1.1984 0.1700 2.1376 V3
3. Calcul du vecteur force global

Nœud 1 : Fx = 0 Fy = 0
Nœud 2 : Fx = 0 Fy = 0
Nœud 3 : Fx = +65 KN Fy = 0
Donc le vecteur force s’écrit :
F = [0 0 0 0 +65 0]T KN
L’indice T signifie la transposée d’une matrice ou d’un vecteur.
4. Appliquer les conditions aux limites et système réduit
Avant l’application des conditions aux limites, le système matriciel KU=F s’écrit :
1.0e+05 *
1.9848 0.4696 -1.7500 0 -0.2348 -0.4696 U1 0
0.4696 0.9391 0 0 -0.4696 -0.9391 V1 0
-1.7500 0 1.8249 -0.2996 -0.0749 0.2996 U2 = 0
0 0 -0.2996 1.1984 0.2996 -1.1984 V2 0
-0.2348 -0.4696 -0.0749 0.2996 0.3097 0.1700 U3 65
-0.4696 -0.9391 0.2996 -1.1984 0.1700 2.1376 V3 0

Et pour appliquer les conditions aux limites, on a :

Nœud 1 : Appuis simple U1 ≠ 0 V1 = 0
Nœud 2 : Appuis double U2 = 0 V2 = 0
Nœud 3 : Libre U3 ≠ 0 V3 ≠ 0
Donc on supprime du système global les lignes et les colonnes d’indice 2, 3 et 4 ce qui nous ramène
au système réduit suivant :

1.0e+05 *
1.9848 -0.2348 -0.4696 U1 0
-0.2348 0.3097 0.1700 U3 = 65
-0.4696 0.1700 2.1376 V3 0

8
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI
5. Calcul du vecteur déplacement
Le calcul du vecteur déplacement revient à résoudre le système après l’application des conditions
aux limites (3 équations, 3 inconnues).
Apres résolution du système, on trouve :
U1 = 0.0002 m
U3 = 0.0024 m
V3 = -0.0001 m

Donc le vecteur déplacement s’écrit :

U = [0.0002 0 0 0 0.0024 -0.0001]T m

6. Calcul de l’effort interne dans la barre N° 03 et les réactions d’appuis au nœud 1.

Pour le calcul des efforts internes dans la barre N°3, il suffit juste de revenir au repère local de la
barre, et d’après le chapitre ‘’élément barre’’, on a :

Fe  Rf e  f e  RT Fe

Fe : vecteur force par rapport au repère global.

fe : vecteur force par rapport au repère local.
R : Matrice de transformation
cos  sin  0 0 
 sin  cos 0 0 
R
 0 0 cos  sin  
 
 0 0 sin  cos 

Donc, pour la barre N°3, fe3  R3T Fe3

Avec :

1.0e+04 *
2.3479 4.6957 -2.3479 -4.6957 0.0002 -43.3333
Fe3=Ke3*Ue3 =
4.6957 9.3915 -4.6957 -9.3915 0 -86.6667
-2.3479 -4.6957 2.3479 4.6957 0.0024 43.3333
-4.6957 -9.3915 4.6957 9.3915 -0.0001 86.6667

9
Université A. Mira, Béjaïa, Année Universitaire 2020-2021
Faculté de Technologie
Département de Génie–Civil
Méthode des Eléments Finis
Enseignant : Mr. L. RAMDANI

Donc, le vecteur force par rapport au repère global ;

Fe3 = [-43.3333 -86.6667 43.3333 86.6667]T

Et pour avoir le vecteur force par rapport au repère local fe3 qui représente aussi les efforts
internes dans la barre :
0.4472 0.8944 0 0 -43.3333 -96.8963
-0.8944 0.4472 0 0 -86.6667 0
fe3=R3T*Fe3 =
0 0 0.4472 0.8944 43.3333 96.8963
0 0 -0.8944 0.4472 86.6667 0

Donc, l’effort interne dans la barre N°3 est un effort de traction Nint3 = 96.8963 KN

Et pour le calcul des réactions d’appuis au nœud 1, il suffit tout simplement de revenir au
vecteur force par rapport au repère global Fe3, avec :

La valeur de la réaction au nœud 1 suivant l’axe x : R1_x = -43.3333 KN

La valeur de la réaction au nœud 1 suivant l’axe y : R1_y = -86.6667KN

10