Научная статья на тему 'МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И ХАРАКТЕРИСТИКИ АСИНХРОННОГО ДВИГАТЕЛЯ ПРИ ПИТАНИИ ОТ ИСТОЧНИКА ТОКА'

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И ХАРАКТЕРИСТИКИ АСИНХРОННОГО ДВИГАТЕЛЯ ПРИ ПИТАНИИ ОТ ИСТОЧНИКА ТОКА Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
112
19
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АСИНХРОННЫЙ ДВИГАТЕЛЬ / ПИТАНИЕ ОТ ИСТОЧНИКА ТОКА / МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / НАСЫЩЕНИЕ МАГНИТОПРОВОДА / ВЫТЕСНЕНИЕ ТОКА / ПЕРЕХОДНЫЙ ПРОЦЕСС / СТАТИЧЕСКИЕ ХАРАКТЕРИСТИКИ

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Маляр В. С., Маляр А. В.

Разработаны методы и математические модели для исследования режимов и характеристик трехфазного асинхронного двигателя с короткозамкнутым ротором при питании обмотки статора от источника тока. Изложены особенности созданных на их основе алгоритмов расчета переходных процессов, установившихся режимов и статических характеристик. Приведено сравнение результатов расчета процессов и характеристик асинхронных двигателей при питании от источника тока и источника напряжения. Исследование установившихся и динамических режимов не может быть осуществлено с достаточной адекватностью на основе известных схем замещения и требует использования динамических параметров, являющихся элементами матрицы Якоби системы уравнений электромеханического равновесия. В математической модели уравнения состояния контуров статора и ротора составлены в неподвижной двухфазной системе координат. Переходные процессы описываются системой дифференциальных уравнений электрического равновесия преобразованных контуров двигателя и уравнением движения ротора, а установившиеся режимы - системой алгебраических уравнений. В основу разработанных алгоритмов положена математическая модель двигателя, в которой учитываются насыщение магнитопровода и явление скин-эффекта в стержнях короткозамкнутой обмотки. Для учета насыщения используются реальные характеристики намагничивания основным магнитным потоком и потоками рассеяния обмоток статора и ротора. На их основе вычисляются дифференциальные индуктивности, которые являются элементами матрицы Якоби систем уравнений, описывающих динамические режимы и статические характеристики. С целью учета скин-эффекта в обмотке короткозамкнутого ротора каждый стержень вместе с короткозамыкающими кольцами разбивается по высоте на несколько элементов. В результате в математической модели рассматриваются эквивалентные обмотки ротора с разными по значению параметрами, между которыми существуют взаимоиндуктивные связи. Решение нелинейной системы алгебраических уравнений электрического равновесия, которой описываются установившиеся режимы, выполняется методом продолжения по параметру. Для расчета статических характеристик используется дифференциальный метод в сочетании с итерационным уточнением методом Ньютона.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Маляр В. С., Маляр А. В.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MATHEMATICAL MODEL AND CHARACTERISTICS OF THE INDUCTION MOTOR WITH A POWER SUPPLY FROM A CURRENT SOURCE

. Methods and mathematical models for studying the modes and characteristics of the three-phase squirrel-cage induction motor with the power supplied to the stator winding from the current source have been developed. The specific features of the algorithms for calculating transients, steady-state modes and static characteristics are discussed. The results of the calculation of the processes and characteristics of induction motors with the power supply from the current source and the voltage source are compared. Steady-state and dynamic modes cannot be studied with a sufficient adequacy based on the known equivalent circuits; this requires using dynamic parameters, which are the elements of the Jacobi matrix of the system of equations of the electromechanical equilibrium. In the mathematical model, the state equations of the stator and rotor circuits are written in the fixed two-phase coordinate system. The transients are described by the system of differential equations of electrical equilibrium of the transformed circuits of the motor and the equation of the rotor motion and the steady-state modes by the system of algebraic equation. The developed algorithms are based on the mathematical model of the motor in which the magnetic path saturation and skin effect in the squirrel-cage bars are taken into consideration. The magnetic path saturation is accounted for by using the real characteristics of magnetizing by the main magnetic flux and leakage fluxes of the stator and rotor windings. Based on them, the differential inductances are calculated, which are the elements of the Jacobi matrix of the system of equations describing the dynamic modes and static characteristic. In order to take into account the skin effect in the squirrel-cage rotor, each bar along with the squirrel-cage rings is divided height-wise into several elements. As a result, the mathematical model considers the equivalent circuits of the rotor with different parameters which are connected by mutual inductance. The non-linear system of algebraic equations of electrical equilibrium describing the steady-state modes is solved by the parameter continuation method. To calculate the static characteristics, the differential method combined with the Newton’s Iterative refinement is used.

Текст научной работы на тему «МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И ХАРАКТЕРИСТИКИ АСИНХРОННОГО ДВИГАТЕЛЯ ПРИ ПИТАНИИ ОТ ИСТОЧНИКА ТОКА»

https://doi.org/10.21122/1029-7448-2021-64-5-421-434 УДК 621.313.3

Математическая модель и характеристики асинхронного двигателя при питании от источника тока

В. С. Маляр1*, А. В. Маляр1*

'-Национальный университет «Львовская политехника» (Львов, Украина)

© Белорусский национальный технический университет, 2021 Bekrusian National Technical University, 2021

Реферат. Разработаны методы и математические модели для исследования режимов и характеристик трехфазного асинхронного двигателя с короткозамкнутым ротором при питании обмотки статора от источника тока. Изложены особенности созданных на их основе алгоритмов расчета переходных процессов, установившихся режимов и статических характеристик. Приведено сравнение результатов расчета процессов и характеристик асинхронных двигателей при питании от источника тока и источника напряжения. Исследование установившихся и динамических режимов не может быть осуществлено с достаточной адекватностью на основе известных схем замещения и требует использования динамических параметров, являющихся элементами матрицы Якоби системы уравнений электромеханического равновесия. В математической модели уравнения состояния контуров статора и ротора составлены в неподвижной двухфазной системе координат. Переходные процессы описываются системой дифференциальных уравнений электрического равновесия преобразованных контуров двигателя и уравнением движения ротора, а установившиеся режимы -системой алгебраических уравнений. В основу разработанных алгоритмов положена математическая модель двигателя, в которой учитываются насыщение магнитопровода и явление скин-эффекта в стержнях короткозамкнутой обмотки. Для учета насыщения используются реальные характеристики намагничивания основным магнитным потоком и потоками рассеяния обмоток статора и ротора. На их основе вычисляются дифференциальные индуктивности, которые являются элементами матрицы Якоби систем уравнений, описывающих динамические режимы и статические характеристики. С целью учета скин-эффекта в обмотке короткозамкнутого ротора каждый стержень вместе с короткозамыкающими кольцами разбивается по высоте на несколько элементов. В результате в математической модели рассматриваются эквивалентные обмотки ротора с разными по значению параметрами, между которыми существуют взаимоиндуктивные связи. Решение нелинейной системы алгебраических уравнений электрического равновесия, которой описываются установившиеся режимы, выполняется методом продолжения по параметру. Для расчета статических характеристик используется дифференциальный метод в сочетании с итерационным уточнением методом Ньютона.

Ключевые слова: асинхронный двигатель, питание от источника тока, математическая модель, насыщение магнитопровода, вытеснение тока, переходный процесс, статические характеристики

Для цитирования: Маляр, В. С. Математическая модель и характеристики асинхронного двигателя при питании от источника тока / В. С. Маляр, А. В. Маляр // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2021. Т. 64, № 5. С. 421-434. https://doi.org/10. 21122/1029-7448-2021-64-5-421-434

Адрес для переписки

Маляр Василий Сафронович

Национальный университет «Львовская

ул. С. Бандеры, 12,

79013, г. Львов, Украина

Тел.: +38 032 258-21-19

vasyl.s.maliar@lpnu.ua

Address for correspondence

Malyar Vasyl S.

Lviv Polytechnic National University

12, Bandera str.,

79013, Lviv, Ukraine

Tel.: +38 032 258-21-19

vasyl.s.maliar@lpnu.ua

Mathematical Model and Characteristics of the Induction Motor with a Power Supply from a Current Source

V. S. Malyar1), A. V. Malyar1)

'-Lviv Polytechnic National University (Lviv, Ukraine)

Abstract. Methods and mathematical models for studying the modes and characteristics of the three-phase squirrel-cage induction motor with the power supplied to the stator winding from the current source have been developed. The specific features of the algorithms for calculating transients, steady-state modes and static characteristics are discussed. The results of the calculation of the processes and characteristics of induction motors with the power supply from the current source and the voltage source are compared. Steady-state and dynamic modes cannot be studied with a sufficient adequacy based on the known equivalent circuits; this requires using dynamic parameters, which are the elements of the Jacobi matrix of the system of equations of the electromechanical equilibrium. In the mathematical model, the state equations of the stator and rotor circuits are written in the fixed two-phase coordinate system. The transients are described by the system of differential equations of electrical equilibrium of the transformed circuits of the motor and the equation of the rotor motion and the steady-state modes by the system of algebraic equation. The developed algorithms are based on the mathematical model of the motor in which the magnetic path saturation and skin effect in the squirrel-cage bars are taken into consideration. The magnetic path saturation is accounted for by using the real characteristics of magnetizing by the main magnetic flux and leakage fluxes of the stator and rotor windings. Based on them, the differential inductances are calculated, which are the elements of the Jacobi matrix of the system of equations describing the dynamic modes and static characteristic. In order to take into account the skin effect in the squirrel-cage rotor, each bar along with the squirrel-cage rings is divided height-wise into several elements. As a result, the mathematical model considers the equivalent circuits of the rotor with different parameters which are connected by mutual inductance. The non-linear system of algebraic equations of electrical equilibrium describing the steady-state modes is solved by the parameter continuation method. To calculate the static characteristics, the differential method combined with the Newton's Iterative refinement is used.

Keywords: induction motor, power supply from a current source, mathematical model, magnetic path saturation, current displacement, transient, static characteristics

For citation: Malyar V. S., Malyar A. V. (2021) Mathematical Model and Characteristics of the Induction Motor with a Power Supply from a Current Source. Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. 64 (5), 421-434. https://doi.org/10.21122/1029-7448-2021-64-5-421-434 (in Russian)

Введение

В современных регулируемых асинхронных электроприводах трехфазные асинхронные двигатели (АД) питаются от преобразователей частоты. Частотно-регулируемый асинхронный электропривод может быть реализован с инвертором напряжения [1, 2] или инвертором тока [3, 4]. Большое число разновидностей схем преобразователей частоты с автономными инверторами тока и напряжения рассматриваются в литературе по преобразовательной технике. Частотные преобразователи со звеном постоянного тока и автономным инвертором управляются током, т. е. по сути двигатель получает питание от источника тока.

Исследованию процессов и характеристик АД при питании от источника напряжения посвящено много научных публикаций, однако поведе-

ние АД при питании от источника тока остается малоисследованным. Анализ электромеханических и электромагнитных процессов при частотном и частотно-токовом способах управления асинхронными электроприводами с автономным инвертором тока представлен в [5, 6]. В основу исследований положены так называемые функциональные схемы систем электропривода и имитационные модели его элементов. Расчеты выполнены с применением известного программного обеспечения, в котором используется линейная модель АД. Однако, несмотря на важность выполненных исследований для определения качественной стороны явлений, динамические режимы требуют корректного учета как насыщения магнитопрово-да АД, так и вытеснения тока в стержнях, оказывающих существенное влияние на активные и реактивные сопротивления контуров [7], а значит, и на их потокосцепления и энергетические показатели. Известны методы анализа [8], использующие линейные математические модели АД, которые дают возможность получить аналитические решения. Однако, согласно [9, 10], неучет изменения активных и индуктивных параметров АД приводит к тому, что результаты расчета существенно отличаются от реальных данных. Эффективность системы как векторного, так и скалярного управления двигателем в значительной степени зависит от точности информации о параметрах электрической машины, используемой в алгоритме управления. Причем динамические процессы описываются системами нелинейных дифференциальных уравнений (ДУ), а для реализации скалярного частотного регулирования АД достаточно воспользоваться уравнениями, описывающими статические режимы, но при этом необходимо учесть изменение параметров вследствие не только изменения частоты, но и насыщения и вытеснения токов в стержнях ротора.

Большинство методов исследования процессов асинхронных электроприводов при питании от источников напряжения и источников тока имеют в своей основе классические [7] либо адаптированные для решаемых задач схемы замещения [11-13] и формулу Клосса для определения электромагнитного момента. Однако использование схем замещения, в которых вытеснение тока в стержнях короткозамкнутого ротора и насыщение магнитопровода учитываются с помощью различных коэффициентов, неприемлемо, поскольку они зависят от протекания динамического процесса и определить их априори с достаточной точностью невозможно.

Задача количественного анализа процессов в АД при питании от источника тока на уровне современных требований может быть решена методами математического моделирования с использованием моделей АД, которые с высокой достоверностью отражают электромагнитные процессы в двигателе. Адекватные по точности и вместе с тем приемлемые с точки зрения вычислительных ресурсов и быстродействия математические модели можно разработать с помощью метода пространственных векторов [14] потокосцеплений, токов и напряжений, позволяющего создать математическую модель на основе теории цепей. При этом процессы описываются системой ДУ, а математическое моделирование динамических режимов АД требует использования динамических электромагнитных параметров [15, 16].

Цель исследований - разработать математические модели, методы и алгоритмы расчета динамических режимов и статических характеристик трехфазных АД при питании от источника тока с учетом насыщения магнитной цепи и вытеснения тока в короткозамкнутой обмотке ротора.

Математическая модель для расчета динамических режимов

В общем случае электромагнитные процессы в динамических режимах АД описываются системой ДУ. В статье задача решается в известной [16] двухфазной системе координат х, у, что позволило существенно упростить систему уравнений электрического равновесия. С целью учета вытеснения тока в стержнях обмотки ротора, а также насыщения магнито-провода ротора потоками рассеивания пазовая часть стержней и коротко-замыкающие кольца разбиваются по высоте на п слоев [17-19]. Таким образом, на роторе получим п короткозамкнутых обмоток, которые преобразуем к эквивалентным двухфазным в соответствии с общей теорией электрических машин.

В системе координат х, у полная система ДУ с учетом принятых исходных положений состоит из уравнений, составленных для эквивалентных контуров:

статора:

Жг

Ж у

уу

Жг

= «оУ уу- Г« + их;

= -«0^ - Чу, + иу;

(1а)

ротора:

Ж VI х

Жг

= у«оУ1 у - ГАх;

Ж V

1 у

Жг

Ж Уп

Жг

Ж Упу

Жг

= - у«оУ1х- гч у;

= У«оУпу - Гп^пх;

= -У«оУпх - Чгу ,

(1б)

где уух, у^, у1х, у1у, ..., упх, упу - потокосцепление преобразованных контуров; юо - частота напряжения питания, рад/с; гу, г1, ..., гп - активное сопро-

тивление контуров статора и ротора; isX, isy, ^, ^, ..., inX, iny, и*х, и*у - ток

и напряжение электрических контуров статора и ротора; 5 = (ю0 -ю)/ю0 -

скольжение ротора; ю - угловая скорость вращения ротора, рад/с.

Для расчета электромеханического переходного процесса необходимо систему (1) дополнить уравнением механического равновесия

Т=У (м -м),

(2)

где р0 - количество пар полюсов; J - момент инерции системы электропривода; Мс - момент сопротивления нагрузки на валу двигателя; Ме - электромагнитный момент АД в ортогональных координатах х, у

Me =1,5 Ро( VsJsy -Vsyisx )-

(3)

Система ДУ, состоящая из уравнений электрического (1) и механического (2) равновесия, с учетом (3) в развернутой форме имеет вид

Lsxsx Lsxsy Lsx1x Lsx1 y Lsxnx Lsxny 0 dijdt

Lsysx Lyy Lylx Lsy1 y Lsynx Lsyny 0 diy/dt

Axsx L1xsy L1x1x Axly L\xnx L\xny 0 dirlJdt

L1ysx 1 -r L ylx L1y1y Llynx L yny 0 dir1yldt

Lnxsx Lnxsy Lnx1x Lnxi1 y Lnxnx Lnxny 0 dirnxldt

Lnysx Lnysy Lny1x Lny1 y Lnynx Lnyny 0 dirny/dt

0 0 0 0 0 0 1 d ю/ dt

»oV™ - rsisx + Ux

-»oVsx - ri + uy

s»oVly - rlix

-s(3oVlx - Vly

s»oVn

-s»oVnx - Vn

a(a(VsAy -Vyix ) -Mc )

, (4)

где а = Ро^; <3 = 1,5Ро.

Дифференциальные индуктивные сопротивления, являющиеся элементами матрицы Якоби системы (4), вычисляются с учетом насыщения магнитопровода с использованием характеристик намагничивания основным магнитным потоком у и потоками рассеивания контуров статора и р°т°ра уаг [16]:

где

Оц ); Vas = Vas (is ); Var = Var (ir X

(5)

^ = V i2+i2; is = V £ + i2; v = V t + ■

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

ской модели АД рассматриваются т = к +1 двухфазных обмоток, между которыми существуют взаимоиндуктивные связи за счет как основного магнитного потока, так и потоков рассеивания. Составляющие 1гх, 7 токов преобразованной обмотки ротора (5) определяются как сумма токов к элементов стержня:

7гх ^ ; 7гу ^ •

]=1 ]=1

При математическом моделировании режимов работы АД, питаемого от источника напряжения, его составляющие их и иу считаются заданными, а токи контуров неизвестными. Интегрируя систему (4) численным методом при заданных значениях напряжения, момента нагрузки на валу АД и момента инерции системы электропривода, можно рассчитать любой динамический режим, в том числе и режим пуска, с учетом насыщения магнитной системы и вытеснения тока в обмотке ротора.

Для расчета переходных процессов при питании АД от источника тока в (1) будем полагать составляющие токов статора и заданными, а напряжений и*х и ибу неизвестными. При этом, не уменьшая общности, примем один из известных токов = 1т, т. е. равным величине изображающего вектора тока статора, а другой 1у = 0. Таким образом, система ДУ для контуров статора (1а) имеет вид:

Ж V *

Ж

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

= ЮоУ *у - Г*1т + их;

ж

(6)

а уравнения для контуров ротора остаются без изменений. Тогда система ДУ для расчета переходного электромеханического процесса при питании АД от источника тока

П1х1х П1х1у П1хпх Ахпу 0

А У1х А у1 у П1упх ^Лупу 0

Ппх1х Ппх1у ^пхпу 0

"Ппу1х 0 ^пу\у 0 Ппупх 0 ^пупу 0 0 1

ЛЛх/Л *ЮоУ1у - ГАх

Жг1 у -*юоУ1х - г71у

Ж7гпх1Ж *ЮоУпу - Гп7пх

Л„у1Л -*ЮоУпх - Гп7пу

Ж ю/ Ж «(СС) -Мс )

(7)

где Ljk - собственные и взаимные дифференциальные индуктивности соответствующих преобразованных контуров.

В результате численного интегрирования системы (7) найдем зависимости от времени токов контуров и угловой скорости ю вращения ротора при известных токах 7кх и 7ку контуров статора. При определении потокосцепле-

ний продольного у5х и поперечного у, контуров статора, входящих в уравнение электромагнитного момента АД, необходимо учитывать, что они зависят от токов всех контуров. Входящие в (7) дифференциальные взаимные индуктивности определяются на основании характеристик намагничивания (5) в соответствии с математической моделью двигателя в координатных осях х, у [16], а активные сопротивления контуров ротора - с учетом вытеснения токов в стержнях ротора [17, 18].

Уравнения (6) позволяют найти расчетные значения напряжений и5х и и5у. В частности, определив производные от потокосцеплений как функции многих переменных по формулам:

d v

dt

sx _ т 1x + + т nx .

_ ^sxlx ,, •••"'" ^sxnx ,, '

d V,

dt

sL _ т rrnx + + т .

-^sylx ^ — ^^synx

(8)

получим

d V и r

Usx -»oVsy + rJm ;

U _ -

sy

dt

- + ®oVs

На рис. 1-4 представлены результаты расчета переходных процессов при пуске двигателя 4Л160Ы6У3 (Р = 15 кВт, и = 220 В, I = 29,9 А, р0 = 3), питаемого от источника номинального напряжения и источника тока величиной = 51т.

1

О -1

м.

о

0,25

0,5

0,75

1,0

30 20 10

О

- 10 -20 -30

Ме

о

b

jilpWM'

0,4

0,8

U 16

2,0

Рис. 1. Изменение электромагнитного момента по времени при пуске двигателя с номинальной нагрузкой, питаемого от: a - источника напряжения; b - источника тока величиной 5Im

Fig. 1. Change in the electromagnetic torque in time when starting a motor with a rated load powered by: from a voltage source (a) and current source of 5Im (b)

a

t

t

1x10 S00 600 «о 200 О

/

/

/

/

у

о

0,2

0,4

о,б о,г

1,0

Рис. 2. Зависимость частоты вращения ротора от времени при пуске двигателя 4A160M6Y3 с номинальной нагрузкой при питании от: а - источника напряжения; b - источника тока

Fig. 2. The dependency of the rotor speed on the time when starting the motor 4A160M6Y3 with a rated load when powered by with rated load powered from: a - voltage source; b - current source

М*

V» Wy

30 20 10 0 10 20 30

Ilk л

1||||| w /

9 V

t -2

0

0,5

1,0

1,5

2,0

Рис. 3. Кривые изменения электромагнитного момента при пуске двигателя без нагрузки, питаемого током 5Im

Fig. 3. Curves of alterations in the electromagnetic moment when starting the motor under no load powered by a current of 5Im

Рис. 4. Зависимости составляющих по осям потокосцеплений контуров при пуске двигателя с номинальной нагрузкой

Fig. 4. Dependencies of the axis components of flux linkages of circuits at the start of the motor under rated load

Математическая модель

для расчета статических характеристик

Для расчета статических характеристик необходимо исключить из ДУ электрического состояния координату времени 7, т. е. перейти от дифференциальных уравнений к конечным. Это можно сделать путем алгебраи-зации уравнений одним из известных в математике методов (конечных разностей, конечных элементов, гармонического баланса, дифференциальных Т-преобразований, точечным методом, сплайн-методом и др.) либо путем преобразования координат, что значительно упрощает решение задачи без существенного снижения точности полученных результатов. Составим систему уравнений электрического равновесия для расчета и исследования статических характеристик АД, питаемого от источника тока.

b

a

t

В установившемся режиме (при постоянном скольжении 5) система ДУ (1) превращается в нелинейную систему конечных уравнений, которая с учетом особенностей питания АД от источника тока (/„ = 1т; 1ьу = 0) имеет вид:

«0^ + и5х = Г51т;

-ю0^5х + иу =0;

ж>0У1у - Щх = 0;

- 5®0^1х- гА у =0; (9)

пу- к?* =0;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

-5®0Упх - Чпу = 0.

Потокосцепление каждого контура при значении скольжения 5 однозначно определяется совокупностью всех токов, два из которых (/„, 15у) заданы, а остальные подлежат определению. Образовав векторы

у = (у*х, Уу, У1х, ,..., Упх, Упу ) ; х = (и5х, и5у , Ах, »1 у, » пх, » пу

5 = ( Г1т ,0,0, ..., 0, )*,

где верхний индекс (*) обозначает транспонирование, запишем систему (9) одним векторным уравнением

г( у, х 5) = 5. (10)

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

Быстрая (квадратичная) сходимость метода требует начального приближения, которое находится в зоне сходимости вычислительного процесса. Для его получения воспользуемся методом дифференцирования по параметру [14]. Для этого введем в систему (10) скалярный параметр в

?(у, х, 5) = в50. (11)

В результате дифференцирования системы (11) по в получим ДУ

Их = 5 0, (12)

а в

где Ж - матрица Якоби системы (10).

Система ДУ (10) в развернутой форме имеет вид

*у1х

пу 1х

*у1у

Л!х1у

"1 у1 у

пу 1 у

пх1у

упх

Л1 упу

1хпу

Жи*хМ Е

Жиу1 Ж Е

ЛЛх/Ж Е

ЖГ1 у/ЖЕ

Ж7гт1Ж Е

Лту!Ж Е

(13)

Интегрируя систему (13) в пределах от е = 0 до е = 1,0 (например, методом Эйлера) за несколько шагов, получим значение вектора х при заданных значениях скольжения * и тока 1т. Начальные условия (значение вектора х) нулевые, поскольку при е = 0 вектор 2 = Е&0 равен нулю. Такой алгоритм позволяет обеспечить сходимость итерационного процесса, поскольку полученное в результате интегрирования решение находится, как правило, в зоне притяжения итерационного метода Ньютона. Уточненное значение вектора х служит начальным условием для расчета статических характеристик (в частности, механической характеристики в диапазоне от * = 1,0 до * = 0). Для вычисления электромагнитного момента АД при питании от источника тока подставим в (3) значения токов: 7*х = 1т; 7!у, = 0. В результате получим

Ме = 1,5 Р0 V *у7*х. (14)

Система (10) конечных уравнений дает возможность исследовать влияние на режим работы АД изменения любой входящей в нее координаты, т. е. рассчитать статическую характеристику. Алгоритм расчета состоит из двух этапов. На первом определяем значения координат вектора х при заданном скольжении (например, * = 1,0) и вектора &0. Согласно методу Ньютона, для каждого значения скольжения приращение Дх(к) вектора х токов на к-м шаге итерации определяется по формуле

ЖМ( к) = -

-2 (х(к)),

(15)

где 2 (х(к)) - вектор невязок системы (10) при х = х(к), заданных значениях * и вектора & .

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

Результаты расчета разработанным методом механических характеристик двигателя 4Л16084У3 (Р =15 кВт, и = 230 В, I = 29,9 А, р = 2) при питании от источника напряжения и источника тока (У = 51т) приведены на

1

х

х

1

*х1 у1х Г1

0

X

*хпупх г

ях

рис. 5, причем для более полной картины расчет выполнен и для генераторного режима (начиная от 5 = -1,0). Как показывают расчеты, механические характеристики АД при питании от источника тока и источника напряжения существенно отличаются.

а Ь

М *

/

— 1

\]

-1,0

-0,5

0

0,5

5 1,0

1,50

0,75 0

-0,75

М *

-1,50:

-1

-0,5

0

0,5

1,0

Рис. 5. Механические характеристики двигателя 4A160S4Y3 при питании от: а - источника напряжения U = Um; b - источника тока I/ = 5Im

Fig. 5. Mechanical characteristics of the motor 4A160S4Y3 with a power supply from: a - voltage source U = Um; b - current source / = 5Im

s

Примеры расчета зависимости пускового электромагнитного момента и активной мощности от величины потребляемого тока двигателя 4Л16084У3 при питании от источника тока приведены на рис. 6.

а Ь

2,500 1,875

1,25 0,25 0

М*

У

3

6

9

12

20 25 10 5 0

P

3

6

9

12

Рис. 6. Зависимость электромагнитного момента (а) и активной мощности (b) от величины тока статора двигателя, питаемого от источника тока, при s = 1,0

Fig. 6. Dependence of the electromagnetic torque (а) and real power (b) on the stator current for the motor powered from a source of current at s = 1,0

ВЫВОДЫ

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

2. Основой разработанных методов является математическая модель асинхронного двигателя высокого уровня адекватности, в которой учиты-

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

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

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

5. Разработанные математические модели позволяют получить высокую точность расчетов при минимальном объеме вычислений и могут быть использованы при разработке частотно-регулируемого асинхронного электропривода.

ЛИТЕРАТУРА

1. Фираго, Б. И. Векторные системы управления электроприводами / Б. И. Фираго, Д. С. Васильев. Минск: Выш. шк., 2016. 159 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2. Соколовский, Г. Г. Электроприводы переменного тока с частотным регулированием / Г. Г. Соколовский. М.: ИЦ «Академия», 2006. 272 с.

3. Leonhard, W. Control of Electrical Drives / W. Leonhard. Berlin: Springer Verlag, 1996. 420 p. https://doi.org/10.1007/978-3-642-97646-9.

4. Мещеряков, В. Н. Система скалярного частотного управления синхронизированной асинхронной машиной с питанием обмоток статора и ротора от автономного инвертора тока / В. Н. Мещеряков, А. М. Башлыков, О. В. Мещерякова // Современные сложные системы управления X THCS'2012: материалы Междунар. науч.-техн. конф. Старый Оскол: ТНТ, 2012. С. 50-52.

5. Косенко, И. А. Сравнение частотно-токового и частотно-напряженческого способов ре-лейно-векторного управления асинхронными электроприводами с автономным инвертором тока / И. А. Косенко // Електротехшка та електроенергетика. 2008. № 1. С. 52-58.

6. Волков, А. В. Электромеханические процессы асинхронного электропривода с автономным инвертором тока при прогнозирующем релейно-векторном регулировании статор-ного напряжения / А. В. Волков, И. А. Косенко // Вюник НТУ «ХП1». Проблеми автома-тизованого електроприводу. Tеорiя i практика. 2008. № 30. С. 164-165.

7. Денник, В. Ф. Характеристики асинхронных двигателей при питании их от источника тока / В. Ф. Денник // Наук. пр. Донец. нац. техн. ун-ту. Сер.: Електротехшка i енергети-ка. 2011. № 11. С. 112-115.

8. Морозов, Д. I. Електромагштш i електромехашчш процеси в асинхронному двигуш при живленш його вщ джерела струму / Д. I. Морозов, I. С. Шевченко, Ю. П. Самчелеев // Вюник НТУ «ХП1». Проблеми автоматизованого електроприводу. Tеорiя i практика. 2013. № 36. С. 112-114.

9. Online-Identification of Electromagnetic Parameters of an Induction Motor / V. K. Ty-tiuk [et al.] // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2020. Т. 63, № 5. С. 423-440. https://doi.org/10.21122/1029-7448-2020-63-5-423-440.

10. Розрахунок статичних характеристик асинхронного частотно-регульованого електропри-вода / В. В. Каневський [та ш.] // Вюник НТУ «ХП1». Проблеми автоматизованого елект-роприводу. Теорiя i практика. 2008. № 30. С. 374-378.

11. Усольцев, А. А. Современный асинхронный электропривод оптико-механических комплексов / А. А. Усольцев. СПб.: СПбГУ ИТМО, 2011. 164 с.

12. Сафарян, В. С. Определение параметров схемы замещения асинхронной машины /

B. С. Сафарян, С. Г. Геворгян // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2015. № 6. С. 20-34.

13. Константинова, С. В. Расчет емкости для работы мини-энергокомплекса на основе асинхронного генератора в автономном режиме / С. В. Константинова, А. Ю. Капустинский, Т. М. Ярошевич // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2021. Т. 64, № 1. С. 40-50. https://doi.org/10.21122/1029-7448-2021-64-1-40-50.

14. Беспалов, В. Я. Математическая модель асинхронного двигателя в обобщенной ортогональной системе координат / В. Я. Беспалов, Ю. А. Мощинский, А. П. Петров // Электричество. 2002. № 8. С. 33-39.

15. Кулагин, Д. О. Математична модель тягового асинхронного двигуна з урахуванням на-сичення магштних кш / Д. О. Кулагин // Науковий вюник НГУ. 2014. № 6. С. 103-110.

16. Фильц, Р. В. Математические основы теории электромеханических преобразователей / Р. В. Фильц. Киев: Наукова думка, 1979. 208 с.

17. Stakhiv, P. Influence of Saturation and Skin Effect on Current Harmonic Spectrum of Asynchronous Motor Powered by Thyristor Voltage Regulator / P. Stakhiv, A. Malyar // Proceedings of the IV International Workshop Computational Problems of Electrical Engineering, June 1-3, 2005. Gdynia, 2005. P. 58-60.

18. Маляр, В. С. Установившиеся режимы и статические характеристики трехфазного асинхронного двигателя при питании от однофазной сети / В. С. Маляр, А. В. Маляр // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2016. Т. 59, № 6.

C. 536-548. https://www.doi.org/10.21122/1029-7448-2016-59-6-536-548.

19. Mathematical Modeling of Processes in Asynchronous Motors with Capacitors Connected in Series / V. Malyar [et al.] // Proceedings of 16th International Conference on Computational Problems of Electrical Engineering (CPEE), Sept. 2-5, 2015. Lviv, 2015. Р. 107-109. https://doi.org/10.1109/cpee.2015.7333350.

Поступила 20.10.2020 Подписана в печать 17.03.2021 Опубликована онлайн 30.09.2021 REFERENCES

1. Figaro B. I., Vasil'ev D. S. (2016) Vector Control Systems of Electric Drives. М^к, Vysheishaya shkola Publ. 159 (in Russian).

2. Sokolovskii G. G. (2006) AC Electric Drives with Frequency Control. Moscow, Akademiya Publ. 272 (in Russian).

3. Leonhard W. (1996) Control of Electrical Drives. Berlin, Springer Veriag Publ. 420. https://doi.org/10.1007/978-3-642-97646-9.

4. Meshcheryakov V. N., Bashlykov А. М., Meshcheryakova О. V. (2012) Scalar Frequency Control System for a Synchronized Asynchronous Machine Powered by the Stator and Rotor Windings from an Autonomous Current Inverter. Sovremennye Slozhnye Sistemy Upravleniya X THCS'2012: materialy Mezhdunar. Nauch.-Tekhn. Konf. = Modern Complex Control Systems X TNCS'2012: Proceedings of the International Scientific and Technical Conference. Stary Оskol, TNT Publ. 50-52 (in Russian).

5. Kosenko I. А. (2008) Comparison of Frequency-Current and Frequency-Stressed Methods of Relay-Vector Control Asynchronous Electric Drives with Autonomous Inverter of Current. Elektroteknika ta Elektroenegetika = Electrical Engineering and Power Engineering, (1), 52-58 (in Russian).

6. Volkov А. V., Kosenko I. А. (2008) Electromechanical Processes of an Asynchronous Electric Drive with an Autonomous Current Inverter with Predictive Relay-Vector Regulation of Stator Voltage. Visnik Natsional'nogo Tekhnichnogo Universitetu "KhPI". Problemi Avtomati-zovanogo Elektroprivodu. Teoriya i Praktika [Bulletin of the National Technical University "Kharkiv Polytechnic Institute". Problems of an Automated Electric Drive. Theory and Practice"], (30), 164-165 (in Russian).

7. Dennik V. F. (2011) Characteristics of Induction Motors Powered from Current Power Supply. Naukovi Pratsi Donets 'kogo Natsional 'nogo Tekhnichnogo Universitetu. Seriya: Elektrotekh-nika i Energetika [Scientific Works of Donetsk National Technical University. Series: Electrical Engineering and Power Engineering], (11), 112-115 (in Russian).

8. Morozov D. I., Shevchenko I. S., Samcheleev Yu. P. (2013) Electromagnetic and Electromechanical Processes in an Induction Motor Powered by a Current Source. Visnik Natsional'nogo Tekhnichnogo Universitetu "KhPI". Problemi Avtomatizovanogo Elektroprivodu. Teoriya i Praktika [Bulletin of the National Technical University "KhPI". Problems of Automated Electric Drive. Theory and Practice], (36), 112-114 (in Ukrainian).

9. Tytiuk V. K., Baranovskaya M. L., Chornyi O. P., Burdilnaya E. V., Kuznetsov V. V., Bo-gatyriov K. N. (2020) Online-Identification of Electromagnetic Parameters of an Induction Motor. Energetika. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob 'edinenii SNG = Energetika. Proceedings of CIS Higher Education Institutions and Power Engineering Associations, 63 (5), 423-440. https://doi.org/10.21122/1029-7448-2020-63-5-423-440.

10. Kanev'skyi V. V., Pustoviy O. V., Sinchuk I. O., Pas'ko O. V. (2008) Calculation of Static Characteristics of Asynchronous Frequency-Controlled Electric Drive. Visnik Natsional'nogo Tekhnichnogo Universitetu "KhPI". Problemi Avtomatizovanogo Elektroprivodu. Teoriya i Praktika [Bulletin of the National Technical University "Kharkiv Polytechnic Institute". Problems of an Automated Electric Drive. Theory and Practice], (30), 374-378 (in Ukrainian).

11. Usoltsev A. A. (2011) Modern Asynchronous Electric Drive of Optical-and-Mechanical Complexes. Saint Petersburg: ITMO University Publ. 164 (in Russian).

12. Safaryan V. S., Gevorgyan S. G. (2015) Ascertainment of the Equivalent Circuit Parameters of the Asynchronous Machine. Energetika. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob'edinenii SNG = Energetika. Proceedings of CIS Higher Education Institutions and Power Engineering Associations, (6), 20-34 (in Russian).

13. Konstantinova S. V., Kapustinskii A. Yu., Yaroshevich T. M. (2021) Calculation of the Capacity for the Operation of a Mini-Energy Complex Based on an Independently Operating Asynchronous Generator. Energetika. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob 'edinenii SNG = Energetika. Proceedings of CIS Higher Education Institutions and Power Engineering Associations, 64 (1), 40-50 https://doi.org/10.21122/1029-7448-2021-64-1-40-50 (in Russian).

14. Bespalov V. Ya., Moshchinskiy Yu. A., Petrov A. P. (2002) Mathematical Model of Asynchronous Motor in Generalized Orthogonal Coordinates System. Elektrichestvo, (8), 33-39 (in Russian).

15. Kulagin D. O. (2014) The Mathematical Model of Asynchronous Traction Motor Taking into Account the Saturation of Magnetic Circuits. Naukovyi Visnyk Natsionalnoho Hirnychoho Universytetu [Scientific Herald of National Mining University], (6), 103-110 (in Ukrainian).

16. Fil'ts R. V. (1979) Mathematical Foundations of the Theory of Electromechanical Transducers. Kiev, Naukova Dumka Publ. 208 (in Russian).

17. Stakhiv P., Malyar A. (2005) Influence of Saturation and Skin Effect on Current Harmonic Spectrum of Asynchronous Motor Powered by Thyristor Voltage Regulator. Proceedings of the IV International Workshop Computational Problems of Electrical Engineering, June 1-3, 2005. Gdynia. 58-60.

18. Malyar V. S., Malyar A. V. (2016) Established Modes and Static Characteristics of Three-Phase Asynchronous Motor Powered with Single Phase Network. Energetika. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob'edinenii SNG = Energetika. Proceedings of CIS Higher Education Institutions and Power Engineering Associations, 59 (6), 536-548. https://www.doi.org/10.21122/1029-7448-2016-59-6-536-548 (in Russian).

19. Malyar V., Hamola O., Maday V., Vasylchyshyn I. (2015) Mathematical Modeling of Processes in Asynchronous Motors with Capacitors Connected in Series. Proceedings of 16th International Conference on Computational Problems of Electrical Engineering (CPEE), September 2-5, 2015. Lviv. 107-109. https://doi.org/10.1109/cpee.2015.7333350.

Received: 20 October 2020 Accepted: 17 March 2021 Published online: 30 September 2021

i Надоели баннеры? Вы всегда можете отключить рекламу.