Algoritmo Matlab assemblaggio Matrici K M

%% Dati sistema - aste assialmente indeformabili - Lumped Mass Matrix + Condensazione Statica - analisi ad un solo telaio

%Estendibile a Simmetric Plan Building(al posto di Kstruct, metto la somma delle matrici laterali dei telai)

%necessaria condensazione statica se sono inclusi i gradi di libertà rotazionali per Kstruct

%massa [m(x)=m]
m=50000; %Kg

%lunghezza
L=3.50; %m

%modulo - sezione
E=2060000000; %KPa
% b=0.3; %m
% h=0.6; %m
I= 0.0006;%A=b*h/12;
zita=0.05;

Nodi

%nodi - o X Y sistema di riferimento assoluto su un incastro
nodo(1).x=0;
nodo(1).y=0;
nodo(1).vinc=1;
nodo(2).x=2*L;
nodo(2).y=0;
nodo(2).vinc=1;
nodo(3).x=0;
nodo(3).y=L;
nodo(4).x=2*L;
nodo(4).y=L;
nodo(5).x=0;
nodo(5).y=2*L;
nodo(6).x=2*L;
nodo(6).y=2*L;
nodo(7).x=0;
nodo(7).y=3*L;
nodo(8).x=2*L;
nodo(8).y=3*L;

Aste

%aste 0 - x(asse trave) - y sistema locale
asta(1).nodoi=1;
asta(1).nodoj=3;
asta(1).massa=0;
asta(1).L=L;
asta(1).E=E;
asta(1).I=I;
%asta(1).A=A;
asta(1).dof=[0 0 1 4]; % ordinati mettendo a sx nodo i e destra nodo j,zero equivale a nodo fermo; il nodo j ha gli stessi spost. del nodo del telaio e quindi assumerà gli stessi dof_telaio
asta(1).trasl=0;

asta(2).nodoi=4;
asta(2).nodoj=2;
asta(2).massa=0;
asta(2).L=L;
asta(2).E=E;
asta(2).I=I;
%asta(2).A=A;
asta(2).dof=[1 5 0 0]; % i fa solo gdl 2 del telaio, j solo 3
asta(2).trasl=0;

asta(3).nodoi=3;
asta(3).nodoj=4;
asta(3).massa=m;
asta(3).L=2*L;
asta(3).E=0.5*E;
asta(3).I=I;
%asta(3).A=A;
asta(3).dof=[0 4 0 5]; %i fa gdl del telaio 1 e 3
asta(3).trasl=1; %perchè trasla come il telaio, ovvero segue il DOF 1 del telaio

asta(4).nodoi=3;
asta(4).nodoj=5;
asta(4).massa=0;
asta(4).L=L;
asta(4).E=0.666*E;
asta(4).I=I;
%asta(4).A=A;
asta(4).dof=[1 4 2 6]; %i fa gdl del telaio 1 e 3
asta(4).trasl=0;

asta(5).nodoi=6;
asta(5).nodoj=4;
asta(5).massa=0;
asta(5).L=L;
asta(5).E=0.666*E;
asta(5).I=I;
%asta(5).A=A;
asta(5).dof=[2 7 1 5]; %i fa gdl del telaio 1 e 3
asta(5).trasl=0;

asta(6).nodoi=5;
asta(6).nodoj=6;
asta(6).massa=m;
asta(6).L=2*L;
asta(6).E=0.333*E;
asta(6).I=I;
%asta(6).A=A;
asta(6).dof=[0 6 0 7]; %i fa gdl del telaio 1 e 3
asta(6).trasl=2;

asta(7).nodoi=5;
asta(7).nodoj=7;
asta(7).massa=0;
asta(7).L=L;
asta(7).E=0.333*E;
asta(7).I=I;
%asta(7).A=A;
asta(7).dof=[2 6 3 8]; %i fa gdl del telaio 1 e 3
asta(7).trasl=0;

asta(8).nodoi=8;
asta(8).nodoj=6;
asta(8).massa=0;
asta(8).L=L;
a