Das Koordinatensystem mit den drei Orthotropieachsen Radial, Transversal, Longitudinal
Holz als typisches orthotropes Material im Ingenieurwesen
Dieses orthotrope Material ist aufgrund seiner inneren Struktur rotationssymmetrisch bezüglich einer Drehung um 180 Grad um eine Achse senkrecht zur Blattebene. Diese Achse, die rote und grüne Linie sind seine Orthotropieachsen. Denn es ist auch symmetrisch gegenüber Drehungen um 180 Grad um die rote und grüne Achse.
Die Orthotropie (von griechisch {{Modul:Vorlage:lang}} Modul:Multilingual:149: attempt to index field 'data' (a nil value) {{Modul:Vorlage:lang}} Modul:Multilingual:149: attempt to index field 'data' (a nil value) „korrekt, senkrecht, gerade stehend“ und {{Modul:Vorlage:lang}} Modul:Multilingual:149: attempt to index field 'data' (a nil value) {{Modul:Vorlage:lang}} Modul:Multilingual:149: attempt to index field 'data' (a nil value) „Weg, Art und Weise“) ist eine spezielle Art der Richtungsabhängigkeit eines Werkstoffs/Materials. Orthotrope Materialien wie im Bild haben die folgenden Eigenschaften:
- Das Kraft-Verformungs-Verhalten ändert sich nicht, wenn das Material um 180 Grad um die Orthotropieachsen gedreht wird.
- Im Bezugssystem parallel zu den Orthotropieachsen gibt es keine Kopplung zwischen Normaldehnungen und Schubverzerrungen.
Ein linear elastisches orthotropes Material besitzt maximal neun Materialparameter.
Ein Material ist isotrop, wenn es richtungsunabhängig dasselbe Kraft-Verformungs-Verhalten hat. Bei anisotropen Materialien ist das Kraft-Verformungs-Verhalten von der Belastungsrichtung abhängig. Die Orthotropie ist ein Spezialfall der Anisotropie und enthält ihrerseits die kubische Anisotropie, transversale Isotropie und Isotropie als Sonderfälle.
Viele Konstruktionswerkstoffe sind orthotrop, z. B. technisches Holz, Gewebe, viele Faser-Kunststoff-Verbunde und Walzbleche mit Textur. Kristalle des rhombischen Kristallsystems sind orthotrop[1]:390, Spezialfälle kommen im tetragonalen[1]:391, hexagonalen[1]:393 und kubischen Kristallsystem vor.
Symmetriegruppe
Die Richtungsabhängigkeit eines Materials zeichnet sich dadurch aus, dass das Kraft-Verformungs-Verhalten unabhängig (invariant) ist gegenüber nur bestimmten Drehungen des Materials: Bei der Orthotropie sind dies alle 180-Grad-Drehungen um die Orthotropieachsen. Diese Drehungen bilden zusammen mit der Punktspiegelung die Symmetriegruppe des orthotropen Materials.[1]:380
Die Invarianz gegenüber diesen Drehungen des Materials veranschaulichen zwei Experimente an einem Teilchen: Im ersten Experiment bringt man am Teilchen eine bestimmte Kraft auf und misst die resultierende Verformung. Im zweiten Experiment dreht man das Material zunächst nacheinander um beliebige Orthotropieachsen – um 180 Grad. Dann bringt man dieselbe Kraft auf wie im ersten Experiment und misst erneut die Verformung. Bei orthotropem Material wird man im zweiten Experiment dieselbe Verformung messen wie im ersten. Und zwar auch bei nicht-linear elastischem Materialverhalten.
Die Abhängigkeit von den Drehungen des Materials erkennt man, wenn man im zweiten Experiment um einen anderen Winkel als 180 Grad dreht. Wenn nicht der Spezialfall transversale Isotropie oder Isotropie vorliegt, wird man nun immer eine andere Verformung messen als im ersten Experiment.
Die angesprochenen Drehungen werden in der Kontinuumsmechanik durch orthogonale Tensoren Q repräsentiert. Eine Symmetriegruppe gR besteht aus denjenigen Transformationen, die die Formänderungsenergie e invariant lassen. Mathematisch wird das mit dem Verzerrungstensor E durch
- $ \mathbf {Q} \in g_{R}\quad \leftrightarrow \quad e(\mathbf {Q\cdot E\cdot Q} ^{\top })=e(\mathbf {E} ) $ für alle E
ausgedrückt.[1]:379 Darin bedeutet „·“ das Matrizenprodukt und das hochgestellte „⊤“ eine Transponierung. Mit Q gehört auch -Q zur Symmetriegruppe, was durch Hinzufügen des negativen Einheitstensors -1, der eine Punktspiegelung repräsentiert, zu gR berücksichtigt wird. Die Symmetriegruppe des orthotropen Materials ist[1]:382
- $ g_{R}=\left\{-\mathbf {1} ,\mathbf {Q} _{1}^{\pi },\mathbf {Q} _{2}^{\pi }\right\} $
Darin steht $ \mathbf {Q} _{i}^{\alpha } $ für den orthogonalen Tensor, der mit dem Winkel α in Radiant um die $ i $-te Orthotropieachse dreht. Die 180-Grad-Drehung um die 3-Achse ist in gR enthalten, denn
- $ \mathbf {Q} _{1}^{\pi }\cdot \mathbf {Q} _{2}^{\pi }\cdot \mathbf {Q} _{3}^{\pi }=\mathbf {1} =(-\mathbf {1} )^{2}\in g_{R} $
Invarianten
In der isotropen Hyperelastizität hängt die Formänderungsenergie von den Hauptinvarianten I1,2,3 des Verzerrungstensors E ab:
- e(E)=e(I1, I2, I3)
Die analoge Darstellung der Anisotropie erfordert, dass ein komplettes System von skalarwertigen Funktionen bekannt ist, die unter allen Transformationen in der Symmetriegruppe gR invariant sind.[1]:380 Bei der Orthotropie bleiben die folgenden Terme invariant:[1]:382
- E11, E22, E33, E232, E132, E122, E12E23E13.
Darin ist Eij := êi·E·êj für i,j=1,2,3 und ê1,2,3 sind die Einheitsvektoren in Richtung der paarweise orthogonalen Orthotropieachsen.
Orthotropie in der Linearen Elastizitätstheorie
Gegeben sind zwei Tensoren zweiter Stufe $ {\boldsymbol {\sigma }} $ und $ {\boldsymbol {\varepsilon }} $ mit 3×3-Koeffizienten $ \sigma _{ij} $ bzw. $ \varepsilon _{ij} $. Der allgemeinste lineare Zusammenhang, den es zwischen diesen Koeffizienten gibt, ist:
- $ f_{C}:\varepsilon _{kl}\rightarrow \sigma _{ij}=\sum _{k,l=1}^{3}C_{ijkl}\varepsilon _{kl} $.
Darin sind $ C_{ijkl} $ 81 Koeffizienten mit denen die neun Komponenten $ \varepsilon _{ij} $ auf neun Komponenten $ \sigma _{ij} $ abgebildet werden. In der linearen Elastizitätstheorie, in der $ {\boldsymbol {\sigma }} $ der symmetrische Spannungstensor und $ {\boldsymbol {\varepsilon }} $ der symmetrische Verzerrungstensor ist, reduziert sich die Anzahl der unabhängigen Tensor-Komponenten auf sechs, so dass nur 36 Koeffizienten $ C_{ijkl} $ unabhängig sind. Diesen Zusammenhang zwischen Spannungen und Verzerrungen kann man nun in Voigt’scher Notation auch als Matrizengleichung schreiben:
- $ \sigma _{i}^{\text{v}}=\sum _{j=1}^{3}C_{ij}^{\text{v}}\varepsilon _{j}^{\text{v}}\Leftrightarrow {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\\\end{bmatrix}}={\begin{bmatrix}C_{1111}&C_{1122}&C_{1133}&C_{1123}&C_{1113}&C_{1112}\\C_{2211}&C_{2222}&C_{2233}&C_{2223}&C_{2213}&C_{2212}\\C_{3311}&C_{3322}&C_{3333}&C_{3323}&C_{3313}&C_{3312}\\C_{2311}&C_{2322}&C_{2333}&C_{2323}&C_{2313}&C_{2312}\\C_{1311}&C_{1322}&C_{1333}&C_{1323}&C_{1313}&C_{1312}\\C_{1211}&C_{1222}&C_{1233}&C_{1223}&C_{1213}&C_{1212}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\\\end{bmatrix}} $.
Die Matrix $ C^{\text{v}} $ mit den 36 unabhängigen Komponenten $ C_{ijkl} $ repräsentiert den Elastizitätstensor des Materials. Im Fall der Hyperelastizität ist diese Matrix symmetrisch, so dass dann nur noch 21 Einträge unabhängig sind.
Materialparameter
Die Koeffizienten $ C_{ijkl} $ des Elastizitätstensors ergeben sich bei orthotroper linearer Elastizität aus nur neun Materialparametern, die in Versuchen an makroskopischen Proben ermittelt werden können:
Formelzeichen |
Bedeutung
|
$ E_{1},E_{2},E_{3} $ |
Elastizitätsmoduln in den Orthotropieachsen
|
$ G_{12},G_{13},G_{23} $ |
Schubmoduln in Ebenen senkrecht zu den Orthotropieachsen
|
$ \nu _{12},\nu _{13},\nu _{23} $ |
Querkontraktionszahlen bei Zug in Richtung einer Orthotropieachse
|
Die Dimension der Elastizitätsmoduln $ E_{1},E_{2},E_{3} $ und Schubmoduln $ G_{12},G_{23},G_{13} $ ist Kraft pro Fläche während die Querkontraktionszahlen $ \nu _{ij} $ dimensionslos sind.
Die Querkontraktionszahlen beschreiben, wie sich eine entlang einer Richtung – z. B. der 1-Richtung – gezogene Materialprobe quer dazu – z. B. in 2-Richtung – kontrahiert. Die entsprechende Querkontraktionszahl wäre dann $ \nu _{12} $. Die Normaldehnung in i-Richtung wird mit $ \varepsilon _{ii} $ bezeichnet. Dann ist für beliebige Werkstoffe die Querkontraktionszahl $ \nu _{ij} $ das negative Verhältnis der Normaldehnung in j-Richtung (Wirkung) zu derjenigen in i-Richtung bei Zug in i-Richtung (Ursache):
- $ \nu _{ij}=-{\frac {\varepsilon _{jj}}{\varepsilon _{ii}}} $
Wegen des Ursache-Wirkungs-Konzepts ist meistens $ \nu _{ij}\nu _{ji}\neq 1 $.
Elastizitätsgesetz für 3D
Ein Material ist linear elastisch orthotrop, wenn eine Orthonormalbasis existiert, so dass das Elastizitätsgesetz dargestellt in Bezug auf diese Basis folgende Form annimmt:
- $ {\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\end{bmatrix}}=\underbrace {\begin{bmatrix}{\frac {1}{E_{1}}}&-{\frac {\nu _{21}}{E_{2}}}&-{\frac {\nu _{31}}{E_{3}}}&&&\\-{\frac {\nu _{12}}{E_{1}}}&{\frac {1}{E_{2}}}&-{\frac {\nu _{32}}{E_{3}}}&&&\\-{\frac {\nu _{13}}{E_{1}}}&-{\frac {\nu _{23}}{E_{2}}}&{\frac {1}{E_{3}}}&&&\\&&&{\frac {1}{G_{23}}}&&\\&&&&{\frac {1}{G_{13}}}&\\&&&&&{\frac {1}{G_{12}}}\end{bmatrix}} _{=:S}{\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\end{bmatrix}} $
Die Matrix S ist die Nachgiebigkeitsmatrix und ihre Symmetrie erfordert:
- $ {\frac {\nu _{21}}{E_{2}}}={\frac {\nu _{12}}{E_{1}}}\,,\quad {\frac {\nu _{31}}{E_{3}}}={\frac {\nu _{13}}{E_{1}}}\,,\quad {\frac {\nu _{32}}{E_{3}}}={\frac {\nu _{23}}{E_{2}}}\,. $
sodass von den zwölf Einträgen nur neun unabhängig sind.
Invertierung der Nachgiebigkeitsmatrix unter Berücksichtigung ihrer Symmetrie liefert die ebenfalls symmetrische Steifigkeitsmatrix
- $ {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\\\end{bmatrix}}={\begin{bmatrix}{\frac {1-\nu _{23}\nu _{32}}{D}}E_{1}&{\frac {\nu _{21}+\nu _{23}\nu _{31}}{D}}E_{1}&{\frac {\nu _{31}+\nu _{32}\nu _{21}}{D}}E_{1}\\{\frac {\nu _{12}+\nu _{13}\nu _{32}}{D}}E_{2}&{\frac {1-\nu _{13}\nu _{31}}{D}}E_{2}&{\frac {\nu _{32}+\nu _{31}\nu _{12}}{D}}E_{2}\\{\frac {\nu _{13}+\nu _{12}\nu _{23}}{D}}E_{3}&{\frac {\nu _{23}+\nu _{21}\nu _{13}}{D}}E_{3}&{\frac {1-\nu _{12}\nu _{21}}{D}}E_{3}\\&&&G_{23}\\&&&&G_{13}\\&&&&&G_{12}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\end{bmatrix}} $
mit
- $ D=1-\nu _{12}\nu _{21}-\nu _{13}\nu _{31}-\nu _{23}\nu _{32}-2\nu _{12}\nu _{23}\nu _{31}\,. $
Die Nachgiebigkeitsmatrix und die Steifigkeitsmatrix sind an denselben Stellen mit von Null verschiedenen Werten besetzt.
Spezialfälle der Orthotropie
In der kubischen Anisotropie sind die Elastizitäts- und Schubmoduln sowie die Querdehnzahlen alle gleich:
- $ E_{1}=E_{2}=E_{3},\quad G_{12}=G_{13}=G_{23},\quad \nu _{12}=\nu _{21}=\nu _{13}=\nu _{31}=\nu _{23}=\nu _{32} $
womit nur drei unabhängige Elastizitätsparameter übrig bleiben. Transversale Isotropie mit fünf unabhängigen Elastizitätsparametern stellt sich ein mit:
- $ {\begin{aligned}E_{2}&=E_{3}\\G_{12}&=G_{13}\\\nu _{12}&=\nu _{13}\\G_{23}&={\frac {E_{2}}{2(1+\nu _{23})}}\end{aligned}} $
In der Isotropie gelten die Identitäten der kubischen Anisotropie und die drei übrig bleibenden unabhängigen Größen sind zusätzlich durch den letzten Zusammenhang in der transversalen Isotropie verbunden, sodass nur noch zwei unabhängige Elastizitätsparameter übrig bleiben.
Ebener Spannungszustand
In dünnwandigen Strukturen aus orthtropem Material sind zwei der Orthotropieachsen oftmals in den Vorzugsrichtungen der Struktur gelegen, wie zum Beispiel bei Holzplatten, und es liegt ein ebener Spannungszustand vor.
Hier ist σ13=σ23=σ33=0 und aus letzterer Identität leitet sich
- $ \varepsilon _{33}=-{\frac {(\nu _{13}+\nu _{12}\nu _{23})\varepsilon _{11}+(\nu _{23}+\nu _{21}\nu _{13})\varepsilon _{22}}{1-\nu _{12}\nu _{21}}}=-{\frac {1}{E_{3}}}(\nu _{31}\sigma _{11}+\nu _{32}\sigma _{22}) $
ab. Das Elastizitätsgesetz vereinfacht sich zu
- $ {\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\2\varepsilon _{12}\end{bmatrix}}={\begin{bmatrix}{\frac {1}{E_{1}}}&-{\frac {\nu _{21}}{E_{2}}}&0\\-{\frac {\nu _{12}}{E_{1}}}&{\frac {1}{E_{2}}}&0\\0&0&{\frac {1}{G_{12}}}\end{bmatrix}}{\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\\\end{bmatrix}} $
bzw.
- $ {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\\\end{bmatrix}}={\begin{bmatrix}{\frac {E_{1}}{1-\nu _{12}\nu _{21}}}&{\frac {\nu _{21}E_{1}}{1-\nu _{12}\nu _{21}}}&0\\{\frac {\nu _{12}E_{2}}{1-\nu _{12}\nu _{21}}}&{\frac {E_{2}}{1-\nu _{12}\nu _{21}}}&0\\0&0&G_{12}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\2\varepsilon _{12}\end{bmatrix}} $
mit jeweils symmetrischer Nachgiebigkeits- bzw. Steifigkeitsmatrix.
In der linearen orthotropen Elastizität für den Ebenen Spannungszustand werden die Schubmoduln $ G_{13},G_{23} $ nicht gebraucht, sodass nur sieben Materialparameter ausreichen, und wenn nur die Spannungen und Verzerrungen in der Ebene interessieren, sind es nur mehr vier Materialparameter.
Ebener Verzerrungszustand
Hier finden die Verzerrungen ausschließlich in der 1-2-Ebene statt, nur die Normalspannung senkrecht zur Ebene darf auftreten. Aus $ \varepsilon _{33}=0 $ leitet sich
- $ \sigma _{33}=\nu _{31}\sigma _{11}+\nu _{32}\sigma _{22}={\frac {E_{3}}{D}}[(\nu _{13}+\nu _{12}\nu _{23})\varepsilon _{11}+(\nu _{23}+\nu _{21}\nu _{13})\varepsilon _{22}] $
mit
- $ D=1-\nu _{12}\nu _{21}-\nu _{13}\nu _{31}-\nu _{23}\nu _{32}-2\nu _{12}\nu _{23}\nu _{31}\,. $
ab. Das Elastizitätsgesetz reduziert sich auf
- $ {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\\\end{bmatrix}}={\begin{bmatrix}{\frac {1-\nu _{23}\nu _{32}}{D}}E_{1}&{\frac {\nu _{21}+\nu _{23}\nu _{31}}{D}}E_{1}&0\\{\frac {\nu _{12}+\nu _{13}\nu _{32}}{D}}E_{2}&{\frac {1-\nu _{13}\nu _{31}}{D}}E_{2}&0\\0&0&G_{12}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\2\varepsilon _{12}\end{bmatrix}} $
bzw.
- $ {\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\2\varepsilon _{12}\end{bmatrix}}={\begin{bmatrix}{\frac {1-\nu _{13}\nu _{31}}{E_{1}}}&-{\frac {\nu _{21}+\nu _{23}\nu _{31}}{E_{2}}}&0\\-{\frac {\nu _{12}+\nu _{13}\nu _{32}}{E_{1}}}&{\frac {1-\nu _{23}\nu _{32}}{E_{2}}}&0\\0&0&{\frac {1}{G_{12}}}\end{bmatrix}}{\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\\\end{bmatrix}} $
mit jeweils symmetrischer Steifigkeits- bzw. Nachgiebigkeitsmatrix.
In der linearen orthotropen Elastizität für den Ebenen Verzerrungszustand werden die Schubmoduln $ G_{13},G_{23} $ nicht gebraucht, sodass nur sieben Materialparameter ausreichen, und wenn nur die Spannungen und Verzerrungen in der Ebene interessieren entfällt zusätzlich E3, sodass nur mehr sechs Materialparameter gebraucht werden.
Stabilitätskriterien
Die Materialparameter können nicht beliebig gewählt werden, sondern müssen gewissen Stabilitätskriterien genügen. Diese folgen aus der Forderung, dass die Steifigkeits- und Nachgiebigkeitsmatrizen positiv definit sein müssen. Dies führt auf die Bedingungen:
- Alle Diagonalelemente der Steifigkeits- und Nachgiebigkeitsmatrix müssen positiv sein (damit sich das Material in Zugrichtung streckt, wenn man daran zieht, und nicht staucht) und
- die Determinante der Steifigkeits- und Nachgiebigkeitsmatrix muss positiv sein (damit es unter Druck komprimiert und nicht expandiert).
Werden an einem realen Werkstoff Materialparameter identifiziert, die diesen Stabilitätskriterien widersprechen, ist Vorsicht geboten. Die Stabilitätskriterien lauten:[2]
- $ {\begin{array}{l}E_{1},E_{2},E_{3},G_{12},G_{13},G_{23}>0\\|\nu _{12}|<{\sqrt {\dfrac {E_{1}}{E_{2}}}}\quad \rightarrow \quad 1-\nu _{12}\nu _{21}>0\\|\nu _{13}|<{\sqrt {\dfrac {E_{1}}{E_{3}}}}\quad \rightarrow \quad 1-\nu _{13}\nu _{31}>0\\|\nu _{23}|<{\sqrt {\dfrac {E_{2}}{E_{3}}}}\quad \rightarrow \quad 1-\nu _{23}\nu _{32}>0\\1-\nu _{12}\nu _{21}-\nu _{13}\nu _{31}-\nu _{23}\nu _{32}-2\nu _{12}\nu _{23}\nu _{31}>0\end{array}} $
Wenn die linke Seite der letzten Ungleichung gegen null geht, setzt das Material einer hydrostatischen Kompression zunehmend Widerstand entgegen. Aus der Symmetrie der Nachgiebigkeitsmatrix folgt ergänzend:
- $ {\begin{array}{l}|\nu _{21}|<{\sqrt {\dfrac {E_{2}}{E_{1}}}}\\|\nu _{31}|<{\sqrt {\dfrac {E_{3}}{E_{1}}}}\\|\nu _{32}|<{\sqrt {\dfrac {E_{3}}{E_{2}}}}\end{array}} $
Herleitung
In der Hyperelastizität ergeben sich die Spannungen aus der Ableitung der Formänderungsenergie nach den Dehnungen. Damit die Spannungen linear in den Dehnungen sind, muss demnach die Formänderungsenergie quadratisch in den Dehnungen sein, denn nur dann ist ihre Ableitung linear. Unter Verwendung der #Invarianten Terme ergibt sich der Ansatz
- $ {\begin{aligned}w({\boldsymbol {\varepsilon }})=&{\frac {a}{2}}\varepsilon _{11}^{2}+{\frac {b}{2}}\varepsilon _{22}^{2}+{\frac {c}{2}}\varepsilon _{33}^{2}+d\varepsilon _{11}\varepsilon _{22}+e\varepsilon _{11}\varepsilon _{33}+f\varepsilon _{22}\varepsilon _{33}\\&+2g\varepsilon _{23}^{2}+2h\varepsilon _{13}^{2}+2p\varepsilon _{12}^{2}\end{aligned}} $
mit neun Parametern a bis p. Um dies nach ε ableiten zu können, müssen die Komponenten εij als Funktion des Tensors ε ausgedrückt werden. Dies gelingt mit der Darstellung des Frobenius-Skalarprodukts ":" als Spur:
- $ \mathbf {A} :\mathbf {B} :=\mathrm {Spur} (\mathbf {A^{\top }\cdot B} ) $
Mit der Abkürzung $ \mathbf {K} _{ij}={\frac {1}{2}}({\hat {e}}_{i}\otimes {\hat {e}}_{j}+{\hat {e}}_{j}\otimes {\hat {e}}_{i}) $ für die symmetrisierten dyadischen Produkte der Orthotropieachsenvektoren ist dann[3]
- $ \mathbf {K} _{ij}:{\boldsymbol {\varepsilon }}={\frac {1}{2}}(\varepsilon _{ij}+\varepsilon _{ji})=\varepsilon _{ij}\quad \rightarrow \quad {\frac {\mathrm {d} \varepsilon _{ij}}{\mathrm {d} {\boldsymbol {\varepsilon }}}}=\mathbf {K} _{ij} $
Aus dem Ansatz der Formänderungsenergie berechnen sich die Spannungen zu
- $ {\begin{aligned}{\boldsymbol {\sigma }}={\frac {\mathrm {d} w}{\mathrm {d} {\boldsymbol {\varepsilon }}}}=&a\varepsilon _{11}\mathbf {K} _{11}+b\varepsilon _{22}\mathbf {K} _{22}+c\varepsilon _{33}\mathbf {K} _{33}\\&+d(\varepsilon _{22}\mathbf {K} _{11}+\varepsilon _{11}\mathbf {K} _{22})+e(\varepsilon _{33}\mathbf {K} _{11}+\varepsilon _{11}\mathbf {K} _{33})+f(\varepsilon _{22}\mathbf {K} _{33}+\varepsilon _{33}\mathbf {K} _{22})\\&+4g\varepsilon _{23}\mathbf {K} _{23}+4h\varepsilon _{13}\mathbf {K} _{13}+4p\varepsilon _{12}\mathbf {K} _{12}\end{aligned}} $
oder in Voigt-Notation im ê1,2,3-System
- $ {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\end{bmatrix}}={\begin{bmatrix}a\varepsilon _{11}+d\varepsilon _{22}+e\varepsilon _{33}\\d\varepsilon _{11}+b\varepsilon _{22}+f\varepsilon _{33}\\e\varepsilon _{11}+f\varepsilon _{22}+c\varepsilon _{33}\\2g\varepsilon _{23}\\2h\varepsilon _{13}\\2p\varepsilon _{12}\end{bmatrix}}={\begin{bmatrix}a&d&e&&&\\d&b&f&&&\\e&f&c&&&\\&&&g&&\\&&&&h&\\&&&&&p\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\end{bmatrix}} $
Im #Elastizitätsgesetz für 3D lassen sich die Parameter direkt ablesen. Ableitung der Spannungen nach den Dehnungen liefert den konstanten und symmetrischen Elastizitätstensor 4. Stufe:
- $ {\begin{aligned}\mathbb {C} :={\frac {\mathrm {d} {\boldsymbol {\sigma }}}{\mathrm {d} {\boldsymbol {\varepsilon }}}}=&a\mathbf {K} _{11}\otimes \mathbf {K} _{11}+b\mathbf {K} _{22}\otimes \mathbf {K} _{22}+c\mathbf {K} _{33}\otimes \mathbf {K} _{33}+d(\mathbf {K} _{11}\otimes \mathbf {K} _{22}+\mathbf {K} _{22}\otimes \mathbf {K} _{11})\\&+e(\mathbf {K} _{11}\otimes \mathbf {K} _{33}+\mathbf {K} _{33}\otimes \mathbf {K} _{11})+f(\mathbf {K} _{22}\otimes \mathbf {K} _{33}+\mathbf {K} _{33}\otimes \mathbf {K} _{22})\\&+4g\mathbf {K} _{23}\otimes \mathbf {K} _{23}+4h\mathbf {K} _{13}\otimes \mathbf {K} _{13}+4p\mathbf {K} _{12}\otimes \mathbf {K} _{12}\end{aligned}} $
Die Tensoren Kii werden Strukturvariable genannt, weil sie die interne Struktur des Materials repräsentieren[1]:387 und mit ihnen auch die invarianten Terme dargestellt werden können. Nicht-linear hyperelastisches Verhalten kann modelliert werden, indem
- die Parameter a bis p durch Funktionen der invarianten Terme ersetzt werden, siehe Hyperelastizität#Orthotrope Hyperelastizität, und/oder
- die invarianten Terme höherer Ordnung im Ansatz zur Formänderungsenergie berücksichtigt werden.[1]:394
Gründe für die Besetztheit der Steifigkeitsmatrix
In diesem Abschnitt wird die Frage geklärt, warum die Steifigkeitsmatrix nur an den entsprechenden Stellen besetzt ist. Im Allgemeinen tauchen in einem linearen Materialgesetz 21 unabhängige Materialkonstanten auf. Im Fall der Orthotropie reduziert sich aber die Zahl der Konstanten auf 9. Warum das so ist, ist nachfolgend dargestellt.
Drehmatrizen bei 180-Grad-Drehungen
Die (linearen) Abbildungen, die 180-Grad-Drehungen um die Orthotropieachsen beschreiben, lassen sich mit Matrizen beschreiben. Wählt man als Bezug eine Basis, deren Basisvektoren sich mit den senkrecht aufeinanderstehenden Drehachsen decken, dann haben diese orthogonalen Matrizen folgende Gestalt
- $ {\begin{aligned}A_{x}={\begin{bmatrix}1&0&0\\0&-1&0\\0&0&-1\\\end{bmatrix}},&&A_{y}={\begin{bmatrix}-1&0&0\\0&1&0\\0&0&-1\end{bmatrix}},&&A_{z}={\begin{bmatrix}-1&0&0\\0&-1&0\\0&0&1\end{bmatrix}},\end{aligned}} $
Diese 3 Matrizen bilden eine echte Untergruppe der Drehgruppe SO(3). Das Produkt dieser drei Matrizen ist die Einheitsmatrix: $ A_{x}A_{y}A_{z}=E $.
Die 3 Matrizen Ax,y,z und zusätzlich die negative Einheitsmatrix -E, die eine Punktspiegelung repräsentiert, bilden die Symmetriegruppe des orthotropen Materials[1]:382. Die Symmetriegruppe eines anisotropen Materials ohne -E ist immer eine echte Untergruppe der Drehgruppe SO(3); SO(3) mit -E ist die Symmetriegruppe eines isotropen Materials.
Symmetriebedingung in Indexschreibweise und Voigt’scher Notation
Gedankenexperiment: Ein Teilchen und dessen Umgebung wird einer bestimmten Deformation unterzogen und damit einem bestimmten Verzerrungstensor $ \varepsilon $. Im einfachsten Fall (der allerdings zur Definition der Orthotropie nicht ausreichend allgemein ist) könnte das Teilchen nur in einer bestimmten Richtung gestreckt werden. Nun ändert man die Streckungsrichtung aktiv. Das heißt, man lässt den materiellen Punkt wie er ist (dreht also das Material nicht) und unterzieht den Punkt aber (derselben) Streckung in anderer Richtung. Man gelangt damit zu einem anderen Verzerrungstensor $ \varepsilon ' $.
Die Änderung der Verzerrungsrichtung kann mit einer Drehmatrix $ A $ beschrieben werden. Es gilt
- $ \varepsilon '=A\,\varepsilon \,A^{-1} $
Mithilfe eines linearen Materialgesetzes $ f_{C} $ lässt sich für gegebenen Verzerrungstensor der zugehörige Spannungstensor ermitteln. Es sei
- $ {\begin{aligned}\sigma &:=f_{C}(\varepsilon )\\\sigma '&:=f_{C}(\varepsilon ')\end{aligned}} $
Im allgemeinen Fall der Anisotropie gilt zwar nicht
- $ \sigma '=A\,\sigma \,A^{-1} $
Aber genau dies fordert man für die oben beschriebene Teilmenge von SO(3) im Fall der Orthotropie: Ein Material heißt orthotrop, wenn für die Funktion $ f_{C} $ folgende Symmetrietransformation für jede der oben genannten (orthogonalen) Drehmatrizen und für beliebige Verzerrungen gilt
- $ {\begin{aligned}Af_{C}(\varepsilon )A^{-1}&=f_{C}(A\varepsilon A^{-1})\Leftrightarrow Af_{C}(\varepsilon )A^{T}=f_{C}(A\varepsilon A^{T})\end{aligned}} $
In Indexschreibweise
- $ {\begin{aligned}\sigma _{mn}'=A_{mo}C_{opjk}\varepsilon _{jk}A_{np}&=C_{mnil}\varepsilon '_{il}=C_{mnil}A_{ij}\varepsilon _{jk}A_{lk}\end{aligned}} $
Nun dieselbe Bedingung in Voigt’scher Notation: Mit der Definition
- $ {\begin{aligned}A_{\sigma }^{\text{v}}:={\begin{bmatrix}A_{11}A_{11}&A_{12}A_{12}&A_{13}A_{13}&A_{12}A_{13}+A_{13}A_{12}&A_{11}A_{13}+A_{13}A_{11}&A_{11}A_{12}+A_{12}A_{11}\\A_{21}A_{21}&A_{22}A_{22}&A_{23}A_{23}&A_{22}A_{23}+A_{23}A_{22}&A_{21}A_{23}+A_{23}A_{21}&A_{21}A_{22}+A_{22}A_{21}\\A_{31}A_{31}&A_{32}A_{32}&A_{33}A_{33}&A_{32}A_{33}+A_{33}A_{32}&A_{31}A_{33}+A_{33}A_{31}&A_{31}A_{32}+A_{32}A_{31}\\A_{21}A_{31}&A_{22}A_{32}&A_{23}A_{33}&A_{22}A_{33}+A_{23}A_{32}&A_{21}A_{33}+A_{23}A_{31}&A_{21}A_{32}+A_{22}A_{31}\\A_{11}A_{31}&A_{12}A_{32}&A_{13}A_{33}&A_{12}A_{33}+A_{13}A_{32}&A_{11}A_{33}+A_{13}A_{31}&A_{11}A_{32}+A_{12}A_{31}\\A_{11}A_{21}&A_{12}A_{22}&A_{13}A_{23}&A_{12}A_{23}+A_{13}A_{22}&A_{11}A_{23}+A_{13}A_{21}&A_{11}A_{22}+A_{12}A_{21}\\\end{bmatrix}}\end{aligned}} $
gilt
- $ {\begin{aligned}{\begin{bmatrix}\sigma '_{11}\\\sigma '_{22}\\\sigma '_{33}\\\sigma '_{23}\\\sigma '_{13}\\\sigma '_{12}\end{bmatrix}}=A_{\sigma }^{\text{v}}{\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\end{bmatrix}}\Leftrightarrow {\sigma '}^{\text{v}}=A_{\sigma }^{\text{v}}{\sigma }^{\text{v}},\qquad \qquad {\begin{bmatrix}\varepsilon '_{11}\\\varepsilon '_{22}\\\varepsilon '_{33}\\\varepsilon '_{23}\\\varepsilon '_{13}\\\varepsilon '_{12}\end{bmatrix}}=A_{\sigma }^{\text{v}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\\varepsilon _{23}\\\varepsilon _{13}\\\varepsilon _{12}\end{bmatrix}}\end{aligned}} $
Mit der neuen Definition
- $ {\begin{aligned}A_{\varepsilon }^{\text{v}}:={\begin{bmatrix}A_{11}A_{11}&A_{12}A_{12}&A_{13}A_{13}&A_{12}A_{13}+A_{13}A_{12}&A_{11}A_{13}+A_{13}A_{11}&A_{11}A_{12}+A_{12}A_{11}\\A_{21}A_{21}&A_{22}A_{22}&A_{23}A_{23}&A_{22}A_{23}+A_{23}A_{22}&A_{21}A_{23}+A_{23}A_{21}&A_{21}A_{22}+A_{22}A_{21}\\A_{31}A_{31}&A_{32}A_{32}&A_{33}A_{33}&A_{32}A_{33}+A_{33}A_{32}&A_{31}A_{33}+A_{33}A_{31}&A_{31}A_{32}+A_{32}A_{31}\\2A_{21}A_{31}&2A_{22}A_{32}&2A_{23}A_{33}&A_{22}A_{33}+A_{23}A_{32}&A_{21}A_{33}+A_{23}A_{31}&A_{21}A_{32}+A_{22}A_{31}\\2A_{11}A_{31}&2A_{12}A_{32}&2A_{13}A_{33}&A_{12}A_{33}+A_{13}A_{32}&A_{11}A_{33}+A_{13}A_{31}&A_{11}A_{32}+A_{12}A_{31}\\2A_{11}A_{21}&2A_{12}A_{22}&2A_{13}A_{23}&A_{12}A_{23}+A_{13}A_{22}&A_{11}A_{23}+A_{13}A_{21}&A_{11}A_{22}+A_{12}A_{21}\\\end{bmatrix}}\end{aligned}} $
ergibt sich
- $ {\begin{aligned}{\begin{bmatrix}\varepsilon '_{11}\\\varepsilon '_{22}\\\varepsilon '_{33}\\2\varepsilon '_{23}\\2\varepsilon '_{13}\\2\varepsilon '_{12}\end{bmatrix}}=A_{\varepsilon }^{\text{v}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{13}\\2\varepsilon _{12}\end{bmatrix}}\Leftrightarrow {\varepsilon '}^{\text{v}}=A_{\varepsilon }^{\text{v}}{\varepsilon }^{\text{v}}\end{aligned}} $
In Voigt’scher Notation erhält man also als Symmetriebedingung
- $ {\begin{aligned}{\sigma '}^{\text{v}}=A_{\sigma }^{\text{v}}{\sigma }^{\text{v}}&=A_{\sigma }^{\text{v}}C^{\text{v}}{\varepsilon }^{\text{v}}=C^{\text{v}}{\varepsilon '}^{\text{v}}=C^{\text{v}}A_{\varepsilon }^{\text{v}}{\varepsilon }^{\text{v}}\end{aligned}} $
Und da dies für beliebige Dehnungen gelten muss, ist die Symmetriebedingung
- $ {\begin{aligned}A_{\sigma }^{\text{v}}C^{\text{v}}=C^{\text{v}}A_{\varepsilon }^{\text{v}}\end{aligned}} $
Spezialfall 180-Grad-Drehungen
Da im Spezialfall der Orthotropie die 3×3-Matrizen $ A $ nur auf der Hauptdiagonalen besetzt sind, vereinfachen sich die Definitionen von oben zu:
$ {\begin{aligned}A_{\sigma }^{\text{v}}=A_{\varepsilon }^{\text{v}}&={\begin{bmatrix}A_{11}A_{11}&0&0&0&0&0\\0&A_{22}A_{22}&0&0&0&0\\0&0&A_{33}A_{33}&0&0&0\\0&0&0&A_{22}A_{33}&0&0\\0&0&0&0&A_{11}A_{33}&0\\0&0&0&0&0&A_{11}A_{22}\\\end{bmatrix}}\\\end{aligned}} $
Die drei 3×3-Matrizen entsprechen also den drei 6x6-Matrizen
- $ {\begin{aligned}A_{x}^{\text{v}}&={\begin{bmatrix}1&0&0&0&0&0\\0&1&0&0&0&0\\0&0&1&0&0&0\\0&0&0&-1&0&0\\0&0&0&0&-1&0\\0&0&0&0&0&1\\\end{bmatrix}},\qquad A_{y}^{\text{v}}={\begin{bmatrix}1&0&0&0&0&0\\0&1&0&0&0&0\\0&0&1&0&0&0\\0&0&0&-1&0&0\\0&0&0&0&1&0\\0&0&0&0&0&-1\\\end{bmatrix}},\qquad A_{z}^{\text{v}}={\begin{bmatrix}1&0&0&0&0&0\\0&1&0&0&0&0\\0&0&1&0&0&0\\0&0&0&1&0&0\\0&0&0&0&-1&0\\0&0&0&0&0&-1\\\end{bmatrix}}\end{aligned}} $
Auswertung der Symmetriebedingungen für den Spezialfall
Die Symmetriebedingung ausgewertet für diese Matrizen ergibt
- $ {\begin{aligned}{\text{wegen }}A_{x}^{\text{v}}C^{\text{v}}=C^{\text{v}}A_{x}^{\text{v}}:\qquad {\begin{bmatrix}C_{11}^{\text{v}}&C_{12}^{\text{v}}&C_{13}^{\text{v}}&C_{14}^{\text{v}}&C_{15}^{\text{v}}&C_{16}^{\text{v}}\\C_{21}^{\text{v}}&C_{22}^{\text{v}}&C_{23}^{\text{v}}&C_{24}^{\text{v}}&C_{25}^{\text{v}}&C_{26}^{\text{v}}\\C_{31}^{\text{v}}&C_{32}^{\text{v}}&C_{33}^{\text{v}}&C_{34}^{\text{v}}&C_{35}^{\text{v}}&C_{36}^{\text{v}}\\-C_{41}^{\text{v}}&-C_{42}^{\text{v}}&-C_{43}^{\text{v}}&-C_{44}^{\text{v}}&-C_{45}^{\text{v}}&-C_{46}^{\text{v}}\\-C_{51}^{\text{v}}&-C_{52}^{\text{v}}&-C_{53}^{\text{v}}&-C_{54}^{\text{v}}&-C_{55}^{\text{v}}&-C_{56}^{\text{v}}\\C_{61}^{\text{v}}&C_{62}^{\text{v}}&C_{63}^{\text{v}}&C_{64}^{\text{v}}&C_{65}^{\text{v}}&C_{66}^{\text{v}}\end{bmatrix}}&={\begin{bmatrix}C_{11}^{\text{v}}&C_{12}^{\text{v}}&C_{13}^{\text{v}}&-C_{14}^{\text{v}}&-C_{15}^{\text{v}}&C_{16}^{\text{v}}\\C_{21}^{\text{v}}&C_{22}^{\text{v}}&C_{23}^{\text{v}}&-C_{24}^{\text{v}}&-C_{25}^{\text{v}}&C_{26}^{\text{v}}\\C_{31}^{\text{v}}&C_{32}^{\text{v}}&C_{33}^{\text{v}}&-C_{34}^{\text{v}}&-C_{35}^{\text{v}}&C_{36}^{\text{v}}\\C_{41}^{\text{v}}&C_{42}^{\text{v}}&C_{43}^{\text{v}}&-C_{44}^{\text{v}}&-C_{45}^{\text{v}}&C_{46}^{\text{v}}\\C_{51}^{\text{v}}&C_{52}^{\text{v}}&C_{53}^{\text{v}}&-C_{54}^{\text{v}}&-C_{55}^{\text{v}}&C_{56}^{\text{v}}\\C_{61}^{\text{v}}&C_{62}^{\text{v}}&C_{63}^{\text{v}}&-C_{64}^{\text{v}}&-C_{65}^{\text{v}}&C_{66}^{\text{v}}\end{bmatrix}}\\{\text{wegen }}A_{y}^{\text{v}}C^{\text{v}}=C^{\text{v}}A_{y}^{\text{v}}:\qquad {\begin{bmatrix}C_{11}^{\text{v}}&C_{12}^{\text{v}}&C_{13}^{\text{v}}&C_{14}^{\text{v}}&C_{15}^{\text{v}}&C_{16}^{\text{v}}\\C_{21}^{\text{v}}&C_{22}^{\text{v}}&C_{23}^{\text{v}}&C_{24}^{\text{v}}&C_{25}^{\text{v}}&C_{26}^{\text{v}}\\C_{31}^{\text{v}}&C_{32}^{\text{v}}&C_{33}^{\text{v}}&C_{34}^{\text{v}}&C_{35}^{\text{v}}&C_{36}^{\text{v}}\\-C_{41}^{\text{v}}&-C_{42}^{\text{v}}&-C_{43}^{\text{v}}&-C_{44}^{\text{v}}&-C_{45}^{\text{v}}&-C_{46}^{\text{v}}\\C_{51}^{\text{v}}&C_{52}^{\text{v}}&C_{53}^{\text{v}}&C_{54}^{\text{v}}&C_{55}^{\text{v}}&C_{56}^{\text{v}}\\-C_{61}^{\text{v}}&-C_{62}^{\text{v}}&-C_{63}^{\text{v}}&-C_{64}^{\text{v}}&-C_{65}^{\text{v}}&-C_{66}^{\text{v}}\end{bmatrix}}&={\begin{bmatrix}C_{11}^{\text{v}}&C_{12}^{\text{v}}&C_{13}^{\text{v}}&-C_{14}^{\text{v}}&C_{15}^{\text{v}}&-C_{16}^{\text{v}}\\C_{21}^{\text{v}}&C_{22}^{\text{v}}&C_{23}^{\text{v}}&-C_{24}^{\text{v}}&C_{25}^{\text{v}}&-C_{26}^{\text{v}}\\C_{31}^{\text{v}}&C_{32}^{\text{v}}&C_{33}^{\text{v}}&-C_{34}^{\text{v}}&C_{35}^{\text{v}}&-C_{36}^{\text{v}}\\C_{41}^{\text{v}}&C_{42}^{\text{v}}&C_{43}^{\text{v}}&-C_{44}^{\text{v}}&C_{45}^{\text{v}}&-C_{46}^{\text{v}}\\C_{51}^{\text{v}}&C_{52}^{\text{v}}&C_{53}^{\text{v}}&-C_{54}^{\text{v}}&C_{55}^{\text{v}}&-C_{56}^{\text{v}}\\C_{61}^{\text{v}}&C_{62}^{\text{v}}&C_{63}^{\text{v}}&-C_{64}^{\text{v}}&C_{65}^{\text{v}}&-C_{66}^{\text{v}}\end{bmatrix}}\\{\text{wegen }}A_{z}^{\text{v}}C^{\text{v}}=C^{\text{v}}A_{z}^{\text{v}}:\qquad {\begin{bmatrix}C_{11}^{\text{v}}&C_{12}^{\text{v}}&C_{13}^{\text{v}}&C_{14}^{\text{v}}&C_{15}^{\text{v}}&C_{16}^{\text{v}}\\C_{21}^{\text{v}}&C_{22}^{\text{v}}&C_{23}^{\text{v}}&C_{24}^{\text{v}}&C_{25}^{\text{v}}&C_{26}^{\text{v}}\\C_{31}^{\text{v}}&C_{32}^{\text{v}}&C_{33}^{\text{v}}&C_{34}^{\text{v}}&C_{35}^{\text{v}}&C_{36}^{\text{v}}\\C_{41}^{\text{v}}&C_{42}^{\text{v}}&C_{43}^{\text{v}}&C_{44}^{\text{v}}&C_{45}^{\text{v}}&C_{46}^{\text{v}}\\-C_{51}^{\text{v}}&-C_{52}^{\text{v}}&-C_{53}^{\text{v}}&-C_{54}^{\text{v}}&-C_{55}^{\text{v}}&-C_{56}^{\text{v}}\\-C_{61}^{\text{v}}&-C_{62}^{\text{v}}&-C_{63}^{\text{v}}&-C_{64}^{\text{v}}&-C_{65}^{\text{v}}&-C_{66}^{\text{v}}\end{bmatrix}}&={\begin{bmatrix}C_{11}^{\text{v}}&C_{12}^{\text{v}}&C_{13}^{\text{v}}&C_{14}^{\text{v}}&-C_{15}^{\text{v}}&-C_{16}^{\text{v}}\\C_{21}^{\text{v}}&C_{22}^{\text{v}}&C_{23}^{\text{v}}&C_{24}^{\text{v}}&-C_{25}^{\text{v}}&-C_{26}^{\text{v}}\\C_{31}^{\text{v}}&C_{32}^{\text{v}}&C_{33}^{\text{v}}&C_{34}^{\text{v}}&-C_{35}^{\text{v}}&-C_{36}^{\text{v}}\\C_{41}^{\text{v}}&C_{42}^{\text{v}}&C_{43}^{\text{v}}&C_{44}^{\text{v}}&-C_{45}^{\text{v}}&-C_{46}^{\text{v}}\\C_{51}^{\text{v}}&C_{52}^{\text{v}}&C_{53}^{\text{v}}&C_{54}^{\text{v}}&-C_{55}^{\text{v}}&-C_{56}^{\text{v}}\\C_{61}^{\text{v}}&C_{62}^{\text{v}}&C_{63}^{\text{v}}&C_{64}^{\text{v}}&-C_{65}^{\text{v}}&-C_{66}^{\text{v}}\end{bmatrix}}\end{aligned}} $
An den letzten 3 Gleichungen erkennt man, dass $ C $ nur folgende Gestalt haben kann
- $ {\begin{aligned}C^{\text{v}}&={\begin{bmatrix}C_{11}^{\text{v}}&C_{12}^{\text{v}}&C_{13}^{\text{v}}&0&0&0\\C_{21}^{\text{v}}&C_{22}^{\text{v}}&C_{23}^{\text{v}}&0&0&0\\C_{31}^{\text{v}}&C_{32}^{\text{v}}&C_{33}^{\text{v}}&0&0&0\\0&0&0&C_{44}^{\text{v}}&0&0\\0&0&0&0&C_{55}^{\text{v}}&0\\0&0&0&0&0&C_{66}^{\text{v}}\end{bmatrix}}\end{aligned}} $
Da diese Voigt’sche Steifigkeitsmatrix außerdem symmetrisch ist, bleibt
- $ {\begin{aligned}C^{\text{v}}&={\begin{bmatrix}C_{11}^{\text{v}}&C_{12}^{\text{v}}&C_{13}^{\text{v}}&0&0&0\\&C_{22}^{\text{v}}&C_{23}^{\text{v}}&0&0&0\\&&C_{33}^{\text{v}}&0&0&0\\&&&C_{44}^{\text{v}}&0&0\\&{\text{sym}}&&&C_{55}^{\text{v}}&0\\&&&&&C_{66}^{\text{v}}\end{bmatrix}}\end{aligned}} $
Zusammenfassung
- Die Orthotropie in der linearen Elastizitätstheorie lässt sich definieren als ein Spezialfall der Anisotropie, bei dem die Steifigkeits- oder Nachgiebigkeitsmatrix eine besonders einfache Form annimmt (9 Konstanten anstelle von 21 Konstanten im allgemeinen Fall).
- Neben der Orthotropie gibt es noch andere Spezialfälle der Anisotropie, z. B. Transversalisotropie, Isotropie etc. Hierbei werden dieselben Symmetriebedingungen angegeben. Nur werden dann andere Untergruppen der Drehgruppe (also andere Matrizen $ A $) betrachtet.
- An der Form des elastischen Gesetzes erkennt man, dass die Kopplung zwischen Zug und Schub für Belastung entlang der Orthotropierichtungen entfällt.
Siehe auch
Weblinks
Einzelnachweise
- ↑ 1,00 1,01 1,02 1,03 1,04 1,05 1,06 1,07 1,08 1,09 1,10 P. Haupt: Continuum Mechanics and Theory of Materials. Springer, 2002, ISBN 978-3-642-07718-0.
- ↑ Holm Altenbach: Kontinuumsmechanik. Einführung in die materialunabhängigen und materialabhängigen Gleichungen. Springer-Verlag, Berlin/Heidelberg 2012, ISBN 978-3-642-24118-5, S. 331, doi:10.1007/978-3-642-24119-2.
- ↑ Die ij-Komponente eines beliebigen Tensors zweiter Stufe T im ê1,2,3-System ist
- $ {\begin{aligned}T_{ij}:=&{\hat {e}}_{i}\cdot \mathbf {T} \cdot {\hat {e}}_{j}={\hat {e}}_{j}\cdot ({\hat {e}}_{i}\cdot \mathbf {T} )\\=&\mathrm {Spur} ({\hat {e}}_{j}\otimes {\hat {e}}_{i}\cdot \mathbf {T} ):=({\hat {e}}_{i}\otimes {\hat {e}}_{j}):\mathbf {T} \end{aligned}} $
Die Fréchet-Ableitung hiervon nach T ist der beschränkte lineare Operator $ {\mathcal {A}} $ der – sofern er existiert – in allen Richtungen H dem Gâteaux-Differential entspricht, also
- $ {\begin{aligned}{\mathcal {A}}(\mathbf {H} )=&\left.{\frac {\mathrm {d} }{\mathrm {d} s}}[({\hat {e}}_{i}\otimes {\hat {e}}_{j}):(\mathbf {T} +s\mathbf {H} )]\right|_{s=0}\\=&({\hat {e}}_{i}\otimes {\hat {e}}_{j}):\mathbf {H} \quad \forall \;\mathbf {H} \end{aligned}} $
Darin ist $ s\in \mathbb {R} $ und der lineare Operator ist das Skalarprodukt mit $ {\hat {e}}_{i}\otimes {\hat {e}}_{j} $. Hier ist $ \mathbf {T} ={\boldsymbol {\varepsilon }} $ ein symmetrischer Tensor, dessen Differential H auch ein symmetrischer Tensor ist. Beim Skalarprodukt mit diesem trägt nur der symmetrische Anteil etwas bei:
- $ ({\hat {e}}_{i}\otimes {\hat {e}}_{j}):\mathbf {H} ={\frac {1}{2}}({\hat {e}}_{i}\otimes {\hat {e}}_{j}+{\hat {e}}_{j}\otimes {\hat {e}}_{i}):\mathbf {H} =\mathbf {K} _{ij}:\mathbf {H} $
Dann wird auch
- $ {\mathcal {A}}={\frac {\partial T_{ij}}{\partial \mathbf {T} }}=\mathbf {K} _{ij} $
geschrieben.
Literatur
- J.F. Nye: Physical Properties of Crystals: Their Representation by Tensors and Matrices. Oxford University Press, 1985, ISBN 978-0-19-851165-6.
- H. Altenbach, J. Altenbach, R. Rikards: Einführung in die Mechanik der Laminat- und Sandwichtragwerke. Deutscher Verlag für Grundstoffindustrie, Stuttgart 1996, ISBN 3-342-00681-1.