MATLAB弹性时程分析法编程

《MATLAB弹性时程分析法编程》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、计算书：课程设计计算书（题一）根据加速度调幅公式：得：=[060011001500210025002900350205015001000600200-700-1300-1700-2000-1800-1500-700-250200-100000];所以经调幅后为=[072.6133.1181.5254.1302.5350.942.4248.1181.512172.624.2-84.7-157.3-205.7-242-217.8-181.5-84.7-30.324.2-12.1000]依次类推可以求出地面运动加速度的差值。因为,,,代入可以算得一、表格第

2、一行数据计算：,t=0.05s代入得起始时刻时：因为所以所以相对加速度=层间剪力=二、表格第二行数据计算：=-250[60.5-6*(-1.8)/0.05-3*(-72)]+240*[3*(-1.8)+(-72)*0.05/2]=-124853N因为所以=-48=-72*0.05-48*0.05/2=-4.8相对加速度==72.6-72=0.6层间剪力==250*0.6=150N依次类推，每次以上一行的数据作为起始值来计算下一行的数值。可以得出所有数据。MATLAB编程“formatshortg;F(21,14)=0;xg=[060011001500

3、210025002900350205015001000600200-700-1300-1700-2000-1800-1500-700-250200-100000];xg1=xg*350/max(xg);xg2=diff(xg1);f(14)=0;f(4)=9000;t=0.05;m=250;c=240;fori=1:21f(3)=xg2(i);f(1)=0.05*i-0.05;f(2)=xg1(i);f(6)=-m*(f(3)-6*f(12)/t-3*f(11))+c*(3*f(12)+f(11)/2*t);f(5)=f(4)+6*m/t^2+3*c

4、/t;f(7)=f(6)/f(5);f(9)=3/t*f(7)-3*f(12)-0.5*f(11)*t;f(8)=6/0.05^2*f(7)-6/0.05*f(12)-3*f(11);f(13)=f(2)+f(11);f(14)=f(13)*m;F(i,:)=f(1,:);f(10)=f(10)+f(7);f(12)=f(12)+f(9);f(11)=f(11)+f(8);enda=max(abs(F(:,13)));b=max(abs(F(:,14)));Fabxlswrite('表格1.xls',F)输出结果：