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: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);if abs(f(10)2&F(1,7)*f(7)0 f(4)=0; else f(4)=9000; endenda=max(abs(F(:,13);b=max(abs(F(:,14);Fabxlswrite(表格2.xls,F)計算書：課程設(shè)計計算書（題二）根據(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 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 依次類推可以求出地面運動加速度的差值。因為 , , 代入可以算得一、表格第一行數(shù)據(jù)計算： , t=0.05s 代入得 起始時刻時： 因為所以所以相對加速度=層間剪力=二、表格第二行數(shù)據(jù)計算：=-250379.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依次類推，每次以上一行的數(shù)據(jù)作為起始值來計算下一行的數(shù)值。可以得出所有數(shù)據(jù)。考慮到此為假定的結(jié)構(gòu)彈塑性恢復(fù)力模型，屈服后結(jié)構(gòu)抗側(cè)移剛度為0，卸載階