mesh2.m

2D MATLAB EFG CODE - SQUARE DOMAIN

%-------------------------------------------------------------------------
% MESH GENERATION PROGRAM FOR 2D BEAM IN BENDING
function[x,conn,numcell,numq] = mesh2(length,height,ndivl,ndivw)

% INPUT DATA
numcell= ndivw*ndivl;
numq = (ndivl+1)*(ndivw+1); 

% SET UP NODAL COORDINATES
for i = 1:(ndivl+1)
   for j = 1:(ndivw+1)
   x(1,((ndivw+1)*(i-1) +j))= (length/ndivl)*(i-1);
   x(2,((ndivw+1)*(i-1) +j))= -(height/ndivw)*(j-1)+height/2;
   end
end

% SET UP CONNECTIVITY ARRAY
for j = 1:ndivl
   for i = 1:ndivw
   elemn = (j-1)*ndivw + i;
   nodet(elemn,1) = elemn + (j-1);
   nodet(elemn,2) = nodet(elemn,1) + 1;
   nodet(elemn,3) = nodet(elemn,2)+ndivw+1;
   nodet(elemn,4) = nodet(elemn,3)-1;
   end
end
conn = nodet';