Анализ системных свойств асинхронного электропривода. Монография
.pdfo |
|
|
|
|
|
|
|
|
1 |
|
o |
|
1 d i1o |
|
|
1 d i2o |
|||||||||||
u1 |
C1u R1 C1i |
|
i1 L1 C1i |
|
|
|
|
|
M12 |
C1i |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
dt |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
C |
|
1 |
|
|
|
|
|
C |
1 |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
d |
|
o |
|
|
|
d |
o |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
1i |
|
|
|
|
|
|
|
1i |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
1 L1 |
|
|
|
|
|
|
i1 |
1 M12 |
|
|
|
|
|
i2 |
; |
|
|
|
(2.81) |
||||||||
|
|
|
d 1 |
|
|
|
|
|
|
d 1 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 d i1o |
|
1 d i2o |
||||||||
o |
|
|
|
|
|
|
|
1 |
o |
|
|
|
|
|
|
||||||||||||
u2 |
C2u |
R2 C2i |
|
i2 M21 C2i |
|
L2 |
C2i |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
dt |
|
|
|
|
|
|
d C |
|
1 |
|
|
|
|
d C |
|
|
|
1 |
o |
|
|
|
|
|
|
|
||||
|
|
|
|
|
o |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
2i |
|
|
|
|
2i |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
2 M21 |
|
d 2 |
i2 2 L2 |
d 2 |
|
|
i1 . |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Перемножив матрицы и векторы в (2.81), получим следующую систему уравнений обобщенной машины:
u |
R i |
|
L |
|
|
di1d |
M |
|
|
|
|
di2d |
L i |
M |
i |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
1d |
11d |
|
1 dt |
|
|
|
|
m dt |
|
1 11q |
|
m |
2q |
|
||||||||||||||||
|
|
|
|
|
|
|
di1q |
|
|
|
|
|
|
|
di2q |
|
|
|
|
|
|
|
|
|
|||||||
u |
R i |
|
L |
|
|
|
M |
m |
|
|
L i |
M |
i |
2d |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
1q |
11q |
1 |
|
dt |
|
|
|
|
|
|
dt |
1 11d |
|
m |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
di2d |
|
|
|
|
|
|
|
di1d |
2 L2i2q Mmi1q |
||||||||||||||||
0 R2i2d |
L2 |
Mm |
|||||||||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
di2q |
|
|
|
|
|
di1q |
|
|
|
|
|
L i |
|
|
|
|
|
||||||||
0 R i |
L |
|
M |
|
|
|
|
|
|
M |
i |
||||||||||||||||||||
|
|
|
|
|
|
m |
|
2 |
2d |
||||||||||||||||||||||
|
|
2 2q |
|
2 |
|
|
|
dt |
|
|
|
|
|
dt |
2 |
|
m 1d |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
;
;(2.82)
;
.
Очевидно, что коэффициенты в уравнениях системы (2.82) являются постоянными величинами.
Пространственная модель обобщенной машины показана на рис. 2.6.
На основании (2.82) можем записать выражения потокосцеплений
фаз статора и ротора соответственно |
|
|
|
1d |
L1i1d Mmi2d ; |
1q L1i1q Mmi2q ; |
(2.83) |
2d |
L2i2d Mmi1d ; |
2q L2i2q Mmi1q . |
(2.84) |
Чтобы получить уравнение электромагнитного момента обобщенной машины, запишем выражение мгновенной потребляемой мощности в осях d, q:
p1 u1di1d u1qi1q . |
(2.85) |
51
Подставим в (2.85) вместо u1d и u1q первое и второе уравнения системы (2.82). Учитывая (2.83) можем записать [29]:
p1 R1 i12d i12q d 1d i1d dt
d 1q |
i |
i |
i |
. |
(2.86) |
|
|||||
dt 1q |
1 1d 1q |
1q 1d |
|
|
|
Первое слагаемое в правой части (2.86) представляет собой электрические потери в обмотке статора, второе и третье – электромагнитную энергию, запасаемую в обмотках статора, четвертое – электромагнитную мощность. Отметим, что электромагнитная энергия, запасаемая в обмотках статора, в формировании электромагнитной мощности не участвует.
Рис. 2.6. Пространственная модель обобщенной электрической машины
Таким образом, выражение электромагнитной мощности можно записать в одной из следующих форм:
pýì 1 1di1q 1qi1d ; pýì 1Mm i1qi2d i1di2q . (2.87)
Для определения электромагнитного момента выражение электромагнитной мощности необходимо поделить на 1 , в результате чего выражение электромагнитного момента обобщенной электрической машины можно записать в одной из следующих форм:
52
Mýo 1di1q 1qi1d ; |
(2.88) |
Mýo Mm i1qi2d i1di2q . |
(2.89) |
Таким путем может быть получена еще одна формула электромагнитного момента обобщенной электрической машины:
Mýo 2qi2d 2di2q , |
(2.90) |
которая после подстановки (2.84) приобретает вид (2.89).
Ранее было показано, что число полюсов обобщенной электрической машины равно единице, поэтому геометрическая частота вращения магнитного поля 1ã равна угловой частоте
питающего напряжения 1. Если многополюсный трехфазный АД
рассматривать с позиций теории обобщенной электрической машины, то справедливо следующее соотношение:
1ã 1 ðï 1ãò , |
(2.91) |
где 1ãò – геометрическая частота вращения магнитного поля многополюсного трехфазного АД.
С учетом (2.91) можем записать следующие выражения электромагнитного момента трехфазного АД:
Mý pï 1di1q 1qi1d |
; |
(2.92) |
||
Mý pï Mm i1qi2d i1di2q |
; |
(2.93) |
||
Mý pï 2qi2d 2di2q . |
|
(2.94) |
||
Запишем уравнение движения обобщенной машины: |
|
|||
J |
d |
pï Mý Mc . |
(2.95) |
|
|
||||
|
dt |
|
|
|
Отметим, что в уравнениях движения (2.95) и электромагнитного момента (2.92) – (2.94) обобщенной машины используются электромагнитный, статический, а также момент инерции, равные по величине аналогичным величинам трехфазного АД.
53
Таким образом, система (2.82) в совокупности с выражением электромагнитного момента (2.93) и уравнением движения (2.95) представляют собой математическую модель обобщенной машины в осях d, q без учета нелинейности магнитопровода и потерь в стали.
Рассмотренное математическое описание процессов электромеханического преобразования энергии в обобщенной машине позволяет определять токи в задачах динамики. Использование преобразований (2.61) – (2.63) и (2.64) – (2.66) позволяет перейти от токов и напряжений обмотки статора обобщенной электрической машины к токам и напряжениям обмотки статора трехфазной асинхронной машины.
Чтобы решить уравнения системы (2.82) численными методами необходимо привести их к нормальной форме. Подобное приведение может быть осуществлено за счет использования совместно с численным методом интегрирования системы дифференциальных уравнений (2.82) метода решения системы линейных алгебраических уравнений относительно производных. При этом систему уравнений (2.82) необходимо преобразовать так, чтобы в их левых частях находились производные искомых функций с соответствующими коэффициентами, а в правых частях – искомые функции.
Другой путь решения системы линейных алгебраических уравнений относительно производных связан с поиском обратной матрицы индуктивностей. Перенесем в уравнениях системы (2.82) производные токов с соответствующими коэффициентами в левую часть, а искомые функции – в правую. Сгруппируем уравнения системы так, чтобы уравнения баланса напряжений фаз статора и ротора, расположенных на одной оси, следовали друг за другом, получим:
L |
di1d |
|
M |
|
|
|
di2d |
|
u |
R i |
|
L i |
|
M |
|
i |
|
; |
||||||||||||||||
|
|
m dt |
|
|
m |
2q |
||||||||||||||||||||||||||||
1 dt |
|
|
1d |
|
|
|
11d |
|
|
|
|
1 11q |
|
|
|
|
|
|||||||||||||||||
|
|
|
di2d |
|
|
|
|
|
|
di1d |
|
|
|
|
|
|
|
|
|
|
L i |
|
|
|
|
|
|
; |
|
|||||
L |
|
M |
m |
|
|
|
R i |
2d |
|
2 |
2q |
M |
|
i |
(2.96) |
|||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
2 |
dt |
|
|
|
|
|
dt |
|
|
2 |
|
|
|
|
2 |
|
|
|
m |
1q |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
di |
|
|
|
di |
2q |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
L |
|
1q |
M |
|
|
|
|
|
|
|
u |
R i |
L i |
|
M |
|
i |
2d |
; |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
1 dt |
|
m dt |
|
1q |
|
|
11q |
|
|
|
|
1 11d |
|
|
m |
|
|||||||||||||||||
|
|
|
di |
|
|
|
|
|
|
di |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
L |
|
2q |
M |
m |
|
|
|
|
1q |
|
R i |
2q |
|
2 |
L i |
2d |
M |
i |
|
. |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
2 |
dt |
|
|
|
|
|
dt |
|
|
2 |
|
|
|
2 |
|
|
m |
1d |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
54
На основании системы (2.96) можем записать матрицу индуктивностей обобщенной машины в виде
|
|
|
L1 |
|
Mm |
|
|
0 |
0 |
|
|
|
|
|
|||
|
|
M |
|
|
|
L |
|
|
0 |
0 |
|
|
|
||||
|
Lo |
m |
|
|
|
|
|
|
Lod |
0 |
(2.97) |
||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
, |
||||
|
|
|
0 |
|
0 |
|
|
|
L1 |
Mm |
|
0 |
Loq |
|
|||
|
|
|
0 |
|
0 |
|
|
|
Mm |
L2 |
|
|
|
|
|||
|
Lod – матрица |
|
|
|
|
|
|
|
|
|
|||||||
где |
индуктивностей |
обобщенной машины по оси d; |
|||||||||||||||
Loq |
– матрица индуктивностей обобщенной машины по оси q; |
0 – |
|||||||||||||||
нулевая матрица размером 2 2. |
|
|
|
|
|
|
|||||||||||
Вид матрицы Lo |
|
позволяет записать обратную матрицу |
1 |
||||||||||||||
|
Lo |
||||||||||||||||
следующим образом [75]: |
|
|
|
Lo 1 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
1 |
|
0 |
|
|
|
|||||||
|
|
|
|
Lo |
|
|
|
d |
|
|
|
|
|
(2.98) |
|||
|
|
|
|
|
|
0 |
|
1 , |
|
||||||||
|
|
|
|
|
|
|
|
|
|
Loq |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Запишем систему (2.96) в векторно-матричной форме: |
|
||||||||||||||||
|
|
|
|
|
|
|
d |
|
|
|
|
1 |
 |
, |
|
|
|
|
|
|
|
|
|
|
io Lo |
|
(2.99) |
||||||||
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
где |
io – вектор токов обобщенной машины; |
 – вектор искомых |
|||||||||||||||
функций системы (2.99).
Отметим, что Lo является матрицей постоянных коэффициентов.
Кроме того, матрицы Lod |
и Loq равны, поэтому можно один раз |
вычислить обратную матрицу, например, Lod 1, а с ее помощью
записать и матрицу Lo 1. Подобный подход позволяет упростить алгоритм численного решения системы дифференциальных уравнений и сократить количество вычислительных операций.
55
Сиспользованием математического описания обобщенной машины
восях d, q без учета насыщения магнитопровода и потерь в стали проведено компьютерное моделирование электромагнитных и электромеханических процессов, а также расчет рабочих характеристик обобщенной машины на основе двигателя АИР80А6У2
всистеме Turbo Pascal 7.0.
Результаты компьютерного моделирования переходных и установившихся процессов и динамическая механическая характеристика обобщенной машины на базе двигателя АИР80А6У2 приводятся на рис. 2.7, 2.8 и 2.11 а.
Моделирование проводилось при условии, что вектор напряжений статора совпадает по направлению с осью d обобщенной машины, т.е.
u1d U1m ; |
u1q 0; |
1 314,15 рад/с. |
(2.100) |
Из рис. 2.7, 2.8 видно, что токи фаз обобщенной машины представляют собой достаточно плавно изменяющиеся функции, которые в установившемся режиме можно рассматривать как сигналы постоянного тока. Токи фаз АД, представленные на рис. 2.9, 2.10, изменяются по более сложным законам, а в установившемся режиме представляют собой гармонические функции. Сопоставить такие функции можно только в том случае, если к величинам одной из систем применить преобразование координат. Однако основным принципом преобразования координат является принцип инвариантности мгновенной мощности, поэтому вывод о корректности выполненных преобразований можно сделать на основании сравнения кривых мгновенной потребляемой мощности. Мгновенная потребляемая мощность трехфазного АД определяется выражением (2.35), а обобщенной электрической машины – выражением (2.85).
56
Рис. 2.7. Прямой пуск обобщенной электрической машины на базе двигателя АИР80А6У2 без нагрузки:
а– фазные токи статора; б – приведенные фазные токи ротора;
в– мгновенная потребляемая мощность, электромагнитный момент и частота вращения ротора
57
Рис. 2.8. Прямой пуск обобщенной электрической машины на базе двигателя АИР80А6У2 с номинальной нагрузкой:
а– фазные токи статора; б – приведенные фазные токи ротора;
в– мгновенная потребляемая мощность, электромагнитный момент и частота вращения ротора
58
Рис. 2.9. Прямой пуск двигателя АИР80А6У2 без нагрузки:
а– фазный ток статора; б – приведенный фазный ток ротора;
в– мгновенная потребляемая мощность, электромагнитный момент и частота вращения ротора
59
Рис. 2.10. Прямой пуск двигателя АИР80А6У2
сноминальной нагрузкой:
а– фазный ток статора; б – приведенный фазный ток ротора;
в– мгновенная потребляемая мощность, электромагнитный момент и частота вращения ротора
60
