Динамическая модель системы распределения мощности между компенсированными инверторами тока

Тип работы:
Реферат
Предмет:
Общие и комплексные проблемы естественных и точных наук


Узнать стоимость новой

Детальная информация о работе

Выдержка из работы

УДК 621. 314
ДИНАМИЧЕСКАЯ МОДЕЛЬ СИСТЕМЫ РАСПРЕДЕЛЕНИЯ МОЩНОСТИ МЕЖДУ КОМПЕНСИРОВАННЫМИ ИНВЕРТОРАМИ ТОКА
Астапович Ю. М., Миргородская Е. Е., Митяшин Н. П., Максимова Н. Н.
ФГБОУ ВПО «Саратовский государственный технический университет имени Гагарина Ю. А. «, Саратов,
Россия (410 012, Саратов, ГСП, ул. Политехническая, 77), e-mail: mityashinnp@mail. ru_
Построение систем электропитания на основе группового включения модулей силовой электроники представляет интерес на предприятиях, имеющих нагрузки широкой гаммы мощностей и требующих питания на частотах, отличных от общепромышленной. Преимущества этого метода состоят в возможности получить необходимые значения мощности источника, используя ограниченный набор силовых модулей, удешевлении резервирования, возможности оперативной адаптации преобразователя к величине и характеру нагрузки. Среди задач теории и практики групповых преобразователей одной из основных является обеспечение высокого качества переходных процессов в источнике при изменении нагрузки и напряжения питающей сети. В статье представлена математическая модель двухмостового группового преобразователя частоты на основе однофазных инверторов тока для питания изменяющейся в широких пределах активно-индуктивной нагрузки. Модель предназначена для исследования динамики автоматической системы распределения нагрузки между составляющими инверторными мостами.
Ключевые слова: инвертор, параллельная работа, мощность, распределение, динамика, дискретные динамические модели.
DYNAMIC MODEL OF POWER DISTRIBUTION BETWEEN COMPENSATED CURRENT INVERTER
Astapovich Y.M., Mirgorodskaya E.E., Mityashin N.P., Maksimova N.N.
Saratov State Technical University named after Gagarin Yu. A., Saratov, Russia (410 012, Saratov, GSP, Street
Polytecheskaya, 77), e-mail: mityashinnp @mail. ru_
Construction of power systems based on a group of power electronicsmodules is of interest to enterprises with a wide range of load capacities that require power at frequencies different from the general industrial. The advantages of this method are possible to get the required power source, using a limited set of power modules, reducing the cost of redundancy, quickly adaptationof inverter to the magnitude and nature of the load. Among the problems of the theory and practice of group one of the main transformers is to provide high-quality source of transients with load and supply voltage. The paper presents a mathematical model of a double-bridge inverter group on the basis of single-phase inverter power supply for varying over a wide range of active-inductive load. The model is designed to study the dynamics of the automatic load sharing between the components of the inverter modules.
Keywords: inverter, parallel operation, power distribution, dynamics, discrete dynamical models. Введение
Увеличение мощности в единице оборудования может быть достигнуто за счет параллельной работы нескольких силовых модулей. В частности, этот прием часто применяется в системах электроснабжения на основе устройств силовой электроники. В настоящей работе решается задача построения линеаризованной модели системы распределения мощности между компенсированными инверторами тока, объединенными по выходу. Такая модель необходима при исследовании динамических процессов в системе при изменении нагрузки или напряжения питающей сети. Цель исследования
Разработка методики моделирования динамических режимов систем распределения мощности между объединёнными по выходу силовыми преобразователями.
Выбор и обоснование метода
Известно, что практическое использование параллельного инвертора тока без использования компенсаторов избыточной реактивной мощности конденсаторной батареи ограничено. Мягкая внешняя характеристика не позволяет использовать его для питания сильно меняющейся нагрузки. Для стабилизации внешней характеристики параллельного инвертора применяют его совместную работу с тиристорно-реакторными или диодно-реакторными компенсаторами [2].
На рис. 1 и 2 показаны схемы однофазных инверторов тока с тиристорно-реакторными (ТРК) и диодно-реакторными компенсаторами (ДРК).
Рисунок 1. Однофазный инвертор тока с тиристорно-реакторным (слева) и с диодно-реакторным (справа) компенсаторами
Отличие работы ТРК и ДРК состоит в характере управления током в реакторе компенсатора. В ТРК длительность токопротекания через реактор задается внешней системой автоматического управления, тогда как в ДРК длительность компенсирующего импульса тока реактора определяется углом запирания в [2].
Строго говоря, оба инвертора являются нелинейными объектами, т. е. непосредственное применение методики построения линейной дискретной динамической модели к ним неправомерно. Однако для каждой из них можно предложить прием аппроксимации процесса компенсации, который приводит к линейной модели. Покажем это для инвертора
На интервале симметрии инвертор с ТРК последовательно принимает три состояния, показанные на рис. 2. Соответствующие эпюры токов реакторов компенсатора и напряжения на коммутирующем конденсаторе приведены на рис. 3.
Т6
Ъл
Тн
тока с ТРК.



Рисунок 2. Состояния инвертора с ТРК на интервале симметрии
На первом подынтервале [0, 61] (рис. 3 слева) реализуется состояние инвертора, показанное на рис. 2а. В начале этого подынтервала включились тиристоры Т1 и Т4, что благодаря показанной на рисунке полярности напряжения и привело к выключению ранее проводивших ток тиристоров Т2 и Т4. Тиристор Т6, включенный на предыдущем интервале,
проводит ток^с^ индуктивности Ьс, который имеет отрицательное значение.
Рисунок 3. Эпюры выходного напряжения и тока реактора Ьс инвертора с ТРК на интервале симметрии при, а & lt- п (слева) иа & gt- п (справа)
Поскольку напряжение и, перезаряжаясь током метает знак, то ток?? са, увеличивающийся по абсолютной величине, пока сохраняется отрицательное значение напряжения, начинает уменьшаться по модулю и в момент времени, соответствующий электрическому углу 0, становится равным нулю.
С этого момента наступает второе состояние (второй подынтервал, рис. 2 б). Оно заканчивается, когда на тиристор Т5 подается управляющий импульс и, поскольку
напряжение u в этот момент в 2 положительное, он открывается. Под действием напряжения u в индуктивности L развивается импульс тока '--La (рис. 3, подынтервал [02, л]).
Если импульсы токов индуктивности L и L не перекрываются во времени, то
вместо двух индуктивностей используют однуиндуктивность Lc. При этом ток этой индуктивности состоит из двух разнополярных колоколообразных импульсов, длительность каждого из которых регулируется моментом подачи импульса на соответствующий тиристор ТРК (на рассматриваемом интервале это Т5, момент времени в 2).
Для регулирования напряжения на выходе инвертора необходимо изменить коэффициент передачи по напряжению. При этом, чем меньше длительность и амплитуда импульсов тока iic, тем большее значение имеет этот коэффициент и, наоборот, при их увеличении коэффициент передачи инвертора по напряжению уменьшается. По определению коэффициент передачи инвертора по напряжению равен
= u/ud,
где U — действующее напряжение на выходе инвертора, Ud — постоянное значение напряжения на его выходе.
В предположении о синусоидальности кривой выходного напряжения длительность каждого из разнополярных импульсов тока? Lc равна 20, где
d = p-a + ?.
Здесь, а — электрический угол управления, равный углу задержки относительно начала соответствующего интервала импульса, открывающего очередной тиристор ТРК (в нашем
случае a = 02, рис. 3) —? — угол запирания рабочего тиристора, соответствующий моменту смены знака выходного напряжения на интервале.
Если угол a задается системой управления ТРК и может считаться известным, то угол? является функцией параметров схемы инвертора и нагрузки.
В [4] показано, что при регулировании по цепи питания инвертора и постоянной нагрузке угол? изменится незначительно (2−3%). Это значение угла может быть найдено из следующего нелинейного уравнения
1
со-z ¦ C---(2(p-a+ ?) + sin2(a-?)) = sin^ + cos j- tg?,
2p
где
z

Определив из этого уравнения угол b при данных значениях параметров нагрузки Z cos j, параметров схемы C и Lc, частоты w и угла а, мы можем найти значения угла 8, а,
следовательно, длительность каждого из подынтервалов а, б, в (рис. 2).
Из второго из этих рисунков имеем длительность первого подынтервала
t = q / w = (28 + a — p)/ w-
длительность второго подынтервала
t2 = (q2 -qi)/ w= (p- 28)/ w-
длительность третьего подынтервала
t3 = (p-q2)/w= (p-a)/w.
При a & gt- п последовательность состояний меняется. В начале интервала имеет место состояние б рис. 2, поскольку импульс управления, последовательный на тиристоре Т6 при a & gt- п, задерживается с предыдущего интервала. Если a & lt- п+в, то при в = a-п возникает состояние, а рис. 2.
Если a & gt- п+в, то этого не происходит, поскольку напряжение на конденсаторе в момент подачи импульса на Т6 имеет запирающую полярность. Поэтому в этом случае инвертор функционирует как параллельный инвертор без компенсатора.
Состояние, а имеет место, пока ток через реактор L сохраняет направление,
проводящее для тиристора Т6. В момент прохождения этого тока через нуль инвертор снова переходит в состояние б и находится в нем до конца интервала.
Из рис. 2 длительность первого подынтервала в случае a & gt- п+в равна
t = q / w = (a-p)/ w-
длительность второго подынтервала
t2 = (q2 -qi)/ w= (2b)/ w-
длительность третьего подынтервала
t3 = (p-q2)/w= (2p -a- 2b)/w.
Случай a = п фактически сводит работу инвертора с ТРК к работе инвертора с ДРК, так как ток в индуктивности Lc начинается в начале интервала (см. рис. 2).
Таким образом, при сделанных предположениях, используя знание моментов перехода инвертора из одного состояния в другое, мы можем построить дискретную динамическую модель инвертора с ТРК при его управлении по цепи питания со стороны звена постоянного тока. Именно такой способ управления потоками энергии через инверторы используется в системе распределения мощности между инверторами тока с ТРК, блок-схема которой приведена на рис. 4.
Динамическая модель системы распределения для инверторов с ТРК строится по методике, описанной в работе [5].
ИППТи
Рисунок 4. Блок-схема системы распределения мощности для двух инверторов с ТРК
Отличие рассматриваемой системы (рис. 4) от системы распределения мощности между некомпенсированными инверторами [3] состоит в отсутствии «смешивания» сигналов обратной связи системы стабилизации выходного напряжения и системы распределения мощности. В результате обе автоматические системы функционируют автономно, что обеспечивает развязку переходных процессов в этих системах. Сказанное позволяет строить динамическую модель системы распределения независимо от системы стабилизации, ограничившись параметрическим влиянием угла, а управления цепью ТРК через силовое воздействие импульсов тока индуктивности Ьс на исследуемые процессы.
Сначала рассмотрим случайа & lt- п. Поскольку на интервале симметрии инвертора в данном случае последовательно возникают три состояния силовой схемы, то мы получаем три соответствующие системы дифференциальных уравнений (ДУ).
ЬХ + X = V-
Ьх + Л X = V- Ьх + Л3 х = V.
Вектор переменных состояния имеет вид
х = column 2Л 2
и1 и 2
2
и 1 1
Ьс
Матрица Ь имеет вид
Ь = (Иа^ Ьп, Ье 2, Сх, С2, Ьл, 4 2, C, 1,4).
Сами системы ДУ отличаются только матрицами при векторе состояния которые
Я = Я, =
ующий вид:
! 0 1 0 0 0 0 0 0 & quot- яг 1 0 1 0 0 0 0 0 0
0 Я{ 2 0 1 0 0 0 0 0 0 ЯГ 2 0 1 0 0 0 0 0
-1 0 0 0 1 0 0 0 0 -1 0 0 0 1 0 0 0 0
0 -1 0 0 0 1 0 0 0 0 -1 0 0 0 1 0 0 0
0 0 -1 0 0 1 0 0 — Я = 0 0 -1 0 0 1 0 0
0 0 0 -1 0 2 1 0 0 0 0 0 -1 0 2 1 0 0
0 0 0 0 -1 -1 0 1 1 0 0 0 0 -1 -1 0 1 0
0 0 0 0 0 0 -1 Я 0 0 0 0 0 0 0 -1 Я 0
0 0 0 0 0 0 -1 0 Я _ 0 0 0 0 0 0 0 0 к
Очевидно, что эти отличия определяются наличием или отсутствием на соответствующем подынтервале тока индуктивности Ьс.
Длительности временных подынтервалов, в течение которых справедливы системы уравнений, определены выше. Поэтому эти значения ц, т2, т3 мы считаем заданными. После записи системы ДУ в форме Коши
X = АД + Ь~1Х- 1 = 1,2,3,
где
А, =-Ь-1 ^,
находим значения расчетного вектора Xв граничные моменты указанных подынтервалов:
Х (г1) = еАт1 хп + (Е — вАт1) Я1−1 X- Х (т1 +т2) = е^2X (г1) + (Е- еАт2) Я2X- Х (т1 + т2 + т3) = еАзТзх (т1 + т2) + (Е- еАзТз) Я3−1 X.
г, е*1 еА1Г2 еА*3
Здесь мы, как и ранее, используем матричные экспоненты е, е, е, позволяющие компактно записать результат решения ДУ аналитическим методом [1]. Вектор ¦ представляет собой вектор входных напряжений разомкнутой системы.
Объединяя записанные решения для подынтервалов и используя метод одного интервала, для значения вектора расчетных переменных в начале следующего интервала получаем равенство
Хп+1 = К • + +t3),
где в данном случае
К = diag (1,1,1,1,1,1,-1,-1,-1),
получаем искомую связь между двумя последовательными значениями и, т. е. дискретную динамическую модель разомкнутой системы
4+1 = ВХп + ^
где
В = КеАТз еАъ еАт1- Е = К (ем еМ2 (Е- ем) Я- + е3 (Е- е4'-2)Я2- + (Е- еАъ) Я3−1).
Для получения модели замкнутой системы распределения дополним последние соотношения уравнениями обратных связей и регуляторов
VI = иа0 + К{ (-/2 X +71 Х2) + (-?2 *5 +71 Хб) —
2 = иа0 + К/(72х -71 Х2) + КА (72X -71×6). В матрично-векторной форме уравнение для вектора V теперь может быть записано так:
V = V, + Сх,
где
~иё 0 ¦ '-- К г 72 К7 0 0 — К7 К7 0 0 0
и& lt- 0 К7 — К7 0 0 К7 — К7 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 — С= 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
Здесь ил 0- текущее значение выпрямленного напряжения, подаваемого на вход ИППТ.
Подставив последнее выражение для вектора V в уравнение разомкнутой системы, получим окончательную дискретную динамическую модель замкнутой системы распределения мощности
/V
хп+х = Пх» + ,
где
ВВ = В + ЕС.
В случаеа & gt- п порядок сомножителей в формулах для матриц В и Е модели системы меняется. Эти матрицы в этом случае имеют вид
В = КеАТз еАт2 еАь- Е = К (еАл еАъ (Е — еАь) Я1 -1 + е2'-3 (Е — еАЛ) Я2−1 + (Е — еАл) Я3−1).
Построенная модель системы распределения была использована при исследовании ее динамики и выборе коэффициентов усиления регуляторов.
Заключение
Разработанная модель служит основой для исследования систем заданного автоматического распределения нагрузки между объединенными по выходу силовыми агрегатами. Информация о динамических свойствах системы и о влиянии отдельных параметров на устойчивость и качество регулирования содержится в матрице D модели и прежде всего в ее спектре.
Список литературы
1. Абгарян К. А. Матричное исчисление с приложениями в теории динамических систем. -М.: Физматлит, 1994. — 544 с.
2. Артюхов И. И. Автономные инверторы тока в системах электропитания / И. И. Артюхов, Н. П. Митяшин, В. А. Серветник // Сарат. политехн. ин-т. — Саратов, 1992. — 152 с.
3. Астапович Ю. М. Автоматическое распределение нагрузки между параллельно работающими агрегатами / Ю. М. Астапович, Н. П. Митяшин, Р. А. Билюков, Н. А. Калистратов // Современные проблемы науки и образования. — 2013. — № 1- URL: http: //www. science-education. ru/107−8571 (дата обращения: 12. 03. 2013).
4. Митяшин Н. П. Многоуровневое моделирование электромеханических комплексов / Н. П. Митяшин, О. Н. Рябов, А. А. Смирнова // Электротехнические комплексы и силовая электроника. Анализ, синтез и управление. Межвуз. науч. сб. СГТУ. — 2002. — С. 74−83.
5. Нугаев Э. К. Дискретные динамические модели структурно симметричных инверторов / Э. К. Нугаев, Н. П. Митяшин, Е. Е. Миргородская, А. А. Щербаков // Управление сложными системами: сб. науч. ст. — Саратов: СГТУ, 2009. — С. 3−9.
Рецензенты:
Игнатьев А. А., д.ф. -м.н., профессор, заведующий кафедрой «Автоматизация, управление, мехатроника» Саратовского государственного технического университета имени Гагарина Ю. А., г. Саратов.
Коломейцев В. А., д.т.н., профессор кафедры «Радиоэлектроника и телекоммуникации» Саратовского государственного технического университета имени Гагарина Ю. А., г. Саратов.

ПоказатьСвернуть
Заполнить форму текущей работой