MATLAB彈性時(shí)程分析法編程

1、計(jì)算書：課程設(shè)計(jì)計(jì)算書（題一）根據(jù)加速度調(diào)幅公式： 得： =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;所以經(jīng)調(diào)幅后為=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 依次類推可以求出

2、地面運(yùn)動(dòng)加速度的差值。因?yàn)?, , 代入可以算得一、表格第一行數(shù)據(jù)計(jì)算： , t=0.05s 代入得 起始時(shí)刻時(shí)： 因?yàn)樗运韵鄬?duì)加速度=層間剪力=二、表格第二行數(shù)據(jù)計(jì)算：=-25060.5-6*(-1.8)/0.05-3*(-72)+240*3*(-1.8)+(-72)*0.05/2=-124853N因?yàn)樗?=-48 =-72*0.05-48*0.05/2=-4.8相對(duì)加速度= =72.6-72=0.6層間剪力= =250*0.6=150N依次類推，每次以上一行的數(shù)據(jù)作為起始值來(lái)計(jì)算下一行的數(shù)值?？梢缘贸鏊袛?shù)據(jù)。MATLAB編程“format short g;F(21,14)=0;xg

3、=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: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/t2+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.052*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