minuskmatrix.m

薄壁结构有限元计算

Ke4=zeros(14,14);
for i=1:14
    Ke4(i,i)=Ke3(i,i);
end
for i=1:7
    Ke4(i,7+i)=Ke3(i,7+i);
end
for i=1:7
    Ke4(7+i,i)=Ke3(7+i,i);
end
Ke4(4,7)=Ke3(4,7);
Ke4(7,4)=Ke3(7,4);
Ke4(11,7)=Ke3(11,7);
Ke4(7,11)=Ke3(7,11);
Ke4(4,14)=Ke3(4,14);
Ke4(14,4)=Ke3(14,4);
Ke4(11,14)=Ke3(11,14);
Ke4(14,11)=Ke3(14,11);
for i=2:6
    Ke4(i,8-i)=Ke3(i,8-i);
end
for i=2:6
    Ke4(7+i,8-i)=Ke3(7+i,8-i);
end
for i=2:6
    Ke4(i,15-i)=Ke3(i,15-i);
end
for i=2:6
    Ke4(7+i,15-i)=Ke3(7+i,15-i);
end
Ke3=Ke4;