assemblestiffnessmatrix.m

是个好东西

function AssembleStiffnessMatrix( ie, k )
% this function assemble stiffness matrix of element
global element K
for i=1:1:2
   for j=1:1:2
      for p=1:1:3
         for q=1:1:3
            m = (i-1)*2+p ;
            n = (j-1)*2+q ;
            M = (element(ie,i)-1)*3+p ;
            N = (element(ie,j)-1)*3+q ;
            K(M,N) = K(M,N) + k(m,n) ;
         end
      end
   end
end
return