clc
clear all
close all
p=40000;
B=3;
Cp=0.42;
eta=0.92;
v=11.5;
rou=1.225;
lamda=6;
i=14;
LB=5.9;
alpha_0=-4.25%需要计算，随机取得一个数据
% p=input('请输入输出功率p=');
% B=input('请输入叶片数B=');
% Cp=input('请输入风能利用系数Cp=');
% eta=input('请输入机电效率eta=');
% v=input('请输入设计风速v=');
% rou=input('请输入空气密度rou=');
% lamda=input('请输入叶尖速比lamda=');
% i=input('请输入叶片分段数i=');
% r=input('请输入机舱半径r=');
% LB=input('请输入叶片长LB=');
D=sqrt(8*p/(Cp*rou*eta*pi*v^3))%风轮直径
R=D/2%风轮半径
r_hub=R-LB
%LB* =R-r_hub;%叶片长度
dbs=LB/i%每段叶片长
for j=1:i
r(j)=r_hub+j*dbs;
lamda_r(j)=lamda*r(j)/R;
theta(j)=atan(lamda_r(j))./3+pi/7
a(j)=sqrt(lamda_r(j).^2+1).*cos(theta(j))
b(j)=sqrt(1+(1-a(j).^2)./lamda_r(j).^2)
end
r%翼型站位
lamda_r%该站位处对应翼型的叶尖速比
theta%中间变量
a%轴向诱导因子
b%周向诱导因子
a=xlsread('4412.xls');
[m,n]=size(a);
dl=max(a(:,4))%最大升阻比的查找
adress=find(a==dl);
cl_adress=adress-3*m;%最大升阻比对应的升力系数位置的查找
cl=a(cl_adress,2)%最大升阻比对应的升力系数
alpha_adress=adress-3*m;%最大升阻比对应的攻角位置的查找
alpha=a(cl_adress,1)%最大升阻比对应的攻角
%升力系数曲线的拟合
% p=polyfit(a(:,1),a(:,2),3)
% y1=polyval(p,a(:,1))
% plot(a(:,1),a(:,2),'k*',a(:,1),y1,'r-')
% plot(a(:,1),a(:,2))
% alpha_0=polyval(p,0)%求零升力攻角
cc=16*pi/9
N=(cc*(LB./r))./(lamda.*(sqrt(lamda^2.*(r./LB).^2+4/9)))%形状参数
L=(r.*N)/(cl*B)%各翼型断面的弦长
L_average=mean(L(:))%叶片的平均弦长
Rz=LB/L_average
K=cl/(alpha-alpha_0)%升力曲线的平均斜率
alpha_m=(L./K).*(1+3/Rz)%校正后的攻角
pha=atan((2/3)*(LB*lamda)./r)%迎风角
theta_i=pha-alpha_m%安装角
a_y=textread('NACA4412.txt');
x3=a_y(:,1);
y3=a_y(:,2);
b=L()
评论1