Galerkin MEF1D

Cours Mthode des lments Finis
Mthode de Galerkin 1D
Berrada and Hadda
Ecole Nationale Suprieure des Arts et Mtiers de Mekns

Universit Moulay Ismail
Octobre 2017
Berrada and Hadda (ENSAM) Mthode des lments Finis Octobre 2017 1 / 32
Sommaire
1 Mthode de Galerkin
2 Principe de la MEF
3 Exemple en 1D
4 Technique dassemblage
5 Application numrique
6 Exercice
Sommaire
3 Exemple en 1D
6 Exercice
Mthode de Galerkin
Sous les hypothses de Lax-Milgram, soit u lunique solution de
(
Trouver u V tel que
(FV )
a(u, v ) = L(v ) v V
Soit Vh un sous-espace de V de dimension finie, on considre le problme

(
Trouver uh Vh tel que
(DS)
a(uh , v ) = L(v ) v Vh
Lobjectif
Calculer une solution approche uh de u dans un sous espace de dimension
finie Vh V
Vh est un ferm de V (car de dimension finie), et V est un Hilbert donc

Vh lest aussi. Do lexistence et lunicit de uh solution de (DS)
Mthode de Galerkin
Un espace Vh est construit laide dun maillage de

h la taille maximale des cellules qui composent le maillage (h>0)
La dimension de Vh reprsente le nombre de degrs de libert (les
inconnues)
Principe
Lorque h 0, la dimension de Vh devient de plus en plus gros et
approchera de mieux en mieux lespace V
Dfinition
Les espaces (Vh )h>0 forment une approximation interne de V si :
1 h > 0, Vh V
2 v V , (vh ) (Vh ) telle que kv vh kV 0
Mthode de Galerkin
Pour u et uh solutions de (FV ) et (DS) respectivement, on peut montrer

que a(u uh , v ) = 0, v Vh .
Pour a(, ) symtrique, il sagit dun produit scalaire sur V . La solution uh

peut tre alors interprt comme la projection orthogonale de u sur Vh au
sens de a(, ).
Lemme de Ca
On a lestimation suivante
M
ku uh kV inf ku v kV
v Vh
Sommaire
3 Exemple en 1D
6 Exercice
Principe de la MEF
Principe de la MEF
Subdiviser (discrtiser) une forme complexe en un grand nombre de
sous-domaines lmentaires de forme gomtrique simple (lments finis)
interconnects en des points appels nuds.
Un maillage Th est constitu

dlments nots K
Etapes de calcul par EF
Etape 1 : Maillage
Crer le maillage en dfinissant les nuds, les lments et les
polynmes dinterpolation
Etape 2 : Calcul
Calcul lmentaire
Pour chaque lment, tablir la matrice lmentaire et le second
membre lmentaire
Assemblage
Assembler les matrices et les vecteurs lmentaires
Conditions aux limites
Modifier le systme global en tenant compte des conditions aux limites
Systme linaire
Rsoudre le systme algbrique obtenu AX = b
Etape 3 : Post traitement
Sommaire
3 Exemple en 1D
6 Exercice
Exemple 1D
Soit le problme suivant

(
u 00 (x ) + u(x ) = f (x ), x ]0, 1[
u(0) = , u(1) =
Formulation variationnelle
(
Trouver u W tel que
R1 0 0
R1
0 (u (x )v (x ) + u(x )v (x )) dx = 0 f (x )v (x ) dx , v W0
o W = {v H 1 (]0, 1[), v (0) = et v (1) = } et W0 = H01 (]0, 1[)

On note Z 1
a(u, v ) = (u 0 (x )v 0 (x ) + u(x )v (x )) dx
0
et Z 1
L(v ) = f (x )v (x ) dx
0
Maillage 1D et espaces polynmiaux
Le maillage du domaine [0, 1] est une dcomposition en segments
successifs
n1
[0, 1] = i=1 Ki , Ki = [xi , xi+1 ]
avec les nuds 0 = x1 < x2 < < xn = 1
Ki
x
0 xi xi+1 1
Les lments du maillage sont dans ce cas les segments Ki
On notera Pk lespace vectoriel des polynmes de degr total k :
En 1D, Pk = vect{1, X , , X k } et dimPk = k + 1 ;
(k+1)(k+2)
En 2D, Pk = vect{X i Y j , 0 i + j k} et dimPk = 2 ;
En 3D, Pk = vect{X i Y j Z l , 0 i + j + l k} et
dimPk = (k+1)(k+2)(k+3)
6 ;
On notera Qk lespace vectoriel des polynmes de degr k par
rapport chaque variable :
En 1D, Qk = Pk ;
En 2D, Qk = vect{X i Y j , 0 i, j k} et dimQk = (k + 1)2 ;
En 3D, Qk = vect{X i Y j Z l , 0 i, j, l k} et dimQk = (k + 1)3 ;
Espace discrtis
On considre lespace P1 des polynmes 1D de degr 1
p P1 a, b R tels que p(x ) = ax + b
dimP1 = 2
On dfinit lespace discrtis de H 1 (]0, 1[)
Vh = {v H 1 (]0, 1[) / v|Ki P1 , i = 1, , n 1}
lespace des fonctions continues P1 par morceaux.
Lespace Vh est de dimension finie
Base nodale de type P1
On considre les fonctions de base (i )1in dfinies par

Les fonctions i sont P1 par morceaux i|Kj P1
La fonction i vaut 1 en xi et 0 sur les autres nuds i (xj ) = ij
1 supp(i ) = Ki1 Ki = [xi1 , xi+1 ]

Fonction chapeau
i|Ki1 (x ) = (x xi1 )/hi1
i|Ki (x ) = (xi+1 x )/hi
i|Kj (x ) = 0, (j 6 {i, i 1})
0 Xi-1 Xi Xi+1
Pn
La ieme composante de v dans cette base est vi = v (xi ) : v = i=1 vi i
La fonction j|Ki est dfinie par sa valeurs aux 2 nuds i = {xi , xi+1 }
Le triplet (Ki , i , P1 ) est appel lment fini de Lagrange 1D de type P1
Un lment fini de Lagrange 1D de type Pk est le triplet (Ki , i , Pk ) tel que
i = {xi + kl hi , l = 0, , k} (k + 1) nuds.
Ecriture matricielle
La fonction v et sa drive peuvent scrire alors

v1
v2
v (x ) = [1 2 n ] et v 0 (x ) = [01 02 0n ] v
| {z } | {z }
T vn BT
| {z }
v
Pour u, v Vh , on a
a(u, v ) = [0,1] u 0 v 0 + uv
R
R
= [0,1] [B T u]T [B T v ] + [T u]T [T v ] L(v ) = fv
R
"[0,1]
Z #
Z Z = fT v
[0,1]
= u T BB T + T v

[0,1] | {z }
[0,1]
LT
| {z } | {z }
S M
Introduction des conditions aux limites
Espace admissible
Wh = {v Vh / v (0) = et v (1) = }
{v = [v1 , v2 , , vn ] Rn / v1 = et vn = }
{v = [v1 , v2 , , vn ] Rn / C v = g}
o
1 0 0
C= g=
0 0 1
Espace associ Wh
Wh0 = {v Vh / C v = 0}
La forme matricielle

Trouver u Wh telle que
u T [S + M]v = LT v v Wh0
Introduction des conditions aux limites
On note L1 = {1, n} la liste des nuds correspondants une condition aux

limites impose et L2 la liste complmentaire (L1 L2 = {1, 2, , n})
On aura alors u(L1 ) = C u = g et v (L1 ) = 0 v Wh0
On introduit la matrice
N =S +M
On trouve alors
T T
u(L1 ) N(L1 , L1 ) N(L1 , L2 ) 0 L(L1 ) 0
=
u(L2 ) N(L2 , L1 ) N(L2 , L2 ) v (L2 ) L(L2 ) v (L2 )
Le systme simplifi devient
N(L2 , L2 )u(L2 ) + N(L2 , L1 )u(L1 ) = L(L2 )
N(L2 , L2 )u(L2 ) = L(L2 ) N(L2 , L1 )g
Sommaire
3 Exemple en 1D
6 Exercice
Technique dassemblage
Les matrices S et M sont calcules lment par lment en profitant du fait que
les fonctions de base sont nulles sauf sur deux lments.
On dcompose le calcul des matrices S et M

n1 Z
X
S +M = [BB T + T ]
i=1 Ki
n1
X
= [Si + Mi ]
i=1
o Si et Mi dsignent les matrices lmentaires

Z
Si = BB T
Ki
Z
Mi = T
Ki
Fonctions de base locales
Dans llment Ki = [xi , xi+1 ], seuls les fonctions i et i+1 sont non nuls
|Ki |Ki
1
Ki
x
0 xi xi+1 1
La fonction vectorielle |Ki est note i

0

0

0

0

0
i
h i h i
i

i = i i = et Bi = [Bi ] [Bi ] =

0i+1

i+1 0
0

0 0

h i xi+1 x 1
On aura alors i = h1i et [Bi ] = h1i
x xi 1
Matrices lmentaires
La matrice Mi est donne par

0 0 0
2
R R

Ki i Ki
i i+1
Mi = 0 [Mi ] [Mi ] = R
0
2i+1
R

Ki i i+1 Ki
0 0 0
La matrice Si est donne par

0 0 0
(0i )2 0i 0i+1
R R
Ki Ki

Si = 0
[Si ] 0 [Si ] = R

0 0 (0i+1 )2
R

Ki i i+1 Ki
0 0 0
Matrices lmentaires
Elment de rfrence K = [0, 1]
On introduit la transformation
Ti : K Ki
t 7 x = xi + thi
On dfinit les fonctions de base locales (1 , 2 ) de llment de

rfrence par
1 : t 7 i (Ti (t)) = 1 t
2 : t 7 i+1 (Ti (t)) = t
Ce qui revient aussi crire
xi+1 x
i|Ki : x 7 1 (Ti1 (x )) =
hi
x xi
i+1|Ki : x 7 2 (Ti1 (x )) =
hi
Matrices lmentaires
Les termes deRMi par changement de variable x = Ti (t)
Le terme Ki
2i
2i (x )dx = K 2i (Ti (t))Ti0 (t)dt
R R
Ki
= hi K 12 (t)dt
R
hi
=
3
de mme pour Ki 2i+1
R
Z Z Z
hi
2i+1 (x )dx = 2i+1 (Ti (t))Ti0 (t)dt = hi 22 (t)dt =
Ki K K 3
R
et pour Ki i i+1
Z Z
hi
i (x )i+1 (x )dx = hi 1 (t)2 (t)dt =
Ki K 6
La matrice [Mi ] est alors
" #
hi 2 1
[Mi ] =
6 1 2
Matrices lmentaires
De mme pour
R les termes de Si
Le terme Ki
(0i )2
(0i )2 (x )dx = K (0i )2 (Ti (t))Ti0 (t)dt
R R
Ki
= K h12 hi dt
R
i
= h1i
de mme pour Ki (0i+1 )2
R
Z Z Z
1 1
(0i+1 )2 (x )dx = (0i+1 )2 (Ti (t))Ti0 (t)dt = 2 hi dt =
Ki K K hi hi
0 0
R
et pour Ki i i+1
Z Z
1 1
0i (x )0i+1 (x )dx = 2 hi dt =
Ki K hi hi
La matrice [Si ] est alors
" #
1 1 1
[Si ] =
hi 1 1
Matrices globales
Les matrices S et M
1
h11

h1 0 0

h11 1
h1 + h2
1
h12 0 0

...
S= ...

1 1 1 1
0 hn2 hn2 + hn1 hn1

1 1
0 0 hn1 hn

h1 h1
3 6 0 0
h1 h1
+ h32 h1
0

6 3 6

...
M= ...

hn2 hn2 hn1 hn1
0 6 3 + 3 6

hn1 hn
0 0 6 3
Second terme
Le second terme lmentaire

R
Ki f i
[Li ] = R

Ki f i+1
Assemblage du second terme L

0
Z n1
X

L= f = Li , Li = [Li ]

[0,1]

i=1
0
Sommaire
3 Exemple en 1D
6 Exercice
Application numrique
On considre un maillage 3 lments de longueurs gales :

Les noeuds alors sont

0
1/3
noeud =

2/3

1
ce qui correspond x1 = 0, x2 = 1/3, x3 = 2/3, x4 = 1

Les lments
1 2
element = 2 3

3 4
ce qui correspond K1 = [x1 , x2 ], K2 = [x2 , x3 ], K3 = [x3 , x4 ]
Application numrique
Calculer les matrices S et M

On pose f = 1, calculer le second terme
Sommaire
3 Exemple en 1D
6 Exercice
Exercice 1D
On considre le problme suivant

(
u 00 (x ) + u(x ) = f (x ), x ]0, 1[
u(0) = , u 0 (1) =
Fin Mthode de Galerkin 1D

Galerkin MEF1D

Transféré par

Droits d'auteur :

Formats disponibles

Galerkin MEF1D

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Galerkin MEF1D

Transféré par

Droits d'auteur :

Formats disponibles

Cours Mthode des lments Finis

Berrada and Hadda

Ecole Nationale Suprieure des Arts et Mtiers de Mekns

Soit Vh un sous-espace de V de dimension finie, on considre le problme

Vh est un ferm de V (car de dimension finie), et V est un Hilbert donc

Un espace Vh est construit laide dun maillage de

Pour u et uh solutions de (FV ) et (DS) respectivement, on peut montrer

Pour a(, ) symtrique, il sagit dun produit scalaire sur V . La solution uh

Un maillage Th est constitu

Soit le problme suivant

o W = {v H 1 (]0, 1[), v (0) = et v (1) = } et W0 = H01 (]0, 1[)

On considre les fonctions de base (i )1in dfinies par

1 supp(i ) = Ki1 Ki = [xi1 , xi+1 ]

La fonction v et sa drive peuvent scrire alors

On note L1 = {1, n} la liste des nuds correspondants une condition aux

Le systme simplifi devient

N(L2 , L2 )u(L2 ) + N(L2 , L1 )u(L1 ) = L(L2 )

N(L2 , L2 )u(L2 ) = L(L2 ) N(L2 , L1 )g

On dcompose le calcul des matrices S et M

o Si et Mi dsignent les matrices lmentaires

La matrice Mi est donne par

La matrice Si est donne par

On dfinit les fonctions de base locales (1 , 2 ) de llment de

Le second terme lmentaire

Assemblage du second terme L

On considre un maillage 3 lments de longueurs gales :

ce qui correspond x1 = 0, x2 = 1/3, x3 = 2/3, x4 = 1

Calculer les matrices S et M

On considre le problme suivant

Vous aimerez peut-être aussi