UNIVERSITE DES FRERES MENTOURI CONSTANTINE

FACULTE DES SCIENCES DE LA TECHNOLOGIE

DEPARTEMENT DE GENI CIVIL

Auteur : DOCTEUR KHEBIZI MOURAD

MAITRE DE CONFERENCES « A »

COURS DE LA METHODE
DES ELEMENTS FINIS
Destiné aux étudiants de :
 Master I : géotechnique

Année universitaire 2019-2020

Cours MEF Introduction

Etabli par : Dr. KHEBIZI Mourad

2
Cours MEF

INTRODUCTION

Ce cours a pour objet d’initier l'étudiant aux notions de base de la méthode des éléments
finis. L’enseignement de ce cours nécessite des connaissances préalables à savoir l’analyse
mathématique, calcul matriciel, résistance des matériaux, méthodes numériques et mécanique
des milieux continus.

Etabli par : Dr. KHEBIZI Mourad

3
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

CHAPITRE 2
METHODE DES ELEMENTS FINIS 1D

2.1. INTRODUCTION

L’étude des structures en barres et poutres par la méthode des déplacements contient une
bonne part d’une analyse par éléments finis et conduit au système linéaire [1] :

[𝐾0 ]{𝑞} = {𝐹}

Où [𝐾0 ]est la matrice de rigidité, {𝑞} le vecteur des 𝑛𝑑 inconnues cinématiques (déplacement)
ou degrés de liberté et {𝐹} le vecteur force.
La méthode des éléments finis 1D consiste à discrétiser les structures par des éléments finis
barres, des éléments finis poutres et des éléments finis portiques. La formulation de chaque type
de ces éléments est présentée dans les sections suivantes de ce chapitre.

2.2. LES ELEMENTS BARRES :

Les éléments finis barres sont utilisés pour discrétiser les structures et les matériaux sollicitées
axialement.

2.2.1. ELEMENT BARRE A DEUX NŒUDS

Nous considérons l’élément barre à deux nœuds de la figure 2.1. Cette barre est sollicitée
par deux forces axiales 𝐹1 et 𝐹2 . Dans ce cas chaque nœud comporte un seul degré de liberté
dans le repère local (une seule inconnue cinématique : 𝑛𝑑 = 1) et la barre comporte deux
inconnues cinématiques (𝑛𝑑 = 2) dans le repère local ( 𝑥, 𝑦).
𝑥, 𝑦: Système d’axes local
𝑋, 𝑌: Système d’axes global

𝑃𝑜𝑖𝑛𝑡

𝑁𝑜𝑒𝑢𝑑

𝛼 𝑁𝑜𝑒𝑢𝑑
𝑋
Figure 2.1 : Elément barre à deux nœuds

Etabli par : Dr. KHEBIZI Mourad

4
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

La matrice de rigidité de cet élément (élément à barre à 2 nœuds) dans le repère local
s’exprime comme suit.
𝐸𝐴 1 −1
[𝐾0 ] = [ ]
𝑙 −1 1
Avec : 𝐸 est le module d’élasticité de barre, 𝐴 la section de la barre et 𝑙 la longueur de la
barre.

NB : La matrice de rigidité de tous les éléments finis doit être symétrique.

Le vecteur des 𝑛𝑑 inconnues cinématiques (déplacements) de l’élément barre à deux nœuds

s’écrit dans le repère local par :
𝑢
{𝑞𝑒 } = {𝑢1 }
2

et le vecteur force de cet élément s’écrit dans le repère local par :

𝐹
{𝐹𝑒 } = { 1 }
𝐹2

2.2.2. ELEMENT BARRE A TROIS NŒUDS

La figure suivante présente un élément barre à trois nœuds, dans lequel chaque nœud
comporte un seul degré de liberté dans le repère local (une seule inconnue cinématique : 𝑛𝑑 =
1) et la barre comporte trois inconnues cinématiques (𝑛𝑑 = 3) dans le repère local ( 𝑥, 𝑦). Ces
trois inconnues cinématiques sont 𝑢1 , 𝑢2 , 𝑢3 .

1 𝑢1 2 𝑢2 3 𝑢3
𝑥
−𝑙 +𝑙
Figure 2.2 : élément barre à trois nœuds

La matrice de rigidité de cet élément (élément à barre à 3 nœuds) dans le repère local
s’exprime par.
7 8 1
−
6 6 6
𝐸𝐴 8 16 8 𝐸𝐴 7 −8 1
[𝐾0 ] = − − = [−8 16 −8]
𝑙 6 6 6 6𝑙
1 −8 7
1 8 7
[ 6 −6 6 ]
N.B : la longueur de la barre égale 2𝑙.

Etabli par : Dr. KHEBIZI Mourad

5
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

Le vecteur des 𝑛𝑑 inconnues cinématiques (déplacement) de l’élément barre à trois nœuds

s’écrit dans le repère local par :
𝑢1
{𝑞𝑒 } = {𝑢2 }
𝑢3
et le vecteur force de cet élément s’écrit dans le repère local par :
𝐹1
{𝐹𝑒 } = {𝐹2 }
𝐹3
Remarque : 𝐾0 symétrique

2.2.3. CALCUL D’UNE STRUCTURE A BARRES :

Pour calculer une structure à barres par la méthode des éléments finis il faut suivre les étapes
suivantes (figure 2.3):
 Nous choisissons un repère global,
 Nous proposons une numérotation globale (numérotation des nœuds et des éléments) ;
 Nous déduisons une connexion des éléments (la connectivité) (voir tableau 2.1).

10 9 8 7 6
9 8 7 6
10

14
12

16

1 5
1 2 3 4 4 𝑋
2 3
Figure 2.3 : Structure constituée par plusieurs éléments barres

Tableau 2.3 : connexion des éléments (la connectivité)

1𝑒𝑟 nœuds 2𝑒𝑚𝑒 nœuds
Elément 1 1 2
Elément 12 2 9
Elément 15 3 7

 Nous tirons les coordonnées de tous les nœuds (0,0), (3,0), (6,0), … … ..
 Nous calculons la longueur de chaque élément :
Longueur de l’élément = 𝑙 = √(𝑥2 − 𝑥1 )2 + (𝑦2 − 𝑦1 )2
 Nous calculons les matrices de rigidité élémentaires dans les repères locaux et dans le
repère global (matrices de rigidité de chaque élément).
 Nous calculons la matrice de rigidité de la structure dans le repère global.
 Nous résolvons le système linéaire [𝐾0 ]{𝑞} = {𝐹}.

Etabli par : Dr. KHEBIZI Mourad

6
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

2.2.4. PASSAGE AU REPERE GLOBAL :

La matrice de passage (matrice de rotation) permet le passage des degrés de liberté globaux aux
degrés de liberté locaux, elle s’écrit comme suit (valable pour l’élément barre à deux nœuds):

[𝑅𝑂] = [𝑐𝑜𝑠𝛼 𝑠𝑖𝑛𝛼 0 0

]
0 0 𝑐𝑜𝑠𝛼 𝑠𝑖𝑛𝛼
Les composantes locales et globales des degrés de liberté sont liées par la transformation.

{𝑞}𝐿 = [𝑅𝑂]{𝑞}𝐺
Pour l’élément barre à deux nœuds, chaque nœud comporte un seul degré de liberté dans le
𝑈1
𝑢1 𝑉
repère local et deux degrés de liberté dans le repère global : {𝑞}𝐿 = {𝑢 } et {𝑞}𝐺 = { 1 }
2 𝑈2
𝑉2
[𝑅𝑂]−1 = [𝑅𝑂]𝑇 , elle permet le passage des degrés de liberté locaux aux degrés de liberté
globaux
Les composantes locales et globales des forces nodales sont liées par la transformation

{𝐹𝑒 }𝐿 = [𝑅𝑂]{𝐹𝑒 }𝐺
𝐹𝑋1
𝐹 𝐹
Avec : {𝐹𝑒 }𝐿 = { 1 } et {𝐹𝑒 }𝐺 = { 𝑌1 }
𝑓2 𝐹𝑋2
𝐹𝑌2
La matrice de rigidité dans le repère global s’exprime par
[𝐾0 ]𝐺 = [𝑅𝑂]𝑇 [𝐾0 ]𝑙 [𝑅𝑂]
𝐸𝐴 +1 −1
Avec [𝐾0 ]𝑙 est la matrice de rigidité dans le repère local : [𝐾0 ]𝑙 = [ ]
𝑙 −1 +1

2.2.5. VECTEUR CHARGE EQUIVALENTE :

Si l’élément support une charge distribué sur sa longueur (figure 2.4), le vecteur force (vecteur
charge équivalente) s’exprime dans le repère local comme suit :
𝑞𝑙/2
{𝐹𝑒 } = { }
𝑞𝑙/2

𝑞(𝑥)
𝑥
1 2
Figure 2.4 : Elément barre à deux nœuds support une charge distribué sur sa longueur

Etabli par : Dr. KHEBIZI Mourad

7
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

2.3. LES ELEMENTS POUTRES EN FLEXION SIMPLE :

Une poutre c’est un corps prismatique obtenu par la translation rigide d’une section
quelconque.
Les éléments poutres sont utilisés pour discrétiser les poutres sollicitées à la flexion simple.

2.3.1. Poutre en flexion simple :

Les théories des poutres sont basées sur plusieurs hypothèses simplificatrices.
Les théories classiques les plus connues et les plus utilisés sont :
- Théorie des poutres de Bernoulli,
- Théorie des poutres de Timoshenko.

2.3.2. Hypothèses
Hypothèse 1 :
Les sections droites avant déformation restent droites après déformation (pas de
gauchissement).
Hypothèse 2 :
- Bernoulli : La normale à une section avant déformation reste normale après
déformation.
Poutres minces (poutres élancées).
- Timoshenko : La normale à une section avant déformation ne reste plus normale à
cette section (mais elle reste droite).
Poutre d’épaisseur modérée.

N.B : Dans ce cours, nous allons étudier uniquement les poutres basées sur la cinématique de
Bernoulli.

2.3.3. Poutre basée sur la cinématique de Bernoulli :

Nous considérons un tronçon d’une poutre sollicitée à la flexion (figure 2.5). Ce tronçon doit
discrétiser par un élément poutre (élément fini) à deux nœuds (figure 2.6), dans e cas chaque
nœud comporte deux degrés de liberté (une translation et une rotation) et l’élément comporte
deux inconnues cinématiques (𝑛𝑑 = 4) dans le repère ( 𝑥, 𝑦). Le vecteur des 𝑛𝑑 inconnues
cinématiques (déplacements) de l’élément poutre à 2 nœuds s’exprime dans le repère local par :
𝑤1
𝜃
{𝑞𝑒 } = {𝑤1 }
2
𝜃2

Etabli par : Dr. KHEBIZI Mourad

8
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

𝑦 𝑞

𝑥
Figure 2.5 : Tronçon d’une poutre sollicité à la flexion simple

1 2
𝑥
𝑤1 𝑤2

𝜃1 𝜃2

Figure 2.6 : degrés de liberté d’un élément poutre à 2 nœuds

La matrice de rigidité de l’élément poutre à deux nœuds (𝑛𝑑 × 𝑛𝑑 ) s’écrit comme suit :

𝑤1 𝜃2 𝑤2 𝜃2
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼
− 3 ⟶ 𝑤1
𝑙3 𝑙2 𝑙 𝑙2
6𝐸𝐼 4𝐸𝐼 6𝐸𝐼 2𝐸𝐼
− 2 ⟶ 𝜃1
[𝐾0 ] = 𝑙2 𝑙 𝑙 𝑙
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼
− 3 − 2 − 2 ⟶ 𝑤2
𝑙 𝑙 𝑙3 𝑙
6𝐸𝐼 2𝐸𝐼 6𝐸𝐼 4𝐸𝐼
− 2 ⟶ 𝜃2
[ 𝑙2 𝑙 𝑙 𝑙 ]
Avec 𝑙 est la longueur de l’élément poutre et 𝐼 son inertie.
Si l’élément supporte une charge transversale uniforme 𝑞 (figure 2.7), on obtient

𝑞𝑙/2
𝑞𝑙 2 /12
{𝐹𝑒 } =
𝑞𝑙/2
{−𝑞𝑙 2 /12}

𝑦 𝑞

𝑥
1 2
Figure 2.7 : Elément supporte une charge transversale uniforme

2.3.4. Poutre basée sur la cinématique de Timoshenko :

Il est demandé aux étudiants de faire un exposé sur les éléments poutres basées sur
l’hypothèse de Timoshenko.

Etabli par : Dr. KHEBIZI Mourad

9
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

2.3.5. ELEMENT PORTIQUE :

Les éléments portique sont utilisés pour discrétiser les poutres (structures) sollicitées à la
flexion composées (chargement axial + un moment fléchissant). La figure 2.8 présente un
élément portique à deux nœuds, chaque nœud comporte trois degrés de liberté dans le repère
local (trois inconnues cinématiques ) et l’élément comporte six inconnues cinématiques (𝑛𝑑 =
6) dans le repère local ( 𝑥, 𝑦). Le vecteur des 𝑛𝑑 inconnues cinématiques (déplacements) de
l’élément portique à deux nœuds s’exprime dans le repère local par :

𝑢1
𝑣1
𝜃
{𝑞𝑒 } = 𝑢1
2
𝑣2
{𝜃2 }

𝑣1 𝑣2

1 (0,0) 2 (𝑙, 0) 𝑥
𝑢1 𝑢2
𝜃1 𝜃2

Figure 2.8: Elément portique à deux degrés de liberté dans le plan.

2.3.5.1. Matrice de rigidité de l’élément portique à deux nœuds (dans le repère local et
dans le plan) :
La matrice de rigidité de l’élément portique s’obtient par une combinaison entre la matrice
de rigidité de l’élément barre et celle de l’élément poutre.

𝐸𝐴 1 −1
⟶ 𝑀𝑒𝑛𝑏𝑟𝑎𝑛𝑎𝑖𝑟𝑒 ⟶ [𝐾0 ] = [ ] ⟶ 𝐸𝑙é𝑚𝑒𝑛𝑡 𝑏𝑎𝑟𝑟𝑒
𝑙 −1 1

12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼

−
𝑙3 𝑙2 𝑙3 𝑙2
[𝐾0 ] → 6𝐸𝐼 4𝐸𝐼 6𝐸𝐼 2𝐸𝐼
− 2
⟶ 𝑓𝑙𝑒𝑥𝑖𝑜𝑛𝑛𝑒𝑙𝑙𝑒 ⟶ [𝐾0 ] = 𝑙2 𝑙 𝑙 𝑙2 ⟶ 𝐸𝑙é𝑚𝑒𝑛𝑡 𝑝𝑜𝑢𝑡𝑟𝑒
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼
− 3 − 2 − 2
𝑙 𝑙 𝑙3 𝑙
6𝐸𝐼 2𝐸𝐼 6𝐸𝐼 4𝐸𝐼
[ [ 𝑙2 − 2
𝑙 𝑙 𝑙 ]

Etabli par : Dr. KHEBIZI Mourad

10
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

[𝐾0 ]𝑝𝑜𝑟𝑡𝑖𝑞𝑢𝑒 = [𝐾0 ]𝑏𝑎𝑟𝑟𝑒 + [𝐾0 ]𝑝𝑜𝑢𝑡𝑟𝑒

N.B : le signe (+) de l’expression précédente ne signifie pas la somme algébrique des deux
matrice mais une combinaison entre les matrice.

La matrice de rigidité ((𝑛𝑑 × (𝑛𝑑 ) de l’élément portique à deux nœuds dans le repère local
s’écrite comme suit :

𝐸𝐴 𝐸𝐴
0 0 − 0 0 → 𝑢1
𝑙 𝑙
𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐼
0 12 6 0 −12 6
𝑙3 𝑙2 𝑙2 𝑙 2 → 𝑣1
𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐴
0 6 2 4 0 −6 2 2
[𝐾0 ] = 𝑙 𝑙 𝑙 𝑙 → 𝜃1
𝐸𝐴 𝐸𝐴
− 0 0 0 0
𝑙 𝑙 → 𝑢2
𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐼
0 −12 3 −6 2 0 12 3 −6 2
𝑙 𝑙 𝑙 𝑙 → 𝑣2
𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐼
[ 0 6 2
𝑙
2
𝑙
0 −6 2
𝑙
4
𝑙 ] → 𝜃2

2.3.5.2. Passage au repère global

La matrice de passage (matrice de rotation) permet le passage des degrés de liberté globaux
(repère global 𝑋, 𝑌 ) aux degrés de liberté locaux (repère local 𝑥, 𝑦 ) elle s’écrit comme
suit (valable pour l’élément portique à deux nœuds dans le plan (figure 2.9 ).

𝛼
𝑋
Figure 2.9 : Elément portique à deux nœuds dans le plan : 𝑋, 𝑌 est le repère global et 𝑥, 𝑦
le repère lovcal

La matrice de rotation du passage local-global est la suivante (pour 1 seul nœud) :

Etabli par : Dr. KHEBIZI Mourad

11
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

cos 𝛼 sin 𝛼 0
[𝑅𝑂] = [−sin 𝛼 cos 𝛼 0]
0 0 1

Les composantes locales et globales des degrés de liberté, pour un seul nœud, sont liées par la
transformation :

𝑢 𝑈
𝑣
{ }= [𝑅𝑂] {𝑉 }
𝜃 𝜃

La matrice du passage local-global pour l’élément portique (à deux nœuds) est la suivante :

cos 𝛼 sin 𝛼 0 ⋮ 0 0 0
−sin 𝛼 cos 𝛼 0 ⋮ 0 0 0
0 0 1 ⋮ 0 0 0
[ℛ𝒪] = ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ⋯
6×6 0 0 0 ⋮ cos 𝛼 sin 𝛼 0
0 0 0 ⋮ − sin 𝛼 cos 𝛼 0
[ 0 0 0 ⋮ 0 0 1]

Les composantes locales et globales des degrés de liberté sont liées par la transformation.

{𝑞}𝐿 = [ℛ𝒪]{𝑞}𝐺
Pour l’élément portique à deux nœuds, chaque nœud comporte trois degrés de liberté dans le
𝑢1 𝑈1
𝑣1 𝑉1
𝜃1 𝜃1
repère local et trois degrés de liberté dans le repère global : {𝑞}𝐿 = 𝑢 et {𝑞}𝐺 =
2 𝑈2
𝑣2 𝑉2
{𝜃2 } { 𝜃2 }
N.B : la matrice [𝑅𝑂]−1 permet le passage des degrés de liberté locaux aux degrés de liberté
globaux.

Les composantes locales et globales des forces nodales de l’élément portique sont liées par la
transformation

{𝐹𝑒 }𝐿 = [𝑅𝑂]{𝐹𝑒 }𝐺
La matrice de rigidité de l’élément portique dans le repère global s’exprime par :

[𝐾0 ]𝑔 = [ℛ𝒪]𝑇 [𝐾0 ]𝑙 [ℛ𝒪]

Avec [𝐾0 ]𝑙 est la matrice de rigidité de l’élément portique dans le repère local.

Etabli par : Dr. KHEBIZI Mourad

12
Cours MEF Chapitre 2 : Méthode des éléments finis 1D

REFERENCES DU CHAPITRE 2

[8].Guenfoud, M (2002-203). Cours de la méthode des éléments finis, destiné aux étudiants
de 5ème année ingénieur promotion 2002-2003, département de génie civil, Université
de Guelma.

Etabli par : Dr. KHEBIZI Mourad

13
Cours MEF Chapitre 3 : Méthode des éléments finis 2D

CHAPITRE 3
METHODE DES ELEMENTS FINIS 2D

3.1. INTRODUCTION :
La méthode des éléments finis 2D est utilisée pour discrétiser les structures (matériaux,
sols,…) obéissent à l’hypothèse de déformation planes ou celle des contraintes planes :
- les sols (déformations planes,….)
- les refends (contraintes planes,….)
- les membranaires (contraintes planes,….)

3.2. Leurs types :

a) Triangulaire b) Quadratique
Figure 3.1 : Eléments finis 2D

3.2.1. Elément rectangle bilinéaire :

On considère, en élasticité plan, un élément fini rectangulaire de dimensions 2𝑎𝑥2𝑏, place
dans le système d’axes (𝑥, 𝑦) indiqué (figure 4.1.a). La théorie enseigne que le champ de
déplacement comporte les deux composantes 𝑢(𝑥, 𝑦) et 𝑣(𝑥, 𝑦) tant dans Ω que sur Γ [1]. Il
faut donc choisir comme degrés de liberté les valeurs de 𝑢 et 𝑣 en certains nœuds.
Si on sélectionne les quatre sommets du rectangle comme nœuds (figure 4.1.b), l’élément
comporte huit inconnues cinématiques [1]
{𝑞}𝑇 = {𝑢1 𝑣1 𝑢2 𝑣2 𝑢3 𝑣3 𝑢4 𝑣4 }

a) Géométrie et axes b) Nœuds et degrés de liberté

Figure 3.1 : Eléments fini rectangulaire d’élasticité plane

Etabli par : Dr. KHEBIZI Mourad

14
Cours MEF Chapitre 3 : Méthode des éléments finis 2D

3.2.1.1. Matrice de rigidité

La matrice de rigidité de l’élément rectangulaire, pour l’état plan de contrainte, s’écrit [1]

𝑐1 𝑐7 ⋮ 𝑐3 −𝑐8 ⋮ 𝑐9 −𝑐7 ⋮ 𝑐5 𝑐8
𝑐2 ⋮ 𝑐8 𝑐4 ⋮ −𝑐7 𝑐10 ⋮ −𝑐8 𝑐6
… … . . ⋮. . … … . . ⋮. . … … . . ⋮. . … …
⋮ 𝑐1 −𝑐7 ⋮ 𝑐5 −𝑐8 ⋮ 𝑐9 𝑐7
𝐸𝑡 ⋮ 𝑐2 ⋮ 𝑐8 𝑐6 ⋮ 𝑐7 𝑐10
[𝐾] = … … . . ⋮. . … … . . ⋮. . … … . . ⋮. . … …
1 − 𝜈2 ⋮ ⋮ 𝑐1 𝑐7 ⋮ 𝑐3 −𝑐8
⋮ ⋮ 𝑐2 ⋮ 𝑐8 𝑐4
… … . . ⋮. . … … . . ⋮. . … … . . ⋮. . … …
⋮ ⋮ ⋮ 𝑐1 −𝑐7
[𝑠𝑦𝑚 ⋮ ⋮ ⋮ 𝑐2 ]

Avec :
𝑡 est l’épaisseur de l’élément.
𝑎 1−𝜈 𝑏
𝛼= 𝜇= 𝛽=
𝑏 2 𝑎

𝜇𝛼 + 𝛽 𝑐1 − 𝛽 𝑐1 − 𝜇𝛼 𝜇+𝜈 𝑐1
𝑐1 = 𝑐3 = 𝑐5 = 𝑐7 = 𝑐9 = −
3 2 2 4 2

𝛼 + 𝜇𝛽 𝑐2 − 𝜇𝛽 𝑐2 − 𝛼 𝜇−𝜈 𝑐1
𝑐2 = 𝑐4 = 𝑐6 = 𝑐8 = 𝑐10 = −
3 2 2 4 2

3.2.1.2. Vecteurs forces

 Force de volume
Soit 𝑏𝑥 et 𝑏𝑦 les composantes constantes de la force de volume (figure 3.2), le vecteur force
dans ce cas s’écrit [1]

𝑏𝑥
𝑏𝑦
……
𝑏𝑥
𝑏𝑦
{𝑓} = 𝑡𝑎𝑏 … . . .
𝑏𝑥
𝑏𝑦
……
𝑏𝑥
{ 𝑏𝑦 }

Etabli par : Dr. KHEBIZI Mourad

15
Cours MEF Chapitre 3 : Méthode des éléments finis 2D

C’est-à-dire une répartition égale à chaque nœud.

Figure 3.2 : Force de volume

 Charge concentrée
Pour une charge concentrée de composantes 𝐹𝑥 et 𝐹𝑦 , agissant au point 𝐻de coordonnées
(figure 3.3) [1]

𝑥𝐻 = 𝜉𝑎 (−1 ≤ 𝜉 ≤ 1) 𝑦𝐻 = 𝜂𝑏 (−1 ≤ 𝜂 ≤ 1)

Le vecteur force (pour le cas d’une charge concentrée) est donne par [1]

(1 − 𝜉)(1 − 𝜂)𝐹𝑥
(1 − 𝜉)(1 − 𝜂)𝐹𝑦
……………………
(1 + 𝜉)(1 − 𝜂)𝐹𝑥
(1 + 𝜉)(1 − 𝜂)𝐹𝑦
1
{𝑓} = …...………………
4 (1 + 𝜉)(1 + 𝜂)𝐹
𝑥
(1 + 𝜉)(1 + 𝜂)𝐹𝑦
……………………
(1 − 𝜉)(1 + 𝜂)𝐹𝑥
{ (1 − 𝜉)(1 + 𝜂)𝐹𝑦 }

Si le point 𝐻 coïncide avec un nœud de l’élément, on vérifie aisément que seul ce nœud est
chargé : on aboutit à une charge purement nodale [1].

Etabli par : Dr. KHEBIZI Mourad

16
Cours MEF Chapitre 3 : Méthode des éléments finis 2D

Figure 3.3 : Charge concentrée au point H

3.2.1.3. Transformation de coordonnées

On désigne par les indices 𝑙 et 𝑔 les grandeurs attachées aux systèmes d’axes local d’étude
d’un élément et global d’étude de la structure. De manière générale, on peut toujours écrire une
relation du genre [1]
{𝑞}𝑙 = [𝑅𝑂]{𝑞}𝑔
Où [𝑅𝑂] est une matrice dite de transformation de coordonnées ou, plus simplement, de
rotation, puisqu’elle opère la rotation de degrés de {𝑞}𝑔 d’un élément fini, exprimés en axes
globaux, à ceux {𝑞}𝑙 définis en axes locaux.

La matrice de rigidité exprimée dans les axes globaux s’écrit donc

[𝐾 ]𝑔 = [𝑅𝑂]𝑇 [𝐾 ]𝑙 [𝑅𝑂]
Le vecteur force défini en axes globaux s’exprime par
{𝑓}𝑙 = [𝑅𝑂]𝑇 {𝑓}𝑔
Considérons un élément fini rectangulaire d’élasticité plane (figure 3.4). La connaissance
des coordonnées des nœuds permet de calculer l’angle 𝛼 orientant les axes locaux (𝑥, 𝑦) par
rapport aux axes globaux (𝑋, 𝑌). En chaque nœud, les déplacements sont 𝑢𝑖 et 𝑣𝑖 (𝑖 = 1,2,3,4).
On a donc simplement
𝑢𝑖 𝑐𝑜𝑠𝛼 𝑠𝑖𝑛𝛼 𝑢𝑖 𝑢𝑖
{𝑣 } = [ ] {𝑣 } = [𝑅𝑂]𝛼 {𝑣 }
𝑖 𝑙 −𝑠𝑖𝑛𝛼 𝑐𝑜𝑠𝛼 𝑖 𝑔 𝑖 𝑔

et, pour tous les nœuds, la relation {𝑞}𝑙 = [𝑅𝑂]{𝑞}𝑔 s’écrit

Etabli par : Dr. KHEBIZI Mourad

17
Cours MEF Chapitre 3 : Méthode des éléments finis 2D

𝑢1 𝑢1
𝑣1 𝑣1
…… [𝑅𝑂]𝛼 ⋮ 0 ⋮ 0 ⋮ 0 … …
𝑢2 … . . ⋮. . … . . ⋮. . … . . ⋮. . … 𝑢2
𝑣2 0 ⋮ [𝑅𝑂]𝛼 ⋮ 0 ⋮ 0 𝑣2
…… = … . . ⋮. . … . . ⋮. . … . . ⋮. . … … … = [𝑅𝑂]{𝑞}𝑔
𝑢3 0 ⋮ 0 ⋮ [𝑅𝑂]𝛼 ⋮ 0 𝑢3
𝑣3 … . . ⋮. . … . . ⋮. . … . . ⋮. . … 𝑣3
…… [ 0 ⋮ 0 ⋮ 0 ⋮ [𝑅𝑂]𝛼 ] … …
𝑢4 𝑢4
{ 𝑣4 }𝑙 { 𝑣4 }𝑔

Figure 3.4 : Rectangle à champ de déplacement bilinéaire en axes locaux et globaux

REFERENCES DU CHAPITRE 3

[1]. Frey, F., Jirousek, J. (1999), Analyse des structures et milieux continus – Méthode des
éléments finis-, Presses polytechnique et universitaire Romandes, Lausanne.

[8]. Guenfoud, M (2002-203). Cours de la méthode des éléments finis, destiné aux étudiants
de 5ème année ingénieur promotion 2002-2003, département de génie civil, Université
de Guelma.

Etabli par : Dr. KHEBIZI Mourad

18
Cours MEF Références)

REFERENCES

[2]. Frey, F., Jirousek, J. (1999), Analyse des structures et milieux continus – Méthode des
éléments finis-, Presses polytechnique et universitaire Romandes, Lausanne.

[3]. Dhatt, G. et Touzot, G. (1984), une présentation de la méthode des éléments finis
(deuxième édition), Maloine S.A, Paris

[4]. KHEBIZI, M. et GUENFOUD M., (2008), Modélisation d'un glissement de terrain lent
par la méthode des éléments finis -Cas du glissement de CILOC à CONSTANTINE –
(2008), séminaire national de génie civil 15-16 décembre 2008, Université de Chlef,
Algérie.

[5]. KHEBIZI, M., GUENFOUD, H. and GUENFOUD, M., (2018), Numerical modelling
of soil-foundation interaction by a new non–linear macro–element, Geomechanics and
Engineering, An International Journal “Tecno-press” , Vol.14, No.4 (2018), 377-386.

[6]. Helwany, S. (2007), Applied soil mechanics with Abaqus applications, John Wiley &
Sons, INC.

[7]. Etemadi, E., Fallahnezhad, E. (2017), Behavior of reinforced concrete interior wide
beam-column connections under lateral loading: A finite element study, International
Journal of Engineering and Technology (IJET), Vol 9 No 3, PP. 2559-2570

[8]. Frey, F. (1980), Analyse des structures et milieux continus – Mécanique des solides-,
Presses polytechnique et universitaire Romandes, Lausanne.

[9]. Guenfoud, M (2002-203). Cours de la méthode des éléments finis, destiné aux étudiants
de 5ème année ingénieur promotion 2002-2003, département de génie civil, Université
de Guelma.

Etabli par : Dr. KHEBIZI Mourad

19