MATLAB Answers

assemble global stiffness matrix

351 views (last 30 days)
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(3,3) and in the end I want it to be like bigK = [ 1 -1 0; -1 2 -1; 0 -1 1]
in which the small matrices will assemble in diagonal way
how I can do it using a for loop ?

  0 Comments

Sign in to comment.

Accepted Answer

John D'Errico
John D'Errico on 12 Oct 2018
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.

  5 Comments

Show 2 older comments
abdelrahman alhammadi
abdelrahman alhammadi on 12 Oct 2018
is there another way of doing it in a for loop in which I only define the small k function and size of big k and then it will do it by itself ?
Torsten
Torsten on 12 Oct 2018
k=[1 -1; -1 1];
n=5;
K=zeros(n+1);
for i=1:n
K(i:i+1,i:i+1)=K(i:i+1,i:i+1)+k;
end
K

Sign in to comment.

More Answers (2)

Torsten
Torsten on 12 Oct 2018
k1 = [1 -1; -1 1];
k2 = k1;
K = zeros(3,3)
K(1:2,1:2) = k1;
K(2:3,2:3) = K(2:3,2:3) + k2;
K

  0 Comments

Sign in to comment.


abdelrahman alhammadi
abdelrahman alhammadi on 12 Oct 2018
what if I have vector a1= [1; -1] and a2= [ 1 ;-1 ] and so on ... the big matrix F= [ 1; 0; 1] if the size was(2) how can I do that in a for loop ?

  1 Comment

John D'Errico
John D'Errico on 12 Oct 2018
Please stop adding multiple answers for every single comment!

Sign in to comment.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!