matlab开发-使用xFoiland ParseCGeometric参数化进行机翼优化

matlab开发-使用xFoiland ParseCGeometric参数化进行机翼优化。包括使用xfoil接口优化翼型所需的所有文件。

%% Airfoil Optimization objective Function
% Authors: Sean Wu & Sheida Hosseini UC Davis MAE Dept. 
% Date:   Mar. 2015

% CC BY-SA Sean Wu 2015

%{ 

This function returns the drag coefficient for a given airfoil as defined
by the PARSEC parameters p=[p1...p11] using Xfoil for the aerodynamic
analysis.

%}

function cd = cdmin（p）

a=parsec（p）;

point=100;
n=（1/point）;
nn=n/5; % to create more points at the LE

% upper surface coordinates
x1=1:-n:0.1;
x2=0.1:-nn:0;

% lower surface coordinates
%xxx=0:n:1;
x3=0:nn:0.1;
x4=0.1:n:1;

xu（1:length（x1））=x1;
xu（length（x1）+1:length（x1）+length（x2））=x2;

xl（1:length（x3））=x3;
xl（length（x3）+1:length（x3）+length（x4））=x4;

X（1:length（xu））=xu;
X（length（xu）+1:length（xu）+length（xl））=xl;

X=X‘;

yu=a（1）*xu.^.5+a（2）*xu.^（1.5）+a（3）*xu.^（2.5）+a（4）*xu.^（3.5）+a（5）*xu.^（4.5）+a（6）*xu.^（5.5）;
yl=-（a（7）*xl.^.5+a（8）*xl.^（1.5）+a（9）*xl.^（2.5）+a（10）*xl.^（3.5）+a（11）*xl.^（4.5）+a（12）*xl.^（5.5））;

Y（1:length（yu）1）=yu;
Y（length（yu）+1:length（yu）+length（yl）1）=yl;


Y=round（Y7）;
X=round（X7）;
coord=[X Y];

% Run Xfoil requires xfoil.m
[ab]=xfoil（coord5E60.5）;

% Plot airfoil evolution comment out if desired
hold on
plot（XY）
        
% objective function
cd=a.CD; 

end

%Sean Wu 2015

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1293  2015-03-22 13:42  airfoil optim\cdmin.m
     文件         606  2015-03-22 13:42  airfoil optim\central.m
     文件        2908  2015-03-22 13:42  airfoil optim\conjugate.m
     文件         534  2015-03-22 13:42  airfoil optim\fwd.m
     文件         959  2015-03-22 13:42  airfoil optim\golden.m
     文件        1691  2015-03-22 13:42  airfoil optim\parsec.m
     文件       11351  2015-03-22 13:42  airfoil optim\xfoil.m
     文件        1337  2015-03-22 13:42  license.txt