CC BY
86
Власов А.Б. и др. Факторный анализ диагностической модели...
УДК 621.313 : 519.237.7
Факторный анализ диагностической модели тепловизионного контроля электрической машины
А.Б. Власов, Е.А. Мухин, Б.Д. Царев
Морская академия МГТУ, кафедра электрооборудования судов
Аннотация. Для диагностики технического состояния асинхронного двигателя на базе количественной термографии необходимо создание диагностической модели. Проведен факторный анализ и оценено влияние различных факторов на диагностическую модель асинхронного двигателя. Показано, что диагностическая модель устойчива к возможным погрешностям при определении различных параметров, используемых в расчетах.
Abstract. The thermovision diagnostic method of electric motors windings needs to create a diagnostic model for estimating influence of different parameters on stator windings heating. The factor analysis has been carried out and influence of different parameters on the diagnostic model of induction motor has been estimated. It has been shown that the diagnostic model is stable to possible errors during determining various parameters used in calculations.
Ключевые слова: тепловизионный контроль, асинхронный двигатель, факторный анализ, диагностика технического состояния, количественная термография, температура обмотки, тепловая модель
Key words: thermovision control, asynchronous motor, factor analysis, operating conditions diagnostic, quantitative thermography, winding temperature, heat model
1. Введение
Задачей количественной термографии является определение диагностических параметров электрооборудования (в частности электрических машин) на основе термограмм, полученных при тепловизионном обследовании. Вычисленные значения параметров позволяют сделать экспертное заключение о техническом состоянии оборудования. Актуальным является решение проблемы определения температур внутренних элементов конструкции электрических машин в процессе эксплуатации, в том числе определение температуры статорной обмотки.
Диагностика технического состояния электрической машины (ЭМ) на базе количественной термографии требует создания диагностической модели, на основе анализа которой можно получить информацию о техническом состоянии ЭМ в процессе эксплуатации под нагрузкой. В работах (Власов, Мухин, 2011; Власов и др., 2010) представлена тепловая диагностическая модель асинхронного двигателя, позволяющая рассчитывать диагностические параметры технического состояния, в том числе температуру обмотки статора машины на основе количественной термографии, а также программа для ЭВМ, позволяющая оперативно реализовывать данную модель.
Исходными данными для диагностической модели являются геометрические параметры двигателя, характеристики конструкционных материалов, температуры различных областей внешней поверхности станины, состояние окружающей среды. При этом значения некоторых из этих параметров не удается определить с достаточной степенью точности по нескольким причинам:
• влияние погрешности экспериментального измерения температуры различных областей оребренной поверхности машины;
• отсутствие необходимых для расчёта параметров в справочной литературе;
• зависимость некоторых параметров от температуры, изменяющейся во время эксплуатации машины. Таким образом, актуальной задачей является оценка влияния различных факторов (параметров),
учитываемых в диагностической модели, на расчет диагностических параметров технического состояния. Чтобы произвести подобную оценку, необходимо определить функциональную зависимость, описывающую взаимосвязь между различными факторами и агрегированной функцией, которой в данном случае является температура обмотки статора ЭМ.
Для нахождения функциональной зависимости используется методика полного факторного эксперимента (ПФЭ) второй степени для четырех параметров (Блохин и др., 1997; Власов, 2004). Агрегированная функция имеет вид:
Y = bo + Ь1Х1б + Ь2Х2б + Ь3Х3б + ЬХ4б + Ь12Х1бХ2б + Ь13Х1бХ3б + Ь14Х1бХ4б + Ь23Х2бХ3б + Ь24Х2бХ4б + /1ч
+ Ь34Х3бХ4б + ЬпХ1б2 + Ь22Х2б' + Ь33Х3б + Ь44%4б ,
где Хб - кодовые (безразмерные) значения параметров, принимающих на концах интервалов значения ±1; i, j изменяются от 1 до 4.
46
Вестник МГТУ, том 16, №1, 2013 г.
стр. 46-51
2. Диагностическая модель электрической машины
Тепловая диагностическая модель ЭМ позволяет определять диагностические параметры технического состояния ЭМ, в том числе температуру обмотки статора, по данным количественной термографии. Исходными данными для расчета являются температуры локальных областей поверхности станины, полученные в результате обработки термограмм.
Согласно предлагаемому алгоритму расчета, по термограммам определяется наиболее горячая точка поверхности станины, которая расположена в центральной области поверхности станины, и находятся температуры в этой точке на конце ребра и у основания ребра. Также по термограммам определяются температуры на краях станины. Перепад температур на ребре (разница между температурой на конце ребра и у основания ребра) используется для расчета коэффициента теплоотдачи а поверхности станины.
3. Обоснование факторов, влияющих на расчёт диагностических параметров
Актуальным является оценка влияния погрешности измерения температуры поверхности станины статора. Температурная чувствительность тепловизора, используемого при диагностических мероприятиях, составляет 0,1 °C. Для оценки влияния погрешности измерения на расчет диагностических параметров технического состояния в качестве факторов для анализа были выбраны следующие параметры: температура поверхности станины в горячей точке у основания ребра (tmax), разница температур на конце ребра и у основания ребра в этой точке (At).
Для расчета температур различных областей поверхности станины статора разработана электротепловая схема конструктивных элементов АД, учитывающая пути распространения тепловых потоков от нагретого сердечника статора, позволяющая рассчитывать значения температуры различных областей поверхности станины статора с относительной погрешностью не более 7 %.
При расчете диагностических параметров АД важным является определение значения теплопроводности станины статора Лст и теплопроводности изоляции паза обмотки статора X При этом в литературных источниках приводятся лишь примерные значения этих параметров. Для исследуемого АД значение теплопроводности станины статора Лст равно 130 Вт/м-К, а значение теплопроводности пазовой изоляции X равно 0,16 Вт/м-К (Сипайлов и др., 1989; Борисенко и др., 1983).
Таким образом, для анализа были выбраны следующие факторы:
• температура поверхности станины в горячей точке у основания ребра (tmax);
• перепад температур на ребре в горячей точке (At);
• теплопроводность станины статора /.С1:
• теплопроводность пазовой изоляции обмотки статора Ли.
4. Факторный анализ
Определим агрегированную функцию, описывающую взаимосвязь исследуемых факторов и температуры статорной обмотки двигателя, с помощью факторного анализа с использованием центрального композиционного ортогонального плана (ЦКОП).
Выбранные базовые значения факторов и шаги варьирования приведены в табл. 1. Например, базовое значение температуры у основания ребра станины в горячей точке Х1 = 57,4 °С. Для проведения факторного анализа значения управляемых переменных (факторов) Х, приводятся к безразмерным величинам по соотношению:
Хб = (X, - Хо) / АХi, (2)
где Хо - базовые значения факторов; АХ - шаги варьирования факторов.
Общее число опытов N, необходимое для описания процесса, определяется соотношением N = 2m + 2т + 1. При количестве факторов m = 4 число опытов N = 25.
Уравнение, которым описывается совокупное влияние четырех факторов, представляется в виде (1). Известно, что коэффициенты полинома (1) при ортогональном планировании независимы друг от друга, и некоторые из них, если доказана их малость относительно остальных, могут не учитываться в расчетах.
Коэффициенты b полинома (1) рассчитываются по формуле:
N N
Ъг Т X* б Y4 Т XI б- (3)
k= 1 k= 1
где N - количество экспериментов, к - номер эксперимента.
Рассмотрим ортогональный центральный композиционный план для учета четырех факторов Х..Х4. Для этого используем матрицу планирования (табл. 3), в которой приведены значения безразмерных переменных в соответствии с данными табл. 1. Рассчитанные по формуле (3) значения коэффициентов полинома bi приведены в табл. 4.
Таким образом, получен следующий полином:
Y = 89,148 + 0,029Х1б + 1,332Х2б - 0,396Х3б - 2,845Х4б + 0,015ХХ2б + 0,033ХХзб +
+ 0,055Х2Х3б + 0,007Х2бХ4б - 0,015Х3бХ4б - 0,053Х1б
2
0,023Х2б2 + 0,044Х3б2 + 0,519Х4б2
(4)
47
Власов А.Б. и др. Факторный анализ диагностической модели...
Таблица 1. Значения переменных при факторном анализе диагностической модели
Фактор г °с Д, °С 2ст, Вт/м-К X, Вт/м-К
Обозначение X, X, X, X4
Базовое значение 57,4 3 130 0,16
Интервал варьирования 0,1 0,2 20 0,03
Верхний уровень X = +1 57,5 3,2 150 0,19
Нижний уровень X = -1 57,3 2,8 110 0,13
Звездная точка X = 1,414 71,54 7,535 158,26 0,21
Звездная точка X = -1,414 43,26 0,465 101,72 0,12
В табл. 3 даны значения температур статорной обмотки Ь. рассчитанные на основе диагностической модели, и значения функции отклика 7ф, рассчитанные по полиному, описываемым соотношением (1). Также в табл. 3 приведен процент отклонения д функции отклика ) ф от соответствующих значений Ь. Значение процента отклонения функции отклика от расчетных значений составляет от 0 до 0,23 %.
После вычисления коэффициентов предполагаемой модели исследуемого процесса оценим их значимость с помощью критерия Стъюдента. Рассчитаем значения величины Yна основе диагностической модели (табл. 2), при соответствующем варьировании фактораX2 (1-0,3; 1; 1+0,3).
Таблица 2. Расчетные данные величины Yпри соответствующем варьировании фактора Х3
т\п Х(Х=1; X2=1; X3=1; X4=1) Y2(Xi=1; X2=1+0,31; X3=1; X4=1) Y3(X=1; X2=1-0,3; X3=1; X4=1) Y ■* ср s2^}
1 87,46 155,14 57,35 99,9 2508,4
2 91,32 58,84 157,54 102,5 2530,1
3 93,8 164,62 61,95 106,7 2762
4 90,37 58,19 156,18 101,5 2494,8
5 85,54 54,73 148,99 96,4 2309,6
6 93,73 164,68 61,95 106,7 2766,5
7 88,12 156,09 57,82 100,6 2532,6
8 92,94 163,54 61,41 105,9 2734,7
9 84,6 53,98 147,55 95,3 2276,0
10 85,61 54,69 148,97 96,4 2309,8
11 90,37 58,28 156,29 101,6 2497,2
12 93,2 163,73 61,54 106,1 2736,6
13 87,2 154,95 57,22 99,7 2506,5
14 84,6 54,07 147,67 95,4 2278,3
15 88,05 156,16 57,82 100,6 2537
16 91,25 58,88 157,55 102,5 2530
17 88,64 79,85 81,92 83,4 21
18 88,6 79,75 81,84 83,3 21
19 90,73 208,49 63,48 120,9 5939,4
20 86,63 62,27 194,93 114,6106 4986,6
21 88,38 79,39 81,42 83,06432 22,2
22 89,25 80,59 82,72 84,18956 20,3
23 85,84 76,44 78,51 80,26418 24,4
24 93,69 85,44 87,53 88,88893 18,3
25 88,78 66,54 69,85 75,05624 143,9
Определяем среднее значение функции отклика Y^ для каждого опыта (табл. 2):
^р = EY/n, (5)
где n - параллельное число опытов.
Определяем дисперсию воспроизводимости эксперимента - среднее значение всех оставшихся на воспроизводимость эксперимента дисперсий функций отклика в параллельных опытах
s2{Y,} = ^2^р)1(п - 1) = 44116,81 / 24 = 2146,2. (6)
Дисперсия ошибок определения коэффициентов полинома будет равна:
s2 = s2(Y,)/mn = 1838,5 / (25-3) = 28,61, где m - количество проведённых опытов.
Рассчитаем значение 1-параметра по формуле (7) для каждого коэффициента и соответствующей ему дисперсии ошибки определения этого коэффициента:
48
Вестник МГТУ, том 16, №1, 2013 г.
стр. 46-51
ti = b/(S1{Yi})0’5.
Для коэффициента Ь0 определяем значение t-параметра:
t0 = b0/(S2{Yi})0’5 = 89,148/(28,61)0’5 = 16,66.
(7)
(8)
Аналогичным образом определяем значения t-параметров всех коэффициентов полинома bi.
Определяем число степеней свободы v = n(m - 1) = 25-(3 - 1) = 50. Критическое значение 4р параметра Стьюдента (Блохин и др., 1997) для коэффициента риска в = 0,001 и степени свободы f = 50 равно 3,51.
При сравнении значения 4р с соответствующим значением t-параметра можно утверждать (табл. 3), что коэффициенты й1, b2, b3, b4, b12, b13, b14, b23, b24 являются незначимыми. В этом случае эффектом взаимодействия учитываемых параметров можно пренебречь и уточненная модель, описывающая исследуемый процесс, примет вид:
Y = Ь0 + ЬцХ! + Й22Х22 + Й33Х32 + Й44Х42. (9)
Максимальное значение процента отклонения функции отклика от расчетных значений при использовании данного полинома составляет 2,63 %.
После уточнения вида модели необходимо проверить её на адекватность исследуемому процессу. Аппроксимирующий полином (9) имеет 5 членов, определяем теоретические значения функции отклика для каждого условия проведения опыта, соответствующего конкретного его номера (табл. 2). Определяем дисперсию адекватности в соответствии с параметром
Лд = (1/5)2 (У5 - 7фзн)2 = 41,47.
Таким образом, дисперсия адекватности 52ад меньше дисперсии воспроизводимости эксперимента s2{ Y,}, что говорит об адекватности модели вида (9).
Для анализа влияния каждого из факторов на значение диагностического параметра технического состояния - температуру обмотки статора двигателя, наглядно отобразим зависимость температуры обмотки от факторов, построив трехмерные поверхности (см. рис.). Для этого необходимо зафиксировать два фактора (были выбраны базовые значения факторов), а два других фактора будут выступать в качестве переменных.
Анализ полученных данных показал, что погрешность измерения температуры tmax в исследуемых пределах незначительно влияет на результат расчета температуры статорной обмотки двигателя. Напротив, увеличение величины перепада температуры на ребре At от 2,8 °C до 3,2 °C приводит к изменению tоб примерно с 87,1 °С до 90 °C, то есть примерно на 3 % (рис. 1а). Это можно объяснить тем, что увеличение разницы температур на ребре связано с возрастающим тепловым потоком из объема и при постоянном значении коэффициента теплоотдачи поверхности станины двигателя увеличивается температура статорной обмотки 4б-
а)
б)
в)
Рис. Зависимость температуры статорной обмотки t^ от различных факторов: а) от температуры у основания ребра и перепада температуры на ребре; б) от температуры у основания ребра и коэффициента теплопроводности пазовой изоляции; в) от температуры у основания ребра и коэффициента теплопроводности станины статора
Повышение коэффициента теплопроводности станины Дст со 110 Вт/м-К до 150 Вт/м-К (рис. в) приводит к изменению температуры обмотки с 89,2 °С до 88 °С (на 1,4 %). Физически это объясняется тем, что при большей теплопроводности станины двигателя теплопередача через металл станины
49
Власов А.Б. и др. Факторный анализ диагностической модели...
возрастает. Естественно, что при фиксированном тепловом потоке увеличение теплопроводности металла станины приводит к увеличению теплоотвода.
При повышении теплопроводности пазовой изоляции Яи с 0,13 Вт/м-К до 0,19 Вт/м-К (рис. б) также происходит снижение температуры обмотки с 92,1 °С до 86,4 °С (на 6,2 %). Это является следствием того, что при большем значении теплопроводности изоляции теплоотдача во внешнюю среду увеличивается (теплоотвод возрастает) и, следовательно, статорная обмотка нагревается в меньшей степени.
Таблица 4. Расчетные значения коэффициентов полинома Ь,
Коэффициенты Полный набор С учетом малости
bo 89,148
Ь\ 0,029 0
Ь2 1,332 0
Ьз -0,396 0
Ь 4 -2,845 0
Ь12 0,015 0
Ь13 0,033 0
Ь14 0 0
Ь23 0,055 0
Ь24 0,007 0
Ь34 -0,015 0
Ьп -0,053
Ь22 -0,023
Ь33 0,044
Ь44 0,519
5. Заключение
Проанализировано воздействие различных факторов на результаты расчетов температуры статорной обмотки ЭМ по диагностической модели АД.
Полученное аналитическое выражение для агрегированной функции Y может быть использовано для расчета важного параметра - температуры статорной обмотки асинхронного двигателя при различных условиях испытаний. Выявлено влияние наиболее значимых факторов в диапазоне их возможных значений в процессе эксплуатации.
Показано, что наибольшее влияние на расчет значения температуры статорной обмотки двигателя оказывает изменение значения теплопроводности пазовой изоляции. Тем не менее, влияние всех факторов на расчет температуры обмотки статора при их изменении в исследуемых пределах мало, что доказывает устойчивость диагностической модели к возможным погрешностям при определении различных параметров, используемых в расчетах.
Литература
Блохин В.Г., Глудкин О.П., Гуров А.И., Ханин М.А. Современный эксперимент: подготовка, проведение, анализ результатов. Подред. О.П. Глудкина. М., Радио и связь, 232 с., 1997.
Борисенко А.И., Костиков О.Н., Яковлев А.И. Охлаждение промышленных электрических машин. М.,
Энергоатомиздат, 296 с., 1983.
Власов А.Б., Мухин Е.А. Методика расчета температуры обмоток электрической машины на основе количественной термографии. Вестник МГТУ, т. 14, № 4, c.671-680, 2011.
Власов А.Б., Мухин Е.А., Царев Б.Д. Расчет температуры обмотки асинхронного двигателя по данным количественной термографии. Свидетельство о государственной регистрации программы для ЭВМ № 2011616586, регистр. 24.08.2011, Федеральная служба по интеллектуальной собственности, патентам и товарным знакам, 2010.
Власов А.Б. Факторный анализ диагностической модели тепловизионного контроля высоковольтных вводов. Вестник МГТУ, т. 7, № 3, c.429-436, 2004.
Сипайлов Г.А., Санников Д.И., Жадан В.А. Тепловые, гидравлические и аэродинамические расчеты в электрических машинах. М., Высш. шк., 239 с., 1989.
50
Таблица 3. Матрица планирования ортогонального центрального композиционного плана 24
№ X06 X16 X26 X36 X46 X16X26 X26X36 X36X46 X16X46 X16X36 X26X46 X16n2 X26n2 X36n2 X46n2 О О О О 5, %
1 1 1 1 1 1 1 1 1 1 1 1 0,2 0,2 0,2 0,2 87,46 87,46 0,00
2 1 1 -1 -1 -1 -1 1 1 -1 -1 1 0,2 0,2 0,2 0,2 91,32 91,18 0,15
3 1 -1 1 -1 -1 -1 -1 1 1 1 -1 0,2 0,2 0,2 0,2 93,80 93,73 0,08
4 1 -1 -1 1 -1 1 -1 -1 1 -1 1 0,2 0,2 0,2 0,2 90,37 90,28 0,10
5 1 -1 -1 -1 1 1 1 -1 -1 1 -1 0,2 0,2 0,2 0,2 85,54 85,55 0,01
6 1 1 1 -1 -1 1 -1 1 -1 -1 -1 0,2 0,2 0,2 0,2 93,73 93,75 0,02
7 1 -1 1 -1 1 -1 -1 -1 -1 1 1 0,2 0,2 0,2 0,2 88,12 88,08 0,04
8 1 -1 1 1 -1 -1 1 -1 1 -1 -1 0,2 0,2 0,2 0,2 92,94 93,01 0,08
9 1 -1 -1 1 1 1 -1 1 -1 -1 -1 0,2 0,2 0,2 0,2 84,60 84,55 0,06
10 1 1 -1 -1 1 -1 1 -1 1 -1 -1 0,2 0,2 0,2 0,2 85,61 85,51 0,12
11 1 1 -1 1 -1 -1 -1 -1 -1 1 1 0,2 0,2 0,2 0,2 90,37 90,38 0,01
12 1 1 1 1 -1 1 1 -1 -1 1 -1 0,2 0,2 0,2 0,2 93,20 93,17 0,04
13 1 -1 1 1 1 -1 1 1 -1 -1 1 0,2 0,2 0,2 0,2 87,20 87,31 0,12
14 1 1 -1 1 1 -1 -1 1 1 1 -1 0,2 0,2 0,2 0,2 84,60 84,64 0,05
15 1 1 1 -1 1 1 -1 -1 1 -1 1 0,2 0,2 0,2 0,2 88,05 88,11 0,06
16 1 -1 -1 -1 -1 1 1 1 1 1 1 0,2 0,2 0,2 0,2 91,25 91,22 0,03
17 1 1,414 0 0 0 0 0 0 0 0 0 1,2 -0,8 -0,8 -0,8 88,64 88,69 0,06
18 1 -1,414 0 0 0 0 0 0 0 0 0 1,2 -0,8 -0,8 -0,8 88,60 88,61 0,01
19 1 0 1,414 0 0 0 0 0 0 0 0 -0,8 1,2 -0,8 -0,8 90,73 90,60 0,15
20 1 0 -1,414 0 0 0 0 0 0 0 0 -0,8 1,2 -0,8 -0,8 86,63 86,83 0,23
21 1 0 0 1,414 0 0 0 0 0 0 0 -0,8 -0,8 1,2 -0,8 88,38 88,29 0,10
22 1 0 0 -1,414 0 0 0 0 0 0 0 -0,8 -0,8 1,2 -0,8 89,25 89,41 0,18
23 1 0 0 0 1,414 0 0 0 0 0 0 -0,8 -0,8 -0,8 1,2 85,84 85,77 0,08
24 1 0 0 0 -1,414 0 0 0 0 0 0 -0,8 -0,8 -0,8 1,2 93,69 93,82 0,14
25 1 0 0 0 0 0 0 0 0 0 0 -0,8 -0,8 -0,8 -0,8 88,78 88,76 0,02
