Loutters的个人资料

【求助，matlab动力学编程出错】 %-------------------------------------------------------- % Ex.2.2 plane 9-bar truss %-------------------------------------------------------- % To compute the dynamic response of node 4 on the vertical direction % % E=2.0e11; A=2.5e-3; density=7860; node_number=6; elment_number=9; nc=[0 0;4 0;4 3;8 0;8 3;12 0]; %node_coordinate en=[1 2;1 3;2 3;2 4;3 4;3 5;4 5;4 6;5 6]; % element_node ed(1:node_number,1:2)=1; % elment_displacement dof=0; for loopi=1:node_number for loopj=1:2 dof=dof+1; ed(loopi,loopj)=dof; end end ek=E*A*[1 0 -1 0;0 0 0 0;-1 0 1 0;0 0 0 0]; % elment_stiffness_matrix em=(density*A)/6*[2 0 1 0; ... %elment_mass_matrix 0 2 0 1; ... 1 0 2 0; ... 0 1 0 2]; k(1:dof,1:dof)=0; %tructural_stiffness_matrix m=k; % structural_mass_matrix, same size with k theta(1:9)=0; el(1:9)=0; e2s(1:4)=0; %index of transform the elment displament number to structural for loopi=1:elment_number for zi=1:2 e2s((zi-1)*2+1)=ed(en(loopi,zi),1); e2s((zi-1)*2+2)=ed(en(loopi,zi),2); end el(loopi)=sqrt((nc(en(loopi,1),1)-nc(en(loopi,2),1))^2+(nc(en(loopi,1),2)-nc(en(loopi,2),2))^2); theta(loopi)=asin((nc(en(loopi,1),2)-nc(en(loopi,2),2))/el(loopi)); lmd=[cos(theta(loopi)) sin(theta(loopi)); -sin(theta(loopi)) cos(theta(loopi))]; t=[lmd zeros(2); zeros(2) lmd]; dk=t'*ek*t/el(loopi); dm=t'*em*t*el(loopi); for jx=1:4 for jy=1:4 k(e2s(jx),e2s(jy))=k(e2s(jx),e2s(jy))+dk(jx,jy); m(e2s(jx),e2s(jy))=m(e2s(jx),e2s(jy))+dm(jx,jy); end end end c=0*k+0*m; nt=1500;dt=0.0001; time=0:dt:nt*dt; q0=zeros(dof,1); dq0=zeros(dof,1); bcdof=zeros(dof,1); fd=zeros(dof,nt); for i=1:nt fd(6,i)=200; end [acc,vel,dsp]=TransResp4(k,c,m,fd,bcdof,nt,dt,q0,dq0); plot(time, dsp(8,:)) xlabel('Time(seconds)') ylabel(' Vertical displ. (m)') —————————————————————————————————————————————— 用到的函数 wilson —————————————————————————————————————————————— function [acc,vel,dsp]=TransResp4(kk,cc,mm,fd,bcdof,nt,dt,q0,dq0) %-------------------------------------------------------------------------- % Purpose: % The function subroutine TransResp4.m calculates transient response of % a structural system using Wilson integration scheme. % Synopsis: % [acc, vel, dsp]=TransResp4(kk,cc,mm,fd,bcdof,nt,dt,q0,dq0) % Variable Description:

【即连续失眠后，昨晚彻底沦落为整夜失眠了】

【求助，徐斌老师的matlab有限元结构动力学分析与工程应用】例子2.4.2，运行不了。请假一下各位高手。跪求指导

【暴暴家乡照】

【一直不想做学霸，可是不知道为什么一直被逼着做学霸】

【本部的有没有，晚上好想去喝酒】

【好忧伤，好忧伤，不知道什么年后哪天来学校】

【发个贴吧，说明我还活着】悲剧的是，我昨天刚回学校，大家竟都放假了。

【今天是在吧里的最后一晚了，不知道什么时候再回来】加油，一切都会好的

【贴吧怎么了，老是看不到回复】 RT

【那些年，那些事】

【切糕体】 “奥巴马宣布，超过万亿的国债将以二百斤切糕的形式归还中国。” “这就是一个保时捷撞翻切糕车之后,保时捷车主弃车逃逸的时代。” “以后谁要给官员送礼一律上切糕。” 中国政府宣布，人民币与切糕挂钩，各国货币与人民币挂钩。德国荷兰政府要求美国将存于美联储的黄金兑换成切糕遭拒。

【乱七八糟】-【只为逗你一笑】

【庆祝啊，不容易】 @公厕女鬼超过你了。嘿嘿。终于到第五页了，不容易啊。

【若果我爱你】-转载如果我爱你　　我就不做那电子邮件　　只是让你对着屏幕　　发出一声无奈的浅笑　　只是让你向着虚无空把手招　　我要做特快专递　　只要思念的信鸽飞抵窗前　　我就会在24小时内飞抵你的床前　　思念可以让爱生，也可以让爱死　　　　如果我爱你　　我就不做那哪怕是最亮的一棵星　　只让你感觉恍若隔世的温情　　只让你做遥不可及的梦　　我要做你床头的一盏灯　　让你于每个夜晚在我温情的抚摸中入梦　　让我在你的枕旁倾听你均匀的呼吸声　　距离可以产生美　　距离也可以让爱疲惫　　　　如果我爱你　　我就不做那令人垂涎欲滴的美味　　让你在飨之时摧残着胃　　让你在品之后因不可多得而失落　　我要做你惯常吃的小菜儿　　让你在咀嚼时谈不上多滋多味　　只是三日不吃便回味　　爱不在美味，只要合口味　　　　如果我爱你　　我就不做那笔挺的西装　　让你在穿时有些累　　我要做你贴身的一件休闲服　　在你站时，让你得体随和　　在你躺时轻松自在　　爱要刻意，便泛味　　　　如果我爱你　　我就不做那装帧精美的精品书　　让你在读时净手　　怕渎了书的圣洁　　让你在放时不敢折痕　　怕亵了书的高雅　　我要做一本平装书　　让你能躺于床上信手拈来　　即使印满你的指纹　　那正是爱打下的烙印　　让你在放时任意折叠　　即使折痕累累　　那也是爱打上的一个幸福结　　爱一旦神圣，凡人便承受不起　　　　如果我爱你　　我就不做那价格不菲的钢琴　　让你在艺术门前止步　　让你在金钱面前却步　　我要做一把廉价的口琴　　让你在烦闷时吹上几口　　烦闷便从琴缝儿中溜走　　让你在伤感时吹上一曲　　伤感便被音符带走　　爱一旦成为奢侈品　　爱就只能当摆设

【经典电视剧】——【大家都说说喜欢的角色】神雕侠侣——陆无双，郭襄倚天屠龙记——小昭射雕英雄传——黄蓉侠客行——阿　绣，侍　剑笑傲江湖——任盈盈天龙八部—— 阿朱

【这不正验证了，我大学做了6年的光棍么】

【刚刚带本科实验，彻底伤了】能动学院，大三某专业。试验总共3组，每组10人，结果你知道我看到了什么么？第一组，9个男生1个女生，第二组，9个男生1个女生，第三组，10个男生0个女生， OH,MY GOD

【又被光棍了，伤不起】

【刚看完80后，大家感觉怎么样】

【不得不说，南校区的浴室水太热了】脱了层皮，还好皮糙肉厚

【十元三毛的爱你会接受吗？】这样的爱你会接受吗？一心一意哦!

【讨论帖】-【学姐和学妹哪个更适合做女朋友】学姐优势在哪，学妹优势在哪呢, 你是喜欢学妹还是学姐?

【好忧伤好忧伤】过几天就23岁生日了，又老了一岁了，伤不起啊，妹纸还没有

【学霸的标准是什么，大家集体讨论一下】想知道自己是不是学霸

【看不到路怕妹纸的帖子，睡不着的有木有】 rt@路帕

【学霸难道就注定单身么】建议学校自习室分为：情侣秀恩爱，单身求邂逅，基友一起走和学霸纯学术四种，以免误伤他人

【今天被发传单了，作为资深的屌丝，伤不起啊】大家都懂的，我就不说什么了

【签上你的生日,看看在帝国理工大学里有没有人与你同月同日生】 RT 12 15

【早安，学弟学妹们】早安上理，早安安徽