N° d’ordre : /2013/DM

RÉPUBLIQUE ALGÉRIENNE DÉMOCRATIQUE ET POPULAIRE

MINISTÈRE DE L’ENSEIGNEMENT SUPÉRIEUR ET DE LA RECHERCHE SCIENTIFIQUE

UNIVERSITÉ DE BATNA
FACULTÉ DE TECHNOLOGIE
DÉPARTEMENT DE MÉCANIQUE
Mémoire présenté Pour l’obtention du diplôme de

MAGISTÈRE
En
Mécanique
Option : CONSTRUCTION MÉCANIQUE

Par

BOUSBIA SALAH Seif Eddine

Analyse vibratoire forcée des structures

par éléments finis

Soutenu le : / / 2013 Devant le jury composé de

Pr OUTTAS TOUFIK Prof. Université de Batna Président

Pr BRIOUA MOURAD Prof. Université de Batna Rapporteur

Pr ZIDANI KAMEL Prof. Université de Batna Examinateur

Pr ZEDIRA HAMMA Prof. Université de Khenchela Examinateur

Année universitaire 2012 / 2013

Dédicaces

A mes chers parents

 Remerciements

Tous mes remerciements vont au Dieu en premier lieu.

Je remercie infiniment mon encadreur Pr. [Link], pour avoir accepté

encadrement dans ce sujet.

Je remercie vivement tout particulièrement Dr. [Link] pour son

énorme aide, sa sympathie, ses conseils et sa disponibilité et surtout sa
patience pendant toute la durée de la réalisation du mémoire.

Je tiens à remercier : les professeurs OUTTAS TOUFIK, ZIDANI

KAMEL et ZEDIRA HAMMA pour nous avoir fait l'honneur d’être
membres du jury. Ainsi que pour avoir consacré une partie de leur temps
précieux pour lire et corriger ce mémoire. On portera beaucoup d’intérêt à
leurs remarques constructives.

Mes sincères remerciements et mes vives reconnaissances tous ceux qui

m’ont aidé de près ou de loin à réaliser ce travail.
 Table des matières

Table des matières

Remerciements
Table des matières
Liste des figures et des tableaux
Nomenclature et Indices
Introduction Générale …………………...…………………………………………….…….1

Chapitre I.
Etudes bibliographiques
I.1. Introduction ……………………………………………..……………………………….3
I.2. Revue bibliographique ………………..…………………………………………………3

Chapitre II :
Formulation de quelques éléments basés sur le champ de déplacement
II.1. Introduction……………………………..……………………………………………….8
II.1.1 Qu'est que c'est la méthode des éléments finis (MEF)………….………………...8
II.1.2 Principes de la MEF…………………………………….………………………...9
II.1.3. Les différents types de modèles d’éléments finis……………….………..……..10
II.2. Le concept d'analyse de rigidité………………….…………………………………....12
II.3. Procédure générale pour la formulation de la matrice de rigidité d’un élément
fini………………………………………………………………………...…………………..14
II.3.1. Etape 1: Identifié le problème ………….…………………………….…………14
II.3.2. étape 2 : Choisir des fonctions de déplacement convenables…………......…….14
II.3.3. Etape 3: Relié les déplacements généraux de l'élément à ses déplacements
nodaux……………………………………………………………...…………………………15
II.3.4. étape 4 : Relations entre déplacements et déformations cas général...………….16
II.3.5. étape 5 : relations entre déformations et contraintes……….……………………17
II.3.6. étape 6 : Relié les charges nodales aux déplacements nodaux ……………..…..17
II.4. Formulation de la matrice de rigidité d’élément fini de poutre..…………...........…19
II.5. Eléments finis rectangulaires pour l'étude de l'élasticité plane……………………..24
II.5.1. l’Elément bilinéaire quadratique (élément rectangulaire à quatre nœuds) (Q4)..24
II.5.2. Elément quadrilatéral quadratique (élément rectangulaire à 8 nœuds) (Q8)…....31
 Table des matières

II.5.3. Présentation de l’Élément T6, Élément Triangulaire à 6 nœuds (ou LST linear
strain triangle)..…………………………………………………………………………….…33
II.6. Assemblage des matrices de rigidité des éléments……………………………….35
II.7. Détermination de matrice de masse des éléments …………………..….………....36
II.7.1. La matrice masse d’un élément de poutre………………………………..…..……….37
II.7.2. La matrice masse d’un élément rectangulaire………………………….…..…………38

Chapitre III :
Analyse dynamique des structures par éléments finis

III.1. Introduction……………………………………………………………..…………….39
III.2. Dynamique d’une poutre par la méthode analytique………………………...…….39
III.2.1. Déformée d’une poutre en flexion dans l’hypothèse de la résistance des
matériaux…………………………………………………………………………………..….39
[Link] différentielle de la déformée en terme de chargement linéaire q(x).....41
[Link] différentielle des vibrations d’une poutre…………..……………..…42
III.2.3.1. Vibrations libres…………………………………………………...…42
III.[Link]. Modes propres de vibration d’une poutre cantilever de section
rectangulaire………...………………………………………………………..…44
III.2.3.2. Vibrations forcée…………………………….…………………...….46
III.3. Etude vibratoire par éléments finis……….…………………………………………47
III.3.1. Élément unidimensionnel (1D) : élément de poutre……………...….……...…47
III.3.1.1. poutre avec des forces verticales et une charge répartie uniforme…..47
III.3.1.2. poutre avec des forces verticales et une charge appliquée au niveau des
nœuds……………………………………………………………………………...… 49
III.3.1.3. Matrice de rigidité d’une poutre ………………………………..……49
III.3.1.4. Poutre avec une force axiale………………………………….………50
III.3.2. Eléments bidimensionnels (2D) : contrainte et déformation plane….…………51

Chapitre IV :
Validation numérique des éléments à champs de déplacement en dynamique
IV.1. Introduction……………………………………………………………...……………55
IV.2. Cas de tests………………………………………………………………………...…..56
 Table des matières

IV.2.1. Valeurs propres en vibration libre d’un solide rectangulaire …………………57

IV.2.1.1. Synoptique du programme de vibration libre d’un solide
rectangulaire……………………………………………………………………………57
IV.2.1.2. Tests et validations……………………….…………………………..….58
IV.2.2. Vibration forcée d’un solide rectangulaire……………………………………..61
IV.2.2.1. Synoptique du programme de vibration forcée d’un solide
rectangulaire……………………………………………………………………………61
IV.2.2.2. Tests et validations………………………………………….………...…62
IV.3. Conclusion………………………………………………………………………….….68

Conclusions générale………………………………….………...…………………………..69

Références bibliographiques.
 Liste des figures et des tableaux

Liste des figures

Fig.I.1 : Trapézoïdale de Fourier de p-élément
Fig.I.2 : Maillages trapézoïdale dans un bar élastique
Fig.I.3 : Les différents types de Maillages d’un cantilever plan
Fig.I.4 : Élément quadrilatère courbé
Fig.I.5 : Élément coque quadrilatère à 8 nœuds de contrainte assumée hybride
Fig.I.6 : Maillage des éléments triangulaires dans la méthode des éléments finis alpha
alternatif (AαFEM)

[Link].1 : Structure treillis (exemple de Pont)..

[Link].2 : Modélisations par éléments finis de structures continues, a Mur de barrage, b plaque
pliée.
[Link].3 : Principes d’analyse par la MEF
[Link].4 : Simple ressort élastique.
[Link].5 : Relation entre force et déplacement pour un simple ressort élastique
[Link].6 : Armature de poutres hyperstatique.
[Link].7.a : Déplacements transversaux et de rotation à chaque extrémité de l'élément de
poutre EP.
[Link].8.b : Position déformé et les forces et des moments en chaque nœud de poutre EP.
[Link].8.a: Elément quadrilatère bilinéaire Q4, a et b : Déplacements nodaux.
[Link].8.b : Elément quadrilatère bilinéaire Q4, c : Forces nodales.
[Link].9 : Déformations dans l’élément fini rectangulaire
[Link].10 : Elément quadratique quadrilatéral Q8
[Link].11: Elément de type T6 et ses 12 degrés de liberté
[Link].12 : deux éléments de poutre 1 et 2 avec L1 et L2

[Link].1 : Géométrie d’une poutre

[Link].2 : Le moment de flexion
[Link].3 : L’équilibre mécanique du tronçon de longueur dx soumis au chargement q(x).
[Link].4 : Vibration libre d’une poutre cantilever
[Link].5 : courbe d’équation de caractéristique
 Liste des figures et des tableaux

[Link].6 : Vibration forcée d’une poutre cantilever

[Link].7 : Mince élément de poutre
[Link].8 : Poutre avec une force axiale
[Link].9 : (a) Mur de cisaillement avec des ouvertures. (b) L'élément typique rectangulaire à
quatre nœuds d'angle

[Link].1: Géométrie et maillage de la poutre.

[Link].2 : Synoptique du programme de vibration libre d’un solide rectangulaire
[Link].3 : Maillage de (3×1) et l’élément Q4
[Link].4: Maillage de (5×1) et l’élément Q4
[Link].5: Maillage de (3×1) et l’élément Q8
[Link].6: Maillage de (5×1) et l’élément Q8
[Link].7: Maillage de (3×1) et l’élément T6
[Link].8: Maillage de (5×1) et l’élément T6
[Link].9 : Synoptique du programme de vibration forcée d’un solide rectangulaire (masse
ponctuelle)
[Link].10 : Géométrie et maillage de la poutre cantilever soumise à une vibration forcée
[Link].11 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au temps
pour l’élément Q8" Test 1
[Link].12 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au temps
pour l’élément Q8" Test 1
[Link].13 : La convergence entre l’élément Q8 et T6
[Link].14 : Géométrie et maillage de solide rectangulaire soumis à une vibration forcée
[Link].15 : Sollicitation harmonique "Les forces par rapport au temps"
[Link].16 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au temps
pour l’élément Q8" Test 2
[Link].17 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au temps
pour l’élément Q8" Test 2

Liste des tableaux

[Link].1: Valeurs propres de la poutre cantilever cas de déformation plane.
[Link].2: Vibration forcée d’un solide rectangulaire "déplacement par rapport au temps"
[Link].3 : Vibration forcée d’un solide rectangulaire avec masse localisée " déplacement
par rapport au temps et forces".
 Nomenclature et Indices

Nomenclature et Indices

A, S : Aire

[A] : Matrice associée à la fonction de déplacement

[B] : Matrice reliant les déformations de l'élément à ses déplacements nodaux

[D] : Matrice d'élasticité

Dx, Dy, Dxy, D1 : Rigidités de flexion d'une plaque

E : Module d'Young

F : Force

{F} : Vecteur des forces nodales

Fx, Fy, Fz: Forces dans les directions x, y et z

G : Module de Coulomb

[H] : Matrice reliant contraintes et déplacements

I : Deuxième moment de l’aire

k : Rigidité

kij : Terme de [K] situé sur la ligne i de la colonne f

[K]Matrice de rigidité

L : Longueur

L1, L2, L3 : Coordonnées barycentriques

Mx, My, Mxy : Moments de flexion interne

Nx, Ny : Fonctions de pondération (chapitre 1 )

t : Épaisseur (t pour « thickness»)

[T] : Matrice de changement de coordonnées ou de base

 Nomenclature et Indices

Tx, Ty, Txz : Moments le long des axes x, y et z

u, v, w : Déplacements le long des axes x, y et z

x, y, z : Système de coordonnées cartésiennes rectangulaires

α1, α2 : Constantes utilisées dans la fonction de déplacement

γxy, γyz, γzx : Déformations de cisaillement

δ : Déplacement

{δ} : Vecteur des déplacements nodaux

Δ : Aire de l'élément

ε : Déformation

εx, εy, εz : Déformations directes

θx, θy, θz : Rotations autour des axes x, y et z

θx, θy, θz : Angles entre les axes des systèmes de coordonnées local et global

ζ : Contrainte

ζx, ζy, ζz : Contraintes directes

ηx, ηy, ηz : Contraintes de cisaillement

ν : Coefficient de Poisson

[ ] : Désigne une matrice

{ } : Désigne un tableau unidimensionnel, vecteur ligne ou colonne

[ e ] { e } : Matrice, vecteur se rapportant à un seul élément

[ ]T : Transposée d'une matrice

{ i}, xi : Désigne des quantités associées au nœud i

(x, y) : Désigne des quantités fonctions de x et de y

[ ]-1 : Inverse u une matrice

 Nomenclature et Indices

Abréviation Signification

MEF Méthode d’éléments finis.

AVMEF Analyse vibratoire par méthode des éléments finis

Q4 Elément rectangulaire à 4 nœuds.

Q8 Elément rectangulaire à 8 nœuds.

E-P Élément poutre.

T6 Élément Triangulaire à 6 nœuds

DDL Degré De Liberté

1D unidimensionnel

2D bidimensionnel
 Introduction générale

Introduction générale
Pour évaluer les caractéristiques vibratoires dans la conception des produits, des pièces
mécaniques et des structures, l’ingénieur a besoin de modèles qui lui permettent de simuler le
comportement dynamique de systèmes physiques complexes. Il peut ainsi prévoir l'influence
de ses décisions au moment de la conception du système.
Les sciences de l'ingénieur (mécanique des solides et des fluides, thermique...)
permettent de décrire le comportement de systèmes physiques grâce à des équations aux
dérivées partielles. La méthode des éléments finis est l'une des méthodes les plus utilisées
aujourd'hui pour résoudre effectivement ces équations. Elle nécessite l'utilisation intensive de
l'ordinateur. C'est une méthode très générale qui s'applique à la majorité des problèmes
rencontrés dans la pratique : problèmes stationnaires ou non stationnaires, linéaires ou non
linéaires, définis dans un domaine géométrique quelconque à une, deux ou trois dimensions.

La méthode des éléments finis consiste à utiliser une approximation simple des
variables inconnues pour transformer les équations aux dérivées partielles en équations
algébriques. Elle fait appel aux trois domaines suivants :

 Sciences de l’ingénieur pour construire les équations aux dérivées partielles.

 Méthodes numériques pour construire et résoudre les équations algébriques.
 Programmation et informatique pour exécuter efficacement les calculs sur ordinateur.

La méthode des éléments finis subit un développement rapide dans plusieurs

directions depuis son apparition en1956, plus exactement après les publications de Turner,
Clough, Martin et Topp [1] qui ont introduit le concept d’élément fini et le travail de Argyris,
et Kesly [2] qui systématisa l’utilisation de la notion d’énergie dans l’analyse des
structures.[22]

En effet l’utilisation des éléments finis a champs de déplacement pour l’analyse

statique, non linéaires et vibratoire connait une certaine extension fantastique qui va de pair
avec l’évolution de l’industrie .Dans le contexte d’avoir de bons résultats dans l’analyse
structurale, l’ingénieur à besoin de plus en plus aux outils numériques de simulation par
éléments finis à la place de outils de prototypage ou expérimentale qui sont très couteux.

Dans ce travail de mémoire une investigation dans le comportement dynamique des

éléments à champs de déplacement Q4 (élément rectangulaire à Quatre nœuds), Q8 (élément

1
 Introduction générale

rectangulaire à huit nœuds),EP (Élément Poutre)et l’élément à champs de déplacement T6

(élément Triangulaire à 6 nœuds) dans l’analyse vibratoire libre et forcée est décrite.

Avec l’outil informatique, on utilise la programmation FORTRAN pour tester les

performances de ces éléments dans les problèmes vibratoires des structures.

Les solutions d'éléments finis obtenus pour ces problèmes sont comparées aux valeurs
analytiques proposées en dynamique des structures.

Les résultats obtenus à l'aide de l’élément T6 sont aussi comparés à celles des
éléments Q4, Q8, et Élément Poutre respectivement.

 L’objectif

La prédiction du comportement de l’élément T6 en analyse vibratoire libre et forcée

des solides rectangulaires par la méthode des éléments finis est l’objectif essentiel de ce
travail.

 Plan de mémoire

Le mémoire est divisé en quatre chapitres :

Le premier est consacré à une analyse bibliographique dont laquelle on donne un

contexte historique sur les travaux antérieures en dynamique.

Le deuxième chapitre est purement consacré à la formulation de quelques éléments

basés sur de champs de déplacements.

Dans le troisième chapitre est consacré au Analyse dynamique des structures par la
méthode d’éléments finis et un petit rappel sur les méthodes analytiques.

Le quatrième chapitre est réservé aux validations de l’élémentT6 en analyse vibratoire

libre et forcée. Pour cela différents tests sont illustrés dans cette étude.

2
Chapitre I : Etudes bibliographiques

CHAPITRE I :
ÉTUDES BIBLIOGRAPHIQUES

I.1. INTRODUCTION
Le domaine d’éléments finis est très vaste et il à connu depuis son apparition un
développement très rapide, chaque jour il y a de nouveaux éléments qui peut résoudre ce
problème ou l’autre.
Dans ce chapitre on fait une analyse bibliographique dont laquelle on donne un
contexte historique sur les travaux antérieurs sur l’utilisation de la méthode des éléments finis
dans l’analyse vibratoire des structures et revue sur les différents éléments utilisés ou
proposés.

I.2. REVUE BIBLIOGRAPHIQUE

En 2003 A.Y.T. Leung, Zhu Bin, Jian jin Zheng, Hao Yang[8], qui utilisent un
élément appelé trapézoïdale de Fourier de p-élément pour l'analyse des vibrations de solides
élastiques à deux dimensions 2D est présentée comme suit :

Fig.I.1 : Trapézoïdale de Fourier de p-élément

Les cas des tests utilisés est la vibration libre longitudinale d’une barre élastique et
d’une poutre cantilever plane et la comparaison est faite avec des éléments : linéaire
Triangulaire T3, linéaire rectangulaire LR,Q8 et Q6 et avec les formules d'intégration
analytique.

3
Chapitre I : Etudes bibliographiques

Fig.I.2 : Maillages trapézoïdale dans un bar élastique

Fig.I.3 : Les différents types de Maillages d’un cantilever plan

Ils montrent que cet élément de Fourier est plus précis que ceux des éléments utilisant
la quadrature de Gauss pour les problèmes de vibration dans le plan, et cet élément est plus
efficace pour les vibrations de haute fréquence.

En 2006 les auteurs Yeon Seok Choo, Noori Choi Byung, Chai Lee [9], ont proposé
des éléments plan quadrilatères et triangulaires élastiques linéaires flexibles, précis et
robustes, avec l'introduire des degrés de liberté supplémentaires de rotation.

4
Chapitre I : Etudes bibliographiques

Ils montrent que les degrés de rotation supplémentaires augmentent la précision et la

qualité de robustesse de l’élément développé.

Les travaux de recherche de A. Houmat [10] 2006, ont été consacrés à l’étude de P-
Version trigonométrique de la méthode des éléments finis basé sur un élément quadrilatère
courbé figure (I.4).Cet élément a été développé et appliqué à l'analyse des vibrations libres
des structures en membranes de forme arbitraire.

Fig.I.4 : Élément quadrilatère courbe

Des nouvelles fonctions de forme ont été exprimées en termes de fonctions

trigonométriques. Le calcul de ces fonctions à des ordres élevés est plus stable par rapport à
ceux de polynômes orthogonaux.
Il a calculé les valeurs de fréquence pour les éléments secteurs ouverts et fermés à
l'aide de maillages à un élément quadrilatère courbe et a montré que la convergence est rapide
et de haute précision par cette méthode. Ces valeurs sont très précises et peuvent servir de
base de comparaison pour les autres techniques de calcul.

Les recherches de Kutlu Darılmaz et de Nahit Kumbasar[11]en 2006, ont été

consacrés à l’étude d’un élément coque quadrilatère à 8 nœuds de contrainte assumée hybride,
avec six dégrée de liberté par nœud, et chaque nœud à trois déplacement et trois rotation :

5
Chapitre I : Etudes bibliographiques

Fig.I.5 : Élément coque quadrilatère à 8 nœuds de contrainte assumée hybride.

Ils ont étudié la performance de l'élément pour différentes cas de structures, mais ils
ont conditionné l’utilisation de l'élément dans l'analyse statique et vibratoire des structures
avec propriétés de matériau (isotrope / orthotrope).

En 2009 N. Nguyen-Thanha, et all [12], ont utilisé la méthode des éléments finis alpha
alternatif (AαFEM) pour l'analyse des vibrations libre et forcée des structures en utilisant les
maillages avec éléments triangulaires (figure.I.6)

Fig.I.6 : Maillage avec éléments triangulaires dans la méthode des éléments finis alpha
alternatif (AαFEM)

6
Chapitre I : Etudes bibliographiques

La méthode proposée améliore considérablement la précision des éléments finis

triangulaires standard et fournit une solution super convergente dans la norme de l'énergie
pour l'analyse statique et dynamique des solides dans les problèmes de deux dimensions.
Ainsi ils ont trouvé que :
 Les résultats numériques de la AαFEM utilisant des éléments triangulaires sont
toujours plus précis que celles de FEM-T3 et même plus précis que ceux de la FEM-
Q4 avec le même nombre de nœuds. Les taux de convergence de la norme de l'énergie
sont asymptotiquement les mêmes que celle de cette technique standard FEM.

 Dans la fréquence propre et l’analyses des vibrations forcées, la AαFEM est toujours
stable et donne des résultats plus précis que la méthode correspondante FEM-T3 et
FEM-Q4;

 La méthode AαFEM est facile à mettre en œuvre dans un programme d'éléments finis
et des maillages triangulaires sont idéals pour les domaines de problèmes complexes.

7
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

CHAPITRE II :
FORMULATION DE QUELQUES ÉLÉMENTS BASÉS SUR LE
CHAMP DE DÉPLACEMENT

Ce chapitre est une synthèse bibliographique des références [3], [4], [5].

II.1. INTRODUCTION

II.1.1 Qu'est que c'est la méthode des éléments finis (MEF)[3] et [5]

La méthode des éléments finis est une technique d'analyse numérique qui permet
d'obtenir des solutions approchées dans une large variété de problèmes d'ingénieries .Tout au
début cette méthode était développée pour étudier les structures complexes. [5]
Un grand nombre de structures utilisées par les ingénieurs sont composées d'une
réunion de parties indépendantes reliées entre elles par des points. Nous appellerons de telles
structures des «structures treillis», les points auxquels les parties indépendantes sont reliées
étant communément appelés nœuds. [3]

Fig. II.1 : Structure treillis (exemple de Pont).

L’analyse de telles structures peut s'effectuer en considérant d'abord le comportement

de chaque partie indépendante puis en assemblant ces parties de telle façon que l'équilibre des
forces et la compatibilité des déplacements soient satisfaits en chaque nœud.

Mais par la suite elle a était étendue et appliquée dans le champ des mécaniques
continues et les ingénieurs ont souvent besoin d'étudier des structures continues, telles que les
poutres épaisses, les plaques sujettes à la flexion, les plaques pliées et les coques.
8
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Dans la méthode des éléments finis, les parties continues sont représentées comme des
structures composées d'éléments indépendants reliés uniquement par les nœuds (fig. II.2).

Fig. II.2 : Modélisations par éléments finis de structures continues,

a Mur de barrage, b plaque pliée.

II.1.2 Principes de la MEF :

 La MEF est basée sur une idée simple : subdiviser (discrétiser) une forme complexe
en un grand nombre de sous domaines élémentaires de forme géométrique simple
(éléments finis) interconnectés en des points appelés nœuds.
 Nous considérons le comportement mécanique de chaque élément séparément, puis
nous assemblons ces éléments de telle façon que l’équilibre des forces et la
compatibilité des déplacements soient satisfaits en chaque nœud.
 La MEF utilise des approximations simples des variables inconnues dans chaque
élément pour transformer les équations aux dérivées partielles en équations
algébriques.
 Les nœuds et les éléments n’ont pas forcement de signification physique particulière,
mais sont basés sur des considérations de précision de l’approximation.

9
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Problème physique Complexe

(Milieu continu)
Sciences de l’ingénieur
pour établir les équations
Equation aux dérivées partielles

Approximation par éléments finis Programmation et

(Discrétisation) informatique pour
exécuter efficacement les
calculs sur l’ordinateur
Système d’équation algébrique
(Forme discrétisée)

Assemblage

Résolution numérique

Solutions Approchées
(numériques)

[Link].3 : Principes d’analyse par la MEF

II.1.3. Les différents types de modèles d’éléments finis: [6]

Les modèles d’élément finis les plus courants sont les suivants :

a) Modèle déplacement
Les éléments finis de cette catégorie sont basés sur une interpolation du champ des
déplacements étendue à tout l’élément, il en résulte que les déplacements sont déterminés de
manière détaillée et unique dans la structure, alors que les contraintes ne peuvent être connues
que par certains moyens et ne sont pas continués aux frontières. Le modèle déplacement est le
plus populaire, le plus connu et le plus développé. Dans ce modèle, les fonctions de
déplacements sont choisies en utilisant le triangle de pascal, de telle sorte que le nombre total

10
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

de constantes indépendantes soit égal au nombre total de degrés de libertés de l’élément, donc
l’invariance géométrique est maintenue.

b) Modèle contrainte
A l’inverse du précédent, un tel élément est construit à partir d’une approximation du
seul champ des contraintes. La connaissance des contraintes est alors détaillée, celle des
déplacements sommaire. Ce modèle est rarement employé, car il conduit en général, à qualité
égale, à un plus grand nombre d’inconnues que le modèle déplacement.

f) Modèle déformation
Pour obtenir une convergence satisfaisante, les éléments dérivant de l’approche
déplacement exigent un grand nombre d’éléments dans le maillage des structures. Par
conséquent, l’approche basée sur le modèle déformation a été développée, pas seulement pour
les éléments courbés mais aussi pour les structures plates. L’approche est basée sur le calcul
des termes exacts qui représentent tous les modes du corps rigide, les autres composantes des
fonctions de déplacements sont basées sur la forme des fonctions de déformations
indépendantes. Au fur et à mesure qu’on vérifie les équations de compatibilité, ce modèle
permet de proposer des champs de déformations en fonction de paramètres propres à la
structure(coefficient de poisson par exemple) tout en vérifiant les équations d’équilibre. Cette
approche mène, généralement, à la représentation des déplacements par un polynôme d’ordre
élevé sans avoir besoin de degrés de liberté additionnels (internes et inutiles). Une
convergence plus rapide est obtenue habituellement quand les résultats sont comparés avec les
éléments du déplacement correspondant, c’est-à-dire, éléments en déplacement qui ont le
même nombre total de degré de liberté.

c) Modèle mixte
Ce modèle est basé sur des approximations indépendantes de deux ou plusieurs
champs (déplacements, déformations, contraintes) étendues à tout l’élément ; en général, ce
modèle conserve les paramètres inconnus de tous les champs comme degrés de liberté ; la
mise au point peut être longue et coûteuse ; il peut toutefois être d’excellente qualité.
Le modèle mixte est parfois utilisé pour corriger certains défauts du modèle déplacement : on
interpole, par exemple, certaines composantes internes de la déformation pour améliorer la
précision au niveau des contraintes ; les paramètres correspondants peuvent alors être
éliminés par condensation, de sorte que l’élément reprend l’aspect d’un modèle déplacement.

11
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

II.2. LE CONCEPT D'ANALYSE DE RIGIDITÉ

Le plus simple composant de structure avec lequel la plupart des ingénieurs sont
familiarisés est le lien articulé à ses deux extrémités qui possède les mêmes caractéristiques
que le ressort élastique (fig. II.4). Pour un tel élément il existe une relation directe entre la
force F dans le ressort et le déplacement δ de son extrémité libre. Cette relation prend la forme
de l'équation II.2.1 :

F = k .δ (II.2.1)

Fig. II.4. - Simple ressort élastique.

La quantité k est appelée la rigidité (ou raideur) du ressort et correspond à la pente de la

courbe: force en fonction du déplacement (fig. II.5). Connaissant les valeurs de la rigidité et
de la force appliquée, l'équation II.2.1 donne le déplacement.

Fig. II.5. - Relation entre force et déplacement pour un simple ressort élastique.
12
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Mais tandis que la valeur d'un seul déplacement est suffisante pour caractériser
complètement l'état de déformation d'un simple ressort (fig. II.5), il est par contre nécessaire
de déterminer les déplacements des joints B, C, D et E lorsqu'on s'intéresse à des structures
plus compliquées comme la structure hyperstatique articulée de la figure II.6

Fig. II.6. - Armature de poutres hyperstatique.

Afin de pouvoir évaluer les contraintes existant dans les différentes parties. Dans toute la
suite, des points comme B, C, D et E seront appelés nœuds. Supposons que, pour une
structure complète, on puisse déterminer une quantité similaire à k dans l'équation II.2.1 ; il
est alors nécessaire d'écrire cette dernière sous la forme matricielle de l'équation II.2.3.

Les quantités {F} et {δ} sont respectivement les vecteurs représentant les charges et les
déplacements nodaux. On a alors :

{F} = [K] {δ} (II.2.3)

Où la quantité [K] est la rigidité de la structure complète.

Ce concept est la base de la méthode d'analyse des rigidités. La quantité [K] est
appelée la matrice de rigidité pour la structure et relie les forces nodales appliquées {F} aux
déplacements nodaux inconnus {δ}. Dans le cas de la structure simple de la figure II.4 pour
laquelle un seul déplacement est possible, cette matrice se réduit à un seul élément et les
vecteurs {F} et {δ} contiennent un terme chacun.

Et nous donnerons des détails sur la manière d'assembler cette matrice pour une
structure complète. Nous allons d'abord examiner la forme prise par la matrice de rigidité

13
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

dans le cas de simples parties articulées. En effet, la bonne compréhension de l'analyse

matricielle des structures comprenant de simples parties rectilignes donne la clé de toute la
; méthode des éléments finis.

II.3. PROCÉDURE GÉNÉRALE POUR LA FORMULATION DE LA

MATRICE DE RIGIDITÉ D’UN ÉLÉMENT FINI

On traite principalement de la première étape de l'analyse des rigidités d'une structure

ou d'un composant, c'est-à-dire la formulation de la matrice de rigidité [Ke] de l'élément.

On présente les sept étapes de base utilisées pour le calcul de la rigidité d'un élément
fini.

II.3.1. ETAPE 1: Identifié le problème :

L'étape I consiste à choisir un système de coordonnées et une numérotation des nœuds

appropriés pour l'élément.

Puisque l'on connait les degrés de liberté de l'élément fini choisi, on peut déterminer le
vecteur de déplacement nodal {δe} et le vecteur de charge nodale {Fe}. La matrice de rigidité
[Ke] pour cet élément est alors définie par l'équation II.3.1.

Fe = Ke  δe (II.3.1)

II.3.2. ETAPE 2 : Choisir des fonctions de déplacement convenables

On choisit une fonction de déplacement qui définit de façon unique l’état de

déplacement de tous les points de 'élément.

Ce modèle peut être représenté de façon commode par un expression polynomiale, et puisque
son but est d'exprimer les déplacements {δ(x,y)} de n'importe quel point en termes de
déplacements nodaux {δe} elle doit contenir un coefficient inconnu pour chaque degré de
liberté de l'élément. L'état des déplacements en chaque point (x,y) de l'élément peut être
décrit sous forme matricielle par l'équation II.3.2.

  x , y   f  x , y   (II.3.2)

14
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

où {α}est le vecteur colonne des coefficients encore inconnus de la fonction polynôme

[f(x,y)].

II.3.3. ETAPE3: Relier les déplacements généraux de l'élément à ses déplacements

nodaux

On exprime maintenant les coefficients de la fonction de déplace ment {α} en fonction

des déplacements nodaux {δe} et représentent dans l'équation II.3.2 on relie les déplacements
des points de l'élément aux déplacements nodaux {δe}.Cela donne par exemple pour le
nœud1:

1e     x1, y1  f  x1, y1  (II.3.3)

En procédant de façon similaire pour tous les autres nœuds on obtient pour n nœuds par
exemple :

1e    f  x1, y1 



 
    f  x , y   
  2     2 2    

e  
e


   
n  f  x n , y n  

e

   
(II.3.4)

Ou

 e    (II.3.5)

Cela donne :

 
    -1  e
(II.3.6)

On remplace {α} par sa valeur dans l'équation II.3.2 et on obtient relation cherchée entre les
déplacements {δ (x, y)} de tous les points de l'élément et les déplacements nodaux :

  x , y   f  x , y   -1 e 

(II.3.7)

15
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

II.3.4. ETAPE 4 : Relations entre déplacements et déformations cas général

On relie maintenant les déformations ε(x, y) du point (x, y) de l’élément aux

déplacements δ(x, y) et donc aussi aux déplacements nodaux {δe}.

Les déformations aux différents points de l'élément peuvent être obtenues à partir de la
fonction de déplacement par une différentiation dont la forme exacte dépend du type de
problème traité. Par exemple, pour un problème d'élasticité plane, les déformations
correspondent aux dérivées premières des déplacements, tandis que pour les problèmes de
flexion, les déformations sont associées à la courbure de l'élément et correspondent aux
dérivées secondes. En général

ε  x, y   {différentielle de   x , y } (II.3.8)

On obtient la forme exacte de cette expression pour chaque classe de problème à partir
de la théorie de l'élasticité. En utilisant l'expression de l'équation II.3.6 pour {δ(x, y)} et en
remarquant que [A]-1et {δe} sont indépendants de x et. y, le vecteur de déformation le
ε(x,y),est donné par :

ε  x, y   différentielle de f  x , y   -1  e 

(II.3.9)

En posant la matrice différentielle de f (x, y) égale à [C], cette équation peut être écrite
de la façon suivante :

ε  x, y   C  -1  e 
(II.3.10)

Où, en général, [C] contiennent des termes en x, et y.

C'est la relation cherchée entre les déformations en chaque point de l'élément et les
déplacements nodaux. Avec :

C  -1   B  (II.3.11)

Cette relation devient :

ε  x, y   B  e 
(II.3.12)
16
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

II.3.5. ETAPE 5 : Relations entre déformations et contraintes

Les contraintes internes à l'élément {ζ(x,y)} sont maintenant reliées aux déformations
le {ε(x,y)}.

Puisque l'on connait déjà une relation entre les déformations internes et les déplacements
nodaux {δe}, on peut relier les contraintes internes {ζ (x, y)} à ces derniers. Dans cette étape,
il est clair que les propriétés élastiques de l'élément devront être considérées. En général,

(x , y )   D ε  x, y 
(II.3.13)

Où [D] est la matrice d'élasticité et contient les propriétés élastiques de l'élément, c'est-à-dire
des quantités telles que le module d'Young' d'élasticité transversale E et le coefficient de

   
Poisson v. Puisque l'on sait à partir de l'équation II.3.12 que ε  x, y    B  
e

 
(x , y )   D  B   e
(II.3.14)

II.3.6. ETAPE 6 : Relié les charges nodales aux déplacements nodaux

Les contraintes internes {ζ(x, y)} sont maintenant remplacées par les charges nodales
{Fe} statiquement équivalentes et ces dernières sont reliées aux déplacements nodaux,
définissant ainsi la matrice de rigidité [Ke] de l'élément,

On utilise le principe des travaux virtuels pour déterminer l'ensemble des charges nodales qui
est statiquement équivalent aux contraintes internes. La condition d'équivalence peut être
exprimée comme suit : durant n'importe quel déplacement virtuel imposé à l'élément, le
travail extérieur total produit par les charges nodales doit être égal au travail interne total des
contraintes. On choisit un ensemble arbitraire de déplacements nodaux représenté par le
vecteur {δ*e} où

 1*e  
 
 *e    2  
  *e
(II.3.15)
 
 *e 
 n 
Le travail extérieur Wext des forces nodales est donné par :

17
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

 F    F   ....   F 

W int  1*e 1
e
2
*e
2
e
n
*e
n
e

   F 
*e  e

(II.3.16)

Si les déplacements arbitrairement imposés produisent des déformations {ε (x, y)*}

aux points de l'élément où les contraintes réelles sont les {ζ (x, y)}, alors le travail interne par
unité de volume est donné par :

    x , y 
T
W int   x , y 
*

(II.3.17)

et le travail interne total est obtenu en intégrant sur le volume total de l'élément c'est-à-dire

    x , y d vol 

 r W int  d vol    r  x , y 
*

(II.3.18)

Maintenant, on sait par l'équation IV que les déformations en chaque point de

l'élément sont exprimées en fonction des déplacements nodaux par {ε(x, y)} = [B] {δe}. D'où,
les {δ*e} étant imposés, les déformations correspondantes s'expriment par :

  x , y     
* *e

(II.3.19)

De plus l'équation II.3.12 relie les contraintes réelles dans l'élément aux déplacements nodaux
par :

  x , y   D   e

(II.3.20)

L'opération finale consiste à égaler travail interne et travail externe produits pendant
les déplacements virtuels {δ*e}. Puisque le principe de base des déplacements virtuels est
valide pour n'importe quel système de déplacements appliqué, on peut choisir ce dernier
comme l'on veut. Dans le cas présent, il est commode de supposer les déplacements nodaux
égaux à l'unité. On obtient alors :

F     D d vol   

e v  e

(II.3.21)

En comparant l'équation II.3.21 avec l'équation l’écrite à nouveau ci-dessous,

18
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

F      
e e e
(II.3.22)

Il est clair que la matrice de rigidité [Ke] de l'élément est donnée par l'expression :

e   v   D d vol 

  (II.3.23)

Donc, pour évaluer la matrice de rigidité de l'élément dans le cas général, il est nécessaire
de calculer la matrice [B] donnée dans l'étape 4 à partir des matrices [A] -1 et [C], et la matrice
[D] donnée dans l'étape 5 ; puis il suffit de faire les multiplications et les intégrations définies
par l'équation II.3.21.

II.4. FORMULATION DE LA MATRICE DE RIGIDITE D’ÉLÉMENT FINI DE

POUTRE : [3], [4]

On considère un Élément de Poutre avec une section transversale uniforme faisant

partie d'une structure continue. Pour un tel élément on peut utiliser le système de coordonnées
et la numérotation de la figureII.7.a

Fig. II.7.a :Déplacements transversaux et de rotation à chaque extrémité de l'élément de

poutre EP.

L’élément de poutre de la figure II.7.b La poutre est chargé par des forces et des moments en
chaque nœud et on suppose qu'elle a une rigidité uniforme à la flexion EI.

19
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Fig. II.7.b :Position déformé et les forces et des moments en chaque nœud de poutre EP.

On associe aux déplacements transversaux et de rotation à chaque extrémité de

l'élément des forces nodales correspondantes, le moment Ty et la force de cisaillement Fz (fig.
II.7.b).

Les vecteurs pour les déplacements et les forcés au nœud 1 peuvent être respectivement écrits
de la façon suivante :

 
w  F 
 
 e 
1

 1
, F e   
1
1

 y 1
 
  1
 (II.4.1)

Les vecteurs complets des déplacements nodaux et des charges nodales pour la poutre 1-2
prennent donc la forme

w  F 

   
1 z1
  e     Fe   
 y 1   y 1 
  
e  
1  

 e
 
 1
, F    
 
  e  w 2 
 2  

  2 
 Fe   Fz 2 
  
  2   y 2 
(II.4.2)

Le déplacement de n'importe quel point de l'élément peut être défini par les deux
composantes w et θy de déplacement et de rotation.

Ainsi le vecteur déplacement est donné par :

20
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

 w 
  x , y    
 y  (II.4.3)

Puisque l'élément possède quatre degrés de liberté (w1, θy1, w2 et θy2), quatre coefficients
inconnus doivent apparaître dans le polynôme représentant le modèle du déplacement.
Supposons que w soit donné par l'équation II.4.4.

w  x  x  x 2  x 3 (II.4.4)
1 2 3 4

Oùα1,α2, α3 et α4 sont des coefficients non encore déterminés. Puisque θy = dw/dx, on obtient:

    2 x  3 x 2
y 2 3 4 (II.4.5)

Ainsi à partir des équations II.4.4 et II.4.5 on obtient le vecteur de déplacements

{δ(x,y)}donné par l'équation II.4.6.

 
 1
 w   1 x x 2 x 3   2 
  x , y        
 y  0 1 2x 3x  3 
 2
 

 4  (II.4.6)

qui définit la matrice [f(x, y)] et le vecteur {α} de l'équation II.3.2.

Pour le simple élément de poutre de la figure II.6, les coordonnées des nœuds sont 0 et
L. L'expression choisie pour définir {δ (x, y)} (équation II.4.4) est :

w    x  x 2  x 3
1 2 3 4
dw
     2 x  3 x 2
y dx 2 3 4 (II.4.7)(II.4.8)

Au nœud 1, x = 0 et donc

w   et   
1 1 y1 2

21
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Au nœud 2, x = L et donc

w     L   L2   L3
2 1 2 3 4
    2 L  3 L2
y2 2 3 4

Cela donne sous forme matricielle :

 1 0 0 0  w
    1 
 1  0 1 0 0  
 2   3 2 3 1  y 1 
  - - -  
3   L2 L L L w 2 
   1   
 4  
2 1 2
-   y 2 
 L3 L2 L3 L2  (II.4.9)

La matrice   -1 correspond à la matrice carrée d'ordre 4 de l'équation II.4.9.

Dans ce cas la seule «déformation» à considérer est la courbure autour de l'axe des y,
Donc, à partir de l'équation II.4.4, le vecteur de déformation est donné par :

-d 2w
  x , y   2
 -2 - 6 x
3 4
dx (II.4.10)

Et sous forme matricielle :

 1 
 
 
ε  x, y  =  0 0 -2 -6x   2 
3 
 4 
(II.4.11)

On a établi à l'étape 3 que {α} = [A]-1{δ}. En remplaçant {α} par sa valeur dans l'équation
II.4.11on obtient alors :

22
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

 1 0 0 0 
 0 1 0 0  w 1 
  
1   y 1 
ε  x, y  =  0 0 -2 -6x    3
 L2 
2 3
  w 
 L L2 L  2 
 2 1 2 1   y 2 
 L3 L2 L3 L2  (II.4.12)

Ce qui donne après multiplication matricielle

w 1 
 
  y 1 
  x , y    L6  12Lx
2 3
4 6x
-
L L2
-
6 12x
+
L3 L3
2 6x

L L2  w 
 2
 y 2 
(II.4.13)

On a ainsi obtenu la matrice [B] de l'équation II.3.13.

 B    
6 12x 4 6x 6 12x 2 6x
 3 - - +  
L
2
L L L2 L3 L3 L L2 (II.4.14)

Donc, dans ce cas particulier, la matrice [D] contient un seul terme correspondant à la rigidité
de flexion EI. Cependant, en général, [D] est d'ordre plus élevé. L'équation II.3.12 devient
donc :

w 1 
 
  y 1 
  x , y     L6  12Lx
2 3
4 6x
-
L L2
-
6 12x
+
L3 L3
2 6x

L L2  w 
 2
 y 2 
(II.4.15)

Toutes les matrices de cette expression ont déjà été calculées explicitement. La matrice [Ke],
tirée de l'équation II.3.23 peut donc être obtenue par

 6 12x 
 2- 3 
 L L 
 4 6x 
L  
-
 L L2      6  12x 4 6x 6 12x 2 6x 
K       L2 L3 - - +   dx
6 12x L L2 L3 L3 L L2
- 2  3 
0

 L L  (II.4.16)
 2 6x 
 - 
 L L2 

23
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

En effectuant les calculs, on obtient l'équation II.4.17et l'intégration donne finalement la

matrice [Ke ] par l'équation II.4.18.

 36 144x 144x 2 24 84x 72x 2 -36 144x 144x 2 12 60x 72x 2 

 L 4 - L5  L 6 -
L3 L 4
 5
L L4
 5 -
L L6
-
L3 L 4
 5 
L 

 24 84x 72x 2 16 48x 36x 2 -24 84x 72x 2 8 36x 36x 2 
L  -  5 -   4 - 5 -  4 
        L 3 L 4 L L 2 L 3 L 4 L 3 L L L 2 L3 L
e
d
0  36 12 60x 72x 2 
2 2
144x 144x 24 84x 72x 36 144x 144x 2
 L4  -  - -   - 5 
L5 L6 L3 L 4 L5 L 4 L5 L6 L3 L 4 L
 2

 12 - 60x  72x 8 36x 36x 2 -12 60x 72x 2 4 24x 36x 2 
 L3 L4 -  4  4 - 5 -  4 
L5 L 2 L3 L L3 L L L 2 L3 L
(II.4.17)

 12 6L -12 6L 
 6L 4L2 -6L 2L2 

   
e

  L3  -12 -6L 12 -6L 

 
 6L 2L2 -6L 4L2  (II.4.18)

II.5.ÉLÉMENTS FINIS POUR L'ÉTUDE DE L'ÉLASTICITÉ PLANE

Les problèmes d’élasticité plane peuvent être divisés en deux groupes distincts, les
problèmes de contraintes plans et les problèmes de déformations plans. Dans le premier cas,
les structures à une petite épaisseur par rapport à ses autres dimensions et, les contraintes
normales au plan sont négligées ; tandis que dans le deuxième cas la structure à une
dimension très grande par rapport à ses autres dimensions et la déformation normale au plan
de la charge est nulle.

On emploie les étapes de base pour calculer la matrice de rigidité des éléments
rectangulaire et triangulaire dans le cas de contraintes et de déformations planes.

II.5.1.ÉlémentRectangulaire à Quatre nœuds Q4, (l’Elément bilinéaire

quadratique) :
On suppose que l'élément rectangulaire Q4 a des côtés de longueur a et b et une
épaisseur t (fig. II.7.a). On utilise le système de numérotation des nœuds de la figure II.7.a.

24
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Puisque pour un problème d'élasticité plane l'élément a deux degrés de liberté en

chaque nœud, chaque élément a huit degrés de liberté. La figure II.7.b montre les huit
déplacements inconnus.

y
a
v2 v4
2 4 2 4
u4
u2
b

u1
x u3
1 1 3
v1 v3
Avec épaisseur t 3

Fig .II.8 (a) Fig. II.8 (b)

Fig. II.8.a : Elément quadrilatère bilinéaire Q4,a et b : Déplacements nodaux.

et la figure II.8.c les forces nodales correspondantes. En utilisant la notation matricielle on

peut écrire les déplacements au nœud 1 de la façon suivante :

F 
F1  Fx 1 
 y1  (II.5.1)

2 fy2 fy4
4
fx4
fx2

fy1
fx3
1 fx1 3
fy3

Fig. II.8(c)

Fig. II.8.b : Elément quadrilatère bilinéaire Q4, c : Forces nodales.

25
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Pour un problème d'élasticité plane, l'état du déplacement en tout point (x, y) de

l'élément peut être représenté par deux composantes, c'est-à-dire :

  x , y   v 
u
  (II.5.2)

Puisque l'élément a huit degrés de liberté, on doit avoir huit coefficients inconnus dans
le polynôme représentant le modèle de déplacement. Deux telles fonctions sont données par
l'équation II.5.3.

u  1   2 x   3 y   4 xy 

v   5   6 x   7 y  8 xy 
(II.5.3)

On peut écrire l'équation II.5.3 sous forme matricielle :

1 
 
 2
 3 
 
u  1 x y xy 0 0 0 0   4 
   
v  0 0 0 0 1 x y xy   5 
 6 
 
 7 
 
 8 (II.5.4)

Ce qui peut être résumé par l'équation générale II.3.2.

  x , y   f  x , y   (II.3.2)

26
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

1 0 0 0 0 0 0 0 
0 0 0 0 1 0 0 0
 
  
 
1 0 b 0 0 0 0 0 
0 0 0 0 1 0 b 0
 
  
    1 a 0 0 0 0 0 0   
e

 
0 0 0 0 1 a 0 0
  
 
1 a b ab 0 0 0 0 
 
 0 0 0 0 1 a b ab 
 
(II.5.5)

Ce qui peut se réduire à :

    
e
(II.3.5)

Où {α} est le vecteur des coefficients inconnus des polynômes donnés par l'équation
II.3.5. Ainsi l'équation II.3.5 définit la matrice [A] pour le cas particulier d'un élément
rectangulaire pour des problèmes d'élasticité plane.

On multiplie alors les deux membres de l'équation II.5.5 par l'inverse de [A] pour
obtenir les coefficients {α} :

      e 
1

(II.3.6)

 1 0 0 0 0 0 0 0
 
 -1 1
0 0 0 0 0 0 
 a a 
1    u1 
   1 1   
 2  
0 0 0 0 0 0
ab b  1
    u 2 
 3  0 1 0 0 0 0 0 0
  
 4   -1 1  2
  0 0 0 0 0 0  u 3 
 5   a a
 
   -1 1   
 6  0 0 0 0 0 0  u 
3
 7   a a
 4
    
4 
-1 1
 8 0 0 0 0 0 0
 a b 
 1 -1 -1 1 
 0 0 0 0 
 ab ab ab ab  (II.5.6)

27
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Définissant ainsi l'inverse de la matrice [A], c'est-à-dire [A]-1. On peut vérifier le résultat en
calculant le produit [A]-1 [A] pour s'assurer qu'il est bien égal à la matrice identité [I].

Il est évident que la relation entre les déformations et les déplacements en tout point
pour un problème d'élasticité plané est indépendante de la forme de l'élément choisi. Ainsi :

 x  u x 

 y  v y 
 xy  u y  v  x 
(II.5.7)

En remplaçant u et v par leur valeur à partir de l'équation II.5.3, on 1obtient l'expression

suivante pour les déformations en tout point de l'élément :


x  1   2 x   3 y   4 xy    2   4 y
x

y   5   6 x   7 y  8 xy    7   8 x
y
 
 xy  1   2 x   3 y   4 xy    5   6 x   7 x   8xy 
y x
 3   4x   6  8 y
(II.5.8)

α2 + α4 bα2 + α4b α2 + α4b

2 4

x
1 α2 α2 α2 3

Fig. II.9 : Déformations dans l’élément fini rectangulaire

Il est d'un intérêt particulier de noter que ces déformations varient linéairement dans
l'élément. En considérant la variation de εx par exemple, elle est indépendante de x et dépend
linéairement de y, variant ainsi dans l'élément comme le montre la figure II.9. De même εy est
indépendant de y et linéairement dépendant de x tandis que γxy varie linéairement à la fois en
x et en y.

28
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

En utilisant ces expressions des déformations dans l'équation II.5.8, on obtient :

 x   2  4 y 
   
  x , y   y   7  8x 
     x     y 
 xy   3 4 6 8 
(II.5.9)

 1 
 
 2
 3 
0 1 0 y 0 0 0 0   
 
=  0 0 0 0 0 0 1 x   4 
 5 
0 0 1 x 0 1 0 y   
6
 
 7 
 
 8 (II.5.10)

Et avec la valeur de {α} donnée par l'équation II.3.6,

  x , y   C    
1 e

(II.3.10)

Ou

  x , y     e
(II.3.12)

Où
  C  
1

(II.3.11)

L’équation II.3.12 relie les déformations en chaque point aux déplacements nodaux. La
matrice [B] est établie à partir de l'équation II.3.11. En prenant la matrice [C] définie par
l'équation II.5.10et [A]-1définie par l'équation II.3.6, on obtient la matrice [B] :

B  
 1 y y 1 y y 
   0 0  0 0 
a ab ab a ab ab
 
 0 
1 x
 0
1 x
 0
x
 0
x 
 b ab b ab ab ab 
 
 1  x 1 y

1 x

y x 1 y

x y 
 b ab a ab b ab ab ab a ab ab ab  (II.5.11)

29
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Les composantes des contraintes et des déformations sont reliées par la matrice [D] où
{σ(x,y)}=[D] {ε(x, y)}. Pour des problèmes d'élasticité plane, [D] est une matrice carrée
d'ordre 3 dont les termes dépendent du problème : contraintes ou déformations planes.

La matrice [D] requise pour l'équation générale II.3.20., c'est-à-dire :

  x , y   D   e

(II.3.20)

Qui donne la relation entre les contraintes en tout point et les déplacements nodaux de la
façon suivante :

Pour l'état de contrainte plane :

 
1  0 
E  
D   
(1- 2 ) 
1 0  (II.5.12)
1- 
0 0 
 2 

Etpour l'état de déformation plane.

 
1-  0 
E  
D    1-
(1- )(1- 2 ) 
0  (II.5.13)
1- 2 
 0 0 
 2 

La matrice [B] étant définie par l'équation II.5.11 et la matrice [D] par l'équation II.5.12. D'où,
dans l'évaluation de [Ke] selon l'équation

e      D d vol 



(II.3.23)

qui pour un élément d'épaisseur constante devient :

e   t     D dx dy



(II.5.14)

Le produit [B]T[D] [B] doit être évalué en premier, et les termes de la matrice qui en résultent
doivent être intégrés sur la surface de l'élément. Donc :

30
 Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

e   t a0 b0 (x , y ) . D .(x , y ).dx . dy


(II.5.15)

0 0 0

1 0 0 
0 0 1  
  1  0  0 1 0 y 0 0 0 0 
a b y 0 x E   
   t  0  0
e
 0  0 0 0 0 0 0 1 x  dx . dy
0  (1- 2 ) 
  1
0 0
  1-  0 0 1 x 0 1 0 y 
0 0 1 0 0 
 2 
0 1 0 
 
0 x y 

(II.5.16)

Tout calcul fait on obtient finalement la matrice de rigidité élémentaire [Ke] :

0 0 0 0 0 0 0 0 
 
0 ab 2
a b 
2
ab 0 0 0  ab
 2 2 
 
0 (1  )ab (1  )a 2b (1  )ab (1  )ab 2

0 0 0
 2 4 2 4 
 
0 ab 2 (1  )a 2b 1  3 (1  )a 3b  (1  )a 2b  ab 2
(1  )a b 
2 2

 ab   0

E .t 2 4 3 2  4 4 8
 e  
  1  2 0 0 0 0 0 0 0 0


 2 
0 (1  )ab (1  )a b2
(1  )ab (1  )ab 
0 0 0
 2 4 2 4 
  ab 2 a 2b 
0  ab 0 0 0 ab 
 4 2 
 3 
0  a 2b (1  )ab 2 (1  )a 2b 2 (1  )ab 2 a b 1  3 (1  )ab  
2
0 a b  
 2 4 8 4 2 3 2  

(II.5.17)

II.5.2. Élément Rectangulaire à 8 nœuds Q8,(Elément quadrilatéral

quadratique).[3] [5] :

L’élément Q8 est un quadrilatère (rectangle) à huit nœuds. Deux degrés de liberté à chaque
nœud (u, v déplacements)a totale c’est 16 DDL, cet élément basé sur le champ de déplacement
dans les directions x et y, sont présentés dans la figure II.10.

En termes de coordonnées généralisées αi, le champ de déplacement est donné par:

31
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

u   1   2 x   3 y   4 x 2   5 xy   6 y 2   7 x 2 y   8 xy 2
v   9   10 x   11 y   12 x 2   13 xy   14 y 2   15 x 2 y   16 xy 2 (II.5.18)

y, v
a
v4
v3
3 7 4 u4
u3
8 6 b

u1 x, u
1 5 2 u2
v1
v2
[Link].10 : Elément quadratique quadrilatéral Q8

La déformation de ce l’élément c’est :

 x   2  2  4 x   5 y  2  7 xy   8 y 2
 y   11   13 x  2  14 y   15 x 2  2  16 xy
(II.5.19)
 xy   3   10   5  2  12  x  2  6   13  y   7 x  2  8   15 xy   16 y
2 2

Pour déterminé la matrice de rigidité [Ke] de l'élément à champs de déplacement Q8 en doit

calcule l’intégrale d’après l’équation (II.3.23) :

e   v   D d vol 

  (II.3.23)

e   t a0 b0 (x , y ) . D .(x , y ).dx . dy



(II.5.15)

La matrice de déformation [B] pour cet élément est :

0 1 0 2x y 0 2 xy y2 0 0 0 0 0 0 0 0 
 
 B   0 0 0 0 0 0 0 0 0 0 0 1 x 2y x2 2xy 
0 0 1 0 x 2y x 2 2xy 0 1 0 2x y 0 2xy y 2 

(II.5.20)

32
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Et [D] c’est la matrice d'élasticité (II.5.12) et (II.5.13).

[Link]ésentation de l’Élément T6, Élément Triangulaire à 6 nœuds (ou

LST linear strain triangle) [4] :

L’élément T6 est représenté dans la figure II.11. Il présente, à part les trois nœuds 1, 2
et 3, des nœuds supplémentaires situés au milieu de chaque côté du triangle. Les degrés de
liberté sont en nombre de 12 (deux degrés de liberté ui, vi pour chaque nœud, i = 1...6).
y, v
v3

3 u3
v6 v5
u5
6 u6 5
v2
v1 v4
u2
u4 2
u1 4
1
x, u

[Link].11 : Elément de type T6 et ses 12 degrés de liberté

En terme de coordonnées généralisées αi, le champ de déplacement est décrit par:

u   1   2 x   3 y   4 x 2   5 y 2   6 xy
(II.5.21)
v   7   8 x   9 y   10 x 2   11 y 2   12 xy

et respectivement, pour le champ de déformations:

u
x    2  2 4 x   6 y
x
v
y    9  2 11 y   12 x (II.5.22)
y
u v
 xy     3  2 5 y   6 x   8  2 10 x   12 y
y x

En utilisant les polynômes de déplacement ci-dessus la matrice [A] de l’équation :

33
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

 e    (II.3.5)

Sera :

1 x 1 y 1 x 12 y 12 x 1 y 1 0 0 0 0 
0 0
 
1 x 2 y 2 x 22 y 22 x 2y2 0 0 0 0 0 0 
1 x 3 y 3 x 32 y 32 x 3y3 0 0 0 0 0 0 

1 x 4 y 4 x 42 y 42 x 4y4 0 0 0 0 0 0 
 
1 x 5 y 5 x 52 y 52 x 5y5 0 0 0 0 0 0 
 
1 x 6 y 6 x 62 y 62 x 6y6 0 0 0 0 0 0 
 A   
0 0 0 0 0 0 1 x 1 y 1 x 12 y 12 x 1 y 1 
 
0 0 0 0 0 0 1 x2 y2 x2 y2 x 2y2
2 2

 
0 0 0 0 0 0 1 x 3 y 3 x 32 y 32 x 3 y 3 
 
0 0 0 0 0 0 1 x 4 y 4 x 42 y 42 x 4 y 4 
0 0 0 0 0 0 1 x 5 y 5 x 52 y 52 x 5 y 5 
 
0 0 0 0 0 0 1 x 6 y 6 x 62 y 62 x 6 y 6 
(II.5.23)

Après avoir les champs de déplacement, la matrice de rigidité de l'élément triangulaire peut
être évaluée par l'expression générale suivant :

K e      B   D  B  dx dy
T

(II.3.23)

La matrice de déformation [B] pour cet élément est :

 0 1 0 2x 0 y 0 0 0 0 0 0 
 B   0 0 0 0 0 0 0 0 1 0 2 y x 
(II.5.24)
0 0 1 0 2y x 0 1 0 2x 0 y 

Et [D] c’est la matrice d'élasticité (II.5.12) et (II.5.13)

Remarque : dans le chapitre quatre le programme d’éléments finis qui va calculer l’intégrale
de l’équation (II.3.23). Il suffit d’entré les valeurs des matrices [B] T, [D] et [B] de chaque
élément pour déterminé la matrice de rigidité élémentaire [Ke].

34
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

II.6. ASSEMBLAGE DES MATRICES DE RIGIDITÉ DES ÉLÉMENTS [4]

Quand la matrice de rigidité [Ke] est déterminée pour un élément typique dans la
structure, les équations pour tous les éléments doivent être assemblés pour obtenir un système
d’équations de la forme :

F    K G .G  (II.6.1)

Qui relie les déplacements et les charges de la structure complète par l’intermédiaire de la
matrice de rigidité.

Dans l’équation (II.6.1), {F} est un vecteur de forces extérieures appliquées aux points
nodaux de la structure. Chaque force extérieure est la somme des forces appliquées aux nœuds
adjacents. Le vecteur {δG} dans (II.6.1) représente les déplacements aux nœuds ; chaque
déplacement dans le vecteur {δG} a la même valeur que les déplacements des nœuds
adjacents. Les termes dans la matrice [KG] sont la somme des termes appropriés dans la
matrice de rigidité des éléments qui ont un déplacement commun.

Pour illustré ce principe considérons l’exemple suivant où nous proposons d’assembler deux
éléments poutres.

[Link].12 : deux éléments de poutre 1 et 2 avec L1 et L2

Remarquons que ces deux éléments ont le nœud 2 en commun. Soient [K1] et [K2] les
matrices de rigidité des éléments 1 et 2.

Élément 1 :

12 6L1 -12 6L1 

 
  6L1 4L12 -6L1 2L12 
   3
1

L1 -12 - 6L1 12 - 6L1  (II.6.2)

 
 6L1 2L12 -6L1 4L12 

35
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

Élément 2 :

12 6L 2 -12 6L 2 
 2
  6L 2 4L 2 -6L 2 2L 2 
2

   3
2

L 2 -12 - 6L 2 12 - 6L 2  (II.6.3)
 
 2 
2 2
6 L 2 L 2 -6 L 2 4 L 2 

Ainsi la matrice de rigidité totale sera :

 12EI 6EI -12EI 6EI 

 L3 0 0 
L12 L13 L12
 1

 6EI 4EI -6EI 2EI 
 L2 0 0 
L1 L12 L1
 1 
 -12EI -6EI 12EI 12EI -6EI 6EI -12EI 6EI 
 L3   2
L12 L13 L 23 L12 L2 L 23 L 2 2 
 K G    1 (II.6.4)
6EI 2EI -6EI 6EI 4EI 4EI -6EI 2EI 
 2
 2  
 L1 L1 L12 L2 L1 L2 L22 L22 
 -12EI -6EI 12EI -6EI 
 0 0 
 L 23 L22 L 23 L22 
 6EI 2EI -6EI 4EI 
 0 0 
 L22 L2 L22 L2 

On remarque que la matrice de rigidité de deux éléments est une matrice symétrique carrée
d’ordre six par six (6x6). Cet ordre est dû au fait que l’élément assemblé possède trois nœuds
donc six degrés de liberté car chaque nœud à deux degrés de liberté.

Pour un assemblage de trois, quatre éléments ou plus d’élément E-P, Q4, Q8, T6 on
procède de la même façon.

II.7. DETERMINATION DE MATRICE DE MASSE DES ÉLÉMENTS : [4]

On construit la matrice de masse élémentaire (d’un élément particulier donné) à partir de

l’expression de l’énergie cinétique :
1
Ec =
2  ( v(x,y,z))2dm (II.7.1)

v(x,y,z) : est le vecteur vitesse en du point courant.

La formulation générale qui détermine la matrice masse [Me] de l’élément est représenté
comme suit:

36
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

T 
M e    A 1  v  (x , y , z ) .  (x , y , z )  dv . A 1  (II.7.2)
     

II.7.1. La matrice masse d’un élément de poutre :

D’après l’équation (II.7.2) la matrice masse [Me] est :
T 
M e    A 1  v  (x ) .  (x ) dv . A 1  (II.7.3)
     

Avec :

 (x )  1 x x 2 x 3  (II.7.4)



Et :

1 0 0 0
0 1 0 0 

 A 1    3 2 3 1  (II.7.5)
   2
L L L2 L 
 2 1 2 1 
 L3 L 2 L3 L 2 

Donc :

 3 2
1 0
L2 L3 
  1 
0 2 1 
1 x 
 L L2 
b h L
 
M e    
  3       1
2  0 0 0 x 2  
x x 2 x 3  dx .dy .dz .  A 1  (II.7.6)
0 0
 L2 L3  x 3 
 
 1 1 
0 0 
 L L2 

L 3L2 L3 L4 
2 20 15 28  1 0 0 0
 2 
 L L3 L 4
L  0
5
 1 0 0 
 12 30 60 105   3 2 3 1 
 bh   (II.7.8)
 L 7L2 4L3 3L4   L2 L L2 L 
 2  
20 15 14   2 1 2 1 
 2 
 L  L3 L4 L5   L3 L2 L3 L2 
 12 20 30 42 

Ainsi :

37
Chapitre II : Formulation de quelques éléments basés sur le champ de déplacement

156 22L 54  13L 

 
bhL  4L2 13L  3L2 
M e   (II.7.9)
  420  156  22L 
 
Symétrique 4L 2 

II.7.2. La matrice masse d’un élément rectangulaire :

D’après l’équation (II.7.2) on a :
T 
M e    A 1  v  (x , y )  .  (x , y )  dv . A 1  (II.7.10)
         

Avec :

1 x y xy 0 0 0 0 
 (x , y )   
1 x y xy 
(II.7.11)
0 0 0 0

T
Pour  A 1  et  A 1  voir le titre II.5.1.
   

Ainsi on obtiendra :

1 1 1 1 
 9 18 36 18 0 0 0 0 
 
1 1 1 1 0 0 0 0 
18 9 18 36 
 
1 1 1 1 0 0 0 0 
 36 18 9 18 
1 1 1 1 
 0 0 0 0 
 M   .h .a.b 18 36 18 9
e  (II.7.12)
   1 1 1 1 
0 0 0 0 
 9 18 36 18 
 1 1 1 1 
0 0 0 0
18 9 18 36 
 
0 0 0 0
1 1 1 1 
 36 18 9 18 
 
0 0 1 1 1 1 
0 0
 18 36 18 9 

On doit appliquer la même procédure pour déterminé la matrice masse des éléments
Q8 et T6.

De la même manière qu’avec les matrices de rigidité, on construit la matrice de masse

globale par assemblage des matrices de masse élémentaires.

38
Chapitre III : Analyse dynamique des structures par éléments finis

CHAPITRE III :
ANALYSE DYNAMIQUE DES STRUCTURES PAR ÉLÉMENTS
FINIS

III.1. INTRODUCTION :

Ce chapitre est consacré au Analyse dynamique des structures par la méthode

d’éléments finis, étude vibratoire sur un élément unidimensionnel, élément poutre. Et sur un
élément bidimensionnel, élément rectangulaire.
Et aussi avec une petite rappelle sur la vibration d’une poutre parla méthode
analytique.

III.2. DYNAMIQUE D’UNE POUTRE PAR LA METHODE

ANALYTIQUE

Ce titre est une synthèse bibliographique des références [13], [14], [15], [16], [17],
[18], [19].

Ces titres sont consacrés à la présentation des théories des vibrations libre et forcée
des poutres, un modèle analytique simplifié. Pour une étude comparative avec les résultats ont
été validés par la méthode des éléments finis.

III.2.1.Déformée d’une poutre en flexion dans l’hypothèse de la résistance

des matériaux.
Une poutre élancée rectiligne d’axe x, de longueur L et de section droite d’aire S,
hauteur h et largeur b vérifiant h, b << L fléchit sous l’action d’un chargement linéique
transversal q(x) et prend une déformée y(x). Pour des chargements modérés, induisant une
déformée telle que le déplacement transversal reste petit devant les dimensions transversales
de la poutre : y(x) << b, h, les sections droites restent droites (ne gauchissent pas) et tournent
simplement l’une par rapport à l’autre.

39
Chapitre III : Analyse dynamique des structures par éléments finis

[Link].1 : Géométrie d’une poutre

M(x) caractérisant le moment de flexion à l’abscisse x résultant du chargement q(x),

écrivons, dans cette hypothèse de flexion faible, l’équilibre mécanique d’un petit tronçon de
longueur dx sous l’action du moment M(x).

[Link].2 : le moment de flexion

En traçant au centre de la section droite terminale la parallèle à la section droite

d’entrée, l’angle a caractérisant la rotation relative des deux sections par rapport à l’état non

fléchi s’écrit sous la forme : soit Le rapport n’est autre que la

déformation d’allongement de sorte que la déformation d’allongement des fibres de la

poutre s’écrit : caractérisant le Module d’Young du matériau constitutif de la

poutre, la contrainte de traction s’écrit :

La force résultante F induite par ces contraintes et le moment de flexion résultant M

sont donnés par :
⁄
∬ ∫
⁄

40
Chapitre III : Analyse dynamique des structures par éléments finis

⁄
∬ ∫
⁄

étant le moment quadratique (couramment appelé moment d’inertie de flexion) de la

section droite par rapport à l’axe de flexion z. F=0 traduit l’absence de force appliquée et la
seconde relation exprime la proportionnalité entre la courbure locale 1/R de la déformée et le
moment de flexion appliqué M et constitue l’équation différentielle de la déformée.
Dans l’hypothèse des petits déplacements envisagés ici, la courbure

⁄
( ( ) )

L’équation différentielle de la déformée se réduit à :

Le signe - provient du fait que la déformée y(x) est repérée dans le référentiel x, y, z alors
que le moment de flexion M(x) est défini dans le trièdre de Frenet : tangente t, normale n et
bi normale r avec t=x, r=-z et n=-y

III.2.2.Équation différentielle de la déformée en terme de chargement

linéaire q(x)

Dans le trièdre de Frenet, écrivons l’équilibre mécanique du tronçon de longueur dx

soumis au chargement linéaire q(x).

[Link].3 :L’équilibre mécanique du tronçon de longueur dx soumis au chargement q(x).

L’équilibre des forces (nullité de la résultante) s’écrit :

41
Chapitre III : Analyse dynamique des structures par éléments finis

L’équilibre des moments par rapport à l’origine du trièdre s’écrit :

⃗⃗ (⃗ ⃗) ⃗⃗ ⃗⃗ soit en projection sur l’axe z :

soit en négligeant le terme du second ordre

En différentiant deux fois l’équation différentielle de la déformée écrite en

terme de moment fléchissant et en tenant compte des deux relations précédentes, l’équation
s’écrit finalement en termes de chargement linéique :

III.2.3. Équation différentielle des vibrations d’une poutre

[Link] LIBRES

[Link].4 : vibration libre d’une poutre cantilever

La masse m(x) étant égale à ρ.S, la force d’inertie linéique q(x) induite par la vibration
sera :

42
Chapitre III : Analyse dynamique des structures par éléments finis

Étant l’accélération induite par la vibration, d’où l’équation différentielle des vibrations
libres:

Compte tenu des hypothèses de flexion faible, la solution de cette équation ne sera
acceptable que dans la mesure où l’amplitude de vibration a reste petite devant l’épaisseur h :
a <<h.

En posant et en séparant les parties temporelles et

spatiales

L’équation différentielle donnant l’amplitude Y(x) de la déformée à la pulsation w s’écrit :

La solution de cette équation s’écrit sous la forme générale :

k1=k, k2=-k, k3=ik et k4= - ik étant les racines de l’équation

⁄
. Avec √

Elle s’écrit donc également sous la forme plus commode :

43
Chapitre III : Analyse dynamique des structures par éléments finis

III.[Link]. Modes propres de vibration d’une poutre cantilever de section rectangulaire

Avec S=bh et , les fréquences propres sont données par :

√
√

et les modes propres par :

Dans le cas d’une poutre encastrée, les conditions aux limites sont :

 Déplacement d’encastrement interdit

 Rotation interdite à l’encastrement

 { } Absence de

moment de flexion à l’extrémité libre

 { } Absence

d’effort tranchant à l’extrémité libre

Les constantes a, b, c, d sont calculées à partir des conditions initiales et des conditions aux
limites. On peut facilement faire apparaitre l'équation caractéristique permettant de trouver
les pulsations propres du système :

| || |

Ce système n’admet une solution a, b, c, d proportionnels à une constante arbitraire près que
si son déterminant est nul, soit l’équation aux pulsations propres définies par les conditions
aux limites :

Tracé l’équation caractéristique et recherche des zéros

44
Chapitre III : Analyse dynamique des structures par éléments finis

[Link].5 : courbe d’équation de caractéristique

Dont les solutions sont : α1 = k1.L =1.875, α2 = k2.L = 4.694, α3 = k3.L =7.855,…

Et les fréquences propres sont :

√ ( ) √

√ √

Donc les pulsations propres de vibration sont (modes fondamentales) :

√ ( ) √

√ ( ) √

√ ( ) √

45
Chapitre III : Analyse dynamique des structures par éléments finis

III.2.3.2. VIBRATIONS FORCÉE

Pour une excitation harmonique : , le problème à résoudre est :

[Link].6 : vibration forcée d’une poutre cantilever

 y 2 (L , t )
 y (0, t )  0  x 2  0

Avec C.L. :  y (0, t ) et  3
 0  y (L , t )  F cos t 
 X  x 3 0

Recherche d’une solution particulière de la forme :

y  X , t  Y  x  cos(t )

D’où

Avec

Dans le cas d’une poutre encastrée, les conditions aux limites sont :

 Déplacement d’encastrement interdit

 Rotation interdite à l’encastrement

 { } Absence de

moment de flexion à l’extrémité libre

46
Chapitre III : Analyse dynamique des structures par éléments finis

 {

}excitation harmonique à l’extrémité libre

a, b, c, d sont solution du système liné[Link] constantes a, b, c,et d sont calculées à

partir des conditions initiales et des conditions aux limites. On peut lourdement faire
apparaitre l'équation caractéristique permettant de trouver les pulsations propres du système.

III.3. ETUDE VIBRATLOIRE PAR ELEMENTS FINIS

Ce titre est une synthèse bibliographique des références [7] et, [4], [6], [20], [21].

III.3.1. ÉLEMENT UNIDIMENSIONNEL (1D) : Élément de poutre

On considéré un élément unidimensionnel (1D) solide, une poutre mince dans la figure
III.7, Les nœuds 1 et 2 d'extrémité sont soumis des forces verticalesF1, F2et des moments M1,
M2qui résulte des translations et des rotations. À chaque nœud il y a deux degrés de liberté.

III.3.1.1. Poutre avec des forces verticales et une charge

répartie uniforme :

L'élément a une longueur L, Est soumise à une charge répartie uniforme q

transversale par unité de longueur. L'équation de l'équilibre pour ce système est donnée par :

La variable continue w est approchée en termes de valeurs nodales discrets, nous avons choisi
d'écrire

{ }

Où = ∂w/∂x au nœud 1 et ainsi de suite. L’équation(III.21) pour être exacte en choisissant

les fonctions de forme cubique:

47
Chapitre III : Analyse dynamique des structures par éléments finis

[Link].7 : Un élément poutre mince

Le dérivé des fonctions de forme est égal à 1 dans un nœud spécifique et égale à zéro à tous
les autres nœuds.

La substitution dans (III.20) et l'application de la méthode de Galerkin’s conduit aux

équations à quatre éléments:

∫ { } { } ∫ { }

Théorème de Green est utilisé pour éviter différencier quatre fois, par exemple :

∫ ∫ ∫

On supposer que EI et q ne sont pas des fonctions de x, (III.23) deviennent :

48
Chapitre III : Analyse dynamique des structures par éléments finis

∫ { } ∫ { }

L’évaluation des intégrales donne :

{ }

[ ] { }

III.3.1.2. Poutre avec des forces verticales et une charge appliquée au

niveau des nœuds :

Le cas ci-dessus est pour une charge répartie uniforme, appliquée sur une poutre. Pour le cas
où la charge est appliquée uniquement au niveau des nœuds, nous avons :

{ } { }
[ ]

Qui représente la relation de rigidité d’élément. Par conséquent, en notation matricielle, nous
avons

[Link] de rigidité d’une poutre :

Si l'élément de la figure III.7a vibrant transversalement il serait soumis à une force

supplémentaire restituée : -ΡA( ), La matrice d'inertie ou de masse :

49
Chapitre III : Analyse dynamique des structures par éléments finis

∫ [ ] { }

L'évaluation des intégraux

[ ] { }

III.3.1.4. Poutre avec une force axiale :

Si l'élément de poutre à la figure III.7 est soumis à une force axiale supplémentaire P
([Link].8), une simple modification de l’équation (III.20), on obtient l'équation différentielle :

Le signe positif correspond à une charge de compression axiale.

[Link].8 : Poutre avec une force axiale

La discrétisation par éléments finis et l'application de la méthode de Galerkin's conduit à une

matrice supplémentaire associé à la force axiale :

∫ { }

L'évaluation de ces intégraux pour les forces P (compression) :

50
Chapitre III : Analyse dynamique des structures par éléments finis

[ ]{ }

III.3.2. ELEMENTS BIDIMENSIONELS (2D) : Contrainte et déformation

plane

L’élément décrit n'est pas d’éléments fini vrais parce qu'il utilisé uniquement pour
résoudre les équations différentielles dans une variable à une seule dimension. Ainsi, le
problème complexe impliquant deux ou trois variables d'espace 2D et 3D.
La situation change radicalement dans analyse des problèmes à deux ou trois
dimensions. Par exemple, considérons le mur plan la figure III.9 (a). Le mur a été subdivisé
en éléments rectangulaires de axb dans laquelle figure III.9 (b) est typique. Ces éléments ont
quatre nœuds aux coins

Les équations à résoudre sont les suivants :

[Link].9 : (a) maillage d’un mur. (b) Élément typique rectangulaire de quatre nœuds aux
coins

51
Chapitre III : Analyse dynamique des structures par éléments finis

(i) l'équilibre

Où σx, σy et τxy sont les seuls contrainte non nuls et Fx, Fy sont les forces du corps, par unité
de volume.

(ii) contrainte-déformation (contraintes planes)

{ } [ ]{ }

Où E est le module d'Young, v est le coefficient de Poisson et εx, εy, et γxy sont les
composants indépendants de faible déformation.

(iii) contrainte - déplacement

{ }

{ }
[ ]

Où u, v sont les composantes du déplacement dans le sens x, y.

En utilisant la notation de la section précédente avec A comme l'opérateur de contrainte-

déplacement et D la matrice contrainte-déformation, ces trois équations deviennent

Où

52
Chapitre III : Analyse dynamique des structures par éléments finis

{ } [ ]

[ ]

En suite :

En écrivant (III.39) :

{ }

{ }

C’est une paire d'équations différentielles partielles dans l'espace avec des variables continues
u et v.

Ceux-ci peuvent être résolus par la discrétisation sur chaque élément à l'aide des fonctions de
forme :

{ }

et

{ }

Dans le cas de l'élément représenté sur la figure rectangulaire III.9 (b) les fonctions de forme
Ni est :

53
Chapitre III : Analyse dynamique des structures par éléments finis

( )( )

( )

( )

Ce résultat appelé le rectangle de linéaire déformation, fonction de Lagrange.

54
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

CHAPITRE IV :
VALIDATION NUMÉRIQUE DES ÉLÉMENTS À CHAMPS DE
DÉPLACEMENT EN ANALYSE VIBRATOIRE

IV.1. INTRODUCTION [4], [22] :

La détermination des modes propres de vibration d’une structure par la technique des
éléments finis conduit à une relation du type suivant :

K X i  = i M  X i  (IV.1)

Avec
K  : La rigidité du système,

M  : La matrice des masses,

X i  : Le vecteur des déplacements de la structure définissant l’i-ème mode propre de
vibration, l’amplitude de déplacement,
λi =  2 : le carré de la pulsation correspondante ou la fréquence angulaire.

Résoudre un problème de valeurs propres consiste à trouver des couples λi{Xi} qui satisfont la
relation (IV.1). A chaque valeur λi correspond un vecteur propre {Xi}.

Mais le plus souvent, l'équation des valeurs propres va décrire une situation physique
telle que les vibrations libres ou forcées d’un solide ou d’un fluide.
Pour un solide élastique en vibration libre l’équation est :

 K e

 -  2 .  M e  . X   0 (IV.2)

Chaque élément de la structure aura une équation de la forme (IV.2). D’où en combinant
l’ensemble de ces équations on obtiendra le système d’équations suivant pour toute la
structure :

 K    .M   . X   0
G
2
G (IV.3)

55
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

Ou {X} et [MG] sont respectivement le vecteur de déplacement et la matrice masse pour tout
la structure.

La structure considérée qui est un corps élastique ayant un nombre infini de degrés de liberté
et de fréquences propres est remplacée en utilisant la méthode des éléments finis, par un
modèle ayant des masses concentrées aux nœuds et un nombre fini de degrés de liberté et de
fréquences propres.

Qui, bien que cette équation n'étant plus sous sa forme standard, elle peut facilement être
convertie en cette forme standard. Dans ce cas la matrice de masse globale M G  peut être
concentrée ou consistante.

K e  X  -  2 . X   0

IV.2. CAS DE TESTS

Dans ce chapitre nous procédons à la validation de l’élément T6 (Triangulaire avec 6
nœuds) basé sur le modèle en déplacement en analyse vibratoire libre et forcée. Pour cela
trois différents tests sont illustrés dans cette validation.

Nominalement ces tests représentent un problème d’un solide rectangulaire élastique

en porte à faux (Poutre cantilever) figure (IV.1), la longueur dans la direction x est 4.0 et
suivant y est de 1.0 avec une section transversale de1.0.

1.0

4.0

[Link].1: Géométrie et maillage de la poutre.

56
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

[Link] propres en vibration libre d’un solide rectangulaire

IV.2.1.1. Synoptique du programme de vibration libre d’un solide rectangulaire

Entrée et initialisation

Former la matrice des masses ponctuelles

Pour tous les éléments

Trouver les coordonnées

Faites pour tous les points

d’intégration

Calculer
BTDB
Additionner dans la matrice de
rigidité

Assemblage globale des mass

ponctuelles.

Assemblage globale des rigidités

Réduire aux problèmes standards des valeurs propres

Extraire et imprimer les valeurs propres

[Link].2 : Synoptique du programme de vibration libre d’un solide rectangulaire

57
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

IV.2.1.2. Tests de validation

Notre poutre avec une rigidité à la flexion de 1/12 soit 0.08333. Le coefficient de
poisson est v=0.3 pour le maillage (3×1), et v=0 pour le maillage (5×1), et le module de
Young est pris E=1 kN/m2. Le problème est en déformation plane avec masse concentrée. La
masse par unité de surface est  =1. Les figures suivant représente les maillages de (3×1) et
(5×1), et les différents éléments Q4, Q8 et T6.

[Link].3 : Maillage de (3×1) et l’élément Q4

[Link].4: Maillage de (5×1) et l’élément Q4

[Link].5: Maillage de (3×1) et l’élément Q8

58
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

[Link].6: Maillage de (5×1) et l’élément Q8

[Link].7: Maillage de (3×1) et l’élément T6

[Link].8: Maillage de (5×1) et l’élément T6

Les résultats des valeurs propres sont illustrés dans le tableau (IV.1) pour le cas de
déformation plane. La fréquence fondamentale ω1 et la fréquence axiale ω2 sont calculées
avec les différents éléments Q4, Q8, poutre et T6. Les résultats obtenus par l’élément T6
montrent que les valeurs des fréquences sont en bon accord avec celles obtenues par l’élément
Q8 et avec celles de la solution analytique. La fréquence fondamentale obtenue par l’élément
Q4 est considérablement plus grande que celle de la solution analytique.

59
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

Fréquences Q4 Q8 Poutre T6

Maillage ω1 0.0805 0.064 0.063 0.064

3×1 ω2 0.353 0.413 0.397 0.414

Maillage ω1 0.068 0.060 0.060 0.060

5×1 ω2 0.391 0.391 0.0393 0.392

Solution ω1 1.875 4  EI / AL4 0.063

exacte
[23] Smith ω2  / 2L E /  0.393

[Link].1: Valeurs propres de la poutre cantilever cas de déformation plane

60
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

[Link] forcée d’un solide rectangulaire

IV.2.2.1. Synoptique du programme de vibration forcée d’un solide rectangulaire

Entrée et initialisation

Chargement harmonique pour une seule entrée

Vider les matrices

Former les masses ponctuelles élémentaires
Pour tous les éléments faites

Trouver la géométrie
Pour tous les points de Gauss

Assembler la rigidité des éléments

Trouver la partie réelle et complexe de K

Assembler la rigidité globale complexe

Résoudre les équations

Isoler les composantes réelle et imaginaire de
la solution

Pour tous les pas de temps

Imprimer les réponses

[Link].9 : Synoptique du programme de vibration forcée d’un solide rectangulaire (masse

ponctuelle)

61
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

IV.2.2.2. Tests de validation :

Pour valider l’emploi de l’élément T6 (triangulaire à 6 nœuds) en vibration forcée on à

proposer deux tests qui ont été traité dans la référence [23] Smith 2004.
 Le premier test :

Pour le premier test figure (IV.4) dans l’analyse vibratoire forcée on utilise la méthode de la
réponse complexe décrite dans la référence [23]. La poutre cantilever est soumise à une force
verticale harmonique cos ωt. Le taux d'amortissement γ est de 0,005 soit 5 % appliqué à tous
les modes du système, le module de Young est E=1 kN/m2, le coefficient de Poisson v=0.3,
fréquence forcée ω=0.3 et la masse par unité de surface est  =1 t/m3. Ici le pas du temps est
pris égal à 1/20 de la période (2π/ω) soit disant 1.047. Le problème est en déformation plane.

1.0

4.0
F(t)=cos ωt

[Link].10 : Géométrie et maillage de la poutre cantilever soumise à une vibration forcée

Les résultats illustrés dans le tableau (IV.3) montrent les déplacements à l’extrémité de
la poutre par rapport au pas de temps imposé (time-step) par l’utilisation des deux éléments
Q8 et T6. Nous pouvons voir que le comportement de l’élément T6 est strictement similaire à
celui de l’élément Q8 en analyse vibratoire forcée.

62
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

Temps Déplacements Déplacements

Q8 T6
0 0 0
1.047 27.21 27.40
2.094 36.68 36.87
3.142 42.55 42.72
4.189 44.26 44.40
5.236 41.64 41.73
6.283 34.94 34.98
7.330 24.82 24.80
8.378 12.27 12.19
9.425 -1.480 -1.607
10.47 -15.09 -15.25
11.52 -27.21 -27.40
12.57 -36.68 -36.87
13.61 -42.55 -42.72
14.66 -44.26 -44.40
15.71 -41.64 -41.73
16.76 -34.94 -34.98
17.80 -24.82 -24.80
18.85 -12.27 -12.19
19.90 1.480 1.607
20.94 15.09 15.25

Tableau .IV.2 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au
temps"

Les résultats des déplacements sont illustrés dans le tableau IV.2 et les figures IV. Ces
résultats montrent bien que le comportement de l’élément T6 est similaire à l’élément Q8.

63
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

[Link].11 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au temps
pour l’élément Q8" Test 1

[Link].12 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au temps
pour l’élément T6" Test 1

64
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

[Link].13 : La convergence entre l’élément Q8 et T6

 Le second test :

Pour le second test figure (IV.5) dans l’analyse vibratoire forcée on utilise la méthode
de superposition modale avec la masse concentrée décrite dans la référence [23]. La poutre
cantilever est soumise à une force verticale harmonique cos ωt. Le taux d'amortissement γ est
de 0,005 soit 5 %, le module de Young est E=1kN/ m2, le coefficient de Poisson v=0.3, la
fréquence forcée ω=0.3. Ici le pas du temps est pris égal à 1/20 de la période (2π/ω) soit
disant 1.047. Le problème est en déformation plane.

1.0

4.0
F(t)=cos ωt

[Link].14 : Géométrie et maillage de solide rectangulaire soumis à une vibration forcée

65
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

Temps Force Déplacements Déplacements

Q8 T6
0 1 0 0
1.047 0.9511 27.32 27.32
2.094 0.8090 36.36 36.37
3.142 0.5878 41.84 41.85
4.189 0.3090 43.22 43.23
-14
5.236 -0.5390.10 40.37 40.39
6.283 -0.3090 33.57 33.59
7.330 -0.5878 23.48 23.50
8.378 -0.8090 11.10 11.11
9.425 -0.9511 -2.373 -2.362
10.47 -1.000 -15.45 -15.61
11.52 -0.9511 -27.32 -27.32
12.57 -0.8090 -36.36 -36.37
13.61 -0.5878 -41.84 -41.85
14.66 -0.3090 -43.22 -43.23
-14
15.71 0.5390.10 -40.37 -40.39
16.76 0.3090 -33.57 -33.59
17.80 0.5878 -23.48 -23.50
18.85 0.8090 -11.10 -11.11
19.90 0.9511 2.373 2.373
20.94 1.000 15.45 15.61

Tableau .IV.3 : Vibration forcée d’un solide rectangulaire avec masse localisée "
déplacement par rapport au temps et forces"

66
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

[Link].15 : Sollicitation harmonique "Les forces par rapport au temps"

[Link].16 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au temps
pour l’élément Q8" Test 2

67
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

[Link].17 : Vibration forcée d’un solide rectangulaire " déplacement par rapport au temps
pour l’élément T6" Test 2

Les résultats illustrés dans le tableau (IV.4) montrent les déplacements à l’extrémité de
la poutre par rapport au pas de temps imposé (time-step) et la force perturbatrice cos ωt par
l’utilisation des deux éléments Q8 et T6. Nous pouvons voir clairement que le comportement
de l’élément T6 est strictement similaire à celui de l’élément Q8 en analyse vibratoire forcée.

IV.3. CONCLUSION
Dans ce chapitre nous avons pu faire programmer l’élément triangulaire T6 et valider
cet élément en analyse vibratoire libre et forcée. A travers plusieurs tests numériques réputés
sévères nous avons constaté que l’élément T6 est similaire à l’élément Q8 dans ces analyses.
La validation de l’élément T6 en analyse vibratoire a été le but essentiel de ce chapitre.

68
Chapitre IV : Validation numérique des éléments à champs de déplacement en analyse vibratoire

Conclusion générale
Dans ce travail une investigation dans l’analyse vibratoire des structures membranaires

par la méthode des éléments finis est faite. Les éléments Q4, Q8 et E-Poutre utilisés dans cette

étude sont basés sur l’approche de déplacement. Nous avons pu utiliser l’élément T6 pour

cette analyse vibratoire.

Nous avons testé la performance de l’élément T6 dans les problèmes vibratoires

forcés et libres des structures et faire la comparaison avec les éléments Q4, Q8 et E Poutre.

Ensuite tous les résultats obtenus à l’aide de ces éléments sont comparés aux solutions

exactes (solution analytique)

L’étude montre que:

 Le comportement des éléments à champs de déplacement Q8 et Poutre sont

similaire à l’élément à champs de déplacementT6en analyse vibratoire forcée

et libre.

 Le comportement de l’élément Q4 s’avère un peu mauvais comparativement

avec les éléments Q8 et T6 suivant le calcul de la pulsation ω2 en vibration

libre.

 L’élément Poutre Donne de meilleur résultat.

 En vibration forcée l’élément Q4 est une mauvaise représentation des solides

rectangulaires.

69
 Références bibliographiques

Références bibliographiques
[1] TURNER, M. R. CLOUGH, H. MARTIN, et L. TOPP. , « Stiffness and deflection
analysis of complex Structure », J. Aearo. sci., 23 N°. 9 SEPT (1956)

[2] ARGYRIS, J. et S. KELSEY. , « Energy Theorems and structural analysis» Butterwoth

Scientific publication London, (1960).

[3] [Link] et [Link] et [Link] et [Link]. «Introduction à la

méthode des éléments finis » Editions EYROLLES .61, bd Saint-Germain, 75005
PARIS (1979)

[4] [Link], [Link]. « Introduction à la méthode des éléments finis pour les
ingénieurs » l'Université d’Oran. Office des publications universitaires, 06 -88.

[5] MARTIN LEVESQUE. « Cours d’éléments Finis en Mécanique des Solides» Ecole
Polytechnique de Montréal

[6] ANTOINE LEGAY. «Calcul des structures par éléments finis» CNAM-Paris, Édition
(2012)

[7] SMITH I.M, GRIFFITH D.V. « Programming the Finite Element Method», Chapter 2,
second Edition. John Wiley & Sons, UK, (1988).

[8] A.Y.T. LEUNG, ZHU BIN, JIANJIN ZHENG, HAO YANG « Analytic trapezoidal
Fourier p-element for vibrating plane problems »Journal of Sound and Vibration. City
University of Hong Kong, Tatchee Avenue, HKSAR, China, (2003).

[9] [Link], [Link], [Link] « Quadrilateral and triangular plane elements with
rotational degrees of freedom based on the hybrid Trefftz method »Korea Advanced
Institute of Science and Technology, Korea,(2006).

[10] [Link] .Article « Free vibration analysis of arbitrarily shaped membranes using

70
 Références bibliographiques

the trigonometric p-version of the finite-element method » Faculty of Engineering,

University of Tlemcen, B.P. 230, Tlemcen 13000, Algeria(2006)

[11] KUTLU DARILMAZ et DE NAHIT KUMBASAR. « An 8-node assumed stress hybrid

element for analysis of shells » Department of Civil Engineering, Istanbul Technical
University, 34469 Maslak, Istanbul, Turkey, 2006.

[12] N. NGUYEN-THANHA, TIMON RABCZUKA, H. NGUYEN-XUANBE,

STÉPHANE PA BORDASC, D. « An alternative alpha finite element method (AαFEM)
for free and forced structural vibration using triangular meshes »Journal of
Computational and Applied Mathematics 2009

[13] L. CHAMPANEY. «Vibration de systèmes continus» Notes du cours du Dynamique des

Constructions. Notes du cours du Dynamique des Constructions, Ecole de ENSMP,
France.

[14] [Link], [Link] et [Link]. « Aide-mémoire Mécanique des

Structures Résistance des Matériaux » DUNOD, Paris, Édition (2008)

[15] SINGIRESU S. RAO« Mechanical Vibrations » Second Edition .PURDUE University.

(2000) S. GRAHAM KELLY. « Mechanical Vibrations theory and applications» the
university of [Link]. Edition 2012

[16] ALAIN PECKER. «dynamique des structures et des ouvrages » Ecole nationale des
ponts et chaussées FRANCE. Édition (2009)

[17] S. GRAHAM KELLY. « Mechanical Vibrations theory and applications» the university
of [Link]. Edition 2012

[18] JEAN-CLAUDE PASCAL. « Vibrations et acoustique 2 » Ecole Nationale Supérieure

d'Ingénieurs du mans - Université du Maine. (2009)

[19] J. BESNIER, R. TARBADAR. «Travaux Pratiques en Comportement Dynamique »

71
 Références bibliographiques

Génie Mécanique. ICAM2 .Nantes.

[20] MAURICE PETYT. « Introduction to finite element vibration analysis » CAMBRIDGE

university PRESS (1998)

[21] JL BATOZ, G DHATT. « Modélisation des structures par éléments finis, Poutre et
Plaque » volume 2, Les Presses de l'Université LAVAL Sainte-Foy, Québec. (1990).

[22] [Link], [Link]. «Une présentation de la méthode des éléments

finis» ,deuxième édition, PARIS (1984).

[23] SMITH I.M, GRIFFITH D.V. « Programming the Finite Element Method », Chapter 9,
John Wiley & Sons, UK, (2004).

72