Introduction à l'élasticité/Relation fondamentale de la dynamique

Introduction

Après avoir traduit les efforts intérieurs par le tenseur des contraintes, nous pouvons revenir sur l’expression de la relation fondamentale de la dynamique et exprimer l'accélération d'un point à partir des contraintes qu’il subit.

Dynamique locale

On s'intéresse, dans cette partie, à un petit volume v(t) que l’on suit dans son déplacement.

Bilan de quantité de mouvement

Nous avons montré que le bilan de quantité de mouvement sur un volume $v(t)$ prenait la forme :

$\underset{v\left(t\right)}{\int }(\rho 𝐚-{𝐟}_v){\mathrm{d}}^{3}V=\underset{\partial v\left(t\right)}{\int }\mathit{\sigma }\left(𝐧\right){\mathrm{d}}^{2}S$

Effectuons le produit scalaire de chaque membre avec un même vecteur u quelconque :

$\underset{v\left(t\right)}{\int }(\rho 𝐚-{𝐟}_v)\cdot 𝐮{\mathrm{d}}^{3}V=\underset{\partial v\left(t\right)}{\int }\mathit{\sigma }\left(𝐧\right)\cdot 𝐮{\mathrm{d}}^{2}S$

Remarquons que le membre de droite peut s'écrire ainsi, en transposant σ :

$\begin{array}{cc}\underset{\partial v\left(t\right)}{\int }\mathit{\sigma }\left(𝐧\right)\cdot 𝐮{\mathrm{d}}^{2}S& =\underset{\partial v\left(t\right)}{\int }𝐧\cdot {\mathit{\sigma }}^{\mathrm{T}}\left(𝐮\right){\mathrm{d}}^{2}S\\ & =\underset{v\left(t\right)}{\int }{\mathrm{d}\mathrm{i}\mathrm{v}}_{𝐱}\left({\mathit{\sigma }}^{\mathrm{T}}\left(𝐮\right)\right){\mathrm{d}}^{3}V\end{array}$

d'après le théorème de flux-divergence (dit également de Stokes ou de Green-Stokes), où div_x est l'opérateur divergence (vectoriel). Nous allons exprimer cette dernière intégrale de manière plus élégante en introduisant la notion de divergence d'un tenseur.

Divergence du tenseur des contraintes

Remarquons que la divergence est une application linéaire : nous allons lui associer un tenseur.

Expression de la relation fondamentale de la dynamique

On a, pour tout volume v(t) et tout vecteur u, la relation suivante :

$\underset{v\left(t\right)}{\int }(\rho 𝐚-{𝐟}_v)\cdot 𝐮{\mathrm{d}}^{3}V=\underset{v\left(t\right)}{\int }(𝐃𝐢{𝐯}_{𝐱}\mathit{\sigma }\cdot 𝐮){\mathrm{d}}^{3}V$

Puisque l’on peut prendre n’importe que volume et n’importe quel vecteur, il en résulte l'égalité suivante :

On peut remarquer qu’il s'agit de trois équations reliant l'accélération au tenseur des contraintes — mais que ce dernier étant de dimension 6, cette relation ne suffit pas à le déterminer. Il faudra donc ajouter des conditions aux limites.

En pratique, on peut être amené à supposer la forme de σ pour effectuer les calculs, quitte à infirmer l'hypothèse faite si l'expérience est en désaccord trop marqué avec le modèle.

Symétrie du tenseur des contraintes

Il est important de remarquer qu'une conséquence de la relation fondamentale de la dynamique est qu'elle nous renseigne sur la symétrie du tenseur des contraintes.

Ce résultat peut se démontrer de plusieurs manières, plus ou moins intuitives. Une démonstration est proposée en annexe et repose sur l’expression du transport du moment dynamique donnée dans un précédent chapitre.

Dynamique matérielle

Lorsqu'un matériau subit de grand déplacements, on ne connaît pas en pratique la position de tous les volumes internes, susceptibles de beaucoup se déformer. Nous allons en fait nous ramener à un domaine fixe en exploitant le fait que la cinématique garde la trace des déformations subies par la pièce.

Cela nous permettra de traiter les propriétés d'une pièce dans son entier, basées sur celles, microscopiques, que nous avons détaillé, mais sans avoir à retourner à ce niveau de détail.

Nous verrons que, dans le cas de l'élasticité, nous serons amenés à ne pas tenir compte des déformations. Cependant, cette partie explique pourquoi cette approximation est légitime et, hors du cadre restrictif des solides élastiques, reste tout à fait générale. En première lecture, les deux sous-parties qui suivent peuvent être survolées.

Dynamique dans la configuration initiale

Rappelons que, dans une intégrale, un changement de variables fait intervenir :

• la transformation (changement du domaine d'intégration) ;
• le jacobien de la transformation (noté J dans la suite).

Moyennant quelques calculs, on peut ainsi montrer que l'équation de la résultante dynamique appliquée à un volume v(t) qu'on suit dans son mouvement :

$\underset{v\left(t\right)}{\int }\rho 𝐚{\mathrm{d}}^{3}V=\underset{v\left(t\right)}{\int }{𝐟}_v{\mathrm{d}}^{3}V+\underset{\partial v\left(t\right)}{\int }\mathit{\sigma }\left({𝐧}_t\right){\mathrm{d}}^{2}S$

est équivalente à la relation suivante (où les notations sont un peu abusives, pour rester simples), qui s'applique sur la configuration initiale :

$\underset{v(0)}{\int }\rho 𝐚J{\mathrm{d}}^{3}V=\underset{v(0)}{\int }{𝐟}_{v}J{\mathrm{d}}^{3}V+\underset{\partial v(0)}{\int }\mathit{\sigma }\cdot {𝐅}^{-\mathrm{T}}\left({𝐧}_0\right)J{\mathrm{d}}^{2}S$

On applique, ici encore, le théorème de Green-Stokes :

${\rho }_0{𝐚}_{𝐩}={𝐟}_{v}J+{𝐃𝐢𝐯}_{𝐩}(J\mathit{\sigma }\cdot {𝐅}^{-\mathrm{T}})$

Avec des notations plus rigoureuses, cette relation s'écrit :

${\rho }_0(𝐩){𝐚}_{𝐩}(𝐩,t)={𝐟}_v(f(𝐩,t),t)J(𝐩,t)+{𝐃𝐢𝐯}_{𝐩}(J(𝐩,t)\mathit{\sigma }(f(𝐩,t),t)\cdot {𝐅}^{-\mathrm{T}})$

Dans la suite, nous allons introduire des notations simplificatrices qui nous permettront de mieux comprendre la situation.

Tenseurs des contraintes de Piola-Kirchhoff

Une propriété importante de ce tenseur, c’est qu’il met en relation les efforts sur une facette à tout instant et ceux subit par la facette dans la configuration initiale :

$\mathit{\pi }\left({𝐧}_0\right){\mathrm{d}}^2{S}_0=\mathit{\sigma }\left({𝐧}_t\right){\mathrm{d}}^2{S}_t$

Par sa définition, on voit également que π n'est, dans le cas général, pas symétrique. On préfère alors souvent utiliser le tenseur suivant :

Ce tenseur, comme le précédent, permet de suivre les efforts intérieurs transportés par la déformation, à partir de la configuration initiale. On aboutit ainsi à la relation fondamentale de la dynamique :

On remarque la présence de F, qui est un terme géométrique, dans cette équation. Cela est lié au fait que l'équation obtenue à la section précédente est obtenue en suivant le volume (on parle d'une description spatiale ou lagrangienne) alors que nous avons ici considéré les déformations en conservant un point de vue fixe (on parle d'une description matérielle ou eulérienne).

Pour des structures élancées, ces termes géométriques ont une grande influence sur la dynamique. Jusqu'ici, la description faite de la dynamique est tout à fait générale : dans la section qui suit, nous faisons les hypothèses de l'élasticité et observons comment cela simplifie grandement le problème.

Dynamique locale infinitésimale

Rappelons les hypothèses de l'élasticité et analysons leurs conséquences : les déformations sont petites et réversibles. On considère alors, au premier ordre :

• que x et p se confondent (de même que leurs dérivées) ;
• que la masse volumique ne varie pas : on confond ρ et ρ₀.

De fait, l'équation obtenue en dynamique matérielle se simplifie beaucoup :

On peut raffiner cette équation en tenant compte des effets dits du second ordre. Nous ne le ferons pas dans le cadre de ce cours et n'utiliserons que la forme ci-dessus pour tous les exemples, exercices et développements.

Remarques et références

Modèle:Bas de page