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:

Hier ist ρ die Dichte, die Beschleunigung des materiellen Punktes, die Schwerebeschleunigung, 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:

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

.

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

Bearbeiten

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:

 

Der Vektor   stellt den Impuls dar, dessen zeitliche Änderung sich aus volumenverteilten und oberflächig eingeleiteten Kräften   bzw.   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

Bearbeiten
 
Spannungen an einem freigeschnittenen Scheibenelement

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   treten auf der Oberfläche mit der Normalen in x-Richtung auf und entsprechend operiert   auf der Oberfläche mit der Normalen in y-Richtung. In der Komponente   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

 

Die Masse   des Scheibenelements ergibt sich aus der Dichte ρ des Materials und dem Volumen  . Division durch dieses Volumen liefert im Grenzwert   und   den lokalen Impulssatz in x- bzw. y-Richtung:

 

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:

 

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

 

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

Der Schnittspannungsvektor

 

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

Bearbeiten

In der lagrangeschen Darstellung lautet der globale Impulssatz

 

die den materiellen Punkten (Partikel) zugeordnete physikalische Größen benutzt, siehe Impulsbilanz. Die Partikel werden durch ihre materiellen Koordinaten   in dem Volumen V des Körpers zu einer festgelegten Zeit t0 im Referenzzustand identifiziert, und auf diese Partikel bezieht sich die Impulsbilanz lokal. Die einem materiellen Punkt zugeordnete Dichte ρ0 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:

 

Die von außen angreifenden, flächenverteilten Kräfte (Spannungen)   sind die mit dem Nennspannungstensor N transformierten Normaleneinheitsvektoren   an der Oberfläche A des Körpers:  . 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:

 

Der Divergenzoperator DIV wird hier groß geschrieben und der Nabla-Operator wird mit einem Index 𝜵0 versehen, weil sie die materiellen Ableitungen nach den materiellen Koordinaten   beinhalten. Es gilt für jedes Tensorfeld  . Die Operatoren in dieser Gleichung sind von den räumlichen Operatoren div bzw. 𝜵 zu unterscheiden, die die räumlichen Ableitungen nach den räumlichen Koordinaten   ausführen, und die in der eulerschen Darstellung benötigt werden. Mit den vorliegenden Ergebnissen kann die Impulsbilanz als verschwindendes Volumenintegral ausgedrückt werden:

 

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

 

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   und zwischen den materiellen und räumlichen Koordinaten braucht nicht unterschieden zu werden, wodurch das eingangs angegebene Bewegungsgesetz entsteht.

Impulssatz in eulerscher Darstellung

Bearbeiten

In der eulerschen Darstellung lautet der globale Impulssatz

 

Die räumlichen Punkte werden durch ihre räumlichen Koordinaten   in dem momentanen Volumen   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:

 

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   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:

 

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  .

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

 

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

 

abgeleitet werden kann. Hier ist die substantielle Zeitableitung der Geschwindigkeit bei festgehaltenem Partikel   zu bilden, das sich zur Zeit t am Ort   befindet und dort die Geschwindigkeit   besitzt:[F 1]

 

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

 

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

 

das eingangs angegebene Bewegungsgesetz entsteht.

Einfluss von Sprungstellen im Impulssatz

Bearbeiten
 
Eine Sprungstelle auf der Fläche as trennt zwei Raumbereiche v+ und v

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 as in Richtung der Sprungstellengeschwindigkeit   und das Volumen v+ weise, siehe Bild rechts.

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

 

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:

 

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 as nicht gehört. Die Summe der drei Gleichungen führt nach Umformungen, wie sie oben bereits angegeben wurden, auf

 

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]

 

abgeleitet werden kann. Wenn die Sprungstelle eine materielle Fläche ist, wie beispielsweise an Materialgrenzen, dann ist   und es folgt:

 

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

Zweites cauchy-eulersches Bewegungsgesetz

Bearbeiten

Das zweite cauchy-eulersche Bewegungsgesetz folgt aus dem 1754 von Leonhard Euler aufgestellten Drallsatz, nach dem die zeitliche Änderung des Drehimpulses   gleich der von außen angreifenden Drehmomente ist:

 

Der Vektor   steht für das von volumenverteilten Kräften ausgehende Drehmoment und der Vektor   für das oberflächig eingeleitete Moment.

Drallsatz am Volumenelement

Bearbeiten
 
Schnittspannungen an einem würfelförmigen Teilkörper

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  , 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   zur Mitte der Schnittfläche am positiven Schnittufer mit Normale in +x-Richtung und die Schnittspannung   wirkt dort auf der Fläche 4L².
  • Das Moment der Schnittspannung am positiven Schnittufer lautet mit dem Kreuzprodukt „ד:  .
  • Am negativen Schnittufer ist der Hebelarm   und die Schnittspannung   operiert auf der gleichen Fläche 4L²:  .
  • Die Momente der Schnittspannungen summieren sich zu  .
  • In den anderen beiden Raumrichtungen ergibt sich entsprechend   und  .
  • 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   und den Trägheitstensor  , der proportional zum Einheitstensor 1 ist. Im Drallsatz   mit der Winkelgeschwindigkeit   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:

 

Die inneren Kräfte in einem Kontinuum sind momentenfrei. Bei kleiner werdendem Teilkörper wird   und somit

 

Hieraus folgt der eingangs erwähnte Satz von der Gleichheit der zugeordneten Schubspannungen.[L 1] Die obige Summe kann mit der Vektorinvariante   von Tensoren koordinatenfrei ausgedrückt werden:

 

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:

 

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

Drehimpulssatz in lagrangescher Darstellung

Bearbeiten

Der Drehimpulssatz lautet in globaler lagrangescher Formulierung:

 

worin die physikalischen Größen zumeist sowohl vom Ort   als auch von der Zeit t abhängen, was hier zwecks kompakter Darstellung unterschlagen wurde. Nur die Dichte ρ0 ist wegen der Massenerhaltung keine Funktion der Zeit und der beliebige Ortsvektor   ist ebenfalls zeitlich fixiert, siehe Drehimpulsbilanz. Die Zeitableitung des ersten Integrals kann wie beim Impulssatz in den Integranden verschoben werden:

 

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

 

Hier wurde die Produktregel  [F 2] und die Definition des Deformationsgradienten   eingesetzt. Die Operatoren 𝜵0 und GRAD bilden den materiellen Vektorgradient mit Ableitungen nach den materiellen Koordinaten X1,2,3, weshalb der Nabla-Operator 𝜵0 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:

 

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   und – wie bei der Herleitung am Volumenelement – die Symmetrie von   abgeleitet werden kann. Die lokale Drehimpulsbilanz in der lagrangeschen Darstellung reduziert sich demnach auf die Forderung

 

Multiplikation von links mit   und von rechts mit   ergibt gleichbedeutend:

 

Der Tensor   ist der zweite piola-kirchhoffsche Spannungstensor, dessen Symmetrie gemäß

 

die Erfüllung der Drehimpulsbilanz sicherstellt. Bei kleinen Verschiebungen stimmen der zweite piola-kirchhoffsche und der cauchysche Spannungstensor näherungsweise überein:  .

Drehimpulssatz in eulerscher Darstellung

Bearbeiten

In globaler eulerscher Formulierung lautet der Drehimpulssatz:

 

worin die physikalischen Größen zumeist Funktionen sowohl vom Ort   als auch von der Zeit t sind, was hier zwecks kompakter Darstellung unterschlagen wurde. Der Vektor   ist beliebig und zeitlich fixiert, und die räumlichen Koordinaten   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:

 

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]

 

Abweichend von der lagrangeschen Darstellung tritt hier der cauchysche Spannungstensor an die Stelle des Nennspannungstensors und wegen   der Einheitstensor an die Stelle des Deformationsgradienten.

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

 

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   geschlossen werden kann. Analog zur lagrangeschen Darstellung reduziert sich die Drehimpulsbilanz in eulerscher Darstellung auf die Forderung nach der Symmetrie des cauchyschen Spannungstensors:

 

Einfluss von Sprungstellen im Drehimpulssatz

Bearbeiten

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

 

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:

 

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

 

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

 

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

Folgerungen aus den Bewegungsgesetzen

Bearbeiten

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

 

Diese Gleichungen werden mit einem Vektorfeld   skalar multipliziert, über das Volumen des Körpers integriert und umgeformt. Es entsteht:

 

Je nach Vektorfeld   ergeben sich verschiedene Aussagen.

Beweis
Skalare Multiplikation des ersten cauchy-eulerschen Bewegungsgesetzes mit dem Vektorfeld   und Integration über das Volumen des Körpers liefert:

 

Der letzte Term auf der rechten Seite wird mit der Produktregel umgeformt:

 

In der materiellen Form wird noch der Deformationsgradient einmultipliziert:

 

Im letzten Schritt wurde   ausgenutzt, und dass im Skalarprodukt mit einem symmetrischen Tensor   nur die symmetrischen Anteile sym(·) etwas beitragen, was auch in der räumlichen Formulierung ausgenutzt wird:  

Das Volumenintegral des Divergenzterms wird mit dem gaußschen Integralsatz in ein Oberflächenintegral umgewandelt:

 

Zusammenführung dieser Ergebnisse resultiert in den angegebenen Gleichungen.

Prinzip von d’Alembert

Bearbeiten

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   werden virtuelle Verschiebungen   eingesetzt, die vom Verschiebungsfeld   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   der Teil der Oberfläche   des Körpers, auf dem Verschiebungsrandbedingungen erklärt sind. Für ein Vektorfeld der virtuellen Verschiebungen   ist dann

 

zu fordern. Auf   können dann keine Oberflächenspannungen vorgegeben werden. Deshalb bezeichnet   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:

 

Die Menge   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   für die Variation des green-lagrangeschen Verzerrungstensors:

 

mit dem virtuellen Deformationsgradient  . In der räumlichen Darstellung bildet sich der virtuelle Verzerrungstensor   aus dem räumlichen virtuellen Verschiebungsgradient  :

 

Bilanz der mechanischen Energie

Bearbeiten

Wenn für   das Geschwindigkeitsfeld eingesetzt wird, folgt die Bilanz der mechanischen Energie:

 

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.

Beweis
Die Zeitableitung der kinetischen Energie ist gleich der Leistung der Trägheitskräfte
 

was die ersten Terme auf den linken Seiten begründet. In der räumlichen Formulierung wurde der reynoldssche Transportsatz und die Massenbilanz angewendet.

Der materielle Gradient der Geschwindigkeit ist die Zeitableitung   des Deformationsgradienten und der symmetrische Anteil des räumlichen Geschwindigkeitsgradienten   ist der Verzerrungsgeschwindigkeitstensor d. Damit schreiben sich die Verformungsleistungen:

 

Die Leistungen der äußeren Kräfte ergeben sich durch Ersetzung des Vektors   durch den Geschwindigkeitsvektor.

Energieerhaltungssatz

Bearbeiten

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

 

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

 

die Verformungsleistung ist. Mit der Abkürzung

 

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

 

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

Bearbeiten

Wird bei kleinen Verformungen, linearer Elastizität und im statischen Fall für das Vektorfeld   das Verschiebungsfeld   eingesetzt, dann ist   der Verschiebungsgradient und  , sodass alle Terme, die   in höherer Ordnung als eins enthalten, vernachlässigt werden können. Es folgt:

 

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]

 

Der Integrand auf der linken Seite ist das Doppelte der Formänderungsenergie   und es entsteht der Satz von Clapeyron[L 8]

 

Fußnoten

Bearbeiten
  1. Die Fréchet-Ableitung einer Funktion   nach   ist der beschränkte lineare Operator   der – sofern er existiert – in alle Richtungen   dem Gâteaux-Differential entspricht, also   gilt. Darin ist   skalar-, vektor- oder tensorwertig aber   und   gleichartig. Dann wird auch   geschrieben.
  2. a b Beweis der Produktregel in kartesischen Koordinaten mit einsteinscher Summenkonvention:
     
    Siehe auch Vektorinvariante und Formelsammlung Tensoralgebra#Kreuzprodukt eines Vektors mit einem Tensor

Literatur

Bearbeiten

Einzelnachweise

Bearbeiten
  1. a b István Szabó: Geschichte der mechanischen Prinzipien. Springer, 2013, ISBN 978-3-0348-5301-9 (eingeschränkte Vorschau in der Google-Buchsuche [abgerufen am 2. Mai 2021]).
  2. P. Haupt: Continuum Mechanics and Theory of Materials. Springer, 2002, ISBN 978-3-642-07718-0, S. 141, doi:10.1007/978-3-662-04775-0.
  3. P. Haupt: Continuum Mechanics and Theory of Materials. Springer, 2002, ISBN 978-3-642-07718-0, S. 144, doi:10.1007/978-3-662-04775-0.
  4. Ralf Greve: Kontinuumsmechanik. Ein Grundkurs für Ingenieure und Physiker. Springer, Berlin u. a. 2003, ISBN 3-540-00760-1, S. 74.
  5. W. H. Müller: Streifzüge durch die Kontinuumstheorie. Springer, 2011, ISBN 978-3-642-19869-4, S. 72.
  6. Ralf Sube: Wörterbuch Physik Englisch: German-English. Routledge, London 2001, ISBN 978-0-415-17338-4 (eingeschränkte Vorschau in der Google-Buchsuche [abgerufen am 17. März 2017]).
  7. 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, hier S. 60, doi:10.1007/978-3-662-39776-3_1 (DOI der englischen Ausgabe).
  8. Martin H. Sadd: Elasticity – Theory, applications and numerics. Elsevier Butterworth-Heinemann, 2005, ISBN 0-12-605811-3, S. 110.