原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
姓名：吴佳侣
班级：机自11-13
学号例： 如图所示有限元模型，弹性模量为 ，厚度为 ，为简化计算取 ，求整体刚度矩阵。
部分分析思路：
各单元信息
单元编号 ① ② ③ ④ 整体编码 1,2,3 2,4,5 5,3,2 3,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/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
K3=transpose(B2)*D*B2*t*mianji;
%%%%%%%%%%%%单元3
bi=0;ci=-a;
bj=a;cj=-a;
bm=-a;cm=0;
mianji=a^2/2;
B3=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
k2=transpose(B3)*D*B3*t*mianji;
%%%%%%%%%%%%单元4
bi=0;ci=a;
bj=-a;cj=-a;
bm=a;cm=0;
mianji=a^2/2;
B4=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
K4=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_2
文档评论（0）