Методика расчета параметров асинхронного двигателя для моделирования режимов его работы в среде Matlab/Simulink
Пантель О. В.
Пантель Олег Валентинович /Pantel Oleg Valentinovich - заведующий лабораторией электромеханики,
кафедра электрических сетей и систем,
Институт ядерной энергии и промышленности,
Севастопольский государственный университет, г. Севастополь
Аннотация: в статье рассматривается методика расчета параметров асинхронного двигателя для моделирования режимов его работы в среде Matlab/Simulink. Приводится порядок и пример расчета параметров.
Ключевые слова: расчет, параметры, асинхронный двигатель, схема, короткое замыкание,
Matlab/Simulink.
Параметрами асинхронного двигателя (АД) являются коэффициенты перед независимыми переменными: напряжением, частотой и моментом сопротивления в уравнениях электромеханического преобразования энергии. Это могут быть индуктивности, активные сопротивления, момент инерции. Точность результатов вычислительного эксперимента на модели определяется точностью задания начальных параметров для конкретного режима работы АД.
При пуске двигателя из-за насыщения его магнитной цепи индуктивные параметры существенно отличаются от их значений при установившемся режиме. Так, например, индуктивности и взаимные индуктивности при пуске на 30-40 % меньше их значений в номинальном режиме.
За счет вытеснения тока в стержнях беличьей клетки активное сопротивление обмотки ротора во время работы АД не остается постоянным. Его значение при пуске двигателя значительно превышает таковое в номинальном режиме работы.
Таким образом, можно утверждать, что применение параметров, соответствующих номинальному режиму при пуске двигателя, приводит к значительным погрешностям в определении ударных значений токов и моментов. В соответствии с рекомендациями [6], при пуске АД для обеспечения более полного совпадения расчетных динамических характеристик двигателя с экспериментальными, необходимо использовать в качестве начальных параметры режима короткого замыкания машины.
Задача определения параметров АД для режима его пуска значительно упрощается при использовании справочной информации [3]. В таблицы [3] включены расчетные значения параметров упрощенной Г-образной схемы замещения АД (рис. 1): главного индуктивного сопротивленияХм; активного сопротивления R" и индуктивного сопротивления рассеяния X" обмотки статора; приведенных к обмотке статора
активного сопротивления R" и индуктивного сопротивления рассеяния X" обмотки ротора.
R2
s
Соотношения параметров схем замещения двигателей единых серий таковы, что переход от Т-образной к упрощенной Г-образной схеме замещения не приводит к заметным погрешностям. Поэтому для
практических расчетов вполне допустимо считать, что R = R, X[ = Х, R = R2, X" = Х'2 .
Для уменьшения погрешности в определении тока намагничивания I, в намагничивающий контур Г-
образной схемы замещения (рис. 1) введено сопротивление Z1 = R + jXx.
Кроме указанных параметров для АД с короткозамкнутым ротором в таблицах [3] даны параметры схемы замещения при коротком замыкании приведенного к обмотке статора активного сопротивления
обмотки ротора с учетом вытеснения тока в стержнях беличьей клетки R2n; активного Якп и индуктивного Хкп сопротивлений короткого замыкания машины.
Все значения параметров схемы замещения (рис. Г.1) приведены в относительных единицах X = X[Ом]- V3 IH /UH, R = R[Om\- V3 Ih /UH, где X[Ом] и R[Om] - соответствующие сопротивления в Ом.
Рассмотрим методику расчета параметров двигателей с короткозамкнутым ротором, для которых известны значения сопротивлений схемы замещения при коротком замыкании. Например, это относится к АД единой серии типа 4А [3], имеющих следующие каталожные данные:
1)
2)
3)
4)
5)
6)
7)
8)
9)
10) 11) 12)
13)
14)
15)
16)
R” [о е.]; 17)
номинальная мощность Рн [кВт]; номинальное линейное напряжение UH [В]; частота сети fH [Гц];
частота вращения магнитного поля (синхронная частота вращения) по [об/мин]; номинальное скольжение sH [%]; критическое скольжение sK [%];
коэффициент мощности в номинальном режиме cos фн;
коэффициент полезного действия в номинальном режиме цн;
кратность пускового тока к;
кратность пускового момента кп;
кратность максимального момента км;
момент инерции J [кг м2];
главное индуктивное сопротивление [о.е.];
активное сопротивление обмотки статора в номинальном режиме R [о.е.];
индуктивное сопротивление рассеяния обмотки статора в номинальном режиме X'x [о.е.]; приведенное к обмотке статора активное сопротивление обмотки ротора в номинальном режиме
приведенное к обмотке статора индуктивное сопротивление рассеяния обмотки ротора в
номинальном режиме X [о.е.];
18) приведенное к обмотке статора активное сопротивление обмотки ротора при коротком замыкании АД R'2 п [о.е.];
19) индуктивное сопротивление короткого замыканияХкп [о.е.].
Порядок расчета параметров
1. Ориентируясь на рекомендуемые значения механических потерь в двигателе ХРмех=(0,01...0,05)Рн, выбираем их величину для конкретной машины.
2. Определяем номинальное значение тока статора АД
I,, =
P -103
V3u„ -
cos срн -Лн
[ А].
(Г1)
3. Определяем активное сопротивление обмотки статора АД
R =
S -1.
[Ом]. (Г2)
4. Определяем приведенное активное сопротивление обмотки ротора АД. Для установившегося номинального режима работы:
R= RU
V3 -1
[Ом]. (Г3)
Для режима прямого пуска:
Rff U
K= -^0^ [Ом]. (Г4)
5. Определяем индуктивности рассеяния обмоток статора и ротора АД. Для установившегося номинального режима работы:
т _ X U
2ж-/н V3 • I„
г _ _XL _U^_
L" _ 2 ж • /я \Г3 • Iя
Для режима прямого пуска:
у
Т _ Tf _ X кп Lis _ Ltr _ '
U„
[ Гя ]
[ Гя ].
[ Гя ].
(Г5)
(Г6)
(Г7)
4Ж- f„ S • 1я
6. Определяем индуктивность цепи намагничивания АД в номинальном режиме
U
х и
L _ и
2ж-f„ 43 • 1я
[ Гя ].
(Г8)
Примечание. Значение индуктивности цепи намагничивания Lm для пускового режима АД необходимо уменьшить на 30.. .40 %.
7. Определяем коэффициент вязкого трения
АР.„
F _■
[ Нм • с]. (Г9)
(2ж • n / 60)2 Здесь номинальная частота вращения АД
пЯ _ no • (1 - s„) [об/мия]. (Г10)
Примечание. Если не заданы номинальное s„ и критическое sK скольжения, они рассчитываются по выражениям
п - п..
п
_ s„ • (км + Jkf-i) -100%
и
s„ _
При отсутствии в каталожных данных значений сопротивлений при коротком замыкании двигателя для определения его параметров в пусковом режиме можно воспользоваться следующим порядком расчета.
1. Задаемся конструктивным коэффициентом с} = (1,01. 1,05), который имеет меньшее значение для машин большей мощности.
2. Ориентируясь на рекомендуемые значения механических потерь в двигателе ЛРмех=(0,01...0,05)Р„, выбираем их величину для конкретной машины.
3. Определяем номинальное значение тока статора АД по выражению (Г.1)
4. Определяем активное сопротивление обмотки статора АД
з (и„/У3)2 • (1 - s,)
2 q • (1 + cj SK ) • км ■ (РЯ +АРмех )
[Ом]
(Г11)
5. Определяем приведенное активное сопротивление обмотки ротора АД
r I • (Ря +АРмех) • К, [Ом] (Г 12)
r 3 (1 - s ) • k2 • I2
v я / t я
6. Определяем индуктивность цепи намагничивания
L
m
1 uh/43__________
2ж/я 1 „U1 -(cosp„)2 - cos^, •sЯ /SK]
[ Гя]
(Г13)
7. Определяем индуктивности рассеяния обмоток статора и ротора АД
L,
L _
1
•V[(U„ Л/3/(kt I„ )]2 - (R + R)2 [r„]
4ж‘ /
J н
8. Определяем индуктивности обмоток статора и ротора АД
Ls « L'r_ Lm + Ls [ Гя]
(Г 15)
(Г 14)
9. Определяем коэффициент вязкого трения по выражению (Г.9)
10. Определяем расчетное значение конструктивного коэффициента
c = 1 + — 1 L
(Г 16)
11. Сравниваем расчетное значение конструктивного коэффициента с ранее принятым.
При необходимости расчет следует повторить, добиваясь минимальной разницы между ними.
Примечание. Для режима пуска двигателя рассчитанное значение индуктивности цепи намагничивания Lm необходимо уменьшить на 30.. .40 %.
Полученные расчетным путем параметры схемы замещения асинхронного двигателя устанавливают в окне задания параметров блока Block Parameters: Asynchronous Machine SI Units.
Пример расчета параметров асинхронного двигателя с короткозамкнутым ротором 4А112МВ6УЗ для номинального режима работы.
Из таблицы Б1 [3]: Рн=4,0 кВт, UH=380 В, fH=50 Гц, n0=1000 об/мин., sH=5,1 %, cos фн=0,81, qH= 0,82, J=0,021 кгм2, Х=2,0, R1=0,077, X=0,073, R20,062, Х'2=0,11.
Порядок расчета.
1. Определяем механические потери в двигателе
АРмех=(0,01...0,05)Рн=0,0354103=140 (Вт).
2. Определяем номинальное значение тока статора АД
I = PH ■ 103
н SUн ■ COSP Лн
4 ■ 103
43 ■ 380 ■ 0,81 ■ 0,82
3. Определяем активное сопротивление обмотки статора АД
R'r U,, 0,077 ■ 380
9,15( А )
Rs =^^ =~’~г' ~~ (Ом)
s 43 ■ 1н 43 ■ 9,15
4. Определяем приведенное активное сопротивление обмотки ротора АД
0,062 ■ 380
R' = R" 'U>
43 ■ i 43 ■ 9,15
= 1,49 (Ом)
5. Определяем индуктивности рассеяния обмоток статора и ротора АД
U 0,073 380
Lls =
l; =
X
2ж ■ fH 43 ■ IH 2 ■ 3,14 ■ 50 43 ■ 9,15 U 0,11 380
Х2
2п ■ f 43 ■ Т 2 ■ 3,14 ■ 50 43 ■ 9,15
6. Определяем индуктивность цепи намагничивания
L =
X,,
U,,
2
380
2л ■ fH 43 ■ IH 2 ■ 3,14 ■ 50 43 ■ 9,15
= 0,0054 ( Гн)
= 0,0084 ( Гн)
= 0,153 ( Гн)
7. Определяем коэффициент вязкого трения
ЛР_
F = ■
140
= 0,014 (Нм ■ с)
(2л ■ пн / 60)2 (2 ■ 3,14 ■ 950/60)2
Здесь номинальная частота вращения АД
П = n0 ■ (1 -sH) = 1000 ■ (1 - 0,051) = 950 (об/мин)
8. Определяем суммарный момент инерции двигателя и исполнительного механизма. Принимаем 3=23д=20,021=0,042 (кгм2).
Литература
1. ЧиликинМ. Г. Общий курс электропривода. М.: Энергоиздат, 1981. 576 с.
2. Герман-Галкин З. Р. Электрические машины: Лабораторные работы на ПК. СПб.: Корона принт, 2003. 256 с.
3. Кравчик А. Э., Шлаф М. М., Афонин В. И., Соболенская Е. А. Асинхронные двигатели серии 4А. М.: Энергоиздат, 1982. 504 с.
4. Елисеева В. А., Шинянский А. В. Справочник по автоматизированному электроприводу. М.: Энергия, 1983. 616 с.
5. Богословский А. П., Певзнер Е. М., Фрейдзон И. Р., Яуре А. Г. Судовые электроприводы. Л.: Судостроение, 1983. 352 с.
6. Копылов И. П. Электрические машины: Учебник для вузов. М.: Энергоатомиздат, 1986. 360 с.