MATLAB弹塑性时程分析法编程

1、体系十质点质皆为30 kN-总忠择恬,- 融（B尼比；口唯迦契记录见表L麒遁谢 _.- 短力模型，妞I位移为2.0*屈服前结构抗则可刖置为P ：一 后抗瞄网妫丸卸嬲段的版移刚度为即小加翱廿算机程户 、二罕遇烈度要求将地I嫌谭一MATLA稀程：format short g;F(21,14)=0;xg=0 600 1100 1500 2100 2500 2900 350 2050 1500 1000600 200 -700 -1300 -1700 -2000 -1800 -1500 -700 -250 200 -100 0 0 0;xg1=xg*2200/max(xg);xg2=diff(xg1)

2、;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/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.05A2*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,:

3、);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)>0f(4)=0;else f(4)=9000;endenda=max(abs(F(:,13);b=max(abs(F(:,14);bxlswrite(' 表格 2.xls',F)计算书：课程设计计算书(题二)根据加速度调幅公式：a。a0,max a(ti)/am2、 am 290Gal 2900(mm/s )得：a02200 a(ti)/2900 22a(ti)/29a(tJ=0 600 1100 1

4、500 2100 2500 2900 350 205015001000600200-700-1300-1700-2000 -1800 -1500 -700 -250 200 -100 0 0 0;所以经调幅后为a0=0 455.2 834.9 1138.5 1593.91897.5 2201.1265.71556.01138.5759 455.4151.8-531.3-986.7-1290.3-1518-1366.2-1138.5 -531.3 -189.8 151.8 -75.9 0 0 0 U2 U1455.2 0 455.2U3 U2 834.9 455.2 379.7 32依次类推可以

5、求生地面运动加速度的差值。c因为 2km ,0.08 ,k 9000kN/m, m 250kN s2/m代人可以算得c 240kN s/m一、表格第一行数据计算：2Kiki 6m/t2 3c/t, t=0.05s代入得Ki623400N / m pmgUa 6U /t 3Ui) c(3UUi t/2)-250 455.2 -113800NKi U Pi U Pi / Ki113800 / 6234000.18mm起始时刻时： U 0 U ' 0 U '' 0因为U '' 6 U /t2 6 Ui'/t 3U i''所以 U1

6、9;'6 ( 0.18) / 0.052432 U ' Ui'' t U '' t/2所以U1'432 0.05/210.8相对加速度=Ug'' U'' 0 0 0(mm/ s2)." .层间剪力 = m(U g U ) 0N二、表格第二行数据计算：U2U10.18U 2'' U1''432U2'U1'10.8U '' 834.9 455.2 379.7(mm/s2) g p m(Ug 6Ui/t 3Ui) c(3Ui Ui t/2)

7、=-250379.7-6*(-10.8)/0.05-3*(-432)+240*3*(-10.8)+(-432)* 0.05/2=-753293N U2 APi /Ki753293/6234001.21mm因为 AU6 AU /t2 6 Ui/t 3Ui所以 AU；1.21 6/0.052 6 ( 10.8)/0.05 3 ( 432)=-312(mm/s2) U2 Ui t AU t/2=-432*0.05-312*0.05/2=-29.4 (mm/ s)相对加速度=Ug U =455.2-432=23.2(mm/s2)层间男力=m(Ug U )=250*23.2=5800N依次类推，每次以上

8、一行的数据作为起始值来计算下一行的 数值。可以得由所有数据。考虑到此为假定的结构弹塑性恢复力模型，屈服后结构抗侧移刚度为0,卸载阶段又恢复为 8000N/m,所以计算时当绝对值4U大于2mmM,要以ki=0N/m代入计算即可。t/s2AugP/(B5Jd/WA?i城如町瓦7M堂)u,伽*)(1®,2455.2咖。渊。0-1L38W也18-10.3IDS834 J3融T-JS3M3-L.a-312-244u/mn和执西/GW2'层间国力川C o000-o. 18-q义-L0.823.25颜输出结果:丁GD-r1C-H丁t7s,英kiAWm)中加Au/na：皿屋"二方/

9、'OLr：山九0u4551nM38m623W0-113793.10 -1132536258-«8. OE702-10l 55217550(X05巧£17黑洸酒 3I0&H8m璐网T6热T二?哪T T/注用泡注中3圈隧邀睡睡蟀蚯Oil3池拈花笳昵 43399000623W0力湖北腐-3.02230C761-1S.1&J453-l.M5bZ238E0.151137.9F103445F.K2tl38:i61山口尊叫©-5J 幽 8068 -4OL31M515出®5Hos T.斯史苏1215农工次E3皑麴乃94614W0用前5巩如七睢1胸

10、口1 Y葡.瞅MEATL岷蚯朗-9.B035H71E0.25IGK. 5517243脸例乃9j614眦-7filOT.3M无蛔娜的-然.2551®-吃短的包E间蛾吸W寸帆比黝:i一帜M脚TB8MM56B137M的潞-3B.30;B5E帕7区升10.35二阻陋组412S,65E172q6：蛔。ntfiLis71T设T曲曳2EJ泯TL族卷11Y4冽貂555.172411m, 24137 BG：i514UKT1蜡题*5-1 =在电把2工，二三垣41.547321 史，丁三三在L.招liyfcl034-常丸电眼期q6UW0-17141.50如、,心眦跳蚯如TlUm巧E %二二二区也575GL

11、62OMW-孤笆阻淤j6UW0-12733814.57-2U725GW65G：5.5S53K5-2,EM轮斑-1上的般0.55455.17Z413E-303.金般06I4W0-12580967. 15PQ4拈134期291.M7466:12.3K67639-122.S2T5281/I5LT闻股9吃，后腹诉q61蛔。-11384504. 33-1913485逸如犷35i 852893-1S3.4DK12E0.65-331.0M4B23-t3 5,1724133061期U-10.34E555J9-1"血血335. ffiCIKE61,95193B15TEE：皿世OlT团由笳E-30144

12、必的j6140-91E3(ffd.79-llSlESLb&222S, W54T9EK,田沌茹-173l(V5-128S. 655172-2ZT,5BG2Bj614W0第刎Bffl,帕T8阻她七1M，蹴刚33956*5-1三，土三二三OLE-1517.24115L724H79q6UW0-2SW122.4S2-t 31TW7702-22西变6“711呱见寸工。猥成0.35:敲51花。22% 5甜比69：i61蛔。5涮1选2-0廊乳斑2K.20E7121TJT3 的 £K210.91512叱(L9-list 甄 OU诋弼E74614W0i32ur. n2.15 / EE 妗 9-E

13、52,. T28114L75536E：5-27.81513220.95-531. rawBsa341.3T93103颗0623400:三二二三：3/52295395r维二7"战63DH32-笈区的13E1祖.附£3濒Q强嫄257305,冏1J2EM3-1&4.735 咂般胸 皿86他帽55h党二加的</ (irco/s'2)/(mn/s)中目对加速度八的，丁2)层间剪力加0000-438. 09702-10. 952175517. 0S5393734271. 343445-744. 9728103一心。. 52367126日9. 50994836223

14、77. 467rg-900. 1356S55-81. 6563829237. 795378959i49. 84474-1302. 44:6107-136.72OS27290. 6573412T2&64. 33531-'1537. 730153一工GIL 72S3335358. 821571589705. 39287-17S1, 935259-290. 21S2187438 0147414109503, 6S53209. 7511193-329. 0240T22马也 2633607118817. 0902-1059. 509844-350. 26S0403495. 6625701

15、12315. 6425-602. 3B30353-391. 8153S23535. 54;79992133886. 999S-203.72S0719-411. 96809554. 8946177138723. 65生4102. 1597946-4:14. 50724:69557. 3322084：139333. 0521393. 7072S11一电。2. 11057055*i5. 4313991136S57- 84：981042. 00e&99-3S6. 217E715510. 974216127742. 5541137. 668851-304. 2257328451. 4619548112fi&5, 48371&66. 514331-226. 6186532376. 59158594:239. 7B9&11310. 74=3939-139. 6846965293. 5C«25673375. £39991