FEM-Course-Matlab/10.四节点八节点四边形单元悬臂梁的Matlab有限元编程/8_nodes/AssembleEleStif.m

function AssembleEleStif( ie, Ke )          % 3级å­<C3A5>程åº<C3A5>，整体定ä½? 
%% å<>•å…ƒåˆšåº¦çŸ©é˜µæ‹¼è£…
global Elements gKA 

ENodes( : , 1 ) = Elements( ie , : ); 
enm = size( ENodes , 1 );      % å<>•å…ƒçš„结点æ?»æ•° 
iLR = 2 * enm;              % å<>•å…ƒçš„æ?»è‡ªç”±åº¦æ•? 
viDof = zeros( 2 * enm , 1 );   % å<>„结点自由度编å<E28093>· 

for il=1:enm 
  inod = ENodes( il ); 
  viDof( 2 * il-1 : 2 * il) = 2 * inod-1 : 2 * inod; 
end
% 拼装：ir为行å±?部编å<E28093>?, irg为对应整体编å<E28093>? 
for ir = 1 : iLR 
  irg = viDof( ir ); 
  % ic为列å±?部编å<E28093>?, icg为对应整体编å<E28093>? 
  for ic = 1 : iLR 
    icg = viDof( ic ); 
    gKA( irg , icg ) = gKA( irg , icg ) + Ke( ir , ic ); 
  end 
end