CC BY
52ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2021. Т. 29. № 4
doi: 10.14498/tech.2021.4.8
УДК 621.311.1
МЕТОД РАСЧЕТА КОЭФФИЦИЕНТА ИСКАЖЕНИЯ СИНУСОИДАЛЬНОСТИ КРИВОЙ НАПРЯЖЕНИЯ, СОЗДАВАЕМОГО ТРЕХФАЗНЫМИ ВЫПРЯМИТЕЛЯМИ
В.С. Осипов
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
E-mail: osipov_4343@mail.ru
Аннотация. Существующие в различных источниках методики расчета представляют выпрямительные установки как источники высших гармоник. При этом действующее значение тока первичной обмотки трансформатора в виде прямоугольников разлагают в ряд Фурье и получают гармонические составляющие тока, кроме основной гармоники высшие 5, 7, 11, 13, 17, 19 и т. д. Амплитуды тока высших гармоник умножаются на соответствующие частотам индуктивные сопротивления питающей сети, результаты возводят в квадрат и суммируют. Очевидно, что все слагаемые — равной величины, так как при увеличении номера гармоники ее амплитуда уменьшается в n раз по сравнению с первой гармоникой, при этом во столько же раз увеличивается частота и индуктивное сопротивление. Коэффициент несинусоидальности Ku определяется как отношение корня квадратного из суммы к величине фазного напряжения. Недостатком определения Ku является некоторая произвольность в количестве учитываемых гармоник. Если принять, например, 9 гармоник, то получим Ku = 7,22 %; если 4 гармоники, то Ku = 4,81 %. При протекании тока первичной обмотки трансформатора амплитуда напряжения на шинах подстанции уменьшается всего на 1,9 В. При этом коэффициент несинусоидальности не может быть равен 7,22 %. Во всех литературных источниках при расчете гармоник тока не дается рекомендаций, на какие номера гармоник необходимо устанавливать на подстанции резонансные фильтры. Не случайно в новом ГОСТ—2013 вычисляют значение коэффициента искажения синусоидальности кривой напряжения K(U)i в процентах как результат i-го наблюдения по формуле (то есть не производится вычисление по гармоникам тока). В данной работе производится другой подход к определению коэффициента искажения синусоидальности кривой напряжения. Выпрямительные установки не являются источниками гармоник, а представляют собой электроприемники с нелинейной характеристикой потребления электрического тока, при этом искажается форма кривой синусоиды питающего напряжения. Эта искаженная синусоида и является источником высших гармоник, которая в данной работе разлагается в ряд Фурье, при этом вычисление интегральных функций коэффициентов ряда Фурье производится в программе Mathcad. Решение произведено для трех схем выпрямителей. Определены гармоники для установки резонансных фильтров. В результате предлагается разработанная методика расчета высших гармоник напряжения и определения коэффициента искажения синусоидальности кривой напряжения на основании разложения в
Осипов Вячеслав Семёнович, доцент кафедры «Электроснабжение промышленных предприятий», кандидат технических наук.
ряд Фурье искаженной кривой синусоиды напряжения при работе трехфазных неуправляемых выпрямителей.
Ключевые слова: выпрямитель, коэффициент, несинусоидальность, ток, напряжение, высшие гармоники, ряд Фурье.
Разработка методики расчета высших гармоник и определения коэффициента искажения синусоидальности кривой напряжения является актуальной проблемой, так как существующие методики позволяют решить эту задачу с большой погрешностью. В цехах современных промышленных предприятий более половины электроэнергии используется в преобразованном виде (на металлургических заводах - более 90 %). Широко применяются вентильные преобразователи (ВП), различного рода преобразователи частоты [1, 2]. Такие нагрузки называют нелинейными, к их числу относятся в первую очередь различного рода ВП.
Эти нагрузки потребляют из сети ток, кривая которого оказывается несинусоидальной, в результате возникают нелинейные искажения кривой напряжения сети или, другими словами, несинусоидальные режимы.
Несинусоидальность напряжения и тока обусловливает дополнительные потери и нагрев, а также ускоренное старение изоляции электрооборудования [3], что подтверждается результатами диагностического обследования различного электроэнергетического оборудования: трансформаторов, кабельных линий, электродвигателей, компенсирующих конденсаторов. Кроме того, ВГ в СЭС промышленных предприятий вызывают нарушение работы и ложные срабатывания устройств релейной защиты и автоматики, приводят к сбоям в работе электронных систем управления и вычислительной техники, создают помехи в аппаратуре телемеханики и связи, искажают показания счетчиков электрической энергии. Прогрессирующее внедрение вентильных преобразователей создает актуальность решения проблемы высших гармоник в электрических сетях путем совершенствования расчета и создания устройств противодействия гармоникам. Коэффициент несинусоидальности является показателем качества электроэнергии.
Для уменьшения влияния высших гармоник применяют подключение преобразовательных устройств через фильтры высших гармоник [3-5] - это сетевые дроссели, линеаторы. Существует несколько схем, при помощи которых можно реализовать активную коррекцию коэффициента мощности, что приводит к уменьшению воздействия высших гармоник. Также применяются пассивные узкополосные и широкополосные фильтры [6-8], которые устанавливаются на шинах ТП. Во всех случаях необходимо знать амплитуду и номера гармоник.
Рассмотрим существующие методики расчета на конкретном примере схемы включения выпрямителя (расчетная схема на рис. 1).
Мощность короткого замыкания системы Бк = 300 мВА. Определяется сопротивление системы, приведенное к напряжению 0,4 кВ.
Индуктивное сопротивление
Сопротивление кабельной линии ААБ (3Х70+1Х35), ХУд = 0,0612 Ом/км, ЯУд = 0,447 Ом/км, длина кабеля 1к = 800 м.
Сопротивления кабеля, приведенные к напряжению 0,4кВ:
300000000
4002
0,00053 Ом.
ЯК = ДУД • 1К-104) = 0,447• 0,8= 0,000572 Ом, = ХУД • 1К •|04] = 0,0612 • 0,8 •('00'] = 0,0000783 Ом.
Рис. 1. Расчетная схема
Сопротивления трансформатора Т1 подстанции мощностью 630 кВА, приведенные к 0,4 кВ. Технические данные: ДРхх = 1,16 кВт, ДРКЗ = 7,6 кВт, Дикз = 5,5 %.
Полное сопротивление, приведенное к 0,4 кВ:
= = М = 0,014 Ом.
т 100 5н 100 630000
Активное сопротивление
= 7600 • 40°2 = 0,00306 Ом,
6300002
=^12Т2 - П- =у10,0142 - 0,003062 = 0,01366 Ом. Суммарные сопротивления питающей сети:
ДПС = 0,000572 + 0,00306 = 0,00363 Ом, ХПС = 0,00053 + 0,0000783 + 0,01366 = 0,0143 Ом.
Трехфазная мостовая схема Ларионова
Расчет параметров трансформатора выпрямительной установки, соединение обмоток «звезда - звезда» (рис. 2).
А В
Рис. 2. Трансформатор
«звезда - звезда» и выпрямитель по схеме Ларионова
Мощность трансформатора Т2 (ТМ100/10) 100 кВа. Технические данные:
ДРХХ = 0,33 кВт, ДРКЗ = 1,97 кВт, Дикз = 4,5 %. Полное сопротивление трансформатора, приведенное к 0,4 кВ:
2 _ ик% ин _ 4,5 4002
100 5н 100 100000
_ 0,072 Ом.
Активное сопротивление
о АРКЗЩ
кт _ -
Б 2
1970•4002 1000002
_ 0,0315 Ом.
Индуктивное сопротивление
2Г2 - В2 0,0722 - 0,03152 _ 0,06474 Ом. Постоянное выпрямленное напряжение иа _ 48 В. Выпрямленный ток 1Л _ 1500 А. Мощность выпрямителя Рв _ 72 кВт. Фазное напряжение вторичной обмотки трансформатора для выпрямителя
[9, 10]:
С
U —-Ud
48
2,34 2,34
= 20,51 В.
Коэффициент трансформации:
ту - _ U10
kt — -
U2,
220 20,51
—10.726.
Ток нагрузки из-за наличия в ней значительной индуктивности сглажен:
I -
' к,'
Действующее значение тока вторичной обмотки трансформатора [9 - 11]:
Id 2 —
( 2
2%
2 jI>
— J-Id — 0.817Id — 0,817-1500 —1225.5 А.
Действующее значение фазного тока первичной обмотки трансформатора [9] при соединении ее в «звезду»:
Id, — !*- — 12255 — 114,26 А. KT 10,726
На рис. 3, а приведена диаграмма трехфазных напряжений. На рис. 3, b -диаграмма токов первичной обмотки трансформатора Т2 фазы А. Эти токи в виде периодической функции могут быть представлены в виде ряда Фурье [14-17]:
a м
f (x) —— + I(an cosnx + bn sinnx), (1)
2 n—1
где коэффициенты ряда a0, an, bn - действительные числа. Коэффициенты ряда определяются по формулам:
1 %
ac — Г j f (x) dx,
%
1 %
an — — jf(x)- cosnxdx,
1 %
bn — — j f (x) - sin nxdx ,
(2)
(3)
(4)
п -1,2,3,... - номер гармоник.
Если функция нечетная (симметричная относительно абсцисс), то достаточно определить только коэффициент Ъп; если функция четная, то определяются коэффициенты а0, ап; если функция аморфная (не четная и не нечетная), то определяются все коэффициенты.
Если - нечетная периодическая функция с периодом 2п, удовлетворяющая условиям разложимости в ряд Фурье, то можно определять
2 %
bn — — j f ( x) - sin nxdx
и разложение в ряд Фурье для нечетной функции имеет вид
f (x) — I bn si
sin nx.
(5)
(6)
1
m=1
Если разложить функцию тока !лх первичной обмотки трансформатора в ряд Фурье [9, 18], то получим кроме первой, т. е. основной гармоники, высшие 5, 7, 11, 13, 17, 19 и т. д. гармонические составляющие (п = 6k ± 1, где k - любое натуральное число, 6 - число пульсов схемы выпрямления за период).
ГОСТ-97 нормирует коэффициенты п-ой гармонической составляющей [1, 2, 19]
К
ип
и (п)
и
•100.
НОМ
(7)
1<И
В)
/ЯЬ-/,
п 2п ^
Рис. 3. Диаграммы напряжений и токов
Коэффициент несинусоидальности напряжения питающей сети определяется по формуле [1, 2, 19]
Ки =
1
40
I и2п -100
п=2
и
(8)
НОМ
где п - номер гармоники; ип - падение напряжения в сопротивлении питающей сети от тока соответствующей гармоники.
В таблице 1 приводятся допустимые значения коэффициентов несинусоидальности.
Таблица 1
Значения коэффициентов гармонических составляющих напряжения Кип %
ГОСТ 32144—2013
Нечетных не кратных т рем Нечетных кратных трем Четных гармонических
Порядок гармониче- % Порядок гармониче- % Порядок гармониче- %
скои составляющей п ской составляющей п ской составляющей п
5 6 3 5 2 2
7 5 9 1,5 4 1
11 3,5 15 0,3 6 0,5
13 3 21 0,2 8 0,5
17 2 >21 0,2 10 0,5
19 1,5 12 0,2
23 1,5 >12 0.2
25 1,5
>25 1,5
В [19] приведен расчет коэффициента искажения синусоидальности кривой напряжения с учетом 5, 7, 11, 13 гармоник:
К = 100% -УЕип2 = 100% - У7,742 + 7,742 + 7,742 + 7,742 = 7 04% (9)
и = иф = 220 = , °, (
что меньше допустимых 8 %, определяемых ГОСТ-97 для напряжения 380 В.
В выражении (9) все слагаемые под корнем равны, так как при увеличении номера гармоники ее амплитуда уменьшается в п раз по сравнению с первой гармоникой, при этом во столько же раз увеличивается частота и индуктивное сопротивление ХПС = Ьш питающей сети. Однако в [19] наблюдаются ошибки, так как в (9) под корнем получены числа путем перемножения сопротивления ХПС на амплитуду тока гармоник. Целесообразнее умножать на действующее значение тока гармоники. Кроме того, для первой гармоники принято действующее, а не амплитудное значение тока первичной обмотки трансформатора.
Для правильного решения необходимо разложение в ряд Фурье функции тока (рис. 3, в). Разложение в ряд Фурье целесообразно выполнять с применением вычислительной техники, например в программе Mathcad [14].
На панели Mathcad размещаются: калькулятор, греческий алфавит, программирование и символьные операции. Затем записывается интеграл по выражению (5), нажимается ^^ (символьные), в следующей строке записывается п = 1 (первая гармоника), затем копируется первое выражение и нажимается (=) на калькуляторе - получаем решение при (х = Ш).
5 ,,„( ( пп \ ( 5пп -п 2281 cos I — I- cos I-
26г_„ . , , , > ^ I 6 1 ^ 6
| (114) sin (пх) dx
п 1 пп
п: = 1
1
—п
5
—п
2б[(П4)sin(пх)dx ^ = 125,703.
тг ~ тг
Таким образом, отношение амплитуды первой гармоники к действующему значению 125,7/114=1,103. В результате получим формулу для определения коэффициента несинусоидальности с учетом действующих значений токов:
к = 100%-1,103 • 1,1 • Хпс •ТУ (10)
" л/2 • иф '
где V - количество учитываемых гармоник.
Подставляя значения в (10), получим Ки = 7,22% при V = 9 и Ки = 4,81 % при V = 4 .
Для проверки получим уменьшение амплитуды фазного напряжения:
Диф = 311,1-0,0722 = 22,46 В.
В соответствии с выражением (1) после суммирования гармоник должны получить функцию, соответствующую рис. 3, в, для чего необходимо учитывать не менее 12 гармоник.
Недостатком определения Ки является некоторая произвольность в количестве учитываемых гармоник. А также во всех источниках при расчете гармоник тока не дается рекомендаций, на какие номера гармоник необходимо устанавливать на подстанции резонансные фильтры.
Не случайно в новом ГОСТ [20] вычисляют значение коэффициента искажения синусоидальности кривой напряжения К^ в процентах как результат 7-го наблюдения по формуле (то есть не производится вычисление по гармоникам тока):
Ки, =
V
40
п=2 -100,
и (1)7
где Вц — действующее значение междуфазного (фазного) напряжения основной частоты для 7-го наблюдения, В.
В данной работе производится другой подход к определению коэффициента искажения синусоидальности кривой напряжения. Выпрямительные установки не являются источниками гармоник, а представляют собой электроприемники с нелинейной характеристикой потребления электрического тока, при этом искажается форма кривой синусоиды питающего напряжения. Эта искаженная синусоида и является источником высших гармоник.
Полное сопротивления питающей сети:
Xпс =Т^ПС^ПС = \/0,003632 + 0,01432 = 0,0148 Ом.
Амплитуда фазного напряжения:
ит1 =л/2 • 220 = 311,1 В.
Величина амплитуды напряжения при протекании тока :
ит2 = ит1 -1,1 • Xпс = 311,1 -114,26 • 0,0148 = 309,4 В.
На рис. 4 показана расчетная графическая функция искаженной формы синусоиды. Это нечетная периодическая функция с периодом 2п, удовлетворяющая условиям разложимости в ряд Фурье по выражению (5). При этом определяются
коэффициенты разложения на отрезках Ъ1п (0-п/6), Ъ2п (п/6-5/6 п), Ъ3п (5/6 п- п), затем определяются значения суммы Ьп= Ьы+ Ъ2п+Ъ3п.
Рис. 4. Расчетная графическая функция
Решение интегральных функций производится в программе Mathcad. Например, решение для Ъ для гармоники п = 5 :
1
—л
6
2 Г
Ъ{: = — I 311.1 • sin ( х)• sin (пх)dx ^
7Г *
311.1-sin (-0.52п + 0.52 ) 311.1 • sin ( 0.52п + 0.52)
2п - 2
2п - 2
п
п := 5.
26
Ъ1 := — I 311.1- sin ( х)• sin (пх) dx ^ 21.44.
тг *
Аналогичные решения - для Ъ2, Ъ3 при варьировании номера гармоник, результаты представлены в табл. 2.
Таблица 2
Коэффициенты ряда Фурье для двенадцати гармоник
2
>
п
№ гарм. 1 2 3 4 5 6 7 8 9 10 11 12
Ъ1п 8,97 16,5 21,44 23,1 21,44 16,98 10,72 3,93 -2,14 -6,5 -8,58 -8,31
Ъ2п 291,56 0 -42,65 0 -42,65 0 -21,32 0 4,26 0 17,06 0
Ъ3п 8,97 -16,5 21,44 -23,1 21,44 -16,98 10,72 -3,93 -2,14 6,5 -8,58 8,31
Ъп 309,9 0 0,234 0 0,234 0 0,12 0 -0,02 0 0,1 0
Наибольшие гармоники - с номерами 3, 5, 7, 11, на эти гармоники можно устанавливать резонансные фильтры. Полученные результаты существенно отличаются от расчетов с использованием гармоник тока.
Коэффициент несинусоидальности в соответствии с выражением (8)
К =У °.2342 + °.2342 + °.121±0.11 • ,00% = 0.166%.
м 220 Коэффициент несинусоидальности предлагается определять как отношение суммы действующих значений гармоник к фазному напряжению
= (0,234 + 0,234 + 0.12 + 0.1)-1.1, ,00% = 0 22%.
42 - 220
Соединение первичной обмотки трансформатора в треугольник (рис. 5) при тех же значениях параметров выпрямителя.
Действующие значения токов фаз первичной обмотки [13] равны:
1й 2
1225.5
1аЬ 1са ^ Т3КТ л/3 -10.728
= 66.03А.
так как у обмотки увеличивается число витков в -у/3 раза при напряжении 380 В.
/в Е
А
/аЬ 3 /Ьс
у Лч Лу
Рис. 5. Первичная обмотка в треугольник
На рис. 6 показана векторная диаграмма для определения тока из сети, например для фазы А:
^ = ^ - ~Г
1Д 1 аЬ са '
Результирующий вектор 1д сдвинут по фазе на 300 относительно 1аЬ. Его величина определяется с учетом фазового сдвига:
!йх = /д = 1аЬ - ^300 + 1сад - ^300 = 2 - 66.03 - ^300 = 114.36А ,
то есть из сети потребляется такой же ток по величине, как при трансформаторе с соединением «звезда - звезда».
'аЬ
Рис. 6. Векторная диаграмма
/
/
с
А
В
С
а
/
А
С
При суммировании токов на графике (рис. 3, с, d, е) если прямоугольные формы не совпадают по фазе, то сила тока будет 7^/2 = 7А/2. Это подтверждается [12] осциллограммами.
Расчетная графическая функция аналогична по форме (рис. 4) с изменением пределов интегрирования: (0-1/3л), (1/3я-2/3я), (2/3п-п), при этом также изменяются амплитуды напряжений.
Величина амплитуды напряжения в пределах (0-1/3п) и (2/3п-п) при протекании тока 7Л:
ит2 =ит1 - 7,1 • 2 пс =311,1 -114,26 • 0,0148 = 309,4В.
Величина амплитуды напряжения в пределах (1/3л-2/3л) при протекании тока 7,1/2:
ипЪ =ипХ -^• 2пс =311,1 -• 0,0148 = 310,25В .
Аналогичный расчет, как и в предыдущем случае, показывает, что номера гармоник будут такими же (см. табл. 2), амплитуда гармоник будет в два раза меньше.
Трехфазная нулевая схема
Фазное напряжение вторичной обмотки трансформатора для выпрямителя [9, 10] (рис. 7):
48
и = и'
1,17 1,17
= 41,02В.
Коэффициент трансформации:
КТ = ию 220
и2ф " 41,02
А В С
= 5,36.
V 3? 3?
Рис. 7. Трехфазная нулевая схема
Действующее значение вторичного тока для трехфазной нулевой схемы с учетом диаграммы (рис. 3) [10, 11]:
1,2 =
2л
(2 ^ —я
0
V
= = 867,05А. л/3
Действующее значение первичного тока
= ^ ^,/2:1500 = .
,10 3КТ 3 • 5,36
Трансформатор Т3 для трехфазной нулевой схемы должен быть мощнее в 1,3 раза, поэтому принят трансформатор Т3 (ТМ160/10). Мощность трансформатора 160 кВа. Технические данные:
ДРХХ=0,51 кВт, ДРКЗ=2,65 кВт, ДиКЗ=4,5 %. Полное сопротивление трансформатора, приведенное к 0,4 кВ:
2т = - Ш = = 0,045Ом.
т 100 5Н 100 160000
Активное сопротивление:
п АРкзи1 2650 • 4002
п =—=- 0,01660м .
Б2Н 1600002
Индуктивное сопротивление:
=у]2Т - Щ =,]0,0452 - 0,01662 = 0,04180м. Полное сопротивление питающей сети:
2ПС =у]ЯПС + ХПС = ^0,003632 + 0,01432 = 0,01480м. Амплитуда фазного напряжения:
ит1 = >/2 • 220 = 311,1В. Величина амплитуды напряжения при протекании тока /,10:
и 20 = и 1 -/,1 • 2ПС = 311,1 -131.52• 0,0148 = 309,2В.
т20 т1 а1 ПС ' ' '
На рис. 8 показана форма искажённой синусоиды, полученная расчётным путём с учетом всех коэффициентов Фурье (2) - (4). Расчет производится в программе Mathcad, как в предыдущем случае.
ит20
/¿Г V
/
!И V
/ 1 / 1 2п
п/6 п \ / юt
Рис. 8. Расчетная графическая функция 3-фазной нулевой схемы
Определение коэффициентов разложения а0, Ь, а производится на отрезках с учетом амплитуды напряжений: (0-п/6), (п/6-5/6 п), (5/6 п- п), (п -2п), затем определяются значения их суммы. Например, для а10:
1 г
:= - f 311.1-sin(x)dx ^ 13.27 .
7Г »
Для других отрезков получено а20 =140,47, а30 = 13,27, а40 =-198,05.
Сумма значений показывает величину постоянного смещения ^ а0 =1,04 .
Определение коэффициентов производится аналогично, например для 1-й гармоники
i
—п
6
1 г
b1: = — I 311.1 - sin ( x)- sin (nx)dx ^ п o
311.1-sin ( 0.52n - Ü.52) 311.1-sin (0.52n + 0.52)
_2n - 2_2n + 2_
п
n := 1,
1
—п
6
1 Г
b1:= - I 311.1 • sin (x)• sin (nx)dx ^ 4.49. n o
Данные расчетов сводятся в таблицу (табл. 3).
Коэффициенты b синусных гармоник
Таблица 3
a10
>
№ гарм. 1 2 3 4 5 6 7 8 9 10 11 12
b1 4,49 8,25 10,72 11,55 10,72 8,49 5,36 1,96 -1,07 -3,25 -4,29 -4,15
b2 145,68 0 -21,31 0 -21,31 0 -10,65 0 2,13 0 8,52 0
bs 4,49 -8,25 10,72 -11,55 10,72 -8,49 5,36 -1,96 -1,07 3,25 -4,29 4,15
b4 155,55 0 0 0 0 0 0 0 0 0 0 0
Sb 310,16 0 0,13 0 0,13 0 0,07 0 -0,01 0 -0,06 0
Аналогично рассчитываются коэффициенты а (табл. 4):
1
—п
16г
a1: = — I 311.1 - sin (x) - cos (nx)dx ^
п 0
5388.41 - cos(0.52n) + 311.1 - n - sin (0.52n) - 6222
п(20n2 - 20) n := 1,
16
a1 := — I 311.1 - sin (x)-cos (nx) dx ^ 12.38.
1Г
Коэффициенты а косинусных гармоник
Таблица 4
№ гарм. 1 2 3 4 5 6 7 8 9 10 11 12
а1 12,38 9,87 6,19 1,97 -2,06 -5,28 -7,22 -7,7 -6,81 -4,9 -2,48 -0,09
а2 0 -85,24 0 -17,05 0 4,87 0 12,18 0 7,75 0 -1,19
а3 -12,38 9,87 -6,19 1,97 2,06 -5,28 7,22 -7,7 6,81 -4,9 2,48 -0,09
а4 0 66,02 0 13,2 0 5,66 0 3,14 0 2 0 1,38
Sa 0 0,52 0 0,09 0 -0,03 0 -0,08 0 -0,05 0 0,01
В трехфазной нулевой схеме спектр наибольших гармоник дополняется четными гармониками с номерами 2, 4, 6, 8, 10.
Определяется коэффициент несинусоидальности:
K,, =
(0,13 + 0,13 + 0,07 + 0,06 + 0,52 + 0,09 + 0,08 + 0,05) -1,1
100% = 0.4%.
л/2 • 220
Для проверки получим уменьшение амплитуды фазного напряжения ДиФ=311,1'0,004=1,24 В, что значительно меньше значений полученных из выражения (10) равных ДиФ =22,46 В. Не может быть коэффициент несинусоидальности 7,22 %, когда при протекании тока первичной обмотки трансформатора амплитуда напряжения уменьшается всего на 1,9 В.
Таким образом, предлагается разработанная методика расчета высших гармоник напряжения и определения коэффициента искажения синусоидальности кривой напряжения на основании разложения в ряд Фурье искаженной кривой синусоиды напряжения при работе трехфазных неуправляемых выпрямителей.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
10.
Жежеленко И.В., Шидловский А.К., Пивняк Г.Г. и др. Электромагнитная совместимость потребителей: моногр. - М.: Машиностроение, 2012. - 351 с.
Жежеленко И.В. Высшие гармоники в системах электроснабжения промпредприятий. -6-е изд., перераб. и доп. - М: Энергоатомиздат, 2010. - 375 с., ил.
Артюхов И.И., Сошинов А.Г., Бочкарева И.И. Электромагнитная совместимость и качество электроэнергии: учеб. пособие. - Волгоград: ИУНЛ ВолгГТУ, 2015. - 124 с. Будко С.И., Малышев С.А., Осипов В.С. Плавное регулирование температуры нагрева печи сопротивления с применением фильтра lineator для подавления высших гармоник // Сборник научных статей по материалам V Международной научно-практической конференции. 11 мая 2021 г., г. Уфа. - Уфа: Изд. НИЦ Вестник науки, 2021. - 236 с. Янченко С.А., Гужов С.В. Работоспособность и качество функционирования электротехнических комплексов и систем в режимах несинусоидальности напряжения: учеб. пособие. - М.: Изд-во МЭИ, 2016. - 44 с.
Зимин Р.Ю. Повышение качества электроэнергии в электротехнических комплексах предприятий нефтедобычи гибридными фильтрокомпенсирующими устройствами: Дис. ... канд. техн. наук. - СПб: Санкт-Петербургский горный университет, 2020. Куско А., Томпсон М. Качество энергии в электрических сетях / пер. с англ. А.Н. Рабод-зян. - М.: Додэка-ХХ1, 2008. - 336 с.: ил.
Козлов В.Г. Электромагнитная совместимость РЭС: Учеб. пособие. - Томск: Томский гос. ун-т систем управления и радиоэлектроники, 2009.
Васильев Б.В., Аристов Е.В., Хузин Р.А. Исследование характеристик трехфазного мостового тиристорного преобразователя: Метод. руководство по курсу «Элементы систем автоматики» Пермь: Пермский гос. техн. ун-т, 2006.
Бар В.И. Основы преобразовательной техники: Курс лекций. - Тольятти: Тольяттинский гос. ун-т, 2002.
1
2
3
4
5
6
7
8
9
11. Гельман М.В., Дудкин М.М., Преображенский К.А. Преобразовательная техника: учеб. пособие. - Челябинск: Изд. центр ЮУрГУ, 2009. - 425 с.
12. Янченко С.А., Гужов С.В. Работоспособность и качество функционирования электротехнических комплексов и систем в режимах несинусоидальности напряжения. - М.: Изд-во МЭИ, 2016.
13. Нестеров С.А. Шестифазные схемы выпрямления // Электроника и информационные технологии. - Саранск: Мордовский гос. ун-т им. Н.П. Огарева, 2009. - № 2 (7).
14. Крючкова И.В. Ряды и преобразование Фурье: метод. указания. - Оренбург: Оренбургский гос. ун-т, 2011. - 59 с.
15. Егорова Г.Ф., Павлова Г.А., Мазуренко И.А. Ряды Фурье: практикум по математическому анализу. - Самара: Самар. гос. техн. ун-т, 2011. - 46 с.
16. Романова Л.Д., Шаркунова Т.А., Елисеева Т.В. Интегральные преобразования: учеб. пособие. - Пенза: Изд-во ПГУ, 2015. - 80 с.
17. Лоссов К.И. Ряды Фурье: учеб. пособие. - M.: МИИГАиК, 2014. - 41 с.
18. Копытин И.И. Методика расчета показателей качества электрической энергии. - Курган: Курганский гос. ун-т, 2011.
19 ГОСТ 13109-97. Электрическая энергия. Совместимость технических средств электромагнитная. Нормы качества электрической энергии в системах электроснабжения общего назначения.
20 ГОСТ 32144-2013. Электрическая энергия. Совместимость технических средств электромагнитная. Нормы качества электрической энергии в системах электроснабжения общего назначения. - М.: Стандартинформ, 2014.
Статья поступила в редакцию 26 июля 2021 г.
METHOD OF CALCULATING THE DISTORTION COEFFICIENT SINUSOIDALITY OF THE VOLTAGE CURVE CREATED BY THREE-PHASE STRAIGHT LINES
V.S. Osipov*
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation E-mail: osipov_4343@mail.ru
Abstract. The existing calculation methods in various sources represent rectifier installations as sources of higher harmonics. In this case, the current value of the primary winding of the transformer in the form of rectangles is decomposed into a Fourier series and harmonic components of the current are obtained, except for the main harmonic higher 5, 7, 11, 13, 17, 19 The amplitudes of the current of the higher harmonics are multiplied by the inductive resistances of the supply network corresponding to the frequencies, the results are squared and summed up. Obviously, all terms are of equal magnitude, since with an increase in the harmonic number, its amplitude decreases by n times compared to the first harmonic, while the frequency and inductive resistance increase by the same number of times. The non-sinusoidal coefficient Ku is defined as the ratio of the square root of the sum to the magnitude of the phase voltage. The disadvantage of determining Ku is some arbitrariness in the number of harmonics taken into account. If we take, for example, 9 harmonics, we get Ku = 7.22%, if 4 harmonics, then Ku = 4.81%. When the current of the primary winding of the transformer flows, the voltage amplitude on the substation tires decreases by only 1.9 V. At the same time, the coefficient of non-sinusoidality cannot be equal to 7.22%%. In all literature sources, when calculating current harmonics, there are no recommendations for which harmonic numbers should be installed at the substation resonant filters. It is no coincidence that in the new GOST-2013, the value of the distortion coefficient of the sinusoidal voltage curve K(U is calculated as a percentage as a result of the i-th observation according to the formula, (that is, no calculation is made by current harmonics). In this paper, another approach is made to determine the distortion coefficient of the sinusoidal voltage curve. Rectifier installations are not sources of harmonics, but are electrical receivers with a nonlinear characteristic of electric current consumption, while the shape of the sinusoid curve of the supply voltage is distorted. This distorted sine wave is the source of the higher harmonics, which in this paper is decomposed into a Fourier series, while the calculation of the integral functions of the coefficients of the Fourier series is performed in the Mathcadprogram. The solution is made for three rectifier circuits. Harmonics for the installation of resonant filters are determined. As a result, a developed method is proposed for calculating the higher voltage harmonics and determining the distortion coefficient of the sinusoidal voltage curve based on the Fourier series expansion of the distorted voltage sine curve during operation of three-phase uncontrolled rectifiers.
Keywords: rectifier, coefficient, non-sinusoidal, current, voltage, higher harmonics, Fourier series.
Vyacheslav S. Osipov (Ph. D. (Techn.), Associate Professor 114
REFERENCES
1. Zhezhelenko I.V., Shidlovsky A.K., Pivnyak G.G., etc. Electromagnetic compatibility of consumers: monogr. M.: Mashinostroenie, 2012. 351 p.
2. ZhezhelenkoI.V. Higher harmonics in power supply systems of industrial enterprises 6th ed., reprint. and add. M: Energoatomizdat, 2010. 375 p.
3. Artyukhov I.I., Soshinov A.G., Bochkareva I.I. Electromagnetic compatibility and power quality: textbook. stipend. Volgograd: IUNL VolgSTU, 2015. 124 p.
4. Budko S.I., Malyshev S.A., Osipov V.S. Smooth regulation of the heating temperature of a resistance furnace using a lineator filter to suppress higher harmonics // Collection of scientific articles based on the materials of the V International Scientific and Practical Conference. May 11, 2021, Ufa. Ufa: Ed. SIC Bulletin of Science, 2021. 236 p.
5. Yanchenko S.A., Guzhov S.V. Efficiency and quality of functioning of electrical complexes and systems in modes of non-sinusoidal voltage: textbook. manual. M.: Publishing House of the MEI, 2016. 44 p.
6. Zimin R.Yu. Improving the quality of electricity in electrical complexes of oil production enterprises with hybrid filter-compensating devices: Dis. ... candidate of Technical Sciences. St. Petersburg: St. Petersburg Mining University, 2020.
7. Cusco A., Thompson M. Energy quality in electric networks / translated from English by A.N. Rabodzyan. M.: Dodeka-XXI, 2008. 336 p.
8. Kozlov V.G. Electromagnetic compatibility of RES: Textbook. Tomsk: Tomsk State University of Control Systems and Radio Electronics, 2009.
9. Vasiliev B.V., Aristov E.V., Khuzin R.A. Investigation of the characteristics of a three-phase bridge thyristor converter: Method. course guide "Elements of automation systems". Perm: Perm State Technical University. 2006.
10. Bar V.I. Fundamentals of transformative technology: A course of lectures. Tolyatti: Togliatti State University, 2002.
11. Gelman M.V., Dudkin M.M., Preobrazhensky K.A. Transformative technique: textbook. stipend. Chelyabinsk: Publishing House of the SUSU Center, 2009. 425 p.
12. Yanchenko S.A., Guzhov S.V. Operability and quality of functioning of electrical complexes and systems in non-sinusoidal voltage modes. M.: Publishing House of MEI, 2016.
13. Nesterov S.A. Six-phase rectification schemes // Electronics and information technology. Saransk: Mordovian State University named after N.P. Ogarev, 2009. № 2 (7).
14. Kryuchkova I.V. Series and Fourier transform: method. instructions. Orenburg: Orenburg State University, 2011. 59 p.
15. Egorova G.F., Pavlova G.A., Mazurenko I.A. Fourier series: a workshop on mathematical analysis. Samara: Samara State Technical University. 2011. 46 p.
16. Romanova L.D., Sharkunova T.A., Eliseeva T.V. Integral transformations: textbook. stipend. Penza: Publishing house of PSU, 2015. 80 p.
17. Lossov K.I. Fourier series: studies. manual. M.: MIIGAiK, 2014. 41 p.
18. Kopytin I.I. Methodology for calculating the quality indicators of electric energy. Kurgan: Kurgan State University, 2011.
19. GOST 13109-97. Electrical energy. Electromagnetic compatibility of technical means. Standards for the quality of electrical energy in general-purpose power supply systems.
20. GOST 32144-2013. Electrical energy. Electromagnetic compatibility of technical means. Standards of quality of electric energy in general-purpose power supply systems. M.: Standartinform, 2014.
