Додаток б

Розрахунок уточнених швидкісних характеристик за допомогою пакету «MatLab

clc

close all

figure(1)

%скоростные характеристики%

%паспортные данные%

w0=62.8; Rrot=0.038; Xrot=0.078; Ro=0.042;Xo=0.107; Ke=1.12; Pn=80; nn=580; Em=320;

Im=155;

Mn=30*10^3*Pn/(pi*nn)

Mk=2.8*Mn

%Расчестные пусковые сопротивления%

R1=0.212; R2=0.150; R3=0.112; R4=0.087; R5=0.062; R6=0.043; R7=0.025;

w=linspace(-w0,w0,1000);

S=(w0-w)/w0;

%приведенные к цепи статора сопротивлдения ротора%

Xr=Xrot*Ke^2;

Rr=Rrot*Ke^2;

Xk=Xo+Xr;

a=Ro/Rr

%естественная скоростная характеристика%

Ske=Rr/sqrt(Ro^2+Xk^2)

M=2*Mk*(1+a*Ske)./(S./Ske+Ske./S+2*a*Ske);

I2=sqrt(M*w0.*S/(3*Rr));

plot(I2,w)

grid on ,

hold on

%реостатная скоростная харакетристика №1%

Rrotsum1=Rr+R1+R2+R3+R4+R5+R6+R7;

Rrsum1=Rrotsum1*Ke^2;

a=Ro/Rrsum1

Skl1=Ske*Rrsum1/Rr

M1=2*Mk*(1+a*Ske)./(S./Skl1+Skl1./S+2*a*Ske);

I21=sqrt(M1*w0.*S/(3*Rrsum1));

plot(I21,w)

grid on ,

hold on

%реостатная скоростная харакетристика №2%

Rrotsum2=Rr+R2+R3+R4+R5+R6+R7;

Rrsum2=Rrotsum2*Ke^2;

a=Ro/Rrsum2

Skl2=Ske*Rrsum2/Rr

M2=2*Mk*(1+a*Ske)./(S./Skl2+Skl2./S+2*a*Ske);

I22=sqrt(M2*w0.*S/(3*Rrsum2));

plot(I22,w)

grid on,

hold on

%реостатная скоростная харакетристика №3%

Rrotsum3=Rr+R3+R4+R5+R6+R7;

Rrsum3=Rrotsum3*Ke^2;

a=Ro/Rrsum3

Skl3=Ske*Rrsum3/Rr

M3=2*Mk*(1+a*Ske)./(S./Skl3+Skl3./S+2*a*Ske);

I23=sqrt(M3*w0.*S/(3*Rrsum3));

plot(I23,w)

grid on,

hold on

%реостатная скоростная харакетристика №4%

Rrotsum4=Rr+R4+R5+R6+R7;

Rrsum4=Rrotsum4*Ke^2;

a=Ro/Rrsum4

Skl4=Ske*Rrsum4/Rr

M4=2*Mk*(1+a*Ske)./(S./Skl4+Skl4./S+2*a*Ske);

I24=sqrt(M4*w0.*S/(3*Rrsum4));

plot(I24,w)

grid on ,

hold on

%реостатная скоростная харакетристика №5%

Rrotsum5=Rr+R5+R6+R7;

Rrsum5=Rrotsum5*Ke^2;

a=Ro/Rrsum5

Skl5=Ske*Rrsum5/Rr

M5=2*Mk*(1+a*Ske)./(S./Skl5+Skl5./S+2*a*Ske);

I25=sqrt(M5*w0.*S/(3*Rrsum5));

plot(I25,w)

grid on,

hold on

%реостатная скоростная харакетристика №6%

Rrotsum6=Rr+R6+R7;

Rrsum6=Rrotsum6*Ke^2;

a=Ro/Rrsum6

Skl6=Ske*Rrsum6/Rr

M6=2*Mk*(1+a*Ske)./(S./Skl6+Skl6./S+2*a*Ske);

I26=sqrt(M6*w0.*S/(3*Rrsum6));

plot(I26,w)

grid on ,

hold on

%реостатная скоростная харакетристика №7%

Rrotsum7=Rr+R7;

Rrsum7=Rrotsum7*Ke^2;

a=Ro/Rrsum7

Skl7=Ske*Rrsum7/Rr

M7=2*Mk*(1+a*Ske)./(S./Skl7+Skl7./S+2*a*Ske);

I27=sqrt(M7*w0.*S/(3*Rrsum7));

plot(I27,w)

grid on ,

hold on

Xrot=0.078; lb=120; le=161.75; Rrdin=0.041; Ern=320; Ke=1.12; w0=62.8

Xr=Xrot*Ke^2;

%построение динамической характеристики при посадке груза%

l_0t=[0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2]

E_ft=[0.13 0.26 0.39 0.52 0.627 0.736 0.825 0.895 0.944 1.0 1.046 1.07 1.1 1.122 1.143 1.163 1.179 1.196 1.2 1.223 1.236 1.25]

l_0=linspace(0.1,2.2)

E_f=interp1(l_0t,E_ft,l_0,'spline')

l0=l_0.*lb;

Ef=E_f.*Em*Ke/sqrt(3);

x0=Ef./l0;

S_=Rrdin*sqrt((le^2-l0.^2)./(l0.^2.*(x0+Xr).^2-le^2*Xr^2))

w=S_.*w0;

l_2=sqrt((le^2-l0.^2)./(1+2*Xr./x0));

l2d=l_2.*Ke;

M=3.*l_2.^2*Rrdin./(w0.*S_);

plot(l2d,-w),

grid on