Учебное пособие 1562
.pdf2.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 |
|||
|
|
|
|
|
Аналогично поступим с потоком Ψ1β, после вычислений имеем:
|
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 |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Аналогично вычисляется поток Ψ2β: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
L |
|
|
i |
|
|
3 |
L |
|
i |
. |
|
|
|
|
|||||
|
|
|
2 |
2 |
2 |
|
|
12 |
|
|
|
|
|||||||||||
|
|
|
|
|
|
2 |
|
|
|
|
1 |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Чтобы получить уравнение момента двигателя в координатах αβ, в формуле (2.1.4) сделаем следующие преобразования: первая квадратная скобка есть ток i2β умноженный на 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 – это отношение числа витков статора и ротора (как и у трансформатора, но в двигателях нужно учитывать еще и распределенность обмоток; в данном случае считаем, что соответствующие коэффициенты одинаковы для статора и ротора). Таким образом:
u’2 = nu2; i’2 = i2/n ; z’2 = 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о, ω – частота питающего напряжения. В виду симметрии, модули токов i1α , i1β равны между собой и на комплексной плоскости в сумме дадут вектор İ1. Аналогично напряжения u1α , u1β на комплексной плоскости образуют вектор Ù1, а токи i2α, i2β – ток İ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′2α и i′2β сдвинуты на один и тот же угол θ относительно токов i1α и i1β, т.е.
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 = Ψ1α; |
|
x3 = Ψ1β; |
x4 = Ψ2α; |
x5 = Ψ2β. |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
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 – максимальная взаимная индуктивность между обмотками статора и ротора. Собствен-