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

第頁(yè)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: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/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)>0f(4)=0;elsef(4)=9000;endenda=max(abs(F(:,13)));b=max(abs(F(:,14)));Fabxlswrite('表格2.xls',F)計(jì)算書：課程設(shè)計(jì)計(jì)算書（題二）根據(jù)加速度調(diào)幅公式：得：=[060011001500210025002900350205015001000600200-700-1300-1700-2000-1800-1500-700-250200-100000];所以經(jīng)調(diào)幅后為=[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]依次類推可以求出地面運(yùn)動(dòng)加速度的差值。因?yàn)?,代入可以算得一、表格第一行數(shù)據(jù)計(jì)算：,t=0.05s代入得起始時(shí)刻時(shí)：因?yàn)樗运韵鄬?duì)加速度=層間剪力=二、表格第二行數(shù)據(jù)計(jì)算：=-250[379.7-6*(-10.8)/0.05-3*(-432)]+240*[3*(-10.8)+(-432)*0.05/2]=-753293N因?yàn)樗?-312=-432*0.05-312*0.05/2=-29.4相對(duì)加速度==455.2-432=23.2層間剪力==250*23.2=5800N依次類推，每次以上一行的數(shù)據(jù)作為起始值來(lái)計(jì)算下一行的數(shù)值?？梢缘贸鏊袛?shù)據(jù)。考慮到此為假定的結(jié)構(gòu)彈塑性恢