matlab单元刚度矩阵,求助：关于有限元三角形单元合成总刚度矩阵怎么处理

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

clc;

clear all;

mu=0.3;E=2e5;t=1;

D=(E/(1-mu^2))*[1 mu 0;mu 1 0;0 0 (1-mu)/2];

%%载入节点和单元信息%%

load Nodefilename.dat;

NodeInformation=Nodefilename;

A1=size(NodeInformation);

NP=A1(1);

load Elementfilename.dat;

ElememtInformation=Elementfilename;

A2=size(ElememtInformation);

NE=A2(1);

%%单元循环%%

K=zeros(2*NP,2*NP);

k=zeros(6,6);

for e=1:NE

i=ElememtInformation(e,7);

j=ElememtInformation(e,8);

m=ElememtInformation(e,9);

xi=NodeInformation(i,2);

yi=NodeInformation(i,3);

xj=NodeInformation(j,2);

yj=NodeInformation(j,3);

xm=NodeInformation(m,2);

ym=NodeInformation(m,3);

a=(xi*(yj-ym)+xj*(ym-yi)+xm*(yi-yj))/2;

bi=yj-ym;

bj=ym-yi;

bm=yi-yj;

ci=xm-xj;

cj=xi-xm;

cm=xj-xi;

B=/(2*a);

k=B.'*D*B*t*a;

K(2*i-1,2*i-1)=K(2*i-1,2*i-1)+k(1,1);

K(2*i-1,2*i)=K(2*i-1,2*i)+k(1,2);

K(2*i-1,2*j-1)=K(2*i-1,2*j-1)+k(1,3);

K(2*i-1,2*j)=K(2*i-1,2*j)+k(1,4);

K(2*i-1,2*m-1)=K(2*i-1,2*m-1)+k(1,5);

K(2*i-1,2*m)=K(2*i-1,2*m)+k(1,6);

K(2*i,2*i-1)=K(2*i,2*i-1)+k(2,1);

K(2*i,2*i)=K(2*i,2*i)+k(2,2);

K(2*i,2*j-1)=K(2*i,2*j-1)+k(2,3);

K(2*i,2*j)=K(2*i,2*j)+k(2,4);

K(2*i,2*m-1)=K(2*i,2*m-1)+k(2,5);

K(2*i,2*m)=K(2*i,2*m)+k(2,6);

K(2*j-1,2*i-1)=K(2*j-1,2*i-1)+k(3,1);

K(2*j-1,2*i)=K(2*j-1,2*i)+k(3,2);

K(2*j-1,2*j-1)=K(2*j-1,2*j-1)+k(3,3);

K(2*j-1,2*j)=K(2*j-1,2*j)+k(3,4);

K(2*j-1,2*m-1)=K(2*j-1,2*m-1)+k(3,5);

K(2*j-1,2*m)=K(2*j-1,2*m)+k(3,6);

K(2*j,2*i-1)=K(2*j,2*i-1)+k(4,1);

K(2*j,2*i)=K(2*j,2*i)+k(4,2);

K(2*j,2*j-1)=K(2*j,2*j-1)+k(4,3);

K(2*j,2*j)=K(2*j,2*j)+k(4,4);

K(2*j,2*m-1)=K(2*j,2*m-1)+k(4,5);

K(2*j,2*m)=K(2*j,2*m)+k(4,6);

K(2*m-1,2*i-1)=K(2*m-1,2*i-1)+k(5,1);

K(2*m-1,2*i)=K(2*m-1,2*i)+k(5,2);

K(2*m-1,2*j-1)=K(2*m-1,2*j-1)+k(5,3);

K(2*m-1,2*j)=K(2*m-1,2*j)+k(5,4);

K(2*m-1,2*m-1)=K(2*m-1,2*m-1)+k(5,5);

K(2*m-1,2*m)=K(2*m-1,2*m)+k(5,6);

K(2*m,2*i-1)=K(2*m,2*i-1)+k(6,1);

K(2*m,2*i)=K(2*m,2*i)+k(6,2);

K(2*m,2*j-1)=K(2*m,2*j-1)+k(6,3);

K(2*m,2*j)=K(2*m,2*j)+k(6,4);

K(2*m,2*m-1)=K(2*m,2*m-1)+k(6,5);

K(2*m,2*m)=K(2*m,2*m)+k(6,6);

end

K;

请指出错误，万分感谢