Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Виноград(Вектор_управ_АД)321стр

.pdf
Скачиваний:
109
Добавлен:
02.04.2015
Размер:
8.03 Mб
Скачать

Для частного случая симметричных процессов удобнее пользоваться следующей формой записи уравнений:

Urs = Rs Irs + ddtΨs + jωk Ψrs ; r

0 = Rr Irr + ddtΨr + j(ωk ω)Ψrr ;

Ψs = Ls Is + Lm Ir

;

 

(3.9)

r

 

r

 

r

;

 

Ψr = Lr Ir + Lm Is

 

 

 

3

 

 

L

r

r

 

M =

 

Z

 

m

Ψ ×I

 

;

2

 

L

 

 

 

p

r

 

s

 

 

 

 

 

r

 

 

 

 

J dω = M Mc . Z p dt

Дифференциальные уравнения электромагнитных процессов обычно упрощают, записывая их относительно каких-либо двух векторных переменных состояния, исключая остальные с помощью уравнений связи. Рассмотрим несколько вариантов более удобной записи уравнений.

3.3. Запись уравнений относительно потокосцеплений статора и ротора

Из уравнений связи (третье и четвертое уравнения системы (3.9)) выразим векторы токов статора и ротора, подставим их в первое, второе и пятое уравнения системы (3.9).

r

 

 

Ls

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψ

=

 

Lm I

 

;

 

 

 

 

 

 

 

 

 

 

 

 

rs

 

 

 

 

rs

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψr

Lm

 

Lr Ir

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

Lm

 

 

 

 

 

 

 

 

 

 

 

L

r

 

 

 

 

 

 

 

 

Ψs

 

 

 

r

 

 

r

 

 

 

 

Ψ

 

 

r

r

 

r

 

 

 

r

L

 

 

(L

 

 

r

 

 

s

rs

 

 

 

 

 

Ψ

 

 

Ψ − L

Ψ )

 

 

L

Ψ

 

 

(L Ψ − L Ψ )

 

Is =

 

 

 

r

r

 

=

r

s

m

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

; I

 

=

 

m

r

 

=

s r

m s

;

 

 

 

Ls

Lm

 

 

L L

L2

r

 

 

 

 

 

 

 

 

L L

 

L L

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lm

Lr

 

 

 

 

s

r

 

m

 

 

 

 

s

m

 

 

s r

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lm

Lr

 

 

 

 

 

31

 

r

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

(Lr Ψs Lm Ψr )

 

 

 

 

dΨs

 

 

 

U s

= Rs

 

+

 

+ jωk Ψs ;

 

 

 

 

Ls Lr

L2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

r

 

 

 

r

 

0 = R

 

(Ls Ψr Lm Ψs )

 

 

 

 

dΨr

+ j(ω

 

 

r

 

+

k

ω)Ψ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ls Lr L2m

 

 

 

 

 

 

 

 

dt

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или в форме Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

R

 

L

 

 

 

 

 

 

R

 

L

 

 

 

 

 

 

 

 

 

 

 

dΨ

s

 

 

 

 

 

 

 

s

r

r

 

 

 

 

s

m

 

r

 

r

 

r

 

 

 

 

 

= −

 

 

 

Ψs +

 

 

 

 

 

Ψr

+U s jωk Ψs ;

 

dt

 

 

 

 

Km

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Km

 

 

 

 

 

 

r

 

 

 

 

 

 

R

L

 

 

 

 

 

R

L

 

 

 

 

 

 

 

 

 

 

 

 

 

dΨ

r

 

 

 

 

m

r

 

 

s

 

r

 

 

 

r

 

 

 

 

=

 

 

 

r

 

 

Ψs

 

 

 

r

 

 

 

Ψr j(ωk

− ω)Ψr

,

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

Km

 

 

 

 

 

 

 

 

 

 

 

Km

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где K

m

= L L L2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

r

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение электромагнитного момента запишется как

 

3

 

Lm

r

(Lr Ψs LmΨr )

 

M =

Z p

 

 

 

 

Ψr ×

 

 

;

2

L

L L

L2

 

 

 

r

 

s r

m

 

M =

3

Z

 

L

m

r

r

 

 

 

 

Ψ

×Ψ

.

2

p Km

 

 

r

s

 

Уравнение движения остается без изменений.

Переходя к проекциям в системе координат (x,y), получим

dΨsx = −

 

Rs Lr

 

Ψ +

 

Rs Lm

 

Ψ +U

x

+ω

Ψ ;

 

 

 

 

 

 

 

 

dt

 

 

 

Km

 

 

sx

 

 

 

 

 

 

 

rx

 

k

sy

 

 

 

 

 

 

 

 

 

Km

 

 

 

 

dΨsy

 

 

 

R L

 

 

 

 

 

 

R L

 

 

 

 

 

= −

s r

Ψ +

 

 

s m

Ψ +U

y

ω

Ψ ;

 

 

 

 

 

 

 

 

dt

 

 

 

Km

 

 

sy

 

 

 

 

 

 

 

ry

 

k

sx

 

 

 

 

 

 

 

 

 

Km

 

 

 

 

dΨrx =

Rr Lm

 

Ψ −

Rr Ls

 

Ψ +(ω

k

ω)Ψ ;

 

 

dt

 

Km

 

 

sx

 

 

 

 

 

 

 

rx

 

 

ry

 

 

 

 

 

Km

 

 

 

 

dΨry

 

R L

 

 

 

R L

 

 

 

 

 

=

 

 

r m

Ψ −

 

 

 

r s

Ψ −(ω

ω)Ψ ;

 

 

 

 

 

 

 

 

dt

 

Km

 

 

sy

 

 

 

 

 

 

 

ry

k

 

 

rx

 

 

 

 

 

Km

 

 

 

 

M =

3

Z

 

Lm

(Ψ Ψ −Ψ Ψ

);

 

p Km

2

 

rx sy ry sx

 

J dω

= M Mc .

 

Z p

dt

 

32

(3.10)

(3.11)

(3.12)

Система уравнений (3.12) представляет собой математическую модель АД, записанную в ортогональной системе координат (x,y), относительно потокосцеплений статора и ротора. Она довольно удобна для вычислений на ЭВМ. Все остальные переменные машины: токи статора, ротора, намагничивания, потокосцепление намагничивания рассчитываются на основе алгебраических уравнений связи этих переменных с потокосцеплениями статора и ротора.

I

sx

=

 

1

 

(L Ψ − L Ψ ); I

sy

=

 

1

 

(L Ψ − L Ψ );

 

 

 

 

 

 

 

 

 

r

sx

m rx

 

 

 

 

 

r

sy

m ry

 

 

 

 

Km

 

 

 

 

 

 

Km

 

 

I

rx

=

 

1

(L Ψ − L Ψ ); I

ry

=

 

1

(L Ψ − L Ψ );

 

 

 

 

 

 

 

 

 

s

rx

m sx

 

 

 

 

s

ry

m sy

 

 

 

 

Km

 

 

 

 

 

Km

 

 

Imx = (Isx + Irx ); Imy = (Isy + Iry ); Ψmx = Lm Imx ; Ψmy = Lm Imy .

Структурная схема асинхронного двигателя, построенная по уравнениям (3.12), изображена на рис.3.5.

3.4. Запись уравнений относительно тока статора и потокосцепления ротора

Из уравнений связи исходной системы уравнений (3.9) выразим ток ротора и потокосцепление статора и подставим их в остальные уравнения.

r

 

 

 

1

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

 

 

r

 

 

 

 

 

r

 

 

 

 

L

m

 

r

 

 

 

r

 

 

 

 

I

r

 

=

 

 

 

 

 

(Ψ − L I

s

);

Ψ = L

s

I

s

+

 

 

 

(Ψ − L I

s

)

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lr

 

 

 

 

r

 

 

 

 

m

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

Lr

 

 

 

r

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2

 

 

r

 

 

 

 

 

L

m

 

r

 

 

 

 

 

 

 

r

 

 

 

L

m

r

 

 

 

 

 

 

 

 

 

 

= L

(1

 

 

 

m

 

)I

s

+

 

 

 

Ψ =σL

I

s

+

 

 

 

Ψ

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

Ls Lr

 

 

 

 

 

 

Lr

 

 

r

 

 

 

 

s

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lr

 

 

 

 

 

 

 

 

 

 

где σ =1

 

 

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

коэффициент рассеяния.

 

 

 

 

 

 

 

 

 

 

L L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

dI

s

 

 

 

L

m

 

dΨ

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

L

m

r

U

 

= R

 

I

 

+σL

 

 

 

 

 

 

+

 

 

 

 

r

 

+ jω

 

 

(σL

I

 

+

 

 

 

 

Ψ );

 

 

 

s dt

Lr

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

s

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

s

 

s

 

 

 

Lr

 

r

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

r

 

 

 

 

 

 

 

r

 

 

 

 

dΨ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

0 =

 

 

 

(Ψ − L

m

I

s

) +

 

 

 

r

 

+ j(ω

k

 

ω)Ψ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lr

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

33

34

Ux

ωk

U y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rs Lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Km

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rs Lr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rr Ls

 

 

 

 

 

 

 

 

Km

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψsx

 

 

 

 

 

 

 

 

 

 

Km

 

Ψ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

R L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r m

 

 

 

 

 

 

 

 

rx

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Km

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Ψsy

 

 

Rr Lm

 

 

 

Ψry

 

 

 

 

 

 

 

 

 

p

 

 

Km

 

 

 

 

 

 

 

Rs Lr

 

 

 

 

 

 

Rr Ls

 

 

 

 

 

 

Km

Km

 

Rs Lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Km

 

 

 

 

 

Lm M

Mc

ωr

3

1

2 Z p Km

Jp

 

ω Z p

Рис.3.5. Структурная схема АД в ортогональных осях (x,y)

 

 

 

dIr

 

 

 

 

 

 

 

r

 

 

 

 

L

m

 

dΨ

r

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

L

m

 

 

 

r

 

 

 

 

σL

 

 

s

 

= −R

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+U

 

jσL

 

ω

 

 

I

 

 

j

 

 

 

 

ω

Ψ ;

 

 

s dt

 

 

 

 

Lr

 

 

dt

 

 

 

 

 

 

 

 

 

Lr

 

 

 

 

 

 

 

 

 

 

s

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

s

 

 

k

 

s

 

 

 

 

 

k

 

 

r

 

 

 

 

 

r

 

 

 

 

R

 

 

 

 

 

 

 

 

 

L

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dΨ

 

 

 

 

r

r

 

 

 

 

 

 

m

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

= −

 

 

Ψ +

 

 

 

 

 

 

 

 

 

I

s

j(ω

k

ω)Ψ

r

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

Lr

 

 

r

 

 

 

 

 

 

 

Lr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из первого уравнения можно исключить производную потокосце-

пления ротора, подставив в него второе уравнение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dIr

 

 

 

 

 

 

 

 

r

 

 

 

 

 

L

m

 

 

 

 

 

 

 

R

r

 

r

 

 

 

 

L

m

R

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

σL

 

 

s

 

= −R

 

 

I

 

 

 

 

 

 

(

 

 

 

 

 

Ψ

 

 

+

 

 

 

 

I

 

 

j(ω

 

 

 

ω)Ψ

 

) +

 

s dt

 

s

s

 

 

 

Lr

 

 

Lr

r

 

 

Lr

 

 

 

s

 

k

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

Lm

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+U s jσLsωk I s

j

 

ωk Ψr

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

m

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

L

m

R

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

L

m

 

 

 

r

 

 

 

r

= −(Rs

+ (

 

 

)

2 Rr )I s

+

 

 

 

 

 

 

 

 

Ψr

jσLsωk I s j

 

 

 

ωΨr

+U s ;

 

 

 

 

 

 

 

 

L2r

 

 

 

 

 

 

 

 

dIr

 

 

 

 

Lr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

L

m

R

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

L

m

 

 

r

 

 

 

 

r

 

 

σL

 

 

s

 

= −K

 

 

R

I

 

 

+

 

 

 

 

 

 

Ψ

 

 

jσL

ω

 

 

I

 

 

 

j

 

 

 

 

ωΨ

 

+U

 

,

s dt

 

 

 

 

 

 

 

 

L2r

 

 

 

 

 

 

 

 

 

Lr

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

s

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

s

 

 

 

k

 

 

s

 

 

 

 

 

 

 

 

r

 

 

 

s

 

где Kr

=1

+ (

 

Lm

)2

Rr

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lr

 

 

Rs

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введя обозначения для постоянных времени статорной и роторной

цепей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T = σLs ;

T =

Lr

, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

Rs

 

r

 

 

 

 

Rr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dI

s

 

 

 

 

 

r

 

 

 

L R

 

 

r

 

 

1

r

 

 

 

r

 

L

 

r

T

 

 

 

= −K

I

 

 

+

 

m r

 

Ψ +

 

 

 

 

U

 

jT

ω

I

 

j

m

 

ωΨ ;

 

 

 

 

 

 

 

 

R

 

 

L R

 

s dt

 

 

 

 

 

r

 

s

 

 

 

L2 R

 

 

 

r

 

 

s

s

k

 

s

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

r s

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

r

s

 

 

 

 

 

r

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

dΨr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −Ψr + Lm Is jTr (ωk ω)Ψr ;

 

 

 

 

 

(3.13)

Tr dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

L

 

 

 

r

 

 

r

 

 

 

J

 

 

dω

 

 

 

 

 

 

 

 

 

M =

 

 

 

 

Z

p

m

 

 

Ψ × I

s

;

 

 

 

 

 

 

 

 

= M M

c

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

Lr

 

 

 

r

 

 

 

 

 

 

Z p

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогичным образом могут быть получены уравнения АД, запи-

санные относительно других пар векторных переменных, например относительно токов статора и ротора, тока статора и потокосцепления намагничивания и т.д.

35

3.5. Уравнения в преобразованных координатах для частных случаев

Переход от уравнений в ортогональной системе координат (x,y), вращающейся с произвольной скоростью ωk , к уравнениям в ортогональной системе координат (α, β) , неподвижной относительно статора, осуществляется с помощью подстановки ωk = 0 и замены ин-

декса «x» на « α», а индекса «y» на « β ».

Для системы уравнений, записанной относительно тока статора и потокосцепления ротора в координатах (x,y) выполним переход в ортогональную систему координат (d,q), ориентированную по вектору

потокосцепления

ротора.

 

В

этом

случае

 

 

 

ωk =ωψ ,

Ψrd = Ψr .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

dId

 

 

= −K

 

I

 

+

Lm Rr

Ψ +T ω I

 

+

 

1

U

 

;

 

 

 

 

 

 

 

 

 

L2r Rs

 

 

 

 

 

 

s dt

 

 

 

 

r

 

d

 

 

 

r

s

 

ψ

q

 

 

Rs

 

 

 

d

 

 

T

 

dIq

= −K

 

I

 

T ω I

 

L

 

ωΨ +

1

 

U

 

;

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L R

 

 

R

 

 

 

 

s dt

 

 

 

r

 

q

 

s

ψ

 

d

 

 

s

 

r

 

s

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

T

dΨr

= −Ψ + L

 

I

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

m

 

d

 

 

Iq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

=ω +ω

 

 

=ω +

L

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

 

 

 

 

 

 

 

 

 

 

 

 

 

Tr

 

 

Ψr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

Lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M =

 

Z

p

Ψ I

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

L

 

 

 

r q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

dω

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= M Mc .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z p

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψrq = 0 ,

(3.14)

Заметим, что в данной системе уравнений существенно упростилась модель роторной цепи, выражение для электромагнитного момента. Количество дифференциальных уравнений уменьшилось на единицу. В установившихся режимах работы двигателя все преобразованные переменные оказываются постоянными величинами. В связи с этим данная система уравнений очень удобна для расчетов процессов в машине и для синтеза векторной системы управления в координатах (d,q). Структурная схема АД в координатах (d,q) изображена на рис.3.6.

36

37

 

 

Lm Rs

 

 

 

 

 

 

 

1 ωr

 

 

 

 

 

 

Mc

 

Z p

 

 

L2 R

 

 

 

 

 

Ud

 

r

 

s

Lm ψr

 

 

 

 

 

 

ω

1

1

I

d

3

L

M

 

Z

p

 

 

 

 

 

 

m

 

 

 

 

RS

Ts P + Kr

 

 

Tr P +1

2 Z p Lr

 

 

J

 

 

 

Ts ω

 

 

 

 

 

 

 

 

 

 

 

 

ψ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

ω

s

 

 

 

 

Ts

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

Tr

 

 

 

 

Uq

 

 

 

 

 

 

 

 

 

 

 

1

1

Iq

 

 

 

 

ω

 

 

 

RS

Ts P + Kr

 

 

 

 

 

 

 

ψ

 

 

 

 

 

 

 

 

 

 

 

 

 

Lm

Lr Rs

Рис.3.6. Структурная схема АД в координатах (d,q)

3.6. Математическое описание АД с учетом насыщения цепи намагничивания

Рассмотренные выше модели АД наиболее широко применяются для анализа и синтеза процессов в электроприводах с векторным управлением. Однако если по условиям работы АД потокосцепления претерпевают существенные изменения, то желателен учет изменения параметров математической модели, связанный с эффектом насыщения магнитной системы. Для АД, питаемых от преобразователей частоты (ПЧ), обычно предполагают, что кратность выходных токов в рабочих режимах привода ограничена (составляет не более 23 значений номинального тока двигателя) и эффектом насыщения цепей рассеяния можно пренебречь. Тем более что, как показывают исследования, эти изменения достаточно слабо влияют на динамические процессы.

Гораздо более существенно сказывается эффект насыщения цепи намагничивания. В частности, в рабочих режимах привода с двигате-

лями серии 4А индуктивность намагничивания Lm может изменяться

до 30 %.

В асинхронном электроприводе регулирование потокосцепления требуется в следующих режимах:

1)при работе привода на скоростях, превышающих номинальную (во 2-й зоне регулирования скорости в режиме постоянства мощности), происходит ослабление поля;

2)при оптимизации энергетических характеристик привода требуется регулирование потока намагничивания в зависимости от нагрузки;

3)при оптимизации динамических характеристик привода, в частности быстродействия процессов "в большом", также требуется регулирование потокосцепления.

Для учета эффекта насыщения используется один из следующих методов: метод статических индуктивностей либо метод динамических индуктивностей. Последний существенно более сложный. Для синтеза систем управления приводом обычно используется более простой метод статических индуктивностей, дающий, тем не менее, достаточно высокую точность в описании динамических процессов. В этом методе нелинейность цепи намагничивания учитывается статической зависимостью между потокосцеплением и током намагничивания, задаваемой таблично или с помощью аналитической аппроксимации:

r

r

 

r

Ψm

 

Ψm = Lm (Im )Im

или

Im =

 

.

Lm (Ψm )

 

 

 

 

 

38

В частности, если характеристика намагничивания задана в табличной форме совокупностью точек (Ψm.n , Im.n ), n =1,..., N , то для

ее аппроксимации удобно использовать следующую гладкую аналитическую функцию:

k

 

Im (Ψm ) = giΨm2i1,i =1,2,...,k .

(3.15)

i=1

Выбор нечетных степеней полинома (3.15) обеспечивает симметрию характеристики относительно начала координат. В связи с возможным разбросом точек исходного массива данных коэффициенты

полиномиальной аппроксимации gi находятся методом наименьших

квадратов. Для удовлетворительной аппроксимации характеристик намагничивания машин нормального исполнения (в диапазоне изме-

нения Ψm = (0...1,3)Ψm nom , включающем зону существенного на-

сыщения) практически достаточно в формуле (3.15) ограничиться k = 4 , т.е.

Im (Ψm ) = g1Ψm + g2Ψm3 + g3Ψm5 + g4Ψm7 .

Применение к этой функции процедуры полиномиальной регрессии приводит к следующему матричному уравнению:

c1

c2

c3

c4

g1

 

d1

 

 

c c c c

g

 

 

d

 

 

,

2

3

4

5

 

2

 

=

2

 

 

c4

c5

c6

 

 

 

 

 

 

 

c3

g3

 

d3

 

 

c4

c5

c6

c7 g4

d4

 

 

 

N

 

 

 

 

 

 

 

 

где c j

= (Ψm2 j )n ; n =1,2,..., N ; j =1,2,....,2k 1;

n=1

N

d j = (Ψm2 j1 )n (Im )n ; k = 4;

n=1

n, N соответственно порядковый номер и общее количество точек

из массива исходных данных.

В результате решения этого уравнения находятся коэффициенты полинома gi и определяется аналитическая зависимость кривой на-

магничивания.

Отметим, что для неявнополюсных машин, к которым относится АД, нелинейный оператор зависит только от амплитуды векторных переменных и не зависит от их направления. Векторы тока и потокосцепления намагничивания остаются при этом сонаправленными.

39

С учетом этого математическая модель АД в преобразованной ортогональной системе координат, вращающейся с произвольной скоро-

стью ωk , примет вид

r

 

 

 

 

r

 

 

 

 

dΨ

s

 

 

 

r

U s

= Rs Is

+

 

 

 

 

 

 

+ jωk Ψs ;

 

 

rdt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

dΨ

 

 

 

 

 

r

0 = R

I

r

+

 

 

 

 

 

 

r

 

+ j(ω

k

ω)Ψ ;

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

dt

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψs = Lσs Is

 

+ Ψm ;

 

 

 

r

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

 

 

 

(3.16)

Ψ = L

 

I

 

 

+ Ψ ;

 

 

rr

 

 

 

σr

 

r

 

 

 

 

 

 

m

 

 

 

 

Ψm = Lm (Im )Im ;

 

 

 

 

 

 

3

 

 

 

L

m

 

r

 

 

r

 

 

M =

 

 

 

Z

p

 

 

 

 

Ψ

r

× I

s

;

 

2

 

 

Lr

 

 

 

 

 

 

 

 

 

 

 

 

J dω = M M c . Z p dt

По аналогии с рассмотренными выше преобразованиями модели (3.9) данную модель также можно преобразовать, записав ее относи-

тельно выбранной пары векторных переменных, например Is и Ψr ,

r

Ψr и т.д.

Ψs и

В общем случае система нелинейных дифференциальных и алгебраических уравнений (3.16) решается итерационно. Однако если предположить, что между шагами расчета ток намагничивания изменяется незначительно, то итерации можно исключить.

Иногда для описания процессов в машине удобно пользоваться не абсолютными, а относительными значениями переменных. Переход к системам относительных (безразмерных) переменных и параметров осуществляется на основе их базисных значений [38].

На рис.3.7 приведена в относительных единицах типовая зависимость Lm (Ψm ) для двигателей серии 4А.

40