Cours1 Amphi
Cours1 Amphi
Cours1 Amphi
H. OUDIN
MMGC SIM herve.oudin@ec-nantes.fr
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
Projet pdagogique
Les TD
Pour finir
Valider votre comprhension des principaux points de cours, Rpondre vos questions sur le thme tudi.
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
Mthodes dapproximation
Formulations variationnelles MEF: les lments finis Exemples dapplication
Formes intgrales
Discrtisation
Mthodes Numriques
Forme matricielle
L (u ) = f ( M , t )
Rsoudre R (u ) = L (u ) f
= 0 sur D
fonction de pondration
R (u ) dV
=0
1re forme intgrale Ne tient pas compte des conditions aux limites du problme
10
u = wi ( M ) qi
i =1
R(
n i =1
w q ) dV = 0 i i
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
11
(u ) =
p 2 p1
Champ inconnu
u = u ( x, y ) z
u=0
12
Formulations variationnelles
MEF: les lments finis Exemples dapplication
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
14
divq + r = 0 q.next = d T =T d
avec
dans
sur 2 sur 1
q =
g r a d T flux de chaleur
1
(divq + r ) T
d TdS +
1
dV = 0
gradT . grad T dV +
r TdV +
i TdS = 0
flux inconnu
15
gradT . grad T dV +
r TdV +
d TdS = 0
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)
16
17
Domaine continu
Domaine discrtis
18
u = u d sur D1 : u = 0
u . u dV
+ : dV
D
f . u dV
D2
T . u dS = 0
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
19
: dV = ?
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
20
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
D = De W = We
[K ] = [K e ]
e
{ F } = { Fe }
e
{U }
[ K ]{
U } = {F }
Systme global
Pour la statique
21
T2 1
Pb de temprature
2 nuds
T ( 0 ) = T1 T (1) = T2
T1 T ( s ) = [1 s ; s ] T2
22
Techniques numriques
Linaire
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
11
23
Techniques numriques
lments triangulaires
24
Techniques numriques
lments ttradriques
Les bases polynomiales sont compltes lments prismatiques bases incompltes lments hexadriques bases incompltes
25
Techniques numriques
Drf
s,t,u
Transformation gomtrique
Drel
x,y,z
{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
Intgration : on montre
f
De
(x, y, z)
dxdydz
Dref
f (s
t, u)
det [J ] dsdtdu
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 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]
27
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
(M )
28
Exemples dapplication
29
30
SNCF
Mannequin HYBRID III 50% dformable
Didier LEVEQUE
31
En biomcanique
Interface os - prothse Articulation du genou ligaments, tendons, cartilages, mnisques
32
Orange II
33
ECN - SNECMA
34
35
A vous de jouer
Soyez acteur de votre formation pour en profiter pleinement