clear all

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%	FEM 
%	Problem 3.7
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%---------- Geometry ------------


% Läs in geometrin
[Edof,Ex,Ey,nelm,ndof,bc]=geom37;
figure(1)
% Rita ut nätet 'meshet'
eldraw2(Ex,Ey,[1 2 1],1:nelm);
hold on

%------ Material properties ----

% Definera material/struktur parametrar
E=2.1*10^11;
A_HEA=2124*10^(-6);
A_RHS=1530*10^(-6);


A=A_RHS; % Change to A_RHS for second solution


ep=[E A];

% Notera att vi först måste skapa en matris K
% med rätt storlek. Fyll med nollor
K=zeros(ndof);


% För varje element: Ta fram elementstyvhetsmatrisen
% genom bar2e. Stoppa in i styvhetsmatrisen med hjälp 
% av assem	
for elnr=1:nelm
   Ke=bar2e(Ex(elnr,:),Ey(elnr,:),ep);
   K=assem(Edof(elnr,:),K,Ke);
end

% Definera lastvektor med 'rätt storlek'	
f=zeros(ndof,1);

% Kraften läggs på i frihetsgrad 10
f(10)=-20*10^3;


% Lös systemet med givna randvillkor bc	
[a,Q]=solveq(K,f,bc);


% För att kunna räkna ut t ex behöver vi 
% förskjutningarna på elementnivå Ed

Ed=extract(Edof,a);

	
displ=[a(9) a(10)];


% Detta kommando räknar ut den deformerade strukturen.
% För att deformationen skall bli tydligare förstorar 
% vi deformationen en faktor 30	
eldisp2(Ex,Ey,Ed,[1 4 0],30);



% Vi kan beräkna (från Hooke's) lag beräkna 
% spänningarna i strukturen
for elnr=1:nelm
   es(elnr)=bar2s(Ex(elnr,:),Ey(elnr,:),ep,Ed(elnr,:));
end


