Advertisement
basyair7

ini.m

May 9th, 2022
576
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. % close all
  2. clear all
  3. clc
  4.  
  5.  
  6. %% DFIG Parameters -> Rotor parameters referred to the stator side
  7.  
  8. f= 50;                      % Stator frequency (Hz)
  9. Ps = 2e6;                   % Rated stator power (W)
  10. n = 1500;                   % Rated rotational speed (rev/min)
  11. Vs = 690;                   % Rated stator voltage (V)
  12. Is = 1760;                  % Rated stator current (A)
  13. Tem = 12732;                % Rated torque (N.m)
  14.  
  15. p = 2;                      % Pole pair
  16. u = 1/3;                    % Stator/rotor turns ratio
  17. Vr = 2070;                  % Rated rotor voltage (non-reached)(V)
  18. smax = 1/3;                 % Maximum slip
  19. Vr_stator = (Vr*smax)*u;    % Rated rotor voltage referred to stator(v)
  20. Rs = 2.6e-3;
  21. Lsi = 0.087e-3;
  22. Lm = 2.5e-3;
  23. Rr = 2.9e-3;
  24. Ls = Lm + Lsi;
  25. Lr = Lm + Lsi;
  26. Vbus = Vr_stator*sqrt(2);
  27. sigma = 1 - Lm^2/(Ls*Lr);
  28. Fs = Vs*sqrt(2/3)/(2*pi*f);
  29.  
  30. J = 127;
  31. D = 1e-3;
  32.  
  33. fsw = 4e3;
  34. Ts = 1/fsw/50;
  35.  
  36. %PI regulators
  37.  
  38. tau_i = (sigma*Lr)/Rr;
  39. tau_n = 0.05/4;
  40. wni = 100*(1/tau_i);
  41. wnn = 1/tau_n;
  42.  
  43. kp_id = (2*wni*sigma*Lr)-Rr;
  44. kp_iq = kp_id;
  45. ki_id = (wni^2)*Lr*sigma;
  46. ki_iq = ki_id;
  47. kp_n = (2*wnn*J)/p;
  48. ki_n = ((wnn^2)*J)/p;
  49.  
  50. % Three blade wind turbine model
  51.  
  52. N = 100;                       % Gearbox ratio
  53. Radio = 42;                    % Radio
  54. ro = 1.225;                    % Air density
  55.  
  56. % Cp and Ct curves
  57.  
  58. beta = 0;                      % Pitch angle
  59. ind2 = 1;
  60.  
  61.     for lambda = 0.1:0.01:11.8
  62.        
  63.         lambdai(ind2) = (1./((1./(lambda-0.02.*beta)+(0.003./(beta^3+1)))));
  64.         Cp(ind2) = 0.73.*(151./lambdai(ind2)-0.58.*beta-0.002.*beta^2.14-13.2).*(exp(-18.4./lambdai(ind2)));
  65.         Ct(ind2) = Cp(ind2)/lambda;
  66.         ind2 = ind2+1;
  67.     end
  68.     tab_lambda = [0.1:0.01:11.8];
  69.    
  70. % Kopt for MPPT
  71. Cp_max = 0.44;
  72. lambda_opt = 7.2;
  73. Kopt = ((0.5*ro*pi*(Radio^5)*Cp_max)/(lambda_opt^3));
  74.  
  75. % Power curve in function of wind speed
  76.  
  77. P = 1.0e+06 *[0,0,0,0,0,0,0,0.0472,0.1097,0.1815,0.2568,0.3418, ...
  78.     0.4437,0.5642,0.7046,0.8667,1.0518,1.2616,1.4976,1.7613,2.0534, ...
  79.     2.3513,2.4024,2.4024,2.4024,2.4024,2.4024, 2.4024];
  80. V = [0.0000, 0.5556,1.1111,1.6667,2.2222,2.7778,3.3333,3.8889,4.4444, ...
  81.     5.0000, 5.5556,6.1111,6.6667,7.2222,7.7778,8.3333,8.8889,9.4444, ...
  82.     10.0000, 10.5556, 11.1111, 11.6667, 12.2222, 12.7778, 13.3333, 13.8889, ...
  83.     14.4444, 15.0000];
  84. figure
  85. subplot(1,2,1)
  86. plot(tab_lambda,Ct,'linewidth',1.5)
  87. xlabel('lambda','fontsize',14)
  88. ylabel('Ct','fontsize',14)
  89. subplot(1,2,2)
  90. plot(V,P,'linewidth',1.5)
  91. grid
  92. xlabel('Wind speed (m/s)','fontsize',14)
  93. ylabel('Power (W)','fontsize',14)
  94.  
  95.  
Advertisement
RAW Paste Data Copied
Advertisement