原创力文档- 1、本文档共5页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
有限元分析大作业程序部分
程序部分
1，计算单元刚度矩阵（输入E,A,I,L）然后输出k（6*6）刚度矩阵：
function k=Beam2D2Node_Stiffness(E,I,A,L)
k=[E*A/L,0,0,-E*A/L,0,0;
0,12*E*I/(L^3),6*E*I/(L^2),0,-12*E*I/(L^3),6*E*I/(L^2);
0,6*E*I/(L^2),4*E*I/L,0,-6*E*I/(L^2),2*E*I/L;
-E*A/L,0,0,E*A/L,0,0;
0,-12*E*I/(L^3),-6*E*I/(L^2),0,12*E*I/(L^3),-6*E*I/(L^2);
0,6*E*I/(L^2),2*E*I/L,0,-6*E*I/(L^2),4*E*I/L];
2，矩阵组装
function y =Beam2D2Node_Assemble(KK,k,i,j)
%以上函数进行单元刚度矩阵的组装
%输入单元刚度矩阵k，单元的节点编号i、j
%输出整体刚度矩阵KK
%-----------------------------------------
DOF(1)=3*i-2;
DOF(2)=3*i-1;
DOF(3)=3*i;
DOF(4)=3*j-2;
DOF(5)=3*j-1;
DOF(6)=3*j;
for n1=1:6
for n2=1:6
KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2);
end
end
y = KK;
输入数据，求k1,k2刚度阵：
E=3E11;
I=6.5E-7;
A=6.8E-4;
L1=1.44;
L2=0.96;
k1=Beam2D2Node_Stiffness(E,I,A,L1);
k2=Beam2D2Node_Stiffness(E,I,A,L2);
k1 =
1.0e+08 *
1.4167 0 0 -1.4167 0 0
0 0.0078 0.0056 0 -0.0078 0.0056
0 0.0056 0.0054 0 -0.0056 0.0027
-1.4167 0 0 1.4167 0 0
0 -0.0078 -0.0056 0 0.0078 -0.0056
0 0.0056 0.0027 0 -0.0056 0.0054
k2 =
1.0e+08 *
2.1250 0 0 -2.1250 0 0
0 0.0264 0.0127 0 -0.0264 0.0127
0 0.0127 0.0081 0 -0.0127 0.0041
-2.1250 0 0 2.1250 0 0
0 -0.0264 -0.0127 0 0.0264 -0.0127
0 0.0127 0.0041 0 -0.0127 0.0081
3，建立整体刚度方程：
T=[0,1,0,0,0,0;-1,0,0,0,0,0;0,0,1,0,0,0;0,0,0,0,1,0;0,0,0,-1,0,0;0,0,0,0,0,1];
k3=T*k2*T;
KK=zeros(12,12);
KK=Beam2D2Node_Assemble(KK,k1,1,2);
KK=Beam2D2Node_Assemble(KK,k3,3,1);
KK=Beam2D2Node_Assemble(KK,k3,4,2)
T =
0 1 0 0 0 0
-1 0 0 0 0 0
0 0 1 0 0 0
0 0 0 0
文档评论（0）