Chapitre 2
: MEF appliquée aux poutres
2.1 Méthode directe pour les structures à éléments discrets:
A. Introduction
La méthode directe est une approche pour des systèmes discrets, basée sur la méthode des
rigidités. Elle est de loin le procédé le plus simple pour introduire les concepts de base de
MEF et présente les avantages suivants :
Application de concept physique (équilibre des forces, conservation d’énergie,
conservation de masse,…) directement à des éléments discrets.
Facile dans son interprétation physique.
Ne demande pas de concept ou de manipulation mathématique sophistiquée.
Son application est limitée à un certain nombre de problèmes pour lesquelles les lois
d’équilibre et de conservation peuvent être facilement exprimées en termes des quantités
physiques que l’on désire obtenir (déplacements).
En général, les systèmes à treillis et les portiques sont constitués d’éléments discrets d’eux-
mêmes selon le sens physique et sont de parfaits exemples pour illustrer la méthode.
B. Elément fini barre (Treillis plan de barres) Le système à treillis est formé d’un
ensemble membrures appelé barres, qui sont sollicités par des efforts agissant le long de leur
axe moyen.
Le schéma statique de la barre impose la présence de rotules aux deux extrémités (nœuds).
Ce sont les premiers éléments présentés par la MEF suivi des éléments poutres assemblés en
ossature.
Figure 2.1 : (a) Structure à treillis représentant une ferme de toiture.
(b)Modèle d’élément barre
1
C. Elément ressort linéaire :
L’élément barre possède des caractéristiques similaires à celles d’un ressort élastique.
Figure 2.2 : Analogie barre- ressort
On considère que chaque élément de structure se comporte comme un ressort élastique c'est-
à-dire que la relation charge déplacement est linéaire.
i k j F
u
Figure2.3 : Déformation d’un ressort élastique
On appelle k la raideur (rigidité) qui correspond à la pente du graphe charge –déplacement.
Pente k
Force F
Pente k
Déplacement u de l’extrémité
Figure2.4 : Relation force /déplacement d’un ressort élastique
Connaissant la valeur de la rigidité et de la charge appliquée on a les relations :
1
F=k . u ↔u= . F
k
D. Formulations en éléments finis :
D.1. Matrice de rigidité élémentaire
2
Figure 2.5 : Ressort équivalent d’une barre à deux rotules
Convention de signes :
Cette même convention est adoptée pour les charges et les déplacements.
-(f,u) +(f,u)
Figure 2.6 : Convention de signes
Equilibre de la barre nous donne:
Noued1 : f 1=−k . ( u2−u 1 )
Noued2 : f 2=k . ( u2−u 1)
Ecriture matricielle pour un élément :
{}[
f1
f2
=
k −k u1
−k k
.
]{ }
u2 signifie que
¿.
Avec
: Vecteur résultant des forces nodales
: Matrice de rigidité élémentaire
: Vecteur résultant des déplacements nodaux
Remarque1:
Une colonne de Ke représente le vecteur charge qui doit être appliqué aux nœuds de l’élément
pour obtenir un état de déformation où le degré de liberté nodal est égal à 1 alors que les
autres sont nuls.
Exemple : u1=0 et u2=1
3
{}[
f1
f2
=
k −k 0
−k k 1 ]{} { }
. =k . −1
1
Donc le produit ¿ . représente la 2ème colonne de Ke
D.2 Matrice de rigidité d’un ensemble d’éléments (Assemblage) :
Il est question dans cette étape de décider comment les matrices de rigidités élémentaires sont
combinées ensembles pour former la matrice d’une structure composée de plusieurs éléments.
Pour la simplicité de l’analyse étudions dans un premier temps le système de deux barres
colinéaires de caractéristiques différentes assimilées aux deux ressorts colinéaires.
Figure 2.7 : Système à deux ressorts colinéaires aux extrémités libres
Pour étudier le système à deux ressorts colinéaires il suffit d’étudier chaque ressort (élément)
tout seul, puis de faire assemblage. Donc la modélisation est la suivante (Figure 2.7).
Figure 2.7 : Système à deux ressorts séparés
4
{ } [ ]{ }
f 1( 1 ) K 1 −K 1 0 u1
{ }[ ]{ }
(1 )
f1 K1 −K 1 u1
Elément 1: = . é quivalent f (1) = − K 1 K 1 0 . u2 (1)
f2
(1 ) −K 1 K1 u2 ⇔ 2
0 0 0 0 0
{ } [ ]{ }
0 0 0 0
{ }[
0
f 1( 2 )
Elément 2 : ( 2) −K
f2
=
K 2 −K 2 u2
2 K2
.
u3
é quivalent
⇔
f
f2
1
]{ }
(2 )
(2 )
= 0 K 2 −K
0 −K 2 K 2
2 . u2
u3
(2)
{ }[ ]{ } { }
( )
f11 K1 −K 1 0 u1 F1
(1)+(2) nous donne : f 2 + f 1 (1 ) (2 ) = −K 1 K 1 + K 2 −K 2 . u2 = F 2
f 2( 2 ) 0 −K 2 K2 u3 F3
Avec fij- Forces internes agissant au nœud i de l’élément j (i,j=1,2)
[ ] {} { }
K1 −K 1 0 u1 F1
¿ −K 1 K 1+ K 2 −K 2 ;= u2 et= F2
0 −K 2 K2 u3 F3
Cette opération est l’opération d’assemblage des matrices de rigidité élémentaires, la matrice
est appelée matrice de rigidité de la structure, le vecteur est le vecteur des inconnues de
déplacements et le vecteur est le vecteur des forces extérieures :
.=¿
Mise en œuvre pratique
La première étape consiste à écrire les deux matrices de rigidité des deux éléments en repérant
les lignes et les colonnes de chaque matrice par les inconnues de déplacements associées :
u1 u 2
[ K
[ K 1 ]= −K1
1
]
−K 1 u1
K 1 u2
u2 u 3
[ K
[ K 2 ]= −K2
2
]
−K 2 u2
K 2 u3
On range ensuite dans la matrice de rigidité de la structure chaque terme des deux matrices
à la ligne et la colonne correspondante :
5
u1 u 2 u3
[ ]
K1 −K 1 0 u1
¿ −K 1 K 1+ K 2 −K 2 u2
0 −K 2 K 2 u3
Le système à résoudre est alors .=¿
La deuxième étape consiste à faire le bilan des déplacements et des forces connus et inconnus.
En prenant un encastrement à l’extrémité gauche et en appliquant une force F à l’extrémité
droite (figure 2.8), on a :
Figure 2.8 : Système à deux ressorts colinéaires à une extrémité encastrée
{} {}
u1 ¿ 0 , connu F 1 inconnu , r é action à l ' encastrement
¿ u2 inconnu ¿ F 2 ¿0
u3 inconnu F3 ¿ F
Si le déplacement est connu en un nœud alors la force est inconnue, si la force est connue
alors le déplacement est inconnu.
La troisième étape est la résolution du système d’équations complet afin de déterminer toutes
les inconnues
6
[ ]{ } { }
K1 −K 1 0 0 F1
.
−K 1 K 1 + K 2 −K 2 u2 = 0
0 −K 2 K2 u3 F
−K 2
F 1= F
K1
1 1 ' K 1+ K 2
u3−u2= F et u 2= F d où u 3= F
K2 K1 K1 K2
E1 S 1 E 2 S2
Avec K 1= et K 2 =
L1 L2
Application :
Soit le système de deux barres colinéaires de caractéristiques différentes assimilées aux deux
ressorts colinéaires, figure 2.9, tel que :
Figure 2.9 : Système à deux ressorts colinéaires sollicités en traction
Le nœud 1 est fixe et les nœuds 2 et 3 supportent deux forces concentrées F2 = 75N et F3 =
50N.
Déterminez les déplacements aux nœuds libres et la réaction à l’appui sachant que les raideurs
des deux ressorts sont k1 = 50N/mm et k2 = 75N/mm.
Barre encastrée et soumise son poids propre
Soit une barre [AB] encastrée en A et soumise son poids propre (figure 2.10).
Soit une charge volumique ⃗f =ρg ⃗x c'est-à-dire ⃗
F =ρgv ⃗x = ρgS . L ⃗x
7
Barre soumise àson poids propre
L
On discrétise la barre en quatre éléments de longueur l= telle qu’elle indique la figure 2.12
4
8
{}
Le champ un solution du problème discrétisé est solution du système d’équations :
ρgLS
R+
8
{} [ ]{ }
F1 2 ρgLS
1 −1 0 0 0 0
F2 8 u2
−1 2 −1 0 0
2 ρgLS = 4 ES
F3 = 0 −1 2 −1 0 u3
8 L
F4 0 0 −1 2 −1 u4
2 ρgLS
F5 0 0 0 −1 1 u5
8
ρgLS
8
La résolution du système peut se faire de la manière suivante :
ρgLS −4 ES
R+ = u2
8 L
2 ρgLS 4 ES
8
=
L
( 2 u2−u 3 )
2 ρgLS 4 ES
8
=
L
( 2 u3−u 4−u2 )
2 ρgLS 4 ES
8
=
L
( 2 u4 −u3−u5 )
ρgLS 4 ES
8
=
L
( u5−u 4 )
Après un calcul tout fait, on trouve :
2
7 ρg L
u2 =
32 E
12 ρg L2
u3 =
32 E
15 ρg L2
u 4=
32 E
2
16 ρg L
u5 =
32 E
Et R=−ρgLS
Élément barre pour le calcul des treillis plans
Les barres composant un treillis plan sont positionnées arbitrairement dans l’espace et font
des angles différents avec le repère global de la structure ( ⃗ X ,⃗
Y ) (Fig. I.7).
On note θ l’angle entre l’axe ⃗
X du repère global et l’axe ⃗x du repère local à la barre. Le vecteur
déplacement d’un point de la barre s’écrit dans le repère local
9
10
