Algoritmo Matlab costruzione spettri NTC08

clear all

close all

clc

%% SPETTRO DI RISPOSTA NTC 08

% Spettro di risposta elastico

%duttilità disponibile

mu=6; %3, 6

% valori del moto al suolo

a_g=0.26*9.81; %m/s^2

F_o=2.36;

T_c_star=0.345; %s

%coefficienti di amplificazione stratigrafica Ss e topografica St

Ss=1;

Cc=1;

St=1;

S=Ss*St;

%Fattore di struttura

K_r=1; %se regolare in pianta, altrimenti 0.8

alfa=2; %valore dell'azione sismica per cui la struttura

diventa labile / valore dell'azione per cui si forma la prima cerniera plastica

q_o=4.5*alfa; %per struttura a telaio, classe di duttilità CDA

% q=K_r*q_o; % SE CHIEDE DI ABBASSARE LO SPETTRO CON MU, PORRE q=mu

q=mu;

zita=5; %smorzamento

eta=sqrt(10/(5+zita));

%Periodi dello spettro

T_c=Cc*T_c_star; %inizio del tratto a velocità costante

T_b=T_c/3; %inizio del tratto ad accelerazione costante

T_d=4*a_g/9.81 + 1.6; %inizio del tratto a spostamento costante

%% Spettro di risposta elastico

i=1;

for T=0:0.001:T_b

T_n(i)=T;

Sa(i)=a_g*S*eta*F_o*(T/T_b + (1/(eta*F_o))*(1- T/T_b));

i=i+1;

end

ii=1;

for T1=T_b:0.001:T_c

T_n1(ii)=T1;

Sa1(ii)=a_g*S*eta*F_o;

ii=ii+1;

end

j=1;

for T2=T_c:0.001:T_d

T_n2(j)=T2;

Sa2(j)=a_g*S*eta*F_o*(T_c/T2);

j=j+1;

end

jj=1;

for T3=T_d:0.001:3

T_n3(jj)=T3;