# É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é massique $\rho$, dans un volume $B$ délimité par une surface $\partial B$. Dans le problème qui nous occupe, cette frontière $\partial B$ 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 $v_{i}$ par rapport au système de référence choisi mais non encore spécifié, ainsi qu'une densité de moment cinétique interne $\ell _{i}$, une force volumique $f_{i}$, un champ de contraintes (ou champ de tensions) $T_{ij}$, un moment de force volumique $M_{i}$ et moment de force surfacique $M_{ij}$. 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 $x_{i}$ et les dérivées spatiales successives ${\tfrac {\partial }{\partial x_{i}}}$, ${\tfrac {\partial ^{2}}{\partial x_{i}\partial x_{j}}}$... par $\partial _{i}$, $\partial _{i}\partial _{j}$... (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
$M=\int _{B}\rho \,\mathrm {d} \tau$
où $\mathrm {d} \tau =\mathrm {d} x_{1}\mathrm {d} x_{2}\mathrm {d} x_{3}$ est un élément de volume.
L'impulsion totale est
$P_{i}=\int _{B}\rho \,v_{i}\,\mathrm {d} \tau$
et le moment cinétique total est
$L_{i}=\int _{B}\left(\varepsilon _{ijk}\,x_{j}\,\rho v_{k}+\ell _{i}\right)\,\mathrm {d} \tau$
où $\varepsilon _{ijk}$ est le symbole de Levi-Civita.
La résultante des forces appliquée au volume $B$ du corps est
$F_{i}=\int _{B}f_{i}\,\mathrm {d} \tau +\int _{\partial B}n_{j}\,T_{ij}\,\mathrm {d} \sigma$
où $\mathrm {d} \sigma$ est un élément de la surface $\partial B$ et $(n_{1},n_{2},n_{3})$ dénote un vecteur unitaire normal à la surface $\partial B$ et pointant hors du volume $B$. De la même, manière, le moment de force résultant appliqué au corps $B$ est
$M_{i}=\int _{B}\left(\varepsilon _{ijk}\,x_{j}\,f_{k}+m_{i}\right)\,\mathrm {d} \tau +\int _{\partial B}\left(\varepsilon _{ijk}\,x_{j}\,n_{q}\,T_{kq}+n_{q}\,M_{iq}\right)\,\mathrm {d} \sigma$
Ces formules sont valables en toute circonstance. Si nous admettons en outre que la surface-frontière $\partial B$ 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 flux-divergence et écrire :
$F_{i}=\int _{B}\left(f_{i}+\partial _{j}T_{ij}\right)\,\mathrm {d} \tau$
$M_{i}=\int _{B}\left(\varepsilon _{ijk}\,x_{j}\,f_{k}+m_{i}+\partial _{q}\left(\varepsilon _{ijk}\,x_{j}\,T_{kq}+M_{iq}\right)\right)\,\mathrm {d} \tau$

## 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, 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 $\partial b(t)$ est régulière, et on oriente la normale unitaire $n(t)$ en un point quelconque de $\partial b(t)$ 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 $\partial b(t)$, on a alors
$\mathrm {D} _{t}Q=\int _{b(t)}\partial _{t}q\,\mathrm {d} \tau +\int _{\partial b(t)}n_{j}qv_{j}\,\mathrm {d} \sigma$
Les notations utilisées ici sont les suivantes : $Q=\textstyle \int _{b}q\,\mathrm {d} \tau$ est une quantité physique arbitraire, et $q$ la densité volumique du champ correspondant ; les opérateurs $\mathrm {D} _{t}$ et $\partial _{t}$ sont les dérivations totales (matérielles) ${\tfrac {\mathrm {d} }{\mathrm {d} t}}$ et partielles (locales) ${\tfrac {\partial }{\partial t}}$ par rapport au temps, respectivement ; $\partial b$ est la surface régulière fermée qui limite le volume $b$, de normale extérieure unitaire $n_{i}(i=1,2,3)$. En utilisant à nouveau le théorème de flux-divergence, il vient
$\mathrm {D} _{t}Q=\int _{b(t)}\left(\partial _{t}q+\partial _{j}(qv_{j})\right)\,\mathrm {d} \tau$
Si $\kappa$ est le taux local de production ou de destruction de la quantité $Q$, c'est-à-dire $\mathrm {D} _{t}Q=\textstyle \int _{b}\kappa \,\mathrm {d} \tau$, alors la formule précédente devient
$\int _{b(t)}\left(\partial _{t}q+\partial _{j}\left(qv_{j}\right)-\kappa \right)\,\mathrm {d} \tau =0$
ou, puisque le volume $b(t)$ est arbitraire et peut donc être pris arbitrairement petit :
$\partial _{t}q+\partial _{j}\left(qv_{j}\right)=\kappa$
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 $\mathrm {D} _{t}q=q\partial _{j}v_{j}=\kappa$.
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) $\mathrm {D} _{t}q=\partial _{t}q+v_{j}\partial _{j}q$.

### Conservation de la masse, ou équation de continuité

Sous forme intégrale, la conservation de la masse s'exprime par
$\mathrm {D} _{t}M=\mathrm {D} _{t}\int _{b(t)}\sigma (x,t)\,\mathrm {d} \tau =\int _{b(t)}\left(\mathrm {D} _{t}\sigma +\rho \partial _{k}v_{k}\right)\,\mathrm {d} \tau =0$
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 : $\mathrm {D} _{t}\rho +\rho \partial _{k}v_{k}=0$.
La forme eulérienne de cette dernière équation est $\partial _{t}\rho +\partial _{k}(\rho v_{k})=0$.
On remarque qu'il s'agit bien ici d'une forme particulière de la loi de conservation générale $\partial _{t}q+\partial _{j}(qv_{j})=\kappa$, où $q=\rho$ et $\kappa =0$. Le fait que $\kappa =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=\rho v_{i}$ et en admettant explicitement que la masse se conserve aussi. On obtient
$\rho \mathrm {D} _{t}v_{i}=f_{i}+\partial _{k}T_{ik}(\xi ,t)$
En effet, l'équation fondamentale de la dynamique exprime le fait que le taux de variation de l'impulsion $\mathrm {D} _{t}P_{i}$ est équilibré par la résultante des forces $F_{i}$ agissant sur le volume. On a donc
$\mathrm {D} _{t}\int _{b}\rho v_{i}\,\mathrm {d} \tau =\int _{b}\left(f_{i}+\partial _{k}T_{ik}\right)\,\mathrm {d} \tau$
Si nous transformons le membre de gauche au moyen du théorème de Reynolds, nous obtenons
$\mathrm {D} _{t}\int _{b(t)}\left[\partial _{t}\left(\rho v_{i}\right)+\partial _{k}\left(\rho v_{i}v_{k}\right)\right]\,\mathrm {d} \tau =\int _{b(t)}\left(f_{i}+\partial _{k}T_{ik}\right)\,\mathrm {d} \tau$
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 $\partial _{t}(\rho v_{i})+\partial _{k}(\rho v_{i}v_{k})=f_{i}+\partial _{k}T_{ik}$.
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 $T_{ik}$ est différentiable. En développant le membre de gauche, on obtient
$\rho \left[\partial _{t}v_{i}+v_{k}\partial _{k}v_{i}\right]+v_{i}\left[\partial _{t}\rho +\partial _{k}\left(\rho v_{k}\right)\right]=f_{i}+\partial _{k}T_{ik}$
Le terme entre crochets dans le membre de gauche représente l'accélération $\mathrm {D} _{t}v_{i}$. 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 $\rho \mathrm {D} _{t}v_{i}=f_{i}+\partial _{k}T_{ik}$.

### 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
$\mathrm {D} _{t}\ell _{i}=m_{i}+\partial _{k}M_{ik}+\varepsilon _{ijk}T_{jk}$
En effet, le taux de variation du moment cinétique $\mathrm {D} _{t}L_{i}$ dans un volume arbitraire $b$ est équilibré par le moment de force résultant $M_{i}$ agissant sur ce volume. Nous avons donc
$\mathrm {D} _{t}\int _{b}\left(\varepsilon _{ijk}x_{j}\rho v_{k}+\ell _{i}\right)\,\mathrm {d} \tau =\int _{b}\left[\varepsilon _{ijk}x_{j}f_{k}+m_{i}+\partial _{q}\left(\varepsilon _{ijk}x_{j}T_{kq}+M_{iq}\right)\right]\mathrm {d} \tau$
Comme précédemment, on transforme le membre de gauche en appliquant le théorème de Reynolds pour obtenir
$\int _{b}\left[\partial _{t}\left(\varepsilon _{ijk}x_{j}\rho v_{k}+\ell _{i}\right)+\partial _{q}\left(\rho \varepsilon _{ijk}x_{j}v_{k}v_{q}+\ell _{i}v_{q}\right)\right]\,\mathrm {d} \tau =\int _{b}\left[\varepsilon _{ijk}x_{j}f_{k}+m_{i}+\partial _{q}\left(\varepsilon _{ijk}x_{j}T_{kq}+M_{iq}\right)\right]\,\mathrm {d} \tau$
En faisant tendre le volume arbitraire $b$ vers zéro, on obtient l'équation exprimant localement la conservation du moment cinétique, à savoir
$\partial _{t}\left(\rho \varepsilon _{ijk}x_{j}v_{k}+\ell _{i}\right)+\partial _{q}\left(\rho \varepsilon _{ijk}x_{j}v_{k}v_{q}+\ell _{i}v_{q}\right)=\varepsilon _{ijk}x_{j}f_{k}+m_{i}+\partial _{q}\left(\varepsilon _{ijk}x_{j}T_{kq}+M_{iq}\right)$
Cette équation peut être réarrangée comme suit :
$\langle \partial _{t}\ell _{i}+v_{q}\partial _{q}\ell _{i}\rangle +\varepsilon _{ijk}x_{j}v_{k}\left[\partial _{t}\rho +\partial _{q}\left(\rho v_{k}\right)\right]+\varepsilon _{ijk}x_{j}\{\rho \partial _{t}v_{k}+\rho v_{q}\partial _{q}v_{k}-f_{k}-\partial _{q}T_{kq}\}=m_{i}+\partial _{q}M_{iq}+\varepsilon _{ijk}\left(T_{kj}-\rho v_{j}v_{k}\right)$
Dans le membre de gauche, les termes entre chevrons représentent la quantité $\mathrm {D} _{t}\ell _{i}$, 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 $\rho \varepsilon _{ijk}v_{j}v_{k}$ est identiquement nul, puisque $v_{j}v_{k}$ est un tenseur symétrique. Admettant que masse et impulsion soient conservées, on obtient donc $\mathrm {D} _{t}\ell _{i}=m_{i}+\partial _{k}M_{ik}+\varepsilon _{ijk}T_{jk}$.

### Symétrie du tenseur des contraintes

Si, en outre, nous supposons que le moment cinétique intrinsèque $\ell _{i}$ d'un point matériel quelconque ne change pas au cours du mouvement de déformation, autrement dit si $\mathrm {D} _{t}\ell _{i}=m_{i}+\partial _{k}M_{ik}$, alors le tenseur des contraintes $T_{jk}$ est symétrique : $\varepsilon _{ijk}T_{jk}=0\equiv T_{jk}=T_{kj}$.

### 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 :
$\partial _{t}\left(\rho v^{2}+U\right)+\partial _{j}\left[\left({\frac {1}{2}}\rho v^{2}+U\right)v_{j}-v_{k}T_{kj}+H_{j}\right]=v_{k}f_{k}+\Gamma$
Le symbole $U$ désigne ici l'énergie interne par unité de volume, $H_{i}$ est le vecteur densité de flux de chaleur, et $\Gamma$ 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 $v^{2}$ représente bien sûr le carré $v_{k}v_{k}$ 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 à la suite d'un séisme, à la suite d'une variation de la longueur du jour, à la suite deà l'application d'un potentiel de marée ou à la suite de 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 $\textstyle Q=\int _{B}\phi \,\mathrm {d} \tau$ est la quantité physique qui se conserve ; elle représente, selon les cas, la masse totale du corps $B$, sa quantité de mouvement totale $P_{i}$, son moment cinétique total $L_{i}$, ou son énergie totale $E$. La forme eulérienne générale d'une loi de conservation est $\partial _{t}\phi +\partial _{j}k_{j}=q$.
Principales lois de conservation de la mécanique des milieux continus :
| $Q$ | $\phi$ | $k_{j}$ | $q$ |
| $M$ | $\rho$ | $\rho v_{j}$ | 0 |
| $P_{i}$ | $\rho v_{i}$ | $\rho v_{i}v_{j}-T_{ij}$ | $f_{i}$ |
| $L_{i}$ | $\ell _{i}$ | $\ell _{i}v_{j}-M_{ij}$ | $m_{i}+\varepsilon _{ijk}T_{kj}$ |
| $E$ | ${\frac {1}{2}}\rho v^{2}+U$ | $\left({\frac {1}{2}}\rho v^{2}+U\right)v_{j}-v_{k}T_{kj}+H_{j}$ | $v_{k}f_{k}+\Gamma$ |

## Hypothèses et lois pour décrire les déformations de la Terre

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 $f_{i}$ se compose d'une force gravifique volumique $f_{i}^{\text{g}}$ causée par l'autogravitation et par l'attraction gravifique de corps extérieurs, d'une force électromagnétique volumique $f_{i}^{\text{em}}$, de forces volumiques inertielles (Coriolis et axifuge) $f_{i}^{\text{in}}$, de forces de marée $f_{i}^{\text{m}}$, et de forces volumiques causées par des charges extérieures $f_{i}^{\text{ch}}$. D'autre part, les forces surfaciques $T_{i}(x,t,n)$ dépendant de l'orientation $n$ qui contribuent en général au tenseur des contraintes $T_{ij}(x,t)$ via les relations de Cauchy $T_{i}(x,t,n)=T_{ij}(x,t)n_{j}(x,t,n)$ sont les forces élastiques qui engendrent un tenseur des contraintes élastiques $T_{ij}^{\text{el}}$, les forces superficielles électromagnétiques qui engendrent un tenseur de Maxwell $T_{ij}^{\text{em}}$, des forces visqueuses qui engendrent un tenseur de Reynolds un tenseur des contraintes élastiques $T_{ij}^{\text{v}}$, des forces de frottement de contact non visqueux qui engendrent un tenseur des forces de friction un tenseur des contraintes élastiques $T_{ij}^{\text{f}}$, et de toute autre force surfacique engendrant un tenseur des contraintes un tenseur des contraintes élastiques $T_{ij}^{\text{s}}$.
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 $\phi ^{g}$ que nous assimilerons ici à l'énergie potentielle gravifique :
$f^{g}=\rho g=-\rho \nabla \phi ^{g}$
Le potentiel $\phi ^{g}$ est déterminé par l'équation de Poisson
$\nabla ^{2}\phi ^{g}=4\pi G\rho$
où $G$ est la constante de gravitation. Il convient de noter que les propriétés du potentiel gravifique impliquent que $\phi ^{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 $\xi$ 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 $\sigma _{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(\xi ,t)$ tandis que le flux variable de ce courant produit un champ d'induction électrique $E(\xi ,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 $f^{\text{em}}=c^{-1}J.B$, avec $J=\sigma _{e}(E+c^{-1}v.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 $(\sigma _{e}=0)$, la force électromagnétique serait nulle $(f^{\text{em}}=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 :
$f^{\text{in}}=\rho \nabla \left({\frac {1}{2}}|\Omega .r|\right)^{2}-2\rho \Omega .v$
où $\Omega$ est le vecteur rotation instantanée, $r$ le vecteur position \mathrm du point matériel $\xi$, 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 $f^{\text{m}}=-\rho \nabla W^{\text{m}}$ dérivant aussi d'un potentiel, le potentiel de marée $W^{\text{m}}$. 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 :
$W^{\text{m}}=\Sigma _{A}GM_{A}d_{A}^{-1}\sum _{n=2}^{\infty }\left({\frac {r}{\mathrm {d} _{A}}}\right)^{n}P_{n}\left(\cos z_{A}\right)$
Ici, $G$ dénote la constante de gravitation $(G=6{,}673\times 10^{-11}{\rm {\;m^{3}\cdot s^{-2}\cdot kg^{-1})}}$, $M_{A}$ est la masse de l'astre générateur de la marée, $d_{A}$ est la distance du centre de masse de la Terre au centre de masse de $A$, $z_{A}$ 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 $P_{n}$ 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
$T_{ik}(\xi ,t)=\left(\kappa -{\frac {2}{3}}\mu \right)\partial _{j}u_{j}\delta _{ik}+\mu (\partial _{i}u_{k}+\partial _{k}u_{i}),$
où $\mathbf {u}$ est le champ des déplacements et où les constantes de matériau $\kappa$ et $\mu$ 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. D'autre part, on peut tenir compte d'une faible anélasticité en faisant dépendre les coefficients $\kappa$ et $\mu$ de la fréquence ; c'est l'approche adoptée notamment par Kanamori et Anderson et par Dziewonski et Anderson.

## 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 : 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(\xi ,t),t]={\text{cste}}$ à un instant donné doivent encore appartenir à la même surface à un autre instant, c'est-à-dire $\sigma F^{\text{milieu 1}}=\sigma F^{\text{milieu 2}}$ ou encore, sous forme eulérienne :
$\partial F^{\text{milieu 1}}+u^{\text{milieu 1}}.\nabla F^{\text{milieu 1}}=\partial F^{\text{milieu 2}}+u^{\text{milieu 2}}.\nabla F^{\text{milieu 2}}$
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.u^{\text{milieu 1}}=n.u^{\text{milieu 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 : $u^{\text{milieu 1}}=u^{\text{milieu 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élépipè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 $V_{1}$ et $V_{2}$ 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 $S_{1}$ de normale unitaire extérieure $n_{1}$, en dessous par une surface $S_{2}$ de normale unitaire extérieure $n_{2}$, et sur les côtés par une surface $S_{3}$ de normale unitaire extérieure $n_{3}$. On applique à ce volume matériel $V=V_{1}+V_{2}$ limité par la surface fermée $S=S_{1}+S_{2}+S_{3}$ l'équation de conservation de l'impulsion sous forme intégrale, soit
$\mathrm {D} _{t}\int \!\!\!\!\!\!\!\subset \!\!\!\!\!\int \!\!\!\!\!\supset \!\!\!\!\!\!\!\int _{V(t)}\rho v_{i}\;\mathrm {d} \tau =\int \!\!\!\!\!\!\!\subset \!\!\!\!\!\int \!\!\!\!\!\supset \!\!\!\!\!\!\!\int _{V(t)}f_{i}\;\mathrm {d} \tau +\int \!\!\!\!\!\!\!\subset \!\!\supset \!\!\!\!\!\!\!\int _{S(t)}T_{ik}n_{k}\;\mathrm {d} \sigma$
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 $T_{ik}^{\text{milieu 1}}n_{k}=T_{ik}^{\text{milieu 2}}n_{k}$, en prenant bonne note du fait que pour $h=0$, nous avons $V=0$, $S_{3}=0$, et $n_{1}=-n_{2}$. Ainsi, les composantes normales du tenseur des contraintes dans leur ensemble doivent rester continues au travers l'interface.