MATLAB弹性时程分析法编程
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
计算书：课程设计计算书（题一）
根据加速度调幅公式：m i a t a a a /)(max ,00*=
)/(29002902s mm Gal a m ==
得：58/)(72900/)(3500i i t a t a a =*=
)(i t a =[0 600 1100 150021002500 2900350 2050
15001000600200 -700 -1300-1700 -2000 -1800-1500 -700-250200 -100 0 0 0];
所以经调幅后为0a =[0 72.6 133.1 181.5 254.1 302.5 350.9 42.4 248.1 181.5 121 72.6 24.2 -84.7 -157.3 -205.7 -242 -217.8 -181.5 -84.7 -30.3 24.2-12.1 0 0 0 ]
6.7206.72''1''2=-=-U U
5.60
6.721.133''2''3=-=-U U
依次类推可以求出地面运动加速度的差值。 因为km c 2=ζ,08.0=ζ ,
m kN k /9000=, m s kN m /2502⋅=
代入可以算得m s kN c /240⋅=
一、表格第一行数据计算：
t c t m k K i i /3/62++=*, t=0.05s
代入得m N K i /623400
=* ）△△2
/3()3/6(''''''''t U U c U t U U m P i i g i *++---=* N 18150-6.72250-=*=
**=i i P U K △△
mm K P U i i 03.0623400/18150
/-=-==**△△ 起始时刻时：0=U 0'=U 0''=U
因为'''2''3/6/6i i U t U t U U -*-*=△△
所以7205.0/)03.0(62''1
-=-*=U △
2/'''''t U t U U i *+*=△△
所以8.12/05.072'1-=*-=U △
相对加速度=)/(0002''''s mm U U g =-=-
层间剪力=N U U m g
0)(''''=- 二、表格第二行数据计算：
03.012-==U U △
72''1''2-==U U △
8.1'1'2-==U U △
)/(5.606.721.1332''s mm U g =-=△
)2/3()3/6(''''''''t U U c U t U U m P i i i i g i *++---=*△△
=-250[60.5-6*(-1.8)/0.05-3*(-72)]+240*[3*(-1.8)+(-72)*0.05/2] =-124853N
mm K P U i i 2.0623400/124853
/2-=-==**△△ 因为'''2''3/6/6i i U t U t U U -*-*=△△
所以)72(305.0/)8.1(605.0/620.02''2
-*--*-*-=U =-48)/(2s mm
2/'''''2t U t U U i *+*=△△
=-72*0.05-48*0.05/2=-4.8)/(s mm
相对加速度='
'''U U g +
=72.6-72=0.6)/(2s mm
层间剪力=)(''''U U m g +
=250*0.6=150N
依次类推，每次以上一行的数据作为起始值来计算下一行的数值。可以得出所有数据。
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*350/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);
end
a=max(abs(F(:,13)));
b=max(abs(F(:,14)));
F
a
b
xlswrite('表格1.xls',F)
输出结果：