ROYAUME DU MAROC

UNIVERSITE MOHAMMED PREMIER

ECOLE NATIONALE DES SCIENCES
APPLIQUEES
OUJDA

DEPARTEMENT MECANIQUE ET MATHEMATIQUES

APPLIQUEES
Filière : GENIE CIVIL
Niveau : 2éme année

Eléments de cours :
MODELISATION NUMERIQUE
CALCUL DE STRUCTURES EN ELEMENTS
FINIS (Chapitre 2)

Prof- BOUSHABA FARID

1
 Plan du cours
Chapitre 1 : Equation des contraintes
1.1. Introduction
1.2. Tenseur de déformation
1.3. Equation d’équilibre
1.4. Loi de comportement
1.5. Conditions aux limites
1.6. Principe des travaux virtuels
1.7. Théorème d’énergie potentiel totale (EPT)

Chapitre 2 : Méthode d’approximation en Eléments finis 1D

2.1. Introduction
2.2. Méthode de RITZ GALERKIN
2.3. Exemple Problème de traction
2.4. Elément poutre
2.5. Principe méthode des éléments finis (MEF)
2.6. Approximation cubique : Elément Poutre
2.7. Calcul de la matrice de rigidité
2.8. Assemblages élément poutre / barre

Chapitre 3 : Elément finis 2D domaine élastique

3.1. état de déformation plane
3.2. Etat de contrainte planes
3.3 Détermination de la matrice de rigidité et vecteur force 2
3.4 Fonction de forme pour élément triangulaire élémentaire
Equation des contraintes Méthode d’approximation en Eléments finis 1D Elément finis 2D domaine élastique

introduction

La méthode des éléments finis consiste à transformer un système aux

dérivés partiels continus en un système discret.

Modèle
Modèle Modèle numérique
d'EDP de Résultat
physique
résolution

3
Equation des contraintes Méthode d’approximation en Elémentsfinis1D Elément finis 2D domaine élastique

Méthode de RITZ GALERKIN

La méthode des éléments finis est basée sur l’approche de RITZ GALERKIN,
le problème se ramène à la recherche des fonctions qui minimisent l’énergie potentielle
totale, La méthode de Ritz-Galerkin consiste à rechercher une solution en
déplacement en construisant une combinaison linéaire avec des fonctions
d’approximation sous la forme :

Trouver la solution approchée du problème consiste à déterminer les scalaires qui

minimisent dans les sous espaces Vn l’énergie potentielle totale

4
:
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Exemple Problème de traction :

Traction ou compression :

On suppose une fonction test d’approximation dans le champ de déplacement :

5
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Exemple Problème de traction :

On remplace le champ de déplacement dans l’équation de déformation et de

travail :

Donc :

On a :

Avec
6
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Exemple Problème de traction :

On remplace l’indice i par m :

On a par définition :

Donc :

La matrice de rigidité est symétrique et les indices i et j sont mués par conséquent :

7
D’ou :
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Elément poutre

Soit un élément de poutre de section A, de longueur L, de module de Young E

et de moment d’inertie I.

Les champs les forces extérieures

sont :

Pour réaliser une sollicitation de flexion

simple on pose :

On impose aussi :

Les moments sur les deux extrémités sont :

8
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Elément poutre

Rappel : Flexion (relation géométrique de la configuration déformée):

On considère la fibre de la poutre après déformation :

Avec I moment quadratique

9
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Elément poutre

Equation d’équilibre locale :

On considère un tronçon élémentaire de poutre de longueur dx

D’après l’équilibre des forces : (*)

D’après l’équilibre des moments :

On remplace dans (*)

Finalement, on obtient l’équation d’équilibre locale :

10
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Principe de la méthode des éléments finis (MEF) :

Dans le cas des chargements simples et des problèmes 1D, la méthode de Ritz
est bien pratique et précise, par contre pour des chargement complexes et des
problèmes de caractères 2 D et 3 D ou les fonctions testes cinématiques admissible
ne se réduit pas à un point au niveau des frontières du domaine, la méthode de Ritz
se trouve heurté et présente des handicaps majeur. La solution à ce type de problème
est les éléments finis.

La méthode des éléments finis est une méthode d’approximation numérique, Elle est
basée sur la méthode de Ritz-Galarkin avec une construction des fonctions testes
continues et dérivablement continues par morceaux.

11
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Approximation linéaire, élément à deux nœuds :

Soit une fonction U(x) définie sur un domaine [0 ; L], on considère deux
nœuds de coordonnées et , et de valeurs nodales et

Le pas du maillage :

12
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Approximation linéaire, élément à deux nœuds :

N est la fonction de forme, on représente les fonctions de formes ou d’interpolation

dans le graphe suivant :

Les fonctions de formes Ni sont des fonctions sous forme de chapeau :

13
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Exemple : Treillis

Un treillis est un assemblage soudé ou articulé d’éléments barres, on suppose que

le chargement s’effectue au niveau des nœuds et le poids propre des barres est
négligeable, le système des structures travaille uniquement en traction-compression

14
,
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Exemple : Treillis

On suppose un élément de barre de nœuds 1-2 lié à un repère local et le système

est relié au repère global ( , ).

15
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Exemple : Treillis

On écrit le système ci- dessous sous forme matricielle :

= =

est la matrice de passage du repère globale au repère local, on note le vecteur

global de composantes :
=

Le vecteur local exprimé dans le repère local porté par les axes du vecteur :

=
16
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Exemple : Treillis

On exprime l’énergie de déformation dans le repère global, on a :

On remplace le vecteur déplacement local en fonction du vecteur déplacement

global, on trouve :

On pose : =

Le calcul des composantes tensorielles de la matrice de rigidité dans le repère global :

=
17
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Exemple : Treillis

Finalement :

Donc la matrice de rigidité globale est :

18
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

=
Approximation cubique : Elément Poutre

Pour déterminer la fonction de forme d’un élément

= poutre qui travaille en flexion
simple, on utilise une approximation cubique de la fonction déplacement V(x)

Pour déterminer les paramètres généraux , on définit les conditions aux limites
statiquement admissible pour un élément poutre ( L, E, A, I )

V1  1 0 0 0 
  0 1 0 0 
 1 
V2  = 1 L L2 L3 
    19
 2  0 1 2 L 3L2 
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Approximation cubique : Elément Poutre

= =

Tout calcul fait :

N1 et N3 représentent la déformée de la poutre encastré pour laquelle on impose

un déplacement unité dans l’autre extrémité.

N2 et N4 représentent la déformée de la poutre encastré sur l’extrémité pour la quelle

on impose une rotation unité à l’autre extrémité.

20
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Calcul de la matrice de rigidité

L’énergie de déformation pour une poutre en flexion simple :

Nivi ; =Ni”vi : D’ou <vi>

Donc : = E I EI

EI

Tout calcul fait :

21
 F
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Cas de chargement réparti : q(x)

On suppose une poutre chargée avec une charge réparti

Les forces nodales équivalentes sur les deux extrémités :

Moments équivalents sur les deux extrémités :

Finalement le vecteur nodal des efforts est :

22
,
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Assemblages élément poutre / barre

Soit deux éléments, un élément barre de caractéristiques ( ) et un élément poutre de
caractéristiques ( )

Le vecteur DDL et le vecteur effort sont :

Pour construire la matrice de rigidité résultante des 2 sollicitations on a :

Elément barre : Elément poutre :

23
Equation des contraintes Méthode d’approximation en Eléments finis1D Elément finis 2D domaine élastique

Assemblages élément poutre / barre

La matrice de rigidité assemblée de l’élément poutre et barre est :

Le vecteur effort assemblé :

24