Théorème de réciprocité de Maxwell-Betti

Le théorème de réciprocité de Maxwell-Betti a été énoncé en 1872 par Enrico Betti. Il énonce que

« Si, à un système élastique quelconque, isostatique ou hyperstatique, on applique successivement un premier système de forces, puis après suppression, un second, la somme des travaux du premier système, pour les déplacements élastiques dus aux forces du second, est égale à la somme des travaux des forces du second système pour les déplacements élastiques dus aux forces du premier. »

— B. de Fontviolant, Résistance des Matériaux analytique et graphique^[1].

Étant donné que les déformations élastiques sont réversibles, l'énoncé porte sur la comparaison de deux états d'un même système, indépendamment du temps ou de toute précédence. Les travaux d'un système de forces, ou chargement, sous les déplacements dus à l'autre chargement, sont appelés travaux réciproques. On montre que cet énoncé vaut aussi bien pour des couples, qui travaillent sous une rotation, que pour des contraintes travaillant sous un certain champ de déformation^[2].

Le théorème de Maxwell-Betti revêt une grande importance en Résistance des Matériaux, surtout dans l'analyse des réseaux de poutres. C'est aussi l'une des sources d'inspiration de la méthode des éléments finis^[3]. G. Colonnetti a établi un énoncé voisin par la suite.

Exemples modifier

Résistance des Matériaux modifier

Démonstration de l'énoncé de Maxwell-Betti pour une poutre console.

Considérons une poutre horizontale, sur laquelle on choisit deux points 1 et 2 arbitrairement, à cela près qu'ils doivent être distincts des points d'appui (ce qui donnerait, cela dit, un résultat exact, mais trivial). Appliquons une force verticale P au point 1 et notons $\Delta _{P2}$ la flèche au point 2 ; puis, ayant supprimé la charge P, appliquons une force Q au point 2, ce qui produit une flèche $\Delta _{Q1}$ au point 1. Le théorème de Maxwell-Betti indique que dans ces conditions :

P\cdot \Delta _{Q1}=Q\cdot \Delta _{P2}.

Système masse-ressort modifier

Système masse-ressort linéaire.

Soient deux masses, de poids respectifs F₁ et F₂, reliées par deux ressorts de raideurs respectives k₁ et k₂ et suspendues au plafond du laboratoire. À l'équilibre, l'allongement des ressorts est donné par:

\underbrace {\begin{pmatrix}k_{1}+k_{2}&-k_{2}\\-k_{2}&k_{2}\end{pmatrix}} _{=:\mathbf {K} }{\begin{pmatrix}u_{1}\\u_{2}\end{pmatrix}}={\begin{pmatrix}F_{1}\\F_{2}\end{pmatrix}}\leftrightarrow {\begin{pmatrix}u_{1}\\u_{2}\end{pmatrix}}={\frac {1}{k_{1}k_{2}}}{\begin{pmatrix}k_{2}&k_{2}\\k_{2}&k_{1}+k_{2}\end{pmatrix}}{\begin{pmatrix}F_{1}\\F_{2}\end{pmatrix}}

Considérons un premier état, où le poids F₁ est négligeable devant F₂. Le déplacement ${\vec {u}}$ sera donné par

{\vec {F}}={\begin{pmatrix}F_{1}\\0\end{pmatrix}}\rightarrow {\vec {u}}={\begin{pmatrix}u_{1}\\u_{2}\end{pmatrix}}={\frac {F_{1}}{k_{1}}}{\begin{pmatrix}1\\1\end{pmatrix}}

Puis un second état, où le poids F₂ est négligeable devant F₁. Le déplacement ${\vec {\tilde {u}}}$ sera donné par :

{\vec {\tilde {F}}}={\begin{pmatrix}0\\{\tilde {F}}_{2}\end{pmatrix}}\rightarrow {\vec {\tilde {u}}}={\begin{pmatrix}{\tilde {u}}_{1}\\{\tilde {u}}_{2}\end{pmatrix}}={\frac {{\tilde {F}}_{2}}{k_{1}}}{\begin{pmatrix}1\\1+{\frac {k_{1}}{k_{2}}}\end{pmatrix}}

Le théorème de Betti exprime que les déplacements, dans les deux configurations précédentes, sont reliés par :

F_{1}{\tilde {u}}_{1}={\frac {F_{1}{\tilde {F}}_{2}}{k_{1}}}={\tilde {F}}_{2}u_{2}

Le travail s'exprime dans chaque cas par ${\vec {F}}=\mathbf {K} {\vec {u}}$ et ${\vec {\tilde {F}}}=\mathbf {K} {\vec {\tilde {u}}}$ . En composant ces deux relations, on obtient :
${\vec {F}}^{\top }{\vec {\tilde {u}}}={\vec {\tilde {u}}}^{\top }{\vec {F}}={\vec {\tilde {u}}}^{\top }\mathbf {K} {\vec {u}}={\vec {\tilde {u}}}^{\top }\mathbf {K} ^{\top }{\vec {u}}=(\mathbf {K} {\vec {\tilde {u}}})^{\top }{\vec {u}}={\vec {\tilde {F}}}^{\top }{\vec {u}}$

La symétrie de la matrice de raideur K est nécessaire pour que les travaux réciproques soient égaux.

En théorie de l'élasticité modifier

Considérons un solide élastique, occupant un volume V et de surface A extérieure, chargé par des forces de contact ${\vec {T}}$ et une densité de forces volumiques ${\vec {f}}$ (par exemple la gravité). Le triplet $\{{\vec {u}},{\boldsymbol {\varepsilon }},{\boldsymbol {\sigma }}\}$ formé du champ de déplacement, du tenseur des déformations et du tenseur des contraintes définit un état élastique du solide, est associé au chargement extérieur $\{{\vec {T}},{\vec {f}}\}$ par la conservation de l'énergie mécanique :

\int _{A}{\vec {T}}\cdot {\vec {u}}\,\mathrm {d} A+\int _{V}{\vec {f}}\cdot {\vec {u}}\,\mathrm {d} V=\int _{V}{\boldsymbol {\sigma }}:{\boldsymbol {\varepsilon }}\,\mathrm {d} V

Le symbole „ $\cdot$ “ désigne le produit scalaire de deux vecteurs et le produit doublement contracté « : » (ou « produit de Frobenius »), est une opération entre deux tenseurs A et B que l'on peut définir comme la trace A : B := tr(A^T · B). Le travail développé par les forces extérieures ${\vec {T}}$ et ${\vec {f}}$ sous le déplacement ${\vec {u}}$ est donc égal au travail de déformation des contraintes σ sous les déformations ε.

Les déplacements et les déformations dépendent de ${\boldsymbol {\varepsilon }}={\tfrac {1}{2}}(\operatorname {grad} ({\vec {u}})+\operatorname {grad} ({\vec {u}})^{\top })$ . Le mot « grad » désigne ici le gradient, et le T en exposant désigne la transposition. Du tenseur de déformation ε, on peut calculer le tenseur des contraintes σ grâce au tenseur d'élasticité symétrique :

\mathbb {C} =\mathbb {C} ^{\top }\quad \rightarrow \quad {\boldsymbol {\sigma }}:=\mathbb {C} :{\boldsymbol {\varepsilon }}={\boldsymbol {\varepsilon }}:\mathbb {C}

Le produit „:“ contracté d'un tenseur de quatrième ordre ( $\mathbb {C}$ ) avec un tenseur de deuxième ordre (ε) donne un tenseur de deuxième ordre (σ). L'isotropie de comportement dictée par la loi de Hooke donne

\mathbb {C} =2\mu \mathbb {I} +\lambda \mathbf {1\otimes 1} \quad \rightarrow \quad {\boldsymbol {\sigma }}=\mathbb {C} :{\boldsymbol {\varepsilon }}={\boldsymbol {\varepsilon }}:\mathbb {C} =2\mu {\boldsymbol {\varepsilon }}+\lambda \operatorname {Sp} ({\boldsymbol {\varepsilon }})\mathbf {1}

Les coefficients $\{\mu ,\lambda \}$ sont les Coefficient de Lamé, $\{\mathbf {1} ,\mathbb {I} \}$ désigne le tenseur identité de deuxième (resp. de quatrième) ordre tous deux symétriques, de sorte qu'on dispose des identités $\mathbf {1} :{\boldsymbol {\varepsilon }}={\boldsymbol {\varepsilon }}:\mathbf {1} =\operatorname {Sp} ({\boldsymbol {\varepsilon }})$ et $\mathbb {I} :{\boldsymbol {\varepsilon }}={\boldsymbol {\varepsilon }}:\mathbb {I} ={\boldsymbol {\varepsilon }}$ . Le théorème de Betti n'exige ni l'isotropie, ni l'homogénéité du matériau.

Soit alors $\{{\vec {\tilde {u}}},{\tilde {\boldsymbol {\varepsilon }}},{\tilde {\boldsymbol {\sigma }}}\}$ un second état d'équilibre élastique du solide relatif au chargement $\{{\vec {\tilde {T}}},{\vec {\tilde {f}}}\}$ . Le théorème de Betti énonce que :

\int _{A}{\vec {T}}\cdot {\vec {\tilde {u}}}\,\mathrm {d} A+\int _{V}{\vec {f}}\cdot {\vec {\tilde {u}}}\,\mathrm {d} V=\int _{V}{\boldsymbol {\sigma }}:{\tilde {\boldsymbol {\varepsilon }}}\,\mathrm {d} V=\int _{V}{\tilde {\boldsymbol {\sigma }}}:{\boldsymbol {\varepsilon }}\,\mathrm {d} V\,.=\int _{A}{\vec {\tilde {T}}}\cdot {\vec {u}}\,\mathrm {d} A+\int _{V}{\vec {\tilde {f}}}\cdot {\vec {u}}\,\mathrm {d} V

La symétrie du tenseur d'élasticité $\mathbb {C}$ , qui résulte de l'hyperélasticité, constitue une hypothèse essentielle^[2].

Théorème du travail maximal modifier

Le théorème de Maxwell-Betti est une conséquence du théorème du travail maximal que nous allons commencer par démontrer.

Soit ${\vec {u}}$ un champ de déplacement intérieur au solide, et soit ${\boldsymbol {\varepsilon }}={\tfrac {1}{2}}(\operatorname {grad} ({\vec {u}})+\operatorname {grad} ({\vec {u}})^{\top })$ le champ de déformation correspondant. Considérons un champ de contrainte ${\tilde {\boldsymbol {\sigma }}}$ , indépendant du champ de déformation précédent, mais satisfaisant à la condition d'équilibre $\operatorname {div} {\tilde {\boldsymbol {\sigma }}}+{\vec {\tilde {f}}}={\vec {0}}$ et compatible avec les conditions de chargement de surface ${\vec {\tilde {T}}}={\tilde {\boldsymbol {\sigma }}}\cdot {\hat {n}}$ , où ${\vec {\tilde {T}}}$ est le vecteur contrainte développé par le chargement en surface et où ${\hat {n}}$ est le vecteur normal unité à la surface (dirigé par convention vers l'extérieur et noté avec un ^). Le travail réciproque des contraintes de surfaces sous le champ de déplacement ${\vec {u}}$ s'exprime avec le tenseur de contrainte :

\int _{A}{\vec {\tilde {T}}}\cdot {\vec {u}}\;\mathrm {d} A=\int _{A}({\tilde {\boldsymbol {\sigma }}}\cdot {\hat {n}})\cdot {\vec {u}}\;\mathrm {d} A=\int _{A}({\tilde {\boldsymbol {\sigma }}}\cdot {\vec {u}})\cdot {\hat {n}}\;\mathrm {d} A

Si la surface est suffisamment régulière (sans points anguleux), cette intégrale de surface peut être transformée en intégrale de volume grâce au théorème de flux-divergence $\textstyle \int _{A}{\vec {v}}\cdot {\hat {n}}\,\mathrm {d} A=\int _{V}\operatorname {div} ({\vec {v}})\,\mathrm {d} V$ :

\int _{A}{\vec {\tilde {T}}}\cdot {\vec {u}}\;\mathrm {d} A=\int _{A}({\tilde {\boldsymbol {\sigma }}}\cdot {\vec {u}})\cdot {\hat {n}}\;\mathrm {d} A=\int _{V}\operatorname {div} ({\tilde {\boldsymbol {\sigma }}}\cdot {\vec {u}})\;\mathrm {d} V

et la formule de Leibnitz donne $\operatorname {div} (\mathbf {T} \cdot {\vec {v}})=\operatorname {div} (\mathbf {T} )\cdot {\vec {v}}+\mathbf {T} ^{\top }:\operatorname {grad} ({\vec {v}})$

\int _{A}{\vec {\tilde {s}}}\cdot {\vec {u}}\;\mathrm {d} A=\int _{V}\operatorname {div} ({\tilde {\boldsymbol {\sigma }}}\cdot {\vec {u}})\;\mathrm {d} V=\int _{V}\operatorname {div} ({\tilde {\boldsymbol {\sigma }}})\cdot {\vec {u}}\;\mathrm {d} V+\int _{V}{\tilde {\boldsymbol {\sigma }}}:\operatorname {grad} ({\vec {u}})\;\mathrm {d} V

En y substituant la conditions d'équilibre $\operatorname {div} {\tilde {\boldsymbol {\sigma }}}=-{\vec {\tilde {f}}}$ et en tenant compte du fait que la composante antisymétrique du gradient de déplacement ne travaille pas dans le produit contracté avec le tenseur (symétrique) des contraintes, c'est-à-dire qu'on a l'identité ${\tilde {\boldsymbol {\sigma }}}:\operatorname {grad} ({\vec {u}})={\tilde {\boldsymbol {\sigma }}}:{\boldsymbol {\varepsilon }}$ , on obtient le théorème de l’énergie élastique :

\int _{A}{\vec {\tilde {T}}}\cdot {\vec {u}}\;\mathrm {d} A+\int _{V}{\vec {\tilde {f}}}\cdot {\vec {u}}\;\mathrm {d} V=\int _{V}{\tilde {\boldsymbol {\sigma }}}:{\boldsymbol {\varepsilon }}\;\mathrm {d} V

Pour que ce résultat soit valable, il faut^[2] que :

le champ de déplacement ${\vec {u}}$ soit admissible, c’est-à-dire
- qu’il soit deux fois différentiable dans l’intérieur du solide élastique linéaire, et
- qu’il soit continu ainsi que son gradient $\operatorname {grad} {\vec {u}}$ dans l'intérieur et à la surface du solide (boule fermée) [V] ;
et que le tenseur des contraintes soit admissible, c’est-à-dire
- qu’il soit continûment dérivable dans l'intérieur du volume V,
- qu’il soit continu, ainsi que sa divergence, sur le volume fermé [V].

Dans l'énoncé du théorème de l'énergie élastique, le champ de contrainte et le champ de déformation sont mutuellement indépendants et ne sont donc pas nécessairement liés par une loi de comportement. Les forces extérieures $\{{\vec {\tilde {T}}},{\vec {\tilde {f}}}\}$ exercent donc le même travail sous le champ de déplacement, que les contraintes internes sous le champ de déformation élastique.

Théorème de réciprocité de Maxwell-Betti modifier

Considérons un second état de déformation correspondant à l'action de forces extérieures $\{{\vec {T}},{\vec {f}}\}$ . Le matériau élastique, en se déformant, développe un champ de contraintes ${\boldsymbol {\sigma }}$ et le déplacement moléculaire ${\vec {\tilde {u}}}$ un champ de déformation ${\tilde {\boldsymbol {\varepsilon }}}={\tfrac {1}{2}}(\operatorname {grad} ({\vec {\tilde {u}}})+\operatorname {grad} ({\vec {\tilde {u}}})^{\top })$ . Comme pour l'état précédent, le théorème de l'énergie élastique s'écrit :

\int _{A}{\vec {T}}\cdot {\vec {\tilde {u}}}\;\mathrm {d} A+\int _{V}{\vec {f}}\cdot {\vec {\tilde {u}}}\;\mathrm {d} V=\int _{V}{\boldsymbol {\sigma }}:{\tilde {\boldsymbol {\varepsilon }}}\;\mathrm {d} V

Cette relation reste valable car $\{{\vec {u}},{\boldsymbol {\varepsilon }},{\boldsymbol {\sigma }}\}$ et $\{{\vec {\tilde {u}}},{\tilde {\boldsymbol {\varepsilon }}},{\tilde {\boldsymbol {\sigma }}}\}$ sont deux états d'équilibre élastique du solide.

Alors, le tenseur d'élasticité étant symétrique, ${\boldsymbol {\sigma }}=\mathbb {C} :{\boldsymbol {\varepsilon }}={\boldsymbol {\varepsilon }}:\mathbb {C}$ et ${\tilde {\boldsymbol {\sigma }}}=\mathbb {C} :{\tilde {\boldsymbol {\varepsilon }}}$ de sorte que :

{\boldsymbol {\sigma }}:{\tilde {\boldsymbol {\varepsilon }}}=({\boldsymbol {\varepsilon }}:\mathbb {C} ):{\tilde {\boldsymbol {\varepsilon }}}={\boldsymbol {\varepsilon }}:(\mathbb {C} :{\tilde {\boldsymbol {\varepsilon }}})={\boldsymbol {\varepsilon }}:{\tilde {\boldsymbol {\sigma }}}={\tilde {\boldsymbol {\sigma }}}:{\boldsymbol {\varepsilon }}

On retrouve bien l'énoncé continu du théorème de Maxwell-Betti :

\int _{A}{\vec {T}}\cdot {\vec {\tilde {u}}}\;\mathrm {d} A+\int _{V}{\vec {f}}\cdot {\vec {\tilde {u}}}\;\mathrm {d} V=\int _{V}{\boldsymbol {\sigma }}:{\tilde {\boldsymbol {\varepsilon }}}\;\mathrm {d} V=\int _{V}{\tilde {\boldsymbol {\sigma }}}:{\boldsymbol {\varepsilon }}\;\mathrm {d} V=\int _{A}{\vec {\tilde {T}}}\cdot {\vec {u}}\;\mathrm {d} A+\int _{V}{\vec {\tilde {f}}}\cdot {\vec {u}}\;\mathrm {d} V

Là encore, quoique le matériau puisse n'être ni homogène, ni isotrope, la symétrie du tenseur d'élasticité est une condition nécessaire essentielle^[2].

Notes modifier

↑ Chap. XIV Théorèmes découlant de l'équation générale de l’Élasticité, t. I, p. 205.
↑ ^{a b c et d} M. E. Gurtin et C. Truesdell et S. Flügge (dir.), Handbuch der Physik, vol. VI2/a : The Linear Theory of Elasticity, Springer, 1972 (ISBN 3-540-05535-5), p. 98 f
↑ (de) Thorsten Kunow, « Verbesserte Berechnung von lokalen Zielgrößen mit der Methode der finiten Elemente unter Verwendung von Grundlösungen » [PDF] (consulté le 25 août 2016)

Voir également modifier

Christian La Borderie, Méthodes Énergétiques

[Bourgault1-1] Chap. XIV Théorèmes découlant de l'équation générale de l’Élasticité, t. I, p. 205.

[gurtin-2] {a b c et d} M. E. Gurtin et C. Truesdell et S. Flügge (dir.), Handbuch der Physik, vol. VI2/a : The Linear Theory of Elasticity, Springer, 1972 (ISBN 3-540-05535-5), p. 98 f

[3] (de) Thorsten Kunow, « Verbesserte Berechnung von lokalen Zielgrößen mit der Methode der finiten Elemente unter Verwendung von Grundlösungen » [PDF] (consulté le 25 août 2016)

[1]

[2]

[3]