MATLAB弹塑性时程分析法编程.doc

MATLAB编程：
format short g;
F(21,14)=0;
xg=[0 600 1100 1500 2100 2500 2900 350 2050 1500 1000 600 200 -700 -1300 -1700 -2000 -1800 -1500 -700 -250 200 -100 0 0 0];
xg1=xg*2200/max(xg);
xg2=diff(xg1);
f(14)=0;f(4)=9000;t=0.05;m=250;c=240;
for i =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);
if abs(f(10))>2&F(1,7)*f(7)>0
f(4)=0;
else f(4)=9000;
end
end
a=max(abs(F(:,13)));
b=max(abs(F(:,14)));
F
a
b
xlswrite('表格2.xls',F)

计算书：课程设计计算书（题二）
根据加速度调幅公式：

得：
=[0 600 1100 1500 2100 2500 2900 350 2050 1500 1000 600 200 -700 -1300 -1700 -2000 -1800 -1500 -700 -250 200 -100 0 0 0];
所以经调幅后为=[0 455.2 834.9 1138.5 1593.9 1897.5 2201.1 265.7 1556.0 1138.5 759 455.4 151.8 -531.3 -986.7 -1290.3 -1518 -1366.2 -1138.5 -531.3 -189.8 151.8 -75.9 0 0 0 ]

依次类推可以求出地面运动加速度的差值。
因为 , ,
,
代入可以算得
一、表格第一行数据计算：
, 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代入计算即可。

输出结果：

专业文档供参考，如有帮助请下载。