Deformationsgradient

Der Deformationsgradient (Formelzeichen: F {\displaystyle \mathbf {F} } ) ist in der Kontinuumsmechanik ein Mittel zur Beschreibung der lokalen Verformung an einem materiellen Punkt eines Körpers. Zur Veranschaulichung kann man sich einen Körper (in Abbildung 1, gelb) vorstellen auf den eine kurze Linie (weil nur lokale Änderungen beschrieben werden, im Bild fett rot) eingeritzt wird. Wird der Körper deformiert (rechts im Bild), wird die eingeritzte Linie nicht nur ihre Lage im Raum ändern, sondern auch gedehnt (oder gestaucht) und verdreht werden. Die Dehnung und Verdrehung beschreibt der Deformationsgradient und ist so ein Maß für die Deformation, daher der Name. Der Anhang Gradient verweist auf die Tatsache, dass lokale Änderungen beschrieben werden. Aus dem Deformationsgradient lassen sich Maße für die lokale Streckung, Verzerrung, Flächen- und Volumenänderung ableiten. Im allgemeinen Fall ist der Deformationsgradient sowohl vom Ort als auch von der Zeit abhängig. Die zeitliche Änderung des Deformationsgradienten gibt Maße für die Änderungsraten der Streckung, Verdrehung, Verzerrung, Flächen- und Volumenänderung. Der Deformationsgradient ist einheitenfrei.

Bei den angesprochenen kurzen Linien handelt es sich mathematisch um Vektoren, die vom Deformationsgradient transformiert werden, wobei die Vektoren im Allgemeinen gedreht und gestreckt werden. Abbildungen von Vektoren leisten Tensoren, siehe Abbildung 2, weswegen der Deformationsgradient ein Tensor ist. Wenn es klar ist, auf welches Koordinatensystem sich der Deformationsgradient bezieht, berechnet er sich wie eine Jacobimatrix und kann dann auch als Matrix notiert werden. Oft bildet der Deformationsgradient die (infinitesimal kleinen) materiellen Linienelemente in der Ausgangs- oder Referenzkonfiguration in die aktuelle oder Momentankonfiguration ab. Ganz allgemein kann eine solche Abbildung aber auch zwischen beliebig anderen zu definierenden Konfigurationen stattfinden.

Die Bewegung eines materiellen Punktes wird mit der Bewegungsfunktion

x → = χ → ( X → , t ) {\displaystyle {\vec {x}}={\vec {\chi }}({\vec {X}},t)}

beschrieben. Der Vektor x → {\displaystyle {\vec {x}}} ist die aktuelle Position des materiellen Punktes X → {\displaystyle {\vec {X}}} zur Zeit t {\displaystyle t} in der Momentankonfiguration. Genauer ist X → {\displaystyle {\vec {X}}} die Position des betrachteten materiellen Punktes in der undeformierten Ausgangs- oder Referenzkonfiguration des Körpers zu einer Zeit t 0 ≤ t {\displaystyle t_{0}\leq t} . Bei festgehaltenem materiellen Punkt X → {\displaystyle {\vec {X}}} beschreibt die Bewegungsfunktion dessen Bahnlinie durch den Raum. Im kartesischen Koordinatensystem mit der Standardbasis { e → 1 , e → 2 , e → 3 {\displaystyle {\vec {e}}_{1},{\vec {e}}_{2},{\vec {e}}_{3}} } hat der Punkt x → {\displaystyle {\vec {x}}} die komponentenweise Darstellung

x → = ∑ i = 1 3 x i e → i = ∑ i = 1 3 χ i ( X → , t ) e → i {\displaystyle {\vec {x}}=\sum _{i=1}^{3}{x}_{i}{\vec {e}}_{i}=\sum _{i=1}^{3}{\chi }_{i}({\vec {X}},t){\vec {e}}_{i}}

und entsprechend gilt X → = ∑ i = 1 3 X i e → i {\displaystyle {\vec {X}}=\sum _{i=1}^{3}{X}_{i}{\vec {e}}_{i}} . Um zu untersuchen, wie sich die aktuelle Position ändert, wenn die Position in der undeformierten Ausgangslage variiert, wird die Ableitung gebildet:

F i j = ∂ χ i ( X → , t ) ∂ X j {\displaystyle F_{ij}={\frac {\partial \chi _{i}({\vec {X}},t)}{\partial {X}_{j}}}}    für    i , j = 1 , 2 , 3 {\displaystyle i,j=1,2,3} .

Die Funktionen F i j {\displaystyle F_{ij}} sind die Komponenten des Deformationsgradienten bezüglich des Basissystems e → i {\displaystyle {\vec {e}}_{i}} .

Um zu einer koordinatenfreien Darstellung zu gelangen, wird das dyadische Produkt ⊗ {\displaystyle {\otimes }} benutzt:

F := G R A D χ → ( X → , t ) = ∂ χ i ∂ X j e → i ⊗ e → j = ( F 11 F 12 F 13 F 21 F 22 F 23 F 31 F 32 F 33 ) e → i ⊗ e → j {\displaystyle \mathbf {F} :=\mathrm {GRAD} \,{\vec {\chi }}({\vec {X}},t)={\frac {\partial \chi _{i}}{\partial X_{j}}}{\vec {e}}_{i}\otimes {\vec {e}}_{j}={\begin{pmatrix}F_{11}&F_{12}&F_{13}\\F_{21}&F_{22}&F_{23}\\F_{31}&F_{32}&F_{33}\end{pmatrix}}_{{\vec {e}}_{i}\otimes {\vec {e}}_{j}}} .

Darin ist F {\displaystyle \mathbf {F} } der Deformationsgradient und G R A D {\displaystyle \mathrm {GRAD} } ist der Operator für den materiellen Gradienten, denn es wird nach den materiellen Koordinaten X j {\displaystyle {X}_{j}} differenziert.

Der Deformationsgradient kann auch mit der Richtungsableitung

F ( X → , t ) ⋅ d X → = d d s χ → ( X → + s d X → , t ) | s = 0 = d x → {\displaystyle \mathbf {F} ({\vec {X}},t)\cdot \mathrm {d} {\vec {X}}=\left.{\frac {\mathrm {d} }{\mathrm {d} s}}{\vec {\chi }}{({\vec {X}}+s\;\mathrm {d} {\vec {X}},t)}\right|_{s=0}=\mathrm {d} {\vec {x}}}

dargestellt werden, was seine Transformationseigenschaften der Linienelemente d X → {\displaystyle \mathrm {d} {\vec {X}}} verdeutlicht, siehe Abbildung 3.

Mathematisch ist das Differential d X → {\displaystyle \mathrm {d} {\vec {X}}} in

d x → = F ⋅ d X → {\displaystyle \mathrm {d} {\vec {x}}=\mathbf {F} \cdot \mathrm {d} {\vec {X}}}

Element des Tangentialraumes T X → V {\displaystyle \mathrm {T} _{\vec {X}}V} im Punkt X → {\displaystyle {\vec {X}}} des Raumes V {\displaystyle V} , den der undeformierte Körper in der Ausgangskonfiguration einnimmt (in Abbildung 4 oben). Das Differential d x → = d χ → {\displaystyle \mathrm {d} {\vec {x}}=\mathrm {d} {\vec {\chi }}} ist entsprechend ein Element des Tangentialraumes T x → v {\displaystyle {\mathrm {T} }_{\vec {x}}v} im Punkt x → {\displaystyle {\vec {x}}} des Raumes v {\displaystyle v} , den der deformierte Körper in der Momentankonfiguration einnimmt (im Bild unten). Damit wird der Deformationsgradient zur Abbildung

F : T X → V ↦ T x → v d X → → d x → = F ( X → , t ) ⋅ d X → {\displaystyle {\begin{array}{llll}\mathbf {F} {:}&\mathrm {T} _{\vec {X}}V&\mapsto &\mathrm {T} _{\vec {x}}v\\&\mathrm {d} {\vec {X}}&\rightarrow &\mathrm {d} {\vec {x}}=\mathbf {F} ({\vec {X}},t)\cdot \mathrm {d} {\vec {X}}\end{array}}} .

Werden jedem materiellen Punkt X → {\displaystyle {\vec {X}}} über eine Referenzkonfiguration konvektive Koordinaten Θ → = ( Θ 1 , Θ 2 , Θ 3 ) {\displaystyle {\vec {\Theta }}=(\Theta _{1},\Theta _{2},\Theta _{3})} zugeordnet, bilden die Tangentenvektoren

G → i := ∂ X → ( Θ → ) ∂ Θ i {\displaystyle {\vec {G}}_{i}:={\dfrac {\partial {\vec {X}}({\vec {\Theta }})}{\partial \Theta _{i}}}}    bzw.    g → i := ∂ χ → ( X → ( Θ → ) ) ∂ Θ i {\displaystyle {\vec {g}}_{i}:={\dfrac {\partial {\vec {\chi }}\left({\vec {X}}({\vec {\Theta }})\right)}{\partial \Theta _{i}}}}

kovariante Basen der Tangentialräume T X → V {\displaystyle \mathrm {T} _{\vec {X}}V} im Punkt X → {\displaystyle {\vec {X}}} bzw. T x → v {\displaystyle \mathrm {T} _{\vec {x}}v} im Punkt x → {\displaystyle {\vec {x}}} (in Abbildung 4 schwarz dargestellt). Die Gradienten der konvektiven Koordinaten

G → i := ∑ j = 1 3 ∂ Θ i ∂ X j e → j {\displaystyle {\vec {G}}^{i}:=\displaystyle \sum _{j=1}^{3}{\dfrac {\partial \Theta _{i}}{\partial X_{j}}}{\vec {e}}_{j}}    bzw.    g → i := ∑ j = 1 3 ∂ Θ i ∂ x j e → j {\displaystyle {\vec {g}}^{i}:=\displaystyle \sum _{j=1}^{3}{\frac {\partial \Theta _{i}}{\partial x_{j}}}{\vec {e}}_{j}}

bilden kontravariante Basen, die zu den kovarianten dual sind (in Abbildung 4 blau dargestellt).

In diesen Basissystemen ausgedrückt bekommt der Deformationsgradient die besonders einfache Form

F = ∑ i = 1 3 g → i ⊗ G → i {\displaystyle \mathbf {F} =\sum _{i=1}^{3}{\vec {g}}_{i}\otimes {\vec {G}}^{i}} .

In dieser Darstellung lässt sich auch sofort mit

F − 1 = ∑ i = 1 3 G → i ⊗ g → i ⇔ g → i = F ⋅ G → i {\displaystyle \mathbf {F} ^{-1}=\sum _{i=1}^{3}{\vec {G}}_{i}\otimes {\vec {g}}^{i}\quad \Leftrightarrow \quad {\vec {g}}_{i}=\mathbf {F} \cdot {\vec {G}}_{i}} .

die Inverse des Deformationsgradienten angeben. Der transponiert inverse Deformationensgradient bildet die kontravarianten Basisvektoren aufeinander ab:

( F T ) − 1 = ∑ i = 1 3 g → i ⊗ G → i ⇔ g → i = ( F T ) − 1 ⋅ G → i {\displaystyle (\mathbf {F} ^{\mathrm {T} })^{-1}=\sum _{i=1}^{3}{\vec {g}}^{i}\otimes {\vec {G}}_{i}\quad \Leftrightarrow \quad {\vec {g}}^{i}=(\mathbf {F} ^{\mathrm {T} })^{-1}\cdot {\vec {G}}^{i}} .

Zumeist wird der Deformationsgradient wie oben in seiner materiellen Darstellung formuliert. Gelegentlich wird aber auch der räumliche Deformationsgradient benutzt. Wegen d e t ( F ) > 0 {\displaystyle \mathrm {det} (\mathbf {F} )>0} kann der Deformationsgradient invertiert und das Ergebnis über X → = χ → − 1 ( x → , t ) {\displaystyle {\vec {X}}={\vec {\chi }}^{-1}({\vec {x}},t)} als Funktion der räumlichen Koordinaten x → {\displaystyle {\vec {x}}} ausgedrückt werden:

f ( x → , t ) = F − 1 ( χ → − 1 ( x → , t ) , t ) = grad ⁡ ( X → , t ) = d ( χ → − 1 ) i d x j e → i ⊗ e → j {\displaystyle \mathbf {f} ({\vec {x}},t)=\mathbf {F} ^{-1}({\vec {\chi }}^{-1}({\vec {x}},t),t)=\operatorname {grad} ({\vec {X}},t)={\frac {\mathrm {d} ({\vec {\chi }}^{-1})_{i}}{\mathrm {d} x_{j}}}{\vec {e}}_{i}\otimes {\vec {e}}_{j}} .

Der räumliche Deformationsgradient f ( x → , t ) {\displaystyle \mathbf {f} ({\vec {x}},t)} bildet dann das Linienelement d x → {\displaystyle \mathrm {d} {\vec {x}}} auf das Linienelement d X → {\displaystyle \mathrm {d} {\vec {X}}} ab:

f : T x → v ↦ T X → V d x → → d X → = f ( x → , t ) ⋅ d x → {\displaystyle {\begin{array}{llll}\mathbf {f} :&\mathrm {T} _{\vec {x}}v&\mapsto &\mathrm {T} _{\vec {X}}V\\&\mathrm {d} {\vec {x}}&\rightarrow &\mathrm {d} {\vec {X}}=\mathbf {f} ({\vec {x}},t)\cdot \mathrm {d} {\vec {x}}\end{array}}} .

Entsprechend hat der räumliche Deformationsgradient in konvektiven Koordinaten die Form

f = ∑ i = 1 3 G → i ⊗ g → i {\displaystyle \mathbf {f} =\sum _{i=1}^{3}{\vec {G}}_{i}\otimes {\vec {g}}^{i}} .

In der Festkörpermechanik treten in vielen Anwendungsbereichen nur kleine Deformationen auf. In diesem Fall erfahren die Gleichungen der Kontinuumsmechanik eine erhebliche Vereinfachung durch geometrische Linearisierung. Dazu werden die Verschiebungen u → ( X → , t ) {\displaystyle {\vec {u}}({\vec {X}},t)} betrachtet, die ein materieller Punkt X → {\displaystyle {\vec {X}}} im Laufe seiner Bewegung erfährt. Weil x → = χ → ( X → , t ) {\displaystyle {\vec {x}}={\vec {\chi }}({\vec {X}},t)} die aktuelle Position des Punktes ist, der in der Ausgangskonfiguration die Position X → {\displaystyle {\vec {X}}} hatte, ist die Verschiebung die Differenz

u → = χ → ( X → , t ) − X → {\displaystyle {\vec {u}}={\vec {\chi }}({\vec {X}},t)-{\vec {X}}} .

Der materielle Gradient der Verschiebungen ist der Tensor

H = G R A D u → = G R A D χ → − G R A D X → = F − I {\displaystyle \mathbf {H} =\mathrm {GRAD} \,{\vec {u}}=\mathrm {GRAD} \,{\vec {\chi }}-\mathrm {GRAD} \,{\vec {X}}=\mathbf {F} -\mathbf {I} } .

und wird Verschiebungsgradient genannt. Wenn L 0 {\displaystyle L_{0}} eine charakteristische Abmessung des Körpers ist, dann wird bei kleinen Verschiebungen sowohl | u → | ≪ L 0 {\displaystyle |{\vec {u}}|\ll L_{0}} als auch ∥ H ∥≪ 1 {\displaystyle \parallel \mathbf {H} \parallel \ll 1} gefordert, so dass alle Terme, die höhere Potenzen von u → {\displaystyle {\vec {u}}} oder H {\displaystyle \mathbf {H} } beinhalten vernachlässigt werden können. In diesem Fall ergeben sich die folgenden Zusammenhänge:

F = I + H F − 1 ≈ I − H d e t ( F ) ≈ 1 + S p ( H ) U ≈ I + 1 2 ( H + H T ) v ≈ I + 1 2 ( H + H T ) R ≈ I + 1 2 ( H − H T ) {\displaystyle {\begin{array}{lcl}\mathbf {F} &=&\mathbf {I} +\mathbf {H} \\\mathbf {F} ^{-1}&\approx &\mathbf {I} -\mathbf {H} \\\mathrm {det} (\mathbf {F} )&\approx &1+\mathrm {Sp} (\mathbf {H} )\\\mathbf {U} &\approx &\mathbf {I} +{\frac {1}{2}}(\mathbf {H} +\mathbf {H} ^{\mathrm {T} })\\\mathbf {v} &\approx &\mathbf {I} +{\frac {1}{2}}(\mathbf {H} +\mathbf {H} ^{\mathrm {T} })\\\mathbf {R} &\approx &\mathbf {I} +{\frac {1}{2}}(\mathbf {H} -\mathbf {H} ^{\mathrm {T} })\end{array}}} .

Die Tensoren U , v {\displaystyle \mathbf {U} ,\mathbf {v} } und R {\displaystyle \mathbf {R} } kommen in der polaren Zerlegung vor, siehe unten.

Der Deformationsgradient F {\displaystyle \mathbf {F} } lässt sich eindeutig "polar" in eine Rotation und eine reine Streckung zerlegen. Durch Anwendung der Polarzerlegung resultiert die Darstellung

F = R ⋅ U = v ⋅ R {\displaystyle \mathbf {F} =\mathbf {R\cdot U} =\mathbf {v\cdot R} } .

Dabei ist R {\displaystyle \mathbf {R} } ein "eigentlich orthogonaler Tensor". Der materielle rechte Strecktensor U {\displaystyle \mathbf {U} } und der räumliche linke Strecktensor v {\displaystyle \mathbf {v} } sind symmetrisch und positiv definit. (Eselsbrücke: U {\displaystyle \mathbf {U} } steht rechts von R {\displaystyle \mathbf {R} } und v {\displaystyle \mathbf {v} } links davon in der polaren Darstellung.)

Anschaulich bedeutet die polare Zerlegung eine Hintereinanderschaltung zweier Transformationen: Im einen Fall eine rotationsfreie Streckung U {\displaystyle \mathbf {U} } mit anschließender Drehung R {\displaystyle \mathbf {R} } und im anderen Fall eine Drehung R {\displaystyle \mathbf {R} } mit anschließender rotationsfreier Streckung v {\displaystyle \mathbf {v} } so wie sie in Abbildung 5 dargestellt sind.

Der rechte Strecktensor berechnet sich gemäß

U = + F T ⋅ F {\displaystyle \mathbf {U} =+{\sqrt {\mathbf {F} ^{\mathrm {T} }\cdot \mathbf {F} }}}

aus der Hauptachsentransformation von F T ⋅ F {\displaystyle \mathbf {F} ^{\mathrm {T} }\cdot \mathbf {F} } , ziehen der Wurzel der Diagonalglieder und Rücktransformation, siehe auch das Beispiel unten. Entsprechend gilt für den linken Strecktensor

v = + F ⋅ F T {\displaystyle \mathbf {v} =+{\sqrt {\mathbf {F\cdot F} ^{\mathrm {T} }}}} .

Der Rotationstensor R {\displaystyle \mathbf {R} } ergibt sich dann aus

R = F ⋅ U − 1 = v − 1 ⋅ F {\displaystyle \mathbf {R} =\mathbf {F\cdot U} ^{-1}=\mathbf {v} ^{-1}\cdot \mathbf {F} } .

Mit Hilfe des Deformationsgradienten können Integrale in der materiellen Darstellung in die räumliche umgerechnet werden. Die zu integrierende Größe sei ein Feld q {\displaystyle q} das skalar-, vektor- oder tensorwertig sein kann und in der materiellen Darstellung q 0 ( X → , t ) {\displaystyle q_{0}({\vec {X}},t)} und der räumlichen

q ( x → , t ) = q 0 ( χ → − 1 ( x → , t ) , t ) {\displaystyle q({\vec {x}},t)=q_{0}\left({\vec {\chi }}^{-1}({\vec {x}},t),t\right)}

vorliege. Dann gelten die Identitäten:

∫ l q ( x → , t ) d x → = ∫ L q 0 ( X → , t ) F ⋅ d X → ∫ a q ( x → , t ) n → d a = ∫ A q 0 ( X → , t ) d e t ( F ) ( F T ) − 1 ⋅ N → d A ∫ v q ( x → , t ) d v = ∫ V q 0 ( X → , t ) d e t ( F ) d V {\displaystyle {\begin{array}{ll}\int _{l}q({\vec {x}},t)\mathrm {d} {\vec {x}}&=\int _{L}q_{0}({\vec {X}},t)\;\mathbf {F} \cdot \mathrm {d} {\vec {X}}\\\int _{a}q({\vec {x}},t){\vec {n}}\;\mathrm {d} a&=\int _{A}q_{0}({\vec {X}},t)\;\mathrm {det} (\mathbf {F} )(\mathbf {F} ^{\mathrm {T} })^{-1}\cdot {\vec {N}}\;\mathrm {d} A\\\int _{v}q({\vec {x}},t)\mathrm {d} v&=\int _{V}q_{0}({\vec {X}},t)\;\mathrm {det} (\mathbf {F} )\;\mathrm {d} V\end{array}}} .

Der Operator d e t ( ⋅ ) {\displaystyle \mathrm {det} (\cdot )} bildet die Determinante und ( ( ⋅ ) T ) − 1 {\displaystyle ((\cdot )^{\mathrm {T} })^{-1}} die transponiert Inverse. L {\displaystyle L} ist eine materielle Linie in der Ausgangskonfiguration und l {\displaystyle l} die zugehörige räumliche in der Momentankonfiguration. Die Oberfläche A {\displaystyle A} des Körpers in der Ausgangskonfiguration hat das Oberflächenelement N → d A {\displaystyle {\vec {N}}\mathrm {d} A} , d. h. die mit dem Flächenstück d A {\displaystyle \mathrm {d} A} multiplizierte Normale N → {\displaystyle {\vec {N}}} des Flächenstücks. Gleiches gilt für das räumliche Flächenelement n → d a {\displaystyle {\vec {n}}\mathrm {d} a} auf der Oberfläche a {\displaystyle a} des Körpers in der Momentankonfiguration. Diese Transformationen sind bei der Zeitableitung der Integrale auf den linken Seiten der Gleichungen nützlich, weil die Gebiete auf den linken Seiten von der Zeit abhängen nicht so aber auf den rechten Seiten.

Die Determinante von F {\displaystyle \mathbf {F} } gibt das lokale Volumenverhältnis J = d e t ( F ) = d v d V > 0 {\displaystyle J=\mathrm {det} (\mathbf {F} )={\frac {\mathrm {d} v}{\mathrm {d} V}}>0} im betrachteten materiellen Punkt bei der Deformation an.

Damit ergibt sich u. a., dass J = d e t ( F ) {\displaystyle J=\mathrm {det} (\mathbf {F} )} positiv sein muss, sonst wäre die Deformation physikalisch nicht möglich (Inversion des materiellen Punktes).

Bleibt bei einer Deformation das Volumen erhalten, also J = 1 {\displaystyle J=1} , liegt Inkompressibilität vor. Bei Gummi- oder Elastomer-Werkstoffen ist dies eine übliche Annahme in der kontinuumsmechanischen Beschreibung und durch das Verhalten dieser Werkstoffklasse annähernd der Fall. Gleiches gilt für die inkompressiblen Flüssigkeiten.

Der Deformationsgradient transformiert neben den Linien-, Flächen- und Volumenelementen auch Tensoren von der Ausgangskonfiguration in die Momentankonfiguration. Diese Transformationen sind für kovariante Tensoren (oftmals Verzerrungstensoren) und kontravariante Tensoren (oftmals Spannungstensoren) unterschiedlich, z. B.:

e = ( F T ) − 1 ⋅ E ⋅ F − 1 S = F ⋅ T ~ ⋅ F T {\displaystyle {\begin{array}{lcl}\mathbf {e} &=&(\mathbf {F} ^{\mathrm {T} })^{-1}\cdot \mathbf {E} \cdot \mathbf {F} ^{-1}\\\mathbf {S} &=&\mathbf {F} \cdot {\tilde {\mathbf {T} }}\cdot \mathbf {F} ^{\mathrm {T} }\end{array}}} .

Der Tensor

e := 1 2 ( I − ( F T ) − 1 ⋅ F − 1 ) {\displaystyle \mathbf {e} :={\tfrac {1}{2}}(\mathbf {I} -(\mathbf {F} ^{\mathrm {T} })^{-1}\cdot \mathbf {F} ^{-1})}

ist der Euler-Almansi-Verzerrungstensor in der Momentankonfiguration,

E := 1 2 ( F T ⋅ F − I ) {\displaystyle \mathbf {E} :={\tfrac {1}{2}}({\mathbf {F} }^{\mathrm {T} }\cdot \mathbf {F} -\mathbf {I} )}

der Green-Lagrange’sche Verzerrungstensor in der Ausgangskonfiguration,

S := d e t ( F ) T {\displaystyle \mathbf {S} :=\mathrm {det} (\mathbf {F} )\mathbf {T} }

der gewichtete Cauchy’sche Spannungstensor, T {\displaystyle \mathbf {T} } der Cauchy’sche Spannungstensor (beide in der Momentankonfiguration) und

T ~ := det ( F ) F − 1 ⋅ T ⋅ ( F T ) − 1 {\displaystyle {\tilde {\mathbf {T} }}:=\det(\mathbf {F} )\mathbf {F} ^{-1}\cdot \mathbf {T} \cdot (\mathbf {F} ^{\mathrm {T} })^{-1}}

der zweite Piola-Kirchhoff’sche Spannungstensor in der Ausgangskonfiguration. Das Skalarprodukt ":" der so einander zugeordneten Tensoren wird von der Transformation nicht verändert, z. B.:

S : e = F ⋅ T ~ ⋅ F T : ( F T ) − 1 ⋅ E ⋅ F − 1 = T ~ : E {\displaystyle \mathbf {S} :\mathbf {e} =\mathbf {F} \cdot {\tilde {\mathbf {T} }}\cdot \mathbf {F} ^{\mathrm {T} }:(\mathbf {F} ^{\mathrm {T} })^{-1}\cdot \mathbf {E} \cdot \mathbf {F} ^{-1}={\tilde {\mathbf {T} }}:\mathbf {E} } .

Die Multiplikative Zerlegung des Deformationsgradienten wird in Materialtheorie angewendet, um die Deformation eines Körpers auf Grund verschiedener Einflüsse zu modellieren. So kann sich ein Körper deformieren, weil er erwärmt wird oder einer äußeren Kraft ausgesetzt wird. Die Deformation kann zusätzlich davon abhängen, wie schnell die Temperatur oder Kraft aufgebracht wird. Die Reaktionen des Materials lassen sich einfacher modellieren, wenn die Phänomene voneinander getrennt betrachtet werden. So kann ein Modell den Einfluss der Temperatur nachbilden und ein anderes Modell die isotherme Verformung durch Kräfte. Die Deformationen aufgrund des einen oder anderen Phänomens können dann anschließend wieder zusammengeführt werden. In der Materialtheorie hat es sich durchgesetzt bei kleinen Verformungen eine additive Zerlegung der Dehnungen und bei großen Verformungen eine multiplikative Zerlegung des Deformationsgradienten zu benutzen.

Seien also a und b zwei Verformungsanteile eines Materials. Für die Modellbildung wird der Deformationsgradient multiplikativ zerlegt:

F = F a ⋅ F b {\displaystyle \mathbf {F} =\mathbf {F} _{a}\cdot \mathbf {F} _{b}} .

Ein Modell beschreibt dann die Entwicklung des Anteils F a {\displaystyle \mathbf {F} _{a}} auf Grund des Einflusses a und ein anderes Modell die Entwicklung des Anteils F b {\displaystyle \mathbf {F} _{b}} auf Grund des Einflusses b. Die Deformationsgradienten F a {\displaystyle \mathbf {F} _{a}} und F b {\displaystyle \mathbf {F} _{b}} erfüllen im Allgemeinen nicht die Kompatibilitätsbedingungen, weswegen es im Allgemeinen kein Bewegungsfeld gibt, aus dem die beiden Anteile per Gradientenbildung abgeleitet werden können.

Weil der Deformationsgradient, wie im vorigen Abschnitt erläutert, Tensoren von einer Konfiguration in die andere transformiert, entspricht die multiplikative Zerlegung der Einführung einer Zwischenkonfiguration. Die Tensoren der Referenzkonfiguration werden mit dem Anteil F b {\displaystyle \mathbf {F} _{b}} und die der Momentankonfiguration mit F a {\displaystyle \mathbf {F} _{a}} in die Zwischenkonfiguration transformiert. Transformation des Green-Lagrange’schen Verzerrungstensors von der Referenzkonfiguration in die Zwischenkonfiguration liefert:

Γ := ( F b T ) − 1 ⋅ E ⋅ F b − 1 = ( F b T ) − 1 ⋅ 1 2 ( ( F a ⋅ F b ) T ⋅ ( F a ⋅ F b ) − I ) ⋅ F b − 1 =: Γ a + Γ b Γ a := 1 2 ( F a T ⋅ F a − I ) Γ b := 1 2 ( I − ( F b T ) − 1 ⋅ F b − 1 ) . {\displaystyle {\begin{array}{lcl}{\boldsymbol {\Gamma }}&:=&(\mathbf {F} _{b}^{\mathrm {T} })^{-1}\cdot \mathbf {E} \cdot \mathbf {F} _{b}^{-1}=(\mathbf {F} _{b}^{\mathrm {T} })^{-1}\cdot {\frac {1}{2}}\left((\mathbf {F} _{a}\cdot \mathbf {F} _{b})^{\mathrm {T} }\cdot (\mathbf {F} _{a}\cdot \mathbf {F} _{b})-\mathbf {I} \right)\cdot \mathbf {F} _{b}^{-1}=:{\boldsymbol {\Gamma }}_{a}+{\boldsymbol {\Gamma }}_{b}\\{\boldsymbol {\Gamma }}_{a}&:=&{\frac {1}{2}}\left(\mathbf {F} _{a}^{\mathrm {T} }\cdot \mathbf {F} _{a}-\mathbf {I} \right)\\{\boldsymbol {\Gamma }}_{b}&:=&{\frac {1}{2}}\left(\mathbf {I} -(\mathbf {F} _{b}^{\mathrm {T} })^{-1}\cdot \mathbf {F} _{b}^{-1}\right).\end{array}}}

Der in die Zwischenkonfiguration transformierte Green-Lagrange’sche Verzerrungstensor zerfällt also in zwei Anteile:

• Ein Anteil Γ a {\displaystyle {\boldsymbol {\Gamma }}_{a}} ist vom Green-Lagrange-Typ und wird mit F a {\displaystyle \mathbf {F} _{a}} gebildet.
• Der andere Anteil ist vom Euler-Almansi-Typ und wird mit F b {\displaystyle \mathbf {F} _{b}} gebildet.

Gleiches gilt, wenn der Euler-Almansi Tensor e {\displaystyle \mathbf {e} } mit F a {\displaystyle \mathbf {F} _{a}} auf die Zwischenkonfiguration transformiert wird:

F a T ⋅ e ⋅ F a = F a T ⋅ 1 2 ( I − ( F a ⋅ ( F b ) T ) − 1 ⋅ ( F a ⋅ F b ) − 1 ) ⋅ F a = Γ a + Γ b {\displaystyle \mathbf {F} _{a}^{\mathrm {T} }\cdot \mathbf {e} \cdot \mathbf {F} _{a}=\mathbf {F} _{a}^{\mathrm {T} }\cdot {\frac {1}{2}}\left(\mathbf {I} -(\mathbf {F} _{a}\cdot (\mathbf {F} _{b})^{\mathrm {T} })^{-1}\cdot (\mathbf {F} _{a}\cdot \mathbf {F} _{b})^{-1}\right)\cdot \mathbf {F} _{a}={\boldsymbol {\Gamma }}_{a}+{\boldsymbol {\Gamma }}_{b}}

In der Zwischenkonfiguration können nun die beiden Phänomene mit den Verzerrungstensoren Γ a {\displaystyle {\boldsymbol {\Gamma }}_{a}} und Γ b {\displaystyle {\boldsymbol {\Gamma }}_{b}} getrennt modelliert werden. Der sich im Modell in der Zwischenkonfiguration ergebende Spannungstensor wird anschließend mit F a {\displaystyle \mathbf {F} _{a}} in die Momentankonfiguration oder mit F b {\displaystyle \mathbf {F} _{b}} in die Referenzkonfiguration transformiert.

Die materielle Zeitableitung des Deformationsgradienten

F ˙ := G R A D χ → ˙ ( X → , t ) = d χ ˙ i d X j e → i ⊗ e → j {\displaystyle {\dot {\mathbf {F} }}:=\mathrm {GRAD} \,{\dot {\vec {\chi }}}({\vec {X}},t)={\frac {\mathrm {d} {\dot {\chi }}_{i}}{\mathrm {d} X_{j}}}{\vec {e}}_{i}\otimes {\vec {e}}_{j}}

ist ein Maß für die Deformationsgeschwindigkeit. Sie hängt über

l := g r a d v → ( x → , t ) = d v i d x j e → i ⊗ e → j = ∑ i , j , k = 1 3 d v i d X k e → i ⊗ e → k ⋅ d X l d x j e → l ⊗ e → j = F ˙ ⋅ F − 1 {\displaystyle \mathbf {l} :=\mathrm {grad} \,{\vec {v}}({\vec {x}},t)={\frac {\mathrm {d} v_{i}}{\mathrm {d} x_{j}}}{\vec {e}}_{i}\otimes {\vec {e}}_{j}=\sum _{i,j,k=1}^{3}{\frac {\mathrm {d} v_{i}}{\mathrm {d} X_{k}}}{\vec {e}}_{i}\otimes {\vec {e}}_{k}\cdot {\frac {\mathrm {d} X_{l}}{\mathrm {d} x_{j}}}{\vec {e}}_{l}\otimes {\vec {e}}_{j}={\dot {\mathbf {F} }}\cdot \mathbf {F} ^{-1}}

mit dem räumlichen Geschwindigkeitsgradient l {\displaystyle \mathbf {l} } des räumlichen Geschwindigkeitsfeldes v → ( x → , t ) {\displaystyle {\vec {v}}({\vec {x}},t)} zusammen. In konvektiven Koordinaten lautet das

F ˙ = ∑ i = 1 3 g → ˙ i ⊗ G → i l = ∑ i = 1 3 g → ˙ i ⊗ g → i = − ∑ i = 1 3 g → i ⊗ g → ˙ i . {\displaystyle {\begin{array}{lcl}{\dot {\mathbf {F} }}&=&\displaystyle \sum _{i=1}^{3}{\dot {\vec {g}}}_{i}\otimes {\vec {G}}^{i}\\\mathbf {l} &=&\displaystyle \sum _{i=1}^{3}{\dot {\vec {g}}}_{i}\otimes {\vec {g}}^{i}=-\sum _{i=1}^{3}{\vec {g}}_{i}\otimes {\dot {\vec {g}}}^{i}.\end{array}}}

Die Berechnung des Deformationsgradienten und seiner polare Zerlegung wird anhand der Scherung eines Quadrates vorgeführt.

Ein Quadrat der Kantenlänge eins wird zu einem flächengleichen Parallelogramm mit Grundseite und Höhe eins verformt, siehe Abbildung 6. Die Punkte des Quadrates haben in der Ausgangskonfiguration die Koordinaten

( X Y ) ∈ [ 0 , 1 ] 2 ⊂ R 2 {\displaystyle {\begin{pmatrix}X\\Y\end{pmatrix}}\in [0,1]^{2}\subset \mathbb {R} ^{2}} .

Die Neigung des Parallelogramms sei

tan ⁡ ( α ) = 5 6 → α ≈ 0,694 7 r a d ≈ 39,805 5 ∘ {\displaystyle \tan(\alpha )={\frac {5}{6}}\quad \rightarrow \quad \alpha \approx 0{,}6947\,\mathrm {rad} \approx 39{,}8055^{\circ }} .

Dann sind die räumlichen Koordinaten der Punkte gegeben durch

( x y ) = χ → ( ( X Y ) , t ) = ( X + tan ⁡ ( α ) Y Y ) = ( X + 5 6 Y Y ) {\displaystyle {\begin{pmatrix}x\\y\end{pmatrix}}={\vec {\chi }}\left({\begin{pmatrix}X\\Y\end{pmatrix}},t\right)={\begin{pmatrix}X+\tan(\alpha )Y\\Y\end{pmatrix}}={\begin{pmatrix}X+{\frac {5}{6}}Y\\Y\end{pmatrix}}} .

Wie üblich wird X {\displaystyle X} mit X 1 {\displaystyle X_{1}} und Y {\displaystyle Y} mit X 2 {\displaystyle X_{2}} identifiziert. Dann bekommt man den Deformationsgradient durch Ableitung:

F = ∑ i , j = 1 2 ∂ x i ∂ X j e → i ⊗ e → j = ( 1 5 6 0 1 ) {\displaystyle \mathbf {F} =\sum _{i,j=1}^{2}{\frac {\partial x_{i}}{\partial X_{j}}}{\vec {e}}_{i}\otimes {\vec {e}}_{j}={\begin{pmatrix}1&{\frac {5}{6}}\\0&1\end{pmatrix}}} .

Die Richtungsableitung liefert über

F ⋅ ( d X d Y ) = ∂ ∂ s χ → ( ( X Y ) + s ( d X d Y ) , t ) | s = 0 = ∂ ∂ s ( ( X + s d X ) + 5 6 ( Y + s d Y ) ( Y + s d Y ) ) | s = 0 = ( d X + 5 6 d Y d Y ) = ( 1 5 6 0 1 ) ( d X d Y ) {\displaystyle {\begin{array}{lcl}\mathbf {F} \cdot {\begin{pmatrix}\mathrm {d} X\\\mathrm {d} Y\end{pmatrix}}&=&{\dfrac {\partial }{\partial s}}{\vec {\chi }}\left.\left({\begin{pmatrix}X\\Y\end{pmatrix}}+s{\begin{pmatrix}\mathrm {d} X\\\mathrm {d} Y\end{pmatrix}},t\right)\right|_{s=0}=\left.{\dfrac {\partial }{\partial s}}{\begin{pmatrix}(X+s\mathrm {d} X)+{\frac {5}{6}}(Y+s\mathrm {d} Y)\\(Y+s\mathrm {d} Y)\end{pmatrix}}\right|_{s=0}\\[3ex]&=&{\begin{pmatrix}\mathrm {d} X+{\frac {5}{6}}\mathrm {d} Y\\\mathrm {d} Y\end{pmatrix}}={\begin{pmatrix}1&{\frac {5}{6}}\\0&1\end{pmatrix}}{\begin{pmatrix}\mathrm {d} X\\\mathrm {d} Y\end{pmatrix}}\end{array}}}

dasselbe Ergebnis. Der Deformationsgradient ist hier vom Ort und der Zeit unabhängig und hat die Determinante eins, was den Erhalt des Flächeninhalts bestätigt. Der rechte Strecktensor U {\displaystyle \mathbf {U} } berechnet sich aus dem rechten Cauchy-Green-Tensor

C = F T ⋅ F = ( 1 0 5 6 1 ) ( 1 5 6 0 1 ) = ( 1 5 6 5 6 61 36 ) {\displaystyle \mathbf {C} =\mathbf {F} ^{\mathrm {T} }\cdot \mathbf {F} ={\begin{pmatrix}1&0\\{\frac {5}{6}}&1\end{pmatrix}}{\begin{pmatrix}1&{\frac {5}{6}}\\0&1\end{pmatrix}}={\begin{pmatrix}1&{\frac {5}{6}}\\{\frac {5}{6}}&{\frac {61}{36}}\end{pmatrix}}}

über Hauptachsentransformation, ziehen der Wurzel der Diagonalglieder und Rücktransformation. Für die Hauptachsentransformation braucht man die Eigenwerte λ 1 , 2 {\displaystyle \lambda _{1,2}} und -vektoren v → 1 , 2 {\displaystyle {\vec {v}}_{1,2}} von C {\displaystyle \mathbf {C} } . Man findet:

λ 1 = 9 4 , v → 1 = 1 13 ( 2 3 ) λ 2 = 4 9 , v → 2 = 1 13 ( 3 − 2 ) {\displaystyle {\begin{array}{lclclcl}\lambda _{1}&=&{\dfrac {9}{4}}&,\quad &{\vec {v}}_{1}&=&{\dfrac {1}{\sqrt {13}}}{\begin{pmatrix}2\\3\end{pmatrix}}\\[3ex]\lambda _{2}&=&{\dfrac {4}{9}}&,\quad &{\vec {v}}_{2}&=&{\dfrac {1}{\sqrt {13}}}{\begin{pmatrix}3\\-2\end{pmatrix}}\end{array}}} .

Mit diesen Eigenwerten und -vektoren erhält man die Hauptachsentransformation

C = ∑ i = 1 2 λ i v → i ⊗ v → i = ( ∑ j = 1 2 v → j ⊗ e → j ) ⋅ ( ∑ i = 1 2 λ i e → i ⊗ e → i ) ⋅ ( ∑ k = 1 2 e → k ⊗ v → k ) = 1 13 ( 2 3 3 − 2 ) ( 9 4 0 0 4 9 ) 1 13 ( 2 3 3 − 2 ) {\displaystyle {\begin{array}{lcl}\mathbf {C} &=&\displaystyle \sum _{i=1}^{2}\lambda _{i}{\vec {v}}_{i}\otimes {\vec {v}}_{i}=\left(\sum _{j=1}^{2}{\vec {v}}_{j}\otimes {\vec {e}}_{j}\right)\cdot \left(\sum _{i=1}^{2}\lambda _{i}{\vec {e}}_{i}\otimes {\vec {e}}_{i}\right)\cdot \left(\sum _{k=1}^{2}{\vec {e}}_{k}\otimes {\vec {v}}_{k}\right)\\[3ex]&=&{\dfrac {1}{\sqrt {13}}}{\begin{pmatrix}2&3\\3&-2\end{pmatrix}}{\begin{pmatrix}{\frac {9}{4}}&0\\0&{\frac {4}{9}}\end{pmatrix}}{\dfrac {1}{\sqrt {13}}}{\begin{pmatrix}2&3\\3&-2\end{pmatrix}}\end{array}}}

und damit den rechten Strecktensor

U = + C = 1 13 ( 2 3 3 − 2 ) ( + 9 4 0 0 + 4 9 ) ( 2 3 3 − 2 ) = 1 13 ( 12 5 5 97 6 ) {\displaystyle \mathbf {U} =+{\sqrt {\mathbf {C} }}={\frac {1}{13}}{\begin{pmatrix}2&3\\3&-2\end{pmatrix}}{\begin{pmatrix}+{\sqrt {\frac {9}{4}}}&0\\0&+{\sqrt {\frac {4}{9}}}\end{pmatrix}}{\begin{pmatrix}2&3\\3&-2\end{pmatrix}}={\frac {1}{13}}{\begin{pmatrix}12&5\\5&{\frac {97}{6}}\end{pmatrix}}} .

Mit seiner Inversen

U − 1 = 1 13 ( 97 6 − 5 − 5 12 ) {\displaystyle \mathbf {U} ^{-1}={\frac {1}{13}}{\begin{pmatrix}{\frac {97}{6}}&-5\\-5&12\end{pmatrix}}}

ergibt sich der Rotationstensor

R = F ⋅ U − 1 = 1 13 ( 12 5 − 5 12 ) = ( cos ⁡ ( β ) − sin ⁡ ( β ) sin ⁡ ( β ) cos ⁡ ( β ) ) {\displaystyle \mathbf {R} =\mathbf {F\cdot U} ^{-1}={\frac {1}{13}}{\begin{pmatrix}12&5\\-5&12\end{pmatrix}}={\begin{pmatrix}\cos(\beta )&-\sin(\beta )\\\sin(\beta )&\cos(\beta )\end{pmatrix}}} ,

siehe Drehmatrix. Der Rotationstensor dreht das im Bild königsblaue Parallelogramm oder hellblaue Quadrat um den Winkel

β = | arcsin ⁡ ( − 5 13 ) | ≈ 0,394 8 r a d ≈ 22,619 9 ∘ {\displaystyle \beta =\left|\arcsin \left(-{\frac {5}{13}}\right)\right|\approx 0{,}3948\,\mathrm {rad} \approx 22{,}6199^{\circ }} .

Den linken Strecktensor kann man nun einfacher aus

v = F ⋅ R T = ( 1 5 6 0 1 ) 1 13 ( 12 − 5 5 12 ) = 1 13 ( 97 6 5 5 12 ) {\displaystyle \mathbf {v} =\mathbf {F\cdot R} ^{\mathrm {T} }={\begin{pmatrix}1&{\frac {5}{6}}\\0&1\end{pmatrix}}{\frac {1}{13}}{\begin{pmatrix}12&-5\\5&12\end{pmatrix}}={\frac {1}{13}}{\begin{pmatrix}{\frac {97}{6}}&5\\5&12\end{pmatrix}}}

ermitteln.

Mathematik:

• Differentialgeometrie
• Jacobimatrix

Formelsammlungen:

• Formelsammlung Tensoralgebra
• Formelsammlung Tensoranalysis

• H. Altenbach: Kontinuumsmechanik. Springer 2012, ISBN 978-3-642-24118-5.
• P. Haupt: Continuum Mechanics and Theory of Materials. Springer Verlag, 2000, ISBN 3-540-66114-X.