Научная статья на тему 'Повышение точности и сокращение времени расчета установившихся режимов электротехнических комплексов повышенной частоты'

Повышение точности и сокращение времени расчета установившихся режимов электротехнических комплексов повышенной частоты Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
53
12
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МЕТОДИКА РАСЧЕТА / ЧАСТОТНЫЙ АНАЛИЗ / ПРЕОБРАЗОВАТЕЛЬ НАПРЯЖЕНИЯ / ПОВЫШЕННАЯ ЧАСТОТА / СОПРОТИВЛЕНИЕ ТРАНЗИСТОРОВ / ЧИСЛЕННОЕ ИНТЕГРИРОВАНИЕ / АНАЛИТИЧЕСКОЕ ВЫРАЖЕНИЕ / ЭДС / РАЗЛОЖЕНИЕ ФУРЬЕ / CALCULATION METHOD / FREQUENCY ANALYSIS / VOLTAGE TRANSFORMER / INCREASED FREQUENCY / RESISTANCE OF TRANSISTORS / NUMERICAL INTEGRATION / ANALYTICAL EXPRESSION / EMF / FOURIER DECOMPOSITION

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Гусенков Алексей Васильевич, Лебедев Владимир Дмитриевич, Шадриков Тимофей Евгеньевич, Танкой Абель, Соколов Александр Михайлович

В современных отечественных и зарубежных работах при моделирования систем переменного тока с несинусоидальной формой питающего напряжения предлагается использовать метод частотного анализа с применением преобразований Фурье. При использовании вышеуказанного метода для расчета установившихся режимов электротехнического комплекса повышенной частоты были выявлены расхождения расчетных и экспериментальных результатов, а также зафиксировано длительное время просчета модели. В целях решения этих проблем актуальной является задача получения обобщенных формульных выражений для величины внутреннего сопротивления IGBT-транзисторов и амплитуды воздействующей ЭДС (напряжения) в составе электротехнического комплекса повышенной частоты при использовании разложения Фурье. Обобщенные выражения для определения внутреннего сопротивления IGBT-транзисторов получены путем анализа физических процессов диффузии зарядов. Амплитуды воздействующей ЭДС (напряжения) определены методом численного интегрирования. Сформированы упрощенные аналитические выражения, пригодные для определения частотной зависимости сопротивления IGBT-транзисторов в открытом состоянии в электротехнических устройствах. Выполнена оптимизация скорости расчета модели электротехнического комплекса повышенной частоты с применением частотного анализа. Полученные обобщенные аналитические выражения позволяют более точно производить вычисление величины внутреннего сопротивления транзистора (разница достигает 70 %). Оптимизация способа расчета в методике позволила сократить время расчета модели электротехнического комплекса повышенной частоты с Nk = 3000 гармоник с 8 часов до 3 минут. Модернизированная методика может быть использована не только для расчета электротехнического комплекса повышенной частоты, но и для анализа любых электрических цепей с компонентами силовой электроники при воздействии непериодических несинусоидальных токов и напряжений.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Гусенков Алексей Васильевич, Лебедев Владимир Дмитриевич, Шадриков Тимофей Евгеньевич, Танкой Абель, Соколов Александр Михайлович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Improving the accuracy and time of calculating steady-state modes of increased frequency electrical systems

The authors of modern Russian and foreign works suggest using the method of frequency analysis with Fourier transforms for modeling alternating current systems with a non-sinusoidal supply voltage. By applying the above-mentioned method to calculating steady-state modes of increased frequency power electrical systems (IFPES), we have found that there are significant differences between the calculated and experimental results, and the model calculation time is long. These problems can be solved by obtaining generalized formula expressions for the internal resistance value of IGBT transistors and the amplitudes of the acting EMF (voltage) in the IFPES applying Fourier decomposition. The generalized expressions for determining the internal resistance of IGBT transistors were obtained by analyzing physical processes of charge diffusion. The amplitude of the acting EMF (voltage) was determined by numerical integration. Simplified analytical expressions suitable for determining the frequency dependence of the resistance of IGBT transistors in the open state in electrical devices have been obtained. The IFPES calculation rate model has been optimized by frequency analysis. The obtained generalized analytical expressions allow making more accurate calculations of the transistor internal resistance value (with a difference of 70 % compared to the previously used value). By optimizing the calculation method we were able to reduce the model calculation from 8 hours to 3 minutes with Nk = 3000 harmonics. The modernized method can be used not only for calculating the IFPES but also for analyzing any electrical circuits with power electronics components exposed to non-periodic and non-sinusoidal currents and voltages.

Текст научной работы на тему «Повышение точности и сокращение времени расчета установившихся режимов электротехнических комплексов повышенной частоты»

ЭЛЕКТРОЭНЕРГЕТИКА

УДК 621.314.228

ПОВЫШЕНИЕ ТОЧНОСТИ И СОКРАЩЕНИЕ ВРЕМЕНИ РАСЧЕТА УСТАНОВИВШИХСЯ РЕЖИМОВ ЭЛЕКТРОТЕХНИЧЕСКИХ КОМПЛЕКСОВ

ПОВЫШЕННОЙ ЧАСТОТЫ

А. ТАНКОЙ1, Т.Е. ШАДРИКОВ1, А.В. ГУСЕНКОВ1, В.Д. ЛЕБЕДЕВ1, А.М. СОКОЛОВ2 1 ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,

г. Иваново, Российская Федерация 2 ФГБОУВО «Ивановский государственный политехнический университет», г. Иваново, Российская Федерация E-mail: Pr3d37@gmail.com

Авторское резюме

Состояние вопроса. В современных отечественных и зарубежных работах при моделирования систем переменного тока с несинусоидальной формой питающего напряжения предлагается использовать метод частотного анализа с применением преобразований Фурье. При использовании вышеуказанного метода для расчета установившихся режимов электротехнического комплекса повышенной частоты были выявлены расхождения расчетных и экспериментальных результатов, а также зафиксировано длительное время просчета модели. В целях решения этих проблем актуальной является задача получения обобщенных формульных выражений для величины внутреннего сопротивления IGBT-транзисторов и амплитуды воздействующей ЭДС (напряжения) в составе электротехнического комплекса повышенной частоты при использовании разложения Фурье.

Материалы и методы. Обобщенные выражения для определения внутреннего сопротивления IGBT-транзисторов получены путем анализа физических процессов диффузии зарядов. Амплитуды воздействующей ЭДС (напряжения) определены методом численного интегрирования. Результаты. Сформированы упрощенные аналитические выражения, пригодные для определения частотной зависимости сопротивления IGBT-транзисторов в открытом состоянии в электротехнических устройствах. Выполнена оптимизация скорости расчета модели электротехнического комплекса повышенной частоты с применением частотного анализа.

Выводы. Полученные обобщенные аналитические выражения позволяют более точно производить вычисление величины внутреннего сопротивления транзистора (разница достигает 70 %). Оптимизация способа расчета в методике позволила сократить время расчета модели электротехнического комплекса повышенной частоты с Nk = 3000 гармоник с 8 часов до 3 минут. Модернизированная методика может быть использована не только для расчета электротехнического комплекса повышенной частоты, но и для анализа любых электрических цепей с компонентами силовой электроники при воздействии непериодических несинусоидальных токов и напряжений.

Ключевые слова: методика расчета, частотный анализ, преобразователь напряжения, повышенная частота, сопротивление транзисторов, численное интегрирование, аналитическое выражение, ЭДС, разложение Фурье

IMPROVING THE ACCURACY AND TIME OF CALCULATING STEADY-STATE MODES OF INCREASED FREQUENCY ELECTRICAL SYSTEMS

A. TANKOI1, T.E. SHADRIKOV1, A. GUSENKOV1, V. LEBEDEV1, A. SOKOLOV2 1 Ivanovo State Power Engineering University, Ivanovo, Russian Federation 2 Ivanovo State Polytechnic University, Ivanovo, Russian Federation E-mail: Pr3d37@gmail.com

Abstract

Background. The authors of modern Russian and foreign works suggest using the method of frequency analysis with Fourier transforms for modeling alternating current systems with a non-sinusoidal supply voltage. By applying the above-mentioned method to calculating steady-state modes of increased frequency power electrical systems (IFPES), we have found that there are significant differences between the calculated and experimental results, and the model calculation time is long. These problems can be solved by obtaining generalized formula expressions for the internal resistance value of IGBT transistors and the amplitudes of the acting EMF (voltage) in the IFPES applying Fourier decomposition.

Materials and methods. The generalized expressions for determining the internal resistance of IGBT transistors were obtained by analyzing physical processes of charge diffusion. The amplitude of the acting EMF (voltage) was determined by numerical integration.

Results. Simplified analytical expressions suitable for determining the frequency dependence of the resistance of IGBT transistors in the open state in electrical devices have been obtained. The IFPES calculation rate model has been optimized by frequency analysis.

Conclusions. The obtained generalized analytical expressions allow making more accurate calculations of the transistor internal resistance value (with a difference of 70 % compared to the previously used value). By optimizing the calculation method we were able to reduce the model calculation from 8 hours to 3 minutes with Nk = 3000 harmonics. The modernized method can be used not only for calculating the IFPES but also for analyzing any electrical circuits with power electronics components exposed to non-periodic and non-sinusoidal currents and voltages.

Key words: calculation method, frequency analysis, voltage transformer, increased frequency, resistance of transistors, numerical integration, analytical expression, EMF, Fourier decomposition

DOI: 10.17588/2072-2672.2019.3.022-031

Введение. Разработанная ранее методика электрического расчета установившегося режима электротехнических комплексов повышенной частоты (ЭТКПЧ) [1], позволяющая вычислять токи и напряжения в различных узлах схемы замещения комплекса, основанная на применении частотного анализа [2], в настоящее время может быть подвергнута доработкам в целях повышения достоверности расчета и значительного увеличения скорости вычислений. Объектами применения методики, например, могут выступать локальные ЭТКПЧ, предназначенные для организации электротепловой обработки (ЭТО) на промышленных предприятиях строительной отрасли. Такой вид электротепловой обработки обладает высокими энергетическими и технико-экономическими характеристиками и может быть успешно реализован на строительных предприятиях, выпускающих готовые (типовые) железобетонные изделия.

Электрическая схема замещения такого комплекса показана на рис. 1, 2.

Выполненные ранее вычисления [1, 2] с использованием указанной выше методики давали завышенные значения высокочастотных составляющих токов и напряжений, особенно в режиме холостого хода ЭТКПЧ. Возможной причиной такого результата является недостаточно точное представление зависимости внутреннего сопротивления транзисторов от частоты гармоник разложения Фурье [4-8].

Другим недостатком методики [1] является большое время, необходимое на вычисление варианта схемы ЭТКПЧ, так как процесс расчета связан с многократным вычислением большого числа определенных интегралов численным методом из-за использования разложения Фурье с числом составляющих, достигающим 3000 и более [1].

С учетом перечисленных недостатков методика расчета установившихся режимов ЭКТПЧ должна быть модернизирована.

Ц 12 и N

Выпрямитель

Полупроводниковый п реобразовател ь

Л 4

Повышающий трансформатор

УГ3

4 у-

220/380 50 Гц

Тр1

Кабельная линия

Понижающий трансформатор

и

Тр2

Нагрузка

"Г"

Реактор

Рис. 1. Электрическая схема ЭТКПЧ: С - емкость на шинах постоянного тока для сглаживания пульсаций; VT1-VT4, VD1-VD4 - силовые транзисторы и диоды полупроводникового преобразователя соответственно; Тр1, Тр2 - повышающий и понижающий силовые трансформаторы соответственно; 1к -длина кабельной линии электропередачи; Цр - реактор

цен я:

цэв я2

11

¡2

'цэв я2

цэн я1

Г 1-0 яо

0, со±

-о яо со±

-о Яо

=^оо

+кэ

и^-отв -*

Отф ЯнП

Рис. 2. Схема замещения ЭТКПЧ: 0к - источник напряжения к-й гармоники; 2вн - внутренне сопротивление источника напряжения; Отн и Отв - емкость обмоток низкого и высокого напряжения силовых трансформаторов соответственно; -5н и -зв - индуктивности рассеивания обмоток низшего и высокого напряжения силовых трансформаторов соответственно; Я1 и Я2 - активные сопротивления обмоток низкого и высокого напряжения силовых трансформаторов соответственно; Яц и -ц - активное сопротивление и индуктивность ветви намагничивания трансформатора; Яр - активное сопротивление реактора; -р - индуктивность реактора; Ян - сопротивление (активное) нагрузки; /К - длина высоковольтной кабельной линии, представленной в виде линии с распределенными параметрами; Я0, -0, О0 - погонные параметры линии (активное сопротивление, индуктивность и емкость на единицу длины)

Материалы и методы. Уточнение влияния внутреннего сопротивления транзистора в преобразователе напряжения. Противо-ЭДС. В [1, 4-8] показано, что, исходя из существующих представлений о диффузионном распространении неосновных носителей заряда от эмиттерного перехода к коллекторному при работе транзистора и явления запаздывания появления тока коллектора относительно тока эмиттера, может быть получена формула (1), описывающая зависимость внутреннего сопротивления транзистора от кратности гармоники:

7 £вн0

а(ю)

= £

вн0<

1 +

(кО

V feJ

(1)

Аит

где 2вн0 =-Г^

- внутреннее сопротивле-

к 0

ние транзистора в открытом состояния при постоянном напряжения; к - 1, 3, 5, ... -кратность гармоники (в кривой ЭДС присутствуют только нечетные гармоники [1]); ^ - частота первой гармоники; 4 - пре-

дельная частота усиления (паспортное значение); |а(ю)| - модуль коэффициента передачи транзистора по току; диг - падение напряжения на открытом транзисторе при постоянном токе ¡ко (паспортные значения).

Формула (1) получена с использованием широко распространенного [3-8] способа вычисления коэффициента передачи транзистора по току а(ю):

а(ю) = ■

1 + I-

1 + У

. 1

(2)

4

где а0 « 1 - коэффициент передачи по току на постоянном напряжении; f - произвольное значение частоты; ю = 2%/; юа = 2%^.

Согласно (1), с увеличением кратности гармоники внутреннее сопротивление транзистора постоянно увеличивается. Существующая методика электрического расчета установившихся режимов работы ЭТКПЧ [1], в основе которой лежит применение формулы (1), согласно экспериментальным данным [1, 11], дает завышенные

с

и

-

р

к

4

ю

2

значения высокочастотных составляющих токов и напряжений. Таким образом, рассмотренная зависимость 2внк от кратности гармоники, по сравнению с тем влиянием, которое определяется соотношением (1), должна быть более выраженной.

Существует выражение [8], позволяющее более точно учесть влияние частоты на величину коэффициента передачи по току. Формулу (2) можно записать в показательной форме:

а(ю) = ■

ап

V

1 +

1

(3)

где р = эг^д — сдвиг (отставание) по фазе

юа

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

тока коллектора от тока эмиттера вследствие относительно медленного процесса диффузионного распространения неосновных носителей от эмиттерного перехода к коллекторному.

Как показано в [5-7] на основе детального анализа процесса диффузии неосновных носителей в базовой области, угол сдвига по фазе при ю = юа должен составлять р = 60о , однако из выражения (3) в этом случае (ю/юа=1) этот угол имеет величину 45о, т.е. формулы (3) и (2) не совсем адекватно отражают фактическую зависимость коэффициента передачи по току от частоты.

Для устранения этого противоречия необходимо ввести поправочный коэффициент т безразмерного характера:

а(ю) = ■

1 +1

ю

У

т

(4)

Тогда угол сдвига по фазе составит

р = эг^д—т.

ю„

(5)

Полагая в этой формуле р = 600 , т.е. р = л/3, а также ю = юа, из соотношения (5) можно определить значение этого безразмерного коэффициента:

т = tg- = 1,73. 3

(6)

Коэффициент передачи по току уже с поправочным коэффициентом может быть записан следующим образом:

а>

(ю)| = ■

1 +

т

1 +

V 'а

т

Применяя совместно формулы (1) и (7), нетрудно получить формулу для определения уточненного значения внутреннего сопротивления транзистора в открытом состоянии:

1 +

к1 V '

'-т

(8)

Если (8) в комплексной форме умножить на величину тока коллектора /«о, то получится следующее выражение:

1к0 ^внк = 0 ^вн0

или

1++А*-1)

а 'а

АУТ (ю) = Лиг|1 +1 у1 | + Евн(ю),

(9)

(10)

где ЛЦ-(ю) - падение напряжения на открытом транзисторе при произвольных значениях частоты напряжения и кратности гармоники; Евн(ю) - противо-ЭДС, отражающая скорректированное значение внутреннего сопротивления (8):

^вн(ю) = Л 0 • гвН0' = 1 Шт ' (11)

'а 'а

Выражения (8) и (10) позволяют предложить новый вариант схемы замещения транзисторного преобразователя напряжения, представленный на рис. 3, в том случае, если для определения внутреннего сопротивления используется выражение (1).

Если 7внк определятся с помощью формулы (8), то противо-ЭДС Евн(ю) необходимо исключить из схемы замещения транзисторного мостового преобразователя напряжения.

Т) Евн(ю)

зн

внк

-0

Рис. 3. Схема замещения транзисторного преобразователя напряжения с использованием противо-ЭДС

Для сравнения результатов расчета по формулам (1) и (8) целесообразно выразить величину внутреннего сопротивления в относительных значениях:

а

0

Vюа

7* = ^вк ^внк у

^вн0

= , 1 +

к1 V '

т

(12)

где т = 1 - для прежней методики расчета с использованием соотношения (1); т = 1,73 - для нового способа определения (8) величины внутреннего сопротивления.

На рис. 4 представлены зависимости, полученные с использованием соотношений (1) и (12).

Анализ полученных зависимостей (рис. 4) показывает, что применение соотношения (8) дает значения величины внутреннего сопротивления на 70 % выше, чем с использованием ранее рассмотренной методики (1).

ое

30 -

К.. '2

1

Г / /.

10 12 14 16

18 ое

Рис. 4. Зависимости величины внутреннего сопротивления транзисторного преобразователя от отношения Ма = 1 - т = 1; 2 -т = 1,73

В , и-2 2000

а)

б)

СН 500и НОИНАЬ

♦ир 1гоо :одар-1220. ошр-

РВЕ£) 17.20кНг : СУС 5В. Имэ

Работа под нагрузкой

В 4000

в)

-4000

В

3000

И 10. Одз АС

СН ИШД

♦1!р 1Ш ОУир -1180.0Шр-р 2360.01) ГВЕЦ 1?. 24кНг : СУС 57. ЭВцэ; ■■■•!••••

Режим холостого хода

Рис. 5. Расчетные (а, б) и экспериментальные осциллограммы (в) напряжения на высоковольтной кабельной линии: а - при использовании в расчетах ранее предложенной формулы (1); б - при использовании формулы (8)

25

20

15

10

4

Для оценки результатов применения уточненной методики был выполнен расчет установившегося режима ЭТКПЧ. На рис. 5 представлены расчетные осциллограммы напряжения высоковольтной кабельной линии в режиме ХХ и под нагрузкой для макета ЭТКПЧ, описанного в [9, 10]. Эти осциллограммы наглядно демонстрируют тот факт, что при использовании нового способа определения 2внк обеспечивается лучшее согласование результатов расчета и опыта. Применение в расчетах ранее предложенной формулы приводит к явно завышенным значениям амплитуды высокочастотных колебаний, т.е. величины внутренних перенапряжений. Это приводит, в свою очередь, к излишним запасам электрической прочности изоляции элементов ЭТКПЧ (высоковольтных силовых трансформаторов [10] и кабельной линии [9]).

Совершенствование методики расчета параметров гармоник при разложении кривой ЭДС в ряд Фурье в целях ускорения вычислений. Другой задачей совершенствования методики электрического расчета установившихся режимов работы ЭТКПЧ является ускорение вычислений. Для ее решения следует обратить внимание на процесс определения параметров гармоник при разложении кривой ЭДС в ряд Фурье. Как показано в [1], наилучшим вариантом представления реальной осциллограммы ЭДС является стилизованная кривая ЭДС, показанная на рис. 5. Согласно [1, 11-15], амплитуда гармоник при разложения исходного сигнала (рис. 6) в ряд Фурье определяется посредством его интегрирования в течение одного периода колебаний 0^7 с использованием формулы

ик = 2 |ЦЦ) • )с№ =

(13)

= 2итГ1('1 + >2 + >3 + ¡4 + >5 +16 + >7 + >8 + У,

где и(И) - аналитическое описание исходного сигнала.

При выполнении интегрирования (13) кривую ЭДС можно рассматривать как функцию, имеющую точки разрыва. В этом случае интеграл от такой функции определяется как сумма интегралов в пределах тех участков, где она остается непрерывной и может быть описана аналитически [14]. Как показано на рис. 6, в рассматриваемом случае можно выделить 9 таких участков.

Рис. 6. Стилизованная осциллограмма и параметры ЭДС транзисторного преобразователя напряжения и(/) «трапеция с паузой»: 7т = (7 -2т) - длительность амплитудного значения при т = 2/1 +12; ит - амплитудное значение ЭДС

Таким образом, интеграл за один период изменения ЭДС во времени необходимо представить, как сумму 9 слагаемых 11^19 (13), каждое из которых соответствует одному из участков, отмеченных на рис. 6, полагая, что в течение паузы (участки 1, 5, 9) ЭДС и^) = 0, на участках 3 и 7 - иф = ит и иф = -ит, а на участках 2, 4, 6, 8 имеет место изменение напряжения по линейному закону. Были составлены выражения для вычисления интегралов формулы (13) на каждом участке:

Ь

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2

1 = 0 • | вт(2л:(2к -1/)= 0;

+* 2 1

2 |

2

2

|+'1+т

ь

вт(2л(2к - 1)/Ц;

/3= | вт(2л(2к - 1)/Ц;

^+*

(14)

|+2 т

>4 =

|+т

1 - t - и + 4 + *

Б1п(2я:(2к - 1/1;

+2^+т

>5= 0 •

| вт(2л(2к - 1)/Ц = 0;

2+2*1+т

3 12

+3 ч+т

( г

>6 =

32

+2 (1+т

-1

^ -1 +

sin(2л(2k -1/);

2

2

2

3 «1+2 т

>7= | (-1^т(2ж(2к -1/1;

32+3*1+т

2 +4(1+2(т

((

8 = |

f - 3, - 32 - 2,

л Л

-1

2+3и+^т

ь+4 -л+2 т

sin(2ж(2k -1/1 1;

>9= 0 • | sin(2ж(2k -1/1 = 0.

2 + 4*1+2 т

Согласно (14),

>1 = >5 = >9 = 0,

(15)

что упрощает использование соотношения (13). Для вычисления определенных интегралов 13 и 17 можно применить известную формулу [10]:

2

| а • sin(bt с-

5

= -а

cos(bt) а

(16)

Ь Ь

= — [COS(bt1) - COS(bt2)].

Интегралы 14 , 16 , ¡8 можно свести к следующему виду:

"2

5 ^ - а)

sin(bt)

о ■ Г ьЬ 2а • Sinl 1

^ =

Ьс

sin(bt1) + Ь^

■ Гbt112 ^

2 • Sin| 1 -1

(17)

Ь2 с

(

sin(bt2) + Ьt2

2 • 12 -1'

Ь2 с

С учетом (16) и (17) выражение (14) может быть вычислено без процедуры численного интегрирования лишь с использованием аналитических выражений.

Например, интеграл ¡3 может быть представлен следующим аналитическим выражением:

1

сс^| (2ж(2к -1)/)| ^ + ^

-^| (2ж(2к -1)/)| ^ + ^ + 7т

(2ж(2к -1)/)

(18)

Интеграл ¡4 может быть выражен формулой

2 к+2+т+

( г* \\

т

>4 =

(2ж(2к -/ + -л + 7

(2ж(2к - ОДНО

sin((2ж(2k -1)/,) д + (2ж(2к -/ I + ^ + Тт

[(2ж(2к -1)/1) + и + Тт^

2 • sin

2 ^ -1

(2ж(2к -1)f1)2(--1)

*п | (2ж(2к -1)/1) 11 + 2^ + Тт) I + (2ж(2к - / Г| + 2-1 + Тп

[(2ж(2к -1)/1) + 2-1 + Тт ^

2 • sin

2

-1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(2ж(2к - 1)/i)2(--1)

(19)

Полученные выражения вида (18)-(19) позволяют при определении амплитуды гармоник ЭДС с помощью формулы (13) полностью отказаться от процедуры численного интегрирования.

Результаты. Несмотря на то, что полученное аналитическое выражение представляет собой набор достаточно громоздких формул (18), их применение в расчетных программах вместо использо-

>

3

2

+

вания численного интегрирования привело к многократному снижению времени вычислении (рис. 7). В качестве тестовой системы использован ПК на базе Intel® Core i3-4005U @ 1.70 GHz с поддержкой Hyper Threading, 8ГБ ОЗУ.

Количество учитываемых гармоник, N

Рис. 7 Длительность вычисления характеристик установившегося режима работы ЭТКПЧ в зависимости от количества учитываемых гармоник: светлый цвет - использование аналитических выражений; темный цвет - использование численного интегрирования

Согласно полученным данным (рис. 7), с увеличением числа гармоник разложения Ык различие в длительности вычислений разными способами быстро нарастает. Причем если при числе гармоник Ык = 100 длительность вычислений оказывается почти одинаковой и не превышает 2 мин, то при использовании рекомендуемого значения Ык= 3000 [1] длительность вычислений по аналитическим выражениям не превышает 3 мин, а при численном интегрировании достигает 8 часов.

Выводы. Полученное уточненное выражение для определения величины внутреннего сопротивления полностью открытого транзистора в зависимости от частоты гармоник разложения кривой ЭДС транзисторного преобразователя напряжения обязательно при выполнении электрического расчета установившихся режимов работы ЭТКПЧ с источником питания на основе транзисторного преобразователя напряжения.

Полученное аналитическое выражение для вычисления интеграла разложения Фурье при определении амплитуды гармоник, входящих в состав кривой ЭДС «трапеция с паузой» преобразователя, позволяет снизить временные затраты на вычисление.

Таким образом, методика электрического расчета установившихся режимов

работы ЭТКПЧ усовершенствована посредством перехода на вычисление аналитических выражений вместо процедуры многократного численного интегрирования, что позволяет значительно уменьшить затраты времени на выполнение расчетов.

Целесообразно дальнейшее развитие и совершенствование методики расчета установившихся режимов ЭТКПЧ, например, в целях учета нелинейных характеристик некоторых элементов схемы замещения ЭТКПЧ.

Список литературы

1. Особенности применения частотного анализа при расчете электрических цепей с транзисторными преобразователями напряжения / А.В. Гусенков, В.Д. Лебедев, А.М. Соколов и др. // Электричество. - 2016. - № 1. - С. 4-12.

2. Осипов Ю.М. Частотный и временной анализ стационарных и переходных характеристик линейных электрических цепей: учеб. пособие по курсам электротехники и ТОЭ. Ч. 2. -СПб.: СПбГИТМО (ТУ), 2002. - 99 с.

3. Розанов Ю.К. Полупроводниковые преобразователи со звеном повышенной частоты. - М.: Энергоатомиздат, 1987. -184 с.

4. Гуртов В.А. Твердотельная электроника: учеб. пособие / ПетрГУ. - Петрозаводск, 2004. - 312 с.

5. Розанов Ю.К. Основы силовой электроники. - М.: Энергоатомиздат, 1992. - 296 с.

6. Comparative analysis of static characteristics of insulated gate bipolar transistors and thyristors with static induction / O.I. Bonomorskii, A.S. Kyuregyan,

A.V. Gorbatyuk, B.V. Ivanov // Russian Electrical Engineering. - 2015. - Т. 86, № 2. - С. 93-97.

7. Моделирование статических и динамических потерь в MOSFET-ключах / В.П. Ба-бенко, В.К. Битюков, В.В. Кузнецов, Д.С. Си-мачков // Российский технологический журнал. -2018. - Т. 6, № 1 (21). - С. 20-39.

8. О возможности и целесообразности физического моделирования электропередачи с нетрадиционными параметрами используемых токов и напряжений / А.В. Гусенков,

B.Д. Лебедев, А.М. Соколов, Т.Е. Шадриков // Энергетик. - 2015. - № 4. - С. 29-32.

9. Выбор рабочего напряжения высоковольтной кабельной электропередачи повышенной частоты / А.В. Гусенков, В.Д. Лебедев, Т.Е. Шадриков, А.М. Соколов // Электротехника. - 2016. - № 10. - С. 50-56.

10. Методика определения емкостных и индуктивных параметров силовых высоковольтных трансформаторов повышенной частоты / Т.Е. Шадриков, А.В. Гусенков, А.А. Дьячков,

А.А. Симакова // Вестник ИГЭУ. - 2016. -Вып. 1. - С. 27-33.

11. Спектральный анализ электродвижущей силы электроустановки повышенной частоты / А.В. Гусенков, В.Д. Лебедев, А.М. Соколов и др. // Электричество. - 2018. -№ 2. - С. 17-26.

12. Стариков А.В., Кузнецов В.В., Ро-кало Д.Ю. Анализ гармонического состава трапецеидального фазного напряжения, формируемого частотным преобразователем // Вестник Самарского государственного технического университета. Сер.: Технические науки. - 2017. - № 3(55). - С. 75-79.

13. Садиков Д.Г., Титов В.Г. Анализ гармонического состава тока и напряжения, потребляемого преобразователями частоты // Вестник Чувашского университета. - 2015. -№ 1. - С. 116-121.

14. Зырянов В.М., Митрофанов Н.А., Соколовский Ю.Б. Исследование гармонического состава напряжения преобразователя частоты // Электрооборудование: эксплуатация и ремонт. - 2016. - № 9. - С. 34-40.

15. Gerard Hurley, William & Gath, Eugene & Breslin, John. (2000). Optimizing the AC resistance of multilayer transformer windings with arbitrary current waveforms // Power Electronics, IEEE Transactions on. 15. 369-376. doi: 10.1109/63.838110.

References

1. Gusenkov, A.V., Lebedev, V.D., Sokolov, A.M., Shadrikov, T.E., Strahov, A.S. Osobennosti primeneniya chastotnogo analiza pri raschete elektricheskikh tsepey s tranzistornymi preobrazovatelyami napryazheniya [Frequency analysis in the calculation of electrical circuits with transistor voltage converters]. Elektrichestvo, 2016, no. 1, pp. 4-12.

2. Osipov, Yu.M. Chastotnyy i vremennoy analiz statsionarnykh i perekhodnykh kharakteristik lineynykh elektricheskikh tsepey. Ch. 2 [Frequency and time analysis of stationary and transient characteristics of linear electric circuits: a study guide for Electrical Engineering and Theory of Electrical Engineering. Part 2]. Saint-Petersburg: SPbGITMO (TU), 2002. 99 p.

3. Rozanov, Yu.K. Poluprovodnikovye preobrazovateli so zvenom povyshennoy chastoty [Semiconductor converters with an increased frequency unit]. Moscow: Energoatomizdat, 1987. 184 p.

4. Gurtov, V.A. Tverdotel'naya elektronika [Solid state electronics]. Petrozavodsk, 2004. 312 p.

5. Rozanov, Yu.K. Osnovy silovoy elektro-niki [Power Electronics Basics]. Moscow: Energoatomizdat, 1992. 296 p.

6. Bonomorskii, O.I., Kyuregyan, A.S., Gor-batyuk, A.V., Ivanov, B.V. Comparative analysis of

static characteristics of insulated gate bipolar transistors and thyristors with static induction. Russian Electrical Engineering, 2015, vol. 86, no. 2, pp. 93-97.

7. Babenko, V.P., Bityukov, V.K., Kuznetsov, V.V., Simachkov, D.S. Modelirovanie staticheskikh i dinamicheskikh poter v MOSFET-klyuchakh [Simulation of static and dynamic losses in MOSFET transistors]. Rossiyskiy tekhnologicheskiy zhurnal, 2018, vol. 6, no. 1(21), pp. 20-39.

8. Gusenkov, A.V., Lebedev, V.D., Sokolov, A.M., Shadrikov, T.E. O vozmozhnosti i tselesoobraznosti fizicheskogo modelirovaniya el-ektroperedachi s netraditsionnymi parametrami ispol'zuemykh tokov i napryazheniy [On possibility and feasibility of physical modeling of power transmission with non-traditional parameters of currents and voltages]. Energetik, 2015, no. 4, pp. 29-32.

9. Gusenkov, A.V., Lebedev, V.D., Shadri-kov, T.E., Sokolov, A.M. Vybor rabochego napryazheniya vysokovol'tnoy kabel'noy el-ektroperedachi povyshennoy chastoty [Selection of operating voltage of high-voltage cable in-creased-frequency power transmission]. Elektrotekhnika, 2016, no. 10, pp. 50-56.

10. Shadrikov, T.E., Gusenkov, A.V., D'yachkov, A.A., Simakova, A.A. Metodika opre-deleniya emkostnykh i induktivnykh parametrov silovykh vysokovol'tnykh transformatorov pov-yshennoy chastoty [Method of determining capacitance and inductance parameters of increased frequency high-voltage power transformers]. Vestnik IGEU, 2016, issue 1, pp. 27-33.

11. Gusenkov, A.V., Lebedev, V.D., Sokolov, A.M., Shadrikov, T.E., Abel', Tankoj. Spektral'nyy analiz elektrodvizhushchey sily el-ektroustanovki povyshennoy chastoty [Spectral analysis of electromotive force of an increased frequency electrical installation]. Elektrichestvo, 2018, no. 2, pp. 17-26.

12. Starikov, A.V., Kuznetsov, V.V., Roka-lo, D.Yu. Analiz garmonicheskogo sostava trapet-seidal'nogo faznogo napryazheniya, formiruemogo chastotnym preobrazovatelem [Analysis of harmonic composition of the trapezoidal phase voltage generated by a frequency converter]. Vestnik Samarskogo gosudarstvennogo tekh-nicheskogo universiteta. Seriya: Tekhnicheskie nauki, 2017, no. 3(55), pp. 75-79.

13. Sadikov, D.G., Titov, V.G. Analiz garmonicheskogo sostava toka i napryazheniya, po-treblyaemogo preobrazovatelyami chastoty [Analysis of harmonic composition of current and voltage consumed by frequency converters]. Vestnik Chu-vashskogo universiteta, 2015, no. 1, pp. 116-121.

14. Zyryanov, V.M., Mitrofanov, N.A., Sokolovskiy, Yu.B. Issledovanie garmonicheskogo sostava napryazheniya preobrazovatelya chastoty [A study of harmonic composition of a frequency converter voltage]. Elektrooborudovanie: eksplu-atatsiya i remont, 2016, no. 9, pp. 34-40.

15. Gerard Hurley, William & Gath, Eugene & Breslin, John. (2000). Optimizing the AC resistance of multilayer transformer windings with

arbitrary current waveforms. Power Electronics, IEEE Transactions on. 15. 369-376. 10.1109/63.838110.

Гусенков Алексей Васильевич,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент, проректор по учебной работе, телефон (4932) 38-57-17, e-mail: avgus@ispu.ru Gusenkov Aleksei Vasilyevich,

Ivanovo State Power Engineering University, Candidate of Engineering Sciences (PhD), Associate Professor, Vice-Rector for Academics, telephone (4932) 38-57-17, e-mail: avgus@ispu.ru

Лебедев Владимир Дмитриевич,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент, зав. кафедрой автоматического управления электроэнергетическими системами, телефон (4932) 26-99-04, e-mail: vd_lebedev@mail.ru Lebedev Vladimir Dmitrievich,

Ivanovo State Power Engineering University, Candidate of Engineering Sciences (PhD), Associate Professor, Head of the Department of Electric Power Systems Automatic Control, telephone (4932) 26-99-04, e-mail: vd_lebedev@mail.ru

Шадриков Тимофей Евгеньевич,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических

наук, доцент кафедры высоковольтной электроэнергетики, электротехники и электрофизики,

телефон (4932) 26-97-23, e-mail: timofeyshadrikov@gmail.com

Shadrikov Timofei Yevgenyevich,

Ivanovo State Power Engineering University,

Candidate of Engineering Sciences (PhD), Associate Professor of the Department of High Voltage Power Engineering, Electrical Engineering and Electrophysics, telephone (4932) 26-97-23, e-mail: timofeyshadrikov@gmail.com

Танкой Абель,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», соискатель, телефон (4932) 26-97-23. Tankoi Abel,

Ivanovo State Power Engineering University, External post-graduate student, telephone (4932) 26-97-23. Соколов Александр Михайлович,

ФГБОУВО «Ивановский государственный политехнический университет», доктор технических наук, доцент, телефон (4932) 26-99-04, e-mail: alex2010fn@yandex.ru Sokolov Aleksandr Mikhailovich,

Ivanovo State Polytechnic University, Doctor of Engineering Sciences (Post-doctoral degree), Associate Professor, telephone (4932) 26-99-04, e-mail: alex2010fn@yandex.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.