nels=6;
nxe=2;
aa=3.0;
bb=3.0;
A=zeros(nel*8,2);
for iel=1:nels
iq=fix((iel-1)/nxe+1);
ip=fix(iel-(iq-1)*nxe);
coord(1:3,1)=aa*(ip-1);
coord(5:7,1)=aa*ip;
coord(4,1)=0.5*(coord(3,1)+coord(5,1));
coord(8,1)=0.5*(coord(7,1)+coord(1,1));
coord(1,2)=-bb*iq;
coord(7:8,2)=-bb*iq;
coord(3:5,2)=-bb*(iq-1);
coord(2,2)=0.5*(coord(1,2)+coord(3,2));
coord(6,2)=0.5*(coord(5,2)+coord(7,2));
A(((iel-1)*8+1):iel*8)=coord
end
nels=6时候数组是48*2,而不是楼主说的29*2.不过不管到底是多大,上边的程序都可以实现 |