Tenseur des contraintes
Le tenseur des contraintes est un tenseur d'ordre 2 utilisé en mécanique des milieux continus pour caractériser l'état de contrainte, c'est-à-dire les efforts intérieurs mis en jeu entre les portions déformées d'un milieu. Le terme a été introduit par Cauchy vers 1822.
Comme les efforts intérieurs sont définis pour chaque surface coupant le milieu (on parle d'ailleurs également d'efforts surfaciques), le tenseur est défini localement, en chaque point du solide. L'état de contrainte du solide est donc représenté par un champ tensoriel. On parle aussi de ce fait de champ de contrainte.
Dans le cadre de l'élasticité linéaire, le champ de contrainte est relié au champ de déformation par la loi de Hooke généralisée, c'est-à-dire que l'on peut écrire l'équation tensorielle Modèle:Mvar où Modèle:Mvar désigne le tenseur des contraintes, Modèle:Mvar le tenseur des déformations et Modèle:Mvar le tenseur des constantes élastiques.
Dans le cadre de la géologie structurale et de la tectonique, on parle fréquemment de tenseur des paléo-contraintes. Il représente la partie anisotrope du tenseur des contraintes, responsable des déformations comme les plis, les failles ou les schistosités. La valeur absolue des termes de la matrice n'est pas accessible, mais il est possible de retrouver l'orientation du triaxe principal, ainsi que le rapport d'intensité entre ces trois axes.
Dans certains cas, il est possible de visualiser ces contraintes par la méthode de photoélasticimétrie.
Construction du tenseur
Dans le cas de certaines sollicitations simples (traction ou compression uniaxiale, flexion, torsion), la contrainte peut être représentée par un simple nombre (scalaire). Dans le cas général, on peut définir une contrainte équivalente (ou effective) qui est aussi un scalaire, mais qui ne définit que très incomplètement l'état de contrainte. Pour le décrire complétement il faut utiliser un tenseur d'ordre 2 ou, dans un référentiel donné, une matrice 3×3.
Prenons une base <math>(\vec{e_1},\vec{e_2},\vec{e_3})</math> et un point Modèle:Math de la pièce. Considérons un cube de matière autour de Modèle:Math, d'arête infinitésimale Modèle:Math, et dont les arêtes sont parallèles aux axes du repère.
Numérotons ses faces :
- les faces Modèle:Mvar et Modèle:Mvar sont les faces normales à <math>\vec{e_i}</math>, en partant du centre du cube, <math>\vec{e_i}</math> pointe vers Modèle:Mvar, la face Modèle:Mvar étant la face opposée.
Nous ne considérons que les faces numérotées positivement. En effet par la troisième loi de Newton (principe d'action-réaction), les forces sur des faces opposées sont elles aussi opposées.
Sur la face j s'exerce un vecteur-force <math>\vec{\mathrm{F}}_j</math> qui a trois composantes :
- <math>\vec{\mathrm{F}}_j=\begin{pmatrix}F_{1j} \\ F_{2j} \\ F_{3j}\end{pmatrix}</math>,
Modèle:Mvar étant la composante selon <math>\vec{e_i}</math> du vecteur-force s'exerçant sur la face Modèle:Mvar.
La surface de chaque facette étant Modèle:Math, on peut définir neuf composantes Modèle:Mvar homogènes à des contraintes :
- <math>\sigma_{ij} = \frac{F_{ij}}{a^2}</math>
On décrit donc l'état de contrainte par le tenseur :
- <math>\mathrm{T}(\mathrm{M})=\begin{pmatrix}
\sigma_{11} & \sigma_{12} & \sigma_{13}\\ \sigma_{21} & \sigma_{22} & \sigma_{23}\\ \sigma_{31} & \sigma_{32} & \sigma_{33}\\ \end{pmatrix}</math> Modèle:Math est un tenseur d'ordre 2, à 3 lignes et 3 colonnes. Il est défini localement pour un point Modèle:Math donné.
En mécanique, on n'utilise pas toujours la notation généralisée <math>(\vec{e_i})</math> pour la base. Si l'on note la base <math>(\vec{x}, \vec{y}, \vec{z})</math>, les composantes du tenseur se notent alors :
- <math>\mathrm{T}(\mathrm{M})=\begin{pmatrix}
\sigma_{xx} & \sigma_{xy} & \sigma_{xz}\\ \sigma_{yx} & \sigma_{yy} & \sigma_{yz}\\ \sigma_{zx} & \sigma_{zy} & \sigma_{zz}\\ \end{pmatrix}</math> Les termes hors diagonale correspondant à du cisaillement, on les note souvent Modèle:Mvar, les composantes du tenseur se notent alors :
- <math>\mathrm{T}(\mathrm{M})=\begin{pmatrix}
\sigma_{xx} & \tau_{xy} & \tau_{xz}\\ \tau_{yx} & \sigma_{yy} & \tau_{yz}\\ \tau_{zx} & \tau_{zy} & \sigma_{zz}\\ \end{pmatrix}</math>
Conventions de signe
L'exposé ci-dessus est le plus naturel en termes de construction mathématique. On remarque que les contraintes normales (les éléments diagonaux du tenseur) sont positives quand le cube considéré subit une extension, et négatives quand il subit une compression. Ce choix est généralement adopté dans les articles et ouvrages de mécanique des milieux continus.
On trouve toutefois aussi la convention de signe contraire (contraintes normales positives en compression, négatives en extension), notamment dans les articles et ouvrages de sciences de la Terre, les roches situées à grande profondeur étant presque toujours en compression. Le signe des éléments non diagonaux du tenseur est également permuté. Dans l'exposé ci-dessus cela correspond à définir <math>F_{ij}</math> comme la valeur algébrique de la force <math>\vec\mathrm F_{ij}</math> exercée sur la face Modèle:Mvar au lieu de Modèle:Mvar.
La convention de signe n'est pas toujours indiquée dans les textes, mais il est facile de voir laquelle a été adoptée à partir des commentaires du texte, notamment quand il est question de pression.
Symétrie du tenseur des contraintes
Le tenseur des contraintes est toujours symétrique, c'est-à-dire que :
- <math>\sigma_{ij} = \sigma_{ji}</math>
ce qui traduit l'équilibre en moment d'un volume infinitésimal.
Le tenseur s'écrit donc :
- <math>\mathrm{T}_{ij} = \begin{pmatrix}
\sigma_{11} & \sigma_{12} & \sigma_{13}\\ \sigma_{12} & \sigma_{22} & \sigma_{23}\\ \sigma_{13} & \sigma_{23} & \sigma_{33}\\ \end{pmatrix} </math>.
Du fait de cette symétrie, on peut écrire le tenseur comme un vecteur, selon la notation de Voigt : en définissant
On peut alors mettre le tenseur sous la forme :
- <math>\mathrm{T}_{i} = \begin{pmatrix}
\sigma_{1} \\ \sigma_{2} \\ \sigma_{3} \\ \sigma_{4} \\ \sigma_{5} \\ \sigma_{6} \\ \end{pmatrix} </math>. Ceci facilite l'écriture de la loi de Hooke généralisée. On voit aussi que l'espace des contraintes est un espace vectoriel à six dimensions.
Contrainte moyenne et pression isostatique
La contrainte moyenne est la moyenne arithmétique des contraintes normales, donc le tiers de la trace du tenseur :
- <math>\bar\sigma= \frac{\mathrm{tr}(\mathrm{T})}{3} = \frac{\sigma_{11} + \sigma_{22} + \sigma_{33}}{3}</math>
Sa valeur ne dépend pas du repère {Modèle:Mvar} utilisé (la trace est un invariant du tenseur).
La pression isostatique Modèle:Mvar est la pression hydrostatique qui, dans le cas d'un matériau isotrope, produirait la même variation de volume (par rapport à l'état de repos) que le tenseur des contraintes Modèle:Math :
- avec la convention usuelle de la mécanique des milieux continus (contraintes normales positives en extension) :
- <math>p=-\bar\sigma</math> ;
- avec la convention contraire (contraintes normales positives en compression) :
- <math>p=\bar\sigma</math>.
La pression isostatique est une généralisation de la notion de pression hydrostatique dans les liquides. En sciences de la Terre, on parle de pression lithostatique (avec souvent la convention des contraintes normales positives en compression).
Principe de la coupure
Une manière simple de déterminer le tenseur des contraintes consiste à employer le « principe de la coupure ». Il s'agit d'une opération de pensée dans laquelle on scie l'objet selon un plan donné.
Supposons un solide se déformant sous l'effet de deux forces extérieures opposées. Si l'on coupe le solide en deux et que l'on sépare les moitiés, alors chaque moitié n'est soumise qu'à une seule force et donc n'est plus déformée mais mise en mouvement. Pour que chaque moitié retrouve sa déformation, il faut exercer une pression sur chacune des faces de la coupure.
Lorsqu'il y a des symétries évidentes à un problème, le choix de plans de coupe judicieux permet de déterminer de manière simple le tenseur des contraintes. C'est ainsi que l'on peut déterminer que dans le cas de la torsion d'un tube, on a un cisaillement pur.
Calcul des vecteurs-contrainte
Considérons le petit élément de volume <math>d\tau</math> délimité par le tétraèdre de sommets Modèle:Math, Modèle:Math, Modèle:Math, Modèle:Math. Les vecteurs normaux aux faces sont donc <math>\vec e_1,\vec e_2,\vec e_3</math> et le vecteur de composantes <math>(1/\mathrm{d}x_1, 1/\mathrm{d}x_2, 1/\mathrm{d}x_3)</math>. La force <math>\vec{\mathrm{F}}</math> s'exerçant sur une face vérifie
- <math>\vec \mathrm{F} = \mathrm{T} \cdot \vec n </math>
où <math>\vec n</math> le vecteur caractéristique de la face, c'est-à-dire le vecteur normal ayant pour norme l'aire de la face.
On a par exemple sur la face Modèle:Math, la relation
- <math>\vec \mathrm{F} = \begin{pmatrix} \mathrm{F}_1 \\ \mathrm{F}_2 \\ \mathrm{F}_3 \end{pmatrix}
= \begin{pmatrix} \sigma_{11} & \sigma_{12} & \sigma_{13}\\ \sigma_{12} & \sigma_{22} & \sigma_{23}\\ \sigma_{13} & \sigma_{23} & \sigma_{33}\\ \end{pmatrix} \cdot \begin{pmatrix} 0\\ 0\\ \frac12\,\mathrm{d}x_1\,\mathrm{d}x_2\\\end{pmatrix}</math>
Contraintes principales
Modèle:Article détaillé Puisque le tenseur des contraintes est une matrice symétrique réelle, il existe (au moins) une base orthonormée de l'espace dans laquelle le tenseur des contraintes est une matrice diagonale :
- <math>\mathrm{T}_{ij} = \begin{pmatrix}
\sigma_1 & 0 & 0 \\ 0 & \sigma_2 & 0 \\ 0 & 0 & \sigma_3 \\ \end{pmatrix}</math>.
Les valeurs Modèle:Math, Modèle:Math et Modèle:Math sont appelées contraintes principales (il ne faut pas les confondre avec les contraintes en notation de Voigt).
On rencontre également souvent la notation :
- <math>\mathrm\overset{=}{\sigma} = \begin{pmatrix}
\sigma_\mathrm I & 0 & 0 \\ 0 & \sigma_\mathrm{II} & 0 \\ 0 & 0 & \sigma_\mathrm{III} \\ \end{pmatrix}_{(\overrightarrow\mathrm{I},\overrightarrow\mathrm{II},\overrightarrow\mathrm{III})}</math> où Modèle:Math, Modèle:Math et Modèle:Math sont les contraintes principales et <math>(\overrightarrow\mathrm{I},\overrightarrow\mathrm{II},\overrightarrow\mathrm{III})</math> un repère principal, c'est-à-dire un repère dans lequel la représentation matricielle du tenseur est diagonale.
Invariants du tenseur des contraintes
Les contraintes principales <math>\sigma_1</math>, <math>\sigma_2</math> et <math>\sigma_3</math> permettent de déterminer les invariants du tenseur, c'est-à-dire les grandeurs caractéristiques qui ne dépendent pas de la base choisie :
- <math>\mathrm{I}_1 = \mathrm{tr}(\mathrm{T}_{ij}) = \sigma_1 + \sigma_{2} + \sigma_{3}\;</math> ;
- <math>\mathrm{I}_2 = \mathrm{tr}[\mathrm{com}(\mathrm{T}_{ij})]
=\frac12[\mathrm{tr}(\mathrm{T}_{ij})^2-\mathrm{tr}(\mathrm{T}_{ij}\mathrm{T}_{ij}^T)] = \sigma_1 \sigma_2 + \sigma_2 \sigma_3 + \sigma_3 \sigma_1</math> ;
- <math>\mathrm{I}_3 = \mathrm{det}(\mathrm{T}_{ij}) = \sigma_1 \sigma_2 \sigma_3\;</math>.
Notons que l'on a :
- <math>p = \frac{\mathrm{I}_1}{3}</math>
Déviateur
Le tenseur des contraintes peut se décomposer en une somme de deux tenseurs : le déviateur Modèle:Math et la pression isostatique Modèle:Math :
- <math>\mathrm{T}_{ij} = \begin{pmatrix}
s_{11} & s_{12} & s_{13} \\ s_{12} & s_{22} & s_{23} \\ s_{13} & s_{23} & s_{33} \\ \end{pmatrix} + \begin{pmatrix} p & 0 & 0 \\ 0 & p & 0 \\ 0 & 0 & p \\ \end{pmatrix} </math>. Cette décomposition simplifie l'expression des énergies de déformation élastique, de changement de volume et de distorsion.
Le déviateur a les mêmes directions principales que le tenseur des contraintes, on a alors dans ce repère :
- <math>s_{ij} = \begin{pmatrix}
s_{1} & 0 & 0 \\ 0 & s_{2} & 0 \\ 0 & 0 & s_{3} \\ \end{pmatrix} </math>. On peut définir des invariants pour le déviateur :
- <math>\mathrm{J}_1 = s_{1} + s_{2} + s_{3} = \mathrm{tr}(s_{ij}) = 0</math> ;
- <math>\mathrm{J}_2 = - s_1 s_2 - s_2 s_3 - s_1 s_3</math> ;
- <math>\mathrm{J}_3 = s_1 s_2 s_3 = \mathrm{det}(s_{ij})</math>.
Notons que l'on a aussi :
- <math>\mathrm{J}_1 = s_{11} + s_{22} + s_{33}</math> ;
- <math>\begin{align}
\mathrm{J}_2 & = \tfrac{1}{2} \sum_{i, j} s_{ij} s_{ji} \\ &= \tfrac{1}{2}(s_1^2 + s_2^2 + s_3^2) \\ & = \tfrac{1}{6} \left (
(\sigma_{11} - \sigma_{22})^2 + (\sigma_{22} - \sigma_{33})^2 + (\sigma_{33} - \sigma_{11})^2
+6(\sigma_{12}^2 + \sigma_{23}^2 + \sigma_{13}^2)
\right ) \\ & = \tfrac{1}{6} \left (
(\sigma_{1} - \sigma_{2})^2 + (\sigma_{2} - \sigma_{3})^2 + (\sigma_{1} - \sigma_{3})^2
\right ) \\ & = \tfrac{1}{3} \mathrm{I}_1^2 - \mathrm{I}_2 \\ \end{align} </math> ;
- <math>\begin{align}
\mathrm{J}_3 & = \tfrac{1}{3} \sum_{i, j, k} s_{ij} s_{jk} s_{ki} \\ &= \tfrac{1}{3}(s_1^3 + s_2^3 + s_3^3) \\ &= s_1s_2s_3 \\ & = \tfrac{2}{27} \mathrm{I}_1^3 - \tfrac{1}{3} \mathrm{I}_1 \mathrm{I}_2 + \mathrm{I}_3 \end{align} </math>. Ces invariants sont utiles pour définir le domaine élastique et la contrainte de comparaison.
Contrainte équivalente
Il est souvent utile de synthétiser le tenseur des contraintes par un scalaire, appelé « contrainte équivalente », Modèle:Mvar :
- <math>\sigma_e = f(\sigma_{ij})</math>
On peut ensuite comparer cette valeur avec les paramètres du matériau établis par l'essai de traction uniaxial, et notamment la limite d'élasticité Re et la résistance à la traction Rm. Les deux contraintes équivalentes les plus utilisées sont celles de Modèle:Pc et de Modèle:Pc.
Cette approche est toutefois insuffisante dans les cas d'endommagement complexe, comme la fatigue pour des contraintes triaxiales.
Triaxialité des contraintes
La triaxialité des contraintes, notée Modèle:Mvar (lettre grecque « êta »), est le rapport entre la contrainte isostatique et la contrainte équivalente de von Mises :
- <math> \eta = \frac{\frac{1}{3} \mathrm{Tr}(\overline{\overline{\sigma}})}{\sigma_{eq}^{VM}} </math>
Ce paramètre est important dans l'étude de l'endommagement et de la mécanique de la rupture. Dans les cas simples, on a :
- cisaillement pur : <math> \eta = 0 </math>
- traction uniaxiale : <math> \eta = 0{,}33 </math>
- équi-bitraction : <math> \eta = 0{,}66 </math>
- compression uniaxiale : <math> \eta = -0{,}33 </math>
- pression hydrostatique : <math> \eta \rightarrow - \infty </math>
Linéarisation des contraintes
La linéarisation des contraintes est une méthode utilisée pour l'étude des enveloppes minces et en particulier des appareils sous pression. Elle consiste à décomposer le tenseur des contraintes en une contrainte de membrane et une contrainte de flexion.
Pour cela, on définit une ligne perpendiculaire à la surface libre, désignée ci-dessous comme l'axe x, et dont la longueur L est l'épaisseur de la paroi :
- contrainte de membrane : <math>\sigma_\mathrm{m} = \frac{1}{\mathrm{L}} \int_0^\mathrm{L} \sigma(x) \cdot \mathrm{d}x</math>
- contrainte de flexion : <math>\sigma_\mathrm{f} = \frac{12}{\mathrm{L}^3} \int_0^\mathrm{L} \sigma(x) \cdot x \cdot \mathrm{d}x</math>
où Modèle:Mvar est une composante du tenseur des contraintes. On linéarise ainsi chaque composante du tenseur.
Voir aussi
Articles connexes
- Tenseur
- Tenseur des constantes élastiques (ou des rigidités)
- Tenseur des souplesses
- Jauge de déformation
- Cercle de Mohr (Christian Otto Mohr)
- Mesures de tenseur de contraintes par diffraction de neutron
Liens externes
- Modèle:Lien brisé
- {{#invoke:Langue|indicationDeLangue}} Modèle:Pdf G. Backus, Continuum Mechanics, Samizdat Press
- {{#invoke:Langue|indicationDeLangue}} Modèle:Pdf B. Kennett, Lecture Notes for a Course on Continuum Mechanics, Samizdat Press
