Zur
Startseite
Aufgabe
Für das skizzierte Fachwerk sind die Stabkräfte zu berechnen.
Gegeben: F.
Berechnung mit dem Programm "Statisch bestimmte ebene Fachwerke"
Mit dem Programm "Statisch bestimmte ebene Fachwerke" wurde folgendes Ergebnis erzielt:
|
|
Kontrollrechnung mit Matlab-Femset
| Knoten-Informationen | ||||||
| Kn.-Koord. | Kn.-Lasten | u-verh. | ||||
| Kn. | x | y | Fx | Fy | ux | uy |
| 1 | 1 | 5 | 0 | 0 | 1 | 1 |
| 2 | 0 | 2 | 0 | 0 | 1 | 1 |
| 3 | 3 | 0 | 0 | 0 | 1 | 1 |
| 4 | 1 | 4 | 0 | 0 | 0 | 0 |
| 5 | 4 | 4 | 0 | -1 | 0 | 0 |
| 6 | 2 | 3 | 0 | 0 | 0 | 0 |
| 7 | 3 | 3 | 0 | 0 | 0 | 0 |
| 8 | 1 | 2 | 0 | 0 | 0 | 0 |
| 9 | 2 | 2 | 0 | 0 | 0 | 0 |
| 10 | 3 | 1 | 0 | 0 | 0 | 0 |
Die Rechnung wird mit Matlab-Femset verifiziert. Matlab-Femset erwartet die Beschreibung des physikalischen Modells mit "Knoten- und Element-Informationen", die in Matrizen zusammengestellt sind (vgl. hierzu die Erläuterungen im Kapitel "Der Stab als finites Element").
Nebenstehend ist die komplette Nummerierung der Stäbe und Knoten zu sehen, wie sie für die Berechnung mit dem Programm "Statisch bestimmte ebene Fachwerke" benutzt wurde. Um die Ergebnisse leichter vergleichen zu können, soll für die Berechnung mit Matlab-Femset die gleiche Nummerierung verwendet werden.
Die Koordinaten der Knoten können aus dem Bild entnommen werden (der Ursprung des nach rechts bzw. oben gerichteten Koordinatensystems liegt im unteren Lager). Für die beiden Kräfte wurden Einheitslasten angenommen. In den beiden Spalten ux und uy zeigt eine 1 an, dass diese Verschiebung verhindert ist (Lager).
| Element-Informationen | ||||||
| El. | Kn. | Kn. | El. | Kn. | Kn. | |
| 1 | 1 | 4 | 8 | 4 | 8 | |
| 2 | 2 | 8 | 9 | 8 | 9 | |
| 3 | 3 | 10 | 10 | 6 | 9 | |
| 4 | 4 | 5 | 11 | 7 | 10 | |
| 5 | 4 | 6 | 12 | 9 | 10 | |
| 6 | 5 | 7 | 13 | 8 | 10 | |
| 7 | 6 | 7 | 14 | 5 | 10 | |
Nebenstehend sind die Element-Informationen zusammengestellt (an welchen Knoten sind die jeweiligen Elemente angeschlossen?):
Im folgenden Matlab-Script findet man diese Matrizen wieder. Nur der farblich abgehobene Teil muss dem aktuellen Problem angepasst werden (das sind genau die Matrizen, die das Modell definieren):
clear all
% Knotenkoordinaten:
% 1 2 3 4 5 6 7 8 9 10
xy = [1 5 ; 0 2 ; 3 0 ; 1 4 ; 4 4 ; 2 3 ; 3 3 ; 1 2 ; 2 2 ; 3 1] ;
% Verhinderte Verschiebungen:
kr = [1 1 ; 1 1 ; 1 1 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 0] ;
% Kraefte:
bk = [0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 -1 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 0] ;
% Koinzidenzmatrix ("An welchen Knoten sind die Staebe angeschlossen?"):
% 1 2 3 4 5 6 7 8 9 10 11 12 13 14
km = [1 4 ; 2 8 ; 3 10 ; 4 5 ; 4 6 ; 5 7 ; 6 7 ; 4 8 ; 8 9 ; 6 9 ; 7 10 ; 9 10 ; 8 10 ; 5 10] ;
% Dehnsteifigkeiten (fuer Statik-Problem nicht relevant, aber formal erforderlich):
nElem = 14 ; ep(1:nElem,1) = 1 ;
drawfw2d (xy , km , kr , bk , 'ek') ; % 2D-Fachwerk zeichnen
% Finite-Elemente-Algorithmus:
[succ uv] = feskal_m (xy , km , ep , kr , bk) ;
if succ == 1 % Nach erfolgreicher Rechnung:
sk = Stabkraefte2D (xy , km , ep , uv) ; % Berechnung der Stabkraefte ...
for i = 1:nElem
disp (['Stab ',num2str(i), ': Stabkraft F = ',num2str(sk(i)),' kN']) ;
end
clf ;
drawfw2dsk (xy , km , kr , bk , sk , 'e') ; % ... und in Graphik-Fenster darstellen
end
Nachfolgend sind links die Ergebnisse der Matlab-Rechnung zu sehen, rechts die von Matlab erzeugte Graphik. Die Ergebnisse bestätigen die oben zu sehenden Stabkräfte:
|
|
Das oben gelistete Matlab-Script und alle verwendeten Functions und DLLs stehen zum Download zur Verfügung: FW4.m, Stabkraefte2D.m, drawfw2d.m, drawfw2dsk.m, Elemat_m.dll, feskal_m.dll.
Knotenrundschnittverfahren
An den 7 Knoten 4, 5, 6, 7, 8, 9, 10 (Knotennummerierung wie für die Matlab-femset-Rechnung oben) werden insgesamt 14 Gleichgewichtsbedingungen formuliert.
Beim Aufschreiben der Gleichgewichtsbedingungen für die Knotenschnitte benötigt man insgesamt drei Winkel aus (nebenstehende Skizze). Die für die Berechnung erforderlichen Winkelfunktionen werden vor dem Aufstellen des Gleichungssystems im Matlab-Script bereitgestellt.
Das nachfolgend gelistete Matlab-Script (Download als FW4Gls.m möglich) baut das Gleichungssystem aus den 14 Gleichgewichtsbedingungen auf und löst es:
clear all
alpha1 = atan(1/1) ; ca1 = cos(alpha1) ; sa1 = sin(alpha1) ;
alpha2 = atan(1/2) ; ca2 = cos(alpha2) ; sa2 = sin(alpha2) ;
alpha3 = atan(3/1) ; ca3 = cos(alpha3) ; sa3 = sin(alpha3) ;
% FS1 FS2 FS3 FS4 FS5 FS6 FS7 FS8 FS9 FS10 FS11 FS12 FS13 FS14
A = [ 0 0 0 1 ca1 0 0 0 0 0 0 0 0 0 ; %4x
1 0 0 0 -sa1 0 0 -1 0 0 0 0 0 0 ; %4y
0 0 0 -1 0 -ca1 0 0 0 0 0 0 0 -ca3 ; %5x
0 0 0 0 0 -sa1 0 0 0 0 0 0 0 -sa3 ; %5y
0 0 0 0 -ca1 0 1 0 0 0 0 0 0 0 ; %6x
0 0 0 0 sa1 0 0 0 0 -1 0 0 0 0 ; %6y
0 0 0 0 0 ca1 -1 0 0 0 0 0 0 0 ; %7x
0 0 0 0 0 sa1 0 0 0 0 -1 0 0 0 ; %7y
0 -1 0 0 0 0 0 0 1 0 0 0 ca2 0 ; %8x
0 0 0 0 0 0 0 1 0 0 0 0 -sa2 0 ; %8y
0 0 0 0 0 0 0 0 -1 0 0 ca1 0 0 ; %9x
0 0 0 0 0 0 0 0 0 1 0 -sa1 0 0 ; %9y
0 0 0 0 0 0 0 0 0 0 0 -ca1 -ca2 ca3 ; %10x
0 0 -1 0 0 0 0 0 0 0 1 sa1 sa2 sa3 ] ; %10y
b = [0 ; 0 ; 0 ; 1 ; 0 ; 0 ; 0 ; 0 ; 0 ; 0 ; 0 ; 0 ; 0 ; 0 ] ; % Rechte Seite
FSi = A \ b ;
for i=1:14
disp(['FS', num2str(i) , ' = ' , num2str(FSi(i))]) ;
end
Nebenstehend sind die berechneten Ergebnisse zu sehen, die die oben beschriebenen Rechnungen bestätigen.
