1
Une (trop ?) brve introduction au calcul par lments finis
F. Feyel , G. Cailletaud
[Link]@[Link] , [Link]@[Link]
F. Feyel, MEALOR : lments finis
Elments finis
Mthode variationnelle d'approximation de la solution de certaines quations aux drives partielles
F. Feyel, MEALOR : lments finis
3
[1] K.J. Bathe. Finite element procedures in engineering analysis. Prentice Hall, Inc., 1982. [2] J.L. Batoz and G. Dhatt. Modlisation des structures par lments finis, I III. Herms, 1991. [3] T. Belytschko, W.K. Liu, and Brian Moran. Nonlinear Finite Elements for Continua and Structures. 2000. [4] J. Besson, G. Cailletaud, J.-L. Chaboche, and S. Forest. Mecanique non linaire des matriaux. Herms, 2001. [5] G.R. Buchanan. Finite element analysis. Schaum s outlines, 1995. [6] P.G. Ciarlet and J.L. Lions. Handbook of Numerical Analysis : Finite Element Methods (Part 1), Numerical Methods for Solids (Part 2). North Holland, 1995. [7] M.A. Crisfield. Nonlinear Finite Element Analysis of Solids and Structures. Wiley, 1991. [8] G. Dhatt and G. Touzot. Une prsentation de la mthode des lments finis. Maloine, 1981. [9] T.J.R. Hughes. The finite element method: Linear static and dynamic finite element analysis. Prentice Hall Inc., 1987. [10] H. Kardestuncer, editor. Finite Element Handbook. Mc Graw Hill, 1987. [11] R. Mc Neal. Finite Element: their design and performance. Marcel Dekker, 1993. [12] O.C. Zienkiewicz and R.L. Taylor. The finite element method, Vol. I-III (Vol.1: The Basis, Vol.2: Solid Mechanics, Vol. 3: Fluid dynamics). Butterworth Heinemann, 2000.
F. Feyel, MEALOR : lments finis
4 +-*/ boucles
Mthodes numriques
Systmes non linaires Systmes diffrentiels Intgration
F. Feyel, MEALOR : lments finis
Systmes non linaires
F. Feyel, MEALOR : lments finis
Point fixe
Solution, erreur : Construction de la suite :
Converge ssi (drive contractante) Ordre 1 Ordre 2
F. Feyel, MEALOR : lments finis
Mthode de Newton
On chercher le zro :
Taylor :
Pour converger l'itration n+1, essayer :
Quasi-Newton : travailler K constant
F. Feyel, MEALOR : lments finis
8 Ordre de la mthode de Newton
Point fixe de
Drive (au voisinage de la solution)
Convergence (inconditionnelle pour les fonctions convexes) quadratique
F. Feyel, MEALOR : lments finis
9 Ordre de la mthode Quasi-Newton
Convergence linaire si
F. Feyel, MEALOR : lments finis
10
Newton : exemple deux variables
Solution
Iter. 1 2 3 4 5 6
x1 0.9166666790843010E+00 0.7959956739925678E+00 0.7862125392351869E+00 0.7861513801371453E+00 0.7861513777574233E+00 0.7861513777574233E+00
x2 0.6666666567325592E+00 0.6190476185745664E+00 0.6180344478073033E+00 0.6180339887499896E+00 0.6180339887498949E+00 0.6180339887498949E+00
f1(x) 0.2500000000000000E+00 0.2847222317424088E+00 0.1682906707772402E-01 0.9673552711856877E-04 0.3741760532709540E-08 0.5019856058607886E-16
f2(x) -0.7500000000000000E+00 0.1736111438108817E+00 0.1456149444031586E-01 0.9570904333691254E-04 0.3741548651090519E-08 -0.7123208273229764E-16
erreur 0.1305153013268777E+00 0.9844296235144490E-02 0.6116147776358272E-04 0.2379721975032112E-08 0.0000000000000000E+00 0.0000000000000000E+00
F. Feyel, MEALOR : lments finis
11
Systmes diffrentiels
F. Feyel, MEALOR : lments finis
12
Intgration Euler :
F. Feyel, MEALOR : lments finis
13
Euler implicite A
RK2
Euler implicite B
Euler point milieu
Euler
F. Feyel, MEALOR : lments finis
14
Runge-Kutta
Euler explicite point(s) intermdiaire(s)
Ordre 2 : rajouter une valuation au point milieu
F. Feyel, MEALOR : lments finis
15
Thta-mthode / Euler implicite
Type A
Newton-Raphson... Type A B
F. Feyel, MEALOR : lments finis
Type B
Rsidu
Jacobien
16
Intgration numrique
F. Feyel, MEALOR : lments finis
17
Choix des points et des poids ? Rectangles Trapzes Simpson
Gauss
F. Feyel, MEALOR : lments finis
18
Intgration de Gauss
Critre : intgration exacte des polynmes d'ordre 2r-1 Exemple l'ordre 2
Soit :
F. Feyel, MEALOR : lments finis
19 Un point d'intgration
5
-1 -1 -0.5 0 0.5 1
F. Feyel, MEALOR : lments finis
20 Deux points d'intgration
5
-1 -1 -0.5 0 0.5 1
F. Feyel, MEALOR : lments finis
Trois points d'intgration
5
21
-1 -1 -0.5 0 0.5 1
F. Feyel, MEALOR : lments finis
22
Intgration de Gauss en 3D
Souvent r1=r2=r3 Rgle spcifique aux triangles, prismes, ttras
F. Feyel, MEALOR : lments finis
23
Variations
F. Feyel, MEALOR : lments finis
24
Fonctionnelles
Fonctionnelle : fonction de fonctions , valeurs dans R Exemple classique : nergie potentielle Variation d'une fonctionnelle Objectif : trouver f qui minimise L(f) Variation (!= drive):
F. Feyel, MEALOR : lments finis
25
Proprits lmentaires
Fonction voisinage : Se comporte comme un oprateur diffrentiel
Commute avec les drives et intgrales, linaire, etc...
F. Feyel, MEALOR : lments finis
26
Euler-Lagrange...
Rsultat fondamental du calcul des variations : trouver f qui minimise
nergie potentielle
Est quivalent rsoudre
Avec et essentielles
F. Feyel, MEALOR : lments finis
naturelles
Essai de traction selon BTCdM-1742-2318
Soit une poutre encastre...
27
Energie potentielle
Equilibre des efforts
Essentielle
F. Feyel, MEALOR : lments finis
Naturelle
28
Autre point de vue formulation forte / faible
Choisir une fonction test admissible Multiplier l'quation et intgrer Intgrer par parties pour obtenir le mme ordre Formuler la forme faible
F. Feyel, MEALOR : lments finis
Application la poutre
Formulation forte :
29
Formulation faible (w : fonction rgulire qui tend vers zro l o les CL essentielles sont appliques) :
F. Feyel, MEALOR : lments finis
30
F. Feyel, MEALOR : lments finis
31
Si B est une forme symtrique, choisissons w=u, et posons
Euler-Lagrange !
tant donn une quation diffrentielle : Forme faible : principe des puissances virtuelles Forme quadratique : principe de l'nergie potentielle
F. Feyel, MEALOR : lments finis
Rsum : forme faible
32
Si u est une solution test qui vrifie :
Alors u est aussi solution de
F. Feyel, MEALOR : lments finis
Rsum : forme quadratique
33
Si u est une solution test qui vrifie :
Alors u est aussi le minimum de la fonctionnelle
F. Feyel, MEALOR : lments finis
Approximation !
Souvent :
34
Soit
avec Rigidit Force
F. Feyel, MEALOR : lments finis
35
Technologie lments finis Mcanique
F. Feyel, MEALOR : lments finis
Approximation lments par lments
Remplacer la formulation continue par une formulation discrte Inconnue Interpolation
36
Gomtrie
forme
Isoparamtrique :
F. Feyel, MEALOR : lments finis
37
Maillage ! Taille de maille : ABSOLUMENT rien voir avec la physique ! Prcision de l'approximation Type d'quations
F. Feyel, MEALOR : lments finis
Iso-paramtrique
Inconnue Gomtrie
38
Sous-paramtrique
Position des noeuds milieux !
(dplacements linaires)
Super-paramtrique
F. Feyel, MEALOR : lments finis
39
Gradient des fonctions de forme
Notation de Voigt symtrique
F. Feyel, MEALOR : lments finis
lment de rfrence
+1
40
-1
+1
-1 Un excellent mapping : les fonctions de forme !
F. Feyel, MEALOR : lments finis
41
F. Feyel, MEALOR : lments finis
42
On connat en fait (simplicit !) :
Or, en particulier :
Soit
F. Feyel, MEALOR : lments finis
43
Statique
F. Feyel, MEALOR : lments finis
44
Problme mcanique
Equilibre Forces imposes
Dplacement impos
F. Feyel, MEALOR : lments finis
45
Principe des puissances virtuelles
Galerkin discrtis
Admissible
F. Feyel, MEALOR : lments finis
Non-linaire !
46
Incrmental !
Newton-Raphson
Consistent
F. Feyel, MEALOR : lments finis
47
Exemple en lasticit :
+ Assemblage !
F. Feyel, MEALOR : lments finis
48
Algorithme gnral
F. Feyel, MEALOR : lments finis
49
Gags et autres bizarreries
F. Feyel, MEALOR : lments finis
50
Patch test
Retrouver une solution uniforme connue Laisser des noeuds libres
F. Feyel, MEALOR : lments finis
51
Modes nergie nulle
liminer les modes nergie nulle autres que rigides Exemple correct : rectangle quatre noeuds
Valeurs propres : { 223.4 90 78 46.36 42 0 0 0 }
F. Feyel, MEALOR : lments finis
52
lment trapzodal
(C. Felippa)
a a
2a
Rgle 1x1 2x2 3x3 4x4 8.77276 8.90944 8.91237 8.91246
Valeurs propres 3.68059 4.09769 4.11571 4.11627 2.26900 3.18565 3.19922 3.19966 0 2.64523 2.66438 2.66496 0 1.54678 1.56155 1.56199 000 000 000 000
F. Feyel, MEALOR : lments finis
53
Positivit du jacobien
Vrifier la position des noeuds milieux Les noeuds milieux au quart du cote produisent une singularit en J est ngatif si la matire sort de l'lment (convexit ?)
F. Feyel, MEALOR : lments finis
Modes parasites
4 3
54
Selon la rgle d'intgration choisie Linaires : 1x1 et 2x2 Quad : 3x3 et 2x2
F. Feyel, MEALOR : lments finis
Verrouillage
+1
55
-L
+L
-1 Incapacit d'un lment reproduire certains champs
F. Feyel, MEALOR : lments finis
56 On obtient :
Donc
F. Feyel, MEALOR : lments finis
