algoritmo Matlab costruzione spettri di risposta elastici

Note preliminari per il funzionamento del programma

Per fare funzionare il programma, è necessario caricarlo interamente in Matlab, creando prima un foglio nel blocco note in cui si carica l'accelerogramma per punti (lo si può scaricare dal sito INGV). Dare a questo file il nome dell'accelerogramma, ad esempio in questo caso elcentro. Quindi, nell'algoritmo, sostituire load elcentro.txt con load ______ il nome dato ed ogni altra istruzione contenente elcentro con il nome dato. Nel caso in cui il file scaricato dell'accelerogramma è disordinato, consiglio di ordinarlo per colonne in Excel e sostituire le virgole con punti, quindi lo si copia nel blocco note, già ordinato e poi proseguire come sopra.

Algoritmo da copiare

clear all
close all
clc
%% COSTRUZIONE SPETTRI DI RISPOSTA PER IL TERREMOTO IN ESAME
% risposta per vari Tn allo stesso sisma del sistema SDOF calcolata con
risoluzione equazione differenziale con metodo numerico interp.della
% forzante
% ------------------- Accelerogramma -------------------
% passo dt assunto costante, i coefficienti dipendono solo dalle caratteristiche
del sistema
load elcentro.txt % accelerogramma
elcentro(:,1)=[]; % tolgo la colonna del tempo, qst istruzione la posso
togliere
acc=elcentro;
n=length(acc);
dt=0.02; % passo d'integrazione=passo di misura
t=0:dt:dt*(n-1);
g=9.807; % controllare se l'accelerogramma fornito è in
m/s2 in caso cambiare
figure(1)
plot(t,acc/g)
title('Accelerogramma El Centro, comp. N-S');
xlabel('t (s)')
ylabel('a(t), g')
grid
zita=0.05;
j=1;
p=-acc; % perché l'accelerazione totale, solo la
parte dovuta al moto del suolo
for Tn=0.01:0.01:max(t)
wn=2*pi/Tn;
k=(wn^2); % con m=1, k è uguale alla
frequenza circolare
wd=wn*sqrt(1-(zita^2));
A=(exp(-zita*wn*dt))*((zita/sqrt(1-zita*zita))*sin (wd*dt) +
cos(wd*dt));