Erreur en relation de comportement 736721 28803452 2008-04-19T22:27:25Z Drébon 94221 /* Notion d'erreur en relation de comportement */ ==Notion d'erreur en relation de comportement== '''L'erreur en relation de comportement''' est une notion utilisée en [[Mécanique des milieux continus|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 <math>\Omega</math> 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 <math>(\partial_1\Omega,\partial_2\Omega)</math> de <math>\partial\Omega</math>, le bord de <math>\Omega</math> telle que sur <math>\partial_1\Omega</math> il n'y ait que des déplacements imposés, et sur <math>\partial_2\Omega</math> que des efforts imposés. ===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 satisfaissant aux conditions aux limites. <math> \underline{u} = \underline{u_d}</math> sur <math>\partial_1\Omega</math> où <math>\underline{u_d}</math> est le déplacement imposé. On note cet espace <math>CA(\Omega)</math> 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 <math>\Omega</math> 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 <math>CA(\Omega)</math>, on le note <math>CA_0(\Omega)</math> ===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 : <math>\underline{\underline{\sigma}}\cdot\underline{n}=\underline{F_d}</math> sur <math>\partial_2\Omega</math> où <math>\underline{F_d}</math> est le champ d'efforts surfaciques imposé. <math>\underline{{div}}(\underline{\underline{\sigma}}) + \underline{f_d} = \rho\underline{\Gamma}</math> où <math>\underline{f_d}</math> est le champ de forces volumiques imposé. Ces conditions définissent l''''espace des champs statiquement admissibles''', noté <math>SA(\Omega)</math>. 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éfinit l'ensemble des '''champs statiquement admissibles à zéro''' comme étant l'ensemble des champs de contraintes statiquement admissibles sur <math>\Omega</math> pour un chargement nul. On note cet ensemble <math>SA_0(\Omega)</math>. Dans le cas général, c'est un '''espace vectoriel de dimension infinie''', qui est la direction de <math>SA(\Omega)</math> ==Recherche d'une solution approchée au problème== La recherche d'une solution approchée consiste à trouver un couple <math>(\underline{u},\underline{\underline{\sigma}})\in CA(\Omega)\times SA(\Omega)</math> ==Erreur en relation de comportement== Il vient naturellement que pour mesurer la qualité de l'approximation, il faut regarder dans quelle mesure le couple <math>(\underline{u},\underline{\underline{\sigma}})</math> vérifie la [[Relation de comportement|relation de comportement]] : <math>\underline{\underline{\sigma}}=\underline{\underline{\underline{\underline{K}}}} \cdot\underline{\underline{\varepsilon(u)}}</math> où <math>\underline{\underline{\underline{\underline{K}}}}</math> est le [[tenseur de Hooke]]. Pour ce faire on définit l''''erreur en relation de comportement''', <math>e</math> par : <math>e=\left\|\underline{\underline{\sigma}}- \underline{\underline{\underline{\underline{K}}}} \cdot\underline{\underline{\varepsilon(u)}}\right\|</math> Ou encore : <math> 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 </math> A partir de l'expression de <math>e^2</math>, on peut facilement aboutir aux [[Théorèmes énergétiques|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 vectoriel|sous-espace]] de dimension finie de <math>CA(\Omega)</math> ou de <math>SA(\Omega)</math> 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 : <math>\underline{\underline{K}}\cdot\underline{u}=\underline{F}</math> où *<math>\underline{\underline{K}}</math> est la [[matrice de rigidité]] *<math>\underline{u}</math> est le déplacement approché *<math>\underline{F}</math> 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 come espace de recherche un espace dont la direction est un espace vectoriel définit par sa base <math>(\varphi_i)</math>, où les <math>\varphi_i</math> sont des fonctions définies sur tout <math>\Omega</math> 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. [http://www.arbredeslemuriens.com/Categorie.php?IDCategorie=AlgoScilab&IDTitre=146 Implémentation rapide en SCILAB de la méthode de Ritz] ===Méthodes aux [[éléments finis]]=== Cette méthode consiste à ne pas définir les fonctions de base définis sur tout <math>\Omega</math>, mais à ''découper'' le domaine en '''éléments''' sur lesquels ont définit des fonctions de bases nulles à l'exterieur de l'élément. Cette méthode à l'avantage d'être facilement programmable et permet de s'affranchir des problèmes d'intégration. {{portail physique}} [[Catégorie:Mécanique des milieux continus]] [[Catégorie:Équations différentielles numériques]]