imported>Koyaanisqatsi01 K (→Einleitung: Tippfehler entfernt) |
imported>Alva2004 (Ableitung präzisiert) |
||
Zeile 1: | Zeile 1: | ||
[[ | [[Datei:Holz coordinateSystem.svg|mini|Das Koordinatensystem mit den drei Orthotropieachsen Radial, Transversal, Longitudinal]] | ||
[[ | [[Datei:Softwood endgrain marked.png|mini|Holz als typisches orthotropes Material im Ingenieurwesen]] | ||
[[Datei:orthotropie.png|mini| | [[Datei:orthotropie.png|mini|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 [[Griechische Sprache|griechisch]] {{lang|grc|ορθός}} {{lang|grc-Latn|orthos}} „korrekt, senkrecht, gerade stehend“ und {{lang|grc|τρόπος}} {{lang|grc-Latn|tropos}} „Weg, Art und Weise“) ist eine spezielle Art der Richtungsabhängigkeit eines Materials. Orthotrope Materialien wie im Bild haben die folgenden Eigenschaften: | Die '''Orthotropie''' (von [[Griechische Sprache|griechisch]] {{lang|grc|ορθός}} {{lang|grc-Latn|orthos}} „korrekt, senkrecht, gerade stehend“ und {{lang|grc|τρόπος}} {{lang|grc-Latn|tropos}} „Weg, Art und Weise“) ist eine spezielle Art der Richtungsabhängigkeit eines [[Werkstoff]]s/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. | # 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. | # Im Bezugssystem parallel zu den Orthotropieachsen gibt es keine Kopplung zwischen [[Dehnung|Normaldehnungen]] und [[Scherung (Mechanik)|Schubverzerrungen]]. | ||
Ein linear elastisches orthotropes Material besitzt maximal neun Materialparameter. | Ein linear elastisches orthotropes Material besitzt maximal neun Materialparameter. | ||
Ein Material ist [[Isotropie|isotrop]], wenn es richtungsunabhängig dasselbe Kraft-Verformungs-Verhalten hat. Bei [[Anisotropie| | Ein Material ist [[Isotropie|isotrop]], wenn es richtungsunabhängig dasselbe Kraft-Verformungs-Verhalten hat. Bei [[Anisotropie|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. | Viele Konstruktionswerkstoffe sind orthotrop, z. B. technisches Holz, Gewebe, viele Faser-Kunststoff-Verbunde und Walzbleche mit Textur. Kristalle des [[Rhombisches Kristallsystem|rhombischen Kristallsystems]] sind orthotrop<ref name="haupt">{{Literatur | ||
|Autor=P. Haupt | |||
|Titel=Continuum Mechanics and Theory of Materials | |||
|Verlag=Springer | |||
|Datum=2002 | |||
|ISBN=978-3-642-07718-0}}</ref>{{rp|390}}, Spezialfälle kommen im [[Tetragonales Kristallsystem|tetragonalen]]<ref name="haupt"/>{{rp|391}}, [[Hexagonales Kristallsystem|hexagonalen]]<ref name="haupt"/>{{rp|393}} und [[Kubische Anisotropie|kubischen Kristallsystem]] vor. | |||
== Symmetriegruppe == | == 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 die Symmetriegruppe des orthotropen Materials<ref | 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.<ref name="haupt"/>{{rp|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 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 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 [[Orthogonaler Tensor|orthogonale Tensoren]] '''Q''' repräsentiert. Eine Symmetriegruppe ''g<sub>R</sub>'' besteht aus denjenigen Transformationen, die die [[Formänderungsenergie]] ''e'' invariant lassen. Mathematisch wird das mit dem [[Verzerrungstensor]] '''E''' durch | |||
:<math>\mathbf Q\in g_R | |||
\quad\leftrightarrow\quad | |||
e(\mathbf{Q\cdot E\cdot Q}^\top)=e(\mathbf E) | |||
</math> für alle '''E''' | |||
ausgedrückt.<ref name="haupt"/>{{rp|379}} Darin bedeutet „·“ das [[Matrizenprodukt]] und das hochgestellte „⊤“ eine [[Transponierte Matrix|Transponierung]]. Mit '''Q''' gehört auch -'''Q''' zur Symmetriegruppe, was durch Hinzufügen des negativen Einheitstensors -'''1''', der eine [[Punktspiegelung]] repräsentiert, zu ''g<sub>R</sub>'' berücksichtigt wird. Die Symmetriegruppe des orthotropen Materials ist<ref name="haupt"/>{{rp|382}} | |||
:<math>g_R=\left\{-\mathbf 1, \mathbf{Q}_1^\pi, \mathbf{Q}_2^\pi\right\} | |||
</math> | |||
Darin steht <math>\mathbf{Q}_i^\alpha</math> für den [[Orthogonaler Tensor|orthogonalen Tensor]], der mit dem Winkel ''α'' in [[Radiant (Einheit)|Radiant]] um die <math>i</math>-te Orthotropieachse dreht. Die 180-Grad-Drehung um die 3-Achse ist in ''g<sub>R</sub>'' enthalten, denn | |||
:<math>\mathbf{Q}_1^\pi\cdot\mathbf{Q}_2^\pi\cdot\mathbf{Q}_3^\pi | |||
=\mathbf 1=(-\mathbf 1)^2\in g_R</math> | |||
== Invarianten == | |||
In der isotropen [[Hyperelastizität]] hängt die Formänderungsenergie von den [[Hauptinvariante]]n I<sub>1,2,3</sub> des Verzerrungstensors '''E''' ab: | |||
:e('''E''')=e(I<sub>1</sub>, I<sub>2</sub>, I<sub>3</sub>) | |||
Die analoge Darstellung der Anisotropie erfordert, dass ein komplettes System von [[Skalar (Physik)|skalarwertigen]] Funktionen bekannt ist, die unter allen Transformationen in der Symmetriegruppe ''g<sub>R</sub>'' invariant sind.<ref name="haupt"/>{{rp|380}} Bei der Orthotropie bleiben die folgenden Terme invariant:<ref name="haupt"/>{{rp|382}} | |||
: E<sub>11</sub>, E<sub>22</sub>, E<sub>33</sub>, E<sub>23</sub><sup>2</sup>, E<sub>13</sub><sup>2</sup>, E<sub>12</sub><sup>2</sup>, E<sub>12</sub>E<sub>23</sub>E<sub>13</sub>. | |||
Darin ist E<sub>ij</sub> := ê<sub>i</sub>·'''E'''·ê<sub>j</sub> für i,j=1,2,3 und ê<sub>1,2,3</sub> sind die [[Einheitsvektor]]en in Richtung der paarweise orthogonalen Orthotropieachsen. | |||
== Orthotropie in der Linearen Elastizitätstheorie == | == Orthotropie in der Linearen Elastizitätstheorie == | ||
Gegeben sind zwei [[Tensor]]en zweiter Stufe <math>\boldsymbol{\sigma}</math> und <math>\boldsymbol{\varepsilon}</math> mit 3×3-Koeffizienten <math>\sigma_{ij}</math> bzw. <math>\varepsilon_{ij}</math>. Der allgemeinste lineare Zusammenhang, den es zwischen diesen Koeffizienten gibt, ist: | Gegeben sind zwei [[Tensor]]en zweiter Stufe <math>\boldsymbol{\sigma}</math> und <math>\boldsymbol{\varepsilon}</math> mit 3×3-Koeffizienten <math>\sigma_{ij}</math> bzw. <math>\varepsilon_{ij}</math>. Der allgemeinste lineare Zusammenhang, den es zwischen diesen Koeffizienten gibt, ist: | ||
:<math> | :<math> | ||
Zeile 59: | Zeile 91: | ||
Die Matrix <math>C^\text{v}</math> mit den 36 unabhängigen Komponenten <math>C_{ijkl}</math> repräsentiert den [[Elastizitätstensor]] des Materials. Im Fall der [[Hyperelastizität]] ist diese Matrix [[Symmetrische Matrix|symmetrisch]], so dass dann nur noch 21 Einträge unabhängig sind. | Die Matrix <math>C^\text{v}</math> mit den 36 unabhängigen Komponenten <math>C_{ijkl}</math> repräsentiert den [[Elastizitätstensor]] des Materials. Im Fall der [[Hyperelastizität]] ist diese Matrix [[Symmetrische Matrix|symmetrisch]], so dass dann nur noch 21 Einträge unabhängig sind. | ||
=== Materialparameter === | === Materialparameter === | ||
Die | Die Koeffizienten <math>C_{ijkl}</math> des Elastizitätstensors ergeben sich bei orthotroper linearer Elastizität aus nur neun Materialparametern, die in Versuchen an [[Unidirektionale Schicht|makroskopischen Proben]] ermittelt werden können: | ||
{| class="wikitable" | {| class="wikitable" | ||
Zeile 109: | Zeile 106: | ||
|} | |} | ||
Die | Die [[Dimension (Größensystem)|Dimension]] der Elastizitätsmoduln <math>E_1,E_2,E_3</math> und Schubmoduln <math>G_{12}, G_{23}, G_{13}</math> ist Kraft pro Fläche während die Querkontraktionszahlen <math>\nu_{ij}</math> 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 <math>\nu_{12}</math>. Die Normaldehnung in ''i''-Richtung wird mit <math>\varepsilon_{ii}</math> bezeichnet. Dann ist für beliebige Werkstoffe die Querkontraktionszahl <math>\nu_{ij}</math> das negative Verhältnis der Normaldehnung in ''j''-Richtung (Wirkung) zu derjenigen in ''i''-Richtung bei Zug in ''i''-Richtung (Ursache): | |||
\nu_{ | |||
\nu_{ | :<math>\nu_{ij} = -\frac{\varepsilon_{jj}}{\varepsilon_{ii}} </math> | ||
</math> | |||
Wegen des Ursache-Wirkungs-Konzepts ist meistens <math>\nu_{ij}\nu_{ji}\ne1</math>. | |||
=== 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: | |||
:<math> | :<math> | ||
Zeile 127: | Zeile 126: | ||
2\varepsilon_{12} | 2\varepsilon_{12} | ||
\end{bmatrix} = | \end{bmatrix} = | ||
\begin{bmatrix} | \underbrace{\begin{bmatrix} | ||
\frac{1}{E_1} & -\frac{\nu_{ | \frac{1}{E_1} & -\frac{\nu_{21}}{E_2} & -\frac{\nu_{31}}{E_3} & & & \\ | ||
& \frac{1}{E_2} & -\frac{\nu_{ | -\frac{\nu_{12}}{E_1} & \frac{1}{E_2} & -\frac{\nu_{32}}{E_3} & & & \\ | ||
& & \frac{1}{E_3} & | -\frac{\nu_{13}}{E_1} & -\frac{\nu_{23}}{E_2} & \frac{1}{E_3} & & & \\ | ||
& & & \frac{1}{G_{23}} & | & & & \frac{1}{G_{23}} & & \\ | ||
& | & & & & \frac{1}{G_{13}} & \\ | ||
& & & & & \frac{1}{G_{12}} | & & & & & \frac{1}{G_{12}} | ||
\end{bmatrix} | \end{bmatrix}}_{=:S} | ||
\begin{bmatrix} | \begin{bmatrix} | ||
\sigma_{11}\\ | \sigma_{11}\\ | ||
Zeile 145: | Zeile 144: | ||
</math> | </math> | ||
Die Matrix S ist die ''Nachgiebigkeitsmatrix'' und ihre Symmetrie erfordert: | |||
:<math> | |||
\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}\,. | |||
</math> | |||
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'' | |||
:<math> | :<math> | ||
Zeile 159: | Zeile 168: | ||
\begin{bmatrix} | \begin{bmatrix} | ||
\frac{1-\nu_{23}\nu_{32}}{D}E_{1}& | \frac{1-\nu_{23}\nu_{32}}{D}E_{1}& | ||
\frac{\nu_{ | \frac{\nu_{21}+\nu_{23}\nu_{31}}{D}E_{1}& | ||
\frac{\nu_{ | \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{1-\nu_{13}\nu_{31}}{D}E_{2}& | ||
\frac{\nu_{ | \frac{\nu_{32}+\nu_{31}\nu_{12}}{D}E_{2} | ||
\\ | \\ | ||
& & | \frac{\nu_{13}+\nu_{12}\nu_{23}}{D}E_{3} & | ||
\frac{1-\nu_{12}\nu_{21}}{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_{23} | ||
\\ | \\ | ||
& | & & & &G_{13} | ||
\\ | \\ | ||
& & & & &G_{12} | & & & & &G_{12} | ||
Zeile 191: | Zeile 201: | ||
-2\nu_{12}\nu_{23}\nu_{31}\,.</math> | -2\nu_{12}\nu_{23}\nu_{31}\,.</math> | ||
Die Nachgiebigkeitsmatrix und Steifigkeitsmatrix sind | Die Nachgiebigkeitsmatrix und die Steifigkeitsmatrix sind an denselben Stellen mit von Null verschiedenen Werten besetzt. | ||
[[ | === Spezialfälle der Orthotropie === | ||
In der [[Kubische Anisotropie|kubischen Anisotropie]] sind die Elastizitäts- und Schubmoduln sowie die Querdehnzahlen alle gleich: | |||
:<math>\ | :<math>E_1=E_2=E_3,\quad G_{12}=G_{13}=G_{23},\quad | ||
E_2 &= | \nu_{12}=\nu_{21}=\nu_{13}=\nu_{31}=\nu_{23}=\nu_{32}</math> | ||
womit nur drei unabhängige Elastizitätsparameter übrig bleiben. [[Transversale Isotropie]] mit fünf unabhängigen Elastizitätsparametern stellt sich ein mit: | |||
:<math>\begin{align} | |||
E_2 &= E_3 | |||
\\ | \\ | ||
G_{12} &= | G_{12} &= G_{13} | ||
\\ | \\ | ||
\nu_{12} &= | \nu_{12} &= \nu_{13} | ||
\\ | \\ | ||
G_{23} &= | G_{23} &= \frac{E_2 }{2 (1+\nu_{23})} | ||
\end{ | \end{align}</math> | ||
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 σ<sub>13</sub>=σ<sub>23</sub>=σ<sub>33</sub>=0 und aus letzterer Identität leitet sich | |||
:<math>\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} ) | |||
</math> | </math> | ||
ab. Das Elastizitätsgesetz vereinfacht sich zu | |||
:<math> | :<math> | ||
Zeile 213: | Zeile 238: | ||
\varepsilon_{11}\\ | \varepsilon_{11}\\ | ||
\varepsilon_{22}\\ | \varepsilon_{22}\\ | ||
2\varepsilon_{12} | |||
\varepsilon_{12} | |||
\end{bmatrix} = | \end{bmatrix} = | ||
\begin{bmatrix} | \begin{bmatrix} | ||
\frac{1}{E_1} & -\frac{\nu_{ | \frac{1}{E_1} & -\frac{\nu_{21}}{E_2} & 0 \\ | ||
-\frac{\nu_{ | -\frac{\nu_{12}}{E_1} & \frac{1}{E_2} & 0 \\ | ||
0 & 0 & \frac{1}{G_{12}} | |||
0 & 0 & \frac{1}{ | |||
\end{bmatrix} | \end{bmatrix} | ||
\begin{bmatrix} | \begin{bmatrix} | ||
Zeile 229: | Zeile 252: | ||
</math> | </math> | ||
bzw. | |||
:<math> | :<math> | ||
\frac{\nu_{12}}{ | \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} | |||
</math> | </math> | ||
mit jeweils symmetrischer Nachgiebigkeits- bzw. Steifigkeitsmatrix. | |||
In der linearen orthotropen Elastizität für den Ebenen Spannungszustand werden die Schubmoduln <math>G_{13}, G_{23}</math> 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 === | === Ebener Verzerrungszustand === | ||
Hier | Hier finden die Verzerrungen ausschließlich in der 1-2-Ebene statt, nur die Normalspannung senkrecht zur Ebene darf auftreten. Aus <math>\varepsilon_{33}=0</math> leitet sich | ||
:<math>\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}] | |||
</math> | |||
mit | |||
:<math> | |||
D=1-\nu_{12}\nu_{21}-\nu_{13}\nu_{31}-\nu_{23}\nu_{32} | |||
-2\nu_{12}\nu_{23}\nu_{31}\,.</math> | |||
ab. Das Elastizitätsgesetz reduziert sich auf | |||
:<math> | :<math> | ||
\begin{bmatrix} | \begin{bmatrix} | ||
\sigma_{11}\\ | |||
\sigma_{22}\\ | |||
\sigma_{12}\\ | |||
\end{bmatrix} | \end{bmatrix} | ||
= | = | ||
\begin{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&0& G_{12} | |||
\end{bmatrix} | \end{bmatrix} | ||
\begin{bmatrix} | \begin{bmatrix} | ||
\varepsilon_{11}\\ | |||
\varepsilon_{22}\\ | |||
2\varepsilon_{12} | |||
\end{bmatrix} | \end{bmatrix} | ||
</math> | </math> | ||
bzw. | |||
:<math> | :<math> | ||
\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} | |||
</math> | |||
mit jeweils symmetrischer Steifigkeits- bzw. Nachgiebigkeitsmatrix. | |||
In der linearen orthotropen Elastizität für den Ebenen Verzerrungszustand werden die Schubmoduln <math>G_{13}, G_{23}</math> nicht gebraucht, sodass nur sieben Materialparameter ausreichen, und wenn nur die Spannungen und Verzerrungen in der Ebene interessieren entfällt zusätzlich ''E''<sub>3</sub>, sodass nur mehr sechs Materialparameter gebraucht werden. | |||
=== Stabilitätskriterien === | === Stabilitätskriterien === | ||
Zeile 279: | Zeile 348: | ||
* die Determinante der Steifigkeits- und Nachgiebigkeitsmatrix muss positiv sein (damit es unter Druck komprimiert und nicht expandiert). | * 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:<ref> | Werden an einem realen Werkstoff Materialparameter identifiziert, die diesen Stabilitätskriterien widersprechen, ist Vorsicht geboten. Die Stabilitätskriterien lauten:<ref name="altenbach">{{Literatur | ||
|Autor=[[Holm Altenbach]] | |||
|Titel=Kontinuumsmechanik | |||
|TitelErg=Einführung in die materialunabhängigen und materialabhängigen Gleichungen | |||
|Verlag=Springer-Verlag | |||
|Ort=Berlin/Heidelberg | |||
|Datum=2012 | |||
|ISBN=978-3-642-24118-5 | |||
|DOI=10.1007/978-3-642-24119-2|Seiten=331}}</ref> | |||
:<math>\begin{array}{l} | :<math>\begin{array}{l} | ||
Zeile 305: | Zeile 382: | ||
\end{array} | \end{array} | ||
</math> | </math> | ||
=== 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 | |||
:<math>\begin{align} | |||
w(\boldsymbol{\varepsilon})=& | |||
\frac a2\varepsilon_{11}^2 + \frac b2\varepsilon_{22}^2 | |||
+ \frac c2\varepsilon_{33}^2 | |||
+ d\varepsilon_{11}\varepsilon_{22} + e\varepsilon_{11}\varepsilon_{33} | |||
+ f\varepsilon_{22}\varepsilon_{33} | |||
\\& | |||
+ 2 g\varepsilon_{23}^2 + 2 h\varepsilon_{13}^2 + 2 p\varepsilon_{12}^2 | |||
\end{align}</math> | |||
mit neun Parametern ''a'' bis ''p''. Um dies nach '''''ε''''' ableiten zu können, müssen die Komponenten ''ε''<sub>ij</sub> als Funktion des Tensors '''''ε''''' ausgedrückt werden. Dies gelingt mit der Darstellung des [[Frobenius-Skalarprodukt#Darstellung als Spur|Frobenius-Skalarprodukts]] ":" als [[Spur (Mathematik)|Spur]]: | |||
:<math>\mathbf{A}:\mathbf{B}:=\mathrm{Spur}(\mathbf{A^\top\cdot B})</math> | |||
Mit der Abkürzung <math>\mathbf{K}_{ij}=\frac12(\hat e_i\otimes\hat e_j+\hat e_j\otimes\hat e_i)</math> für die symmetrisierten [[Dyadisches Produkt|dyadischen Produkte]] der Orthotropieachsenvektoren ist dann<ref name="Frechet">Die ij-Komponente eines beliebigen Tensors zweiter Stufe '''T''' im ê<sub>1,2,3</sub>-System ist | |||
:<math>\begin{align}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{align}</math> | |||
Die [[Fréchet-Ableitung]] hiervon nach '''T''' ist der beschränkte lineare Operator <math>\mathcal{A}</math> der – sofern er existiert – in allen Richtungen '''H''' dem [[Gâteaux-Differential]] entspricht, also | |||
:<math>\begin{align}\mathcal{A}(\mathbf{H}) | |||
=&\left.\frac{\mathrm d}{\mathrm ds} | |||
[(\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{align}</math> | |||
Darin ist <math>s\in\R</math> und der lineare Operator ist das Skalarprodukt mit <math>\hat e_i\otimes\hat e_j</math>. Hier ist <math>\mathbf{T}=\boldsymbol\varepsilon</math> ein symmetrischer Tensor, dessen Differential '''H''' auch ein symmetrischer Tensor ist. Beim Skalarprodukt mit diesem trägt nur der symmetrische Anteil etwas bei: | |||
:<math>(\hat e_i\otimes\hat e_j):\mathbf{H} | |||
=\frac12(\hat e_i\otimes\hat e_j+\hat e_j\otimes\hat e_i):\mathbf H | |||
=\mathbf{K}_{ij}:\mathbf H</math> | |||
Dann wird auch | |||
:<math>\mathcal{A} | |||
=\frac{\partial T_{ij}}{\partial\mathbf{T}} | |||
=\mathbf{K}_{ij} | |||
</math> | |||
geschrieben. | |||
</ref> | |||
:<math>\mathbf{K}_{ij}:\boldsymbol{\varepsilon} | |||
=\frac12(\varepsilon_{ij}+\varepsilon_{ji})=\varepsilon_{ij} | |||
\quad\rightarrow\quad | |||
\frac{\mathrm d\varepsilon_{ij}}{\mathrm d\boldsymbol{\varepsilon}} | |||
=\mathbf{K}_{ij} | |||
</math> | |||
Aus dem Ansatz der Formänderungsenergie berechnen sich die Spannungen zu | |||
:<math>\begin{align} | |||
\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} ) | |||
\\& | |||
+ 4 g \varepsilon_{23} \mathbf{K}_{23} | |||
+ 4 h \varepsilon_{13} \mathbf{K}_{13} | |||
+ 4 p \varepsilon_{12} \mathbf{K}_{12} | |||
\end{align}</math> | |||
oder in Voigt-Notation im ê<sub>1,2,3</sub>-System | |||
:<math>\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}\\ | |||
2 g\varepsilon_{23}\\ | |||
2 h\varepsilon_{13}\\ | |||
2 p\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} | |||
</math> | |||
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: | |||
:<math>\begin{align} | |||
\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} ) | |||
\\& | |||
+ 4 g \mathbf{K}_{23}\otimes\mathbf{K}_{23} | |||
+ 4 h \mathbf{K}_{13}\otimes\mathbf{K}_{13} | |||
+ 4 p \mathbf{K}_{12}\otimes\mathbf{K}_{12} | |||
\end{align}</math> | |||
Die Tensoren '''K'''<sub>ii</sub> werden ''Strukturvariable'' genannt, weil sie die interne Struktur des Materials repräsentieren<ref name="haupt"/>{{rp|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.<ref name="haupt"/>{{rp|394}} | |||
== Gründe für die Besetztheit der Steifigkeitsmatrix == | == Gründe für die Besetztheit der Steifigkeitsmatrix == | ||
Zeile 310: | Zeile 517: | ||
=== Drehmatrizen bei 180-Grad-Drehungen === | === Drehmatrizen bei 180-Grad-Drehungen === | ||
Die (linearen) Abbildungen, die 180-Grad-Drehungen um die | 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 [[Orthogonale Matrix|orthogonalen Matrizen]] folgende Gestalt | ||
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 | |||
[[ | |||
::<math> | ::<math> | ||
Zeile 324: | Zeile 527: | ||
0 & 0 & -1 \\ | 0 & 0 & -1 \\ | ||
\end{bmatrix}, | \end{bmatrix}, | ||
A_y | &&A_y = | ||
\begin{bmatrix} | \begin{bmatrix} | ||
-1 & 0 & 0\\ | -1 & 0 & 0\\ | ||
Zeile 330: | Zeile 533: | ||
0 & 0 & -1 | 0 & 0 & -1 | ||
\end{bmatrix}, | \end{bmatrix}, | ||
A_z | &&A_z = | ||
\begin{bmatrix} | \begin{bmatrix} | ||
-1 & 0 & 0\\ | -1 & 0 & 0\\ | ||
Zeile 339: | Zeile 542: | ||
</math> | </math> | ||
Diese 3 Matrizen (und zusätzlich die [[ | Diese 3 Matrizen bilden eine echte Untergruppe der [[Drehgruppe]] SO(3). Das Produkt dieser drei Matrizen ist die [[Einheitsmatrix]]: <math>A_xA_yA_z=E</math>. | ||
Untergruppe | |||
Die 3 Matrizen ''A''<sub>x,y,z</sub> und zusätzlich die negative Einheitsmatrix -''E'', die eine [[Punktspiegelung]] repräsentiert, bilden die ''Symmetriegruppe'' des orthotropen Materials<ref name="haupt"/>{{rp|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 === | === Symmetriebedingung in Indexschreibweise und Voigt’scher Notation === | ||
Zeile 517: | Zeile 721: | ||
=== Spezialfall 180-Grad-Drehungen === | === Spezialfall 180-Grad-Drehungen === | ||
Da im Spezialfall der Orthotropie die 3×3- | Da im Spezialfall der Orthotropie die 3×3-Matrizen <math>A</math> nur auf der [[Hauptdiagonale]]n besetzt sind, vereinfachen sich die Definitionen von oben zu: | ||
der [[Hauptdiagonale]]n besetzt sind, vereinfachen sich die Definitionen | <math> | ||
von oben zu | |||
\begin{align} | \begin{align} | ||
A^{\text{v}}_\sigma= | A^{\text{v}}_\sigma= | ||
Zeile 554: | Zeile 756: | ||
\end{bmatrix},\qquad | \end{bmatrix},\qquad | ||
A^{\text{v}}_y | A^{\text{v}}_y | ||
= | |||
\begin{bmatrix} | \begin{bmatrix} | ||
1 & 0 & 0 & 0 & 0 & 0 \\ | 1 & 0 & 0 & 0 & 0 & 0 \\ | ||
Zeile 564: | Zeile 766: | ||
\end{bmatrix},\qquad | \end{bmatrix},\qquad | ||
A^{\text{v}}_z | A^{\text{v}}_z | ||
= | |||
\begin{bmatrix} | \begin{bmatrix} | ||
1 & 0 & 0 & 0 & 0 & 0 \\ | 1 & 0 & 0 & 0 & 0 & 0 \\ | ||
Zeile 676: | Zeile 878: | ||
== Siehe auch == | == Siehe auch == | ||
* [[Materialmodell]] | |||
*[[Materialmodell | * [[Spezielle lineare Gruppe]] | ||
*[[Spezielle lineare Gruppe]] | |||
== Weblinks == | == Weblinks == | ||
{{ | {{Wiktionary|orthotrop}} | ||
== Einzelnachweise == | == Einzelnachweise == | ||
<references /> | <references /> | ||
== Literatur == | |||
* {{Literatur | |||
|Autor=J.F. Nye | |||
|Titel=Physical Properties of Crystals: Their Representation by Tensors and Matrices | |||
|Verlag=[[Oxford University Press]] | |||
|Datum=1985 | |||
|ISBN=978-0-19-851165-6}} | |||
* {{Literatur | |||
|Autor=[[Holm Altenbach|H. Altenbach]], J. Altenbach, R. Rikards | |||
|Titel=Einführung in die Mechanik der Laminat- und Sandwichtragwerke | |||
|Verlag=Deutscher Verlag für Grundstoffindustrie | |||
|Ort=Stuttgart | |||
|Datum=1996 | |||
|ISBN=3-342-00681-1}} | |||
[[Kategorie:Werkstoffeigenschaft]] | [[Kategorie:Werkstoffeigenschaft]] | ||
[[Kategorie:Kontinuumsmechanik]] | [[Kategorie:Kontinuumsmechanik]] |
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:
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.
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
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
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
In der isotropen Hyperelastizität hängt die Formänderungsenergie von den Hauptinvarianten I1,2,3 des Verzerrungstensors E ab:
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
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.
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:
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:
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.
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):
Wegen des Ursache-Wirkungs-Konzepts ist meistens $ \nu _{ij}\nu _{ji}\neq 1 $.
Ein Material ist linear elastisch orthotrop, wenn eine Orthonormalbasis existiert, so dass das Elastizitätsgesetz dargestellt in Bezug auf diese Basis folgende Form annimmt:
Die Matrix S ist die Nachgiebigkeitsmatrix und ihre Symmetrie erfordert:
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
mit
Die Nachgiebigkeitsmatrix und die Steifigkeitsmatrix sind an denselben Stellen mit von Null verschiedenen Werten besetzt.
In der kubischen Anisotropie sind die Elastizitäts- und Schubmoduln sowie die Querdehnzahlen alle gleich:
womit nur drei unabhängige Elastizitätsparameter übrig bleiben. Transversale Isotropie mit fünf unabhängigen Elastizitätsparametern stellt sich ein mit:
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.
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
ab. Das Elastizitätsgesetz vereinfacht sich zu
bzw.
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.
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
mit
ab. Das Elastizitätsgesetz reduziert sich auf
bzw.
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.
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:
Werden an einem realen Werkstoff Materialparameter identifiziert, die diesen Stabilitätskriterien widersprechen, ist Vorsicht geboten. Die Stabilitätskriterien lauten:[2]
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:
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
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:
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]
Aus dem Ansatz der Formänderungsenergie berechnen sich die Spannungen zu
oder in Voigt-Notation im ê1,2,3-System
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:
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
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.
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
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.
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
Mithilfe eines linearen Materialgesetzes $ f_{C} $ lässt sich für gegebenen Verzerrungstensor der zugehörige Spannungstensor ermitteln. Es sei
Im allgemeinen Fall der Anisotropie gilt zwar nicht
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
In Indexschreibweise
Nun dieselbe Bedingung in Voigt’scher Notation: Mit der Definition
gilt
Mit der neuen Definition
ergibt sich
In Voigt’scher Notation erhält man also als Symmetriebedingung
Und da dies für beliebige Dehnungen gelten muss, ist die Symmetriebedingung
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
Die Symmetriebedingung ausgewertet für diese Matrizen ergibt
An den letzten 3 Gleichungen erkennt man, dass $ C $ nur folgende Gestalt haben kann
Da diese Voigt’sche Steifigkeitsmatrix außerdem symmetrisch ist, bleibt