Ecole MECANO « Mécanique des Nano-objets » Les Escandilles Autrans - 15 au 19 mars 2010 INTRODUCTION AUX ELEMENTS FINIS INTRODUCTION Modélisation physique et couplages Exemples d’utilisation LA METHODE DES ELEMENTS FINIS Équilibre, comportement et conditions aux limites Formulation faible Discrétisation, assemblage, résolution Analyse de la solution Exemple COMPLEMENTS Éléments finis isoparamétriques Non linéarités (matériau, géométrie, contact, …) Cinématiques particulières (poutres, plaques, …) R. Fortunier Ecole des Mines de Saint-Etienne Centre SMS, UMR CNRS 5146 [email protected]Poutre sous son poids - analytique - traité au tableau en parallèle fil rouge !
37
Embed
INTRODUCTION AUX ELEMENTS FINIS - im2np.fr · Éléments finis isoparamétriques Non linéarités (matériau, géométrie, contact, …) Cinématiques particulières (poutres, plaques,
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
INTRODUCTION AUXELEMENTS FINIS
INTRODUCTIONModélisation physique et couplages
Exemples d’utilisation
LA METHODE DES ELEMENTS FINISÉquilibre, comportement et conditions aux limites
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
NON-LINEARITES• λ, D, S ou L non constantes (changements de phase, élasticité non linéaire, plasticité, …)• conditions aux limites (rayonnement, contacts)• géométrie (en mécanique)• précipitation, réactions chimiques, diode, …
COUPLAGES MULTI-PHYSIQUES• λ, D, S ou L quantités fonctions d’inconnues d’une autre physique• couplage par le terme volumique (mise en forme, …)
LES ELEMENTS FINIS EN MECANIQUE :Statique, Linéaire, Petits déplacements, Petites déformations
EFFET DU TEMPS• loi de comportement avec effet du temps (fluage, …)• thermique ou diffusion non stationnaire• dynamique en mécanique
INTRODUCTIONModélisation physique et couplages
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
TRANSPORT
boitier de direction (PSA, ESI)
Roue de moto (ESI)
INTRODUCTIONExemples d’utilisation
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
train à grande vitesse italien (ESI)
TRANSPORT
INTRODUCTIONExemples d’utilisation
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
fermeture d’une capsule (Danone, ESI)
AGRO-ALIMENTAIRE
INTRODUCTIONExemples d’utilisation
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
AGRO-ALIMENTAIRE
fabrication d’une bouteille en verre (ESI)
INTRODUCTIONExemples d’utilisation
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
croissance cristalline(Cl. Maurice)
SCIENCE DES MATERIAUX
nitruration (P. Duranton et al., ESI)
INTRODUCTIONExemples d’utilisation
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
croissance cristalline (B. Serre, Zset)
GB movement
Zener pinning (G. Couturier)
SCIENCE DES MATERIAUX
INTRODUCTIONExemples d’utilisation
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
flexion 4P d’une éprouvette en cuivre (S. Villert)
SCIENCE DES MATERIAUX
INTRODUCTIONExemples d’utilisation
Injection de résine dans un renfort (G. Pacquaut, Zset) Frittage (J. Bruchon, CIMLIB)
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
EQUATIONS D’EQUILIBRE
LOI DE COMPORTEMENT ET CONDITIONS AUX LIMITES
FORMULATION FAIBLE
DISCRETISATION
ASSEMBLAGE
RESOLUTION
ANALYSE DE LA SOLUTION
EXEMPLE : ACTIONNEMENT D’UN MICRO INTERRUPTEUR
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
EQUATIONS D’EQUILIBRE
0ft dv ds Ω, ΩAA ΩΩ
A =+⊂∀ ∫∫∂
( ) 0 dv . Ω, ΩAΩ
A =+∇⊂∀ ∫ fσ
1,2,3)(j 0 f x
σ , ou . , j
3
1i i
ij ==+∂
∂Ω∈∀=+∇Ω∈∀ ∑
=
x0fσx
σ.n = tn
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
LOI DE COMPORTEMENT ET CONDITIONS AUX LIMITES
= : εLσ
( )∇+∇= 21 uuε t
∇= : uLσ
( ) 0fuLx : . , =+∇∇Ω∈∀ ( )
( ) T
UΩ pour . :
Ω pour pour : .
∂∈=∇=∂∈=
Ω∈=+∇∇
xTnuLtxUu
x0fuL
x
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
FORMULATION FAIBLE
( ) 0fuL : . , =+∇∇
fonctions test v(x)
( ) 0 dv . dv : . . ,E v =+∇∇∈∀ ∫∫ΩΩ
fvuLvv
0 dv :: dv . dS .:E pour tout que telE Trouver vu
=∇∇−+∈∈
∫∫∫ΩΩΩ∂
uLvfvTvvu
t
T
( ) ( ) U
1U
1
sur /H sur /H
Ω∂=Ω∈=Ω∂=Ω∈=
0vvUuu
v
u
EE
Conditions aux limites« naturelles »
Conditions aux limites« essentielles »
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
DISCRETISATION
∑∑==
==Ω∈∀
ee n
1JJ
eJ
n
1II
eI
e )(N )(et )(N )( , ee vxxvuxxux
∑∑==
∇=∇∇=∇Ω∈∀
ee n
1J
eJ
eJ
n
1I
eI
eI
e )(N )(et )(N )( , vxxvuxxux tt
x
y
Nie(x,y)
1
0 nœud i
nœud j
nœud k
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
nœud n1 / n1 x y z2 / n2 x y z…n / nn x y z
elements m1 / imat ne i1 i2 … ine …2 / imat ne i1 i2 … ine ……m / imat ne i1 i2 … ine …
numéro absolu (global)
numéro introduit (global)
numéro introduit du nœud numéro 2de l’élément 1 (local) MATRICES DE LOCALISATION Ae
DISCRETISATION
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
ASSEMBLAGE
∑∑==
==n
1LL
eJL
eJ
n
1KK
eIK
eI A et A vvuu
3numérotation locale
de l’élément 191
219
numérotation localede l’élément 6
1
2
3
6
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
RESOLUTION
∑∑
∑ ∫∫∑
∑ ∫∫∫∫
= =
= ΩΩ∩Ω∂=
= ΩΩ∩Ω∂ΩΩ∂
=
⎟⎟⎟
⎠
⎞
⎜⎜⎜
⎝
⎛+=
+=+
m
1e
eI
n
1I
eI
m
1e
eI
eI
n
1I
eI
m
1e
.
dv N dS N .
dv . dS . dv . dS .
e
e
Fv
fTv
fvTvfvTv
eeT
eeTT
∑∑∑
∑∑∑ ∫
∑∫∫
= = =
= = = Ω
= ΩΩ
=
⎟⎟⎟
⎠
⎞
⎜⎜⎜
⎝
⎛∇∇=
∇∇=∇∇
m
1e
n
1I
n
1J
eI
eIJ
eJ
m
1e
n
1I
n
1J
eI
eeteJ
m
1e
e e
e e
. .
. dv N . . N .
dv :: dv ::
uKv
uLv
uLvuLv
e
e
tt
0 dv :: dv . dS .:E pour tout que telE Trouver vu
=∇∇−+∈∈
∫∫∫ΩΩΩ∂
uLvfvTvvu
t
T
∑∑∑= ==
=><==m
1e
n
1I
eI
eIKK
n
1KKK
e
A avec . . FFFvFv [ ] ∑∑= =
=><=m
1e
n
1JI,
eJL
eIJ
eIKKL
e
A A avec . . v KKuKv
vecteur des sollicitations matrice de rigidité
tenseur des contraintes
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
[ ] ( ) 0 - .
:E pour tout que telE Trouver nvnu=><
∈><∈uKFv
vu
UJ
nvUI
nusur est J noeud le si / sur est I noeud le si /
Ω∂=><=Ω∂==
0vvUuu
EE
nœud I=KIJ uJ FI
uI = U
+P +PU
(pénalisation)
RESOLUTION
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
ANALYSE DE LA SOLUTION
VALEURS NODALESvecteur u contenant le « vecteur déplacement » uI de chaque nœud I.réactions concentrées aux nœuds : R = F – [K].u (résidu)
valeurs relativement précisestracé d’isovaleurs sur la structure
VALEURS DANS LES ELEMENTS (aux points d’intégration)Tenseur des contraintes et des déformations (gradient de déplacement)
Interpolation aux nœuds (moyennation)valeurs moins précises aux noeudstracé d’isovaleurs sur la structure
LA METHODE DESELEMENTS FINIS
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
Numerical simulations (Zset)
Parameters :• Cantilever thickness e• Contact thickness ec• Constant initial gap depth g0=1,7µm• Without stopper• (Hole position)• Influence of the meshActuation voltage :• 30V, 50V, 80VElectrostatic force :• The current actuation gap g is updated locally
(initial value g0) to evaluate the loading force
22
20
212
SVFg
ε=
contact
Cantilever
Subtrate
e
ec
g0 = 1,7µm
EXEMPLEMICRO-INTERRUPTEUR
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
Mechanical model Optimization
3
23
Initial gap(µm)
Actuated voltage (V)
Cantilever thickness e (µm)
Contact thickness ec(µm)
1.7 30 2 1
50 3 1.2
80 4 1.4
x
4
16
25
11 .
...
.
zy Notations :
example : e2 ec1 30V
Results :Cantilever final shape (along X axis)Mean contact pressure
EXEMPLEMICRO-INTERRUPTEUR
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
24
7
6
5
4
3
2
1
Case number Nb elts (thickness) Nodes Nb Elements Nb
3 1168 725
5 2766 1960
7 6631 5180
10 50112 43570
3 (quadratic) 4170 725
5 (quadratic) 10177 1960
7 (quadratic) 24966 5180
Influence of the mesh (case e4 ec0.7 50V)
EXEMPLEMICRO-INTERRUPTEUR
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
25
Influence of the mesh
case 4
case 3
case 1
EXEMPLEMICRO-INTERRUPTEUR
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010
Influence of the mesh - maximum deflection (U3)
.
At this point, end of the actuated area
Cases 5, 6, 7(quadratic elements)
Cases 1 to 4(linear elements)
EXEMPLEMICRO-INTERRUPTEUR
Ecole MECANO « Mécanique des Nano-objets »Les Escandilles Autrans - 15 au 19 mars 2010