该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
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;
请指出错误,万分感谢