CC BY
65Транспорт
УДК 621.314.5
DOI: 10.12737/article 59b11cbccce9f4.34179572
И.Ю.Бутарев
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ТРЕХФАЗНОГО КОРРЕКТОРА КОЭФФИЦИЕНТА МОЩНОСТИ НА ОСНОВЕ ИНВЕРТОРА НАПРЯЖЕНИЯ
Проведено моделирование электромагнитных процессов трехфазного корректора коэффициента мощности на основе инвертора напряжения. Разработана модель трехфазного корректора коэффициента мощности. Проведен бифуркационный анализ динамических процессов в типовом коррек-
торе коэффициента мощности на основе инвертора напряжения.
Ключевые слова: корректор коэффициента мощности, PFC, инвертор напряжения, моделирование, силовая электроника, бифуркация, нелинейная динамика.
I.Yu. Butarev
SIMULATOR OR THREE-PHASE CORRECTOR OF POWER FACTOR BASED ON VOLTAGE INVERTER
The model of a three-phase corrector of a power factor is developed. A bifurcation analysis of dynamic processes in a standard corrector of a power factor on the basis of a voltage inventor is carried out. Three-phase correctors of a power factor are intended for the compensation of a negative impact of such loads upon a mains supply. A three-phase corrector of a power factor on the basis of a voltage inverter is a system complex enough, and for the correct control a definition of instant values of voltage and current at each phase of input sine-wave voltage and output voltage is required. The presence of a feedback in the corrector of a power factor may cause nonlinear effects (bifurcations), which may result in non-designed modes of operation and converter failure. That is why it is re-
quired that a simulator taking into account the non-designed mode appearance should be developed. The presented simulator of electromagnetic processes of the three-phase corrector of a power factor on the basis of a voltage inverter allowed defining a possibility of undesirable mode appearance at rated operation modes. The example of a simulation with the aid of Math-worksMatlab program is shown, bifurcation diagrams are formed and the analysis of the impact of load resistance in the designed mode of corrector operation is carried out. The recommendations for the choice of three-phase corrector basic parameters at engineering designing are shown.
Key words: power factor corrector, PFC, voltage inverter, simulation, power electronics, bifurcation, nonlinear dynamics.
Для повышения коэффициента мощности, устранения высших гармоник в сети применяют пассивные и активные корректоры [2; 3] коэффициента мощности. Активный корректор коэффициента мощности (APFC - Active Power Factor Correction) имеет в своем составе накопительный дроссель, ключевые транзисторы и управляющую микросхему. АККМ по сути является импульсным преобразователем, имеющим за счет высокой частоты работы на выходе практически правильную синусоиду. Для электросети общего пользования блок питания, оснащенный модулем АККМ, представляется в виде обычного активного сопротивления. При использовании АККМ коэффициент мощ-
ности блоков питания находится, как правило, в диапазоне 0,95.. .0,99 (95... 99%).
Наличие нелинейных элементов в составе АККМ может привести к возникновению автоколебаний, бифуркаций, хаотическому режиму работы и выходу устройства из строя [4; 5]. Необходимо оценить условия и режимы работы, способствующие возникновению бифуркаций.
Трехфазные корректоры коэффициента мощности призваны скомпенсировать негативное влияние таких нагрузок на питающую сеть. Топологий для реализации силовой части корректоров коэффициента мощности достаточно много. Один из вариантов - это использование схемы трехфазного корректора коэффициента мощности на основе инвертора напряжения [1].
Трехфазный корректор коэффициента мощности на основе инвертора напряжения является достаточно сложной системой, и для корректного управления требуется определение мгновенных значений напряжения и тока на каждой фазе входного синусоидального напряжения и напряжения выхода. Схема трехфазного корректора мощности на основе инвертора представлена на рис. 1.
Наличие обратной связи в корректоре коэффициента мощности может вызывать нелинейные эффекты (бифуркации) [6], которые могут привести к непроектным режимам работы и выходу из строя преобразователя. Поэтому требуется разработать математическую модель, учитывающую возникновение непроектных режимов.
Исследуемая модель реализована в виде компьютерной программы в Math-
works Matlab в блочной модели (рис. 2) корректора коэффициента мощности в Matlab Sim Power Systems. Блочная модель позволяет исследовать поведение корректора коэффициента мощности с учетом обратной связи и моделировать ситуации, при которых возможно проявление нелинейности и хаоса. К примеру, на рис. 3 представлен результат моделирования, при котором получен переход из проектного режима в хаотический.
Для ускорения расчетов и возможности анализа нелинейной динамики корректора можно создать самостоятельно математическую модель. Математическая модель данного преобразователя представляет собой систему дифференциальных уравнений, построенных в базисе коммутационно-разрывных функций.
Рис. 1. Принципиальная схема трехфазного корректора коэффициента мощности на основе инвертора напряжения: Яи, Я2, Яз - активное сопротивление дросселей соответствующих фаз; Ьи, Ь3 - индуктивность дросселей соответствующих фаз; С - емкость конденсатора; Ян - сопротивление нагрузки; и¡, и2, и3 - входное трехфазное напряжение; УТ1-УТ6 - силовые транзисторные ключи IGBT; УВ1-УБ6 - антипараллельные диоды
ключей IGBT
При создании математической модели представленного трехфазного корректора коэффициента мощности можно условно обозначить 7 возможных вариантов схем замещения, которые могут присутствовать на тактовом интервале работы. По топологии 7 возможных вариантов могут быть сведены к 3 (нагрузка рассоеди-
нена с источником (происходит разряд конденсатора на нагрузку), включены два ключа сверху и ключ снизу, включены ключ сверху и два ключа снизу). При этом на каждом тактовом интервале возможно до 3 коммутаций, каждая из которых изменяет топологию схемы. Каждая из схем замещения описывается системой диффе-
ренциальных уравнений, которая в матричной форме имеет вид
f = А,Х + В„ (1)
где Аг - матрица постоянных коэффициентов на г-м участке размером 3*3; Вг -вектор вынуждающих воздействий на г-м участке; Х=[/1,/2,ис] - вектор переменных состояния.
Рис. 2. Блочная модель трехфазного корректора коэффициента мощности в Matlab Simulink
Рис. 3. Проявление бифуркации в трехфазном корректоре коэффициента мощности
Рассмотрим первую схему замещения, представленную на рис. 4. Эта схема является частным случаем схемы, представленной на рис. 1 , при открытых силовых ключах VT1,VT4,VT6. В данном слу-
чае открыты верхний ключ в первой стойке и два нижних ключа во второй и третьей стойках.
Система уравнений (1) для этой схемы замещения может быть представлена с
с11
помощью законов Кирхгофа:
-и, + М + Ь^-Л. + ис + и2 - ад - Ь2 = 0;
Ж
ад
-и2 + К12 + + из - К/3 - Ь3 = 0;
= 1с + Ь;
/ + !г + /3 = 0.
(2)
Приведем ее к форме Коши, поочередно находя зависимость производных переменных состояния от самих переменных.
Пусть V = -
—2 ^ ЬЬ ^ -2-3 дем систему (2) к форме Коши:
Приве-
Ж/1
ж
= (- + - )уЦ - ^и2 - -2уи3 - (Ь + - )уис -
-(— +—+ —уЩ//г + ьуЩЬ; ад
— = --УЦ + (- +—)уи2 - + ь^ис + &
+-УЯ/ - (— + — )УЯ2/2 + 1уЯ3/3;
=-1 /,--1- ис.
С СЯ с
(3)
ж
Система уравнений (3) приведена к форме (1), следовательно, можно выделить Ay - матрицу постоянных коэффициентов для соответствующей рис. 2 схемы заме-
А =
щения трехфазного корректора коэффициента мощности, B7 - вектор вынуждающих воздействий для соответствующей схемы замещения.
- [((-2 + -3 К + -2К >] (-К - -2К )У -(Ь2 + -3 )У
(-3К - АК )у -[((-, + -3 )К2 + АЯ3 )У ]
1
С Б,
0
(- + -) ---1уих + (- + -)уи2 - -уиз
-V
-(С^}
;(4)
0
.(5)
Если рассмотреть и решить две другие типовые схемы замещения (включены два ключа сверху и ключ снизу, отключение источника от нагрузки), то можно получить аналогичные (4) и (5) матрицы,
причем матрица A будет незначительно отличаться некоторыми элементами, а вектор (5) будет одинаковым для всех схем замещения.
1
Рис. 4. Схема замещения трехфазного корректора коэффициента мощности при открытых силовых ключах
УТ1,УТ3,УГ 5 или УТ2,УТ4,УТ6
Рассмотрим вектор (5). В элементах матрицы присутствуют выражения, содержащие параметр фазных напряжений трехфазного источника и, и2 и из. Целесообразно для упрощения расчетных
выражений прибегнуть к разложению матрицы на три составные части, содержащие одно из фазных напряжений в своих элементах соответственно.
B = в* + B* + B*;
[(L2 + -[ L3V]U2 -[ L2V ]U3
в* = -[ LV ] Ui ; в2 = [(L + L3)v ]U2 ; вЗ = -[ LV ]U3
0 0 0
Входные напряжения каждой из фаз можно представить как
U = Em sinCt+ф) = Em (sinct- cos^ + cosct- sin^); U2 = Em sin(ct+ ф2) = Em (sin cot- cos ф2 + cos cot- sin ф2); U3 = Em sin(ct+^) = Em (sinct-cos^ + cosct-sin^). Расчет моментов коммутации и перехода корректора из одной схемы замещения в другую происходит в расчете системы разностных функций:
Фм (ZM ) = a2(P2UA Xai(U3-PUc )-Рг I A ) -Usaw;
al(U3 -PUc)-РзIВ)-Usaw;
ai(U3 -PUc)-P3IC) Usaw•
Исследование нелинейной динамики
-РкБ (2кВ ) = а2(Р2иВ () =а2(-Р2иС Здесь Цл, ис - значение фазных напряжений источника в конкретный момент времени коммутации; /л, /в, 1с - значение фазных токов источника в конкретный момент времени коммутации; ис -напряжение конденсатора; ^а^^ - развертывающее пилообразное напряжение; а и Р - коэффициенты пропорционального регулятора и обратной связи соответственно.
Чтобы найти время коммутаций на тактовом интервале работы, нужно знать значения переменных состояния в эти моменты времени. Для расчета переменных состояния в моменты коммутаций верно следующее уравнение:
X (t) = eA(t-to) X о+e
A(t-о) f c>-A(t-to)
Bdr
Зная значения переменных состояния в любой момент времени для каждой из схем замещения и моменты коммутации этих схем, можно составить кусочно-непрерывную математическую модель электромагнитных процессов в корректоре коэффициента мощности на основе инвертора напряжения в базисе коммутационно-разрывных функций.
трехфазного корректора коэффициента мощности проводилось на основе специализированных алгоритмов, представленных в [6].
В ходе расчетов выбран следующий набор параметров: ит=311 В; w=314 Гц; Rн=30 Ом; &в=40 кГц; Яи= Яь2= Яьэ=1 Ом; Ь1= Ь2= Ьз=0,8 мГн; С=25 мкФ; ир=10 В; изн=7 В; а1=50; а2=0,5; р1=0,015; Р2=1; Р3=0,0015.
На рис. 5 представлена двухпарамет-рическая карта динамических режимов, на которой отмечены области существования различных режимов в пространстве двух параметров: сопротивления нагрузки Ян и амплитуды входного напряжения и^.Области существования различных динамических режимов отмечены символами Пгу (/ - да-цикл, характерный для данной области, j - номер области на карте динамических режимов). В частности, область Пду представляет собой первую область существования основного (проектного) режима с частотой н (1-цикл). Область Пх соответствует недетерминированным ре-
о
жимам функционирования преобразователя (т ^ го).
Как видно из рис. 5, при номинальной амплитуде входного напряжения ит =310 Вв системе возможно появление непроектных режимов при падении сопротивления нагрузки Ян ниже 22 Ом. С ростом входного напряжения растет и граница перехода из проектного 1 -циклового режима. В случае низкого сопротивления
нагрузки возможно проявление различных вариаций непроектных режимов. При сопротивлении нагрузки 16 Ом и амплитудном напряжении 310 В корректор будет работать в трехцикловом режиме, то есть частота колебаний будет в три раза выше номинальной, что недопустимо. Поэтому рекомендуется не допускать работу корректора на низкоомную нагрузку.
Рис. 5. Двухпараметрическая карта динамических режимов трехфазного корректора коэффициента мощности
Неправильный подбор параметров для корректора коэффициента мощности может привести к снижению его КПД, выходу из нормального режима работы, перегреву и поломке. Особенно большой угрозе подвержены конденсаторы фильтров в данных устройствах.
Представленная математическая модель электромагнитных процессов трехфазного корректора коэффициента мощности на основе инвертора напряжения позволила выявить возможность возникновения нежелательных режимов при номинальном режиме работы
трехфазного корректора коэффициента мощности. Получена система дифференциальных уравнений, определяющих электромагнитные процессы в трехфазном корректоре коэффициента мощности, приведен пример моделирования с помощью программы Mathworks Matlab, построены бифуркационные диаграммы и проведен анализ влияния сопротивления нагрузки в проектном режиме работы корректора. Даны рекомендации для выбора основных параметров трехфазных корректоров при инженерном проектировании.
Работа выполнена в соответствие с заданиемМинобрнауки № 8.1729.2017/4.6
СПИСОК ЛИТЕРАТУРЫ
1. Giauoris, D. Control of Fast Scale Bifurcations in Power-Factor Correction Converters / D. Giauoris, S. Banejjee, B. Zahavi, V. Pickert // IEEE Transactions on Circuits and Systems - II: Express Briefs. - 2007. - Vol. 54. - № 9. - P. 805-809.
2. Barbosa, P.M. Three-Phase Power Factor Correction Circuits for Low-Cost Distributed Power Systems / P.M. Barbosa // Faculty of the Virginia Polytechnic Institute and State University. - 2002. - P. 245.
3. Prasad, A.R. An active power factor correction technique for three-phase diode rectifiers / A.R. Prasad, P.D. Ziogas, S. Manias // IEEE Trans. Power Electron. - 1991.- Vol. 6. - P. 83-92.
4. Naik, R. Third harmonic modulated power electronics interface with 3-phase utility to provide a regulated DC output and to minimize line-current harmonics /R. Naik, M. Rastogi, N. Mohan // Conf. Rec. IEEE-IAS Annu. Meeting. - 1992. -P. 689-694.
5. Tognoli, M. A DSP based control for a symmetrical three-phase two-switch PFC-power supply for variable output voltage / M. Tognoli, A.C. Rufer // IEEE PESC'96. - 1996. - P. 15881594.
6. Андриянов, А.И. Алгоритмы для бифуркационного анализа обобщенной модели преобразователей постоянного напряжения / А.И. Андриянов, Н.И. Булохов // Справочник. Инженерный журнал. - 2013. - № 10. - С. 30-39.
1. Giauoris, D. Control of Fast Scale Bifurcations in Power-Factor Correction Converters / D. Giauoris, S. Baneijee, B. Zahavi, V. Pickert // IEEE Transactions on Circuits and Systems - II: Express Briefs. - 2007. - Vol. 54. - № 9. - P. 805809.
2. Barbosa, P.M. Three-Phase Power Factor Correction Circuits for Low-Cost Distributed Power Systems / P.M. Barbosa // Faculty of the Virginia Polytechnic Institute and State University. - 2002. - P. 245.
3. Prasad, A.R. An active power factor correction technique for three-phase diode rectifiers / A.R. Prasad, P.D. Ziogas, S. Manias // IEEE Trans. Power Electron. - 1991.- Vol. 6. - P. 83-92.
4. Naik, R. Third harmonic modulated power electronics interface with 3-phase utility to
provide a regulated DC output and to minimize line-current harmonics /R. Naik, M. Rastogi, N. Mohan // Conf. Rec. IEEE-IAS Annu. Meeting. -1992. - P. 689-694.
5. Tognoli, M. A DSP based control for a symmetrical three-phase two-switch PFC-power supply for variable output voltage / M. Tognoli, A.C. Rufer // IEEE PESC'96. - 1996. - P. 15881594.
6. Andriyanov, A.I. Algorithms for bifurcation analysis of generalized model of constant voltage converters / A.I. Andriyanov, N.I. Bulokhov //
Reference Book. Engineering Journal. - 2013. -№ 10. - pp. 30-39.
Статья поступила в редколлегию 31.03.17. Рецензент: д.т.н., профессор Брянского государственного технического университета
Потапов Л.А.
Сведения об авторах:
Бутарев Игорь Юрьевич, науч. сотрудник Брянского государственного технического университета, е-таП: igorbutarev@gmail.com.
Butarev Igor Yurievich, Research Assistant, Bryansk State Technical University, e-mail: igorbu-tarev@gmail.com.
