% FEMSET-Berechnung: Ebenes Fachwerk mit 10 Staeben und 7 Knoten einschliesslich
% graphischer Darstellung und Berechnung der Stabkraefte

clear all

%       1       2        3           4          5        6             7        <-- Knotennummern
xy = [0 1000 ; 0 0 ; 1000 1000 ; 1000 500 ; 1000 0 ; 2000 1000  ;  2000 500] ;  % Knotenkoordinaten
kr = [1   1  ; 1 1 ;   0    0  ;   0   0  ;   0  0 ;   0    0   ;    0   0 ] ;  % Geom. Randbednggn.
bk = [0   0  ; 0 0 ;   0    0  ;   0   0  ;   0  0 ;   0 -10000 ; 10000  0 ] ;  % Knotenlasten

%      1     2     3     4     5     6     7     8     9     10                 <-- Elementnummern
km = [6 7 ; 3 6 ; 3 7 ; 5 7 ; 3 4 ; 4 5 ; 1 3 ; 1 4 ; 2 4 ; 2 5] ;              % Topologie
ep(1:10,1) = 210000*300 ;                                                       % Dehnsteifigkeiten

clf ;
drawfw2d (xy , km , kr , bk , 'ek') ;                  % 2D-Fachwerk zeichnen

% Finite-Elemente-Algorithmus:
[succ uv] = femalg_m (xy , km , ep , kr , bk) 

if succ == 1                                           % Nach erfolgreicher Rechnung: 
    drawfw2duv (xy , km , kr , bk , uv) ;              % Verformtes System zeichnen
    sk = Stabkraefte2D (xy , km , ep , uv)             % Berechnung der Stabkraefte ...
    figure ;                                           % ... und zeichnen des Systems in ..
    drawfw2dsk (xy , km , kr , bk , sk , 'e') ;        % ... in ein neues Graphik-Fenster
end