Computing Assembled Stiffness Matrix - Code

Question 1

for in = 1:nnode;

nodei = ?

Answer 1

lnods(ielem, in);

Question 2

for il = 1:ndofn;
ie = ?
ig = ?

Answer 2

(in - 1) * ndofn +il;

(nodei - 1) * ndofn + il;

Question 3

for jn = 1:nnode;

nodej = ?

Answer 3

nodej = lnods(ielem, jn);

Question 4

for jl = 1:ndofn;
je = ?
jg = ?
gstif(ig, jg) = ?

Answer 4

je = (jn - 1) * ndofn + jl;
jg = (nodej - 1) * ndofn + jl;
gstif(ig, jg) = gstif(ig, jg) + estif(ie, je);