% Aufgabe 18-15, Berechnungsmodell (2 Elemente, 4 Knoten):
clear all
xy = [0 ; 1200 ; 1200 ; 2500] ;  % Knotenkoordinaten (Doppelknoten am Gelenk)
km = [1 2 ; 3 4] ;               % Koinzidenzmatrix
ep = [3e9 0  0   ;               % 1. Spalte: Biegesteifigkeit,
      5e9 0 -1.4 ] ;             % 2. und 3. Spalte: Linienlastintensitaet
kr = [ 0    0    ; 0 0 ;   0  0  ; 1  0 ] ;  % Verhinderte Verformungen
bk = [ 0 -220000 ; 0 0 ; -800 0  ; 0  0 ] ;  % Einzelkraefte und -momente
sc = [50    0    ; 0 0 ;  300 0  ; 0 5e6] ;  % Federn, Drehfeder
li = [ 0    0    ; 1 0 ;   1  0  ; 0  0 ] ;  % Gleichheit von Verformungen 
 
% Finite-Elemente-Algorithmus (Berechnung der Knotenverformungen):
[succ uv] = femalg_m (xy , km , ep , kr , bk , sc , li) 