CC BY
67
ТЕХНИЧЕСКИЕ НАУКИ
УДК 621.313.333.2:621.3.016.1
ПОВЫШЕНИЕ ТОЧНОСТИ РАСЧЕТА ВРАЩАЮЩЕГО МОМЕНТА АСИНХРОННОГО ДВИГАТЕЛЯ
С УЧЕТОМ ПОВЕРХНОСТНОГО ЭФФЕКТА В СТЕРЖНЯХ РОТОРА М.Р. ВИНОКУРОВ
(Донской государственный технический университет),
А.А. МОИСЕЕНКО
(Омский государственный технический университет),
Н.Ю. МАСЛОВЦЕВА
(Рудненский индустриальный институт)
Проведена количественная оценка точности определения механического момента асинхронного двигателя при использовании Т-образной схемы замещения с учетом поверхностного эффекта в стержнях ротора по сравнению с Г-образной схемой замещения.
Ключевые слова: механический момент, схема замещения, асинхронный двигатель, скольжение, вихревые токи, ферромагнитная среда.
Введение. При расчете вращающего момента асинхронного двигателя (АД), как правило, используется Г-образная схема замещения. Кроме того, отмечается [1], что переход от Т-образной схемы замещения к упрощенной Г-образной приводит к погрешности 3-5% в определении силы тока ротора. Если учесть, что момент АД пропорционален квадрату тока ротора, то погрешность определения момента АД оказывается весьма существенной. Известно также, что на величину вращающего момента оказывают влияние вихревые токи в стержнях ротора и нелинейность магнитных характеристик статора [2].
Постановка задачи и метод исследования. Рассмотрим результаты расчета вращающего момента для Г-образной схемы замещения. Как известно [1], если пренебречь действием тока намагничивания в поперечной ветви (т.е. принять его равным нулю) и не учитывать действие вихревых токов, то схема замещения Ад имеет вид, показанный на рис.1, где R1 - активное сопротивление обмотки статора; 7Х1 - индуктивное сопротивление рассеяния обмотки статора; R'2 - приведенное активное сопротивление обмотки ротора; 7Х'2 - приведенное индуктивное сопротивление рассеяния обмотки ротора; Я'2 -—- - приведенное сопротивление, учитывающее активную мощность, эквивалентную механическому моменту АД; - - скольжение.
Рис.1. Г-образная схема замещения асинхронного двигателя
Вращающий момент в этом случае определяется известным выражением [1]:
1
Р
ю
ті 22Г R2
ю
mR2
и,2
Ю5
Л + Л
(1)
+ (Хі + Х 2)2
где Рмех - механическая мощность АД; ю - угловая скорость АД; т - поправочный коэффициент; 12Г - приведенный ток в нагрузке Г-образной схемы замещения; и1 - действующее значение напряжения фазы статорной обмотки.
Расчет вращающего момента по приведенному соотношению (1) дает завышенные результаты по двум причинам. Одна заключается в том, что в этой схеме замещения не учитывается влияние на ток ротора тока намагничивания в поперечной ветви, который создает падение напряжения на сопротивлениях обмотки статора, что приводит к уменьшению мощности на приведенном сопротивлении ротора Л'2. Другая обусловлена тем, что не учитывается влияние вихревых токов в стержнях ротора.
Если рассмотреть вращающий момент, получаемый с использованием Т-образной схемы замещения, т.е. учесть ветвь намагничивания (рис. 2), то для расчета приведенного комплексного
значения тока ротора 12 можно получить выражение вида:
Л Х о
где
12= и ,
2 1 М +
Л Л
М = Л1Л0 - X Х0 + Л - Х1Х2 + Л,— - Х0 Х2;
5 5
л Л
N = л х1 + л х0 +^ х1 + л Х2 + х0 + л Х2.
5 5
(2)
Выражение для 1'2 в источнике [1] приводится в несколько иной форме. На наш взгляд, целесообразно использовать выражение (2) для удобства сравнения с выражением для момента, определяемого по Г-образной схеме замещения.
„1-5
Рис.2. Т-образная схема замещения асинхронного двигателя
С учетом зависимости (2) вращающий момент АД для Т-образной схемы замещения АД можно рассчитать следующим образом:
МТСЗ =
Р..
ті22 Л / 5 _ и2 т( +Х02)
ю
ю ю (М2 + N2)
(Л2 + Хо2)
т
Л2/ 5 = п\ — и2/ 5,
ю
(3)
П =
М2 + N ’
где П - вещественный коэффициент для данного типа АД, определяемый параметрами эквивалентной схемы замещения и скольжения.
2
Выражение (3) позволяет выполнять расчет вращающего момента в реальном времени с помощью микроконтроллера, при этом требуется измерять всего 2 рабочих параметра - напряжение питания и частоту вращения двигателя.
Проведем сравнение расчетных данных вращающего момента, полученных по формулам (1) и (3) (т.е. для Г- и Т-образной схем замещения АД) для асинхронного двигателя типа
4А80А2У3 со следующими паспортными данными:
- номинальное напряжение и1, В.......................... 3x220/380
- номинальная мощность Р2ном, кВт.......................1,5
- номинальное скольжение 5ном...........................0,042
- критическое скольжение 5кр............................0,37
- номинальная частота напряжения питания /1ном, Гц......50
- номинальная частота вращения ротора п2ном, мин-1 ..... 2835
Параметры схемы замещения, Ом: Х=3,73, R1=7,82, Х2=4,21, R2=2,91, Х0=133, R0=11,8. Эквивалентные параметры АД приведены в справочнике [3].
Расчетные значения вращающего момента получены для значений скольжения в диапазоне 0,01<5<5кр и приведены в табл.1.
Таблица 1
Вращающий момент, Н-м 5 5н 5кр
0,01 0,02 0,03 0,042 0,1 0,2 0,37
Мгсз 1,49 2,8 3,94 5,33 8,49 9,55 7,41
Мтсз 1,4 2,63 3,71 5,04 8,1 9,19 7,15
Относительные значения разностей вращающих моментов (относительная погрешность),
вычисленных по Г- и Т-образной схемам замещения, определяются следующим образом:
М - М 5 = ГСЗ ТСЗ -100%.
Значения 5 в зависимости от скольжения 5 приведены в табл.2. График зависимости 5=/5) приведен на рис.3.
________________________________________________________________________________________Таблица 2
Погрешность 5 5н 5кр
0,01 0,02 0,03 0,042 0,1 0,2 0,37
5, % 6,27 6,01 5,77 5,7 4,57 3,8 3,46
0 0.05 0,10 0.15 0.20 0;25 0;30 С>
Рис.3. Зависимость 5 от скольжения 5
Как видно из приведенных расчетных данных, использование Т-образной схемы замещения позволяет существенно увеличить точность расчета вращающего момента.
Однако расчет вращающего момента по выражению (3) основан на том, что значения самих параметров схемы замещения считаются постоянными при изменении скольжения.
В установившемся режиме работы в широком диапазоне изменения скольжения рекомендуется учитывать роторные вихревые токи, оказывающие влияние на параметры Т-образной схемы замещения [1].
Для оценки влияния вихревых токов при изменении частоты тока в стержнях ротора (изменение скольжения я) на параметры схемы замещения ротора рассмотрим зависимости полного, активного и реактивного сопротивлений стержня ротора от скольжения для АД типа 4А80А2У3. Данные зависимости получим исходя из численного расчета электромагнитного поля, создаваемого переменным током в стержне ротора. В исследовании [4] проведен численный расчет электромагнитного поля для различных форм пазов ротора, при этом не учитывается влияние ферромагнитной среды статора, что снижает точность расчета параметров стержней ротора.
Выполним расчет электромагнитного поля с учетом влияния ферромагнитной среды статора, основываясь на следующих условиях:
- решение задачи электромагнитного поля ротора осуществим в предположении, что ферромагнитные среды ротора и статора линейны;
- электропроводностью ферромагнетика ротора и статора пренебрегаем;
- материал стержней ротора - алюминий, у=37-106 См/м;
- геометрия и размеры пазов ротора и статора для данного типа двигателя (рис.4);
- действующее значение тока в стержне ротора 10 А;
- для расчета электромагнитного поля и параметров стержня ротора используем программу ЕІсиґ,
- расчет выполняем в диапазоне скольжений ях<я<3яном, скольжение моделируем через задание соответствующей частоты переменного тока стержня ротора.
С учетом вращения ротора относительно статора, взаимное расположение пазов ротора и статора в различные
моменты времени изменяется, при этом можно выделить два крайних случая:
- статор представлен сплошным металлическим телом (ц>>цо, у~0);
- паз статора расположен симметрично над пазом
ротора.
Для учета ферромагнитной среды статора необходимо выполнить расчет электромагнитного поля и параметров стержня ротора для двух указанных
крайних случаев, а затем, используя расчетные данные, вычислить средние значения этих пара-
метров.
На рис.5 показаны расчетные силовые линии индукции магнитного поля, создаваемого переменным током в стержне ротора для двух крайних случаев.
В табл.3 приведены результаты расчета параметров алюминиевого стержня с учетом действия поверхностного эффекта на различных частотах тока в стержне ротора для случая, когда статор представлен сплошным металлическим телом (рис.5, а).
В табл.4 приведены результаты расчета параметров алюминиевого стержня с учетом действия поверхностного эффекта на различных частотах тока в стержне ротора для случая, когда паз статора расположен симметрично над пазом ротора (рис.5, б).
[ 8'5 ]
\ / ^
1 1 ^ 1 6,8
1 і
сО-1 ^ І ^ і ґ^\
Ісз- | ^ 1
1 6,1 1
4* 7 »1\
\ 1
33\-і ,
Рис.4. Геометрия и размеры пазов
ротора и статора
асинхронного двигателя 4А80А2У3
а) б)
Рис.5. Расчетная картина магнитного поля, создаваемого переменным током в стержне ротора двигателя 4А80А2У3: а - статор представлен сплошным металлическим телом; б- паз статора расположен над пазом ротора
Таблица 3
5 /, Гц 2, Ом R, Ом X, Ом L, Гн
0,12 6 6,766Е-04 4,295Е-04 5,228Е-04 1,387Е-05
0,1 5 6,116Е-04 4,293Е-04 4,356Е-04 1,387Е-05
0,08 4 5,528Е-04 4,291Е-04 3,485Е-04 1,387Е-05
0,06 3 5,024Е-04 4,290Е-04 2,614Е-04 1,387Е-05
0,042 2,1 4,663Е-04 4,289Е-04 1,830Е-04 1,387Е-05
0,03 1,5 4,483Е-04 4,289Е-04 1,307Е-04 1,387Е-05
0,02 1 4,376Е-04 4,289Е-04 8,713Е-05 1,387Е-05
0,0032 0,16 4,291Е-04 4,288Е-04 1,394Е-05 1,387Е-05
Таблица 4
5 /, гц 2, Ом R, Ом X, Ом L, Гн
0,12 6 6,542Е-04 4,295Е-04 4,935Е-04 1,309Е-05
0,1 5 5,945Е-04 4,293Е-04 4,113Е-04 1,309Е-05
0,08 4 5,407Е-04 4,291Е-04 3,290Е-04 1,309Е-05
0,06 3 4,949Е-04 4,290Е-04 2,468Е-04 1,309Е-05
0,042 2,1 4,624Е-04 4,289Е-04 1,727Е-04 1,309Е-05
0,03 1,5 4,463Е-04 4,289Е-04 1,234Е-04 1,309Е-05
0,02 1 4,367Е-04 4,289Е-04 8,225Е-05 1,309Е-05
0,0032 0,16 4,290Е-04 4,288Е-04 1,316Е-05 1,309Е-05
На основе расчетных данных, приведенных в табл.3 и 4 для двух крайних случаев расположения пазов ротора относительно пазов статора, рассчитаны средние значения полного, активного и реактивного сопротивлений стержня на единицу длины для различных значений скольжения. Результаты расчета приведены в табл.5.
_____________________________________________________________________________________Таблица 5
5 / Гц 7, Ом Я, Ом X, Ом L, Гн
0,12 6 6,654Е-04 4,295Е-04 5,08Е-04 1,350Е-05
0,1 5 6,031Е-04 4,293Е-04 4,23Е-04 1,350Е-05
0,08 4 5,468Е-04 4,291Е-04 3,39Е-04 1,350Е-05
0,06 3 4,987Е-04 4,290Е-04 2,54Е-04 1,350Е-05
0,042 2,1 4,644Е-04 4,289Е-04 1,78Е-04 1,350Е-05
0,03 1,5 4,473Е-04 4,289Е-04 1,27Е-04 1,350Е-05
0,02 1 4,372Е-04 4,289Е-04 8,47Е-05 1,350Е-05
0,0032 0,16 4,291Е-04 4,288Е-04 1,36Е-05 1,350Е-05
Примем значения сопротивлений стержня ротора при номинальном режиме (5ном=0,042 для выбранного типа двигателя) за 100% и рассчитаем относительные изменения этих значений сопротивлений в зависимости от скольжения по выражениям:
7 - 7
Еу = -І----Н -100 %,
Я - Я
=-
Яы
X. - X,
ЕХ =-
хн
Результаты расчета сведены в табл.6.
Таблица 6
5 Е-, % Ея, % Ех, %
0,12 43,28 0,14 185,39
0,1 29,87 0,09 137,64
0,08 17,74 0,05 90,45
0,06 7,39 0,02 42,70
0,042 0,00 0,00 0,00
0,03 -3,68 0,00 -28,65
0,02 -5,86 0,00 -52,42
0,0032 -7,60 -0,02 -92,36
Из табл.6 следует, что изменение активного сопротивления ротора под действием вихревых токов незначительно и при 5=35ном составляет 0,14%. Изменение реактивного сопротивления стержней ротора при увеличении скольжения существенно и объясняется увеличением частоты тока, протекающего в стержнях.
Для учета изменения реактивных сопротивлений стержней ротора при расчете вращающего момента определим значения приведенного сопротивления Х'2 для эквивалентной схемы замещения при выбранном значении скольжения по формуле:
8 Х
Хп _ VI . \тг л,
. — У1 + ---- ,
2, 2 2 шо
8 Х
где X 2у^ - приращение приведенного реактивного сопротивления ротора в зависимости от скольжения.
Уточненный расчет вращающего момента выполнен по выражению (3) с учетом изменения реактивных сопротивлений стержней ротора.
Результаты расчета значений вращающего момента сведены в табл.7:
- для Г-образной схемы замещения (МГСЗ);
- Т-образной схемы замещения (МТСЗ);
- Т-образной схемы замещения с учетом изменения Х2 (МТСЗ,ПЭ).
Таблица 7
5 Мгсз, Н-м Мтсз, Н-м Мгсз.п.э, Н-м
0,12 10,27 9,82 8,34
0,1 9,43 9,00 8,24
0,08 8,34 7,94 7,64
0,06 6,93 6,57 6,51
0,042 5,33 5,04 5,04
0,03 4,06 3,82 3,83
0,02 2,85 2,68 2,69
0,0032 0,50 0,47 0,47
По данным табл.7 построены механические характеристики АД для перечисленных выше случаев (рис.6).
ГІ2, мин"1
3000
2950
2900
2850
2800
2750
2700
2650
2600 -О
Рис.6. Расчетные механические характеристики двигателя 4А80А2У3
Полученные результаты могут быть использованы при реализации микропроцессорного устройства бесконтактного контроля вращающего момента АД [5] для уточнения косвенного метода измерения момента с использованием Т-образной схемы замещения АД в широком диапазоне изменения скольжения.
Выводы. 1. Использование Т-образной схемы замещения позволяет повысить точность расчета вращающего момента асинхронного двигателя для номинального режима на 5% по сравнению с расчетом по Г-образной схеме замещения.
2. При понижении оборотов асинхронного двигателя ниже номинальных (5>5ном) погрешность расчета вращающего момента по Г-образной схеме замещения понижается незначительно (см. рис.3).
3. Влияние изменения реактивного сопротивления ротора Х'2 существенно возрастает при
?>2?
ном
4. При 5<35ном влиянием вихревых токов на изменение активных сопротивлений стержней ротора можно пренебречь.
Библиографический список
1. Копылов И.П. Электрические машины: учеб. для вузов / И.П. Копылов. - М.: Высш. шк.: Логос, 2000. - 607 с.
2. Иванов-Смоленский А.В. Электрические машины: учеб. для вузов. В 2 т. Т.1. - 2-е изд., перераб. и доп. / А.В. Иванов-Смоленский. - М.: Изд-во МЭИ, 2004. - 656 с.
3. Асинхронные двигатели серии 4А: справочник / А.Э. Кравчик [и др.]. - М.: Энергоатом-издат, 1982. - 504 с.
4. Попов А.П. Исследование поверхностного эффекта в проводниках произвольного сечения с помощью программы Е1си: / А.П. Попов, А.О. Чугулёв // Омский науч. вестн. - 2009. -№3(83). - С.183-187.
5. Попов А.П. Моделирование микропроцессорного устройства контроля вращающего момента асинхронного двигателя с использованием его эквивалентной схемы замещения и импульсного датчика оборотов / А.П. Попов, М.Р. Винокуров, А.А. Моисеенко // Вестн. Донск. гос. техн. ун-та. - 2011. - Т.11, №1. - С.77-82.
Материал поступил в редакцию 10.05.2011.
References
1. Kopy'lov I.P. E'lektricheskie mashiny': ucheb. dlya vuzov / I.P. Kopy'lov. - M.: Vy'ssh. shk.: Logos, 2000. - 607 s. - In Russian.
2. Ivanov-Smolenskij A.V. E'lektricheskie mashiny': ucheb. dlya vuzov. V 2 t. T.1. - 2-e izd., pererab. i dop. / A.V. Ivanov-Smolenskij. - M.: Izd-vo ME'I, 2004. - 656 s. - In Russian.
3. Asinxronny'e dvigateli serii 4A: spravochnik / A.E'. Kravchik [i dr.]. - M.: E'nergoatomizdat, 1982. - 504 s. - In Russian.
4. Popov A.P. Issledovanie poverxnostnogo e'ffekta v provodnikax proizvol'nogo secheniya s pomoshh'yu programmy' Elcut / A.P. Popov, A.O. Chugulyov // Omskij nauch. vestn. - 2009. -#3(83). - S.183-187. - In Russian.
5. Popov A.P. Modelirovanie mikroprocessornogo ustrojstva kontrolya vrashhayushhego momenta asinxronnogo dvigatelya s ispol'zovaniem ego e'kvivalentnoj sxemy' zameshheniya i impul'snogo datchika oborotov / A.P. Popov, M.R. Vinokurov, A.A. Moiseenko // Vestn. Donsk. gos. texn. un-ta. -2011. - T.11, #1. - S.77-82. - In Russian.
ACCURACY INCREASE OF CALCULATING INDUCTION MOTOR TORQUE BASED ON SURFACE EFFECT IN ROTOR BARS
M.R. VINOKUROV
(Don State Technical University),
A.A. MOISEYENKO
(Omsk State Technical University),
N.Y. MASLOVTSEVA
(Rudny Industrial Institute)
A quantitative assessment of the accuracy for determining the mechanical moment of the induction motor when using T equivalent circuit with skin effect in rotor bars compared with L equivalent circuit has been conducted. Keywords: mechanical torque, equivalent circuit, induction motor, slip, eddy currents, ferromagnetic medium.
