90 lines
2.9 KiB
Matlab
90 lines
2.9 KiB
Matlab
function [coordinates, nodes, nP, nE, dxv, dyv,dzv,zm,vm,xrao,yrao,zrao] = GenerateNode_final_new(dx, dy,dz ,nx,ny,nz,NTG)
|
|
% coordinates xy cordinate，即每个节点的x,y坐标
|
|
% nodes connection，即每个子区域对应的八个节点数
|
|
% nP numbers of point，即节点数量
|
|
% nE numbers of element，即子区域数量
|
|
nP = (nx + 1) * (ny + 1)*(nz+1);%计算节点数量
|
|
nE = nx * ny*nz;%子区域数量计算
|
|
% dx=dx*ones(nx,1);%此处是为了暂时妥协
|
|
% dy=dy*ones(nx,1);
|
|
% dz=dz*ones(nx,1);
|
|
dxv = zeros(nE, 1);%每个子区域x方向上尺寸矩阵
|
|
dyv = zeros(nE, 1);%每个子区域y方向上尺寸矩阵
|
|
dzv = zeros(nE, 1);
|
|
zm=zeros(nE, 1);%基质网格计算深度（以下平面为基准计算得到的高度，值为正数）
|
|
%给每个子区域的x方向和y方向赋尺寸
|
|
vm=zeros(nE, 1);%基质网格体积
|
|
for k=1:nz
|
|
for i = 1 : nx
|
|
for j = 1 : ny
|
|
dxv(i+(j-1)*nx+nx*ny*(k-1)) = dx(i);
|
|
dyv(i+(j-1)*nx+nx*ny*(k-1)) = dy(j);
|
|
dzv(i+(j-1)*nx+nx*ny*(k-1))=dz(k);
|
|
vm(i+(j-1)*nx+nx*ny*(k-1))=dx(i)*dy(j)*dz(k)*NTG(i+(j-1)*nx+nx*ny*(k-1));
|
|
|
|
end
|
|
end
|
|
end
|
|
coordinates = zeros(nP, 3);%初始化每个节点的x,y坐标矩阵
|
|
nodes = zeros(nE, 8);%初始化每个子区域对应的八个节点数的矩阵
|
|
%计算每个节点的坐标
|
|
x0 = zeros(nx+1,1);
|
|
y0 = zeros(ny+1,1);
|
|
z0=zeros(nz+1,1);
|
|
sumx = 0;
|
|
sumy = 0;
|
|
sumz=0;
|
|
for i = 1 : nx
|
|
sumx = sumx+dx(i);
|
|
x0(i+1) = sumx;
|
|
end
|
|
for i = 1 : ny
|
|
sumy = sumy+dy(i);
|
|
y0(i+1) = sumy;
|
|
end
|
|
for i = 1 : nz
|
|
sumz = sumz+dz(i);
|
|
z0(i+1) = sumz;
|
|
end
|
|
for k=1:nz+1
|
|
for j = 1 : ny + 1
|
|
for i = 1 : nx + 1
|
|
coordinates(i + (j - 1) * (nx + 1)+(nx+1)*(ny+1)*(k-1), 1) = x0(i);%(i,j）节点的x坐标赋值
|
|
coordinates(i + (j - 1) * (nx + 1)+(nx+1)*(ny+1)*(k-1), 2) = y0(j);%(i,j）节点的y坐标赋值
|
|
coordinates(i + (j - 1) * (nx + 1)+(nx+1)*(ny+1)*(k-1), 3) = z0(k);
|
|
end
|
|
end
|
|
end
|
|
%计算每个子区域的八个顶点(节点）对应的节点编号
|
|
for k=1:nz
|
|
for j = 1 : ny
|
|
for i = 1 : nx
|
|
nodes(i + (j - 1) * nx+nx*ny*(k-1), 1) = i + (j - 1) * (nx + 1)+(k-1)*(nx+1)*(ny+1);
|
|
nodes(i + (j - 1) * nx+nx*ny*(k-1), 2) = nodes(i + (j - 1) * nx+nx*ny*(k-1), 1) + 1;
|
|
nodes(i + (j - 1) * nx+nx*ny*(k-1), 3) = nodes(i + (j - 1) * nx+nx*ny*(k-1), 1) + nx + 1;
|
|
nodes(i + (j - 1) * nx+nx*ny*(k-1), 4) = nodes(i + (j - 1) * nx+nx*ny*(k-1), 1) + nx + 2;
|
|
nodes(i + (j - 1) * nx+nx*ny*(k-1), 5) = nodes(i + (j - 1) * nx+nx*ny*(k-1), 1)+(nx+1)*(ny+1);
|
|
nodes(i + (j - 1) * nx+nx*ny*(k-1), 6) = nodes(i + (j - 1) * nx+nx*ny*(k-1), 5) + 1;
|
|
nodes(i + (j - 1) * nx+nx*ny*(k-1), 7) = nodes(i + (j - 1) * nx+nx*ny*(k-1), 5) + nx + 1;
|
|
nodes(i + (j - 1) * nx+nx*ny*(k-1), 8) = nodes(i + (j - 1) * nx+nx*ny*(k-1), 5) + nx + 2;
|
|
zm(i+(j-1)*nx+nx*ny*(k-1))=sum(coordinates(nodes(i + (j - 1) * nx+nx*ny*(k-1), :),3))/8;
|
|
end
|
|
end
|
|
end
|
|
|
|
xrao=[0];
|
|
for i=1:nx
|
|
xxiang=sum(dx(1,1:i));
|
|
xrao=[xrao xxiang];
|
|
end
|
|
yrao=[0];
|
|
for i=1:ny
|
|
yxiang=sum(dy(1,1:i));
|
|
yrao=[yrao yxiang];
|
|
end
|
|
zrao=[0];
|
|
for i=1:nz
|
|
zxiang=sum(dz(1,1:i));
|
|
zrao=[zrao zxiang];
|
|
end
|
|
end |