bem_bing_liu.m

叶片元素计算法 基于叶素-动量理论（BEM）. 采用叶素-动量理论可以计算风轮旋转面中的轴向诱导。 风力发电用

clear
clc

v_0=[5:1:25]
rho=1.225
R=30.5
B=3
theta=-4
rpm_max=21.92
omega=(2*pi/60)*rpm_max



     
[a,a_prime,c_r,c_n,phi,r]= function_bem(v_0, theta,omega);



for i=1: length(v_0)
     T(i)=0
     P(i)=0
    for j=1:(length(r)-1)

         T(i) = T(i) + 1/2*rho*B*c_r(j)*c_n(j,i)*v_0(i)^2*(1-a(j,i))^2/(sin(phi(j,i)*pi/180))^2 *(r(j+1)-r(j)); % (equation 8.21)
         P(i) = P(i) + 4*pi*rho*omega^2*v_0(i)*a_prime(j,i)*(1-a(j,i))*r(j)^3*(r(j+1)-r(j)); % (equation 6.29)
    end
         lambda(i) = omega*R/v_0(i);
         c_p(i) = P(i)/(1/2*rho*v_0(i)^3*pi*R^2); % (equation 6.20)
end
 


 figure(1)
 plot(v_0,T,'r'), xlabel('v_0   (m/s)'), ylabel('T  (N)'),title('---------------')
 figure(2)
 plot(v_0,P,'b'), xlabel('v_0   (m/s)'), ylabel('P  (W)') ,title('--------------')
 figure(3)
 plot(lambda,c_p,'k'), xlabel('\lambda'), ylabel('c_p'),title('--------------')
  figure(4)
 plot(v_0,c_p,'k'), xlabel('\V0'), ylabel('c_p'),title('--------------')