CC BY
450
УДК 62-83-523
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ АСИНХРОННОГО ДВИГАТЕЛЯ В НЕПОЛНОФАЗНОМ РЕЖИМЕ РАБОТЫ
Г.И. Однокопылов, А.Д. Брагин
Томский политехнический университет E-mail: OGIz@yandex.ru
Рассмотрен процесс формирования математической модели асинхронного двигателя в составе электроприводов опасных производственных объектов в неполнофазном аварийном режиме работы с учетом полей прямой и обратной последовательности, что позволяет учитывать влияние эллиптического поля в неполнофазных аварийных режимах работы. Процесс переключения в аварийный двухфазный режим работы выполнен с использованием матрицы отказов двигателя. Приведены результаты моделирования переходных процессов при обрыве фазы статора двигателя. Выполнено сравнение временных диаграмм, полученных экспериментально, по частоте вращения и току при помощи автоматизированного испытательного стенда и при помощи имитационной модели в среде Matlab Simulink. Показана адекватность предложенной математической модели, позволяющей исследовать асинхронный электродвигатель в аварийном двухфазном режиме работы трехфазного электродвигателя.
Ключевые слова:
Математическая модель асинхронного двигателя, неполнофазный режим работы, поле прямой и обратной последовательности, матрица отказов двигателя, обрыв фазы двигателя.
Исследование аварийных состояний асинхронного двигателя в составе исполнительных электроприводов опасных производственных объектов предполагает рассмотрение работы двигателя в неполнофазных режимах работы. Обеспечение отказоустойчивого управления со свойством живучести асинхронного электропривода предполагает разработку алгоритмов восстановления работоспособности с полным или частичным восстановлением работоспособности с круговым вращающимся полем.
При рассмотрении двухфазного режима работы трехфазного асинхронного двигателя можно выделить два различных режима:
• работа в аварийном режиме с эллиптическим полем в нерегулируемом электроприводе;
• работа в аварийном режиме с алгоритмическим формированием кругового вращающегося поля в частотно-регулируемом электроприводе [1]. Целью настоящей статьи является формирование математической модели асинхронного двигателя, учитывающей влияние эллиптического поля в неполнофазных режимах работы.
При рассмотрении аварийного двухфазного режима работы двигателя с эллиптическим полем в воздушном зазоре его можно разложить на поля прямой и обратной последовательности, используя метод симметричных составляющих. Обобщенная модель электрической машины учитывает лишь поле прямой последовательности. Существует два варианта представления модели асинхронного двигателя, работающего с эллиптическим полем -модель обобщенной электрической машины с двумя статорами и двумя роторами, модель обобщенной электрической машины с двумя статорами и одним ротором [2]. Первая модель - более проста в описании, однако дает не точные результаты моделирования переходных процессов, вторая же модель имеет более сложные выражения для результирующего электромагнитного момента, но дает
более точные результаты расчетов в динамических режимах работы.
В связи с вышеперечисленными особенностями математического описания обобщенной электрической машины с эллиптическим полем для формирования математической модели выбрана модель с двумя статорами и двумя ротором (рис. 1).
°^ М| х°
Рис. 1. Схема электрической машины с двумя статорами и двумя роторами
При исследовании несимметричных машин с учетом пространственных гармоник или при наличии несинусоидальных питающих напряжений приводить машину к двухфазной неправомерно, так как поля в зазоре двухфазной и многофазной машин в этих условиях отличаются друг от друга [2]. Так как асинхронный двигатель при обрыве одной из фаз питающей сети является несимметричным - в качестве математической модели может быть принята обобщенная электрическая машина в трехфазной системе координат.
При составлении уравнений и рассмотрении переходных процессов асинхронных машин использованы общепринятые допущения и ограничения, связанные с понятием «идеализированная машина»: машина не насыщена, потерями в стали пренебрегаем; фазные обмотки симметричны и сдвинуты на угол 120 электрических градусов; намагничивающие силы обмоток и магнитные поля распределены синусоидально, вдоль окружности воздушного зазора; воздушный зазор равномерен; ротор симметричен. Реальная распределенная обмотка заменяется сосредоточенной, оси фаз которой сдвинуты на фазный угол, а намагничивающая сила ее равна намагничивающей силе реальной обмотки.
В качестве основы были взяты уравнения модели асинхронного двигателя в заторможенной системе координат [1]. Для упрощения расчетов уравнения для полей прямой и обратной последовательности представлены отдельными системами.
При определении напряжений, подаваемых на обмотки машин, отвечающих за поля прямой и обратной последовательности, учитывалось, что при постоянном питающем напряжении и 1=иА=ив=ис=сопв1 напряжения прямой и обратной последовательности зависят от скольжения [3].
Уравнения обобщенной электрической машины имеющей два статора и один ротор для поля прямой последовательности имеют вид:
А Т,А
Уравнения для обратной последовательности:
А Т л
А А Т,
и1В — КА* + ;
А А Т
и,с — ЯАс +- 1С ■
А
АТ I-
0 = КАа +~!Т+(т^-т, )ю/^
АТ
о = я/, + —^+(Тс -Та)® ^^3; А
АТ
о = яа + —+(Та—т, )®/Тз. А
Потокосцепления для всех фаз:
Т1А — АбКа — 2^тЧь — 2АтЬ + +АтКа — 2 Ал^й — 2 АЛ ■
Т1В — — 2Ат\л + -^18 — 2Ал^С —
— 2 АтКа + Ал^Ь — . АЛ ■
Т1С — — 2 Аяг']А — 2Ая^В + А?АС —
— 2АтКа — 2Аи^й + АЛ ■
Т1а — АтКл — 2Ат/1В — . —т А1С +
+АяКа — 2 Аи^й — 2 'Ая*! ■
Т1Ь — — 2 Ат\л + -тА1В — .-тА1С —
— 2 АтКа + —яг1А — ^ 'АЛ ■
Т1с — — 2АтА1Л — ^-тА1В + Ая^С —
— 2АтКа — 2Аи^й + А;*, •
и2 А ЯАА +
и2 В ЯлА2 В +
и 2 С Я!А2С +
А/
АТ
Ы Т 2;
А/
АТ2
А/
0 - ЯгА2а +
о — я/,.,, +
0 — ЯАс +
А Т., А/
А Т
а Т21
А/
АТ2 ____2
А/
+ (Т2, -Т.с )®/л/3; + (Т.с -Т.а )®/>/3; + (Т 2а -Т 2, )т/£.
Потокосцепления для всех фаз:
Т2А — АБ/2Л — 2 —тЧв — . Цп^Т +
Т ■ —1Т ■ —1Т ■ ■
+ т'.а 2 тА2Ь . тА2с ■
Т2 В — — 2 Ат'/2А + — 2 Ат^.^ —
— 2 Ат/2а + АтА2Ь — . АиАс ■
Т2С — — 2АиАл — 2-т^2В + А; А2С — + 2 -тАа — 2 -т ^2Ь + Ат ■ 2с ■
Т2а — ^тЧл — 2 Ат Чв — . Ат А2С +
+LR^2a — 2^т^Ы — . Ат А2с ■
Т24 — — 2 ^т^А + АтЬВ — . Ат А2С —
— 2 ^тЧа + — 2 Ат А2с ■
Т2с — — 2 ^тЧл — 2 Ат Чв + Ат А2С —
1
1
2 ’Ат Ьа 2 АтА2Ь + ■ 2с •
Электромагнитные моменты Ы{ и М2, создаваемые полями прямой и обратной последовательности, определяются по следующим выражениям:
(/1А/1с + ■ 1ВА1а + ^С^Ь ) —
_—(А1 Аг1Ь + ^В^с + /1С/1а ) _
(/2АА2с + А2вЧа + ■2сЧь ) —
_—(/2 АА2Ь + ■ 2ВА2с + А2СА2а ) _
Уравнение движения электропривода, учитывающее моменты прямой и обратной последовательности:
м л, м Л Ша
М. + М2 - м = —------.
12 с т.
р Ш
Итоговые выражения для уравнений поля прямой последовательности:
= [ аг
[Ц]-[[Я ] + [ Ь2]~31 [і^,
для поля обратной последовательности:
Ж [/2
Ж
= [ АГ
[и2] -|[ Я ]+[ А2] [І 2
[/і]=
[иі] =
ии
ик
0
0
0
[І2] =
[и2] =
~и2А " " я, 0 0 0 0 0 "
и2 В 0 я, 0 0 0 0
и2 С ; [ я ] = 0 0 я, 0 0 0
0 0 0 0 Яг 0 0
0 0 0 0 0 Яг 0
0 0 0 0 0 0 Яг _
[ А] =
а5 - Ат /2 - А» /2 А»
- Ат/2 А - 1т /2 -А» /2
-Ат /2 -А» /2
Ат - Ат /2
А -Ат /2 -Ат /2 Ат
- Ат /2 - Ат /2 А - а» /2 - а» /2
-Ат/2 Ат -а»/2 -а»/2 А -А,, /2
Ат -Ат /2 -А» /2 А
- Ат /2 - Ат /2
- Ат /2 - Ат /2
[ а2] =
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 а - а 0 ь -ь
- а 0 а -ь 0 ь
а - а 0 ь - ь 0
а = 3 = 2 А», ь = ( АЯ + і 2 А» )
трицу отказов, позволяющую переходить в аварийный двухфазный режим работы:
"М~
[ М ] =
где [£х], [г2] - матрица неизвестных (токи статора и ротора прямой и обратной последовательностей); [Ці], [и2] - матрица напряжений прямой и обратной последовательностей; [Л] - матрица сопротивлений; И, И - матрицы индуктивностей; [Ь1]~1 -обратная матрица индуктивностей статора:
иі
Мв
мс 0 0 0
При этом нулевому состоянию переменных МА, Мв, Мс будет соответствовать номинальный трехфазный режим работы, а единичному - обрыв соответствующей фазы обмотки статора. Тогда итоговое выражение для матрицы активных сопротивлений примет следующий вид:
[ Ям ] = (1 + к06р[ М ])[ Я],
где кобр - коэффициент, имитирующий введение высокоомного сопротивления в цепь статора АД.
Окончательные выражения для уравнений полей прямой и обратной последовательности примут следующий вид соответственно:
Жіі]
Ж
Ж [І2]
ж
= [ Аіґ П^]-|[ Ям ]
[ а2] ^1 [/і]
= [АГ][и2]-| [Я
Имитация обрыва фазы статора осуществляется путем приравнивания коэффициента МА, Мв или Мс, соответствующего отказавшей фазе, нулю. При этом на выводы обмотки машины, отвечающей за поля обратной последовательности, подается напряжение, зависящее от текущего значения скольжения, и она подключается встречно с машиной прямой последовательности.
Расчет модели проводился в среде МаШЬ 81ши-Нпк. Адекватность предложенной математической модели была проверена путем сравнения расчетных и экспериментальных данных, полученных при помощи автоматизированного испытательного стенда (рис. 2). Результаты сравнения переходных процессов, рассчитанных на модели и полученных экспериментально, приведены на рис. 3, 4.
При работе в трехфазном режиме поля обратной последовательности в воздушном зазоре отсутствуют и соответствующие уравнения в расчетах не участвуют.
Для имитации обрыва фазы статора необходимо в исходные выражения для уравнений поля прямой и обратной последовательности ввести ма-
Рис. 2. Автоматизированный испытательный стенд для исследования аварийных и неполнофазных режимов работы асинхронного двигателя: 1 - резистивная нагрузка; 2 - компьютерная измерительная система М1С-300; 3 - источник питания; 4 - блок формирования аварийных режимов работы; 5 - асинхронный двигатель АИР63А2; 6 - машина постоянного тока серии 1ПИ 12.11
Переходные процессы при пуске двигателя, полученные
Переходные процессы при пуске двигателя, полученные с помощью
А, 1д, 1б, 1с
■? -‘-А; ^Б?
А, 1а, 1б, 1с
0 0.2 0.4 0.6 0.8
1 1.2 1, с
1 1.2 1, с
Рис. 3. Переходные процессы при пуске двигателя, полученные экспериментально и при помощи математической модели: а) по скорости; б) по току
Переходные процессы при обрыве фазы Переходные процессы при обрыве фазы
статора двигателя, полученные статора двигателя, полученные с
экспериментально помощью математической модели
рад/с, га рад/с, га
0.6 0.8 1 1.2 1.4 1.6 1.8 1, с 0.6 0.8
а
1 1.2 1.4 1.6 1.8 1 с
А, 1а, 1б, 1с 3
2 1 0 -1 -2 ■3
0.6 0.
А, 1д, 1б, 1с
1 1.2 1.4 1.6 1.8 1, с
б
Рис. 4. Переходные процессы при обрыве фазы статора двигателя, полученные экспериментально и при помощи математической модели: а) по скорости; б) по току
б
Разница между частотой вращения двигателя и амплитудой тока, полученных экспериментально и на математической модели, в трехфазном установившемся режиме работы составляет 2,8 и 9,6 % соответственно. В аварийном двухфазном режиме работы эта разница составляет 2,3 % для частоты вращения и 2,5 % для тока для установившихся значений. Отклонения в переходных режимах работы вызваны принятыми в ходе составления математической модели допущениями, а также несимметрией питающей сети при проведении эксперимента.
Таким образом, разработанная математическая модель позволяет исследовать асинхронный электродвигатель в аварийном двухфазном режиме работы с учетом поля обратной последовательности, а также в двухфазном режиме работы с алгоритмическим поддержанием кругового вращающегося поля, когда поле обратной последовательности в воздушном зазоре отсутствует.
Выводы
1. Разработана математическая модель асинхронного двигателя для исследования неполнофазного аварийного режима работы с учетом полей прямой и обратной последовательности с использованием матрицы отказов двигателя.
2. На основе разработанной имитационной модели в среде Matlab 8тиНпк и проведенных подтверждающих экспериментов для случая обрыва фазы статора асинхронного двигателя при эллиптическом вращающемся поле расхождение результатов по току составило 2,3 %, по частоте вращения - 2,5 % в установившемся режиме, в трехфазном установившемся режиме работы расхождение результатов по току составило 2,8 %, по частоте вращения, 6 %, что подтверждает адекватность разработанной математической модели.
^ИСОК ЛИТЕРАТУРЫ 2.
1. Однокопылов Г.И., Однокопылов И.Г. Обеспечение живучести электродвигателей переменного тока: монография. - Томск: Изд-во ТПУ, 2005. - 187 с.
Копылов И.П. Математическое моделирование электрических машин. - 3-е изд., перераб. и доп. - М.: Высшая школа, 2001. - 327 с. Иванов-Смоленский А.В. Электрические машины: в 2-х т. Т. 1. - 2-е изд., перераб. и доп. - М.: Изд-во МЭИ, 2004. - 652 с.
Поступила 20.05.2013 г.
UDC 62-83-523
MATHEMATICAL MODEL OF INDUCTION MOTOR IN PHASE FAILURE MODE
G.I. Odnokopylov, A.D. Bragin Tomsk Polytechnic University
The paper considers the formation process of mathematical model of induction motor in electric drives of hazardous production facilities in phase failure emergency mode considering direct and reverse sequence fields that allows considering the impact of the elliptic field in phase failure emergency modes. The switching process into emergency two-phase operation was made using motor failure matrix. The paper introduces the results of transient modeling at loss of motor stator phase. The timing charts of speed and current obtained experimentally were compared using an automated test bench and a simulation model in the Matlab Simulink. Adequacy of the proposed mathematical model which allows investigating the induction motor in emergency two-phase operation of three-phase motor is proved.
Key words:
Mathematical model of induction motor, phase failure mode, field of direct and negative sequence, the matrix of motor failures, loss of motor phase.
REFERENCES
1. Odnokopylov G.I., Odnokopylov I.G. Obespechenie zhivuchesti elektrodvigateley peremennogo toka (AC motor survivability support). Tomsk, TPU Publ., 2005. 187 p.
2. Kopylov I.P. Matematicheskoe modelirovanie elektricheskikh mashin (Mathematical modeling of electrical machines). Moscow, Vysshaya shkola, 2001. 2d Iss., 327 p.
3. Ivanov-Smolenskiy A.V. Elektricheskie mashiny (Electrical machines). Moscow, MEI Publ., 2004. 2dIss., 1, 652 p.
