求整体刚度矩阵matlab程序(共7页).doc

姓名：吴佳侣班级：机自11-13学号：03111323例： 如图所示有限元模型，弹性模量为 ，厚度为 ，为简化计算取 ，求整体刚度矩阵。部分分析思路：各单元信息单元编号整体编码1,2,32,4,55,3,23,5,6局部编码以整编码体表示的单元刚度矩阵整体刚度矩阵程序：%求整体刚度矩阵clcclear allsyms E u t a%定义变量%写出单元刚度矩阵,单元面积为a2/2%单元1bi=0;ci=a;bj=-a;cj=-a;bm=a;cm=0;mianji=a2/2;B1=1/2/mianji*bi 0 bj 0 bm 0 0 ci 0 cj 0 cm ci bi cj bj cm bm;