Creating Global Stiffness matrix

59 views (last 30 days)
Alpha Boy
Alpha Boy on 4 Jul 2020
Commented: Ruairi Kavanagh on 23 Apr 2023
Hello,
I am quite noob in Matlab. I am trying to add two matrix to create a global stiffness matrx.
k1 =
10150000 -5000000 -3750000 2000000 -6400000 3000000
-5000000 12400000 3000000 -10000000 2000000 -2400000
-3750000 3000000 3750000 0 0 -3000000
2000000 -10000000 0 10000000 -2000000 0
-6400000 2000000 0 -2000000 6400000 0
3000000 -2400000 -3000000 0 0 2400000
k2 =
10150000 -5000000 -3750000 2000000 -6400000 3000000
-5000000 12400000 3000000 -10000000 2000000 -2400000
-3750000 3000000 3750000 0 0 -3000000
2000000 -10000000 0 10000000 -2000000 0
-6400000 2000000 0 -2000000 6400000 0
3000000 -2400000 -3000000 0 0 2400000
I will find a global stiffness matrix(K) (8X8). How can I find this matrix?
  5 Comments
Daniyal Raja
Daniyal Raja on 24 Nov 2022
Moved: the cyclist on 24 Nov 2022
sem question my g

Sign in to comment.

Answers (1)

the cyclist
the cyclist on 5 Jul 2020
If I've understood correctly, then
K = zeros(8);
vec1 = [1 2 3 4 7 8];
K(vec1,vec1) = K(vec1,vec1) + k1;
vec2 = [5 6 7 8 3 4];
K(vec2,vec2) = K(vec2,vec2) + k2;
gives what you want.
Because K starts out as all zeros, the addition with k1 could have been written
K(vec1,vec1) = k1;
but I thought it might be easier to see what was going on this way, and also that k1 and k2 are being treated symmetically.
  4 Comments
the cyclist
the cyclist on 5 Jul 2020
I'm assuming you've defined k1 and k2, such as
k1 =[
10150000 -5000000 -3750000 2000000 -6400000 3000000
-5000000 12400000 3000000 -10000000 2000000 -2400000
-3750000 3000000 3750000 0 0 -3000000
2000000 -10000000 0 10000000 -2000000 0
-6400000 2000000 0 -2000000 6400000 0
3000000 -2400000 -3000000 0 0 2400000];
k2 = [
10150000 -5000000 -3750000 2000000 -6400000 3000000
-5000000 12400000 3000000 -10000000 2000000 -2400000
-3750000 3000000 3750000 0 0 -3000000
2000000 -10000000 0 10000000 -2000000 0
-6400000 2000000 0 -2000000 6400000 0
3000000 -2400000 -3000000 0 0 2400000];
Ruairi Kavanagh
Ruairi Kavanagh on 23 Apr 2023
how would i would i change this code if i had 8 matrices fo 4x4 that needed to fit into a 6x6 matrix

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!