MATLAB弹塑性时程分析法编程

1、 - - r M复加型血町位吸皿m 6抗働務刚度肖;激阶艮師训噺渥辟 爾r度设萨郵I讚桂楝滸豪 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 xg1=xg*2200/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

2、(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/tA2+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

3、(11)+f(8); if abs(f(10)2 else f(4)=9000; end end a=max(abs(F(:,13); b=max(abs(F(:,14); xIswriteC 表格 2.xls,F) 计算书：课程设计计算书 （题二） 根据加速度调幅公式: a。 a0,maxa(ti ) / am am 290Gal 2900(mm/s2) a(ti) =06001100 1500 2100 25002900 350 205015001000600 200 - -700 -1300-1700 -2000 -1800-1500 -700-250 200 - 100 0 0 0;

4、所以经调幅后为a0 =0 455.2 834.91138.5 1593.9 1897.52201.1265.7 1556.0 1138.5759 455.4 151.8-531.3-986.7 -1290.3 -1518-1366.2- 1138.5 2200 a(ti)/290022a(ti)/29 a。 得: -531.3-189.8 151.8-75.900 0 u2 u1455.2 IIII U3 U2834.9 0455.2 455.2379.7 依次类推可以求出地面运动加速度的差值。 因为 c 2 Jkm 0.08 k 9000kN/m, m 250kN s2/m 代入可以算得c

5、240kN s/m 、表格第一行数据计算: Ki ki 6m/t 2 3c/t , t=0.05s 代入得 Ki 623400N / m 0 Pi (u； 6u/t 3UJ c(3u U- t -250 455.2 -113800N Ki U Pi U Pi /Ki 113800 / 623400 0.18mm II U 0 因为 U U /t2 6 Ui /t 3U i II 所以 U1 U Ui ( 0.18) / 0.052 II U t /2 432 所以 U1 432 0.05/2 10.8 起始时刻时： 2 0(mm/ s2) 相对加速度 =U； U 0 0 IIII 层间剪力 =

6、m（U ； U ） 0N 二、表格第二行数据计算： U2 H U2 U1 II U1 0.18 432 U1 U2 U； 834.9 10.8 455.2 379.7(m m/ s2) AAHIn P m( U； 6Ui/t 3Ui ) c(3Ui 5 t/2) =-250379.7-6*( -10.8)/0.05 -3*( -432)+240*3*( -10.8)+(-432)* 0.05/2 =-753293N au2 ap 因为AU /K 753293/6234001.21mm 6 AU /t2 6 Ui/t 3Ui 所以AU? 1.21 6/0.052 6 ( 10.8)/0.05 3

7、 ( 432) =-312(mm/s2) U2 Ui t AU t/2 =-432*0.05 -312*0.05/2= -29.4 (mm/s) IIII 相对加速度=Ug U 2 =455.2-432=23.2 (mm/s ) .11II 层间剪力=m(Ug U ) =250*23.2=5800N 依次类推，每次以上一行的数据作为起始值来计算下一行的 数值。可以得出所有数据。 t/s Au/mi 血7庶2) (1 S.2 蕊 230W 朋綁0 -IL3bOO O-IS -ICiS 张 TSOW Sll %6伍9迟黑 fU 15K. 172414 9 6I4W0 -11 曲曲IM 3 -lei

8、SmBS 4:7.12b5rat -41.54TSI&B -6-. 3981 (U5 【137.迪側-snsiQ粗 0 6I4W0 -123U?141.gQ -20Ll?Tti2Ulit 刖出5昭出4 -20l LhiT2?fcli *Sd. 02457966 鼻5 币El妙迦丁 -溺.(49刃5 4 6UW0 -12魏eu.曲-姒 7激肿朗 305. 355665 -2,53015S9S -1O2D2IS5 0l55 4砥H如弗占3点瑁出 0 6I4W0 -1258036?.T5 -20l4?tlE449 四h永网曲 12.3溺陨9 7里裁毗31 IB 田用门3形-迎,巧昵讯 i eiiwo

9、 -lliS4504.囹-ISl34S2c8E5 648.或I03Tf35i B92899 14S,& 0.S5 -931.0144023 -mii31139 9 S：4W -10546053.19 -1613；33：431 155.5SC11526 fJ.S5193ET5 -leiTiTZSlE Ol? T址20喘 坟工卿27聘 0 亂灿0 -aiE3(174. i9 -1131EaOfi&2 醐.亂弭慚 T*L 07856 -1 血 5曲竝 0.75 -W,W2如用跚 0 GLWO -5fiWefl.0U -&.13WH IK 1296Mt 迟翱坯T5 -1 證叩航2 U =1517,2U

10、J?g 151.72413? t 6UW0 -2900122.22叫肌刚神02-2迅则矇召34,711(2602 -2呢破也32 &J5 -n&,5l7Z41 2 球58Ca的 ) SL4W0 前劑氐5酿-O-aS-MLSEE -ZM,20S7m T1 詞師92開 唆9151205 .OlS -1137.B1034 hOLaWT 9 bUWO 132tHV. HT 2.155：册439 -&53.7422811 1l75536E35 -20T.3151022 Ol95 -531.M4403a 341.37&S103 m 趣威 2364117.2fil 3.码2295 595 -M.lT37oe

11、2 iXOSld册眾-2说砧內L3耳 1 1 -1391 E55ir24 311.3793103 m 633W0 25阖術朗4.12M3 P眺nm -5,训EZ捉-2乩跖巧LT? .fctJiM .55L3O8 2) 0 -435. 03702 -744. 972S103 -900. 1S56555 -1302, 4电&10干 -1537, 730153 -1761. 995259 20玄 T511193 -1059. 509844 -602.3630353 -203.72&0719 102. 1597946 393. 7072611 1042- ooe&9g 1437- 658S51 166

12、6. 614331 1310, 7 43939 1577- 59721 1201- 490502 62S. 74S2206 132. 5745124 i/ /(mm/s)相对加速度层间勇力/K 0 0 0 -10. 952175517.05539379 -10. 5236712689.5099483S -81. 6563829237.7S537S9 -136. 720527290,6573412 -207. 7253335353.S215T15 4271. 34S445 22377. 4=8709 59148. 84474 72&64, 33531 89705. 39237 -290. 21S

13、2187433. 0147414=109503. GS53 -329* 0240722Tb. 2G3360T113617. 0302 -350. 2680403495.6625701 -391, Slb3&2353b.54799E2 -ill. 96809554.894:6177 12515. 6425 133S86. 999S 138723. 65i4 -4=14, 50724=69557. 332208135333. 0521 -4=02* 11C57055i5. 431399113&357. 84=98 -3&G. 2176715510-974216 127743. 554 -304. 2257328451.咗619543 -226. 6186532376. 9591585 -139,