script

prog='rxKZAD'

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

tip='4MTKF(H)200L6';  %Тип двигателя
pp=3;                  %Число пар полюсов
u1n=220;           %Номинальное фазное напряжение,В
f1n=50;               %Номинальная частота,Гц
Pn=22000        %Номинальная мощность, Вт
nn=935;            %Номинальная скорость,об/мин
I1n=51;             %Номинальный ток статора,А 
kpdn=0.83;       %Номинальный кпд
cosfi1n=0.79;  %Номинальный cos fi1n
Mk=760;           %Критический момент, Нм
Mp=706;          %Пусковой момент,Нм
I1p=275;          %Пусковой ток статора, A

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

pi=3.1416;
won=2*pi*f1n/pp; %Номинальная синхронная скорость,рад/с
wn=nn/9.55;      %Номинальная скорость,рад/с
Mn=Pn/wn;        %Номинальный момент,Нм
sn=(won-wn)/won; %Номинальное скольжение
Mko=Mk/Mn;       %Критический момент, о.е.
Mpo=Mp/Mn;       %Пусковой момент, о.е.
I1po=I1p/I1n;    %Пусковой ток статора, о.е.

%Приведенное активное сопротивление ротора r2p

%Если предположить, что полные потери состоят из постоянных и переменных потерь,
%и постоянные примерно равны 1/3 полных потерь,
%а механические потери составляют половину постоянных потерь, 
%то механические потери dРмех определяются из уравнения

 dPmex=Pn*(1/kpdn-1)/6;
 
 Pmex=Pn+dPmex;
 
 r2p= Pmex*sn/(3*(1-sn)*(I1n*cosfi1n)^2);

%Активное сопротивление фазной обмотки статора r1 

dP2n=3*(I1n*cosfi1n)^2*r2p;
Pemn=dP2n/sn;
Memn=Pemn/won;
a1 = 1;
sk=sn*(Mko+sqrt((Mko^2-1)+2*a1*sn*(Mko-1))/(1-2*sn*a1*(Mko-1))); 
a=(Memn*(sn/sk+sk/sn)-2*Mk)/(2*sk*(Mk-Memn)) ;                            
r1=r2p*a ;                                                                                                      

%Индуктивные сопротивления рассеяния статора и приведенное ротора
    xk=r2p*sqrt(1-a^2*sk^2)/sk;
    x1=xk/2;
    x2=xk/2;

%Ток намагничивания в номинальной точке по формуле  В.А. Шубенко
sinfi1n=sqrt(1-cosfi1n^2);
Imn=I1n*(sinfi1n-cosfi1n*sn/sk);

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

%Номинальный электромагнитный момент,
Memn=2*Mk*(1+a*sk)/(sn/sk+sk/sn+2*a*sk);

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

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

Imn
xmn
 x1
x2
r1
r2p









