1

Sommaire

I/ Introduction : ....................................................................................................................................... 3
1. Structure : .................................................................................................................................... 3
2. Discrétisation du domaine :......................................................................................................... 3
3. Hypothèse :.................................................................................................................................. 4
4. Travail demandé : ........................................................................................................................ 4
II/ Matrice de Rigidité Globale : .............................................................................................................. 5
1. Matrice de rigidité de la poutre horizontale : Elément fini 3 ..................................................... 5

2. Matrice de rigidité de la poutre verticale : Elément fini 1 .......................................................... 6

3. Matrice de rigidité de la poutre verticale : Elément fini 2 .......................................................... 7
4. Matrice de rigidité de la poutre verticale : Elément fini 4 .......................................................... 7
5. Matrice de rigidité de la poutre verticale : Elément fini 5 ......................................................... 8
6. Assemblage :................................................................................................................................ 9
III/ Vecteurs forces et déplacements : .................................................................................................. 12

2
 I/ Introduction :

1. Structure :
La structure étudiée est constituée de 5 poutres considérées comme des
éléments finis (Figure 1). Il s’agit du Mi i-projet numéro 11.

Figure 1: Figure Mini-projet numéro 11

2. Discrétisation du domaine :
Il s’agit des éléments finis ci-après:
Tableau 1 : Tableau des connections entre éléments

Elément Nœud i Nœud j Longueur

1 1 3 L1 = 2.5 m

2 2 4 L2 = 2 m

3 3 4 L3 = 2.5 m

4 3 5 L4 = 2 m

5 4 6 L5 = 1.5 m

3
  La structure est soumise à des sollicitations connues :
- Forces : F1  20 KN et F2  10 KN ;
- Moments : M 1  20 KN .m et M 2  15 KN .m

 Ainsi que des sollicitations inconnues au niveau des appuis :

- Encastrement : Nœud 1 ;
- Appuis 2 et 5 .

 Les conditions aux limites :

- Appuis si ples aux œuds 2 et ;
- E ast e e t au iveau du œud 1 .

C’est u e st u tu e situ e da s le pla . C’est u p o l e bidimensionnel.

3. Hypothèse :
On suppose que les 5 poutres ont les mêmes caractéristiques à savoir :
- E : Module d’You g ;
- S : Section ;
- I : Moment quadratique.

4. Travail demandé :

 Définir la matrice de rigidité globale de la structure ;

 Définir les vecteurs forces et déplacements ;
 Calculer le vecteur déplacement ;
 Cal ule les a tio s d’appuis.

4
II/ Matrice de Rigidité Globale :
Pour définir la matrice de rigidité globale de la structure étudiée, on doit
déterminer les matrices de rigidité élémentaires des 5 poutres « éléments finis ».
Pour un élément finis quelconque e , de longueur L e , situé dans le plan et faisant
un angle  ave l’axe ho izo tal, et travaillant en traction et en flexion simple, la
matrice de rigidité élémentaire est définie comme suit :

 
 A e c 2  12B e s 2 A e  12B e  cs  6B e L e s 
 A e c 2  12B e s 2   A e  12B e  c s  6 B e L es

 A e  12B e  cs A e s  12B e c
2 2
6B e L e c  A e  12B e  c s 
 A e s  12B e c
2 2
 6 BeLec 

 2 B e L2e 
K e    62 B e L e s 2  6 BeLec
2
6B e L e c 4 BeL 6 Be Les

e


  A e c  12B e s   A e  12B e  c s 6 Be Les A e c  12B e s
2 2
A e  12B e  cs 6 BeLes 
  A  12B  cs
 e e 
 A e s 2  12B e c 2   6 BeLec A e  12B e  c s A e s 2  12B e c 2  6 B e L e c

  6 BeLes 6 BeLec 2 B e L2e 6 Be Les  6 BeLec 4 B e L2e 

Avec :

 ES
A e  L
 e

 EI
B e  L3
 e

 c  cos 
s  sin 

cs  cos  sin 




Appliquons ce résultat pour les 5 poutres de la structure étudiée.

1. Matrice de rigidité de la poutre horizontale : Elément fini 3

  0
c  1

Pour cette poutre : s  0
cs  0

L e  L 3  2.5m

5
La matrice de cette poutre est donnée par :
 ES ES 
 L 0 0  0 0 
L3
 3

 0 12EI 6EI
0 
12EI 6EI 
 L33 L23 L33 L23 
 6EI 4EI 6EI 2EI 
 0 0  2 
 L23 
K 
3

ES
L3
ES
L3 L3

  0 0 0 0 
 L3 L3 
 12EI 6EI 12EI 6EI 
 0   0  
 L33 L23 L33 L23 
 0 6EI 2EI
0
6EI
 2
4EI 
 L23 L3 L3 L3 

2. Matrice de rigidité de la poutre verticale : Elément fini 1

 
  2

c  9
Pour cette poutre s  1
cs  0

L e  L 1  2.5m


La matrice de cette poutre est donnée par :

 12EI 6EI 12EI 6EI 
 L3 0   0 
L21 L31 L21 
 1

 0 ES ES
0 0  0 
 L1 L1 
 6EI 4EI 6EI 2EI 
 2 0 0 
K    12LEI
1 1 L1
6EI
L21
12EI
L1 
6EI 
 3 0 0 
 L1 L21 L31 L21 
 0 
ES
0 0
ES
0 
 L1 L1
 
  6EI 0
2EI 6EI
0
4EI 
 L21 L1 L21 L 1 

6
 3. Matrice de rigidité de la poutre verticale : Elément fini 2

 
  2

c  9
Pour cette poutre s  1
cs  0

L e  L 2  2m


La matrice de cette poutre est donnée par :

 12EI 6EI 12EI 6EI 

 L3 0   0 
L22 L32 L22 
 2

 0 ES ES
0 0  0 
 L2 L1 
 6EI 4EI 6EI 2EI 
 2 0 0 
K    12LEI
2 2 L2
6EI
L22
12EI
L2 
6EI 
 3 0 0 
 L2 L22 L32 L22 
 0 
ES
0 0
ES
0 
 L2 L1
 
  6EI 0
2EI 6EI
0
4EI 
 L22 L2 L22 L 2 

4. Matrice de rigidité de la poutre verticale : Elément fini 4

 
  2

c  9
Pour cette poutre s  1
cs  0

L e  L 4  2m


7
 La matrice de cette poutre est donnée par :

 12EI 6EI 12EI 6EI 

 L3 0   0 
L24 L34 L24 
 4

 0 ES ES
0 0  0 
 L4 L4 
 6EI 4EI 6EI 2EI 
 2 0 0 
K    12LEI
4 4 L4
6EI
L24
12EI
L4 
6EI 
 3 0 0 
 L4 L24 L34 L24 
 0 
ES
0 0
ES
0 
 L4 L1
 
  6EI 0
2EI 6EI
0
4EI 
 L 4 L 4 
2
L4 L24

5. Matrice de rigidité de la poutre verticale : Elément fini 5

 
  2

c  9
Pour cette poutre s  1
cs  0

L e  L 5  1.5m


La matrice de cette poutre est donnée par :

 12EI 6EI 12EI 6EI 
 L3 0   0 
L25 L35 L25 
 5

 0 ES ES
0 0  0 
 L5 L5 
 6EI 4EI 6EI 2EI 
 2 0 0 
K    12LEI
5 5 L5
6EI
L25
12EI
L5 
6EI 
 3 0 0 
 L5 L25 L35 L25 
 ES ES 
 0  0 0 0 
L5 L5
 
  6EI 0
2EI 6EI
0
4EI 
 L25 L5 L25 L 5 

8
 6. Assemblage :

La Matrice de rigidité globale K  est une matrice 18 18 . En effet la structure
contient 6 œuds et ha ue œud à 3 degrés de liberté ( 2 déplacements et 1
rotation). Vu sa taille, on va la subdiviser en 4 matrices K 1  , K 2  , K 3  et
K 4 

K   K 3  K 4  

 1  2  

K  K  

Les 4 matrices sont comme suit :

 
 U1 V1 1 U2 V2 2 U3 V3 3 
  EI   EI   EI   EI  
  12 3  0   6 2  0 0 0   12 3  0   6 2  
  L1   L1   L1   L1  
  ES   ES  
 0   0 0 0 0 0    0 
  L1   L1  
  EI   EI   EI   EI  
  6 2  0   4  0 0 0  6  0   2  
  L1   L 
 1
 L2 
 1
 L 
 1 
  EI   EI  
 0 0 0  12 3  0   6 2  0 0 0 
  L2   L2  
 
K 1   0 0 0

0  ES 

0 0 0 0 
  L2  
  EI   EI  
 0 0 0   6 2  0  4  0 0 0 
  L2   L2  
        
  12 

EI
3 
0 EI
 6 2  0 0 0
EI
 12 3 
ES EI
 12 3  0  0  0 EI EI
 6 2  0  6 2  
 L1   L1   L1 L 3 L4   L1 L4  
 
 0  ES   ES EI ES   
   0 0 0 0 0  0  0   12 3  
EI
 0  6 2  0  
  L1   L1 L 3 L4   L3  
 
         EI  
  6 EI
  L21 
0   2
 L 
EI
 0 0 0  6
 L2
EI
 0  6
EI
  0  6
EI
 0   4
EI
 4
EI
 4 
 1  1 L24  
 L23 

 L
 1 L3 L 4  
  
 Pr emier quart de la matrice globale 

Figure 2 : Premier quart de la matrice globale

9
  
 U4 V4 4 U5 V5 5 U 6 V 6 6 
 
 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
 EI   EI  
  12 3  0   6 2  0 0 0 0 0 0
 L2   L2  
  ES  
 0    0 0 0 0 0 0 0
  L2  
K 2   
 EI   EI  
  6 2  0  2  0 0 0 0 0 0
  L2   L2  
  ES   EI   EI  
    0 0   12 3  0   6 2  0 0 0 
  L3   L4   L4  
  EI   EI   ES  
 0   12 3   6 2  0    0 0 0 0
  L3   L3   L4  
  EI   EI   EI   EI  
 0   6 2   2   6 2  0  2  0 0 0
    L3 
3 
L 
4
L 
4
L 
 
 Deuxième quart de la matrice globale 
Figure 3: Deuxième quart de la matrice globale

 
 U1 V1 1 U2 V2 2 U3 V3 3 
  EI      
0 0 0   12 3  0  6 EI2    ES  0 0 
  L2   L2   L3  
  ES   EI   EI  
0 0 0 0    0 0   12 3    6 2  
  L2   L3   L3 
  EI   EI   EI   EI  
0 
0 0  6 2   0  2    0 
 6 L2   2  
  L 2   L 2   3   L3  
   EI  
EI 
K 3   0 0 0 0 0 0   12 3  0  6 2  
  L 4   L4  
 
0 0 0 0 0 0 0 
ES
0 
 L4 
  EI   EI  
0 0 0 0 0 0   6 2  0  2  
  L4   L4  
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 
   
 Troisième quart de la matrice globale 
Figure 4 : Troisième quart de la matrice globale

10
11
  
 U4 V4 4 U5 V5 5 U6 V6 6 
 EI ES EI   EI EI   EI   
 12 3   12 3  0  0  0  6 2  0  6 2  0 0 0   12 3  0   6 EI2 
 L 2 L 3 L 5   L 2 L 5   L 5   L 5 
  ES EI ES     ES  
 0  0  0   12 3  
EI
 0  6 2  0  0 0 0 0    0 
  L2 L 3 L5   L3   L5  
  EI     EI   EI   EI  
 6 2  0  6 2  EI EI
 0  6 2  0   4 4
EI
 4 
EI
0 0 0  6 2  0  2  
  L 2 L 5   L3   L2 L3 L5   L5   L5  
  EI   EI  
 0 0 0  12 3  0  6 2  0 0 0 
  L4   L4  
 
K 4    0 0 0
 ES 
0   0 0 0 0 
  L4  
 
   
 0 0 0
EI
 6 2  0  4 
EI
0 0 0 
  L4   L4  
 
       
 EI
  12 3  0 EI
 6 2  0 0 0
EI
 12 3  0 EI
 6 2  
  L5   L5   L5   L5  
 
   


0  
ES
 0 0 0 0 0 
ES
 0 

 L5   L5 
 
  EI   EI   EI   EI  

  6 2  0  2  0 0 0  6 2  0  4 
L5  
 
L5   L5 
5
L   
   
 Q uatrièmequart de la matrice globale 
Figure 5: Quatrième quart de la matrice globale

III/ Vecteurs forces et déplacements :

Les vecteurs forces F et déplacements U sont définis comme suit :
 R 1 X  Inconnue   U 1  0
   
 R 1Y  Inconnue   V1  0 
M 1Z  Inconnue   1  0 
   
F2 X  0  U 2 
 R  Inconnue  V  0 
 2Y   2 
M 2 Z  0   2 
   
F3 X  F1  20KN  U 3 
F3 Y  0   V3 
   
M 3 Z  0   3
F    U   

F4 X  0  U
 4 
F4 Y  0   V4 
   
M 4 Z  M 1  20KN .m   4 
F  0  U 
 5X   5 
 R 5 Y  Inconnue   V5  0 
   
M 5 Z  M 2  15KN .m   5 
F6 X  F2  10KN  U 6 
   
F6 T  0   V6 
M  0   
 6Z   6 

12
Les lig es , , , et o tie e t des i o us ui doive t t e eti s si l’i souhaite soud e le
problème sur Matlab.

On remplace les valeurs numériques de la longueur de chaque élément sur Matlab pour obtenir la

matrice globale de rigidité K  . Par souci de clarté, on la subdivisera en deux K  , K 

1  2 
.

On rappelle que = �. � et = �. �

On ne peut résoudre le système sans supprimer les lignes et colonnes correspondants aux réactions
inconnues en outre la ligne et colonne 1,2,3,5 et 14.

On extrait alors la matrice B (13x13) après avoir opéré les commandes suivantes :

A=K

A=A(4:18,4:18)

A(2,:)=[]

A(:,2)=[]

B=A

B(10,:)=[]

B(:,10)=[]

Encore une fois , pour visualiser la matrice B , on la subdivise en deux :

13
14
On peut à présent résoudre le système réduit {�} ∗ {�} = {�} en utilisant la commande : u=B\F sur
Matlab (Elimination de Gauss) avec F comme suit :

O o tie t ai si u e at i e des d pla e e ts * u’o agrandit par la suite puisque

U1=V1=�1 =V2=V5 =0.

15
 U1  0 
V  0 
 1 
1  0 
 
U 2  (25 * (10752 * a^3 + 139200 * a^2 * b + 276350 * a * b^2 + 9375 * b^3))/(2* a * b * (5376 * a^2 + 58400 * a * b + 65625 * b^2)) 
V2  0 
 
 2  (125 * (5376 * a^2 + 21520 * a * b + 1875 * b^2))/(4* a * (5376 * a^2 + 58400 * a * b + 65625 * b^2)) 
U  -(250* (224 * a - 1875 * b))/(5376* a^2 + 58400 * a * b + 65625 * b^2) 
 3 
V  -(15000* (8 * a + 25 * b))/(5376* a^2 + 58400 * a * b + 65625 * b^2)3 
 
 3  (25 * (16128 * a^2 + 44800 * a * b + 46875 * b^2))/(4* a * (5376 * a^2 + 58400 * a * b + 65625 * b^2)) 
U   U 4  (75 * (1792 * a^2 + 18720 * a * b + 28125 * b^2))/(b* (5376 * a^2 + 58400 * a * b + 65625 * b^2)) 

V  (27000 * (8 * a + 25 * b))/(5376* a^2 + 58400 * a * b + 65625 * b^2) 
 4 
 4  (125 * (5376 * a^2 + 21520 * a * b + 1875 * b^2))/(4* a * (5376 * a^2 + 58400 * a * b + 65625 * b^2)) 
 
U 5  - (5 * (167552 * a^2 + 737300 * a * b + 1021875 * b^2))/(2* a * (5376 * a^2 + 58400 * a * b + 65625 * b^2)) 
V5  0 
 
 5  (5 * (209664 * a^2 + 1625600 * a * b + 1809375 * b^2))/(4* a * (5376 * a^2 + 58400 * a * b + 65625 * b^2)) 
U  (15 * (71680 * a^3 + 646656 * a^2 * b + 937400 * a * b^2 + 346875 * b^3))/(8* a * b * (5376 * a^2 + 58400 * a * b + 65625 * b^2))
 6 
V6  (27000 * (8 * a + 25 * b))/(5376* a^2 + 58400 * a * b + 65625 * b^2) 
  (5 * (43008 * a^2 + 6200 * a * b - 271875 * b^2))/(2* a * (5376 * a^2 + 58400 * a * b + 65625 * b^2)) 
 6 
 

16
Les a tio s d’appuis so t e suite al ul es e appli ua t l’ uatio F=K*U.

17
18