MINISTÈRE DE L’ENSEIGNEMENT SUPÉRIEUR ET DE LA

RECHERCHE SCIENTIFIQUE
UNIVERSITÉ Abderrahmane MIRA de BEJAIA
FACULTÉ DES SCIENCES EXACTES
Département de Mathématiques

Mémoire présenté pour l’obtention du diplôme de Master

en Mathématiques

Option : Analyse et Probabilités

Par

MEHENAOUI Naima
KHALFOUNE Samia

THÈME

Etude variationnelle du système d’élasticité linéaire

Soutenu publiquement, le 25 /05/2017 devant le jury composé de :

Mme H. BECHIR M. C. A Université A. Mira de Béjaia Présidente

Mme S. TAS Prof Université A. Mira de Béjaia Promotrice
Mr B. KERAI M. A. A Université A. Mira de Béjaia Examinateur
1
 Remerciements

Tout d’abord, nous remercions Dieu, le tout puissant, de nous avoir aidé à accomplir ce
modeste travail.
Nous adressons nos chaleureux remerciements et témoignons notre profonde recon-
naissance à notre promotrice Mme S. TAS pour les conseils qu’elle n’a pas cessé de nous
prodiguer, ses orientations, ses judicieuses remarques et sa disponibilité.
Nous remercions vivement Mme H. Bechir pour l’honneur qu’elle nous fait en présidant
le jury de soutenance.
Nous tenons à remercier également Mr B. Kerai d’avoir bien voulu examiner ce travail
ainsi que pour son aide et ses conseils forts utiles.
En…n, nous témoignons notre gratitude à tous ceux qui nous ont aidé de près ou de
loin à réaliser ce mémoire.

i
 Dédicaces

Je dédie ce mémoire à :
Ma mère, pour son soutien, ses sacri…ces, ses précieux conseils, son assistance et sa
présence dans ma vie .
Mon père, pour les longues années de sacri…ce et de privation a…n de m’aider à avancer
dans la vie.
Un grand merci donc à mes parents pour les valeurs inculquées et leur soutien perma-
nent.
Ma soeur qui ne cesse d’être pour moi un exemple de persévérance, de courage et de
générosité.
Mon adorable frère Said.
Toute ma grande famille.
Tous mes amis.
Mes professeurs.
Toute ma promotion de 2017.
KHALFOUNE Samia

ii
 Dédicaces

Je dédie ce modeste travail à :

Ma mère qui a tant oeuvré pour ma réussite.

Mon père, pour son soutien, ses sacri…ces, ses précieux conseils et sa présence dans
ma vie.
Ma grand mère
Mes deux frères.
Tous mes amis et camarades de promotion 2017.

MEHENAOUI Naima

iii
 Table des matières

Introduction 2

Historique 4

1 Préliminaires et outils de base 5

1.1 Tenseurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.1 Pourquoi parle-t-on de tenseur ? . . . . . . . . . . . . . . . . . . . . 5
1.1.2 Dé…nition d’un tenseur . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.3 Tenseur des déformations . . . . . . . . . . . . . . . . . . . . . . . . 6
1.1.4 Tenseur des contraintes . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1.5 Tenseur d’élasticité . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2 Quelques rappels sur les espaces de Sobolev . . . . . . . . . . . . . . . . . 10
1.2.1 Espace L2 ( ) et (L2 ( ))n . . . . . . . . . . . . . . . . . . . . . . . 10
1.2.2 Espace de Sobolev H 1 ( ) . . . . . . . . . . . . . . . . . . . . . . . 11
1.2.3 Espace de Sobolev H m ( ) . . . . . . . . . . . . . . . . . . . . . . . 12
1.2.4 Espace de Sobolev H s (Rn ) . . . . . . . . . . . . . . . . . . . . . . . 13
1.2.5 Espace de Sobolev H01 ( ) . . . . . . . . . . . . . . . . . . . . . . . 15
1.2.6 Théorème de trace et formule de Green . . . . . . . . . . . . . . . . 16
1.2.7 Un résultat de compacité . . . . . . . . . . . . . . . . . . . . . . . . 17
1.3 Problèmes aux limites elliptiques linéaires . . . . . . . . . . . . . . . . . . 17
1.3.1 Formulation variationnelle . . . . . . . . . . . . . . . . . . . . . . . 17
1.3.2 Théorème de Lax-Milgram . . . . . . . . . . . . . . . . . . . . . . . 18

iv
 Table des matières

2 Modélisation du problème d’élasticité 19

2.1 Equations générales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19

3 Formulation variationnelle du problème de l’élasticité linéaire 23

3.1 Problème de Dirichlet homogène . . . . . . . . . . . . . . . . . . . . . . . . 23
3.2 Problème mixte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

Conclusion 48

Bibliographie 49

v
 Introduction

La modélisation mathématique a acquis une importance considérable, ces dernières

décennies, dans de nombreux domaines notamment dans les applications industrielles et
les sciences de l’ingénieur. Elle consiste à transformer un phénomène physique en des
modèles abstraits accessibles à l’analyse et au calcul. La modélisation d’un problème réel
utilise les lois de la physique (mécanique, thermodynamique, électromagnétisme, acous-
tique, etc.). Ces lois sont généralement, écrites sous forme de bilans qui se traduisent par
des Equations aux Dérivées Partielles (EDP).
Dans ce mémoire, nous nous intéressons à l’analyse mathématique de l’élasticité linéaire
qui est un système d’équations aux dérivées partielles de type elliptique correspondant à
un modèle physique stationnaire, c’est-à-dire indépendant du temps. Ce système a pour
objectif d’étudier le comportement mécanique des solides déformables, considérés comme
milieux continus, dont le matériau obéit à une loi constitutive linéaire réversible. En
outre, les déformations et les déplacements provoqués par des actions extérieures sont
considérés très faibles. Soulignons aussi que dans ce mémoire, nous nous attachons prin-
cipalement à étudier les e¤ets mécaniques (déformations, contraintes) se produisant dans
les solides en déformation, en nous préoccupant assez peu des phénomènes thermiques qui
peuvent se produire simultanément. De manière plus précise, nous allons supposer que la
température n’in‡ue pas sur les relations entre contraintes et déformations dans le solide.
Dans ces conditions, nous décrivons les déformations des solides élastiques au moyen de
la seule équation (mécanique) de conservation de la quantité de mouvement car elle est
tout à fait découplée de l’équation (thermique) de conservation de l’énergie.

1
 Introduction

Nous exposons, en premier lieu, les dé…nitions de di¤érents types de tenseur qui sont
primordiales pour la compréhension du problème d’élasticité ainsi que quelques rappels sur
les espaces de Sobolev et la formulation variationnelles de problèmes aux limites elliptiques
linéaires. Ensuite, nous montrons que le problème aux limites étudié est bien posé pour
cette équation aux dérivées partielles elliptiques, c’est-à-dire qu’elle admet une solution
unique qui dépend continûment des données. Nous montrons aussi que ces solutions
minimisent une énergie correspondante. Nous considérons d’abord des conditions aux
limites de type Dirichlet, puis des conditions mixtes. L’inégalité de Korn est l’ingrédient
essentiel dans la démonstration de la coercivité de la forme bilinéaire associée au problème.
Pour cela, nous nous appuyons sur l’approche variationnelle et nous faisons appel aux
espaces de Sobolev le cadre le plus approprié pour l’étude des EDP en général.

2
 Historique

C’est au début du XIX siècle qu’est née la théorie de l’élasticité, qui a rapidement attiré
l’attention de nombreux savants et mathématiciens.
Il y a un siècle et demi ou même moins, des traités tels que ceux de Saint-Venant, de
Love ou de Todhunter pouvaient contenir l’histoire du développement de cette discipline.
Dans cette dernière décennie, l’ouvrage d’Amy Dahan a permis de rafraîchir les mémoires
sur une partie de l’histoire de l’élasticité en France qui avait un lien étroit avec Cauchy.
Tout au long du XIXe siècle, la théorie de l’élasticité a connu une période d’intense
développement.
L’histoire de l’élasticité des solides est étroitement liée au développement des moyens
mathématiques. Ce lien a suscité la mise au point de nouveaux outils mathématiques.
Le développement de la théorie de l’élasticité est principalement dû à Navier, Poisson,
Cauchy, Lamé, Kirchho¤, Maxwell, Clebsch et Thomson.
La théorie de l’élasticité peut se réduire à l’établissement des relations fondamentales
liant ses di¤érents paramètres physiques et à l’étude des conditions aux limites à l’aide des
équations aux dérivées partielles, comme l’a démontré Cauchy. La théorie de l’élasticité
des corps n’est pas une théorie mathématique, mais une branche de la mécanique qui, à
son tour, représente une partie de la physique. Ce qui nous amène à tenir compte des
di¤érentes propriétés d’ordre pratique.
L’histoire des théories et des expériences, qui ont donné lieu au développement de
la théorie de l’élasticité des solides, présente un très grand intérêt. Plusieurs théories
physiques ont été utiles pour saisir la réalité du phénomène de l’élasticité. En e¤et, à

3
 Historique

travers les siècles, la théorie de l’élasticité des solides est passée par plusieurs étapes qui
ont contribué à son développement et à son accomplissement.
Dans les années 1820, plusieurs savants français vont aborder le domaine de l’élasticité
des solides. Il s’agit en premier lieu d’écrire les équations aux dérivées partielles qui dé…nis-
sent le mouvement des points d’un corps solide, c’est-à-dire son processus de déformation,
quand ce corps est soumis à l’action d’un système de forces extérieures ou intérieures, ou
quand il est en état de mouvement interne. La dé…nition de la notion de pression, on dit
aujourd’hui contrainte, se pose. On peut considérer que la théorie moderne de l’élasticité
est née en 1821, quand Navier donna pour la première fois, dans son célèbre mémoire, les
équations générales de l’équilibre et du mouvement d’un corps solide élastique. La plupart
des problèmes que pose cette théorie ont déjà été partiellement résolus ou discutés dans
des cas particuliers aux XVIIe et XVIIIe siècles ; de Galilée à Euler en passant par Hooke
et les Bernoulli, plusieurs résultats d’un grand intérêt avaient été obtenus. Seulement,
aucune théorie mathématique rigoureuse n’a été établie. Au XIX siècle, la loi de Hooke
fut le point de départ de plusieurs recherches qui débouchèrent sur la naissance d’une
vraie théorie mathématique de l’élasticité.

4
 1
CHAPITRE

Préliminaires et outils de
base

1.1 Tenseurs

1.1.1 Pourquoi parle-t-on de tenseur ?

Comment exprimer mathématiquement les contraintes qui s’exercent en chaque point

d’un milieu continu (solide, ‡uide, etc.) ? Les scalaires permettent certes d’exprimer des
pressions et les vecteurs des forces, mais ces entités mathématiques s’avèrent insu¢ santes.
Une contrainte en un point d’un milieu continu n’est pas une simple pression ni même
représentable par un unique vecteur force.
Le mot tenseur vient donc bien de l’idée d’exprimer, entre autres, des sortes de ten-
sions, des contraintes, etc. Mais la notion de tenseur se révèle d’une grande généralité,
avec des applications dans bien d’autres domaines que la seule mécanique des milieux
continus (électromagnétisme, relativité, géométrie di¤érentielle).

1.1.2 Dé…nition d’un tenseur

La notation tensorielle est une nomenclature, c’est-à-dire une façon systématique d’appeler
les variables et les champs d’intérêt physique et mathématique. Elle permet de faire

5
 1.1. Tenseurs

simplement des opérations sur nos équations de physique mathématique qui autrement
seraient beaucoup plus lourde.
Une grandeur bij à deux indices qui est appliquée à un vecteur cj , donne un autre
vecteur ai selon la loi
ai = bij cj

est un tenseur du second ordre. Un vecteur est encore appelé tenseur du premier
ordre et un scalaire, tenseur d’ordre 0. Plus généralement, et par récurrence, une
grandeur bi1 :::ik j1 :::jl à (k + l) indices qui, appliquée à un tenseur d’ordre l, cj1 :::jl , donne
un tenseur d’ordre k, ai1 :::ik selon la loi

ai1 :::ik = bi1 :::ik j1 :::jl cj1 :::jl

est un tenseur d’ordre (k + l):

Le plus simple des tenseurs du second ordre est le tenseur de Kronecker, dé…ni par
8
< 1 si i = j
ij =
: 0 si i 6= j

C’est bien un tenseur, car il applique un vecteur ai sur lui-même :

ai = ij aj

Remarque 1.1.1 Un scalaire, un vecteur, une forme linéaire, une forme bilinéaire, une
application linéaire sont des tenseurs particuliers.

1.1.3 Tenseur des déformations

Les déformations d’un objet sont mesurées à partir d’une position initiale, qui est générale-
ment la position de repos dans laquelle aucune force n’est appliquée à l’objet.
Le champ vectoriel déplacement, généralement noté u, est, par dé…nition pour chaque
point, le vecteur reliant sa position au repos à sa position actuelle dans la con…guration
déformée.

6
 1.1. Tenseurs

On considère les premières dérivées spatiales de ce champ u à travers la matrice A qui

est son gradient : 0 1
@ux1 @ux1 @ux1
: :
B @x1 @x2 @xn C
B @ux2 @ux2 @ux2 C
B : : C
B @x1 @x2 @xn C
B C
A = grad u = B : : : : : C
B C
B C
B : : : : : C
@ A
@uxn @uxn @uxn
@x1 @x2
: : @xn

A peut être écrite sous la forme d’une somme de deux matrices n n:

1 1
A = grad u = (grad u + grad uT ) + (grad u grad uT )
|2 {z } |2 {z }
E R

E est une matrice symétrique tandis que R est antisymétrique et est liée au rotationnel
rot u du champ de déplacement.

Dé…nition 1.1.1 La matrice E représente les déformations du matériau et dé…nit le

tenseur des déformations noté ".
On a donc
0 1
@ux1 1 @ux1 @u 1 @ux1 @uxn
( + @xx12 ) : : ( + )
B @x1 2 @x2 2 @xn @x1 C
B 1 @ux2 @ux1 @ux2 1 @ux2 @uxn C
B ( + ) : : ( + ) C
B 2 @x1 @x2 @x2 2 @xn @x2 C
B C
"=B : : : : : C
B C
B C
B : : : : : C
@ A
1 @uxn @ux1 1 @uxn @ux2 @uxn
(
2 @x1
+ @xn
) (
2 @x2
+ @xn
) : : @xn

Les termes de la matrice peuvent s’écrire :

1 @ui @uj
(")ij = ( + ); 1 i; j n
2 @xj @xi

1.1.4 Tenseur des contraintes

Après la description des déformations, nous introduisons le tenseur des contraintes qui va
décrire en chaque point la répartition des forces internes.

7
 1.1. Tenseurs

Les contraintes représentent les e¤orts de cohésion dans un solide qui permettent à la
matière de résister aux sollicitations. Si nous prenons un cube de matériel et le soumettons
à une charge arbitraire, nous pouvons mesurer la contrainte sur plusieurs directions. Ces
mesures forment un tenseur de second degré, le tenseur des contraintes.

Vecteur contrainte et tenseur des contraintes

Le vecteur contrainte caractérise les e¤orts de contact exercés à travers un élément de

~ sur une partie D du milieu continu : le vecteur contrainte est
surface dS de normale N
dé…ni par
~
~ ) = lim df
T~ (N ; df~ = T~ (N
~ )dD:
dS!D dS

Pour connaître l’état de contrainte en un point donné, il faut connaître les vecteurs
~.
contraintes associés à toutes les facettes, c’est-à-dire à tout vecteur unitaire N
~ à T~
Il existe donc une application linéaire, le tenseur des contraintes, faisant passer N

T~ = N
~ (1.1)

Dé…nition 1.1.2 Le tenseur des contraintes est une application linéaire de l’espace vec-
toriel de Rn dans lui-même. Si l’on choisit une base orthonormée ~ei cette application

8
 1.1. Tenseurs

linéaire est représentée par une matrice d’éléments ij (i; j = 1; :::n) et la relation (1.1)
donne la relation matricielle
2 3 2 32 3
T : : N1
6 1 7 6 11 12 1n
76 7
6 7 6 76 7
6 T2 7 6 21 22 : : 2n 7 6 N2 7
6 7 6 76 7
6 7 6 76 7
6 : 7=6 : : : : : 76 : 7
6 7 6 76 7
6 7 6 76 7
6:7 6 : : : : : 76 : 7
4 5 4 54 5
Tn n1 n2 : : nn Nn

Le tenseur est un tenseur réel symétrique d’ordre 2 représentable par une matrice n n
symétrique.

1.1.5 Tenseur d’élasticité

Le comportement élastique est caractérisé par une relation linéaire entre contraintes et
déformations. Cette relation s’écrit

ij = Eijkl "kl

ou inversement
"ij = ijkl kl

où Eijkl et ijkl sont les composantes de deux applications E et , inverses l’une

de l’autre, de l’espace des tenseurs symétriques dans lui-même. Ce sont les tenseurs
d’élasticité. Souvent E est appelé tenseur de rigidité et tenseur de complaisance.
Compte-tenu de la symétrie des tenseurs des contraintes et des déformations, on doit avoir
par exemple, pour E
Eijkl = Ejikl et Eijlk = Eijkl

Si on ajoute de plus l’hypothèse thermodynamique (le tenseur d’élasticité est symétrique)

alors on a
Eijkl = Eklij

9
 1.2. Quelques rappels sur les espaces de Sobolev

1.2 Quelques rappels sur les espaces de Sobolev

1.2.1 Espace L2 ( ) et (L2 ( ))n

Dé…nition 1.2.1 L’espace L2 ( ) des fonctions mesurables et de carré intégrables sur

à valeurs dans R; est un espace vectoriel muni du produit scalaire et de la norme associée:
pour tout f; g 2 L2 ( ) :
Z Z
1
(f; g)L2 ( ) = f (x)g(x)dx; kf kL2 ( ) =( jf (x)j 2 dx) 2 :

Proposition 1.2.1 ([6]) L’espace L2 ( ) muni du produit scalaire (:; :)L2 ( ) et de la norme
k:kL2 ( ) est un espace de Hilbert.

Dé…nition 1.2.2 Soit f 2 F(0; Rn ); 1une fonction dé…nie sur à valeurs vectorielles
f (x)
B 1 C
B C
B : C
B C
B C
dans Rn (n 2): On note f = B : C.
B C
B C
B : C
@ A
fn (x)
L’espace (L2 ( ))n = L2 ( ) ::: L2 ( ) est l’espace des fonction f 2 F( ; Rn ) dont les
| {z }
n fois
composantes (fi )i=1;:::n sont dans L2 ( ):

Pour faire de (L2 ( ))n un espace de Hilbert, on le munit du produit scalaire cartésien:

X
n
(f; g)(L2 ( ))n = (fi ; gi )L2 ( )
i=1

associé à la norme:
q
kf k(L2 ( ))n = kf1 k 2 L2 ( ) + ::: + kfn k 2 L2 ( ) :

10
 1.2. Quelques rappels sur les espaces de Sobolev

1.2.2 Espace de Sobolev H 1 ( )

On note D( ) l’espace des fonctions de classe C 1 à support compact dans : Notons que
les fonctions de D( ) s’annulent, ainsi que toutes leurs dérivées, sur le bord de :

Théorème 1.2.1 ([6])

L’espace D( ) est dense dans L2 ( ), c’est-à-dire pour tout f 2 L2 ( ); il existe une suite
fn 2 D( ) telle que
lim kf fn kL2 ( ) =0
n!1

Soit un ouvert de Rn ; n 2 N : Soit u 2 L2 ( ); u s’identi…e à une distribution sur ;

@u
notée encore u (car L2 ( ) L1 loc ( ) D0 ( )). On note ; j = 1; :::n les dérivées
@xj
partielles premières de u au sens des distributions.

Dé…nition 1.2.3 On appelle "espace de Sobolev" d’ordre 1 sur ; l’espace

@u
H 1( ) = u 2 L2 ( ) : 2 L2 ( ) 8j = 1; :::n :
@xj

On munit H 1 ( ) du produit scalaire, pour u; v 2 H 1 ( ):

Z X n Z
@u @v
(u; v)H 1 ( ) = u(x)v(x)dx + (x) (x)dx
j=1
@xj @xj

et de la norme correspondante
Z n Z
X
2 @u 1
kukH 1 ( ) =( ju(x)j dx + (x) 2 dx) 2
j=1
@xj
Z n Z
X 1
2
=( ju(x)j dx + jruj 2 dx) 2 :
j=1

En physique ou en mécanique, l’espace de Sobolev est souvent appelé espace d’énergie

…nie, au sens où il est constitué des fonctions d’énergie …nie (c’est-à-dire de norme
kukH 1 ( ) …nie). Les fonctions d’énergie …nie peuvent éventuellement être “singulières”, ce
qui a un sens physique possible (concentration ou explosion locale de certaines grandeurs).

Proposition 1.2.2 ([6]) L’espace H 1 ( ); muni de la norme k:kH 1 ( ) ; est un Hilbert.

11
 1.2. Quelques rappels sur les espaces de Sobolev

1.2.3 Espace de Sobolev H m ( )

On peut aisément généraliser la dé…nition de l’espace de Sobolev H 1 ( ) aux fonctions qui

sont m fois dérivables au sens faible, m 2 N:

Dé…nition 1.2.4 Pour un entier m 2 N, l’espace de Sobolev H m ( ) est dé…ni par

H m ( ) = u 2 L2 ( ):8 2 NN avec j j m, @ v 2 L2 ( )

où la dérivée partielle @ v est à prendre au sens faible et

@j j X
n
@ = , j j= j
@x1 1 ::::@xnn j=1

On munit H m ( ) du produit scalaire : pour u; v 2 H m ( )

Z X
(u; v)H m ( ) = @ u(x)@ v(x) dx
j j m

et de la norme Z X
kukH m ( ) =( @ u(x)@ u(x) dx)1=2
j j m

Proposition 1.2.3 ([6]) L’espace de Sobolev H m ( ), munit de la normek:kH m ( ) ; est un

espace de Hilbert.

Les fonctions de H m ( ) ne sont pas toujours continues ou régulières (cela dépend de m

et de la dimension n). Mais si m est su¢ samment grand, alors toute fonction de H m ( )
est continue.
Pour = Rn , il est possible de dé…nir les espaces de Sobolev en utilisant la transformée
de Fourier.

12
 1.2. Quelques rappels sur les espaces de Sobolev

1.2.4 Espace de Sobolev H s (Rn )

Rappels sur la transformée de Fourier

Dé…nition 1.2.5 Soit u 2 S(Rn )

S(Rn ) = u 2 C 1 (Rn ), 8 , 2 Nn , lim x D u(x) = 0 .

jxj!1

On dé…nit la transformée de Fourier de u 2 S(Rn ) par

Z
1
u^( ) = F (u( )) = n u(x) exp( ix: ) dx; 8 2 Rn .
(2 ) 2 Rn

Propriétés
1) 8u 2 S(Rn ) : u^ 2 S(Rn ):
2) F : S(Rn ) ! S(Rn ) est linéaire, bijective et bicontinue. Son inverse, noté F; est
dé…ni par : Z
1
Fv(x) = n v( )eix d (TFI)
(2 ) 2
Rn

3) 8u 2 S(Rn ); 8 2 Nn ; on a :

F(D u) = ij j
u^ (R1)

F(x u) = ij j D u^ (R2)

Dé…nition 1.2.6 On dé…nit la transformée de Fourier de toute distributions temperée

T 2 S 0 (Rn ) par :
hFT; 'i = hT; F'i ; 8' 2 S(Rn )

De même, la transformée de Fourier inverse de T 2 S 0 (Rn ) est donnée par :

FT; ' = T; F' ; 8' 2 S(Rn )

F FT = FFT = T; 8T 2 S 0 (Rn )

Propriétés
1) F : S 0 (Rn ) ! S 0 (Rn ) est linéaire, bijective et bicontinue.

13
 1.2. Quelques rappels sur les espaces de Sobolev

2) Les relations (R1) et (R2) sont véri…ées au sens des distributions, 8T 2 S 0 (Rn ):
3) Si u 2 S 0 (Rn )\L1 (Rn ) alors les expressions de Fu et Fu sont données explicitement
par (TF) et (TFI). De plus, est une fonction continue bornée sur Rn :

Théorème 1.2.2 (de Plancherel)

Si u 2 L2 (Rn ) alors u^ 2 L2 (Rn ): De plus, l’application F : L2 (Rn ) ! L2 (Rn ) est une

isométrie. l’identité :
kukL2 (Rn ) = k^
ukL2 (Rn )

Théorème 1.2.3 L’espace H m (Rn ) peut être dé…ni par :

m
H m (Rn ) = u 2 L2 (Rn )=(1 + j j2 ) 2 u^ 2 L2 (Rn ):

et la norme k:kH m (Rn ) est équivalente à la norme :

Z
1
u 7 ! ( (1 + j j 2 )m j^
u( )j 2 d ) 2
Rn

Dé…nition 1.2.7 Pour s 2 R, l’espace H s (Rn ) est dé…nit par :

n 0 s
o
H s (Rn ) = u 2 S (Rn ) : ((1 + j j2 ) 2 u^) 2 L2 (Rn ) .

On munit H s (Rn ) du produit scalaire :

s s
(u , v)s = ((1 + j j2 ) 2 u^ , (1 + j j2 ) 2 v^)L2 ((Rn )
Z
(u , v)s = (1 + j j2 )s u^( )^ v( ) d
Rn

et de la norme associée
Z
1
kuks = ( (1 + j j2 )s j^
u( )j2 d ) 2 .
Rn

14
 1.2. Quelques rappels sur les espaces de Sobolev

1.2.5 Espace de Sobolev H01 ( )

Dé…nition 1.2.8 Soit D( ) l’espace des fonctions de classe C 1 à support compact dans
. L’espace de Sobolev H01 ( ) est l’adhérence de D( ) dans H 1 ( ).

H01 ( ) est en fait le sous-espace de H 1 ( ) constitué des fonctions qui s’annulent sur le
bord (@ ) puisque tel est le cas des fonctions de D( ). En général, H01 ( ) est strictement
plus petit que H 1 ( ) car D( ) est un sous-espace strict de D( ): Une exception impor-
tante est le cas où = Rn : en e¤et, dans ce cas = Rn = et D(Rn ) est dense dans
H 1 (Rn ), donc, on a H01 (Rn ) = H 1 (Rn ). Cette exception se comprend aisément puisque
l’espace entier Rn n’a pas de bord.

Proposition 1.2.4 ([6]) L’espace de Sobolev H01 ( ), munit de la norme k:kH 1 ( ) ; est un
espace de Hilbert.

Théorème 1.2.4 ([6]) (Inégalité de Poincaré)

Soit un ouvert de Rn borné dans au moins une direction de l’espace. Il existe une
constante c > 0; telle que, pour toute fonction v 2 H01 ( ) :
Z Z
jv(x)j dx c jrv(x)j2 dx:
2

15
 1.2. Quelques rappels sur les espaces de Sobolev

1.2.6 Théorème de trace et formule de Green

Théorème 1.2.5 ([6])

Soit un ouvert borné de Rn , de classe C 1 . Alors D( ) est dense dans H 1 ( ) et
l’application restriction 0 dé…nie par

0 : (D( ); k:kH 1 ( ) ) ! (C 0 ( ); k:kL2 ( ) )

v7 ! 0 (v) = vj

est linéaire et continue (pour D( ) muni de la norme k:kH 1 ( ) et C 0 ( ) muni de la norme

k:kL2 ( ) ), et 0 se prolonge par continuité de H 1 ( ) dans L2 ( ) .
On dé…nit ainsi :

0 : H 1 ( ) ! L2 ( )

v7 ! 0 (v) = vj

et la continuité s’écrit :

9c > 0; 8v 2 H 1 ( ); k 0 (v)kL2 ( ) c kvkH 1 ( ) :

0 est appelée "application trace".

Théorème 1.2.6 ([6]) (Formule de Green)

Soit un ouvert borné de Rn ; de classe C 1 par morceaux. Alors si u et v sont des fonctions
de H 1 ( ); on a :
Z Z Z
@v @u
u( )dx = ( )vdx + uvNj d ; j = 1; :::; n
@xj @xj

où Nj est la j-ème composante de la normale extérieure:

16
 1.3. Problèmes aux limites elliptiques linéaires

Théorème 1.2.7 ([1]) (Formule de la divergence)

Soit un domaine de Rn , N sa normale extérieure. Soit u 2 (C 1 ( ))n un champ de
vecteur dé…ni sur : On a : Z Z
div udx = u:N d

1.2.7 Un résultat de compacité

Rappelons tout d’abord que, dans un espace de Hilbert de dimension in…nie, il n’est pas
vrai que, de toute suite bornée, on puisse extraire une sous-suite convergente au contraire
de ce qui se passe en dimension …nie.

Théorème 1.2.8 ([6]) (de Rellich)

Si est un ouvert de Rn borné et de classe C 1 alors de toute suite bornée de H 1 ( ); on
peut extraire une sous-suite convergente dans L2 ( ) (on dit que l’injection canonique de
H 1 ( ) dans L2 ( ) est compacte).

Remarque 1.2.1 Si on remplace H 1 ( ) par H01 ( ) alors non seulement le Théorème

(1.2.4) de Rellich reste vrai, mais en plus il n’est pas nécessaire de supposer que l’ouvert
est régulier.

1.3 Problèmes aux limites elliptiques linéaires

1.3.1 Formulation variationnelle

Nous appelons "formulation variationnelle" tout problème du type suivant :

8
< Trouver u 2 H tel que :
: a(u; v) = hl; vi , 8v 2 H:

où :
i) H est un espace de Hilbert (réel) dont le produit scalaire est noté (:; :)H et la norme
associée k:kH .

17
 1.3. Problèmes aux limites elliptiques linéaires

ii) a (:; :) est une forme bilinéaire sur H H; qui de plus est continue sur H H c’est
à dire :
9Ca > 0; 8u; v 2 H : ja (u; v)j Ca kukH kvkH .

iii) l est une forme linéaire sur H: De plus, l est continue sur H c’est à dire :

9Cl > 0; 8v 2 H : jl (v)j Cl kvkH :

Dé…nition 1.3.1 On dit qu’une forme bilinéaire a (:; :) : H H ! R est coercive si :

9 > 0 tel que 8v 2 H a (v; v) kvk 2 H

1.3.2 Théorème de Lax-Milgram

Théorème 1.3.1 ([1]) (Lax-Milgram)

Soit H un espace de Hilbert (réel). Soit a (:; :) une forme bilinéaire sur H H; continue
et coercive.
i) Pour tout l 2 H 0 ; il existe un unique u 2 H; solution de

a(u; v) = hl; vi ; 8v 2 H:

De plus, il existe une constante C > 0; indépendante de l telle que

kukH C klkH 0 (1.2)

ii) Si de plus a (:; :) est symétrique alors l’unique solution u de (1.2) est caractérisée par

J(u) = minJ(v)
v2H

où
1
J(v) = a(v; v) hl; vi
2

18
 2
CHAPITRE

Modélisation du
problème d’élasticité

Lorsqu’un solide est soumis à des forces extérieures, il peut se déplacer et ou se déformer.
Un solide est dit élastique si, conformément à l’intuition, il reprend spontanément sa forme
de départ lorsque les contraintes extérieures sont supprimées. Ainsi, les déformations des
solides élastiques sont par dé…nition réversibles.
Un système d’élasticité linéaire est représenté par un ensemble d’équations aux dérivées
partielles. Ces équations modélisent les déformations d’un matériau solide élastique sous
l’hypothèse de petits déplacements et de petites déformations, cette hypothèse permet
d’obtenir des équations linéaires ; d’où le nom d’élasticité linéaire. De plus, dans la théorie
de l’élasticité linéaire, les déformations et les contraintes sont supposées liées linéairement.
La linéarité est entendue au sens d’une relation de proportionnalité entre la sollicitation
et la réponse du milieu.

2.1 Equations générales

Le domaine Rn (n = 2 ou 3) qui représente un solide soumis à une force volumique
f (x). Le déplacement d’un point matériel M0 (dans une con…guration de référence) sous
!
l’action des forces extérieures f est dé…ni par u = M0 M , où M est la position actuelle

19
 2.1. Equations générales

du point M0 . Le problème consiste donc à calculer les composantes du déplacement u(x)

que subit un corps élastique soumis à diverses sollicitations.

Nous écrivons l’équilibre d’une partie quelconque D de : Cela fait intervenir le tenseur
des contraintes ; dé…ni implicitement en exprimant que la force exercée sur un élément
de surface dS, de normale unitaire N est

T dS = N dS

L’équation d’équilibre s’écrit

Z Z
(u):N dS + f dx = 0
@D D

et en utilisant la formule de la divergence, on obtient la forme ponctuelle

div( (u)) + f = 0 (1)

où div( (u)) est le vecteur de composantes

X
n
@ ij (u)
(div( (u)))j =
j=1
@xj

Nous relions maintenant le tenseur des contraintes aux tenseur des déformations.
Une telle relation est encore une loi de comportement. Nous nous restreindrons ici

20
 2.1. Equations générales

à l’élasticité linéaire. Dans ce cadre, la relation entre ces deux tenseurs est linéaire, donc
s’exprime elle-même à l’aide d’un tenseur d’ordre 4, le tenseur d’élasticité :

(u)i;j = Eijkl "kl (u) 1 i; j n

En pratique, cette loi est souvent trop générale et il est possible de la simpli…er en
considérant un matériau homogène (il a les mêmes propriétés en tout point) et isotrope
(en un point donné, il a les mêmes propriétés dans toutes les directions). La loi linéaire
de comportement élastique reliant le tenseur des contraintes au tenseur des déformations
est alors :
(u) = 2 "(u) + trace("(u)) I: (2)

Cette loi linéaire porte le nom de la loi de Hooke, où I est le tenseur identité d’ordre
2, et sont des coe¢ cients indépendants appelés les constantes de Lamé et qui sont
spéci…ques à chaque matériau. En fait, grâce aux tests mécaniques permettant de mesurer
les coe¢ cients E et propres à chaque matériau, les constantes et sont données par
les relations
E E
= et =
2(1 + ) (1 2 )(1 + )
où E est le module de Young qui correspond à la rigidité du matériau, alors que
est le coe¢ cient de Poisson et il mesure le degré de compressibilité du matériau.
Donc, dans le cas isotrope, le tenseur d’élasticité E est dé…ni par

Eijkl = (Iik Ijl + Iil Ijk ) + Iij Ilk :

Le problème consiste donc à déterminer un champ de déplacement u véri…ant l’équation

d’équilibre (1) et la loi de comportement (2): Bien entendu, ces équations doivent être
complétées par des conditions aux limites convenables sur le bord de . Celles-ci peu-
vent être de deux types :
Déplacement imposé : Sur une partie du bord 0, le déplacement est une fonction
donnée
u = u0 sur 0: (3)

Une telle condition aux limites s’appelle une condition essentielle.

21
 2.1. Equations générales

Contrainte imposée : Sur une partie du bord N, la contrainte normale est imposée

(u):N = gN sur N: (4)

Cette condition aux limites s’appelle une condition naturelle.

Finalement, déterminer le déplacement u revient à résoudre le système d’équations
aux dérivées partielles (1) et (2), avec les conditions aux limites (3) et (4).
Ainsi, le problème aux limites considéré est :
8
>
> div(2 "(u) + trace("(u)) I) = f dans
>
<
(P ) u = u0 sur 0
>
>
>
: (u):N = gN sur N

22
 3
CHAPITRE

Formulation
variationnelle du
problème de l’élasticité
linéaire

3.1 Problème de Dirichlet homogène

Soit un ouvert borné et connexe de Rn (n = 2 ou n = 3 dans les applications) de frontière
; C 1 par morceaux. On considère le problème suivant : étant donnée f 2 (L2 ( ))n ;
trouver u solution de
8
< div(2 "(u) + trace("(u)) I) = f dans
(P1 )
: u = 0 sur @
Les équations du problème (P1 ) décrivent les petits déplacements u à partir de l’état
naturel d’un solide élastique homogène et isotrope soumis à une densité volumique f dans
; les déplacements u étant …xés nuls sur le bord de :

23
 3.1. Problème de Dirichlet homogène

Munissons les espaces (L2 ( ))n et (H 1 ( ))n des normes hilbertiennes suivantes
X
n
1
kvkL2 ( ) = ( kvi k2L2 ( ) ) 2
i=1
X
n
1
kvkH 1 ( ) =( kvi k2H 1 ( ) ) 2
i=1

En vue d’établir l’existence et l’unicité de la solution de (P1 ) et a…n d’appliquer le

théorème de Lax-Milgram, il faut véri…er la coercivité de la forme bilinéaire. Pour cela,
on a besoin d’un lemme qui est un cas particulièrement simple de l’inégalité de Korn
qu’on verra un peu plus loin.

Lemme 3.1.1 Soit un ouvert de Rn : Pour toute fonction v 2 (H 1 0 ( ))n , on a

p
krvk(L2 ( ))n 2 k"(v)kL2 ( ) (1.1)

Démonstration.
Soit v 2 (D( ))n : on a
Z n Z
2 1X @vi @vj 2
j"(v)j dx = ( + ) dx
4 i;j=1 @xj @xi
n Z
1X
2 2
@vi @vj @vi @vj
= ( + +2 )dx
4 i;j=1 @xj @xi @xj @xi

puisque "(u)ij est symétrique, donc

Z n Z
1X
2
2 @vi @vi @vj
j"(v)j dx = ( + )dx
2 i;j=1 @xj @xj @xi

Intégrons par partie deux fois

Z Z Z
@vi @vj @ 2 vj @vi @vj
dx = vi dx = dx
@xj @xi @xj @xi @xi @xj

Finalement
Z n Z
1X
2
2 @vi @vi @vj
j"(v)j dx = ( + )dx
2 i;j=1 @xj @xi @xj
Z X n 2 Xn 2
1 @vi @vi
= + dx
2 i;j=1
@xj i=1
@xi
Z
1
jrvj2 dx
2

24
 3.1. Problème de Dirichlet homogène

Donc
p
2 k"(v)kL2 ( ) krvk(L2 ( ))n 8v 2 v 2 (D( ))n

Par densité de v 2 D( ) dans H01 ( ); on en déduit (1.1).

Démontrons maintenant le théorème 3.1.1.

Théorème 3.1.1 Soit un ouvert borné de Rn et = @ . Soit f 2 (L2 ( ))n . Il existe

une unique solution faible u 2 (H 1 0 ( ))n de (P1 ):

Démonstration.
1ere Etape : Recherche de la formulation variationnelle.
En utilisant le fait que
X
n
@ ij
div = ( )1 i;j n
j=1
@xj
@ui @uj
ij = ( + ) + (div(u))
@xj @xi
et trace("(u)) = div(u); on déduit les équations de (P1 ); pour 1 i n
Xn
@ @ui @uj
( + ) + (div(u)) ij = fi (1.2)
j=1
@xj @xj @xi

Pour trouver la formulation variationnelle, on multiplie chaque équation (1.2) par une
fonction test vi 2 H01 ( ) et on intègre sur ; on obtient:
X n Z Z
@ @ui @uj
( + ) + (div(u)) ij vi dx = fi vi dx
j=1
@xj @xj @xi

En appliquant la formule de Green, on trouve

n Z
X n Z
X
@ui @uj @vi @ui @uj
( + ) dx ( + ) v i Nj d
j=1
@xj @xi @xj j=1
@xj @xi

n Z
X n Z
X Z
@vi
+ (div(u)) dx (div(u))vi Nj d = fi vi dx
j=1
@xj j=1

Or, vi = 0 sur , donc

Xn Z Xn Z Z
@ui @uj @vi @vi
( + ) dx + (div(u)) dx = fi vi dx
j=1
@x j @x i @x j j=1
@xj

25
 3.1. Problème de Dirichlet homogène

On somme alors ces équations pour i allant de 1 à n, a…n de faire apparaître la

divergence de la fonction v = (v1 ; :::; vn ) et de simpli…er la première intégrale car

X n
@ui @uj @vi 1 X @ui @uj @vi
n
@vj
( + ) = ( + )( + ) = 2"(u):"(v)
i;j=1
@x j @x i @x j 2 i;j=1
@x j @x i @x j @x i

On obtient alors la formulation variationnelle :

8
>
< Trouver u 2 (H01 ( ))n tel que 8v 2 (H01 ( ))n
(F V ) R R R
>
: 2 "(u):"(v)dx + div(u) div(v)dx = f:vdx

2eme Etape : Résolution du problème varaitionnel.

On pose
Z Z
a(u; v) = 2 "(u):"(v)dx + div(u) div(v)dx
Z
L(v) = f:vdx

a(:; :) est une forme bilinéaire continue sur (H01 ( ))n (H01 ( ))n :
L est une forme linéaire continue sur (H01 ( ))n :
Montrons la coercivité en trois étapes.
Premièrement, montrons que
Z Z Z
2 2
2 j"(v)j dx + jdiv(v)j dx j"(v)j 2 dx (1.3)

avec = min(2 ; (2 + n )):

Pour cela, on utilise une inégalité algébrique: si on note
X
n
A:B = aij bij
i;j=1

le produit scalaire usuel des matrices symétriques, on peut décomposer toute matrice
réelle symétrique A sous la forme

1 1
A = Ad + Ah avec Ad = A tr A Id et Ah = tr A Id
n n

26
 3.1. Problème de Dirichlet homogène

de telle maniére que

Ad :Ah = 0 et jAj 2 = Ad 2
+ Ah 2 :

On a alors

2 jAj 2 + (tr A) 2 = 2 Ad + A h 2
+ (tr A) 2

=2 Ad 2
+ Ah 2
+ (tr A) 2

=2 Ad 2
+ (2 + n ) Ah 2

min(2 ; 2 + n ) Ad 2
+ Ah 2

jAj 2

Ce qui donne le résultat pour A = "(v):

Deuxièment, on utilise le lemme (3.1.1) qui donne une constante C > 0 telle que, pour
tout v 2 (H01 ( ))n Z Z
2
j"(v)j dx C jrvj 2 dx: (1.4)

Troisièment, on utilise l’inégalité de Poincaré (composante par composante) qui donne

une constante Ca > 0 telle que, pour tout v 2 (H01 ( ))n
Z Z
jvj dx Ca jrvj 2 dx:
2
(1.5)

Au total, ces trois inégalités (1.3),(1.4) et (1.5) conduisent à la coercivité

Z Z
2
a(v; v) = 2 j"(v)j dx + jdiv(v)j 2 dx
Z
j"(v)j 2 dx
Z
C jrvj 2 dx
Z
C Ca jvj 2 dx

Cp kvk2H 1 ( )

27
 3.1. Problème de Dirichlet homogène

Ainsi, d’après le théorème de Lax-Milgram, il existe une unique solution u 2 (H 1 ( ))n

solution faible de la formulation variationnelle (F V ):
3eme Etape : Equivalence avec le problème aux limites.
La troisième étape consiste à véri…er qu’en résolvant la formulation variationnelle
(F V ); on a bien résolu le problème aux limites (P1 ):
On a
Z Z Z
2 "(u):"(v)dx + div(u) div(v)dx = f:vdx , 8v 2 (H01 ( ))n

ou encore
X n Z X n Z n Z
X
@ui @uj @vi @vi
( + ) dx+ (div(u)) ij dx = fi vi dx , 8vi 2 H01 ( )
i;j=1
@xj @xi @xj i;j=1
@xj i=1

En appliquant la formule de Green pour vi 2 H01 ( ), on obtient :

X n Z Xn Z
@ @ui @uj
( + ) + (div(u)) ij vi dx = fi vi dx
i;j=1
@x j @x j @x i i=1

ou encore
Z Z
div(2 "(u) + div(u) I):vdx = f:vdx , 8v 2 (H01 ( ))n

Puisque div(u) = tr("(u)), on a

Z Z
div(2 "(u) + tr("(u)) I):vdx = f:vdx , 8v 2 (H01 ( ))n :

Comme D( ) est dense dans H01 ( ); alors en particulier pour ' 2 (D( ))n ; on a
Z Z
div(2 "(u) + tr("(u)) I):'dx = f:'dx

Or, D( ) est dense dans L2 ( ); donc :

Z Z
div(2 "(u) + tr("(u)) I):'dx = f:'dx 8' 2 (L2 ( ))n

Par suite
div(2 "(u) + tr("(u)) I) = f dans (L2 ( ))n

28
 3.1. Problème de Dirichlet homogène

ou encore

div(2 "(u) + tr("(u)) I) = f presque partout dans :

D’autre part : u 2 (H01 ( ))n ) u = 0 sur :

On a donc bien retrouvé l’équation et la condition aux limites (P1 ):
Ainsi la solution faible u véri…e bien (P1 ):

Proposition 3.1.1 Soit un ouvert borné de Rn , et f 2 (L2 ( ))n . L’application qui

à f 2 (L2 ( ))n fait correspondre la solution unique u 2 (H01 ( ))n de la formulation
variationnelle de (P1 ) est linéaire et continue de (L2 ( ))n dans (H01 ( ))n : En particulier,
il existe une constante C > 0 telle que, pour tout f 2 (L2 ( ))n ;

kukH 1 ( ) C kf kL2 ( ) (1.6)

Démonstration. La linéarité de cette application est évidente. La continuité est une

conséquence du théorème de Lax-Milgram. En e¤et, on a
Z
a(u; u) = f:udx Cp kuk 2 H 1 ( )

En appliquant l’inégalité de Cauchy-Schwarz, on obtient

kukH 1 ( ) C kf kL2 ( )

D’où le résultat.

Remarque 3.1.1 L’inégalité (1.6) est ce qu’on appelle une estimation d’énergie. Elle
garantit que l’énergie de la solution est contrôlée par celle de la donnée. Les estimations
d’énergie sont très naturelles d’un point de vue physique et très utiles d’un point de vue
mathématique.

L’analyse du problème aux limites (P1 ) avec une condition aux limites de Dirichlet sur
tout le bord @ est un peu trompeuse par sa simplicité. En e¤et, dès que l’on introduit
une autre condition aux limites (par exemple, de Neumann) sur une partie du bord,

29
 3.2. Problème mixte

la démonstration de la coercivité de la formulation variationnelle devient beaucoup plus

di¢ cile car on doit remplacer l’inégalité élémentaire du Lemme 3.1.1 par sa généralisation,
nettement plus technique, dite inégalité de Korn. Rappelons qu’on ne peut pas, en
général, se contenter d’une condition aux limites de Dirichlet sur l’intégralité du bord
@ car elle s’interprète comme le fait que le solide est …xé et immobile sur son bord. En
pratique, tout le bord n’est pas bloqué et souvent une partie du bord est libre de bouger, ou
bien des forces surfaciques sont appliquées sur une autre partie. Les conditions aux limites
de Neumann s’écrivent (u):N = h sur N où h est une fonction à valeurs vectorielles.
La condition de Neumann s’interprète en disant que h est une force appliquée sur le bord
(plus précisément, h est une densité de forces surfacique, homogène à une pression, tandis
que f est une densité de forces volumique). Si h = 0, aucune force ne s’applique et le
bord peut bouger sans restriction : on dit que le bord est libre.
Nous allons maintenant aborder le cas général du problème d’élasticité linéaire pour un
matériau anisotrope, et nous allons considérer le système d’élasticité avec les conditions
aux limites mixtes ( un mélange de Dirichlet et de Neumann).

3.2 Problème mixte

Soit un ouvert borné et connexe de Rn (n = 2 ou 3), de frontière ; C 1 par morceaux.
Soit 0 une partie de mesure super…cielle strictement positive, et soit N son complémen-
taire dans . Considérons le problème en déplacement u suivant :
8
>
> div( (u)) = f dans
>
<
(P2 ) u = u0 sur 0
>
>
>
: (u):N = h sur N

où la frontière du domaine est donnée par = 0[ N et le vecteur N représente la

normale unitaire extérieure sur l’interface N. Ce système d’équations est celui que satis-
fait, en petit déplacement, un corps occupant dans la con…guration initiale un domaine ,
et soumis à l’action des forces volumiques f s’exerçant dans , et des forces surfaciques h

30
 3.2. Problème mixte

s’appliquant sur la partie N de la frontière de . Le corps subit un déplacement imposé

u0 sur sa frontière 0.

Sans perte de généralité, nous traitons dans cette section, le cas homogène u0 = 0, les
autres cas s’en déduisent aisément.
On considère une relation linéaire plus générale entre le tenseur des contraints ij et
le tenseur des déformations linéaire "ij :
X
ij (u) = Eijkl "kl (u)
k;l

Munissons les espaces (L2 ( ))n et (H 1 ( ))n des normes hilbertiennes suivantes
X
n
1
kvkL2 ( ) =( kvi k2L2 ( ) ) 2
i=1
Xn
1
kvkH 1 ( ) =( kvi k2H 1 ( ) ) 2
i=1

On pose
V = v 2 (H 1 ( ))n : v = 0 sur 0 = (H 1 0 ( ))n .

V est un espace de Hilbert. En e¤et, l’application

0 : H 1 ( ) ! L2 ( 0)

est linéaire et continue. Or, V = ker( 0 ). Donc V est un sous espace fermé de
(H 1 ( ))n :D’où : V est un Hilbert.

31
 3.2. Problème mixte

Théorème 3.2.1 Soit un ouvert de Rn borné, connexe de frontière , C 1 par morceaux,

et = 0 [ N avec mes ( 0) > 0. La semi-norme jvjH 1 ( ) = krvk(L2 ( ))n dé…nit une
norme sur H 1 0 ( ) équivalente à la norme induite par k:kH 1 ( ) : On note

8v 2 H 1 0 ( ) : kvkH 1 ( ) = krvk(L2 ( ))n

0

et
(u; v)H 1 ( ) = (ru; rv)L2 ( )
0

dé…nit un produit scalaire sur H 1 0 ( ):

Démonstration. On véri…e que jvjH 1 ( ) est une norme sur H 1 0 ( ) : c’est une semi-
norme et
@v
jvjH 1 ( ) =0) = 0 pour tout i = 1; :::n
@xi
D’où v est constante (presque partout) sur qui est connexe. Puisque mes( 0) > 0; on
déduit lorsque v 2 H 1 0 ( ) que v = 0 presque partout dans :
Montrons qu’il existe deux constante C1 > 0 et C2 > 0 telle que

C1 kvkH 1 ( ) < krvk(L2 ( ))n < C2 kvkH 1 ( )

Trivialement C2 = 1: Il reste à véri…er l’existence de C1 : Supposons le contraire : il existe

une suite de fonctions (un ) de H 1 0 ( ) telle que :
1
kun kH 1 ( ) > kruk(L2 ( ))n
n
On pose :
un
vn =
kun kH 1 ( )

telle que : 8
< kvn kH 1 ( =1
)
: krvk 21
(L ( ))n <
n
Donc, on peut extraire une sous suite noté encore (vn ) qui converge faiblement vers v dans
H 1 ( ): Sachant que est borné et régulier, le théorème de Rellich indique que l’injection

32
 3.2. Problème mixte

canonique de H 1 ( ) dans L2 ( ) est compacte. Donc (vn ) converge fortement vers v dans
L2 ( ): En particulier (vn ) est une suite de Cauchy dans L2 ( ) :

kvn vm kL2 ( ) = 0 quand n; m ! 1

1
On en déduit avec krvk(L2 ( ))n < que
n

kvn vk k 2 H 1 ( ) = kvn vk k 2 L2 ( ) + kr(vn vk )k 2 ! 0+0=0

n;m!1

et (vn ) est une suite de Cauchy dans H 1 0 ( ): H 1 0 ( ) est complet pour la norme k:kH 1 ( ) (car
fermé dans H 1 ( ) qui est complet), et (vn ) converge vers un élément w 2 H 1 0 ( ) tel que
rw = 0: D’où w = 0 puisque mes( 0) 6= 0: Ce qui n’est pas compatible avec kwkH 1 ( ) = 1;
d’où l’existence de C1 > 0:
Si la fonction f 2 (L2 ( ))n et h 2 (L2 ( ))n , et si les conditions de symétrie suivantes:

Eijkl = Ejikl = Eijlk 1 i; j; k; l n

sont satisfaites alors le tenseur des contraintes est symétrique.

Pour trouver la formulation variationnelle du problème (P2 ), on multiplie la première
équation de (P2 ) par une fonction test v qui s’annule sur 0; c’est à dire v 2 V et on
intègre sur ; on obtient
Z Z
div( (u)):vdx = f:vdx

XZ @ ij (u)
Z
vi dx = f:vdx (2.1)
i;j
@xj

On suppose u 2 (H 2 ( ))n alors 8 v 2 (H 1 ( ))n , on a d’après la formule de Green

XZ @vi X Z @ ij (u) XZ
( ij (u) dx = ( )vi dx + ( ij (u)Nj vj )) d
i;j
@xj i;j
@xj i;j

Ainsi, l’équation (2.1) devient :

XZ @vi
Z XZ
( ij (u) dx = f:vdx + ( ij (u)Nj vj )) d ; 8v 2 V
i;j
@x j i;j

33
 3.2. Problème mixte

Or, v = 0 sur 0 ce qui implique

XZ @vi
Z XZ
( ij (u) ) dx = f:vdx + ( ij (u)Nj vj )) d ; 8v 2 V
i;j
@xj i;j N

En tenant compte de la condition au limites (u):N = h sur N , on obtient

XZ @vi
Z Z
( ij (u) dx = f:vdx + h:v d 8v 2 V
i;j
@xj N

Grâce à la symétrie de , ona

XZ @vi XZ
( ij (u) )dx = ( ij (u)"ij (v) dx
i;j
@x j i;j

Ainsi, la formulation variationnelle est

8
>
< Trouver u 2 V; tel que 8v 2 V
(F V ) P R R R
>
: ( ij (u)"ij (v) dx = f:vdx + N h:v d
i;j

On pose :

XZ
a(u; v) = ( ij (u)"ij (v)) dx
i;j

et Z Z
L(v) = f:v dx + h:v d
N

Nous allons maintenant énoncer et démontrer un résultat d’existence et d’unicité de la

solution du problème variationnel (FV) de l’élasticité linéaire, pour cela on doit démontrer
un résultat très important qui est "l’inégalité de Korn".
On démontre l’inégalité de Korn en s’appuyant sur le lemme de Lions suivant.

Lemme 3.2.1 (Lions)

Soit un ouvert régulier de Rn et v un élément de H 1
( ) tel que

@v 1
2H ( ); 8i = 1; :::; n:
@xi

alors v 2 L2 ( ):

34
 3.2. Problème mixte

Démonstration. On introduit l’espace

1 @v 1
X( ) = v2H ( ): 2H ( ); 8i
@xi

qui est un Hilbert pour la norme

Xn
@v
2
1
(kvk2H 1( )+ )2 :
i=1
@xi H 1( )

On montre en plusieurs étapes que X( ) = L2 ( ).

1) Dans le cas où = Rn ; le lemme de Lions se déduit rapidement d’un calcul de trans-
1
formé de Fourier. En e¤et une distribution tempérée v est dans H ( ) si et seulement
1
si (1 + j j 2 ) 2 v^ 2 L2 ( ):
Par conséquent,
Z Z
1 @v 1 v ( )j 2
j^ i j^2
v ( )j 2
v2H ( ); 2H ( ), d < +1 et d < +1
@xi 1 + j j2 1 + j j2

En sommant toutes les intégrales, on obtient

Z
v ( )j 2 d < +1
j^

c’est-à-dire v^ 2 L2 ( ): D’où v 2 L2 ( ) (d’après le théorème de Plancherel).

2) Il su¢ t de montrer que (2.2) pour le demi-espace
n 0 0
o
= (x ; xn ) : x 2 Rn 1 ; xn > 0 :

En e¤et, soient 0; 1 ; :::; N tels que

X
N

0 2 D( ) ; 1 2 D( ) ; i = 1; :::; N; i =1
i=1

et supp( i ) carte locale dé…nissant : Notons que

' 2 D( ); v ! 'v l’application de X( ) à lui-même.

On peut écrire aussi

X
N
v= iv
i=1

35
 3.2. Problème mixte

On peut considérer 0v comme un élément de X(Rn ); en utilisant la première étape

0v 2 L2 ( ): On aura ce résultat si on montre que

iv 2 L2 ( )

Si on suppose que est une variété continuellement di¤érentiable de dimension n 1,

0 0
puis l’image de iv est dans l’espace X( ) avec = (x ; xn ) : x 2 Rn 1 ; xn > 0 ; ce
qui achève la preuve de cette étape.
3) On introduit
d'
H01 (0; 1; L2 (Rn 1 )) = 'j'; 2 L2 (0; 1; L2 (Rn 1 )); '(x0 ; 0) = 0
dxn
où x0 = (x1 ; x2 ; :::; xn 1 ):

1
H (0; 1; L2 (Rn 1 )) dual de H01 (0; 1; L2 (Rn 1 ))

où L2 (Rn 1 ) est identi…é à son dual,

dv
Y( )= vjv; 2 H 1 (0; 1; L2 (Rn 1 ))
dxn
On montre que Y ( ) est dense dans X( ). En e¤et, soit m une suite régularisante dans
D(Rn 1 ). Pour v 2 X( ); on dé…nit

vm = v(x0 ) m;

ou plus précisément Z
vm = v(x0 y 0 ; xn ) m (y
0
)dy 0 :
R2
Alors comme m ! 1; on a vm ! v dans X( ). Notez que vm est en particulier contenu
dans Y ( ).
4) D( ) est dense dans X( ). Il su¢ t de montrer la densité dans Y ( ): Comme une
simple conséquence du théorème de Hahn-Banach, il su¢ t de montrer que pour toute
fonctionnelle M qui est nule sur D( ); est nule sur un espace entier Y ( ). Soit v 7 ! M (v)
une forme linéaire continue dans Y ( ) et donc de la forme
Z 1
dv
M (v) = (f; v) + g; dxn
0 dxn
f; g 2 H01 (0; 1; L2 (Rn 1 ))

36
 3.2. Problème mixte

On suppose que M = 0 sur D( ): Si f~ ; g~ indique la continuité de f et g en 0 pour xn < 0;

la supposition est équivalente à
d~
g
f~ =0
dxn
d~
g
Donc 2 H01 (0; 1; L2 (Rn 1 )), mais alors
dxn
Z 1 Z 1
dv dg
g; dxn = ; v dxn 8v 2 Y ( )
0 dxn 0 dxn

Donc M = 0 sur Y ( ):
5) Pour v 2 D( ) on pose
8
< v(x) si xn > 0
P v(x) =
: a v(x0 ; x ) + a v(x0 ; 2x )
1 n 2 n

a2
où a1 + a2 = 1; a1 + = 1: On véri…e que v 7 ! P v est continue dans D( ); avec
2
la topologie induite par X( ) dans X(Rn ): En utilisant cela, et du point précédent,
nous pouvons étendre P à une application linéaire continu de X( ) ! X(Rn ) et tel
que P restreint à équivaut à v. Ensuite, pour v 2 X( ); P v 2 X(Rn ), donc comme
conséquence de l’étape 1), P v 2 L2 (Rn ) donc

v 2 L2 ( )

Il reste à montrer la continuité de P:

8
> @u
@ < xn > 0
Pv = @xn
@xn > @v 0 @v 0
: a1 (x ; xn ) 2a2 (x ; 2xn ) xn < 0
@xn @xn
@v
On pose = w, on introduit
@xn
8
< w(x) si xn > 0
Qw =
: a w(x0 ; x ) 2a w(x0 ; 2x ) x < 0
1 n 2 n n

On a besoin de montrer que P (resp Q) est continu dans D( ) avec la topologie induite
1 1
par H ( ) dans H (Rn ); donc par transposition, t P (resp t Q) est continu en tant que

37
 3.2. Problème mixte

fonctionnelle H 1 (Rn ) dans H01 ( ). Maintenant, pour H 1 (Rn ), nous avons

t 1 xn
P '(x) = '(x) + a1 '(x0 ; xn ) + a2 '(x0 ; );
2 2
t xn
Q'(x) = '(x) a1 '(x0 ; xn ) a2 '(x0 ; )
2

Alors
t
P '(x0 ; 0) =t Q'(x0 ; 0) = 0

d’où le résultat.
On adopte dans l’énoncé et la démonstration du théorème de l’inégalité de Korn la
convention d’Einstein de sommation des indices répétés.

Théorème 3.2.2 (L’inégalité de Korn)

Soit un ouvert borné et régulier. Il existe une constante C>0 (dépendant de ) tel que:
Z Z
"ij (v) "ij (v)dx + vi vi dx C kvk2H , 8v 2 (H 1 ( ))n

Notons que l’inégalité de Korn est loin d’être une trivialité car le premier membre
de l’inégalité ne fait intervenir que certaines combinaisons linéaires des dérivées partielles
@vi @v1 @v2 1 @v1 @v2
premières (en dimension n = 2; , et ( + )) tandis que le second
@xj @x1 @x2 2 @x1 @x2
@vi
membre de l’inégalité fait intervenir toutes les dérivées : On s’attend donc plutôt à ce
@xj
que "(v) contienne moins d’informations que ru, et il est étonnant de voir que l’on peut
néanmoins contrôler dens L2 dans toutes les dérivées partielles de v à l’aide de ses seules
combinaisons linéaires, et de v elle-même. La raison profonde est que, contrairement aux
apparences, "(v) contient bien toute l’information de rv, puisque grâce à l’identité (2:3),
on peut récupérer toutes les dérivées secondes de v à partir de "(v):
Démonstration.
Soit E l’espace des fonction v 2 (L2 ( ))n telles que:

"ij (v) 2 (L2 ( ))n ; 8i; j:

E est un espace de Hilbert pour la norme

Z Z
( "ij (v) "ij (v)dx + vi vi dx)1=2 :

38
 3.2. Problème mixte

En utilisant la dé…nition de "ij nous constatons que

@ 2 vi @ @ @
= "ik (v) + "ij (v) "jk (v) (2.3)
@xj @xk @xj @xk @xi

@
Si v 2 E alors "ij 2 L2 ( ): Par conséquent, "ij (v) 2 H 1
( ): Donc, en utilisant
@xk
l’égalité établie, on a
@ 2 vi 1
2H ( ) 8i; j; k
@xj @xk
@vi
En appliquant le Lemme 3.2.1 à ; on voit que
@xk
@vi
2 L2 ( ); 8i; k:
@xk

Ainsi
v 2 (H 1 ( ))n :

Nous avons l’égalié algébrique E = (H 1 ( ))n et puisque l’injection

(H 1 ( ))n E

est continue, on applique le théorème du graphe fermé pour obtenir un isomorphisme des
espaces de Banach.

Remarque 3.2.1 L’interprétation mécanique de l’inégalité de Korn est la suivante :

l’énergie élastique, proportionnelle à la norme du tenseur des déformations "(u) dans
L2 ( ) contrôle la norme du déplacement u dans (H 1 ( ))n , à l’addition près de la norme
de u dans L2 ( ).

Pour montrer que l’inégalite de Korn entraine la coercivité du problème variationnel

de l’élasticité linéarisée, nous allons avoir besoin du lemme suivant.

Lemme 3.2.2 Soit v 2 D0 ( ) une distribution telle que "(v) = 0: Alors v est un déplace-
ment rigide in…nitésimal

v(x) = a ^ x + b avec a; b 2 Rn :

39
 3.2. Problème mixte

@v
Démonstration. En e¤et, si "(v) = 0; on déduit de l’identité (2.3) que =0
@xj @xk
pour tout j; k.

Par conséquent, v est une fonction a¢ ne de la forme

v(x) = Ax + b où A 2 Mn et b 2 Rn :

Comme rv = A, on déduit que : AT + A = 0: Donc A est antisymétrique.

On peut par conséquent représenter son action sur un vecteur x comme un produit
vectoriel par un vecteur a:
P
La loi de comportement ij = Eijkl "kl (u) correspond à un matériau anisotrope.
k;l
Lorsqu’en outre ce matériau est non homogène, les coe¢ cients d’élasticité Eijkl dépendent
de x ; on suppose alors que les fonctions Eijkl appartiennent à L1 ( ); 1 i; j; k; l n
et que la propriété d’ellipticité est uniforme, c’est-à-dire qu’il existe une constante >0
indépendante de x: On a alors

Théorème 3.2.3 (Existence et unicité de la solution)

Si le tenseur d’élasticité E 2 (L1 ( ))n , véri…e la propriété de symétrie

Eijkl (x) = Ejikl (x) = Eijlk (x); 1 i; j; k; l n; 8x 2

et si la condition d’ellipticité est satisfaite

X X
9 0 tel que Eijkl "ij "kl ( "kl )2
i;j;k;l k;l

pour tout tenseur symétrique ", alors le problème variationnel (FV)

8
< Trouver u 2 V tel que
(F V )
: a(u; v) = L(v); 8v 2 V:

admet une solution unique dans V .

Si de plus, on a
Eijkl = Eklij 1 i; j; k; l n

40
 3.2. Problème mixte

alors la solution de (F V ) correspond à l’unique solution du problème de minimisation

min J(v)
v2V

où
1
J(v) = a(v; v) L(v)
2

Démonstration.
a(:; :) est une forme bilinéaire sur V V et L est une forme linéaire sur V:
1) Montrons la coercivité

XZ
a(u; u) = ( ij (u)"ij (u)) dx
i;j 1

XZ X
= ( Eijkl "kl (u)) "ij (u) dx
i;j k;l
Z X
= ( Eijkl "kl (u)) "ij (u) dx
i;j;k;l
Z X
( "kl (u))2 dx
k;l
X
k"ij k2L2 ( )
i;j

X
Pour simpli…er la notation, on met : k"(u)k = k"ij k2L2 ( ) : Donc :
i;j

a(u; u) k"(u)k

Pour montrer que a(:; :) est coercive, il su¢ t de trouver Cp > 0 tel que

k"(u)k Cp kuk2V

On raisonne par l’absurde. On suppose donc qu’il existe une suite (un ) de V telle que
8
< k"(u )k ! 0 quand n ! 1
n
: et ku k 1 = 1 n H ( )

41
 3.2. Problème mixte

Comme la suite (un ) est bornée dans (H 1 ( ))n ; on peut extraire une sous suite notée
(un ) encore telle que
un ! u faiblement dans (H 1 ( ))n :

La trace étant continue, on a aussi u = 0 sur 0, c’est-à-dire u 2 V: De plus :

k"(un )k ! "(u) faiblement dans (L2 ( ))n :

Par ailleurs, la première partie de notre hypothèse est que

k"(un )k ! 0 fortement dans (L2 ( ))n :

On déduit donc que "(u) = 0:

Par le lemme 3.2.2, il vient que v est un déplacement rigide in…nitésimal. Or, un
déplacement in…nitésimal ne peut s’annuler que
a) sur une droite si a 6= 0
b) nulle part si a = 0 et b 6= 0
c) partout si a = b = 0:
Comme il s’annule sur une partie du bord qui est de mesure positive, on est néces-
sairement dans la situation c), c’est à dire que u = 0:
On utilise le théorème de Rellich, pour déduire que un ! 0 fortement dans (L2 ( ))n :
Par l’inégalité de Korn, il vient que

run ! 0 fortement dans (L2 ( ))n ;

d’où …nalement un ! 0 fortement dans (H 1 ( ))n : Ceci contredit la seconde partie de

notre hypothèse.
La forme bilinéaire a(:; :) est donc coercive.

42
 3.2. Problème mixte

2) Montrons la continuité de a(:; :):

XZ X @vi
ja(u; v)j = ( Eijkl "kl (u) dx
i;j k;l
@xj
XZ @vi
jEijkl "kl (u)j dx
i;j;k;l
@xj
XZ @vi
kEk1 "kl (u) dx
i;j;k;l
@xj

C( ) kukV kvkV

La forme bilinéaire a(:; :) est donc continue sur V .

3) Montrons la continuité de L:
Z Z
jL(v)j = f:vdx + h:vd

Z Z
f:vdx + h:vd

En utilisant l’inégalité de Cauchy-Schwarz

jL(v)j kf kL2 ( ) kvkL2 ( ) + khkL2 ( ) kvkL2 ( ) :

Or, d’après le théorème 1.2.2 l’application trace 0 : H 1 ( ) ! L2 ( ) est continue et

kvkL2 ( ) c kvkH 1 ( )

On a donc

jL(v)j kf kL2 ( ) kvkH 1 ( ) + c khkL2 ( ) kvkH 1 ( )

(kf kL2 ( ) + c khkL2 ( ) ) kvkH 1 ( )

donc L est une forme linéaire continue sur V .

On peut aisément conclure en utilisant le théorème de Lax-Milgram qu’il existe une unique
solution faible u 2 V . De plus, si le tenseur E est tel que

Eijkl (x) = Eklij (x); 8 x 2

43
 3.2. Problème mixte

alors la forme bilinéaire a(:; :) est symétrique. En e¤et

XZ
a(u; v) = ( ij (u)"ij (v)) dx
i;j 1

XZ X
= ( Eijkl "kl (u)) "ij (v) dx
i;j k;l
XZ
= (Eijkl "kl (u)) "ij (v) dx
i;j
XZ X
= ( Eklij "ij (v)) "kl (u) dx
k;l i;j
XZ
= ( ij (v))"kl (u) dx
k;l

= a(v; u)

Dans ce cas, l’équation variationnelle a(v; u) = L(v); 8 v 2 V correspond à la condition

d’optimalité de minimisation
1
min ( a(v; v) L(v)). (2.4)
v2V 2
Finalement, on montre que la solution unique de (FV) est bien une solution du problème
aux limites (P2 ):
En e¤et, grâce à la symétrie de , on a
XZ XZ @vi
a(u; v) = ( ij (u)"ij (v)) dx = ij (u) dx
i;j i;j
@xj

Soit ' 2 (D( ))n : Alors,

XZ @'i
Z Z
a(u; ') = ij (u) dx = f:' dx + h:' d
i;j
@xj
N

et comme la fonction test ' est nulle sur N , alors

XZ @'i
Z
ij (u) dx = f:' dx (2.5)
i;j
@xj

En appliquant la formule de Green, on a

XZ @'i X Z @ ij (u) Z
ij (u) dx = 'i dx = div( (u)):' dx
i;j
@xj i;j
@xj

44
 3.2. Problème mixte

Ainsi, (2.5) devient Z Z

div( (u)):' dx = f:' dx

ce qui implique
div( (u)) = f dans (L2 ( ))n :

Comme f 2 (L2 ( ))n alors

div( (u)) = f presque partout dans :

Soit v 2 (H 1 0 ( ))n , d’après la formule de Green, on a

X Z @ ij (u) XZ @vi XZ
vi dx = ij (u) dx ij (u)Ni vj d :
i;j
@xj i;j
@xj i;j

Or
X@ ij (u)
dx = div( (u))
i;j
@xj

Donc Z XZ @vi XZ
div( (u)):vdx = ij (u) dx ij (u)Ni vj d :
i;j
@xj i;j

Par suite
Z XZ Z Z
a(u; v) = div( (u)):vdx + ij (u)Ni vj d = f:v dx + h:v d :
i;j
N

Comme div( (u)) = f presque partout dans ; alors

XZ Z
ij (u)Ni vj d = h:v d :
i;j
XZ XZ ZN
ij (u)Ni vj d + ij (u)Ni vj d = h:v d :
i;j i;j
N 0
Z Z ZN
(u):N:vd + (u):N:vd = h:v d :
N 0 N

Comme v = 0 sur 0, on a
Z Z
(u):N:vd = h:v d :
N N

45
 3.2. Problème mixte

D’où :
h = (u):N preque partout sur N

En…n, comme u 2 V = (H 1 0 ( ))n , et que est un ouvert régulier, la trace de u est

bien dé…nie et u = 0 presque partout sur 0: Ainsi, la solution faible u 2 V est bien une
solution du problème (P2 ):

Remarque 3.2.2 L’interprétation de l’égalité (2.4) est la suivante. La solution de la

formulation variationnelle (F V ) réalise le minimum d’une énergie mécanique. En termes
mécaniques, J(v) est la somme de l’énergie de déformation
Z
1 X
( ( ij (u)"ij (u)) dx)
2 i;j 1

et de l’énergie potentielle des forces extérieures (ou travail des forces extérieures au
signe près) Z Z
f:udx h:ud :
N

Proposition 3.2.1 L’application de (L2 ( ))n (L2 ( ))n dans (H 1 ( ))n qui à f fait
correspondre u, unique solution faible de problème aux limites (P2 ) est linéaire continue.

Démonstration.
Il est clair que (f; g) ! (L2 ( ))n (L2 ( ))n est linéaire. Pour montrer qu’elle est continue
de (L2 ( ))n (L2 ( ))n dans (H 1 ( ))n ; on prend v = u dans la formulation variationnelle
(F V ) et on utilise la coercivité de la forme bilinéaire, on obtient :
XZ Z Z
a(u; u) = ( ij (u)"ij (u) dx = f:udx + h:ud Cp kuk2H 1 ( )
i;j N

En utilisant l’inégalité de Cauchy-Schwarz

1
kuk2H 1 ( ) (kf kL2 ( ) kukL2 ( ) + khkL2 ( ) kukL2 ( ) ):
Cp

Grâce au théorème de trace, on a

kuk2H 1 ( ) C(kf kL2 ( ) + khkL2 ( ) ) kukH 1 ( ) :

46
 3.2. Problème mixte

Donc
kuk2H 1 ( ) C(kf kL2 ( ) + khkL2 ( ) ) kukH 1 ( ) :

Ce qui prouve la continuité.

Remarque 3.2.3 La formulation variationnelle n’est rien d’autre que le principe des
travaux virtuels en mécanique. Poursuivant cette analogie, l’espace V est l’espace des
déplacements v cinématiquement admissibles, et l’espace des tenseurs symétriques
2 (L2 ( ))n , tels que div = f dans et :N = h sur N , est celui des tenseurs de
contraintes statiquement admissibles.

47
 Conclusion

Dans ce mémoire, nous avons montré que le problème d’élasticité linéaire est bien posé,
c’est-à-dire, nous avons montré l’existence et l’unicité de la solution ainsi que la stabilité
de celle-ci par rapport à des perturbations des données, en appliquant l’approche vari-
ationnelle. En premier lieu, nous avons établi la formulation variationnelle, puis nous
avons résolu la formulation variationnelle en utilisant le théorème de Lax-Milgram, a…n
d’appliquer ce théorème, la seule hypothèse non triviale à véri…er porte sur la coercivité
de la forme bilinéaire a(:; :), nous avons réglé le problème en démontrant l’inégalité de
Korn. En…n, nous avons interprété la formulation variationnelle pour véri…er qu’on a bien
résolu le problème aux limites de l’élasticité.
L’approche variationnelle s’avère riche et puissante, elle permet de construire la méth-
ode des éléments …nis, fournissant ainsi un moyen d’approcher la solution qui, bien sou-
vent, n’est pas calculable explicitement.
L’élasticité linéaire ne représente qu’un comportement simpli…é très particulier des
solides. Dans la nature, la plupart des problèmes rencontrés obéissent à des lois de
comportement non linéaires, c’est-à-dire le tenseur des contrainte est une fonction
non linéaire du tenseur des déformation "(u) dans le cadre de l’élasticité non linéaire
en petits déplacements. Il en résulte que les équations correspondantes de l’équilibre
sont non linéaires, contrairement aux équations rencontrées dans l’élasticité linéaire. Les
phénomènes mécaniques non linéaires sont un domaine très actif actuellement de la mé-
canique des solides, en liaison avec la recherche des matériaux nouveaux (polymères,
matériaux composites, . . . ) et de l’étude de leurs propriétés mécaniques.

48
 Bibliographie

[1] ALLAIRE G., Analyse numérique et optimisation, Editions de l’Ecole polytechnique,

Palaiseau (2005).

[2] BONY J.-M., Cours d’analyse. Théorie des distributions et analyse de Fourier, Edi-
tions de l’Ecole Polytechnique, Palaiseau (2001).

[3] Brézis H., Analyse fonctionnelle, Masson, Paris (1983).

[4] Lions J. L. , Problèmes aux limites dans les équations aux dérivées partielles, Presses
de l’Univ. de Montréal (1965).

[5] RAVIART P.-A.,THOMAS J.-M., Introduction à l’analyse numérique des équations

aux dérivées partielles, Masson, Paris (1983).

[6] Cours des espaces de Sobolev et de formulation variationnelle du master 1 et du master

2.

49