结构化网格生成程序

对去边四边形区域，求解泊松方程，实现结构画网格划分，用于CAE分析。matlab程序。

function [xy_Lxy_Rxy_Uxy_D] = linear_boundary（x_maxx_miny_maxy_minnum_node1num_node2n1n2）
% 构造线性矩形边界
% 主要用到线性插值函数linspace
% INPUT:
%      x_max | x_min | y_max | y_min 分别是矩形外边界的范围，注意要包含翼型
%      num_node1 矩形边界x方向节点数，num_node2为y向节点数
%      n1 为矩形下边界节点数，n2为矩形左边界节点数
% OUTPUT：
%      依次为  【左，右，上，下】边界坐标，每个均为2列（包含xy 方向）的列向量
%% 导入翼型坐标并处理
xy0 = importdata（‘CH10 （smoothed）.txt‘）; % 该文件由profili软件得到，去除里面的文字，只保留点坐标
% plot（xy0（:1）xy0（:2））;
num_xy0 = size（xy01）; % 应该为偶数上下面点数相同

a = xy0（ 2 : num_xy0/2  : ）; % 处理数据格式
b = xy0（ num_xy0/2+1 : end :）;
b = flipud（b）;
xy0 = [b;a];                   % 得到逆时针排序的翼型点坐标
% plot（xy0（:1）xy0（:2））;
num_xy0 = size（xy01）;
clear a b;
%% 构造边界
% 构造下边界xy_D
% num_node1 = 89;
n_temp = （num_node1 - num_xy0）/2 ;
a_temp = zeros（n_temp+12）;
b_temp = a_temp;
a_temp（:1） = linspace（x_max100n_temp+1）;
a_temp（:2） = 0;
b_temp（:1） = linspace（100x_maxn_temp+1）;
b_temp（:2） = 0;
xy_D = [a_temp;xy0（2:end-1:）;b_temp]; 
clear a_temp b_temp xy0;                   % 得到xy_D
% 构造上边界xy_U
xy_U = xy_D;
% n1 = 39;
% n2 = 10;
n3 = num_node1 - n1 - n2;
xy_U（1:n11） = linspace（x_maxx_minn1）;
xy_U（1:n12） = y_min;

xy_U（n1:n1+n22） = linspace（y_miny_maxn2+1）;
xy_U（n1:n1+n21） = x_min;

xy_U（n1+n2:end1） = linspace（x_minx_maxn3+1）;
xy_U（n1+n2:end2） = y_max;
% 构造左边界xy_L
% num_node2 = 6;
xy_L = zeros（num_node22）;
xy_L（:2） = linspace（y_min0num_node2）;
xy_L（:1） = x_max;
% 构造右边界
xy_R = xy_L;
xy_R（:2） = linspace（y_max0num_node2）;
xy_R（:1） = x_max;
clear num_node1 num_node2 ;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1533  2010-05-08 15:33  划分网格\CH10 （smoothed）.txt

     文件       1771  2010-05-09 00:07  划分网格\linear_boundary.m

     文件       1301  2010-05-09 10:43  划分网格\Main_proc.m

     文件        127  2010-05-09 09:00  划分网格\plot_boundary.m

     文件       2081  2010-05-09 00:12  划分网格\solve_4_3_13.m

     目录          0  2014-04-11 17:10  划分网格

----------- ---------  ---------- -----  ----

                 6813                    6