La Méthode des Éléments Finis: TD1
Ibrahim ALAME
01/10/2024
1 Problème de Dirichlet
On pose Ω =]0, 1[, V = H01 (Ω) et on considère le problème de Dirichlet
d2 u
− 2 + 4u = 4 sur Ω,
dx (1)
u(0) = u(1) = 0
1. Calculer explicitement la solution du problème (??). (La solution du problème étant connue
explicitement, la suite de l’exercice n’a évidement qu’un intérêt pédagogique)
2. Ramener l’étude du problème à un problème variationnel de type : Chercher u ∈ V solution
de
∀v ∈ V, a(u, v) = ℓ(v) (2)
où a(·, ·) est une forme bilinéaire et ℓ est une forme linéaire sur V que l’on déterminera.
3. On admet que le problème (??) a une solution u ∈ V et une seule. Pour construire une ap-
proximation uh de u, nous allons choisir un sous-espace de V constitué de fonctions continues
1
affines par intervalles. De façon plus précise, soit n un entier naturel et h = n+1 , à ce pas
h, nous associons les points xi = ih, 0 ≤ i ≤ n + 1 qui subdivisent l’intervalle Ω = [0, 1] en
n + 1 intervalles Ki = [xi , xi+1 ], 0 ≤ i ≤ n de longueur h. On choisit alors pour sous-espace
de dimension finie de V l’espace
Vh = {v ∈ V ; v|Ki ∈ P1 , 0 ≤ i ≤ n} = {v ∈ C 0 (Ω); v(0) = v(1) = 0; v|Ki ∈ P1 , 0 ≤ i ≤ n}
(3)
où P1 désigne l’espace des fonctions polynômes de degré inférieur ou égal à 1.
Déterminer la dimension de Vh et montrer que la suite des fonctions φi ∈ Vh , 1 ≤ i ≤ n
définies par φi (xj ) = δij , 1 ≤ i ≤ n constitue une base de Vh .
4. Donner l’expression analytique de φi (x) en fonction de i et x et tracer sa courbe représentative.
5. Pour tout i = 1, · · · , n on désigne par ui la valeur de la solution approchée uh au point xi :
ui = uh (xi ), 1≤i≤n (4)
Pn
Montrer que uh = j=1 uj φj , et que (uj )1≤i≤n sont solution du système linéaire
n
X
a(φj , φi )uj = ℓ(φi ), 1≤i≤n
j=1
1
6. Calculer a(φi , φj )en distinguant les trois cas : |i − j| > 1, |i − j| = 1 et i − j = 0.
7. Donner l’écriture matricielle du système linéaire.
8. Pour n = 4, Calculer et tracer une solution approchée du problème (??).
2 Problème de Neumann
On pose V = H 1 (0, 1) et on considère le problème de Neumann
d du
− (1 + x) + u = 1 sur ]0, 1[,
dx dx
(5)
′
u (0) = u′ (1) = 0
1. Montrer que le problème (??) se ramène à un problème variationnel de type : Chercher u ∈ V
solution de
∀v ∈ V, a(u, v) = ℓ(v) (6)
où a(·, ·) est une forme bilinéaire et ℓ est une forme linéaire sur V que l’on déterminera.
2. On subdivisent l’intervalle [0, 1] en n + 1 intervalles Ki = [xi , xi+1 ], 0 ≤ i ≤ n de longueur
1
h = n+1 . On choisit alors d’approcher V par le sous-espace de dimension finie de Vh suivant :
Vh = {v ∈ C 0 ([0, 1]); v|Ki ∈ P1 , 0 ≤ i ≤ n} (7)
où P1 désigne l’espace des fonctions polynômes de degré inférieur ou égal à 1. Les fonctions
φi , 0 ≤ i ≤ n + 1, définies par
1 − |x−x i|
si x ∈ [xi−1 , xi+1 ]
φ (x) = h 1≤i≤n
i
0 sinon
1 − hx si x ∈ [0, h]
φ0 (x) = (8)
0 sinon
1 − 1−x
si x ∈ [1 − h, 1]
φn+1 (x) = h
0 sinon
constituent une base de Vh . Les coordonnées dans cette base d’une fonction uh ∈ Vh sont les
Pn+1
nombres ui = uh (xi ), 0 ≤ i ≤ n + 1. On a donc uh = i=0 ui φi .
y
φ0 φj φn+1
x
0=x0 xj 1=xn+1
Montrer que (uj )0≤i≤n+1 sont solution du système linéaire
n+1
X
a(φj , φi )uj = ℓ(φi ), 0≤i≤n+1
j=0
3. Calculer a(φi , φj ) et ℓ(φi ).