Cauchy-eulersche Bewegungsgesetze

Die cauchy-eulerschen Bewegungsgesetze von Augustin-Louis Cauchy und Leonhard Euler sind die lokalen Formen des Impuls- und Drallsatzes in der Kontinuumsmechanik. Es sind Bewegungsgleichungen, die, wenn sie lokal, d. h. in jedem Punkt eines Körpers erfüllt sind, sicherstellen, dass die Bewegung des Körpers als Ganzes – inklusive Verformungen – dem Impuls- bzw. Drallsatz gehorcht.

Das erste cauchy-eulersche Bewegungsgesetz korrespondiert mit dem Impulssatz und lautet im geometrisch linearen Fall an einem materiellen Punkt des Körpers:

$ \rho {\vec {a}}=\rho {\vec {k}}+\operatorname {div} \;{\boldsymbol {\sigma }} $

Hier ist ρ die Dichte, $ {\vec {a}} $ die Beschleunigung des materiellen Punktes, $ {\vec {k}} $ die Schwerebeschleunigung, $ {\boldsymbol {\sigma }} $ der cauchysche Spannungstensor und div der Divergenzoperator. Die spezifische Impulsänderung bestimmt sich demnach aus der spezifischen Schwerkraft und dem Antrieb durch einen Spannungsanstieg. Alle Variablen in der Gleichung sind im Allgemeinen sowohl vom Ort als auch von der Zeit abhängig.

Das zweite cauchy-eulersche Bewegungsgesetz entspricht dem lokal formulierten Drallsatz, der sich auf die Forderung nach der Symmetrie des cauchyschen Spannungstensors reduziert:

$ {\boldsymbol {\sigma }}={\boldsymbol {\sigma }}^{\top } $

Das Superskript „⊤“ markiert die Transposition. Die Symmetrie entspricht dem Satz von der Gleichheit der zugeordneten Schubspannungen.^{[L 1]}

σ_ij = σ_ji, i,j = 1,2,3

Bei großen Verschiebungen können beide Bewegungsgesetze in lagrangescher Betrachtungsweise materiell oder eulerscher Betrachtungsweise räumlich formuliert werden. Die Struktur der Gleichungen bleibt dabei erhalten, aber es kommt zu Modifikationen in den Abhängigkeiten oder im Spannungstensor.

Die cauchy-eulerschen Bewegungsgesetze sind die Basis für die eulerschen Gleichungen der Strömungsmechanik, der Navier-Stokes- und der Navier-Cauchy-Gleichungen. Eine der Grundgleichungen der Verschiebungsmethode in der Finite-Elemente-Methode ist das #Prinzip von d’Alembert in der lagrangeschen Fassung, das eine aus den cauchy-eulerschen Gesetzen folgende Aussage ist.

Für Begriffsklärung empfiehlt sich die Lektüre des Artikels zur Kontinuumsmechanik. Die verwendeten Operatoren und Rechenregeln sind in den Formelsammlungen zur Tensoralgebra und Tensoranalysis aufgeführt.

Erstes cauchy-eulersches Bewegungsgesetz

Das erste cauchy-eulersche Bewegungsgesetz folgt aus dem 1687 von Isaac Newton formulierten und nach ihm benannten zweiten newtonschen Gesetz, das dem Impulssatz entspricht, demgemäß die Änderung des Impulses mit der Zeit gleich der auf einen Körper wirkenden äußeren Kräfte ist:

$ {\frac {\mathrm {d} {\vec {p}}}{\mathrm {d} t}}={\vec {F}}_{v}+{\vec {F}}_{a} $

Der Vektor $ {\vec {p}} $ stellt den Impuls dar, dessen zeitliche Änderung sich aus volumenverteilten und oberflächig eingeleiteten Kräften $ {\vec {F}}_{v} $ bzw. $ {\vec {F}}_{a} $ ergibt. In dem die Kontinuumsmechanik den Körper als Punktmenge idealisiert, wird aus der obigen Gleichung eine Integralgleichung, in der der spezifische Impuls, die spezifische Schwerebeschleunigung und die oberflächig wirkenden Kräfte über das Volumen bzw. über die Oberfläche integriert werden. Bei kleinen Verformungen kann das erste cauchy-eulersche Bewegungsgesetz am Volumenelement hergeleitet werden.

Impulssatz am Volumenelement

Der zweidimensionale Fall im ebenen Spannungszustand lässt sich leichter veranschaulichen und soll daher vorangestellt werden. Dazu wird eine ebene Scheibe der Dicke h betrachtet, die durch in der Ebene wirkende Kräfte belastet wird, siehe oberen Bildteil. Aus dieser Scheibe wird gedanklich ein rechteckiges Stück (gelb) herausgeschnitten, parallel zu dessen Kanten ein kartesisches Koordinatensystem definiert wird, in dem es die Breite dx und Höhe dy hat. Nach dem Schnittprinzip entstehen an den Schnittflächen Schnittspannungen, die an die Stelle des weggeschnittenen Teils treten, siehe Schnittreaktion. Bei einem (infinitesimal) kleinen Scheibenelement können die Schnittspannungen als über die Fläche konstant angenommen werden. Die Schnittspannungen $ {\vec {t}}_{x}=\sigma _{xx}{\hat {e}}_{x}+\sigma _{xy}{\hat {e}}_{y} $ treten auf der Oberfläche mit der Normalen in x-Richtung auf und entsprechend operiert $ {\vec {t}}_{y}=\sigma _{yx}{\hat {e}}_{x}+\sigma _{yy}{\hat {e}}_{y} $ auf der Oberfläche mit der Normalen in y-Richtung. In der Komponente $ \sigma _{ij} $ bezieht sich der erste Index also auf die Flächennormale und der zweite Index auf die Wirkrichtung. Nach Voraussetzung gibt es keine Spannungen senkrecht zur Scheibenebene. An den Flächen, deren Normalen in positive Koordinatenrichtung weisen, ist das positive Schnittufer und die Spannungen wirken in positiver Richtung. An den Flächen, deren Normalen in negative Koordinatenrichtung weisen, ist das negative Schnittufer und die Spannungen wirken in negativer Richtung, siehe Bild. Sie sind im Gleichgewicht mit den Schnittspannungen an den benachbarten, weggeschnittenen Teilen des Körpers.

Das zweite newtonsche Gesetz besagt, dass die an dem Scheibenelement angreifenden Spannungen – multipliziert mit ihrer Wirkfläche – das Scheibenelement beschleunigen. An dem Scheibenelement führt das unter Berücksichtigung der Schwerebeschleunigung in x- und y-Richtung auf

$ {\begin{array}{ll}\mathrm {d} m\,a_{x}=\mathrm {d} m\,k_{x}\!\!\!\!\!\!&+\,\sigma _{xx}(x+\mathrm {d} x,y)h\,\mathrm {d} y+\sigma _{yx}(x,y+\mathrm {d} y)h\,\mathrm {d} x\\&-\,\sigma _{xx}(x,y)h\,\mathrm {d} y-\sigma _{yx}(x,y)h\,\mathrm {d} x\\\mathrm {d} m\,a_{y}=\mathrm {d} m\,k_{y}\!\!\!\!\!\!&+\,\sigma _{xy}(x+\mathrm {d} x,y)h\,\mathrm {d} y+\sigma _{yy}(x,y+\mathrm {d} y)h\,\mathrm {d} x\\&-\,\sigma _{xy}(x,y)h\,\mathrm {d} y-\sigma _{yy}(x,y)h\,\mathrm {d} x\end{array}} $

Die Masse $ \mathrm {d} m=\rho h\,\mathrm {d} x\,\mathrm {d} y $ des Scheibenelements ergibt sich aus der Dichte ρ des Materials und dem Volumen $ h\,\mathrm {d} x\,\mathrm {d} y $. Division durch dieses Volumen liefert im Grenzwert $ \mathrm {d} x\rightarrow 0 $ und $ \mathrm {d} y\rightarrow 0 $ den lokalen Impulssatz in x- bzw. y-Richtung:

$ \left.{\begin{array}{rcl}\rho a_{x}&=&\rho k_{x}+{\frac {\partial \sigma _{xx}}{\partial x}}+{\frac {\partial \sigma _{yx}}{\partial y}}\\\rho a_{y}&=&\rho k_{y}+{\frac {\partial \sigma _{xy}}{\partial x}}+{\frac {\partial \sigma _{yy}}{\partial y}}\end{array}}\right\}\leftrightarrow \quad \rho a_{i}=\rho k_{i}+\sum _{j=1}^{2}{\frac {\partial \sigma _{ji}}{\partial x_{j}}}\,,\quad i=1,2 $

wenn – wie üblich – die Koordinaten nach dem Schema x→1, y→2, z→3 durchnummeriert werden. In drei Dimensionen resultieren die gleichen Differentialgleichungen analog, nur wird von eins bis drei summiert:

$ \rho a_{i}=\rho k_{i}+\sum _{j=1}^{3}{\frac {\partial \sigma _{ji}}{\partial x_{j}}}\,,\quad i=1,2,3 $

Multiplikation dieser Gleichungen mit dem Basisvektor ê_i der Standardbasis und Addition der resultierenden drei Gleichungen mündet in der Vektorgleichung

$ {\begin{array}{cccccccl}\underbrace {\displaystyle \sum _{i=1}^{3}\rho a_{i}{\hat {e}}_{i}} \!\!\!\!&=&\!\!\!\!\underbrace {\displaystyle \sum _{i=1}^{3}\rho k_{i}{\hat {e}}_{i}} \!\!\!\!&+&\!\!\!\!\underbrace {\displaystyle \sum _{k=1}^{3}{\hat {e}}_{k}{\frac {\partial }{\partial x_{k}}}} \!\!\!\!&\cdot &\!\!\!\!\underbrace {\displaystyle \sum _{i,j=1}^{3}\sigma _{ji}{\hat {e}}_{j}\otimes {\hat {e}}_{i}} \\\rho {\vec {a}}\!\!\!\!&=&\!\!\!\!\rho {\vec {k}}\!\!\!\!&+&\!\!\!\!\nabla \!\!\!\!&\cdot \!\!\!\!&{\boldsymbol {\sigma }}&=\;\rho {\vec {k}}+\operatorname {div} \;{\boldsymbol {\sigma }}\end{array}} $

Der Nabla-Operator „$ \nabla $“ liefert im Skalarprodukt die Divergenz div des cauchyschen Spannungstensors $ {\boldsymbol {\sigma }} $, der eine Summe von Dyaden ist, die mit dem dyadischen Produkt „$ \otimes $“ der Basisvektoren und seinen Komponenten $ \sigma _{ji} $ gebildet werden. Die Vektorgleichung ist die koordinatenfreie Version des lokalen Impulssatzes, die in beliebigen Koordinaten eines Inertialsystems gilt.

Der Schnittspannungsvektor

$ {\vec {t}}_{x}=\sigma _{xx}{\hat {e}}_{x}+\sigma _{xy}{\hat {e}}_{y}+\sigma _{xz}{\hat {e}}_{z}={\boldsymbol {\sigma }}^{\top }\cdot {\hat {e}}_{x}={\hat {e}}_{x}\cdot {\boldsymbol {\sigma }} $

in der Schnittfläche mit Normalenvektor in x-Richtung ist im cauchyschen Spannungstensor zeilenweise eingetragen, was sinngemäß auch für Schnittspannungsvektoren in y- und z-Richtung gilt.

Impulssatz in lagrangescher Darstellung

In der lagrangeschen Darstellung lautet der globale Impulssatz

$ {\frac {\mathrm {d} }{\mathrm {d} t}}\int _{V}\rho _{0}({\vec {X}})\,{\dot {\vec {\chi }}}({\vec {X}},t)\,\mathrm {d} V=\int _{V}\rho _{0}({\vec {X}})\,{\vec {k}}_{0}({\vec {X}},t)\,\mathrm {d} V+\int _{A}{\vec {t}}_{0}({\vec {X}},t)\,\mathrm {d} A $

die den materiellen Punkten (Partikel) zugeordnete physikalische Größen benutzt, siehe Impulsbilanz. Die Partikel werden durch ihre materiellen Koordinaten $ {\vec {X}}\in V $ in dem Volumen V des Körpers zu einer festgelegten Zeit t₀ im Referenzzustand identifiziert, und auf diese Partikel bezieht sich die Impulsbilanz lokal. Die einem materiellen Punkt zugeordnete Dichte ρ₀ ist auf Grund der Massenbilanz keine Funktion der Zeit. Der aufgesetzte Punkt steht hier wie im Folgenden für die substantielle Ableitung, also für die Zeitableitung bei festgehaltenem Partikel, denn die Gesetze der Mechanik beziehen sich auf die Partikel und nicht auf die Raumpunkte.

Bei den Integralen oben ist das Integrationsgebiet materiell festgelegt, sodass es sich also mit dem Körper mitbewegt, ohne dass neue Partikel zum Gebiet hinzukommen oder wegfallen. Dies wird durch die Großschreibung V bzw. A der Integrationsgebiete symbolisiert.

Weil das Referenzvolumen V somit nicht von der Zeit abhängt, kann die Zeitableitung des Integrals in den Integranden verschoben werden:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \frac{\mathrm{d}}{\mathrm{d}t}\int_V\rho_0(\vec{X})\,\dot{\vec{\chi}}(\vec{X},t)\,\mathrm{d}V =\int_V\rho_0(\vec{X})\,\ddot{\vec{\chi}}(\vec{X},t)\,\mathrm{d}V

Die von außen angreifenden, flächenverteilten Kräfte (Spannungen) Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{t}_0 sind die mit dem Nennspannungstensor N transformierten Normaleneinheitsvektoren Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{N} an der Oberfläche A des Körpers: Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{t}_0 =\mathbf{N}^\top\cdot\vec{N}=\mathbf P\cdot\vec N . Darin ist P = N^⊤ der erste-Piola-Kirchhoff Tensor und ^⊤ bedeutet die Transponierung. Das Oberflächenintegral der Oberflächenspannungen wird mit dem gaußschen Integralsatz in ein Volumenintegral umgewandelt:

$ {\begin{aligned}\int _{A}{\vec {t}}_{0}\,\mathrm {d} A=&\int _{A}\mathbf {N} ^{\top }\cdot {\vec {N}}\,\mathrm {d} A=\int _{V}\nabla _{0}\cdot \mathbf {N} \,\mathrm {d} V\\=&\int _{A}\mathbf {P} \cdot {\vec {N}}\,\mathrm {d} A=\int _{V}\operatorname {DIV} (\mathbf {P} )\,\mathrm {d} V\end{aligned}} $

Der Divergenzoperator DIV wird hier groß geschrieben und der Nabla-Operator wird mit einem Index 𝜵₀ versehen, weil sie die materiellen Ableitungen nach den materiellen Koordinaten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{X} beinhalten. Es gilt für jedes Tensorfeld Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf T(\vec X)\colon\;\operatorname{DIV}\mathbf T=\nabla_0\cdot(\mathbf T^\top) . Die Operatoren in dieser Gleichung sind von den räumlichen Operatoren div bzw. 𝜵 zu unterscheiden, die die räumlichen Ableitungen nach den räumnlichen Koordinaten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{x} ausführen, und die in der eulerschen Darstellung benötigt werden. Mit den vorliegenden Ergebnissen kann die Impulsbilanz als verschwindendes Volumenintegral ausgedrückt werden:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \int_V (\rho_0\ddot{\vec{\chi}}-\rho_0\vec{k}_0-\nabla_0\cdot\mathbf{N})\,\mathrm{d}V =\int_V (\rho_0\ddot{\vec{\chi}}-\rho_0\vec{k}_0-\operatorname{DIV}\mathbf{P}) \,\mathrm{d}V =\vec{0}

Diese Gleichung gilt für jeden Körper und jeden seiner Teilkörper, sodass – Stetigkeit des Integranden vorausgesetzt – auf das erste cauchy-eulersche Bewegungsgesetz in der lagrangeschen Darstellung

$ \rho _{0}({\vec {X}})\,{\ddot {\vec {\chi }}}({\vec {X}},t)=\rho _{0}({\vec {X}})\,{\vec {k}}_{0}({\vec {X}},t)+\nabla _{0}\cdot \mathbf {N} ({\vec {X}},t)=\rho _{0}({\vec {X}})\,{\vec {k}}_{0}({\vec {X}},t)+\operatorname {DIV} \mathbf {P} ({\vec {X}},t) $

geschlossen werden kann. Das Vorkommen der materiellen Koordinaten und des Nennspannungstensors N bzw. des ersten-Piola-Kirchhoff’schen Spannungstensors P an Stelle des cauchyschen Spannungstensors berücksichtigt die Formänderung des bei der Betrachtung am Volumenelement oben herausgeschnittenen Teilkörpers bei großen Deformationen. Bei kleinen Verschiebungen ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf{N}\approx\mathbf{P}\approx\boldsymbol{\sigma} und zwischen den materiellen und räumlichen Koordinaten braucht nicht unterschieden zu werden, wodurch das eingangs angegebene Bewegungsgesetz entsteht.

Impulssatz in eulerscher Darstellung

In der eulerschen Darstellung lautet der globale Impulssatz

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \frac{\mathrm{d}}{\mathrm{d}t} \int_v\rho(\vec{x},t)\vec{v}(\vec{x},t)\,\mathrm{d}v =\int_v\rho(\vec{x},t)\vec{k}(\vec{x},t)\,\mathrm{d}v+\int_a\vec{t}(\vec{x},t)\,\mathrm{d}a

Die räumlichen Punkte werden durch ihre räumlichen Koordinaten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{x}\in v in dem momentanen Volumen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): v zur Zeit t identifiziert, siehe Impulsbilanz. Anders als in der lagrangeschen Darstellung sind die Integrationsgrenzen als Oberflächen des Körpers von der Zeit abhängig, was bei der Berechnung der Impulsänderung zu berücksichtigen ist. Nach dem reynoldsschen Transportsatz gilt:

$ {\begin{aligned}{\frac {\mathrm {d} }{\mathrm {d} t}}\int _{v}\rho ({\vec {x}},t){\vec {v}}({\vec {x}},t)\,\mathrm {d} v=&\int _{v}{\frac {\mathrm {d} }{\mathrm {d} t}}(\rho {\vec {v}})\,\mathrm {d} v+\int _{a}\rho {\vec {v}}({\vec {v}}\cdot \mathrm {d} {\vec {a}})=\int _{v}\left[{\frac {\mathrm {d} }{\mathrm {d} t}}(\rho {\vec {v}})+\operatorname {div} ({\vec {v}})\rho {\vec {v}}\right]\,\mathrm {d} v\\=&\int _{v}{\big (}\underbrace {{\dot {\rho }}{\vec {v}}+\operatorname {div} ({\vec {v}})\rho {\vec {v}}} _{={\vec {0}}}+\rho {\dot {\vec {v}}}{\big )}\,\mathrm {d} v=\int _{v}\rho {\dot {\vec {v}}}\,\mathrm {d} v\end{aligned}} $

Der aufgesetzte Punkt steht für die substantielle Ableitung und in der ersten Zeile wurde das Oberflächenintegral mit dem gaußschen Integralsatz in ein Volumenintegral überführt. Der unterklammerte Term trägt auf Grund der lokalen Massenbilanz Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \dot{\rho}+\rho\operatorname{div}\vec{v}=0 in der eulerschen Darstellung nichts bei.

Das Oberflächenintegral der von außen angreifenden Spannungen wird wie in der lagrangeschen Darstellung mit dem gaußschen Integralsatz in ein Volumenintegral überführt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \int_a\vec{t}\,\mathrm{d}a =\int_a\boldsymbol{\sigma}^\top\cdot\vec{n}\,\mathrm{d}a =\int_v\nabla\cdot\boldsymbol{\sigma}\,\mathrm{d}v =\int_v\operatorname{div}(\boldsymbol{\sigma})\,\mathrm{d}v

denn der cauchysche Spannungstensor σ ist wegen des zweiten cauchy-eulerschen Bewegungsgesetzes unten symmetrisch. Der Nabla-Operator 𝜵 und der Divergenzoperator div beinhalten die räumlichen Ableitungen nach den räumlichen Koordinaten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec x .

Mit den vorliegenden Ergebnissen kann die Impulsbilanz als verschwindendes Volumenintegral ausgedrückt werden:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \int_v(\rho\dot{\vec{v}}-\rho\vec{k}-\nabla\cdot\boldsymbol{\sigma})\,\mathrm{d}v=\vec{0}

Diese Gleichung gilt für jedes Volumen, sodass – Stetigkeit des Integranden vorausgesetzt – das erste cauchy-eulersche Bewegungsgesetz in der eulerschen Darstellung

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \rho(\vec{x},t)\,\dot{\vec{v}}(\vec{x},t) =\rho(\vec{x},t)\,\vec{k}(\vec{x},t)+\nabla\cdot\boldsymbol{\sigma}(\vec{x},t)

abgeleitet werden kann. Hier ist die substantielle Zeitableitung der Geschwindigkeit bei festgehaltenem Partikel $ {\vec {X}} $ zu bilden, das sich zur Zeit t am Ort Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{x}=\vec{\chi}(\vec{X},t) befindet und dort die Geschwindigkeit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{v}(\vec{x},t)=\dot{\vec{\chi}}(\vec{X},t) besitzt:^{[F 1]}

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \dot{\vec{v}}(\vec{x},t) :=&\left.\frac{\mathrm{d}}{\mathrm{d}t}\vec{v}(\vec{\chi}(\vec{X},t),t) \right|_{\vec X\;\text{fest}} =\frac{\partial}{\partial\vec{x}}\vec{v}(\vec{x},t)\cdot\dot{\vec{\chi}}(\vec{X},t) +\frac{\partial}{\partial t}\vec{v}(\vec{x},t) \\=& \operatorname{grad}(\vec{v})\cdot\vec{v}+\frac{\partial\vec{v}}{\partial t} \end{align}

Der räumliche Operator grad berechnet den räumlichen Geschwindigkeitsgradienten mit Ableitungen nach den räumlichen Koordinaten x_1,2,3. Der konvektive Anteil

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \operatorname{grad}(\vec{v})\cdot\vec{v} :=(\nabla\otimes\vec{v})^\top\cdot\vec{v} =\vec{v}\cdot(\nabla\otimes\vec{v}) =(\vec v\cdot\nabla)\vec v

in der substantiellen Beschleunigung berücksichtigt das Hindurchfließen des Materials durch das bei der Betrachtung am Volumenelement oben festgehaltene Volumen V bei großen Verschiebungen. Bei kleinen Verschiebungen kann der quadratische konvektive Anteil vernachlässigt werden, sodass mit

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{a}=\dot{\vec{v}}=\frac{\partial\vec{v}}{\partial t}

das eingangs angegebene Bewegungsgesetz entsteht.

Einfluss von Sprungstellen im Impulssatz

Die verlangte örtliche Stetigkeit der Integranden wird unter realen Verhältnissen verletzt, wenn beispielsweise Dichtesprünge an Materialgrenzen oder Stoßwellen auftreten. Solche flächigen Sprungstellen können jedoch berücksichtigt werden, wenn die Fläche selbst örtlich stetig differenzierbar ist und so in jedem ihrer Punkte einen Normalenvektor besitzt. Die Fläche – im Folgenden Sprungstelle genannt – muss keine materielle Fläche sein, kann sich also mit einer anderen Geschwindigkeit bewegen als die Masse selbst. Durch diese Fläche wird die Masse in zwei Stücke v⁺ und v⁻ geteilt und es wird vereinbart, dass der Normalenvektor der Sprungstelle a_s in Richtung der Sprungstellengeschwindigkeit $ {\vec {v}}_{s} $ und das Volumen v⁺ weise, siehe Bild rechts.

Dann lautet das Reynolds-Transport-Theorem mit Sprungstelle:^{[L 2]}

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \frac{\mathrm{d}}{\mathrm{d}t}\int_{v}\rho\vec{v}\,\mathrm{d}v =& \int_{v}\frac{\partial(\rho\vec{v})}{\partial t}\,\mathrm{d}v +\int_{a}\rho\vec{v}\;(\vec{v}\cdot\,\mathrm{d}\vec{a}) +\int_{a_s}[[\rho\vec{v}\;((\vec{v}_s-\vec{v})\cdot\vec{n})]]\mathrm{d}a_s \\=& \int_{v}\rho\dot{\vec{v}}\,\mathrm{d}v +\int_{a_s}[[\rho\vec{v}\;((\vec{v}_s-\vec{v})\cdot\vec{n})]]\mathrm{d}a_s \end{align}

Der zweite Term mit der Sprungklammer [[...]] kommt neu hinzu. Die Integrale über die von außen angreifenden Kräfte werden getrennt für die Volumina v⁺ und v⁻ berechnet:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \left(\int_{v}\rho\vec k\;\mathrm{d}v+\int_{a}\vec t\,\mathrm{d}a\right)^+ =& \int_{v^+}\rho\vec k\;\mathrm{d}v^++\int_{a^+}\boldsymbol{\sigma}^\top\cdot\vec{n}\,\mathrm{d}a^+ -\int_{a_s}\boldsymbol{\sigma}^\top\cdot\vec{n}\,\mathrm{d}a_s \\ \left(\int_{v}\rho\vec k\;\mathrm{d}v+\int_{a}\vec t\,\mathrm{d}a\right)^- =& \int_{v^-}\rho\vec k\;\mathrm{d}v^-+\int_{a^-}\boldsymbol{\sigma}^\top\cdot\vec{n}\,\mathrm{d}a^- +\int_{a_s}\boldsymbol{\sigma}^\top\cdot\vec{n}\,\mathrm{d}a_s \end{align}

Die Normale soll immer nach außen gerichtet sein und geht daher auf der Sprungstelle einmal mit positivem und einmal mit negativem Vorzeichen ein. Die Vereinigung der Oberflächen a⁺ und a⁻ ergibt die Oberfläche a des gesamten Volumens v, zu dessen Oberfläche die innere Fläche a_s nicht gehört. Die Summe der drei Gleichungen führt nach Umformungen, wie sie oben bereits angegeben wurden, auf

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \int_v(\rho\dot{\vec{v}}-\rho\vec{k}-\nabla\cdot\boldsymbol{\sigma})\,\mathrm{d}v =\int_{a_s}[[\rho\vec{v}\;((\vec{v}-\vec{v}_s)\cdot\vec{n})-\boldsymbol{\sigma}^\top\cdot\vec{n}]]\mathrm{d}a_s

Jenseits der Sprungstelle verschwindet die rechte Seite und die lokale Impulsbilanz ohne Sprungstelle folgt. An der (flächigen) Sprungstelle ist dv=0 und die linke Seite kann vernachlässigt werden, sodass bei Stetigkeit des Integranden mit der Sprungklammer in der Fläche^{[L 3]}

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): [[\rho\vec{v}\;((\vec{v}_s-\vec{v})\cdot\vec{n})+\boldsymbol{\sigma}^\top\cdot\vec{n}]] =[[\rho\vec{v}\otimes(\vec{v}_s-\vec{v})+\boldsymbol{\sigma}^\top]]\cdot\vec{n} =\vec0

abgeleitet werden kann. Wenn die Sprungstelle eine materielle Fläche ist, wie beispielsweise an Materialgrenzen, dann ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{v}_s=\vec{v} und es folgt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): [[\boldsymbol{\sigma}^\top]]\cdot\vec{n}=\vec0 \quad\Leftrightarrow\quad {\boldsymbol{\sigma}^+}^\top\cdot\vec{n} =\vec{t}^+ ={\boldsymbol{\sigma}^-}^\top\cdot\vec{n} =\vec{t}^-

Die Schnittspannungen auf beiden Seiten einer materiellen Sprungstelle müssen gleich sein.^{[L 4]}

Zweites cauchy-eulersches Bewegungsgesetz

Das zweite cauchy-eulersche Bewegungsgesetz folgt aus dem 1754 von Leonhard Euler aufgestellten Drallsatz, nach dem die zeitliche Änderung des Drehimpulses Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{L} gleich der von außen angreifenden Drehmomente ist:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \dot{\vec{L}} =\vec{M}_v+\vec{M}_a

Der Vektor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{M}_v steht für das von volumenverteilten Kräften ausgehende Drehmoment und der Vektor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{M}_a für das oberflächig eingeleitete Moment.

Drallsatz am Volumenelement

Es wird ein belasteter Körper betrachtet, aus dem gedanklich ein würfel­förmiger Teilkörper (im Bild gelb) herausgeschnitten wird, der die Kantenlänge 2L hat und in dessen Massenmittelpunkt ein zu den Würfelkanten parallel ausgerichtetes kartesisches Koordinatensystem gelegt wird. An den Würfelflächen entstehen dem Schnittprinzip zufolge Schnittspannungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{t} , die an die Stelle des weggeschnittenen Teilkörpers treten und die nach dem cauchyschen Fundamentaltheorem die mit dem cauchyschen Spannungstensor transformierten Normalenvektoren an die Schnittfläche sind. Bei infinitesimal kleinem Würfel können die Schnittspannungen als über die Fläche konstant angenommen werden und zu einer Resultierenden aufintegriert werden, die den Würfel aus Symmetriegründen in den Flächenmitten belasten. Für die in der Würfelmitte angreifenden Momente gilt:

• Vom Schwerpunkt des Würfels weist der Vektor $ L{\hat {e}}_{x} $ zur Mitte der Schnittfläche am positiven Schnittufer mit Normale in +x-Richtung und die Schnittspannung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{t}_{+x}=\hat{e}_x\cdot\boldsymbol{\sigma}_{+x} wirkt dort auf der Fläche 4L².
• Das Moment der Schnittspannung am positiven Schnittufer lautet mit dem Kreuzprodukt „ד: Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{M}_{+x}=4L^2(L\hat{e}_x)\times(\hat{e}_x\cdot\boldsymbol{\sigma}_{+x}) .
• Am negativen Schnittufer ist der Hebelarm Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): -L\hat{e}_x und die Schnittspannung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{t}_{-x}=-\hat{e}_x\cdot\boldsymbol{\sigma}_{-x} operiert auf der gleichen Fläche 4L²: Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{M}_{-x}=4 L^2(-L\hat{e}_x)\times(-\hat{e}_x\cdot\boldsymbol{\sigma}_{-x}) .
• Die Momente der Schnittspannungen summieren sich zu $ {\vec {M}}_{x}=4L^{3}{\hat {e}}_{x}\times [{\hat {e}}_{x}\cdot ({\boldsymbol {\sigma }}_{+x}+{\boldsymbol {\sigma }}_{-x})] $.
• In den anderen beiden Raumrichtungen ergibt sich entsprechend Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{M}_y=4L^3\hat{e}_y\times[\hat{e}_y\cdot(\boldsymbol{\sigma}_{+y}+\boldsymbol{\sigma}_{-y})] und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{M}_z=4L^3\hat{e}_z\times[\hat{e}_z\cdot(\boldsymbol{\sigma}_{+z}+\boldsymbol{\sigma}_{-z})] .
• Im infinitesimal kleinen Würfel kann von ortsunabhängiger Dichte ρ und ortsunabhängigem Schwerefeld ausgegangen werden, das daher in der Würfelmitte kein Moment verursacht.

Bei homogener Dichte hat der Würfel die Masse Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): m=\rho(2L)^3=8\rho L^3 und den Trägheitstensor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \boldsymbol{\Theta}=\tfrac{m}{6}(2L)^2\mathbf1=\tfrac{16}{3}\rho L^5\mathbf1 , der proportional zum Einheitstensor 1 ist. Im Drallsatz $ {\vec {M}}={\dot {\vec {L}}}=({\boldsymbol {\Theta }}\cdot {\vec {\omega }}{\dot {)\,}} $ mit der Winkelgeschwindigkeit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec\omega ist die Massenträgheit demnach von fünfter Ordnung in den Abmessungen des ausgeschnittenen Teilkörpers während die Momente nur von dritter Ordnung sind, und das gilt auch bei einem nicht würfelförmigen Quader mit unterschiedlichen Dimensionen in x-, y- und z-Richtung, siehe Trägheitstensor eines Quaders. Bei kleiner werdendem Teilkörper geht die Massenträgheit schneller gegen null als die Momente, woraus das Boltzmann-Axiom resultiert:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec M =4L^3\sum_{k=1}^3 \hat{e}_k\times[\hat{e}_k\cdot (\boldsymbol{\sigma}_{+k}+\boldsymbol{\sigma}_{-k})] =\vec0

Die inneren Kräfte in einem Kontinuum sind momentenfrei. Bei kleiner werdendem Teilkörper wird Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \boldsymbol{\sigma}_{+k}\approx\boldsymbol{\sigma}_{-k}\approx\boldsymbol{\sigma} und somit

$ {\begin{aligned}{\vec {0}}=&\sum _{k=1}^{3}{\hat {e}}_{k}\times ({\hat {e}}_{k}\cdot {\boldsymbol {\sigma }})=\sum _{k=1}^{3}{\hat {e}}_{k}\times \left({\hat {e}}_{k}\cdot \sum _{i,j=1}^{3}\sigma _{ij}{\hat {e}}_{i}\otimes {\hat {e}}_{j}\right)\\=&\sum _{i,j=1}^{3}\sigma _{ij}{\hat {e}}_{i}\times {\hat {e}}_{j}=(\sigma _{12}-\sigma _{21}){\hat {e}}_{3}+(\sigma _{23}-\sigma _{32}){\hat {e}}_{1}+(\sigma _{31}-\sigma _{13}){\hat {e}}_{2}\end{aligned}} $

Hieraus folgt der eingangs erwähnte Satz von der Gleichheit der zugeordneten Schubspannungen.^{[L 1]} Die obige Summe kann mit der Vektorinvariante Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{\mathrm i} von Tensoren koordinatenfrei ausgedrückt werden:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{0} = \sum_{i,j=1}^3\sigma_{ij}\hat{e}_i\times\hat{e}_j = \vec{\mathrm i}\left(\sum_{i,j=1}^3\sigma_{ij}\hat{e}_i\otimes\hat{e}_j\right) = \vec{\mathrm i}(\boldsymbol{\sigma})

Denn in der Vektorinvariante ist das dyadische Produkt ⊗ durch das Kreuzprodukt × ausgetauscht. Nur der schiefsymmetrischer Anteil des Tensors trägt zu seiner Vektorinvariante bei, die hier verschwindet, sodass die Symmetrie des cauchyschen Spannungstensors folgt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \boldsymbol{\sigma}=\boldsymbol{\sigma}^\top

Diese Tensorgleichung, die in beliebigen Vektorraumbasen eines Inertialsystems gilt, ist die koordinatenfreie Version des lokalen Drallsatzes.

Drehimpulssatz in lagrangescher Darstellung

Der Drehimpulssatz lautet in globaler lagrangescher Formulierung:

$ {\frac {\mathrm {d} }{\mathrm {d} t}}\int _{V}({\vec {\chi }}-{\vec {c}})\times \rho _{0}{\dot {\vec {\chi }}}\,\mathrm {d} V=\int _{V}({\vec {\chi }}-{\vec {c}})\times \rho _{0}{\vec {k}}_{0}\,\mathrm {d} V+\int _{A}({\vec {\chi }}-{\vec {c}})\times {\vec {t}}_{0}\,\mathrm {d} A $

worin die physikalischen Größen zumeist sowohl vom Ort Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{X} als auch von der Zeit t abhängen, was hier zwecks kompakter Darstellung unterschlagen wurde. Nur die Dichte ρ₀ ist wegen der Massenerhaltung keine Funktion der Zeit und der beliebige Ortsvektor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{c} ist ebenfalls zeitlich fixiert, siehe Drehimpulsbilanz. Die Zeitableitung des ersten Integrals kann wie beim Impulssatz in den Integranden verschoben werden:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \frac{\mathrm{d}}{\mathrm{d}t} \int_V (\vec{\chi}-\vec{c})\times\rho_0\dot{\vec{\chi}}\,\mathrm{d}V = \int_V[ \underbrace{\dot{\vec{\chi}}\times\rho_0\dot{\vec{\chi}}}_{=\vec{0}} + (\vec{\chi}-\vec{c})\times\rho_0\ddot{\vec{\chi}}]\,\mathrm{d}V = \int_V(\vec{\chi}-\vec{c})\times\rho_0\ddot{\vec{\chi}}\,\mathrm{d}V

Das Oberflächenintegral wird wie gehabt mit dem gaußschen Integralsatz in ein Volumenintegral umgeschrieben:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \int_A (\vec{\chi}-\vec{c})\times\vec{t}_0\,\mathrm{d}A =& \int_A (\vec{\chi}-\vec{c})\times (\mathbf{N}^\top\cdot\vec{N})\,\mathrm{d}A = \int_A [(\vec{\chi}-\vec{c})\times\mathbf{N}^\top]\cdot\vec{N}\,\mathrm{d}A \\=& -\int_A[\mathbf{N}\times(\vec{\chi}-\vec{c})]^\top\cdot\vec{N}\,\mathrm{d}A = -\int_V\nabla_0\cdot[\mathbf{N}\times(\vec{\chi}-\vec{c})]\,\mathrm{d}V \\=& -\int_V\left\{(\nabla_0\cdot\mathbf{N})\times(\vec{\chi}-\vec{c}) -\vec{\mathrm i}\big([\nabla_0\otimes(\vec{\chi}-\vec{c})]^\top\cdot\mathbf{N}\big) \right\}\,\mathrm{d}V \\=& \int_V\left\{(\vec{\chi}-\vec{c})\times(\nabla_0\cdot\mathbf{N}) +\vec{\mathrm i}\big(\mathbf{F}\cdot\mathbf{N}\big) \right\}\,\mathrm{d}V \end{align}

Hier wurde die Produktregel $ \nabla \cdot (\mathbf {T} \times {\vec {f}})=(\nabla \cdot \mathbf {T} )\times {\vec {f}}-{\vec {\mathrm {i} }}\left((\nabla \otimes {\vec {f}})^{\top }\cdot \mathbf {T} \right) $^{[F 2]} und die Definition des Deformationsgradienten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf{F}=\operatorname{GRAD}\vec{\chi}=(\nabla_0\otimes\vec\chi)^\top eingesetzt. Die Operatoren 𝜵₀ und GRAD bilden den materiellen Vektorgradient mit Ableitungen nach den materiellen Koordinaten X_1,2,3, weshalb der Nabla-Operator 𝜵₀ einen Index 0 und der Operator GRAD hier in Abgrenzung zum räumlichen Gradienten grad groß geschrieben wird.

Mit den vorliegenden Ergebnissen kann die Drehimpulsbilanz als verschwindendes Volumenintegral ausgedrückt werden:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \int_V[(\vec{\chi}-\vec{c})\times(\underbrace{\rho_0\ddot{\vec{\chi}} -\rho_0\vec{k}_0-\nabla_0\cdot\mathbf{N}}_{=\vec{0}}) -\vec{\mathrm i}\big(\mathbf{F}\cdot\mathbf{N}\big) ]\,\mathrm{d}V =-\int_V\vec{\mathrm i}\big(\mathbf{F}\cdot\mathbf{N}\big)\,\mathrm{d}V =\vec{0}

Der unterklammerte Term trägt wegen der lokalen Impulsbilanz nichts bei. Das letzte Integral gilt für jeden beliebigen Teilkörper, sodass bei stetigem Integrand Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{\mathrm i}\big(\mathbf{F}\cdot\mathbf{N}\big)=\vec{0} und – wie bei der Herleitung am Volumenelement – die Symmetrie von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf{F\cdot N} abgeleitet werden kann. Die lokale Drehimpulsbilanz in der lagrangeschen Darstellung reduziert sich demnach auf die Forderung

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf{F\cdot N}=(\mathbf{F\cdot N})^\top=\mathbf{N^\top\cdot F^\top}

Multiplikation von links mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf{F}^{-1} und von rechts mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf{F}^{\top-1} ergibt gleichbedeutend:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathbf{F}^{-1}\cdot\mathbf{F\cdot N\cdot F}^{\top-1} =\mathbf{F}^{-1}\cdot\mathbf{N^\top\cdot F^\top\cdot F}^{\top-1} \;\leftrightarrow\; \mathbf{N\cdot F}^{\top-1} =(\mathbf{N\cdot F}^{\top-1})^\top

Der Tensor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \tilde{\mathbf{T}}:=\mathbf{N\cdot F}^{\top-1} ist der zweite piola-kirchhoffsche Spannungstensor, dessen Symmetrie gemäß

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \tilde{\mathbf{T}}=\tilde{\mathbf{T}}^\top

die Erfüllung der Drehimpulsbilanz sicherstellt. Bei kleinen Verschiebungen stimmen der zweite piola-kirchhoffsche und der cauchysche Spannungstensor näherungsweise überein: Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \tilde{\mathbf{T}}\approx\boldsymbol{\sigma} .

Drehimpulssatz in eulerscher Darstellung

In globaler eulerscher Formulierung lautet der Drehimpulssatz:

$ {\frac {\mathrm {d} }{\mathrm {d} t}}\int _{v}({\vec {x}}-{\vec {c}})\times \rho {\vec {v}}\,\mathrm {d} v=\int _{v}({\vec {x}}-{\vec {c}})\times \rho {\vec {k}}\,\mathrm {d} v+\int _{a}({\vec {x}}-{\vec {c}})\times {\vec {t}}\,\mathrm {d} a $

worin die physikalischen Größen zumeist Funktionen sowohl vom Ort Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{x} als auch von der Zeit t sind, was hier zwecks kompakter Darstellung unterschlagen wurde. Der Vektor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{c} ist beliebig und zeitlich fixiert, und die räumlichen Koordinaten Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{x} stellen Integrationsvariable dar, die daher auch nicht von der Zeit abhängen.

Das erste Integral wird wie bei der Impulsbilanz mit dem reynoldsschen Transportsatz berechnet:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \frac{\mathrm{d}}{\mathrm{d}t} \int_v (\vec{x}-\vec{c})\times\rho\vec{v}\,\mathrm{d}v =& \int_v\left\{ \frac{\mathrm{d}}{\mathrm{d}t}[(\vec{x}-\vec{c})\times\rho\vec{v}] +\operatorname{div}(\vec{v}) (\vec{x}-\vec{c})\times\rho\vec{v} \right\}\,\mathrm{d}v \\=& \int_v (\vec{x}-\vec{c})\times [\underbrace{\dot{\rho}\,\vec{v}+\operatorname{div}(\vec{v})\rho\,\vec{v}}_{=\vec0} +\rho\,\dot{\vec{v}}] \,\mathrm{d}v = \int_v (\vec{x}-\vec{c})\times\rho\,\dot{\vec{v}}\,\mathrm{d}v \end{align}

Der unterklammerte Term trägt aufgrund der Massenbilanz nichts bei. Das Oberflächenintegral in der Drehimpulsbilanz wird analog zur lagrangeschen Darstellung mit dem gaußschen Integralsatz in ein Volumenintegral umgeschrieben:^{[F 2]}

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \int_a(\vec{x}-\vec{c})\times\vec{t}\,\mathrm{d}a =& \int_a(\vec{x}-\vec{c})\times(\boldsymbol{\sigma}^\top\cdot\vec{n})\,\mathrm da = \int_a[(\vec{x}-\vec{c})\times\boldsymbol\sigma^\top]\cdot\vec{n}\,\mathrm{d}a \\=& -\int_a[\boldsymbol\sigma\times(\vec{x}-\vec{c})]^\top\cdot\vec{n}\,\mathrm{d}a = -\int_v\nabla\cdot[\boldsymbol\sigma\times(\vec{x}-\vec{c})]\,\mathrm{d}v \\=& -\int_v\left\{(\nabla\cdot\boldsymbol{\sigma})\times(\vec{x}-\vec{c}) -\vec{\operatorname{i}}\left([\nabla\otimes(\vec{x}-\vec{c})]^\top \cdot\boldsymbol{\sigma}\right) \right\}\,\mathrm{d}v \\=& \int_v\left\{(\vec{x}-\vec{c})\times(\nabla\cdot\boldsymbol{\sigma}) +\vec{\operatorname{i}}(\mathbf1\cdot\boldsymbol{\sigma}) \right\}\,\mathrm{d}v \\=& \int_v\left[(\vec{x}-\vec{c})\times(\nabla\cdot\boldsymbol{\sigma}) +\vec{\operatorname{i}}(\boldsymbol{\sigma}) \right]\,\mathrm{d}v \end{align}

Abweichend von der lagrangeschen Darstellung tritt hier der cauchysche Spannungstensor an die Stelle des Nennspannungstensors und wegen $ (\nabla \otimes {\vec {x}})^{\top }=\operatorname {grad} {\vec {x}}=\mathbf {1} $ der Einheitstensor an die Stelle des Deformationsgradienten.

Mit den vorliegenden Ergebnissen kann die Drehimpulsbilanz als verschwindendes Volumenintegral ausgedrückt werden:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \int_v[ (\vec{x}-\vec{c})\times(\underbrace{\rho\dot{\vec{v}} -\rho\vec{k}-\nabla\cdot\boldsymbol{\sigma}}_{=\vec{0}}) -\vec{\operatorname{i}}(\boldsymbol{\sigma}) ]\,\mathrm{d}v =-\int_v\vec{\operatorname{i}}(\boldsymbol{\sigma})\,\mathrm{d}v =\vec{0}

Der unterklammerte Term trägt wegen der lokalen Impulsbilanz nichts bei und das letzte Integral gilt für jedes beliebige Volumen, sodass bei stetigem Integrand auf Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec{\operatorname{i}}(\boldsymbol{\sigma})=\vec{0} geschlossen werden kann. Analog zur lagrangeschen Darstellung reduziert sich die Drehimpulsbilanz in eulerscher Darstellung auf die Forderung nach der Symmetrie des cauchyschen Spannungstensors:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \boldsymbol{\sigma}=\boldsymbol{\sigma}^\top

Einfluss von Sprungstellen im Drehimpulssatz

Analog zum ersten cauchy-eulerschen Bewegungsgesetz lautet das Reynolds-Transport-Theorem mit Sprungstelle hier:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \frac{\mathrm{d}}{\mathrm{d}t}\int_{v}(\vec x-\vec{c})\times\rho\vec{v}\,\mathrm{d}v = \int_{v}(\vec x-\vec{c})\times\rho\dot{\vec{v}}\,\mathrm{d}v +\int_{a_s}[[(\vec x-\vec{c})\times\rho\vec{v}\;((\vec{v}_s-\vec{v})\cdot\vec{n})]]\mathrm{d}a_s

Der zweite Term mit der Sprungklammer [[...]] kommt neu hinzu. Die Integrale über die von außen angreifenden Kräfte werden getrennt für v⁺ und v⁻ berechnet:

$ {\begin{array}{l}\displaystyle \left(\int _{v}({\vec {x}}-{\vec {c}})\times \rho {\vec {k}}\;\mathrm {d} v+\int _{a}({\vec {x}}-{\vec {c}})\times {\vec {t}}\,\mathrm {d} a\right)^{+}=\ldots \\\displaystyle \ldots =\int _{v^{+}}({\vec {x}}-{\vec {c}})\times \rho {\vec {k}}\;\mathrm {d} v^{+}+\int _{a^{+}}({\vec {x}}-{\vec {c}})\times {\boldsymbol {\sigma }}^{\top }\cdot {\vec {n}}\,\mathrm {d} a^{+}-\int _{a_{s}}({\vec {x}}-{\vec {c}})\times {\boldsymbol {\sigma }}^{\top }\cdot {\vec {n}}\,\mathrm {d} a_{s}\\\displaystyle \left(\int _{v}({\vec {x}}-{\vec {c}})\times \rho {\vec {k}}\;\mathrm {d} v+\int _{a}({\vec {x}}-{\vec {c}})\times {\vec {t}}\,\mathrm {d} a\right)^{-}=\ldots \\\displaystyle \ldots =\int _{v^{-}}({\vec {x}}-{\vec {c}})\times \rho {\vec {k}}\;\mathrm {d} v^{-}+\int _{a^{-}}({\vec {x}}-{\vec {c}})\times {\boldsymbol {\sigma }}^{\top }\cdot {\vec {n}}\,\mathrm {d} a^{-}+\int _{a_{s}}({\vec {x}}-{\vec {c}})\times {\boldsymbol {\sigma }}^{\top }\cdot {\vec {n}}\,\mathrm {d} a_{s}\end{array}} $

Die Summe der drei Gleichungen führt nach Umformungen, wie sie oben bereits angegeben wurden, auf

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{array}{l} \displaystyle\int_v[ (\vec{x}-\vec{c})\times(\rho\dot{\vec{v}}-\rho\vec{k}-\nabla\cdot\boldsymbol{\sigma}) -\vec{\operatorname{i}}(\boldsymbol{\sigma}) ]\,\mathrm{d}v =\ldots \\ \qquad\qquad\ldots= \int_{a_s}[[(\vec x-\vec{c})\times\rho\vec{v}\;((\vec{v}-\vec{v}_s)\cdot\vec{n}) -(\vec x-\vec{c})\times\boldsymbol{\sigma}^\top\cdot\vec{n}]]\mathrm{d}a_s \end{array}

Jenseits der Sprungstelle verschwindet die rechte Seite und die Symmetrie des Spannungstensors folgt wie oben. An der (flächigen) Sprungstelle ist dv=0 und die linke Seite kann vernachlässigt werden, sodass bei Stetigkeit des Integranden mit der Sprungklammer in der Fläche

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): (\vec x-\vec{c})\times[[\rho\vec{v}\;((\vec{v}_s-\vec{v})\cdot\vec{n}) +\boldsymbol{\sigma}^\top\cdot\vec{n}]] =(\vec x-\vec{c})\times[[\rho\vec{v}\otimes(\vec{v}_s-\vec{v}) +\boldsymbol{\sigma}^\top]]\cdot\vec{n} =\vec0

abgeleitet werden kann, was wegen der Sprungbedingung im ersten cauchy-eulerschen Bewegungsgesetz identisch erfüllt ist.

Folgerungen aus den Bewegungsgesetzen

Aus den Bewegungsgesetzen können weitere, materialunabhängige, zu Prinzipien äquivalente Gleichungen gefolgert werden. Das erste cauchy-eulersche Bewegungsgesetz lautet:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \text{materiell:}\quad & \rho_0\,\ddot{\vec\chi}\!\!\!\!\!\!\!\!\!\!&=&\;\rho_0\,\vec{k}_0+\nabla_0\cdot\mathbf{N} \\ \text{räumlich:}\quad & \rho\,\dot{\vec v}\!\!\!\!\!\!\!\!\!\!&=&\;\rho\,\vec{k}+\nabla\cdot\boldsymbol{\sigma} \end{align}

Diese Gleichungen werden mit einem Vektorfeld Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec q skalar multipliziert, über das Volumen des Körpers integriert und umgeformt. Es entsteht:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \text{materiell:} & \int_V\rho_0\,\ddot{\vec{\chi}}\cdot\vec{q}\,\mathrm{d}V +\int_V\tilde{\mathbf{T}}: \operatorname{sym}(\mathbf{F}^\top\cdot\operatorname{GRAD}\vec q)\,\mathrm{d}V=\ldots \\&\qquad\ldots= \int_V\rho_0\,\vec{k}_0\cdot\vec{q}\,\mathrm{d}V +\int_A\vec{t}_0\cdot\vec{q}\,\mathrm{d}A \\ \text{räumlich:}& \int_v\rho\,\dot{\vec{v}}\cdot\vec{q}\,\mathrm{d}v +\int_v\boldsymbol{\sigma}:\operatorname{sym(grad}(\vec q))\,\mathrm{d}v =\ldots \\&\qquad\ldots= \int_v\rho\,\vec{k}\cdot\vec{q}\,\mathrm{d}v+\int_a\vec{t}\cdot\vec{q}\,\mathrm{d}a \end{align}

Je nach Vektorfeld $ {\vec {q}} $ ergeben sich verschiedene Aussagen.

Prinzip von d’Alembert

Das Prinzip von d’Alembert hat eine grundlegende Bedeutung für die Lösung von Anfangsrandwertaufgaben der Kontinuumsmechanik, insbesondere der Verschiebungsmethode in der Finite-Elemente-Methode. Für das Vektorfeld Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec q werden virtuelle Verschiebungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \delta\vec{u} eingesetzt, die vom Verschiebungsfeld $ {\vec {u}} $ unabhängige, gedachte, weitgehend beliebige, differenzielle Verschiebungen sind und die mit den geometrischen Bindungen des Körpers verträglich sind. Die virtuellen Verschiebungen müssen verschwinden, wo immer Verschiebungsrandbedingungen des Körpers vorgegeben sind. Sei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {A}^{u} der Teil der Oberfläche Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): A des Körpers, auf dem Verschiebungsrandbedingungen erklärt sind. Für ein Vektorfeld der virtuellen Verschiebungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \delta\vec{u}_0 ist dann

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \text{materiell}\quad&\delta\vec{u}_0(\vec{X})=\vec{0}\quad\text{für alle}\;\vec{X}{\in}{A}^{u} \\ \text{räumlich}\quad&\delta\vec{u}(\vec{x})=\vec{0}\quad\text{für alle}\;\vec{x}{\in}{a}^{u} \end{align}

zu fordern. Auf $ {A}^{u} $ können dann keine Oberflächenspannungen vorgegeben werden. Deshalb bezeichnet Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {A}^{\sigma}=A\setminus A^{u} den Teil der Oberfläche des Körpers, auf dem Oberflächenspannungen wirken (können) was entsprechend auch in der räumlichen Formulierung definiert wird. So entsteht:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{array}{llcl} \text{materiell:} &\displaystyle \int_V\rho_0\,\ddot{\vec{\chi}}\cdot\delta\vec{u}_0\,\mathrm{d}V +\int_V\tilde{\mathbf{T}}:\delta\mathbf{E}\,\mathrm{d}V \!\!\!\!\!\!&=&\!\!\!\!\!\!\displaystyle \int_V\rho_0\,\vec{k}_0\cdot\delta\vec{u}_0\,\mathrm{d}V +\int_{A^\sigma}\vec{t}_0\cdot\delta\vec{u}_0\,\mathrm{d}A \\&\text{für alle}\;\delta\vec{u}_0\in\mathcal{V}_0 \\ \text{räumlich:}&\displaystyle \int_v\rho\,\dot{\vec{v}}\cdot\delta\vec{u}\,\mathrm{d}v +\int_v\boldsymbol{\sigma}:\delta\boldsymbol{\varepsilon}\,\mathrm{d}v \!\!\!\!\!\!&=&\!\!\!\!\!\!\displaystyle \int_v\rho\,\vec{k}\cdot\delta\vec{u}\,\mathrm{d}v +\int_{a^\sigma}\vec{t}\cdot\delta\vec{u}\,\mathrm{d}a \\&\text{für alle}\;\delta\vec{u}\in\mathcal{V} \end{array}

Die Menge Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \mathcal{V}_{(0)} enthält alle zulässigen, materiellen bzw. räumlichen, virtuellen Verschiebungsfelder. Auf der linken Seite steht die virtuelle Arbeit der Trägheitskräfte und die virtuelle Deformationsarbeit und auf der rechten Seite die virtuelle Arbeit der äußeren Kräfte (volumen- und oberflächenverteilt.)

In der materiellen Darstellung stehen die virtuelle Verzerrungen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \delta\mathbf{E} für die Variation des green-lagrangeschen Verzerrungstensors:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \delta\mathbf{E} :=\frac{1}{2}(\mathbf{F}^\top\cdot\delta\mathbf{F}+\delta\mathbf{F}^\top\cdot\mathbf{F}) =\operatorname{sym}(\mathbf{F}^\top\cdot\delta\mathbf{F})

mit dem virtuellen Deformationsgradient $ \delta \mathbf {F} :=\operatorname {GRAD} \delta {\vec {u}} $. In der räumlichen Darstellung bildet sich der virtuelle Verzerrungstensor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \delta\boldsymbol\varepsilon aus dem räumlichen virtuellen Verschiebungsgradient Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \delta\mathbf{h}:=\operatorname{grad}\delta\vec u :

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \delta\boldsymbol{\varepsilon}:=\operatorname{sym\,grad}\delta\vec{u} =\frac{1}{2}(\delta\mathbf{h}+\delta\mathbf{h}^\top)

Bilanz der mechanischen Energie

Wenn für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec q das Geschwindigkeitsfeld eingesetzt wird, folgt die Bilanz der mechanischen Energie:

$ {\begin{array}{llcl}{\text{materiell:}}&{\frac {\mathrm {d} }{\mathrm {d} t}}\int _{V}{\frac {\rho _{0}}{2}}\,{\dot {\vec {\chi }}}\cdot {\dot {\vec {\chi }}}\,\mathrm {d} V&=&\int _{V}\rho _{0}\,{\vec {k}}_{0}\cdot {\dot {\vec {\chi }}}\,\mathrm {d} V+\int _{A}{\vec {t}}_{0}\cdot {\dot {\vec {\chi }}}\,\mathrm {d} A-\int _{V}{\tilde {\mathbf {T} }}:{\dot {\mathbf {E} }}\,\mathrm {d} V\\{\text{räumlich:}}&{\frac {\mathrm {d} }{\mathrm {d} t}}\int _{v}{\frac {\rho }{2}}{\vec {v}}\cdot {\vec {v}}\,\mathrm {d} v&=&\int _{v}\rho \,{\vec {k}}\cdot {\vec {v}}\,\mathrm {d} v+\int _{a}{\vec {t}}\cdot {\vec {v}}\,\mathrm {d} a-\int _{v}{\boldsymbol {\sigma }}:\mathbf {d} \,\mathrm {d} v\end{array}} $

Auf der linken Seite steht die zeitliche Änderung der kinetischen Energie und auf der rechten Seite steht die Leistung der äußeren Kräfte (volumen- und flächenverteilt) abzüglich der Verformungsleistung. Dieser Satz wird auch Arbeitssatz^{[L 5]} oder „Satz von der geleisteten Arbeit“ (englisch Theorem of work expended^{[L 6]}) genannt.

Energieerhaltungssatz

In einem konservativen System gibt es eine skalarwertige Funktion W_a, die potentielle Energie, deren negative Zeitableitung gemäß

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): -\frac{\mathrm{d}}{\mathrm{d}t}W_a :=\int_V\rho_0\,\vec{k}_0\cdot\dot{\vec{\chi}}\,\mathrm{d}V +\int_A\vec{t}_0\cdot\dot{\vec{\chi}}\,\mathrm{d}A = \int_v\rho\,\vec{k}\cdot\vec{v}\,\mathrm{d}v+\int_a\vec{t}\cdot\vec{v}\,\mathrm{d}a

die Leistung der äußeren Kräfte ist, und eine Formänderungsenergie W_i, deren Zeitableitung

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \frac{\mathrm{d}}{\mathrm{d}t}W_i :=\int_V\tilde{\mathbf{T}}:\dot{\mathbf{E}}\,\mathrm{d}V =\int_v\boldsymbol{\sigma}:\mathbf{d}\,\mathrm{d}v

die Verformungsleistung ist. Mit der Abkürzung

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): K :=\int_V\frac{\rho_0}{2}\,\dot{\vec{\chi}}\cdot\dot{\vec{\chi}}\,\mathrm{d}V =\int_v\frac{\rho}{2}\vec{v}\cdot\vec{v}\,\mathrm{d}v

für die kinetische Energie schreibt sich die Bilanz der mechanischen Energie:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \frac{\mathrm{d}}{\mathrm{d}t}K+\frac{\mathrm{d}}{\mathrm{d}t}W_i+\frac{\mathrm{d}}{\mathrm{d}t}W_a =\frac{\mathrm{d}}{\mathrm{d}t}(K+W_i+W_a)=0 \quad\Leftrightarrow K+W_i+W_a=E

Die mechanische Gesamtenergie E, bestehend aus der kinetischen Energie, der Formänderungsenergie und der potentiellen Energie, ist mithin in einem konservativen System zeitlich konstant, was als Energieerhaltungssatz bekannt ist.

Satz von Clapeyron

Wird bei kleinen Verformungen, linearer Elastizität und im statischen Fall für das Vektorfeld Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \vec q das Verschiebungsfeld $ {\vec {u}}({\vec {X}},t):={\vec {\chi }}({\vec {X}},t)-{\vec {X}} $ eingesetzt, dann ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \operatorname{grad}\vec q=\operatorname{grad}\vec u=\mathbf{H} der Verschiebungsgradient und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \|\mathbf{H}\|\ll1 , sodass alle Terme, die Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \|\mathbf{H}\| in höherer Ordnung als eins enthalten, vernachlässigt werden können. Es folgt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \begin{align} \operatorname{GRAD}\vec u=\mathbf{H}=\operatorname{GRAD}\vec{\chi}-\mathbf1 =\mathbf{F-1} \\\rightarrow\quad \operatorname{sym}(\mathbf{F}^\top\cdot\operatorname{GRAD}\vec u) = \operatorname{sym}(\mathbf{(1+H^\top)\cdot H}) \approx \operatorname{sym}(\mathbf{H}) = \frac12(\mathbf{H+H^\top})=\boldsymbol{\varepsilon} \end{align}

Der symmetrische Anteil des Verschiebungsgradienten ist der linearisierte Verzerrungstensor. Der zweite piola-kirchhoffsche Spannungstensor geht bei kleinen Verformungen in den cauchyschen Spannungstensor über und es resultiert der Arbeitssatz^{[L 7]}

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \int_v\boldsymbol{\sigma}:\boldsymbol{\varepsilon}\,\mathrm{d}v =\int_v\rho\,\vec{k}\cdot\vec{u}\,\mathrm{d}v+\int_a\vec{t}\cdot\vec{u}\,\mathrm{d}a

Der Integrand auf der linken Seite ist das Doppelte der Formänderungsenergie $ {\boldsymbol {\sigma }}:{\boldsymbol {\varepsilon }}=2w $ und es entsteht der Satz von Clapeyron^{[L 8]}

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): 2\int_v w\,\mathrm{d}v =\int_v\rho\,\vec{k}\cdot\vec{u}\,\mathrm{d}v+\int_a\vec{t}\cdot\vec{u}\,\mathrm{d}a

Fußnoten

Literatur

• Holm Altenbach: Kontinuumsmechanik. Springer, 2012, ISBN 978-3-642-24118-5. 
• Ralf Greve: Kontinuumsmechanik. Ein Grundkurs für Ingenieure und Physiker. Springer, Berlin u. a. 2003, ISBN 3-540-00760-1. 
• Morton E. Gurtin: The Linear Theory of Elasticity. In: C. Truesdell (Hrsg.): Festkörpermechanik : Teil 2 (= S. Flügge [Hrsg.]: Handbuch der Physik. Band 6a, Teilbd. 2). Springer, Berlin/Heidelberg 1972, ISBN 3-540-05535-5, S. 1–295, doi:10.1007/978-3-662-39776-3_1 (DOI der englischen Ausgabe). 
• P. Haupt: Continuum Mechanics and Theory of Materials. Springer, 2002, ISBN 978-3-642-07718-0, doi:10.1007/978-3-662-04775-0. 

Einzelnachweise