Formulation avec conditions de Dirichlet et de Neumann

On considère ΩR3 tel que Ω=ΩDΩN. On se donne g:ΩNR3 la charge normale(densité de force) sur ΩN. Le problème mélangeant condition de Dirichlet, imposant le déplacement, et de Neumann, imposant les contraintes normales, en élasticité linéaire s’exprime sous la forme suivante

Problème mixte
σ(d)+f=0 in Ωσ(d)=λ(d)I+μ(d+dT) in Ωd=0 on ΩDσ(d)\normal=g sur ΩN
si g=0 sur ΩN on dira que la surface est laissée libre (de force ou de charge).

Le choix naturel d’espace fonctionel est X[H10,ΩD(Ω)]3 pour la solution et les fonctions tests avec

X[H10,ΩD(Ω)]3={v[H1(Ω)]3 s.t. v=0  on ΩD}

La formulation variationnelle est obtenue en faisant le produit scalaire de la première équation par la fonction test v. Observez de plus que:

Intégration par partie (Théorème de la divergence)
Ωσ(d)  v=Ωσ(d) : vΩv(σ(d)\normal)

et

Symmetry of σ(d)
σ(d):v=σ(d):ε(v)

et nous obtenons

Ωσ(d) : ε(v)Ωv(σ(d)\normal)=Ωfv

La formulation faible s’écrit alors

Formulation faible du problème mixte

On cherche dH10,ΩD(Ω)]3 telle que

a(d,v)=Ωfv+ΩNg\normalv[H10,ΩD(Ω)]3

Avec

a(d,v)=Ωσ(d) : ε(v)=Ωλdv+Ω2με(d):ε(v)
En mécanique du continuum, la fonction de test v joue le rôle d’un déplacement virtuel et la formulation faible précédente exprime le principe du travail virtuel.

1. Problème bien posé

Le problème faible avec condition aux limites mixtes est bien posé grâce à l’inégalité de Korn.

2. Finite element approximation

On considère une approximation du problème mixte par élément fini H1-conforme à partir d’une famille de maillages affines, géométriquement conformes {Th}h>0 et d’un élément fini de Lagrange de degré k1 décrit par {ˆK,ˆP,ˆΣ}. On introduit Ensuite

Xkh={vh[C0(Ω)]3|KTh,vhφgeoK[ˆP]3; vh=0 on ΩN}
Problème discret

Trouver dXkh telle que

a(dh,vh)=Ωfvh+ΩNg\normalvhXkh
Proposition: Estimation a priori pour le problème mixte

Soit d une solution du problème faible, on a tout d’abord que lim. De plus si \disp{d} \in [H^{l+1}(\Omega)]^3 \cap X^k_h pour l\in \{1\ldots k\} alors il existe c telle que

\forall h, \quad \|\disp{d}-\disp{d}_h\|_{1,\Omega} \leq c h^l |\disp{d}|_{l+1,\Omega}.