% Aufgabe 19-8, Matlab-Femset: Polygon aus biege- und dehnsteifen geraden Traegern
clear all
R = 52 ; E = 210000 ; I = 5*20^3/12 ; A = 5*20 ; F = 1000 ;

n  = 3 ;                                 % n Elemente ---> n+1 Knoten
xy = zeros (n+1,2) ; km = zeros (n  ,2) ;
ep = zeros (n  ,4) ; kr = zeros (n+1,3) ; bk = zeros (n+1,3) ;

dphi = (pi/2)/n ; phi  = 0 ;
for i=1:n+1
    xy(i,:) = [R*cos(phi) , R*sin(phi)] ; % Knotenkoordinaten
    phi = phi + dphi ;
end

for i=1:n
    km(i,:) = [i , i+1] ;                 % Koinzidenzmatrix (Knoten der Elemete)
    ep(i,:) = [E*I , E*A , 0 , 0] ;       % Elementparameter
end

kr(1,:)   = [0   1 1] ;                   % Verhinderte Verschiebungen unten (Knoten 1)
kr(n+1,:) = [1   0 1] ;                   % Verhinderte Verschiebungen links (Knoten n+1)
bk(1,:)   = [F/2 0 0] ;                   % Belastung am Knoten 1
 
% Finite-Elemente-Algorithmus:
[succ uv] = femalg_m (xy , km , ep , kr , bk)