Cours ECN : Mthode des Elments Finis
Mthode des Elments Finis
MEF UE-35 : MEEFI
Projet pdagogique - site WEB Mthodes dapproximation Formulations variationnelles MEF: les lments finis Exemples dapplication
H. OUDIN
MMGC SIM [Link]@[Link]
Cours ECN : Mthode des Elments Finis
Objectifs de cet enseignement
Prsenter les principes de base de la MEF
Problmes lmentaires Formulations variationnelles Mthodes numriques comprendre gnraliser appliquer
Parcours pdagogique (site Web, poly)
Treillis Portiques Mthodes variationnelles EDP Mthodes numriques MEFLAB
Cours ECN : Mthode des Elments Finis
Objectifs de cet enseignement
Prsenter la notion de Modle
Hypothses de modlisation Comment formuler un problme de physique pour pouvoir le traiter numriquement Hypothses de discrtisation Comment le traiter numriquement.
Utiliser un code de calcul industriel
Aborder les problmes danalyse et de validation de modles via des exemples simples.
Cours ECN : Mthode des Elments Finis
Projet pdagogique
Ltudiant est lacteur principal de sa formation
Le parcours pdagogique est organis en thmes
Travail en Autonomie
Supports pdagogiques sur le WEB + le poly
Vous pouvez travailler chez vous
Activits proposes
Comprendre Apprendre Appliquer Valider
Vido & Prsentations PowerPoint (site) Polycopi + exercices corrigs (site) + QCM (site) Exercices (site) MEFLAB (site) Maple Exercices traits en TD
Cours ECN : Mthode des Elments Finis
Projet pdagogique Site WEB
Le menu donne accs aux documents en ligne
Objectifs : [Link] Travail en autonomie rgulier Pouvoir ragir en TD sur vos difficults
Cours ECN : Mthode des Elments Finis
Projet pdagogique
Les TD
Pour finir
Valider votre comprhension des principaux points de cours, Rpondre vos questions sur le thme tudi.
Travail personnel avant les TD
Confrence (2*1h)
Grgory LEGRAIN Nicolas CHEVAUGEON
lerreur de discrtisation XFEM
Evaluation
Note individuelle (coef 7) Note collective (coef 4) DS sans documents Projet pondr par votre TA
Cours ECN : Mthode des Elments Finis
Mthode des Elments Finis
Projet pdagogique - site WEB
Mthodes dapproximation
Formulations variationnelles MEF: les lments finis Exemples dapplication
Cours ECN : Mthode des Elments Finis
Mthodes dapproximation : gnralits
Systme physique continu
Mise en quations formulation mathmatique du problme (EDP) Formes diffrentielles Problme aux limites Rsidus pondrs Formulation mathmatique du problme (PTV) Forme Variationnelle Discrtisation du milieu Mthodes des lments finis
Systme physique discret Formulation mathmatique du problme
(q. de Lagrange)
Formes intgrales
Discrtisation
Mthodes Numriques
Forme matricielle
Cours ECN : Mthode des Elments Finis
Rsidus Pondrs : Formulation
Modle math. pos sur un domaine continu Systme d'quations diffrentielles : "EDP"
M D
Conditions aux limites M D C(u ) = e( M , t )
L (u ) = f ( M , t )
Si u solution approche R(u) : rsidu (erreur commise)
(M , t)
Rsoudre R (u ) = L (u ) f
= 0 sur D
fonction de pondration
R (u ) dV
=0
Annulation du Rsidu pondre sur le domaine
1re forme intgrale Ne tient pas compte des conditions aux limites du problme
Cours ECN : Mthode des Elments Finis
10
Rsidus Pondrs : Approximation
Soit une approximation n paramtres:
u = wi ( M ) qi
i =1
R(
n i =1
w q ) dV = 0 i i
Fcts de forme Une quation n inconnues
Comment construire un systme matriciel ?
i de 1 n
Nombre fini de Fcts de pondration
i (M)
R(
n w (M) q j j j =1
) dV = 0
Systme matriciel
Attention Fcts de forme doivent vrifier toutes les CL
en pratique cest impossible pour un Pb de lingnieur
Cours ECN : Mthode des Elments Finis
11
Rsidus Pondrs : Exemple
coulement d'un fluide visqueux incompressible dans une conduite
(u ) =
p 2 p1
Coefficient de viscosit cinmatique du fluide. le champ des vitesses
Champ inconnu
u = u ( x, y ) z
Soit lapproximation 1 paramtre u ( x , y ) = ( x 2 a 2 )( y 2 a 2 ) q elle vrifie les conditions aux limites M
p 2 p1 collocation a2 p p q = 0, 31 2 2 1 Galerkin a p p1 0, 2947 a 2 2 Solution de rfrence au centre q = 0, 25
u=0
Cours ECN : Mthode des Elments Finis
12
Mthode des Elments Finis
Projet pdagogique - site WEB Mthodes dapproximation
Formulations variationnelles
MEF: les lments finis Exemples dapplication
Cours ECN : Mthode des Elments Finis
13
Formulation variationnelle
EDP R(u) = 0 dans D CL sur D Forme intgrale 1 D R(u) = 0 Formulation CL sur D forte
Objectif : transformer la Forme intgrale 1 Pour faire apparatre les CL intgration par parties Forme intgrale 1 D R(u) = 0 CL sur D
TH d'Ostrogradsky
D g ( , u ) + 2 h ( , u ) = 0 Formulation CL1 sur 1 faible PTV C L 2 su r 2
CL sur les flux : drives spatiales de u
Cours ECN : Mthode des Elments Finis
14
Formulation variationnelle : Exemple
EDP
Champ inconnu : T temprature
Conduction thermique dans un Four
2 Four Rsistance
divq + r = 0 [Link] = d T =T d
avec
dans
Condition de flux Condition sur T
o Flux nul
sur 2 sur 1
o Pice chauffer Flux nul
q =
g r a d T flux de chaleur
1
Annulation de lerreur pondre T
(divq + r ) T
d TdS +
1
dV = 0
gradT . grad T dV +
r TdV +
i TdS = 0
flux inconnu
Cest la forme variationnelle du problme
Cours ECN : Mthode des Elments Finis
15
Formulation variationnelle : Exemple
Choix T = 0 sur 1 Champ virtuel thermiquement admissible
TTh admissible
gradT . grad T dV +
r TdV +
d TdS = 0
Nous obtenons une quation 1 champ T
Il faut satisfaire la condition :
T = Td
sur 1
La Formulation variationnelle est directement utilisable dans la Mthode des lments Finis . Rsultat MEFLAB doptimisation des rsistances pour que la temprature dans la pice soit proche de la consigne fixe (cest un projet EF). Ce que nous venons de prsenter pour un Pb de conduction thermique Peut tre fait pour dautres Pb de physique (cours en ligne, poly, exo de cours)
Cours ECN : Mthode des Elments Finis
16
Mthode des Elments Finis
Projet pdagogique - site WEB Mthodes dapproximation Mthodes variationnelles
MEF: les lments finis
Exemples dapplication
Cours ECN : Mthode des Elments Finis
17
Mthode des Elments Finis : MEF
Ides de base
Point de dpart : Formulation Variationnelle Approximation de la solution par sous-domaines : lments finis forme simple approximation sur des variables physiques
Forces nodales
Charge rpartie Dplacements imposs
Domaine continu
Domaine discrtis
Cours ECN : Mthode des Elments Finis
18
MEF : Approximation lments finis
Formulation Variationnelle PTV en Mcanique
uCA
u = u d sur D1 : u = 0
u . u dV
+ : dV
D
f . u dV
D2
T . u dS = 0
Efforts donns sur D2
Approximation lments Finis
Pour chaque lment :
D = De W = W e
{u ( M )} = [ N ( M ) ]{U e } { u } = [ N ( M ) ]{ U e }
(Galerkin)
Mmes familles de fonctions pour u et u
Cours ECN : Mthode des Elments Finis
19
MEF : Ecriture matricielle
Pour les efforts internes
: dV = ?
Rappel : Notation matricielle T { } =< xx , yy , zz , 2 xy , 2 xz , 2 yz >
{ } =< xx , yy , zz , xy , xz , yz > { ( M ) } = [ L ] {u ( M )} Oprateur gradient en petites dformations { ( M ) } = [ D ( M ) ]{ ( M ) } Loi de comportement
T = { } { }
Approximation EF { ( M ) } = [ L ] [ N ( M ) ] {U e } = [ B ( M ) ] {U e } { ( M ) } = [ D ( M ) ][ B ( M ) ]{U e }
De
: dV = { U e } [ K e ]{U e }
T
avec
[ K e ] = [ B ( M )]T [ D ( M )] [ B ( M )]
De
dVe
Matrice raideur lmentaire
Cours ECN : Mthode des Elments Finis
20
MEF : Ecriture matricielle
Pour les efforts externes
De
f . u dV e =
De
{ u }T { f } dV e
Approximation EF
{ u } = [ N ( M ) ] {U e }
T De
= { U e }
{ N ( M )} { f } dV e
T
Assemblage
Vecteur force gnralise lmentaire
D = De W = We
[K ] = [K e ]
e
{ F } = { Fe }
e
On dfini un vecteur global
{U }
[ K ]{
U } = {F }
Systme global
Pour la statique
Dmarche utilise pour ltude des treillis et des portiques
Cours ECN : Mthode des Elments Finis
21
MEF : Techniques numriques
Approximation nodale Exemple 1D
T1 0
T(s)
T2 1
Pb de temprature
2 nuds
approximation 2 paramtres : T = a0+ a1 s
Identification aux nuds :
Fonctions dinterpolation
T ( 0 ) = T1 T (1) = T2
T1 T ( s ) = [1 s ; s ] T2
Variables nodales signification physique
Exemple : approximation utilisant 3 lments
Cours ECN : Mthode des Elments Finis
22
Techniques numriques
Linaire
lments une dimension
Base polynomiale
(1 x )
Quadratique
(1 x x2 )
Cubique
(1 x x2 x3 )
Type Lagrange
Interpolation
1 1
N 1N2 1
1
N11 N
N2
N2
N3 N
3 N
N1
N4
2
1 s
s s 0 0
Type Hermite 2 variables par nud exemple : lment poutre v et
11
Cours ECN : Mthode des Elments Finis
23
Techniques numriques
lments triangulaires
lments deux dimensions
Les bases polynomiales sont compltes lments quadrilatraux
Les bases polynomiales sont incompltes lments toriques
zo
symtrie cylindrique
Cours ECN : Mthode des Elments Finis
24
Techniques numriques
lments ttradriques
lments trois dimensions
Les bases polynomiales sont compltes lments prismatiques bases incompltes lments hexadriques bases incompltes
Cours ECN : Mthode des Elments Finis
25
Techniques numriques
Drf
s,t,u
Transformation gomtrique
Drel
x,y,z
x =< N g ( s , t , u ) > {x n } y =< N g ( s , t , u ) > {y n } z =< N ( s , t , u ) > {z } g n
nuds
{x n }, {y n }, {z n }
Drivation : on montre
x 1 y = [J ] z s t u
J matrice jacobienne de la transformation Fct de Ng et
{xn },{yn },{zn }
==> matrices [B(s,t,u)]e
Intgration : on montre
f
De
(x, y, z)
dxdydz
Dref
f (s
t, u)
det [J ] dsdtdu
Cours ECN : Mthode des Elments Finis
26
Techniques numriques
Intgration numrique
Dref
dv
f ( )
i i =1
Npi
Calcul des matrices lmentaires [M e ] = < N ( ) > T < N ( ) > det [J ] dv ref
[K e ] = [ B ( )] T [ D ] [ B ( )]
Dref
Dref
det [J ] dv ref
Ng et {x n }, {y n }, {z n }
Pour chaque lment
Pour chaque point d intgration Calcul de [J] et [J]-1 au point d intgration Construction de [D] et [B] Calcul de [B]T [D] [B] det[J] i Calcul de [N]T [N] det[J] i Accumuler dans [K] et [M]
Cours ECN : Mthode des Elments Finis
27
Bilan : Dmarche lments finis
Domaine continu
Discrtisation gomtrique Construction de lapproximation nodale
{u } = [N e ]{u e }
Calcul des matrices lmentaires Assemblage
2 E de =
De
: dV = {u n }T
De
D B dV {u n }
[K ] = [K e ]
{u n }T [K e ]{u n }
Prise en Compte des Conditions aux limites et Rsolution de lquation matricielle valuation des grandeurs lmentaires
[K ]{U } = {FD }+ {FI }
Rsolution U { } dplacemen ts nodaux {FI } efforts de liaisons
(M )
} = [D ( M )]{ ( M ) } = [D ( M )][B ( M )]{u n }
Cours ECN : Mthode des Elments Finis
28
Mthode des Elments Finis
Projet pdagogique - site WEB Mthodes dapproximation Mthodes variationnelles MEF: les lments finis
Exemples dapplication
Cours ECN : Mthode des Elments Finis
29
Refroidissement d'une jante en alliage d'aluminium
Le remplissage, le refroidissement, le transfert de chaleur de la pice au moule, et la solidification sont modliss.
Laboratoire de mtallurgie physique de l'EPFL Ecole Polytechnique Fdrale de Lausanne
Cours ECN : Mthode des Elments Finis
30
SNCF
Mannequin HYBRID III 50% dformable
Modle lments finis du sige impact
Didier LEVEQUE
Cours ECN : Mthode des Elments Finis
31
En biomcanique
Interface os - prothse Articulation du genou ligaments, tendons, cartilages, mnisques
Jean ROYER - MMGC
Cours ECN : Mthode des Elments Finis
32
Orange II
Mat 45 m Voiles 1000 m2 110 Pieds ( 37,80m ) 30 Tonnes Carbone-Nomex
Design : Gilles Ollier
Gilles MARCKMANN Laurent GORNET - MMGC
Cours ECN : Mthode des Elments Finis
33
ECN - SNECMA
Modlisation de la perte dune aube dans un racteur
Laurent STAINIER - MMGC
Cours ECN : Mthode des Elments Finis
34
lments finis et Level Set method
Voir le site de FEDKIW : Stanford
milie MARCHANDISE JF. REMACLE : UCL Nicolas CHEVAUGEON : MMGC
Cours ECN : Mthode des Elments Finis
35
A vous de jouer
Soyez acteur de votre formation pour en profiter pleinement
