Voir la source de Mécanique des milieux continus/Description de l’évolution du milieu continu

{{Chapitre
  | idfaculté = sciences de l'ingénieur
  | numéro    = 1
  | précédent = [[../|Sommaire]]
  | suivant   = [[../Grandeurs admettant une densité par rapport au volume/]]
  | niveau    = 14
}}

== Choix des coordonnées  ==
<math>\epsilon</math> : Espace affine (espace des points)<br />
E : Espace vectoriel associé (espace des vecteurs)<br />
[[Fichier:Milieu continu 1.jpg|center|Milieu continu à deux instant t0 et t1]]
<br />
On définit le domaine <math>D_0</math> occupé par le milieu continu à l’instant 0 avec une répartition continue de matière à l’intérieur.<br />
Bord de <math>D_0</math> : surface <math>\Sigma_0</math><br />
Bort de <math>D_t</math> : surface <math>\Sigma_t</math><br />
Dans l’espace des points <math>\epsilon</math>, on peut choisir une origine O. Les points <math>M_0</math> de <math>D_0</math> sont déterminés par les composantes <math>X_1</math>, <math>X_2</math> et <math>X_3</math> du vecteur <math>\textstyle{\vec X = \overrightarrow{OM_0}}</math> dans la base <math>\vec{e_1}</math>, <math>\vec{e_2}</math> et <math>\vec{e_3}</math> choisie orthonormée en général.<br />
À l’instant t, le point <math>M_0</math> de <math>D_0</math> est devenu un point M de <math>D_t</math> déterminé par les composantes <math>x_1</math>, <math>x_2</math> et <math>x_3</math> du vecteur <math>\textstyle{\vec x = \overrightarrow{OM}}</math>.<br />
<div style="text-align: center;"><math>\vec X = X_1 \vec{e_1} + X_2 \vec{e_2} + X_3 \vec{e_3} \qquad ; \qquad \vec x = x_1 \vec{e_1} + x_2 \vec{e_2} + x_3 \vec{e_3}</math></div><br />
Connaître le mouvement d’un milieu continu, c’est d’être capable d’exprimer la position du point M en fonction de sa position initiale <math>M_0</math> et du temps t.<br />
<math>x_1=\phi_1(X_1,X_2,X_3,t)</math><br />
<math>x_2=\phi_2(X_1,X_2,X_3,t)\, \qquad ou \qquad \vec x = \vec \phi(\vec X,t)</math><br />
<math>x_3=\phi_3(X_1,X_2,X_3,t)</math><br />
Connaître le mouvement, c’est connaître la fonction <math>\phi</math> qui permet de faire la transformation entre les deux états étudiés (<math>D_0</math> et <math>D_t</math>).<br />
{{Définition
  | titre   = Hypothèses
  | contenu =
* Un point M est issu d’un seul point <math>M_0</math>.
On appelle trajectoire la courbe suivie par le point <math>M_0</math> au cours du temps.<br />
* On suppose que les trajectoires ne se mélangent pas.
* On suppose que la fonction <math>\phi</math> est continue et dérivable au moins 2 fois par rapport au temps. (→ permet d'introduire la vitesse et l'accélération en termes de fonction continu)
* On suppose que la fonction <math>\phi</math> est inversible par rapport à X. À tout instant t, <math>\vec x = \vec \phi(\vec X,t)</math> peut s’inverser en <math>\vec X =\vec \psi(\vec x,t)</math>.
}}
<br />
Connaître le mouvement peut se regarder d’une deuxième façon : c’est d’être capable de dire où se trouvait à l’instant 0 la particule qui se trouve en x à l’instant t.<br />
Les 4 variables (<math>X_1</math>, <math>X_2</math>, <math>X_3</math>, t) sont dites variables de Lagrange. C’est le choix de variable que l’on fait quand on exprime le mouvement de <math>\vec x =\vec \phi(\vec X,t)</math>.<br />
Les 4 variables (<math>x_1</math>, <math>x_2</math>, <math>x_3</math>, t) sont dites variables d'Euler. C’est le choix de variable que l’ont fait quand on veut exprimer le mouvement de <math>\vec X =\vec \psi(\vec x,t)</math><br />.
<br />
<u>Exemple :</u> La vitesse d’un point M peut être exprimée en fonction de <math>M_0</math> et de t (variables de Lagrange) ou en fonction de M et de t (variables d’Euler).<br />
Un solide rigide est un milieu continu particulier. Il garde toujours la même forme. Son champ des vitesses s’exprime facilement en variables d’Euler.<br />
<u>Rappel :</u> <math>\overrightarrow{V(M/T_0)} = \overrightarrow{V(G/T_0)} + \overrightarrow{\Omega(S/T_0)} \wedge \overrightarrow{GM}</math> avec <math>T_0=(O;\vec e_1;\vec e_2;\vec e_3)</math> (Pour un solide indéformable).<br />
<math>\overrightarrow{V(G/T_0)}</math> et <math>\overrightarrow{\Omega(S/T_0)}</math> : ne dépendent que du temps.<br />
<math>\overrightarrow{GM}</math> : Position actuelle à l’instant t.<br />
Notation : 
{|class="wikitable"
|-
| <math>y_1(t)</math>
| <math>y_2(t)</math>
| <math>y_3(t)</math>
| coordonnée du centre d’inertie
|-
| <math>\omega_1(t)</math>
| <math>\omega_2(t)</math>
| <math>\omega_3(t)</math>
| composantes du vecteur rotation instantané
|-
| <math>\frac{dy_1(t)}{dt}</math>
| <math>\frac{dy_2(t)}{dt}</math>
| <math>\frac{dy_3(t)}{dt}</math>
| composantes de la vitesse de G
|}
<br />
La vitesse du point M a pour composantes :<br />
<div style="text-align: center;"><math>  \begin{matrix} 
V_1(x_1,x_2,x_3,t) = \frac{dy_1(t)}{dt} + \omega_2(t) \left[ x_3-y_3(t) \right] - \omega_3(t) \left[ x_2-y_2(t) \right]\\
V_2(x_1,x_2,x_3,t) = \frac{dy_2(t)}{dt} + \omega_3(t) \left[ x_1-y_1(t) \right] - \omega_1(t) \left[ x_3-y_3(t) \right]\\
V_3(x_1,x_2,x_3,t) = \frac{dy_3(t)}{dt} + \omega_1(t) \left[ x_2-y_2(t) \right] - \omega_2(t) \left[ x_1-y_1(t) \right]
\end{matrix}
</math></div>
<br />
Ou encore <div style="text-align: center;"><math>\vec V(\vec x,t) = \frac{d\vec y(t)}{dt} + \vec \omega \wedge \left[\vec x - \vec y(t) \right]</math></div><br />
Comme on prend toujours le même repère <math>T_0</math>, on ne le précise pas mais on précise bien les 4 variables <math>x_1</math>, <math>x_2</math>, <math>x_3</math> et t.

== Vitesse ==
=== Définition ===
La vitesse d’un point, c’est la dérivée de sa position par rapport au temps.

=== Expression de la vitesse en variable Lagrange  ===
En dérivant la position <math>\vec x =\vec \phi(\vec X,t)</math> par rapport au temps on obtient un vecteur <math>\frac{\partial \vec \phi(\vec X,t)}{\partial t}</math> dépendant de <math>\vec X</math> et de t.<br />
Ce n’est pas très pratique d’exprimer la vitesse en fonction de la position initiale et du temps. En effet quand on souhaite connaître une vitesse, c’est celle d’un point qu’on observe <math>\vec x</math> et non d’un point dont sa position initiale était <math>\vec X</math>.

=== Expression de la vitesse en variable d’Euler  ===
Le vecteur <math>\frac{\partial \vec \phi(\vec X,t)}{\partial t}</math> peut être ré exprimé en varaiable d’Euler puisque <math>\vec X =\vec \psi(\vec x,t)</math>. En faisant le remplacement de <math>\vec X</math> par <math>\vec \psi(\vec x,t)</math> on obtient :<br />
<div style="text-align: center;"><math>\vec V(\vec X,t) = \frac{\partial \vec \phi(\vec X,t)}{\partial t} \qquad -> \qquad \vec V(\vec x,t)= \frac {\partial \vec \phi(\vec \psi(\vec x,t),t)}{\partial t}</math></div>

=== Méthode pratique pour exprimer la vitesse en variable d’Euler  ===

# Dérivée partiellement par rapport au temps la position <math>\vec x =\vec \phi(\vec X,t)</math>. On obtient un vecteur <math>\frac{\partial \vec \phi(\vec X,t)}{\partial t}</math>.
# Inverser <math>\vec x =\vec \phi(\vec X,t)</math> en <math>\vec X =\vec \psi(\vec x,t)</math>.
# Remplacer <math>\vec X</math> par <math>\vec \psi(\vec x,t)</math> dans <math>\frac{\partial \vec \phi(\vec X,t)}{\partial t}</math>.
<br />
{{Définition
  | titre   = Exemple - Dilatation d'une boule
  | contenu =
[[Fichier:Chap2 vitesse 1.JPG|right|Représentation d'une boule qui se dilate]]
Facteur de dilatation <math>\lambda (t)</math> supposé dérivable et positif.<br />
Le vecteur <math>\vec x</math> est définit par <math>\vec x(\vec X,t)=\lambda (t) \vec X</math>.<br />
À l'instant t=0, <math>\vec x</math> est en <math>\vec X</math> (position initiale) ce qui implique que <math>\lambda (0)=1</math>.<br />
# <math>\frac{\partial \vec \phi(\vec X,t)}{\partial t} = \frac {d \lambda (t)}{dt} \vec X</math>
# <math>\vec x(\vec X,t)= \vec \phi(\vec X,t) =\lambda (t) \vec X \quad \to \quad \vec X = \frac {1}{\lambda (t)} \vec x </math>
# <math>\frac{\partial \vec \phi(\vec X,t)}{\partial t} = \frac {d \lambda (t)}{dt} \vec X = \frac {d \lambda (t)}{dt} \frac {1}{\lambda (t)} \vec x </math>
<br />
La vitesse est donc égale à : <math>\vec V(\vec x,t)= \frac {d \lambda (t)}{dt} \frac {1}{\lambda (t)} \vec x </math> (vitesse colinéaire à <math>\vec x \to</math> déformation radiale)
}}

== Accélération ==
=== Définition ===
L’accélération, c’est la dérivée seconde de la position par rapport au temps.

=== Calcul des composantes de l’accélération  ===
On recherche les composantes <math>\gamma^1</math>, <math>\gamma^2</math> et <math>\gamma^3</math> du vecteur accélération <math>\vec \gamma</math> en dérivant par rapport au temps les composantes de la vitesse.<br />
Par exemple <math>\gamma^1(x^1,x^2,x^3,t) </math> s’obtient en dérivant <math>V^1(x^1,x^2,x^3,t) </math> par rapport au temps en n’oubliant pas que <math>x^1</math> ,<math>x^2</math> et <math>x^3</math> dépendent du temps (<math>\vec x=\vec \phi(\vec X,t)</math>). On dit que x dépend implicitement du temps si on dérive <math>V^1</math> par rapport au temps.

== Dérivée particulaire d’une fonction  ==
=== Définition ===
Soit <math>h(x^1,x^2,x^3,t) </math> une fonction de la position actuelle <math>\vec x</math> et du temps t, on appelle dérivée particulaire de la fonction h la fonction obtenue en dérivant par rapport au temps l’expression de la fonction h sans oublier que la position x dépend du temps (<math>\vec x=\vec \phi(\vec X,t)</math> ; on dit que <math>\vec x</math> dépend du temps de façon implicite).

<u>Remarque :</u> La dérivée partielle <math>\frac{\partial h}{\partial t}</math> ne tient compte que de la dépendance explicite par rapport au temps. Cela signifie qu’elle ne prendrait pas compte que <math>\vec x</math> dépend implicitement du temps t.

<u>Notation :</u> La dérivée particulaire sera notée <math>\frac{dh}{dt}</math>. C’est la dérivée partielle + la contribution de la dépendance implicite de <math>x^1</math>,<math>x^2</math> et <math>x^3</math> par rapport au temps.

=== Calcul de la dérivée particulaire d’une fonction  ===
D’après le théorème de la dérivée d’une fonction composée, on peut écrire :

<div style="text-align: center;"><math>\underbrace{\frac {dh}{dt}}_{\text{derivee particulaire}} = \underbrace{\frac{\partial h}{\partial t}}_{\text{derivee partielle}}+\underbrace{\frac{\partial h}{\partial x^1} \frac{\partial x^1}{\partial t} + \frac{\partial h}{\partial x^2} \frac{\partial x^2}{\partial t} + \frac{\partial h}{\partial x^3} \frac{\partial x^3}{\partial t}}_{\text{dependance implicite de}\ x\ \text{par rapport au temps}}</math></div>

== Dérivée particulaire d’un vecteur  ==
=== Définition ===
=== Composantes de <math>\frac{d\vec w}{dt}</math>  ===

La différenciation du vecteur <math>\vec w</math> se fait, selon les coordonnées choisies :

<div style="text-align: center;"><math>d\vec w = \frac{\partial \vec w}{\partial t} dt + \frac{\partial \vec w}{\partial x_1} dx_1 + \frac{\partial \vec w}{\partial x_2} dx_2 + \frac{\partial \vec w}{\partial x_3} dx_3</math></div>

En divisant par <math>dt</math>, on trouve :

<div style="text-align: center;"><math>\frac{d\vec w}{dt} = \frac{\partial \vec w}{\partial t} + \frac{dx_1}{dt}\frac{\partial \vec w}{\partial x_1} + \frac{dx_2}{dt} \frac{\partial \vec w}{\partial x_2}+ \frac{dx_3}{dt} \frac{\partial \vec w}{\partial x_3}</math></div>

Soit :

<div style="text-align: center;"><math>\frac{d\vec w}{dt} = \frac{\partial \vec w}{\partial t} + (\vec v \cdot \overrightarrow{\textrm{grad}}) \vec w</math></div>

=== Première expression de l’accélération  ===
=== Deuxième expression de l’accélération, une expression complètement vectorielle ===

== Ligne de courant et trajectoires  ==
Un mouvement de milieu continu est donné.

Si l’on fixe une particule , la trajectoire au sens de la cinématique est exactement la trajectoire vue comme trajectoire du champ eulérien des vitesses au sens des systèmes dynamiques, le paramètre étant ici le temps .

Si l’on fixe un instant , la position de à cet instant est . La trajectoire de au sens des systèmes dynamiques, pour le champ eulérien des vitesses à l’instant , est une courbe qui peut être sans rapport avec le mouvement de

C’est ainsi qu’apparait la notion de ligne de courant (on devrait préciser : à l’instant ). Pour ces courbes, étant fixé, le paramètre n’a aucun rapport avec le temps.
Voici un exemple (le détail des calculs n’a pas d’importance pour l’instant).

 est un carré du plan, le mouvement est défini pour  par

De la relation on déduit le champ des vitesses

La figure suivante représente les pour à avec un pas . Pour , on obtient le carré initial (en bas), et l’on suit le mouvement d’un réseau de 36 points, où sont représentés les vecteurs . 
Les lignes de courant à un instant fixé sont les courbes d’équation

La figure suivante représente des lignes de courant à l’instant . 
Le mouvement n’est pas stationnaire, les trajectoires sont distinctes des lignes de courant. Elles sont rectilignes.

 Dans le cas où le champ des vitesses à  fixé n’admet pas de solution exacte facilement exprimable, ou si le champ n’est pas défini analytiquement mais en des points donnés d’un maillage numérique ou expérimental, il faut procéder à une résolution numérique et donc commencer par discrétiser les équations... Dans ce cas, le résultat n’est généralement pas une expression analytique ou paramétrique des lignes de courant, mais juste leur tracé comme sur la figure ci-dessous qui représente quelques lignes de courant (les courbes en rouge) du champ de vitesse à un instant  d’un jet plan en impact sur une plaque perpendiculaire au jet (les vecteurs vitesses - un peu flous - sont colorées en fonction de leur norme).

 
=== Définition de la trajectoire d’une particule qui était en <math>\vec X</math> à l’instant t=0 ===
=== Ligne de courant à l’instant t fixé passant par un point <math>\vec x_0</math>  ===

== Notation pour les dérivées  ==

{{Bas de page
  | idfaculté = sciences de l'ingénieur
  | précédent = [[../|Sommaire]]
  | suivant   = [[../Grandeurs admettant une densité par rapport au volume/]]
}}