地震波位移、速度和加速度反应谱绘制.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

最近下载