Navier-Cauchy-Gleichungen

Die Navier-Cauchy-, Navier- oder Navier-Lamé-Gleichungen (nach Claude Louis Marie Henri Navier, Augustin-Louis Cauchy und Gabriel Lamé) sind ein mathematisches Modell der Bewegung – inklusive Deformation – von elastischen Festkörpern. Bei der Herleitung der Modellgleichungen wird sowohl geometrische- als auch physikalische Linearität (lineare Elastizität) vorausgesetzt. Die Gleichungen lauten koordinatenunabhängig vektoriell

$\begin{array}{rl}\varrho \ddot{\overrightarrow{u}}=& G[\mathrm{\Delta }\overrightarrow{u}+\frac{1}{1-2\nu }\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})]+\varrho \overrightarrow{k}\\ =& G[\mathrm{\Delta }\overrightarrow{u}+\frac{1}{1-2\nu }\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{(}\mathrm{d}\mathrm{i}\mathrm{v}(\overrightarrow{u}))]+\varrho \overrightarrow{k}\end{array}$

oder in kartesischen Koordinaten

$\varrho \frac{{\partial }^2{u}_i}{\partial t^2}=G\sum _{k=1}^3[\frac{{\partial }^2{u}_i}{\partial x_k^2}+\frac{1}{1-2\nu }\frac{{\partial }^2{u}_k}{\partial x_i\partial x_k}]+\varrho k_i,i=1,2,3$

Es handelt sich um ein partielles Differentialgleichungssystem zweiter Ordnung in drei unbekannten Verschiebungen $\overrightarrow{u}(\overrightarrow{x},t)$, die im Allgemeinen sowohl vom Ort $\overrightarrow{x}$ als auch von der Zeit t abhängen. Verschiebungen sind die Wege, die die Partikel eines Körpers bei einer Bewegung – inklusive Deformation – zurücklegen. Die Materialparameter ϱ, G und ν sind die Dichte, der Schubmodul bzw. die Querkontraktionszahl, 𝜵, Δ = 𝜵², grad und div der Nabla-, Laplace-, Gradienten- bzw. Divergenzoperator und $\varrho \overrightarrow{k}$ repräsentiert eine volumenverteilte Kraft, wie die Schwerkraft eine ist.

Jedes Material im festen Aggregatzustand hat einen mehr oder weniger ausgeprägten linear-elastischen Bereich, zumindest bei kleinen und langsamen Verformungen, die bei vielen Anwendungen, vor allem im technischen Bereich, vorliegen.

Historisches

Claude Louis Marie Henri Navier leitete diese, nach ihm benannte Gleichung 1821 aus einem molekularen Modell ab, das auf Materialien mit identischen ersten und zweiten Lamé-Konstanten beschränkt ist. Die allgemeinere, hier vorgestellte Gleichung mit zwei verschiedenen Elastizitätskonstanten, erschien erstmals in einer Arbeit von Cauchy 1828.

Herleitung

Ausgangspunkt ist das erste Cauchy-Euler’sche Bewegungsgesetz bei kleinen Verschiebungen

$\varrho \ddot{\overrightarrow{u}}=\mathrm{\nabla }\cdot \mathit{\sigma }+\varrho \overrightarrow{k}$

das der Impulsbilanz entspricht. Zusätzlich zu den eingangs beschriebenen Variablen tritt hier der infolge des Drallsatzes symmetrische Spannungstensor σ auf. Dessen Abhängigkeit von den Verschiebungen ergibt sich mit dem linearisierten Verzerrungstensor

$\mathit{\epsilon }=\frac{1}{2}(\mathrm{\nabla }\otimes \overrightarrow{u}+\mathrm{\nabla }\otimes {\overrightarrow{u}}^{\mathrm{\top }})=\frac{1}{2}(\mathrm{grad}(\overrightarrow{u}{)}^{\mathrm{\top }}+\mathrm{grad}(\overrightarrow{u}))$

aus dem Hooke’schen Gesetz:

$\begin{array}{rl}\mathit{\sigma }=& 2G[\mathit{\epsilon }+\frac{\nu }{1-2\nu }\mathrm{S}\mathrm{p}(\mathit{\epsilon })\mathbf{𝟏}]\\ =& G[\mathrm{\nabla }\otimes \overrightarrow{u}+\mathrm{\nabla }\otimes {\overrightarrow{u}}^{\mathrm{\top }}+\frac{\nu }{1-2\nu }\mathrm{S}\mathrm{p}(\mathrm{\nabla }\otimes \overrightarrow{u}+\mathrm{\nabla }\otimes {\overrightarrow{u}}^{\mathrm{\top }})\mathbf{𝟏}]\\ =& G[\mathrm{\nabla }\otimes \overrightarrow{u}+\mathrm{\nabla }\otimes {\overrightarrow{u}}^{\mathrm{\top }}+\frac{2\nu }{1-2\nu }(\mathrm{\nabla }\cdot \overrightarrow{u})\mathbf{𝟏}]\end{array}$

Das Superskript ^⊤ steht für die Transposition des dyadischen Produkts ⊗, 1 für den Einheitstensor und der Operator Sp extrahiert die von der Transposition unbeeinflusste Spur, die bei dem Gradient eines Vektorfeldes gleich der Divergenz desselben ist. Die im ersten Cauchy-Eulerschen Bewegungsgesetz auftretende Divergenz wird bereitgestellt:

$\begin{array}{rl}\mathrm{\nabla }\cdot \mathit{\sigma }=& G[\mathrm{\Delta }\overrightarrow{u}+\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})+\frac{2\nu }{1-2\nu }\mathrm{\nabla }\cdot ((\mathrm{\nabla }\cdot \overrightarrow{u})\mathbf{𝟏}\left)\right]\\ =& G[\mathrm{\Delta }\overrightarrow{u}+\frac{1}{1-2\nu }\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})]\end{array}$

In Kombination mit dem obigen Bewegungsgesetz $(\varrho \ddot{\overrightarrow{u}}=\mathrm{\nabla }\cdot \mathit{\sigma }+\varrho \overrightarrow{k})$ führt das auf die Navier-Cauchy-Gleichungen:

$\begin{array}{rl}\varrho \ddot{\overrightarrow{u}}=& G[\mathrm{\Delta }\overrightarrow{u}+\frac{1}{1-2\nu }\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})]+\varrho \overrightarrow{k}\\ =& \mu \mathrm{\Delta }\overrightarrow{u}+(\lambda +\mu )\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})+\varrho \overrightarrow{k}\end{array}$

In der letzten Gleichung wurden alternativ die erste und zweite Lamé-Konstante λ und μ eingesetzt. Gelegentlich ist es bequem noch die Identität $\mathrm{\Delta }\overrightarrow{u}=\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})-\mathrm{\nabla }\times (\mathrm{\nabla }\times \overrightarrow{u})$ auszunutzen:

$\varrho \ddot{\overrightarrow{u}}=(\lambda +2\mu )\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})-\mu \mathrm{\nabla }\times (\mathrm{\nabla }\times \overrightarrow{u})+\varrho \overrightarrow{k}$

Das Kreuzprodukt × mit dem Nabla-Operator bildet die Rotation eines Vektorfeldes.

Koordinaten der Navier-Cauchy Gleichungen

Die Navier-Cauchy Gleichungen sind in kartesischen-, zylinder- und Kugelkoordinaten bekannt.

Kartesische Koordinaten

$\begin{array}{rl}\varrho {\ddot{u}}_x=& \mu \mathrm{\Delta }{u}_x+(\lambda +\mu )\frac{\partial }{\partial x}(\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z})+\varrho k_x\\ \varrho {\ddot{u}}_y=& \mu \mathrm{\Delta }{u}_y+(\lambda +\mu )\frac{\partial }{\partial y}(\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z})+\varrho k_y\\ \varrho {\ddot{u}}_z=& \mu \mathrm{\Delta }{u}_z+(\lambda +\mu )\frac{\partial }{\partial z}(\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z})+\varrho k_z\end{array}$

mit $\mathrm{\Delta }f=\frac{{\partial }^{2}f}{\partial x^2}+\frac{{\partial }^{2}f}{\partial y^2}+\frac{{\partial }^{2}f}{\partial z^2}$.

Zylinderkoordinaten

$\begin{array}{rl}\varrho {\ddot{u}}_{\rho }=& \mu (\mathrm{\Delta }{u}_{\rho }-\frac{u_{\rho }}{{\rho }^2}-\frac{2}{{\rho }^2}\frac{\partial u_{\phi }}{\partial \phi })\\ & +(\lambda +\mu )\frac{\partial }{\partial \rho }(\frac{1}{\rho }\frac{\partial }{\partial \rho }(\rho u_{\rho })+\frac{1}{\rho }\frac{\partial u_{\phi }}{\partial \phi }+\frac{\partial u_z}{\partial z})+\varrho k_{\rho }\\ \varrho {\ddot{u}}_{\phi }=& \mu (\mathrm{\Delta }{u}_{\phi }-\frac{u_{\phi }}{{\rho }^2}+\frac{2}{{\rho }^2}\frac{\partial u_{\rho }}{\partial \phi })\\ & +(\lambda +\mu )\frac{1}{\rho }\frac{\partial }{\partial \phi }(\frac{1}{\rho }\frac{\partial }{\partial \rho }(\rho u_{\rho })+\frac{1}{\rho }\frac{\partial u_{\phi }}{\partial \phi }+\frac{\partial u_z}{\partial z})+\varrho k_{\phi }\\ \varrho {\ddot{u}}_z=& \mu \mathrm{\Delta }{u}_z\\ & +(\lambda +\mu )\frac{\partial }{\partial z}(\frac{1}{\rho }\frac{\partial }{\partial \rho }(\rho u_{\rho })+\frac{1}{\rho }\frac{\partial u_{\phi }}{\partial \phi }+\frac{\partial u_z}{\partial z})+\varrho k_z\end{array}$

mit $\mathrm{\Delta }f=\frac{1}{\rho }\frac{\partial }{\partial \rho }\left(\rho \frac{\partial f}{\partial \rho }\right)+\frac{1}{{\rho }^2}\frac{{\partial }^2{u}_{\phi }}{\partial {\phi }^2}+\frac{{\partial }^2{u}_z}{\partial z^2}$.

Kugelkoordinaten

Die Formulierung benutzt den Sinus und Cosinus sin bzw. cos und den Tangens tan.

$\begin{array}{rl}\varrho {\ddot{u}}_r=& \mu (\mathrm{\Delta }{u}_r-\frac{2u_r}{r^2}-\frac{2}{r^2\sin (\vartheta )}\frac{\partial u_{\phi }}{\partial \phi }-\frac{2u_{\vartheta }}{r^2\tan (\vartheta )}-\frac{2}{r^2}\frac{\partial u_{\vartheta }}{\partial \vartheta })+\varrho k_r\\ & +(\lambda +\mu )\frac{\partial }{\partial r}(\frac{1}{r^2}\frac{\partial }{\partial r}(r^2{u}_r)+\frac{1}{r\sin (\vartheta )}\frac{\partial }{\partial \vartheta }(u_{\vartheta }\sin (\vartheta ))+\frac{1}{r\sin (\vartheta )}\frac{\partial u_{\phi }}{\partial \phi })\\ \varrho {\ddot{u}}_{\vartheta }=& \mu (\mathrm{\Delta }{u}_{\vartheta }+\frac{2}{r^2}\frac{\partial u_r}{\partial \vartheta }-\frac{u_{\vartheta }}{r^2{\sin }^2(\vartheta )}-2\frac{\cos (\vartheta )}{r^2{\sin }^2(\vartheta )}\frac{\partial u_{\phi }}{\partial \phi })+\varrho k_{\vartheta }\\ & +\frac{\lambda +\mu }{r}\frac{\partial }{\partial \vartheta }(\frac{1}{r^2}\frac{\partial }{\partial r}(r^2{u}_r)+\frac{1}{r\sin (\vartheta )}\frac{\partial }{\partial \vartheta }(u_{\vartheta }\sin (\vartheta ))+\frac{1}{r\sin (\vartheta )}\frac{\partial u_{\phi }}{\partial \phi })\\ \varrho {\ddot{u}}_{\phi }=& \mu (\mathrm{\Delta }{u}_{\phi }-\frac{u_{\phi }}{r^2{\sin }^2(\vartheta )}+\frac{2}{r^2\sin (\vartheta )}\frac{\partial u_r}{\partial \phi }+\frac{2\cos (\vartheta )}{r^2{\sin }^2(\vartheta )}\frac{\partial u_{\vartheta }}{\partial \phi })+\varrho k_{\phi }\\ & +\frac{\lambda +\mu }{r\sin (\vartheta )}\frac{\partial }{\partial \phi }(\frac{1}{r^2}\frac{\partial }{\partial r}(r^2{u}_r)+\frac{1}{r\sin (\vartheta )}\frac{\partial }{\partial \vartheta }(u_{\vartheta }\sin (\vartheta ))+\frac{1}{r\sin (\vartheta )}\frac{\partial u_{\phi }}{\partial \phi })\end{array}$

mit $\mathrm{\Delta }f=\frac{{\partial }^{2}f}{\partial {\rho }^2}+\frac{2}{r}\frac{\partial f}{\partial \rho }+\frac{1}{r^2\sin (\vartheta )}\frac{\partial }{\partial \vartheta }(\sin (\vartheta )\frac{\partial f}{\partial \vartheta })+\frac{1}{r^2{\sin }^2(\vartheta )}\frac{{\partial }^{2}f}{\partial {\phi }^2}$.

Randbedingungen

Im konkreten Berechnungsfall der Navier-Cauchy-Gleichungen sind Randbedingungen zu definieren. Als geometrische oder Dirichlet-Randbedingungen werden in den Auflagern die Verschiebung vorgegeben, oftmals ganz unterdrückt. Die dynamischen oder Neumann-Randbedingungen entsprechen flächenverteilten Kräften $\overrightarrow{t}$ (Vektoren mit der Dimension Kraft pro Fläche), die auf Oberflächen des Körpers wirken.

Lösungsmethoden

Für einfache Fälle, siehe das Beispiel unten, gerade Stäbe und ebene Scheiben können analytische Lösungen angegeben werden. Bei unregelmäßig geformten Körpern bietet sich als numerisches Werkzeug die Verschiebungsmethode in der Finite-Elemente-Methode an.

Spezialfälle

Ebene Probleme

Analytische Lösungen für die dreidimensionalen Navier-Cauchy-Gleichungen sind schwer aufzufinden. Die meisten erlangten Lösungen basieren auf 2-dimensionalen Modellen:

Ebener Verzerrungszustand

Hier verschwindenden alle Verzerrungen in Dickenrichtung was bei sehr dicken oder langen zylindrischen Körpern vorkommt.

Ebener Spannungszustand

Bei ihm verschwindenden Spannungen in Dickenrichtung, was in dünnen Flächenträgern vorkommt. Im Inneren von Körpern kann dieser Zustand bei ortsabhängigen Spannungen nur näherungsweise erfüllt sein.

Verallgemeinerter ebener Spannungszustand

Wegen dieses Mangels des ebenen Spannungszustands wurde eine Theorie entwickelt, die mit über die Dicke gemittelten Größen arbeitet.

Antiebener Verzerrungszustand

Dieses Modell setzt ausschließlich Verschiebungen in Dickenrichtung voraus.

Axialsymmetrischer Verzerrungszustand.

Dies ist die Formulierung in #Zylinderkoordinaten, wo eine Abhängigkeit von der Umfangsrichtung φ entfällt.

Der ebene Verzerrungszustand und der ebene Spannungszustand sind die fundamentalen ebenen Theorien und liefern sehr ähnliche Feldgleichungen. Diese können auf nur eine Gleichung in einer unbekannten, die Airy’sche Spannungsfunktion zurückgeführt werden.

Die weitreichenden Eigenschaften komplexwertiger Funktionen können mit dem Ansatz U := u+iv ausgenutzt werden, wo i die imaginäre Einheit ist. Im statischen Fall führt dieser Ansatz auf

$2\mu \frac{{\partial }^{2}U}{\partial \overline{z}\partial z}+(\lambda +\mu )\frac{\partial }{\partial \overline{z}}(\frac{\partial U}{\partial z}+\frac{\partial \overline{U}}{\partial z})=0$

Darin stellt der Überstrich den konjugiert komplexen Wert dar. Die allgemeine Lösung dieser Gleichung lässt sich mit zwei zu bestimmenden Funktionen Φ und Ψ darstellen:

$2\mu U=\gamma \mathrm{\Phi }(z)-z\overline{{\mathrm{\Phi }}^{\prime }(z)}-\overline{\mathrm{\Psi }(z)}$

Der Parameter γ hängt von der Querkontraktionszahl ab:

$\gamma =\{\begin{array}{ll}3-4\nu & \text{im ebenen Spannungszustand}\\ \frac{3-\nu }{1+\nu }& \text{im ebenen Verzerrungszustand}\end{array}$

Siehe Airysche Spannungsfunktion#Darstellung mit komplexen Funktionen.

Harmonische Schwerkraft

Im Gleichgewicht schreiben sich die Navier-Cauchy-Gleichungen

$\overrightarrow{0}=(\lambda +2\mu )\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})-\mu \mathrm{\nabla }\times (\mathrm{\nabla }\times \overrightarrow{u})+\varrho \overrightarrow{k}$

Die Divergenz (𝜵·) und Rotation (𝜵×) dieser Gleichung liefern:

$\begin{array}{rl}(\lambda +2\mu )\mathrm{\Delta }(\mathrm{\nabla }\cdot \overrightarrow{u})=(\lambda +2\mu )\mathrm{\nabla }\cdot (\mathrm{\Delta }\overrightarrow{u})=& -\mathrm{\nabla }\cdot (\varrho \overrightarrow{k})\\ \mu \mathrm{\nabla }\times (\mathrm{\Delta }\overrightarrow{u})=& -\mathrm{\nabla }\times (\varrho \overrightarrow{k})\end{array}$

Wenn die Schwerkraft $\varrho \overrightarrow{k}$ sowohl divergenz- als auch rotationsfrei ist, dann resultiert

$\mathrm{\nabla }\cdot (\mathrm{\Delta }\overrightarrow{u})=0\text{und}\mathrm{\nabla }\times (\mathrm{\Delta }\overrightarrow{u})=\overrightarrow{0}$

Ein Vektorfeld, dessen Divergenz und Rotation verschwinden, ist harmonisch, so dass im Gleichgewicht von $\mathrm{\nabla }\cdot (\varrho \overrightarrow{k})=0$ und $\mathrm{rot}(\varrho \overrightarrow{k})=\overrightarrow{0}$ auf

$\mathrm{\Delta }\mathrm{\Delta }\overrightarrow{u}=\overrightarrow{0}$

geschlossen werden kann. Letzteres ist die sogenannte biharmonische Differentialgleichung.

Inkompressibilität

Bei Inkompressibilität verschwindet die Spur des Verzerrungstensors, denn sie gibt die Volumendehnung an:

$\mathrm{Sp}(\mathit{\epsilon })=\frac{1}{2}\mathrm{Sp}(\mathrm{\nabla }\otimes \overrightarrow{u}+\mathrm{\nabla }\otimes {\overrightarrow{u}}^{\mathrm{\top }})=\mathrm{Sp}(\mathrm{\nabla }\otimes \overrightarrow{u})=\mathrm{\nabla }\cdot \overrightarrow{u}=0$

Bei Inkompressibilität ist der Kugel-Anteil des Spannungstensors unbestimmt und wird zum Drucktensor zusammengefasst:

$\mathit{\sigma }=-p\mathbf{𝟏}+2G\mathit{\epsilon }=-p\mathbf{𝟏}+G\mathrm{\nabla }\otimes \overrightarrow{u}+G\mathrm{\nabla }\otimes {\overrightarrow{u}}^{\mathrm{\top }}$

Der Skalar p ist der Druck, der sich erst im konkreten Berechnungsfall aus den Randbedingungen und Naturgesetzen ergibt. Für die Divergenz des Spannungstensors hat dies die Konsequenz (siehe die obigen Anmerkungen):

$\begin{array}{rl}\mathrm{\nabla }\cdot \mathit{\sigma }=& -\mathrm{\nabla }\cdot (p\mathbf{𝟏})+G\mathrm{\nabla }\cdot (\mathrm{\nabla }\otimes \overrightarrow{u})+G\mathrm{\nabla }\cdot (\mathrm{\nabla }\otimes \overrightarrow{u}{)}^{\mathrm{\top }}\\ =& -\mathrm{\nabla }p+G\mathrm{\Delta }\overrightarrow{u}+G\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})\\ =& -\mathrm{\nabla }p+G\mathrm{\Delta }\overrightarrow{u}\end{array}$

Das erste Cauchy-Euler’sche Bewegungsgesetz schreibt sich dann

$\varrho \ddot{\overrightarrow{u}}=\mathrm{\nabla }\cdot \mathit{\sigma }+\varrho \overrightarrow{k}=-\mathrm{\nabla }p+G\mathrm{\Delta }\overrightarrow{u}+\varrho \overrightarrow{k}.$

Zu diesen drei Gleichungen in den vier Unbekannten $\{p,\overrightarrow{u}\}$ wird noch $\mathrm{\nabla }\cdot \overrightarrow{u}=0$ zum Abschluss benötigt.

Wellengleichungen

Division der Navier-Cauchy-Gleichung mit Rotation durch die Dichte ergibt bei vernachlässigbarer Schwerkraft:

$\begin{array}{rl}\ddot{\overrightarrow{u}}=& \frac{\lambda +2\mu }{\varrho }\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})-\frac{\mu }{\varrho }\mathrm{\nabla }\times (\mathrm{\nabla }\times \overrightarrow{u})\\ =:& c_l^2\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{u})-c_t^2\mathrm{\nabla }\times (\mathrm{\nabla }\times \overrightarrow{u})\end{array}$

Die Faktoren c_l und c_t haben die Dimension einer Geschwindigkeit. Gemäß dem Helmholtz-Theorem lässt sich jedes reale Verschiebungsfeld eindeutig in einen divergenzfreien und einen rotationsfreien Anteil zerlegen:

$\overrightarrow{u}={\overrightarrow{u}}_l+{\overrightarrow{u}}_t$ mit $\mathrm{\nabla }\times {\overrightarrow{u}}_l=\overrightarrow{0}$ und $\mathrm{\nabla }\cdot {\overrightarrow{u}}_t=0$

Für den rotationsfreien Anteil gibt es ein Skalarpotential, dessen Gradientenfeld er ist, und für den divergenzfreien Anteil existiert ein Vektorfeld, dessen Rotation er ist:

$\begin{array}{rl}{\overrightarrow{u}}_l=& \mathrm{\nabla }\phi ,{\overrightarrow{u}}_t=\mathrm{\nabla }\times \overrightarrow{a}\\ \overrightarrow{u}=& \mathrm{\nabla }\phi +\mathrm{\nabla }\times \overrightarrow{a}\end{array}$

Dies in die Navier-Cauchy-Gleichung eingesetzt

$\ddot{\overrightarrow{u}}=\frac{{\partial }^2}{\partial t^2}\mathrm{\nabla }\phi +\frac{{\partial }^2}{\partial t^2}\mathrm{\nabla }\times \overrightarrow{a}=c_l^2\mathrm{\nabla }(\mathrm{\nabla }\cdot \mathrm{\nabla }\phi )-c_t^2\mathrm{\nabla }\times (\mathrm{\nabla }\times (\mathrm{\nabla }\times \overrightarrow{a}))$

zeigt mit $\mathrm{\nabla }\times (\mathrm{\nabla }\times \overrightarrow{f})=\mathrm{\nabla }(\mathrm{\nabla }\cdot \overrightarrow{f})-\mathrm{\Delta }\overrightarrow{f}$ und Δ = 𝜵²:

$\begin{array}{rl}\frac{{\partial }^2}{\partial t^2}\mathrm{\nabla }\phi +\frac{{\partial }^2}{\partial t^2}\mathrm{\nabla }\times \overrightarrow{a}=& c_l^2\mathrm{\nabla }(\mathrm{\Delta }\phi )-c_t^2\mathrm{\nabla }(\mathrm{\nabla }\cdot (\mathrm{\nabla }\times \overrightarrow{a}))+c_t^2\mathrm{\Delta }\mathrm{\nabla }\times \overrightarrow{a}\\ =& c_l^2\mathrm{\nabla }\mathrm{\Delta }\phi +c_t^2\mathrm{\nabla }\times \mathrm{\Delta }\overrightarrow{a}\end{array}$

oder umgestellt mithilfe vom Satz von Schwarz, der die Reihenfolge der Ableitungen freigibt:

$\mathrm{\nabla }(\frac{{\partial }^2\phi }{\partial t^2}-c_l^2\mathrm{\Delta }\phi )+\mathrm{\nabla }\times (\frac{{\partial }^2\overrightarrow{a}}{\partial t^2}-c_t^2\mathrm{\Delta }\overrightarrow{a})=\overrightarrow{0}$

Diese Gleichung wird gewiss erfüllt, wenn die in den Klammern stehenden Terme verschwinden, die Wellengleichungen darstellen:

$\begin{array}{rl}\frac{{\partial }^2\phi }{\partial t^2}-c_l^2\mathrm{\Delta }\phi =& 0\\ \frac{{\partial }^2\overrightarrow{a}}{\partial t^2}-c_t^2\mathrm{\Delta }\overrightarrow{a}=& \overrightarrow{0}\end{array}$

Die obere Gleichung beschreibt Longitudinalwellen, die sich mit der Geschwindigkeit

$c_l=\sqrt{\frac{\lambda +2\mu }{\varrho }}$

ausbreiten und die untere Transversalwellen, die sich mit der Geschwindigkeit

$c_t=\sqrt{\frac{\mu }{\varrho }}$

fortpflanzen. Wegen $c_l>c_t$ werden Longitudinalwellen als P-Wellen (Primärwellen) und die Transversalwellen als S-Wellen (Sekundärwellen) bezeichnet, denn diese treffen später ein.

Beispiel

Datei:HorizontalSwingingRod.gif

Longitudinalwellen eines elastischen Stabes in den ersten beiden Moden

Bei der Longitudinalwelle des geraden Stabes, der in 1-Richtung liegt (im Bild waagerecht), bewegen sich alle Querschnittsflächen parallel zur 1-Richtung und Schubverzerrungen treten nicht auf.

Unter Vernachlässigung der Schwerebeschleunigung lautet die Navier-Cauchy-Gleichung in 1-Richtung:

$\begin{array}{rl}\rho \frac{{\partial }^2{u}_1}{\partial t^2}=& G[\frac{{\partial }^2{u}_1}{\partial x_1^2}+\frac{{\partial }^2{u}_1}{\partial x_2^2}+\frac{{\partial }^2{u}_1}{\partial x_3^2}\\ & +\frac{1}{1-2\nu }(\frac{{\partial }^2{u}_1}{\partial x_1\partial x_1}+\frac{{\partial }^2{u}_2}{\partial x_1\partial x_2}+\frac{{\partial }^2{u}_3}{\partial x_1\partial x_3})]\end{array}$

Aus der Querdehnung in k-Richtung (k=2, 3)

${\epsilon }_{kk}=-\nu {\epsilon }_{11}\to \frac{\partial u_k}{\partial x_k}=-\nu \frac{\partial u_1}{\partial x_1}\to \frac{{\partial }^2{u}_k}{\partial x_1\partial x_k}=-\nu \frac{{\partial }^2{u}_1}{\partial x_1^2}$

und der Abwesenheit von Schubverzerrungen

${\epsilon }_{1k}=\frac{1}{2}(\frac{\partial u_1}{\partial x_k}+\frac{\partial u_k}{\partial x_1})=0\to \frac{{\partial }^2{u}_1}{\partial x_k^2}=-\frac{{\partial }^2{u}_k}{\partial x_1\partial x_k}=\nu \frac{{\partial }^2{u}_1}{\partial x_1^2}$

folgt mit dem Elastizitätsmodul $E=2G(1+\nu )$:

$\rho \frac{{\partial }^2{u}_1}{\partial t^2}=G[\frac{{\partial }^2{u}_1}{\partial x_1^2}+2\nu \frac{{\partial }^2{u}_1}{\partial x_1^2}+\frac{1}{1-2\nu }\frac{{\partial }^2{u}_1}{\partial x_1^2}-\frac{2\nu }{1-2\nu }\frac{{\partial }^2{u}_1}{\partial x_1^2}]=E\frac{{\partial }^2{u}_1}{\partial x_1^2}$

Senkrecht zur Stabachse werden die Navier-Cauchy-Gleichungen verletzt, weswegen diese Näherung nur für schlanke Stäbe gültig ist.

Mit $u=u_1,x=x_1$ und der Wellenausbreitungsgeschwindigkeit $c=\sqrt{E/\rho }$ ergibt sich die Schwingungsgleichung für den geraden Stab:

$\frac{{\partial }^{2}u}{\partial t^2}=c^2\frac{{\partial }^{2}u}{\partial x^2}$

Der Produktansatz $u(x,t)=aT(t)U(x)+C$ mit freien Parametern a und C, die in der Schwingungsgleichung heraus fallen und der Anpassung an Randbedingungen dienen, sowie zwei noch zu bestimmenden Funktionen T und U ergibt:

$\ddot{T}U=c^{2}T{U}^{″}\to \frac{\ddot{T}}{T}=c^2\frac{U^{″}}{U}$

Der Strich ( )' gibt wie üblich die Ableitung nach der x-Koordnate wieder. Weil die Funktionen auf der linken Seite der letzten Gleichung nur von der Zeit und die auf der rechten Seite nur von der x-Koordinate abhängen, sind die Brüche Konstanten:

$\begin{array}{rlrl}\frac{\ddot{T}}{T}=& -{\omega }^2\to & T(t)=& \sin (\omega t+\alpha )\\ \frac{U^{″}}{U}=& -{\lambda }^2\to & U(x)=& u\sin (\lambda x+\beta )\\ -{\omega }^2=& -c^2{\lambda }^2\to & \omega =& \pm c\lambda \end{array}$.

Die Amplitude der Funktion T und der Faktor a werden der Amplitude u der Funktion U zugeschlagen. Die Amplitude u, die Verschiebung C, die Kreisfrequenzen ω und λ sowie die Phasenwinkel α und β müssen an die Anfangs- und Randbedingungen angepasst werden. Bei fester Einspannung bei x₀ ist

$U(x_0)=u\sin (\lambda x_0+\beta )=0\to \beta =-\lambda x_0$

Andere Werte für β sind zwar möglich, führen aber auf gleichwertige Lösungen und Translationen werden mit dem Parameter C realisiert. An einem freien Ende bei x = x₀ wird die Normalkraft N = E A U′(x₀) vorgegeben, wo der Faktor E der Elastizitätsmodul und A die Querschnittsfläche des Stabes ist. So wird mit der Kraft die Ableitung der Funktion U am freien Ende festgelegt:

$U^{\prime }(x_0)=u\lambda \cos (\lambda x_0+\beta )=\frac{N}{EA}$

Im konkreten Fall hier, wird anfänglich maximale Auslenkung mit

$T(0)=\sin (\alpha )=1\to \alpha =\frac{\pi }{2}\to T(t)=\cos (\omega t)$,

feste Einspannung in C = x = 0, ein unbelastetes freies Ende bei x = L und anfängliche Auslenkung R am freien Ende angenommen:

$\begin{array}{rl}U(0)=& 0\to \beta =0\\ U^{\prime }(L)=& u\lambda \cos (\lambda L)\equiv 0\to \lambda =\frac{(2n-1)\pi }{2L}=\frac{\omega }{c}\\ U(L)=& R=u\sin \left(\frac{(2n-1)\pi }{2}\right)=(-1{)}^{n+1}u\end{array}$

Der Zähler n = 1,2,… beziffert den Schwingungsmode. Die finale Form der Bewegungsfunktion ist somit:

$u(x,t)=(-1{)}^{n+1}R\cos \left(\frac{(2n-1)\pi }{2}\frac{ct}{L}\right)\sin \left(\frac{(2n-1)\pi }{2}\frac{x}{L}\right)$

Das Bild zeigt die mit den Parametern aus der Tabelle berechnete Lösung.

Parameter	Länge L	Endverschiebung R	Mode n	Wellengeschw. c
Einheit	mm	mm	-	mm/s
Wert	100	10	1,2	1

Wegen $\frac{2\cdot 2-1}{2\cdot 1-1}=3$ schwingt die zweite Mode dreimal schneller als die erste.

Siehe auch

• Eulersche Gleichungen (Strömungsmechanik)
• Navier-Stokes-Gleichungen
• Formelsammlung Tensoranalysis
• Formelsammlung Tensoralgebra

Weblinks

• Lamé-Navier’sche Gleichungen. Universität Stuttgart, abgerufen am 30. Oktober 2015. 

Literatur

• Holm Altenbach: Kontinuumsmechanik. Einführung in die materialunabhängigen und materialabhängigen Gleichungen. 2. Auflage. Springer Vieweg, Berlin u. a. 2012, ISBN 978-3-642-24118-5. 
• P. Haupt: Continuum Mechanics and Theory of Materials. Springer, 2000, ISBN 3-540-66114-X. 
• Ralf Greve: Kontinuumsmechanik. Springer, 2003, ISBN 3-540-00760-1. 
• W. H. Müller: Streifzüge durch die Kontinuumstheorie. Springer, 2011, ISBN 978-3-642-19869-4, S. 106 ff. (Hier taucht in einem Beispiel der Name Lamé-Naviersche Gleichungen auf.). 
• Martin H. Sadd: Elasticity - Theory, applications and numerics. Elsevier Butterworth-Heinemann, 2005, ISBN 0-12-605811-3, S. 440 f. (die Winkel φ und ϑ haben dort die Bedeutung von Θ und ϕ). 
• M. E. Gurtin: The Linear Theory of Elasticity. In: S. Flügge (Hrsg.): Handbuch der Physik. Band VI2/a, Bandherausgeber C. Truesdell. Springer, 1972, ISBN 3-540-05535-5.