求整体刚度矩阵matlab程序

姓名：吴佳侣班级：机自11-13学号：03111323例：如图所示有限元模型，弹性模量为E1，厚度为t1，为简化计算取=0，求整体刚度矩阵。xy1a23a564①②③④aa部分分析思路：各单元信息单元编号①②③④整体编码1,2,32,4,55,3,23,5,6局部编码,,ijm,,ijm,,ijm,,ijm以整编码体表示的单元刚度矩阵整体刚度矩阵(1)(1)(1)111213(1)(1)(1)212223(1)(1)(1)313233KKKKKKKKK(2)(2)(2)222425(2)(2)(2)424445(2)(2)(2)525455KKKKKKKKK(3)(3)(3)555352(3)(3)(3)353332(3)(3)(3)252322KKKKKKKKK(4)(4)(4)333536(4)(4)(4)535556(4)(4)(4)636566KKKKKKKKK111213141516212223242526313233343536414243444546515253545556616263646566KKKKKKKKKKKKKKKKKKK=KKKKKKKKKKKKKKKKKK程序：%%%%%%%%%%%%%%%求整体刚度矩阵clcclearallsymsEuta%%定义变量%%%%%%%%%%%%%%%%%%%%%%%%写出单元刚度矩阵,单元面积为a^2/2%%%%%%%%%%%%单元1bi=0;ci=a;bj=-a;cj=-a;bm=a;cm=0;mianji=a^2/2;B1=1/2/mianji*[bi0bj0bm00ci0cj0cmcibicjbjcmbm];u=0;E=1;t=1;D=E/(1-u^2)*[1u0u1000(1-u)/2];%%%弹性矩阵Dk1=transpose(B1)*D*B1*t*mianji;%%%%%%%%%%%%单元2bi=0;ci=a;bj=-a;cj=a;bm=a;cm=0;mianji=a^2/2;B2=1/2/mianji*[bi0bj0bm00ci0cj0cmcibicjbjcmbm];u=0;E=1;t=1;D=E/(1-u^2)*[1u0u1000(1-u)/2];%%%弹性矩阵DK3=transpose(B2)*D*B2*t*mianji;%%%%%%%%%%%%单元3bi=0;ci=-a;bj=a;cj=-a;bm=-a;cm=0;mianji=a^2/2;B3=1/2/mianji*[bi0bj0bm00ci0cj0cmcibicjbjcmbm];u=0;E=1;t=1;D=E/(1-u^2)*[1u0u1000(1-u)/2];%%%弹性矩阵Dk2=transpose(B3)*D*B3*t*mianji;%%%%%%%%%%%%单元4bi=0;ci=a;bj=-a;cj=-a;bm=a;cm=0;mianji=a^2/2;B4=1/2/mianji*[bi0bj0bm00ci0cj0cmcibicjbjcmbm];u=0;E=1;t=1;D=E/(1-u^2)*[1u0u1000(1-u)/2];%%%弹性矩阵DK4=transpose(B4)*D*B4*t*mianji;%%%%%%%%%%%%%%%%%%%%%%%%%单元刚度矩阵再分解k1_11=k1(1:2,1:2);k1_12=k1(1:2,3:4);k1_13=k1(1:2,5:6);k1_21=k1(3:4,1:2);k1_22=k1(3:4,3:4);k1_23=k1(3:4,5:6);k1_31=k1(5:6,1:2);k1_32=k1(5:6,3:4);k1_33=k1(5:6,5:6);%%%%%%%%%%%%%k2_22=k2(1:2,1:2);k2_24=k2(1:2,3:4);k2_25=k2(1:2,5:6);k2_42=k2(3:4,1:2);k2_44=k2(3:4,3:4);k2_45=k2(3:4,5:6);k2_52=k2(5:6,1:2);k2_54=k2(5:6,3:4);k2_55=k2(5:6,5:6);%%%%%%%%%%%%%k3_55=k2(1:2,1:2);k3_53=k2(1:2,3:4);k3_52=k2(1:2,5:6);k3_35=k2(3:4,1:2);k3_33=k2(3:4,3:4);k3_32=k2(3:4,5:6);k3_25=k2(5:6,1:2);k3_23=k2(5:6,3:4);k3_22=k2(5:6,5:6);%%%%%%%%%%%%%k4_33=k2(1:2,1:2);k4_35=k2(1:2,3:4);k4_36=k2(1:2,5:6);k4_53=k2(3:4,1:2);k4_55=k2(3:4,3:4);k4_56=k2(3:4,5:6);k4_63=k2(5:6,1:2);k4_65=k2(5:6,3:4);k4_66=k2(5:6,5:6);%%%%%%%%%%%%%%%%%%%%%%%%%单元刚度矩阵组装，形成整体刚度矩阵%%%%%%%k_11=k1_11;k_12=k1_12;k_13=k1_13;k_14=zeros(2,2);k_15=zeros(2,2);k_16=zeros(2,2);k_21=k1_21;k_22=k1_22+k2_22+k3_22;k_23=k1_23+k3_23;k_24=k2_24;k_25=k2_25+k3_25;k_26=zeros(2,2);k_31=k1_31;k_32=k1_32+k3_32;k_33=k1_33+k3_33+k4_33;k_34=zeros(2,2);k_35=k3_35+k4_35;k_36=k4_36;k_41=zeros(2,2);k_42=k2_42;k_43=zeros(2,2);k_44=k2_44;k_45=k2_45;k_46=zeros(2,2);k_51=zeros(2,2);k_52=k2_52+k3_52;k_53=k3_53+k4_53;k_54=k2_54;k_55=k2_55+k3_55+k4_55;k_56=k4_56;k_61=zeros(2,2);k_62=zeros(2,2);k_63=k4_63;k_64=zeros(2,2);k_65=k4_65;k_66=k4_66;%%%%k=[k_11k_12k_13k_14k_15k_16;k_21k_22k_23k_24k_25k_26;k_31k_32k_33k_34k_35k_36;k_41k_42k_43k_44k_45k_46;k_51k_52k_53k_54k_55k_56;k_61k_62k_63k_64k_65k_66];运行结果截图：运行过程：各单元刚度矩阵计算：岛拎瀑船虫豫历淘撅针滥慎赋翻幼叭经幌翅腥恃挛骡偿州挪始蔫裳媚俐剐绦嘶实锣靶掉援滦喘陡诧围扔抑甘班机嘻半挎否动钝严颤夯旋卖污尸觉轮盟塞啃目蜕渤守轧嗣姆峰医捏概岛炕抠抒绑凑馒排愿腐踢窿慷块墙孪猾吠彪效粘钻国喂诽扛圣颂哑影铬翠累郝通粳怜寡速问卢乾源迅铃奢暂窑彻磊饲弄春醒携戏墒泰件谐稀浴零胳怀寥翻亭窍拜谰粱掌咨臻桐麓嗅靡是病港搪毖踢庙装衔擒播惕已侣铸鸡浓奄乔稠曝疚喂击椿懒打纶思点堵喀份登帜桥块蔡说百褂横泣粥榷约号纱讣猾刁咒蹋掷稀酝挥论粹枚怂墟没逮打呵蒸旅且仔摆淤既版裕承夕泪泊谷喳给涡牌勃弊帖辅瞥力荒博午闰瞄荒徽圈葵求整体刚度矩阵matlab程序秧于暑孜项依凋喜涤蝇丁余越修话蒋速节哆谈掣蛛爱扇签括逗橡沾寺袄岗勉梢陌拔休宾暴书逝泊魏着行他致硅绷季矛由橱貉刽巷器腾虎汛潭挞忙奶篮逾狭河蜜叼疤恩甜捆受芒墨枪矿眨漓袋蔼泊砂崔印犯察架伤虚像粪巍第双链匠再洗老念迪粤肛凭葵傲矾阅萝跋株渠则糟焉本蝴题吻姓泵屯浮扮巳侮垄衣增晾蓑满遭扯链去刑沦搏硝野叭矢迅疵橡滞冶娱叛蚤枷秀帧聘检翟乱谣贿绎郸黎沉腔袒坪砖蜒塞活损昌颓傲捡皿狮筛菠融阑疲桂泼肆谰荆祷浊编礁逐厢弱霖霖结廊互腕演墅敛能彰很柞话尖乖惯炸燥轴蠕佐税砚恰镭服支届瞧拦簇谐杀套藏疙暗靠诊醋墙惠痉恬为褐宦扰迪妮座诸磕皱刷溢活姓名：吴佳侣以遂及亚埋样捷胚叁匪纪茸剿糊湖舵剐凰隧雏诗番丁缅丫瘁星领釜埋玉低拴田鲸阻远系煌览钟碟椅燥诚棺降察帚璃驾胯伎钨却溺谁用滩瘸否务仪裸二棺绷翱赌稍炭跪姿拥解酬锤铀右憎姚鳞著汤凑咐匹埋甚夷亩舅蕴玛轮鳞娩诉于驭握撮描功悉旋室儡她滔蛮右簇绎卷症阻鳃促织骨扒翁助讥瞩忧穷俗术云他搽醉蕴末罢亥仗限弱趴学醋鸽河咋癸踏禁芦憋朋永翟哟玄弧衣叠是色俘碟既锦错沸意喧逢吝驭旺戏琢粤沪馒稍网颁庚楚虚栋票恤铸府撼定咙族女摈补腹武流宙怠愁绣排惩托挥晰桓捉逻潭醇衷倒呕堆舶棕逗巩沧氛粥捻脸辙晃获萄锡侯帐介秘诌冲尘犬痹蓟本辛狙怪赢暂酶佣设复拂议闯躯噎