Glättung der Flächenschnittgrößen in RFEM

Im Programm bestehen verschiedene Möglichkeiten für die Anzeige der ungeglätteten oder geglätteten Ergebnisse. Sie sind im RFEM-Handbuch beschrieben.

• Online-Handbuch RFEM 5

Auch in einem früheren Fachbeitrag sind Glättungsvarianten zur Auswertung von Flächenschnittgrößen und -spannungen vorgestellt.

• KB 001571 | Verlauf der Schnittgrößen/Spannungen bei Flächen - Glättungsvarianten

Folgender Beitrag erläutert anhand eines Beispiels, wie die FE-Schnittgrößen bei der Glättung für die Werteermittlung herangezogen werden. Bei diesem Prozess spielt neben der Position eines finiten Elements auch die Plattentheorie – Mindlin oder Kirchhoff – eine Rolle.

Beispiel: Platte nach Mindlin und Kirchhoff

1 Beispiel

Eine Platte mit den Abmessungen 4 m ⋅ 3 m wird als zweidimensionales Modell erzeugt. Die Platte ist beidseits an den kürzeren Randlinien gelenkig gelagert. Die Verdrehung der längeren Randlinien ist durch eine Einspannung behindert. Um Querdehnungeffekte zu unterdrücken, wird die Querdehnungszahl des Materials zu null gesetzt. Für einen Vergleich der Schnittgrößenermittlung nach den Plattenbiegetheorien von Kirchhoff und Mindlin wird eine Plattendicke von 80 cm gewählt. Damit ist ein Verhältnis d/L = 0.2 im Grenzbereich zwischen den beiden Biegetheorien gegeben (siehe FAQ 003158).

• FAQ | Ist es ausreichend die Berechnung mit Flächenelementen durchzuführen oder sollen besser Volumenelemente zum Einsatz kommen?

Die Maschenweite des FE-Netzes wird der Einfachheit halber mit 1 m festgelegt. Somit liegen 4 ⋅ 3 finite Elemente vor.

Auf die Platte wirkt eine Last von p = 5 kN/m². Das Eigengewicht wird im Lastfall nicht angesetzt.

Nicht geglättete Schnittgrößen

Die Ergebnisse des Rechenkerns für die Schnittgrößen v_x und m_x werden für einen Längsschnitt betrachtet, der in der Mitte der Platte verläuft. Der Verlauf der elementweisen, nicht geglätteten Schnittgrößen kann mit der Anzeigeoption "Nicht durchlaufend" dargestellt werden. Nach den Biegetheorien von Mindlin und Kirchhoff liegen die im Bild 02 gezeigten Verläufe vor.

2 Nicht geglättete Ergebnisse nach Mindlin (links) und Kirchhoff (rechts)

Standardglättung für innenliegende Knoten

Für FE-Knoten, die innerhalb einer Fläche liegen, wird zunächst das arithmetische Mittel aus den Knotenwerten der benachbarten finiten Elemente gebildet. Für diesen Ansatz müssen die Elemente in einer Fläche und auf der gleichen Seite einer möglichen internen Linie liegen. Der Knoten darf kein benutzerdefinierter Knoten in der Fläche sein.

Für die Unstetigkeitsstelle bei x = 1.00 m sind dies intern folgende Werte (das Ergebnis dieses Schritts ist für die Ausgabe nicht verfügbar).

• Mindlin:

  $${\mathrm{v}}_{\mathrm{x}}\left(1\right) = \frac{7.5 · 2.5}{2} = 5$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(1\right) = \frac{8.772 · 5.351}{2} = 7.061$$

  Standardglättung für innenliegende Knoten nach Mindlin
• Kirchhoff:

  $${\mathrm{v}}_{\mathrm{x}}\left(1\right) = \frac{7.785 · 2.595}{2} = 5.19$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(1\right) = \frac{7.954 · 7.904}{2} = 7.929$$

  Standardglättung für innenliegende Knoten nach Kirchhoff

Glättung für Randknoten

An den Rändern einer Fläche (Glättungsoption "Durchlaufend innerhalb einer Fläche") oder eines Modells ("Durchlaufend gesamt") liegen keine benachbarten finiten Elemente vor, die für eine Mittelung der Knotenwerte dienen könnten. Daher wird ein anderer Ansatz gewählt, der in zwei Schritten erfolgt.

Im ersten Schritt werden die gemittelten Werte für diejenigen Knoten berechnet, die sich nicht am Rande einer Fläche oder des Modells befinden. Die Werte der Knoten an den Rändern werden so berechnet, dass die ursprünglichen Werte in den Mitten der finiten Elemente erhalten bleiben. Im zweiten Schritt werden dann die gemittelten Werte für die Randknoten bestimmt.

Für den Randknoten bei x = 0.00 m ergeben sich folgende Werte.

• Mindlin:

  $${\mathrm{v}}_{\mathrm{x}}\left(0\right) = 7.5 · (7.5 - 5) = 10$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}(0.5) = \frac{5.351 · 2.982}{2} = 4.167$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(0\right) = 4.167 · (4.167 - 7.061) = 1.272$$

  Glättung für Randknoten nach Mindlin
• Kirchhoff:

  $${\mathrm{v}}_{\mathrm{x}}\left(0\right) = 7.785 · (7.785 - 5.19) = 10.381$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}(0.5) = \frac{7.954 · 0.379}{2} = 4.167$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(0\right) = 4.167 · (4.167 - 7.929) = 0.404$$

  Glättung für Randknoten nach Kirchhoff

Querkräfte

Nach Mindlin werden die Querkräfte als erste Ableitung der Durchbiegung berechnet.

Diese Werte werden vom Rechenkern bestimmt und direkt verwendet. Die Glättung der Querkräfte erfolgt dann wie oben beschrieben je nach Lage des FE-Knotens in der Fläche.

Bei der Biegetheorie nach Kirchhoff werden die Querkräfte als dritte Ableitung der Durchbiegung berechnet.

Diese Werteermittlung über die dritte Ableitung ist mit einem erheblichen Genauigkeitsverlust verbunden. Aus diesem Grund werden die Querkräfte des Rechenkerns nicht verwendet, sondern durch einen verbesserten Ansatz aus den Ableitungen der Momente bestimmt.

Die in den Gleichungen dargestellten Momente m_x, m_y und m_xy repräsentieren die geglätteten Werte, die nach den oben beschriebenen Verfahren bestimmt werden. Damit ergeben sich genauere Werte für die Querkräfte als vonseiten des Rechenkerns.

Momente

Während die Momentenwerte in den Integrationspunkten den theoretischen Werten entsprechen, so führt die Extrapolation bei der Glättung über eine hyperbolische Paraboloidfunktion zu einem Verlust an Genauigkeit. Ein hyperbolisches Paraboloid drückt die Momentenverteilung im Element nur angenähert aus. Aus diesem Grund kommt ein verbesserter Algorithmus zur Anwendung, der die Extrapolation durch eine weiterentwickelte Querkraft-Integralmethode ersetzt. Da die Querkraftverteilung im Flächenelement der Form eines hyperbolischen Paraboloids entsprechend angenommen wird (Fläche zweiter Ordnung), stellt das Integral dieser Fläche eine Fläche dritter Ordnung dar, die den Momentenverlauf mit einer größeren Genauigkeit abbildet.

Dieser Ansatz korrespondiert mit den oben genannten Gleichungen zur Ermittlung der Querkräfte über die Ableitungen der Momente. Die Momente werden dann mit folgenden Gleichungen bestimmt.

Die Integrationskonstanten m_x,0 und m_y,0 werden aus den Bedingungen des Wertes in Elementmitte berechnet.

Die Werte der Schnittgrößen werden nach den oben beschriebenen Verfahren bestimmt, die Integrale durch eine numerische Integration. Die so ermittelten Momente werden dann nach dem Verfahren für innenliegende Knoten oder Randknoten geglättet.

Bei der Biegetheorie nach Mindlin ergeben sich im Beispiel folgende Werte.

• Erstes Element:

  $${\mathrm{v}}_{\mathrm{x}}\left(\mathrm{x}\right) = 10 - 5 \mathrm{x}$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(\mathrm{x}\right) = \int {\mathrm{v}}_{\mathrm{x}}\mathrm{dx} {\mathrm{m}}_{\mathrm{x},0} = 10 \mathrm{x} - \frac{5}{2} {\mathrm{x}}^2 {\mathrm{m}}_{\mathrm{x},0}$$

  
  

  $${\mathrm{m}}_{\mathrm{x},\mathrm{C}} = {\mathrm{m}}_{\mathrm{x}}(0.5) = 4.167 = {\int }_0^{1}10 \mathrm{x} - \frac{5}{2} {\mathrm{x}}^2 {\mathrm{m}}_{\mathrm{x},0} = {\left[\frac{10}{2} {\mathrm{x}}^2 - \frac{5}{6} {\mathrm{x}}^3 {\mathrm{m}}_{\mathrm{x},0}\mathrm{x}\right]}_0^1 = 4.167 {\mathrm{m}}_{\mathrm{x},0}$$

  
  

  $${\mathrm{m}}_{\mathrm{x},0} = {\mathrm{m}}_{\mathrm{x}}\left(0\right) = 4.167 - 4.167 = 0$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}{\left(1\right)}^{-} = 10 - \frac{5}{2} 0 = 7.5$$

  Erstes Element bei der Biegetheorie nach Mindlin
• Zweites Element:

  $${\mathrm{v}}_{\mathrm{x}}\left(\mathrm{x}\right) = 5 - 5 \mathrm{x}$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(\mathrm{x}\right) = \int {\mathrm{v}}_{\mathrm{x}}\mathrm{dx} + {\mathrm{m}}_{\mathrm{x},0} = 5 \mathrm{x} - \frac{5}{2} {\mathrm{x}}^2 + {\mathrm{m}}_{\mathrm{x},0}$$

  
  

  $${\mathrm{m}}_{\mathrm{x},\mathrm{C}} = {\mathrm{m}}_{\mathrm{x}}(1.5) = 9.167 = {\int }_0^{1}5 \mathrm{x} - \frac{5}{2} {\mathrm{x}}^2 + {\mathrm{m}}_{\mathrm{x},0} = {\left[\frac{5}{2} {\mathrm{x}}^2 - \frac{5}{6} {\mathrm{x}}^3 + {\mathrm{m}}_{\mathrm{x},0}\mathrm{x}\right]}_0^1 = 1.667 + {\mathrm{m}}_{\mathrm{x},0}$$

  
  

  $${\mathrm{m}}_{\mathrm{x},0} = {\mathrm{m}}_{\mathrm{x}}(1{)}^{+} = 9.167 - 1.667 = 7.5$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(1\right) = \frac{{\mathrm{m}}_{\mathrm{x}}(1{)}^{-} + {\mathrm{m}}_{\mathrm{x}}(1{)}^{+}}{2} = \frac{7.5 + 7.5}{2} = 7.5$$

  Zweites Element bei der Biegetheorie nach Mindlin

3 Geglättete Schnittgrößenverläufe v-x und m-x nach Mindlin

Bei der Biegetheorie nach Kirchhoff ergeben sich folgende Werte.

• Erstes Element:

  $${\mathrm{v}}_{\mathrm{x}}\left(\mathrm{x}\right) = 10.038 - 5.025 \mathrm{x}$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(\mathrm{x}\right) = \int {\mathrm{v}}_{\mathrm{x}}\mathrm{dx} {\mathrm{m}}_{\mathrm{x},0} = 10.038 \mathrm{x} - \frac{5.025}{2} {\mathrm{x}}^2 {\mathrm{m}}_{\mathrm{x},0}$$

  
  

  $${\mathrm{m}}_{\mathrm{x},\mathrm{C}} = {\mathrm{m}}_{\mathrm{x}}(0.5) = 4.167 = {\int }_0^{1}10.038 \mathrm{x} - \frac{5.025}{2} {\mathrm{x}}^2 {\mathrm{m}}_{\mathrm{x},0} = {\left[\frac{10.038}{2} {\mathrm{x}}^2 - \frac{5.025}{6} {\mathrm{x}}^3 {\mathrm{m}}_{\mathrm{x},0}\mathrm{x}\right]}_0^1 = 4.182 {\mathrm{m}}_{\mathrm{x},0}$$

  
  

  $${\mathrm{m}}_{\mathrm{x},0} = {\mathrm{m}}_{\mathrm{x}}\left(0\right) = 4.167 - 4.182 = - 0.015 {\mathrm{m}}_{\mathrm{x}}{\left(1\right)}^{-} = 10.038 - \frac{5.025}{2} - 0.015 = 7.511$$

  Erstes Element bei der Biegetheorie nach Kirchhoff
• Zweites Element:

  $${\mathrm{v}}_{\mathrm{x}}\left(\mathrm{x}\right) = 5.013 - 5.013 \mathrm{x}$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(\mathrm{x}\right) = \int {\mathrm{v}}_{\mathrm{x}}\mathrm{dx} + {\mathrm{m}}_{\mathrm{x},0} = 5.013 \mathrm{x} - \frac{5.013}{2} {\mathrm{x}}^2 + {\mathrm{m}}_{\mathrm{x},0}$$

  
  

  $${\mathrm{m}}_{\mathrm{x},\mathrm{C}} = {\mathrm{m}}_{\mathrm{x}}(1.5) = 9.167 = {\int }_0^{1}5.013 \mathrm{x} - \frac{5.013}{2} {\mathrm{x}}^2 + {\mathrm{m}}_{\mathrm{x},0} = {\left[\frac{5.013}{2} {\mathrm{x}}^2 - \frac{5.013}{6} {\mathrm{x}}^3 + {\mathrm{m}}_{\mathrm{x},0}\mathrm{x}\right]}_0^1 = 1.671 + {\mathrm{m}}_{\mathrm{x},0}$$

  
  

  $${\mathrm{m}}_{\mathrm{x},0} = {\mathrm{m}}_{\mathrm{x}}(1{)}^{+} = 9.167 - 1.671 = 7.496$$

  
  

  $${\mathrm{m}}_{\mathrm{x}}\left(1\right) = \frac{{\mathrm{m}}_{\mathrm{x}}(1{)}^{-} + {\mathrm{m}}_{\mathrm{x}}(1{)}^{+}}{2} = \frac{7.511 + 7.496}{2} = 7.503$$

  Zweites Element bei der Biegetheorie nach Kirchhoff

4 Geglättete Schnittgrößenverläufe v-x und m-x nach Kirchhoff

Beim Glättungsprozess laufen die oben beschriebenen Schritte sukzessive im Programm ab. Die geglätteten Werte lassen sich grafisch mit der Anzeigeoption "Durchlaufend innerhalb Flächen" darstellen.

Knowledge Base

KB 001606 | Glättung der Flächenschnittgrößen in RFEM