Combining smaller matrices into a bigger one via a loop

 clc;
Ln=5;%input('Length between node''Ln');
Bn=2.5;%input('Beam between node''Bn');
h=0.02;%input('Plate thickness between node''h');
M=7600;%input('Mass of the Helicopter''M');
Lweb=0.07;%input('Length of web of the stiffner''Lweb');
Lf=0.07;%input('Length of flange of the stiffner''Lf');
Tweb=0.015;%('web thickness''Tweb');
Tf=0.015;%('flange thickness''Tf');
Ss=0.4;%input('Stiffners Spacing''Ss');
ymax=0.005;%input('Maximum Displacement of materials''ymax');
E=207*10^9;
Es=69*10^9;
G=69*10^12;
g=9.81;
% Area of Hellideck
Ah=Bn*Ln;
Aweb=2*Tweb;
Af=Lf;
%As=Aweb+Af;
% Total length of stiffner
LTf=Aweb+Af
% Number of longtidunal Spacing
L1=LTf+Ss;
Nl=Ln/L1
Nv=Bn/L1
Np=Nl*12
Nl1=2*Nl
%Nl2=hm.Nl1
%Nl3=int(Nl2)
Nv1=4*Nv;
%Nv2=int(Nv1)
K11=zeros(Np,Np);
%K1=zeros(Nv1,Nl1)
Iy=(Ss*h^3)/12;
Iz=(h*Ss^3)/12;
F=M*g;
Q=0
K=zeros(30,60)
for L=L1:L1:Ln
    q=F/L;
   for B=L1:L1:Bn
       As=h*B;
       A=(As*E)/L;
       I=(B*h^3)/12;
       Asy=B*h;
       Qy=(12*E*Iz)/(G*Asy.*L^2);
       By=(E.*I)/((1+Qy).*L^3);
       Qz=(12*E*Iy)/(G*Asy.*L^2);
       Bz=(E*Iy)/((1+Qz).*L^3);
R=[1,0,0,0,0,0,-1,0,0,0,0,0;0,cos(Q),sin(Q),0,cos(Q),sin(Q),0,cos(Q),sin(Q),0,cos(Q),sin(Q);0,-sin(Q),cos(Q),0,-sin(Q),cos(Q),0,-sin(Q),cos(Q),0,-sin(Q),cos(Q);0,0,0,-1,0,0,0,0,0,1,0,0;0,-cos(Q),-sin(Q),0,-cos(Q),-sin(Q),0,-cos(Q),-sin(Q),0,-cos(Q),-sin(Q);0,sin(Q),-cos(Q),0,sin(Q),-cos(Q),0,sin(Q),-cos(Q),0,sin(Q),-cos(Q);-1,0,0,0,0,0,1,0,0,0,0,0;0,cos(Q),sin(Q),0,cos(Q),sin(Q),0,cos(Q),sin(Q),0,cos(Q),sin(Q);0,-sin(Q),cos(Q),0,-sin(Q),cos(Q),0,-sin(Q),cos(Q),0,-sin(Q),cos(Q);0,0,0,-1,0,0,0,0,0,1,0,0;0,-cos(Q),-sin(Q),0,-cos(Q),-sin(Q),0,-cos(Q),-sin(Q),0,-cos(Q),-sin(Q);0,sin(Q),-cos(Q),0,sin(Q),-cos(Q),0,sin(Q),-cos(Q),0,sin(Q),-cos(Q)];    
K1=[A,0,0,0,0,0,-A,0,0,0,0,0;0,12.*By,0,0,0,6.*L.*By,0,-12.*By,0,0,6.*L.*By,0;0,0,12.*Bz,-6.*L.*Bz,0,0,0,0,-12.*Bz,-6.*L.*Bz,0,0;0,0,-6.*L.*Bz,(4+Qz).*L.^2.*Bz,0,0,0,0,6.*L.*Bz,(2-Qz).*L.^2.*Bz,0,0;0,6.*L.*By,0,0,(4+Qy).*L^2.*By,0,0,-6.*L.*By,0,0,(2-Qy).*L.^2.*By,0;-A,0,0,0,0,0,A,0,0,0,0,0;0,-12.*By,0,0,-6.*L.*By,0,0,12.*By,0,0,-6.*L.*By,0;0,0,-12.*Bz,6.*L.*Bz,0,0,0,0,12.*Bz,6.*L.*Bz,0,0;0,0,-6.*L.*Bz,(2-Qz).*L.^2.*Bz,0,0,0,6.*L.*Bz,(4+Qz).*L.^2.*Bz,0,0,0;0,6.*L.*By,0,0,(2-Qy).*L.^2.*By,0,0,-6.*L.*By,0,0,(4+Qy).*L.^2.*By,0;0,0,0,0,0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0,0,0,0,0];
R1=R.';
KG=K1.*R.*R1;
f=[F;q.*(L/2);q.*(L^2/12);q.*(L/2);-q.*(L^2/12);0;F;q.*(L/2);-q.*(L^2/12);q.*(L/2);q.*(L^2/12);0];
%k11=sum(:1,) :12:Np
% applying downbary conditions

 KG(12,:)=[];
KG(:,12)=[];
KG(11,:)=[];
KG(:,11)=[];
KG(10,:)=[];
KG(:,10)=[];
KG(8,:)=[];
KG(:,8)=[];
KG(7,:)=[];
KG(:,7)=[];
KG(6,:)=[];
KG(:,6)=[];
k=KG

 f(12,:)=[];
f(11,:)=[];
f(10,:)=[];
f(8,:)=[];
f(7,:)=[];
f(6,:)=[];

   end

 end

 %%The KGs are displayed but id loved to have them arranged into a bigger
%%matrix such that once the inner loops finish running all KG from it at
%%every time should be concanated vertically for every set of B untill the
%%L is exhausted.the final matrix size that KG should for is a 30 by 60
%%matrix of combined KG

Hi experts,

I'm quite new to matlab and I have a problem.im working on an FEA project and I'm required to write a simple code that combines the matrices into the K zeros matrix in a get in a certain order. The smaller matrices for each KG and the code prints them out from the two for loops which are 50 in number correctly.but it is expected to arrange or cantecate those KG matrix vertically once the inner for loop runs and keep that progress until the entire outer loops finishes. I have done Alot of internet search for a few weeks and nothing is working .please help me. here is the code

[Image Analyst]

Try this: http://www.mathworks.com/matlabcentral/answers/13205#answer_18099

[Walter Roberson]

Please attach your .m file. When you posted it, the spacing got messed up terribly.

[rees adah]

Ok...please how do I do that?I am quite new on this platform and this is my first time of asking for assistance.or can I send it to your box so you could help me upload?I just copied the entire m file and pasted it there before the questions.theres nothing more to keep away

[rees adah]

Thanks Walter for editing the question. My problem is how to fill up the K zero matrix with the sets of KG that are output from the loop which are a 6 by 6 matrix such that for every inner loop they are cancanated vertically into the K zero matrix.i was thinking if one could use multidimensional arrays to store up all the KG smaller matrices and then out of the loops do a fresh arrangement buh none of them worked

[rees adah]

Please can I get some help?I was wondering if I could combine all the subsequent matrix into a multidimensional array on different pages and store them there for use later...is there a way I could do that via the inner loop?it didn't work for me

Answer 1

[rees adah]

Please how do I store all the output of the KG matrices(50 of (6* 6)) from the loops into a multidimensional k array?I really need need plz

[Walter Roberson]

MultidimensionalArray(:, :, SetNumber) = TheKGArray

[rees adah]

Thanks if this is what you mean K(:,:,B)=KG This assigns only the last KG array to all the pages of the multidimensional array rather than the individual KG.it still didn't work

[rees adah]

Can someone please help!?

[Walter Roberson]

K(:, :, L-L1+1, B-L1+1) = KG;

[rees adah]

I get the error "subscript indices must be real positive integers or logical"...I really dont know what else to do...can't I store all of the k values into a multidimensional array?

[rees adah]

Please could you help check or run to see it? It's not storing into the K multidimensional array..any help will be appreciated