- Equations fondamentales de la mecanique des milieux continus
-
Équations fondamentales de la mécanique des milieux continus
Les équations fondamentales de la mécanique des milieux continus regroupent les énoncés mathématiques les plus couramment utilisés dans le cadre de la mécanique des milieux continus, principalement applicables à la déformation des solides et à la mécanique des fluides.
Quelques quantités physiques intégrales du système
En général, on étudie les déformations et la dynamique de la Terre en termes d'un modèle de continuum caractérisé par une distribution interne de masse non spécifiée pour le moment, autrement dit par une densité[1] massique ρ, dans un volume B délimité par une surface . Dans le problème qui nous occupe, cette frontière n'est pas donnée a priori, mais doit être déterminée en supposant qu'elle correspond à la forme d'équilibre du corps déformable. En outre, nous supposerons pour le moment qu'il existe un champ de vitesse vi par rapport au système de référence choisi mais non encore spécifié, ainsi qu'une densité de moment cinétique interne , une force volumique fi, un champ de contraintes (ou champ de tensions) Tij, un moment de force volumique Mi et moment de force surfacique Mij. Toutes ces quantités peuvent être des fonctions explicites ou implicites du temps t. Sauf indication différente, nous utiliserons un système de référence cartésien dont l'origine (0, 0, 0) se trouve au point O. Nous dénoterons les coordonnées spatiales par xi et les dérivées spatiales successives , , ... par , , ... (avec i,j = 1,2,3). En règle générale, nous adoptons ici la notation indicielle pour les vecteurs et les tenseurs, ainsi que la règle de sommation d'Einstein sur les indices muets répétés.
La masse totale du modèle est
où dτ = dx1dx2dx3 est un élément de volume.
L'impulsion totale est
et le moment cinétique total est
où est le symbole de Levi-Civita.
La résultante des forces appliquée au volume B du corps est
où dσ est un élément de la surface et (n1,n2,n3) dénote un vecteur unitaire normal à la surface et pointant hors du volume B. De la même, manière, le moment de force résultant appliqué au corps B est
. Ces formules sont valables en toute circonstance. Si nous admettons en outre que la surface-frontière de B est régulière et que les champs des tensions et des moments de force surfaciques sont différentiables dans B, on peut appliquer le théorème de la divergence de Gauss et écrire :
.Lois de conservation
Les équations régissant un milieu continu peuvent être déduites, soit en utilisant un « principe de correspondance » pour adapter les lois de la mécanique des points discrets de Newton aux points matériels de la mécanique des milieux continus,[2] soit en exprimant de façon directe le fait que certains attributs physiques, tels que la masse, l'impulsion, le moment cinétique, l'énergie cinétique, l'énergie interne, l'énergie totale, etc., contenus dans un volume arbitraire b(t) se déformant au cours du temps t, ne peuvent pas changer de manière arbitraire, mais que leurs changements précis sont régis par des lois de conservation. Ces dernières peuvent s'obtenir en exprimant le taux de changement dans le temps dans un volume arbitraire de matière au moyen du théorème de transport de Reynolds.
Soit b(t) une partie quelconque se déformant au cours du temps du volume matériel B(t). On suppose que sa frontière est régulière, et on oriente la normale unitaire n(t) en un point quelconque de vers l'extérieur. Pour une quantité physique arbitraire q associée à un point matériel appartenant au volume b(t) ou à sa frontière , on a alors
. Les notations utilisées ici sont les suivantes : est une quantité physique arbitraire, et q la densité volumique du champ correspondant ; les opérateurs Dt et sont les dérivations totales (matérielles) et partielles (locales) par rapport au temps, respectivement ; est la surface régulière fermée qui limite le volume b, de normale extérieure unitaire ni(i = 1,2,3). En utilisant à nouveau le théorème de la divergence de Gauss, il vient
. Si κ est le taux local de production ou de destruction de la quantité Q, c'est-à-dire , alors la formule précédente devient
ou, puisque le volume b(t) est arbitraire et peut donc être pris arbitrairement petit :
. Ceci est la forme eulérienne générale d'une loi de conservation en physique des milieux continus. La forme lagrangienne équivalente, valable pour chaque point matériel individuel, est .
On déduit cette forme lagrangienne aisément de la forme eulérienne à l'aide de l'identité (voir l'article Concepts de base en théorie des milieux continus) .
Conservation de la masse, ou équation de continuité
Sous forme intégrale, la conservation de la masse s'exprime par[3]
. Cette formule fournit l'équation de continuité, exprimant la conservation de la masse sous forme différentielle, en faisant tendre le volume arbitraire b(t) tendre vers un point : .
La forme eulérienne de cette dernière équation est .
On remarque qu'il s'agit bien ici d'une forme particulière de la loi de conservation générale , où q = ρ et κ = 0. Le fait que κ = 0 indique bien qu'il n'y a ni production ni destruction de masse, laquelle est donc conservée.
Conservation de l'impulsion, ou équation de mouvement
De manière similaire, à peine plus compliquée, on peut établir l'équation de conservation de l'impulsion régissant le mouvement d'un point matériel, en posant q = ρvi et en admettant explicitement que la masse se conserve aussi. On obtient
. En effet, l'équation fondamentale de la dynamique exprime le fait que le taux de variation de l'impulsion DtPi est équilibré par la résultante des forces Fi agissant sur le volume. On a donc
. Si nous transformons le membre de gauche au moyen du théorème de Reynolds, nous obtenons
. Dès lors, en faisant tendre le volume arbitraire b vers zéro, on obtient la forme eulérienne générale de l'équation de mouvement, à savoir .
Il convient de noter que dans cette équation, on ne présuppose aucunément qu'il y ait conservation de masse, mais seulement que le tenseur des contraintes Tik est différentiable. En développant le membre de gauche, on obtient
. Le terme entre crochets dans le membre de gauche représente l'accélération Dtvi. Le terme entre accolades dans ce même membre de gauche représente simplement l'équation de conservation de la masse sous forme eulérienne. Ainsi donc, si nous supposons que la masse reste conservée — ce qui est pratiquement toujours le cas en géophysique interne — l'équation de conservation de la quantité de mouvement s'écrit .
Conservation du moment cinétique, ou équation de rotation
L'équation exprimant la conservation du moment cinétique, qui est primordiale dans les études traitant de la rotation d'un corps déformable comme la Terre, peut être obtenue de manière semblable à celle exprimant la conservation de l'impulsion. En supposant qu'il y ait conservation de la masse et conservation de la quantité de mouvement, on obtient
. En effet, le taux de variation du moment cinétique DtLi dans un volume arbitraire b est équilibré par le moment de force résultant Mi agissant sur ce volume. Nous avons donc
. Comme précédemment, on transforme le membre de gauche en appliquant le théorème de Reynolds pour obtenir
. En faisant tendre le volume arbitraire b vers zéro, on obtient l'équation exprimant localement la conservation du moment cinétique, à savoir
. Cette équation peut être réarrangée comme suit :
. Dans le membre de gauche, les termes entre chevrons représentent la quantité , les termes entre crochets s'annulent si la masse est conservée, et les termes entre accolades s'annulent si à la fois la masse et l'impulsion sont conservées. Dans le membre à droite du signe =, le terme est identiquement nul, puisque vjvk est un tenseur symétrique. Admettant que masse et impulsion soient conservées, on obtient donc .
Symétrie du tenseur des contraintes
Si, en outre, nous supposons que le moment cinétique intrinsèque d'un point matériel quelconque ne change pas au cours du mouvement de déformation, autrement dit si , alors le tenseur des contraintes Tjk est symétrique : .
Conservation de l'énergie
L'équation exprimant la conservation de l'énergie peut s'obtenir par la même méthode que celle employée pour la masse, la quantité de mouvement et la quantité de rotation. On trouve :
. Le symbole U désigne ici l'énergie interne par unité de volume, Hi est le vecteur densité de flux de chaleur, et Γ est le taux de production ou de destruction intrinsèque d'énergie par unité de volume, par réactions chimiques, par radioactivité, ou par tout autre processus interne. Le symbole v2 représente bien sûr le carré vkvk de la vitesse, c'est-à-dire deux fois l'énergie cinétique par unité de masse.
Tableau récapitulatif des principales lois de conservation
Le tableau ci-dessous résume les principales lois de conservation générales, sous forme eulérienne, utiles pour résoudre les problèmes posés par la mécanique des milieux continus, en particulier lorsqu'il s'agit d'étudier les déformations globales de la Terre suite à un séisme, suite à une variation de la longueur du jour, suite à l'application d'un potentiel de marée ou suite à l'application d'une distribution de masse sur la surface externe. Les notations employées dans le tableau sont celles de ce texte. Rappelons que l'intégrale est la quantité physique qui se conserve ; elle représente, selon les cas, la masse totale du corps B, sa quantité de mouvement totale Pi, son moment cinétique total Li, ou son énergie totale E. La forme eulérienne générale d'une loi de conservation est .
Principales lois de conservation de la mécanique des milieux continus Q φ kj q M ρ ρvj 0 Pi ρvi ρvivj − Tij fi Li E vkfk + Γ Hypothèses et lois pour décrire les déformations de la Terre
Article connexe : Lois de déformation.Dans les problèmes de géodynamique théorique, on considère en général que la résultante des forces externes par unité de volume fi se compose d'une force gravifique volumique causée par l'autogravitation et par l'attraction gravifique de corps extérieurs, d'une force électromagnétique volumique , de forces volumiques inertielles (Coriolis et axifuge) , de forces de marée , et de forces volumiques causées par des charges extérieures . D'autre part, les forces surfaciques Ti(x,t,n) dépendant de l'orientation n qui contribuent en général au tenseur des contraintes Tij(x,t) via les relations de Cauchy Ti(x,t,n) = Tij(x,t)nj(x,t,n) sont les forces élastiques qui engendrent un tenseur des contraintes élastiques , les forces superficielles électromagnétiques qui engendrent un tenseur de Maxwell , des forces visqueuses qui engendrent un tenseur de Reynolds un tenseur des contraintes élastiques , des forces de frottement de contact non visqueux qui engendrent un tenseur des forces de friction un tenseur des contraintes élastiques , et de toute autre force surfacique engendrant un tenseur des contraintes un tenseur des contraintes élastiques .
La force de gravité par unité de masse g, plus simplement appelée accélération gravifique ou gravité, dérive d'un potentiel gravifique φg que nous assimilerons ici à l'énergie potentielle gravifique[4] :
. Le potentiel φg est déterminé par l'équation de Poisson
, où G est la constante de gravitation. Il convient de noter que les propriétés du potentiel gravifique impliquent que φg et g sont des fonctions continues partout dans l'espace, et en particulier sur les frontières entre deux milieux continus différents.
En ce qui concerne le champ de force produit par des interactions électromagnétiques, on admet d'habitude que les points matériels ne transportent pas de charge électrique, c'est-à-dire qu'ils sont électriquement neutres. De la sorte, on suppose qu'aucune force de Lorentz ne s'applique directement à un point massique ξ se mouvant dans le champ magnétique de la Terre. Cependant, le matériau terrestre conduit l'électricité à des degrés divers, autrement dit possède une conductivité électrique σe finie non nulle. Dès lors, lorsque ce matériau se déplace à une vitesse v non nulle, il interagit avec le champ géomagnétique interne et engendre une densité de courant électrique induit J qui modifie le champ magnétique initial en donnant naissance à un champ d'induction magnétique variable B(ξ,t) tandis que le flux variable de ce courant produit un champ d'induction électrique E(ξ,t). Le résultat de toutes ces interactions du matériau terrestre en mouvement avec le champ géomagnétique B est la naissance d'un champ de force électromagnétique induit fem = c − 1J.B, avec J = σe(E + c − 1v.B), lequel modifie à son tour le mouvement. On a utilisé ici des unités électrostatiques de Gauss, c étant la vitesse de la lumière dans le vide. Il convient de remarquer que dans un diélectrique parfait (σe = 0), la force électromagnétique serait nulle (fem = 0). Sur une interface électromagnétique de normale extérieure n, les quantités n.E,n.H et n.B sont continues. Lorsqu'on étudie les déformations élastiques de la Terre, notamment la propagation des ondes sismiques, les oscillations libres, les déformations de marée, etc., on néglige le plus souvent de considérer de tels phénomènes magnétohydrodynamiques ou magnétoélastiques, en admettant simplement que les forces de rappel élastiques sont beaucoup plus importantes que les forces électromagnétiques. Il faut cependant garder à l'esprit que ces forces existent et sont essentielles lorsqu'on considère la génération du champ géomagnétique interne dans le noyau externe de la Terre par effet de dynamo autoexcitée. Elles peuvent éventuellement devenir importantes lorsqu'on étudie les déformations et oscillations du noyau, mais sont négligeables en pratique lorsqu'on considère des déformations globales de la Terre de faible amplitude.
Les forces d'inertie peuvent se formuler ainsi :
où Ω est le vecteur rotation instantanée, r le vecteur position \mathrm du point matériel \xi</math>, et v la vitesse instantanée de ce point. Le premier terme représente la force axifuge et le second terme représente la force de Coriolis par unité de volume.
Les déformations dues aux marées sont produites par des forces dérivant aussi d'un potentiel, le potentiel de marée Wm. Dans le cas des marées terrestres, compte tenu de la sensibilité des appareils de détection actuellement sur le marché, il est nécessaire de considérer les termes quadrupolaires (n = 2) et termes octupolaires (n = 3) pour la Lune et le Soleil dans le développement général du potentiel de marée en multipôles :
. Ici, G dénote la constante de gravitation (G = 6,673.10 − 11m3s − 2kg − 1), MA est la masse de l'astre générateur de la marée, dA est la distance du centre de masse de la Terre au centre de masse de A, zA est la distance zénithale de A et r est la distance du centre de masse de la Terre au point où l'on mesure la marée. Le symbole Pn désigne, comme d'habitude, le polynôme de Legendre de degré n. La première sommation dans l'expression ci-dessus s'étend en théorie sur tous les astres A pouvant engendrer un effet de marée sur la Terre ; en pratique, elle se limite à la Lune, au Soleil et, dans des calculs très précis, à Vénus et Jupiter.
Loi de Hooke
Dans de nombreuses applications géophysiques, particulièrement en sismologie, la Terre est traitée comme un matériau linéairement élastique et isotrope, ce qui signifie que le tenseur lagrangien des tensions élastiques est exprimé par la loi de Hooke. Dans un référentiel cartésien, on a donc
où est le champ des déplacements et où les constantes de matériau κ et μ sont les modules de compression et de cisaillement, respectivement. Dans les problèmes qui concernent des durées s'échelonnant en gros entre une fraction de seconde et quelques heures, c'est-à-dire des laps de temps typiques des périodes des ondes sismiques, des oscillations libres et des marées habituelles, le module de compression à considérer est le module de compression adiabatique, évalué pour des conditions isentropiques. Pour le module de cisaillement, Léon Brillouin a montré en 1940 qu'il n'y avait pas lieu de distinguer entre un cisaillement effectué sous des conditions isentropiques ou des conditions isothermes.[5] D'autre part, on peut tenir compte d'une faible anélasticité en faisant dépendre les coefficients κ et μ de la fréquence ; c'est l'approche adoptée notamment par Kanamori et Anderson[6] et par Dziewonski et Anderson[7].
Conditions aux frontières
Afin de pouvoir étudier théoriquement les déformations globales de la Terre (ou d'une planète tellurique), il est nécessaire d'adjoindre aux équations différentielles qui décrivent le mouvement de déformation les conditions qui s'appliquent au centre et sur la surface externe de la Terre, ainsi que des conditions adéquates qui s'appliquent aux diverses frontières internes entre des milieux continus différents. Ainsi, à côté de la continuité du potentiel gravifique et de la force de gravité, il faut tenir compte de deux types généraux de conditions à appliquer lorsqu'on rencontre une interface simple[8] : des conditions cinématiques et des conditions dynamiques.
Conditions cinématiques
Les conditions cinématiques s'obtiennent aisément en exprimant mathématiquement le fait que les points matériels appartenant à une surface-frontière F[x(ξ,t),t] = cste à un instant donné doivent encore appartenir à la même surface à un autre instant, c'est-à-dire σFmilieu 1 = σFmilieu 2 ou encore, sous forme eulérienne :
. Comme nous nous limitons à des surfaces régulières, la fonction F admet des dérivées partielles continues par rapport au temps et par rapport aux coordonnées spatiales. Il s'ensuit que la composante normale du champ de déplacement relative à l'interface doit être continue, autrement dit
n.umilieu 1 = n.umilieu 2. Ces conditions n'excluent pas la possibilité de glissement d'un milieu par rapport à l'autre, et on appellera une frontière sur laquelle un tel glissement est possible une « interface glissante » (en anglais : slip boundary). Par contre, si nous voulons exclure la possibilité de glissement, nous devons imposer la continuité à travers la frontière à la fois des composantes normale et tangentielle du déplacement : umilieu 1 = umilieu 2.
Une frontière de ce type est appelée « interface soudée » ou « interface non-glissante » (en anglais : slip-free boundary).
Conditions dynamiques
Pour établir les conditions dynamiques, on considère un petit parallélipipède droit (une « boîte d'allumettes ») de hauteur h. On suppose cette « boîte » partitionnée par une surface-frontière en deux volumes V1 et V2 remplis de matériau du milieu continu 1 et de matériau du milieu continu 2, respectivement. La « boîte » est limitée au-dessus par une surface S1 de normale unitaire extérieure n1, en-dessous par une surface S2 de normale unitaire extérieure n2, et sur les côtés par une surface S3 de normale unitaire extérieure n3. On applique à ce volume matériel V = V1 + V2 limité par la surface fermée S = S1 + S2 + S3 l'équation de conservation de l'impulsion sous forme intégrale, soit
. En supposant que les forces volumiques sont continues à travers la frontière, et en laissant tendre h vers zéro, nous obtenons les conditions d'interface dynamiques , en prenant bonne note du fait que pour h = 0, nous avons V = 0, S3 = 0, et n1 = − n2. Ainsi, les composantes normales du tenseur des contraintes dans leur ensemble doivent rester continues au travers l'interface.
Notes
- ↑ En géodésie et en géophysique, le terme densité dénote toujours la densité absolue, c'est-à-dire la masse volumique, et non pas la densité relative, c'est-à-dire le rapport de la masse volumique du milieu à la masse volumique de l'eau.
- ↑ Toutefois, en mécanique des milieux continus des forces internes agissant sur chaque point matériel jouent un rôle essentiel et sont généralement exprimées par la divergence du tenseur des contraintes. Par contre, en mécanique des points de Newton, les forces internes agissant sur un point quelconque se réduisent à un couple de forces opposées de même intensité, et sont donc éliminées des équations.
- ↑ On indique, si nécessaire, qu'une quantité quelconque q variable dans l'espace et dans le temps se rapporte à un point spatial (description eulérienne) en mettant explicitement l'argument x. Si q se rapporte à un point matériel (description lagrangienne), on met explicitement l'argument ξ.
- ↑ Ceci est la convention utilisée le plus souvent en physique. En géodésie et en mécanique céleste, on assimile en général la fonction potentielle (ou potentiel) V au travail à fournir contre les forces de gravité, de sorte que l'on a V = − φ. En géodésie, les cotes potentielles du champ de gravité vont décroissantes avec la distance au centre, en physique elles vont croissantes.
- ↑ L. Brillouin, Influence de la température sur l'élasticité d'un solide, Mémorial des Sciences mathématiques, fascicule 99, p. 20–23, Gauthier-Villars, Paris, 1940.
- ↑ H. Kanamori & D.L. Anderson, Importance of physical dispersion in surface waves and free oscillations, Reviews of Geophysics and Space Physics, n°15, p. 105–112, 1977.
- ↑ A.M. Dziewonski & D.L. Anderson, Preliminary reference Earth model, Physics of the Earth and Planetary Interiors, n°25, p. 297–356, 1981.
- ↑ On se restreint ici aux hypothèses usuelles faites en sismologie et en géodynamique théorique, à savoir qu'on considère uniquement des interfaces déformables simples. Une telle interface simple est une surface régulière, c'est-à-dire une surface qui en chacun de ses points possède une seule et unique normale. Des points matériels faisant initialement partie d'une telle surface resteront continûment sur cette surface pendant qu'elle se meut et se déforme. Ainsi, sur une interface simple, la formation de cavités ou l'interpénétration (donnant lieu à un mélange de particules) de deux milieux dans une portion finie de volume sont exclues de la description.
Voir aussi
Bibliographie
- R. Aris (1962). Vectors, Tensors, and the Basic Equations of Fluid Mechanics, Prentice-Hall, Englewood Cliffs (Corrected reprint, Dover Publications, 1962). ISBN 0-486-66110-5.
- L. Brillouin (1960). Les tenseurs en mécanique et en élasticité, Masson, Paris.
- C. Denis (1993). Global deformations and evolution of the Earth, Acta Geodaetica, Geophysica et Montanistica Hungarica, vol. 28, 15–131. Akadémiai Kiadó, Budapest.
- H. Jeffreys (1969). Cartesian Tensors, Cambridge University Press, Cambridge. ISBN 0-521-09191-8.
- L. Sedov (1971). A course in continuum mechanics, volume 1. Basic equations and analytical techniques, Wolters-Noordhoff Publishing, Groningen. ISBN 90-01-79680-X.
Liens internes
Liens externes
- Portail de la géodésie et de la géophysique
- Portail de la physique
- Portail de l’astronomie
Catégories : Milieux continus | Loi en physique
Wikimedia Foundation. 2010.