计算书：课程设计计算书（题一）

根据加速度调幅公式：a0?a0,max?a(ti)/am

am?290Gal?2900(mm/s2)

得：a0?350?a(ti)/2900?7a(ti)/58

a(ti)=[0

600 1100 150021002500 2900350 2050

15001000600200 -700 -1300-1700 -2000 -1800-1500

-700-250200 -100 0 0 0];

所以经调幅后为a0=[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 ]

''U2?U1''?72.6?0?72.6

''''U3?U2?133.1?72.6?60.5

依次类推可以求出地面运动加速度的差值。

c因为??2km,??0.08 ,

k?9000kN/m, m?250kN?s2/m

代入可以算得c?240kN?s/m 一、表格第一行数据计算：

Ki??ki?6m/t2?3c/t, t=0.05s

?N/m 代入得Ki?623400?''''''''△P??m(△U?6U/t?3U)?c(3U?U） igii?t/2?-250?72.6?-18150N

?Ki?△U?△Pi

??△U?△P/K/623400??0.03mm ii??18150'''起始时刻时：U?0U?0U?0

''2'''△U?6?△U/t?6?U/t?3U因为ii

''2△U?6?(?0.03)/0.05??72 所以1△U'?Ui''?t?△U''?t/2

'△U所以1??72?0.05/2??1.8 ''''2相对加速度=Ug?U?0?0?0(mm/s)

''''m(U?U)?0N 层间剪力=g二、表格第二行数据计算：

U2?△U1??0.03

''U2?△U1''??72

'U2?△U1'??1.8

''△Ug?133.1?72.6?60.5(mm/s2)

?''''''''△P??m(△U?6U/t?3U)?c(3U?Uigiiii?t/2)

=-250[60.5-6*(-1.8)/0.05-3*(-72)]+240*[3*(-1.8)+(-72)*0.05/2] =-124853N

??△U2?△P/K/623400??0.2mm ii??124853''2'''因为△U?6?△U/t?6?Ui/t?3Ui

''2U??0.20?6/0.05?6?(?1.8)/0.05?3?(?72) 所以22=-48(mm/s)

'△U2?Ui''?t?△U''?t/2

=-72*0.05-48*0.05/2=-4.8(mm/s) 相对加速度=Ug?U

2 =72.6-72=0.6(mm/s) ''''层间剪力=m(Ug?U)

'''' =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)

输出结果：