风力机叶片设计MATLAB程序

对于选定翼型，利用MATLAB计算其升力系数，基于叶素栋梁理论计算风力机性能，优化设计长度归一化的风力机叶片模型。

clear;
syms sita;
alfa=3*pi/180;
y = 0.5*cos（sita）;
% pu=0.043;qu=0.081;au=0.61;bu=1.6;cu=0.48;du=1.1;
pu=0.01;qu=0.05;au=1;bu=1;cu=1;du=0.6;
% zu = pu*（1 + cos（sita））^au *（1 - cos（sita））^bu +qu*（1 + cos（sita））^cu*（1 - cos（sita））^du;
zu = pu*（sin（sita））^au *（sin（sita））^bu +qu*（sin（sita））^cu*（1 - cos（sita））^du;
dzu=diff（zu）;
% pl=0.022;ql=0.130;al=2.6;bl=1;cl=0.7;dl=1.7;
pl=0.022;ql=0.130;al=2.6;bl=1;cl=0.7;dl=1.7;
zl = pl*（1 + cos（sita））^al* （1 -cos（sita））^bl-ql*（1 +cos（sita））^cl*（1 -cos（sita））^dl;
dzl=diff（zl）;
Vu=（sin（sita）+2*zu）^2*（sin（sita-alfa）+sin（alfa））^2/（4*dzu^2*sin（sita）^2+sin（sita）^4）;
Vl=（sin（sita）+2*zl）^2*（sin（sita-alfa）+sin（alfa））^2/（4*dzl^2*sin（sita）^2+sin（sita）^4）;
Cpu=1-Vu;
Cpl=1-Vl;
dCLu=Cpu*（dzu*sin（alfa）-sin（sita）*cos（alfa）/2）;
dCLl=Cpl*（dzl*sin（alfa）-sin（sita）*cos（alfa）/2）;
% CLu=integral（@（sita）dCLu0pi）;CLl=integral（@（sita）dCLlpi2*pi）;
% CL=CLu+CLl;
CLu=0;CLl=0;
% for sita=0:0.001:pi
%     CLu=CLu+0.001*subs（dCLu）;
%     CLl=CLl+0.001*subs（dCLl）;
% end
% CL=CLu+CLl;
ezplot（yzu）hold on;
grid on;
% ezplot（sitaCpl）hold on;
% ezplot（sitadCLu）hold on;
% ezplot（sitadCLl）hold on;

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

     文件       1154  2017-12-14 15:27  m文件\airfoil_752.m

     文件       5981  2017-12-22 11:12  m文件\blade_752_20171219.m

     文件       2832  2017-12-15 10:37  m文件\blade_752_performance.m

     目录          0  2018-05-22 14:54  m文件

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

                 9967                    4