地震波位移、速度和加速度反应谱绘制.doc

想预览更多内容,点击预览全文

申明敬告:

本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。如果您已付费下载过本站文档,您可以点击这里二次下载

文档介绍

地震波位移、速度和加速度反应谱绘制

问题:绘制EI-Centro-NS地震波位移、速度和加速度的反应谱

答:读取EI-Centro-NS地震波的加速度数据,使用Matlab进行编程,利用Matlab内部函数lsim求解,绘制出4种不同阻尼比的情况下的反应谱图,4种阻尼比分别为0、0.02、0.05、1。

Matlab程序如下:

%%%%%%画出加速度时程曲线%%%%%%

clear

format compact

t=xlsread('C:\Users\shaoqun123\Desktop\el-centro地震波南北向数据.xlsx','Sheet1','A1:A2675');

data=xlsread('C:\Users\shaoqun123\Desktop\el-centro地震波南北向数据.xlsx','Sheet1','B1:B2675');

%% 计算反应谱

kesi=[ 0 0.02 0.05 0.1];

T=0.1:0.1:5;

Sd=zeros(4,length(T));

Sv=zeros(4,length(T));

Sa=zeros(4,length(T));

for i=1:length(kesi) for j=1:length(T) w=2*pi./T(j); A=[0 1;-w^2 -2*kesi(i)*w]; B=[0;1]; C=[1 0; 0 1]; D=[0 ;0]; [y,Z]=lsim(A,B,C,D,data,t); d=Z(:,1); v=Z(:,2); a=-2*kesi(i)*w*v-w^2*d; Sd (i,j)=max(abs(d)); Sv(i,j)=max(abs(v)); Sa(i,j)=max(abs

您可能关注的文档

最近下载