MATLAB: Assemble global stiffness matrix

#assemble

I am trying to assemble a matrix in which I would like to put small matrices into a big one but I don't know how to do that,

I know that I have to define the big matrix first and then show put the small ones inside

for example I have k1=[1 -1; -1 1] , K2=[1 -1; -1 1] and the define bigK in which bigK=zeros(5,5) and in the end I want it to be like bigK = [ 1 0 -1 0 0; -1 0 1 0 0; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0]

in which the small matrices will assemble in diagonal way

how I can do it using a for loop ?

Best Answer

result=blkdiag(k1,K2)

Related Solutions

MATLAB: Assemble global stiffness matrix

k1 = [120 -120;-120 120];
k = zeros(5,5);
k([1 3], [1 3]) = k1;

Result

k =
   120     0  -120     0     0
     0     0     0     0     0
  -120     0   120     0     0
     0     0     0     0     0
     0     0     0     0     0

MATLAB: Assemble global stiffness matrix

K = zeros(10,10);
k = [1 -1; -1 1];
for n = 1:9;
  rind = n + [0 1];
  cind = rind;
  K(rind, cind) = K(rind, cind) + k;
end

There are better ways to do this in terms of pure efficiency. But the above will suffice. First learn how to write at least decent code that solves the task at hand. Only when you decide if the code is a problem do you want to worry about optimizing it for speed.

Best Answer

Related Solutions

MATLAB: Assemble global stiffness matrix

MATLAB: Assemble global stiffness matrix

Related Question