Chapitre 1
Chapitre 1
Chapitre 1
Chapitre 1 :
Chargée de la matière : Dr. Z.BENADLA
contact : [email protected]
I. Introduction
• Les logiciels de simulation par éléments finis sont de plus en plus
utilisés dans l'enseignement des Sciences de l'Ingénieur.
Méthode des éléments finis, M1 VOA– Semestre 2
2019‐2020 Dr. Z. BENADLA 7
II. Méthode des éléments finis
II.1 Définition
Méthode des éléments finis, M1 VOA– Semestre 2
C’est un outil de simulation numérique le plus utilisé et le plus efficace, utilisé par les
Ingénieurs de tout secteur.
Elle s’applique aux problèmes de mécanique des milieux continus (MMC).
La vaste utilisation de la MEF dans tous les domaines est due à la similarité des
différents problèmes de l’Ingénieur (équations différentielles, système d’équations
d’équilibre…etc.
2019‐2020 Dr. Z. BENADLA 8
II. Méthode des éléments finis…
II.2 Principe
Méthode des éléments finis, M1 VOA– Semestre 2
Modèle physique
Modèle EF Elément fini
2019‐2020 Dr. Z. BENADLA 9
II. Méthode des éléments finis…
II.2 Principe
Méthode des éléments finis, M1 VOA– Semestre 2
2019‐2020 Dr. Z. BENADLA 10
II. Méthode des éléments finis…
II.2 Principe…
Méthode des éléments finis, M1 VOA– Semestre 2
Théorie des poutres (RDM) Théorie des plaques
Théorie MMC
2019‐2020 Dr. Z. BENADLA 11
II. Méthode des éléments finis…
II.3 La discrétisation
Méthode des éléments finis, M1 VOA– Semestre 2
Exprimer résultats Données du départ
la théorie retenue est issue de la mécanique des milieux continus par une
transformation nommée discrétisation.
Toutes les théories de la MMC peuvent être discrétisées, 3D, 2D (contraintes planes,
déformations planes, axisymétrie...) ou théories des poutres, plaques ou coques.
La discrétisation de la MMC repose sur deux modifications :
• Une hypothèse cinématique supplémentaire sur l'allure du champ de déplacement,
qui est recherché sous la forme d'une combinaison linéaire (inconnue) de fonctions
de forme connues, définies plus tard ;
• Une transformation des équations d'équilibre du milieu continu, destinée à les
rendre mathématiquement compatibles avec l'hypothèse supplémentaire introduite.
2019‐2020 Dr. Z. BENADLA 12
II. Méthode des éléments finis…
II.3 La discrétisation…
Méthode des éléments finis, M1 VOA– Semestre 2
les lois utilisées au cours d'une simulation par éléments finis ne sont pas exactement les
lois de la MMC : il s'agit de lois discrétisées, donc différentes, qui conduisent ainsi à des
solutions différentes.
résolution approchée des
les éléments finis
équations aux dérivées partielles.
2019‐2020 Dr. Z. BENADLA 13
II. Méthode des éléments finis…
II.4 Le Maillage
Méthode des éléments finis, M1 VOA– Semestre 2
Dans une simulation par éléments finis, le produit est modélisé par un maillage,
dont la définition s'appuie sur deux objets géométriques.
1. Un découpage du domaine géométrique occupé par le produit en sous‐
domaines de formes simples (3D, 2D ou 1D (lignes)), selon la théorie dans
laquelle on se place
2019‐2020 Dr. Z. BENADLA 14
II. Méthode des éléments finis…
II.4 Le Maillage…
Méthode des éléments finis, M1 VOA– Semestre 2
Quelques types d'éléments bidimensionnels courants
2019‐2020 Dr. Z. BENADLA 16
II. Méthode des éléments finis…
II.5 Définition des nœuds et des éléments
Méthode des éléments finis, M1 VOA– Semestre 2
Un nœud est un point numéroté, défini par ses coordonnées dans un repère choisi
(cartésien, cylindrique ou sphérique).
Un élément est défini par ses connectivités.
Changement physique Changement mécanique
2019‐2020 Dr. Z. BENADLA 17
II. Méthode des éléments finis…
II.5 Définition des nœuds et des éléments…
Méthode des éléments finis, M1 VOA– Semestre 2
Sa forme
Son nombre
2019‐2020 Dr. Z. BENADLA 18
II. Méthode des éléments finis…
II.5 Définition des nœuds et des éléments…
Méthode des éléments finis, M1 VOA– Semestre 2
(c)
2019‐2020 Dr. Z. BENADLA 19
III. Fonctions de forme
III.1 Définition
Méthode des éléments finis, M1 VOA– Semestre 2
2019‐2020 Dr. Z. BENADLA 20
III. Fonctions de forme…
III.2 Barre soumise à la traction i j
Méthode des éléments finis, M1 VOA– Semestre 2
x = 0 →u(x)= u1
x = L →u(x)= u2
Et
Fonctions de forme
2019‐2020 Dr. Z. BENADLA 21
III. Fonctions de forme…
III.2 Barre soumise à la flexion i j
Méthode des éléments finis, M1 VOA– Semestre 2
La variation n’est pas linéaire entre les deux v(x) = α1 +α2 x +α3 x2 +α4 x3
nœuds , elle est de la forme : θ(x)= du(x)/dx = α2 +2α3 x +3α4 x2
x = 0 →v(x)= v1 et θ(x)= θ1
x = L →v(x)= v2 et θ(x)= θ2
2019‐2020 Dr. Z. BENADLA 22
IV.Avantages et inconvénients de la MEF
Avantages :
Méthode des éléments finis, M1 VOA– Semestre 2
•Simples,
•Adaptées au calcul
Inconvénients :
•Approchées
C’est la méthode numérique de référence mais
pas l’unique!
2019‐2020 Dr. Z. BENADLA 23
V. Analyse par éléments finis
A quoi sert l’analyse par EF?
Méthode des éléments finis, M1 VOA– Semestre 2
Avantages
Gains de coûts
Gains de temps… réduire le temps préalable à la mise sur le marché !
Création de modèles plus sûrs et de meilleure qualité
2019‐2020 Dr. Z. BENADLA 24
VI. Algorithme de la MEF
Structure à n nœuds et m éléments
Pour chaque élément Pour la structure
Méthode des éléments finis, M1 VOA– Semestre 2
Construction de la matrice de Assemblage des rigidités [Ke]
rigidité [ke] (repère local)
Assemblage de {F} en prenant en
Construction du vecteur de charges compte les éventuelles charges
{fe} (repère local) nodales
Calculer la matrice de passage [] Résolution du système [K]{U}={F}
liant les 2 repères local et global après introduction des conditions
aux limites (d’appuis)
Calculer la matrice de rigidité
exprimée en repère global [Ke]
Calcul du vecteur de charges {Fe}
(repère global)
Calcul des efforts interne en repère
local 25
2019‐2020 Dr. Z. BENADLA