地震工程作业

作业 1绘制 1940 El Centro，N-S 重量地震动的绝对加速度、相对速度和相对位移反应谱.地震动：在 PEER Ground Motion Database 自行下载经典的 1940 El Centro，N—S 重量。 要求：在此模板内完成，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：(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，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（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(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)*(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）; for i=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）＊Displ...