Orthotrope Materialgesetze

Im Regelfall ist bei den orthotropen Materialien das hookesche Gesetz mit $\mathrm\nu\;=\;\frac{\mathrm E}{2\;\mathrm G}\;\;-\;1\;\leq\;0,5$ und $-1\;<\;\mathrm\nu\;<\;0,5$ nicht mehr gültig.

deformazione trasversale

Die folgenden Materialparameter beziehen sich auf zweidimensionale Steifigkeiten und, soweit nicht anders benannt, auf das Material Holz. Es wird ein lokales Achsensystem wie in Bild 01 zugrunde gelegt.

• E_x = Steifigkeit in lokaler x-Richtung der Fläche
• E_y = Steifigkeit in lokaler y-Richtung der Fläche
• G_xz = Schubsteifigkeit in lokaler x-Richtung der Fläche (Dickenrichtung der Platte)
• G_yz = Schubsteifigkeit in lokaler y-Richtung der Fläche (Dickenrichtung der Platte)
• G_xy = Schubsteifkgeit in Scheibenebene
• ν_xy = Querdehnung in x-Richtung
• ν_yx = Querdehnung in y-Richtung

Die Spannungen in Bild 02 stehen mit den hier genannten Steifigkeiten in Beziehung.

tensioni tangenziali

Das Materialgesetz unterliegt folgenden Regeln.

Gleichung 1:

Gleichung 2:
${\mathrm\nu}_\mathrm{xy}\;\geqslant\;{\mathrm\nu}_\mathrm{yx}$

Gleichung 3:

Gleichung 4 (Steifigkeiten in Scheibenebene ESZ):

Das Verhältnis der Dehnungen in obigen Gleichungen unterstreicht die Beziehungen in Bild 01.

Die Steifigkeiten in Scheibenebene berechnen sich wie folgt.

Gleichung 5:

Querdehnung ν

Wie bereits in Bild 01 erläutert, resultieren aus dem weicheren Materialverhalten in der jeweiligen Richtung veränderte Verformungen und Spannungen in dieser Richtung.

Verhältnis der Dehnungen:

Gleichung 6:

Gleichung 7:

Für ${\mathrm\sigma}_\mathrm x\;\neq\;0,\;{\mathrm\sigma}_\mathrm y\;=\;0,\;{\mathrm\sigma}_\mathrm{xy}\;=\;0$ ergeben sich mit dem hookeschen Gesetz folgende Gleichungen.

Gleichung 8:

Gleichung 9:

Gleichung 10:

Gleichung 11:

Gleichung 12:

Gleichung 13:

Steifigkeitsmatrix

Berechnung der globalen Steifigkeitsmatrix der Platte.

Gleichung 14:

Biegeanteile:

Gleichung 15:
${\mathrm D}_{11}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{11}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm x\;\mathrm d^3}{12\;\left(1\;-\;{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}\right)}$

Gleichung 16:
${\mathrm D}_{12}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{12}\;\widehat{=\;}\mathrm{sgn}\left({\mathrm\nu}_\mathrm{xy}\right)\;\sqrt{{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}\;{\mathrm D}_{11}\;{\mathrm D}_{22}}$

Gleichung 17:
${\mathrm D}_{22}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{22}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm y\;\mathrm d^3}{12\;\left(1\;-\;{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}\right)}$

Gleichung 18:
${\mathrm D}_{33}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{33}\;\widehat{=\;}{\mathrm G}_\mathrm{xy}\;\frac{\mathrm d^3}{12}$

Scheibenanteile:

Gleichung 19:
${\mathrm D}_{66}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{11}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm x\;\mathrm d}{1\;-\;{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}}$

Gleichung 20:
${\mathrm D}_{67}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{12}\;\widehat{=\;}\mathrm{sgn}\left({\mathrm\nu}_\mathrm{xy}\right)\;\sqrt{{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}\;{\mathrm D}_{66}\;{\mathrm D}_{77}}$

Gleichung 21:
${\mathrm D}_{77}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{22}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm y\;\mathrm d}{1\;-\;{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}}$

Gleichung 22:
${\mathrm D}_{88}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{33}\;\widehat{=\;}{\mathrm G}_\mathrm{xy}\;\mathrm d$

Schubanteile:

Gleichung 23:

Gleichung 24:

Voraussetzung für diese Gleichungen ist natürlich immer, dass die Steifigkeitsmatrix positiv definiert ist, also alle Eigenwerte der Matrix positiv sind.

Aus diesem Grund überprüft RFEM unter anderem die Definition der Querdehnung nach folgender Gleichung.

Gleichung 25:

Beispiel

Mit folgendem Beispiel (Bild 03) wird das orthotrope Materialverhalten erläutert. Hierfür wird dem orthotropen Material ein isotropes Material gegenübergestellt. Zusätzlich wird die Steifigkeit der orthotropen Platte einmal mit der hohen Steifigkeit in x-Richtung und einmal in y-Richtung definiert.

• 3D Viewer | Google
• 3D Viewer | Babylon
• Realtà virtuale

Sistema

System:

• Plattendicke 200 mm
• Material C24
• Orthotrope Steifigkeiten
  
  $\begin{array}{l}{\mathrm{E}}_{\mathrm{x}} = 1.100,0 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{E}}_{\mathrm{y}} = 37,0 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{G}}_{\mathrm{y}} = 6,9 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{G}}_{\mathrm{x}} = 69,0 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{G}}_{\mathrm{xy}} = 69,0 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{\nu }}_{\mathrm{xy}} = 2,52\\ {\mathrm{\nu }}_{\mathrm{yx}} = 0,085\end{array}$

  Formula 29
• Isotrope Steifigkeiten
  
  $\begin{array}{l}\mathrm{E} = 1.100 \mathrm{kN}/{\mathrm{cm}}^2\\ \mathrm{G} = 500 \mathrm{kN}/{\mathrm{cm}}^2\\ \mathrm{\nu } = 0,1\end{array}$

  Formula 30
• Dimension b = 2,0 m, l = 4,0 m
• Last 20 kN/m²
• FE-Netzgröße 50 cm

Das System ist unverschieblich in vertikaler z-Richtung gelagert. Die Lagerbedingungen in x- und y-Richtung wurden so gewählt, dass keine Zwängungen entstehen.

Die Berechnung erfolgt nach Theorie I. Ordnung mit linear elastischem Materialverhalten und Lagerbedingungen.

Das hookesche Gesetz ergibt mit den gegebenen Werten folgende Querdehnung.

Gleichung 26:
$\mathrm\nu\;=\;\left(\sqrt{{\mathrm E}_\mathrm x\;\cdot\;{\mathrm E}_\mathrm y}\right)\;/\;\left(2\;\mathrm G\right)\;-\;1\;=\;6,97$

Diese hohe Querdehnung ist mit dem gewählten Materialmodel nicht möglich. Mit den Gleichungen aus [1] können die Werte jedoch angepasst werden.

Gleichung 27:

Gleichung 28:

Gleichung 29:

Gleichung 30:

• 3D Viewer | Google
• 3D Viewer | Babylon
• Realtà virtuale

Superfici con rigidezze

Ergebnisse
Wie zu erwarten, ergeben sich die größten Verformungen mit der Orientierung der Steifigkeiten in y-Richtung (Bild 06). Die Lagerreaktion und das Moment der isotropen Platte sind in Bild 05 dargestellt.

• 3D Viewer | Google
• 3D Viewer | Babylon
• Realtà virtuale

Risultati materiale isotropo

Da die Platte mit der hohen Steifigkeit in y-Richtung (E_y = 1.100 kN/cm²) die hohe Tragfähigkeit in dieser Richtung aufweist, sind dort auch die Lagerreaktionen höher (125,4 kNm zu 58,3 kNm).

• 3D Viewer | Google
• 3D Viewer | Babylon
• Realtà virtuale

Deformazioni e reazioni vincolari di placche ortotropiche

Die maximalen Biegemomente ergeben sich für die orthotropen Platten analog für m_x mit der Steifigkeit in x-Richtung und für m_y mit der hohen Steifigkeit in y-Richtung.

Für die Platte mit der hohen Steifigkeit in y-Richtung ist das maximale Biegemoment m_y nahezu in der Mitte der Platte (Bild 07).

Biegemomente der orthotropen Platte mx (links) und my (rechts)

Variation der Querdehnung

Die Querdehnung nach den Dehnungsbeziehungen kann die in der Tabelle aufgeführten maximalen und minimalen Werte erreichen.

L'analisi con spettro di risposta ASCE 7-16 in RFEM specificatamente alla minuto 52:25 fornirà una visione dettagliata del flusso di lavoro in RFEM e RF-DYNAM Pro per l'applicazione della matrice di rigidezza geometrica per tenere conto degli effetti P-Delta secondo ASCE 7.

Die eingangs eingeführte Platte mit der hohen Steifigkeit (E_x = 11.000) wird hierfür mit diesen hohen Querdehnungen definiert. Die anderen Steifigkeiten der Platte bleiben dabei jedoch unverändert.

In Bild 08 sind die Ergebnisse der Variation von ν_xy = 5,44 bis -5,44 aufgeführt.

Für ν_xy = 5,44 sind die Lagerreaktionen qualitativ identisch zum isotropen Materialverhalten. Das Biegemoment vergrößert sich von m_x = 18,1 kNm/m (isotrope Platte) zu m_x = 34,9 kNm/m (orthotrope Platte).

Gegenüber der orthotropen Platte mit üblicheren Querdehnungen (ν_xy = 2,5) verkleinert sich das Biegemoment leicht.

Ergebnisse der Variation gemäß Tabelle (links: νxy = 5,44, Centro: νxy = 0, a destra: νxy = -5,44)

Mit ν_xy = 0 verlagert sich der hohe Ausschlag der Lagerreaktion am freien Ende der Platte zu einem konstanten Wert von 43 kN/m.

Das Moment m_x vergrößert sich auf 38,1 kNm/m. Im Vergleich zum vorherigen Ergebnis (ν_xy =5,44) zeigt sich der Einfluss der Querdehnung. Für ν =0 wird keine Verformung oder Verzerrung aufgrund der Querdehnung hervorgerufen.

Für ν_xy = -5,44 ergibt sich ein Durchschlag am freien Plattenende und es werden die Lagerreaktionen negativ. Das maximale Moment ergibt sich in der Mitte der Platte mit 59,5 kNm/m.

Die Platte verhält sich jetzt mehr wie eine einaxial gespannte Platte ohne das dritte Lager in seiner Längsrichtung.

Anhand Bild 01 und den dort aufgeführten Beziehungen lässt sich dieses Verhalten erklären.

Aufgrund der hohen negativen Querdehnung (ν_xy = -5,44) wird die Platte am freien Rand komplett überdrückt und kann sich daher hier nicht verformen.

Der Einfluss der Orthotropie in y-Richtung liegt hierbei nahe null (E_y ≈ 0).

Zusammenfassung

Mit dem orthotropen Materialmodell in RFEM ist eine nahezu beliebige Definition von Materialparametern möglich. Anhand der Variation der Querdehnungen können sehr unterschiedliche Ergebnisse erreicht werden. Eine Querdehnung nach Modifikation der Werte nach [1] ergibt Werte nahe der Lösung für einen Einfeldträger.

Gleichung 31:

Zu hohe negative Querdehnungen bilden ein geändertes statisches System ab, welches nicht mehr der Modellierung entspricht.