- Erreur en relation de comportement
-
L'erreur en relation de comportement est une notion utilisée en mécanique des milieux continus pour quantifier l'erreur commise lors de la résolution approchée d'un problème régi par des équations d'équilibre, des conditions aux limites et une relation de comportement lorsque les champs approchés sont admissibles.
Lorsque l'on cherche la solution à un problème du type trouver le champ de déplacement et de contrainte d'un solide soumis à un chargement donné et avec certains déplacements imposés, la méthode naturelle consiste à partir des équations locales d'équilibre et à les intégrer. Cependant selon le type de domaine et/ou de chargement cela peut s'avérer impossible.
On est alors amené à chercher des solutions approchées. On se place dans l'hypothèse des petites perturbations, dans le cadre de l'élasticité linéaire.
Soit Ω l'espace occupé par le solide.
De plus on suppose que le problème est bien posé, c’est-à-dire qu'il existe une partition à deux éléments
de 
, le bord de Ω telle que sur 
il n'y ait que des déplacements imposés, et sur 
que des efforts imposés.Sommaire
Champs de déplacement cinématiquement admissibles
Champs de déplacement cinématiquement admissibles
Le champ de déplacement solution vérifie certaines conditions aux limites et possède une certaine régularité issue de la structure du milieu. Ces conditions définissent l'espace des champs cinématiquement admissibles, qui est l'ensemble des champs de déplacement suffisamment réguliers satisfaisant aux conditions aux limites.
sur où 
est le déplacement imposé.On note cet espace CA(Ω)
Dans le cas général, cet espace est un espace affine de dimension infinie
Champs de déplacement cinématiquement admissibles à zéro
On définit l'espace des champs cinématiquement admissibles à zéro comme étant l'ensemble des champs de déplacement suffisamment réguliers sur Ω satisfaisant à des conditions aux limites nulles
Dans le cas général, cet espace est un espace vectoriel de dimension infinie, c'est la direction de CA(Ω), on le note CA0(Ω)
Champs de contraintes statiquement admissibles
Champs de contraintes statiquement admissibles
Le champ de contraintes solution vérifie lui aussi certaines conditions aux limites, possède aussi une certaine régularité et satisfait aux équations d'équilibre locales du milieu :
sur où 
est le champ d'efforts surfaciques imposé. 
où 
est le champ de forces volumiques imposé.Ces conditions définissent l'espace des champs statiquement admissibles, noté SA(Ω).
Dans le cas général, cet espace est un espace affine de dimension infinie.
Champs de contraintes statiquement admissibles à zéro
De même que pour les champs cinématiquement admissibles, ont défini l'ensemble des champs statiquement admissibles à zéro comme étant l'ensemble des champs de contraintes statiquement admissibles sur Ω pour un chargement nul.
On note cet ensemble SA0(Ω). Dans le cas général, c'est un espace vectoriel de dimension infinie, qui est la direction de SA(Ω)
Recherche d'une solution approchée au problème
La recherche d'une solution approchée consiste à trouver un couple
Erreur en relation de comportement
Il vient naturellement que pour mesurer la qualité de l'approximation, il faut regarder dans quelle mesure le couple
vérifie la relation de comportement :
où 
est le tenseur de Hooke.Pour ce faire on définit l'erreur en relation de comportement, e par :
Ou encore :
![e^2 = \frac{1}{2} \int_\Omega {Tr}\left[\left(\underline{\underline{\sigma}} - \underline{\underline{\underline{\underline{K}}}}\cdot\underline{\underline{\varepsilon(u)}}\right)
\cdot\underline{\underline{\underline{\underline{K}}}}^{-1}
\cdot\left(\underline{\underline{\sigma}}- \underline{\underline{\underline{\underline{K}}}}
\cdot\underline{\underline{\varepsilon(u)}}\right)\right]
d\Omega](8/6487096f6b7b75608828169fc058b196.png)
A partir de l'expression de e2, on peut facilement aboutir aux théorèmes énergétiques.
Application aux méthodes approchées
Ces théorèmes sont utiles pour la recherche de solutions approchées. Pour ce faire on se place dans un sous-espace de dimension finie de CA(Ω) ou de SA(Ω) dont on se donne une base.
La résolution du problème de minimisation de l'énergie potentielle se ramène alors à la résolution du système linéaire suivant :
où
est la matrice de rigidité
est le déplacement approché
est le vecteur des efforts extérieurs généralisés
Les différences entre les méthodes approchées découlant de cette approche portent essentiellement sur le choix du sous-espace de recherche.
Méthode de Ritz
Cette méthode est la plus intuitive, elle consiste à choisir comme espace de recherche un espace dont la direction est un espace vectoriel définit par sa base (φi), où les φi sont des fonctions définies sur tout Ω et sont des fonctions qu'on pense assez proches de la solution exacte.
Cette méthode est pertinente si on a déjà une idée de la solution, mais n'est pas facilement programmable, et pose toujours un problème d'intégration si le domaine n'est pas simple.
Implémentation rapide en SCILAB de la méthode de Ritz
Méthode aux éléments finis
Cette méthode consiste à ne pas définir les fonctions de base définis sur tout Ω, mais à découper le domaine en éléments sur lesquels ont définit des fonctions de bases nulles à l'extérieur de l'élément.
Cette méthode à l'avantage d'être facilement programmable et permet de s'affranchir des problèmes d'intégration.
Catégories :- Mécanique des milieux continus
- Équations différentielles numériques
Wikimedia Foundation. 2010.
