Algoritmo Matlab assemblaggio Matrici K M portale1piano

Dati sistema - aste assialmente indeformabili

close allclear allclc

Massa

%massa [m(x)=m]m=100; %Kg/m, massa fondamentale, dato del problema

Lunghezza

%lunghezzaL=2; %lunghezza fondamentale

Modulo e sezione

%modulo - sezioneE=33000000; %KPa, per C32/40b=0.3; %mh=0.6; %mI= b*(h^3)/12;A=b*h;

Nodi

%nodi - o X Y sistema di riferimento assoluto su un incastronodo(1).x=0;nodo(1).y=0;nodo(2).x=0;nodo(2).y=L;nodo(3).x=2*L;nodo(3).y=L;nodo(4).x=2*L;nodo(4).y=0;

Aste

%aste 0 - x(asse trave) - y sistema localeasta(1).nodoi=1;asta(1).nodoj=2;asta(1).massa=m;asta(1).L=L;asta(1).E=E;asta(1).I=I;asta(1).A=A;asta(1).dof=[0 0 1 2]; % 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(2).nodoi=2;asta(2).nodoj=3;asta(2).massa=1.5*m;asta(2).L=2*L;asta(2).E=E;asta(2).I=I;asta(2).A=A;asta(2).dof=[0 2 0 3]; % i fa solo gdl 2 del telaio, j solo 3

asta(3).nodoi=3;asta(3).nodoj=4;asta(3).massa=m;asta(3).L=L;asta(3).E=E;asta(3).I=I;asta(3).A=A;asta(3).dof=[1 3 0 0]; %i fa gdl del telaio 1 e 3

Dati telaio

%dati telaion_nodi=length(nodo);