% Výpočet deformace konstrukce jeřábu zatížené zvedaným břemenem

% uzly
nd = [0 0; 10 0; 0 20; 10 40; ...
      0 60; 10 80; 0 80; -30 80; ...
      30 80; 50 80; 70 80; 50 90; ...
      10 90; 0 90; -20 90];

% elementy
el = [1 2; 1 3; 3 5; 5 7; ...
      2 4; 4 6; 2 3; 3 4; ...
      4 5; 5 6; 8 7; 7 6; ...
      6 9; 9 10; 10 11; 11 12; ...
      12 13; 13 14; 14 15; 15 8; ...
      8 14; 6 14; 13 9; 9 12; ...
      7 14; 6 13; 10 12; 7 13 ];

% materiálové parametry
mat = 5e4 * ones(1,length(el));

% působící síla - uzel 11, Fx = 0, Fy = -10
force = [11 0 -20];

% nulové posuvy na krajích
fix = [1 1; 1 2; 2 2 ];

% vlastní výpočet
u = mkp(el, nd, mat, force, fix)