资源描述：
《地震工程作业.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、作业1绘制1940ElCentro，N-S分量地震动的绝对加速度、相对速度和相对位移反应谱。地震动：在PEERGroundMotionDatabase自行下载经典的1940ElCentro，N-S分量。要求：在此模板内完成，A4纸打印。自编程序与软件（Bispec或Seismosigna等）计算反应谱进行对比。提交自编写程序。Matlab程序：clearfid=fopen('E:EarthquakeElcentro.txt');[Accelerate,count]=fscanf(fid,'%g');%count读入的记录的量Accelerate=9.8*Accel
2、erate';%单位统一为m和stime=0:0.02:(count-1)*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
3、,length(TA));MVel=zeros(3,length(TA));MAcc=zeros(3,length(TA));j=1;forDamp=[0.00,0.05,0.1]t=1;forT=0.0:0.02:4Frcy=2*pi/T;DamFrcy=Frcy*sqrt(1-Damp*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(1-Damp*Damp)+c);A(1,2)=e_t*s/DamFrcy;A(
4、2,1)=-Frcy*e_t*s/sqrt(1-Damp*Damp);A(2,2)=e_t*(-s*Damp/sqrt(1-Damp*Damp)+c);d_f=(2*Damp^2-1)/(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)*
5、(c-Damp/sqrt(1-Damp^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);fori=1:(count-1)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)*Dis
6、place(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);endMDis(j,t)=max(abs(Displace));MVel(j,t)=max(abs(Velocity));ifT==0.0MAcc(j,t)=max(abs(Accelerate));elseMAcc(j,t)=max(abs(AbsAcce));endDisplace=zeros(1,co
7、unt);Velocity=zeros(1,count);AbsAcce=zeros(1,count);t=t+1;endj=j+1;endcloseallfigure%绘制位移反应谱plot(TA,MDis(1,:),'-b',TA,MDis(2,:),'-r',TA,MDis(3,:),':k')title('Displacement')xlabel('Tn(s)')ylabel('Displacement(m)')legend('ζ=0','ζ=0.02','ζ=0.05')gridfigure%绘制速度反应谱plot(TA,MVel(1,:)
