CC BY
40
Электротехнические комплексы и системы
65
ЭЛЕКТРОТЕХНИЧЕСКИЕ КОМПЛЕКСЫ И СИСТЕМЫ
УДК 612
А.В.Губенков
МОДЕЛЬ АСИНХРОННОГО ЭЛЕКТРОДВИГАТЕЛЯ С КАБЕЛЕМ И УСТРОЙСТВОМ КОММУТАЦИИ В СТАТОРНОЙ ЦЕПИ
В сетях электроснабжения с электродвига-тельной нагрузкой имеют место переходные процессы, которые могут изменять питающее напряжение электродвигателей. При этом возможно возникновение неблагоприятных условий работы коммутационной аппаратуры, что может привести к ее отключению и что, в свою очередь, может привести к распаду всего технологического процесса.
В этой связи необходим анализ условий работы коммутационной аппаратуры, а также тех последствий, которые возникают при внезапном прекращении электроснабжения одного из потребителей в сети с электродвигательной нагрузкой.
Рассмотрим модель асинхронного двигателя с кабелем и устройством коммутации (пускателем) в статорной цепи. На рис. 1 отображена схема рассматриваемой модели.
ного асинхронного электродвигателя с кабелем в статорной цепи. Модель имеет следующий вид:
+і'
Рис. 1- Схема АД с кабелем и устройством коммутации в статорной цепи
Коммутирующее устройство и асинхронный двигатель получают питание через протяжённый кабельный участок. Коммутирующее устройство имеет произвольное расположение на кабельном участке и своими главными контактами делит его на 2 части с сопротивлениями Кк1, Ьк1 и Як2, Ьн.
В [1] описана математическая модель одиноч-
Ls \+—'
V Ls ;
DWsa — U а +
+ kr —— • (— Rrira — P®VrP )— (k + Rs У sa,
D^sp= Up + (1)
kr — • (— Rrirp + p°Wra)—(Rk + Rs Xp,
L
s
DWra — —Rrira — p°4/rp,
DWrP —— Rrirp + p°4/ra ■
где Rk, Lk - активное и индуктивное сопротивление участка кабеля; Rs- активное сопротивление статора; L’s- переходная индуктивность; Rr- активное сопротивление ротора; Ua, Uр - мгновенные значения напряжений вторичной обмотки трансформатора по неподвижным осям а и Р; isa isp ira irp - мгновенные значения токов статора (индекс s) и ротора (индекс r) по осям а и Р; p -число пар полюсов; о — угловая частота вращения ротора; Щр, \^rca Щр — потокосцепления статора (индекс s) и ротора (индекс r) по осям а и р. Знак D перед переменными в уравнениях обозначает первую производную по времени этой переменной.
Для того, чтобы включить в (1) модель коммутационного устройства, необходимо определиться со следующими условиями и ограничениями:
1. Расположение коммутационного устройства на кабельном участке может быть произвольным. На рис. 1 показаны точки расположения пускателя. После вторичной обмотки трансформатора и участком кабельной сети, точка А; непосредственно на участке кабельной сети, точка B; рядом с асинхронным двигателем, после участка кабельной сети точка С;
2. Влиянием тока управления коммутационного устройства можно пренебречь, т. к. пусковой и рабочий ток коммутационного аппарата на 3-4 порядка меньше соответственно пускового и рабочего тока асинхронной машины.
Ua
Ub
3. Асинхронный двигатель и коммутационный аппарат получают питание через один кабельный участок, следовательно, начиная с момента включения и заканчивая моментом отключения коммутационного аппарата, они электрически связаны и влияют на работу друг друга. Напряжение на зажимах коммутационного аппарата является функцией напряжения на зажимах асинхронного двигателя.
Рассмотрим условия Кирхгофа (рис. 1) для напряжений контура: вторичная обмотка трансформатора - кабельный участок - коммутирующее устройство - статорная обмотка электродвигателя [1].
Щ а ~ Щка + Щкуа + Щ.ва, Щр = Щкр + ЩкуР + ир.
(2)
соответственно падение
напряжения в кабельном участке (индекс к) и падение напряжения в статоре АД (индекс 5) по неподвижным осям а и р.
Падение напряжение на всём кабельном участке складывается из падений напряжения на первом и втором кабеле ик=иы+ик2. и^а, и^р - падение напряжения на сопротивлении межконтакт-ного промежутка коммутирующего устройства по неподвижным осям а и р.
Динамические процессы в кабельной сети описываются [1]
ика= + КкКа,
(3)
Ж йївр
икр = Ьк ~~Т ^ Кк^вр.
Ж
где Як,Ьк - активное и индуктивное сопротивление общего участка кабеля.
Напряжение, прикладываемое к обмотке контактора коммутирующего устройства, выражается формулой:
и^-Ща — Щк,
(4)
Разность берётся по абсолютной величине т.к. модели последних пускателей работают на выпрямленном напряжении.
Коммутационное устройство также влияет на работу асинхронной машины. Рассмотрим случай, когда при снижении напряжения на зажимах обмотки электромагнита коммутирующего аппарата, по абсолютной величине, ниже итп происходит отключение коммутирующего аппарата и как следствие отключение подачи напряжения на статор асинхронной машины, и, наоборот, при превышении итЫ - включение электромеханической системы.
Будем различать процессы включения и отключения коммутационного устройства. При включении аппарата подвижные и неподвижные контакты соударяются, в результате упругой деформации контакты размыкаются, но под действием контактной пружины замыкаются вновь, процесс носит затухающий характер. Данное яв-
ление называется вибрацией контактов (дребезг контактов). При отключении аппарата возникает электрическая дуга в межконтактном промежутке.
В момент включения аппарата в силовой цепи происходит резкое уменьшение сопротивления коммутирующего органа 2к и быстрый рост тока нагрузки.
При отключении аппарата, в межконтактном промежутке возникает электрическая дуга, которая имеет нелинейное электрическое сопротивление, зависящее в основном от величины протекающего в ней тока.
Дуговые процессы в электрических аппаратах до сих пор не имеют строгого математического описания [4] поэтому необходимо ввести функцию, имитирующую сопротивление электрической дуги при отключении коммутирующего устройства и определяющую сопротивление межкон-тактного промежутка коммутирующего органа 2к.
В общем случае данная функция имеет сложный нелинейный характер, который зависит от сопротивления электрической дуги во время коммутации, а также восстанавливающейся электрической прочности межконтактного промежутка и восстанавливающегося напряжения. Косвенный учёт этих трёх факторов будем производить с помощью функции, которая не имеет разрыва в момент коммутации.
аг^( - (к - 0.005)-107 + 12 (5)
ж ' 2
где ( - время независимая переменная; (к- момент времени, в который начинает происходить переключение аппарата; коэффициент в уравнении 107- определяет скорость изменения функции от одного состояния до другого (угол наклона);
0.005 - промежуток времени, за который происходит коммутация. Общий вид функции представлен на рис. 2.
Рис. 2- Функция изменения сопротивления коммутирующего аппарата
Учитывая, что комплексное сопротивление межконтактного промежутка содержит активную и индуктивную составляющую
/2 2 1ку =дЯку + ) , функцию, описывающую
процесс коммутации, можно представить двумя системами уравнений:
4.95x10
Яку (хкі) —
0.001, если хк( < хтіп, (
10е
аеЩ(і - ік - 0.005) ■ 10
(6)
если, Хкі — хтп .
0.001, если Хк^ < хтп,
( 7
acrtg(і - ік - 0.005) ■ 10
20 ■
если, хкі - хт
(7)
где хи - расстояние между якорем контактора и сердечником магнитной системы (отсчитывается от положения якоря при минимальном зазоре); хтіп - минимальное расстояние, при котором главные контакты контактора остаются замкнутыми.
Падение напряжения на сопротивлении меж-контактного промежутка можно представить в виде:
ТТ т Жі,а , П ;
и куа ~ Ьку , + Яку‘,а,
икф — тку"
Жі
Ж
Ф
Ж
(8)
Якуі,ф.
Сопротивление межконтактного промежутка обычно считается чисто активным, а его индуктивность - близкой к нулю [2].
Дополним существующую модель (1) асинхронного двигателя с кабелем в статорной цепи, моделью коммутирующего устройства. Для этого подставим уравнение (2) и (8) в уравнение модели одиночного асинхронного двигателя. После преобразований, аналогичных [1], получим:
1 +
тк + тку
Б/,
: и а + кг
тк + тку
■( Кгіга рт//гф) ^к + Яку + Я,)-
1 +
тк + тку
Ь,
Б/,ф = и ф+ к,
Ьк + Ьку (9)
Я-і
Ггф + р®//га ■ — Я і —
)-(Як
+ Яку + Я, },ф
Ь,
Я,),
Гга р®/гф,
(В/г
Б/Ф — ~Яг^гф + рт/та
где Ьк=Ьк1+ Ьк2 и Як=Як1+ Як2.
Уравнение (9) представляет собой модель одиночного асинхронного двигателя с кабелем в статорной цепи и произвольным расположением коммутирующего органа на кабельном участке.
Рассмотрим динамические процессы, происходящие в электромагните при его включении. В общем случае уравнение баланса напряжений на обмотке электромагнита при его включении на напряжение ик( имеет вид [2]:
Ж/кі
икі — ікіЯкі +"
Ж
(10)
где к - ток, протекающий по обмотке контактора КУ; Ккг активное сопротивление обмотки контактора; цц - суммарное потокосцепление обмотки электромагнита контактора. Представим потокосцепление как цы^к^а, где Ьа - динамическое зна-
чение индуктивности электромагнита, и подставим в (10), тогда:
Ж ЖЬж
икі — ікіЯкі + ЬЖ ~ЖГ + ікі Ж, Жі Жі
(11)
Индуктивность электромагнита является величиной, зависящей от геометрических размеров магнитопровода или же его магнитного сопротивления. При включении контактора, т.е. замыкании магнитной системы, силовые линии магнитной индукции сгущаются, тем самым увеличивается потокосцепление обмотки электромагнита с магнитной системой, как следствие индуктивность увеличивается, а при отпадании якоря (отключении) происходит обратный процесс. Именно поэтому в динамических процессах коммутационных устройств вводится понятие динамической индуктивности электромагнита Ьа
Подставив (3) в (4) по оси а получим:
иа - Ьк1 ■
Жі
■- Як1і,
— ікіЯкі + Ж/к- , (12)
Жі
преобразовав которое, будем иметь:
и а -
Ьк1
,1",
Б/,
кг Ьк1 (Ягі,
га + рю/гР)
Ь,
(13)
Это уравнение описывает влияние переходных процессов в асинхронной машине с кабельным участком на динамические процессы в электромагните.
Модель не учитывает изменение положения якоря хкі в зависимости от подводимого к обмотке контактора напряжения икі.
Рассмотрим дифференциальные уравнения электромагнита [3]:
Ж¥1й — икі + ікі (Ы,хЫ )Які, Жхкі
Ж
■ — V,
(14)
м
где: V - скорость движения якоря; Рет- электромагнитная сила; Рт- сила веса якоря, в зависимости от рассматриваемой конструкции контактора может действовать в ту же сторону, что и электромагнитная сила, или в противоположную ей сторону; Рт},- сила механизма противодействующая электромагнитной; М- масса якоря.
Ток контактора к и электромагнитная сила Рет есть функции, зависящие от текущего положения якоря хк, и текущего потокосцепления ЦЦ.
В приложении 12 [3] рассматривается программа расчёта «БС» статических характеристик электромагнита методом участков. Программа написана на языке программирования РшскВаБЮ. Существующая в «БС» подпрограмма расчёта координат точек кривой зависимости потокосцеп-ления от тока «Л.Р», определяет набор статиче-
+
п
+
ж
,
+
ских характеристик, которые выводятся в файлы данных. Далее эти данные импортируются в оболочку MathCad, где и производится двумерная сплайн-интерполяции статических характеристик электромагнита. В результате последней определяются функции зависимости тока контактора к и электромагнитной силы Рет, от текущего положения якорях хк, и текущего потокосцепления Щр
Недостатком данного метода расчёта, является погрешность, которая возникает при сплайн-интерполяции характеристик. И как пишет сам автор: «кривые имеют подозрительно вычурный характер, что подтверждает желательность увеличения числа точек статических характеристик и увеличения значения максимальной МДС, для которой производится расчёт потокосцепле-ний...».
Избавиться от этого недостатка можно путём нахождения заданного значения к и Рет как функции хк от и щ в одном программном блоке, и не использовать средства MathCad. Для этого необходимо расширить подпрограмму расчёта координат точек кривой зависимости потокосцепления от тока «Би» приложения 12 [3], ещё одним общим циклом подекадного снижения МДС и нахождения заданного потокосцепления при заданном магнитном зазоре.
Общий вид функции тока обмотки к, построенный при помощи двумерной сплайн-интерполяции статических характеристик, представлен на рис. 3.
Так как динамические процессы в коммутирующем устройстве и асинхронном двигателе влияют друг на друга, то уравнения (6), (7), (9) и (14) следует рассматривать как общую систему уравнений.
Дополнительно произведём замену уравнения баланса напряжений на обмотке контактора в системе уравнений (14), уравнением (13). Дополним
Рис. 3- Функция изменения тока обмотки контактора от суммарного потокосцепления и перемещения якоря
модель асинхронного двигателя уравнением движения Бт. Получим общую систему уравнений (15) модели рассматриваемой на рис.1:
(( Т Т \
1 +
Ьк + Ьку
Ь,
Б/,а — иа + кг
Ьк + Ьку
Ь,
( (га рт/гф і,а)к + яку + Я,)' Б/га — -Ягіга - рт/гф,
Ьк + Ьку
Ь,
Б/ф — ир + кг-
( Яг^р + рт//га) і,фФЯк + яку + Я,)' Б/гф — -Ягігф + рт/га,
Бт — ■
Б/кі =
J
иа- Ьк1Б/га- кг^ * Ь, Ь,
іягіга+ рт/гфР і,аяк1 - ікі (/кі,хкі )які
Бхкі — v,
— рет -/кі, хкі ) + рт ртИ М ’
0001, если хкі < хтіп> аг^(і - ік - 0.005~) ■ 107
Бу —-
яку —
Ьку —
10
если хк — хтіп, 0.0001, если хкі < хтіп,
15)
20 ■
7
агсі§ (і - ік - 0.005) ■ 10
если хкі — хтіп.
Коэффициенты 10 и 20 в системах уравнений для активной и индуктивной составляющих межкон-тактного промежутка определяют величину сопротивления последнего.
В полученной системе уравнений первые 5 уравнений описывают динамику в асинхронной машине, последующие 3 уравнения позволяют получить динамические характеристики электромагнита, последние 2 функции Яку и Ьку в системе
уравнений являются связующим звеном между
Ток обмотки
Рис. 4 - Ток обмотки
X
X
Ж
71
с
Изменение суммарного потокосцепления обмотки контактора
Скорость движения якоря
-1 .................................................................................................У................................................................................................................................
1
1. с
Рис. 7 - Скорость движения якоря
ра
динамикой в асинхронной машине и процессами в электромагните контактора. Полученная модель позволяет рассчитать динамические характеристики модели с учётом изменения положения якоря хкі в зависимости от минимально допустимого напряжения питания обмотки контактора итіп.
По полученной математической модели (15) была реализована программная модель в среде Бе1рЫ 6. В качестве асинхронного двигателя был взят двигатель типа ДКВ355ЬВ4, мощностью 315кВт, напряжение питания ин=1140 В. В качестве коммутирующего устройства - пускатель се-
рии ПВИ-250 с контактором КТУ-4Б.
Рассмотрим режим повторного включения асинхронного двигателя. Для этого в момент времени 1=0.4с. увеличим сопротивление обмотки контактора до таких значений, при которых якорь опадёт от магнитной системы, а в момент времени 1=0.62с. - установим номинальное сопротивление обмотки контактора.
Ниже представлены полученные динамические характеристики. На рисунках 4-8 представлены характеристики контактора.
Характеристики асинхронного двигателя представлены на рисунках 9-11.
Таким образом, полученная математическая модель (система уравнений (15)) позволяет оценить динамические характеристики комплекса состоящего из асинхронного двигателя с кабелем в статорной цепи и устройством коммутации, которое имеет произвольное расположение на кабельном участке. Особое внимание было уделено работе коммутирующего устройства (пускателю), динамические характеристики которого зависят от
Скорость
Момент
1. с
Рис. 10 - Момент двигателя
двигателя
динамики асинхронного двигателя, а также от расположения на кабельном участке.
70
В.Г. Каширских, А.Г. Анисимов
СПИСОК ЛИТЕРАТУРЫ
1. Ещин Е.К. Моделирование электромеханических процессов многодвигательных электроприводов горных машин. - Кемерово: Кузбасский гос. техн. ун-т, 1999. - 115с.
2. И.С. Таев, Б.К. Буль, А.Г. Годжелло и др., Основы теории электрических аппаратов. Учеб. Для вузов. М.: Высш шк. 1987. - 352с.
3. Переходные процессы в электрических машинах и аппаратах и вопросы их проектирования: Учеб. пособие для вузов/ Гольдберг О.Д., Буль О.Б., Свириденко И.С., Хелемская С.П.; Под ред. О.Д. Гольдберга -М.: Высш. шк., 2001. 512с.: ил.
4. Программирование и применение ЭВМ в расчётах электрических аппаратов: Учеб. пособие для
вузов по спец. «Электр. аппараты» / Никитенко А.Г., Гринченко В.П., Иванченко А.Н. - М.: Высш шк.
1990. - 231с.: ил.
□ Автор статьи:
Губенков Александр Вячеславович - аспирант каф. вычислительной техники и информационных технологий
УДК 622.271.4:621.879-83:178.1
В.Г. Каширских, А.Г. Анисимов
ОЦЕНКА ПАРАМЕТРОВ ДВИГАТЕЛЯ ПОСТОЯННОГО ТОКА С ПОМОЩЬЮ МЕТОДА НАИМЕНЬШИХ КВАДРАТОВ
Знание текущих значений параметров двигателей постоянного тока (ДПТ) в процессе их работы в составе электроприводов карьерных экскаваторов и буровых станков необходимо для создания современных систем управления, диагностики и защиты.
Параметры конкретных ДПТ (активные и индуктивные сопротивления, а также постоянные времени обмоток возбуждения и якоря) могут значительно отличаться от каталожных данных, которые являются усредненными для серии двигателей, а также могут изменяться после некоторых видов ремонта. Кроме того, значения параметров зависят от режимов нагружения и теплового состояния ДПТ. Непосредственное измерение значений параметров при работе ДПТ невозможно, поэтому целесообразно использование наблюдающих устройств. В качестве исходной математической модели ДПТ используем уравнения [1]:
- ' Я Т Шв.. ив — 'в яв ^ Тв .
ш
Шя- + кФа,
и я — 'я Яя 2 + Тя 2'
Шг
(1)
где ив, 1в - напряжение и ток обмотки возбуждения; Яв, Ье - активное сопротивление и индуктивность обмотки возбуждения; и„ ¡я - напряжение и ток обмотки якоря; Яя£, Тя£ - суммарное активное сопротивление и индуктивность цепи якоря; к -конструктивный коэффициент; Ф - магнитный поток; а - частота вращения якоря.
Преобразуем уравнения (1):
1 ( ^
я
ив Тв
в \ 1(
я
я 2
и я — Тя 2 —^— кФа
и запишем их в виде:
в
'в — а1' ив — а2 ' ~Г ш
I я — аз • и я — а4 •
1
а5 -а, 1
(2)
где а, ——; а2 ——; а3 — —
1 2 3 я
я 2
Т 2 кФ
а4 — ——; а5 — •
Я я 2 Я я 2
Для применения метода наименьших квадратов представим уравнения (2) в матричной форме [2]:
%1 — Х1 - А1 + ув; 1
^ (3)
^2 — Х2 - А2 +уя,)
где
— 'в; ^2 — 'я ;
X: —
и в ; х 2 — и я я — а
_ в Ш _ я
