UNIVERSITE BATNA 2

Faculté de Technologie
Département de Génie Civil

Résumé du cours de
Méthode des Éléments Finis

Mr YAHIAOUI DJARIR
Année universitaire 2019-2020
 PRINCIPES DE LA METHODE
DES ELEMENTS FINIS

Introduction

La Méthode des éléments finis est une technique numérique pour trouver des solutions
approximatives d'équations différentielles partielles. Il résulte de la nécessité de résoudre des problèmes
complexes dans le domaine élasticité ou plastique et d'analyse structurelle en génie civil, mécanique et
aérospatial. Dans une simulation structurale, La Méthode des éléments finis aide à produire des
visualisations de la rigidité et la résistance de la structure. Cela contribue également à minimiser le poids
du matériau et son coût des structures. Et permet une visualisation détaillée et indique la répartition des
contraintes et des déformations dans les éléments structurale.

1.1 Du problème physique au modèle mathématique ou numérique

Nous partons d'un problème physique, le cadre précis de l'étude est défini par les hypothèses
simplificatrices qui permettent de définir un modèle mathématique. La difficulté pour l'ingénieur est de
savoir choisir parmi les lois de la physique celles dont les équations traduiront avec la précision voulue la
réalité du problème physique. Un bon choix doit donner une réponse acceptable pour des efforts de mise
en œuvre non prohibitifs. Les différentes étapes d'analyse d'un problème physique s'organisent suivant le
processus schématisé par la figure suivante :
 ÉLEMENT BARRE

Introduction

Nous nous intéressons ici aux structures composées de barres droites sollicitées uniquement en
traction ou compression. Ce sont des "treillis", structures formées de barres liées entre elles par des liaisons
de type "rotule", "pivot" ou "articulation" et chargée par des forces s’exerçant uniquement au niveau des
liaisons.

Figure 3-1 Treillis de barres articulées sous une charge ponctuelle

Pour dériver la matrice de rigidité d'un élément barres élastiques linéaires en utilisant les étapes
générales décrites au chapitre 2. Nous allons inclure une introduction à la fois sur un système local et global
des coordonnées. Et discuterons également de la transformation d'un vecteur du système de coordonnées
local vers un système de coordonnées globales.

3.1 Dérivation matrice de la rigidité pour un élément barre dans les coordonnées locales

En considérant un élément barre élastique d’une section transversale constante comme illustré dans
la figure 3-2. La barre est soumise à des forces de traction T dirigées vers le long de l'axe local de la barre
et appliquées aux nœuds 1 et 2.

A est la section constante de l’élément barre, E le module d'élasticité, et L la longueur initiale. Les degrés
de liberté nodaux sont des déplacements axiaux locaux représentés par u1 et u2 comme indiqués dans la
figure 3-2.
 Figure 3-2 Barre soumise à des forces de traction T; Les déplacements nodaux positifs et les forces sont
tous dans la direction x locale

D’après la loi de Hooke [Équation. (a)] et la relation déformation-déplacement [Équation. (b)], nous
écrivons

𝜎𝑥 = 𝐸𝜀𝑥 (a)

𝑑𝑢
𝜀𝑥 = (b)
𝑑𝑥

D’après l'équilibre de la force, nous avons

𝐴𝜎𝑥 = 𝑇 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 (c)

Pour une barre avec des charges appliquées uniquement aux extrémités. En utilisant l'équation (b)
dans l'équation (a) et ensuite l'équation (a) dans l'équation (c) et en différenciant par rapport à x, on obtient
l'équation différentielle régissant la barre avec un comportement élastique linéaire comme suit

𝑑 𝑑𝑢
( 𝐴𝐸 )=0 (d)
𝑑𝑥 𝑑𝑥

u est la fonction de déplacement axial sur le long de l'élément dans la direction x et A et E seront supposés
constants sur toute la longueur de la barre.

Les hypothèses suivantes sont utilisées pour dériver la matrice de la rigidité des éléments barres:

1. La barre ne peut supporter la force de cisaillement ou le moment de flexion, c'est-à-dire

f1y =0, f2y = 0, m1 = 0 et m2 = 0.

2. Tous les effets des déplacements transversaux sont ignorés.

3. La loi de Hooke s’applique ; C'est-à-dire que la contrainte axial x est liée à la déformation axiale x
par 𝜎𝑥 = 𝐸𝜀𝑥 .

4. Aucune charge intermédiaire appliquée.

Étape 1 Sélectionner du type d'élément

Représentez la barre avec numérotation les nœuds à chaque extrémité (Figure 3-2).
Étape 2 Sélectionner une fonction de déplacement

Parce qu’une fonction linéaire avec des points aux extrémités a un chemin unique. Par conséquent,

Supposons une variation de déplacement linéaire sur le long de l'axe x de la barre parce qu’une fonction
linéaire avec des points aux extrémités a un chemin unique. Ces points d'extrémité sont les valeurs nodales
u1 et u2.

𝑢 = 𝑎1 + 𝑎2 𝑥 (3.1.1)

Avec le nombre total des coefficients ai est toujours égal au nombre total de degrés de liberté associés à
l'élément. Ici, le nombre total de degrés de liberté est deux déplacements axiaux à chacun des deux nœuds
de l'élément. En utilisant la même procédure que dans la section 2.2 pour l'élément de ressort, nous
exprimons l'équation. (3.1.1) comme:

𝑢2 − 𝑢1
𝑢=( ) 𝑥 + 𝑢1 (3.1.2)
𝐿

Figure 3-3 Déplacement u tracé sur la longueur de l'élément

La raison pour laquelle nous convertissons la fonction de forme de déplacement de l'équation. (3.1.1) à
l'équation (3.1.2) est qu'elle nous permet d'exprimer la déformation en terme de déplacements nodaux en
utilisant la relation de déformation/ déplacement donnée par l'équation. (3.1.5) et ensuite relier les forces
nodales aux déplacements nodaux dans l'étape 4.

Sous forme matricielle, l'équation (3.1.2) devient

𝑢1
𝑢 = [𝑁1 𝑁2 ] {𝑢 } (3.1.3)
2

Avec des fonctions de forme données par

𝑥 𝑥
𝑁1 = 1 − 𝑁2 = (3.1.4)
𝐿 𝐿

La fonction de déplacement linéaire u (Équation. (3.1.2)) est tracée sur la longueur de l'élément barre, est
illustrée dans la figure 3-3.

Étape 3 Définir les relations déformation- déplacement et déformation- contrainte

La relation déformation - déplacement est

 𝑑𝑢 𝑢2 − 𝑢1
𝜀𝑥 = = (3.1.5)
𝑑𝑥 𝐿

Les équations (3.1.3) et (3.1.4) ont été utilisées pour obtenir l'équation (3.1.5), et la relation contrainte -
déformation est

𝜎𝑥 = 𝐸𝜀𝑥 (3.1.6)

Etape 4 Dérivation de la matrice de rigidité des éléments barre

La matrice de rigidité de l'élément barre est dérivée comme suit.

𝐴𝜎𝑥 = 𝑇 (3.1.7)

En utilisant les équations. (3.1.5) et (3.1.6) dans l'équation (3.1.7), nous obtenons

𝑢2 − 𝑢1
𝑇 = 𝐴𝐸 ( )
𝐿
(3.1.8)

En outre, par la convention de signe de force nodale de la figure 3-2,

𝑓1𝑥 = 𝑇 (3.1.9)

Lorsque nous substituons l'équation (3.1.8), l'équation (3.1.9) devient

𝐴𝐸
𝑓1𝑥 = (𝑢2 − 𝑢1 ) (3.1.10)
𝐿

De même,

𝑓2𝑥 = 𝑇 (3.1.11)

On remplace l'équation (3.1.8) dans l'équation (3.1.11)

𝐴𝐸
𝑓2𝑥 = (𝑢2 − 𝑢1 )
𝐿
(3.1.12)

Exprimant les équations. (3.1.10) et (3.1.12) sous forme matricielle, nous avons

𝑓1𝑥 𝐴𝐸 1 −1 𝑢1
{ }= [ ]{ } (3.1.13)
𝑓2𝑥 𝐿 −1 1 𝑢2

{𝑓} = 𝐾 {𝑑 }, nous avons, de l'équation (3.1.13),

[𝑘 ] =
𝐴𝐸 1 −1
[ ] (3.1.14)
𝐿 −1 1
L'équation (3.1.14) représente la matrice de rigidité d'un élément barre dans les coordonnées locales. Dans
l'équation (3.1.14), AE / L pour un élément barre est similaire à la constante de raideur k pour un élément
ressort.

Étape 5 Assemblages des matrice des éléments barres pour obtenir la matrice de rigidité globale

Assemblez les matrices de rigidité globale et les forces à l'aide de la méthode de rigidité directe. Cette étape
s'applique aux structures composées de plusieurs éléments.

[𝐾 ] = ∑𝑁
𝑒=1[𝑘
(𝑒)
] et {𝐹} = ∑𝑁
𝑒=1{𝑓
(𝑒)
} (3.1.15)

Toutes les matrices de rigidité locale des éléments barres [k e] doivent être transformées en matrice de
rigidité globale [k] (à moins que les axes locaux ne coïncident avec les axes globaux).

Étape 6 Calcul des déplacements nodaux

Déterminer les déplacements en imposant des conditions aux limites et en résolvant simultanément un
système d'équations, {𝐹 } = [𝐾] {𝑑 },

Étape 7 Calcul des forces interne

Enfin, déterminer les déformations et contraintes dans chaque élément barre par substitution des
déplacements dans les équations. (3.1.5) et (3.1.6).
 ÉLEMENT POUTRE

Introduction

Dans ce chapitre on développe la matrice de rigidité pour un élément de poutre travaillent en flexion,
la plupart des structures soit des bâtiments, des ponts et de beaucoup d'autres structures contiennent des
éléments en poutre. Une poutre droite (Figure 1.1) est un solide engendré par une surface plane A, constante
ou non, dont le centre de gravité G décrit le segment G1G2, le plan qui contient A restant normal à G1G2.
De plus, les dimensions de A sont faibles (sans être négligeables) par rapport à G1G2.

G1G2 est la fibre moyenne de la poutre.

A est la section droite de la poutre.

Figure 4-1. Elément de poutre.

4.1 Rigidités des Poutres

Dans cette section, la dérivation de la matrice de rigidité pour un élément de poutre est exposée.
Une poutre est un élément structurel généralement long et mince soumis à une charge transversale qui
produit des effets de flexion importants par rapport à des effets de torsion ou axiaux. Cette déformation de
flexion est mesurée comme un déplacement transversal et une rotation. Dans ce sens, les degrés de liberté
considérés par nœud sont un déplacement transversal et une rotation.

Soit l'élément de poutre illustré dans la figure 4-1. La poutre est de longueur L avec des coordonnées locales
axiales x et des coordonnées locales transversales y. Les déplacements nodaux transversaux dans le repère
local sont donnés par vi et les rotations par i. Les forces nodales locales sont données par fiy et les moments
de flexion par mi comme montré dans la figure 4–1 et 4-2. Nous négligeons tout d'abord tous les effets des
forces axiales.
 Figure 4–1 Elément poutre avec des déplacements, des rotations, des forces, et des moments nodaux
positifs.

Figure 4–2 Conventions de signe pour les forces de cisaillement et les moments de flexion.

(A) poutre non déformée sous charge w (x) (B) poutre déformée en raison du chargement appliqué

(C) élément différentiel de poutre.

Figure 4-3 Poutre sous charge distribuée

Dans tous les nœuds, les conventions de signes suivantes sont utilisées :

1. Les moments sont positifs dans le sens inverse des aiguilles d'une montre.

2. Les rotations sont positives dans le sens inverse des aiguilles d'une montre.

3. Les forces sont positives dans la direction positive.

4. Les déplacements sont positifs dans la direction positive de y.

 La figure 4-2 indique les conventions de signes utilisées dans la théorie de la poutre pour les forces de
cisaillement positives V et les moments de flexion m.

- Matrice de rigidité des poutres basée sur la théorie de la poutre Euler-Bernoulli

(Compte tenu uniquement des déformations de flexion)

L'équation différentielle régissant le comportement élastique linéaire de la poutre est basée sur :

- Au cours de la déformation, les sections droites restent perpendiculaires à la courbe moyenne.

- Les sections droites restent planes selon Navier-Bernoulli (pas de gauchissement).

L'hypothèse de Bernoulli permet de négliger le cisaillement dans le cas de la flexion : le risque de

rupture est alors dû à l'extension des fibres situées à l'extérieur de la poutre, et la flèche est due au moment
fléchissant. Cette hypothèse n'est pas valable pour les poutres courtes, car ces dernières sont hors des limites
de validité du modèle de poutre, à savoir que la dimension des sections doit être petite devant la longueur
de la courbe moyenne. Ceci est illustré dans la figure 4-3, où un plan par la ligne verticale a-c (figure 4-3
(a)) est perpendiculaire à l'axe longitudinal x avant la flexion, et ce même plan par a’-c’ (tournant à travers
l'angle  dans la Figure 4- 3 (b)) reste perpendiculaire à l'axe x courbé après la flexion. Cela se produit en
pratique uniquement lorsqu’un moment constant existe dans la poutre.

L'équation différentielle est dérivée comme suit. Considérons la poutre représentée sur la figure 4-3
soumise à un chargement distribué w(x) (force / longueur). D’après l'équilibre de la force et du moment
d'un élément différentiel de la poutre, montré à la figure 4-3 (c), nous avons

∑ 𝐹𝑦 = 0: 𝑉 − (𝑉 + 𝑑𝑉 ) − 𝑥 (𝑥 )𝑑𝑥 = 0 (4.1.1a)

En simplifiant l'équation (4.1.1a), on obtient

𝑑𝑉
−𝑤𝑑𝑥 − 𝑑𝑉 = 0 𝑜𝑢 𝑤 = − (4.1.1a)
𝑑𝑥

𝑑𝑥 𝑑𝑀
∑ 𝑀2 = 0: − 𝑣𝑑𝑥 + 𝑑𝑀 + 𝑤 (𝑥 )𝑑𝑥 =0 𝑜𝑢 𝑉 = (4.1.1c)
2 𝑑𝑥

La forme finale de l'équation (4.1.1c), reliant la force de cisaillement au moment de flexion, est obtenue
en divisant l'équation gauche par dx puis en prenant la limite de l'équation lorsque dx= 0. Le terme w(x)
disparaît alors.

En outre, la courbure κ de la poutre est liée au moment par

1 𝑀
𝜅= = (4.1.1d)
𝜌 𝐸𝐼
 ρ est le rayon de la courbure comme illustré dans la figure 4-4b, v est la fonction de déplacement
transversal dans la direction y (voir figure 4-4a), E est le module d'élasticité et I est le moment d'inertie
suivent l'axe z (où l'axe z est perpendiculaire aux axes x et y). (I = bh3 / 12 pour une section transversale
rectangulaire de base b et de hauteur h indiquée sur la figure 4-4c).

La courbure pour des petites rotations  = dv / dx est donnée par

𝑑2 𝑣
𝜅= (4.1.1d)
𝑑𝑥 2

En remplaçant l'équation (4.1.1e) dans (4.1.1d), nous obtenons

𝑑2 𝑣 𝑀
= (4.1.1f)
𝑑𝑥 2 𝐸𝐼

On résoudra l'équation (4.1.1f) pour M et en remplaçant ce résultat en (4.1.1c) et (4.1.1b), on obtient

𝑑2 𝑑2 𝑣
2
𝐸𝐼 = −𝑤(𝑥) (4.1.1g)
𝑑𝑥 𝑑𝑥 2

Pour EI constante et seulement des forces et des moments nodaux, l'équation (4.1.1g) devient

𝑑4 𝑣
𝐸𝐼 =0 (4.1.1g)
𝑑𝑥 4

a) Partie de corbure de la poutre b)Rayon de corbure à v(x)

C) Section rectangulaire .

Figure 4-4 .Déformation d’une poutre en flexion

Nous allons maintenant suivre les étapes décrites au chapitre 1 pour développer la matrice de rigidité pour
un élément de la poutre, puis illustrer des solutions complètes pour les poutres.

Étape 1 Sélectionner le type d'élément

Représenter la poutre avec numérotation des nœuds à chaque extrémité et en général numéroter l'élément
(Figure 4-1).

Étape 2 Sélectionner une fonction de déplacement

Supposons que la variation de déplacement transversal à travers la longueur de l'élément est égale à :

𝑣(𝑥 ) = 𝑎1𝑥 3 + 𝑎2 𝑥 2 + 𝑎3 𝑥 + 𝑎4 (4.1.2)

La fonction de déplacement cubique Eq. (4.1.2) est appropriée car il existe quatre degrés de liberté (un
déplacement transversal vi et une petite rotation i à chaque nœud).

En utilisant la même procédure que celle décrite à la section 2.2, nous exprimons v en fonction des degrés
nodaux de liberté v1, v2, 1 et 2 comme suit :

𝑣(0) = 𝑣1 = 𝑎4
𝑑𝑣(0)
= 𝜃1 = 𝑎4
𝑑𝑥
(4.1.3)
𝑣(𝐿) = 𝑣2 = 𝑎1𝐿3 + 𝑎2 𝐿2 + 𝑎3 𝐿 + 𝑎4
𝑑𝑣(0)
= 𝜃1 = 3𝑎1 𝐿2 + 2𝑎2 𝐿 + 𝑎3
𝑑𝑥

= dv/dx pour une petite rotation  supposée. La résolution des équations (4.1.3) pour a1 à a4 en termes
de degrés de liberté nodaux et en remplaçant dans l'équation (4.1.2), nous avons

2 1 −3 1
𝑣= (𝑣1 − 𝑣2) + (𝜃1 + 𝜃2) + (𝑣1 − 𝑣2 ) − (2𝜃1 + 𝜃2 )𝑥 2 + 𝜃1𝑥 + 𝑣1 (4.1.4)
𝐿3 𝐿2 𝐿2 𝐿

Sous forme matricielle, nous exprimons l'équation (4.1.4) comme

𝑣 = [𝑁]{𝑑 } (4.1.5)

Avec

𝑣1
𝜃
{𝑑 } = {𝑣1 } (4.1.6a)
2
𝜃2

[𝑁] = [𝑁1 𝑁2 𝑁3 𝑁4 ] (4.1.6b)

 1 1
𝑁1 = 3
(2𝑥 3 − 3𝑥 2 𝐿 + 𝐿3 ) 𝑁2 = (𝑥 3 𝐿 − 2𝑥 2 𝐿2 + 𝑥𝐿3 )
𝐿 𝐿3
1 1 (4.1.7)
𝑁3 = 3
(−2𝑥 3 + 3𝑥 2𝐿 ) 𝑁4 = 3
(𝑥 3𝐿 − 𝑥 2 𝐿2 )
𝐿 𝐿

N1, N2, N3 et N4 sont appelés les fonctions de forme pour un élément de poutre. Ces fonctions de
forme cubique sont connues sous le nom des fonctions d’interpolation. Pour l'élément de poutre, N1 = 1
lorsqu'il est évalué au nœud 1 et N1 = 0 lorsqu'il est évalué au nœud 2. Parce que N2 est associé à 1, nous
avons, à partir des deux équations (4.1.7), (dN2/ dx)= 1 lorsqu'il est évalué au nœud 1. Les fonctions de
forme N3 et N4 ont des résultats analogues pour le nœud 2.

Étape 3 Définir les relations déformation-déplacement et déformation-contrainte

Soit la relation de déformation axial / déplacement suivante soit valide

𝑑𝑢
𝜀𝑥 (𝑥, 𝑦) = (4.1.8)
𝑑𝑥

Figure 4-5 segment de la poutre (a) avant déformation et (b) après déformation ; (C) Angle de rotation de
la section transversale ABCD

u est la fonction de déplacement axial. À partir de la configuration déformée de la poutre représentée sur
la figure 4-5, on relie le déplacement axial au déplacement transversal par

𝑑𝑣
𝑢 = −𝑦 (4.1.9)
𝑑𝑥

Nous rappelon l'hypothèse fondamentale de poutre, que les sections transversales de la poutre (la section
transversale ABCD) soit planes avant la déformation par flexion restent planes après la déformation et, en
général, avec une petite rotation (dv /dx) . En utilisant l'équation (4.1.9) dans l'équation (4.1.8), nous
obtenons
 𝑑2 𝑢
𝜀𝑥 (𝑥, 𝑦) = −𝑦 (4.1.10a)
𝑑𝑥 2

Utilisant la loi de Hooke (x = E x) et remplaçant l’Eq. (4.1.1f) pour d2v = dx2 dans l'équation (4.1.10a),
on obtient la formule de flexion de la poutre ou la contrainte de flexion comme suit :

−𝑀𝑦
𝜎𝑥 = (4.1.10b)
𝐼

À partir de la théorie de poutre, le moment de flexion et la force de cisaillement sont liés à la fonction de
déplacement transversal comme suit :

𝑑2 𝑣 𝑑3 𝑣
𝑚(𝑥 ) = 𝐸𝐼 𝑉 = 𝐸𝐼 (4.1.11)
𝑑𝑥 2 𝑑𝑥 3

Étape 4 Dérivée la matrice de rigidité

Tout d'abord, dériver la matrice de rigidité de l’élément poutre en utilisant l’approche d'équilibre
direct. Nous relions les forces de cisaillement et les moments de flexion de la poutre (figures 4-1 et 4-2),
ainsi que les équations. (4.1.4) et (4.1.11), pour obtenir

𝑑3 𝑣(0) 𝐸𝐼
𝑓1𝑦 = 𝑣 = 𝐸𝐼 = (12𝑣1 + 6𝐿𝜃1 − 12𝑣2 + 6𝐿𝜃2)
𝑑𝑥 3 𝐿3

𝑑2 𝑣(0) 𝐸𝐼
𝑚1 = −𝑚 = −𝐸𝐼 = (6𝐿𝑣1 + 4𝐿2 𝜃1 − 6𝐿𝑣2 + 2𝐿2 𝜃2) (4.1.12)
𝑑𝑥 2 𝐿3

𝑑3 𝑣(𝐿) 𝐸𝐼
𝑓2𝑦 = −𝑉 = −𝐸𝐼 = (−12𝑣1 − 6𝐿𝜃1 + 12𝑣2 + 6𝐿𝜃2)
𝑑𝑥 3 𝐿3

𝑑2 𝑣(𝐿) 𝐸𝐼
𝑚2 = 𝑚 = 𝐸𝐼 = (6𝐿𝑣1 + 2𝐿2 𝜃1 − 6𝐿𝑣2 + 4𝐿2 𝜃2 )
𝑑𝑥 2 𝐿3

Les équations (4.1.12) relient les forces nodales aux déplacements nodaux. Sous forme matricielle, les
équations (4.1.12) deviennent

𝑓1𝑦 12 6𝐿 −12 6𝐿 𝑣1
𝑚 2
−6𝐿 2𝐿2 ] {𝜃1}
[ 1 ] = 3 [ 6𝐿 4𝐿
𝐸𝐼
(4.1.13)
𝑓2𝑦 𝐿 −12 −6𝐿 12 −6𝐿 𝑣2
𝑚2 6𝐿 2𝐿2 −6𝐿 24 𝜃2

La matrice de rigidité est alors égale

12 6𝐿 −12 6𝐿
2
[𝑘 ] = 3 [ 6𝐿 4𝐿
𝐸𝐼 −6𝐿 2𝐿2 ] (4 .1.14)
𝐿 −12 −6𝐿 12 −6𝐿
6𝐿 2𝐿2 −6𝐿 24
L'équation (4.1.13) indique que [k] relie les forces transversales et les moments de flexion aux
déplacements transversaux et aux rotations, alors que les effets axiaux ont été négligés.

Matrice de rigidité des poutres basée sur la théorie de la poutre Timoshenko (Y compris la
déformation transversale du cisaillement)
La théorie de poutre avec déformation due au cisaillement est dérivée comme suit. la déformation
par cisaillement (déformation due à la force du cisaillement V) est maintenant incluse. En se référant à la
figure 4-6, on observe une section d'une poutre de longueur différentielle dx avec la section transversale
supposée rester plane mais non plus perpendiculaire à l'axe neutre (axe x) en raison de l'inclusion de la
force de cisaillement entraînant une rotation indiquée par le terme b. La courbure totale de la poutre dans
un point x se compose de deux parties, une provoquée par la flexion et une par la force de cisaillement, de
sorte que la rotation de la courbure au point x est donnée par

𝑑𝑣
= 𝜃 (𝑥 ) + 𝛽(𝑥) (4.1.15a)
𝑑𝑥

La rotation due au moment de flexion et due à la force transversale de cisaillement est donnée par (x) et
(x).

Nous supposons que la courbure est linéaire et que la courbure angulaire (pente) est très faible.

La relation entre le moment de flexion et la déformation de flexion (courbure) est comme suit :

𝑑𝜃(𝑥)
𝑀(𝑥 ) = 𝐸𝐼 (4.1.15b)
𝑑𝑥

Et la relation entre la force de cisaillement et la déformation due à cisaillement est donnée par

𝑉 (𝑥 ) = 𝑘𝑠 𝐴𝐺𝛽(𝑥) (4.1.15c)

La différence entre dv /dx et  représenté la déformation du cisaillement zy (=) de la poutre

𝑑𝑣
𝛾𝑦𝑧 = = −𝜃 (4.1.15d)
𝑑𝑥

Soit l'élément différentiel dans la Figure 4-3 (c) et les équations. (4.1.1b) et (4.1.1c) obtenues en sommant
les forces transversales puis en additionnant les moments de flexion. Nous remplaçons maintenant l'éq.
(4.1.15c) pour V et l’Eq. (4.1.15b) pour M dans les équations. (4.1.1b) et (4.1.1c) avec  de l'équation
(4.1.15a) pour obtenir les deux équations différentielles suiventes :

𝑑 𝑑𝑣
𝑘𝑠 𝐴𝐺 − 𝜃 = −𝑤 (4.1.15e)
𝑑𝑥 𝑑𝑥

𝑑 𝑑𝜃 𝑑𝑣
𝐸𝐼 + 𝑘𝑠 𝐴𝐺 −𝜃 = 0 (4.1.15f)
𝑑𝑥 𝑑𝑥 𝑑𝑥
Pour dériver la matrice de rigidité pour l'élément de poutre, y compris la déformation transversale du
cisaillement, nous supposons que le déplacement transversal doit être donné par la fonction cubique dans
l'équation. (4.1.2). Nous choisissons la déformation de cisaillement transversale  consistant avec le
polynôme cubique pour v(x), de sorte que  est une constante donnée par

𝛾=𝑐 (4.1.15g)

En utilisant la fonction de déplacement cubique pour v, la relation de rotation donnée par l'équation
(4.1.15a), et la déformation de cisaillement donnée par l’Eq. (4.1.15g), ainsi que la relation de moment de
flexion -courbure, Eq. (4.1.15b) et la relation de déformation de cisaillement- force de cisaillement Eq.
(4.1.15c), dans la relation de force de cisaillement – moment de flexion Eq. (4.1.1c), nous obtenons

𝑐 = 6𝑎1 𝑔 (4.1.15h)

g= EI / ks AG et ks A est la surface de cisaillement. Les surfaces de cisaillement, As, varient selon les
formes des sections. Par exemple, pour une forme rectangulaire As est égale à 0.83 fois la section
transversale A, pour une section transversale circulaire, on prend 0,9 fois la section transversale.

Utiliser les équations. (4.1.2), (4.1.15a), (4.1.15g) et (4.1.15h) permettent à  d'être exprimé comme un
polynôme en x comme suit:

𝜃 = 𝑎3 + 2𝑎2 𝑥 + (3𝑥 2 + 6𝑔)𝑎1 (4.1.15i)

Utiliser les équations. (4.1.2) et (4.1.15i), nous pouvons exprimer les coefficients a 1 à a4 en termes de
déplacements nodaux v1 et v2 et les rotations 1 et 2 de la poutre aux extrémités x = 0 et x = L comme
précédemment. Les expressions pour a1 à a4 sont ensuite données comme suit :

2𝑣1 +𝐿𝜃1 −2𝑣2 +𝐿𝜃2

𝑎1 =
𝐿(𝐿2 +12𝑔)

−3𝐿𝑣1 −(2𝐿2 +6𝑔)𝜃1 +3𝐿𝑣2 +(−𝐿2 +6𝑔)𝜃2

𝑎2 = (4.1.15j)
𝐿(𝐿2 +12𝑔)

−12𝑔𝑣1 +(𝐿3 +6𝑔𝐿)𝜃1 +12𝑔𝑣2 −6𝑔𝐿𝜃2

𝑎3 =
𝐿(𝐿2 +12𝑔)

𝑎4 = 𝑣1

En substituant dans l'équation (4.1.2), nous obtenons

2𝑣1 +𝐿𝜃1 −2𝑣2 +𝐿𝜃2 3𝐿𝑣1 +(2𝐿2 +6𝑔)𝜃1 −3𝐿𝑣2 −(−𝐿3 +6𝑔)𝜃2
𝑣= 𝑥3 − 𝑥2 −
𝐿(𝐿2 +12𝑔) 𝐿(𝐿2 +12𝑔)

12𝑔𝑣1 −(𝐿3 +6𝑔𝐿)𝜃1 −12𝑔𝑣2 −6𝑔𝐿𝜃2

𝑥 + 𝑣1 (4.1.15k)
𝐿(𝐿2 +12𝑔)
D'une manière similaire à l'étape 4 (pour dériver la matrice de rigidité pour l'élément de poutre sans prendre
la déformation du cisaillement), nous avons

𝐸𝐼(12𝑣1 +6𝐿𝜃1 −12𝑣2 +6𝐿𝜃2 )

𝑓1𝑦 = 𝑣(0) = 6𝐸𝐼𝑎1 =
𝐿(𝐿2 +12𝑔)

𝐸𝐼(−6𝐿𝑣1 −(4𝐿2 +12𝑔)𝜃1 +6𝐿𝑣2 +(−2𝐿2 +12𝑔)𝜃2 )

𝑚1 = −𝑚(0) = −2𝐸𝐼𝑎2 = (4.1.15l)
𝐿(𝐿2 +12𝑔)

𝐸𝐼(−12𝑣1 −6𝐿𝜃1 +12𝑣2 +6𝐿𝜃2 )

𝑓2𝑦 = −𝑉(𝐿) =
𝐿(𝐿2 +12𝑔)

𝐸𝐼[6𝐿𝑣1 +(2𝐿2 −12𝑔)𝜃1 −6𝐿𝑣2 +(4𝐿2 +12𝑔)𝜃2

𝑚2 = 𝑚(𝑙) =
𝐿(𝐿2 +12𝑔)

Sous forme matricielle l’Eq. (4 . 1. 15 l) devient

𝑓1𝑦 12 6𝐿 −12 6𝐿 𝑣1
2 2
𝑚 𝐸𝐼 6𝐿 (4𝐿 + 12𝑔) −6𝐿 (2𝐿 − 12𝑔) 𝜃1
[ 1] = 2 [ ] {𝑣 } (4.1.15m)
𝑓2𝑦 𝐿(𝐿 +12𝑔) −12 −6𝐿 12 −6𝐿 2
2 2 𝜃
𝑚2 6𝐿 (2𝐿 − 12𝑔) −6𝐿 (4𝐿 + 12𝑔) 2

La matrice de rigidité, y compris la déformation de flexion et cisaillement, est alors donnée par

12 6𝐿 −12 6𝐿
2 2
𝐸𝐼 6𝐿 (4𝐿 + 12𝑔) −6𝐿 (2𝐿 − 12𝑔)
[𝑘 ] = 2 [ ] (4 .1.15n)
𝐿(𝐿 +12𝑔) −12 −6𝐿 12 −6𝐿
2
6𝐿 (2𝐿 − 12𝑔) −6𝐿 (4𝐿2 + 12𝑔)

Dans l'équation (4.1.15n) rappelez-vous que g représente le terme de cisaillement transversal, et si on fixe
g = 0, on obtient l'équation (4.1.14) pour la matrice de rigidité de la poutre en négligeant la déformation du
cisaillement transversal. Pour voir plus facilement l'effet du facteur de correction de cisaillement, nous
définissons le terme de correction de cisaillement non dimensionnel
𝜑=12EI/(ksAGL2)=12g=L2 et réécrivons la matrice de rigidité comme

12 6𝐿 −12 6𝐿
2 2
𝐸𝐼 6𝐿 (4𝐿 + 𝜑) −6𝐿 (2𝐿 − 𝜑)
[𝑘 ] = 2 [ ] (4 .1.15o)
𝐿(𝐿 +𝜑) −12 −6𝐿 12 −6𝐿
6𝐿 (2𝐿2 − 𝜑) −6𝐿 (4𝐿2 + 𝜑)

4.2 Étape 5 Assemblage des équations des éléments pour obtenir les équations globales et introduire
des conditions aux limites
 Considérons la poutre dans la figure 4-6 comme exemple pour illustrer la procédure d'assemblage
des matrices de rigidité des éléments de poutre. Supposons que EI est constant dans toute la poutre. Une
force de 1000 N et un moment de 1000 N-mm sont appliqués au milieu de la poutre. L'extrémité gauche
est un support fixe et l'extrémité droite est un support de rotule simple.

Figure 4-6 Poutre soumise à force et moment.

Les charges sont supposées être appliquées uniquement aux nœuds. En utilisant l'équation (4.1.14), les
matrices de rigidité globale pour les deux éléments sont maintenant données par

12 6𝐿 −12 6𝐿
2
[𝑘 (1) ] = 3 [ 6𝐿 4𝐿
𝐸𝐼 −6𝐿 2𝐿2 ] (4 .2.1)
𝐿 −12 −6𝐿 12 −6𝐿
6𝐿 2𝐿2 −6𝐿 24

12 6𝐿 −12 6𝐿
2
[𝑘 (2) ] = 3 [ 6𝐿 4𝐿
𝐸𝐼 −6𝐿 2𝐿2 ] (4 .2.2)
𝐿 −12 −6𝐿 12 −6𝐿
6𝐿 2𝐿2 −6𝐿 24

La matrice de rigidité globale peut être assemblée en utilisant la méthode de rigidité directe. Lorsque la
matrice de rigidité globale a été assemblée, les forces nodales globales externes sont liées aux déplacements
nodaux globaux. Grâce à la superposition directe et aux équations. (4.2.1) et (4.2.2), les équations de
régulation pour la poutre sont donc données par

𝐹1𝑦 12 6𝐿 −12 6𝐿 0 0
𝑀1 6𝐿 4𝐿2 −6𝐿 2𝐿2 0 0 𝑣1
𝐹2𝑦 𝐸𝐼 −12 −6𝐿 12 + 12 −6𝐿 + 6𝐿 −12 6𝐿 𝜃
= 3 {𝑣1 } (4.2.3)
𝑀2 𝐿 6𝐿 2𝐿2 −6𝐿 + 6𝐿 4𝐿2 + 4𝐿2 −6𝐿 2𝐿 2
𝐹3𝑦 0 0 −12 −6𝐿 12 −6𝐿 𝜃2
[ 𝑀3 ] [ 0 0 6𝐿 2𝐿2 −6𝐿 4𝐿 2 ]

Compte tenu des conditions aux limites, avec le nœud 1 fixe et un appui simple au noeud 3, nous avons

𝜃1 = 0 𝑣1 = 0 𝑣3 = 0 (4.2.4)
En considérant les troisième, quatrième et sixième équations dans la matrice (4.2.3) correspondant aux
linges avec des degrés de liberté inconnus et utilisant des équations. (4.2.4), et remplaçant F2y = 1000 N,
M2 = 1000 N-MM et M3 = 0. nous obtenons

−1000 24 0 6𝐿 𝑣2
𝐸𝐼
{ 1000 } = 3 [ 0
𝐿
8𝐿2 2𝐿2 ] {𝜃2 } (4.2.5)
0 6𝐿 2𝐿2 4𝐿2 𝜃3
Nous allons maintenant résoudre l'équation. (4.2.5) pour trouver le déplacement nodal inconnu v2 et les
rotations nodales inconnues 2 et 3.

4.3 Vecteur Charges exprimé en repère local

Figure 4.7 –Elément de poutre plane chargé uniformément.

Considérant une poutre plane chargée uniformément, le vecteur de charge {fe} se réduit à celui des forces
nodales de surface (4.18) soit

2 3
𝐹1𝑦 1 − 3𝑥2 + 2𝑥2
𝐿 𝐿
𝑀1 2𝑥2 𝑥3
0 𝐹2𝑦 𝑥 − + 0
𝐿
{𝑓𝑒 } = {𝑓𝑒𝑠 } = ∫𝑠 [𝑁]𝑇 . {𝑓 𝑠 }. 𝑑𝑆𝑒 = ∫0 [𝑁]𝑇 . {
𝐿 𝐿 𝐿2
−𝑞 } 𝑑𝑥 = = ∫0 2 3 . {−𝑞 } 𝑑𝑥 =
𝑒 𝑦 𝑀2 3𝑥 2𝑥 𝑦
2 + 3
𝐹3𝑦 𝐿 𝐿
𝑥2 𝑥3
[ 𝑀3 ] { − +
𝐿 𝐿3 }
𝑞𝑦 𝐿
− 2
𝑞𝑦 𝐿2
− 12
𝑞𝑦 𝐿
(6.3.4)
− 2
𝑞𝑦 𝐿2
{ 12 }

D’une manière générale et dans le cas de forces orientées dans le sens y négatif, le vecteur de charge d’un
élément de poutre plane soumise à un système de charges transversales est égale à l’inverse des réactions
et moments d’encastrement de la poutre bi-encastrée subissant le même chargement. En d’autres termes, ce
vecteur traduit les actions nodales équivalentes au chargement appliqué sur la poutre (cf. tableau 6.1 des
charges nodales équivalentes).