CC BY
258
Мониторинг, контрольи диагностика электрооборудования
УДК 621.313 https://doi.org/10.18503/2311-8318-2018-3(40)-70-78
Мугалимов Р.Г. \ Мугалимова А.Р. 2, Калугин Ю.А. \ Одинцов К.Э. 1
1 Магнитогорский государственный технический университет им. Г.И. Носова
2 OOO «МГТУ - Энергосбережение +»
Методика диагностики и идентификации неисправностей обмоток асинхронного
двигателя в режиме его функционирования
В статье выполнен обзор методов диагностики неисправностей асинхронных двигателей. На основе обзора источников информации показано, что отсутсвуют методики, аппаратные и программные комплексы, позволяющие определять неисправности обмоток асинхронных двигателей в режиме их функционирования. В работе предлагается методика, позволяющая определять место неисправностей обмоток АД, а также предсказывать развитие неисправностей. В данной работе дается научное обоснование методики диагностирования развития и идентификации неисправностей обмоток АД, основанное на регистрации реальных фазных напряжений питания, потребляемых токов и моделировании на их основе результирующего вращающегося магнитного поля двигателя, а также энергетические характеристики и последующего их сравнения с результирующим вращающимся магнитным полем двигателя с заведомо исправными и неисправными обмотками по признакам, определяющим искажения вращающегося магнитного поля статора, и идентификаторами неисправностей обмоток. Сформулирована постановка задачи. Задача заключается в разработке диагностических признаков неисправностей обмоток статора и ротора АД по осциллограммам фазных напряжений и токов, зарегистрированных при работе в режимах нормального функционирования электропривода. Критерием решения задачи является получение идентификаторов (признаков) витковых коротких замыканий и обрыва проводников i-x фаз обмотки статора, обрыва стержней короткозамкнутого ротора и ухудшения магнитной системы двигателя. В качестве идентификаторов приняты: координаты вектора магнитного потока статора на комплексной плоскости, определяющие форму его годографа; действующие значения фазных токов, полной, активной, реактивной мощностей, величины КПД и cos ср. Задача решается методом математического моделирования результирующего вращающегося магнитного поля и его составляющих по величинам фазных токов статора, токов намагничивания и токов ротора, а также величин КПД и cos ф асинхронного двигателя, работающего в симметричных и несимметричных режимах. Приводится пример решения задачи с приме-нием предложенной методики и разработанного программного комплекса.
Ключевые слова: асинхронный двигатель, диагностика, моделирование, идентификация, программный комплекс.
Введение
При эксплуатации асинхронных электродвигателей (АД) они аварийно отключаются от источника электропитания, например, из-за перегрузок по току или вращающему моменту, повреждения обмоток статора или ротора, разрушения подшипников или из-за комбинации перечисленных причин. Аварийное отключение электропривода из-за повреждения обмоток или разрушения подшипников АД часто ведет к нарушению реализуемого технологического процесса, простою оборудования, утрате материальных ресурсов и значительным экономическим потерям. Кроме того, как правило, при повреждении обмоток электродвигатель выходит из строя и направляется на капитальный ремонт, что увеличивает эксплуатационные расходы предприятия. В итоге снижается конкурентная способность выпускаемой продукции. Для снижения числа аварийных ситуаций в асинхронных электроприводах применяются различные методы контроля и диагностики электрических и механических параметров, состояний обмоток и подшипников, а также режимов нагрева электродвигателей.
В данной работе предлагается методика и даётся её научное обоснование. Разработана программа диагно-
© МугалимовР.Г., Мугалимова А.Р., Калугин Ю.А., Одинцов К.Э., 2018
обмотки двигателя, неисправности обмоток двигателя,
стирования и идентификации неисправностей обмоток асинхронного двигателя путем регистрации реальных сигналов фазных напряжений, токов и на их основе моделирования формы и координат вращающегося магнитного поля, действующих значений фазных токов, мощностей, КПД и cos ф в режиме функционирования двигателя, затем сравнения полученных параметров с аналогичными параметрами заведомо исправного двигателя и с параметрами-идентификаторами тестовых неисправностей обмоток статора и ротора.
Научно-техническая проблема и её актуальность
По экспертным оценкам в реальном секторе экономики РФ эксплуатируется до 180 млн единиц трехфазных асинхронных электродвигателей [1]. По оценкам специалистов электроремонтных предприятий, цехов и участков, ежегодно около 1,3% (до 1,36 млн единиц) парка АД выходит их строя по причине разрушения обмоток, из-за несвоевременного выявления раз -рушающих факторов и поступает на капитальный ремонт. Разумеется, диагностировать техническое состояние каждого АД не представляется возможным как по методическим, инструментальном, кадровым, так и по экономическим соображениям. Однако наиболее значимые и мощные асинхронные электроприводы, применяемые для реализации непрерывных технологических
процессов, например в металлургической, нефтегазовой, химической и других отраслях производства, а также в технологиях, обеспечивающих жизнедеятельность людей, должны диагностироваться.
Известны методы и средства диагностики технического состояния АД, основанные на измерении величин или токов, или температуры, или вибрации, или акустических характеристик электрических машин. Например, методика диагностики, основанная на непрерывном контроле температуры конкретной точки обмотки статора, сравнение её с температурой окружающей среды, определения превышения температуры над окружающей средой и выработки сигнала на отключение двигателя от источника электропитания [2, 3]. Однако эта методика не идентифицирует причину перегрева обмотки, т.е. не различает из-за какого фактора обмотка перегревается: из-за технологической перегрузки, из-за витковых замыканий или обрыва фаз обмотки статора (ротора) или из-за износа подшипников, или других причин.
Известна также методика, основанная на измерении фазных токов и его спектрального анализа, позволяющая диагностировать состояние подшипниковых узлов [4-7]. Сущность методики заключается в периодической регистрации и выполнении спектрального анализа потребляемых токов и получении спектральной характеристики токов АД. Она сравнивается с аналогичными спектральными характеристиками токов АД, полученными при работе с заведомо исправными подшипниками.
Идентифицировать неисправности обмоток статора (ротора) позволяет методика, основанная на регистрации токов АД и исследовании их спектра с применением вейвлет-преобразований [8]. Известна методика оценки параметров АД путем косвенного измерения сопротивления обмотки статора в режиме функционирования, основанная на формировании сигнала, пропорционального падению напряжения в меди обмотки статора [9]. Однако эта методика позволяет контролировать изменение величины сопротивления в функции температуры и не позволяет диагностировать вид неисправности обмоток АД. Известны методика и прибор для контроля энергетических характеристик асинхронных электроприводов в режиме их функционирования [10]. Однако они не выявляют неисправности электромагнитной системы АД.
Приведенные методики диагностирования АД по потребляемым токам основаны на применении специальных компьютерных программ и требуют достаточно сложного математического обеспечения.
В данной работе дается научное обоснование методики диагностирования развития и идентификации неисправностей обмоток АД, основанное на регистрации реальных фазных напряжений питания, потребляемых токов и моделировании на их основе результирующего вращающегося магнитного поля двигателя, а также его энергетических характеристик и последующего их сравнения с результирующим вращающимся магнитным полем и энергетическими характеристиками двигателя с заведомо исправными и неисправными обмотками по признакам, определяющим искажения вращающегося магнитного поля статора, и идентификаторами неисправностей обмоток.
Постановка задачи
Задача заключается в разработке диагностических признаков неисправностей обмоток статора и ротора АД по осциллограммам фазных напряжений и токов, зарегистрированных при работе в режимах нормального функционирования электропривода.
Критерий решения задачи
Критерием решения задачи является получение идентификаторов (признаков) витковых коротких замыканий и обрыва проводников i-x фаз обмотки статора, обмотки фазного ротора, обрыва стержней короткоза-мкнутого ротора и ухудшения магнитной системы двигателя. В качестве идентификаторов приняты: координаты вектора магнитного потока статора на комплексной плоскости, определяющие форму его годографа; действующие значения фазных токов, полной, активной, реактивной мощностей, величины КПД и cos ф.
Условия, допущения и ограничения
1. Осциллографирование фазных напряжений и токов осуществляется с использованием штатных датчиков напряжений и токов в установившихся режимах при холостом ходе без отсоединения вала АД от вала рабочей машины и при работе с нагрузкой.
2. Продолжительность осциллографирования электрических сигналов от 180 до 300 с. При использовании аналоговых датчиков их сигналы оцифровываются.
3. Используются паспортные (каталожные) данные, параметры магнитной системы тестируемого АД и его трехфазная Т-образная электрическая схема замещения.
4. При исследовании учитывается первая гармоника синусоидального трехфазного напряжения, модуль, фаза и частота которого могут изменятся в пределах, определяемых ГОСТ 9630-80.
5. Трехфазный источник питания обладает большой мощностью, его внутреннее сопротивление бесконечно маленькое.
Задача решается методом математического моделирования результирующего вращающегося магнитного поля и его составляющих по величинам фазных токов статора, токов намагничивания и токов ротора, а также фазных мощностей, величин КПД и cos ф асинхронного двигателя, работающего в симметричных и несимметричных режимах.
Научное обоснование метода
Из классической теории АД известно, что при подаче на его трехфазную обмотку симметричного трехфазного напряжения U1 = U1-e'vi, U2=U2-e7'(vi"120°), U3=U3-e"~vi+12C|0) при идентичных эквивалентных фазных сопротивлениях Z3l=Z32=Z33=Z3-eVz3 в обмотке протекают одинаковые по модулю, но сдвинутые по начальной фазе на 120° токи:
J _ J . jl ).
J _ J jl"Vz -120°) .
2 _ 2 '
J _ J jl-Vz +120°)
13 — 13 t:
В установившихся режимах работы АД симметричная трехфазная система фазных токов статора, при расположении в пространстве магнитных осей фаз обмотки статора со сдвигом на 120 электрических градусов относительно друг друга, возбуждает круговое вращающееся магнитное поле статора. При нарушении симметрии источника питания или симметрии параметров обмоток АД, т.е. при 7Э, ф7зз, что имеет место при неисправностях электромагнитной системы, трехфазная система токов статора становится несимметричной, а круговое вращающееся магнитное поле статора искажается в эллипсоидное вращающееся поле, которое характеризуется координатами - большим и маленьким диаметрами и углом наклона оси большого диаметра относительно вещественной оси комплексной плоскости.
Идея предлагаемого метода диагностирования и идентификации неисправностей обмоток АД основана на формировании годографа вектора вращающегося магнитного поля и его координат, полученных методом моделирования по реальным фазным токам, и сравнении с годографами векторов вращающихся магнитных полей, их координат и энергетических показателей заведомо исправного АД и этого же двигателя с моделируемыми тестовыми неисправностями, число которых Ы=2к, где к=2л; А - число аргументов, определяющих место и вид неисправности обмотки.
Обоснование методики диагностирования неисправностей обмоток АД основано на анализе параметров электрической Т-образной схемы замещения (рис. 1), и их связи с реальными параметрами обмоток статора, ротора и параметрами магнитной системы электрической машины.
На рис. 1 введены обозначения:
и - напряжение /-й фазы источника питания;
7Ч =км 'Кц1+]'к22/Хц1 - комплексное сопротивление фазы обмотки статора;
■м,=кт,Ят1+]-К,Хт1 - комплексное сопротивление фазы ветви намагничивания (главное сопротивление двигателя);
7г=к/ Яг2 +/-Хг2 - комплексное сопротивление фазы обмотки ротора, приведенное к обмотке статора;
Янг=Я„--(1/^-1) - активное сопротивление, пропорциональное механической нагрузке АД, а/ - скольжение АД;
км, к„, км/ - коэффициенты, характеризующие неисправности, соответственно в обмотках статора, ротора и ветви намагничивания.
7 .
и
1
А р
7'
^ п
я„
Рис. 1. Т-образная схема замещения 1-ой фазы АД
Установим связи между величинами сопротивлений схемы замещения с реальными физическими параметрами обмотки статора с учетом возможных её неисправностей. Активное сопротивление Я; фазы в Ом исправной обмотки статора определяется известным соотношением [11]
I
1ср
Я; = w1 — 3
уе-«1 • с, • q-10
(1)
где w1 - число эффективных витков фазы обмотки статора; /1ср - средняя длина витка, мм; у0 - удельная проводимость материала проводника (меди) при рабочей температуре 9, См/мкм; а1 - число параллельных ветвей обмотки статора; с1 - число элементарных проводников в одном эффективном проводнике; q - площадь поперечного сечения элементарного проводника.
Удельная проводимость материала проводника (меди) в Ом при рабочей температуре 0 определяется по известной формуле
7-с1+че" 20°С)],
I V I 20
где у20 - удельная электрическая проводимость материала проводника при 20°С, См/мкм; а - температурный коэффициент сопротивления, °С-1, для меди и
алюминия а=0,004°С-1; 0 - температура обмотки АД.
Из выражения (1) видно, что активное сопротивление Я; фазы обмотки статора зависит от электрической проводимости материала, из которого изготовлена обмотка, числа эффективных витков, средней длины витка, числа элементарных проводников в одном эффективном проводнике, числа параллельных ветвей и от температуры АД. При витковых замыканиях обмотки активное сопротивление Я; уменьшается, а при обрыве элементарных проводников или параллельных ветвей -увеличивается. При математическом моделировании тестовых неисправностей эти признаки неисправностей обмотки статора учитываются коэффициентом
к, =-
к,.
к, • к,.
(2)
где
kwl^= (1-^ ^ к*= (1-* ^ ка,/ = (1-*);
- число
короткозамкнутых витков обмотки статора; си - число разомкнутых (оборванных) элементарных проводников; а1г- - число разомкнутых параллельных ветвей.
При решении поставленной задачи, используя выражение (2), представляется возможным моделировать тестовые короткие замыкания витков, обрывы элементарных проводников и параллельных ветвей обмотки статора. Поэтому при моделировании неисправностей обмоток статора активное сопротивление обмотки изменяется по формуле
Я/ = к/ • я;.
Индуктивное сопротивление исправной обмотки статора Х;, Ом, обусловленное потоками рассеяния магнитного поля, определяется известным соотношением [11]
w
1
т
7
п
х = 1,58 • / -/1 1 р ■ 41-108
(3)
где / - частота тока, Гц; /1 - длина сердечника статора, мм; w1 - число эффективных витков фазы обмотки статора; р - число пар полюсов АД; 41 - число пазов на полюс и фазу; ^ - коэффициент проводимости рассеяния обмотки статора, для АД средней и большой мощности ^=3,5^4,3.
Из приведенного отношения (3) видно, что индуктивное сопротивление обмотки статора при витковых замыканиях обмотки, т.е. уменьшении числа эффективных витков, уменьшается в квадратичной зависимости. Поэтому при моделировании коротких замыканий витков, т.е. при кыф\ изменение индуктивного сопротивления обмотки статора в Ом учитывается выражением вида
X . = к2. ■ X.
(4)
С учетом изложенного комплексное сопротивление 1-й фазы обмотки статора в Ом, при моделировании неисправностей обмотки, представляется векторной суммой
1 .= к ■ 1к2. ■ Хл,
31 31 1 J w■,l 1 '
(5)
где Я]Х1 - параметры обмотки статора, определяемые по каталожными данным или из справочников по АД.
При решении задачи коэффициент к3/ определяется отношением (2) и задается исследователем в зависимости от моделирования вида неисправностей: замыкания витков, числа разомкнутых элементарных проводников, числа разомкнутых параллельных ветвей.
Величина главного индуктивного сопротивления АД, Ом, определяется отношением [11]
х _ Е. _ ке • ^
Хт=г=~г~
(6)
где ке - коэффициент обмотки статора по ЭДС; ке=0,95^0,97; и1н - номинальное фазное напряжение АД, В; 1т - намагничивающий ток 1-й фазы обмотки статора, А.
Намагничивающий ток (А), на один полюс обмотки статора определяется выражением [11]:
!_■ =
2,22 • ^ • р т • w1
(7)
где ^ - МДС обмотки статора или суммарное магнитное напряжение участков магнитной системы АД (воздушный зазор, зубцы статора, зубцы ротора, спинка ротора, спинка статора); т1 - число фаз двигателя.
После подстановки выражения (7) в (6) получим:
Х_ =-
к. -а.
■ т • w1
2,22 • ^ • р
(8)
Из анализа выражения (8) следует, что величина главного индуктивного сопротивления АД, при неизменных геометрических параметрах магнитной систе-
мы, прямо пропорциональна числу витков обмотки статора. Поэтому при решении поставленной задачи величина главного индуктивного сопротивления 1-й фазы АД (Ом) определяется произведением
Х_. = к • Х_
(9)
где Хт - находится по каталожным или справочным данным исследуемого АД.
Комплекс главного сопротивления, 1т, двигателя содержит вещественную составляющую Ят. Эта величина пропорциональна магнитным потерям в электротехнической стали - потерям активной мощности в спинке и зубцах ферромагнитного сердечника статора. При заданных частоте тока, технологии изготовления сердечника, геометрических размерах сердечника и числа пазов статора магнитные потери пропорциональны квадрату магнитной индукции, В, или квадрату фазного напряжения, т.е. РС=Б2 =и\и. При среднем коэффициенте насыщения магнитной системы, равного отношению Р/Рб=1,35^1,77, где - магнитное напряжение в воздушном зазоре АД и номинальном напряжении питания, магнитные потери можно считать неизменными. Поэтому при решении поставленной задачи вещественную составляющую комплекса главного сопротивления АД принято считать не зависимой от числа витков обмотки статора. С учетом изложенного комплексное сопротивление ветви намагничивания 1-й фазы АД в Ом выражаются векторной суммой
1т/ - КтI + 7 ' кж ' ХтI,
(10)
Р„
к ■ Р
кс - коэффициент потерь
где Кт т 1т тК
активной мощности в стали АД; кс=0,0Н0,025.
Электрическая Т-образная схема замещения АД (рис. 1) содержит ветвь обмотки ротора. Комплексное сопротивление ротора, приведенное к обмотке статора, имеет вид:
1г = + JX2,
где К2,Х2 - активное и реактивное сопротивления ротора, приведенные к обмотке статора.
Асинхронный двигатель может изготавливаться с короткозамкнутой обмоткой или с фазной обмоткой ротора. Вещественная составляющая сопротивления обмотки ротора К2 зависит от электрической проводимости материала, числа пазов, числа эффективных витков w2, длины витка /2, числа параллельных ветвей а2, числа элементарных проводников с2, если обмотка фазная, и от числа короткозамкнутых стержней, геометрических размеров сердечника и его пазов, если обмотка короткозамкнутая.
При витковых замыканиях фазной обмотки ротора его активное сопротивление уменьшается с увеличением числа замкнутых витков и увеличивается при обрыве элементарных проводников и размыкании параллельных ветвей. Сопротивление (Ом) определяется произведением
- кг 2 • ^
(11)
где
к,. =-
к
к,- к
(12)
где К21= (1-^ (1-С^ ка2/ = (1-2); ^ - число короткозамкнутых витков обмотки ротора; с2/ - число разомкнутых (оборванных) элементарных проводников; а2/ - число разомкнутых параллельных ветвей. При решении задачи кг2 определяется условиями моделирования неисправностей роторной обмотки АД, т.е. величинами коэффициентов к2, кС2, каг
Для АД с короткозамкнутым ротором возможны обрывы стержней короткозамкнутой обмотки. С увеличением числа оборванных стержней активное сопротивление Я2 увеличивается. При этом величина Яг2 (Ом) определяется отношением
Я - Я2
я . — —,
(13)
где кс2= ^1- ^) - доля исправных стержней; с2/ - число
разомкнутых стержней короткозамкнутого ротора.
Индуктивное сопротивление обмотки ротора, приведенное к обмотке статора, Х2 (Ом) пропорционально квадрату числа витков (стержней) и частоте тока ротора, определяемое скольжением /:
X . = к2 ■ Х2 • • ..
г I w^l 2 /
(14)
При работе АД в режиме реального холостого хода скольжение =5хх=0,005-0,009, а при работе АД с нагрузкой в диапазоне (0,75^1,0)-Р2н, скольжение •=5Н=0,025-0,03. В справочниках по АД, или в каталогах, Х2 определяется для номинального скольжения. Поэтому при моделировании магнитного поля АД при работе без нагрузки, т.е. при холостом ходе, индуктивное сопротивление X. (Ом) снижается до величины, определяемой выражением:
К, = Я2
—• ,
V • .
(18)
При моделировании результирующего магнитного поля учитывается, что при реальном холостом ходе двигатель преодолевает механическую мощность трения в подшипниках и мощность, затрачиваемую на вентиляцию. Из опыта эксплуатации АД механическая мощность реального холостого хода находится в диа-пазоне (0,005-0,009) -Р2н.
Из изложенного и с учетом рис. 1 эквивалентное сопротивление 1-й фазы АД в Ом выражается равенством
7, = 7„
•( 7. + К,)
+ 7. + ЯН1 •
(19)
Поэтому из приведенного анализа величин сопротивлений электрической Т-образной схемы замещения и их связи с реальными физическими параметрами обмоток статора, ротора и параметрами магнитной системы электрической машины трехфазный асинхронный двигатель при работе в симметричных и не симметричных (аварийных) режимах может быть представлен трехфазной электрической схемой замещения (рис. 2).
Для моделирования тестовых режимов работы двигателя расчет приведенной электрической схемы (см. рис. 2), осуществляется по известным выражениям, следующим из закона Ома:
г _ и - иШ .
1 _ 1
г _ и 2 - Цщ ,
2 7 '
э2
г _ У 3 - У Ш ,
з г-г ;
Х. - Х2хх - Кг I ' Х2 ' ~,
2 •„
а при работе с номинальной нагрузкой
Х.= Х 2 = к2. • Х 2.
г, 2хх w,l 2
(15)
(16)
При решении поставленной задачи комплексное сопротивление 1-й фазы обмотки ротора в Ом определяется векторной суммой
7' . = к 2 • К + 1 ■ к2. ■ Х'.
г/ г 2 2 Л 2
(17)
При моделировании коэффициент кг2 задается исследователем.
Для учета механической мощности, приложенной к валу АД, последовательно с сопротивлением 7. включено активное сопротивление Я^ (см. рис. 1). Величина этого сопротивления в Ом пропорциональна механической мощности и определятся известной формулой [11]
У, • — + и2 ■ — + и. ■ —
1 г-г 2 г-г 3 Г7
У _ 7 э1_£э2_7 э3
- 111
-+-+-
7э\ 7э2 7 э3
Величины эквивалентных сопротивлений .-х фаз обмоток АД определяются по формуле (19).
СГ ^ [Ь
Рис. 2. Электрическая схема замещения трехфазного АД
с 2
Моделирование результирующего вращающегося магнитного потока АД основано на известном принципе векторного сложения магнитных потоков, возбужденных токами фаз обмоток статора, оси которых смещены в пространстве относительно друг друга на угол 120 электрических градусов, токами ротора и токами намагничивания [12]. При этом учитывается, что вектор магнитного потока Ф определяется потоком вектора магнитной индукции В через площадь ¡¡:
«5 = ^BdS.
Магнитная индукция является функцией от вектора напряженности магнитного поля Н, относительных магнитных проницаемостей вакуума д0 и ферромагнитного вещества дг:
В = •Н.
Зависимость В(Н) для ферромагнитного вещества нелинейная. На практике она представляется так называемой кривой намагничивания [12].
Напряженность магнитного поля связана с возбуждающими его токами 11, известным законом полного тока:
Yjr ■ W = §Hdl.
При решении поставленной задачи моделирование результирующего вращающегося магнитного поля и его составляющих осуществляется возбуждающими их токами ветвей схемы замещения 1-х фаз АД (см. рис. 1). Токи фаз 1-х ветвей схемы замещения АД определяются по закону Ома:
1. = ^;
31 1.
i =
i =
и.
Z.,
и...
Z „ + R
и = и. -1 • 1..
тп / / 3/
Формирование признаков неисправностей обмоток АД осуществляется компьютерным моделированием по программе исследования.
Программа моделирования числа Ы, тестовых неисправностей обмоток АД формируется с учетом числа аргументов, влияющих на неисправности. Определим число аргументов, вызывающих несимметричный режим работы АД. Для этого все аргументы разделим на две группы.
Группа № 1 - аргументы, определяющие место неисправности: обмотка статора з; обмотка ротора г; фаза статора 3!, фаза статора з2, фаза статора з3; фаза ротора г1, фаза ротора г2, фаза ротора г3; замыкание витков обмотки статора w1l; обрыв элементарного проводника обмотки статора с11; размыкание (обрыв) параллельной
ветви обмотки статора aXi; замыкание витков обмотка ротора w2i; оообрыв элементарного проводника обмотки ротора c2i; размыкание (обрыв) параллельной ветви обмотки ротора a2i. Группа № 1 содержит число аргументов Aj = 14.
Группа № 2 - аргументы, определяющие развития неисправностей в обмотках АД: AwXi, Ac1i, Aa1i, Aw2i, Ac2i, Aa2i. Эта группа аргументов содержит число A2=6. Причем аргументы при моделировании тестовых неисправностей могут задаваться исследователем с шагом: Aw1i=(0,1-0,25)^; Дс1г=(1-7); Aab=(R4); Aw2i=(0,R0,25)w2; Ac2i=(R7); Aa2i=(R4). Максимальные значения аргументов c1, c2, a1, a2 определяются обмоточными данными АД. Величина шага аргументов определяет точность диагностирования развития неисправностей обмоток АД.
Таким образом, место неисправностей обмоток АД определяется числом тестов NAl=2Al=216=65536. Минимальное число тестов, определяющих развитие неисправностей обмоток АД при Aw1i =w1, Aw2i =w2, cb=1, c2i=1, a1i=1, a2i=1, равно NA2=2A2=26=64. Максимальное число тестов, определяющих развитие неисправностей, зависит от заданного шага изменения аргументов A2. Аргументы cXi, c2i, a1i, a2i - целые числа.
После формирования тестовых годографов векторов вращающихся магнитных полей и энергетических показателей АД выполняется сравнение реального годографа и энергетических показателей АД, полученных по зарегистрированным фазным напряжениям и токам. Если координаты годографа вектора вращающегося магнитного поля и энергетические показатели отличаются от годографа вектора магнитного поля и энергетических показателей исправного двигателя не более чем на 3-5%, то двигатель будет считаться исправным. Если показатели режимов работы АД отличаются больше чем на 5%, то выполняется сравнение реального годографа и энергетических показателей АД, полученных по зарегистрированным фазным напряжениям и токам с аналогичными параметрами тестовых неисправностей АД. В результате сравнения параметров выявляется номер тестовой неисправности обмоток АД. Признаки неисправностей обмоток при моделировании симметричных и несимметричных режимов работы АД формируются по алгоритмам [13-15].
Пример реализации методики диагностики
и идентификации неисправностей обмоток двигателя
Моделируются неисправности обмотки статора двигателя с номинальными данными: Р2н=30 кВт, UlH=220 В;/1=50 Гц; п2н=2945 об/мин; ^=3000 об/мин;
=0,905; cos ф =0,9; обмоточные данные, o.e.: г,=0,003, xs=0,073, rm=0,928, xm=3,8, r=0,018, xr=0,11 [16]. При осциллограф ировании зарегистрированы фазные токи ij=98-sin(a>-/-23°), i2=91-sin(ra-/+100°), i3=89-sin(ra-M44°).
Для исправного двигателя вычислены параметры: IlH = 55,8 A; ZH=3,942 Ом; Rs=0,012 Ом, X,=0,288 Ом, Rm=3,658 Ом, Xm=14,98 Ом, Rr=0,07 Ом,ХГ=0,434 Ом.
Программа исследования неисправностей обмоток статора и ротора АД определяются таблицами истинности, табл. 1 и 2, комбинациями их функций f и fr, а также приращениями коэффициентов ksi, kri.
Таблица 1
Программа формирования неисправностей обмотки статора
№ Вход Выход
п/п ¿2 fs
0 0 0 0 J0=S3'S2'S1
1 0 0 1 f =£3 S2 -Si
2 0 1 0 f2=S3'S2'S1
3 0 1 1 f3=S3'S2'S\
4 1 0 0 J4=S3 'S2 'S1
5 1 0 1 f5=S3'S2'S1
6 1 1 0 J6=S3'S2'S1
7 1 1 1 j7=S3 *S2 -S1
Таблица 2 Программа формирования неисправностей обмотки ротора
№ Вход Выход
п/п r3 r2 r1 fr
0 0 0 0 J II 1
1 0 0 1 1f 11 г
2 0 1 0 f2=r'r2 r1
3 0 1 1 j II г
4 1 0 0 f4=r3'r2'r1
5 1 0 1 J5=r3-r2-r1
6 1 1 0 f6=r3-r2'r1
7 1 1 1 J7=r3 r2 r1
В табл. 1, 2 строки под номером 0 соответствуют заведомо исправному асинхронному двигателю.
Решение задачи осуществляется по алгоритму:
1. Задаются исходные данные двигателя и оцифрованные сигналы инструментально зарегистрированных фазных напряжений и токов.
2. Вычисляются параметры Т-образной схемы замещения АД.
3. Задаются коэффициенты ksi и kn-, характеризующие неисправности обмоток АД.
4. Для заданных неисправностей обмоток вычисляются токи ветвей электрической схемы (см. рис. 2), проверяются законы Кирхгофа, баланс мощностей, определяются действующие значения фазных токов, мощностей, КПД, cos ф; создается таблица тестовых данных.
5. Для комбинаций строк табл. 1, 2 и коэффициентов ksi и kri формируются годографы векторов результирующих вращающихся магнитных полей двигателя и составляются серии энергетических данных (координаты главных осей годографа, действующие значения фазных токов, мощностей, КПД, cos ф).
6. Главные координаты вращающихся магнитных полей, соответствующие коэффициентам, характеризующим неисправности обмоток АД, строки 1-7 табл. 1 и 2, сравниваются с соответствующими координатами нулевой строки - координатами магнитных полей, КПД и cos ф заведомо исправного двигателя.
7. Формируется результирующее вращающееся магнитное поле и энергетические параметры двигателя по инструментально зарегистрированным фазным токам.
8. Осуществляется сравнение параметров реального результирующего вращающегося магнитного поля и энергетических характеристик с параметрами вращающихся магнитных полей и энергетических характеристик, полученных моделированием для заданных неисправностей обмоток АД.
9. Определяются совпадающие признаки параметров неисправностей реальной обмотки АД с параметрами, полученнымиметодом моделирования, и идентифицируется номер комбинации неисправностей.
С помощью разработанной методики и программ -ного комплекса при заданных А^1=0,25, Ас1 = 1, Да1 = 1, Д^2=1, Ас2=1, Ая2=1 смоделированы годограф вектора вращающегося магнитного поля исправного двигателя (рис. 3, сплошная линия) и годограф вектора вращающихся магнитного поля, сформированного по осциллограммам токов (рис. 3, пунктирная линия). Сравнительный анализ годографов векторов магнитных полей показывает следующие идентификаторы (показатели): круговое вращающееся магнитное поле для зарегистрированных фазных токов формирует эллипсоидный годограф с большим диаметром, равным 1,16 ед. и маленьким диаметром, равным 0,94 ед. от диаметра кругового вращающегося поля исправного двигателя; ось большого диаметра расположена под углом ф=155°; действующие значения фазных токов, А: ^=69,6, /2=64,5, /3=63,1.
Идентификаторы (показатели) исправного двигателя: годограф круговой с диаметром равным 1,0; угол ф=0°; действующие значения фазных токов, А: 11=61,84, /2=61,84, /3=61,84.
Сравнение идентификаторов годографов исправного и испытуемого АД показывает, что отличительные признаки соответствуют строке под номером один табл. 1 и строке под номером нуль для табл. 2, что соответствует короткому замыканию витков в первой фазе обмотки статора. Идентификаторы развития неисправностей свидетельствуют о замыкании витков в пределах (0,2-0,25)^. Точность предсказания может быть повышена изменением шага Д^ и предварительной обработкой тестовых данных.
100
50
-50
-100
-100
-50
0
50
100
Рис. 3. Годографы вращающихся магнитных полей: для исправного и для испытуемого двигателей
0
Заключение
1. Разработана методика диагностирования и идентификации неисправностей обмоток АД в режиме нормального функционирования, основанная на построение годографов векторов вращающихся магнитных полей и моделирования энергетических характеристик двигателя, работающего в симметричных и несимметричных режимах.
2. Разработаны алгоритм и программный комплекс, реализующие методику диагностирования и идентификации неисправностей обмоток АД.
3. Методика и программный комплекс позволяют диагностировать и идентифицировать неисправности обмоток АД с точностью не ниже 86,6%.
4. Методика и программный комплекс рекомендуются к применению в службах главного энергетика, электрикам цехов промышленных предприятий для периодического мониторинга функционирующих асинхронных электроприводов.
список литературы
1. Электротехника: учеб. пособие для вузов. Э45 / под ред. B.C. Панитюшина. Изд. 2-е, перераб. и доп. М.: Высш. шк., 1976. 560 с.
2. ГОСТ 8865-93. Системы электрической изоляции. Оценка неисправности и классификации.
3. ГОСТ 7217-87. Машины электрические вращающиеся. Двигатели асинхронные. Методы испытаний.
4. Барков A.B., Борисов A.A. Современные возможности диагностирования машины с электроприводом по току двигателя // Металлургические процессы и оборудование. 2013. № 1 (31). C. 61-65.
5. Obaid R.R., Habelter T.G., Stack J.R. Stator current analysis for bearing damage detection in induction motors // The 4th IEEE International symposium on diagnostics for electrical machines, power electronics and drives, SDEMPED 2003. Proceedings. New Jersey. 2003. P. 182-187.
6. Silva J.L.H., Cardoso A.J.M. Bearing failures diagnosis in three-phase induction motors by extended Park's vector approach // The 31st Annual Conference of IEEE Industrial Electronics Society (IECON). 2005. P. 2591-2596.
7. Onel I.Y., Dalci K.B., Senol I. Detection of outer raceway bearing defects in small induction motors using stator current analysis // Sadhana-Academy Proceedings in Engineering Sciences. 2005. Vol. 30 (6). P. 713-722.
8. Купцов В.В., Петушков М.Ю., Сарваров А.С. Современные методы диагностирования асинхронных двигателей и их развитие: монография. Магнитогорск: ГОУ ВПО «МГТУ», 2010. 247 с.
9. Брякин И.В., Бочкарев И.В., Келебаев К.К. Диагностика параметров электрических машин переменного тока // Электротехнические системы и копмлексы. 2017. №4 (37). С. 38-44.
10. Малацион А.С., Малацион Н.В. Контроль энергетических характеристик асинхронных электроприводов с преобразователями частоты в условиях производства // Труды IX Международной (XX Всероссийской) конференции по автоматизированному электроприводу. АЭП-2016. Пермь, 3-7 октябрь 2016. Пермь: Изд-во Перм. нац. иссл. политехи. ун-та 2016. С. 68-70.
11. Гурин Я.С., Кузнецов В.И. Проектирование серий электрических машин. М.: Энергия, 1978. 480 с.
12. Бессонов Л.А. Теоретические основы электротехники. М.: Высшая шк., 1964. 480 с.
13. Bentley J.L. Multidimensional binary search trees used for associative searching // Communications of the ACM. 1975. Vol. 18 (9). P. 509-517.
14. Breiman L., Friedman J.H., Olshen R.A., Stone C.J. Classification and regression trees. Monterey, CA: Wadsworth & Brooks/Cole Advanced Books & Software. 1984. P. 368.
15. Cortes C., Vapnik V. Support-vector networks // Machine Learning. 1995. Vol. 20. P. 273-297.
16. Асинхронные двигатели серии 4A: справочник / Крав-чикА.Э, Шлаф М.М., Афонин В.И., Саболенская Е.А. М.: Энергоиздат, 1982. 504 с.
Поступила в редакцию 19 мая 2018 г.
Information in English
Method of Inductions Motor Malfunction Diagnostics and Identification During Operation
Rif G. Mugalimov
D.Sc. (Eng.), Associate Professor, Professor, Department of Industrial Power Supply, Nosov Magnitogorsk State Technical University, Magnitogorsk, Russia. E-mail: Energosberegenie@rambler.ru.
Aliya R. Mugalimova
Ph.D. (Eng.), Engineer, Limited Liability Company MSTU - Energy Saving +, Magnitogorsk, Russia. E-mail: Energos-beregenie@rambler.ru.
Yurij A. Kalugin
Undergraduate student, Nosov Magnitogorsk State Technical Unversity, Magnitogorsk, Russia. E-mail: py-ro74rus@gmail. com.
Konstantin E. odintsov
Ph.D. (Eng.), Associate Professor, Department of Industrial Power Supply, Nosov Magnitogorsk State Technical University, Magnitogorsk, Russia.
The article contains an overview of methods for diagnostics of malfunctions in induction motors (IMs). Based on the review of information sources, it is shown that there are no methods, hardware or software complexes that make it possible to determine the malfunctions of windings of induction motors during their operation. In this paper, the technique that allows us
to determine the location of faults in the windings of an IM, as well as to predict the development of faults, is purposed. In this paper, scientific justification for the procedure of diagnostics for the development and identification of malfunctions in the windings of AD is given. The procedure is based on recording of real phase supply voltages, consumed currents, and simulating the
resultant rotating magnetic field of the motor on their basis, as well as energy characteristics and their subsequent comparison with the resulting rotating magnetic field of the engine with obviously serviceable and faulty windings on the grounds determining the distortions of the rotating magnetic field of the stator, and by the fault indicators of the windings. The investigation problem is formulated. The problem consists in developing diagnostic signs of malfunctions in the stator winding and rotor windings from the oscillograms of phase voltages and currents registered during normal operation of the electric drive. The criterion of problem solving is obtaining the identifiers (signs) of the short-circuiting malfunctions and conductor breaking of stator winding i-th phases, short-circuited rotor rod breaking and deteriorating of the motor magnetic system. The coordinates of the vector of the stator magnetic flux on the complex plane, which determine the shape of its hodograph; effective phase currents values, full, active and reactive powers, energy conversation efficiency and cos 9 values are taken as identifiers. The problem is solved by mathematical modeling of the resultant rotating magnetic field and its components in terms of the magnitude of the stator phase currents, magnetizing currents and rotor currents, as well as the values of the energy conversation efficiency and cos 9 of the asynchronous motor operating in symmetric and asymmetric modes. An example of the solution of the problem with the application of the proposed method and the developed software package is given.
Keywords: asynchronous motor, diagnostics, motor windings, motor winding malfunctions, simulation, identification, software complex.
References
1. Elektrotekhnika [Electrical engineering]. Textbook for universities. Ed. V.S. Panityushin. Moscow, 1976. 560 p. (In Russian)
2. GOST 8865-93. Electrical insulation systems. Malfunction evaluation and classification. (In Russian)
3. GOST 7217-87. Electric rotating machines. The engines are asynchronous. Test methods. (In Russian)
4. Barkov A. V., Borisov A. A. Possibilities of diagnostics for machines with electric drive using motor current. Metallur-gicheskie protsessy i oborudovanie [Metallurgical processes and equipment]. 2013, no. 1 (31), pp. 61-65. (In Russian)
5. Obaid R.R., Habelter T.G., Stack J.R. Stator current analysis for bearing damage detection in induction motors. The 4th IEEE International symposium on diagnostics for electrical
machines, power electronics and drives, SDEMPED 2003. Proceedings. New Jersey, 2003, pp. 182-187.
6. Silva J.L.H., Cardoso A.J.M. Bearing failures diagnosis in three-phase induction motors by extended Park's vector approach. The 31st Annual Conference of IEEE Industrial Electronics Society (IECON). 2005, pp. 2591-2596.
7. Onel I.Y., Dalci K.B., Senol I. Detection of outer raceway bearing defects in small induction motors using stator current analysis. Sadhana-Academy Proceedings in Engineering Sciences. 2005, vol. 30(6), pp. 713-722.
8. Kuptsov V.V., Petushkov M.Yu., Sarvarov A.S. Sovremen-nye metody diagnostirovaniya asinkhronnykh dvigatelej i ikh razvitie [Modern methods of asynchronous motor diagnostics]. Magnitogorsk: Nosov Magnitogorsk State Technical University, 2010. 247 p. (In Russian)
9. Bryakin I.V., Bochkarev I.V., Kelebaev K.K. Diagnostics of AC electrical machines. Elektrotekhnichskie sistemy i kop-mleksy [Electrotechnical systems and complexes]. 2017, no. 4(37), pp. 38-44. (In Russian)
10. Malatsion A.S., Malatsion N.V. Control of power characteristics of induction motors equipped with frequency converter in manufacturing environment. Trudy IX Mezhdunarodnoj (XX Vserossijskoj) konferentsii po avtomatizirovannomu eh-lektroprivodu. AEHP-2016 [Proceedings of IX (XX All Russian) conference on automatic electric drive]. Perm, 3-7 October, 2016. Perm, 2016, pp. 68-70. (In Russian)
11. Gurin YA.S., Kuznetsov V.I., Proektirovanie serij ehlek-tricheskikh mashin [Design of electrical machines]. Moscow: Energiya, 1978. 480 p. (In Russian)
12. Bessonov L.A. Teoreticheskie osnovy ehlektrotekhniki [Fundamentals of electrical engineering]. Moscow: Vysshaya shkola, 1964. 480 p. (In Russian)
13. Bentley J.L. Multidimensional binary search trees used for associative searching. Communications of the ACM. 1975, vol. 18 (9), pp. 509-517.
14. Breiman L., Friedman J.H., Olshen R.A., Stone C.J. Classification and regression trees. Monterey, CA: Wadsworth & Brooks/Cole Advanced Books & Software. 1984. P. 368.
15. Cortes C., Vapnik V. Support-vector networks. Machine Learning. 1995, vol. 20, pp. 273-297.
16. Kravchik A.Eh., SHlaf M.M., Afonin V.I., Sabolen-skaya E.A. Asinkhronnye dvigateli serii 4A [Asynchronous motors of 4A series]. Moscow: Energoizdat, 1982. 504 p.
Мугалимов Р.Г., Мугалимова А.Р., Калугин Ю.А., Одинцов К.Э. Методика диагностики и идентификации неисправностей обмоток асинхронного двигателя в режиме его функционирования // Электротехнические системы и комплексы. 2018. № 3(40). С. 70-78. Шр8:/Мого^/10.18503/2311-8318-2018-3(40)-70-78
Mugalimov R.G., Mugalimova A.R., Kalugin Yu.A., Odintsov K.E. Method of Inductions Motor Malfunction Diagnostics and Identification During Operation. El-ektrotekhnicheskie sistemy i kompleksy [Electrotechnical Systems and Complexes], 2018, no. 3(40), pp. 70-78. (In Russian). https://doi.org/10.18503/2311-8318-2018-3(40)-70-78
