Problème de l'obstacle - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Introduction

Le problème de l'obstacle est un problème classique de mécanique. Pour visualiser ce problème, il faut imaginer une membrane recouvrant un objet appelé alors obstacle (tel un film cellophane recouvrant un rôti de boeuf !). En effet, ce problème consiste à trouver une courbe solution u qui a une position très précise par rapport à l'obstacle et qui en plus vérifie une propriété de minimisation de longueur.

Pour mieux appréhender le problème, considérons-le en dimension 1. La membrane est alors un élastique recouvrant un objet. Cet élastique qui se trouve toujours au-dessus de l'objet tend à minimiser sa longueur. De plus, pour de petites variations, minimiser sa longueur revient à minimiser son énergie, en effet paramétrant l'élastique sur l'intervalle [ − 1,1] par

\begin{cases}           x(t) =  t \\            y(t) =  u(t) \end{cases}

pour t \in [-1,1] .

la longueur de l'élastique entre − 1 et 1 est :

.

Nous voudrions minimiser cette longueur, c'est-à-dire trouver

\min_{u \in K} \int_{-1}^{1} \sqrt{1 + u'(t)^2} \mathrm dt .

Or, si u' est suffisamment petit, nous avons \sqrt{1 + u'(t)^2} \sim 1 + \tfrac{1}{2}u'(t)^2 .

Minimiser la longueur de l'élastique revient donc à trouver

{\min_{u \in K} { \int_{-1}^{1} \left(   1 + {\frac{1}{2}} u'(t)^2  \right) \mathrm dt}} ,

ce qui revient à minimiser

{\frac{1}{2}} { \int_{-1}^{1} u'(t)^2 \mathrm dt} = {\frac{1}{2}} a(u,u) = {\frac{1}{2}} E(u) .

Approche théorique

Existence et unicité

Le problème de l'obstacle peut être vu comme une application du théorème de Stampacchia, en considérant la proposition suivante:

Proposition —  L'espace de Hilbert H_0^1{(-1,1)}, le convexe fermé non vide

K = \{ u \in H_0^1(-1,1) \text{ tel que } u \geq \chi \text{ presque partout } \}

\chi \in C^1 et χ < 0 aux bords, la forme bilinéaire, continue, coercive et symétrique

{{refnec| \begin{array}{cccc}              a : & H \times H & \to & \mathbb{R} \\                  & (u,v) & \mapsto & \displaystyle\int_{-1}^1 u'(t)v'(t) \mathrm dt                 \end{array} }}

et la forme linéaire continue

{{refnec|}}

vérifient les hypothèses du théorème de Stampacchia.

(Les démonstrations se font aisément en utilisant les inégalités de Hölder et de Poincaré.)

Le théorème de Stampacchia s'applique donc, et ainsi il existe un unique

u \in K = \{ u \in H_0^1(-1,1) \text{ tel que } u \geq \chi \text{ presque partout } \}

tel que :

a(u,v-u)  \  \geq  \  L(v-u) \quad \forall  \  v \in K.

Ce qui est équivalent à dire qu'il existe un unique u \in K tel que :

  {\int_{-1}^1 {u'(t)(v-u)'(t) \mathrm dt}} \  \geq  \  {\int_{-1}^1 {g(t)(v-u)(t) \mathrm dt}} \quad \forall  \  v \in K.

De plus, a étant symétrique, alors u est caractérisé par la propriété :

\begin{cases}      u \in K \\      \displaystyle{\frac{1}{2}} {\int_{-1}^1 {(u'(t))^2 \mathrm dt}} - {\int_{-1}^1 {g(t)u(t) \mathrm dt}} = {\min_{v \in K}\left(\frac{1}{2}{\int_{-1}^1 {(v'(t))^2 \mathrm dt}} - {\int_{-1}^1 {g(t)v(t) \mathrm dt}}\right)}   \end{cases}

Propriétés de la solution

Démontrons à présent quelques propriétés vérifées par la solution u, en utilisant la dérivée seconde de u. Cependant u\in H^1_0 , u' est alors une dérivée faible et se pose alors le problème de définir la dérivée seconde. C'est pourquoi il faut utiliser une théorie plus générale : la théorie des distributions.

Proposition —  Avec les notations du théorème de Stampacchia,

  • - u''-g \geq 0 au sens des distributions sur Ω = ] − 1,1[
  • u'' − g = 0 au sens des distributions sur  \omega=\{t \mid u(t) > \chi (t)\}.

Comparaison d'énergies

Proposition —  Soit T0 la tangente à χ au point d'abscisse t_0\in [-1,0] telle que T0( − 1) = 0. Soit T1 la tangente à χ au point d'abscisse t_1\in [0,1] telle que T1(1) = 0. Définissons deux fonctions dans K :

 u(t) = \left\{ \begin{array}{ll} T_0(t) \  & t\in [-1,t_0],\\  \chi (t) \  & t\in [t_0,t_1],\\ T_1(t) \  & t\in [t_1,1]. \end{array} \right.

 u_1(t) = \left\{ \begin{array}{ll} T_0(t) \  & t\in [-1,0],\\ T_1(t) \  & t\in [0,1]. \end{array} \right.

  • L'énergie de u, \textstyle E(u)= { \int_{-1}^1 u'(t)^2 \mathrm dt} est inférieure à celle de u1, \textstyle E(u_1)= { \int_{-1}^1 u_1'(t)^2 \mathrm dt} .
  • En considérant u1 telle qu'elle a été définie précédemment et u2 un triangle centré (voir figure 2), l'énergie de u2 est supérieure à celle de u1.
  • L'énergie d'un triangle décentré u3 est supérieure à celle d'un triangle centré u2 (voir figure 2).
  • Soit u_*~ la solution tordue (voir figure 2). Alors l'énergie de u_*~ est supérieure à celle de u.

Conclusion

Fig. 3: Solution dans le cas d'un obstacle concave

La solution semble donc être (voir Fig. 3)

 u(t) = \begin{cases} T_0(t) \  & t\in [-1,t_0],\\  \chi (t) \  & t\in [t_0,t_1],\\ T_1(t) \  & t\in [t_1,1]. \end{cases} ,

et on vérifie facilement que c'est le cas car cette solution vérifie le théorème de Stampacchia qui garantit son unicité :

 { \int_{-1}^{1} u'(v-u)'}  =  4{\int_{t_0}^{t_1} (v - \chi)} \geq 0.
Page générée en 0.372 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales | Partenaire: HD-Numérique
Version anglaise | Version allemande | Version espagnole | Version portugaise