matlab求刚度,求整体刚度矩阵matlab程序.doc

姓名：吴佳侣

班级：机自11-13

学号例： 如图所示有限元模型，弹性模量为 ，厚度为 ，为简化计算取 ，求整体刚度矩阵。

部分分析思路：

各单元信息

单元编号①②③④整体编码1,2,32,4,55,3,23,5,6局部编码

以整编码体表示的单元刚度矩阵

整体刚度矩阵

程序：

%%%%%%%%%%%%%%%求整体刚度矩阵

clc

clear all

syms E u t a%%定义变量

%%%%%%%%%%%%%%%%%%%%%%%%写出单元刚度矩阵,单元面积为a^2/2

%%%%%%%%%%%%单元1

bi=0;ci=a;

bj=-a;cj=-a;

bm=a;cm=0;

mianji=a^2/2;

B1=1/2/mianji*[bi 0 bj 0 bm 0

0 ci 0 cj 0 cm

ci bi cj bj cm bm];

u=0;E=1;t=1;

D=E/(1-u^2)*[1 u 0

u 1 0

0 0 (1-u)/2]; %%%弹性矩阵D

k1=transpose(B1)*D*B1*t*mianji;

%%%%%%%%%%%%单元2

bi=0;ci=a;

bj=-a;cj=a;

bm=a;cm=0;

mianji=a^2/2;

B2=1