function AssembleEleStif( ie, Ke )          % 3级子程序，整体定�? 
%% 单元刚度矩阵拼装
global Elements gKA 

ENodes( : , 1 ) = Elements( ie , : ); 
enm = size( ENodes , 1 );      % 单元的结点�?�数 
iLR = 2 * enm;              % 单元的�?�自由度�? 
viDof = zeros( 2 * enm , 1 );   % 各结点自由度编号 

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