matlab弹塑性时程分析法编程

MATLAB 编程：at 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/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);if abs(f(10))>2else f(4)=9000;endenda=max(abs(F(:,13)));b=max(abs(F(:,14)));Fabxlswrite( 表格 2.xls ,F)计算书：课程设计计算书(题二)根据加速度调幅公式： miata/)(mx,0)/(2902sGl得： 9/(/)0 ii tataa=[0 600 1100 1500 2100 2500 2900 350 )(ita2050 1500 1000 600 200 -700 -1300 -1700 -2000 -1800 -1500 -700 -250 200 -100 0 0 0];所以经调幅后为 =[0 455.2 834.9 1138.5 1593.9 0a1897.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 ]2.4502.45 1 2U7.39.9.83 2 3依次类推可以求出地面运动加速度的差值。因为 , ,kmc208., Nk/90mskN/252代入可以算得 c/40一、表格第一行数据计算：, t=0.05s tctmkKii /3/62代入得 Ni /40 )△△ 2/3()/6( tUctUP iigi 1380-2.450-iiPUK△△ mii 18.06234/180/ △△起始时刻时： 0U 0 U因为 2 3/6/6iitt△△所以 4205./)18.(2 1 △ / tUtUi △△所以 8.102/05.43 1 △相对加速度= )/(2 smg层间剪力= NUmg0)( 二、表格第二行数据计算： 18.012U△ 432 1 2△ 8.0 1 2△ )/(7.392.459.32 smUg△ )2/3(/6( tUctUmPiiiigi △△=-250[379.7-6*(-10.8)/0.05-3*(-432)]+240*[3*(-10.8)+(-432)* 0.05/2]=-753293N mKPUii 21.62340/759/2 △△因为 2 3/6/6iiUttU△△所以 )432(05./)8.1(605./621.2 2 △=-312 )/(sm2/ 2tUti △△=-432*0.05-312*0.05/2=-29.4 )/(sm相对加速度= g=455.2-432=23.2 )/(2s层间剪力= )( Umg=250*23.2=5800N依次类推，每次以上一行的数据作为起始值来计算下一行的数值。可以得出所有数据。考虑到此为假定的结构弹塑性恢复力模型，屈服后结构抗侧移刚度为 0，卸载阶段又恢复为 8000N/m,所以计算时当绝对值△U 大于 2mm 时，要以 ki=0N/m 代入计算即可。输出结果：