cb

"lvrt of doubly fed induction matchine by using crowbar" simulation in simulink of matlab that initialization condutions will be uploaded in another file

clear
clc
NumberOfUnits=6;

% Asynchronous machine data
Pnom=1.5e6/0.9*NumberOfUnits; % Nominal power （VA）
Vnom=575;       % Line-Line voltage （Vrms）
Fnom=60;        %Hz
Rs=0.00706;     %pu
Lls=0.171;      %pu
Rr=0.005;       %pu
Llr=0.156;      %pu
Lm=2.9;         %pu
H=5.04;         %Inertia constant （s）
F=0.01;         %Friction factor （pu）
p=3;            %Number of pairs of poles


Ls=Lls+Lm;
Lr=Llr+Lm;
sig=1-（Lm^2/（Ls*Lr））;
%PWM data
PWM_freq=27*Fnom;   %Not used with average model

%DC link
Vdc_nom=1200;       %Volts
C_DClink=10000e-6*NumberOfUnits; %Farads

%Choke data
R_RL=0.30/100;  %pu
L_RL=0.30;      %pu

%Turbine data
Pmec=1.5e6*NumberOfUnits;  % Nominal power （W）
power_C=0.73;   %pu
speed_C=1.2;    %pu

%Tracking characteristic speeds: [speed_A ... speed_D]
speed_ABCD=[0.7 0.71 1.2 1.21];
speed_A=speed_ABCD（1）; %pu of synchronous speed
speed_B=speed_ABCD（2）; %pu of synchronous speed
speed_C=speed_ABCD（3）; %pu of synchronous speed
speed_D=speed_ABCD（4）;

power_A=0; %pu of mechanical power
power_B=power_C*（speed_B/speed_C）^3; %pu of mechanical power
power_D=1;



Kp_grid_side_cur_reg=2.5;
Ki_grid_side_cur_reg=500;
Kp_dc_reg=0.002;
Ki_dc_reg=0.05;

Kp_Q=0.05;
Ki_Q=5;
Kp_rotor_side_cur_reg=0.3;
Ki_rotor_side_cur_reg=8;

%Imax_grid_conv=Pmax_pu/Vnom_pu=Pmax_pu/1;
Pmax=500;
Imax_grid_conv=Pmax;




%Ts_Control=1.000e-04;
%Ts_Power=5.0e-06;
%B=0;

%%%%%??????????????//Pnom=Pmec;
%%
	%Wind turbine typical characteristic data
    c2=116;c3=0.4;c4=5;c5=21;
    %Computation of c1 and c6 for lambda_nom et cpmax
    Pelec_base=Pnom;
    wind_base=11;
    P_wind_base=power_C;
    speed_nom=speed_C;
    pitch_angle=0;
    cp_nom=0.48;
    lambda_nom=8.1;
    %cp_nom=0.4382;
    %lambda_nom=6.325;
    cp_max=cp_nom;

    lambda_i_nom=1/（1/lambda_nom-0.035）;
    k_nom=-（c2*c5/lambda_i_nom-c4*c5-c2）*exp（-c5/lambda_i_nom）/lambda_nom^2;
    c1=cp_max/（（c2/lambda_i_nom-c4）*exp（-c5/lambda_i_nom）+k_nom*lambda_nom）;
    c6=k_nom*c1;
    

	wr_elec_pu=1e-3:0.01:1.5;
    wr_mec_pu=wr_elec_pu/speed_nom;
    wind_max=ceil（（1.2/P_wind_base）^（1/3）*wind_base）;
    wind_min=round（（（0.1）/P_wind_base）^（1/3）*wind_base）;
    wind1=wind_min:（wind_base-wind_min）/5:wind_base;
    wind2=wind_base+（wind_base-wind_min）/5:（wind_base-wind_min）/5:wind_max;
    wind_pu=[wind1 wind2]/wind_base;
    
   	fig_num=figure;
	clf（fig_num）; 
    for k=1:length（wind_pu）
        for i=1:length（wr_mec_pu）;
            lambda_pu=wr_mec_pu（i）/wind_pu（k）;
            lambda=lambda_pu*lambda_nom;
            lambda_i=1/（1/（lambda+0.08*pitch_angle）-0.035/（pitch_angle^3+1））;
            cp=c1*（c2/lambda_i-c3*pitch_angle-c4）*exp（-c5/lambda_i）+c6*lambda;
            cp_pu=cp/cp_nom;
            Pwind_pu（ki）=wind_pu（k）^3*cp_pu*P_wind_base;
        end
    end
    for k=1:length（wind_pu）
        plot（wr_elec_puPwind_pu（k:））
        hold on
    end
    
    plot（[speed_Dspeed_D0

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

     文件       5093  2010-09-27 10:02  Initialisation.m

     文件     836552  2010-11-12 17:29  IMC_cb.mdl

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

               841645                    2