CC BY
31
УДК 621.313
DOI: 10.14529/power160211
ТЕОРЕТИЧЕСКИЕ АСПЕКТЫ ПОСТРОЕНИЯ ТЕРМОДИНАМИЧЕСКИХ МОДЕЛЕЙ ЭЛЕКТРОДВИГАТЕЛЕЙ
В.П. Метельков
Уральский федеральный университет им. первого Президента России Б.Н. Ельцина, г. Екатеринбург
Получено обобщенное математическое описание термодинамической модели электрической машины на основе тепловых схем для произвольного количества узлов и охлаждающих сред, учитывающее влияние температуры узлов на величину греющих потерь, выделяющихся в этих узлах. Показано, что введение учета зависимости мощности потерь в узлах термодинамической модели от температуры этих узлов из-за температурной зависимости сопротивлений обмоток не изменяет свойство линейности модели. Аналитически доказано, что при отсутствии зависимости мощности потерь в узлах термодинамической модели от температуры этих узлов изменение во времени свободных составляющих термодинамических процессов в этой модели всегда имеет апериодический затухающий характер, однако при наличии температурной зависимости мощности потерь это положение не всегда справедливо. Аналитически определены условия, при которых в случае наличия температурной зависимости мощности потерь происходит неограниченное увеличение температуры элементов электрической машины без выхода на установившийся уровень.
Ключевые слова: термодинамическая модель, электрический двигатель, тепловая схема, дифференциальные уравнения, собственные числа.
Введение
В настоящее время в экономике России используются десятки миллионов электродвигателей. В большинстве промышленных установок с электроприводом мощность двигателей завышена на 20-40 %, что говорит о существовании возможностей для снижения стоимости электропривода и потребления электроэнергии им за счет более рационального выбора мощности двигателей [1]. С учетом того факта, что классические косвенные методы оценки теплового состояния двигателей не всегда дают адекватный результат [2], возрастает роль численных методов расчета на основе термодинамических моделей (ТДМ). Значение математического моделирования для оценки теплового состояния двигателей повышается также из-за того, что функционирование механической части электропривода может иметь свои особенности, которые часто затруднительно учесть аналитически [3, 4].
Существует два подхода к численному моделированию тепловых процессов в электрических машинах - моделирование на основе тепловых схем (ТС) с сосредоточенными параметрами и на основе использования пространственно-временной дискретизации объекта исследования (FEM- и CFD-моделирование). При всех преимуществах FEM и CFD, в задачах термодинамического моделирования электрических машин, эти подходы часто требуют построения трехмерных моделей [5, 6], что обусловливает очень высокие требования к вычислительным ресурсам и длительное время решения таких задач. Термодинамические модели на основе ТС с сосредоточенными параметрами
предъявляют значительно меньшие требования к вычислительным ресурсам, обеспечивают сравнительно быстрый расчет, что позволяет использовать их при решении задач, связанных с перебором значительного количества вариантов, например при анализе чувствительности модели к изменениям параметров [5, 7, 8]. Также ТДМ на основе ТС позволяют, в некоторых случаях, провести исследования аналитическим путем с целью выяснения влияния особенностей модели на температуру элементов машины в установившихся и переходных режимах. Это обусловливает актуальность создания и исследования таких ТДМ электрических машин.
В настоящей работе ставится задача получения математического описания обобщенной термодинамической модели электродвигателя на основе метода тепловых схем и исследования некоторых свойств и особенностей функционирования этой модели.
1. Математическое описание
термодинамической модели
В литературе приводятся различные формы математического описания моделей электрических машин, пригодные для исследования термодинамических процессов и/или установившегося теплового режима. Так, в [9] рассматривается термодинамическое уравнение для одной греющейся массы, в [10, 11] рассматриваются ТС с несколькими узлами с точки зрения исследования установившегося теплового режима. В [12, 13] получено математическое описание термодинамических процессов для системы, состоящей из двух и трех
греющихся масс, а в [14] предложена система уравнений для расчета тепловых переходных процессов, состоящая из восьми уравнений для статора и семи - для ротора электрической машины. В [15-18] дано математическое описание термодинамических процессов в виде дифференциального уравнения в матричной форме для произвольного количества греющихся масс. Однако во всех упомянутых случаях математическое описание получено для одной температуры охлаждающей среды. В то же время существуют электрические двигатели с использованием жидкостного охлаждения, например асинхронные короткозамк-нутые взрывозащищенные двигатели серии ЭКВ [19], в которых применяется способ охлаждения двигателя ICW37 - водяное охлаждение корпуса. Выпускаются с водяным охлаждением также тяговые электродвигатели. Некоторые крупные двигатели выпускаются с двумя контурами охлаждения - водяным и воздушным. Например, высоковольтные двигатели АЗ-3000-10-1000УХЛ4, АСЗ-3150-10-375УХЛ4, АН32-18-66-12Т3 со степенью защиты 1Р44 имеют воздушно-водяное охлаждение ICW37A97 и ICW37A91, двигатели серии АЗМ и АЗМС (1Р44) используют способ охлаждения ICW37А71 [20]. В подобных случаях возникает вопрос о необходимости учета в математическом описании не единственной температуры охлаждающей среды.
Также в некоторых случаях приходим к ТС с несколькими температурами охлаждающей среды и при необходимости учета постепенного повышения температуры охлаждающего потока внутри машины, как при аксиальной, так и при радиальной системе вентиляции [21, 22], что актуально для крупных машин как защищенного, так и закрытого исполнения. В [22] получены уравнения теплового баланса электродвигателя для произвольного количества температур охлаждающей среды, но эти уравнения приведены в форме, пригодной лишь для исследования установившегося режима.
Получим обобщенное математическое описание термодинамических процессов в электрическом двигателе для произвольного количества греющихся элементов (узлов ТДМ) при произвольном количестве охлаждающих сред. Обозначим количество узлов модели - п, а количество охлаждающих сред - т.
Рассмотрим тепловой баланс для /'-го узла. Количество тепла, выделяющееся за элементарный промежуток времени dt в /-м узле, определяется суммарной мощностью потерь (Д Р^) в нем:
Qli = Д Ъ. (1)
Количество тепла, отданное /-м узлом другим узлам, определяется разностью их температур:
^ = Г}=1 (ег - (2)
(1*1)
где - - номер узла, с которым происходит тепло-
обмен, ]<е {1,п}, 0^ и 0, - температуры узлов, между которыми происходит теплообмен; -тепловая проводимость между узлами / и
Количество тепла, отданное /-м узлом в охлаждающие среды:
Qзl = 5Х=1(0г - еок)К1к^, (3)
где k - номер охлаждающей среды; 0Л - температура ^й охлаждающей среды; %Ык - тепловая проводимость между /-м узлом и ^й охлаждающей средой.
Количество тепла, затраченное на изменение температуры /-го узла, определяется его теплоемкостью С/:
= С^. (4)
Уравнение теплового баланса
Qlí = Q2i + Qзi + Qлi, с учетом (1)-(4), получим дифференциальное уравнение, описывающее изменение температуры /-го узла:
+ ИГ=1(01 - 0ок)^о1к. (5)
После преобразований из (5) можно записать:
= - 0 ( Е"' = 1 Ч/ + ЕГ=1 ^о¿к I + V и*1) )
' ] = 1 0 ^у + Т=1 0ок)^о1к. (6)
Учтем влияние температуры узлов термодинамической модели на мощность греющих потерь в них. Для потерь в обмотках машины используем известную линейную зависимость сопротивления от температуры. Температурная зависимость других видов потерь может быть приведена к линейному виду путем линеаризации в пределах актуального диапазона изменения температуры этих узлов. При объединении нескольких элементов конструкции электрической машины в один узел ТДМ в нем может оказаться более одного источника греющих потерь, каждый из которых имеет свою зависимость мощности потерь от температуры. Поэтому представим суммарную мощность потерь в /-м узле следующим образом:
ДР1 = Е?=1ДР5и(1) [1 + к.01(1) (01 - е^.од)], (7) где 5 - количество источников потерь в /-м узле; ДР^эд - мощность потерь в 1-м источнике /-го узла; кВ1(1) - температурный коэффициент для 1-го источника потерь в /-м узле; 05щ1) - температура для 1-го источника потерь в /-м узле, относительно которой рассчитывается мощность этого источника.
Для тех узлов, где нет зависимости мощности потерь от температуры, принимаем ^¿(¿)= 0.
Запишем (6) с учетом (7):
= ДР1 - 0 I Е"'=1 Ь] + Т=1 К1к -
\ ц*1)
—Е?=1 ^еод) + '"=1 0(8)
где ДР[ = — ^еод^.од] + Ик = 1 0ок^о1к.
Как видим из (8), введение учета зависимости мощностей потерь в узлах модели от температуры при линейной связи этих величин не нарушает линейность дифференциального уравнения для /-го узла (здесь и далее принимаем допущение об отсутствии температурной зависимости теплоем-костей и тепловых проводимостей). Это дает возможность записать (8) в векторно-матричной форме, для чего введем следующие обозначения: АР -матрица-столбец мощностей потерь в узлах модели размерностью пх1; 0 - матрица-столбец температур узлов размерностью пх1; С - диагональная матрица теплоемкостей узлов размерностью пхп; Л - симметричная (= Аji) общая матрица тепловых проводимостей размерностью пхп.
АР =
А p; - i01! iCl
А P2 ; e = 02 ; c =
А pn - -0n- r
n
^21 - "^22"" ^2 n
^nl ^n2 ■■'
Л =
Диагональные элементы в матрице проводимостей Л определяются следующим выражением:
^И = 1 ^о1к + ТГ] = 1 ^Ц —
и*о
— £?=1 АР3ц(1) кщц. (9)
С учетом введенных обозначений уравнение (8) в векторно-матричной форме запишем следующим образом:
С— = АР' + Л0. (10)
dt
Уравнению (10) соответствует векторно-мат-ричная структурная схема, показанная на рис. 1.
Рис. 1. Структурная схема термодинамической модели двигателя при произвольном количестве базовых температур
2. Исследование решения уравнения термодинамической модели
Исследуем свойства решения уравнения (10) в случае, когда привод работает при неизменных параметрах внешней среды, при постоянстве нагрузки и скорости. При этих условиях мощность потерь в узлах модели будет зависеть только от температуры этих узлов. Для выяснения характера процессов в рассматриваемой ТДМ представим уравнение (10) в следующем виде:
% = AB + F, (П)
где А = С_1Л; Р=С1АР .
Запишем общее решение уравнения (11) [23,
24]:
0 = eAtQinitial + £ G(t — s)F(s)ds, (12) где 0initiai - вектор начальных температур узлов модели; G(t — s) = eA(t~s)- матрица Грина размерностью nxn.
Взяв интеграл в уравнении (12), получим:
0 = (Е — eM)®steady + eM®initial, (13)
где 0steady = —A~1F = —Л_1АР - вектор установившихся значений температур, Е - единичная диагональная матрица размерностью nxn.
Характер изменения во времени температуры узлов будет определяться характером изменения свободных составляющих решения уравнения (11), так как вектор внешних воздействий на систему в рассматриваемом случае постоянен (изменение мощности потерь из-за изменения температуры учтено в диагональных элементах матрицы A). Это означает, что вид процессов будет определяться характером и знаком собственных чисел матрицы A [23]. Из практики известно, что свободные составляющие термодинамических переходных процессов в электрических машинах носят апериодический затухающий характер. Выясним, всегда ли справедливо это утверждение. Свободные составляющие переходных процессов могут иметь указанный характер только в том случае, когда все собственные числа матрицы системы являются вещественными и отрицательными.
Сформулируем теорему:
Матрица A = С^Л, где С и Л - симметричные матрицы порядка n, имеет вещественные собственные числа y, ....
Симметричные матрицы имеют вещественные собственные числа, однако матрица A, являясь произведением двух симметричных матриц С1 и Л, не обязательно сама должна быть симметричной, поэтому неочевидно, что ее собственные числа вещественные. Докажем эту теорему, опираясь на подход, предложенный Ф.Р. Гантмахером [24]. Введем вспомогательную матрицу
V = VC (14)
в смысле V V = С и V Vх = E, где E - единичная диагональная матрица.
Правомерность использования (14) в данном случае следует из того, что для положительно определенной матрицы C существует положительно
определенная матрица
Vc, такая, что (Vc) = С [25]. Из физического смысла ясно, что матрица С положительно определена, так как ее диагональные элементы определяются теплоемкостями. Также отметим, что если « - собственные числа
матрицы С, то собственными числами матрицы VC будут числа Jx(c) [25]. Поскольку матрица С диагональная, ее собственные числа равны ее диагональным элементам, которые ненулевые. Это зна-
чит, что матрица C невырожденная, но тогда и VC будет невырожденная, так как все ^Х(С)^0. Поэтому для матрицы V существует и обратная матрица Vх. С учетом (14) и перечисленных выше обстоятельств, A = С4 Л можем записать в виде
А = (VC)"1(VC)"1A(VC)"1VC, (15)
то есть
А = V~1SV, (16)
где
S = V-iAV~i. (17)
Для диагональной матрицы V (поскольку С диагональная) и симметричной матрицы Л, матрица S также является симметричной матрицей. Поскольку S вещественная и симметричная, то она имеет вещественные собственные числа [23]. Выражение (16) показывает, что матрицы A и S связаны преобразованием подобия. Матрицы, обладающие свойством подобия, имеют один и тот же спектр собственных чисел [25], откуда следует вещественность собственных чисел матрицы A. Теорема доказана.
Исследуем вопрос о знаке собственных чисел матрицы A. Рассмотрим сначала процессы в ТДМ без учета зависимости мощности потерь в узлах от температуры. При kei = 0 модуль каждого диагонального элемента матрицы тепловых проводимо-стей Л в соответствии с (9) равен:
М = EE=i К* + YTj=i ,(i = 1.....n), (18)
(J*i)
Из (18) следует, что для матрицы Л выполняется условие диагонального преобладания по строкам. Учитывая, что все элементы в каждой строке матрицы A получены делением элементов соответствующей строки матрицы Л на одно и то же значение теплоемкости (ap = rkiJC1), для элементов матрицы A также выполняется условие диагонального преобладания по строкам:
Ы = EkU^ + E"=i aij,(i = 1.....n). (19)
(J*i)
В соответствии с теоремой Гершгорина [24], каждое собственное число матрицы А всегда расположено в области локализации (в круге Гершго-рина), определяемой соотношением:
xl<D=i^,(i = 1.....п). (20)
(J*i)
Перепишем соотношение (20) с учетом знака диагональных элементов ай следующим образом:
х + Е
k = l aoik + YTj=i 1
U*i)
<Enj=
(21)
(J*i)
Поскольку 2Г=1 ао1к > 0 и Е"=1 а^ > 0, ясно,
а*1)
что условие (21) может выполняться только при отрицательных значениях %^А)к. Причем должно быть
Х(Л)<-ЕГ=1 С^к- (22)
Таким образом, доказано, что при отсутствии зависимости мощности потерь в узлах ТДМ от температуры этих узлов, изменение во време-
ни свободных составляющих термодинамических процессов в этой модели всегда имеет апериодический затухающий характер.
Отметим, что сформулированное выше положение относится не только к термодинамическим процессам в электрических машинах, но имеет и более общий характер. Оно справедливо для любых термодинамических систем с симметричными матрицами тепловых проводимостей. Иногда в тепловых схемах электрических машин используют однонаправленные тепловые проводимости, например, в тепловых схемах двигателей защищенного исполнения [22]. Однако такие тепловые проводимости определяют перетоки тепла по направлению к охлаждающим средам и поэтому входят только в диагональные элементы матрицы Л, не изменяя ее симметричного характера.
3. Влияние температурного изменения
сопротивлений обмоток на свойства модели
Теперь рассмотрим свойства системы с учетом зависимости мощности потерь в узлах от температуры. В этом случае модуль каждого диагонального элемента матрицы Л в соответствии с (9) определяется следующим образом:
М =
ЕГ=1 ^oifc + Е
п i 7 = 1 Ч/
0>Q
Е!=1ДР5и(1) &ег(о|. (23)
С учетом (23) запишем выражение для матрицы А, определяющее область локализации ее собственных чисел:
х + ЕГ=1 aoik+Enj=
1 "¿7
-Cf'ELi^
(j*i) rsti(i) kei(i) I
< En7 =
1 "¿7-
(24)
(J*i)
Анализ неравенства (24) показывает, что при Ег=1 ДРбЫ(1) кщ1)ф0 оно может выполняться как при отрицательных, так и при положительных значениях х^А) (включая и Х(Л)= 0). В любом случае для этого должно выполняться неравенство
Х(А) < СI 1(ESl=lДPst.i(l) кещ) - ЕГ=1 ^огл). (25)
Из (25) видно, что при Ef=1 APstm квт > > ЕГ=1 ^оik может оказаться, что Х(Лр > 0. Это соответствует ситуации неограниченного увеличения температуры 1-го узла без выхода на установившийся уровень.
В таблице приведены значения собственных чисел матрицы A = С1 Л для 6-массовой ТДМ двигателя 4A100L4 (модель включает в качестве своих узлов обмотку статора с разделением на лобовые и пазовую части, сталь статора, внутренний воздух, ротор и корпус).Обозначено: M/MN - относительная величина момента двигателя. Учитывалось влияние изменения температуры на мощность потерь в обмотках статора и ротора.
На рис. 2 и 3 показаны графики изменения превышения температуры (т) лобовых частей об-
Величины собственных чисел матрицы А
Двигатель 4A100L4
M/MN 0 0,5 1,0 1,5 2,0
-361,7809 -361,7808 -361,7807 -361,7805 -361,7802
-10,2185 -10,1960 -10,1285 -10,0160 -9,8586
w100, 1/c -2,5035 -2,4870 -2,4383 -2,3600 -2,2567
-1,0839 -1,0786 -1,0621 -1,0327 -0,9878
-0,0594 -0,0554 -0,0427 -0,0196 +0,0165
-0,2781 -0,2729 -0,2578 -0,2334 -0,2005
т, °С
5000
-5000 / ------------1-------------1-------------1-------------1--------------------------
-10000 -------г --{--- ---}--
-15000-'-'-'-'-1-'-
0 100 200 300 400 500 600 700
t, С
Рис. 2. Графики t(¿) и т (t) при MIMn = 1,0
мотки статора двигателя 4A100L4 и второй производной этой величины (т) при разных величинах нагрузки. Рис. 2 соответствует ситуации, когда все собственные числа матрицы А отрицательны, а рис. 3 - ситуации, при которой одно из собственных чисел имеет положительное значение (см. таблицу) и процесс нагрева двигателя приобретает характер «теплового саморазгона», о чем свидетельствует переход графика второй производной в положительную область.
Выводы
1. Получено обобщенное математическое описание ТДМ электрической машины для произвольного количества узлов и охлаждающих сред, учитывающее влияние температуры узлов на величину греющих потерь, выделяющихся в этих узлах.
2. Показано, что введение учета зависимости мощности потерь в узлах ТДМ от температуры этих узлов при линейной (или линеаризованной) связи этих величин не изменяет свойство линейности модели.
3. Аналитически доказано, что при отсутствии зависимости мощности потерь в узлах ТДМ от температуры этих узлов изменение во времени свободных составляющих термодинамических процессов в этой модели всегда имеет апериодический затухающий характер, однако при наличии температурной зависимости мощности потерь это положение не всегда справедливо. Аналитически определены условия, при которых в случае нали-
т, °С
Рис. 3. Графики т(0 и x(t) при M/Mn = 2,0
чия температурной зависимости мощности потерь происходит неограниченное увеличение температуры элементов электрической машины без выхода на установившийся уровень.
Литература
1. Онищенко, Г.Б. Проблемы и перспективы развития электропривода / Г.Б Онищенко, М.Г. Юньков // Труды VIII Международной (XIX Всероссийской) конференции по автоматизированному электроприводу АЭП-2014: в 2 т. Т. 1. -Саранск: Изд-во Мордов. ун-та, 2014. - С. 5-9.
2. Зюзев, А.М. К оценке теплового состояния асинхронного двигателя в повторно-кратковременном режиме / А.М. Зюзев, В.П. Метельков // Электротехника. - 2014. - № 9.- С. 19-23. DOI: 10.3103/S10683 71214090090
3. Программный комплекс для исследования эксплуатационных режимов электроприводов буровых установок / А.М. Зюзев, В.М. Липанов, В.П. Метельков и др. // Электротехника. - 2003. -№ 7. - С. 25-31.
4. Зюзев, А.М. Математическая модель спус-ко-подъемного агрегата буровой установки / А.М. Зюзев, В.П. Метельков, В.Н. Радченко // Изв. вузов. Горный журнал. - 2003. - № 5. - С. 22-30.
5. Boglietti, A. TEFC Induction Motors Thermal Models: A Parameter Sensitivity Analysis / A. Boglietti, A. Cavagnino, D.A. Staton // IEEE Transactions on Industry Applications. - 2005. - Vol. 41, issue 3. -pp. 756-763. DOI: 10.1109/TIA.2005.847311
6. Chin, Y.K. Transient Thermal Analysis using both Lumped-Circuit Approach and Finite Element Method of a Permanent Magnet Traction Motor / Y.K. Chin, D.A. Staton. // IEEE Africon. - Gaborone, Botswana, 2004. - pp. 1027-1035. DOI: 10.1109/AFRIC0N. 2004.140684 7
7. Evolution and Modern Approaches for Thermal Analysis of electrical machines / A. Boglietti,
A. Cavagnino, D. Staton et al. // IEEE Transactions on Industrial Electronics. - 2009. - Vol. 56, no. 3. -pp. 871-882. DOI: 10.1109/TIE.2008.2011622
8. Tang, W.H. A Simplified Transformer Thermal Model Based on Thermal-Electric Analogy / W.H. Tang, Q.H. Wu, Z.J. Richardson // IEEE Transactions on Power Delivery. - 2004. - Vol. 19, no. 3. -pp. 1112-1119. DOI: 10.1109/TP WRD.2003.822968
9. Копылов, И.П. Проектирование электрических машин /Б.К. Клоков, В.П. Морозкин, Б.Ф. Токарев. - М.: Высш. шк., 2005. - 767 с.
10. Филиппов, И.Ф. Теплообмен в электрических машинах / И. Ф. Филиппов. - Л.: Энергоатом-издат, 1986. - 256 с.
11. Гольдберг, О.Д. Проектирование электрических машин / О.Д. Гольдберг, И. С. Свириденко. -М.: Высш. шк., 2006. - 430 с.
12. Федоров, М.М. Динамические тепловые модели узлов электрических машин / М.М. Федоров // Електромашинобудування та електрообла-днання: Мiжвiд. наук.-техн. зб.- 1999. - Вип. 53. -C. 70-73.
13. Электромеханические и тепловые режимы асинхронных двигателей в системах частотного управления / Р.Т. Шрейнер, А.В. Костылев,
B.К. Кривовяз, С.И. Шилин. - Екатеринбург: ГОУ ВПО «Рос. гос. проф.-пед. ун-т», 2008. - 361 с.
14. Расчетно-теоретическое определение изменения теплового состояния обдуваемых асинхронных двигателей в режиме короткого замыкания и после него /А.Н. Бурковский, О.И. Пилипенко, О.А. Рыбалко, А.А. Мельник // Взрывозащищенное электрооборудование: сб. науч. тр. УкрНИИВЭ. -Донецк, 2007. - С. 183-193.
15. Воробьев, В.А. Нагревание и тепловая защита асинхронных двигателей средней мощности
при периодической нагрузке: дис. ... канд. техн. наук. - Свердловск, 1982. - 181 с.
16. Петрушин, В.С. Тепловые расчеты нестационарных режимов работы асинхронных двигателей регулируемых электроприводов / В.С. Петрушин, А.М. Якимец, В.Л. Кобрин //Електротехтка i електромеханiка. - 2003. - № 4. - С. 65-68.
17. Петрушин, В. С. Особенности тепловых расчетов неустановившихся режимов работы регулируемых асинхронных двигателей / В. С. Петрушин, А.М. Якимец, О.В. Каленик. // Електромашинобудування та електрообладнання: мiжвiд. наук.-техн. зб. - 2008. - Вип. 71. - С. 47-51.
18. Вареник, Е.А. Тепловые переходные процессы в элементах конструкции асинхронных двигателей при неподвижном роторе / Е.А. Вареник, М.М. Федоров, В.Е. Михайлов // Взрывозащищенное электрооборудование: сб. науч. тр. УкрНИИВЭ. -Донецк, 2009. - С. 146-154.
19. Савин, Б.В. Исследования и разработки взрывозащищенных асинхронных двигателей мощностью 0,25-2000 кВт /Б.В. Савин, И.Г. Ширнин // Электротехническая промышленность. Сер. Электрические машины. - 1977. - №10 (80). -С. 7-8.
20. Продукция ООО «ПО «Ленинградский электромашиностроительный завод». Технический каталог [Диск]. - Русэлпром, 2009. - 160 с.
21. Борисенко, А.И. Охлаждение промышленных электрических машин / А.И. Борисенко, О.Н. Костиков, А.И. Яковлев - М.: Энергоатом-издат, 1983. - 296 с.
22. Сипайлов, Г.А. Тепловые гидравлические и аэродинамические расчеты в электрических машинах /Г.А. Сипайлов, Д.И. Санников, ВА. Жадан. -М.: Высш. шк., 1989. - 239 с.
23. Корн, Г. Справочник по математике для научных работников и инженеров / Г. Корн, Т. Корн. - М.: Наука, 1984. - 832 с.
24. Гантмахер, Ф.Р. Теория матриц / Ф.Р. Гант-махер. - М.: Наука, 1966. - 576 с.
25. Ланкастер, П. Теория матриц / П. Ланкастер. - М.: Наука, 1973. - 280 с.
Метельков Владимир Павлович, канд. техн. наук, доцент кафедры «Электропривод и автоматизация промышленных установок» Уральского федерального университета им. первого Президента России Б.Н. Ельцина, г. Екатеринбург; [email protected].
Поступила в редакцию 4 мая 2016 г.
DOI: 10.14529/power160211
THEORETICAL ASPECTS OF DEVELOPING THERMODYNAMIC MODELS OF ELECTRIC MOTORS
V.P. Metelkov, [email protected]
Ural Federal University named after the first President of Russia B.N. Yeltsin, Ekaterinburg, Russian Federation
The paper provides a generalized mathematical description of the thermodynamic model for an electrical machine, which is based on thermal circuits for any number of nodes and cooling media. This description considers an influence of temperature of nodes on the thermal losses at these nodes. Considering the dependence of power losses in the nodes of the thermodynamic model on the temperature of these nodes due to temperature dependency of winding resistance does not change the model linearity. The analysis proves that if there is no dependence of the power losses in the nodes of the thermodynamic model on the temperature of these nodes, intime change of transient components of thermodynamic processes of this model is always of an acyclic damped nature. However, if there is a temperature dependence of the power losses, this provision is not always relevant. The authors have analytically defined the conditions of non-limited temperature rise of electric machine components without switch to steady-state mode (provided there is a temperature dependence of the power losses).
Keywords: thermodynamic model, electric motor, thermal circuit, differential equations, eigenvalues.
References
1. Onishchenko G.B., Yunkov M.G. [Power Drive Systems Development Problems and Prospects]. Trudy VIII Mezhdunarodnoy (XIX Vserossiyskoy) konferentsii po avtomatizirovannomu elektroprivodu [Proc. of the VIII Int. (XIX All_Russian) Conf. on the Automatic Electric Drive]. Saransk, 2014, vol. 1, pp. 5-9.
2. Zyuzev A.M., Metelkov V.P. Toward the Evaluation of the Thermal State of an Induction Motor in the Recursive Short_Term Mode. Russian Electrical Engineering, 2014, vol. 85, no. 9, pp. 554-558. DOI: 10.3103/S1068371214090090
3. Zyuzev A.M., Lipanov V.M., Metelkov V.P., Polyakov V.N., Popov A.S., Shut'ko V.F. The Software Package for the Study of Operating Modes of Electric Drives for Drilling Rigs. Russian Electrical Engineering, 2003, no. 7, pp. 25-31.
4. Zyuzev A.M., Metelkov V.P., Radchenko V.N. Mathematical Model of Lowering and Lifting Unit for Drilling Rig, Izvestiya vuzov. Gornyy zhurnal [Proceedings of the universities. Mining Journal], 2003, no. 5, pp. 22-30.
5. Boglietti A., Cavagnino A., Staton D.A. TEFC Induction Motors Thermal Models: A Parameter Sensitivity Analysis. IEEE Transactions on Industry Applications, May-June 2005, vol. 41, issue 3, pp. 756-763. D0I:10.1109/TIA.2005.847311
6. Chin Y.K., Staton D.A. Transient Thermal Analysis Using Both Lumped-Circuit Approach and Finite Element Method of a Permanent Magnet Traction Motor. IEEE Africon, Gaborone, Botswana, September 2004, pp. 1027-1035. DOI: 10.1109/AFRIC0N.2004.1406847
7. Boglietti A., Cavagnino A., Staton D., Shanel M., Mueller M., Mejuto C. Evolution and Modern Approaches for Thermal Analysis of Electrical Machines. IEEE Transactions on Industrial Electronics, March 2009, vol. 56, no. 3, pp. 871-882. DOI: 10.1109/TIE.2008.2011622
8. Tang W.H., Wu Q.H., Richardson Z.J. A Simplified Transformer Thermal Model Based on Thermal-Electric Analogy. IEEE Transactions on Power Delivery, July 2004, vol. 19, no. 3, pp. 1112-1119. DOI: 10.1109/TPWRD.2003.822968
9. Kopylov I.P., Klokov B.K., Morozkin V.P., Tokarev B.F. Proektirovanie elektricheskikh mashin [Design of Electrical Machine]. Moscow, Vyssh. Shk. Publ., 2005. 767 p.
10. Filippov I.F. Teploobmen v elektricheskikhmashinakh [Heat Transfer in Electrical Machines]. Leningrad, Energoatomizdat, 1986. 256 p.
11. Goldberg O.D., Sviridenko I.S. Proektirovanie elektricheskikh mashin [Design of Electrical Machine] Moscow, Vyssh. Shk., 2006. 430 p.
12. Fedorov M.M. [Dynamic Models of Thermal Units of Electrical Machines]. Electrical Machinery and Electrical Equipment: Collection of Scientific Papers, 1999, issue 53, pp. 70-73.
13. Shreiner R.T., Kostylev A.V., Krivovyaz V.K., Shilin S.I. Elektromekhanicheskie i teplovye rezhimy asinkhronnykh dvigateley v sistemakh chastotnogo upravleniya [Electromechanical and Heat Modes of Induction Motors in Frequency Control Systems]. Ekaterinburg, Russian State Professional Pedagogical University Press, 2008. 361 p.
14. Burkovskiy A.N., Pilipenko O.I., Rybalko O.A., Melnik A.A. Computational and Theoretical Determination of the Change of the Thermal State Windswept Asynchronous Motors in Short Circuit Mode and After. Explosion-proof Electrical Equipment: Collection of Scientific Papers UkrNIIVE. Donetsk, 2007, pp.183-193.
15. Vorobyov V.A. Nagrevanie i teplovaya zashchita asinkhronnykh dvigateley sredney moshchnosti pri periodicheskoy nagruzke: dis. ... cand. tekhn. nauk [Heating and Thermal Protection of Asynchronous Motors of Average Power Under Periodic Load: Tesis for a Cand.'s Degree]. Sverdlovsk, 1982. 181 p.
16. Petrushin V.S., Yakimets A.M., Kobrin V.L.Calculation of Thermal Non-stationary Modes of Operation of Asynchronous Motors for Variable Speed Drives, Electrotechnics and Electrical Engineering, 2003, no. 4, pp. 65-68.
17. Petrushin V.S., Yakimets A.M., Kalenik O.V. Features Thermal Calculations of Transient Modes, of Regulated Asynchronous Motors. Electrical Machinery and Electrical Equipment: Collection of Scientific Papers, 2008, issue 71, pp. 47-51.
18. Varenik E.A., Fedorov M.M., Mikhailov V.E. [Thermal Transients in Structural Elements of Induction-Motors with a Fixed Rotor]. Explosion-proof Electrical Equipment: Collection of Scientific Papers UkrNIIVE. Donetsk, 2009, pp.146-154.
19. Savin B.V., Shirnin I.G. Савин, Б.В. Research and Development of Explosion-proof Asynchronous Motors 0,25-2000 kW. Electrotechnical Industry. A Series of Electrical Machines, 1977, no. 10 (80), pp. 7-8.
20. Production of "Leningrad Electric Machine Plant". Technical Catalog. Ruselprom, 2009. 160 p. [Disk].
21. Borisenko A.I., Kostikov O.N., Yakovlev A.I. Okhlazhdenie promyshlennykh elektricheskikh mashin [Cooling of Industrial Electric Machines]. Moscow, Energoatomizdat, 1983. 296 p.
22. Sipaylov G.A., Sannikov D.I., Zhadan V.A. Teplovye, gidravlicheskie i aerodinamicheskie raschety v elektricheskikh mashinakh [Thermal, Hydraulic and Aerodynamic Computations in Electrical Machines]. Moscow, Vyssh. Shk. Publ., 1989. 239 p.
23. Korn G., Korn T. Mathematical Handbook for Scientists and Engineers. Moscow, Nauka, 1984. 832 p.
24. Gantmaher F.R. Teoriya matrits [Matrix Theory]. Moscow, Nauka, 1966. 576 p.
25. Lankaster P. Teoriya matrits [Matrix Theory]. Moscow, Nauka, 1973. 280 p.
Received 4 May 2016
ОБРАЗЕЦ ЦИТИРОВАНИЯ
FOR CITATION
Метельков, В.П. Теоретические аспекты построения термодинамических моделей электродвигателей / В.П. Метельков // Вестник ЮУрГУ. Серия «Энергетика». - 2016. - Т. 16, № 2. - С. 82-89. DOI: 10.14529/power160211
Metelkov V.P. Theoretical Aspects of Developing Thermodynamic Models of Electric Motors. Bulletin of the South Ural State University. Ser. Power Engineering, 2016, vol. 16, no. 2, pp. 82-89. (in Russ.) DOI: 10.14529/power160211
