Université A. Belkaïd ‐ Tlemcen    Méthode des éléments finis (GV841)
Faculté de Technologie    Filière : Travaux publics
Département de Génie Civil    Master1 voies et ouvrages d’art
 
Chapitre 4 : Eléments finis barres 
 
Exercice 1 : 
Soit la poutre ci‐dessous de longueur « 3L », avec « L=1,00  m » et de section « A=900 cm² ».  Elle est 
soumise  à  une  charge  concentrée  « F=2MN »,  appliquée  au  nœud  « 2 ».  Cette  poutre  est  de  rigidité 
variable :  les  tronçons    entre  les  nœuds  « 1 »  et  « 3 »  ont  une  rigidité  E1=10 000  MPa.  Alors  que  le 
tronçon « 3‐4 » a une rigidité E2=30 000 MPa. 
1/ Ecrire la matrice de rigidité de cette poutre.  
2/ Déterminer les déplacements aux nœuds. 

 
Discrétisation de cette poutre en des éléments et des nœuds.
Nœuds    Eléments 
N°  X  Y    N°  NO  NE  Longueur  Module de Young 
1  0  0    1  1  2  L  E1 
2  L  0    2  2  3  L  E1 
3  2L  0    3  3  4  L  E2 
4  3L  0             

1/ La matrice de rigidité de la poutre :

a/ Matrices de rigidité locales :
La poutre est soumise à un effort de traction‐compression dont la forme de la matrice de rigidité est :  

Elément 1 [MN/ml] Elément 2 [MN/ml] Elément 3 [MN/ml]

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 1/10

 b/ Assemblage:
   
   

  u1  u2  u3  u4      u1  u2  u3  u4  

                       

        0  0  u1   9  ‐9 0  0  u1
          0  u2     ‐9  18  ‐9  0  u2
[K]=  [K]= 102 x 
  0        u3   0  ‐9  36  ‐27 u3
 
  0  0      u4   0  0  ‐27  27 u4
                       

2/ Déplacements aux nœuds : 
  u1  u2  u3  u4             

                           

  9  ‐9 0  0  u1    F1        F1    0   

    ‐9  18  ‐9  0  u2    Avec    F2  =  F=2MN   
= F2 
102 x  
           
  0  ‐9  36  ‐27  u3  F3  F3  0 

  0  0  ‐27  27  u4    F4        F4    0   

                           

 
 
Conditions aux limites 
Nous  avons  un  appui  double  au  nœud  1  et  un  appui  simple  au  nœud  4  alors    les  2  DDL  sont  nuls :  
(u1=0 et u4=0). On élimine alors les 1ère et 4ème lignes et colonnes, on obtient alors le système suivant : 
9 9 0 0 0
9 18 9 0 2   18 9 2 10    ,
10  
0 9 36 27 0 9 36 0 ,
0 0 27 27 0
 
 
 
 
 
 
 
 
 
 
 

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 2/10

Exercice 2 : 
 
Soit le système treillis constitué de 3 barres de mêmes rigidités 
« E=200 KN/mm² » et de sections égales « A=3000 mm² », (figure 
ci‐contre). 
1/  Déterminer  la  matrice  de  rigidité  de  chaque  barre  de  ce 
système dans le repère local ; 
2/Exprimer  la  rigidité  de  chaque  barre  dans  un  repère  global 
(X,Y) placé au nœud « 1 » 
3/ Ecrire la matrice de rigidité du système (Assemblage). 
4/ Déterminer les déplacements aux nœuds. 
5/ Calculer les efforts dans les barres. 
1/ Matrices de rigidité locales : 
Définitions des nœuds et les éléments : 
Nœuds    Eléments 
N°  X  Y    Longueur 
N°  NO  NE 
(m) 
1  0  0    1  1  2  3√2 4,24 
2  3    3    2  2  3  3√2 4,24 
3  6  0    3  1  3  6 
 
La poutre est soumise à un effort de traction‐compression dont la forme de la matrice de rigidité est :  

Elément 1 [KN/ml] Elément 2 [KN/ml] Elément 3 [KN/ml]

√ √

 
2/La matrice de rigidité globale élémentaire : 
Les  3  éléments  du  système  sont  orientés  d’une  manière  différente,  pour  pouvoir  calculer  les 
déplacements au niveau des nœuds, on doit considérer un repère commun (X,Y). 
Avant de faire l’assemblage, on doit écrire la matrice de chaque élément dans le repère global (X,Y), ce 
qu’on appelle « matrice de rigidité globale élémentaire ». 
La matrice de rigidité de chaque élément transformée dans le repère (X,Y) est donnée par la relation :  
 
[Ke]=[]T[k][] 
 
 
 

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 3/10

Avec [] étant la matrice de transformation : 
0 0
Dans le plan (X,Y), elle est égale à :   ce qui nous donne : 
0 0

 
Et   et   et i et j étant les nœuds n’origine et d’extrémité des éléments. 

Elément 1 [KN/ml] Elément 2 [KN/ml] Elément 3 [KN/ml]

6 3
0,707   0,707 1 et 0
, 4,24

 
√ √
   

On peut écrire les matrices sous la forme :  

√ 1 1
K  ; avec  10  
1 1

De même : 

√ 1 1
K  ; avec  10  
1 1

1 0
Et  K  ; avec  10  
0 0

 
3/ Matrice de rigidité su système (Assemblage) : 
La matrice de rigidité su système est obtenue par la somme des différentes matrices de rigidité 
élémentaires. Sa taille est de 2x nombre de nœuds (6x6) : 

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 4/10

 √ √ √ √

√ √ √ √
, , , ,
√ √ √ √ , , , ,
√ , , , , ,
√ √ √ √ , , , , ,
  
√
, , , ,
√ √ √ √ , , , ,
√ √ √ √

4/Les déplacements aux nœuds : 
Les conditions aux limites : 
Nous avons un appui double au nœud 1 (Ux=0 et Uy=0) et un appui simple au nœud 3 (Uy=0), alors  les 
3  DDL  correspondent  aux  1ère,  2ème  et  6ème  lignes  et  colonnes.  Ces  3  lignes  et  colonnes  seront 
supprimées, telles que : 
, , , ,
, , , ,
, , , , ,
 
, , , , ,
, , , ,
, , , ,

0
0
Le vecteur forces est donné par :  100  
200
0
0
 
On obtient alors le système à résoudre suivant : 
, , 100 ,
, , 200    ,  
, , , 0 ,

5/ Les efforts dans les barres : 
Pour calculer les efforts dans les barres, il faut revenir sur les axes locaux (x,y) : 
L’effort local dans une barre quelconque est donné par :   N ij 
EA
u j  ui 
L

En fonction des coordonnées globales :   

Ce qui donne :   

Elément 1 [KN] Elément 2 [KN] Elément 3 [KN]

     

,   ,

 
   

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 5/10

Exercice 3 : 
 
Soit  le  système  montré  en  figure  ci‐contre.  Celui‐ci  est  composé  d’une 
poutre de caractéristiques « E et A » et d’un ressort de rigidité « k ». Une 
charge d’intensité « P=2 KN » est appliquée au nœud « 2 ». 
1/Quel est le nombre de degrés de liberté pour ce système ? 
2/  Ecrire  la  matrice  de  rigidité  locale  de  chaque  élément  de  cette 
structure. 
3/ Ecrire la matrice globale élémentaire pour chaque élément. 
4/ Donner la matrice de rigidité globale de la structure 
5/ Déterminer le déplacement au nœud 2. 
6/ Calculer les efforts normaux dans les éléments. 
Données : E = 10 000 MPa, A= 25 mm², L=1,00 m et k = 100 KN/ml. 
 
1/ Le nombre de degrés de liberté pour ce système : 
Localement, nous avons un 1DDL par nœud, ce qui nous fait 2DDL par élément. 
Les éléments ne sont pas orientés dans le même repère, donc on doit se référer 
par rapport un repère global (X,Y) qu’on placera n’importe où, par exemple : 
 
Dans le repère global, nous avons 2DDL par nœud ce qui nous fait 6DDL pour tout 
le système. 
 
2/ La matrice de rigidité locale de chaque élément de cette structure : 
Elément 1-2 [KN/ml] Elément 2-3 [KN/ml]

3/ La matrice globale élémentaire pour chaque élément : 
Elément 2 [KN/ml]
Elément 1 [KN/ml]
1 0
0 ; 1  0,707
√2

  ,  

 
4/ La matrice de rigidité globale de la structure 
0
 
0

0 0
Avec :  K  ; avec  100  
0 1

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 6/10

 1 1
Et  K  ; avec  88,39  
1 1

0 0 0 0
0
0 100 0 100
0 0 88,39 88,39 88,39 88,39
 
0 100 88,39 188,39 88,39 88,39
88,39 88,39 88,39 88,39
0
88,39 88,39 88,39 88,39
5/ Le déplacement au nœud 2 : 
Les 2DDL des nœuds 1 et 3 sont bloqués, donc on supprime les lignes et colonnes suivantes : (1ère, 2ème, 
5ème et 6ème)  ce qui reste comme système : 
88,39 88,39 0 2
88,39 188,39 20    2  

6/ Les efforts normaux dans les éléments : 
Pour calculer les efforts dans les barres, il faut revenir sur les axes locaux (x,y) : 

Elément 1 [KN] Elément 2 [KN]

   

 
 
   

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 7/10

 Exercice 4 : 
 
Soit  le  système  treillis  3D  montré  en  figure.  Il  est  soumis  à 
une charge concentrée au nœud 5 d’intensité 4KN. Les barres 
de la structure ont la même section « A = 25 mm² ». La barre 
« 1‐5 » a une rigidité égale à « 2E » alors que les autres barres 
possèdent la même rigidité « E=10 000 MPa ». Ce système est 
encastré à sa base. 
1/  Quel  est  le  nombre  de  degrés  de  liberté  de  ce  système ? 
Qu’elle serait alors la taille des matrices de rigidités (locales et globales) ? 
2/ Ecrire les matrices de rigidités locales. 
3/ Exprimer la rigidité élémentaire dans le repère global (X,Y,Z). 
4/ Ecrire la matrice de rigidité de la structure. 
5/ Déterminer le déplacement au nœud 5. 
 
1/  Le  nombre  de  degrés  de  liberté  de  ce  système  et  la  taille  des  matrices  de  rigidités  (locales  et 
globales) : 
 
Le système étant globalement dans le repère (X, Y, Z), alors chaque nœud possède 3DDL (Ux, Uy et Uz). 
Nous  avons  5  nœuds  ce  qui  donne  15  DDL  au  total.  Maintenant,  localement,  il  s’agit  d’un  système 
treillis,  alors  les  barres  sont  soumises  à  l’effort  de  traction‐compression,  ce  qui  veut  dire  1DDL  pour 
chaque nœud. 
La taille des matrices de rigidités locales est toujours (2x2) alors que la taille de la matrice globale est 
de (15x15). 
 
2/ Les matrices de rigidités locales : 
Avant de définir les matrices de rigidités, nous allons définir les connectivités des éléments : 
 
Nœuds    Eléments 
N°  X  Y  Z    N°  NO  NE  Rigidité  Longueur (m) 
1  0  0  0    I 1  5  2EA  5 
2  4  0  1    II 2  5  EA  √10 3,16 
3  4  0  ‐1    III 3  5  EA  √10 3,16 
4  6  0  0    IV 4  5  EA  √13 3,61 
5  4  3  0             
 
Elément 1-5 [KN/ml] Elément 2-5 [KN/ml] Elément 3-5 [KN/ml] Elément 4-5 [KN/ml]
2EA 1 1 EA 1 1 EA 1 1 EA 1 1
K K K K
L 1 1 L 1 1 L 1 1 L 1 1
1 1 1 1 1 1 1 1
K 100 K 79,11 K 79,11 K 69,25  
1 1 1 1 1 1 1 1
 

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 8/10

 3/ La rigidité élémentaire dans le repère global (X,Y,Z) : 
La matrice de rigidité de chaque élément transformée dans le repère (X,Y,Z) est toujours donnée par la 
relation :  
[Ke]=[]T[k][] 
Avec [] étant la matrice de transformation : 
0 0 0
Dans le repère (X,Y,Z), elle est égale à :   ce qui nous donne : 
0 0 0

Elément I [KN/ml] Elément II [KN/ml]

0,8 ;  0,6 ;  0  0 ;  0,95 ; 0,316
, ,

0,64 0,48 0 0,64 0,48 0 0 0 0 0 0 0

0,48 0,36 0 0,48 0,36 0 0 0,9 0,3 0 0,9 0,3
0 0 0 0 0 0 0 0,3 0,1 0 0,3 0,1
K 100   K 79,11  
0,64 0,48 0 0,64 0,48 0 0 0 0 0 0 0
0,48 0,36 0 0,48 0,36 0 0 0,9 0,3 0 0,9 0,3
0 0 0 0 0 0 0 0,3 0,1 0 0,3 0,1

Elément III [KN/ml] Elément IV [KN/ml]

0 ;  0,95 ; 0,316 0,55;  0,83 ;  0
, , , ,

0 0 0 0 0 0 0,3 0,46 0 0,3 0,46 0

0 0,9 0,3 0 0,9 0,3 0,46 0,69 0 0,46 0,69 0
0 0,3 0,1 0 0,3 0,1 0 0 0 0 0 0
K 79,11   K 69,25  
0 0 0 0 0 0 0,3 0,46 0 0,3 0,46 0
0 0,9 0,3 0 0,9 0,3 0,46 0,69 0 0,46 0,69 0
0 0,3 0,1 0 0,3 0,1 0 0 0 0 0 0

 
4/ La matrice de rigidité de la structure :  
0 0 0
0 0 0
0 0 0  
0 0 0

0,64 0,48 0
Avec :  K  ; avec  100 0,48 0,36 0  
0 0 0
0 0 0
K  ; avec  79,11 0 0,9 0,3  
0 0,3 0.1
0 0 0
K  ; avec  79,11 0 0,9 0,3  
0 0,3 0.1
0,3 0,46 0
K  ;avec  69,25 0,46 0,69 0   
0 0 0

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 9/10

 0,64 0,48 0 0 0 0 0 0 0 0 0 0 0,64 0,48 0
100 0,48 0,36 0 0 0 0 0 0 0 0 0 0 100 0,48 0,36 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 79,11 0 0,9 0,3 0 0 0 0 0 0 79,11 0 0,9 0,3
0 0 0 0 0,3 0.1 0 0 0 0 0 0 0 0,3 0.1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 79,11 0 0,9 0,3 0 0 0 79,11 0 0,9 0,3  
0 0 0 0 0 0 0 0,3 0.1 0 0 0 0 0,3 0.1
0 0 0 0 0 0 0 0 0 0,3 0,46 0 0,3 0,46 0
0 0 0 0 0 0 0 0 0 69,25 0,46 0,69 0 69,25 0,46 0,69 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0,64 0,48 0 0 0 0 0 0 0 0,3 0,46 0 0,64 0,48 0 0 0 0 0 0 0 0,3 0,46 0
100 0,48 0,36 0 79,11 0 0,9 0,3 79,11 0 0,9 0,3 69,25 0,46 0,69 0 100 0,48 0,36 0 79,11 0 0,9 0,3 79,11 0 0,9 0,3 69,25 0,46 0,69 0
0 0 0 0 0,3 0.1 0 0,3 0.1 0 0 0 0 0 0 0 0,3 0.1 0 0,3 0.1 0 0 0

, , , ,
, , , ,

, , , ,  
, , , ,
, , , ,
, , , ,

, , , ,
, , , , , , , ,
, , , , ,
 
5/ Déterminer le déplacement au nœud 5 : 
Le système est encastré à la base, c’est‐à‐dire les nœuds 1n 2, 3 et 4 sont bloqués ; ce qui permet d’éliminer les lignes et colonnes correspondantes : 
Nœud 1 : Lignes et colonnes 1, 2 et 3   ; Nœud 2 : Lignes et colonnes 4, 5 et 6   ; ; Nœud 3 : Lignes et colonnes 7, 8 et 9   ; Nœud 4 : Lignes et colonnes 10, 11 
et 12 ; Le système à résoudre est donc : 
, , ,
, ,    ,  
,

Dr. Zahira BENADLA (GC, Faculté de Technologie, Tlemcen) 10/10