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 - 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 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