Cours MEF

Introduction à la
Méthode des Éléments Finis
20h
Cursus ingénieur par alternance
2ième Année itii

Objectifs :
Aborder les notions fondamentales de la méthode des éléments finis à partir de l'étude des
structures treillis et portiques.
Présenter sa généralisation aux problèmes généraux de la physique.

Application :
A partir de ces notions l'étudiant doit pouvoir utiliser les outils numériques tels que MEFLAB
et/ou RDM6.
L’application à des problèmes de l’ingénieur plus complexes, en utilisant un code de type
industriel seront abordées en Troisième année.

Déroulement : sur 5 séances de 4h

Séance 1 : Étude des treillis - L’élément barre 2h cours + TD1
Séance 2 : Erreur d’approximation « colonne » 2h cours + TD2
Séance 3 : L’élément poutre 2h cours + TD3
Séance 4 : Étude des portiques 2h cours + TD4
Séance 5 : Généralisation 2h cours + 2h DS

Supports :
Tous les supports utiles sont sur le site : [Link]
Vous avez le polycopié de cours à conserver pour l'an prochain
Ci-dessous les textes de TD que nous traiterons cette année.
 Étude des treillis
TD 1 : Exercice Tr-1
Thème : Mise en œuvre des techniques éléments finis, changement de base, assemblage, résolution, calcul des efforts aux
appuis, calcul des contraintes, et calcul des efforts aux nœuds.

Intéressons-nous à la réponse statique du treillis plan représenté par la figure ci-

contre. Toutes les barres ont les mêmes caractéristiques (E, S). A
Modélisation
Définissez vos vecteurs globaux. A
Calcul de la matrice raideur
Exprimer les matrices raideurs de chaque élément
En déduire l'expression de la matrice raideur assemblée complète. G
Extraire la matrice raideur réduite. mg
Résolution statique - Post traitement
Déterminer la déformée statique.
Calculer les efforts aux appuis, et vérifier l'équilibre global de la structure.
Calculer les efforts sur chaque élément, et vérifier l'équilibre du nœud chargé.
Écrire les équations d'équilibre de l’élément à 45°, retrouver les résultats précédents

TD 2 : Exercice Tr-5
Thème : Modes rigides et utilisation des symétries du problème.
Intéressons-nous à la réponse statique d’un carré de coté A , chargé par deux forces G
égales et opposées appliqués à deux sommets opposés. A F
Modélisation. A
Que pensez-vous de la matrice raideur de cette structure ?
Proposer des conditions aux limites pour résoudre ce problème.
G
Proposer un modèle qui tienne compte des symétries. F
Ce modèle résout-il les problèmes de modes rigides ?
Calcul de la matrice raideur
Exprimer les matrices raideurs de chaque élément
En déduire l'expression de la matrice raideur assemblée.
Résolution statique - Post traitement
Déterminer la déformée statique.
Calculer les efforts aux appuis, et vérifier l'équilibre global de la structure.
Calculer les efforts sur chaque élément, et vérifier l'équilibre des nœuds de la structure.
 Étude des portiques
TD 3 : Étude d’une poutre sous son poids propre
Objectifs : mise en œuvre de la méthode des éléments finis, et illustrer la notion d'erreur liée à l'approximation.
Nous cherchons la réponse statique de la poutre sur appuis représentée par la G
figure ci contre. g
Modèle à 1 élément.
Définissez vos vecteurs globaux : {U } {FI } (bilan inconnues – équations)
Déterminer la matrice raideur.
le vecteur force généralisé associé au poids propre.
En déduire le système réduit et calculer les déplacements nodaux, représenter la déformée.
Calculer la flèche au centre de la poutre, et comparer à la solution analytique
5 ρ gS A 4
v ( A / 2) = −
384 EI
Calculer les efforts aux appuis, et vérifier l'équilibre global de la structure.
Calculer les efforts sur l'élément et tracer les diagrammes de l'effort tranchant et du moment de flexion.
Comparer à la solution analytique.
Modèle à 2 éléments.
Déterminer la matrice raideur assemblée complète.
Déterminer le vecteur force généralisé associé au poids propre de la structure.
Prise en compte de la symétrie
Préciser le nouveau maillage et les conditions de symétrie.
Écrivez le système réduit des équations, calculez les déplacements nodaux, comparez à la solution analytique.
Calculer les efforts aux nœuds, comparez à la solution analytique.
Calculer les efforts sur l'élément et tracer les diagrammes de l'effort tranchant et du moment de flexion.
Comparer à la solution analytique.

TD 4 : Extrait du DS 2010.
Intéressons-nous à la réponse statique de la structure représentée sur la figure 1 ci-dessous.

G
yo 1 A 2
A 3

G A /2
xo
4
Figure1 : définition du système mécanique étudié
Cette structure est constituée d'une poutre encastrée à une extrémité, soutenue par un appui simple en son centre et à
l'autre extrémité par une barre. Le chargement est représenté par une charge répartie uniformément sur la moitié droite
de la poutre. Il en résulte que la poutre ne travaille qu'en flexion, et le barre uniquement à l'effort normal.
Les caractéristiques mécaniques sont : ES raideur de la barre
EI raideur de la poutre
Questions de cours
• Un modèle éléments finis nous donnera la solution exacte en tout point ? (oui – non – pourquoi ?)
• Les valeurs des efforts aux appuis seront exactes ? (oui – non)
• Les déplacements nodaux seront exacts ? (oui – non)

Modèle éléments finis

On utilise le modèle éléments finis le plus simple possible avec 2 éléments poutre et 1 élément barre. Vous utiliserez la
numérotation nodale indiquée sur la figure.

{U } = {v1 θ1 v2 θ 2 v3 θ3 v4 }
T
On posera

• Préciser les conditions aux limites du problème, et définissez le vecteur « des inconnues efforts » associé à ces
conditions aux limites.
• Exprimer la matrice raideur de chaque élément sur ses variables nodales, en déduire la matrice raideur
assemblée.
S A2
ES / A
On posera α = =
EI / A3 I
• Exprimer le vecteur associé aux charges extérieures.
• Pour α = 3 , résoudre les équations du système réduit et représenter la déformée de la structure.
• Calculer les efforts aux appuis, représenter sur une figure ces efforts, les équations d’équilibre sont-elles
vérifiées ?

• Donner l’expression du moment flexion sur chaque élément M f 1 ( x ) et M f 2 ( x ) , avec x la coordonnée

locale sur l’élément considéré
• Tracer le diagramme du moment de flexion sur la structure, qu’en concluez-vous ?
• Quelle est l’erreur maximale commise sur le calcul du moment de flexion, où se situe-t-elle ?
• Notre modèle est-il cependant suffisant pour dimensionner la structure ?
(Justifier votre réponse)
La structure étudiée est en fait une passerelle surplombant une voie rapide, dessinez le problème réel (section de la
double voie) en indiquant ses caractéristiques, et son chargement.