ini.m

% close all
clear all
clc


%% DFIG Parameters -> Rotor parameters referred to the stator side

f= 50;                      % Stator frequency (Hz)
Ps = 2e6;                   % Rated stator power (W)
n = 1500;                   % Rated rotational speed (rev/min)
Vs = 690;                   % Rated stator voltage (V)
Is = 1760;                  % Rated stator current (A)
Tem = 12732;                % Rated torque (N.m)

p = 2;                      % Pole pair
u = 1/3;                    % Stator/rotor turns ratio
Vr = 2070;                  % Rated rotor voltage (non-reached)(V)
smax = 1/3;                 % Maximum slip
Vr_stator = (Vr*smax)*u;    % Rated rotor voltage referred to stator(v)
Rs = 2.6e-3;
Lsi = 0.087e-3;
Lm = 2.5e-3;
Rr = 2.9e-3;
Ls = Lm + Lsi;
Lr = Lm + Lsi;
Vbus = Vr_stator*sqrt(2);
sigma = 1 - Lm^2/(Ls*Lr);
Fs = Vs*sqrt(2/3)/(2*pi*f);

J = 127;
D = 1e-3;

fsw = 4e3;
Ts = 1/fsw/50;

%PI regulators

tau_i = (sigma*Lr)/Rr;
tau_n = 0.05/4;
wni = 100*(1/tau_i);
wnn = 1/tau_n;

kp_id = (2*wni*sigma*Lr)-Rr;
kp_iq = kp_id;
ki_id = (wni^2)*Lr*sigma;
ki_iq = ki_id;
kp_n = (2*wnn*J)/p;
ki_n = ((wnn^2)*J)/p;

% Three blade wind turbine model

N = 100;                       % Gearbox ratio
Radio = 42;                    % Radio
ro = 1.225;                    % Air density

% Cp and Ct curves

beta = 0;                      % Pitch angle
ind2 = 1;

    for lambda = 0.1:0.01:11.8

        lambdai(ind2) = (1./((1./(lambda-0.02.*beta)+(0.003./(beta^3+1)))));
        Cp(ind2) = 0.73.*(151./lambdai(ind2)-0.58.*beta-0.002.*beta^2.14-13.2).*(exp(-18.4./lambdai(ind2)));
        Ct(ind2) = Cp(ind2)/lambda;
        ind2 = ind2+1;
    end
    tab_lambda = [0.1:0.01:11.8];

% Kopt for MPPT
Cp_max = 0.44;
lambda_opt = 7.2;
Kopt = ((0.5*ro*pi*(Radio^5)*Cp_max)/(lambda_opt^3));

% Power curve in function of wind speed

P = 1.0e+06 *[0,0,0,0,0,0,0,0.0472,0.1097,0.1815,0.2568,0.3418, ...
    0.4437,0.5642,0.7046,0.8667,1.0518,1.2616,1.4976,1.7613,2.0534, ...
    2.3513,2.4024,2.4024,2.4024,2.4024,2.4024, 2.4024];
V = [0.0000, 0.5556,1.1111,1.6667,2.2222,2.7778,3.3333,3.8889,4.4444, ...
    5.0000, 5.5556,6.1111,6.6667,7.2222,7.7778,8.3333,8.8889,9.4444, ...
    10.0000, 10.5556, 11.1111, 11.6667, 12.2222, 12.7778, 13.3333, 13.8889, ...
    14.4444, 15.0000];
figure
subplot(1,2,1)
plot(tab_lambda,Ct,'linewidth',1.5)
xlabel('lambda','fontsize',14)
ylabel('Ct','fontsize',14)
subplot(1,2,2)
plot(V,P,'linewidth',1.5)
grid
xlabel('Wind speed (m/s)','fontsize',14)
ylabel('Power (W)','fontsize',14)