Calcul des Structures par la Méthode des

Éléments Finis
Deuxième Année DUT Génie Mécanique - Maroc

Mohamed Diouri
Année Universitaire 2025-2026

“La précision de l’analyse numérique dépend de la qualité de la discrétisation géométrique et

de la pertinence du choix des fonctions d’approximation.”

1
TABLE DES MATIÈRES S3 DUT Génie Mécanique

Table des matières

1 Introduction Générale 4
1.1 Applications Pratiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

2 Rappel de Résistance des Matériaux 4

2.1 Hypothèses Fondamentales de la Calcul des Structures . . . . . . . . . . . . . . 4
2.2 Efforts Intérieurs et Déformations Unitaires . . . . . . . . . . . . . . . . . . . . 4
2.3 Équations Différentielles du Poutre . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.4 Méthodes Classiques de Résolution . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.4.1 Méthode des Différentielles . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.4.2 Méthode des Efforts Intérieurs . . . . . . . . . . . . . . . . . . . . . . . . 5

3 Fonctions d’Approximation : Lagrange et Hermite 6

3.1 Principe de l’Approximation dans le Méthode des Éléments Finis . . . . . . . . 6
3.2 Fonctions de Lagrange . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.2.1 Élément Linéaire à 2 Nœuds . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.2.2 Élément Quadratique à 3 Nœuds . . . . . . . . . . . . . . . . . . . . . . 7
3.3 Fonctions d’Hermite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.3.1 Élément de Poutre Cubique à 2 Nœuds . . . . . . . . . . . . . . . . . . . 7
3.4 Dérivées et Déformations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.4.1 Dérivées des Fonctions de Lagrange . . . . . . . . . . . . . . . . . . . . . 7
3.4.2 Dérivées des Fonctions d’Hermite . . . . . . . . . . . . . . . . . . . . . . 8
3.5 Critères de Choix des Fonctions d’Approximation . . . . . . . . . . . . . . . . . 8

4 Introduction à la Méthode des Éléments Finis 8

4.1 Principe Général du Méthode des Éléments Finis . . . . . . . . . . . . . . . . . 8
4.2 Formulation Variationnelle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
4.3 Équation Matricielle du Méthode des Éléments Finis . . . . . . . . . . . . . . . 9

5 Changement de Repère en Méthode des Éléments Finis : Concepts Généraux 9

5.1 Nécessité du Changement de Repère . . . . . . . . . . . . . . . . . . . . . . . . 9
5.2 Principe de la Transformation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
5.3 Cas Plan vs Cas Spatial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

6 Éléments de Barre : Formulation et Changement de Repère 10

6.1 Modélisation de l’Élément de Barre . . . . . . . . . . . . . . . . . . . . . . . . . 10
6.2 Formulation dans le Repère Local . . . . . . . . . . . . . . . . . . . . . . . . . . 10
6.3 Changement de Repère : Barre Plan . . . . . . . . . . . . . . . . . . . . . . . . 11
6.4 Matrice de Rigidité Globale : Barre . . . . . . . . . . . . . . . . . . . . . . . . . 12
6.5 Cas Particuliers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
6.5.1 Barre Horizontale (α = 0) . . . . . . . . . . . . . . . . . . . . . . . . . . 12
6.5.2 Barre Verticale (α = 90) . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
6.5.3 Barre à 45° (α = 45) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

7 Éléments de Poutre : Formulation et Changement de Repère 12

7.1 Modélisation de l’Élément de Poutre . . . . . . . . . . . . . . . . . . . . . . . . 12
7.2 Formulation dans le Repère Local . . . . . . . . . . . . . . . . . . . . . . . . . . 13
7.2.1 Extension Axiale (Lagrange Linéaire) . . . . . . . . . . . . . . . . . . . . 13
7.2.2 Flexion (Hermite Cubique) . . . . . . . . . . . . . . . . . . . . . . . . . 13
7.3 Matrice de Rigidité Locale : Poutre . . . . . . . . . . . . . . . . . . . . . . . . . 13

2
TABLE DES MATIÈRES S3 DUT Génie Mécanique

7.4 Changement de Repère : Poutre Plan . . . . . . . . . . . . . . . . . . . . . . . . 14

7.5 Matrice de Rigidité Globale : Poutre . . . . . . . . . . . . . . . . . . . . . . . . 14

8 Combinaison Barre-Poutre : Modélisation Hybride 15

8.1 Nécessité des Modèles Mixtes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
8.2 Problématique de Compatibilité des DDL . . . . . . . . . . . . . . . . . . . . . 15
8.3 Stratégie de Modélisation Hybride . . . . . . . . . . . . . . . . . . . . . . . . . . 15
8.3.1 Numérotation Unifiée des DDL . . . . . . . . . . . . . . . . . . . . . . . 15
8.3.2 Matrices de Transformation Adaptées . . . . . . . . . . . . . . . . . . . . 16
8.3.3 Matrice de Rigidité Élémentaire Hybride . . . . . . . . . . . . . . . . . . 16
8.4 Assemblage de la Matrice Globale . . . . . . . . . . . . . . . . . . . . . . . . . . 16
8.5 Exemple Numérique : Portique avec Contreventement . . . . . . . . . . . . . . . 16
8.5.1 Calcul des Matrices Locales . . . . . . . . . . . . . . . . . . . . . . . . . 17
8.5.2 Transformation et Assemblage . . . . . . . . . . . . . . . . . . . . . . . . 17
8.5.3 Système Global et Résolution . . . . . . . . . . . . . . . . . . . . . . . . 17
8.6 Vérification de la Conservation d’Énergie . . . . . . . . . . . . . . . . . . . . . . 17

9 Synthèse et Perspectives 18
9.1 Bilan du Cours . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
9.2 Limites de l’Approche et Extensions . . . . . . . . . . . . . . . . . . . . . . . . 18
9.3 Applications Professionnelles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
9.4 Recommandations pour la Suite . . . . . . . . . . . . . . . . . . . . . . . . . . . 18

10 Annexes 19
10.1 Formules Synthétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
10.1.1 Éléments de Barre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
10.1.2 Éléments de Poutre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
10.2 Exercices d’Application . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
10.2.1 Exercice 1 : Barre Incliné . . . . . . . . . . . . . . . . . . . . . . . . . . 19
10.2.2 Exercice 2 : Poutre Cantilever . . . . . . . . . . . . . . . . . . . . . . . . 19
10.2.3 Exercice 3 : Portique Simple . . . . . . . . . . . . . . . . . . . . . . . . . 19
10.2.4 Exercice 4 : Structure Hybride . . . . . . . . . . . . . . . . . . . . . . . 20

11 Introduction 20

12 Repères local et global 20

13 Élément barre (barre axiale en 2D) 20

13.1 Degrés de liberté . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
13.2 Matrice de transformation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
13.3 Matrice de raideur locale (barre 2 noeuds) . . . . . . . . . . . . . . . . . . . . . 21

14 Élément poutre 2D (poutre d’Euler–Bernoulli + axiale) 21

14.1 Degrés de liberté . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
14.2 Matrice de raideur locale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
14.3 Matrice de transformation pour la poutre 2D . . . . . . . . . . . . . . . . . . . . 22
14.4 Raideur globale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
14.5 Remarques pratiques et implémentation numérique . . . . . . . . . . . . . . . . 22

15 Exemple compact (formule explicite pour la barre) 23

16 Conclusion 23

3
2 RAPPEL DE RÉSISTANCE DES MATÉRIAUX S3 DUT Génie Mécanique

1 Introduction Générale
Ce cours de calcul des structures s’adresse aux étudiants de deuxième année DUT Génie
Mécanique et constitue une transition essentielle entre la Calcul des Structures classique et les
méthodes numériques modernes. Il vise trois objectifs principaux :
1. Rappeler et consolider les concepts fondamentaux de la Calcul des Structures
2. Introduire les outils mathématiques nécessaires au Méthode des Éléments Finis
(fonctions d’approximation)
3. Maîtriser les éléments de base du Méthode des Éléments Finis pour les structures
filaires (barres et poutres)

1.1 Applications Pratiques

Les concepts développés s’appliquent dans de nombreux domaines du génie mécanique :

— Conception de châssis automobiles — Calcul de ponts et ouvrages d’art

— Dimensionnement d’ossatures métal-
— Optimisation de machines-outils
liques
— Analyse de structures aéronautiques — Simulation de treuils et palans

2 Rappel de Résistance des Matériaux

2.1 Hypothèses Fondamentales de la Calcul des Structures
La Calcul des Structures classique repose sur un ensemble d’hypothèses simplificatrices qui
permettent de passer du comportement local du matériau à celui de la structure :
1. Homogénéité et isotropie : Propriétés constantes dans tout le volume
2. Linéarité : Relation σ = Eε valable dans la limite d’élasticité
3. Petites déformations : |ε| ≪ 1, rotations négligeables
4. Sections planes : Hypothèse de Saint-Venant pour les poutres
5. Fibres normales : Hypothèse de Bernoulli-Navier (pas de déformation de cisaillement)

2.2 Efforts Intérieurs et Déformations Unitaires

Effort intérieur Symbole Déformation Relation constitutive

Effort normal N εx = ∆L
L
N = EAεx
Moment de flexion y My κz = ρ1z My = EIz κz
Moment de flexion z Mz κy = ρ1y Mz = EIy κy
Moment de torsion Mt γxy = rθ
L
Mt = GJ Lθ
Ty
Effort tranchant y Ty γxz = GA s
Ty = GAs γxz
Tz
Effort tranchant z Tz γyz = GA s
Tz = GAs γyz

Table 1 – Relations effort-déformation en Calcul des Structures

4
2 RAPPEL DE RÉSISTANCE DES MATÉRIAUX S3 DUT Génie Mécanique

2.3 Équations Différentielles du Poutre

Pour une poutre soumise à des charges volumiques q = (qx , qy , qz ) :

dN
+ qx =0 (1a)
dx
dTy d2 Mz
− + qy =0 (1b)
dx dx2
dTz d2 My
+ + qz =0 (1c)
dx dx2
dMt
+ mt =0 (1d)
dx
En combinant les équations (1b) et (1c) avec l’hypothèse de Bernoulli, on obtient les équa-
tions classiques :

d4 v d4 w
EIz = qy (x), EIy = qz (x) (2)
dx4 dx4

2.4 Méthodes Classiques de Résolution

2.4.1 Méthode des Différentielles
Résolution des équations différentielles d’ordre 4 avec 4 conditions aux limites :

∫ x
1
2
v(x) = A1 + A2 x + A3 x + A4 x + 3
(x − ξ)qy (ξ)dξ
EIz 0
+ termes particuliers selon qy (x) (3)

2.4.2 Méthode des Efforts Intérieurs

Pour les structures isostatiques ou hyperstatiques :
1. Établir les équilibres statiques aux nœuds
2. Écrire les relations de continuité des déplacements
3. Appliquer les conditions aux appuis
4. Résoudre le système algébrique

N (x) P M (x) T (x) Q

RA , M A RB RC

Figure 1 – Poutre hyperstatique - Méthode des efforts intérieurs

5
3 FONCTIONS D’APPROXIMATION : LAGRANGE ET HERMITE
S3 DUT Génie Mécanique

3 Fonctions d’Approximation : Lagrange et Hermite

3.1 Principe de l’Approximation dans le Méthode des Éléments Fi-
nis
Le principe fondamental du Méthode des Éléments Finis consiste à approximer les champs
de déplacements continus par des fonctions simples définies sur des sous-domaines discrets
(éléments finis). Cette approximation s’écrit :
∑
n
u(x) ≈ uh (x) = Ni (x)ui = N(x)ue (4)
i=1

Où :
— N(x) : Matrice des fonctions d’approximation
— ue : Vecteur des déplacements nodaux
— n : Nombre de nœuds de l’élément

3.2 Fonctions de Lagrange

Les fonctions de Lagrange sont polynomiales et assurent la continuité des déplacements (C 0 )
entre éléments adjacents.

3.2.1 Élément Linéaire à 2 Nœuds

Pour un élément de longueur L avec nœuds aux extrémités x1 = 0, x2 = L :

L−x
N1 (x) = 1 − ξ = (5)
L
x
N2 (x) = ξ = (6)
L
Où ξ = x/L est la coordonnée réduite (0 ≤ ξ ≤ 1).
Propriétés :
— N1 (0) = 1, N1 (L) = 0, N2 (0) = 0, N2 (L) = 1
— N1 (x) + N2 (x) = 1 (partition de l’unité)
— Continuité C 0 : Déplacements continus, dérivées discontinues

Ni
1
N1 (x) = 1 − x
L

x
N2 (x) = L
0 x
x1 = 0 x2 1= L

Figure 2 – Fonctions de Lagrange linéaires - Élément à 2 nœuds

6
3 FONCTIONS D’APPROXIMATION : LAGRANGE ET HERMITE
S3 DUT Génie Mécanique

3.2.2 Élément Quadratique à 3 Nœuds

Pour un élément avec nœuds aux positions x1 = 0, x2 = L/2, x3 = L :

N1 (x) = 2ξ 2 − 3ξ + 1 (7)
N2 (x) = −4ξ 2 + 4ξ (8)
N3 (x) = 2ξ 2 − ξ (9)

Ces fonctions permettent une approximation parabolique plus précise que l’approximation
linéaire.

3.3 Fonctions d’Hermite

Les fonctions d’Hermite sont utilisées pour les problèmes où les dérivées des déplacements
(rotations, courbures) sont des degrés de liberté. Elles assurent la continuité C 1 (déplacements
et dérivées continus).

3.3.1 Élément de Poutre Cubique à 2 Nœuds

Pour un élément de poutre de longueur L avec 2 DDL par nœud (déplacement transversal
v et rotation θ = dv/dx), l’approximation s’écrit :

v(x) = N1 (x)v1 + N2 (x)Lθ1 + N3 (x)v2 + N4 (x)Lθ2 (10)

Les fonctions hermitiennes cubiques sont :

N1 (x) = 1 − 3ξ 2 + 2ξ 3 (11a)
N2 (x) = ξ − 2ξ 2 + ξ 3 (11b)
N3 (x) = 3ξ 2 − 2ξ 3 (11c)
N4 (x) = −ξ 2 + ξ 3 (11d)

Propriétés essentielles :

Fonction Ni (0) Ni′ (0) Ni (L) Ni′ (L)

N1 (x) 1 0 0 0
N2 (x) 0 1/L 0 1/L
N3 (x) 0 0 1 0
N4 (x) 0 0 0 1/L

Table 2 – Propriétés des fonctions hermitiennes aux nœuds

3.4 Dérivées et Déformations

3.4.1 Dérivées des Fonctions de Lagrange
Pour l’élément linéaire :
dN1 1 dN2 1
=− , = (12)
dx L dx L
duh u 2 − u1
εx = = = constante (13)
dx L

7
4 INTRODUCTION À LA MÉTHODE DES ÉLÉMENTS FINIS S3 DUT Génie Mécanique

Ni

vh (x) (approximation)
N1 (x) N3 (x)
x
v1 , θ1 Élément poutre hermitien v2 , θ2

Figure 3 – Fonction hermitienne cubique - Déplacement transversal

3.4.2 Dérivées des Fonctions d’Hermite

Pour la flexion (élément cubique) :
dN1 1
= (−6ξ + 6ξ 2 ) (14a)
dx L
dN2 1
= (1 − 4ξ + 3ξ 2 ) (14b)
dx L
d 2 N1 6
= 2 (−1 + 2ξ) (14c)
dx2 L
d 2 N2 2
= 2 (−2 + 3ξ) (14d)
dx2 L
La courbure (déformation de flexion) est :
d 2 vh 1
κz = − 2
= − 2 [6(−1 + 2ξ)v1 + · · · ] (15)
dx L

3.5 Critères de Choix des Fonctions d’Approximation

Critère Lagrange Linéaire Lagrange Quadratique Hermite Cubique

Continuité C0 C0 C1
Précision Moyenne Bonne Excellente
DDL par nœud 1 1 2
Applications Barres, membranes Barres complexes Poutres, plaques
Complexité Faible Moyenne Élevée
Coût calculatoire Faible Moyen Élevé

Table 3 – Comparaison des fonctions d’approximation

4 Introduction à la Méthode des Éléments Finis

4.1 Principe Général du Méthode des Éléments Finis
La méthode des éléments finis transforme un problème continu en un problème matriciel
discret par trois étapes principales :
1. Discrétisation géométrique : Division du domaine Ω en sous-domaines Ωe
2. Approximation locale : Expression des solutions sur chaque élément
3. Assemblage global : Combinaison des contributions élémentaires

8
5 CHANGEMENT DE REPÈRE EN MÉTHODE DES ÉLÉMENTS FINIS : CONCEPTS
GÉNÉRAUX S3 DUT Génie Mécanique

Élément fini e
Domaine continu Ω

Ni (x)
Maillage grossier (3 éléments)

Figure 4 – Principe de discrétisation en Méthode des Éléments Finis

4.2 Formulation Variationnelle

Le problème mécanique s’énonce par le principe des travaux virtuels :
∫ ∫ ∫
δε : σ dΩ = δ · b dΩ + δ · t dS (16)
Ω Ω ∂Ωt

Avec la relation constitutive σ = Dε (tenseur d’élasticité).

4.3 Équation Matricielle du Méthode des Éléments Finis

En introduisant l’approximation h = Ne et εh = Be, l’équation faible devient :
∫ (∫ ∫ )
T T T T T
δe B DBe dΩ = δe N b dΩ + N t dS (17)
Ωe Ωe ∂Ωe

On obtient le système élémentaire :

K e e =e (18)

Où :
∫
Ke = BT DB dΩ (19a)
∫ Ωe
∫
T
e = N b dΩ + NT t dS (19b)
Ωe ∂Ωe

L’assemblage global donne :

K= (20)

5 Changement de Repère en Méthode des Éléments Fi-

nis : Concepts Généraux
5.1 Nécessité du Changement de Repère
Dans les structures réelles, les éléments ne sont généralement pas alignés avec les axes du
repère global de l’ouvrage. Le changement de repère est nécessaire pour :
1. Formuler les relations constitutives dans le repère local de l’élément
2. Respecter l’anisotropie géométrique de chaque élément
3. Assembler correctement les contributions élémentaires dans le repère global
4. Appliquer les conditions aux limites dans le repère global

9
6 ÉLÉMENTS DE BARRE : FORMULATION ET CHANGEMENT
S3 DUT
DE REPÈRE
Génie Mécanique

5.2 Principe de la Transformation

Le principe fondamental du changement de repère repose sur la conservation de l’énergie
mécanique. Pour un déplacement virtuel δ̆ :

δW = δggT = δllT (21)

Avec la transformation g = Tl, on obtient :

g= TTl , K g = TT K l T (22)

Où T est la matrice de transformation (matrice de rotation).

5.3 Cas Plan vs Cas Spatial

Aspect Cas Plan (2D) Cas Spatial (3D)

Degrés de liberté 2 translations 3 translations + 3 rotations
Matrice de rotation 2×2 3 × 3 ou 6 × 6
Angle de rotation α (1 paramètre) 3 angles d’Euler
Complexité Faible Élevée
Applications Portiques plans, treillis 2D Structures 3D complexes

Table 4 – Comparaison plan/spatial

Ce cours se concentre sur le cas plan, plus accessible pour le niveau DUT.

6 Éléments de Barre : Formulation et Changement de

Repère
6.1 Modélisation de l’Élément de Barre
L’élément de barre est le plus simple en Méthode des Éléments Finis. Il modélise un membre
axialement porteur sans résistance à la flexion.
Caractéristiques :
— 2 nœuds aux extrémités
— 2 DDL par nœud en 2D : uxi , uyi et uxj , uyj
— 4 DDL totaux par élément
∆L
— Déformation uniquement axiale : εx = L

6.2 Formulation dans le Repère Local

Dans le repère local (xl , yl ) aligné avec l’axe de la barre :
Degrés de liberté locaux : uli , ulj (déplacements axiaux)
Approximation : ( xl ) xl
ul (xl ) = 1 − u l i + u lj (23)
L L
Déformation axiale :
dul u l − u li
εl = = j (24)
dxl L

10
6 ÉLÉMENTS DE BARRE : FORMULATION ET CHANGEMENT
S3 DUT
DE REPÈRE
Génie Mécanique

Yg

u xj , u yj
jxl
yl

u xi , u yi α L sin α

Xg
i L cos α

Figure 5 – Élément de barre incliné - Repères local et global

Énergie de déformation :
∫ L
1 1 EA
Ue = EAε2l dxl = (ulj − uli )2 (25)
2 0 2 L

Matrice de rigidité locale :

[ ]
EA 1 −1
Kbl = (26)
L −1 1

6.3 Changement de Repère : Barre Plan

Pour passer du repère global (Xg , Yg ) au repère local (xl , yl ), on utilise l’angle α entre l’axe
de la barre et l’axe Xg .
Matrice de rotation élémentaire :
[ ]
cos α sin α
R= (27)
− sin α cos α

Transformation des déplacements :

{ } { }
u li u
= R xi (28)
v li u yi

Puisque la barre ne résiste qu’aux déplacements axiaux, on projette :

uli = cos α uxi + sin α uyi (29)

Matrice de transformation complète (4×2) :

[ ]
cos α sin α 0 0
Tb = (30)
0 0 cos α sin α

Vérification de l’orthogonalité :
[ ]
1 0
Tb TTb = (31)
0 1

11
7 ÉLÉMENTS DE POUTRE : FORMULATION ET CHANGEMENT
S3 DUT
DE Génie
REPÈREMécanique

6.4 Matrice de Rigidité Globale : Barre

En appliquant la transformation :

Kbg = TTb Kbl Tb (32)

Développement :
[ ] [ ]
EA cos α cos α sin α 0 0
Kbg = [1 − 1]
L 0 0 0 cos α sin α
[ ] [ ]
EA 0 0 0 cos α sin α
+ [1 − 1]
L cos α cos α sin α 0 0
Résultat final :
 
cos2 α cos α sin α − cos2 α − cos α sin α
EA 
 cos α sin α sin2 α − cos α sin α − sin2 α  
Kbg = (33)
L  − cos α 2
− cos α sin α 2
cos α cos α sin α 
− cos α sin α − sin2 α cos α sin α sin2 α

6.5 Cas Particuliers

6.5.1 Barre Horizontale (α = 0)
 
1 0 −1 0
EA  0 0 0
Kbg = 0  (34)
L −1 0 1 0
0 0 0 0

6.5.2 Barre Verticale (α = 90)

 
0 0 0 0
EA  0 −1
Kbg = 0 1  (35)
L 0 0 0 0
0 −1 0 1

6.5.3 Barre à 45° (α = 45)

 
1 1 −1 −1
EA 
1 1 −1 −1
Kbg =  (36)
2L −1 −1 1 1
−1 −1 1 1

7 Éléments de Poutre : Formulation et Changement de

Repère
7.1 Modélisation de l’Élément de Poutre
L’élément de poutre modélise la flexion et l’extension axiale. Il est plus complexe que l’élé-
ment de barre.
Caractéristiques :

12
7 ÉLÉMENTS DE POUTRE : FORMULATION ET CHANGEMENT
S3 DUT
DE Génie
REPÈREMécanique

— 2 nœuds aux extrémités

— 3 DDL par nœud en 2D : uxi , uyi , θzi
— 6 DDL totaux par élément
— Déformations : εx (axiale) et κz (flexion)

Yg

yl jxul xj , uyj , θzj

θ zi
u yi α
u xi
zl Xg
i

Figure 6 – Élément de poutre incliné - Repères et DDL

7.2 Formulation dans le Repère Local

Dans le repère local (xl , yl , zl ) :
Degrés de liberté : {uli , vli , θli , ulj , vlj , θlj }

7.2.1 Extension Axiale (Lagrange Linéaire)

( xl ) xl
ul (xl ) = 1 − u l i + u lj (37)
L L
ulj −uli
Déformation : εl = L

7.2.2 Flexion (Hermite Cubique)

vl (xl ) = N1 (xl )vli + N2 (xl )Lθli + N3 (xl )vlj + N4 (xl )Lθlj (38)
Avec les fonctions (11a)-(11d). La courbure est :
d 2 vl
κ lz = − (39)
dx2l

7.3 Matrice de Rigidité Locale : Poutre

La matrice de rigidité locale découle de l’énergie totale :
∫ L
1 EA 1
Ue = Uaxial + Uf lexion = (ulj − uli )2 + EIz κ2lz dxl (40)
2 L 2 0
Contribution axiale :
 
1 0 0 −1 0 0
0 0 0 0 0 0
 
EA 
0 0 0 0 0 0 
Kpl,axial = (41)
L 
−1 0 0 1 0 0 
0 0 0 0 0 0
0 0 0 0 0 0

13
7 ÉLÉMENTS DE POUTRE : FORMULATION ET CHANGEMENT
S3 DUT
DE Génie
REPÈREMécanique

Contribution flexion (calcul intégral) :

 
0 0 0 0 0 0
0 12 6L 0 −12 6L 
 
p EIz 
 0 6L 4L2 0 −6L 2L2 

Kl,f lexion = 3  (42)
L 0 0 0 0 0 0 
0 −12 −6L 0 12 −6L
0 6L 2L2 0 −6L 4L2

Matrice totale :
Kpl = Kpl,axial + Kpl,f lexion (43)

7.4 Changement de Repère : Poutre Plan

La transformation pour la poutre doit traiter les translations et les rotations.
Matrice de rotation pour les translations :
 
cos α sin α 0
Rtrans = − sin α cos α 0 (44)
0 0 1

Transformation des déplacements nodaux :

   
 u li   u xi 
vl = Rtrans uyi (45)
 i  
θ li θ zi

Matrice de transformation complète (6×6) :

 
c s 0 0 0 0
−s c 0 0 0 0
[ ]  
Rtrans 03×3 0 0 1 0 0 0
Tp = =
0
 (46)
03×3 Rtrans  0 0 c s 0
0 0 0 −s c 0
0 0 0 0 0 1

Où c = cos α, s = sin α.

7.5 Matrice de Rigidité Globale : Poutre

Kpg = TTp Kpl Tp (47)

Cette opération donne une matrice 6×6 symétrique et définie positive. Le calcul explicite
est complexe mais peut être structuré par blocs.
Cas particulier : Poutre horizontale (α = 0)

Kpg = Kpl (48)

Cas particulier : Poutre verticale (α = 90)
La matrice se transforme selon la rotation de 90° des axes.

14
8 COMBINAISON BARRE-POUTRE : MODÉLISATION HYBRIDE
S3 DUT Génie Mécanique

8 Combinaison Barre-Poutre : Modélisation Hybride

8.1 Nécessité des Modèles Mixtes
Dans la pratique, la plupart des structures réelles combinent des éléments purement axiaux
(barres de treillis) et des éléments fléchissants (poutres de portique). Cette combinaison est
indispensable pour :
— Représenter fidèlement l’architecture des structures industrielles
— Optimiser le calcul en utilisant l’élément adapté à chaque fonction
— Modéliser les contreventements et les éléments porteurs

D C

E
Éléments barre (contreventement)

A Éléments poutre (portique) B

P = 50 kN

Figure 7 – Structure hybride : Portique avec contreventement en treillis

8.2 Problématique de Compatibilité des DDL

La principale difficulté technique réside dans la compatibilité des degrés de liberté aux
nœuds communs :

Élément DDL Nœud i DDL Nœud j Total Interprétation physique

Barre uxi , uyi u xj , u yj 4 Translations uniquement
Poutre u xi , u y i , θ z i u xj , u y j , θ z j 6 Translations + Rotations

Table 5 – Comparaison des DDL - Barre vs Poutre

Aux nœuds de connexion :

— Les translations ux , uy doivent être identiques pour les deux types
— La rotation θz n’existe pas pour la barre
— Les efforts se cumulent selon les conventions de signe

8.3 Stratégie de Modélisation Hybride

8.3.1 Numérotation Unifiée des DDL
Pour gérer cette hétérogénéité, on utilise une numérotation globale qui assigne :
— 3 DDL par nœud pour les nœuds de poutres : {ux , uy , θz }
— 2 DDL par nœud pour les nœuds de barres uniquement : {ux , uy }
— DDL conditionnels aux interfaces : θz = 0 pour les barres

15
8 COMBINAISON BARRE-POUTRE : MODÉLISATION HYBRIDE
S3 DUT Génie Mécanique

Nœud Connexion DDL Attribués Index Global Commentaire

A Poutre uniquement u xA , u yA , θ z A 1, 2, 3 DDL complets
B Poutre + Barre u xB , u y B , θ z B 4, 5, 6 Rotation active
C Poutre uniquement uxC , uyC , θzC 7, 8, 9 DDL complets
E Barre uniquement u xE , u yE 10, 11 Pas de rotation

Table 6 – Numérotation des DDL dans une structure hybride

8.3.2 Matrices de Transformation Adaptées

Pour les éléments de poutre (6 DDL) : Utilisation de Tp (6×6)
Pour les éléments de barre (4 DDL) : Utilisation de Tb (4×2) avec sélection des
translations uniquement :
[ ]
ubarre = I2×2 02×1 unud (49)
Où unud = {ux , uy , θz }T représente les 3 DDL du nœud.

8.3.3 Matrice de Rigidité Élémentaire Hybride

Pour une barre connectée à un nœud de poutre, la matrice devient :
[ ]T [ ]
Kbhybride = I2×2 02×1 Kbg I2×2 02×1 (50)
Cette matrice 6×6 (pour un nœud complet) a des termes nuls dans les lignes/colonnes
correspondant à θz .

8.4 Assemblage de la Matrice Globale

L’assemblage suit le schéma classique mais avec une gestion spéciale des interfaces :
∑
Kg = Keg avec Keg de taille adaptée (51)
e

Exemple pour la structure de la figure 7 :

— Élément AB (poutre) : Contribution 6×6 aux DDL 1-6
— Élément BE (barre) : Contribution 6×6 aux DDL 4-5, 10-11 (avec θzB , θzE = 0)
— Élément BC (poutre) : Contribution 6×6 aux DDL 4-6, 7-9
La matrice globale résultante est 11×11 (9 DDL poutres + 2 DDL barre).

8.5 Exemple Numérique : Portique avec Contreventement

Considérons le portique renforcé de la figure 7 avec les données suivantes :
Paramètres géométriques :
√
— LAB = LBC = 6 m, LBE = 6 2 m
— αAB = 0, αBC = 90, αBE = 45
Propriétés des sections :
— Poutres : Ap = 20 cm², Ip = 800 cm�, E = 200 GPa
— Barre : Ab = 10 cm²
Charge appliquée : P = 50 kN verticale au nœud B

16
8 COMBINAISON BARRE-POUTRE : MODÉLISATION HYBRIDE
S3 DUT Génie Mécanique

8.5.1 Calcul des Matrices Locales

Élément AB (poutre horizontale) :
 
3.33 0 0 −3.33 0 0
 0 2.22 0.30 0 −2.22 0.30 
 
 0
8 0.30 0.080 0 −0.30 0.040 
p
Kl,AB = 10   N/m (52)
 −3.33 0 0 3.33 0 0  
 0 −2.22 −0.30 0 2.22 −0.30
0 0.30 0.040 0 −0.30 0.080
Élément BE (barre diagonale) :
[ ]
1 −1
Kbl,BE = 2.36 × 10 8
(53)
−1 1

8.5.2 Transformation et Assemblage

Transformation poutre AB (α = 0) : Kpg,AB = Kpl,AB
Transformation poutre BC (α = 90) :
 
0 1 0
R90◦ = −1 0 0 (54)
0 0 1
Transformation barre BE (α = 45) :
 
0.5 0.5 0 −0.5 −0.5 0
 0.5 0.5 0 −0.5 −0.5 0
 

8 0 0 0 0 0 0
Kg,BE = 1.18 × 10 
b  (55)
−0.5 −0.5 0 0.5 0.5 0 
−0.5 −0.5 0 0.5 0.5 0
0 0 0 0 0 0

8.5.3 Système Global et Résolution

Après assemblage et application des conditions aux limites (uxA = uyA = θzA = 0, uxC =
uyC = 0), le système réduit est 5×5.
Vecteur de charge réduit : r = {−50000, 0, 0, 0, 0}T N
Résultats typiques :
— Déplacement vertical nœud B : uyB ≈ −2.1 mm
— Rotation nœud B : θzB ≈ −0.12 rad
— Déformation barre BE : εBE ≈ 0.18%

8.6 Vérification de la Conservation d’Énergie

La transformation hybride conserve l’énergie totale :

uTg Kg ug = uTl Kl ul (56)

Cette propriété est vérifiée par :
TT K l T = K g ⇒ TT T = I (57)
Pour les structures hybrides, cette propriété s’étend aux interfaces en conservant la compa-
tibilité cinématique.

17
9 SYNTHÈSE ET PERSPECTIVES S3 DUT Génie Mécanique

9 Synthèse et Perspectives
9.1 Bilan du Cours
Ce cours a permis de maîtriser les concepts fondamentaux du calcul des structures par
Méthode des Éléments Finis au niveau DUT :
1. Bases théoriques de la Calcul des 4. Éléments de barre : Modélisation et
Structures et passage au numérique transformation
2. Fonctions d’approximation de La-
5. Éléments de poutre : Flexion et cou-
grange et Hermite
plage axial
3. Principes du Méthode des Élé-
ments Finis et formulation variation- 6. Structures hybrides : Combinaison et
nelle compatibilité

9.2 Limites de l’Approche et Extensions

Modélisation développée Extensions pour les niveaux

supérieurs
Éléments 1D linéaires Éléments 2D/3D (membranes,
coques, solides)
Petites déformations Grandes déformations géomé-
triques
Matériau linéaire élastique Plasticité, hyperélasticité
Statique linéaire Dynamique, instabilités
Structures filaires Structures continues
Changement de repère plan Transformations spatiales com-
plexes

Table 7 – Limites et perspectives du Méthode des Éléments Finis

9.3 Applications Professionnelles

Les compétences acquises s’appliquent directement dans :
— Bureaux d’études : Dimensionnement de structures métalliques
— Industrie automobile : Calcul de châssis et suspensions
— Génie civil : Analyse de bâtiments et ponts
— Aéronautique : Modélisation de fuselages et ailes
— Machines-outils : Optimisation des structures porteuses

9.4 Recommandations pour la Suite

1. Pratique logicielle : Maîtriser un code Méthode des Éléments Finis commercial (Robot
Structural, SAP2000)
2. Vérification analytique : Toujours comparer avec les solutions Calcul des Structures
classiques
3. Analyse de convergence : Étudier l’effet du raffinement du maillage
4. Projets réels : Appliquer ces concepts à des études industrielles concrètes

18
10 ANNEXES S3 DUT Génie Mécanique

10 Annexes
10.1 Formules Synthétiques
10.1.1 Éléments de Barre
Matrice locale : [ ]
EA 1 −1
Kbl = (58)
L −1 1
Matrice globale :  
c2 cs −c2 −cs
EA 
 cs2 s
2
−cs −s2 

Kbg = (59)
L  −c −cs c2 cs 
−cs −s2 cs s2

10.1.2 Éléments de Poutre

Matrice locale flexion :
 
12 6L −12 6L
EI  6L 4L2 −6L 2L2 
Kpl,f lexion = 3   (60)
L −12 −6L 12 −6L
6L 2L2 −6L 4L2

10.2 Exercices d’Application

10.2.1 Exercice 1 : Barre Incliné
Une barre d’acier (E = 200 GPa, A = 5 cm²) de longueur L = 3 m est inclinée à α = 30.
Une force Fx = 20 kN est appliquée à l’extrémité.
1. Écrivez la matrice de rigidité globale
2. Calculez les déplacements aux extrémités
3. Vérifiez l’allongement axial

10.2.2 Exercice 2 : Poutre Cantilever

Une poutre en béton armé (E = 30 GPa, I = 2000 cm�) de 4 m de long, encastrée d’une
part, soumise à P = 10 kN à l’extrémité libre.
1. Comparez la solution Méthode des Éléments Finis (1 élément) avec la solution analytique
2. Étudiez l’effet d’un maillage plus fin
3. Analysez la distribution des moments fléchissants

10.2.3 Exercice 3 : Portique Simple

Un portique rectangulaire 3 m × 2 m avec une charge horizontale de 15 kN au linteau.
1. Modélisez avec 2 éléments poutres
2. Calculez les déplacements et rotations
3. Vérifiez l’équilibre des efforts aux nœuds

19
13 ÉLÉMENT BARRE (BARRE AXIALE EN 2D) S3 DUT Génie Mécanique

10.2.4 Exercice 4 : Structure Hybride

Portique avec une diagonale de contreventement. Comparez :
1. Le déplacement sans contreventement
2. Le déplacement avec contreventement (barre)
3. L’effet de rigidification apporté par la diagonale
Changement de repère pour l’élément réel : barre et poutre (repère plan)

11 Introduction
Ce chapitre explique la procédure et les matrices nécessaires pour ramener la raideur d’un
élément fini défini dans son repère local vers le repère global (repère plan 2D). On donne les
matrices de transformation usuelles pour l’élément barre (axial) et pour l’élément poutre 2D
(axial + flexion), ainsi que la formule générale
K (g) = T T K (l) T
où K (l) est la matrice de raideur dans le repère local et K (g) sa forme dans le repère global.

Packages recommandés
% dans le préambule de [Link]
\usepackage{amsmath,amssymb,mathtools}

12 Repères local et global

Soit un élément raide posé entre deux nœuds 1 et 2. On note :
— (X, Y ) le repère global.
— (x, y) le repère local de l’élément tel que l’axe x est dirigé de 1 vers 2 (axe longitudinal
de l’élément).
— L la longueur de l’élément, θ l’angle entre l’axe local x et l’axe global X (mesuré positi-
vement à partir de X vers Y ).
Les coordonnées du vecteur directeur unité sont
c = cos θ, s = sin θ.

13 Élément barre (barre axiale en 2D)

13.1 Degrés de liberté
Pour une barre 2D placée dans un plan, on peut conserver pour chaque nœud les degrés de
liberté globaux (uX , uY ) (déplacements en X et Y ). Le comportement axial ne dépend que de
la projection des déplacements le long de l’axe de l’élément.
Ordre des degrés de liberté local choisi (vectoriel) :
{u1 , u2 }local
où ui est le déplacement axial local au nœud i. Dans le repère global on a l’ordre usuel
{u1X , u1Y , u2X , u2Y }.

20
14 ÉLÉMENT POUTRE 2D (POUTRE D’EULER–BERNOULLIS3
+ DUT
AXIALE)
Génie Mécanique

13.2 Matrice de transformation

La transformation qui relie les déplacements locaux aux déplacements globaux s’écrit :
 
( ) u 1X
u1  u1Y 
= Tbarre 
u2X 

u2
u2Y
avec ( )
c s 0 0
Tbarre = .
0 0 c s

13.3 Matrice de raideur locale (barre 2 noeuds)

Dans le repère local, la matrice de raideur axiale (2x2) est
( )
(l) EA 1 −1
Kbarre = ,
L −1 1

où E est le module d’Young et A la section.

La matrice globale (4x4) s’obtient par :
(g) T (l)
Kbarre = Tbarre Kbarre Tbarre .

Si l’on effectue le calcul, on obtient la forme connue

 2 
c cs −c2 −cs
(g) EA  cs2 s
2
−cs −s2 
.
Kbarre =
L  −c −cs c2 cs 
−cs −s2 cs s2

14 Élément poutre 2D (poutre d’Euler–Bernoulli + axiale)

14.1 Degrés de liberté
Pour la poutre 2D (élément beam-axial combiné) on admet généralement pour chaque nœud
trois degrés de liberté :

{u, v, φ} où u = dépl. axial, v = dépl. transverse, φ = rotation (en plan).

L’ordre nodal global conventionnel pour l’élément à deux nœuds est :

{u1X , v1Y , φ1 , u2X , v2Y , φ2 }.

14.2 Matrice de raideur locale

Dans le repère local (axes x longitudinal et y transverse), la matrice de raideur combinée
′ ′
(6×6)pourl( lmentpoutre−−axial(Euler−−Bernoulli)s
) critenblocs(sparantlapartaxialeetlapartf lexion)
(l) K a 0
Kpoutre = , où Ka est la contribution axiale (2×2) et Kb la contribution flexion/transverse
0 Kb
(4 × 4) :

21
14 ÉLÉMENT POUTRE 2D (POUTRE D’EULER–BERNOULLIS3
+ DUT
AXIALE)
Génie Mécanique

Axial : ( )
EA 1 −1
Ka = .
L −1 1
Flexion (Euler–Bernoulli) : en notation usuelle (nœuds 1 et 2) :
 
12 6L −12 6L
EI  6L 4L 2
−6L 2L2 
Kb = 3  ,
L −12 −6L 12 −6L
6L 2L2 −6L 4L2

avec I le moment d’inertie.

Remarque : l’ordre dans Kb suit la séquence des degrés de liberté transverses/rotations :
{v1 , φ1 , v2 , φ2 }.

14.3 Matrice de transformation pour la poutre 2D

On construit la matrice de rotation R (6×6)[Link]
pour les déplacements et un coefficient scalaire pour les rotations (les rotations en plan se trans-
forment de la même façon que les directions tangentielles) :
 
R2 0 0
Tpoutre =  0 R2 0  ,
0 0 R2

mais pour expliciter l’ordre {u1 , v1 , φ1 , u2 , v2 , φ2 } on écrit explicitement (bloc par bloc) :
 
c s 0 0 0 0
−s c 0 0 0 0
 
 0 0 1 0 0 0
Tpoutre =   0 0 0 c s 0 .

 
 0 0 0 −s c 0
0 0 0 0 0 1

Explication : les déplacements axiaux/transverses se transforment par la rotation 2D [c s; −s c],

tandis que la rotation scalaire locale φ est la même variable numérique (dans le cas plan) —
on conserve donc le coefficient 1 pour la composante rotation.

14.4 Raideur globale

La matrice globale s’obtient par
(g) T (l)
Kpoutre = Tpoutre Kpoutre Tpoutre .

Le calcul donne une matrice 6×6compltequicouplelescomposantesaxialesetf lexionnellessuivantladirectio

en partition 3×3deblocs2×2 et effectuer la transformation bloc par bloc.

14.5 Remarques pratiques et implémentation numérique

— Dans une routine d’assemblage d’un code EF, calculez d’abord K (l) (dans l’ordre local
choisi), construisez T à partir de c = cos θ, s = sin θ, puis effectuez K (g) = T T K (l) T .
— Vérifiez la symétrie et la positivité définie de K (g) (numériquement) pour détecter des
erreurs de construction.

22
17 MATRICES DE RAIDEUR DANS LE REPÈRE GLOBAL S3 DUT Génie Mécanique

— Attention au sens positif de θ et à l’ordonnancement des degrés de liberté : une inversion

d’ordre conduit à des signes erronés ou à des couplages inversés.
— Si vous utilisez des éléments avec des DOF supplémentaires (par ex. cisaillement ou
Timoshenko), la taille des blocs et la forme de K (l) et T changent mais le principe K (g) =
T T K (l) T reste identique.

15 Exemple compact (formule explicite pour la barre)

(g)
Pour rappel : la matrice 4×4Kbarre explicitée plus haut est pratique pour implémentation
directe :  2 
c cs −c2 −cs
(g) EA  cs2 s
2
−cs −s2 
.
Kbarre = 
L −c −cs c2 cs 
−cs −s2 cs s2

16 Conclusion
Ce chapitre fournit les formules standards pour ramener les matrices de raideur locales (barre
et poutre) au repère global. Dans un code EF, c’est l’étape clé avant l’assemblage globale. Les
mêmes principes s’appliquent en 3D en remplaçant les matrices de rotation 2D par des matrices
de rotation 3D adaptées.

17 Matrices de raideur dans le repère global

17.1 Barre 2D (axiale)
La matrice de raideur globale d’un élément barre incliné dans le plan (X, Y ) s’écrit :
 
c2 cs −c2 −cs
 
(g) EA 
cs s2 −cs −s2 
,
Kbarre =
L  −c 2
−cs c 2
cs 

−cs −s 2
cs s 2

avec c = cos θ et s = sin θ.

17.2 Poutre 2D (axiale + flexion)

Pour l’élément poutre plan (Euler–Bernoulli), avec les degrés de liberté

{u1 , v1 , φ1 , u2 , v2 , φ2 },

la matrice de raideur globale complète est :

 
K11 K12 K13 K14 K15 K16
K12 K22 K23 K24 K25 K26 
 
K13 K23 K33 K34 K35 K36 
Kpoutre =  ,
(g)
K14 K24 K34 K44 K45 K46 
 
K15 K25 K35 K45 K55 K56 
K16 K26 K36 K46 K56 K66

23
17 MATRICES DE RAIDEUR DANS LE REPÈRE GLOBAL S3 DUT Génie Mécanique

où les coefficients s’écrivent explicitement :

( ) ( )
E AL2 c2 + 12Is2 Ecs AL2 − 12I 6EIs
K11 = 3
, K12 = 3
, K13 = − 2 ,
L L L
( 2 2 )
E AL s + 12Ic2 6EIc 4EI
K22 = 3
, K23 = 2
, K33 = ,
L L L
K14 = −K11 , K15 = −K12 , K16 = K13 ,
K24 = −K12 , K25 = −K22 , K26 = K23 ,

K34 = −K13 , K35 = −K23 , K36 = 2EI

L
,
K44 = K11 , K45 = K12 , K46 = −K13 ,
K55 = K22 , K56 = −K23 , K66 = K33 .
(g) (g)
Cette forme explicite illustre la symétrie de la matrice Kpoutre = (Kpoutre )T .

24