CC BY
91
ЭЛЕКТРОТЕХНИКА. ЭНЕРГЕТИКА ОМСКИЙ НАУЧНЫЙ ВЕСТНИК № 3 (113) 2012
*
Ю. Б. Сениченков. — СПб. : БХВ-Петербург, 2007. — 352 с. — ISBN 5-94157-580-7.
2.Brenan, K. E. Numerical solution of initial-value problems in differential-algebraic equations / K. E. Brenan, S .L. Campbell, L. R. Petzold. — Philadelphia : SIAM, 1996. — 256 p.
3.Ковалёв, Ю. З. Разработка алгоритмов исследования динамики обобщенного электромеханического преобразователя энергии на ЭЦВМ : дис. ... д-ра техн. наук : 05.09.01 / Юрий Захарович Ковалёв. — Москва, 1980. — 415 с.
4.Ковалёв, В. З. Моделирование электротехнических комплексов и систем как совокупности взаимодействующих подсистем различной физической природы : дис. ... д-ра техн. наук : 05.09.03 / В. З. Ковалёв. — Омск, 2000. — 370 с.
5.Rosenbrock, H. H. Some general implicit processes for the numerical solution of differential equations / H. H Rosenbrock // Computer J. — 1963 — vol. 5, pp. 329 — 330.
6.Артемьев, С. С. Численные методы решения задачи Коши для систем обыкновенных и стохастических дифференциальных уравнений : моногр. / С. С. Артемьев ; под ред. Г. А. Михайлова. — Новосибирск : [б. и.], 1993. — 156 c.
7.Хайрер, Э. Решение обыкновенных дифференциальных уравнений. Жесткие и дифференциально-алгебраические задачи / Э. Хайрер, Г. Ваннер : пер. с англ. — М. : Мир, 1999. — 685с. - ISBN 5-03-003117-0.
САВЧЕНКО Антон Анатольевич, ассистент кафедры энергетики Академического института прикладной энергетики, г. Нижневартовск.
Адрес для переписки: e-mail: gruzo-work@yandex.ru КОВАЛЁВ Александр Юрьевич, кандидат технических наук, директор Нижневартовского филиала Омского государственного технического университета.
Адрес для переписки: e-mail: gruzo-work@yandex.ru
Статья поступила в редакцию 04.06.2012 г.
© А. А. Савченко, А. Ю. Ковалёв
уДК 621318 М. В. СЕМЕНЯК
В. К. ФЁДОРОВ
Омский государственный технический университет
СТАТИСТИЧЕСКИЕ МОДЕЛИ НЕСИНУСОИДАЛЬНОСТИ НАПРЯЖЕНИЯ
В данной статье отражены канонические гармоники напряжения. Здесь определили коэффициенты несинусоидальности и пульсаций. Кроме того, в статье рассмотрены неканонические гармоники и энтропия.
Ключевые слова: канонические гармоники, математическое ожидание, дисперсия, неканонические гармоники, энтропия.
Тенденция к широкому внедрению управляемых преобразователей УП определяется потребностями увеличения экономической эффективности производства. В то же время проявляется отрицательное влияние применения УП на электрические параметры СЭОУ (наличие высших гармоник в кривых напряжения и тока).
Искажение синусоидальности кривой напряжения характеризуется коэффициентом гармоник
ТР1
, (1)
щ
где и, и — соответственно действующее значение у-ой и основных гармоник напряжения.
В соответствие с Г0СТ-13109-2007 на качество электроэнергии коэффициент гармоник не должен превышать 5%. Но как показали многочисленные экспериментальные исследования [1-3] в электрических сетях со значительным удельным весом статических преобразователей несинусоидальность напряжения значительно превышает допустимые 5 %.
Канонические гармоники. В установившемся режиме a = const УВП может быть представлен как некоторый генератор, дающий на входе постоянную и переменную составляющую ЭДС. Постоянная составляющая определяется выражением [4]
ЛЬ
(2)
П П
----+а
2 т
где Е^0) — максимальная выпрямленная ЭДС, равная
(3)
я т
где Ефм— амплитудное значение фазной ЭДС, т — число фаз преобразователя, ю0 — круговая частота питающего напряжения.
Переменная составляющая ЭДС Е(а), играющая роль помехи, имеет гармонические колебания, зависящие от угла а
00
^(а) = X Емк (а) +фх ) (4)
к=1
где Емк(а) — амплитуда К-й гармоники, равная
(5)
(кт) -1
Фк — фазовый сдвиг К-й гармоники относительно момента открывания вентиля
ФК = аісід
ктїда
(6)
При расчётах спектрального состава первичного линейного (переменного) и выпрямленного токов вводятся два допущения:
а)трёхфазная система линейных напряжений симметрична;
б) в цепи постоянного тока имеется сглаживающий реактор.
Результирующие коэффициенты несинусои-дальности и коэффициент пульсаций определим как
Кнс = ^(Канс)2 + (ККнс)2, Кп=^(Кап)2+(Ккп)2,
(7)
(8)
По-по
где д = 0,05-0,1 — относительная среднеквадратическая погрешность определения корреляционной функции, п0 — число пересечений функций 1^) линии математического ожидания в единицу времени.
Для исследуемых процессов п0=0,2с-1, при д, = 00,5 получаем Т=200 с. Шаг квантования, опре-
енты несинусоидальности и пульсаций, обусловленные действием канонических и нормальных гармоник напряжения.
При случайном изменении суммарной нагрузки первичный (выпрямленный) ток представляется случайным процессом. Параметры гармонических ЭДС УВП остаются при этом детерминированными величинами. Для совершенствования методик расчёта и изменения параметров несинусоидальных режимов, а также для выбора рациональных методов и средств снижения уровней гармоник в электрических сетях важное значение имеет установление закономерностей, характеризующих изменение во времени гармоник тока и напряжения и эквивалентных их значений при случайном изменении суммарной нагрузки УВП. Известно, что случайный процесс изменения гармоник тока и напряжения в промышленных электрических сетях является стационарным в широком смысле [1, 5], это позволяет анализировать стохастические процессы в рамках спектрально-корреляционной теории случайных процессов.
Ниже рассматриваются корреляционные функции и спектральные плотности указанных процессов для отдельных гармоник и их эквивалентных значений, а также соответствующие частотные характеристики и законы распределения.
Регистограммы процессов изменения отдельных гармоник тока во времени отдельно строились с помощью анализатора гармоник. При обработке регистограмм процесс І (і) задавался цифровой последовательностью с шагом квантования Ді, связанным с минимальной длительностью реализаций соотношением
Т
Л* = -
п
где п — число точек дискретной последовательности.
Значение Т определялось по выражению
Т = -
Рис. 1. Возникновение анормальных гармоник в линейном токе и выпрямленном напряжении:
1 — источник энергии, 2 — опорное напряжение,
3 — вход системы управления, 4 — обратная связь,
5 — управляющие импульсы
делённый по теореме Котельникова, Дt=0,2 с (максимальная частота в спектре огибающей для прокатных станов / =2,5+5 Гц).
м ' ^'
Число точек дискретной последовательности
п = — = 10000. М
(9)
Числовые характеристики процессов 7 (і) находились с помощью ЭВМ по следующим выражениям [6]: математическое ожидание —
к=1
у/к
(10)
дисперсия
І(Л*-яі2(/у))2 к=1_____________
(11)
где с2Ь — среднеквадратическое отклонение процесса.
Нормированная корреляционная функция
п-т к=1
(12)
Нормированная спектральная плотность п-т
(13)
л-1
где Гц.
Г
Неканонические гармоники и энтропия. В [7]
показано, что в установившемся режиме работы управляемого выпрямителя при разбросе углов управления от установочных значений в системе импльсно-фазового управленая (ИФУ) в линейном токе выпрямителя наряду с каноническими появляются неканонические (анормальные) гармоники (рис. 1).
Отклонение моментов возникновения управляющих импульсов от установочных значений зачастую имеет случайный характер из-за разброса параметров отдельных элементов системы ИФУ, колебаний напряжения и нагрузки.
Для таких случайных процессов момент появления tJK, любого п-го импульса к-ой реализации процесса может быть представлен в виде:
і м = п* Т+Е ,
(14)
к
а
ОМСКИЙ НАУЧНЫЙ ВЕСТНИК № 3 (113) 2012 ЭЛЕКТРОТЕХНИКА. ЭНЕРГЕТИКА
239
ЭЛЕКТРОТЕХНИКА. ЭНЕРГЕТИКА ОМСКИЙ НАУЧНЫЙ ВЕСТНИК № 3 (113) 2012
Чі
240
где T — длина тактового интервала; £,n — случайная величина с нулевым средним значением.
Отметим, что среднее значение tn любого n-го импульса равно установочному значению момента возникновения этого имцульса. При таких условиях работы системы ИФУ линейный ток выпрямителя можно считать стационарным случайным процессом, спектральный состав которого представляет собой ряд гармоник со случайными амплитудными и начальными фазами. Коммутация тока с одного вентиля на другой не нарушает стационарности линейного тока управляемого выпрямителя.
Представим линейный ток выпрямителя случайной стационарной функцией времени вида
о 00
т= К*™ cos лев t+isn sin nat),
л=2
где icn и in образуют множество взаимно независимых случайных величин, каждая из которых имеет нулевое среднее значение. Закон распределения коэффициентов in и in предполагается таким, что вероятность того, что ряд (15) сходится, равна единице. Покажем, что коэффициенты in и in распределены по нормальному закону.
Для того, чтобы в этом убедиться, фиксируем произвольные значения n и рассматриваем заданное множество С в плоскости i , i , которому при-
cn sn' 1 J 1
надлежит пара коэффициентов (icn, isn).
Если i(t) — стационарный случайный процесс, то произвольный сдвиг во времени на величину т должен оставить все вероятностные характеристики этого процесса без изменения. Поскольку
Приняв для удобства
получим
PAl )=h.(i2 ), PAl )=h,(l2 ),
1' cn’ 1' cn>' 2' sn’ 2' sn'1
P(i2 +i2 )=h.(l2 )*h,(l2 ).
' cn sn’ 1' cn’ 2' sn’
Причём (20) выполняется тождественно относительно i и i . Предположив i =0, найдём
сп йп 1 ^ sn ' ^
P(i"cn)=h1(i'2cn)*h2(0).
Аналогично
P(l2sn)=h1(0)*h2(i'2sn).
(21)
(22)
Таким образом, обе функции Л пропорциональны функции Р, поэтому
P(i2 +l2 )=к*(г2 )*P(l2 ),
cn sn sn sn
(23)
где к — коэффициент пропорциональности. Откуда
1пР(г'2 + г'2 ) = 1пк + 1пР(г'2 )+1пР(г'2 ) (24)
' сп йп ' сп’ ' йп’
для всех гсп и гп. Это справедливо лишь тогда, когда 1пР(г2) есть линейная функция своего аргумента, то
lnP(i2 )=R*i2 +R',
cncn
(25)
где R и R' — постоянные. Таким образом, получаем
icncosna>-(t-x) +
+isn sin лю ■ (t - х) = (icn cos лют -- zsn sin лсот) - COS ЛЮ t -- (icn sin лют+isn cos лют) • sin ЛЮ t = = i’cn cos nat + i’sn sin neat,
(16)
то действие переноса на любую величину состоит в замене коэффициентов (і , і ) новыми (і’ , і’ ),
11 і V Сп> snl ' сп' sn’1
определёнными описанным выше способом.
Из (16) видно, что точка (і’, і'п) связана с точкой (і і) поворотом около начала координат на угол пют. Следовательно, если точка (і, і) принадлежит множеству С, то точка (і’сп, і'п) принадлежит множеству Ст, представляющему собой множество С, повёрнутое на угол пют [8]
(17)
для всех т и всех С.
Очевидно, что для удовлетворения (17) совместная плотность вероятности коэффициентов (і, і) должна зависеть только от радиального расстояния от начала координат
P(l , l )=P(l2 +i2 ), (18)
cn sn cn sn
Р(ІСП'І*п) = Р{І2сп +і2зп) = г-е-К^™+‘2™\ (26)
где г — постоянная.
Из (26) видно, что ісп и і — независимые переменные и нормально распределённые косинусные и синусные коэффициенты высших гармоник.
Теперь вернёмся к (15) и представим его как
'(0= Е!п8Іп(па)( + ¥п)' (27)
Л=1
где іп — амплитуда анормальной п-й гармоники тока; у — начальная фаза этой гармоники;
i„=^i2cn+isn,
1 J
где гп и у — независимые случайные величины в силу независимости г и г .
•' сп йп
Найдём совместную плотность вероятности величин г и у .
пп
В общем случае
P(l ,w ) = l *P(l ,l ),
n n n cn sn
(28)
где для удобства берётся квадрат радиуса (г2сп + г2вп) вместо самого радиуса и где Р — функция, которую надо отыскать. Пусть Р1(гсп) и Р2(г'йп) — соответствующие плотности вероятности для гсп и по отдельности. В силу независимости этих величин
P(l , l )=Pl(i )*P,(l )
' cn sn' 1' cn’ 2' sn’
(19)
где ln - якобиан преобразования от переменных lcn, ln к переменным in, Yn. Поэтому
P(in<Wn) = r-in-e R'{2"
(29)
Это релеевский закон распределения. Закон распределения для ln
есть
2л 2
(30)
Для уп
со
Р{Чп)= ІР{Іп'Чп№п=С°ПЗІ. (31)
—00
Понятно, что п-я гармоника тока распределена по релеевскому закону; начальная фаза п-й гармоники подчинена равномерному распределению на отрезке (п,т). Поскольку доказательство проведено для произвольного значения п, то оно справедливо для любого п.
5. Литвак, В. В. Исследование несимметрии и несинусои-дальности напряжения на Западно-Сибирском металлургическом заводе / В. В. Литвак, В. В. Прокопчик, В. К. Фёдоров // Кибернетика электроэнергетических систем : сб. науч. тр.— Вып. 4. - Томск : ТГУ, 1975. - С. 35-40.
6. Румшиский, Л. З. Математическая обработка результатов эксперимента / Л. З. Румшиский. — М. : Наука, 1971. — 157 с.
7. Красовский, А. А. Изменение энтропии непрерывных динамических систем / А. А. Красовский // Техническая кибернетика. - 1963. - Вып. 5. - С. 3-11.
8. Лэннинг, Д. Случайные процессы в задачах автоматического управления / Д. Лэннинг, Д. Бэттин. — М. : ИИЛ, 1958. — 349 с.
Библиографический список
1. Жежеленко, И. В. Высшие гармоники в системах электроснабжения металлургических заводов / И. В. Жежеленко // Электричество. — 1979.- № 11. — С. 37-43.
2. Жежеленко, И. В. Показатели качества электроэнергии на промышленных предприятиях / И. В. Жежеленко. - М. : Энергия, 1977. - 126 с.
3. Исследование несинусоидальных кривых токов и напряжений в заводских системах электроснабжения : отчёт о НИР / Р. И. Борисов, В. К. Фёдоров. - Томск : ТПИ, 1976. - 68 с. -Гос. регистрация № 75060767.
4. Булгаков,А.А.Новаятеорияуправляемыхвыпрямителей/ А. А. Булгаков. - М. :Наука, 1970. - 211 с.
СЕМЕНЯК Мария Владимировна, аспирантка кафедры электроснабжения промышленных предприятий Омского государственного технического университета (ОмГТУ), ассистент кафедры электротехники и электрификации сельского хозяйства Омского государственного аграрного университета им. П. А. Столыпина.
Адрес для переписки: e-mail: mary_semenyak@mail.ru ФЁДОРОВ Владимир Кузьмич, доктор технических наук, профессор кафедры электроснабжения промышленных предприятий ОмГТУ.
Статья поступила в редакцию 25 04.2012 г.
© М. В. Семеняк, В. К. Фёдоров
УДК 621317 В. Ю. СЫСОЛЯТИН
Омский государственный технический университет
ЦИФРОВОЕ
УСТРОЙСТВО ЗАРЯДА-РАЗРЯДА ХИМИЧЕСКИХ ИСТОЧНИКОВ ТОКА
Рассматривается цифровое устройство, предназначенное для автоматизированного заряда аккумуляторов в разнополярном импульсном режиме в формате: за-ряд—пауза—разряд—пауза. Данное устройство позволяет устанавливать различные интервалы времени длительности каждой из фаз с одновременной индикацией их значений на жидкокристаллическом дисплее. В статье рассматриваются функциональная и принципиальная электрические схемы устройства, а также алгоритм функционирования программы.
Ключевые слова: заряд-разряд химического источника тока, блок-схема программы, принципиальная электрическая схема, микропроцессорное устройство.
В процессе эксплуатации различных ХИТ (химических источников тока) важное значение имеет технология заряда кислотных аккумуляторов. В настоящее время в научных публикациях существует большой объем противоречивой информации, описывающей использование различных режимов заряда аккумуляторных батарей, а также их влияние на изменение характеристик аккумуляторов.
Кроме того, существующие простейшие технологии не учитывают всех особенностей физико-химических процессов, протекающих в аккумуляторе, и поэтому не являются оптимальными. Оптимальные технологии заряда достаточно трудоёмки и без
автоматизации трудновыполнимы, поскольку у рядового пользователя зачастую нет времени, чтобы методично следить за многочасовыми процессами заряда и выполнять все процедуры.
Целью данной работы является:
—разработка структурной и принципиальной схемы цифрового устройства для широтно-импульсного заряда-разряда аккумуляторов;
—разработка алгоритма управляющей программы, обеспечивающей возможность изменения временных интервалов режима заряда и разряда в широких пределах, с целью изучения данного процесса.
ОМСКИЙ НАУЧНЫЙ ВЕСТНИК № 3 (113) 2012 ЭЛЕКТРОТЕХНИКА. ЭНЕРГЕТИКА
