Algoritmo Matlab spettri Risposta metodoNewmark

clear all

close all

clc

%% COSTRUZIONE DELLO SPETTRO DI PROGETTO ELASTICO - METODO DI NEWMARK E HALL

%84.1% probabilità di non superamento

% Parametri del suolo

g=9.81; %m/s^2

ag0=0.26*g; %accelerazione massima al suolo, misurata in

g perchè nell'accelerogramma è in g

vg0=1.2192*ag0/(g); %perchè si misura in cm/s * 1/g (g da trasformare in

cm/sec^2)

ug0=6*(vg0^2)/(ag0); %cosi ottengo lo spostamento del suolo in cm

perchè g è in m/s*s *100= cm/s*s

% Coefficienti di amplificazione della pseudoaccelerazione - pseudovelocità e

spostamento

zita=5; %espresso in %

alfa_a= 4.38-1.04*log(zita);

alfa_v= 3.38-0.67*log(zita);

alfa_d= 2.73-0.45*log(zita);

% Valori amplificati

Aa=alfa_a*ag0;

Va=alfa_v*vg0;

Da=alfa_d*ug0;

% Periodi dello spettro - dal diagramma tetralogaritmico - in secondi

Ta=1/33;

Tb=1/8;

Tc=(2*pi*Va)/(Aa); % dall'intersezione

Td=(2*pi*Da)/(Va);

Te=10;

Tf=33;

%% Spettro della pseudoaccelerazione spettrale di progetto --- scala

logaritmica

% %dati ricavati dal diagramma tetralogaritmico

%

% %periodi Tn<Ta

% Tn1=0:Ta/500:Ta;

% for i=1:length (Tn1)

% A1(i)=ag0; %per periodi inferiori a Ta non amplifica -

struttura rigida

% end

%

% %periodi Ta<Tn<Tb

% Tn2_=Ta:((Tb-Ta)/500):Tb;

% ii=1;

% for Tn2=Ta:(Tb-Ta)/500:Tb

% A2(ii)= Tn2*(Aa-ag0)/(Tb-Ta) - ((Ta*(Aa-ag0)/(Tb-Ta))-ag0); %retta tra

Tb e Tc - tratto di transizione

% ii=ii+1;

% end

%

% %periodi Tb<Tn<Tc

% Tn3=Tb:((Tc-Tb)/500):Tc;

% for j=1:length(Tn3)

% A3(j)=Aa; %tratto ad accelerazione costante

% end

%

% %periodi Tc<Tn<Td

% Tn4_=Tc:(Td-Tc)/500:Td;

% jj=1;

% for Tn4=Tc:((Td-Tc)/500):Td

% A4(jj)=Tn4*((2*pi/Td)*Va-Aa)/(Td-Tc) - ((Tc*((2*pi/Td)*Va-Aa)/(Td-Tc))-

Aa); %tratto a velocità costante

% jj=jj+1;

% end

%

% %periodi Td<Tn<Te

% Tn5_=Td:(Te-Td)/500:Te;

% k=1;

% for Tn5=Td:(Te-Td)/500:Te

% A5(k)=Tn5*((4*pi*pi/(Te^2))*Da - (2*pi/Td)*Va)/(Te-Td)- (Td*(((4*pi*pi/

(Te^2))*Da - (2*pi/Td)*Va)/(Te-Td)) - (2*pi/Td)*Va); %tratto a spostamento

amplificato

% k=k+1;

% end

%

% %periodi Te<Tn<Tf

% Tn6_=Te:(Tf-Te)/50:Tf;

% kk=1;

% for Tn6=Te:(Tf-Te)/50:Tf

% A6(kk)= Tn6*((((4*pi*pi/(Tf^2))*ug0 - (4*pi*pi/(Te^2))*Da))/(Tf-Te)) -

(Te*(((4*pi*pi/(Tf^2))*ug0 - (4*pi*pi/(Te^2))*Da)/(Tf-Te)) - (4*pi*pi/

(Te^2))*Da); %tratto a spostamento no amplificato

% kk=kk+1;

% end

%

% figure(2)

% plot (Tn1,A1/(g))

% hold on

% plot (Tn2_,A2/(g))

% plot (Tn3,A3/(g))

% plot (Tn4_,A4/(g))

% plo