Klassische Laminattheorie

Die klassische Laminattheorie ist ein Verfahren zur Berechnung der Scheiben- und Plattensteifigkeit so wie der Spannungen eines ebenen Mehrschichtenverbunds. Die Schichten des Verbunds bestehen für gewöhnlich aus orthotropen Faser-Kunststoff-Verbunden bzw. unidirektionalen Schichten. Es können aber auch isotrope Schichten behandelt werden. Von besonderer Bedeutung ist die Berechnung der Kopplung zwischen Scheibenbelastungen und Plattenverformungen. Bei vorhandener Kopplung kann sich ein Bauteil z.B. unter Zug krümmen, man spricht umgangssprachlich dann von Verzug.

Die klassische Laminattheorie ist die Grundlage einer Vielzahl von Berechnungsprogrammen für faserverstärkte Kunststoffe. Für eine Festigkeitsberechnung ist die Ermittlung der Schichtspannungen eines Mehrschichtenverbunds nach der klassischen Laminattheorie unerlässlich. Die klassische Laminattheorie wird auch im deutschsprachigen Raum, in Anlehnung an den englischen Begriff classical lamination theory, als CLT abgekürzt.

Annahmen

Die klassische Laminattheorie beruht in großen Zügen auf der Kirchhoff'schen Plattentheorie. Sie gilt nur für einen infinitesimalen ungestörten Ausschnitt. Es werden also kleine Laminatrandeffekte oder Lasteinleitungsprobleme berücksichtigt. Im Einzelnen gelten die folgenden Annahmen:

• Das Elastizitätsgesetz der Einzelschichten ist ideal linear elastisch.
• Das Laminat ist dünn (Dicke $t$ ist klein gegenüber den restlichen Abmessungen)
• Die Laminatdicke $t$ ist konstant
• Die Theorie I. Ordnung ist gültig (kleine Verformungen)
• Die Bernoullische Annahmen sind gültig (ebene Querschnitte, schubstarr in Dickenrichtung)
• Der Spannungszustand ist aufgrund der Dünnwandigkeit eben ($\sigma _{z}=\tau _{zx}=\tau _{zy}=0$).
• Die Schichten sind ideal miteinander verklebt.
• Das Laminat liegt in der $x,y$-Ebene.

Berechnungsablauf

Vor der Berechnung müssen folgende Größen der Schichten festgelegt werden:

• Steifigkeitsmatrix $Q_{ij}$ jeder UD-Schicht (Elastizitätsgesetz der Unidirektionale Schicht)
• Schichtwinkel $\alpha$ der unidirektionalen Schichten. Bei isotropen Schichten ist der Schichtwinkel beliebig.
• Schichtdicken $t$
• Schichtreihenfolge

Der Konstrukteur versucht die obige Größen so zu wählen, daß durch die äußeren Belastungen des Laminats eine möglichst günstige Beanspruchung jeder UD-Schicht entsteht. Des weiteren muss er dafür sorgen, daß er die geforderten Steifigkeitsanforderungen erfüllt. Der nachfolgende Berechnungsgang wird daher oft iterativ durchlaufen.

1. Transformation der $N$ Scheibensteifigkeitsmatrizen der UD-Schichten ins globale $x,y$-System $Q_{k}\rightarrow {\hat {Q}}_{k}$
2. Berechnung der Scheibensteifigkeitsmatrix $A$, Plattensteifigkeitsmatrix $D$ und Koppelsteifigkeitsmatrix $B$. Zusammensetzen der Scheiben-Platten-Steifigkeitsmatrix
3. Invertieren der Scheiben-Platten-Steifigkeitsmatrix
4. Berechnen der globalen Dehnungen ${\hat {\epsilon }}$ und Krümmungen ${\hat {\kappa }}$
5. Transformieren der globalen Dehnung in die Dehnungen jeder UD-Schicht im Schichtkoordinatensystem
6. Berechnung der Spannungen in jeder UD-Schicht im Schichtkoordinatensystem mit Hilfe der Scheibensteifigkeitsmatrix der UD-Schicht $Q_{ij}$
7. Es folgt für gewöhnlich eine Festigkeitsanalyse mittels Bruchkriterien für Faserkunststoffverbunde.

Als Nebenprodukt erhält man aus der Inversen der Scheiben-Platten-Steifigkeitsmatrix die Ingenieurskonstanten des geschichteten Verbunds.

Äußere Lasten

In der klassischen Laminattheorie wird nicht mit Spannungen, als äußere Lasten, gerechnet sondern mit deren Flüssen. Ein Kraft- oder Momentenfluss ist eine breitenbezogene Größe.

Bei isotropen Scheiben und Platten führen Scheibenlasten nur zu Scheibenverformungen (Dehnungen und Schiebungen). Bei geschichteten Laminaten können Scheibenlasten auch zu Plattenverformungen führen (Krümmungen und Drillungen). Daher ist die folgene Unterscheidung in Scheiben- und Plattenlasten notwendig.

Scheibenlasten

Scheibenlasten sind Normalspannungen und Schub in der Laminatebene bzw. deren Flüsse. $z$ entspricht der Koordinatenrichtung normal zur Laminatebene.

${\begin{bmatrix}n_{x}\\n_{y}\\n_{xy}\end{bmatrix}}=\int _{t}^{}{\begin{bmatrix}\sigma _{x}\\\sigma _{y}\\\tau _{xy}\end{bmatrix}}\,\mathrm {d} z$

Plattenlasten

Die Plattenlasten bestehen aus den Biegemomentflüssen sowie dem Drillmomentenfluss. Das Drillmoment kann als Torsionsmoment interpretiert werden.

${\begin{bmatrix}m_{x}\\m_{y}\\m_{xy}\end{bmatrix}}=-\int _{t}^{}{\begin{bmatrix}\sigma _{x}\\\sigma _{y}\\\tau _{xy}\end{bmatrix}}\cdot z\,\mathrm {d} z$

Scheiben-Platten-Steifigkeitsmatrix

Die Scheiben-Platten-Steifigkeitsmatrix beschreibt das elastische Verhalten des gesamten Laminats. Sie setzt sich aus Drei Untermatrizen zusammen, der Scheibensteifigkeitsmatrix A, der Plattensteifigkeitsmatrix D und der Koppelmatrix B, die die beiden ersten Matrizen koppelt.

Scheibensteifigkeitsmatrix $A_{ij}$

Schematische Darstellung der Dehnungs-Schiebungs-Kopplung durch die A16 oder A26 Terme

Die Scheibensteifigkeitsmatrix ergibt sich aus der Parallelschaltung der Scheibensteifigkeiten ${\hat {Q}}_{ij}$ aller $N$ Einzelschichten (unidirektionale Schicht) unter Gewichtung der Schichtdicke der Einzelschicht.

$A_{ij}=\sum _{k=1}^{N}{\hat {Q}}_{ij,k}\cdot t_{k}$

Plattensteifigkeitsmatrix $D_{ij}$

Schematische Darstellung der Biege-Drill-Kopplung durch die D16 oder D26 Terme

Die Plattensteifigkeitsmatrix ergibt sich aus der Parallelschaltung der Biegesteifigkeiten aller Einzelschichten (unidirektionale Schicht) zuzüglich deren Steineranteil.

${\begin{matrix}D_{ij}=\sum _{k=1}^{N}{\hat {Q}}_{ij,k}\cdot t_{k}\cdot (&\underbrace {\frac {t_{k}^{2}}{12}} &+&\underbrace {(z_{k}-{\frac {t_{k}}{2}})^{2})} \\{}&{}{\rm {Biegesteifigkeit}}&{}&{}{\rm {Steineranteil}}\end{matrix}}$

Koppelsteifigkeitsmatrix $B_{ij}$

Schematische Darstellung der Dehnungs-Krümmungs-Kopplung durch die B11,B12 oder B22 Terme

Die Koppelsteifigkeitsmatrix setzt sich aus den Scheibensteifigkeiten der Einzelschichten, gewichtet mit deren statischem Moment, zusammen. Daraus ergibt sich, daß für symmetrisch geschichtete Laminate die Kopplung zwischen Scheibe und Platte verschwindet.
Die Kopplung zwischen Dehnungen und Krümmungen wird bei Bimetallschaltern ausgenutzt.

${\begin{matrix}B_{ij}=-\sum _{k=1}^{N}{\hat {Q}}_{ij,k}\cdot t_{k}\cdot &(\underbrace {z_{k}-{\frac {t_{k}}{2}})} \\{}&{}{\rm {stat.Moment}}\end{matrix}}$

Interpretation

Zusammengesetzt schreibt sich das Elastizitätsgesetz des Scheiben-Platten-Elements, bezogen auf die Neutralebene (Index 0), als

${\begin{bmatrix}n_{x}\\n_{y}\\n_{xy}\\m_{x}\\m_{y}\\m_{xy}\end{bmatrix}}={\begin{bmatrix}A_{11}&A_{12}&A_{16}&B_{11}&B_{12}&B_{16}\\A_{12}&A_{22}&A_{16}&B_{12}&B_{22}&B_{26}\\A_{16}&A_{26}&A_{66}&B_{16}&B_{26}&B_{66}\\B_{11}&B_{12}&B_{16}&D_{11}&D_{12}&D_{16}\\B_{12}&B_{22}&B_{16}&D_{12}&D_{22}&D_{26}\\B_{16}&B_{26}&B_{66}&D_{16}&D_{26}&D_{66}\\\end{bmatrix}}{\begin{bmatrix}\epsilon _{x}\\\epsilon _{y}\\\gamma _{xy}\\\kappa _{x}\\\kappa _{y}\\\kappa _{xy}\end{bmatrix}}_{0}$

Ausgewählte Laminate und deren Eingenschaften des Elastizitätsgesetzes:

• UD-Schicht: orthotrop als Scheibe ($A_{1,6},A_{2,6}=0$), orthotrop als Platte ($D_{1,6}=D_{2,6}=0$), keine Scheiben-Platten-Kopplung ($B_{ij}=0$)
• UD-Schicht außerhalb der Symmetrieachsen: anisotrop als Scheibe ($A_{1,6},A_{2,6}\neq 0$), anisotrop als Platte ($D_{1,6},D_{2,6}\neq 0$), keine Scheiben-Platten-Kopplung ($B_{ij}=0$)
• ausgeglichener Winkelverbund oder Kreuzverbund: orthotrop als Scheibe ($A_{1,6},A_{2,6}=0$), orthotrop als Platte ($D_{1,6},D_{2,6}=0$), Scheiben-Platten-Kopplung ($B_{11},B_{22}\neq 0$)
• ausgeglichener Winkelverbund oder Kreuzverbund, symmetrisch geschichtet: orthotrop als Scheibe ($A_{1,6},A_{2,6}=0$), orthotrop als Platte ($D_{1,6},D_{2,6}=0$), keine Scheiben-Platten-Kopplung ($B_{ij}=0$)

Die Platten-Orthotropie wird aufgrund der fehlenden Wölbungs-Drillungs-Kopplung auch Wölbungsorthotropie genannt.

Ingenieurskonstanten

Die Ingeneurskonstanten ${\hat {E}}_{x},{\hat {E}}_{y},{\hat {G}}_{xy},\nu _{xy}$ erhält man aus der inversen Scheiben-Platten-Steifigkeitsmatrix. ${\begin{bmatrix}{\bar {A}}&{\bar {B}}\\{\bar {B}}&{\bar {D}}\end{bmatrix}}={\begin{bmatrix}A&B\\B&D\end{bmatrix}}^{-1}$

• ${\hat {E}}_{x}={\frac {1}{{\bar {A}}_{1,1}}}$

• ${\hat {E}}_{y}={\frac {1}{{\bar {A}}_{2,2}}}$

• ${\hat {G}}_{xy}={\frac {1}{{\bar {A}}_{6,6}}}$

• ${\hat {\nu }}_{yx}=-{\frac {{\bar {A}}(A_{12})}{{\bar {A}}_{1,1}}}$

Grenzen der Theorie

Die Gültigkeitsgrenzen der CLT ergeben sich hauptsächlich aus den Annahmen der Kirchhoffschen Plattentheorie (siehe: Annahmen). Für dicke und gedrungene Platten nimmt die Qualität der Lösung ab. Auch sehr schubweiche Schichten unter Querkraftbelastung sollten nicht mit der CLT berechnet werden. Grund ist der steigende Anteil der Schubabsenkung gegenüber der Biegeabsenkung der Platte.

Bei reiner Scheibenbelastung beeinflusst die Scheibendicke die Qualität der Lösung nicht. Bei sehr dicken Einzelschichten kann es zu Abweichungen kommen, da bei der Scheibenbelastung die Verformung durch den interlaminaren Schub nicht berücksichtigt wird. Für die Praxis ist dies jedoch kaum von Bedeutung, da Laminate gewöhnlich aus dünnen Einzelschichten bestehen.

Anwendung in der Finite-Elemente-Methode

Die Ingenieurskonstanten können direkt in der Finite-Elemente-Methode in Kombination mit Volumenelementen benutzt werden, um das globale Verhalten eines geschichteten Verbunds zu simulieren. Sind die Schichtspannungen von Interesse müssen diese, wie in der manuellen CLT-Rechnung, aus den globalen Dehnungen berechnet werden. Die in der FEM berechneten Spannungen sind die globalen Spannungen des Verbunds und entsprechen nicht den Schichtspannungen. Die globale Spannung darf nicht direkt in die Schichtspannungen transformiert werden.

Bei der Formulierung über die Ingenieurskonstanten, unter Verwendung von 3-dimensionalen Volumenelementen, werden weder Schiebe-Dehn-Kopplungen (Anisotropie als Scheibe) noch Biege-Dehn-Kopplungen berücksichtigt. Auch die Anisotropie als Platte kann nicht abgebildet werden. Es existieren spezielle Schalenelemente, die diese Kopplungen berücksichtigen können.

Für den in der Praxis häufigen Fall der symmetrisch geschichteten, ausgeglichenen Laminate (Orthotropie) ist die Modellierung über Volumenelemente und Ingenieurskonstanten jeodch zulässig.

Berechnungsprogramme

Es existieren eine Vielzahl von CLT-Programmen, die zum Teil konstenfrei erhältlich sind. Sie basieren zumeist auf Tabellenkalkulationsprogrammen. Teilweise ist eine Datenbank mit Faserhalbzeugen und Matrixsystemen angeschlossen.

• Lamicens (kostenfrei) [1]
• ESAComp (kommerziell)
• Compositor (kommerziell)
• ...

Literatur

J. Wiedemann: Leichtbau, Band 1: Elemente. Springer-Verlag, Berlin 1986. ISBN 3-5-40164049

H. Altenbach, J. Altenbach, R. Rikard: Einführung in die Mechanik der Laminattragwerke und Sandwichstragwerke. Deutscher Verlag für Grundstoffindustrie, 1996. ISBN 3-3-42006811