Transversale Isotropie

Transversale Isotropie

Bildhafte Erklärung der Transversalen Isotropie.
Der Werkstoff ist rotationssymmetrisch bezüglich der 1-Achse, die senkrecht auf der isotropen 2-3-Ebene steht.
Ein so orientierter Rundstab aus diesem Material kann um seine Längsachse gedreht werden, ohne dass sich seine Eigenschaften ändern.

Die transversale Isotropie (von lateinisch transversus „quer“ sowie altgr. {{Modul:Vorlage:lang}} Modul:Multilingual:149: attempt to index field 'data' (a nil value) isos „gleich“ und {{Modul:Vorlage:lang}} Modul:Multilingual:149: attempt to index field 'data' (a nil value) tropos „Drehung, Richtung“) ist eine spezielle Art der Richtungsabhängigkeit eines Werkstoffs. Transversal isotrope Materialien haben die drei Eigenschaften:

  1. Es gibt eine Vorzugsrichtung, die 1-Richtung im Bild, in der das Kraft-Verformungs-Verhalten des Materials anders ist als senkrecht dazu.
  2. Senkrecht zur Vorzugsrichtung, in 2- und 3-Richtung, sind die Materialeigenschaften unabhängig von der Richtung (isotrope Ebene) und
  3. in einem Bezugssystem parallel zur Vorzugsrichtung gibt es keine Kopplung zwischen Normaldehnungen und Schubverzerrungen.

In Ebenen, die nicht senkrecht zur Vorzugsrichtung sind, ist das Kraft-Verformungs-Verhalten des Materials richtungsabhängig.

Den speziellen Fall, dass ein Material (an einem Teilchen) unabhängig von der Belastungsrichtung jeweils dasselbe Kraft-Verformungs-Verhalten zeigt, bezeichnet man als Isotropie. Den allgemeinen Fall, dass das Kraft-Verformungs-Verhalten von der Belastungsrichtung abhängt, bezeichnet man dagegen als Anisotropie. Die transversale Isotropie ist ein Sonderfall der Orthotropie und Anisotropie und enthält ihrerseits die Isotropie als Spezialfall.

Ein linear elastisches transversal isotropes Material besitzt maximal fünf Materialparameter.

Unidirektional verstärkte Kunststoffe sind im ungeschädigten Zustand in guter Näherung transversal isotrop.

Bedeutung in der Konstruktion

Unidirektional verstärkte Kunststoffe sind im ungeschädigten Zustand in guter Näherung transversal isotrop. Sie haben eine hohe Festigkeit in Richtung der Fasern und sind senkrecht dazu nachgiebiger. In der Konstruktion werden transversal isotrope Werkstoffe gerne eingesetzt, denn sie gestatten die Werkstoffeigenschaften an die Belastung anzupassen. Unter anderem die geringe Dichte bei hoher Festigkeit in Belastungsrichtung haben zu einer starken Zunahme der Nutzung der faserverstärkten Kunststoffe geführt. Durch Schädigung verlieren diese Werkstoffe im Allgemeinen ihre transversale Isotropie.

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 transversalen Isotropie sind dies beliebige Drehungen um die Vorzugsrichtung oder 180-Grad-Drehungen senkrecht zur Vorzugsrichtung. Diese Drehungen bilden die Symmetriegruppe des transversal isotropen Materials[1].

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 beliebig parallel zur Vorzugsrichtung oder um 180 Grad senkrecht dazu. Dann bringt man dieselbe Kraft auf wie im ersten Experiment und misst erneut die Verformung. Bei transversal isotropem 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 senkrecht zur Vorzugsrichtung dreht. Wenn nicht der Spezialfall der 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 w invariant lassen. Mathematisch wird das mit dem Verzerrungstensor E durch

QgRw(QEQ)=w(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 transversal isotropen Materials ist[1]:382

gR={1,Q1φ,Q2π}

Darin steht Qiα für den orthogonalen Tensor, der mit dem Winkel α in Radiant um die i-te Orthotropieachse dreht; der Winkel φ ist beliebig. Die Orthotropieachsen sind zum einen die Vorzugsrichtung (Faserrichtung) ê1 und ê2,3 sind dazu senkrecht in der isotropen Ebene und zueinander senkrecht.

Invarianten

In der isotropen Hyperelastizität hängt die Formänderungsenergie von den Hauptinvarianten I1,2,3 des Verzerrungstensors E ab:

w(E)=w(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 transversalen Isotropie bleiben die folgenden Terme invariant:[1]:382

E22+E33, E11, E22E33-E232, E122+E132, det(E)

Darin ist Eij := êi·E·êj für i,j=1,2,3, wobei ê1,2,3 die Einheitsvektoren in Richtung der paarweise orthogonalen Orthotropieachsen sind und det bildet die Determinante.

Transversal isotrope Elastizität

Ein transversal isotroper linear elastischer Werkstoff zeichnet sich dadurch aus, dass in seiner Steifigkeits- oder Nachgiebigkeitsmatrix die Koppelterme nicht besetzt sind. Schubspannungen in Ebenen parallel oder senkrecht zur Vorzugsrichtung führen nicht zu Normaldehnungen. In einem solchen Material existiert eine Orthonormalbasis ê1,2,3 in der die Spannungs-Dehnungs-Beziehung

[ε11ε22ε332ε232ε312ε12]=[1E1ν21E2ν31E3ν12E11E2ν32E3ν13E1ν23E21E31G231G131G12]=:S[σ11σ22σ33σ23σ31σ12]

mit der gezeigten Nachgiebigkeitsmatrix S zwischen den Spannungen σij und den Dehnungen εij vorliegt. Die Dimensionen der Elastizitätsmoduln E1,E2,E3 und Schubmoduln G12,G23,G13 sind Kraft pro Fläche während die Querkontraktionszahlen νij dimensionslos sind. Die Indizes der Querkontraktionszahlen sind sorgfältig definiert durch das negative Verhältnis der Normaldehnung εjj in j-Richtung (Wirkung) zu derjenigen εii in i-Richtung bei Zug in i-Richtung (Ursache):

νij=εjjεii

Wegen des Ursache-Wirkungs-Konzepts ist meistens νijνji1.

Materialparameter

Die zwölf in der obigen Nachgiebigkeitsmatrix vorkommenden Kennwerte ergeben sich bei transversal isotroper, linearer Elastizität aus nur fünf Materialparametern, die in Versuchen an makroskopischen Proben ermittelt werden können:

Formelzeichen Bedeutung
E Elastizitätsmodul in Vorzugsrichtung
E Elastizitätsmodul senkrecht zur Vorzugsrichtung
ν Querkontraktionszahl bei Zug in Vorzugsrichtung
G Schubmodul in Ebenen parallel zur Vorzugsrichtung
G Schubmodul in der isotropen Ebene

Aufgrund der transversalen Isotropie sind die folgenden Ausdrücke im 1-2-3-System identisch:[2]

E1=EE2=E3=EG12=G13=GG23=Gν12=ν13=νν21=ν31ν23=ν32

Aus thermodynamischen Gründen (vergleiche Cauchy-Elastizität und Hyperelastizität) ist die Nachgiebigkeitsmatrix symmetrisch und legt so

ν21E2=ν12E1ν21=ν31=EEν

fest. Die Querkontraktionszahl in der Ebene senkrecht zur Vorzugsrichtung ist schließlich durch die Isotropieannahme gebunden:

G=E2(1+ν23)ν23=ν32=E2G1.

So sind alle zwölf Kennwerte auf die fünf Materialparameter zurückgeführt. Isotropie stellt sich mit

E1=E2=EG12=G23=G=E2(1+ν)

als Spezialfall ein.

Spannungs-Dehnungs-Beziehung

Damit lautet das Elastizitätsgesetz bei transversal isotroper, linearer Elastizität:

[ε11ε22ε332ε232ε312ε12]=[1E1ν12E1ν12E10001E2ν23E20001E20001G2300sym1G1201G12][σ11σ22σ33σ23σ31σ12]

zwischen den Spannungen σij und den Dehnungen εij. Durch Invertierung der Nachgiebigkeitsmatrix erhält man die Steifigkeitsmatrix:

[σ11σ22σ33σ23σ13σ12]=[C11112ν12(λ+G23)2ν12(λ+G23)000λ+2G23λ000λ+2G23000G2300symG120G12][ε11ε22ε332ε232ε132ε12]

mit

λ=ν12ν21+ν23(1ν232ν12ν21)(1+ν23)E2C1111=1ν231ν232ν12ν21E1.

Diese für kleine Dehnungen in voigtscher Notation geschriebene lineare Matrizengleichung zwischen Spannungen und Dehnungen lässt sich mit Hyperelastizität auf nichtlinear elastisches transversal isotropes Verhalten verallgemeinern.

Ebener Spannungszustand

In dünnwandigen Strukturen aus transversal isotropem Material ist die Vorzugsrichtung oftmals in den Vorzugsrichtungen der Struktur gelegen, wie zum Beispiel bei der unidirektionalen Schicht aus der Faser-Kunststoff-Verbunde bestehen, und es liegt ein ebener Spannungszustand vor.

Hier ist σ132333=0 und aus letzterer Identität leitet sich

ε33=ν12(1+ν23)ε11+(ν23+ν21ν12)ε221ν12ν21=1E3(ν31σ11+ν32σ22)

ab. Das Elastizitätsgesetz vereinfacht sich zu

[ε11ε222ε12]=[1E1ν21E20ν12E11E20001G12][σ11σ22σ12]

bzw.

[σ11σ22σ12]=[E11ν12ν21ν21E11ν12ν210ν12E21ν12ν21E21ν12ν21000G12][ε11ε222ε12]

mit jeweils symmetrischer Nachgiebigkeits- bzw. Steifigkeitsmatrix. Die Elastizitätsgesetze sind dieselben wie bei der Orthotropie.

In der linearen transversal isotropen Elastizität für den Ebenen Spannungszustand senkrecht zur isotropen Ebene werden alle Materialparameter gebraucht; nur wenn ausschließlich die Spannungen und Verzerrungen in der Ebene interessieren entfällt der Schubmodul G23, sodass nur vier Materialparameter ausreichen.

Ebener Verzerrungszustand parallel zur Vorzugsrichtung

Beim ebenen Verzerrungszustand parallel zur Vorzugsrichtung finden die Spannungen und Verzerrungen ausschließlich in der 1-2-Ebene statt, nur die Normalspannung senkrecht zur Ebene darf auftreten. Das ist näherungsweise bei einer dicken unidirektionalen Schicht, die flächig belastet wird, der Fall. Aus ε33=0 leitet sich

σ33=ν31σ11+ν32σ22=E3D[ν12(1+ν23)ε11+(ν23+ν12ν21)ε22]

mit D=(1+ν23)(1ν232ν12ν21) ab. Das Elastizitätsgesetz reduziert sich auf

[σ11σ22σ12]=[1ν232DE1ν21(1+ν23)DE10ν12(1+ν23)DE21ν21ν12DE2000G12][ε11ε222ε12]

bzw.

[ε11ε222ε12]=[1ν12ν21E1(1+ν23)ν21E20(1+ν23)ν12E11ν232E20001G12][σ11σ22σ12]

mit jeweils symmetrischer Steifigkeits- bzw. Nachgiebigkeitsmatrix.

In der linearen transversal isotropen Elastizität für den Ebenen Verzerrungszustand senkrecht zur isotropen Ebene werden alle fünf Materialparameter gebraucht.

Ebener Verzerrungszustand in der isotropen Ebene

Bei einem prismatischen Körper in Vorzugsrichtung, der nur geringfügig gestaucht oder gestreckt wird, liegt in guter Näherung ein ebener Verzerrungszustand in der isotropen Ebene vor. Dann ergibt sich

σ11=ν12(σ22+σ33)=ν12E2(ε22+ε33)1ν232ν12ν21

und das Elastizitätsgesetz

[σ22σ33σ23]=[(1ν12ν21)E2D(ν23+ν12ν21)E2D0(ν23+ν12ν21)E2D(1ν12ν21)E2D000G12][ε22ε332ε23]

mit D=(1+ν23)(1ν232ν12ν21) bzw.

[ε22ε332ε23]=[1ν12ν21E2ν23+ν12ν21E20ν23+ν12ν21E21ν12ν21E20001G12][σ22σ33σ23]

mit jeweils symmetrischer Steifigkeits- bzw. Nachgiebigkeitsmatrix. Hier wird der Elastizitätsmodul in Vorzugsrichtung E1 nicht gebraucht, weswegen nur vier Materialparameter zu bestimmen sind.

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:[3]

E1,E2,G12,G23>0|ν23|<1|ν12|<E1E21ν12ν21>01ν232ν12ν21>0

Wenn die linke Seite der letzten Ungleichung gegen null geht, setzt das Material einer hydrostatischen Kompression zunehmend Widerstand entgegen. Aus der Symmetrie-Beziehung folgt ergänzend:

|ν21|<E2E1.

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. Mit den #Invarianten lässt sich der Ansatz

w(ε)=a2ε112+b2(ε22+ε33)2+cε11(ε22+ε33)+(db)(ε22ε33ε232)+2e(ε122+ε132)

mit fünf Parametern a bis e machen. Um dies nach ε ableiten zu können, müssen die Komponenten εij als Funktion des Tensors ε ausgedrückt werden. Das gelingt mit der Darstellung des Frobenius-Skalarprodukts ":" als Spur:

A:B:=Spur(AB)

Mit der Abkürzung Kij=12(e^ie^j+e^je^i) für die symmetrisierten dyadischen Produkte der Orthotropieachsenvektoren ist dann[4]

Kij:ε=12(εij+εji)=εijdεijdε=Kij

Aus dem Ansatz der Formänderungsenergie berechnen sich die Spannungen zu

σ=dwdε=aε11K11+b(ε22+ε33)(K22+K33)+c(ε11(K22+K33)+(ε22+ε33)K11)+(db)(ε33K22+ε22K332ε23K23)+4e(ε12K12+ε13K13)

oder in Voigt-Notation im ê1,2,3-System

[σ11σ22σ33σ23σ13σ12]=[aε11+cε22+cε33cε11+bε22+dε33cε11+dε22+bε33(bd)ε232eε132eε12]=[acccbdcdbbd2ee][ε11ε22ε332ε232ε132ε12]

In der #Spannungs-Dehnungs-Beziehung lassen sich die Parameter direkt ablesen, womit sich die 44-Komponente ergibt zu

bd2=E22(1+ν23)=G23

Ableitung der Spannungen nach den Dehnungen liefert den konstanten und symmetrischen Elastizitätstensor 4. Stufe:

C:=dσdε=aK11K11+b(K22+K33)(K22+K33)+c((K22+K33)K11+K11(K22+K33))+(db)(K22K33+K33K222K23K23)+4e(K12K12+K13K13)

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

  1. die Parameter a bis e durch Funktionen der invarianten Terme ersetzt werden, siehe Hyperelastizität#Orthotrope Hyperelastizität, und/oder
  2. der invariante Term höherer Ordnung im Ansatz zur Formänderungsenergie berücksichtigt wird.[1]:394

Beispiel

Zug einer transversal isotropen Materialprobe mit einer Kraft F in einem Winkel α zur Vorzugsrichtung
Beispiel für die Richtungsabhängigkeit des Elastizitätsmoduls und der Querdehnzahlen bei transversaler Isotropie

Ein transversal isotroper linear elastischer Werkstoff habe die Kennwerte

E1=2000MPaE2=1000MPaG12=700MPaG23=350MPaν12=0,25

Die Stabilitätskriterien werden erfüllt:

E1,E2,G12,G23>0|ν23|=ν23=E22G231=0,4285<1|ν12|=0,25<E1E2=1,4142|ν21|=ν12E2E1=0,125<E2E1=0,70711ν12ν21=0,96875>01ν232=0,8163>0(1ν232ν12ν21)(1+ν23)=0,7270>0.

Wird eine Probe dieses Materials wie im oberen Bild in einem Winkel α zur Vorzugsrichtung einaxial belastet, würde man den Elastizitätsmodul und die Querdehnzahlen, wie im unteren Bild gezeigt, messen. Bei Isotropie wären die Kurven konzentrische Kreise.

Siehe auch

Weblinks

Wiktionary: transversal – Bedeutungserklärungen, Wortherkunft, Synonyme, Übersetzungen
Wiktionary: isotrop – Bedeutungserklärungen, Wortherkunft, Synonyme, Übersetzungen

Einzelnachweise

  1. 1,0 1,1 1,2 1,3 1,4 1,5 1,6 P. Haupt: Continuum Mechanics and Theory of Materials. Springer, 2002, ISBN 978-3-642-07718-0, doi:10.1007/978-3-662-04775-0.
  2. Helmut Schürmann: Konstruieren mit Faser-Kunststoff-Verbunden. 2. Auflage. Springer, 2008, S. 182 f.
  3. H. Altenbach: Kontinuumsmechanik: Einführung in die materialunabhängigen und materialabhängigen Gleichungen. Springer, 2012, ISBN 3-642-24119-0.
  4. Die Fréchet-Ableitung einer Funktion f nach x ist der beschränkte lineare Operator A der – sofern er existiert – in alle Richtungen h dem Gâteaux-Differential entspricht, also
    A(h)=ddsf(x+sh)|s=0=lims0f(x+sh)f(x)sh
    gilt. Darin ist sR,f,xundh skalar-, vektor- oder tensorwertig aber x und h gleichartig. Dann wird auch
    A=fx
    geschrieben.

Literatur

  • 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.