2019-2020

C.C. Éléments-Finis final

ESIX MéSN2 - Master Mécanique 1
27 Novembre 2019 – durée : 90 minutes

Chaque candidat doit en début d’épreuve porter son nom dans le coin des copies. Il devra en outre
porter son numéro de place sur chacune de ses copies ou intercalaires ou pièces annexés. Le candidat
devra en outre éteindre son téléphone portable. Aucun document n’est autrorisé.

1 Questions de cours
On considère un espace de Hilbert séparable noté V et on considère un sous-espace affine fermé Vadm ⊂
V , on lui associe le sous espace vectoriel fermé V0 :
∀u, v ∈ Vadm , u − v ∈ V0 .
On considère a une forme bilinéaire définie sur V × V et L une forme bilinéaire définie sur V . On
considère le problème :
(
Trouver u ∈ Vadm tel que
(1)
a(u, u∗ ) = L(u∗ ) ∀u∗ ∈ V0
1. Sous quelles conditions satisfaites par la forme bilinéaire a et la forme linéaire L est ce que le
problème (1) est bien posé ?
2. Soit V h ⊂ V un sous-espace vectoriel fermé de V Ecrire le problème de Galerkin associé à (1)
dans V h .
3. Quel est le lien entre la solution de (1) et la solution du problème de Galerkin associé ?

2 Problème de Galerkin
On considère le problème variationnel suivant

 Trouver
 u(x) ∈ V0 = H01 ([0, 1]) tel que
Z Z
(2) 0 0

 u (x)v (x)dx = v(x)dx ∀v ∈ V0
[0,1] [0,1]
0
où u désigne la dérivée (au sens généralisé) de u. On admettra que problème (2) est bien posé. On
considère les fonctions suivantes
ϕ1 (x) = x(1 − x)
ϕ2 (x) = sin(πx)

1
 2019-2020 3 PROBLÈME 2D

1. Rappeler la définition de l’espace de Sobolev H01 ([0, 1]).

2. Montrer que les fonctions ϕ1 et ϕ2 sont bien admissibles dans Vadm .

3. Trouver les solutions de Galerkin du problème (2) dans les sous-espaces
(a) V1h = V ect{ϕ1 }
(b) V2h = V ect{ϕ2 }

4. Écrire (sans justifier) le problème aux limites équivallent à (2)

3 Problème 2D
Soit un domaine Ω ⊂ R2 . On considère un problème aux limites sur une fonction scalaire inconnue u

 −div(A∇u) = f
 Ω
(3) u = ud ΓD

(A∇u).n = 0 ΓN


où f est une fonction donnée de carré intégrable sur Ω, et u,n désigne la dérivée normale de u sur ΓN
et A est une matrice symétrique définie positive. Le bord de Ω est la réunion disjointe de ΓD et ΓN . Sur
ΓN , n désigne la normale unitaire extérieure.
1. Écrire la formulation variationnelle associée au problème (3) : on désignera par Vadm l’espace
affine des solutions admissibles et par V0 son espace vectoriel associé.

2. Donner un exemple non-nul de solution admissible du problème (3) dans le cas où ud = 0 et

Ω = {(x, y) ∈ R2 /x ∈ [0, a], y ∈ [0, b]}

ΓD = {(x, y) ∈ ∂Ω/x = 0 ou x = a}

avec a, b, h trois réels positifs non-nuls.

3. On considère dans toute la suite un maillage quadrangulaire à 4 noeuds sommets, conforme de Ω
: Soit V h l’espace vectoriel des fonctions continues polynômiales de degré 1 sur chaque élément
El de la discrétisation de Ω précédente:
m
[
Ω= El
l=1

et dont une base, constituée par les fonctions ”chapeau” ou fonctions de forme φi , est complètement
définie par les conditions

φi (Nj ) = δij , ∀i, j = 1, 2, . . . , n,

où δij est le symbole de Kronecker et les Nj désignent les n noeuds du maillage. Il s’agit donc
de l’élément isoparamétrique quadrangulaire à 4 noeuds (QUA4). On choisit de travailler sur un
élément de référence dont vous rappelerez la géométrie et dont les paramètres sont les variables r
et s.

2
 2019-2020 3 PROBLÈME 2D

4. Explicitez les expressions des fonctions de forme polynomiales de degré 1 Φi sur l’élément de
référence définies de façon analogue:

Φi (j) = δij ∀i, j = 1, 2, 3, 4

où on a noté 1,2,3,4, les quatres noeuds sommets définissant l’élément de référence.
5. Sur un élément El dont les sommets sont les noeuds {Nl1 , Nl2 , Nl3 , N4 }, qu’on notera pour
simplifier {1, 2, 3, 4}, de coordonnées :
   
x1 y1
x2  y2 
x̂ = x3  , ŷ = y3  .
  

x4 y4

Explicitez, la décomposition de u ∈ V h , en tant que fonction des variables r et s en fonction des

fonctions Φ1 , Φ2 , Φ3 , Φ4 et du vecteur û composé des valeurs de u aux noeuds {1, 2, 3, 4} :
 
u1
u2 
û = 
u3  .


u4

6. Déterminer la matrice Br telle que

∇r,s,t u = Br û,
∇r,s,t u est le gradient de u par rapport aux variables r, s, t.
7. Sur un élément El , trouver J telle que

∇r,s,t u = J∇x,y,z u.

où ∇x,y,z u est le gradient de u par rapport aux variables x, y, z .

8. Déduire, toujours sur El , en fonction des fonctions de forme Φi des coordonnées x̂, ŷ et ẑ les
matrices B et Bi telles que
∇u = B[û]
u = Bi [û]

9. En déduire les matrices élémentaires Kl et Ml telles que pour tout u, u∗ ∈ V h :

Z
∇u∇u∗ = û> Kl û∗
El
Z
uu∗ = û> Ml û∗
El

Préciser comment le calcul numérique des intégrales est effectuée.

3
 2019-2020 3 PROBLÈME 2D

10. Décrire l’assemblage des matrices Kl et Ml formant respectivement les matrices de ’rigidité’ K
et de ’masse’ M telles que
Z
∇u∇u∗ = û> K û∗
Ω
Z
uu∗ = û> M û∗
Ω

où û désigne désormais le vecteur colonne composé des valeurs de u en tous les noeuds du mail-
lage.

11. Déterminer R en fonction de fb et de M tel que ∀u∗ ∈ V h :

Z
f u∗ = R> û∗
Ω

où fb désigne le vecteur de Rn dont les composantes sont les valeurs de f aux n noeuds du
maillage.
12. Décrire précisément, en fonction de K, M , R, û, le problème (variationnel) à résoudre correspon-
h
dant à une solution éléments-finis de (3) dans Vadm = V h ∩ Vadm .
13. Soit LD la liste des noeuds du maillage appartenant au bord ΓD et soit L2 la liste des noeuds du
maillage n’appartenant pas à ΓD , autrement dit, L2 est la liste complémentaire de LD :

L2 ∪ LD = {1, 2, . . . , n}.

En introduisant les décompositions :

U1 = û[LD], U2 = û[L2],
R1 = R[LD], R2 = R[L2]

K11 = K[LD, LD] K12 = K[LD, L2]

K21 = K[L2, LD] K22 = K[L2, L2]

Déduire du problème variationnel, le système linéaire à résoudre pour obtenir la solution élément-
finis.
14. Proposer une solution analytique au problème (3), pouvant servir de référence pour valider un
code de calcul numérique.
15. Proposer une mesure de l’erreur de discrétisation par la méthode des éléments finis.