作业 1绘制 1940 El Centro,NS 重量地震动得绝对加速度、相对速度与相对位移反应谱。地震动:在 PEER Ground Motion Database 自行下载经典得 1940 El Centro,NS 重量。 要求:在此模板内完成,A4 纸打印。自编程序与软件(Bispec 或 Seismosigna 等)计算反应谱进行对比。提交自编写程序。Matlab 程序:clearfid = fopen('E:\Earthquake\El centro、txt');[Accelerate,count] = fscanf(fid,'%g'); %count 读入得记录得量Accelerate=9、8*Accelerate'; %单位统一为 m 与 stime=0:0、02:(count1)*0、02; %单位 sDisplace=zeros(1,count); %相对位移Velocity=zeros(1,count); %相对速度AbsAcce=zeros(1,count); %绝对加速度DampA=[0、00,0、02,0、05]; %三个阻尼比TA=0、0:0、02:4; %TA=0、000001:0、02:4;结构周期Dt=0、02; %地震记录得步长%记录计算得到得反应,MDis 为某阻尼时最大相对位移%MVel 为某阻尼时最大相对速度,MAcc 某阻尼时最大绝对加速度MDis=zeros(3,length(TA));MVel=zeros(3,length(TA));MAcc=zeros(3,length(TA));j=1; for Damp=[0、00,0、05,0、1] t=1; for T=0、0:0、02:4 Frcy=2*pi/T ; DamFrcy=Frcy*sqrt(1Damp*Damp); e_t=exp(Damp*Frcy*Dt); s=sin(DamFrcy*Dt); c=cos(DamFrcy*Dt); A=zeros(2,2); A(1,1)=e_t*(s*Damp/sqrt(1Damp*Damp)+c); A(1,2)=e_t*s/DamFrcy; A(2,1)=Frcy*e_t*s/sqrt(1Damp*Damp); A(2,2)=e_t*(s*Damp/sqrt(1Damp*Damp)+c); d_f=(2*Damp^21)/(Frcy^2*Dt); d_3t=Damp/(Frcy^3*Dt); B=zeros(2,2); B(1,1)=e_t*((d_f+Damp/Frcy)*s/DamFrcy+(2*d_3t+1/Frcy^2)*c)2*d_3t; B(1,2)=e_t*(d_f*s/DamFrcy+2*d_3t*c)1/Frcy^2+2*d_3t; B(2,1)=e_t*((d_f+Damp/Frcy)*(cDamp/sqrt(1Damp^2)*s)(2*d_3t+1/Frcy^2)*(DamFrcy*s+Damp*Frcy*c))+1/(Frcy^2*Dt); B(2,2)=e_t*(1/(Frcy^2*Dt)*c+s*Damp/(Frcy*DamFrcy*Dt))1/(Frcy^2*Dt); for i=1:(count1) Displace(i+1)=A(1,1)*Displace(i)+A(1,2)*Velocity(i)+B(1,1)*Accelerate(i)+B(1,2)*Accelerate(i+1); Velocity(i+1)=A(2,1)*Displace(i)+A(2,2)*Velocity(i)+B(2,1)*Accelerate(i)+B(2,2)*Accelerate(i+1); AbsAcce(i+1)=2*Damp*Frcy*Velocity(i+1)Frcy^2*Displace(i+1); end MDis(j,t)=max(abs(Displace)); MVel(j,t)=max(abs(Velocity)); if T==0、0 MAcc(j,t)=max(abs(Accelerate));else MAc...
