Théorèmes énergétiques - Définition

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

Introduction

Les théorèmes énergétiques permettant de poser un problème de mécanique sous forme d'un problème de minimisation, et donc d'utiliser toutes les méthodes existantes pour ce genre de problèmes. Il y a deux théorèmes énergétiques, l'un est beaucoup plus connu que l'autre, il s'agit du théorème de l'énergie potentielle. Le second théorème est similaire au premier, la seule différence étant que dans un cas on parle de l'énergie potentielle et dans l'autre de l'énergie complémentaire.

Théorème de l'énergie potentielle

Le champ de déplacement solution d'un problème statique de mécanique des milieux continus minimise l'énergie potentielle du milieu.

Exemples

ExemplethmEN1.svg

Soit un solide sur un plan sans frottement, relié à un ressort de raideur 'k' et soumis à une force 'F'.

On repère la position de référence du solide par 'x=0', et on suppose que cette position coïncide avec la position de repos (ou longueur à vide) du ressort. On note 'T' la tension du ressort.

Résolution en utilisant le théorème de l'énergie potentielle

L'énergie potentielle du système s'écrit :

 E_p = \frac{1}{2} k x^2 - F x

La minimisation de cette énergie donne :

0=\frac{\partial E_p}{\partial x}=k x -F

 x=\frac{F}{k}

Résolution par l'énergie complémentaire

L'énergie complémentaire du système s'écrit :

 E_c= \frac{1}{2k}T^2

Ici comme l'on se doit de chercher 'T' dans l'espace des champs admissibles, on a du fait de l'effort 'F' imposé, 'T' = '-F', et comme on a un problème à 1 ddl, on a la solution.

Théorème de l'énergie complémentaire

Le champ de déplacement solution d'un problème statique de mécanique des milieux continus minimise l'énergie complémentaire du milieu

Démonstrations

À partir de l'erreur en relation de comportement

Établissement de l'expression de l'énergie potentielle et de l'énergie complémentaire

Soit e l'erreur en relation de comportement. Partons de l'expression de e2. On a :

 e^2\left(\underline{u},\underline{\underline{\sigma}}\right) = \frac{1}{2} \int_\Omega \operatorname{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] \mathrm d\Omega =\frac{1}{2}\left[ \int_\Omega \operatorname{Tr}\left(\underline{\underline{\sigma}}\cdot \underline{\underline{\underline{\underline{K}}}}^{-1} \cdot\underline{\underline{\sigma}}\right)\mathrm d\Omega +  \int_\Omega \operatorname{Tr}\left(\underline{\underline{\underline{\underline{K}}}}\cdot  \underline{\underline{\varepsilon}} \cdot\underline{\underline{\varepsilon}}\right) \mathrm d\Omega - 2 \int_\Omega \operatorname{Tr}\left(\underline{\underline{\sigma}}\cdot \underline{\underline{\varepsilon}}\right)\mathrm d\Omega \right]

On applique alors la formule de Green :

 \int_\Omega \operatorname{Tr}\left(\underline{\underline{\sigma}} \cdot\underline{\underline{\varepsilon}} \right)\mathrm d\Omega =\int_{\partial\Omega} \left(\underline{\underline{\sigma}}\cdot\underline{n}\right)\cdot\underline{u}\, \mathrm dS -\int_\Omega  \underline\operatorname{div} \left(\underline{\underline{\sigma}}\right)\cdot\underline{u}\, \mathrm d\Omega

On peut séparer l'intégrale sur le bord de Ω en deux intégrales : une sur \partial_1\Omega , partie de \partial\Omega sur laquelle on a imposé le déplacement et une sur \partial_2\Omega , partie de \partial\Omega sur laquelle on a imposé l'effort. On rappelle que dans ce cas on a :

\partial\Omega=\partial_1\Omega\cup\partial_2\Omega et \partial_1\Omega\cap\partial_2\Omega=\empty

On obtient donc, en utilisant les équations d'équilibre et les conditions aux limites :

 \int_\Omega \operatorname{Tr}\left(\underline{\underline{\sigma}} \cdot\underline{\underline{\varepsilon}} \right)\mathrm d\Omega =\int_{\partial_1\Omega} \left(\underline{\underline{\sigma}}\cdot\underline{n}\right)\cdot\underline{u_d}\, \mathrm dS +\int_{\partial_2\Omega} \underline{F_d}\cdot\underline{u}\, \mathrm dS +\int_\Omega \underline{f_d}\cdot\underline{u}\,\mathrm d\Omega

Soit

 e^2\left(\underline{u},\underline{\underline{\sigma}}\right) = E_p(\underline{u}) + E_c(\underline{\underline{\sigma}})

Avec :

E_p(\underline{u})=\frac{1}{2} \int_\Omega \operatorname{Tr}\left(\underline{\underline{\underline{\underline{K}}}}\cdot\underline{\underline{\varepsilon}} ^2\right) \mathrm d\Omega - \left(\int_\Omega \underline{f_d}\cdot\underline{u}\,\mathrm d\Omega + \int_{\partial_2\Omega} \underline{F_d}\cdot\underline{u}\, \mathrm dS\right)
E_c(\underline{\underline{\sigma}})= \frac{1}{2} \int_\Omega \operatorname{Tr}\left(\underline{\underline{\sigma}}\cdot\underline{\underline{\underline{\underline{K}}}}^{-1} \cdot\underline{\underline{\sigma}}\right)\mathrm d\Omega - \int_{\partial_1\Omega} \left(\underline{\underline{\sigma}}\cdot\underline{n}\right)\cdot\underline{u_d}\, \mathrm dS

Étude des différents termes de ces expressions

Énergie potentielle
\frac{1}{2} \int_\Omega \operatorname{Tr}\left(\underline{\underline{\underline{\underline{K}}}}\cdot\underline{\underline{\varepsilon}} ^2\right) \mathrm d\Omega

est l'énergie de déformation exprimée en déplacement. C'est l'analogue de l'énergie de déformation élastique d'un ressort \tfrac{1}{2} k (l-l_0)^2 . C'est une forme bilinéaire symétrique continue et coercive.

\left(\int_\Omega \underline{f_d}\cdot\underline{u}\,\mathrm d\Omega + \int_{\partial_2\Omega} \underline{F_d}\cdot\underline{u}\, \mathrm dS\right)

est le travail des efforts imposés (donc connues) dans le champ de déplacement inconnu. C'est une forme linéaire continue.

Énergie complémentaire
\frac{1}{2} \int_\Omega \operatorname{Tr}\left(\underline{\underline{\sigma}}\cdot\underline{\underline{\underline{\underline{K}}}}^{-1} \cdot\underline{\underline{\sigma}}\right)\mathrm d\Omega

est l'énergie de déformation exprimée en contrainte, c'est l'analogue de l'énergie de déformation élastique d'un ressort exprimée en fonction de la force \tfrac{1}{2}\tfrac{F^2}{k} . C'est une forme bilinéaire symétrique continue et coercive.

\int_{\partial_1\Omega} \left(\underline{\underline{\sigma}}\cdot\underline{n}\right)\cdot\underline{u_d}\, \mathrm dS

est le travail des efforts inconnus dans le champ de déplacement imposé (donc connu) Ep(u) ne dépendant que de u et des forces imposées. C'est une forme linéaire continue.

Démonstration des théorèmes énergétiques

Nous savons que la solution exacte est telle que e^2(\underline{u},\underline{\underline{\sigma}})=0 , donc elle réalise un minimum de cette fonction. De plus comme  e^2\left(\underline{u},\underline{\underline{\sigma}}\right) est la somme de deux fonctions de variables indépendantes, le couple (\underline{u},\underline{\underline{\sigma}}) est tel que \underline{u} réalise le minimum de E_p(\underline{u}) et \underline{\underline{\sigma}} réalise le minimum de E_c(\underline{\underline{\sigma}}) .

De plus nous sommes dans les hypothèses d'application du théorème de Stampacchia qui garantit l'existence et l'unicité des solutions de chacun de ces problèmes de minimisation.

De plus le déplacement et la contrainte étant liés par la relation de comportement, on peut choisir de minimiser l'une ou l'autre de ces énergies pour obtenir le couple solution.

À partir de l'équation d'équilibre local

On peut faire des formulations variationnelles en déplacement v CA0 (ci dessous) ou contrainte to CA 0 et obtenir a(u,v)=l(v) ou a(sigma, to)=l(to) équations à résoudre.

Énergie potentielle, formulation en déplacement

on cherche \underline{u} , cinématiquement admissible, solution de

\underline{div}(\underline{\underline{\sigma}})+\underline{f}_d=\underline{0}

vérifiant la relation de comportement  \underline{\underline{\sigma}}=\underline{\underline{\underline{\underline{K}}}}.\underline{\underline{\varepsilon}}(\underline{u})

On multiplie par un champ cinématiquement admissible à zéro \underline{u}^\star :

\underline{div}(\underline{\underline{\sigma}}).\underline{u}^\star+\underline{f}_d.\underline{u}^\star=0,\ \forall \underline{u}^\star \ CA0


Un intégration donne :

\int_{\Omega}\underline{div}(\underline{\underline{\sigma}}).\underline{u}^\star d\Omega+\int_{\Omega}\underline{f}_d.\underline{u}^\star d\Omega=0

La formule de Green permet d'arriver à :

\int_{\Omega}\underline{div}(\underline{\underline{\sigma}}).\underline{u}^\star d\Omega=\int_{\partial\Omega}(\underline{\underline{\sigma}}.\underline{n}).\underline{u}^\star d\Gamma-\int_{\Omega}\sigma:\underline{\underline{\varepsilon}}(\underline{u}^\star)d\Omega

Soit:

\int_{\Omega}\underline{\underline{\sigma}}:\underline{\underline{\varepsilon}}(\underline{u}^\star) d\Omega =  \int_{\Omega}\underline{f}_d.\underline{u}^\star d\Omega + \int_{\partial\Omega}(\underline{\underline{\sigma}}.\underline{n}).\underline{u}^\star=0

On décompose le bord du domaine en deux: lieux à efforts imposés {\partial\Omega}_F  : \partial\underline{\underline{\sigma}}.\underline{n}=\underline{F}_d , lieu à déplacements imposés {\partial\Omega}_u  : \underline{u}=\underline{u}_d et \underline{u}^\star=0

On obtient alors :

\int_{\Omega}\underline{f}_d.\underline{u}^\star d\Omega + \int_{{\partial\Omega}_F}\underline{F}_d.\underline{u}^\star d\Gamma -\int_{\Omega}\underline{\underline{\sigma}}:\underline{\underline{\varepsilon}}(\underline{u}^\star)=0

En utilisant la relation de comportement, on aboutit au problème suivant :

trouver  \underline{u} cinématiquement admissible tel que :

 a(\underline{u},\underline{u}^\star)=l(\underline{u}^\star),\ \forall \underline{u}^\star \ CA0

Avec a(\underline{u},\underline{v})=\int_{\Omega} \underline{\underline{\underline{\underline{K}}}}.\underline{\underline{\varepsilon}}(\underline{u})\underline{\underline{\varepsilon}}(\underline{v})d\Omega et  l(\underline{v})=\int_\Omega \underline{f}_d.\underline{v}d\Omega + \int_{{\partial\Omega}_F} \underline{F}_d.\underline{v}d\Gamma

On peut alors vérifier que l'on est dans le cadre d'application des théorèmes de Lax-Milgram et de Stampacchia, qui assurent d'une part, l'existence et l'unicité de la solution, d'autre par l'équivalence de ce problème au problème de minimisation de la fonctionnelle I(\underline{u}) définie par :

 I(\underline{u})=\frac{1}{2} a(\underline{u},\underline{u}) - l(\underline{u})

Cette fonctionnelle n'est autre que l'énergie potentielle.

Énergie complémentaire, formulation en contrainte

L'établissement du théorème de l'énergie complémentaire se fait de la même façon, mais en posant une formulation variationnelle en contrainte, pour cela il faut partir de la relation de comportement et la multiplier par un champ de contraintes statiquement admissible à zéro.

Page générée en 0.154 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