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

Дискретно-полевые модели электрических машин. Часть I II

.pdf
Скачиваний:
21
Добавлен:
15.11.2022
Размер:
16.43 Mб
Скачать

Рис. 7.10. Распределение тангенциальной составляющей магнитной индукции (огибающая) в исследуемой области

Рис. 7.11. Распределение тангенциальной составляющей магнитной индукции (огибающая) вдоль радиальной координаты

211

свидетельствует об их идентичности (см. рис. 7.3, 7.8–7.11). Разность значений радиальной индукции в первой точке на рис. 7.3 и 7.9 объясняется дискретностью задания токовой нагрузки при решении двумерной задачи и в виде бегущей волны – в одномерном варианте.

8. МОДЕЛИРОВАНИЕ РАБОЧЕГО РЕЖИМА КОРОТКОЗАМКНУТОГО АСИНХРОННОГО ДВИГАТЕЛЯ

8.1. СТАЦИОНАРНЫЙ РЕЖИМ КОРОТКОЗАМКНУТОГО АСИНХРОННОГО ДВИГАТЕЛЯ НА ОСНОВЕ ОДНОМЕРНОЙ МОДЕЛИ

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

Наряду с ранее принятыми допущениями, используемыми в режиме идеального холостого хода, введём дополнительные.

1.Короткозамкнутая обмотка ротора равномерно распределена

ввоздушном зазоре. Влияние короткозамыкающих колец учитывается соответствующим увеличением сопротивлений стержней:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Rр = Rст +

RK

 

= RCT

1

+

RK

 

 

 

 

 

= RстK р , (8.1)

 

 

 

 

 

 

 

 

 

πp

 

Rст

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

πp

 

 

2sin

 

 

 

 

 

 

 

 

 

2sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z 2

 

 

 

 

 

 

 

 

 

Z 2

 

212

где Kр – коэффициент, учитывающий влияние короткозамыкающих колец.

Удельная электропроводность материала стержня

 

γр

=

l ст

.

(8.2)

 

 

 

RрS ст

 

Если обмотка ротора равномерно распределена в зазоре машины, то расчётная электропроводность материала ротора определяется как

γ

 

=

γрS стZ 2

.

(8.3)

э

 

 

 

πδDр

 

2.Электропроводность и магнитная проницаемость материала ротора не зависят от пространственных координат ϕ и Z.

3.Обмотка ротора симметрична.

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

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot

 

rot A

= J ст + J р ,

(8.4)

µ

 

 

 

 

 

 

 

 

 

 

где J р – плотность тока ротора.

Выполняя операции rot, преобразуя полученное выражение подобно режиму холостого хода и проектируя уравнения на координатную ось Z, получим

1

 

1AZ

 

+

1

 

 

1

AZ

 

= −J сZ J рZ.

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.5)

 

 

 

µ

 

2

 

 

 

µ

φ

R R

R

 

R

φ

 

 

 

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

213

1

 

1

A

=

1

 

 

1

 

A

+

1

 

 

 

1

 

A

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

R

 

 

 

 

, (8.6)

 

 

 

µ

 

R

 

µ

 

 

 

 

 

 

R R

R

 

R

 

R

я

R R

µ

 

R

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

первый член которой представляет распределение магнитной индукции в ярме машины, а второй– потоки рассеяния ротора. Здесь и в дальнейшем индексы единственных составляющих плотности тока и векторногопотенциалабудутопущены.

Первую составляющую (8.6) преобразуем подобно тому, как это делалось для режима холостого хода. В результате преобразований будем иметь

1

 

1

A

= −qA .

 

 

 

 

 

R

 

 

 

 

 

(8.7)

 

 

 

 

µ

 

R R

R

я

 

 

 

 

 

 

 

 

 

 

 

 

Потоки рассеяния в пазах ротора пропорциональны току ротора. Поэтому вторую составляющую (8.6) можно записать в виде

1

 

1

A

 

= Kσ 2

 

 

 

 

 

 

R

 

 

 

 

 

 

J р .

(8.8)

 

 

 

 

µ

 

R R

R

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент Kσ 2 в этом выражении можно рассматривать как отношение индуктивного сопротивления рассеяния ротора к сопротивлению намагничивающего контура X σ 2 Xµ . Подставляя выраже-

ния (8.7) и (8.8) в (8.5), получим

1

 

2A

qA = −µ

0 J ст µ

0 (1+ Kσ 2 ) J р .

(8.9)

 

 

R02 φ2

 

 

 

Плотность равномерно распределённого в зазоре тока ротора выражается через векторный потенциал в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

J = γ

э

+V × B+

grad φ .

(8.10)

 

 

 

 

 

t

 

 

214

Для задачи в плоскопараллельном приближении при отсутствии в проводящей среде свободных электрических зарядов grad φ = 0 .

Учитывая, что плотность тока и векторный потенциал имеют по од-

ной аксиальной составляющей JZ и АZ,

а магнитная индукция – две

составляющие B R и Bφ , плотность тока ротора с учётом принятых

допущений может быть записана в виде

 

 

 

 

A

ω

 

A

 

 

J р = γэ

 

,

(8.11)

t

 

φ

 

 

 

где ω – частота вращения ротора.

Подставляя это выражение в (8.9) и выполняя промежуточные

преобразования, в окончательном виде будем иметь

 

 

 

1

 

2A

− µ

γ

(1+ K

 

)ω

A

µ

γ

(1+

K

)

AqA= −

µ

. (8.12)

 

 

 

 

 

R02 φ2 0

э

 

σ 2

 

φ

0

э

 

σ 2

t

0 J с

 

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

Покажем, что представленная одномерная модель с учётом вышеуказанных допущений полностью соответствует процессам асинхронной машины.

Для стационарного режима будем считать, что при гармоническом изменении величин во времени ток статорной обмотки создаёт бегущую волну:

J

= J

exp j (ω

0

t

pφ) .

(8.13)

c

с.м

 

 

 

 

Тогда и векторный потенциал также представляет бегущее магнитное поле

 

exp j (ω

 

t

pφ) .

(8.14)

A = A

0

м

 

 

 

 

215

Подставляя эти выражения в уравнение (8.12) и выполняя промежуточные преобразования, получим

p2

+ jµ

γ

 

(1

+ K

 

 

)(

ω

 

 

pω)+

q

 

 

=

µ

 

 

.

(8.15)

 

 

 

σ 2

0

 

с.м

 

0

 

э

 

 

 

 

 

 

 

 

 

 

Aм

 

0 J

 

 

R02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Частота вращения ротора в функции скольжения s

 

 

 

 

 

 

 

 

 

 

 

 

 

ω =

ω 0

(1 s ) ,

 

 

 

 

 

 

 

 

 

(8.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω 0

 

pω= ω

0s .

 

 

 

 

 

 

 

 

 

(8.17)

В этом случае выражение векторного потенциала из (8.15)

 

 

 

 

=

 

 

 

 

 

 

 

 

 

µ0 J с.м

 

 

 

 

 

 

 

.

 

 

 

(8.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aм

 

 

 

 

 

µ

 

γ

 

(1 +

K

 

)ω

 

s

 

q

 

 

 

 

 

 

 

 

 

 

 

 

0

э

σ 2

0

 

 

 

 

 

 

 

 

 

 

α

2 1 + j

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как было показано выше (7.40), выражение 1 +

 

q

представляет

α

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

коэффициент насыщения магнитной цепи машины K . Следовательно, векторный потенциал может быть представлен в виде

 

=

 

 

µ0 J с.м

 

 

 

 

 

 

 

.

(8.19)

 

 

µ

γ

(1 + K

 

)ω

 

s

Aм

 

+ j

σ 2

0

 

 

 

 

α 2 K µ 1

0

э

 

 

 

 

 

 

 

 

 

 

 

 

α

2

K µ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Размерность выражения

 

 

 

 

 

µ0ω 0

 

Ом с 1

2

 

 

м

с

м= Ом м

 

α 2 K µ

совпадает с размерностью удельного сопротивления. Обозначим это отношение как

X 20 =

µ0ω 0

.

(8.20)

 

 

α 2 K µ

 

216

Рассмотрим физический смысл этой величины. Индуктивное сопротивление намагничивающего контура, отнесённое к обмотке ротора, определяется как отношение ЭДС обмотки ротора е2 к намагничивающему току I0:

X µ 2 =

e2

.

(8.21)

 

 

I 0

 

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

X µ 2 =

E

 

2lδ W 2 K об2

=

E

 

lδ

W 2 K об2 m2

.

(8.22)

 

 

 

 

 

 

 

 

 

 

J 0 δ2 pτ/ m2

 

J 0

 

pτ δ

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

= X µ 2

 

 

pτ δ

 

 

.

 

(8.23)

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

J 0

 

 

 

 

 

m2

 

 

 

 

 

δ W 2 K об2

 

 

С другой стороны, индуктивное сопротивление намагничивающего контура, приведённое к обмотке ротора, выраженное через параметры машины, записывается в виде [4]

 

 

 

 

 

 

4m2 f

 

 

2

 

2

 

 

 

X µ 2 =

 

µ0τ lδ

 

W 2 K

об2

.

(8.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

δK µ

p

 

 

 

Подставляя это выражение в (8.23) и выполняя преобразования,

в окончательном виде получим

 

 

 

 

 

 

 

 

 

 

 

 

E

=

4 f µ

0τ2W 2 K об2

.

 

(8.25)

 

 

 

J 0

π

 

 

 

K µ

 

 

 

 

 

 

 

 

 

 

 

 

Для короткозамкнутой обмотки ротора W2

= 1/2, Kоб2 = 1 и от-

ношение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

=

f µ0 τ2

=

ω0µ0

 

= X 20 .

(8.26)

 

J 0

 

π2 K µ

α 2 K µ

 

 

 

 

 

 

 

 

 

217

Таким образом, X 20 можно рассматривать как удельное индуктивное сопротивление намагничивающего контура, приведённое к обмотке ротора. В этом случае приведённое ниже произведение представляет собой сумму индуктивныхсопротивлений

 

 

 

 

 

X 20 (1 + Kσ

2 ) = X 20 + Xσ

2 ,

 

 

 

 

 

 

(8.27)

и векторный потенциал записывается как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

µ0 J с.м

 

 

 

 

 

 

 

.

 

 

 

(8.28)

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aм

 

α

 

 

 

1 + jγ s (

X

20 + X

σ 2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

K µ

 

 

 

 

э

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом допущений напряжённость электрического поля опре-

деляется выражением

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aм

 

 

Aм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eм = −

 

 

t

ω

∂ϕ

 

 

 

 

= − jω

0 sAм .

 

 

 

 

 

(8.29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Умножая левую и правую части уравнения (8.28) на – jω

0·γ э·s, полу-

чим плотность тока ротора

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J р.м = γ E sм=

 

 

 

jγ эs X 20 J с.м

 

 

 

.

 

 

 

(8.30)

 

 

 

 

 

1 + jγ эs ( X 20+ X σ

2 )

 

 

 

Разделив правую часть этого выражения на γ эs

и приведя его к обще-

му знаменателю, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

+ j

X

 

 

 

+ jX

 

 

 

 

 

= − j

X

 

 

 

(8.31)

 

 

 

 

 

 

 

 

 

 

 

 

J р.м γ эs

 

 

20 J р.м

 

 

 

 

σ

2 J

р.м

 

 

 

 

20 J с.м

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

+ jX

 

 

 

 

= − j

X

20

 

 

 

j

X

20

 

 

=

 

.

(8.32)

 

 

 

 

 

 

 

 

 

 

J р.м γ эs

 

 

σ 2 J

р.м

 

 

 

 

 

 

J с.м

 

 

 

J

р.м

 

E м

 

 

Отсюда плотность тока ротора

218

Рис. 8.1. Схема замещения роторной цепи асинхронной машины

J р.м =

 

E м

.

(8.33)

 

 

1

+ jX σ 2

 

 

 

γ эs

 

 

 

 

 

 

Полученное выражение аналогично выражению тока ротора асинхронной машины

I 2s

=

 

E

2s

.

(8.34)

R2

+ j X 2s

 

 

 

 

 

 

 

S

 

 

 

Разница заключается в том, что в выражении (8.33) использованы распределённые (дифференциальные) параметры вместо интегральных величин выражения (8.34).

Полученному уравнению магнитного поля (8.12) соответствует роторная цепь Т-образной схемы замещения асинхронной машины

(рис. 8.1).

Таким образом, решая краевую задачу с принятыми допущениями, получают решение, соответствующее реальным условиям асинхронной короткозамкнутой машины с насыщением магнитопровода и рассеянием ротора.

Электромагнитные процессы асинхронной машины могут быть описаны с использованием схемы замещения машины с полым немагнитным ротором [37].

Для осуществления перехода к схеме замещения машины с полым немагнитным роторам запишем выражение напряжённости электрического поля, используя (8.29)

 

= − jω

s

=

j X 20 s J с.м

.

(8.35)

 

E

0

Aм

1 + j X 20sγ σ

 

219

Плотность тока ротора

 

 

 

 

 

 

= γ

 

= j X 20γ s J с.м .

 

(8.36)

J

 

 

E

1 + j X 20γ sσ

 

 

Полученное выражение может быть преобразовано к виду

 

 

1

 

 

 

 

 

 

 

J

 

 

 

+ j X 20σ = − j X 20 J с.м .

(8.37)

 

 

 

γs

 

 

 

 

Умножив и разделив левую часть уравнения на σ, получим

 

 

 

1

 

 

 

 

 

σJ

 

 

 

 

+ j X 20 = − j X 20 J

с.м .

(8.38)

γσs

 

 

 

 

 

 

 

Если обозначить σJ= J2′′м , то полученное уравнение может быть представлено в виде

J2′′м

1

= − j X 20 (Jс.м + J2′′м ) .

(8.39)

 

 

γσs

 

Этому уравнению соответствует схема замещения вторичной цепи асинхронной машины с полым немагнитным ротором, индуктивное сопротивление ротора которой равно нулю.

Преобразуем уравнение Кирхгофа статорной обмотки асинхронной машины. ЭДС намагничивающего контура Т-образной схемы замещения

 

= − jX

 

(J

+ J ′′

 

) = −

jX 20

(σJ

с.м

+ σJ ′′

) =

20

 

 

 

E1m

 

с.м

 

σ

 

 

(8.40)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jX 20

(σJс.м + J2′′м ).

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в уравнение σ =1 + Kσ2 =

X 20 + X σ2

, будем иметь

 

 

 

 

 

 

 

 

 

 

X 20

 

 

220