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

Учебное пособие 1562

.pdf
Скачиваний:
10
Добавлен:
30.04.2022
Размер:
1.39 Mб
Скачать

2.1.1 Переход от трехфазной модели к двухфазной

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

Рис. 2.1.1

Однако, такой переход возможен, если токи в обмотках изменяются со сдвигом по фазе на 2/3π, даже если и не по синусоидальному закону, и магнитная система машины линейна.

На рис.2.1.1 приведены условные изображения обмоток трехфазного двигателя и двухфазного. Полагаем, что ось совпадает с осью А. На рис. 2.1.2 векторы магнитных индукций катушек трехфазной и двухфаз-

ной систем. Найдем проекции BA, BB, Bc на оси

и :

B

BA

 

1

BB BC

,

2

 

 

 

 

 

 

 

 

 

 

 

B

3

BB

BC .

 

2

 

 

 

 

 

 

Если подставить значения BA, BB, BС в первое уравнение (синусоидальные функции), то получим B= 3/2 Bm sin t. Чтобы максимальное значение B было равно Bm, следует взять коэффициент 2/3, тогда:

B

2

BA

1

BB

BC

,

 

3

3

(2.1.5)

 

 

 

 

 

 

 

 

B

 

1

 

BB

 

BC .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично можно преобразовать потоки сцепления, токи и т.д.

В системе (2.1.1) уравнения статора преобразуем к осям , . Для этого умножим первое на 2/3 , второе и третье на –1/3 и сложим все уравнения:

d

2

 

1

 

1

 

2

 

1

 

1

 

2

 

1

 

1

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

3 1A

3

1B

3

1C

R1 3 i1A

3 i1B

3 i1C

3 u1A

3 u1B

3 u1C

 

 

учитывая (2.1.3) получим:

Рис. 2.1.2

 

 

 

 

Рис. 2.1.3

 

d 1

R

i

 

u .

(2.1.6a)

 

 

1

 

dt

1

 

 

 

 

 

 

 

 

 

Поступая аналогично, запишем вторую формулу:

d 1

R

i

 

u .

(2.1.6b)

 

1

dt

1

 

 

 

 

 

 

 

 

Система уравнений (2.1.6) описывает статор эквивалентной двухфазной машины. Для ротора запись несколько сложнее. Катушки ротора вращаются относительно статора (рис. 1.1), поэтому векторы B2 и B2 ротора

смещаются относительно неподвижных координат

и

статора:

B2

2

B2a cos

B2 b cos(

 

2

)

B2c cos(

 

2

) ,

3

3

3

 

 

 

 

 

(2.1.7)

B2

2

B2a sin

B2 b sin(

2

 

)

B2c sin(

 

2

 

) ,

3

3

 

3

 

 

 

 

 

 

 

 

 

т.е. если в неподвижных катушках 2

 

и 2 индукции будут изменяться по

закону (2.1.7) , то получится такое же поле, что и от вращающегося ротора. Что бы получить уравнения для ротора, умножим четвертое уравнение в (2.1.1) на 2/3cosγ, пятое на 2/3cos(γ+2π/3), шестое на 2/3cos(γ-2π/3) и

сложим их:

2

 

d 2a

cos

d 2 b

cos(

2

)

 

d 2c

cos(

 

 

2

)

3

 

 

dt

dt

3

 

dt

 

 

3

 

 

 

 

 

 

 

 

 

 

 

R

2

 

i2a cos

i2 b cos(

2

)

i2c cos(

2

 

) 0

 

 

 

 

 

 

 

 

2 3

 

3

3

 

Т.к. γ являтся функцией времени, то воспользуемся соотношением:

d

cos

d

cos

sin

d ,

dt

dt

dt

 

учтем, что

d

dt

-скорость вращения ротора, тогда

d 2

 

R 2 i2

0.

dt

2

 

 

 

Аналогично для второго уравнения, умножая на синусы соответствующих выражений:

d 2

 

R2 i2 0.

 

2

dt

 

 

Чтобы привести потоки сцепления к осям αβ, умножим первые три строки в системе (2.1.2) на соответствующие коэффициенты и сложим:

2

 

 

 

1

 

1

 

 

 

L1

2

i1A

1

i1B

1

i1C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 1A

3

 

1B 3

 

1C

3

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

L12 [

2

i 2a

cos

2

i 2 b

cos(

2

)

 

 

2

i 2c

cos(

2

)

1

i 2a

cos(

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

3

3

 

 

3

3

3

3

1

i2 b

cos

1

i2c

cos(

2

 

)

 

1

i 2a

cos(

2

)

 

1

i 2 b

cos(

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

3

3

 

 

3

3

 

3

3

 

 

1

i2c

cos

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сложим все косинусы при каждом из токов i2a, i2b, i2c, тогда:

 

L1i1 L12 i2a cos

i2 b cos(

2

) i

2c cos(

2

)

,

1

 

 

3

3

 

 

 

 

 

 

 

 

следовательно

 

L1

3

Lm i1

3

L12 i

.

 

 

 

1

2

2

2

 

 

 

 

 

Аналогично поступим с потоком Ψ, после вычислений имеем:

 

L1

3

Lm i1

3

L12 i

.

 

 

 

1

2

2

2

 

 

 

 

 

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

 

2

 

cos

 

cos(

2

 

)

 

 

 

 

cos(

2

) .

2

 

2a

2 b

 

 

 

 

 

 

 

 

 

2c

 

3

 

 

 

 

 

 

3

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После преобразований, учитывая, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos cos(

 

2

 

 

)

cos(

 

 

 

2

 

) 0,

 

 

 

 

 

3

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2 i2

3

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

L12 i1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично вычисляется поток Ψ:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

i

 

 

3

L

 

i

.

 

 

 

 

 

 

 

2

2

2

 

 

12

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Чтобы получить уравнение момента двигателя в координатах αβ, в формуле (2.1.4) сделаем следующие преобразования: первая квадратная скобка есть ток iумноженный на 3/2, во второй и третьей скобках синусы представим как разность или сумма некоторого угла и ±2/3π. В результате получим:

 

3

2

 

.

Mд

 

L12 i1 i

2 i1 i

2

2

 

 

 

 

2.1.2 Сводка формул двухфазной модели

Напряжения фаз А, В, С связаны с напряжениями фаз α и β следующими формулами:

u

2

uA

1

u B uC ,

3

3

 

 

 

 

 

 

u

 

1

 

u B

 

uC .

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

Токи фаз α и β преобразуются в токи фаз А, В, С:

iA

i ,

 

 

 

 

i B

 

1

i

 

3

 

i ,

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iC

 

1

i

 

 

3

i .

2

 

 

 

 

 

 

2

 

 

Уравнения для статора:

d 1

R

i

 

u ,

 

1

dt

1

 

 

 

 

 

 

ddt1 R1i1 u .

Уравнения для короткозамкнутого ротора:

d 2

 

R 2 i2

dt

2

 

 

d 2

 

R2 i2

dt

2

 

 

Связь потоков сцепления с токами:

 

L1

3

Lm

i1

1

2

 

 

 

 

 

 

 

L1

3

Lm

i1

1

2

 

 

 

 

 

0,

0.

23 L12 i2,

23 L12 i2,

 

L2 i2

3

 

 

 

,

 

 

 

 

 

 

 

2

2

L12 i1

 

 

 

 

 

 

 

 

 

L

i

 

3

L

 

i

.

2

2

 

12

2

 

2

 

 

 

1

 

 

 

 

 

 

 

 

 

Уравнения динамики двигателя:

d

J dt Mд Mн

Выражение для момента на валу:

 

3

2

 

.

Mд

 

L12 i1 i

2 i1 i

2

2

 

 

 

 

В данных выводах нет приведения к статорной обмотке. В большинстве книг уравнения записываются сразу для приведенных значений, что порождает массу недоразумений. Как известно, приведение к статорной обмотке осуществляется для того, чтобы токи, напряжения, потоки сцепления ротора были того же порядка, что и статорные. Коэффициент приведения n – это отношение числа витков статора и ротора (как и у трансформатора, но в двигателях нужно учитывать еще и распределенность обмоток; в данном случае считаем, что соответствующие коэффициенты одинаковы для статора и ротора). Таким образом:

u2 = nu2; i2 = i2/n ; z2 = n2z2; Ψ2 = nΨ2.

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

d

 

'2

 

 

 

 

 

'2

 

R'2 i'2

 

 

 

0,

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

'2

 

 

 

 

 

'2

 

R'2 i'2

 

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

 

3

Lm i1

 

3

 

L12 n i'

,

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

 

3

Lm i1

 

3

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

2

 

L12 n i'2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

3

 

 

 

 

 

 

,

 

 

'2

 

L2 n

i'2

 

 

 

 

L12 n i1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'

 

 

L

 

n2 i'

 

 

 

3

L

 

n

i

 

.

 

 

2

 

2

2

2

12

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

Mд

 

 

 

 

 

 

L12 n i1 i'2

 

 

 

i1 i'

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Считается, что приведенные значения индуктивности ротора L2΄ = L2n2 и взаимной индуктивности между статором и ротором L12΄ = L12n примерно равны максимальной взаимной индуктивности Lm, только L2΄ отличается от Lm на индуктивность рассеивания:

L'

L'

 

 

3

L

 

,

 

 

L'12

Lm .

 

 

 

 

 

 

2

 

2

2

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В приведенной системе уравнений обозначим

LM = 3/2Lm, тогда

 

d

1

 

R

 

 

i

 

 

 

u

,

 

 

 

 

dt

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

1

 

R

 

i

 

 

 

u

,

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

'2

 

 

 

 

 

 

 

'

2 R'2 i'2

0,

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

'2

 

 

 

 

 

 

 

'

 

R'

 

i'

 

0. ,

 

dt

 

 

 

 

 

 

 

2

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

L1

 

 

 

 

LM i1

 

 

LM i'2 ,

1

 

 

L1

 

 

 

 

 

LM i1

 

 

LM i'2 ,

'2

 

(L'2

 

 

 

LM )i'2

 

LM i1 ,

'2 (L'2 LM )i'2 LM i1,

J

d

 

M

 

M

 

,

 

 

dt

д

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mд

3

LM

i1 i'

2 i1 i'

.

 

 

2

2

 

 

 

 

 

 

 

 

 

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

p

u i u i ,

а мощность трехфазного двигателя:

pABC uA iA

 

 

u Bi B

uC iC

u i

 

1

 

u

 

3

u

1

i

 

3

i

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

u

3

u

1

i

3

i

3

 

(u i

 

u i ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

2

2

2

 

 

 

 

 

 

 

т.е. как и в формуле для момента, мощность трехфазного двигателя в 1.5 раза больше мощности двухфазного двигателя. Поэтому можно рекомендовать полученную систему уравнений, но помнить, что потребляемая двигателем мощность должна быть в 1.5 раза больше, чем рассчитываемая, а механическая мощность на валу ωМ определяется верно. Такая модель физически не верна, т.к. в двигателе не выполняется закон сохранения энергии: потребляемая мощность меньше выдаваемой. Второй способ: уменьшить момент на валу в 1.5 раза, момент сопротивления в 1.5 раза и момент инерции в 1.5 раза. Механическая мощность так же уменьшится в 1.5 раза и к.п.д. двухфазного двигателя будет таким же, что и трехфазного, т.е. двухфазный двигатель составляет 2/3 трехфазного.

Поэтому, и в том и другом случае, если двухфазная модель включается в более сложную схему в виде части (макроэлемента), следует в формулах перехода от двухфазной системы координат к трехфазной вводить коэффициент 3/2.

Врассмотренной модели двигатель имел только одну пару полюсов р

=1. Если число пар полюсов р = 2, то поворот ротора на 180о, по отношению к взаимной индуктивности, равносилен повороту на 360о, т.е. полному повороту. Поэтому в трехфазной модели перед углом γ следует поставить множитель p. Если обозначить скорость ротора Ω, то

d , dt

во всех же формулах, где были производные косинусов или синусов, появится множитель р.

d

'2

 

 

p

'2

R'2 i'2

 

 

0,

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

'2

 

 

p

'

 

R'

 

i'

 

 

0. ,

 

dt

 

2

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

d

 

 

M

 

M

,

 

 

 

 

 

 

 

 

 

д

 

 

 

 

 

 

 

dt

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mд

 

 

 

3p

LM

 

i1 i'2

 

 

i

1 i'

.

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

В трехфазной модели вместо угла γ должно быть рγ, а в формуле для момента вместо L12 должно быть pL12.

2.1.3 Статика асинхронного двигателя

Один из самых сложных вопросов математического моделирования - это определение параметров модели. Данные, приводимые в справочниках или даже заводские, очень часто не сходятся между собой. Поэтому следует проверять хотя бы установившееся состояние.

Будем считать, что в двухфазной модели асинхронного двигателя в установившемся состоянии напряжения и токи изменяются по гармоническому закону, например, для токов статора:

 

 

.

 

 

.

i1 Jm 2 I1 e j t , i1

Jm 2 I1 e j t ,

где Jm – мнимая часть комплексного тока, причем, сдвиг во времени между токами 90о, ω – частота питающего напряжения. В виду симметрии, модули токов i, iравны между собой и на комплексной плоскости в сумме дадут вектор İ1. Аналогично напряжения u, uна комплексной плоскости образуют вектор Ù1, а токи i, i– ток İ2. Таким образом, подставив токи и напряжения в уравнения:

d

1

R

i

 

u

,

 

 

1

dt

1

 

 

 

 

 

 

 

 

d

1

R

i

 

u

,

 

 

1

 

 

dt

1

 

 

 

 

 

 

 

 

1

 

L1i1

 

LM i'2 ,

1

 

L1i1

 

 

LM

i'2 ,

на комплексной плоскости получим:

.

.

.

.

R1 I1

j L1 I1

j LM I'2

U1 ,

.

.

.

.

R1 I1

j L1 I1

j LM I'2

U1

сложим правые и левые части, тогда:

 

 

.

 

 

 

 

.

 

 

 

 

 

 

.

 

 

.

 

 

 

 

R1 I1 j L1 I1 j LM I'2

 

U1 .

 

 

Аналогично поступим с уравнениями

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

'2

 

p

 

'2

 

R'2 i'2

0,

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

'2

 

p

 

'

 

 

R'

 

i'

 

0. ,

 

 

 

 

 

dt

 

 

2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'2

 

 

L'2 i'2

 

LM i1 ,

 

 

 

 

 

'

2

 

 

L'

i'

2

 

L

M

i

1

,

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

.

 

 

 

 

 

.

 

 

 

 

 

 

.

 

.

 

 

R2'

I'2

j L'2 I'2

 

 

j LM I1

 

p L'2 I'2

p LM I1

0

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

.

 

 

 

 

.

 

 

 

 

 

 

.

 

.

 

 

R2'

I'2

j L'2 I'2

 

j LM I1

 

p L'2 I'2

 

p LM I1

0,

 

сложив эти два уравнения, получим:

 

.

 

 

.

 

 

 

.

 

.

 

.

 

 

 

. .

R'

I'

2

j L'

I'

2

j L

M

I1

p L'

(I'

2

I'

2

) p L

M

(I1 I1 ) 0 .

2

 

2

 

 

 

2

 

 

 

 

Учтем, что мгновенное значение uα = Um sin ωt, а uβ = -Umcos ωt, то в комплексной форме

.

.

. . .

,

U

U, U

jU , U1 U U

аналогично и для токов, поэтому

.

.

.

.

 

.

 

.

 

. .

I1

I1

jI1

I'

2

I'

2

jI'

2

 

 

 

 

 

 

 

Кроме того, скольжение

1s,

/p

следовательно

R' .

 

 

.

 

 

 

.

2

I'

 

j L'

I'

 

j L

 

I1 0.

 

2

2

M

s

 

2

 

 

 

 

 

 

 

 

 

 

 

Чтобы получить уравнение для момента, учтем, что в установившемся состоянии токи i′и i′сдвинуты на один и тот же угол θ относительно токов iи i, т.е.

i1

 

Im1 sin(

t

)

 

i1

 

Im1 cos(

t

)

,

i'

I'

 

 

 

 

 

i'

I'

 

 

 

 

 

sin(

t

 

)

 

cos(

t

 

)

2

m2

 

 

 

 

2

m2

 

 

 

 

тогда

M

 

3p

L

 

i

 

i'

 

i

 

i'

 

3pL

I

I'

sin .

д

 

M

1

2

1

2

 

2

 

 

 

 

 

 

M 1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример: Асинхронный короткозамкнутый двигатель имеет следующие номинальные параметры: мощность на валу P2 = 5.5 кВт; напряжение

питания U = 220 В;

к.п.д. двигателя η = 0.855; скорость n = 1445

об/мин; cosφ = 0.85;

число пар полюсов р = 2;

кроме того,

известны ко-

эффициенты рассеивания статорной ζ1 = 1.015 и роторной

ζ2 = 1.03

обмоток; частота сети f = 50 Гц.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определить:

 

 

 

частоту сети в рад./с ω = 2πf ;

 

 

 

ω

= 100π;

 

скорость двигателя в рад./c

 

Ω = πn/30;

Ω = 151.3

 

скольжение s = 1 - pΩ/ω;

 

 

 

s

= 0.037

 

 

номинальный момент в Н. м

 

Mc = P2/Ω;

 

 

Mc = 36.35

 

номинальный ток фазы трехфазного

 

 

 

 

 

 

 

двигателя в А I1 = P2/(3ηU1 cosφ);

 

 

I1 = 11.47

 

Для расчета параметров двигателя

R1, R2, L1, L2, Lm комплексные уравне-

ния статики

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

.

 

 

 

 

.

 

.

 

 

R1 I1 j L1 I1 j LM I'2

U1 ,

 

 

 

R' .

 

 

.

 

 

 

 

.

0 ,

 

 

2

 

I'

 

j L'

I'

 

j L

 

I1

 

 

 

 

 

2

2

M

 

 

 

s

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

запишем в проекциях на действительную и мнимую оси:

R I

1a

L I

1b

L

I'

U

,

1

1

 

M 2b

1

 

R I

1b

L I

L

I'

0,

 

1

1

1a

 

M 2a

 

 

R2'

 

'

'

 

'

 

 

 

 

I

2a

L2 I

2 b

LM I1b

0,

s

 

 

 

 

 

 

 

 

R2'

 

'

'

 

'

 

 

 

 

I

2 b

L2 I

2a

LM I1a

0,

s

 

 

 

 

 

 

 

 

кроме того, введем уравнение выходной мощности:

3 (1 s) L

M

(I

1b

I'

I

1a

I'

) P .

 

 

2a

 

2 b

2

Это уравнение получается из уравнения для момента

P

M

д

3 (1 s) L

M

I

I'

sin(

1 2

) ,

2

 

 

1

2

 

 

где φ1, φ2 – фазы токов статора и ротора, представив синус разности через синусы и косинусы углов, а затем через отношения сторон треугольников комплексных токов. Ниже приведена страница расчета параметров асинхронного двигателя в программе MathCAD:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

p

 

 

 

 

 

 

2

 

 

 

 

 

U1

220

P2

5500

kw

.85

0.855

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1445

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

151.320046

 

 

s

1

 

p

 

 

 

 

 

 

s

0.036667

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I1

 

 

 

 

 

P2

 

 

 

 

 

 

 

I1

11.466575

 

 

 

Mc

 

P2

 

 

Mc

36.346804

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 kw U1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1

1.32

 

 

 

R2

.957

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

.169

 

 

 

 

L2

.1715

 

 

Lm

.1665

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

L1

2

 

 

L2

 

 

 

 

Xm

 

 

 

 

Lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lm

 

 

Lm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I1a

 

I1 kw

 

 

 

 

 

 

 

I1b

 

 

 

I1

1

 

 

 

 

kw2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I2a

9

 

 

 

 

 

 

 

 

 

 

 

 

 

I2b

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Given

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1 I1a

 

 

 

1 Xm I1b

 

 

Xm I2b

 

 

U1

 

 

R1 I1b

 

 

 

 

 

1 Xm I1a

 

 

 

 

Xm I2a

 

0

 

 

 

 

 

 

 

 

 

 

 

R2

 

I2b

 

 

 

 

2 Xm I2a

 

 

 

 

Xm I1a

 

 

0

 

R2

I2a

 

 

 

2 Xm I2b

 

 

Xm I1b

 

0

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 1

 

 

s ) Xm ( I1b I2a

 

 

 

I1a I2b)

 

 

P2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

Find( I2a I2b Xm R1 R2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9.566496

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.174982

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

33.935754

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.833975

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.762229

 

 

 

 

 

 

 

 

 

 

 

 

Lm

 

 

 

A2

 

L1

 

 

1 Lm

 

 

 

 

 

 

 

 

 

L2

 

 

 

2 Lm

 

R1

A3

R2

A4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1

1.833975

 

 

R2

0.762229

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lm

0.108021

 

 

L1

0.109643

 

 

 

 

L2 0.111265

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

2.1.4 Пуск асинхронного двигателя

Для решения в программе MathCAD следует всю систему уравнений привести к стандартной форме:

dx

f (x, t) .

 

dt

 

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

 

 

 

 

 

 

 

x1 = ω;

x2 = Ψ;

 

x3 = Ψ;

x4 = Ψ;

x5 = Ψ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

o

100

 

U1

310

 

p

2

 

P2

5500

Mc

36.347

J

.20188

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1

1.833975 R2

.762229

L1

0.109643

L2 0.111265

Lm

0.108021

 

 

 

 

 

0

 

A

1.5 p

Lm

B

 

L1

Lm

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

J

 

Lm L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

B1 1 x3

 

 

 

B1 2 x5 B2 1 x2

 

 

 

B2 2 x4 B1 1

x2

 

 

 

B1 2

x4 B2 1 x3

 

 

 

B2 2

x5

 

 

Mc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

D( t x)

 

 

 

 

 

 

R1 B1 1 x2

 

 

B1 2 x4

 

 

 

U1 sin(

o t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1 B1 1 x3

 

 

 

B1 2 x5

U1 c os( o t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1 x5 p R2 B2 1 x2 B2 2 x4 x1 x4 p R2 B2 1 x3 B2 2 x5

Z rkfixed( x013000D)

Решение ведется на интервале времени от 0 до 1 с, число точек при фиксированном шаге равно 3000.

Результаты расчетов:

i 1

3000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mgi

A J

 

 

 

 

B1 2 Zi 6 B2 1 Zi 3

 

 

 

B2 2 Zi 5

 

B1 1 Zi 3

 

 

 

B1 2 Zi 5 B2 1 Zi 4

 

 

 

B2 2 Zi 6

 

B1 1 Zi 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mg3000 36.346 Z3000 2 151.27

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Mg – момент двигателя в Н. м;

 

 

 

Zi,2 – скорость в рад./с

i

1

3000

 

 

 

 

 

 

 

 

 

 

Mg i

Графики переходного процесса:

Zi 5

B1 1 Zi 3

B1 2 Zi 5

B2 1 Zi 4

B2 2 Zi 6

A J B1 1

Zi 4

B1 2

Zi 6

B2 1 Zi 3

B2 2

 

 

2 00

 

 

 

 

 

 

 

 

 

 

 

 

1 50

 

 

 

 

 

 

 

 

 

 

 

 

Zi 2

 

 

 

 

 

 

 

 

 

 

 

 

1 00

 

 

 

 

 

 

 

 

 

 

 

 

Mgi

 

 

 

 

 

 

 

 

 

 

 

 

5 0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0 .2

 

0 .4

 

0 .6

0 .8

 

1

 

 

 

 

 

 

 

 

Zi 1

 

 

 

 

Рис. 1.1.4

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

Для этого представим уравнения асинхронного двигателя в виде:

 

 

 

 

di

1

 

 

 

 

 

 

di'

 

 

 

 

 

u ,

 

L

 

 

 

 

 

L

 

 

 

 

 

2

 

 

 

R i

 

 

 

 

dt

 

 

 

 

 

dt

 

 

 

1

 

 

 

 

 

 

M

 

 

1 1

 

 

 

 

 

 

 

 

di

1

 

 

 

 

 

 

di

'

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

L

 

 

 

 

 

2

 

 

 

R

 

i

 

 

u ,

 

 

 

dt

 

 

 

 

 

dt

 

 

 

 

 

1

 

 

 

 

 

 

M

 

 

1

 

1

 

 

 

 

L'

di2'

 

 

 

L

 

 

di1

 

R'

i'

e

 

0 ,

dt

 

 

 

dt

 

 

2

 

 

 

 

M

 

 

 

2 2

 

 

X

 

 

 

 

 

 

di'

 

 

 

 

 

 

 

di

1

 

 

 

 

 

 

 

 

 

 

 

 

L'

 

 

 

2

 

 

L

 

 

 

 

 

 

 

 

R

'

i'

e

 

 

0 ,

 

 

dt

 

 

 

 

dt

 

 

 

 

2

 

 

 

 

 

 

 

M

 

 

2 2

 

Y

 

e

X

 

 

 

 

 

 

p

 

(L'

2

i'

L

 

 

i

1

) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

M

 

 

 

 

e

Y

 

 

 

 

 

 

p

 

(L'

i'

L

 

 

i

) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

M 1

 

 

 

 

J

d

 

 

 

Md

 

 

 

MC ,

 

 

 

 

 

.

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

d

 

15.

p (i

 

 

i'

i

 

 

i'

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

1

 

2

 

 

 

 

На схеме обмотки статора 1α и 1β индуктивно связаны посредством одинаковых трансформаторов К1 и К2 с обмотками ротора 2α и 2β. Нелинейные источники ЭДС ex и ey выдают напряжение равное произведению скорости на поток сцепления. Уравнения динамики представлено ветвью с индуктивностью J и нелинейным зависимым источником Md. В текстовой части токи записаны непосредственно с указанием узлов.

В программе Simulink (MATLAB) в пакете Power System - моделирование электротехнических и энергетических систем - встроена модель асинхронного двигателя. Модель двухфазная с преобразованием трехфазных напряжений в двухфазные и двухфазных токов в трехфазные. На вход модели подается трехфазная система напряжений. Вводятся следующие параметры двигателя (обозначения программы):

Pn

- номинальная мощность (В.А),

Vn - фазное напряжение (В),

fn

- частота (Гц),

Rs

- сопротивление статора (Ом),

Lls

- индуктивность рассеяния статора (Гн),

R'r

- сопротивление ротора приведенное к статору (Ом) ,

L'lr - приведенная индуктивность рассеяния ротора (Гн), Lm - взаимная индуктивность (Гн),

J - суммарный момент инерции (кг.м2) ,

F - суммарный вязкий коэффициент трения (Н.м.с). p - число пар полюсов.

В последней строке вводятся начальные условия: s - начальное скольжение,

θe - электрический угол (град), isa,isb,isc - начальные токи фаз статора (А), pha,phb,phc - фазовые углы (град.).

Входными сигналами модели двигателя являются: напряжения фаз АВС (В) и момент сопротивления на валу двигателя (Н.м) - Tm.

Выходные величины выводятся с помощью блока Machines Measurement:

ira, irb, irc

- токи фаз a,b,c ротора (А),

irq, ird

 

- токи двухфазной модели ротора (А),

phirq, phird

- магнитные потоки ротора по осям q,d (Вб),

vrq, vrd

- напряжения на фазах ротора двухфазной модели (В),

isa, isb, isc

- токи фаз a,b,c статора (А),

isd, isq

 

- токи статора двухфазной модели (А),

phisq, phisd

- магнитные потоки статора по осям q,d (Вб),

vsq, vsd

- напряжения на фазах статора двухфазной модели (В),

wm

- скорость ротора (рад/с),

Te

- электромагнитный момент (Нм),

thetam

- угол поворота ротора (рад).

Переходный процесс полностью совпадает с рассчитанным ранее.

Трехфазная модель значительно сложнее, т.к. из-за изменений взаимных индуктивностей с углом поворота ротора, нельзя использовать обычную взаимную индуктивность, принятую в программе MicroCap. Вместо трансформаторов в схеме стоят нелинейные источники А1, В1, С1 и А2, В2, С2, моделирующие изменение взаимной индуктивности с углом пово-

рота, и линейные, которые дают максимальную взаимную индуктивность Lm. Сопротивления R1a, R1b, R1c и R2a, R2b, R2c – активные сопротивления статора и ротора (они указаны разными, т.к. они присутствуют в текстовом описании). Интегрирование скорости осуществляется лапласовским источником LF равным 1/s. Переменные х1, х2, х3 и y1, y2, y3 являются промежуточными, они несколько ускоряют вычисления.

Текстовое описание схемы:

. d efi n e p 2

. define y 2*Pi/3

. define Fi 2*v(out)

. define x1 c os (Fi)

. define x2 c os (Fi+y)

.define x3 cos(Fi-y)

.define y1 sin(Fi)

.define y2 sin(Fi+y)

.define y3 sin(Fi-y)

.define A2 i(R2a)*x1+i(R2b)*x2+i(R2c)*x3

.define B2 i(R2a)*x3+i(R2b)*x1+i(R2c)*x2

.define C2 i(R2a)*x2+i(R2b)*x3+i(R2c)*x1

.define A1 i(R1a)*x1+i(R1b)*x3+i(R1c)*x2

.define B1 i(R1a)*x2+i(R1b)*x1+i(R1c)*x3

.define C1 i(R1a)*x3+i(R1b)*x2+i(R1c)*x1

.define Mom - p*Lo12*(i(R1a)*(i(R2a)*y1+

+i(R2b)*y2+i(R2c)*y3)+i(R1b)*(i(R2a)*y3+ +i(R2b)*y1+i(R2c)*y2)+i(R1c)*(i(R2a)*y2+ +i(R2b)*y3+i(R2c)*y1))

.MODEL 3PHASEA SIN (F=50 A=310)

.MODEL 3PHASEB SIN (F=50 A=310 PH=-y)

.MODEL 3PHASEC SIN (F=50 A=310 PH=y)

Результаты моделирования лучше сравнить сразу по двум моделям: трехфазной и двухфазной. Как видно из рис. 2.1.5 различие довольно значительное.

Рис. 1.1.5

2.2 Моделирование синхронного генератора

2.2.1 Переход от осей А, В, С к вращающимся осям d, q

Пусть вектор индукции В создается тремя катушками А, В, С (рис.2.

2.1 а)

B BA BB CC ,

в каждой из которых индукция пропорциональна току, следовательно, если токи синусоидальны во времени, то синусоидальными будут и векторы индукции:

BA Bm sin( t); BB

Bm sin( t

2

) ; BC

Bm sin( t

2

) . (2.2.1)

3

3

 

 

 

 

 

Рис. 2.2.1

Найдем проекцию векторов BA ,BB ,BC на некоторую плоскость (рис.2.2.2)

B BA cos BB cos(

2

) BC cos(

2

)

3

3

 

 

 

Рис. 1.2.2

Рис. 2.2.3

(во втором слагаемом знак у аргумента cos изменен на противоположный, т.к. функция cos – четная). Подставляя значения ВA, BB, BC и учитывая,

что sinx cosy = = ½[sin(x+y) + sin(x-y)], получим:

B Bm sin

t cos

Bm sin( t

 

2

) cos(

 

2

)

 

Bm sin( t

2

 

)

cos(

2

)

3

3

3

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bm[sin( t

)

sin(

 

t

)

sin(

t

 

 

 

4

)

sin(

t

 

)

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin(

t

 

4

 

)

 

sin(

t

)]

 

3

Bm sin( t

),

 

 

 

 

 

3

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если угол

является функцией времени

=

 

t +

o, то индукция не ме-

няется во времени, причем максимальное значение равно 1.5Bm, т.е. суммарный вектор индукции вращается с угловой скоростью . Точно такое же поле можно получить и от двух катушек (рис.1.1b), как неподвижных –

координаты ,

, так и вращающихся вместе с ротором – координаты d,q.

В координатах

d,q

вектор индукции неподвижен. Аналогично предыду-

щему выводу, найдем (рис.2.2.3) (угол

отсчитывается от вектора индук-

ции фазы А):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

B

cos

B

 

cos(

2

 

)

B

 

 

cos(

2

) ,

(2.2.2a)

 

B

 

 

 

 

C

 

 

 

d

 

A

 

 

 

3

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

учитывая, что ось q

сдвинута на /2

относительно оси d , то cos( + /2) =

-sin , следовательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

[B sin

B

 

sin(

2

)

B

 

 

sin(

 

2

)] .

(2.2.2b)

B

 

 

 

 

C

 

 

q

 

A

 

 

 

3

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что при подстановке (2.2.1) в (2.2.2) также получится коэффициент 3/2. Преобразование (2.2.2) можно перенести на токи iA, iB, iC , потоки сцепления A, B, C , напряжения и т.д., однако принято в формулах преобразований учитывать коэффициент 2/3, чтобы получалось, что максимальное значение, например, тока равно

Imax = Im,

тогда

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

id

 

 

[iA cos

iB cos(

 

 

 

)

 

iC cos(

 

 

 

 

)]

 

 

3

 

3

 

3

 

 

 

i

 

 

2

[i

 

sin

i

 

sin(

 

2

)

i

 

sin(

 

2

)] ,

 

(2.2.3)

q

3

A

B

3

C

3

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iA Im sin(

 

 

t);

iB

Im sin( t

 

 

 

2

 

) ;

iC

I m sin( t

2

) .

 

 

 

 

 

3

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обратный переход от осей d,q к A,B,C получаем из рис.2.2.3

BA

Bd cos

Bq sin

 

 

 

 

 

 

 

2

 

 

 

2

(2.2.4)

BB

Bd cos(

)

Bq sin(

)

3

3

 

 

 

 

 

 

 

BC

Bd cos(

 

2

 

)

Bq sin(

2

 

)

3

 

3

 

 

 

 

 

 

 

 

Аналогично можно записать преобразования токов, потоков сцепления и т.д.

2.2.2 Уравнения неявнополюсной синхронной машины

Уравнения неявнополюсной синхронной машины (без демпферной обмотки):

uA

r

iA

d

A

 

 

 

;

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

uB

r

iB

d

B

;

(2.2.5)

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

uC

r

iC

d

C

 

 

;

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

r

i

 

 

d

f

 

;

 

f

f

 

 

 

 

 

 

f

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где r – активное сопротивление статорной обмотки, одинаковое по всем фазам;

rf - активное сопротивление роторной обмотки. Для потока сцепления, например, фазы А

A = LАА iA + LAB iB + LAC iC + LAf if ,

где LАА - собственная индуктивность фазы обмотки статора; LАВ , LAC - взаимная индуктивность между обмотками статора; LAf – взаимная индуктивность между обмоткой А статора и ротора. Для остальных потоков сцепления запись аналогичная. Если учесть, что взаимные индуктивности статорных обмоток равны половине от максимальной взаимной индуктивности Lm между катушками фаз А, В, С статора, т.к. они расположены под углом 120о (примерно такая же максимальная взаимная индуктивность будет и между катушками статора и ротора, приведенного к статору (см. далее)), то

 

 

 

LAf

= Lafcosγ.

 

 

 

 

Кроме того

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iA + iB + iC = 0,

 

 

 

следовательно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

1

L

 

i

 

L

 

cos

i

,

A

AA

 

 

m

A

af

 

2

 

 

 

 

 

f ,

 

 

 

 

 

 

 

 

 

 

 

 

где γ – угол между катушкой фазы А статора и ротора, Laf – максимальная взаимная индуктивность между обмотками статора и ротора. Собствен-