www.pudn.com > zhichixiangliangjijiemian.rar > MESH2.M

% MESH GENERATION PROGRAM FOR 2D BEAM IN BENDING 
function[x,conn,numcell,numq] = mesh2(length,height,ndivl,ndivw) 
ndivw=20; 
ndivl=8; 
length=2; 
height=5; 
 
% 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)-length/2; 
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'; 
conn