预览加载中,请您耐心等待几秒...
1/7
2/7
3/7
4/7
5/7
6/7
7/7

在线预览结束,喜欢就下载吧,查找使用更方便

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

MATLAB编程: formatshortg; F(21,14)=0; xg=[060011001500210025002900350205015001000600200-700-1300-1700-2000-1800-1500-700-250200-100000]; xg1=xg*2200/max(xg); xg2=diff(xg1); f(14)=0;f(4)=9000;t=0.05;m=250;c=240; fori=1:21 f(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/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); ifabs(f(10))>2&F(1,7)*f(7)>0 f(4)=0; elsef(4)=9000; end end a=max(abs(F(:,13))); b=max(abs(F(:,14))); F a b xlswrite('表格2.xls',F) 计算书:课程设计计算书(题二) 根据加速度调幅公式: 得: =[060011001500210025002900350205015001000600200-700-1300-1700-2000-1800-1500-700-250200-100000]; 所以经调幅后为=[0455.2834.91138.51593.91897.52201.1265.71556.01138.5759455.4151.8-531.3-986.7-1290.3-1518-1366.2-1138.5-531.3-189.8151.8-75.9000] 依次类推可以求出地面运动加速度的差值。 因为,, , 代入可以算得 一、表格第一行数据计算: ,t=0.05s 代入得 起始时刻时: 因为 所以 所以 相对加速度= 层间剪力= 二、表格第二行数据计算: =-250[379.7-6*(-10.8)/0.05-3*(-432)]+240*[3*(-10.8)+(-432)*0.05/2] =-753293N 因为 所以 =-312 =-432*0.05-312*0.05/2=-29.4 相对加速度= =455.2-432=23.2 层间剪力= =250*23.2=5800N 依次类推,每次以上一行的数据作为起始值来计算下一行的数值。可以得出所有数据。 考虑到此为假定的结构弹塑性恢复力模型,屈服后结构抗侧移刚度为0,卸载阶段又恢复为8000N/m,所以计算时当绝对值△U大于2mm时,要以ki=0N/m代入计算即可。 输出结果: