（1）网格离散程序 grid_generation_juxing.m
clc
clear
clf
%%%%%%% 区域几何尺寸及网格划分参数
H=0.06; %区域总高
L=0.08; %区域总长
Nx=4; %水平方向的网格数量，选择能被 5 整除的数
Ny=4; %竖直方向的网格数量
theta=0; %网格平面旋转角度
%%%%%%% 区域几何尺寸及网格划分参数
%%%%%%% 总单元数和结点数
E=Nx*Ny; %总单元数
Nz=(2*Nx+1)*(2*Ny+1); % 二次单元结点总数
Nd=(Nx+1)*(Ny+1); % 线性单元结点总数
%%%%%%% 总单元数和结点数
%%%%%%% 单元间距
Dx=L/Nx/2; %水平方向网格间距
Dy=H/Ny/2; %竖直方向网格间距
%%%%%%% 单元间距
%%%%%%% 结点分布拓扑
AAA=zeros(Ny*2+1,Nx*2+1);
for i=1:2:2*Nx+1
AAA(1,i)=(i+1)/2;
end
for i=1:Nx
AAA(1,2*i)=(Nx+1)*(Ny+1)+i;
end
for i=1:2*Nx+1
AAA(2,i)=(Nx+1)*(Ny+1)+Nx+i;
评论0