script

prog='FR AD r1r2x1x2'; 

%          Расчет сопротивлений статора и ротора ADFR
%                  по каталожным данным

%         Каталожные данные двигателя

%Тип двигателя
tip='4MTF(H)200L6' 
pp=3;                  %Число пар полюсов
u1n=220;           %Номинальное фазное напряжение,В
f1n=50;               %Номинальная частота,Гц
Pn=22000;        %Номинальная мощность, Вт
nn=935;            %Номинальная скорость,об/мин  935
I1n=55;             %Номинальный ток статора,А 
kpdn=0.83;       %Номинальный кпд
cosfi1n=0.79;  %Номинальный cos fi1n
Mk=638;            %Критический момент, Нм
I2n=60;             %Номинальный ток ротора, А
E20=235;         %Напряжение на кольцах неподвижного разомкнутого ротора, В 

%Обмоточные данные двигателя [Яуре]

%r1=0.235;         %Невыключаемое активное сопротивление статора,Ом
%r2=0.073;         %Невыключаемое активное сопротивление ротора,Ом
%x1=0.27;           %Невыключаемое индуктивное сопротивление статора,Ом
%x2=0.57/2.37;  %Невыключаемое индуктивное сопротивление ротора,Ом 0.232      

%Расчетные данные двигателя

pi=3.1416;
won=2*pi*f1n/pp; %Номинальная синхронная скорость,рад/с
wn=nn/9.55;      %Номинальная скорость,рад/с
Mn=Pn/wn;        %Номинальный момент,Нм
sn=(won-wn)/won; %Номинальное скольжение
Mko=Mk/Mn;       %Критический момент, о.е.
ke=0.95*u1n/(E20/1.73) ;                                                              % 1.5386

%Приведенный ток ротора
I2np=I2n/ke;                                                                                     %  38.99 A
I2np1=I1n*cosfi1n ;                                                                        % 43.45 A

%Активное сопротивление ротора
%r2 = r2кат20 гр•(1+а(меди)•(t перегр гр – 20гр)) = 0,073•(1+0,0043•(140–20)) = 0,11 Ом,

 r2=E20*sn/(1.73*I2n)                                                          %   0.1473 Ом

%Приведенное среднее сопротивление ротора
 r2p=r2*ke^2;                                                  %=0.153*1.5386^2=0.362 Ом
 
%Потери в роторе
 dP2n=3*I2np^2*r2p;
 
%Электромагнитная мощность номинальная
Pemn=dP2n/sn

%Номинальный электромагнитный момент,
Memn=Pemn/won

%Активное сопротивление статора r1
a1 = 1;
sk=sn*(Mko+sqrt((Mko^2-1)+2*a1*sn*(Mko-1))/(1-2*sn*a1*(Mko-1))); % 0,4159
a=(Memn*(sn/sk+sk/sn)-2*Mk)/(2*sk*(Mk-Memn))                                    % a=0.75
r1=r2p*a 

%Индуктивные сопротивления рассеяния статора и приведенное ротора
    xk=r2p*sqrt(1-a^2*sk^2)/sk;                                                                 %0.8
    x1=xk/2                                                                                                    %0.4
    x2p=xk/2;                                                                                                  %0.4
    x2=x2p/ke^2                                                                                            %0.169
    
%Ток намагничивания в номинальной точке по формуле  В.А. Шубенко
sinfi1n=sqrt(1-cosfi1n^2) ;
Imn=I1n*(sinfi1n-cosfi1n*sn/sk)                                                            %26.92

 %Индуктивное сопротивление контура намагничивания
 
 xmn=u1n/Imn-x1                                                                                      %7.77
 
 %Номинальная ЭДС статора
E1n=Imn*xmn;                                                                                 %209.2/220=0.951*u1n

%Номинальный момент на валу 
Mn=Pn/wn;                                                                                                     %224.7

%Потери момента в номинальном режиме, Нм
dMxn=Memn-Mn                                                                                 %17.66/242.36=0.073*Memn

