Meine Merkliste

Home
Lexikon
Klassische_Laminattheorie

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.

Weiteres empfehlenswertes Fachwissen

Empfindlichkeitsprüfung von Laborwaagen

Lassen Sie sich von statischen Aufladungen nicht die Wägegenauigkeit stören.

Sicherer Wägebereich zur Sicherstellung genauer Resultate

Inhaltsverzeichnis

1 Annahmen
2 Berechnungsablauf
3 Äußere Lasten
- 3.1 Scheibenlasten
- 3.2 Plattenlasten
4 Scheiben-Platten-Steifigkeitsmatrix
5 Grenzen der Theorie
6 Anwendung in der Finite-Elemente-Methode
7 Berechnungsprogramme
8 Literatur

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 keine 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 ( $σ z = τ z x = τ z y = 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 i j$ jeder UD-Schicht (Elastizitätsgesetz der Unidirektionale Schicht)
Schichtwinkel $α$ 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.

Transformation der $N$ Scheibensteifigkeitsmatrizen der UD-Schichten ins globale $x, y$ -System $Q_k\rightarrow \hat Q_k$
Berechnung der Scheibensteifigkeitsmatrix $A$ , Plattensteifigkeitsmatrix $D$ und Koppelsteifigkeitsmatrix $B$ . Zusammensetzen der Scheiben-Platten-Steifigkeitsmatrix
Invertieren der Scheiben-Platten-Steifigkeitsmatrix
Berechnen der globalen Dehnungen $\hat \epsilon$ und Krümmungen $\hat \kappa$
Transformieren der globalen Dehnung in die Dehnungen jeder UD-Schicht im Schichtkoordinatensystem
Berechnung der Spannungen in jeder UD-Schicht im Schichtkoordinatensystem mit Hilfe der Scheibensteifigkeitsmatrix der UD-Schicht $Q i j$
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 folgende 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 i j$

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 i j$

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^2_{k}} {12}} & + & \underbrace{(z_k-\frac {t_k} {2})^2)} \\ {} & {}\rm{Biegesteifigkeit} & {} & {}\rm{Steineranteil} \end{matrix}$

Koppelsteifigkeitsmatrix $B i j$

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_{26}&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 Eigenschaften 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 i j = 0$ )
UD-Schicht außerhalb der Symmetrieachsen: anisotrop als Scheibe ( $A_{1,6},A_{2,6}\ne 0$ ), anisotrop als Platte ( $D_{1,6},D_{2,6}\ne 0$ ), keine Scheiben-Platten-Kopplung ( $B i j = 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} \ne 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 i j = 0$ )

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

Ingenieurskonstanten

Die Ingenieurskonstanten $\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_{1,2}}{\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 jedoch zulässig.

Berechnungsprogramme

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

LamiCens (kostenfrei) [1]
ESAComp (kommerziell)
Compositor (kommerziell)
AlfaLam (lineares Werkstoffgesetz), AlfaLam.nl (nichtlineares Werkstoffgesetz), beide Programme sind kostenfrei [2]
...

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

Kategorie: Verbundwerkstoff

Dieser Artikel basiert auf dem Artikel Klassische_Laminattheorie aus der freien Enzyklopädie Wikipedia und steht unter der GNU-Lizenz für freie Dokumentation. In der Wikipedia ist eine Liste der Autoren verfügbar.

Klassische Laminattheorie

Weiteres empfehlenswertes Fachwissen

Empfindlichkeitsprüfung von Laborwaagen

Lassen Sie sich von statischen Aufladungen nicht die Wägegenauigkeit stören.

Sicherer Wägebereich zur Sicherstellung genauer Resultate

Inhaltsverzeichnis

Annahmen

Berechnungsablauf

Äußere Lasten

Scheibenlasten

Plattenlasten

Scheiben-Platten-Steifigkeitsmatrix

Scheibensteifigkeitsmatrix $A i j$

Plattensteifigkeitsmatrix $D i j$

Koppelsteifigkeitsmatrix $B i j$

Interpretation

Ingenieurskonstanten

Grenzen der Theorie

Anwendung in der Finite-Elemente-Methode

Berechnungsprogramme

Literatur

Über chemie.de

Über LUMITOS

Werben bei LUMITOS

Klassische Laminattheorie

Weiteres empfehlenswertes Fachwissen

Empfindlichkeitsprüfung von Laborwaagen

Lassen Sie sich von statischen Aufladungen nicht die Wägegenauigkeit stören.

Sicherer Wägebereich zur Sicherstellung genauer Resultate

Inhaltsverzeichnis

Annahmen

Berechnungsablauf

Äußere Lasten

Scheibenlasten

Plattenlasten

Scheiben-Platten-Steifigkeitsmatrix

Scheibensteifigkeitsmatrix Aij

Plattensteifigkeitsmatrix Dij

Koppelsteifigkeitsmatrix Bij

Interpretation

Ingenieurskonstanten

Grenzen der Theorie

Anwendung in der Finite-Elemente-Methode

Berechnungsprogramme

Literatur

Über chemie.de

Über LUMITOS

Werben bei LUMITOS

Scheibensteifigkeitsmatrix $A i j$

Plattensteifigkeitsmatrix $D i j$

Koppelsteifigkeitsmatrix $B i j$