Mthode Des Elements Finis

Projet tutoré
Présenté par :
Bounejmate Sarah & Bouyahya Salma
Pour obtenir le diplôme de
Licence fondamentale : Sciences Mathématiques et

Applications
' $
Résolution d’équations aux dérivées

partielles par la méthode des éléments finis
& %
soutenue le : 25 Juin 2015

devant le jury composé de :
Mme.Touria Ghemires Président Professeur de l’enseignement supérieur à la FSR.
Mr.Abdellah Alla Encadrant Professeur de l’enseignement supérieur à la FSR.
Mme.Nadia Raissi Examinateur Professeur de l’enseignement supérieur à la FSR.
Année universitaire : 2014/2015

Table des figures
1 Géométrie à étudier (réf.[4]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2 Maillage de la géométrie à étudier (réf.[4]) . . . . . . . . . . . . . . . . . . . . 9
1.1 Diagramme expliquant la relation entre la dérivée au sens classique et la dé-

rivée au sens des distributions . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.1 Exemple d’éléments finis en 1 2 et 3 dimensions (réf.[9]) . . . . . . . . . . . . 33

3.2 Maillage en éléments finis (quadrilatères) d’un avion (réf.[7]) . . . . . . . . . 33
3.3 Fonction de base ϕi (réf.[1]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.1 Graphe représentant la solution exacte et la solution approchée par la mé-

thode des éléments finis de -u"(x)=f(x) (à gauche) et -u"(x)+u(x)=f(x) (à droite)
pour n=3 et u(x)=x(x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
pour n=10 et u(x)=x(x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
pour n=100 et u(x)=x(x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
4.4 Graphe représentant respectivement l’erreur pour -u"(x)=f(x) et -u"(x)+u(x)=f(x)
avec u(x)=x(x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
2
pour n=3 et u(x)=x(x-1) e(x−1/2) . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2
pour n=10 et u(x)=x(x-1) e(x−1/2) . . . . . . . . . . . . . . . . . . . . . . . . . 46
2
pour n=100 et u(x)=x(x-1) e(x−1/2) . . . . . . . . . . . . . . . . . . . . . . . . 47
4.8 Graphe représentant l’erreur respectivement pour -u"(x)=f(x) et -u"(x)+u(x)=f(x)
2
avec u(x)=x(x-1) e(x−1/2) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
pour n=3 et u(x)=x4 (x-1) ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
pour n=10 et u(x)=x4 (x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
pour n=100 et u(x)=x4 (x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
thode des éléments finis de -u"(x)+u(x)=f(x) pour n=100 et u(x)=x4 (x-1) . . . 49
2
TABLE DES FIGURES
4.13 Graphe représentant l’erreur respectivement pour -u"(x)=f(x) et -u"(x)+u(x)=f(x)

avec u(x)=x4 (x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3
Liste des tableaux
4.1 Tableau représentant l’erreur respectivement pour -u"(x)=f(x) et -u"(x)+u(x)=f(x)

avec u(x)=x(x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
2
avec u(x)=x(x-1) e(x−1/2) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
avec u(x)=x4 (x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
4
Table des matières
Remerciement 7
Introduction 8
1 L’espace de Sobolev H1 11
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.2 Rappel sur les espaces de Banach et de Hilbert . . . . . . . . . . . . . . . . . 11
1.3 Les distributions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.4 L’espace de Sobolev H1 (Ω) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.5 Notion de trace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.6 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2 Formulation variationnelle 20
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.2 Formulations variationnelles des problèmes modèles . . . . . . . . . . . . . . 20
2.2.1 Problème de Dirichlet homogène . . . . . . . . . . . . . . . . . . . . . 20
2.2.2 Problème de Dirichlet non homogène . . . . . . . . . . . . . . . . . . 21
2.2.3 Problème de Neumann . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.3 Théorème de Lax-Milgram . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.3.1 Application aux problèmes modèles . . . . . . . . . . . . . . . . . . . 27
2.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3 Méthode des éléments finis 30

3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.2 Approximation interne générale . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.3 La méthode des éléments finis . . . . . . . . . . . . . . . . . . . . . . . . . . 32
3.3.1 La méthode des éléments finis en dimension un . . . . . . . . . . . . . 33
3.3.2 Élément fini P1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
4 Applications 39
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.2 Algorithme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.3 Code Matlab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.4 Applications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
4.4.1 Mise en œuvre de la méthode des éléments finis en dimension 1 pour
u(x)=x(x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
2
u(x)=x(1-x) e(x−1/2) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
u(x)=x4 (x-1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
4.5 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
Conclusion et perspectives 50
5
TABLE DES MATIÈRES
Bibliographie 52
6
Remerciement
Avant tout, nous rendons grâce à Dieu le tout puissant pour l’énergie, l’esprit et l’endu-
rance sans lesquelles ce projet ne serait pas achevé.
Nous adressons nos sincères remerciements et nos profondes gratitudes à notre profes-
seur tuteur, Mr.Abdellah Alla, qui n’a pas cessé de nous pousser à aller de l’avant, pour
ses innombrables conseils et consignes, les documents qu’il a mis à notre disposition, son
temps, et la confiance qu’il nous a témoigné.
Nous tenons aussi à remercier Mme.Touria Ghemires, qui nous a honoré en acceptant
d’être président du jury.
Nous remercions également Mme.Nadia Raissi, qui nous a fait honneur en acceptant de
faire partie du jury en tant qu’examinateur.
Sans oublier de remercier toute l’équipe pédagogique du département de Mathématiques

de la Faculté des Sciences de Rabat, ainsi que l’ensemble des intervenants responsables de
la formation du cycle licence qui font de leur mieux pour que les étudiants aient une forma-
tion solide et parfaite.
Nous serons à jamais reconnaissantes à nos parents, nos proches, et nos amis pour leur
soutien inconditionné, et la confiance qu’ils nous ont accordé. Nous leur remercions pour
tout ce qu’ils ont bien voulu faire de nous.
7
Introduction
Depuis leur introduction en 1950, les équations aux dérivées partielles (en abrégé EDPs)
ont servi comme un outil de modélisation pour représenter analytiquement le compor-
tement dynamique de certains systèmes physiques. Aujourd’hui, l’application des EDPs
n’épargne presque aucun domaine : biologie, dynamique des populations, traitement de
signal, etc... d’où l’intérêt porté à leur résolution.
Certaines EDPs peuvent être résolues analytiquement. Toutefois, un nombre important
d’EDPs n’admettent pas de solutions analytiques. C’est dans cette optique que les recherches
se sont penchées sur les méthodes numériques pour arriver à approcher les solutions de ces
équations. Notons que malgré ces efforts indéniables, il n’existe pas de méthode universelle
pour la résolution numérique des EDPs. L’algorithme de résolution dépend très étroitement
du type du problème posé. C’est pour celà que nous allons restreindre notre champs d’étude
aux équations de type elliptique.
Les équations aux dérivées partielles sont regroupées suivant trois grandes catégories :
les EDPs paraboliques, hyperboliques, et elliptiques.
Pour mieux comprendre ce que celà veut dire, nous allons voir le cas des équations aux
dérivées partielles du deuxième ordre portant sur des fonctions réelles de deux variables
réelles u(x, y). Une telle équation s’écrit :
∂2 u ∂2 u ∂2 u ∂u ∂u
a 2 +b +c 2 +d +e + f u = g. (1)
∂x ∂x∂y ∂y ∂x ∂y
Pour simplifier, nous supposons que les coefficients a, b, c, d, e, f sont constants et g une
fonction continue.
On dit que l’équation (1) est elliptique si b2 − 4ac < 0, parabolique si b2 − 4ac = 0, et hyper-
bolique si b2 − 4ac > 0. L’origine de ce vocabulaire est bien sûr la classification des coniques
du plan. En effet, il est bien connu que l’équation du deuxième degré :
ax2 + bxy + cy 2 + dx + ey + f = 0,
définit une courbe plane qui est (sauf dans certains cas dégénérés) une ellipse si b2 − 4ac < 0,
une parabole si b2 − 4ac = 0, et une hyperbole si b2 − 4ac > 0 (réf. [6]).
Il existe plusieurs techniques permettant de résoudre les équations aux dérivées par-
tielles. On pense par exemple à la méthode des différences finies, Cette méthode travaille
directement sur une discrétisation de l’équation via une discrétisation du domaine. Mais
nous allons nous intéresser à une autre méthode dont l’esprit est un peu différent et plus
abstrait, c’est la méthode des éléments finis qui est basée essentiellement sur la discréti-
sation de l’espace des solutions cherchées, plutôt que le domaine sur lequel l’équation est
posée (bien qu’en définitive dans la méthode des éléments finis on verra que pour discréti-
ser l’espace, on utilisera une discrétisation du domaine). A part la méthode des différences
finies et des éléments finis, il existe bien d’autres méthodes d’approximation comme la mé-
thode du volumes finis, la méthode spectrale, etc... Mais, la méthode des éléments finis reste
8
INTRODUCTION
la plus largement répandue. Cette popularité est sans doute fondée. En effet, la méthode des
éléments finis est très générale, elle possède une base mathématique rigoureuse qui est fort
utile, elle permet de faire le traitement de plusieurs géométries complexes, et enfin, une
fois un algorithme de résolution est écrit, il peut être appliqué à n’importe quel problème
en changeant seulement les inputs. L’objectif de notre mémoire serait donc, d’introduire les
bases de la résolution des EDPs par la méthode des éléments finis (réf[3]).
Le principe de la méthode des éléments finis est simple, on transforme notre système
d’EDP continu en un système discret. La géométrie considérée initiale étant complexe est
subdivisée en un nombre fini de sous domaines simples appelés éléments dont l’assemblage
permet d’obtenir la géométrie initiale. Chaque élément est lié a son voisin par un nœud.
Considérons à titre d’exemple la géométrie suivante et supposons que l’on veuille calculer
sa surface :
Figure 1 – Géométrie à étudier (réf.[4])
Cette géométrie n’étant pas usuelle, on ne sait pas comment calculer directement sa sur-
face. La méthode des éléments finis propose qu’on la subdivise en un nombre fini d’objets
simples (carreaux, rectangles ..) pour lesquels le calcul de la surface ne pose pas de pro-
blèmes.
Figure 2 – Maillage de la géométrie à étudier (réf.[4])
9
INTRODUCTION
Pour arriver à notre objectif nous allons suivre le cheminement suivant : Dans un premier
temps, nous allons introduire quelques outils primordiaux de l’analyse mathématique. Ces
outils nous seront essentiels pour l’étude des équations aux dérivées partielles. Il s’agit, en
particulier, de la théorie des distributions et des espaces de Sobolev. Ensuite, dans le chapitre
2 nous allons parler de la formulation variationnelle, du théorème de Lax-Milgram. Nous
consacrons le troisième chapitre à la méthode des éléments finis que l’on va élaborer en
dimension 1, pour finir avec la simulation de celle ci sur Matlab pour quelques exemples
tests.
10
1 L’espace de Sobolev H1
1.1 Introduction
La majorité des phénomènes sont modélisés par le biais des équations aux dérivées par-
tielles. Il s’est avéré que l’espace des fonctions continument dérivables suffisamment de fois
pour donner un sens à l’EDP ne soit pas le cadre adéquat pour faire son étude défaut de
complétude d’où la nécessité d’introduire un nouvel espace complet appelé espace de Sobo-
lev H1 que nous allons définir dans ce qui suit.
L’objectif de ce chapitre est de donner les outils de base nécessaires à l’utilisation de

la méthode des éléments finis (expliquée dans le chapitre 3). Parmi les outils de base, on
retrouve les notions de distribution, des espaces de Hilbert et de Sobolev, et la notion
de trace. L’étude de ces notions pourrait faire l’objet d’un livre (et même de plusieurs) et il
va de soi que nous nous contenterons d’un survol assez rapide mais relativement complet.
Nous omettrons toutefois quelques détails qui ont bien sûr leur importance mais qui ne sont
pas essentiels pour une bonne compréhension de la méthode des éléments finis. Pour une
étude approfondie de ces notions nous renvoyons à ces deux livres [8] et [6] sur lesquels est
essentiellement basé notre rapport.
1.2 Rappel sur les espaces de Banach et de Hilbert

Espace de Banach :
Soit {un }n∈N une suite dans un espace normé (V , k.k). On dit que cette suite est de Cauchy
si :
lim kun − um k = 0.
n,m→+∞
L’espace normé (V , k.k) est dit de Banach si toute suite de Cauchy dans V converge vers un
élément de V (pour la norme k.k). En d’autres termes, un espace de Banach est un espace
normé complet.
Exemple :
Soit Ω un ouvert borné de Rn , l’espace V = C 0 (Ω)muni de la norme du maximum, est un
espace de Banach :
kvk∞ = max |v(x)|.
x∈Ω
Par contre, même si l’application k.kp : C 0 (Ω) −→ R (1 ≤ p ≤ ∞), définie par :

Z ! p1
kvkp = |v|p , ∀v ∈ v.
Ω
11
CHAPITRE 1. L’ESPACE DE SOBOLEV H1
est une norme sur C 0 (Ω) , l’espace (C 0 (Ω), k.kp ) n’est pas un espace de Banach.
Définition 1.2.1. Un produit scalaire dans un espace linéaire V est une forme (. , .) : V × V −→ R
avec les propriétés suivantes :
1. (u, u) ≥ 0, ∀u ∈ V et (u, u) = 0 ⇔ u = 0.
2. (u, v) = (v, u), ∀u, v ∈ V .
3. (αu + βv, w) = α(u, w) + β(v, w), ∀u, v, w ∈ V α, β ∈ R.
Le couple (V , (. , .)) est un espace linéaire muni d’un produit scalaire.
Théorème 1.2.1. (Inégalité de Cauchy-Schwarz). Soit (V , (. , .)) un espace linéaire muni d’un
produit scalaire. Alors on a :
p p
|(u, v)| ≤ (u, u) (v, v), ∀u, v ∈ V .
Espace de Hilbert :
Théorème 1.2.2. Soit (V , (. , .)) un espace linéaire muni d’un produit scalaire. Posons kvk =
p
(v, v) pour tout v ∈ V . Alors le couple (V , k.k) est un espace normé et
|(u, v)| ≤ kuk.kvk. ∀u, v ∈ V .
Un espace de Hilbert est un espace de Banach dont la norme est induite par un produit
scalaire.
Exemple :
L’espace (Rn , k.k2 ) est un espace de Hilbert.
 n 1
X  2 p
kxk2 =  xi2  = (x, x), ∀x = (xi )ni=1 ∈ Rn ,
i=1
où (. , .) est le produit scalaire dit euclidien.

Soit Ω un ouvert borné de Rn , les espaces L2 (Ω) et H1 (Ω) qu’on va introduire dans les
paragraphes suivants sont aussi des espaces de Hilbert.
1.3 Les distributions

Soient Ω un ouvert de Rn (n ≥ 1) et f une fonction f : Ω −→ R. On appelle support de f
l’ensemble :
suppf = {x ∈ Ω tel que : f (x) , 0}.
On dit que la fonction f est a support compact s’il existe un sous-ensemble compact K ⊂ Ω
tel que supp f ⊂ K. On définit ainsi l’espace D(Ω) par
D(Ω) = {f ∈ C ∞ (Ω) tel que : f est a support compact}.
12
On peut maintenant définir l’espace des distributions sur Ω.

Soit T une application linéaire de D(Ω) dans R. On note par < T , ϕ > la valeur de T cor-
respondante à l’élément ϕ ∈ D(Ω), c’est-à-dire < T , ϕ >= T (ϕ). On dit que T est continue
si : D E D E
lim T , ϕn = T , ϕ ,
n→+∞
où {ϕk }k∈N est une suite arbitraire de D(Ω) qui converge vers ϕ dans D(Ω).
On appelle distribution sur D(Ω) toute application T : D(Ω) −→ R linéaire et continue.
L’espace des distributions sur Ω est l’espace D0 (Ω), dual de D(Ω), c’est-a-dire, l’espace des
applications T : D(Ω) −→ R linéaires et continues (au sens spécifié ci-dessus).
Exemple :
Soit a un point de Ω. Le delta de Dirac relatif au point a, noté par δa , est la distribution
définie par : D E
δa , ϕ = ϕ(a), ∀ϕ ∈ D(Ω).
δa est linéaire. En effet,
D E
δa , ϕ + ψ = (ϕ + ψ)(a),
= ϕ(a) + ψ(a),
D E D E
= δa , ϕ + δa , ψ , ∀ϕ, ψ ∈ D(Ω).
Elle est aussi continue. En effet,

D E
lim δa , ϕn = lim ϕn (a),
n→+∞ n→+∞
= ϕ(a),
D E
= δa , ϕ , ∀ϕn ∈ D(Ω) tel que lim δn = δ.
n→+∞
On dit qu’une suite de distributions {Tn } converge dans D0 (Ω) vers une distribution T
si : D E D E
lim Tn , ϕ = T , ϕ , ∀ϕ ∈ D(Ω).
n→+∞
Considérons l’espace L2 (Ω) des fonctions à carré sommable défini par :

Z
2
L (Ω) = {f : Ω −→ R tel que |f (x)|2 dx < ∞}.
Ω
Cet espace est un espace de Hilbert.

Z
(f , g)L2 (Ω) = f (x)g(x) dx.
Ω
La norme de L2 (Ω) associée à ce produit scalaire est donc donnée par :

Z ! 21
21
2
kf kL2 (Ω) = (f , f )L2 (Ω) = |f (x)| dx .
Ω
13
A toute fonction f ∈ L2 (Ω) on associe une distribution Tf ∈ D0 (Ω) définie par :

D E Z
Tf , ϕ = f (x)ϕ(x) dx, ∀ϕ ∈ D(Ω). (1.1)
Ω
On a le résultat suivant :
Lemme 1.3.1. L’espace D(Ω) est dense dans L2 (Ω). En d’autres mots, toute fonction f ∈ L2 (Ω)
peut être approchée (pour la norme (k.kL2 (Ω) ) par des fonctions de D(Ω).
Soit T ∈ D0 (Ω) avec Ω un ouvert de Rn . La dérivée partielle de T par rapport à xi avec

i = (1, . . . , n) au sens des distributions est une nouvelle distribution définie de la manière
suivante :
D ∂T E D ∂ϕ E
,ϕ = − T , , ∀ϕ ∈ D(Ω).
∂xi ∂xi
De manière analogue on peut définir les dérivées successives. Pour chaque multi-indice
α = (α1 , . . . , αn ) ∈ Nn , on définit :
D E D E
D α T , ϕ = (−1)|α| T , D α ϕ ∀ϕ ∈ D(Ω).
Remarque 1.3.1. la dérivation au sens des distributions est une généralisation de la dérivée clas-
sique pour les fonctions, au sens suivant : si une fonction f est de classe C 1 , alors la dérivée de sa
distribution associée Tf coïncide avec la distribution Tf0 associée à la dérivée au sens classique f 0
de f . Cette relation est résumé dans le diagramme suivant (réf.[2]) :
Figure 1.1 – Diagramme expliquant la relation entre la dérivée au sens classique et la dérivée
au sens des distributions
1.4 L’espace de Sobolev H1(Ω)

Soit Ω un ouvert de Rn , et v une fonction de L2 (Ω), en général les fonctions intégrables
∂v
ou les fonctions de carré intégrable ne sont pas forcement dérivables, et donc ∂x n’aurait pas
i
de sens (au sens usuel). Mais pourtant, il va falloir considérer leurs dérivées pour que notre
EDP ait un sens. C’est ici ou l’on va se servir de la notion de dérivée faible(ou dérivée au sens
∂v
des distributions). En effet, ∂x a un sens au sens des distributions donnée par définition par
i
D ∂v E D ∂ϕ E
, ϕ = − v, pour tout ϕ ∈ D(Ω).
∂xi ∂xi
14
Définition 1.4.1. On appelle espace de Sobolev H1 (Ω) d’ordre 1 sur Ω (ouvert de Rn ), l’espace :
∂v
H1 (Ω) = { v ∈ L2 (Ω) tel que ∈ L2 (Ω), ∀i = 1, . . . , n },
∂xi
= { v ∈ L2 (Ω) tel que ∇(v) ∈ (L2 (Ω))n }.
On lui associe le produit scalaire :

Z  n

X ∂u ∂v 
(u, v)H1 (Ω) = uv +  dx.

Ω
 ∂xi ∂xi 
i=1
et la norme associée est notée :

Z  n !2  1/2
  X ∂v  
kvkH1 (Ω) = (v, v)1/2
1
H (Ω)
=  v 2 +  dx .
∂x
   
Ω i=1i
Définissons maintenant un autre espace de Sobolev qui est H10 (Ω), et qui nous sera très
utile pour les problèmes avec conditions aux limites.
Définition 1.4.2. L’espace de Sobolev H10 (Ω) est défini comme étant l’adhérence de D(Ω) dans
H1 (Ω).
Théorème 1.4.1. (Inégalité de Poincaré). Si Ω est borné, alors il existe une constante C =
C(Ω) > 0 telle que :
∀v ∈ H10 (Ω), kvkL2 (Ω) ≤ C(Ω)k∇vkL2 (Ω) . (1.2)
Démonstration :
On va montrer le résultat pour v ∈ D(Ω), puis on va conclure par densité de D(Ω) dans
1
H0 (Ω).
0 0
Puisque Ω est borné, alors on peut supposer que Ω ⊂ Ω × [a, b], avec Ω un ouvert de Rn-1 ,
0
on note x = (x , xn ) ∈ Rn-1 × R. On peut écrire que :
n o
Ω ⊂ (x0 , xn ) ∈ Rn−1 × R a ≤ xn ≤ b , a, b ∈ R .
Soit v ∈ D(Ω) et soit ṽ son prolongement par 0 en dehors de Ω, on a ṽ ∈ D(Rn ) et ṽ(x0 , a) = 0.

D’où : Z xn
∂ṽ 0
0
ṽ(x , xn ) = (x , t) dt ∀(x0 , xn ) ∈ Rn .
a ∂xn
On a aussi :
Z xn 2
0 2 ∂ṽ 0
|ṽ(x , xn )| = (x , t) dt
a ∂xn
Z xn !2
∂ṽ 0
≤ (x , t) dt .
a ∂xn
15
D’où par l’inégalité de Cauchy-Schwarz :
∂ṽ 0 2
xn
Z Z xn
0 2
|ṽ(x , xn )| ≤ (x , t) dt × |1| dt,
a ∂xn a
∂ṽ 0 2
Z xn
≤ (xn − a) × (x , t) dt,
a ∂xn
∂ṽ 0 2
Z +∞
≤ (xn − a) × (x , t) dt.
−∞ ∂xn
Ainsi : Z Z 2
0 2 0 ∂ṽ
|ṽ(x , xn )| dx ≤ (xn − a) × (x) dx.
Rn−1 Rn ∂xn
D’où en intégrant en xn , et ṽ(x0 , .) étant nul à l’extérieur de ]a, b[, on obtient :
Z Z b Z 2
0 2 ∂ṽ
|ṽ(x , xn )| dx ≤ (xn − a) dxn × (x) dx,
Rn a Rn ∂xn
Z 2
1 ∂ṽ
≤ (b − a)2 × (x) dx.
2 Rn ∂xn
Ainsi :
2
1 ∂v
kv(x)k2L2 (Ω) ≤ (b − a)2 × (x) .
2 ∂xn L2 (Ω)
b−a ∂v
=⇒ kvkL2 (Ω) ≤ √ × ∀v ∈ D(Ω).
2 ∂xn L2 (Ω)
L’inégalité de Poincaré étant vérifiée pour toute fonction de D(Ω), et on sait que toute
fonction de H10 (Ω) est une limite d’une suite de fonctions de D(Ω) (car D(Ω) est dense dans
H10 (Ω)), on en déduit alors le résultat.
Remarque 1.4.1. — Ce résultat nous montre en particulier que pour Ω borné, H10 (Ω) est un
sous-espace de H1 (Ω). Cependant, la fonction constante égale à 1 appartient à H1 (Ω) et
n’appartient pas à H10 (Ω), sinon elle vérifierait l’égalité ci-dessus ce qui est absurde car son
gradient est nul.
— H10 (Ω) est un espace de Hilbert.
Lemme 1.4.1. Si Ω est borné, la semi-norme
|v|H1 (Ω) = k∇vkL2 (Ω)

 n Z !2 1/2
X ∂v 
=  dx . (1.3)
Ω ∂xii=1
est une norme sur l’espace H10 (Ω) équivalente à la norme k.kH1 (Ω) .
16
Démonstration :
On se limite ici au cas de l’espace H10 (Ω). D’après (1.2), on tire que si |v|H1 (Ω) = 0 avec
v ∈ H10 (Ω), alors kvkL2 (Ω) = 0, ce qui implique que v = 0 puisque k.kL2 (Ω) est une norme. Il en
résulte donc que |.|H1 (Ω) est une norme sur H1 (Ω).
Soit v ∈ H10 (Ω). Puisque :

Z n !2  1/2


v 2 +
X ∂v  
kvkH1 (Ω) =   dx ,
∂x
  
Ω i=1 i
1
= kvk2L2 (Ω) + |v|2H1 (Ω) 2 ,
il est clair que :

|v|H1 (Ω) ≤ kvkH1 (Ω) .
D’autre part, d’après (1.2), on a :
1
kvkH1 (Ω) = kvk2L2 (Ω) + |v|2H1 (Ω) 2
1
≤ (1 + C 2 (Ω))|v|2H1 (Ω) 2
1
= (1 + C 2 (Ω) 2 |v|H1 (Ω) .
Ce qui complète la démonstration.
1.5 Notion de trace

Lors de la résolution des équations aux dérivées partielles, il nous faudra introduire des
conditions aux limites, ce qui nous amène à parler de la restriction à la frontière Γ = ∂Ω,
d’une fonction de H1 (Ω).
Une fonction v ∈ H1 (Ω) n’est pas nécessairement continue, et définir sa trace v|Γ (i.e sa ’va-
leur’) sur le bord Γ = ∂Ω pose un problème.
Les fonctions de H1 (Ω) ne sont pas nécessairement continues. Pour définir leur valeur
au bord on va les approcher par des fonctions régulières.
On admet que les fonctions de D(Rn ) sont denses dans H1 (Rn ). La démonstration de ce
résultat se fait par troncature et régularisation. D’autre part si Ω est un ouvert de classe C 1
par morceaux, c’est à dire que localement sa frontière est le graphe d’une fonction C 1 , ou
bien si Ω = Rn alors D(Ω) est dense dans H1 (Ω). Nous allons montrer que le fait d’associer
une valeur au bord pour une fonction est une opération continue pour la topologie de H1 .
Plus précisément, nous avons le théorème suivant :
Théorème 1.5.1. (théorème de trace) Soit Ω un ouvert borné de Rn , de frontière Γ de classe C 1

par morceaux.
Alors :

1. D Ω est dense dans H1 (Ω).
17
2. l’application restriction γ0 définie par :


0
D Ω −→ C (Γ )


γ0 : 
 v 7−→ γ0 (v) = v|Γ


est linéaire continue (pour D Ω muni de la norme k.kH1 (Ω) et C 0 (Γ ) muni de la norme k.kL2 (Γ ) )
et γ0 se prolonge par continuité de H1 (Ω) dans L2 (Γ ).
On définit ainsi : 
H1 (Ω) −→ L2 (Γ )


γ0 : 
 v 7−→ γ0 (v) = v|Γ

Et la continuité s’écrit :
∃c > 0, ∀v ∈ H1 (Ω) ; kγ0 (v)kL2 (Γ ) ≤ ckvkH1 (Ω) .
Remarque 1.5.1. le noyau de γ0 est en fait H10 (Ω).
Pour la démonstration de ce théorème nous aurons besoin du lemme suivant :

Lemme 1.5.1. Pour toute fonction v de D Rn+ , nous avons l’inégalité suivante :
kv(. , 0)kL2 (Rn-1 ) ≤ kvkH1 (Rn+ ) . (1.4)
Démonstration
:
Soit D Rn+ , alors on peut écrire :
Z ∞ !
0 2 ∂ 0 2
|v(x , 0)| = − |v(x , xn )| dxn ,
0 ∂xn
Z∞
∂v 0
= −2 v(x0 , xn ) (x , xn ) dxn ,
0 ∂xn
et en utilisant l’inégalité de Cauchy-Schwartz et l’inégalité suivant :
2ab ≤ a2 + b2 .
On obtient :
Z ∞ !1/2 Z ∞ 2 !1/2
0 2 0 2 ∂v 0
|v(x , 0)| ≤ 2 v(x , xn ) dxn (x , xn ) dxn ,
0 0 ∂xn
Z ∞ Z ∞ 2
0 2 ∂v 0
≤ v(x , xn ) dxn + (x , xn ) dxn .
0 0 ∂xn
Par intégration en x0 , on en déduit :

Z Z 2!
0 2 0 2∂v
|v(x , 0)| dx ≤ |v(x)| + (x) dx.
Rn-1 Rn+ ∂xn
18
D’où :
kv(. , 0)kL2 (Rn-1 ) ≤ kvkH1 (Rn+ ) .
Revenons maintenant à la démonstration du théorème (1.5.1).
Démonstration :
Pour démontrer le théorème (1.5.1), on se restreint au cas Ω = Rn+ , qui est le modèle le
plus simple d’ouverts de Rn de frontière C ∞ .
Rn+ := {x = (x0 , xn ); xn > 0}.
Dans ce cas, la frontière de Ω est l’hyperplan :
Γ = {x = (x0 , 0); x0 ∈ Rn-1 } Rn-1 .
Le cas d’un domaine général Ω est traité par cartes locales : on se place sur un morceau de
∂Ω et on le redresse grâce à un difféomorphisme en une partie de Rn sur lequel on applique
le résultat désormais connu. Puis on recolle les morceaux. n
Pour
revenir au cas Ω = R+ , on
va utiliser l’inégalité (1.4) du lemme (1.5.1). Comme D Rn+ est dense dans H1 (Rn+ ) cette
inégalité se prolonge aux fonctions de ce dernier espace et on obtient, pour ces fonctions,
une trace dans L2 du bord. A noter que ce n’est pas le cas pour les fonctions de L2 (Rn+ ) : pour
celles-ci on ne peut pas définir de valeur au bord (réf.[5]).
Corollaire 1.5.1. H10 (Ω) est le sous-espace de H1 (Ω) constitué des fonctions qui s’annulent sur
le bord Γ = ∂Ω et on a :
H10 (Ω) = { v ∈ H1 (Ω) tel que γ0 (v) = 0L2 (Γ ) }.
Démonstration :
Comme les fonctions de H10 sont des limites de fonctions de D(Ω) pour la topologie de
H10 et que celles-ci s’annulent sur ∂Ω, la continuité de γ0 implique que la fonction limite a
une trace nulle. L’inclusion dans un sens est donc prouvée. L’autre sens est admis.
Le résultat suivant permet de déduire la formule de Green.

Corollaire 1.5.2. (intégration par parties en dimension n) Soit Ω un ouvert borné de Rn , de
frontière Γ de classe C 1 par morceaux. Alors si u et v sont des fonctions de H1 (Ω), on a :
Z Z Z
∂u ∂v
v dx = − u v dx + uvni ds, i = 1, . . . , n
Ω ∂xi Ω ∂xi Γ
où ni est la i-ème composante de la normale sortante au domaine Ω.
1.6 Conclusion
Nous avons donc préparé le terrain, nous disposons maintenant de tous les outils néces-
saires pour élaborer la méthode des éléments finis. On a commencé par l’espace de Sobolev
pour assurer la complétude de l’espace dans lequel la solution approchée est définie, la no-
tion de dérivée au sens des distributions pour que l’EDP ait un sens, et enfin le théorème de
trace pour pouvoir donner des valeurs à la solution(qui n’est pas forcément continue) aux
bords.
19
2 Formulation variationnelle
2.1 Introduction
Nous abordons dans ce chapitre l’étude de certains problèmes aux limites elliptiques
linéaires intervenant de manière courante en mécanique et en physique. L’objet principal
est de trouver une formulation variationnelle de ces problèmes qui nous permet d’élaborer
une démonstration aisée de l’existence et de l’unicité des solutions. L’approche que nous
allons suivre est appelée une approche variationnelle.
Au cours de ce chapitre, l’exemple modèle d’EDP que nous allons traiter sera le Lapla-
cien, pour lequel nous étudions le problème aux limites suivant :

−∆u = f dans Ω,


 (2.1)
u = 0
 sur Γ .
où nous imposons des conditions aux limites de type Dirichlet. Dans (2.1), Ω est un ouvert
de l’espace Rn , Γ = ∂Ω est son bord (ou frontière), f est un second membre (une donnée du
problème), et u l’inconnue.
Le principe de l’approche variationnelle pour la résolution des équations aux déri-
vées partielles, est de remplacer l’équation du problème (2.1) par une formulation équi-
valente, dite variationnelle, obtenue en intégrant l’équation multipliée par une fonction
quelconque, qu’on appelle fonction test.
2.2 Formulations variationnelles des problèmes modèles

2.2.1 Problème de Dirichlet homogène
On appelle problème de Dirichlet une équation de Laplace avec des conditions aux li-
mites de type Dirichlet c’est à dire on spécifie les valeurs de la solution sur la frontière.
Pour tout Ω ouvert de Rn et de frontière Γ := ∂Ω le problème de Dirichlet s’énonce de la
façon suivante :
étant donné une fonction f ∈ L2 (Ω), trouver une fonction u définie dans Ω et solution de :



 (2.2)
u = 0
 sur Γ .
Supposons qu’il existe une solution u ∈ H2 (Ω) de (2.2) (une solution assez régulière).
Soit v∈ H10 (Ω), en multipliant la première équation de (2.2) par la fonction test v et en inté-
grant sur Ω, on obtient :
Z Z
−(∆u)v dx = f v dx.
Ω Ω
20
CHAPITRE 2. FORMULATION VARIATIONNELLE
Et en utilisant la formule de Green, on a :

Z Z Z
∂u
∇u.∇v dx − v dσ = f v dx.
Ω Γ ∂n Ω
comme v=0 sur le bord, on obtient ce que l’on appelle la formulation variationnelle du
problème (2.2) : 


Trouver u ∈ZH10 (Ω), tel queZ:
(2.3)

1
∀v ∈ H0 (Ω) ∇u.∇v dx = f v dx.



Ω Ω
Réciproquement, si u ∈ H10 (Ω)

est solution de (2.3), alors puisque
D(Ω) ⊂ H10 (Ω), on a pour tout ϕ ∈ D(Ω) :
Z n Z Z
X ∂u ∂ϕ
∇u.∇ϕ dx = dx = f ϕ dx.
Ω i=1 Ω ∂xi ∂xi Ω
cette égalité s’écrit encore :

n D
X ∂u ∂ϕ E D E
, = f ,ϕ ,
∂xi ∂xi
i=1
et par définition de la dérivation au sens des distributions il vient :

n D 2
X ∂ u E D E
− 2
, ϕ = f , ϕ ,
i=1
∂xi
c’est à dire :
−∆u = f .
La formulation variationnelle est donc équivalente pour les solutions régulières, à la formu-
lation classique du problème (2.2) de Dirichlet.
2.2.2 Problème de Dirichlet non homogène

Soit Ω un ouvert de Rn de frontière Γ := ∂Ω , on considère le problème suivant : étant
donné une fonction f ∈ L2 (Ω), trouver une fonction u définie dans Ω et solution de :



 (2.4)
u = d
 sur Γ .
Pour se ramener au problème de Dirichlet homogène, on va construire un relèvement, c’est-

à-dire une fonction u0 ∈ H1 (Ω) telle que γ0 (u0 ) = d où γ0 est l’application trace.
Si d ∈ Im(γ0 ) = γ0 (H1 (Ω)) qui est un sous-espace de L2 (Γ ), on sait d’après le théorème (1.5.1)
21
de trace qu’il existe u0 ∈ H1 (Ω) telle que d = γ0 (u0 ).

On cherche donc u sous la forme u=ũ+u0 avec ũ∈ H10 (Ω) et u0 ∈ H1 (Ω) telles que γ0 (u0 ) = d.
Soit v ∈ H10 (Ω), on multiplie la première équation de (2.4) par v et on intègre sur Ω on ob-
tient :
Z Z
−(∆u)v dx = f v dx.
Ω Ω
Comme u = ũ + u0 , on trouve :
Z Z Z
− ∆ũ.v dx − ∆u˜0 .v dx = f v dx.
Ω Ω Ω
D’où par la formule de Green, on obtient :

Z Z Z Z Z
∂ũ ∂ u0
∇ũ.∇v dx − v dσ + ∇u0 .∇v dx − v dσ = f v dx.
Ω Γ ∂n Ω Γ ∂n Ω
Et puisque ũ ∈ H10 (Ω), alors ũ = 0 sur Γ , et on se ramène donc au problème suivant :




 Trouver ũ ∈ H10Z(Ω), tel que Z:
Z
(2.5)
 ∇ũ.∇v dx = f v dx − ∇u0 .∇v dx ∀v ∈ H10 (Ω).



Ω Ω Ω
C’est la formulation variationnelle du problème (2.4) de Dirichlet non homogène.
2.2.3 Problème de Neumann

Toujours avec les mêmes hypothèses sur Ω, on cherche maintenant une fonction u telle
que :

−∆u + u = f dans Ω,



 ∂u (2.6)

 =0 sur Γ .
∂n

Supposons qu’il existe une solution u ∈ H2 (Ω) de (2.6) (u suffisamment régulière).

Soit v∈ H1 (Ω), en multipliant la première équation de (2.6) par la fonction test v et en inté-
grant sur Ω, on obtient :
Z Z Z
−(∆u)v dx + uv dx = f v dx,
Ω Ω Ω
et en utilisant la formule de Green, on a :

Z Z Z Z
∂u
∇u.∇v dx − v dσ + uv dx = f v dx,
Ω Γ ∂n Ω Ω
∂u
comme ∂n
= 0 sur le bord Γ , on obtient ce que l’on appelle la formulation variationnelle du
22
problème (2.6) de Neumann :


1


Trouver u ∈ HZ
0 (Ω), tel que :Z Z
(2.7)

1
∀v ∈ H (Ω) ∇u.∇v dx + uv dx = f v dx.



Ω Ω Ω
Réciproquement, si u∈ H2 (Ω) est solution du problème (2.7), alors en particulier on a en

prenant v ∈ D(Ω) dans (2.7), et on utilisant la formule de Green :
Z Z Z
−(∆u)v dx + uv dx = f v dx.
Ω Ω Ω
On a donc −∆u + u = f , au sens des distributions et comme f et u sont dans L2 (Ω), cette
égalité à lieu dans L2 (Ω).
Reprenons maintenant pour une fonction test v dans H1 (Ω) et intégrons par parties l’équa-
tion du problème (2.7), on obtient :
Z Z Z Z
∂u
−(∆u)v dx + uv dx − v dσ = f v dx.
Ω Ω Γ ∂n Ω
Z
∂u
D’après l’égalité obtenue précédemment, v = 0 pour tout v∈ H1 (Ω), on retrouve alors
Γ ∂n
la formulation variationnelle.
2.3 Théorème de Lax-Milgram

Cette section est consacrée au théorème de Lax-Milgram qui sera l’outil essentiel per-
mettant de démontrer des résultats d’existence et d’unicité de solutions de la formulation
variationnelle dans un espace de Hilbert. D E
On note V un espace de Hilbert réel de produit scalaire . , . et de norme k.k.
Nous considérons une formulation variationnelle du type :

 Trouver u ∈ V , tel que :


 (2.8)
a(u, v) = L(v) ∀v ∈ V .

On se donne :
1. L(.) est une forme linéaire continue sur V , c’est-à-dire que v −→ L(v) est linéaire de
V dans R, et il existe une constante C > 0 telle que :
|L(v)| ≤ Ckvk ∀v ∈ V .
2. a(., .) est une forme bilinéaire sur V , c’est-à-dire que u −→ a(u, v) est une forme li-
néaire de V dans R pour tout v ∈ V , et v −→ a(u, v) est une forme linéaire de V dans
R pour tout u ∈ V .
3. a(., .) est continue sur V × V , c’est-à-dire, qu’il existe une constante M > 0 telle que :
|a(u, v)| ≤ MkukV kvkV ∀u, v ∈ V .
23
4. a(., .) est V-elliptique (cœrcive), c’est-à-dire, qu’il existe une constante strictement
positive α telle que :
|a(v, v)| ≥ αkvk2V ∀v ∈ V .
Théorème 2.3.1. (Lax-Milgram) soit V un espace de Hilbert, et soit L une forme linéaire conti-
nue sur V, et soit a(. , .) une forme bilinéaire continue sur V ×V . Si a(. , .) est V-elliptique (cœrcive),
alors il existe une unique solution u du problème variationnel (2.8).
Démonstration :
On considère l’application suivante :
v 7−→ a(u, v).
Elle est linéaire continue donc d’après le théorème de Riesz on a :
∃A : V −→ V tel que a(u, v) = hAu, vi .
D’autre part, on considère l’application :
v 7−→ L(v).
Elle est aussi linéaire et continue, donc d’après le théorème de Riesz, on a :
∃τ : V 0 −→ V tel que L(v) = hτL, vi .
Alors :
hAu, vi = hτL, vi ∀v ∈ V ,
d’où :
Au = τL.
Montrons maintenant l’existence et l’unicité de la solution u.
On considère l’application suivante :
T : v 7−→ T (v) = v − ρ(Av − τL) avec ρ > 0.
Montrons que T est contractante.

Soit v1 , v2 ∈ V , on a :
kT (v1 ) − T (v2 )k2 = kv1 − v2 − ρ(Av1 − τL + ρ(Av2 − τL))k2

= kv1 − v2 − ρA(v1 − v2 ))k2 .
On pose, v = v1 − v2 , on obtient :
kT (v)k2 = kv − ρAvk2
D E
= kvk2 − 2ρ v, Av + ρ2 kAvk2 .
On sait que a(. , .) est V-elliptique, alors :

D E
v, Av = a(v, v) ≥ αkvk2 .
24
D E
=⇒ −2ρ v, Av ≤ −2ραkvk2 .
Puisque a(. , .) est continue, il existe donc une constante M > 0 telle que :
kAvk2 = a(v, v) ≤ M||v||2 .
D’où :
D E
kT (v)k2 ≤ kvk2 − 2ρ v, Av + ρ2 Mkvk2
≤ kvk2 − 2ραkvk2 + ρ2 Mkvk2 .
Ainsi :
kT (v)k2 ≤ (1 − 2ρα + ρ2 M) kvk2
T est une contraction si :
1 − 2ρα + ρ2 M < 1.
Ce qui implique que :
ρ(ρM − 2α) < 0
=⇒ ρM − 2α < 0
=⇒ ρ < 2α
M.
2α
Donc T est une contraction si ρ ∈]0, M [.
D’après le théorème de point fixe, il existe une unique solution u dans V qui vérifie :
T (u) = u.
Or
T (u) = u − ρ(Au − τL).
C’est-à-dire :
u = u − ρ(Au − τL).
D’où il existe une unique solution u telle que :
Au = τL.
Théorème 2.3.2. Sous les mêmes hypothèses du théorème (2.3.1) et si de plus la forme bilinéaire
a(., .) est symétrique, c’est-à-dire :
a(u, v) = a(v, u) ∀u, v ∈ V ,
nous avons :
Le problème variationnel (2.8) est équivalent au Problème de minimisation suivant :

Trouver une fonction u ∈ V telle que :



1 (2.9)
J(u) = min J(v) = min a(v, v) − L(v)



v∈V v∈V 2
25
Démonstration :
En premier lieu, on démontre que si u est l’unique solution du problème variationnel
(2.8), alors u minimise forcement la fonctionnelle J sur tout espace V, soit donc w ∈ V quel-
conque.
On va montrer que J(u + w) ≥ J(u). On a :
1
J(u + w) = a(u + w, u + w) − L(u + v).
2
Puisque a(., .) est bilinéaire, et L(.) est linéaire, nous avons :
1
J(u + w) = [a(u, u) + a(u, w) + a(w, u) + a(w, w)] − L(u) − l(w),
2
et la symétrie de a(., .) nous donne :
1 1

J(u + w) = a(u, u) − L(u) + (a(u, w) − L(w)) + a(w, w).
2 2
Dans le terme de droite, on reconnait dans la première parenthèse J(u), tandis que l’ex-
pression à l’intérieur de la deuxième parenthèse est nulle puisque u est solution du pro-
blème variationnel (2.8).
Et puisque a(., .) est V-elliptique alors le dernier terme est toujours positif.
On a donc :
1
J(u + w) = J(u) + a(w, w) ≥ J(u).
2
Donc J(u) est certainement inférieure ou égale à J(u +w). Puisque ce raisonnement est valide
quel que soit w ∈ V , u minimise bien la fonctionnelle J sur l’espace V.
Inversement, si u minimise J, on démontre que u est aussi une solution du problème
variationnel (2.8).
Considérons pour ce faire, la fonction de la variable réelle t définie par :
g(t) = J(u + tv),

1
= a(u + tv, u + tv) − L(u + tv).
2
En développant, on trouve :
1
g(t) = J(u) + t(a(u, v) − L(v)) + t 2 a(v, v).
2
On pose :
1
P (t) = t(a(u, v) − L(v)) + t 2 a(v, v).
2
C’est un polynôme de degré 2 en t.
Alors :
g(t) = J(u) + P (t).
Or :
J(u) ≤ J(u + tv) = g(t).
26
Donc :
P (t) ≥ 0,
et ceci n’est vrai, que si le discriminant D de P (t) est inférieur ou égal à 0.
C’est-à-dire :
D = (a(u, v) − L(v))2 ≤ 0.
Ceci implique que :
D = 0.
Ainsi :
a(u, v) = L(v).
D’où u est solution de (2.8).
2.3.1 Application aux problèmes modèles

• Problème de Dirichlet
On rappelle qu’il s’agit de trouver u ∈ H10 (Ω) tel que :
Z Z
1
∀v ∈ H0 (Ω) ∇u.∇v dx = f v dx.
Ω Ω
On applique donc le théorème (2.3.1) ci-dessus avec :

Z Z
a(u, v) = ∇u.∇v dx et L(v) = f v dx.
Ω Ω
Ceci sur l’espace V = H10 (Ω) muni de la norme :
kvkV = |v|H1 (Ω) = k∇vkL2 (Ω) .

0
• La forme bilinéaire a(., .) est continue sur H10 (Ω)×H10 (Ω), en effet, à l’aide de l’inégalité
de Cauchy-Schwarz, on a :
Z
|a(u, v)| = ∇u.∇v dx ,
ZΩ
≤ |∇u.∇v| dx,
Ω
Z !1/2 Z !1/2
2 2
≤ |∇u| dx |∇v| dx ,
Ω Ω
= k∇ukL2 (Ω) k∇vkL2 (Ω) ,
= kukV kvkV .
D’où :
|a(u, v)| ≤ kukV kvkV .
• D’autre part, a(., .) est V-elliptique, en effet :
27
Z
a(v, v) = (∇v)2 dx,
ZΩ
≤ |∇v|2 dx,
Ω
= k∇vkL2 (Ω) ,
= kvkV .
D’où :
a(v, v) = kvkV .
• Et enfin, L(.) est continue sur H10 (Ω), en effet à l’aide des inégalités de Cauchy-Schwarz
et de Poincaré, on a :
Z
|L(v)| = f v dx ,
Ω
Z !1/2 Z !1/2
|L(v)| ≤ |f |2 dx . |v|2 dx ,
Ω Ω
= kf kL2 (Ω) .kvkL2 (Ω) ,
≤ C(Ω)kf kL2 (Ω) .k∇vkL2 (Ω) .
On pose :
M = C(Ω)kf kL2 (Ω) .
D’où :
|L(v)| ≤ MkvkV .
Le théorème (2.3.1) de Lax-Milgram nous assure donc l’existence et l’unicité de la solution
u∈ H10 (Ω) du problème (2.2) de Dirichlet.
De plus, puisque la forme bilinéaire a(., .) est symétrique, le théorème (2.9), nous assure que
cette fonction u minimise la fonctionnelle quadratique suivante :
Z Z
1 2
J(v) = (∇v) dx − f v dx, sur l ’espace H10 (Ω).
2 Ω Ω
• Problème de Neumann
Rappelons qu’il s’agit de trouver u∈ H1 (Ω) tel que :
Z Z Z
1
∀v ∈ H (Ω) ∇u.∇v dx + uv dx = f v dx.
Ω Ω Ω
On applique donc le théorème (2.3.1) de Lax-Milgram, avec :

Z Z Z
a(u, v) = ∇u.∇v dx + uv dx, et L(v) = f v dx.
Ω Ω Ω
Ceci sur l’espace V = H1 (Ω) muni de la norme : kvkV = kvkH1 (Ω) .

• La forme bilinéaire a(., .) est continue sur H1 (Ω)×H1 (Ω), en effet à l’aide de l’inégalité
de Cauchy-Schwarz, on a :
28
Z Z
a(u, v) = ∇u.∇v dx + uv dx ≤ kukH1 (Ω) kvkH1 (Ω) .
Ω Ω
D’où
|a(u, v)| ≤ kukV kvkV .
• D’autre part a(., .) est V-elliptique. En effet,
Z Z
2
a(v, v) = (∇v) dx + (v)2 dx,
Ω Ω
= kvk2V .
D’où :
a(v, v) ≥ kvk2H1 (Ω) .
• Enfin L(.) est continue sur H1 (Ω), en effet à l’aide de l’inégalité de Poincaré, on a :
Z !1/2 Z !1/2
2 2
|L(v)| ≤ |f | dx . |v| dx ,
Ω Ω
= kf kL2 (Ω) .kvkL2 (Ω) ,
≤ C(Ω)kf kL2 (Ω) .kvkV .
En posant :
M = C(Ω)kf kL2 (Ω) .
On obtient :
|L(v)| ≤ MkvkV .
Le théorème (2.3.1) de Lax-Milgram nous assure donc l’existence et l’unicité de la solution
u ∈ H1 (Ω) du problème (2.6) de Neumann.
De plus, puisque la forme bilinéaire a(. , .) est symétrique, alors d’après le théorème (2.9), la
solution u minimise la fonctionnelle quadratique :
Z Z ! Z
1 2 2
J(v) = (∇v) dx + v dx − f v dx sur l ’espace H1 (Ω).
2 Ω Ω Ω
2.4 Conclusion
En se plaçant dans un cadre fonctionnel adéquat ; l’espace de Sobolev d’ordre 1, les mo-
dèles posés sous forme d’EDP avec des conditions aux limites, peuvent être écrits sous une
forme plus générale, dite forme faible ou variationnelle. Cette formulation permet de poser
un cadre conceptuel général qui nous fournit des outils d’analyse puissants, à savoir le théo-
rème de Lax-Milgram, ce théorème à la fois simple est puissant, nous assure l’existence et
l’unicité de la solution cherchée.
Lors de la transformation du système d’EDPs sous forme variationnelle, on suit la dé-
marche suivante : On choisit l’espace de Sobolev des fonctions tests en fonction du problème
considéré, On multiplie l’EDP par une fonction test, et enfin on intègre par parties.
29
3 Méthode des éléments finis
3.1 Introduction
Dans ce chapitre nous présentons la méthode des éléments finis qui est la méthode nu-
mérique de référence pour le calcul des solutions de problèmes aux limites elliptiques. Le
principe de cette méthode est directement issu de l’approche variationnelle que nous avons
vu dans le chapitre précédent. L’idée de base est de remplacer l’espace de Hilbert V , sur
lequel est posée la formulation variationnelle, par un sous-espace Vh de dimension finie. Le
problème approché posé sur Vh se ramène à la simple résolution d’un système linéaire, dont
la matrice est appelée matrice de rigidité. Par ailleurs, on peut choisir le mode de construc-
tion de Vh de manière à ce que le sous-espace Vh soit une bonne approximation de V, et que
la solution uh dans Vh de la formulation variationnelle, soit proche de la solution exacte u
dans V .
Remarque 3.1.1. Si le problème modèle est sous la forme : −∆u + u = f , on obtient au plus de la
matrice de rigidité, une autre matrice dite de masse.
Le plan de ce chapitre est le suivant. Dans la section (3.1) nous détaillons le processus
d’approximation variationnelle interne, et dans la section (3.2) nous présentons la méthode
des éléments finis en une dimension d’espace.
3.2 Approximation interne générale

Nous considérons à nouveau le cadre général de la formulation variationnelle introduit
dans le Chapitre 2.
Étant donné un espace de Hilbert V , une forme bilinéaire continue et cœrcive a(u, v), et une
forme linéaire continue L(v), on considère la formulation variationnelle suivante :

Trouver u ∈ V tel que :


 (3.1)
a(u, v) = L(v), ∀vh ∈ Vh ,

dont on sait qu’elle admet une unique solution par le Théorème (2.3.1). L’approximation
interne de (3.1) consiste à remplacer l’espace de Hilbert V par un sous-espace de dimension
finie Vh , c’est-à-dire à chercher la solution de :

Trouver uh ∈ Vh tel que :


 (3.2)
a(uh , vh ) = L(vh ), ∀vh ∈ Vh ,

La résolution de l’approximation interne (3.2) est facile comme le montre le lemme suivant.
30
CHAPITRE 3. MÉTHODE DES ÉLÉMENTS FINIS
Lemme 3.2.1. Soit V un espace de Hilbert réel, et Vh un sous-espace de dimension finie. Soit
a(u, v) une forme bilinéaire continue et cœrcive sur V , et L(v) une forme linéaire continue sur V .
Alors l’approximation interne (3.2) admet une unique solution. Par ailleurs cette solution peut
s’obtenir en résolvant un système linéaire de matrice définie positive et symétrique si a(u, v) est
symétrique.
Démonstration :
L’existence et l’unicité de uh ∈ Vh , solution de (3.2), découle du théorème (2.3.1) de Lax-
Milgram appliqué à Vh . Pour mettre le problème sous une forme plus simple, on introduit
Nh
X
une base (ϕj )1≤j≤Nh de Vh . Si uh = uj ϕj , on pose U = (u1 , . . . , uNh ) le vecteur dans RNh des
j=1
coordonnées de uh . Le problème (3.2) est équivalent à :
 
XNh 
Trouver U ∈ RNh tel que a  uj ϕj , ϕi  = L(ϕi ), ∀i = 1, . . . , Nh .
 
 
j=1
Ce qui s’écrit sous la forme d’un système linéaire :

AU = b. (3.3)
avec, pour 1 ≤ i, j ≤ Nh , (A)ij = a(ϕj , ϕi ), et (b)i = L(ϕi ).
Dans les applications mécaniques la matrice A est appelée matrice de rigidité. La
proposition suivante permet de montrer que la matrice A est inversible.
Proposition 3.2.1. La matrice A est définie positive c’est-à-dire qu’elle vérifie la propriété sui-
vante :
∀ U ∈ RNh − {0}, hAU , U i > 0.
Avant de donner la démonstration de cette proposition, expliquons pourquoi elle im-
plique l’inversibilité de la matrice A. En effet, si A n’est pas inversible, alors il existe un
vecteur U , 0 dans RNh tel que AU = 0. Ceci implique en particulier que hAU , U i = 0 avec
U , 0, ce qui contredit le fait que A est définie positive.
Démonstration :
Soit U = (u1 , . . . , uNh ) , 0 dans RNh . Le vecteur AU a pour coordonnées
   
X Nh Nh
X  X Nh Nh
X 
A1j uj , . . . , AN j uj  =  uj a(ϕj , ϕ1 ), . . . , uj a(ϕj , ϕNh ) .
   

   
j=1 j=1 j=1 j=1
On a donc
 
Nh X
X Nh X Nh Nh
X  Z1
hAU , U i = ui uj a(ϕj , ϕi ) = a  uj ϕj , ui ϕi  = a(uh , uh ) = |uh0 |2 dx,
 
i=1 j=1

j=1 i=1
 0
Nh
X
où uh (x) est une fonction de Vh dénie par uh (x) = uj ϕj (x).
j=1
Comme les (u1 , . . . , uNh ) sont tous non nuls, on vérifie aisément que uh est non nulle, ce qui
R1
implique que 0 |uh0 |2 > 0.
31
Nous allons maintenant comparer l’erreur commise en remplaçant l’espace V par son
sous-espace Vh . Plus précisément, nous allons majorer la différence ku − uh k où u est la solu-
tion dans V de (3.1) et uh elle dans Vh de (3.2).
Le lemme suivant, dû à Jean Céa, montre que la distance entre la solution exacte u et la solu-
tion approchée uh est majorée uniformément par rapport au sous-espace Vh par la distance
entre u et Vh .
Lemme 3.2.2. (lemme de Céa) Soit u une solution du problème (3.1) et uh une solution du
problème (3.2). Alors il existe une constante C > 0 indépendante de h tel que :
ku − uh kV ≤ C inf ku − vh kV .
vh ∈Vh
Démonstration :
Soit vh ∈ Vh quelconque, et wh = vh − uh ∈ Vh ⊂ V .
Comme wh est dans V et vh c’est une fonction test valide dans (3.1) et (3.2).
On a donc :
a(uh , wh ) = L(wh ).
Et
a(u, wh ) = L(wh ).
D’où :
a(u − uh , wh ) = 0.
Par cœrcitive et continuité de a(. , .), on a :
αku − uh k2 ≤ a(u − uh , u − uh ),
= a(u − uh , u − vh + wh ),
= a(u − uh , u − vh ) + a(u − uh , wh ),
≤ Mku − uh k ku − vh k.
Donc
M
ku − uh k ≤ ku − vh k.
α
M
et ceci pour tout vh , d’oú le résultat, avec C = α.
3.3 La méthode des éléments finis

Le principe de la méthode des éléments finis est de construire des espaces d’approxi-
mation interne Vh , des espaces fonctionnels usuels H1 (Ω), H10 (Ω), . . . , dont la définition est
basée sur la notion géométrique de maillage du domaine Ω. Un maillage est un recouvre-
ment de l’espace en volumes élémentaires très simples : triangles, tétraèdres, parallélépi-
pèdes.
Le paramètre h de Vh correspond à la taille maximale des mailles ou cellules qui com-
posent le maillage. Typiquement une base de Vh sera constituée de fonctions dont le support
est localisé sur une ou quelques mailles.
Ceci aura deux conséquences importantes : d’une part, quand h −→ 0, l’espace Vh sera de
plus en plus "grand" en terme de dimension et approchera de mieux en mieux l’espace V
32
Figure 3.1 – Exemple d’éléments finis en 1 2 et 3 dimensions (réf.[9])
tout entier, et d’autre part, la matrice de rigidité A sera creuse, c’est-à-dire que la plupart de
ses coefficients seront nuls (ce qui limitera le coût de la résolution numérique). La méthode
des éléments finis est une des méthodes les plus efficace et les plus populaire pour résoudre
numériquement des problèmes aux limites.
Figure 3.2 – Maillage en éléments finis (quadrilatères) d’un avion (réf.[7])
3.3.1 La méthode des éléments finis en dimension un

Nous choisissons le domaine Ω =]0, 1[. En dimension 1 un maillage est simplement
constitué d’une collection de points (xj )0≤j≤N +1 tels que
0 = x0 ≤ x1 ≤ x2 ≤ . . . ≤ xN ≤ xN +1 = 1.
Le maillage sera dit uniforme si les points xj sont équidistants, c’est-à-dire que :
1
xj = jh avec h= , 0 ≤ j ≤ N + 1.
N +1
Les points xj sont aussi appelés les sommets (ou nœuds) du maillage.
Nous considérons le problème modèle suivant :


−u”(x) = f (x) dans ]0, 1[,


 (3.4)
u(0) = u(1) = 0.

dont nous savons qu’il admet une solution unique dans H10 (Ω) si f ∈ L2 (Ω).
33
Dans ce qui suit on notera Pk l’ensemble des polynômes à coefficients réels d’une
variable réelle de degré inférieur ou égal à k :
 

 X k 

j
 
Pk =  p(x) = a x , a ∈ .
 
 j j R

 
 j=0 
Remarque 3.3.1. Généralement, le sous-espace Vh de V est noté

Vhk = {vh ∈ C 0 (Ω) tel que vh|K ∈ Pk (K)}
Sa dimension est égale à :
dim Vhk = (k + 1)(N + 1) − N
3.3.2 Élément fini P1

Pour k=1 et en notant Vh1 = Vh on a :
Vh = {v ∈ C 0 ([0, 1]), v|[xj ,xj+1 ] ∈ P1 , 0 ≤ j ≤ n, v(0) = v(1) = 0}. (3.5)

La méthode des éléments finis P1 est alors simplement une méthode d’approximation va-
riationnelle interne appliquée à l’espaces Vh définis par (3.5).
On peut représenter les fonctions de Vh , affines par morceaux, à l’aide d’une fonctions
de base très simple. Introduisons la "fonction chapeau" ϕj définie par :
 x−xj−1


 xj −xj−1 si x ∈ [xj−1 , xj ]

 xj+1 −x

ϕj (x) = 
 xj+1 −xj si x ∈ [xj , xj+1 ]


0

sinon
Lemme 3.3.1. L’espace Vh , défini par (3.5), est un sous-espace de H1 (]0, 1[) de dimension n + 2,
et toute fonction vh ∈ Vh est définie de manière unique par ses valeurs aux sommets (xj )0≤j≤N +1 :
N
X +1
vh (x) = vh (xj )ϕj (x) ∀x ∈ [0, 1].
j=0
Démonstration :
Les fonctions continues et de classe C 1 par morceaux appartiennent H1 (Ω). Donc Vh
est bien un sous-espaces de H1 (]0, 1[). Le reste de la preuve est immédiat en remarquant que
ϕj (xi ) = δij , où δij est le symbole de Kronecker qui vaut 1 si i = j et 0 sinon.
Corollaire 3.3.1. Les {ϕj }0≤j≤N +1 constituent une base de Vh .
Démonstration :
Puisque la famille {ϕj }0≤i≤N +1 est de dimension finie, et sa dimension est égale à la
dimension de Vh , on montre seulement que cette famille est libre.
Soient donc αj des scalaires dans R, avec j ∈ {0, . . . , N + 1}.
On suppose que :
N
X +1
αj ϕj (xi ) = 0 avec i ∈ {0, . . . , N + 1}.
j=0
34
• Pour i = 0, on a :
α0 ϕ0 (x0 ) + α1 ϕ1 (x0 ) + . . . + αN +1 ϕN +1 (x0 ) = 0 =⇒ α0 ϕ0 (x0 ) = 0,
d’où :
α0 = 0.
• Pour i = 1, on a :
α0 ϕ0 (x1 ) + α1 ϕ1 (x1 ) + . . . + αN +1 ϕN +1 (x1 ) = 0 =⇒ α1 ϕ1 (x1 ) = 0,
d’où :
α1 = 0.
..
.
• Pour i = N + 1, on a :
α0 ϕ0 (xN +1 ) + α1 ϕ1 (xN +1 ) + . . . + αN +1 ϕN +1 (xN +1 ) = 0 =⇒ αN +1 ϕN +1 (xN +1 ) = 0,
d’où :
αN +1 = 0.
Donc : ∀j ∈ {0, . . . , N + 1}, αj = 0.
Ainsi, {ϕj }0≤i≤N +1 est une famille libre.
Figure 3.3 – Fonction de base ϕi (réf.[1])
Décrivons la résolution pratique du problème de Dirichlet (3.4) par la méthode des élé-
ments finis P1 . La formulation variationnelle (3.2) de l’approximation interne devient ici :

Trouver uh ∈ Vh tel que :


(3.6)
 u (x)v 0 (x) dx = 1 f (x)vh (x) dx ∀vh ∈ Vh .
R1 R

 0
0 h h 0
On décompose uh sur la base des (ϕj )1≤j≤n et on prend vh = ϕi , ce qui donne :
N
X Z 1 Z 1
uh (xj ) ϕj0 (x)ϕi0 (x) dx = f (x)ϕi (x) dx.
j=1 0 0
35
R1
En notant U = (uh (xj ))1≤j≤N , b = ( 0 f (x)ϕi (x) dx)1≤i,j≤n , et en introduisant la matrice de
rigidité :
Z1 !
0 0
A= ϕj (x)ϕi (x) dx .
0 1≤i,j≤N
la formulation variationnelle dans Vh revient à résoudre dans Rn le système linéaire

suivant :
AU = b.
Comme les fonctions de base ϕj ont un "petit" support, l’intersection des supports de ϕj et
ϕi est souvent vide et la plupart des coefficients de A sont nuls. Les coefficients non nuls se
calculent facilement :
Z1 Z xi
0 0 (−1) 1 1
Ai,i−1 = a(ϕi−1 , ϕi ) = ϕi (x)ϕi−1 (x) dx = dx = − .
0 xi−1 h h h
Z 1 Z xi Z xi +1
1 1 2
Ai,i = a(ϕi , ϕi ) = (ϕi0 (x))2 dx
= 2
dx + dx = .
0 xi−1 h xi (−h)2 h
Z1 Z xi +1
0 0 1 (−1) 1
Ai,i+1 = a(ϕi+1 , ϕi ) = ϕi+1 (x)ϕi (x) dx = dx = − .
0 xi h h h
D’où : 


− 1h si j = i − 1 ou j = i + 1
2

Aij = 
 h si j = i (3.7)

0

si j , {i − 1, i, i + 1}.
Finalement, la matrice de rigidité est une matrice tridiagonale :
 2 −1 0 . . . 0 
 
−1 2 −1 0
 

 .
.

0 −1 2 −1 0 .
 
1  

A =  (3.8)

h  . . . . . . 
 . . . 0 
 . 
 .
 . 0 −1 2 −1

0 . . . 0 −1 2
 
Remarque 3.3.2. Si le problème modèle est sous la forme : −u”+u = f , on a au plus de la matrice
de rigidité (3.8), la matrice de masse qui s’écrit de la forme suivante :
4 1 0 . . . 0
 
1 4 1 0
 

 .. 
h 0 1 4 1 0 . 

M =  .. .. .. (3.9)
6 

. . . 0
 . 
 .
 . 0 1 4 1

0 ... 0 1 4
 
36
Pour obtenir le second membre b, il faut calculer l’intégrale :

Z xi +1
b= f (x)ϕi (x) dx pour tout 1 ≤ i ≤ N .
xi −1
L’évaluation exacte du second membre b peut être difficile ou impossible si la fonction f est
compliquée. En pratique on a recours à des formules de quadrature (ou formules d’intégra-
tion numérique) qui donnent une approximation des intégrales définissant b. Par exemple,
on peut utiliser la formule du point milieu :
Z xi+1
1 xi + 1
ψ(x) dx ≈ ψ( ).
xi+1 − xi xi xi
ou la formule des trapèzes :
Z xi+1
1 1
ψ(x) dx ≈ (ψ(xi+1 ) + ψ(xi )) .
xi+1 − xi xi 2
ou bien encore la formule de Simpson :
Z xi+1
1 1 1 2 x + xi
ψ(x) dx ≈ ψ(xi+1 ) + ψ(xi ) + ψ( i+1 ).
xi+1 − xi xi 6 6 3 2
On utilise la formule du trapèze, et on trouve :
Z xi
x − xi−1
f (x)ϕi (x)dx = i (f (xi )ϕi (xi ) + f (xi−1 )ϕi−1 (xi−1 )) ,
xi−1 2
h
= f (xi ).
2
et de la même manière on trouve :
Z xi+1
(xi+1 − xi )
f (x)ϕi (x)dx = f (xi+1 )ϕi (xi+1 ) + f (xi )ϕj (xj ) ,
xi 2
h
= f (xi ).
2
D’où
b = h × f (xi ).
Ainsi, le système à résoudre s’écrit :
 2 −1 0 . . . 0    
 
  u1   f (x1 ) 

−1 2 −1 0

    
 .   u2   f (x2 ) 
1  0 −1 2 −1 0
 ..   .   . 
  .   . 
2
 . . . . . .  . = . 
h   . . . 0   .   . 
 .   ..   .. 
 .
 . 0 −1 2 −1 u  f (x )
   
 N N
0 . . . 0 −1 2

On remarque que la matrice A est tridiagonale, symétrique et définie positive, alors le

système ci-dessus admet une unique solution.
37
3.4 Conclusion
La méthode des éléments finis 1D consiste donc à :
• Choisir N points entre 0 et 1 et choisir les fonctions ϕi .
• Construire la matrice de rigidité A.
• Déterminer le vecteur b (avec une méthode d’intégration).
• Résoudre le système linéaire AU = b où U désigne le vecteur des inconnus.
38
4 Applications
4.1 Introduction
Le but de ce chapitre est la mise en œuvre de la méthode des éléments finis en dimension
1 pour les deux problèmes modèles suivants :

−u”(x) + u(x) = f (x) sur ]0, 1[,



 u(0) = u(1) = 0.


−u”(x) = f (x) sur ]0, 1[,



 u(0) = u(1) = 0.

Nous allons ensuite pour les deux problèmes, comparer la solution obtenue par approxi-
mation, et la solution exacte. Nous utilisons le logiciel de calcul scientifique Matlab pour
illustrer la méthode des éléments finis.
4.2 Algorithme
Etape 1 : Discrétisation du domaine selon le degré de liberté fourni par l’utilisateur
(Maillage).
Etape 2 : Calcul de la matrice de masse et de rigidité.
Etape 3 : Calcul du second membre.
Etape 4 : Prise en compte des conditions aux limites.
Etape 5 : Résolution du système linéaire.
Etape 6 : Représentation graphique de la solution approchée.
4.3 Code Matlab

—— On cherche à résoudre : −u”(x) + u(x) = f (x) sur [0, 1] ——
—— La fonction u est donnée, deux fois dérivable et f continue ——
—— On fait d’abord le grand ménage ——
39
CHAPITRE 4. APPLICATIONS
clc
close all
clear all
dbstop if error ;
— On demande à l’utilisateur de saisir le nombre de nœuds désirés —

nddl=input(’donner le nombre de degré de liberté nddl : ’) ;
————– On saisit notre second membre f(x)=-u"(x)+u(x) ————–

f=@(x)(NOTRE SECOND MEMBRE f) ;
————– On saisit notre solution exacte ————–

sln = inline(’ON SAISIT LE u’) ;
————– On discrétise notre intervalle [0,1] ————–

xx=linspace(0,1,nddl) ;
————– On définit le pas de chaque maille ————–

for i=1 :nddl-1
h(i)=xx(i+1)-xx(i) ;
end
h;
————– On définit les conditions aux limites ————–

x0=0 ;
xL=1 ;
ua=0 ;
ub=0 ;
————– On définit notre matrice de rigidité ————–

A=zeros(nddl-2) ;
for i=1 :nddl-2
for j=1 :nddl-2
if(abs(i-j)>1)
A(i,j)=0 ;
elseif(abs(i-j)==0)
A(i,j)=(1./h(i)) + (1./h(i+1)) + (1./3)*(h(i)+h(i+1)) ;
else
A(i,j)=(-1./h(i)) + h(i)./6 ;
end
end
end
40
—— On calcule le second membre du problème approché ——
Bb=[] ;
for i=1 :nddl-2
Bb=[Bb 0.5*(h(i)+h(i+1))*f(xx(i+1))] ;
B=[Bb ] ;
end
————– On calcule la solution approchée ————–
u=[ua ; A B’ ; ub] ;
————– On dessine la solution exacte avec un ————–

————– graphe de couleur noire et d’épaisseur 3 ————–
figure(140) ;
plot(0 :0.01 :1,sln(0 :0.01 :1),’-.K’,’LineWidth’,3) ;
hold on ;
————– On dessine la solution approchée en marron ————–
plot(xx,u,’o-’,’Color’,[ 8.7059e-001 4.9020e-001 0],’LineWidth’,2) ;
————– On définit la légende et on nomme les axes ————–
Texte...legend(’bf sol exacte’,’bf sol approchée pour n=40’) ;

xlabel(’fontnameTimesfontsize14 bf x’,’Color’,[ 8.7059e-001 4.9020e-001 0]) ;
ylabel(’fontnameTimes fontsize14 bf u(x)’,’Color’,[ 8.7059e-001 4.9020e-001 0]) ;
grid ;
————– On définit les couleurs de l’interface graphique ————–
set(gcf,’Color’,[ 9.5294e-001 8.7059e-001 7.3333e-001]) ;

whitebg([ 9.5294e-001 8.7059e-001 7.3333e-001]) ;
41
4.4 Applications
Dans ce paragraphe, nous traitons les deux problèmes -u”(x)+u(x) = f (x) et −u”(x) = f (x)
sur [0,1] pour chacun des exemples suivants :
Exemple 1 :

u(x) = x(x − 1) sur ]0, 1[



u(0) = u(1) = 0

Exemple 2 :
 2
u(x) = x(1 − x)e(x−1/2) sur ]0, 1[



u(0) = u(1) = 0

Exemple 3 :

u(x) = x4 (1 − x) sur ]0, 1[



u(0) = u(1) = 0

Nous traçons la solution exacte et approchée pour chaque exemple, et regardons la

convergence pour chacune à partir du graphe.
42
4.4.1 Mise en œuvre de la méthode des éléments finis en dimension 1

pour u(x)=x(x-1)
Pour n=3 :
Figure 4.1 – Graphe représentant la solution exacte et la solution approchée par la mé-
thode des éléments finis de -u"(x)=f(x) (à gauche) et -u"(x)+u(x)=f(x) (à droite) pour n=3 et
u(x)=x(x-1)
Pour n=10 :
u(x)=x(x-1)
43
Pour n=100 :
Figure 4.3 – Graphe représentant la solution exacte et la solution approchée par la méthode
des éléments finis de -u"(x)=f(x) (à gauche) et -u"(x)+u(x)=f(x) (à droite) pour n=100 et
u(x)=x(x-1)
On sait que l’erreur vaut :
Z xn+1 !1/2
2
kErrkL2 = |u(x) − uh | dx .
x0
avec
n+1
X
uh = ϕi (x)ui .
i=0
Donc : 1/2
Z
n 2
 xn+1 X 
kErrkL2 =  u(x) − ϕi (x)ui dx .
 
 x 
0 0
Et on a :  x−xi−1

 x −xi−1 si x ∈ [xi−1 , xi ]
 xii+1 −x


ϕi (x) = 
 xi+1 −xi si x ∈ [xi , xi+1 ]

0

sinon.
Donc : Z n
xn+1 X n Z
X xi+1
ϕi (x)ui dx = ϕi (x)ui + ϕi+1 (x)ui+1 dx.
x0 i=0 i=0 xi
On se sert de cette formule pour calculer l’erreur pour ces exemples et on a :

————– On calcule l’erreur sur Matlab ————–
d=0 ;
for i=1 :nddl-2
g=@(x)abs(x.*(x-1)-u(i).*((x-xx(i+1))/(xx(i)-xx(i+1)))-u(i+1).*((x-xx(i))/(xx(i+1)-xx(i))))2 ;
d=d+quad(g,xx(i),xx(i+1)) ;
end
d=d1/2 ;
44
Figure 4.4 – Graphe représentant respectivement l’erreur pour -u"(x)=f(x) et -

u"(x)+u(x)=f(x) avec u(x)=x(x-1)
n 3 10 100
kErr1kL2 3.2275.10−2 2.1251.10−3 1.8534.10−5
kErr2kL2 2.9270.10−2 1.8448−3 1.6243.10−5
Table 4.1 – Tableau représentant l’erreur respectivement pour -u"(x)=f(x) et -

u"(x)+u(x)=f(x) avec u(x)=x(x-1)
Remarque :
Dans tout ce qui suit, Err1 représente l’erreur pour -u"(x)=f(x), et Err2 l’erreur pour -u"(x)+u(x)=f(x).
Commentaire : La solution u(x) est est régulière dans cet exemple, s’annule en 0 et 1 et
ne représente aucune singularité. A partir de n=10 la solution approchée converge vers la
solution exacte.
45

2
pour u(x)=x(1-x) e(x−1/2)
Pour n=3 :
2
u(x)=x(x-1) e(x−1/2)
Pour n=10 :
2
u(x)=x(x-1) e(x−1/2)
46
Pour n=100 :
2
u(x)=x(x-1) e(x−1/2)
Figure 4.8 – Graphe représentant l’erreur respectivement pour -u"(x)=f(x) et -

2
u"(x)+u(x)=f(x) avec u(x)=x(x-1) e(x−1/2)
n 3 10 100
kErr1kL2 6.0688.10−2 5.0346.10−3 4.4660.10−5
kErr2kL2 5.5407.10−2 4.4860.10−3 4.0162.10−5

2
u"(x)+u(x)=f(x) avec u(x)=x(x-1) e(x−1/2)
Commentaire : Il s’agit d’une fonction gaussienne centrée en 1/2, elle converge moins
rapidement que le premier exemple mais pour n=100 il est clair que la solution approchée
et exacte coïncident.
47

pour u(x)=x4 (x-1)
Pour n=3 :
u(x)=x4 (x-1) )
Pour n=10 :
u(x)=x4 (x-1)
48
Pour n=100 :
u(x)=x4 (x-1)
des éléments finis de -u"(x)+u(x)=f(x) pour n=100 et u(x)=x4 (x-1)
À première vue, il semble que la solution approchée coïncide en tous les points avec
la solution exacte. Mais, en agrandissant l’image on s’aperçoit que même pour un nombre
d’itérations très grand (ici n=100), la solution approchée ne vaut pas la solution exacte mais
elles s’approchent de plus en plus l’une de l’autre au fil des itérations avec une erreur à ne
pas oublier quoiqu’elle soit petite.
49
Figure 4.13 – Graphe représentant l’erreur respectivement pour -u"(x)=f(x) et -

u"(x)+u(x)=f(x) avec u(x)=x4 (x-1)
n 3 10 100
kErr1kL2 3.3002.10−2 3.9726.10−3 3.9657.10−5
kErr2kL2 2.9590.10−2 3.6269.10−3 3.7001.10−5

u"(x)+u(x)=f(x) avec u(x)=x4 (x-1)
Commentaire : Contrairement aux deux premiers exemples, dans ce cas u(x) représente
une singularité : gradient fort au voisinage de 1. On remarque que pour n petit la solu-
tion approchée est loin de la solution exacte et que celles-ci ne coïncident qu’à partir d’un
nombre considérable d’itérations n=100.
4.5 Conclusion
Au fils des itérations, on remarque que la solution approchée converge vers la solution
exacte.
La solution approchée converge plus rapidement pour -u"(x)+u(x)=f(x) que pour -u"(x)=u(x).
La solution exacte et la solution approchée ont une convergence de même ordre.
50
Conclusion et perspectives
La méthode des éléments finis est largement utilisée pour l’analyse des problèmes d’in-
génierie, non seulement en raison de son aptitude à résoudre une variété de problèmes phy-
siques formulés en termes d’équations aux dérivées partielles, mais aussi en raison de sa
capacité à gérer des géométries complexes et des conditions aux limites. Cependant cette
méthode représente quelques points faibles, qu’il serait faux de sous-estimer :
• La méthode des éléments finis s’appuie sur une formulation variationnelle de l’équa-
tion aux dérivées partielles, mais Il n’existe pas forcement de formulation variation-
nelle pour toute équation aux dérivées partielles.
• La méthode n’est pas bien adaptée à la résolution numérique d’équations non li-
néaires.
• Complexité de mise en œuvre de la méthode.
• Grand coût en temps de calcul et mémoire.
Il existe bien d’autres méthodes numériques de résolution d’équations aux dérivées par-
tielles comme les méthodes de volumes finis, d’éléments finis de frontière (ou méthode in-
tégrale),spectrale, de Fourier, etc.
Parlons par exemple de la méthode de volumes finis, c’est la méthode de choix pour les
équations de conservation non linéaires. Cette méthode consiste à intégrer, sur des volumes
élémentaires, les équations écrites sous forme intégrale. Elle est plus avantageuse que la
méthode des éléments finis puisque sa mise en œuvre est simple avec des volumes élémen-
taires rectangles, et elle permet de traiter des géométries plus complexes avec des volumes
de forme quelconque. C’est une méthode particulièrement bien adaptée aux équations de la
mécanique des fluides.
51
Bibliographie
[1] ALEXANDRE ERN, Aide mémoire : Éléments finis, Dunod, 2005.
[2] ALFIO QUARTERONI : Introduction à la méthode des éléments finis, École Polytechnique
fédération de Lausanne, 2013-2014.
[3] ANDRÉ FORTIN, ANDRÉ GARON :Les éléments finis : de la théorie à la pratique, École
Polytechnique de Montréal,1997-2006.
[4] Dr. H. JERRY QI : Finite element analysis, 2006.
[5] EMMANUEL MAITRE : Notes du Cours Éléments Finis - ENSIMAG 2A - Version du 21
février 2011.
[6] GRÉGOIRE ALLAIRE : Analyse numérique et optimisation. Éditions de l’école
polytechnique-Octobre 2012, 91128 Palaiseau Cedex.Méthode numérique.
[7] MARI VUORINEN et PEBTTI VARPASUO et JUKKA KäHKöNEN, Reaction-time
res- ponse of a large commercial aircraft using finite element method, Mai 2011,
www.kolumbus.fi (webographie)
[8] RAVIART, P.A et J.M.THOMAS : Introduction à l’analyse numérique des équations aux
dérivées partielles. Masson, Paris,1983.
[9] The Finite Element Method Theory, www.illustrations.marin.ntnu.no (webographie).
52

Mthode Des Elements Finis

Transféré par

Droits d'auteur :

Formats disponibles

Mthode Des Elements Finis

Transféré par

Informations du document

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

Mthode Des Elements Finis

Transféré par

Droits d'auteur :

Formats disponibles

Projet tutoré

Pour obtenir le diplôme de

Licence fondamentale : Sciences Mathématiques et

Résolution d’équations aux dérivées

soutenue le : 25 Juin 2015

Année universitaire : 2014/2015

1.1 Diagramme expliquant la relation entre la dérivée au sens classique et la dé-

3.1 Exemple d’éléments finis en 1 2 et 3 dimensions (réf.[9]) . . . . . . . . . . . . 33

4.1 Graphe représentant la solution exacte et la solution approchée par la mé-

4.13 Graphe représentant l’erreur respectivement pour -u"(x)=f(x) et -u"(x)+u(x)=f(x)

4.1 Tableau représentant l’erreur respectivement pour -u"(x)=f(x) et -u"(x)+u(x)=f(x)

3 Méthode des éléments finis 30

Sans oublier de remercier toute l’équipe pédagogique du département de Mathématiques

Figure 1 – Géométrie à étudier (réf.[4])

Figure 2 – Maillage de la géométrie à étudier (réf.[4])

L’objectif de ce chapitre est de donner les outils de base nécessaires à l’utilisation de

1.2 Rappel sur les espaces de Banach et de Hilbert

Par contre, même si l’application k.kp : C 0 (Ω) −→ R (1 ≤ p ≤ ∞), définie par :

|(u, v)| ≤ kuk.kvk. ∀u, v ∈ V .

où (. , .) est le produit scalaire dit euclidien.

1.3 Les distributions

D(Ω) = {f ∈ C ∞ (Ω) tel que : f est a support compact}.

On peut maintenant définir l’espace des distributions sur Ω.

Elle est aussi continue. En effet,

Considérons l’espace L2 (Ω) des fonctions à carré sommable défini par :

Cet espace est un espace de Hilbert.

La norme de L2 (Ω) associée à ce produit scalaire est donc donnée par :

A toute fonction f ∈ L2 (Ω) on associe une distribution Tf ∈ D0 (Ω) définie par :

Soit T ∈ D0 (Ω) avec Ω un ouvert de Rn . La dérivée partielle de T par rapport à xi avec

1.4 L’espace de Sobolev H1(Ω)

On lui associe le produit scalaire :

et la norme associée est notée :

Soit v ∈ D(Ω) et soit ṽ son prolongement par 0 en dehors de Ω, on a ṽ ∈ D(Rn ) et ṽ(x0 , a) = 0.

D’où par l’inégalité de Cauchy-Schwarz :

Lemme 1.4.1. Si Ω est borné, la semi-norme

|v|H1 (Ω) = k∇vkL2 (Ω)

il est clair que :

Ce qui complète la démonstration.

1.5 Notion de trace

Théorème 1.5.1. (théorème de trace) Soit Ω un ouvert borné de Rn , de frontière Γ de classe C 1

2. l’application restriction γ0 définie par :

∃c > 0, ∀v ∈ H1 (Ω) ; kγ0 (v)kL2 (Γ ) ≤ ckvkH1 (Ω) .

Remarque 1.5.1. le noyau de γ0 est en fait H10 (Ω).

Pour la démonstration de ce théorème nous aurons besoin du lemme suivant :

kv(. , 0)kL2 (Rn-1 ) ≤ kvkH1 (Rn+ ) . (1.4)

et en utilisant l’inégalité de Cauchy-Schwartz et l’inégalité suivant :

Par intégration en x0 , on en déduit :

Le résultat suivant permet de déduire la formule de Green.

où ni est la i-ème composante de la normale sortante au domaine Ω.

2.2 Formulations variationnelles des problèmes modèles

Et en utilisant la formule de Green, on a :

Réciproquement, si u ∈ H10 (Ω)

cette égalité s’écrit encore :

et par définition de la dérivation au sens des distributions il vient :

2.2.2 Problème de Dirichlet non homogène

Pour se ramener au problème de Dirichlet homogène, on va construire un relèvement, c’est-

de trace qu’il existe u0 ∈ H1 (Ω) telle que d = γ0 (u0 ).

D’où par la formule de Green, on obtient :

Et puisque ũ ∈ H10 (Ω), alors ũ = 0 sur Γ , et on se ramène donc au problème suivant :

Texte...legend(’bf sol exacte’,’bf sol approchée pour n=40’) ;