HAM7246

‫الجمـهوريــة الجزائريــة الديمقراطيــة الشعبيــة‬
REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE
Ministère de l’enseignement supérieur ‫وزارة التعليــم العالــي و البحــث العلمـي‬

et de la recherche scientifique ‫جامعـــــــةاألخوة منتوري قسنطينـــــــــة‬
Université Frères Mentouri Constantine1 ‫كلـيــة علوم التكنولوجية‬
Faculté des Sciences de la Technologie ‫قسم الهندسة الميكانيكية‬
Département Génie Mécanique
N0 d’ordre :87/DS/2018
06/GM/2018
THESE.
En vue de l’obtention du
Doctorat en Sciences
Spécialité : Génie Mécanique
Option : Construction mécanique
HAMRITFareh
Thème
ANALYSE DYNAMIQUE DES STRUCTURES MECANIQUES

PAR LA METHODE DES ELEMENTS FINIS
Soutenue publiquementle : ……………… /…………………/………………….........
Devant le Jury composé de:

Président : Mr. BOUCHOUCHA Ali Professeur à l’Université Frères Mentouri Constantine
Rapporteur :Mr. NECIB Brahim Professeur à l’Université Frères Mentouri Constantine
Examinateurs :
Mr. MERABET Abderrazak MCA à l’Université Frères Mentouri Constantine
Mr. KHEMISSA Mohamed Professeur à l’Université M’sila
Mr. ZAOUI Moussa MCA à l’Université M’sila
Membre invité :Mr. DRISS Zied Professeur à l’Université de Sfax (Tunisie)
Année Universitaire 2017-2018

Résume
Résumé
Les structures mécaniques discrètes ou continues ontune grande importance dans le domaine
d’application de l’industrie moderne contemporaine. Cependant, durant leur fonctionnement ces
structures sont souvent soumisesà des sollicitations externes considérables ou à desfortes amplitudes
de vibrations qui peuvent leur causer de larges déformations et des contraintes internes pouvant leur
provoquer des fissurations internes oumême leur destruction totale.Afin d’éviter ces types de problèmes,
le concept de l'analyse statique et dynamique de ces structures est recommandé, et en raison de la complexité
de leur forme et de leur dimension, la méthode des éléments finis est la plus utilisée.Cette dernière est
actuellement reconnue comme étant une technique très puissante pour l’analyse statique et dynamique
des structures discrètes ou continues de forme compliquée appliquée dans le domaine de la mécanique,
l’aéronautique,le génie civil,le maritime ou la robotique. En conséquence, le calcul et le
dimensionnement de ces systèmes mécaniques par la méthode des éléments finis joue actuellement
un rôle important au service de l’industrie pour d’éventuels dimensionnement et prédiction de leur
durée de vie.Notre travail consiste à l’analyse statique et dynamique des systèmes mécaniques discrets
et continus bidimensionnels utilisant la méthode des éléments finis en se basant sur les éléments
principaux barres, poutres et plaques , sous l’effet des excitations extérieures avec différentes
conditions aux limites. Les structures discrètes considérées sont bidimensionnelles en charpente
métallique interconnectées aux nœuds par soudage, rivetage ou boulonnés sous différentes conditions
aux limites. Leurs éléments sont modélisés comparativement par des éléments barres et des éléments
poutres, tandis que pour les structures continues, les éléments sont des plaques minces rectangulaires
avec différentes conditions aux limites.Les forces d'excitation sont basées sur des forces périodiques,
aléatoires ou impulsives et une solution numérique par développement d’un programme pour décrire
le comportement de ces structures est réalisée. Les matrices de masse et de rigidité de toutes les
structures sont déterminées respectivement par assemblage des éléments barres, poutre et plaque en
se basant sur l’énergie cinétique et de déformation pour chaque élément. Les déplacements, les
réactions aux nœuds et les forces axiales dans tous les éléments ainsi que les contraintes transversales
et les modes propres des structures sous différentes conditions limites ont été aussi calculées et de
bons résultats ont été obtenus comparés à ceux obtenus utilisant d’autres logiciels déjà existants. En
fait, l’analyse utilisant la méthode éléments finis permettra au bon dimensionnement et au design
adéquat des structures mécaniques industrielles complexes en fonction de différentes conditions aux
limites, de leur chargement interne et de leur niveau vibratoire.
Mots clés : Structures mécaniques, rigidité, déformation, contrainte, modes de vibration, éléments finis,
élément barre, élément poutre, élément plaque, assemblage, énergie de déformation, forces internes.
Résume
Abstract:
The mechanicaldiscrete or continuumstructures are actually of great importance in the application

field of contemporary modern industry. However, during their life time these structures are often
subjected to considerable external stresses or to high amplitudes of vibrations which can cause them
large deformations and internal stresses which can cause them internal cracking or even their total
destruction. In order to avoid these types of problems, the concept of static and dynamic analysis of
these structures is recommended, and due to the complexity of their shape and size, the finite element
method is the most used. The latter is currently recognized as a very powerful technique for the static
and dynamic analysis of discrete or continuous structures of complicated form applied in the field of
mechanics, aeronautics, civil engineering, maritime or robotics. Consequently, the calculation and
dimensioning of these mechanical systems by the finite element method plays an important role at the
service of the industry for possible sizing and prediction of their lifetime. Our work consists of static
and dynamic analysis of two-dimensional discrete and continuous mechanical systems using the
finite element method based on the main elements of bars, beams and plates, under the effect of
external excitations with different boundary conditions. The discrete structures considered are two-
dimensional in metallic framework interconnected to the nodes by welding, riveting or bolted under
various boundary conditions. Their elements are modeled comparatively by bar elements and beam
elements, while for continuous structures the elements are rectangular thin plates with different
boundary conditions. The excitation forces are based on periodic, random or impulsive forces and a
numerical solution by development of a program to describe the behavior of these structures is
realized. The mass and stiffness matrices of all the structures are determined respectively by
assembling the bars, beam and plate elements based on the kinetic and deformation energy for each
element. The displacements, the node reactions and the axial forces in all the elements as well as the
transverse stresses and the eigenvalues of the structures under different boundary conditions were
also calculated and good results were obtained compared to those obtained using other software
Already existing. In fact, analysis using the finite element method will allow the proper dimensioning
and design of complex industrial mechanical structures according to different boundary conditions,
their internal loading and their vibratory level.
Keywords:Mechanical structures, stiffness, deformation, stress, vibration modes, finite elements, bar
element, beam element, plate element, assembly, deformation energy, internal forces.
‫‪Résume‬‬
‫خالصة‪:‬‬
‫الهياكل الميكانيكية المنفصلة أوالمتصلةذات أهمية كبيرة في مجال تطبيق الصناعة الحديثة المعاصرة‪.‬ومع ذلك‪ ،‬غالبا ما تتعرض هذه‬
‫الهياكل إلى الضغوط الخارجية الكبيرة أو لسعة عالية من االهتزازات التي يمكن أن تسبب لهم تشوهات كبيرة والضغوط الداخلية‬
‫التي يمكن أن تسبب تكسير الداخلي أو حتى تدميرها الكلي‪ .‬من أجل تجنب هذه األنواع من المشاكل‪ ،‬فمن المستحسن مفهوم التحليل‬
‫والدينامكية من هذه الهياكل‪ ،‬ويرجع ذلك إلى تعقيد الشكل والحجم وطريقة العناصر المحدودة هو األكثر استخداما‪ .‬واعترفت هذه‬
‫األخيرة حاليا على تقنية قوية جدا لتحليل والدينامكية هياكل منفصلة أو متصلة من شكل معقد المطبقة في مجال الميكانيكا‪ ،‬علم‬
‫الطيران والفضاء‪ ،‬الهندسة المدنية‪ ،‬البحرية والروبوتكا‪ .‬ونتيجة لذلك‪ ،‬وحساب وأبعاد هذه األنظمة الميكانيكية باستخدام طريقة‬
‫العناصر المحدودة يلعب دورا هاما في خدمة هذه الصناعة إلمكانية التحجيم والتنبؤ بها حياتهم‪ .‬يتكون عملنا التحليل والدينامكية‬
‫األنظمة الميكانيكية المنفصلة والمتصلة ثنائية األبعاد باستخدام طريقة العناصر المحددة استنادا إلى العناصر الرئيسية من الحانات‪،‬‬
‫والحزم وألواح‪ ،‬تحت تأثير اإلثارات الخارجية مع الشروط الحدية المختلفة‪ .‬الهياكل المنفصلة تعتبر هي ثنائية األبعاد في إطار‬
‫معدني مترابطة إلى العقد عن طريق اللحام‪ ،‬التثبيت أو انسحب تحت شروط الحدود المختلفة‪ .‬وعلى غرار عناصرها نسبيا من قبل‬
‫عناصر شريط وعناصر شعاع‪ ،‬أما بالنسبة للهياكل مستمرة العناصر هي لوحات رقيقة مستطيلة مع الشروط الحدية المختلفة‪ .‬وتقوم‬
‫قوات اإلثارة على قوات الدورية‪ ،‬عشوائية أو التسرع والحل العددي بواسطة وضع برنامج لوصف يتحقق سلوك هذه الهياكل‪ .‬يتم‬
‫تحديد المصفوفات الكتلة وصالبة من جميع الهياكل على التوالي عن طريق تجميع القضبان‪ ،‬عناصر شعاع ولوحة على أساس الطاقة‬
‫الحركية وتشوه لكل عنصر‪ .‬كما تم حساب عمليات النزوح‪ ،‬ردود فعل عقدة والقوى المحورية في جميع العناصر فضال عن الضغوط‬
‫عرضية القيم الذاتيةوالهياكل تحت شروط الحدود المختلفة وتم الحصول على نتائج جيدة مقارنة بتلك التي تم الحصول عليها باستخدام‬
‫برامج أخرى إذا الموجودة‪ .‬في الواقع‪ ،‬فإن التحليل باستخدام طريقة العناصر المحدودة يسمح لألبعاد السليم وتصميم الهياكل‬
‫الميكانيكية الصناعية المعقدة وفقا لشروط الحدود المختلفة‪ ،‬تحميلها الداخلي ومستوى اهتزازي‪.‬‬
‫الكلمات المفتاحية‪ :‬الهياكل الميكانيكية‪ ،‬وصالبة‪ ،‬وتشوه‪ ،‬واإلجهاد‪ ،‬وسائط االهتزاز‪ ،‬والعناصر المحدودة‪ ،‬عنصر شريط‪ ،‬عنصر شعاع‪،‬‬
‫عنصر لوحة‪ ،‬والتجمع‪ ،‬والطاقة تشوه‪ ،‬القوى الداخلية‪.‬‬
Remerciements
Je tiens à remercier Monsieur NECIB Brahim Professeur à l’université de

Constantine1,qui n’a jamais compté les heures afin que mon travail puisse avancer le
plus rapidement possible, également avec un grand dévouement, ils m’a consacré tout
leur temps à suivre de prés l’évolution de cette thèse à orienter les différentes étapes et
à palier toutes les difficultés auxquelles j’avais à faire face. Merci pour
l’enrichissement scientifique que vous m’avez apporté et pour la confiance que vous
m’avez accordée.
Je remercie Monsieur BOUCHOUCHA Ali, Professeur à l’université des frères

Mentouri Constantine, pour m’avoir fait l’honneur de présider le jury de cette thèse.
Monsieur MERABET Abderrazak, Professeur à l’université frères Mentouri

Constantine, et Monsieur KHEMISSA Mohamed Professeur à l’université de M’sila,
et Monsieur ZAOUI Moussa, Maître de conférences à l’université de M’sila, pour
l’intérêt qu’ils ont témoigné pour mon travail en acceptant de rapporter ce manuscrit.
J’exprime aussi mes vives reconnaissances à Monsieur ZIED DRISS, Professeur à

l’Ecole Nationale d’Ingénieurs de Sfax (Tunis) , pour l’intérêt qu’il a porté à ce travail
en acceptant de faire partie du jury de cette thèse.
Enfin, mes sincères remerciements s’adressent à tous les enseignants du l’université

Mohamed BOUDIAF (M’sila), et plus particulièrement à tous les membres du
laboratoire des Systèmes Electromécaniques (LASEM), et à tous ceux qui m’ont aidé
de près ou de loin à l’élaboration de ce travail.
Dédicaces
Je dédie ce modeste travail :
A mes très chers parents
A ma très chère femme Fatima;
A mes chers enfants; Sohaib, mouncif, tamim
A mes très chers frères et mes sœurs
A mes très amis B.Bakri, A.Righi, M .Meddah
A ma grande famille;
Fareh HAMRIT
Table des matières
1
Analyse dynamique des structures mécaniques par MEF
Table des matières
Table des matières
Résume ………………………………………………………………………………………… I
Remerciement………………………………………………………………………………….. II
Dédicace………………………………………………………………………………………… III
Table des matières …………………………………………………………………………… 1

Index des notations…………………………………………………………………………… 6
Liste des Figures et tableaux ………………………………………………………………… 9
Introduction ……………………………………………………………………………… 12
Chapitre 1 : Introduction générale et Etat de l’art
1 .1 Introduction ……………………………………………………………………………… 18
1.2 La méthode des éléments finis (MEF) …………………………………………………... 18
1.3 Historique…………………………………………………………………………………. 20
1.4 Problématique…………………………………………………………………………….. 22
1.5 Applications et état de l’art……………………………………………………………… 23
1.5.1 la construction mécanique…………………………………………………….. 24
1.5.2 Aéronautique et l’aérospatiale………………………………………………… 26
1.5.3 le Génie Civil…………………………………………………………………… 27
1.5.4 la Robotique et la Biomécanique……………………………………………… 28
Chapitre 2 : Techniques de modélisation utilisées en éléments finis
2.1 Introduction………………………………………………………………………………. 31
2.2 Principes de formulation de la méthode des éléments finis…………………………….. 31
2.2.1 Principe de la méthode des forces-déplacements……………………………………… 31
2.2.2 Principe de la méthode variationnelle…………………………………………………. 32

2
Table des matières
2.2.Champd’application………………………………………………………………………. 33
2.3 Techniques de modélisations utilisées en éléments finis………………………………… 34
2.3.1 Méthode de maillage (nœuds et éléments) …………………………………………….. 34
2.3.2Etapes de modélisation par la MEF…………………………………………………… 35
2.4 Classification des éléments de la MEF…………………………………………………… 37
2.4.1 Introduction à l’état plan de contrainte et de déformation…………………………... 37
2.4.2 Connaissances requises pour la réalisation des logiciels MEF……………………..... 38

2.4.3 Connaissances requises pour la réalisation des logiciels MEF……………………...... 38
2.5Discrétisation des structures par éléments finis……………………................................. 39
2.5.1 Types d’éléments…………………………………………………………………………. 39
2.5.2 Assemblage………………………………………………………………………………. 41
2.5.3 Elimination directe de la bande de matrice……………………………………………….. 43
2.6 Types de numérotation…………………………………………………………………….. 44
2.6.1 Numérotation verticale…………………………………………………………………… 44
2.6.2 Numérotation horizontale………………………………………………………………… 44
2.7 Calcul d’erreurs dans l’analyse par MEF…………………………………………………… 45
2.7.1 Erreurs usuelles…………………………………………………………………………… 45
2.7.2 Erreurs détectées (Vérification du modèle) ……………………………………………….. 46
2.8 Critique des résultats obtenus par MEF……………………………………………………... 47
Chapitre 3: Détermination de matrice de rigidité des éléments des structures mécaniques
3.1 Introduction………………………………………………………………………………. 51
3.2 Détermination de la matrice de rigidité de l’élément de barre…………………………. 51
3.2.1Matrice de rigidité en coordonnées locales…………………………………………….. 51
3.2.2Matrice de rigidité en coordonnées globales…………………………………………… 54
3.3Détermination de la matrice de rigidité pour un élément poutre……………………….. 58
3.3.1 Matrice de rigidité en coordonnées locales……………………………………………. 58
3.3.2Equation de rigidité de l’élément poutre en coordonnées locales………………………. 60

3
Table des matières
3.3.3Matrice de rigidité en coordonnées globales…………………………………………… 62

3.3.4Fonction de déplacement………………………………………………………………… 64
3.3.5Transformation des coordonnées locales aux coordonnées globales………………….. 67
3.3.6 Détermination de la force axiale………………………………………………………… 68
3.4 Détermination de la matrice de rigidité pour un élément plaque……………………… 68
Chapitre 4 : Mise en équation de l’analyse dynamique des structures mécaniques

4.1 Introduction………………………………………………………………………………… 76
4.2 Analyse dynamique d’un élément barre………………………………………………….. 77
4.3 Détermination de la matrice masse d’un élément barre…………………………………. 80
4.4 Vibration d’un élément poutre…………………………………………………………….. 81
4.5 Détermination des modes propres pour un élément poutre……………………………… 84
4.6 vibration d’un élément plaque……………………………………………………………... 85
4.7 Détermination des modes propres d’un élément poutre et matrice masse…………….... 88
4.8 Détermination de la matrice masse……………................................................................... 88
Chapitre 5: Analyse dynamique des structures mécaniques discrètes utilisant la MEF

5.1 Introduction………………………………………………………………………………… 92
5.2 Structures en éléments barre et poutre…………………………………………………… 92
5.3 Caractéristique de la structure en élément poutre………………………………………. 93
5.4 Les déplacements nodaux suivant l’axe ox……………………………………………….. 93
5.5 Les déplacements nodaux suivant l’axe oy………………………………………………. 94
5.6. Les déplacements angulaires……………………………………………………………… 96
5.7. Les réactions……………………………………………………………………………….. 96
5.8. Les forces axiales…………………………………………………………………………… 97
5.9. Les modes propres………………………………………………………………………… 98
Chapitre 6: Analyse dynamique de structures continues planes utilisant la MEF

6.1 Introduction……………………………………………………………………………….. 102
6.2 Structures en éléments plaques…………………………………………………………… 102
6.3 Premier Cas : encastrement d'un côté et libre d'autre coté…………………………….. 106

4
Table des matières
6.3.1 les modes propres de la structure………………………………………………………… 106

6.3.2 Déplacements nodaux……………………………………………………………………. 108
6.3.3 Réactions aux appuis……………………………………………………………………... 110
6.4 Deuxième Cas : encastrement de deux cotés……………………………………………… 112
6.4.1 les modes propres de la structure……………………………………………………….. 112
6.4.2 déplacements nodaux……………………………………………………………………. 112
6.4.3 Réactions aux appuis ……………………………………………………………………. 117
6.5 Les déplacements nodaux suivant l’axe ox………………………………………………... 120
6.6 Les déplacements nodaux suivant l’axe oy………………………………………............... 120
6.7 Les contraintes normales dans les éléments………………………………………............. 121
6.8 Les contraintes tangentielles dans les éléments………………………………………........ 122
6.9 Les modes propres suivant l’axe ox………………………………………............................ 122
6.10 Les modes propres suivant l’axe oy……………………………………….......................... 123
6.11 Conclusion………………………………………................................................................... 124
Conclusion générale ….................................…………………………………............................ 125
Perspectives……………………………………….................................…….............................. 128
Références bibliographiques………………………………………............................................. 130

5
Index des notations
Index des notations
 ij : Tenseur de contrainte
 ij : Tenseur de déformation
 : Tenseur des contraintes N / m  2
ε: Les déformations internes de la structure
 : Dérivée partielle
E : Module de «Young » N / m 2 
 : Coefficient de poisson
 : Densité du matériau
𝐸𝑐 : Energie cinétique
𝐸𝑝 : Energie potentielle
𝑋̅, 𝑌̅ : Coordonnées locales
𝑋, 𝑌 : Coordonnées globales
ds : Surface différentielle m  2
dv : Volume différentiel m 
3

n : Normale
u1 et u 2 : Déplacements nodaux m
f i x  : Fonctions de forme

f i x  : Dérivée première de la fonction de forme

6
Index des notations
U : Énergie de déformation totale N.m
Ti : Vecteur de contrainte
F : Vecteur de force
F (m) : Force de volume
f ( p) : Force de surface
K  : Matrice de rigidité
K ij : Coefficients de la matrice de rigidité
M  : Matrice de masse
mij : Coefficients de la matrice de masse
W : Le travail des forces volumiques et surfaciques
 T : Transposé d’un vecteur
tr  : Trace 
X , Y , Z : Vecteurs unitaires
a : Dimension de la plaque suivant l’axe x m
b : Dimension de la plaque suivant l’axe y m
e:Epaisseur de la plaque [m]

U
Fi  : Forces généralisées N 
qi
q i : Déplacements généralisées m

7
Index des notations

f i (x) : Dérivées secondes des fonctions de forme
Fxi , Fy j : Forces aux nœuds pour chaque élément N 
𝜔 : La fréquence [rad/s]
{x} : Vecteur des déplacements aux nœuds
{𝑥̈ } : Vecteur des accélérations aux nœuds
{F(t)} : Vecteur des forces aux nœuds
𝜆 : cos 𝜃
: sin 𝜃
𝜕𝑣
𝜃= : Rotation de la section (rad)
𝜕𝑥
𝐸𝐴 : Rigidité longitudinale (N)
𝐸𝐼: Rigidité de flexion (𝑁. 𝑚2 )
𝑣1
𝜃
{𝑣1 }: Tenseur de déplacement d’un élément poutre
2
𝜃2
𝐹1
𝑀
{ 1 } : Tenseur forces d’un élément poutre.
𝐹
𝑀2

8
Listes des figures et tableaux
Listes des figures
Figure 1.1 : Discrétisation d’une structure en nœuds et éléments

Figure 1.2 : Principales étapes pour une analyse par la Méthode des Eléments Finis
Figure 1.3 : Représentation des nœuds et les éléments
Figure 1.4 : Forces et déplacements de quatre nœuds.
Figure 1.5 : Assemblage des nœuds.
Figure 1.6 : Structures en élément barre ou poutre.
Figure 1.7 : Structures en élément poutre et plaque.
Figure 1.8 : Structure en élément coque.
Figure 1.9 : Structure en élément poutre dans le domaine génie civil
Figure 1.10 : Structure en élément poutre dans le domaine hydraulique.
Figure 1.11 : Structureen élément poutre d’un bras de robot.
Figure 1.12 : Structure en élément poutre d’un robot
Figure 1.13 : Structure en élément poutre d’un châssis d’automobile
Figure 2.1 : Système local et global de coordonnées
Figure 2.2: Forces nodales dans le système local et global d’axes de coordonnées
Figure 2.3: Système local et global de coordonnées
Figure 2.4: Elément rectangulaire à 4 nœuds
Figure 2.5: a) Simulation de la flexion. b) Mode de déformation correcte en flexion
pure. c) Mode de déformation de l’élément bilinéaire quadratique
Figure 2.6: Fonction de forme f2 de l’élément bilinéaire quadrilatéral
Figure 2.7: Elément quadratique quadrilatéral
Figure 2.8: Fonctions de forme
Figure 3.1: Les forces axiales pour un élément barre
Figure 3.2: Elément poutre en flexion
Figure 3.3: Système à un degré de liberté
Figure 3.4: Elément rectangulaire de huit degrés de liberté.

9
Figure 4.1 : Structures en élément poutre avec des différentes conditions aux limites
Figure 4.2 : Comparaison des déplacements suivant l’axe ox
Figure 4.3 : Comparaison des déplacements suivant l’axe oy
Figure 4.4 : Comparaison des déplacements angulaire
Figure 4.5 : Comparaison des déplacements angulaire
Figure 4.6 : Comparaison des forces axiales
Figure 4.7 : Les Modes propres de trois structures
Figure 5.1 : Elément plaque rectangulaire encastré-libre
Figure 5.2 : Elément plaque rectangulaire encastré-encastré
Figure 5.3 : Mode propre 1 d’une plaque encastrée-libre
Figure 5.8 : Première mode de la structure A, (a)Ansys, (b) Abaqus.
Figure 5.9 : Deuxième mode de la structure A, (a)Ansys, (b) Abaqus.
Figure 5.10:Troisième mode de la structure A, (a)Ansys, (b) Abaqus.
Figure 5.11: Quatrième mode de la structure A, (a)Ansys, (b) Abaqus.
Figure 5.12: Cinquième mode de la structure A, (a)Ansys, (b) Abaqus.
Figure 5.13: Les déplacements nodaux de Cinque premiers Modes (premier cas)
Figure 5.14: Les réactions aux appuis de Cinque premiers Modes (premier cas)
Figure 5.15: Première mode de la structure B, (a)Ansys, (b) Abaqus.
Figure 5.16: Deuxième mode de la structure B, (a)Ansys, (b) Abaqus.
Figure 5.17: Troisième mode de la structure B, (a)Ansys, (b) Abaqus.
Figure 5.18: Quatrième mode de la structure B, (a)Ansys, (b) Abaqus.
Figure 5.19: Cinquième mode de la structure B, (a)Ansys, (b) Abaqus.
Figure 5.20: Les déplacements nodaux de Cinque premiers Modes (deuxième cas)
Figure 5.21: Les réactions aux appuis de Cinque premiers Modes (deuxième cas)

10
Figure 5.22: Comparaison des déplacements nodaux suivant l’axe ox

Figure 5.23: Comparaison des déplacements nodaux suivant l’axe oy
Figure 5.24: Comparaison des contraintes normales
Figure 5.25: Comparaison des contraintes tangentielles
Figure 5.26: Comparaison des modes propres suivant l’axe ox
Figure 5.27: Comparaison des modes propres suivant l’axe oy
Liste des tableaux

Tableau 2.1 Type des champs d’application de la MEF
Tableau 2.2 Types d’éléments finis
Tableau 2.3 Bande d’une matrice à deux dimensions
Tableau 2.4 Vérification des erreurs
Tableau.5.1 Fréquences des modes propres des structures (Hz)

11
Introduction

12
Introduction
INTRODUCTION
Les structures mécaniques discrètes ou continu ont une grande importance dans
le domaine de technologie tel que les structures discrètes composées d’éléments barres
ou d’éléments poutres assemblées par soudures ou rivetage en des points
appelés«nœuds»,etsoumises àdes forces extérieures, sous l’effet de ces forces, la
structure peut se déformer et des contraintes internes dans chaque élément peuvent se
manifester, ces structures sont caractérisées par un nombre finis de paramètres inconnus
tel que les déplacements et les réactions des nœuds et les forces axiales dans les
éléments. Ces types des structures supportent des vibrations propres et forcées.
L’analyse statique et dynamique de ces structures est une discipline très importante, qui
s’est développée pour répondre à des besoins de construction, initialement dans le
domaine du Génie mécanique, Génie-Civil et aérodynamique. Elle repose sur
l’utilisation de modèles simplifiés, qui vont permettre l’analyse des structures de façon
rapide. Ces modèles exploitent une caractéristique essentielle des structures qui sont des
solides déformables tridimensionnel.
– les structures minces dont une dimension (l’épaisseur) est très petite devant les deux
autres, et qui sont appelées plaques ou coques
– les structures élancées dont une dimension (la longueur) est très grande devant les
deux autres, et qui sont appelées poutre ou arc
Ces structures constituent aujourd’hui l’immense majorité des structures industrielles, et
dans tous les domaines : aéronautique, automobile, construction ferroviaire et navale,
génie civil, etc. Ces deux types de structures sont aussi souvent combinés entre elles,
par exemple en renforçant des plaques par des poutres.
Bien que pour l’étude du comportement dynamique des systèmes continus, la
forme faible discrète s’appuyant sur la méthode de Galerkin soit une approche à
première vue plus puissante que la formulation forte basé sur l’écriture directe des
13
Introduction
équations différentielles du mouvement, elle laisse rapidement entrevoir ses limites en

présence de géométries complexes, son inconvénient principal résidant dans le fait que
les fonctions de forme qui caractérisent les champs des déplacements réels et virtuels
approchés sont définies globalement sur l’ensemble du domaine étudié.
La méthode des éléments finis contourne cette difficulté en subdivisant le corps
déformable de forme compliquée en un assemblage de sous-domaines ou éléments
simples et en construisant les solutions approchées sous la forme d’une combinaison de
fonctions de forme à support compact. Dérivée en tant que cas particulier d’application
de la méthode de Galerkin, elle constitue un procédé numérique général et systématique
pour la discrétisation des structures et connait depuis plusieurs années un succès
considérable, notamment grâce au développement de nombreux codes d’élément finis,
dont le domaine d’application s’étend largement au l’élastodynamique linéaire [1] .
La méthode des éléments finis est actuellement reconnue comme étant une
technique très puissante pour l’analyse statique et dynamique des structures discrètes
ou continues, dans le domaine d’engineering. Elle est utilisée aujourd’hui dans de
nombreux secteurs de l’industrie tel que, la mécanique, le génie civil, l’aérospatial et la
robotique.
Notre travail est consacré à l’utilisation de cette méthode pour l’analyse statique et
dynamique des structures en élément barre, poutre et plaque sous l’effet des excitations
extérieures avec des différentes conditions aux limites. La compréhension de cette
méthode donne nécessité à l’élaboration de certaines connaissances scientifiques telles
que , la théorie d’élasticité, la mécanique des milieux continus, la résistance des
matériaux, la dynamique des structures, de plus les notions de base de mathématique
appliquée. Si la structure présente un comportement trop complexe ou un système
continu défini par un nombre infini de paramètres, il devient très difficile d’analyser ou
de trouver la solution analytiquement. Cependant, la méthode des éléments finis pousse
la possibilité de trouver la solution la plus parfaite possible, en remplaçant le système
14
Introduction
continu par un système discret, caractérisé par un nombre fini de paramètres.

Les principales étapes de la méthode des éléments finis sont :
1- La géométrie est décomposée en domaines de forme géométrique simple. Les
éléments sont reliés entre eux en des points appelés (nœuds). Elémentairement il peut
être un élément barre ou poutre.
2- Le champ de déplacements dans chaque élément est défini en fonction des
déplacements des nœuds de l’élément ainsi que l’énergie de déformation et sa matrice
de rigidité.
3- La matrice de rigidité globale est construite à partir des matrices de rigidité
élémentaires.
4- Après la mise en place des conditions aux limites et des charges, on calcul les
déplacements inconnus dans chaque nœud, les déformations et les contraintes dans
chaque élément.
Cette thèse est constituée de six chapitres, une introduction générale et une
conclusion. Le premier chapitre est réservé pour une introduction générale et l’état de
l’art , dont les domaines d’application de ces types des structures compliquées dans
l’industrie, on trouve aussi dans ce chapitre une représentation des différentes
structures mécaniques de plus en plus complexes .
Le deuxième chapitre, une description détaillé des techniques utilisant les
éléments finis pour la modélisation des différentes structures mécaniques, tel que
Cette méthode a été appliquée pour l’analyse statique et dynamique des structures
mécanique, on a cité quelques étapes pour l’analyse par cette méthode, ainsi la
discrétisation des structures par éléments finis, une description détaillé pour
l’assemblage par la méthode pour déterminer la matrice globale de toute la structure,
même pour les numérotations des nœuds pour minimisée la largueur de la bande de la
matrice.
15
Introduction
Le troisième chapitre est consacré à les différentes étapes pour la détermination

des matrices de rigidité pour chaque type de structure, en va rassembler les matrices
élémentaires de chaque élément pour obtenir la matrice globale, dans ce chapitre une
présentation pour l’application de la méthode des éléments finis sur des différentes
structures on utilisant l’élément barre, poutre et plaque, ce chapitre aussi réservé pour la
détermination des déplacements nodaux, les forces axiales dans les éléments, et les
réactions aux appuis.
Dans le quatrième chapitre la mise en équation de l’analyse dynamique des
structures par la méthode des éléments finis, pour l’analyse dynamique des structures
mécanique en élément barre, poutre et plaque. Tel que les vibrations des systèmes
mécaniques basées sur la détermination de la matrice masse de la structure, pour
trouver les modes propres de chaque structure, avec l’utilisation de théorème de
Lagrange qui liée à l’énergie totale de la déformation.
Le cinquième chapitre consiste à l’application de la méthode des éléments finis

sur structures mécaniques discrètes (élément barre et élément poutre), pour l’analyse
statique et dynamique de ces structures, avec des différentes conditions aux limites.
Pour voir l’effet de la condition sur les déplacements nodaux ou les forces axiales et les
réactions aux appuis, d’un côté, et d’autre coté sur les modes propres des structures. On
a élaborer dance ce chapitre un programme de calcul des déplacements, des réactions
dans chaque nœud et des forces axiales dans chaque élément pour l’analyse statique et
les modes propres pour le cas dynamique en comparant les résultats de ces calculs avec
les résultats du logiciel « SAP 2000 », et validant par le logiciel Abaqus et Ansys.
Le sixième chapitre réservé pour l’application de la méthode des éléments finis sur les
structures mécaniques continus (élément plaque), tel que la détermination des
déplacements nodaux et les réactions aux appuis pour l’état statique, et les modes
propres de la structure avec des différentes conditions aux limites.
16
Chapitre 1 :Introduction généraleet Etat de l’art
Chapitre 1 :
Introduction générale
etEtat de l’art
17
1.1 Introduction
La modélisation numérique est devenue incontournable dans le monde de la

conceptionindustrielle et de la recherche scientifique. Les problèmes à résoudre
étant de plus en pluscomplexes, il est nécessaire d’adopter une approche adaptée
au problème traité. En effet, la méthode des éléments finis (MEF) – ou (FEM
« en anglais »)est actuellement la méthode la plus connue et la plus utilisée dans
le domaine de la simulation numérique.
Dans les domaines d’application pour l’analyse statique ou dynamique des
structures mécaniques, la méthode deséléments finis (MEF) est la méthode de
référence mais n’est pas optimale pour les structures très complexes car elle
nécessite de mailler un grand volume d’air pour obtenir une bonneprécision.
Le principe de cette méthode des éléments finis consiste à restreindre la
détermination du champ considérer par un nombre fini de points du milieu,
appelés « nœuds » et l’assemblage de tous ces points qui représente le maillage
de tout le corps (Fig. I.1).
1.2 La méthode des éléments finis (MEF)
La méthode des éléments finis (MEF) est une méthode d’approximation

numérique de solutions de problèmes aux limites statiques ou dynamiques des
milieux continus ou fluides en se basant sur les équations différentielles ou aux
dérivées partielles avec des conditions aux limites déterminées.C’est une
méthode de discrétisation destinée à chercher une solution approchée à un
problème physique mais sous une forme variationnelle, c’est-à-dire à chercher
une solution approchée du problème. Elle consiste à décomposer la structure en
élément finis limités par des frontières, et à écrire la fonctionnelle globale à
minimiser comme une somme de fonctionnelles relatives à chaque élément[2].
18
L’appellation « éléments finis » vient de la décomposition du domaine d’étude

en éléments qui sont souvent représentés par un maillage des frontières et
assemblage de tous les éléments (Figure 1.1).
Figure 1.1 : Décomposition d’une structure en maillage par éléments

En effet sur les frontières des éléments et éventuellement dans leur domaine,
sont placés des points appelée nœud et les coefficients de la représentation sont
exprimés en fonction des valeurs nodales des champs inconnues, valeurs qui
sont appelées degrés de liberté, et qui constituent les inconnues du problème
discrétisé.
L’assemblage des éléments adjacents est une opération nécessaire pour
exprimer la fonctionnelle globale comme une somme des fonctionnelles
élémentaires. Cet assemblage s’effectue en écrivant l’égalité des inconnues
nodales aux nœuds communs des éléments. L’assemblage une fois effectué ne
laisse subsister que des inconnues nodales indépendantes qui constituent les
degrés de liberté définitifs du problème. Les fonctions d’interpolation sont
choisies en principe de façon à réaliser la continuité des champs inconnus à la
traversée des frontières communes des éléments (conformité). La frontière du
domaine globale est approchée par la réunion des frontières élémentaires des
éléments concernés. Donc la méthode se présente comme une méthode de
résolution des équations différentielles basée essentiellement sur les différentes
formulations énergétique.
19
1.3 Historique
L'idée fondamentale derrière la méthode des éléments finis remonte loin en

arrière. Les grecs par exemple avaient reconnu que l'on peut approcher la
solution d'un problème complexe en le divisant en problèmes plus simples. On
peut par exemple approcher le périmètre d'un cercle en calculant le périmètre
d'un polygone à n côtes, pourvu que n soit suffisamment grand. Il suffit alors de
connaître la longueur d'un segment de droite, problème beaucoup plus simple
que celui de la longueur d'un arc de cercle.
L'application à la solution des équations aux dérivées partielles est
évidemment plus récente et est intimement liée au développement de
l'informatique. Courant [42] a introduit le concept de formulation variationnelle,
qui est à la base de toute méthode d'éléments finis. Pour la méthode de Ritz [63],
on part d'un problème posé dans un espace de dimension infinie. On approche
ensuite la solution du problème initial en cherchant une solution dans une suite
croissante de sous-espaces de dimension finie. Ces problèmes approchés sont en
général beaucoup plus facile à résoudre. On peut de plus espérer que la solution
du problème en dimension infinie peut être obtenue par un passage à la limite.
Le choix des fonctions de base constituant ces espaces de dimension finie est
délicat et initialement on les construisait globalement sur le domaine. C'est
Courant qui eut l'idée d'introduire des fonctions à support local qui simplifient
grandement leur construction. La théorie derrière la méthode des éléments finis
a pris une forme plus rigoureuse avec les travaux de Strang et Fix [67].
L’origine de la méthode des éléments finis peut se trouver dans les travaux
de Fermat et Bernouilli (1743) avec le calcul des variations, puis il faut attendre
le début du XX èmesiècle avec les progrès en analyse avec la méthode de
Galerkin se basant sur des théorèmes de projection dans les espaces de Hilbert.
20
En 1943 Robert Courant introduit le principe variationnel avec des fonctions de

base à support locaux ouvrant la voie à une division d’un domaine considéré en
« éléments ». Cependant ce n’est qu’avec le développement des ordinateurs que
ces travaux trouvent leurs applications avec les travaux pionniers de
Zienckiewiz et Argyris qui définiront la méthode en 1960. Ce qui amène le
succès de la méthode et sa puissance est l’apport du calcul matriciel de grande
dimension ce qui permettra à la méthode connaître un développement fulgurant
accompagné par les progrès de l’informatique. En conséquence, la méthode des
élémentsfinis est une méthode puissante basée sur une théorie mathématique
rigoureuse.
Aujourd’hui, les élémentsfinis sont un outil majeur, incontournable en

mécanique de solide, de fluides des interactions et des structures et applicable
dans de nombreux domaines impliquant des problèmes d’équations aux dérivées
partielles comme par exemple en mathématiques financières ou
l’électromagnétisme. Aussi de nombreux codes industriels (solveurs) existent et
sont généralement couplés à des logiciels de CAO ou Computer Aided Design
(CAD) ; par exemple : Ansys, Abaqus, Robot, LS-dyna, Feap, Code-Aster,
Cast3M et bien d’autres.
Par ailleurs, les bases théoriques de la M.E.F reposent d’une part sur la
formation énergétique de la mécanique des structures dont la formulation des
théorèmes énergétiques de l’élasticité a été effectuée au siècle dernier, en 1819
NAVIER définit une méthode d’étude des systèmes hyperstatiques basée sur
l’application des conditions d’équilibre et de compatibilité par MAX WELL en
1864 et CASTIGLIANO en 1878 [3]. D’autres part les méthodes
d’approximation, c’est au début du 20ème siècle qu’ont été acquis des résultats
fondamentaux dans le domaine des méthodes d’approximation sous l’impulsion
de RITZ en 1902 et de GARLEKIN en 1915 puis en 1943 COURANT établit les
21
bases de la M.E.F.
En 1954 DENKE systématise la méthode des forces. En fin de 1955

AGURIS présente une approche unifiée des méthodes de déplacement et des
forces, puis l’année suivante TURNER et CLOUTH publient une présentation
systématique de la méthode des déplacements, ces deux publications sont
particulièrement importantes et représentent véritablement le début de la M.E.F
comme technique de calcul des structures complexes. Dès lors, on laisse au
développement de nouveaux éléments tels que membranes, plaques, coques. Par
ailleurs le domaine d’application de la M.E.F limité au début à la statique
linéaire s’étend maintenant à la dynamique.
Cette historique perçut incomplète si l’on omet de mentionner le

développement de programmes généraux d’analyse (encore appelés codes
généraux), à partir des années 60, ce phénomène a été particulièrement
important parce qu’il a véritablement abouti à faire rentrer la M.E.F dans la
pratique industrielle.
Actuellement, la méthode des éléments finis est reconnue comme étant une
technique très puissante pour l’analyse des contraintes et des déformations à
l’état statique ou dynamique des structures discrètes ou continues, dans de
nombreux domaines d’industrie modernes afin leur prévoir un bon concept et un
bon dimensionnement afin d’éviter leur rupture et leur désastre.
1.4 Problématique
Il est certain que d’une part, la M.E.F se prête bien à la programmation sur
ordinateur (les procédures numériques peuvent être rendues automatiques et
modulaires) et que d’autre part, elle se caractérise par son universalité et son
adoptabilité au traitement des problèmes les plus divers. Pour illustrer le
fonctionnement de la méthode des éléments finis et pour justifier l'introduction
d'un certain nombre d'outils mathématiques, la convention de sommation par
22
rapport aux indices répétés (convention d’Einstein) est utilisée. Tandis que la
convention des dérivées partielles est désignée par des indices précédés d’une
virgule. En conséquence, le calcul tensoriel joue un rôle important dans
l’utilisation de cette méthode. La MEF est maintenant reconnue comme l'une
des principales méthodes de résolution des équations aux dérivées partielles
(EDP) dans des géométries quelconques, que ce soit en dimension un, deux ou
trois. On trouve même des méthodes d'éléments finis en dimension 4, soit en
espace-temps...
En effet, la MEF est une technique numérique de modélisation appliquée

pour l’analyse statique et dynamique des structures mécaniques discrètes ou
continues.Si la structure possède un nombre important des paramètres inconnus,
la méthode des éléments finis discrétisée la structure en des parties
infinitésimales pour faciliter la solution.Actuellement, la méthode des éléments
finis reconnue comme une technique de résolution des équations aux drivées
partielles dans des géométries quelconque.
1.5 Applications et état de l’art
Les premières applications véritables de la méthode des éléments finis se

retrouvent en 1956 en mécanique des structures quand un groupe de chercheurs
(Turner, Clough, Martin et Topp [37]) de Boeing utilisent cette méthode pour
calculer la voilure d'un avion. La méthode des éléments finis est maintenant
reconnue comme l'une des principales méthodes de résolution des équations aux
dérivées partielles (EDP) dans des géométries quelconques, que ce soit en
dimension un, deux ou trois. On trouve même des méthodes d'éléments finis en
dimension 4, soit en espace-temps...
Les applications de la méthode des éléments finis sont très variées ettrès
nombreuses dans plusieurs domaines d’ingénierie et diverses disciplines
23
scientifiques et physiques, que ce soit en construction mécanique,charpente

métallique, mécanique des solides, élasticité, plasticité, fracture mécanique,
mécanique des fluides, aéronautiques et aérospatiales, génie civil, ainsi qu’en
d’autres problèmes thermiques, électro- magnétiques, chimiques,
biomécaniquesetbien d’autres.
En effet, ce champ d’application est très vaste carelle afait preuve deson
efficacité de calcul dans différents cas de problèmes simples ainsi que des
problèmes de calculs de grande complexité à l’état statique et l’état
dynamique.Il couvre toutes les applications de la mécanique des structures
(statique linéaire, plasticité, flambage, matériaux composites, dynamique, chocs,
frottements...) et aussi de la mécanique des fluides, la rhéologie, les échanges
thermiques et les calculs électromagnétiques[15].
Le but de l'analyse des structures mécaniques utilisant la MEF consiste à
déterminer les effets des charges statiques et dynamiques sur les pièces, les
assemblages, les mécanismes et les structures afin d'éviter tout risque de
défaillance lors de leur utilisation.En réalité, la conception pilotée par la
simulation intervient au début du processus de développement de produits. Elle
aide à détecter les défauts dès le début, permet de réduire les coûts, accélère la
conception et améliore la qualité du produit.
En fait, Les applications sont tout aussi nombreuses et variées et les
ingénieurs de diverses disciplines utilisent la MEF, que ce soit en mécanique des
fluides ou des solides, mais aussi pour les problèmes thermiques, électro-
magnétiques, chimiques, etc. On retrouve aussi des applications en physique, et
notamment en astrophysique. Les principaux domaines suivants utilisant cette
méthode sont:
1.5.1 laconstruction mécanique
La méthode des éléments finis constitue aujourd'hui un outil de base

couramment utilisé dans l'industrie mécanique. Elle consiste à remplacer et à
24
enrichir l’expérimentale qui est très couteuse et à mieux comprendre les

principes physiques des systèmes mécaniques de géométrie compliquée ou de
grandes dimensions afin d’optimiser les process et/ou les produits existants ainsi
qu’à concevoir de nouveaux produits.L’application de la MEF dans l’industrie
mécanique se résume essentiellement en l’analyse et la modélisation des pistons,
des engrenages, des châssis et carrosseries de véhicules, des structures discrètes
et de la charpente métallique, etc.(figures 1.2).
(a) Structure mécanique (b) Structure discrète en charpente métallique
(c) Structure de carrosserie de voiture (d) Structure continue - engrenage
Figure 1.2 : Applications de la MEF en Mécanique
25
1.5.2 Aéronautique et l’aérospatiale
La MEF joue aussi un rôle important dans le calcul et le dimensionnement

des structures de l’aéronautique et de l’aérospatiale vue leur dimension et leur
complexité afin d’augmenter la vitesse des avions (du subsonique au
supersonique et sonique) et de choisir les matériaux convenables de ces
structures pour stabiliser l’écoulement aérodynamique autour des ailes et de
minimiser les éventuels tourbillons. En fait, on retrouve les premières
applications véritables de cette méthode en 1956 en mécanique des structures où
un groupe de chercheurs (Turner, Clough, Martin et Topp [37]) de Boeing ont
utilisés cette méthode pour calculer la voilure d'un avion.
Aujourd’hui, cette méthode est souvent utilisée pour l’analyse aérodynamique

des écoulements autour des ailes et la détermination des contraintes élastiques et
aéro élastiques des les structures internes des avions (fuselage, ailes et hélices
d’hélicoptères) afin d’éviter leurs fractures et leurs désastres(figure 1.3).
(a) Structure interne d’un avion (b) Structure interne d’une aile d’avion
26
(c)Modélisation aérodynamique d’un avion (d) Structure de la station aérospatiale
Figure 1.3 : Applications de la MEF en Aéronautique et Aérospatiale
1.5.3 le Génie Civil
La MEF est souvent utilisée en génie civil où les structures des bâtiments,
des ponts et chalets sont souvent constitués en plaques, poutres, treillis et
charpente métallique assemblées par rivetage, soudage ou coulage. Elle est
utilisée pour l’analyse et le désigne de ces structures pour leur prévention contre
le séisme et les vibrations de grandes amplitudes afin d’augmenter leur durée de
vie et éviter leur destruction.
(a) Structuret en charpente métallique (a) Structure d’un pont en arc
27
(c) Structure en construction bâtiments (d) Structure discrète d’une grue
Figure 1.4 : Applications de la MEF en génie civil
1.5.4 la Robotique et la Biomécanique
Les systèmes robotiques et biomécaniques sont souvent constitués

d’éléments assemblés en articulation pour pouvoir atteindre les différents
mouvements en trois dimensions. Leur champ d’application est réalisé dans
plusieurs domaines industriels et biomécaniques et peut être analysé par la
méthode des éléments finis afin de déterminer le mouvement, le déplacement et
la rotation tridimensionnelle des différentes parties constituants le robot. Le
design et la conception des bras d’un robot est basée sur des systèmes et des
matériaux intelligents à savoir combien il fléchit s’il soulève un poids bien
déterminé. Il est extrêmement important de concevoir des bras de robot pour être
très rigide, comme la flexion réduit considérablement la précision de position.
En conséquence, le bras doit être aussi léger que possible pour économiser de
l'argent sur les actionneurs « matériaux intelligents » qui sont utilisés pour leur
commande à distance. C'est-à-dire économiser la puissance des moteurs car les
brasde robot plus lourds ont besoin de moteursplus forts.
28
(a) Mouvement tridimensionnel d’un robot (b) Robot spatial

gratte ciel
(c) Mouvements biomécaniques (d) Mouvement d’un bras
(e) Mouvements d’un jambe
Figure 1.5 : Applications de la MEF pour le mouvement des robots

1.6 Conclusion
En conclusion, les structures considérées dans les applications sont souvent

constitués d’éléments poutres et plaques assemblées en maillage important,
chose qui nous a permis de choisir et de considérer dans notre travail une
analyse des structures discrètes et continues utilisant les éléments barre, poutre
et plaque.
29
Chapitre 2 :Techniques de modélisation utilisées en éléments finis
Chapitre 2 :
Techniques de modélisation
utilisées en éléments finis
30
2.1 Introduction
Les techniques habituelles de modélisation utilisées en éléments finis

reposent essentiellement sur la méthode des forces – déplacements dans un
comportement linéaire, c'est-à-dire les déplacements en tous points sont
proportionnels aux efforts extérieurs. Cela impose que le matériau ait un
comportement linéaire (élastique) et que les déformations soient
proportionnelles aux déplacements (petits déplacements).Le calcul d’une
structure par la méthode des éléments finis exige, pour être efficace, d’accorder
un plus grand soin à la phase de modélisation. La précision des résultats dépend
de la qualité du modèle et également du maillage réalisé (nœuds et éléments).
2.2 Principes de formulation de la méthode des éléments finis
2.2.1 Principe de la méthode des forces-déplacements
La méthode des forces-déplacements repose essentiellement sur l’équation

d’élasticité décrite mathématiquement par la loi :
 Forces   Matrice de   Déplacemen ts 

     généralisé 
 généralisé es    Rigidité de  
 aux noeuds   la Structure   aux noeuds 

  
Où le vecteur des forces est supposés inconnues et les déplacements sont

connues. Inversement, si les forces aux nœuds sont connues et les déplacements
doivent être déterminés la relation s’exprime comme suit :
 Déplacemen ts   Matrice de   Forces 
     généralisé es 
 généralisé    Camplaisance  
 aux noeuds  de la Structure   aux noeuds 
  
La matrice de complaisance est l’inverse de la matrice de rigidité.Et si la

matrice de rigidité est de grande dimension, son inverse devient plus difficile.
Ceci est fonction des dimensions de la structure, de son maillage et de
l’assemblage de tous ses éléments et de leurs conditions aux limites. La méthode
31
des éléments finis a donc connu un développement intense à partir des

annéescinquante sous l’impulsion et la disponibilité des premiers ordinateurs
[1]. Dés lors on a assisté au développement de nouveaux éléments
tridimensionnels tels que : treillis, tétraèdres, membranes, plaques orthotropes,
coques…etc.
2.2.2 Principe de la méthode variationnelle

Le principe variationel consiste à mettre en place une formulation de la
méthode des éléments finis en se basant sur des outils de mathématiques
appliqués. Il s'agit de mettre en place, à l'aide des principes hérités de la
formulation variationnelle ou formulation faible, un algorithme discret
mathématique permettant de rechercher une solution approchée d’une équation
aux dérivées partielles sur un domaine compact avec des conditions aux limites
ou aux bornes et/ou dans l'intérieur du compact [8]. Il s'agit donc avant tout de la
résolution d'un problème discrétisé où, grâce à la formulation variationnelle, les
solutions du problème vérifient des conditions d'existence plus faibles que celles
des solutions du problème de départ. Comme de nombreuses autres méthodes
numériques, outre l'algorithme de résolution en soi, se posent les questions de
qualité de la discrétisation telles que: l’existence et l’unicité de la solution, sa
stabilité et sa convergence, et en fin la minimisation de l’erreur (entre la solution
discrète du problème et sa solution unique initiale).
La MEF peut donc aboutir à des systèmes intégro-différentiels basés sur

des équations aux dérivées partielles linéaires dont on cherche une
approximation par des fonctions affinées. Cette méthode permet donc de
résoudre de manière discrète une EDP dont on cherche une solution approchée
« suffisamment » fiable. Elle comporte des conditions aux bords permettant
d'assurer existence et unicité d'une solution. La discrétisation passe par la
définition d'un espace de fonctions tests approprié, sur lequel les solutions de la
32
(re)formulation variationnelle de l'équation est exacte. Cela nécessite la

définition d'un maillage du domaine en fragments quelconques : les éléments
finis. Ces fragments peuvent par exemple être de forme quelconque mais
doivent former un pavage de l'espace considéré. Dans ce cadre, et moyennant
une attention toute particulière au problème des bords du domaine, on obtient
une formulation algébrique dite discrétisation du problème initial. La solution de
ce problème algébrique, si elle existe elle est unique.En conclusion, les
équations développées sont généralisées à partir des éléments finis.Leur
résolution est basée sur les conditions aux limites entre les éléments et leurs
conditions aux limites vues que le système considéré existent réellement et qui
est d’un ressort d’ingénierie.
2.2.3 Champ d’application

La méthode des éléments finis est une méthode de calcul numérique qui,
ayant un profond caractère plus physique qu’abstrait et a été inventée plutôt par
les ingénieurs que par les mathématiciens.Cette méthode a été appliquée pour la
première fois dans des problèmes liés à l’analyse des contraintes et depuis, elle a
été étendue dans d’autres problèmes liés au milieu continu, fluide ou thermique.
Dans toutes les applications utilisant cette méthode, l’analyse recherche à
calculer l’état de quantité et de qualité dutype de champ considéré, comme par
exemple :
Application Type de champ

Analyse des contraintes Champ des contraintes ou champ des déplacements
Analyse thermique Champ de température ou flux de chaleur
Ecoulement des fluides Fonction de courant ou fonction du potentiel de
vitesse
Tableau 2.1 : Type des champs d’application de la MEF
33
2.3 Techniques de modélisations utilisées en éléments finis
2.3.1 Méthode de maillage (nœuds et éléments)

La MEF pourrait être définie sous la forme de petits éléments assemblés
ensemble en des points appelés nœuds. Suivant la structure à analyser discrète
ou continue, elle est divisée en plusieurs éléments (petites pièces) et assemblées
par soudures ou rivetage pour les structures discrètes en des points appelés
«nœuds» ou en des points de contact pour les corps continus. Ces éléments sont
ensuite reconnectés par l’intermédiaire des nœuds pour reconstruire le corps
entier (fig.2.1).
Éléments Nœuds
Figure2.1Discrétisation d’une structure en nœuds et éléments
Le comportement de chaque élément est décrit par un set d’équations

algébriques. Dans l’analyse des contraintes ces équations sont des équations
d’équilibre des nœuds. Du fait que le nombre de ces équations est très grand
(centaines ou milliers), l’utilisation d’un ordinateur est absolument obligatoire.
Autrement dit, dans un élément, une quantité de champ (ex. le champ de
déplacement) est interpolé à partir des valeurs existantes dans les nœuds. En
connectant les éléments ensemble, la quantité de champ devient interpolée sur
l’entier de la structure. Les meilleures valeurs de la quantité de champ dans les
nœuds sont celles qui minimisent certaines fonctions (telle que l’énergie totale).
Le processus de minimisation génère un set d’équations algébriques simultanées
pour les différentes valeurs de la quantité de champ dans les nœuds.
34
Ce set d’équations est décrit sous forme matricielle par :

F  K   (2.1)
Où :{F} = Tenseur des Forces (Forces nodales) .

[K] = Matrice des constates (Matrice de rigidité) .
 = Tenseur de déplacement (déplacement aux nœuds).
2.3.2 Etapes de modélisation par la MEF

Pour faire une analyse par la MEF on doit prendre soin de la modélisation :
choisir une numérotation convenable des nœuds afin de minimiser la matrice de
rigidité du corps considéré ainsi que le type d’élément finis approprié et
convenable qui puissent représenter de façon adéquate l’action physique du
corps. Ainsi les principales étapes pour une analyse utilisant la méthode des
éléments finis sont représentées comme suit :
Un problème à résoudre
Anticiper le comportement
physique. Pensez comment les Solutions
résultats offerts par la MEF L’analyse par éléments
OUI expérimentalesou
seront vérifiéspour voir s’ils finis est-elle nécessaire ? NON
analytiques
sont raisonnables
Concevez un modèle STOP

initial à l’aide des
éléments finis
Preprocessing:
Preparer le modèle
Penser à réviser le modèle en

Procession : améliorant le modèle déjà
Résoudre les équations du existant
modèle proposé
NON
Post procession : Les résultats sont-ils raisonnables ?

Afficher les résultats obtenus Les erreurs estimées sont petites ?
suite au roulage
OUI
LOGICIEL
STOP
Figure2.2 :Principales étapes pour une analyse par la M E F
35
Il est souhaitable de ne pas utiliser des éléments « déformés » ou des

éléments grossiers pour représenter des variations considérables d’une certaine
quantité de champ. A l’autre extrême un sur-raffinage pourrait conduire à une
perte de temps pour l’analyste ainsi qu’à un sur chargement de la mémoire de
l’ordinateur. Cependant, même si un grand nombre d’éléments est utilisé dans la
discrétisation il y a une erreur dénommée erreur de discrétisation qui existe du
fait que la structure physique et le modèle mathématique ont une infinité de
degrés de liberté (qui sont au fait les déplacements pour une infinité de points de
la structure), tandis que le modèle avec éléments finis a un nombre fini de
degrés de liberté.
A part les erreurs introduites par l’analyste lors de la discrétisation,
l’ordinateur introduit des erreurs numériques par l’arrondissage ou par le
tronquée des nombres qui sont introduits dans les matrices et qui servent à la
résolution des équations. Les logiciels d’analyse par éléments finis sont devenus
sur une large échelle des instruments de calcul facile à utiliser et peuvent
afficher les résultats sous une forme très attractive. Même le plus inhabile
utilisateur peut offrir une réponse quelle qu’elle soit. Mais une carte joliment
colorée des contraintes et des déformations peut être obtenue pour n’importe
quel modèle bon ou mauvais. Il se peut que la plupart des analyses par MEF
soient si défectueuses qu’elles ne peuvent pas être dignes de confiance.
Un utilisateur responsable doit comprendre suffisamment bien la nature
physique du problème et le comportement des éléments finis afin de préparer un
modèle convenable et de bien évaluer la qualité des résultats. La responsabilité
des résultats obtenus revient à l’ingénieur qui utilise le logiciel et non pas au
vendeur de ce logiciel, même si les résultats sont affectés par les erreurs du
programme. A présent il existe un grand nombre de logiciels pour l’analyse par
la méthode des éléments finis tels que : ANSYS, ABAQUS, COSMOS-M,
PATRAN, IDEAS, SAP2000 et autres.
36
En résumé, pour qu’un problème mécanique soit résolu utilisant la MEF,

il faut que les données d’entrée (procession) soient introduites dans l’ordinateur
utilisant un programme choisi, puis la post procession (pré procession) sera
introduite en tenant compte des charges, des contraintes ou des déplacements
appliquées et enfin attendre après un certain temps l’obtention des données de
sortie (voir schéma ci-dessous).
Procession des données d’entrée(Post procession)

- Chargements
- Contraintes
-Déplacements
- Propriétés des matériaux
-Températures
Données d’entréeDonnées de sortie

(Pré procession) - Génération du réseau
- Coordonnées nodalesd’éléments finis
- Supports (appuis etc.)
- Blocages (conditions à la limite)
2.4 Classification des éléments de la MEF
2.4.1 Introduction à l’état plan de contrainte et de déformation

En cas où la variation du déplacement ou de la contrainte sont négligeables
tout le long de l’axe z (la direction normale au plan d’analyse), on considère que
le corps est un problème plan (état plan de contrainte ou état plan de
déformation). Si par contre les déplacements et les contraintes peuvent varier
dans toutes les directions x, y ou z la structure en cause peut être appelée tri
dimensionnel (ou Solide 3D). Un cas spécial de solide ayant symétrie axiale
(ex : une cloche) s’appelle de façon usuelle « Solide de révolution ». Les
chargements à leur tour peuvent ou non être distribués de façon axiale
37
symétrique. Une plaque plane qui supporte des chargements dans son plan est un
problème plan. Par contre, si la plaque est chargée par des forces qui n’agissent
pas dans son plan, cela représente un problème de flexion de plaque ou, plus
simplement un problème de plaque (en anglais : PLATE PROBLEM ). Par
contre, si la plaque est courbe elle devient une coque (SHELL). Les réservoirs,
par exemple, peuvent être considérés dans l’analyse par la MEF, des coques[3].
En conclusion les éléments finis peuvent être divisés en plusieurs
catégories en fonction de la structure : éléments plans, éléments solides 3D,
éléments solides à symétrie axiale, éléments de plaque, éléments de coque. On
trouve de même des éléments de barre articulée (en anglais : TRUSS), éléments
de poutre (en anglais : BEAM), éléments de fondation élastique etc.
2.4.2 Connaissances requises pour la réalisation des logiciels MEF
La MEF a un caractère pluridisciplinaire. Pour pouvoir réaliser des logiciels

qui puissent résoudre un certain type de problèmes dans le domaine du génie
mécanique, il s’impose de maitriser les disciplines suivantes :
- la mécanique des structures à l’état statique, dynamique, vibratoire et
la résistance des matériaux ;
- l’analyse numérique c'est-à-dire des procédés et algorithmes de calcul,
et de graphes sur l’ordinateur ;
- la programmation linéaire (C++, Pascal etc.)
Ainsi les familles des grands logiciels, avec de multiples facilités, réalisés
par des compagnies spécialisées qui sont généralement utilisés par des collectifs
de recherche, font partie, entre autres, NASTRAN, ANSYS, COSMOS,
ALGOR, IMAGES3D.. etc.
2.4.3 Connaissances nécessaires à l’utilisateur de la MEF
Un utilisateur de la MEF est mis dans la situation de résoudre un certain
38
problème. On doit mentionner dès le début que le logiciel appliqué au problème

respectif ne le résout pas. Il ne fait que résoudre un modèlecréé par l’utilisateur.
Les résultats peuvent être confirmés ou pas, en fonction du modèle choisi par
l’utilisateur. La modélisation est une activité de simplification de la structure en
l’encadrant ses différentes portions dans une des catégories suivantes : barres,
plaques, blocs massifs, en tenant compte des chargements, appuis etc. La
modélisation correcte (la plus proche de la réalité) est un problème d’expérience,
d’inspiration et moins de la connaissance des fondements théoriques de la
méthode.
2.5 Discrétisation des structures par éléments finis
2.5.1 Types d’éléments

Bien que pour l’étude du comportement élastodynamique des système
continus, la forme faible discrète s’appuyant sur la méthode de Galerkin, soit
une approche à première vue plus puissante que la formulation forte basé sur
l’écriture directe des équations différentielles du mouvement [1]. Elle laisse
rapidement entrevoir ses limites en présence de géométries complexes, son
inconvénient principal résidant dans le fait que les fonctions de forme qui
caractérisent les champs des déplacements réels et virtuels approchés sont
définies globalement sur l’ensemble du domaine étudié. La méthode des
éléments finis contourne cette difficulté en subdivisant le corps déformable de
forme compliquée en un assemblage de sous-domaines ou éléments simples et
en construisant les solutions approchées sous la forme d’une combinaison de
fonctions de forme à support compact. Dérivée en tant que cas particulier
d’application de la méthode de Galerkin, elle constitue un procédé numérique
général et systématique pour la discrétisation spatiale des structures et connait
depuis plusieurs années un succès considérable, notamment grâce au
développement de nombreux codes d’élément finis, dont le domaine
39
d’application s’étend largement au-delà de l’élastodynamique linéaire. Les

structures mécaniques courantes comportant souvent des zones ou des portions
dont l’épaisseur est relativement faible en comparaison des deux autres
dimensions géométriques. Comme il n’est pas économique au niveau
algorithmique (consommation exagérée de temps de calcul et mauvaise gestion
de la mémoire) de modéliser de telles parties structurelles par des éléments finis
solides, il est nécessaire de développer des éléments finis spécifiques
La Méthode des Eléments Finis a développé une série de types d’éléments
finis qui, pour le début, peuvent être classifiés en :
- éléments finis unidimensionnels (généralement des barres);
- éléments finis bidimensionnels (plaques et mêmes volumes) ;
- éléments finis tridimensionnels (blocs massifs).
paraboliques
Eléments Linéaires Cubiques
(quadratiques)
Unidimensionnels
bidimensionnels
tridimensionnels
autres types Masse Ressort Contact
Tableau 2.2 Types d’éléments finis
Les éléments finis sont générés par des points qui ne sont que des nœuds
de la structure. Il existe des éléments ayant un degré supérieur à ceux cubiques
40
(qui sont les plus performants) mais le plus couramment sont utilisés les
éléments linéaires et paraboliques.
Certains éléments finis ont des nœuds intérieurs pour améliorer la
précision, mais la méthode des éléments finis ne travaille pas avec ces nœuds.
Ils sont générés et ensuite condensés dans la phase de calcul des matrices de
rigidité des éléments.
2.5.2 Assemblage
Dans les deux cas d’élément poutre ou élément barre, et après la
formulation des équations de rigidité de tous les éléments on passe à la
formulation de l’équation de rigidité globale de la structure. La procédure
consiste à superposer tous les éléments agissant aux nouds et par conséquent
l’assemblage des vecteurs forces généralisées, des matrices de rigidité et des
vecteurs déplacements généralisées de tous les éléments de la structure, pour
éclaircir la procédure d’assemblage on considère une partie de la structure
représenté sur la (figure.2.3)
Figure2.3 : Représentation des nœuds et les éléments
La première étape consiste à éliminer les conditions aux limites existantes

aux nœuds, et essayer de faire apparaître tous les degrés de liberté et les
composantes des forces existantes aussi aux nœuds comme représentés sur la
(figure2.4).
41
Figure2.4 :Forces et déplacements de quatre

nœuds.
Après la formulation de l’équation de rigidité de chaque élément, on
effectue l’assemblage tel que les déplacements nodaux pour chaque élément sont
les mêmes que les déplacements aux points communs.
Figure2.5:Assemblage des nœuds
Les forces aux nœuds pour chaque élément doivent être en équilibre avec
les forces de chaque nœud commun. Alors la superposition des efforts se traduit
42
par les équations suivantes :
𝐹𝑥1 = 𝑃𝑥1 (é𝑙é𝑚 1) + 𝑃𝑥1 (é𝑙é𝑚 2) + 𝑃𝑥1 (é𝑙é𝑚 3)

Fy1 = Py1 (élém 1) + Py1 (élém 2) + Py1 (élém 3)
Fx2 = Px2 (élém 1) + Px2 (élém 2) + Px2 (élém 3)
L’assemblage est lié essentiellement à la numérotation de la structure et la

manière de la numérotation joue un rôle très important sur la grandeur de largeur
de la bande de la structure.
2.5.3 Elimination directe de la bande de matrice

Par une numérotation judicieuse des degrés de liberté la structure de la
matrice de rigidité globale est à bande diagonale, ce qui permet l’utilisation
d’algorithmes d’inversion appropriée d’une manière économique. Dans ce cas
on cherche à réduire en générale la largeur de bande moyenne ou maximale.
Mais les méthodes modernes de stockage en mémoire d’ordinateur ont modifié
largement cette exigence. Pour réduire la largeur de bande, la numérotation
devra être telle que tous les éléments couplés aient des degrés de liberté dont la
numérotation soit aussi voisine que possible [4]. Il existe quelques règles assez
générales, mais le cas de domaines multiplement connexes est encore d’une
optimisation délicate.
Pour les structures à deux dimensions, et si les points nodaux sont bien
numérote, tel que la largeur de la bande de la matrice de rigidité peut être plus
étroites. Et donc les problèmes de stockage peuvent être réduites, par le stockage
des coefficients dans la bande de l’ordinateur, et les équations peuvent être
résolues en utilisant la méthode d’élimination directe aux lieux des itérations, la
largeur de la bande pour une structure assemblée d’éléments finis peut être
43
trouvé par l’équation générale suivante

L = 2[mn + (n -1)] +1
m : La différence maximale entre deux nœuds d’un même élément fini.

n : Nombre de degrés de liberté pour chaque nœud.
2.6 Types de numérotation

2.6.1 Numérotation verticale
Dans cet exemple tel que la numérotation des points nodaux est verticale, la
largeur de la bande de la matrice pour une structure à deux dimensions, c’est -à-
direchaque nœud possède deux degrés de liberté égale:
n=2m=4
L = 2[4´ 2 + (2 -1)] +1 = 19
2.6.2 Numérotation horizontale
Dans cet exemple tel que la numérotation des points nodaux est
horizontale, la largeur de la bande de la matrice pour une structure à deux
dimensions, c’est -à- dire chaque nœud possède deux degrés de liberté égale :
n=2m=5
44
L = 2[5´ 2 + (2 -1)] +1 = 23
Donc la largeur de la bande est minimale dans le premier système de
numérotation c.-à-d. la numérotation verticale.
Tableau 2.3Bande d’une matrice à deux dimensions
2.7 Calcul d’erreurs dans l’analyse par MEF

2.7.1 Erreurs usuelle
Une matrice de rigidité globale [K] peut être singulière au cas où :
- Les propriétés du matériau (modules d’élasticité) sont nulles.
- Unou plusieurs nœuds ne sont pas connectés à un élément quelconque.
- Le manque de supports, ou bien les supports sont insuffisants pour
prévenir des mouvements de corps rigide.
- Un mécanisme est créé du fait que le modèle (une partie du modèle) est
insuffisamment bloqué.
- Il y a de grandes différences en ce qui concerne les rigidités.
- Une partie de la structure a flambé.
- Dans une analyse non-linéaire les supports ou les connexions ont atteint
45
unerigidité nulle, donc une partie de la structure est insuffisamment

bloquée.
Si une matrice [K] est singulière, elle déclenche un message d’erreur et arrête
le processus d’analyse numérique. Si l’exécution s’arrête ou l’exécution
continue mais les résultats sont bizarres il est clair que quelque chose est louche
et dans ce cas une recherche de l’erreur est absolument nécessaire. Il est plus
dangereux s’il existe des erreurs qui puissent conduire à des résultats plausibles,
mais pas précis, exacts. Dans cette catégorie on peut énumérer :
- une utilisation erronée des types d’éléments, par exemple, les éléments
de coque sont utilisés à la place des éléments de type “solide” ou bien
des éléments de symétrie axiale à la place des éléments plans.
- les appuis (supports) sont mal disposés en direction, type ou location.
- les chargements sont mal disposés en direction, type, location ou
grandeur.
- autres données peuvent être incorrectes, par exemple une utilisation
incorrecte des unités de mesure (vitesse angulaire  en rot/sec au lieu
de rad/sec).
- un élément peut être défini plusieurs fois, la duplication est très
difficile à détecter à cause du fait qu’il est très difficile d’observer
l’erreur lors du traçage.
- une connexion peut être sans aucune signification du point de vue
physique.
-
2.7.2 Erreurs détectées : Vérification du modèle
Dans ce cas des vérifications sont faites par l’analyste (l’utilisateur) ou
automatiquement faites par le logiciel. Et dans les deux cas plusieurs
contestations seront observées (Tableau 3).
46
Vérification faite par l’analyste Vérification automatique faite par le

(l’utilisateur) logiciel
vérification manuelle de la location un certain nœud n’est pas connecté à un
de certains nœuds. élément quelconque.
vérification de l’introduction des les nœuds sont tout près l’un de l’autre
propriétés des matériaux. ou coïncident mais ne sont pas connectés.
vérification de l’introduction de les éléments partagent un nœud mais
l’aire des sections transversales n’utilisent pas le même set de degrés de
pour les poutres. liberté que celui-ci.
affichage de la frontière du modèle, le coefficient de Poisson n’est pas dans
au cas où une partie du contour a l’intervalle (0 ; 0,5) ; un test analogue
une forme de “fissure” dans un peut détecter des propriétés impossibles
matériau, quelques nœuds sont soit d’un matériau orthotrope.
adjacents soit coïncident.
Tableau 2.4Vérification des erreurs
2.8 Critique des résultats obtenus par MEF

Au cours d’une post-procession les déplacements sont ceux qui doivent
être examinés premièrement. D’habitude les logiciels tracent seulement des
lignes droites entre les nœuds, de telle manière que les formes courbes
supposées être des éléments déformés des poutres, les bords des éléments de
plaque ou de coque, ne sont pas visibles.
Au cas où une symétrie est anticipée cet aspect serait visible dans le champ
de déplacements. Les déformations produites par une variation de température
sont tracés par le logiciel seulement si la configuration non-déformée existe à
une température de 0o.
Les réactions dans les supports doivent être examinées afin de vérifier si elles
satisfaisant les équations d’équilibre statique, par exemple de vérifier si la
somme des composantes des réactions sur la direction Oxéquilibre les forces
appliquées sur cette direction. A noter que toutes les réactions doivent être
rapportées au même système de coordonnées.
Avant d’examiner les résultats qui concernent les contraintes, avant de faire
47
une évaluation, on doit se poser les questions suivantes :

 est-ce qu’elles les coordonnées aux systèmes globaux ou locaux
d’axes ?
 au cas où le système d’axes est bien défini, comment ces axes sont-
elles orientées ?
 les contraintes résultantes sont-elles reportées (par exemple celles dues
aux moments résultants) ?
 dans les éléments de poutre, de plaque ou de coque, les contraintes
peuvent rester valides à la surface supérieure, moyenne ou inférieure ?
 est-ce que les contraintes sont concentrées dans les nœuds ?
 ceci serait incorrect au cas où le système d’axes de coordonnées ne
s’adapte pas ou s’il existe des discontinuités d’épaisseurs ou des
discontinuités des propriétés du matériau ;
Quelques caractéristiques qui puissent définir un champ de contrainte précis,
exact sont présentées ci-dessous :
 le contour des contraintes doit être normal à un plan de symétrie
réfléchissant ;
 sur une frontière libre (sans chargements), une des contraintes
principales doit être nulle ;
 sur une frontière chargée par une pression p, une des contraintes
principales doit être ˝- p ˝ ;
 les trajectoires des contraintes principales doivent être normales ou
tangentes aux frontières libres, à celles sollicitées par une pression p ;
 dans les problèmes de symétrie axiale, les contraintes normales -
radiales ou circonférentielles doivent être égales sur l’axe de
révolution ;
 le rapprochement, l’exactitude des contours des contraintes est un
autre guide ; au cas où les contours des contraintes ont des incréments
48
égaux entre eux, les éléments qui présentent une dispersion des
contraintes doivent être plus raffinés que ceux qui ne présentent pas de
dispersion ;
 au cas où certains phénomènes tels que le flambement, les
déformations dans le domaine plastique, l’action variable des efforts
n’ont pas été pris en compte, une nouvelle analyse par MEF devrait
être faite.
49
Chapitre 3 :Détermination de la matrice de rigidité des éléments des structures
mécaniques
Chapitre 3 :
Détermination de la matrice de
rigidité des éléments des
structures mécaniques
50
mécaniques
3.1 Introduction
Dans ce chapitre on va présenter les différentes étapes pour la
détermination des matrices de rigidité pour différents typesd’éléments de
structures mécaniques. Cette analyse sera réalisée en se basant sur l’énergie de
déformation (U) dans l’élément du corps mécanique considéré et le théorème de
Castigliano pour déterminer les vecteurs de forces{F} appliquées aux nœuds
« i »de l’élément.
En conséquence, la matrice de rigidité de l’élément [ki] pourrait être
déterminée. Et la matrice globale [K] de tout le système mécanique pourrait être
obtenue par assemblage de toutes les matrices des éléments élémentaires.Dans
ce contexte, l’élément barre, poutre et plaque sera considéré et
analyséséparément à l’état statique de l’élément et de toute la structure
mécanique.
3.2 Détermination de la matrice de rigidité de l’élément de barre

3.2.1 Matrice de rigidité en coordonnées locales
On considère un élément de barre uniforme, prismatique et élastique, de
longueur L, de module élastique E et d’aire de la section transversale A (fig. 3.1)
en coordonnées locales (xy). Un nœud est localisé à chacune des extrémités de
la barre et les seuls déplacements qui sont permis sur l’élément sont ceux axiaux
suivant l’axe (x).
On déplace d’abord le premier nœud, ensuite le deuxième et dans chaque

cas, on calcule les forces qui doivent être appliquées dans les nœuds pour
maintenir le même état de déplacement (δ). Ces forces sont faciles à déterminer
à partir des formules élémentaires de la résistance des matériaux et données :
EA
F  (3.1)
L
51
mécaniques
En introduisant la théorie des éléments finis, on note que le vecteur de force au
nœud « i » (i =1,2) associée au déplacement du nœud « j » (j = 1,2) suivant l’axe
(x) serait donnait par (Fij). Pour les deux cas les forces sont données par :
EA EA
F11 = F21=  u1 et F12= F22 =  u2 (3.2)
L L
Avec ui le vecteur de déplacement au nœud « i ».
A,E
F11 1 2 F21
u1 L
A,E
F12 1 2 F22
L u2
Figure 3.1 : Elément barre
A une position quelconque (x) le long de barre, la fonction de déplacement u(x)

une fonction linéaire décrite sous forme polynomiale du 1ierdegré donnée par :
𝑢(𝑥) = 𝑎1 + 𝑎2 𝑥
Où 𝑎1 et 𝑎2 sont des constantes à déterminer par les conditions aux limites aux
points nodaux 1 et 2. C'est-à-dire :
Pour x = 0 :
u(0) = u1 = 𝑎2 ;
Et pour x = 𝑙 :
𝑢 (𝑥) = 𝑢 (𝑙) = 𝑢2 = 𝑎1 + 𝑎2 𝑙
Donc:
𝑎2 = u1 ,
𝑢2 − 𝑢1
𝑎2 =
𝑙
Après substitution des coefficients 𝑎𝑖 dans l’équation (3.3), et après
réarrangement on obtient la forme de fonction de déplacement suivante :
52
mécaniques
𝑢(𝑥) = 𝑓1 (𝑥)𝑢1 + 𝑓2 (𝑥)𝑢2 (3.3)
Avec : fi(x) sont les fonctions de formes données :
𝑥 𝑥
𝑓1 (𝑥) = 1 − et 𝑓2 (𝑥) =
𝑙 𝑙
Par ailleurs, en utilisant la méthode de l’énergie pour déterminer la matrice de

rigidité del’élément barre à section droite uniforme, la déformation par unité de
longueur est donnée par la relation suivante :
𝑢(𝑥 + ∆𝑥) − 𝑢(𝑥) 𝜕𝑢 (3.4)

𝜀 = lim =
∆𝑥→0 ∆𝑥 𝜕𝑥
Où ∆𝑥est une longueur infinitésimale fournie par la présence de la force axiale

de la barre, et la déformation axiale de cet élément barre est obtenue par la
substitution des équationsci-dessous dans l’équation (2.3) :
  f ( x)u1  f ( x)u2 (3.5)
Et la force axiale de cet élément barre suivant les coordonnées locales

est donnée par la relation suivante :
u
S  A  EA  EA
x
S  EA( f ( x)u1  f ( x)u2 ) (3.6)
Connaissant que l’énergie de déformation de l’élément barre est
donnée par :
2
l
2
U=  2E
V
dV = 
0
2E
Adx
l l
S2 EA
U  dx   ( f ( x)u1  f ( x)u 2 ) 2 dx (3.7)
0
2 EA 2 O
On passe maintenant à l’application du théorème de Castigliano qui

stipule que :
53
mécaniques
U
Fi =
 ui
En conséquence, les forces axiales aux nœuds 1 et 2 suivant (x) seront

données par :
𝜕𝑈 𝐸𝐴 𝑙 ′
𝐹𝑥1 = = ∫ (𝑓 (𝑥)𝑢1 + 𝑓2′ (𝑥)𝑢2 ) (𝑓1′ (𝑥)) 𝑑𝑥
𝜕𝑢1 2 0 1
l l
 ( EA f1f1dx)u1  ( EA f1f 2dx)u 2 (3.8)
0 0
𝑙 𝑙
𝜕𝑈
𝐹𝑥2 = = (𝐸𝐴 ∫ 𝑓1 𝑓2 𝑑𝑥 ) 𝑢1 + (𝐸𝐴 ∫ 𝑓2′ 𝑓2′ 𝑑𝑥) 𝑢2
′ ′
𝜕𝑢2 0 0
En forme matricielle, l’équation qui exprime la relation entre les charges

nodales est les déplacements nodaux sera :
𝐹𝑥1 𝐾 𝐾12 𝑢1 (3.9)

{ } = [ 11 ]{ }
𝐹𝑥2 𝐾21 𝐾22 𝑢2
Où symboliquement :
Fx   K u (3.10)
Tel que [k] est la matrice de rigidité, et les coefficients de rigidité sont
déterminent par la relation suivante :
𝑙
𝐾𝑖𝑗 = 𝐸𝐴 ∫0 𝑓𝑖′ (𝑥)𝑓𝑗′ (𝑥)𝑑𝑥 (3.11)
Et pour i = 1à2 et j = 1à 2 , et la matrice de rigidité ij k est une matrice

symétrique, et après substitution de l’équation (2.3) dans l’équation
(2.10),l’équation de rigidité pour un élément barre orientée suivant les
coordonnées locales prend la forme finale suivante :
𝐹𝑥1 𝐸𝐴 1 −1 𝑢1
{ }= 𝑙 [ ]{ } (3.12)
𝐹𝑥2 −1 1 𝑢2
3..2.2Matrice de rigidité en coordonnées globales
54
mécaniques
3.2.2.1 Transformation des coordonnées du système local au système
global
Dans la figure 2.1 est représenté un exemple typique de barre ainsi que les
deux repères de coordonnées, local xOy et global XOY. Les déplacements dans
le système local sont u1 , v1 et u2, v2 et dans le système global U1, V1 et U2, V2.
L’angle  est considéré l’angle entre l’axe X et la direction positive de la barre
1-2.
Y
y
2’
déformée
1’ 
x
v2 V2
2 u2
v1 V1
U2 
u1 initiale
1 X
U1

Figure3.1 : Système local et global de coordonnées

Les déplacements dans le système local xOy peut être exprimé en fonction des
déplacements globaux.
𝑢1 = 𝑈1 cos 𝜃 + 𝑉1 sin 𝜃
(3.13)
𝑢2 = 𝑈2 cos 𝜃 + 𝑉2 sin 𝜃
ou sous forme matricielle :

𝑈1
𝑢1 cos 𝜃 sin 𝜃 0 0 𝑉
{𝑢 } = [ ] . { 1 } ⇒ {𝑢} = {𝑇}{𝑈} (3.14)
2 0 0 cos 𝜃 sin 𝜃 𝑈2
𝑉2
Où [T] s’appella matrice de transformation du système local dans le système

celui global. Aussi,
X 2  X1 Y Y
cos   , sin   2 1 et L X 2  X 1 2  Y2  Y1 2 ,
L L
Avec L est la longueur de l’élément de barre.
55
mécaniques
3.2.2.2 Transformation des forces du système local en système global

(3.1)par
On considère la barre articulée de la figure 2.2 soumise dans le système local
les forces 𝐹𝑥1 et 𝐹𝑥2 appliquées dans les nœuds 1 et 2 dans le système local
d’axes de coordonnées.
Y
y
Fy2, v2 x
𝐹𝑥2 𝑢2
2 Fx2 ,u2
Fx1 , u1 1 
X
𝐹𝑥1 𝑢1 Fy1, v1
Figure3.2 Forces nodales dans le système local et global d’axes de coordonnées

Les composantes de ces forces dans le système global seront :
Fx1  Fx1 cos  , Fy1  Fy1 cos  (3.15)
Fx 2  Fx 2 sin  , Fy 2  Fy 2 sin 
Sous forme matricielle, les relations 3.3 peuvent être écrites sous la forme :
𝐹𝑥1 𝑐𝑜𝑠𝜃 0
𝐹𝑦1
= [ 𝑠𝑖𝑛𝜃 0 ] . {𝐹̅𝑥1 } (3.16)
𝐹𝑥2 0 𝑐𝑜𝑠𝜃 𝐹̅𝑥2
{𝐹𝑦2 } 0 𝑠𝑖𝑛𝜃
3.2.2.3 Détermination de la matrice de rigidité d’un élément barre dans

les coordonnées globales
cos 𝜃 0
sin 𝜃 0 𝐸𝐴 1 −1 cos 𝜃 sin 𝜃 0 0
[𝐾] = [ ]. 𝐿 .[ ].[ ]=
0 cos 𝜃 −1 1 0 0 cos 𝜃 sin 𝜃
0 sin 𝜃
56
mécaniques
𝜆2 𝜆µ −𝜆2 −𝜆µ
2
𝐸𝐴 𝜆µ µ −𝜆µ −µ2
. (3.17)
𝐿
−𝜆 2
−𝜆µ 𝜆2 𝜆µ
[−𝜆µ −µ2 𝜆µ µ2 ]
Ou et 𝜇 sont les cosinus directeurs de la barre = cos 𝜃 et 𝜇 = sin 𝜃

La matrice de rigidité dans le système global d’axes est symétrique, singulière
et, comme on peut facilement remarquer, les éléments de la diagonale principale
sont positifs.
𝑢1
𝑢1  𝜇 00 𝑣1 (3.18)
{ }=[ ]{ }
𝑢0 0 0 𝜇 𝑢2
𝑣2
La définition de l’énergie c’est un demi la multiplication de la force nodale par

le déplacement correspondant :
1
𝑈 = 2 (𝐹𝑥1 𝑢1 + 𝐹𝑥2 𝑢2 )
(3.19)
𝑇
1 𝐹𝑥1 𝑢
= { } { 1}
2 𝐹𝑥2 𝑢2
On peut écrire l’équation de rigidité (3.11) pour une barre orientée
suivant l’axe locale X comme suite :
𝐹𝑥1 𝐸𝐴 1 −1 𝑢1
{ }= 𝐿 [ ]{ } (3.20)
𝐹𝑥2 −1 1 𝑢2
Après substitution des deux équations (3.18) et (3.20) dans l’équation

(3.19)on obtient:
2
𝑢1 𝑇  𝜇 −2 −𝜇 𝑢1
𝐸𝐴 𝑣 𝜇 𝜇2 −𝜇 −𝜇 2 𝑣1
𝑈 = 2𝐿 {𝑢1 } { }
2 𝜇 𝑢2
− 𝜇 (3.21)
2 − 2
𝑣2
[−𝜇 −𝜇 2 𝜇 𝜇 2 ] 𝑣2
Pour obtenir l’équation de rigidité d’une barre orientée arbitrairement on
utilisant la théorème de CASTIGLIANO :
57
mécaniques
𝜕𝑈
𝜕𝑢1
𝐹𝑥1 𝜕𝑈 2 𝜇 −2 −𝜇 𝑢1
𝐹𝑦1 𝜕𝑣1 𝐸𝐴 𝜇 𝜇 2 −𝜇 −𝜇 2 𝑣1
= = 𝐿 {𝑢 } (3.22)
𝐹𝑥2 𝜕𝑈
− 2 − 𝜇  2  𝜇 2
{𝐹𝑦2 } 𝜕𝑢2 2 𝑣
𝜕𝑈
[−𝜇 −𝜇 2
𝜇 𝜇 ] 2
{ 𝜕𝑣2 }
Où symboliquement :
{F}= [K]{q} (3.23)
[𝐾] : est la matrice de rigidité d’un élément barre, qui on peut écrire sous
forme simple :
𝐸𝐴 𝐾0 −𝐾0
[𝐾] = [ ]
𝐿 −𝐾0 𝐾0
2 𝜇
[𝐾0 ] = [ ]
𝜇 𝜇 2
La matrice de rigidité dans le système global d’axes est symétrique, singulière
et, comme on peut facilement remarquer, les éléments de la diagonale principale
sont positifs.
3.3 Détermination de la matrice de rigidité pour un élément poutre

3.3.1 Matrice de rigidité en coordonnées locales
Un élément poutre de longueur L, de section A, de module d’élasticité E et de
moment d’inertie I, est considéré. Aux nœuds 1 et 2 sont appliqués des forces
transversales Yi et des moments de flexion Mi ; avec vi et θi sont les déflexions
transversales et les angles de rotations aux nœuds 1 et 2 respectivement – Fig
3.3
58
mécaniques
Y
M2 , θ2
M1 , θ1
θ 2
1
Y2 , v2 X
Y1 , v1
x
L
Figure 3.3 : Elément poutre en coordonnées locales
La déflexion de l’élément poutre à n’importe quelle position x ; v(x) est donnée

par l’équation différentielle de l’élément non chargé, c'est-à-dire :
 4 v( x)  4 v( x)
EI =0  =0 (3.24)
 x4  x4
Donc : v(x) = a1 + a2 x + a3 x2 + a4 x3 (3.25)

Où ai sont les constantes d’intégrations.
Tandis que θ(x) représente l’angle de rotation ou bien la tangente de v par
rapport à x. autrement dit pour des petites rotations
 v ( x) (3.26)
θ(x) = = a2 + 2 a3 x + 3 a4 x2
x
En appliquant les conditions aux limites à x =0 et x = L, on aura et forme

matricielle :
 v1  1 0 0 0  a1 
  0 0  a 
 1  1 0  2
  =   (3.27)
v 2  1 L L2 L3  a3 
 2  0 1 2 L 3L 
 a 4 
Les constantes d’intégration seront obtenues juste par inversion de l’équation

(3.27) :
59
mécaniques
 a1   L3 0 0 0   v1 
a     
 2 1  0 L3 0 0   1 (3.28)
  =  
 a 3  L  3L  2 L 3L  L2 
3 2
 v 2 
a 4   2 L 2 L   2 
Ou bien : a = T  q (3.29)

Remplaçons ces constantes dans l’équation de forme on aura :
3x 2 2x 2 3x 2 x2 2x3 x3 2x3 x3
v( x ) = v1 + x θ1 - v 1 - θ1 + v2 - θ2 + v 1 + θ1 - v 2 + θ2
L2 L L2 L L3 L2 L3 L2
Ou bien : v( x ) = f1( x ) v1 + f2( x ) θ1 + f3( x ) v2 + f4( x ) θ2 (3.29)
Avec : fi( x ) sont nommées les fonctions de formes données par :
x x
f1 ( x ) = 1 – 3 ( ) 2 + 2 ( ) 3
L L
x2 x3
f2 ( x ) = x – 2 +3 2
L L
(3.30)
x
f3( x ) = 3- 2 ( ) 3
L
x2 x3
f4 ( x ) = – + 2
L L
3 .3.2 Equation de rigidité de l’élément poutre en coordonnées locales
En se basant sur l’énergie de déformation dans la poutre et en utilisant le
théorème de Castigliano ; l’énergie de déformation dans la poutre est donnée
par :
 2v 2
L L
M2 EI
U= 0 2EI dx = 2 (
0 x 2
) dx (3.31)
 2v
Où : 2 = f1"(x) v1 + f2" (x) θ1 + f3" (x) v2 + f4" (x) θ2
x
Avec les dérivées secondes des fonctions de formes sont données par :
x 6
f1"( x ) = 12 3
- 2
L L
4 6x
f2"( x ) = - - 2
L L
60
mécaniques
6 x (3.32)
f3"( x ) = 2
- 12 3
L L
2 x
f4"( x ) = - -6 2
L L
Et les forces aux nœuds sont données :
Fi = kijqj ; i,j =1 – 4
Où qi sont les degrés de liberté et kij sont éléments du tenseur (matrice) de
rigidité ; on aura :
U U  2v   2v
L
EI
F1 = Y1 =
 q1
=
 v1
=
2  2(
0
) (
 x 2 v1 x 2
) dx
L
= EI  ( f1" ( x) v1  f 2 " ( x)1  f 3" ( x) v2  f 4 " ( x) 2 ) f1" ( x) dx (3.33)
0
L
= EI  ( f1" ( x) f1" ( x) v1  f 2 " ( x) f1" ( x)1  f 3" ( x) f1" ( x) v2  f 4 " ( x) f1" ( x) 2 ) dx
0
Utilisant l’équation (3.32) :

F1 = Y1 = k11 v1 + k12θ1 + k13 v2 k14θ2 (3.34)
Par identification des formules (3.32) et (3.34), on peut écrire :
L
k11 = EI  f1" ( x) f1" ( x) dx
0
L
k12 = EI  f1" ( x) f 2" ( x) dx
0
L
k13 = EI  f1" ( x) f 3" ( x) dx
0
L
k14 = EI  f1" ( x) f 4 " ( x) dx
0
De la même façon :
U U U
F2 = M1 = , F3 = Y2 = et F4 = M2 =
 1 v 2  2
En conséquence, on peut déterminer tous les éléments de la matrice de rigidité

en fonction des dérivées secondes des fonctions de formes comme suit :
61
mécaniques
L
kij = EI  f i " ( x) f j " ( x) dx (3.35)
0
Pour i=1 et j=1 on peut trouver k11 :

L
x 6 2 EI
k11 = EI  (12 3
 2 ) dx = 12 3
0 L L L
Même chose pour les autres éléments, on trouve l’équation de rigidité

ainsi que la matrice de rigidité kij qui sont :
 12 6 12 6 
 L2 
L L2 L 
 Y1     v1 
 
M   
 1 EI  6 4 
6
2 
  1
  =  L L  
 Y2  L  
v 2 
  (3.36)
M 2   12 6 12 6  2 
 L2 
L L2
 
L
 
 
 6 6 
 L 2  4 
L 
En conséquence la matrice de rigidité pour l’élément poutre serait:

 12 6 12 6 
 L2 
L L2 L 
 
 
 6 
kij = k  =
EI 6
 L 4  2 
L
L  
  (3.37)
 12 6 12 6
 L2 
L L2
 
L
 
 
 6 6 
 L 2  4 
L 
3.3.3Matrice de rigidité en coordonnées globales
Dans la figure 2.3 est représenté un exemple typique de poutre ainsi que les
deux repères de coordonnées, local xOy et global XOY. Les déplacements dans
le système local sont u, v et, et dans le système global𝑢, 𝑣et 𝜃. L’angle  est
considéré l’angle entre l’axe X et la direction positive de la poutre. Les
déplacements dans le système local xOy peut être exprimé en fonction des
déplacements globaux [14].
62
mécaniques
Figure 3.4 Changement Système local et global de coordonnées
.
Dans ce cas l’axe de la poutre est orienté est orienté suivant (x) comme présenté
sur la figure 3.4 ci-dessous, tout en considérant une force axiale (Xi) dans
chaque nœud et dans la direction x. Dans ce cas l’élément barre est associé avec
l’élément poutre en coordonnées locales. En conséquence la matrice de rigidité
de cet élément peut être obtenu en combinant les deux sous matrices trouvées à
partir de l’élément barre et l’élément poutre trouvés aux paravents.
On part de l’expression de la matrice de rigidité pour un élément de poutre dans

un système local d’axes de coordonnées :
 EA EA 
 L 0 0  0 0 
L
 12 EI 6 EI 12 EI 6 EI 
 0 0  
 L3
L 2
L3 L2 
 6 EI 4 EI 6 EI 2 EI 
 0 
L 
0
K    EA L2 L
EA
L2 (3.38)
 0 0 0 0 
 L L 
 12 EI 6 EI 12 EI 6 EI 
 0  3
 2
0  2 
 L L L3 L 
 0 6 EI 2 EI 6 EI 4 EI 
0 
 L2 L L2 L 
63
mécaniques
3.3.4 Fonction de déplacement
La fonction de déplacement d’un élément poutre, l'équation d'équilibre est

donnée par la formule suivante [3]:
 4v (3.39)
o
x 4
Et le comportement de la poutre est décrit par une fonction polynomiale du
3 degré, qui est la fonction de déplacement et la solution de cette équation [15].:
𝑣(𝑥) = 𝑎1 + 𝑎2 𝑥 + 𝑎3 𝑥 2 + 𝑎4 𝑥 3
𝜕𝑣 (3.40)
𝜃(𝑥) = = 𝑎2 + 2𝑎3 𝑥 + 3𝑎4 𝑥 2
𝜕𝑥
Les 𝑎𝑖 sont les coefficients de déplacement, le vecteur de coefficient de

déplacement {𝑎𝑖 }Peut être obtenu on utilisant les conditions aux limites aux
points nodaux. Les équations des déplacements aux nœuds peuvent être écrites
sous la forme matricielle suivante :
𝑎1
𝑣(𝑥) 1 𝑥 𝑥2
𝑥 ] . { 𝑎2 }
3
{ }=[ 𝑎3 (3.41)
𝜃(𝑥) 0 1 2𝑥 3𝑥 2
𝑎4
Ainsi pour l'élément de poutre de la (figure 2.3)ou les cordonnés des nœuds sont
O et L on aura :
Au nœud (1) x = 0 :
𝑣 = 𝑣1 = 𝑎1
𝜕𝑣
𝜃= = 𝜃1 = 𝑎2 + 2𝑎3 𝑥 + 3𝑎4 𝑥 2 = 𝑎2
𝜕𝑥
Au nœud (2) x = L :
𝑣 = 𝑣2 = 𝑎1 + 𝑎2 𝐿 + 𝑎3 𝐿2 + 𝑎4 𝐿2
𝜕𝑣
𝜃= = 𝜃2 = 𝑎2 + 2𝑎3 𝐿 + 3𝑎4 𝐿2
𝜕𝑥
En exprimant (2.26)sous forme matricielle on aura:
𝑣1 1 0 0 0 𝑎1
𝜃 0 1 0 0 ] . { 𝑎2 }
{𝑣1 } = [ 𝑎3
2 0 𝑙 𝑙2 𝑙3
𝜃2 0 1 2𝑙 3𝑙 2 𝑎4
Ou bien :
64
mécaniques
{𝑞(𝑥)} = [𝐴]{𝑎}
La forme inverse de cette équation donne les expressions des constantes 𝑎𝑖 et

après substitution des coefficients (𝑎𝑖 ) dans l’équation et après réarrangement,
on obtient la forme finale de la fonction de déplacement [7].
2
2 x2 3x 2
V  x   v1  x1  2
3x
v1  1  2 v2 
L L L (3.42)
2 3 3 3
x 3x x 2x x3
 2  3 v1  2 1  3 v2  2  2
L L L L L
V ( x)  f1 ( x)v1  f 2 ( x)1  f 3 ( x)v2  f 4 ( x) 2
On peut avoir les fonctions de forme qui sont [12].
 x
2
x
3
 f 1  x   1  3   2 
 L L
  x   x3 
2
 f 2  x   x  2    2  (3.43)
  L  L 
 2 3
 x x
f
 3  x   3    2 
 L L
  x   x3 
 f 4  x       2 
  L  L 
On passe maintenant à l’application du théorème de Castigliano.

∂U
𝐹𝑖 = (3.44)
∂qi
𝐹𝑖 : Les forces généralisées

𝑄𝑖 : Les déplacements nodaux généralisés
Pour un élément poutre à section droite uniforme, l’énergie de
déformation est donnée par [18].
𝑙 2
𝑀2 𝐸𝐼 𝑙 𝜕 2 𝑣 (3.45)
𝑈=∫ 𝑑𝑥 = ∫ ( 2 ) 𝑑𝑥
0 2𝐸𝐼 2 0 𝜕𝑥
On combinant les deux équations (3.45) et (3.44) on aura:
𝐸𝐼 𝜕 2 𝑣 𝜕 𝜕 2 𝑣 (3.46)
𝐹𝑖 = ( ) ( ) 𝑑𝑥
2 𝜕𝑥 2 𝜕𝑞𝑖 𝜕𝑥 2
65
mécaniques
On dérivant l’équation de déplacement (3.43):
𝜕 2𝑣
= 𝑓1′′ (𝑥)𝑣1 + 𝑓2′′ (𝑥)𝜃1 + 𝑓3′′ (𝑥)𝑣2 + 𝑓4′′ (𝑥)𝜃2
𝜕𝑥 2
On peut calculer maintenant les forces et les moments en fonction des

déplacements. [6].
𝜕𝑈 𝜕𝑈
𝐹𝑦1 = 𝑀𝑙 =
𝜕𝑣1 𝜕𝜃1
𝜕𝑈 𝜕𝑈
𝐹𝑦2 = 𝑀2 =
𝜕𝑣2 𝜕𝜃2
𝐹𝑖 = 𝐾𝑖1 𝑣1 + 𝐾𝑖2 𝜃1 + 𝐾𝑖3 𝑣2 + 𝐾𝑖4 𝜃2
Avec :
𝑙
𝐾𝑖𝑗 = 𝐸𝐼 ∫0 𝑓𝑖′′ (𝑥) + 𝑓𝑗′′ (𝑥) 𝑑𝑥
(3.47)
𝑙
6 𝑥 2 𝐸𝐼
= 𝐸𝐼 ∫ (− 2 + 12 3 ) 𝑑𝑥 = 12 3
0 𝐿 𝐿 𝐿
𝑙
𝐾11 = 𝐸𝐼 ∫ 𝑓1′′ (𝑥) + 𝑓1′′ (𝑥) 𝑑𝑥
0
𝑙
𝐾12 = 𝐸𝐼 ∫ 𝑓1′′ (𝑥) + 𝑓2′′ (𝑥) 𝑑𝑥 (3.48)
0
𝑙
𝐾13 = 𝐸𝐼 ∫ 𝑓1′′ (𝑥) + 𝑓3′′ (𝑥) 𝑑𝑥
0
𝑙
𝐾14 = 𝐸𝐼 ∫ 𝑓1′′ (𝑥) + 𝑓4′′ (𝑥) 𝑑𝑥
0
.
Les 𝑓𝑖′′ (𝑥)𝑒𝑡𝑓𝑗′′ (𝑥)sont les dérivées secondes des fonctions de forme exprimées
dans les équations (2.28),et après toutes les procédures de dérivation et de
substitution, on aboutit à l’équation de rigidité en coordonnées locales :
12 6 12
− 𝐿26
𝐹𝑦1 𝐿2𝐿
6 6 𝐿 𝑣1
𝑀 4 −𝐿 𝜃
𝐸𝐴
{ 1 } = 𝐿 𝐿12 2 {𝑣1 }
𝐹𝑦2 − 𝐿2 − 6
12
−𝐿
6 2 (3.49)
𝑀2 𝐿
𝐿2 𝜃2
−𝐿 4 ]
6 6
[ 𝐿 2
Ou symboliquement :
{𝐹} = [𝐾]{𝑞} (3.50)
[𝐾] :est la matrice de rigidité d’une poutre en coordonnées locales.
66
mécaniques
3.3.5 Transformation des coordonnées locales aux coordonnées globales

Les six forces nodales et moments de flexion dans un système de coordonnées
locales sont transformés en coordonnées globales, par l’équation suivante :
 FX 1   cos  sin  0 0 0 0   FX 1 
   sin   
 FX 2   cos  0 0 0 0  FX 2 
F 
 Y1   0 0 1 0 0 0  FY 1 

     (3.51)
M 1   0 0 0 cos  sin  0  M 1 
F   0 0 0  sin  cos  0  FY 2 
 Y2     
  0 1 
M2 
M 2  
0 0 0 0

L’équation (3.51)peut être sous forme symbolique :

{F} = [T]{F} (3.52)
On peut aussi transformer les six degrés de libertés du système de

coordonnées locales au système de coordonnées globales:
{q} = [T]{q}
La transposée de la matrice de transformation [T] est la même avec sa matrice

inverse:
[T]−1 = T T
On peut écrire l’équation (3.52)sous la forme inverse:

{F} = [T]−1 {F} = [T]T {F}
On substituant cette équation dans équation force déplacement (3.52),

on obtient:
{F} = [K]{q}
[T]{F} = [K][T]{q}
Comme la matrice [T] est orthogonale, on peut alors déduire que :

{F} = [T]T [K][T]{q} (3.53)
Après toutes les substitutions pour la matrice [K] et la matrice de transformation
67
mécaniques
[T] dans l’équation (3.53),on obtient la forme finale de l’équation de rigidité
pour un élément poutre dans système de coordonnées globales.
12
R2 + L2 μ2
Fx1 12
(R − L2 ) μ
12
Rμ2 + L2 2
Fx2 6 6
Fy1 −Lμ 
EI
= L L 4 12
M1 2 12 2 12
−R − L2 μ (−R + 2 ) μ
6
μ R2 + L2 μ2 (3.54)
Fy2 L L
12
12
Rμ2 + 2 2
12 12 2
(−R + L2 ) μ −Rμ − 2  − L 
2 6 (R − ) μ L
{ M2 } L
L2
6
6 6 2
6
μ −L 4
[ −Lμ  L ]
L
3.3.6 Détermination de la force axiale

La force axiale dans un élément poutre ou barre orientée arbitrairement dans
le plan d’un angle 𝜃 avec l’horizontal est donnée par l’équation (3.22)
EA 2
Fx2 = [ (μ2 − μ1 ) + μ(v2 − v1 )]
L
EA
Fy2 = [μ(μ2 − μ1 ) + μ2 (v2 − v1 )]
L
Et la force axiale défini comme suite:
S = Fx2  + Fy2 μ
EI 2
= ( + μ2 )[(μ2 − μ1 ) + μ(v2 − v1 )]
L
EI
= [(μ2 − μ1 ) + μ(v2 − v1 )]
L
Ou bien sous forme matricielle :
EI μj − ui
Sn = [n μn ] { v − v } (3.55)
L j i
3.4 Détermination de la matrice de rigidité pour un élément plaque

(élément rectangulaire à 4 nœuds)
3.4.1 Définition d’une plaque mince
Une plaque est un solide délimité par deux plans parallèles , tel que le plan
moyen, ou plan médian situé à équidistance entre les faces (c'est l'équivalent de
la courbe moyenne des poutres) , la plaque est mince dont l’épaisseur est petit
68
mécaniques
comparée par les deux autres dimensions .
3.4.2 Fonctions de déplacement
L’élément rectangulaire à 4 nœuds est un quadrilatère (rectangle) à quatre
nœuds. Ses
degrés de liberté sont présentés dans la (figure 3.5)
y, v
a a
v3 v4
3 4
u4
u3
b
x, u
u1 b
u2
1 v1 2 v2
Figure3.5 Elément rectangulaire à 4 nœuds

En termes de coordonnées généralisées i, le champ de déplacement est donné
par:
u  1   2 x   3 y   4 xy
(3.56)
v   5   6 x   7 y   8 xy
Pour représenter la déformation d’un élément rectangulaire à n’importe quel

point défini par (x, y) sur la plaque, on choisit une simple fonction polynomiale
pour définir les fonctions de déplacement comme suite.
u ( x, y)  c1  c 2 x  c3 y  c 4 xy
v( x, y)  c1  c 2 x  c3 y  c 4 xy
Le terme “bilinéaire” vient du fait que la forme des expressions pour u et v est le
produit de deux polynômes linéaires (c1 + c2x)(c3 + c4y) où ci sont des
constantes. Il y déformation sera:
69
mécaniques
u
x   2  4 y
x
v
 yz   7  8 x (3.57)
y
u v
 xy    3  6   4 x  8 y
y x
Le champ de déformation montre que la quantité x est indépendante de x, ce qui

veut dire que l’élément rectangulaire à 4 nœuds ne peut pas modeler de façon
exacte une poutre encastrée sollicitée par une force transversale, là où la
déformation axiale varie de façon linéaire par rapport à l’axe Ox. (figure. 3.6).
En plus, l’élément rectangulaire à 4 nœuds ne peut pas modeler un état de
flexion pure, malgré son habileté de représenter une déformation x qui varie
linéairement par rapport à l’axe Oy.
F
M M M M
(a) (b) (c)

Figure.3.6 a) Simulation de la flexion. b) Mode de déformation correcte en
flexion pure. c) Mode de déformation de l’élément bilinéaire quadratique
3.4.3 Fonctions de forme

Tous les côtés de l’élément rectangulaire à 4 nœuds se déforment sous
forme de lignes droites. En conséquence les angles droits dans l’intérieur de
l’élément ne restent pas droits suite à la déformation, donc les déformations de
cisaillement (glissements) apparaissent partout dans l’intérieur de l’élément
exceptant au long de l’ace Oy. Comme on le sait bien, dans la flexion pure xy =
0; en conclusion, l’élément rectangulaire à 4 nœuds utilisé dans la modélisation
de la flexion, introduit un effet parasite de déformation de cisaillement.
70
mécaniques
Les fonctions de forme ont les expressions suivantes :
 x  y
f 1 ( x, y )  1  1  
 a  b
x y
f 2 ( x , y )  1  
a b (3.58)
xy
f 3 ( x, y ) 
ab
 x y
f 4 ( x , y )  1  
 a b
y
4 3
x
1
2
1
Figure.3.7 Fonction de forme f2 de l’élément bilinéaire quadrilatéral

Le champ de déformation sera :
 u1 
v 
 
1
 1 y 1 y y y
  a 1  b  0 1   0 0  0  u2 

 x     a b ab ab  v(3.59)

   1 x x x 1  x   
 2
 y    0  1   0  0 0 1   u 3 
   b a ab ab b  a  v 
 xy   1  x  y   
3
1 y x 1 y x y 1 x
 1     1    1   1    u
 b  a  a b ab a  b  ab ab b a ab   4 
v 4 
 
 
Ou bien sous forme matricielle :

{𝜀} = [𝐴]{𝑞}
3.4.4 Elément quadrilatéral quadratique (élément rectangulaire à 8 nœuds)

L’élément élément rectangulaire à 8 nœuds est présenté dans la (figure 3.7). En
terme de coordonnées généralisées i est donné par:
71
mécaniques
u   1   2 x   3 y   4 x 2   5 xy   6 y 2   7 x 2 y   8 xy 2
(3.60)
v   9   10 x   11 y   12 x 2   13 xy   14 y 2   15 x 2 y   16 xy 2
y, v
a a
v3 v4
3 7 4
u4
u3
b
8 6 x, u
u1 b
u2
1 v1 v2 2
Figure.3.8 Elément quadratique quadrilatéral
y y
4 3 4 3
6
x x
1 1 2
2
1 1
Figure.3.9 Fonctions de forme

3.4.5 Détermination de la matrice de rigidité d’un élément plaque par la
méthode de l’énergie de déformation
On utilisant l’énergie de déformation d’un élément plaque pour déterminer la

matrice de rigidité élémentaire, tel que l’énergie de déformation par rapport à
chaque degré de liberté donnée par le théorème de Castigliano comme suit.
a b t
U        dz d y dx
1 T
(3.61)
2000
D’après le théorème de Castigliano, les forces égale la dérivée de l’énergie par

rapport les déplacements nodaux, en supposant que l’épaisseur e est constante,
donc la formule (3.61)devienne comme suite :
72
mécaniques
e T 
b a
U q   ACAdx dy q (3.62)
2  0 0 
L’équation (3.61)peut être résolue ainsi:

U
Fi 
qi
Après avoir résolu le système, on obtient la forme matricielle suivante :
F   K q
où [k] est maintenant la matrice qui donne l’expression de la matrice de rigidité
d’un élément plaque,
K   e  AT C A dx dy (3.63)

aire
Les coefficients de rigidité pour un élément en état plaque peuvent être obtenus
en utilisant l'équation de l'énergie de déformation (3.62),les fonctions de
déplacement, et le premier théorème de Castigliano:
K   
K11 K 21  (3.64)
 K12 K 22 
En partant de l’équation (3.63)on aura :
 f i  f j  1   f i  f j 
K11 ij  E
  x  x      e dA
1  2 A  2  y  y 
 f i  f j  1   f i  f j 
K 22 ij  E
  y  y      e dA (3.65)
1  2 A  2  x  x 
  f i  f j 1   f i  f j 
K 21 ij  E
   x  y     e dA
1  2 A  2  y  x 
Donc la relation constitutive pour un élément plaque entre les forces et les
déplacements peuvent être écrits sous la forme:
73
mécaniques
𝐶1
𝐹𝑥1
𝐶2 𝐶3 𝑢1
𝐹𝑦1 𝑣1
𝐶4 −𝐶5 𝐶1
𝐹𝑥2 𝑢2
𝐹𝑦2 𝐶5 𝐶6 −𝐶2 𝐶3
𝑣2 (3.66)
= −𝐶1 −𝐶2 𝐶7 −𝐶5 𝐶1 𝑢3
𝐹𝑥3
𝐹𝑦3
2 −𝐶3 𝐶5 𝐶8 𝐶2 𝐶3 𝑣3
−𝐶2 2 −𝐶1 𝑢4
𝐹𝑥4 𝐶7 𝐶5 𝐶2 𝐶4 −𝐶5 𝐶1
2 { 𝑣4 }
𝐹
{ 𝑦4 } −𝐶 𝐶8 −𝐶3
−𝐶2 𝐶3
[ 5 𝐶2 𝐶5 𝐶6 ]
2
Alors on peut avoir la matrice de rigidité, avec les coefficients de rigidité Kij
 b 1  a  E e
C1    
 3a 6 b  1  2
  1   E e
C2    
4 8  1  2
 a 1  b  E e
C3    
 3b 6 a  1  2
 b 1  a  E e
C4     
 3a 12 b  1 
2
(3.67)
  1   Ee
C5    
4  1 
2
8
 a 1  b  E e
C6    
 6b 6 a  1  2
 b 1  a  E e
C7    
 6a 6 b  1  2
 a 1  b  E e
C8     
 3b 12 a  1 
2
74
Chapitre 4 :
Mise en équation de l’analyse

dynamique des structures
mécaniques par la MEF
75
4.1 Introduction
Toutes les structures mécaniques soumises à des forces extérieures variables
sont susceptibles d’entrer en vibration. Ces forces d’excitation fonction du
temps sont périodiques, impulsionnelles ou aléatoires [23]. Les vibrations des
systèmes mécaniques ont le plus souvent des effets nocifs qui se ramènent
pratiquement à trois types :
- Les résonances entrainant des bruits excessifs ou des ruptures par amplitude
de déformation excessive ;
- La fatigue entrainant la rupture des matériaux ou des détériorations
d’équipements ;
- Et enfin les instabilités de toutes natures, (instabilités de combustion, de
pilotage, d’asservissement, instabilités par flambage dynamique, etc.).
Dans tous les cas, des calculs pour décrire le comportement dynamique des
structures mécaniques on commencera par exposer la théorie des structures
discrétisées et continu pour déterminer la matrice de rigidité pour une analyse
statique. Puis on donnera quelques théorèmes utiles pour déterminer les modes
propres pour une vibration libre. On terminera par quelques formules donnant la
réponse d’une structure linéaires à des excitations forcées, pour une vibration
forcée.
Les structures mécaniques font l’objet de trois types de calculs classiques
destinés à prévoir leur comportement vibratoire :
- Le calcul des vibrations libres des structures linéaires conservatives. En effet
les forces d’amortissement, le plus souvent inaccessibles au calcul, sont en
général estimées ou mesurées dans leur aspect généralisé modal, puis
introduites après coup.
- Le calcul des réponses des structures dissipatives passives à des excitations de
type périodique, impulsionnelles, ou aléatoire.
76
- Le calcul des structures actives, afin de déterminer les domaines d’instabilité

des paramètres qui définissent les structures et les actions extérieures.
4.2 Analyse dynamique d’un élément barre
Considérons une structure élastique linéaire conservative occupant un domaine
𝑜 de frontière o dans l’état naturel, et supposons la soumise à une densité
volumique f de force donnée et à une densité surfacique F de force donnée sur
o. si les forces dépendent du temps le déplacement U résultant en un point
𝑀𝑂 𝑑𝑒 ̅̅̅̅
𝑂 dépendra de MO et du temps t [21] :
𝑈 = 𝑈(𝑀𝑂 , 𝑡), 𝑓 = 𝑓(𝑀𝑂 , 𝑡), 𝐹 = 𝐹(𝑀𝑂 , 𝑡) (4.1)
MO et t étant les variables de Lagrange.

La structure est alors soumise à une densité volumique de forces d’inertie et si
on appelle :
𝜌 = 𝜌(𝑀𝑂 ) (4.2)
la masse spécifique au point MO, cette densité s’écrit −𝜌𝑈̈à un instant t

quelconque, ou
𝜕 2𝑈
𝑈̈ = (4.3)
𝜕𝑡 2
La discrétisation, par exemple par la méthode des éléments finis en

déplacement, s’obtient dans chaque élément e par l’interpolation :
𝑈 = 𝑇(𝑀𝑂 ). 𝑞(𝑡) (4.4)
Ou q(t) est la colonne des déplacements (généralisés) nodaux.

La densité volumique de forces d’inertie s’écrit donc :
−𝜌𝑈̈ = −𝜌(𝑀𝑂 )𝑇(𝑀𝑂 )𝑞̈ (𝑡) (4.5)
Et donne un travail virtuel :
̅̈ 𝛿𝑈𝑑 = −𝛿𝑞 ̅̅̅̅̅̅̅̅
̅̅̅ [∫ 𝜌(𝑀𝑂 )𝑇(𝑀 ̅̅̅ 𝑀𝑒 = ̅̅̅̅
∫ − 𝜌𝑈  𝑂) 𝑇(𝑀𝑂 )𝑑] 𝑞̈ (𝑡) = −𝛿𝑞 𝑀𝑒 𝑞̈ (𝑡), 𝑀𝑒 (4.6)
𝑒 𝑒
Me est la matrice de masse de l’élément e rapportée aux déplacements nodaux q.

Finalement, après assemblage, l’équation discrétisée de la structure s’écrira:
77
[
mq  k q  F (t ) (4.7)
̅,
𝑎𝑣𝑒𝑐 𝑀 = 𝑀 ̅
𝐾=𝐾
Les matrices M et K sont les matrices de masse et de rigidité de la structure

globale discrétisée
Alors la force de déplacement d’après la loi de Hooke donnée par la formule
suivante [7]:
𝑑𝑁 𝜕2 𝑈
= 𝜌𝐴 (4.8)
𝑑𝑥 𝜕𝑡 2
Après discrétisation l’énergie de déformation du système s’écrit:

L’énergie potentielle:
1 1 𝑙 𝜕𝑈
𝐸𝑝 = ∫𝛺 𝜎𝜀𝑑𝛺 = ∫0 𝐸𝐴 ( )2 𝑑𝑥 (4.9)
2 2 𝜕𝑥
1 1 𝑙 𝜕𝑈
𝐸𝑝 = ∫𝛺 𝜎 = ∫0 𝐸𝐴 ( )2 𝑑𝑥 (4.10)
2 2 𝜕𝑥
Et l’énergie cinétique :
2 2
1 𝜕 𝑈 2 𝑙 𝜕 𝑈 2 (4.11)
𝐸𝐶 = 2 ∫𝛺 𝜌( ) 𝑑𝛺 = ∫0 𝜌𝐴( ) 𝑑𝑥
𝜕𝑡2 𝜕𝑡2
Les deux formes quadratiques Ep et Ec sont semi-positives dans le cas général.

Le principe d’Hamilton donne classiquement l’équation de Lagrange :
𝑑 𝜕𝐸𝐶 𝜕𝐸𝑝 U (4.12)
− = Fq  -
𝑑𝑡 𝜕𝑞̇ 𝜕𝑞 q
Les déplacements généralisées 𝑞 , au nombre de 6, sont les modes de corps

solide, ou déplacement d’ensemble. Si on s’intéresse aux mouvements libres de
la structure, 𝐹(𝑡) = 0 et l’équation (1), après élimination des réactions d’appuis,
devient:
mq k q  0 (4.13)
78
Dont les solutions non nulles sont harmoniques du type :

𝑞(𝑡) = 𝑞𝑒 𝑗𝑤𝑡
L’équation (3.13) devient alors l’équation aux vibrations libre classiques
d’un système linéaire conservatif discrétisé :
(4.14)
[−𝑤 2 𝑀 + 𝐾]𝑞 = 0
Lorsque les déplacements d’ensemble ont été éliminés, et si les
déplacements q sont indépendants, les deux matrices M et K sont définies
positives, de même que les formes quadratiques correspondantes Ec et Ep.
Les déplacements et efforts imposées nuls aux extrémités :
𝑈(0, 𝑡) = 0 𝑒𝑡 𝑈(𝑙, 𝑡) = 0
(4.15)
𝜕𝑈 𝜕𝑈
(0, 𝑡) = 0 𝑒𝑡 (𝑙, 𝑡) = 0
𝜕𝑥 𝜕𝑥
Apres une séparation des variables on peut obtenir:

𝑈(𝑥) = 𝑀(𝑥). 𝑇(𝑥) (4.16)
Pour assurer la stabilité de la solution, les deux fonctions sont donc égales a une
constante :
1 𝑑2𝑀 1 𝑑2𝑇 𝜔2
= =− (4.17)
𝑀 𝑑𝑥 2 𝐶 2 𝑑𝑡 2 𝐶2
On a donc :
𝑑2𝑀 𝜔 2
2
+( ) 𝑀 =0 (4.18)
{ 𝑑𝑥 2 𝐶
𝑑 𝑇
+ 𝜔2 𝑇 = 0
𝑑𝑡 2
𝜔𝑥 𝜔𝑥
𝑀(𝑥) = 𝐴. 𝑠𝑖𝑛 + 𝐵. 𝑐𝑜𝑠
{ 𝑐 𝑐
(4.19)
𝑇(𝑡) = 𝐶. 𝑠𝑖𝑛𝜔𝑡 + 𝐷. 𝑐𝑜𝑠𝜔𝑡
Tel que les constantes A, B, C et D calculées d’après les conditions initiales
et les conditions aux limites. Alors une vibration linéaire d’un élément barre
s’écrit par un système d’équation différentielle de la forme suivante:
𝑀𝑋̈ + 𝐾𝑋 = 𝐹(𝑡)
79
Fx1 u1 Fx2 u2
A B
e
x dx
L
Figure 4.1. les forces axiales pour un élément barre

L’équation forces-déplacements d’un élément barre s’écrit donc comme suite :
𝜌𝐴𝐿
𝐹𝑥1 − 𝑢̈ 1 𝐸𝐴 1 −1 𝑢1
2
{ }= 𝐿 [ ]{ }
𝜌𝐴𝐿
𝐹𝑥2 − 2 𝑢̈ 2 −1 1 𝑢2 (4.20)
Ou bien sous la forme suivante :

𝐹𝑥1 𝐸𝐴 1 −1 𝑢1 𝜌𝐴𝐿 1 0 𝑢̈ 1
{ }= 𝐿 [ ] {𝑢 } + 2 [ ]{ } (4.21)
𝐹𝑥2 −1 1 2 0 1 𝑢̈ 2
Tel que la déplacement est sinusoïdal de la forme suivante :
𝑢 = 𝑎 sin 𝜔𝑡
Donc :
𝐹𝑥1 𝐸𝐴 1 −1 𝜌𝐴𝐿 1 0 𝑢1
{ }=[𝐿 [ ] − 𝜔2 2 [ ]] {𝑢 }
𝐹𝑥2 −1 1 0 1 2 (4.22)
4.3 Détermination de la matrice masse d’un élément barre

Pour déterminer la matrice masse d’un élément barre on utilisant l’équation de
Lagrange à partir de [26]:
𝑑 𝜕𝐸 𝜕𝐸𝑃
( 𝑐) −
𝑑𝑡 𝜕𝑞̇ 𝜕𝑞𝑖
= 𝐹𝑥𝑖 (4.23)
𝑖
𝑚𝑖 𝑞̈ 𝑖
Avec : 𝐸𝑐 = , et 𝑚𝑖 𝑞̈ 𝑖 = 𝐹𝑥𝑖
2
𝑑 𝜌𝐴 𝑙 𝜕 𝐸𝐴 𝑙 𝜕
𝐹𝑥1 = 𝑑𝑡 ( 2 ) ∫0 2{𝑢̇ (𝑥)} 𝜕 𝑢̇ (𝑢̇ (𝑥)𝑑𝑥 + ∫0 𝑢(𝑥) 𝜕𝑢 (𝑢(𝑥))𝑑𝑥 (4.24)
1 2 1
𝑙 𝑙
𝑑
= (𝜌𝐴) ∫ (𝑓1 (𝑥)𝑢̇ 1 + 𝑓2 (𝑥)𝑢̇ 2 )(𝑓1 (𝑥)𝑑𝑥 + 𝐸𝐴 ∫ (𝑓1′ (𝑥)𝑢1 + 𝑓2′ (𝑥)𝑢2 )(𝑓1′ (𝑥)𝑑𝑥
𝑑𝑡 0 0
𝑙 𝑙
𝑢̈
= [𝜌𝐴 ∫ 𝑓1 (𝑥)𝑓1 (𝑥)𝑑𝑥 𝜌𝐴 ∫ 𝑓1 (𝑥)𝑓2 (𝑥)𝑑𝑥] { 1 }
0 0
𝑢̈ 2
80
𝑙 𝑙 𝑢1
+ [𝐸𝐴 ∫0 𝑓1′ (𝑥)𝑓1′ (𝑥)𝑑𝑥 𝐸𝐴 ∫0 𝑓1′ (𝑥)𝑓2′ (𝑥)𝑑𝑥] {𝑢 } (4.25)
2
Sous forme implicite :

𝑢
𝐹𝑥2 = [𝑚21 𝑚22 ] {𝑢̈ 1 } + [𝑘21 𝑘22 ] {𝑢1 } (4.26)
𝑢̈ 2 2
Donc la matrice de rigidité et de masse pour un élément barre soumis des forces
vibratoire s’écrit sous une forme générale comme suite:
𝑙
𝑚𝑖𝑗 = 𝜌𝐴 ∫0 𝑓𝑖 (𝑥)𝑓𝑗 (𝑥)𝑑𝑥
𝑙
(4.27)
𝑘𝑖𝑗 = 𝐸𝐴 ∫ 𝑓𝑖′ (𝑥)𝑓𝑗′ (𝑥)𝑑𝑥
0
4.4 Vibration d’un élément poutre

La solution d’équilibre statique d’une structure élastique linéaire en
élément poutre rendait minimale son énergie potentielle totale dans la classe des
déplacements cinématiquement admissibles, soit une structure d’un élément
poutre occupant un domaine O , de frontière O, soumise à des forces
extérieures de densité volumique , tel que v(x,t) est la flèche de la section A
𝜕𝑣
d’abscisse x a l’instant t, avec un angle (𝜃 = ) entre la section A et l’axe oy ,
𝜕𝑥
Ou 𝑀𝑧 (𝑥, 𝑡) est le moment fléchissant en un point quelconque de la poutre,

causé par l’effet combiné du système de forces extérieures. L’intégration est
faite sur la longueur totale de la poutre.
D’après le théorème de castigliano, on considère la poutre montrée sur la figure
4.2 , cette poutre a fléchi sous l’action de la force Fi qui est appliquée
graduellement. Le potentiel U acquis par la poutre est égal au travail externe fait
par la force Fi .
81
Y F
M1 , θ1, u1, v1 M2 , θ2, u2, v2

X
A B
Figure.4.2 élément poutre en flexion
Nous avons les équations d’équilibre de la poutre, tel que les forces sont
parallèles a l’axe oy , et les moments autour de l’axe oz .
L’équation d’équilibre pour les forces [22]:
𝜕𝑇 𝜕2 𝑉
= 𝜌𝐴 (4.28)
𝜕𝑥 𝜕𝑡 2
L’équation d’équilibre pour les moments[22] :
𝜕𝑀 𝜕2𝜃 (4.28)
𝑇+ = 𝜌𝐼 2
𝜕𝑥 𝜕𝑡
Avec la loi de comportement dynamique d’une structure en élément poutre en
flexion :
𝑑𝑣
𝑑𝜃 2
𝑑 𝑉
𝜃 =
𝑑𝑥
𝑀 = 𝐸𝐼 = 𝐸𝐼 2 et[ 𝑀 𝜕 2𝑉
(4.29)
𝑑𝑥 𝜕𝑥
= 2
𝐸𝐼 𝜕𝑥
Avec I est le moment quadratique de flexion :

𝐼 = 𝐴𝑖𝑥2
On peut aussi déterminer L’équation de mouvement appliquée sur la poutre, On
a:
𝜕2 𝑉 𝑑4𝑉 𝑑4𝑉
𝜌𝐴 2 + 𝐸𝐼 4 = 𝜌. 𝐼 2 2 (4.30)
𝜕𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑡
𝜕4 𝑉
La solution de cette équation dépend duYcouple de l’inertie(−𝜌𝐼. 𝑑𝑥 ) qui
𝜕𝑥 2 𝜕𝑡 2
𝜕𝑀
est négligeable devant l’effort tranchant (𝑇 + ) 𝑑𝑥 , cette expression devient
𝜕𝑥
X
égale à :
82
𝜕2 𝑉 𝜕4 𝑉
𝜌𝐴 + 𝐸𝐼 =0 (4.31)
𝜕𝑡 2 𝜕𝑥 4
Après une substitution des paramètres d’un côté on peut obtenir les suivant:
𝜆 𝜆 𝜆 𝜆
𝑉(𝑥) = 𝑎. 𝑐𝑜𝑠ℎ 𝑥 + 𝑏. 𝑠𝑖𝑛ℎ 𝑥 + 𝑐. 𝑐𝑜𝑠 𝑥 + 𝑑. 𝑠𝑖𝑛 𝑥
𝑙 𝑙 𝑙 𝑙
𝜆 𝜆 𝜆 𝜆 𝜆 𝜆 𝜆 𝜆 (4.32)
𝜃(𝑥) = −𝑎 𝑠𝑖𝑛ℎ 𝑥 + 𝑏 𝑐𝑜𝑠ℎ 𝑥 − 𝑐 𝑠𝑖𝑛 𝑥 + 𝑑 𝑐𝑜𝑠
𝑙 𝑙 𝑙 𝑙 𝑙 𝑙 𝑙 𝑙
𝜆2 𝜆 𝜆 𝜆 𝜆
𝑀(𝑥) = 𝐸𝐼 2 (𝑎. 𝑐𝑜𝑠ℎ 𝑥 + 𝑏. 𝑠𝑖𝑛ℎ 𝑥 + 𝑐. 𝑐𝑜𝑠 𝑥 + 𝑑. 𝑠𝑖𝑛 𝑥)
𝑙 𝑙 𝑙 𝑙 𝑙
3
𝜆 𝜆 𝜆 𝜆 𝜆
𝑇(𝑥) = −𝐸𝐼 3 (𝑎. 𝑠𝑖𝑛ℎ 𝑥 − 𝑏. 𝑐𝑜𝑠ℎ 𝑥 + 𝑐. 𝑠𝑖𝑛 𝑥 − 𝑑. 𝑐𝑜𝑠 𝑥)
𝑙 𝑙 𝑙 𝑙 𝑙
Peut alors écrire ces équations sous forme matricielle :
1 0 1 0
𝑉(0) 0
𝜆
0
𝜆 𝑎
𝑙 𝑙 (4.33)
𝜃(0) 𝑏
{ }= 𝜆2 { }
𝑀(0) 𝜆2 0 𝐸𝐼 2 0 𝑐
𝐸𝐼 2 𝜆3
𝑙
𝑇(0) 𝑙
𝐸𝐼 𝜆3 𝑑
[ 0 𝑙3 0 𝐸𝐼 3 ]
𝑙
Après quelques transformations, et l’utilisation de la pulsation ω, ainsi, la

solution harmonique, on peut écrire la solution générale comme suite :
𝜋2 𝐸𝐼
𝜔𝑛 = 𝑛2 √𝜌.𝑠 (4.34)
𝑙2
𝑥
𝑉𝑛 (𝑥, 𝑡) = sin 𝑛𝜋 (𝛼𝑛 cos 𝜔𝑡 + 𝛽 sin 𝜔𝑡)
𝑙
On utilisant l’énergie cinétique d’un élément poutre en vibration de flexion:

𝜌𝐴 1
𝐸𝑐 =
2
∫0 𝑣̇ (𝑥)2 𝑑𝑥 (4.35)
Et l’énergie potentielle :
𝐸𝐼 𝐿
𝐸𝑝 =
2
∫0 𝑣 ′′ (𝑥)2 𝑑𝑥 (4.36)
On utilisant la méthode de Lagrange pour dériver l’énergie potentielle par

rapport les déplacements, on trouvant la forme suivante :
mq k q  F (t ) (4.37)
83
Les équations de mouvements résultantes sont de la forme explicite suivante :

u1 
 Fx1   12 2   
   L2    
   12  v1 
 Fy1  12 2  
    
   L2 L2   
M   6 6   
 EI   L  
1 
 
4   1
L
     
F  L   12  2 12

6
 
12 2
  u 
 x2
  L2 L2 L L2   2
   12 12 6 12 12 2   
F   L2   2 2
L

L
 
L
 
L2
    (4.38)
 y 2   6 6 6 6  v 2 
      2    4  
M 2 
  
 L L L L 
  
 2 
 
 156 2  u1 
   
 156 156 2
 v1 
AL   22 L 22 L 4 L2  1 
2    
420  54 2  54  13L 156 2  u 2 
  54 542 13L  156 1562  v 
   2
 13L  13L  3L2 22 L  22 L 4 L2   2 
4.5 Détermination des modes propres pour un élément poutre

Le système d’équations différentielles d’une vibration libre pour un élément
poutre, sous forme matricielle suivante :
(4.39)
mq k q  0
Ainsi que, les déplacements nodaux sont sinusoïdal de la forme :
q( x)  a sin t (4.40)
On obtient comme précédemment lorsque on substituant l’équation (4.39) dans
l’équation (4.40)
k    mq( x)  0
2
Puis, on multipliant cette dernière par [𝑚]−1 , on obtient :
k m
1 1

  2 mm .q( x)  0 Avec [𝐴] = [𝑘][𝑚]−1 (4.41)
La réponse du système d’équation pour trouver la solution doit le déterminant

est nul :
A  I  0
La solution prend la forme :
(4.42)
84
𝑎11 − 𝜆 𝑎12 … . . 𝑎1𝑛

| 𝑎21 𝑎22 − 𝜆 … . . 𝑎2𝑛 | = 0
𝑎𝑛1 𝑎𝑛2 𝑎𝑛𝑛 − 𝜆
Soit, après développement :
𝜆𝑛 + 𝑎1 𝜆𝑛−1 + 𝑎2 𝜆𝑛−2 … … … 𝑎𝑛−1 𝜆 + 𝑎𝑛 = 0 (4.43)
Cette équation (3.43) exprime les pulsations propres d’une poutre, tels que les
scalaires 𝜆 = 𝜔2 sont les valeurs propres de la matrice[𝐴].
4.6 vibration d’un élément plaque

le système d’équations différentielles d’une vibration dissipatif pour un élément
plaque, et l’équation de mouvement pour un oscillateur harmonique amorti
soumis à une force extérieure F(t) s’écrit sous la forme matricielle suivante [16]:
mxt   cxt   kxt   F (t )
c k
F(t)
Figure.4.3 système à un degré de liberté
Ce système est constitué de n équation indépendantes, l’équation caractéristique

associée à cette équation du mouvement [16] :
(4.44)
𝑚𝑟 2 + 𝑐𝑟 + 𝑘 = 0
La solution générale du système s’obtient par la résolution de l’équation

caractéristique, qui devient comme suite:
c c 2  4km
r1, 2    (4.45)
2m 2m
85
Ou, d’une autre manière :

𝑟1,2 = −𝛼𝑤0 ± 𝑤√𝛼 2 − 1 (4.46)
Après une substitution des paramètres, on peut écrire l’équation du mouvement
sous la forme :
x  2w0x  w0 2 x  0 (4.47)
Et cherchons maintenant la solution du type :
𝑥 = 𝐴𝑒 𝑟1𝑡 + 𝐵𝑒 𝑟2𝑡 (4.48)
Tel que les constantes A et B sont déterminés par les conditions initiales.
La solution générale s’écrit donc comme suite :
𝑥(𝑡) = (𝑐1 + 𝑐2 𝑡)𝑒 −𝛼𝑤0𝑡 (4.49)
Energie d’une vibration pour un élément plaque :

L’énergie totale d’une vibration pour un élément poutre d’après le théorème de
Lagrange est la somme de l’énergie potentielle et l’énergie cinétique [24] :
l’énergie potentielle :
𝑛 𝑛
1 (4.50)
𝐸𝑝 = ∑ ∑ 𝐾𝑖𝑗 𝑥𝑖 𝑥𝑗
2
𝑖=1 𝑗=1
Avec 𝐾𝑖𝑗 = 𝐾𝑗𝑖
- l’énergie cinétique :
1
𝐸𝑐 = ∑𝑛𝑖=1 ∑𝑛𝑗=1 𝑚𝑖𝑗 𝑥𝑖° 𝑥𝑗° avec 𝑚𝑖𝑗 = 𝑚𝑗𝑖
2
(4.51)
𝑥𝑖° = 𝑑𝑥𝑖 ⁄𝑑𝑡
- l’énergie dissipative:
1
𝐸𝑑 = ∑𝑛𝑖=1 ∑𝑛𝑗=1 𝐶𝑖𝑗 𝑥𝑖° 𝑥𝑗° Avec:𝐶𝑖𝑗 = 𝐶𝑗𝑖 (4.52)
2
Par définition les équations de Lagrange s’écrivent :
86
𝑑 𝜕𝐸 𝜕𝐸 𝜕𝐸𝑑 𝜕𝐸𝑐
( 𝑐 ) − 𝜕𝑥𝑐 +
𝑑𝑡 𝜕𝑥 ° 𝜕𝑥𝑖°
+
𝜕𝑥𝑖
(4.53)
𝑖 𝑖
Tel que les fi représentent les forces généralisées, on dérivant l’énergie potentiel
𝐸𝑝 , pour obtenir le travail global des forces extérieures s’écrit :
𝑑𝑤 = ∑ 𝑓𝑖 𝑑𝑥𝑖
𝑖=1
Pour une vibration libre, nous avons :

𝑑 𝜕𝐸 𝜕𝐸 𝜕𝐸𝑝
( 𝑐 ) − 𝜕𝑥𝑐 +
𝑑𝑡 𝜕𝑥̇ 𝜕𝑥
=0 (4.54)
𝑖 𝑖
Le plus souvent, Ec est une forme quadratique des 𝑥̇ 𝑖 et :

𝑑 𝜕𝐸 𝜕𝐸𝑝
( 𝑐) +
𝑑𝑡 𝜕𝑥̇ 𝜕𝑥
=0 (4.55)
𝑖
Pour un système libre dissipatif :

𝑑 𝜕𝐸 𝜕𝐸𝑑 𝜕𝐸𝑝
( 𝑐) −
𝑑𝑡 𝜕𝑥̇ 𝑥̇ 𝑖
+
𝜕𝑥
=0 (4.56)
𝑖
Détermination de l’équation caractéristique du structure:
nous avons n équations linéaires caractérisant l’état vibratoire du système, et

après une utilisation des équation de Lagrange, on obtient :
∑𝑛𝑗=1(𝑚𝑖𝑗 𝑥̈ 𝑖 + 𝐶𝑖𝑗 𝑥̇ 𝑥 + 𝑘𝑖𝑗 𝑥) = 𝑓𝑖 (𝑖 = 1,2, … … … … , 𝑛) (4.57)

Comme il s’agit d’un système libre conservatif. 𝑓𝑖 = 0, λ𝑖𝑗 = 0
𝑛
(4.58)
∑(𝑚𝑖𝑗 𝑥̈ 𝑖 + 𝑘𝑖𝑗 𝑥) = 0 (𝑖 = 1,2, … … … … , 𝑛)
𝑗=1
On cherchant les solutions de type :
𝑥𝑖 = 𝑎𝑖 ejωt (4.59)
|−𝑚𝑖𝑗 𝜔2 + 𝐾𝑖𝑗 | = 0
87
4.7 Détermination des modes propres d’un élément poutre et matrice masse
Il s’agit que la détermination des modes propres d’une vibration d’un élément
plaque (rectangle) comme montré dans la (figure 4.4), basée essentiellement sur
la détermination de matrice masse de la structure sachant que la matrice de
rigidité a été déterminer dans le deuxième chapitre. L'élément a une longueur a,
la largeur b, et l’épaisseur constante t, chacun des quatre coins possède deux
degrés de liberté, les déplacements u et v dans la direction de x et de y,
respectivement les quatre coins s'appellent habituellement les points nodaux.
ainsi cet élément possède huit masses nodales (quatre paires de mx et my) et huit
déplacements nodaux ou degrés de liberté nodaux (quatre paires de u et de v).
y, v
m3,v3
m4,v4
m1,v1
m2,v2
Figure.3.4Elément rectangulaire de huit degrés de liberté.
4.8 Détermination de la matrice masse
La détermination des modes propres d’une structure par la technique des

éléments finis nécessite le remplace de la force d’inertie due à la masse de la
structure par des forces d’inerties équivalentes appliquées aux nœuds.
La forme fléchie d’une structure sujette à des vibrations en mode propre est
donnée par :
{𝑓(𝑥, 𝑦, 𝑧)}. 𝑓(𝑡) = [Φ(x, y, z){α}sinωn t (4.60)
88
Les forces d’inerties par unité de volume dues aux vibrations de la structure ont
les composantes suivantes :
𝜕 2𝑢
𝑋 = −𝜌 2
𝜕 𝑡
𝜕2 v
𝑌 = −𝜌 𝜕2 𝑡 (4.61)
𝜕 2𝑤
𝑍 = −𝜌
𝜕 2𝑡
Où 𝜌 la masse volumique.le vecteur de charge répartie {q(x, y, z)}, sera donc

comme suite :
𝜕2 𝑢
𝜌
𝜕2 𝑡
𝜕2 𝑣 (4.62)
{q(x, y, z)} = 𝜌
𝜕2 𝑡
𝜕2 𝑤
{𝜌 𝜕2 𝑡 }
Sachant que :
𝑢(𝑥, 𝑦, 𝑧)
{f(x, y, z)}𝑓(𝑡) = { 𝑣(𝑥, 𝑦, 𝑧) } 𝑓(𝑡) (4.63)
𝑤(𝑥, 𝑦, 𝑧)
En combinant les deux équations (4.60) et (4.62) :
𝜕2 (4.64)
{q(x, y, z)} = −𝜌 {[Φ(x, y, z){α}sin𝜔𝑛 𝑡]}
𝜕2𝑡
= −𝜔𝑛2 𝜌{[Φ(x, y, z){α}sin𝜔𝑛 𝑡]}
Tel que:{α} = [𝐴1−1 ]{𝛿}

Ainsi, que :
{q(x, y, z)} = 𝜔𝑛2 𝜌[Φ(x, y, z)[𝐴1−1 ]{𝛿}sin𝜔𝑛 𝑡] (4.65)
En substituant l’équation (4.65) dans (4.60), on obtient :
𝑄𝑐,𝑟 = {𝑑𝛿}𝑇 [𝐴1−1 ]𝑇 . ∫𝑣 [Φ(x, y, z)]𝑇 . 𝜔𝑛2 𝜌[Φ(x, y, z)][𝐴1−1 ]{𝛿}sin𝜔𝑛 𝑡𝑑𝑣 (4.66)
89
Le travail virtuel des forces d’inerties équivalentes appliquées aux nœuds est:
𝑄𝑐,𝑟 = {𝑑𝛿}𝑇 . {𝑁}
En égalisant les deux travaux obtenus on aura:
{𝑁} = 𝜔𝑛2 𝜌[𝐴1−1 ]𝑇 [∫ [Φ(x, y, z)]𝑇 . [Φ(x, y, z)]𝑑𝑣] . [𝐴1−1 ]{𝛿}sin𝜔𝑛 𝑡

𝑣
{𝑁} = 𝑤𝑛2 [𝑀𝑒 ]{𝛿}sin𝜔𝑛 𝑡
Donc:
[𝑀𝑒 ] = 𝜌[𝐴1−1 ]𝑇 [∫𝑣 [Φ(x, y, z)]𝑇 . [Φ(x, y, z)]𝑑𝑣] . [𝐴1−1 ] (4.67)
[𝑀𝑒 ] la matrice masse élémentaire pour un élément plaque.
1 𝑥𝑦𝑥𝑦 0 0 0 0
[Φ(x, y)] = [ ]
0 0 0 0 1 𝑥𝑦𝑥𝑦
La matrice masse d’un élément rectangulaire comme suite:
 1 1 1 1 
 9 0 0 0 0 
18 36 18
 1 1 1 1 
 0 0 0 0 
 18 9 18 36 
 1 1 1 1 
 36 0 0 0 0 
18 9 18 (4.68)
 1 1 1 1 
 0 0 0 0 
M e    .a.b.h  18 36 18 9
1 1 1

1 
 0 0 0 0
9 18 36 18 
 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 
 18 36 18 9 
90
Chapitre 5 :
Analyse dynamique des

structures mécaniques
discrètes utilisant la MEF
91
Chapitre 5 : Analyse dynamique des structures mécaniques discrètes par EF
5.1 Introduction
La méthode des éléments finis est basé sur la forme géométrique de la

structure (élément barre, élément poutre, élément plaque ou élément coque…).
La discrétisation de ces structures fournit une matrice élémentaire de rigidité en
fonction de l'énergie totale de la déformation, de contrainte, pour fournir une
matrice de rigidité globale. Après la mise en œuvre des conditions aux limites et
des charges, nous calculons l'inconnu dans tous les nœuds déplacements et les
contraintes (force axiale) dans tous les éléments. Les forces sont appliquées à la
structure afin de déterminer la matrice de masse. L'analyse dynamique de
flexion sélectionne les meilleures conditions pour limiter la structure.
Dans ce chapitre, un programme de calcul établi par langage PASCAL a

été développé. Les déplacements, les forces et les réactions aux nœuds ainsi que
les forces axiales dans chaque élément et les modes propres des toutes les
structures ont été analysés sous différentes charges appliquées et les conditions
aux limites. Les résultats obtenus sont comparés à ceux trouvés en utilisant des
programmes existants, tels que "SAP 2000" et "Abaqus". Une bonne
comparaison a été observée.
5.2 structures en éléments barre et poutre
Dans notre étude, nous avons choisi trois exemples d’une structure en élément
poutre avec des différentes conditions aux limites. D'une part, afin de voir l'effet
de la condition aux sur les résultats de réactions ou les forces axiales et même
sur les modes propres de la toute la structure.
D'autre part, en voir l’effet de la structure au même sur notre résultats, tel que la
forme et caractéristiques physiques et mécaniques du matériau,
92
(a)
(b)
(c)
Figure.5.1structures en élément poutre avec des différentes conditions aux

limites
5.3 Caractéristique de la structure en élément poutre
La structure est un élément de poutre métallique et supporté uniquement à

ces extrémités. Ces caractéristiques géométriques ont une section carrée S =
0,0244 m2, une longueur L = 2,85 m et une hauteur H = 0,45 m. Les propriétés
mécaniques du matériau sont donnés par le module d'Young, E = 200.10 9 (N /
m2) et la densité ρ = 7849 kg.m-3, La force appliqué F= 10000,00 N
5.4 Les déplacements nodaux suivant l’axe ox
On peut remarquer d’après la figure (5.2), que le déplacement nodal suivant

l’axe ox, sera maximal dans le nœud 2 et le nœud 17 dans le premier cas de
condition aux limites figure (5.1.a), tel que 𝑢2 = 0.000727 𝑚 et 𝑢17 =
0.00075 𝑚 , ce qui indique que la structure supporte une déplacement dans les
deux nœuds,
93
Displacement to the axis X (m )

0,0008 SAP 2000 program
Developed program
0,0007
0,0006
0,0005
0,0004
0,0003
0,0002
0,0001
0,0000
-0,0001
0 2 4 6 8 10 12 14 16 18 20
Node
Figure.5.2 Comparaison des déplacements suivant l’axe ox
D’après les résultats obtenus on peut remarquer que les déplacements

nodaux suivant l’axe oxseront minimales dans la structure 𝑁 0 1 , qui permet de
choisir la meilleure condition aux limites, on peut remarquer aussi dans ce cas la
symétrie de déplacement nodaux par rapport les nœuds symétriques, on peut
observer aussi que dans les deux nœuds 2 et 17 les déplacements seront
maximales.
94
5.5 Les déplacements nodaux suivant l’axe oy
Dispalcement to the axis Y (n) (m)

0,0002
Sap 2000 program
0,0000
developed program
-0,0002
-0,0004
-0,0006
-0,0008
-0,0010
-0,0012
-0,0014
-0,0016
0 2 4 6 8 10 12 14 16 18 20
Node
Figure.5.3Comparaison des déplacements suivant l’axeoy
D’après les résultats obtenus on peut remarquer que les déplacements

nodaux suivant l’axe oyseront minimales dans la structure 𝑁 0 1 , qui permet de
choisir la meilleure condition aux limites, on peut remarquer aussi dans ce cas la
symétrie de déplacement nodaux suivant l’axe oy par rapport les nœuds
symétriques, on peut observer aussi que le déplacement maximal dans le nœud
10 et nul dans les deux nœuds 2 et 18. Alors Les déplacements les plus grands
situés au milieu de la structure.
95
5.6. Les déplacements angulaires
Sap 2000 progam

Developed program
0,0002
Angle of rotation rad)

0,0001
0,0000
-0,0001
-0,0002
0 2 4 6 8 10 12 14 16 18 20
Nodes
Figure.5.4Comparaison des déplacements angulaire
Dans le cas des déplacements angulaire, on peut remarquer la diminution

des valeurs des déplacements selon arrangement des nœuds, et nul dans le nœud
1 et maximale dans le nœud 18 puisque librement appui, ce qui explique le
renforcement de la structure à ce niveau.
5.7. Les réactions
Sap 2000 program

15030 Devloped program
15020
Reactions (m)
15010
15000
14990
14980
14970
0 2 4 6 8 10 12 14 16 18
Node
Figure.5.5comparaison des déplacements angulaire
96
Dans le calcul des réactions aux appuis, on peut constater que nous avons
deux appuis sur le nœud 1 et 17, alors que la réaction suivant l’axe ox, est
minimale comparer avec la réaction suivant l’axe oy. Ce qui aura un impact très
important sur les modes propres de la structure.
5.8. Les forces axiales
30000
Sap 2000 program
Developed program
20000
The axial strengths (N)
10000
-10000
-20000
-30000
0 10 20 30 40
Element
Figure.5.6Comparaison des forces axiales
En ce qui concerne les forces axiales dans chaque élément de la structure,

on peut remarquer que, Il ya quelques éléments supporte des contraintes
maximale, comme dans les éléments 17, 20, 22, 25 et 27, Cela nécessite le
renforcement sur élément, sachant que nous avons choisi la meilleure condition
aux limite pour que aussi les forces axiales soit minimale. on peut voir aussi à
partir de la figure (Fig.5.6) qu’elles sont les forces axiales les plus grandes, donc
si la structure subit des chargements pareils le problème de sécurité ne se pose
pas, par contre si la structure soumise à des chargement plus importants dans ce
cas on peut jouer sur les conditions aux limites.
97
5.9. Les modes propres
Structure Structure Structure (b) Structure (c)

(a)
Fréquence (Hz) 23,83 170,47 23,83
De premier mode
Fréquence (Hz) 91,51 352,54 133,83
De deuxième
mode
Fréquence (Hz) 165,36 375,52 185,77
De troisième mode
Fréquence (Hz) 302,19 552,45 326,51
De quatrième
mode
Tableau.5.1Fréquences des modes propres des structures (Hz)
Mode propre 1 de la structure (a)
98
Mode propre 1 de la structure (b)
Mode propre 1 de la structure (c)
99
Figure.5.7 les Modes propres de trois structures
100
Chapitre 6 :
Analyse dynamique de
structures continues planes
utilisant la MEF
101
Chapitre 6 : Analyse dynamique de structures continues planes utilisant EF
6.1 Introduction
Apres une étude statique, en utilisant les résultats obtenus (chapitre 3), pour une
analyse dynamique de la structure continu (élément plaque, en acier la plaque
est mince, dont l’épaisseur est petit par rapport aux autres dimensions). Notre
analyse basée essentiellement sur la détermination des modes propres (les
Cinque premier modes), qu’ils dominent le comportement de la réponse
dynamique de toute la structure. Et aussi la détermination des déplacements
nodaux de chaque mode
6.2 Structures en éléments plaques

Nous considérons la plaque rectangulaire représentée sur la figure
6.1.encastrée à l’une de ses extrémités et libre sur le reste de son pourtour. Les
caractéristiques du matériau élastique sont les suivants, la longueur est 2m et de
largeur de 1m, le module de Young 𝐸 = 2.1011 𝑁⁄𝑚2 , et de masse
volumique ρ = 7800kg/m3, avec un module de poisson µ = 0.3, et d’épaisseur =
0,05 m, Pour la résolution par éléments finis, on utilise un maillage régulier
composé de 16 éléments et 25 nœuds, pour notre analyse et pour une validation
par logiciel Abaqus on utilise un maillage plus précis de 200 éléments et 231
nœuds, et pour confirmer la validité de nos résultats on a utilisé un maillage
tétraédrique pour logiciel Ansys 17.0 d’un côté pour voire l’effet du maillage
sur les résultats obtenus et d’autre côté voire aussi les différences entre les types
de maillage. Et les éléments avec cisaillement transverse. Ses éléments utilisent
une approximation quadratique du déplacement hors plan et des rotations.
Chaque élément dispose de huit nœuds et donc de 58 degrés de liberté, avec les
conditions aux limites on obtenir les résultats, des déplacements nodaux, des
réactions aux appuis, des forces axiales dans les éléments et les modes propres
de toute la structure.
102
Figure.6.1 élément plaque rectangulaire encastré-libre
Figure.6.2 élément plaque rectangulaire encastré-encastré
Afin de l’analyse dynamique de la structure en élément plaque, nous obtenons la

réponse dynamique telle que les fréquences et les modes propres de la structure.
Le premier mode est à une fréquence qu’il s’agit du premier mode de flexion
de la plaque.
103
Figure.6.3 Mode propre 1 d’une plaque encastrée-libre
Le deuxième mode est à une fréquence beaucoup plus élevée. Il s’agit du

premier mode de torsion. On l’appelle mode de torsion par comparaison avec le
Phénomène de torsion des poutres. Il s’agit bien évidemment d’un mouvement
de flexion de la plaque. [28]
Dans la suite, nous étudierons une vibration de flexion puis une vibration de
torsion, ces deux modes auront, on s’en doute, une importance capitale pour ces
deux types de vibration. Les autres modes de la plaque sont les suivants :
Le troisième mode c’est le deuxième mode de flexion.
104
Le quatrième mode c’est le deuxième mode de torsion.
Figure.6.6Mode propre 5 d’une plaque encastrée-libre
Le cinquième mode c’est le troisième mode de flexion.
Figure.6.7Mode propre 5 d’une plaque encastrée-libre
105
Dans notre étude nous sommes limités les Cinque premiers modes de la
structure puisqu’ils présentent la plupart de la réponse dynamique de la
structure.
6.3 Premier Cas : encastrement d'un côté et libre d'autre coté

6.3.1 les modes propres de la structure
(b)
(a)
Figure.6.8Première mode de la structure A, (a)Ansys, (b) Abaqus.
(b)
(a)
Figure. 6.9Deuxième mode de la structure A, (a)Ansys, (b) Abaqus.
106
(b)
(a)
Figure6.10 Troisième mode de la structure A, (a)Ansys, (b) Abaqus.
(b)
(a)
Figure. 6.11Quatrième mode de la structure A, (a)Ansys, (b) Abaqus.
(b)
(a)
Figure. 6.12Cinquième mode de la structure A, (a)Ansys, (b) Abaqus.
107
6.3.2Déplacements nodaux
(a) Déplacements nodaux de la structure pour Mode propre 1( premier cas)
(b)Déplacements nodaux de la structure pour Mode propre 2(premier cas)
108
( c)Déplacements nodaux de la structure pour Mode propre 3 ( premier cas)
(d)Déplacements nodaux de la structure pour Mode propre 5 (premier cas)
(e)Déplacements nodaux de la structure pour Mode propre 5 ( premier cas)
Figure.6.13 Les déplacements nodaux de Cinque premiers Modes(premier cas)
109
6.3.3 Réactions aux appuis
(a) Les réactions aux appuis de la structure pour Mode propre 1 de la structure premier cas
(b)Les réactions aux appuis de la structure pour Mode propre 2 de la structure premier cas
110
( c)Les réactions aux appuis de la structure pour Mode propre 3 de la structure premier cas
(d)Les réactions aux appuis de la structure pour Mode propre 5 de la structure premier cas
111
(e)Les réactions aux appuis de la structure pour Mode propre 5 de la structure premier cas
Figure. 6.14 Les réactions aux appuis de Cinque premiers Modes (premier cas)
6.4 Deuxième Cas : encastrement de deux cotés
6.4.1 les modes propres de la structure
(b)
(a)
Figure.6.15Première mode de la structure B, (a)Ansys, (b) Abaqus.
112
(b)
(a)
Figure.6.16Deuxième mode de la structure B, (a)Ansys, (b) Abaqus.
(b)
(a)
Figure. 6.17Troisième mode de la structure B, (a)Ansys, (b) Abaqus.
(b)
(a)
Figure. 6.18Quatrième mode de la structure B, (a)Ansys, (b) Abaqus.
113
(b)
(a)
Figure. 6.19Cinquième mode de la structure B, (a)Ansys, (b) Abaqus.
6.4.2 déplacements nodaux
(a) Déplacements nodaux de la structure pour Mode propre 1(deuxième cas)

deuxièmecas
114
( b)Déplacements nodaux de la structure pour Mode propre 3(deuxième cas)
( c)Déplacements nodaux de la structure pour Mode propre 2 (deuxième cas)
115
(d)Déplacements nodaux de la structure pour Mode propre 5 (deuxième cas)
(e)Déplacements nodaux de la structure pour Mode propre 5 (deuxième cas)
Figure.6.20 Les déplacements nodaux de Cinque premiers Modes (deuxième cas)
116
6.4.3 Réactions aux appuis
(a) Les réactions aux appuis de la structure pour Mode propre 1 (deuxième cas)
deuxièmecas
(b)Les réactions aux appuis de la structure pour Mode propre 2 (deuxième cas)
117
( c)Les réactions aux appuis de la structure pour Mode propre 3 (deuxième cas)
(d)Les réactions aux appuis de la structure pour Mode propre 5 (deuxième cas)
118
( e)Les réactions aux appuis de la structure pour Mode propre 5 (deuxième cas)
Figure.6.21 Les réactions aux appuis de Cinque premiers Modes (deuxième cas)
119
6.5 Les déplacements nodaux suivant l’axe ox
0,0010
Programme
0,0008
Déplacement suivant l'axe X (u)
Sap2000
0,0006
0,0004
0,0002
0,0000
-0,0002
-0,0004
-0,0006
-0,0008
-0,0010
0 5 10 15 20 25
noeuds
Figure. 6.22 Comparaison des déplacements nodaux suivant l’axe ox

On remarque que les déplacements nodaux suivant l’axe ox est maximale
dans les nœuds 25 est nuls dans les nœuds 1, 2, 3, 5, 5 et 6 dépend la condition
au limite de la structure, tel que une meilleure structure qui subit les
déplacements nodale, Et identifier la valeur maximale du déplacement. Et son
impact sur la structure, Ce que nous savons la meilleure structure avec les
conditions aux limites spécifiques.
6.6 Les déplacements nodaux suivant l’axe oy
0,010 Programme
Déplacement suivant l'axe Y (v)
Sap2000
0,005
0,000
-0,005
-0,010
-0,015
0 5 10 15 20 25
noeud
Figure.6.23 Comparaison des déplacements nodaux suivant l’axe oy

On remarque que les déplacements nodaux suivant l’axe oy est maximale
120
dans les nœuds 11 et 15 est nuls dans les nœuds 1, 2, 3, 5, 5 et 6 dépend la

condition au limite de la structure, tel que une meilleure structure qui subit les
déplacements nodale, Et identifier la valeur maximale du déplacement. Et son
impact sur la structure, Ce que nous savons la meilleure structure avec les
conditions aux limites spécifiques.
6.7 Les contraintes normales dans les éléments
résultat de programme dévloppé

Résultat de logicielle Sap2000
9.00E+008
8.50E+008
8.00E+008
7.50E+008
X
7.00E+008
6.50E+008
6.00E+008
0 2 4 6 8 10 12 14 16 18
Element
Figure.6.24 Comparaison des contraintes normales
On remarque que les contraintes normales pour cette structure sont

symétriques, Ce qui permet de dire que cette structure avec cette condition aux
limites, donnez-nous les meilleurs résultats par rapport à les contraintes
normales.
121
6.8 Les contraintes tangentielles dans les éléments
(contraine de programme dévloppé)

6,00E+008
logicielle Sap2000
4,00E+008
2,00E+008
0,00E+000
XY
-2,00E+008
-4,00E+008
-6,00E+008
0 2 4 6 8 10 12 14 16 18
Element
Figure.6.25 Comparaison des contraintes tangentielles

La même chose remarquée sur les contraintes tangentielles, soit sur la
symétrie de la distribution de ces contraintes, soit sur la valeur maximale Et son
impact sur la structure.
6.9 Les modes propres suivant l’axe ox
Mode N:1 Programme

Mode N:2 Sap2000
0,06
0,04
0,02
Les Modes
0,00
-0,02
-0,04
-0,06
-0,08
0 2 4 6 8 10
noeuds
Figure.6.26 Comparaison des modes propres suivant l’axe ox

En fait, cette étude est basée principalement sur la détermination des
122
modes propres de la structure. Et nous visons toujours le premier mode de la

structure.
Donc, chaque fois qu'ils les valeurs des modes propres sont minimales la
structure être en bonne condition qui peut résister à des forces vibratoires.
6.10 Les modes propres suivant l’axe oy
Mode N:1
0.06 Mode N:2
0.04
0.02
Les modes
0.00
-0.02
-0.04
0 5 10 15 20 25
noeuds
Figure.6.27 Comparaison des modes propres suivant l’axe oy

En fait, cette étude est basée principalement sur la détermination des
modes propres de la structure. Et nous visons toujours le premier mode de la
structure.
Donc, chaque fois qu'ils les valeurs des modes propres sont minimales la
structure être en bonne condition qui peut résister à des forces vibratoires
123
6.11Conclusion
Le but de cette analyse vibratoire d’une structure tridimensionnelle (élément plaque),
est la détermination des fréquences propres de la structure (vibration libre), dans les
différentes conditions aux limites.
On a déterminé les cinq premiers modes des deux structures avec des différentes
conditions aux limites telles que les quatre premiers modes dominent le comportement
vibratoire presque 90 pour cent de la structure. Et le cinquième mode juste pour voir la
réponse vibratoire après les quatre premiers modes.
Concernant les déplacements nodaux on peut dire que la structure subie des
déplacements maximaux au milieu ce qui explique le renforcement de la plus part des
structures, dans les deux structures,
Pour les forces axiales et les réactions aux appuis les valeurs maximales se concentrent
autour des éléments les plus prêt de l’encastrement par contre elles se situent au milieu
de structure dans le cas d’encastrement de deux cotés.
124
Conclusion générale
125
Conclusions et perspectives
Conclusion générale
Dans notre étude, nous nous sommes intéressés à l’étude du comportement

statique et dynamique des structures mécaniques par la méthode des éléments
finis, avec des différents types de conditions aux limites. Notre but étant la
détermination des déplacements pour chaque nœud et les forces axiales pour
chaque élément et les réactions aux appuis de structure. Pour se faire, nous
avons commencé notre travail par l’analyse statique des différentes structures en
élément barre, poutre et plaque avec des différentes conditions aux limites. Pour
voir l’effet de la condition aux limites sur le comportement statique et
dynamique de la structure, et l’introduction des propriétés mécaniques du
matériau de la structure permet de choisir la meilleure structure subie des forces
statique ou vibratoires. Une nécessité de développer un programme de calcul.
Notre programme en langage (Pascal), et les résultats obtenus comparées avec
les résultats de logiciel SAP2000, et une validation par Abaqus. Dans un
premier temps, nous avons étudié une nouvelle démarche qui consiste à
discrétiser la structure globale au plusieurs structures élémentaires, basée sur les
théories de la méthode des éléments finis, tel que la détermination des matrices
élémentaires pour chaque élément et avec l’assemblage on peut avoir la matrice
globales de toute la structure.
Cette matrice de rigidité nous a permis de calculer les déplacements dans
chaque élément et les forces axiales dans les éléments et les réactions sur les
appuis.
Dans un deuxième temps, nous nous sommes intéressés au développement

d’une nouvelle approche qui permet d’étudier le comportement dynamique de
ces structures, tel que la détermination de matrice masse très important pour
notre calcul de modes propres de la structure, Les comparaisons effectuées entre
125
nos résultats numériques avec ceux tirés de logiciel Sap 2000 ont montré une
bonne concordance, ce qui prouve la validité de la méthode d’analyse. Les
résultats obtenus montrent que nous pouvons choisir la meilleure structure avec
la meilleure condition aux appuis, On pourrait dire aussi que nous obtenons des
meilleurs résultats au niveau des modes propres, puisque notre étude statique
était bonne.Par ailleurs, pour une compréhension plus concrète de la méthode
des éléments finis pour l’analyse statique et dynamique des structures
mécaniques basée sur comprendre les principes de la mécanique de milieu
continu, comme nous avons conclu que la méthode est très puissante pour ce
type d’analyse.
Notre programme de calcul nous a permis de calculer numériquement le

comportement statique et dynamique des structures mécanique. Et nous avons
aussi étudié l’influence des paramètres géométriques de la structure, l’effet des
conditions aux limites sur les résultats.
126
Perspectives
127
Perspectives
En perspectives, nous proposons de développer notre programme, Il s’agit en
fait de tout un processus qui commence par la construction de la géométrie de la
structure à l’aide d’un logiciel de dessin tels que SOLIDWORKS ou bien
AUTOCAD.
Tel que l’utilisation de ces logiciels permettant de traiter des structures à

géométries complexes. Comme application, nous avons étudié différents types
des structures avec des géométries plus ou moins compliquées, tels que ces
structures en élément barre , poutre et plaque Se compose de formes simples et
parfois difficile. Comme on peut être dans l'avenir que nous utilisons cette
méthode pour l’analyse des structures dans des formes plus complexes. Surtout
quand nous utilisons des structures avec des matériaux intelligents, Lorsque les
formes de ces structures nécessitent des calculs très précis et les erreurs les plus
bas possibles. les turbines à pales droites, à pales incurvées et à pales inclinées.
Par ailleurs, nous avons pu appliquer cette même méthode sur des structures à
géométries plus ou moins compliquées, et dans de nombreux domaines de
l’industrie surtout la mécanique , le génie civil, l’aéronautique et même la
biomécanique.
128
References bibliographiques
128
Résume
REFERENCES BIBLIOGRAPHIQUES
[1] A. Fortin, Les éléments finis de la théorie à la pratique, Ecole

polytechnique de Montréal, 2011.
[2] A. El Ghazi, S. El Hajji, Guide d’utilisation du logiciel MATLAB,

département mathématiques et informatique, Université Mohammed V,
Maroc, 2004.
[3] Pierre Spiteri, Présentation générale de la méthode des éléments finis,

Ecole nationale supérieur de Toulouse, article [AF505].
[4] R.A. Sauer, J.C. Mergel, 2014, A geometrically exact finite beam element
formulation for thin film adhesion and debonding, Finite Elements in
Analysis and Design 86, 120-135.
[5] G. Alotta, G. Failla, M. Zingales, 2014, Finite element method for a

nonlocal Timoshenko beam model, Finite Elements in Analysis and
Design, 89, 77-92.
[6] A. Ghoneim, 2015, A mesh free interface-finite element method for

modeling isothermal solutal melting and solidification in binary systems,
Finite Elements in Analysis and Design, 95, 20-41.
[7] J.G. Kim, J.K. Lee, H.J. Yoon, 2015, Free vibration analysis for shells of
revolution based on p-version mixed finite element formulation, Finite
Elements in Analysis and Design, 95,12-19
[8] R. Brighenti, S. Bottoli, 2014, A novel finite element formulation for

beams with composite cross-section, International Journal of Mechanical
Sciences, 89, 112-122
129
Résume
[9] M. Donà, A. Palmeri, M. Lombardo, A. Cicirello, 2015, An efficient two-

node finite element formulation of multi-damaged beams including shear
deformation and rotatory inertia, Computers and Structures, 147, 96-106.
[10] J.G. Teng, D. Fernando, T. Yu, 2015, Finite element modelling of

debonding failures in steel beams flexurally strengthened with CFRP
laminates, Engineering Structures 86, 213-224
[11] D.Y. Gao, J. Machalová, H. Netuk, 2015, Mixed finite element

solutions to contact problems of nonlinear Gao beam on elastic
foundation, Nonlinear Analysis: Real World Applications, 22, 537-550.
[12] M. Lalanne, P. Berthier, J. Der Hagobian, Mécanique des vibrations

linéaires, MASSON Edition, 1995.
[13] T.Y. Yang, Finite element structural analysis, Prentice-Hall, 1986.
[14] B. Necib, Théorie générale de la mécanique des milieux continus,

Université de Constantine, 2010.
[15] B. Necib, Les méthodes des éléments finis (MEF), Université de

Constantine, 2011.
[16] P. Trompette, Mécanique des structures par la méthode des éléments finis
(Statique et dynamique), 1992.
[17] M. Abdi, K. Mohamed, Numerical Simulation and Active Vibration

Control of Piezoelectric Smart Structure. International Review of
Mechanical Engineering, 3, 175-181, 2009.
[18] O.C. Zienkiewiez, The finite element method in engineering analysis,

Paris, 1991.
[19] A. Lebeid, B.Necib, 2011, Analysis and numerical modelling of ceramic

piezoelectric beam, International Review of Mechanical Engineering, 3,
130
Résume
454-456.
[20] N. Brie, Mise en œuvre d'un calcul de modes propres d'une structure, EDF
R&D, 2014.
[21] Antoine Kgay, calcul des structures par éléments finis, 2012, Cnam-Paris.
[22] Yves Debard,, RDM-Eléments finis, 2009, institut universitaire de

technologie de Nans.
[23] ANDRI Fortin, les éléments finis de la théorie à la pratique, 2011,

Université Laval.
[24] P.Ciarlet, Exercices d’analyse numérique matricielle et d’optimisation,

1988, Masson
[25] Daniel Choi, Méthode des éléments finis, avril 2010, université de Caen,
France.
[26] Michel KERN, Introduction a la méthode des éléments finis, 2005, Ecole
nationale supérieure des Mins de France.
[27] ERIC Goncalves da silva, Méthode et analyse numérique, janvier 2011.
[28] Géradin M. et Rixen D., Théorie des Vibrations - Application à la

dynamique des structures, Masson, 1997.
[29] Jean Michel Génevaux, Dynamique des structures, V8, Janvier 2016, Ecole
nationale d’ingénieurs du Mans.
[30] Sylvain Drapier, Dynamique des Solides et desStructures, octobre 2016,

Ecole Nationale Supérieure des Mines de Saint-Etienne.
[31] Alexei PANTCHICHKINE, IntroductionAuxsystèmesdynamiquesetàla

modélisation, 2008 , Institut Fourier, Université Grenoble-1
[32] Sylvain Drapier, Dynamique des Solides et des Structures, 4iem edition
131
Résume
octobre 2009, Ecole Nationale Supérieure des Mines de Saint-Etienne
[33] Agassant, J. F., P. Avenas, J. Ph. Sergent et P. Carreau:

PolymerProcessing, Principles andModeling. Hanser Publishers, Munich,
1991.
[34] Arnold, D., F. Brezzi et M. Fortin: A Stable Finite Element for the Stokes
Equations. Calcolo,21 :337{344, 1984.
[35] Bartle, R.G.: The elements of integration. Wiley, New York, 1966.
[36] Bathe, K. J.: Finite Element Procedures. Prentice Hall, 1996.
[37] Bathe, K. J. and T. Sussman: A Finite Element Formulation for Nonlinear

Incompressible Elasticand inelastic Analysis. Computers and Structures,
26(1/2) :357{409, 1987.
[38] Bonet, J. and D. R. Wood: Nonlinear continuum mechanics for _nite

element analysis. Cambridge,1997.
[39] Brezzi, F. and M. Fortin: Mixed and Hybrid Finite Element Methods.
Springer-Verlag, 1991.
[40] Burden, R.L. et J.D. Faires: NumericalAnalysis. Brooks/Cole, Pacific

Grove, septième édition,2001, ISBN 0-534-38216-9.
[41] Ciarlet, P.G.: The Finite Element Method for Elliptic problems. North-
Holland, Amsterdam,1986.
[42] Courant, R. L.: Variational Methods for the Solution of Problems of

Equilibrium and Vibration.Bulletin of the American Mathematical
Society, 49:1{23, 1943.
[43] Céa, J.: Approximation variationnelle des problémes aux limites. Ann.
132
Résume
Inst. Fourier (Grenoble),14 :345{444, 1964.
[44] Delfour, M.C. and J. P. Zol_ezio: Shapes and Geometries : Analysis,

Différential Calculus, andOptimization. Advances in Design and Control.
SIAM, Philadelphia, 2001.
[45] Dhatt, G. et G. Touzot: Une présentation de la méthode des éléments finis.

Les presses del'Université Laval, Québec, 1981.
[46] Duvaut, G.: Mécanique des milieux continus. Collection mathématiques

appliquées pour lamaîtrise. Masson, 1990.
[47] Duvaut, G. et J.L. Lions: Les inéquations en mécanique et en physique.

Dunod, Paris, 1972.
[48] Engelman, M. et M. Jamnia: Transient Flow Past a Circular Cylinder : a

Benchmark Solution.Int. J. Numer. Meth. Fluids, 11 :985{1000,
1990.381-382 Bibliographie
[49] Fortin, A.: Analyse numérique pour ingénieurs. Presses Internationales

Polytechnique, troisième édition, 2008.
[50] Fortin, A., M. Jardak, R. Pierre et J. J. Gervais: Localization of Hopf

Bifurcations in FluidFlow Problems. Int. J. Numer. Meth. Fluids, 24
:1185{1210, 1997.
[51] Fortin, M. and A. Fortin: Experiments with Several Elements for Viscous
Incompressible Flows.Int. J. Num. Meth. Fluids, 5 :911{928, 1985.
[52] Fortin, M. et R. Glowinski: Résolution Numérique de Problèmes aux

Limites par des Méthodesde Lagrangien Augmenté. Dunod, 1983.
[53] Gasquet, C. et P. Witomski: Analyse de Fourier et applications. Masson,

Paris, 1990.
[54] George, A. et J. Liu: Computer Solution of Large Sparse Positive Definite
133
Résume
Systems. Prentice-Hall, 1981.
[55] Jean, R.: Mesure et intégration. Les presses de l'université du Québec,

Montréal, 1975.
[56] Lucquin, B.: Equations aux dérivées partielles et leurs approximations.

Mathématiques à l'université. Ellipses, Paris, 2004.
[57] The MathWorks Inc.: Matlab User's guide, 1996.
[58] Medvinsky, A. B., S. V. Petrovskii, I. A. Tikhonova, H. Malchow et B. L.

Li: SpatiotemporalComplexity of Plankton and Fish Dynamics. SIAM
review, 44(3), 2002.
[59] Ozell, B. et C. Pic: Manuel d'utilisation du logiciel configurable de

visualisation scientifiqueVU. CERCA, Montréal, 1998.
[60] Philippin, G.: Cours d'analyse à l'usage des ingénieurs, tome 2. Presses de
l'Université deMontréal, Montréal, 1994.
[61] Raviart, P.A. et J.M. Thomas: Introduction à l'analyse numérique des

équations aux dérivéespartielles. Masson, Paris, 1983.
[62] Reddy, J.N.: Applied Functional Analysis and Variational Methods in

Engineering. McGrawHill, New York, 1986.
[63] Reddy, J.N.: An Introduction to the Finite Element Method. McGraw-

Hill, New York,deuxiémeédition, 1993.
[64] Ritz, W: UebereineneueMethodezurLosunggewisser Variations probleme

der mathemati-schenPhysik. J. ReineAngew. Math., 135 :1{61, 1908.
[65] Roache, P.J.: Code veri_cation by the method of manufactured solutions.

ASME, 124 :4{10,mars 2002.
[66] Schwartz, L.: Théorie des distributions. Dunod, Paris, 1965.
134
Résume
[67] Strang, G. et G. Fix: An Analysis of the Finite Element Method. Prentice-

Hall, EnglewoodCliés, 1973.
[68] Swokowski: Analyse. De Boeck Université, Bruxelles, cinquième édition,

1993.Bibliographie 383
[69] Turner, M.J., R.W. Clough, H.C. Martin et L.C. Topp: Stifness and
Deflection Analysis ofComplex Structures. Journal of the Aeronautical
Sciences, 23 :805{882, 1956.
[70] Zienkiewicz, O. C. et R. L. Taylor: The finite Element Method : basic

formulation and linearproblems, tome 1. McGraw-Hill, 1989.
135

HAM7246

Transféré par

Droits d'auteur :

Formats disponibles

HAM7246

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

HAM7246

Transféré par

Droits d'auteur :

Formats disponibles

‫الجمـهوريــة الجزائريــة الديمقراطيــة الشعبيــة‬

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

Ministère de l’enseignement supérieur ‫وزارة التعليــم العالــي و البحــث العلمـي‬

ANALYSE DYNAMIQUE DES STRUCTURES MECANIQUES

Soutenue publiquementle : ……………… /…………………/………………….........

Devant le Jury composé de:

Année Universitaire 2017-2018

The mechanicaldiscrete or continuumstructures are actually of great importance in the application

Je tiens à remercier Monsieur NECIB Brahim Professeur à l’université de

Je remercie Monsieur BOUCHOUCHA Ali, Professeur à l’université des frères

Monsieur MERABET Abderrazak, Professeur à l’université frères Mentouri

J’exprime aussi mes vives reconnaissances à Monsieur ZIED DRISS, Professeur à

Enfin, mes sincères remerciements s’adressent à tous les enseignants du l’université

A mes très chers parents

A ma très chère femme Fatima;

A mes chers enfants; Sohaib, mouncif, tamim

A mes très chers frères et mes sœurs

A mes très amis B.Bakri, A.Righi, M .Meddah

Table des matières

Table des matières …………………………………………………………………………… 1

Chapitre 1 : Introduction générale et Etat de l’art

Chapitre 2 : Techniques de modélisation utilisées en éléments finis

Analyse dynamique des structures mécaniques par MEF

2.4.2 Connaissances requises pour la réalisation des logiciels MEF……………………..... 38

2.5.1 Types d’éléments…………………………………………………………………………. 39

2.5.3 Elimination directe de la bande de matrice……………………………………………….. 43

2.6 Types de numérotation…………………………………………………………………….. 44

2.6.1 Numérotation verticale…………………………………………………………………… 44

2.6.2 Numérotation horizontale………………………………………………………………… 44

2.7 Calcul d’erreurs dans l’analyse par MEF…………………………………………………… 45

2.7.1 Erreurs usuelles…………………………………………………………………………… 45

2.7.2 Erreurs détectées (Vérification du modèle) ……………………………………………….. 46

2.8 Critique des résultats obtenus par MEF……………………………………………………... 47

Chapitre 3: Détermination de matrice de rigidité des éléments des structures mécaniques

3.2 Détermination de la matrice de rigidité de l’élément de barre…………………………. 51

3.2.1Matrice de rigidité en coordonnées locales…………………………………………….. 51

3.2.2Matrice de rigidité en coordonnées globales…………………………………………… 54

3.3Détermination de la matrice de rigidité pour un élément poutre……………………….. 58

3.3.1 Matrice de rigidité en coordonnées locales……………………………………………. 58

3.3.2Equation de rigidité de l’élément poutre en coordonnées locales………………………. 60

Analyse dynamique des structures mécaniques par MEF

3.3.3Matrice de rigidité en coordonnées globales…………………………………………… 62

Chapitre 4 : Mise en équation de l’analyse dynamique des structures mécaniques

Chapitre 5: Analyse dynamique des structures mécaniques discrètes utilisant la MEF

Chapitre 6: Analyse dynamique de structures continues planes utilisant la MEF

Analyse dynamique des structures mécaniques par MEF

6.3.1 les modes propres de la structure………………………………………………………… 106

Analyse dynamique des structures mécaniques par MEF

Index des notations

 : Tenseur des contraintes N / m  2

ε: Les déformations internes de la structure

𝑋̅, 𝑌̅ : Coordonnées locales

u1 et u 2 : Déplacements nodaux m

Analyse dynamique des structures mécaniques par MEF

U : Énergie de déformation totale N.m

F (m) : Force de volume

K ij : Coefficients de la matrice de rigidité

mij : Coefficients de la matrice de masse

W : Le travail des forces volumiques et surfaciques

 T : Transposé d’un vecteur